심해 다중 채널 합성 개구 소나(SAS) 마이크로 내비게이션 위상 오차 보정 서보 엔진

1. 물리적 배경: 합성 개구 소나(SAS)와 가간섭성 붕괴

합성 개구 소나(Synthetic Aperture Sonar, SAS)는 잠수정이 이동하며 방출한 연속 음향 핑(Ping)의 반사 신호를 코히런트(Coherent) 적분하여, 물리적 안테나 길이의 100배가 넘는 가상 개구면을 형성함으로써 수백 미터 전방에서도 수 센티미터($< 3\text{ cm}$)급 해저 등방향 분해능을 달성합니다.

그러나 수중 음속 $c \approx 1,500\,\text{m/s}$의 저속 전파 특성으로 인해 핑 간격 동안 플랫폼의 위치 오차 $\Delta r > \lambda / 8 \approx 1.8\,\text{mm}$가 발생하면 위상 상쇄 간섭이 발생하여 영상이 완전히 흐려집니다.

2. 수학적 모델링: DPCA 중복 위상 중심 교차 상관

다중 수신 배열 소나(Multi-Receiver SAS)에서는 연속된 두 핑 사이에서 플랫폼의 전진 속도 $v$와 펄스 반복 주기(PRI)에 의해 공간적으로 일치하는 위상 중심(Phase Center) 쌍이 존재합니다. 이 중복 소나 수신 채널 신호 $s_1(t)$와 $s_2(t)$의 복소 교차 상관을 수행합니다:

피크 상관 시점의 복소 위상각 $\Delta \phi$로부터 파장 대비 시선 방향 변위 편차 $\Delta r$을 직접 유도합니다:

3. 10,000Hz 0-Alloc C# 마이크로 내비게이션 위상 서보 엔진

심해 UUV 임베디드 DSP 및 에지 컴퓨팅 보드에서 핑 레이트 동기화 및 0.1ms 초저지연 위상 회전을 수행하는 C# 무할당 구조체를 구현합니다.

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.OceanSensing
{
    public readonly ref struct SasPingChannelSample
    {
        public readonly double InPhaseReal;
        public readonly double QuadratureImag;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public SasPingChannelSample(double real, double imag)
        {
            InPhaseReal = real;
            QuadratureImag = imag;
        }
    }

    public ref struct SasMicronavPhaseServo
    {
        private readonly double _wavelengthMeters;
        private double _accumulatedPhaseDriftRad;
        private double _crossCorrReal;
        private double _crossCorrImag;

        public SasMicronavPhaseServo(double acousticFreqHz, double soundSpeedMps = 1500.0)
        {
            _wavelengthMeters = soundSpeedMps / acousticFreqHz;
            _accumulatedPhaseDriftRad = 0.0;
            _crossCorrReal = 0.0;
            _crossCorrImag = 0.0;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public void AccumulateCorrelator(in SasPingChannelSample s1, in SasPingChannelSample s2)
        {
            // 복소 켤레 곱: s1 * s2*
            _crossCorrReal += s1.InPhaseReal * s2.InPhaseReal + s1.QuadratureImag * s2.QuadratureImag;
            _crossCorrImag += s1.QuadratureImag * s2.InPhaseReal - s1.InPhaseReal * s2.QuadratureImag;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public void SolvePhaseCorrection(Span correctionOut)
        {
            double phaseError = Math.Atan2(_crossCorrImag, _crossCorrReal);
            double rangeDisplacementMeters = (_wavelengthMeters / (4.0 * Math.PI)) * phaseError;

            _accumulatedPhaseDriftRad += phaseError;

            correctionOut[0] = phaseError;                // 라디안 위상 오차
            correctionOut[1] = rangeDisplacementMeters;   // 미터 단위 시선 궤적 오차
            correctionOut[2] = _accumulatedPhaseDriftRad; // 누적 보정 위상

            // 다음 프레임을 위한 상관기 리셋
            _crossCorrReal = 0.0;
            _crossCorrImag = 0.0;
        }
    }
}

4. 교육 및 국방/해양방재 연계