13.1.1.1 의사거리(Pseudorange) 측정 방정식 유도 및 오차 항 분석

PX4-Autopilot 생태계가 3차원 절대 공간의 위치를 자율적으로 인식하고 궤적을 제어하기 위한 가장 근본적인 첫 시작점은, 기체 보드에 장착된 위성 항법(GNSS) RF 수신기가 각 단일 위성으로부터 발송된 원시 디지털 코드를 복조하여 그 거리를 역산해 내는 계측 과정이다. 이 최초의 스칼라 거리 관측치를 의사거리(Pseudorange) 라 엄밀히 칭한다. 이는 우주 공간의 진정한 기하학적 물리 직선거리(True Geometric Range)가 아니라, 송수신 하드웨어 클록(Clock)의 비동기 편향과 대기권 통과 지연 굴절로 인해 수학적으로 심각하게 오염된 일종의 ’가짜 오차 거리(Pseudo-)’를 내포하고 있음을 뜻한다.

본 절에서는 단일 엣지 수신기가 이 의사거리를 물리적으로 어떤 지배 방정식하에 도출하는지 측정 뼈대 원리를 밝히고, 이 관측치에 끈적하게 융합된 핵심적인 오차 항(Error Terms)들의 스칼라 역학을 학술적으로 분해 분석한다.

1. 시간 지연(Time Delay: \Delta t) 계측 역학

지구 상공 궤도를 도는 GPS 위성은 주기적으로 자신의 궤도 좌표 변수(Ephemeris)와 현재 타임스탬프의 정밀 발사 시간(Time of Transmission, t_{Tx})을 PRN(Pseudo-Random Noise) 코드 형태로 디지털 변조하여 광속의 전파로 우주 공간에 하향 흩뿌린다. 지상관제소나 공중에 떠 있는 로버(Rover) 기체 수신기는 안테나 렌즈에 맺힌 이 백색 잡음과 같은 신호를 칩셋에서 상관 획득(Acquisition/Correlation)한 후, 수신기 자체의 내부 수정 발진기 시간(Time of Reception, t_{Rx})과 들어온 PRN 코드를 정밀 수학 위상 상관 연산하여 둘 간의 기계적 도달 시간차 통계 \Delta t를 구하게 된다.

\Delta t = t_{Rx} - t_{Tx}

이 측정 도출된 도달 시간차 스칼라에 진공 광속 상수 c (정확히 2.99792458 \times 10^8 m/s)를 비례하여 곱해주면 가장 원초적인 1차원 의사거리 관측치 스파이크 P_{raw}를 얻게 된다.

P_{raw} = c \cdot \Delta t = c \cdot (t_{Rx} - t_{Tx})

2. 측정 방정식(Measurement Equation) 유도 및 에러 분리

하지만 우주의 위성 하드웨어 시간 t_{Tx}와 지상 수신기 하드웨어 시간 t_{Rx}는 결코 절대 우주 진리 시간(True GPS System Time, t_{true})과 완벽히 오차 없이 동기화되어 있지 않다.

• 수신기 하드웨어 수신 시간: t_{Rx} = t_{true, Rx} + dt_u (수신기 저가 쿼츠 시계의 하드웨어 오차 편향 dt_u 수반 지연)
• 위성 하드웨어 송신 시간: t_{Tx} = t_{true, Tx} + dt^s (위성 탑재 원자 시계의 상대성 오차 편향 dt^s 수반 지연)

이를 앞선 물리적 시간차 방정식에 대입하면 수학적으로 다음과 같이 행렬화 전개된다.
P = c \cdot (t_{true, Rx} - t_{true, Tx}) + c \cdot (dt_u - dt^s)

여기서 진리 도달 시간차 (t_{true, Rx} - t_{true, Tx})에 절대 광속 c를 곱한 값이 바로 완전한 우주 진공 상태를 가정할 경우의 기하학적 실제 통달 거리(Geometric True Euclidean Range, \rho) 역학이 된다.
그러나 실제 전파는 뚫기 힘든 대류권과 전리층이라는 불규칙한 지구 대기 플라즈마/수증기 매질을 겪으면서 굴절률(Refractive Index) 변화에 따른 엄청난 추가적 파동 통과 시간 지연(d_{ion}, d_{trop})을 가혹하게 마주하게 된다. 여기에 최종적으로 도심 지형지물 렌즈 반사에 의한 다중 경로 왜곡 스펙트럼과 하드웨어 내부 아날로그 열잡음(\epsilon_P)까지 마지막으로 합산 적층하면, 마침내 완전한 측위 공학적 의사거리 지배 방정식이 참혹하게 도출된다.

P = \rho + c \cdot (dt_u - dt^s) + d_{ion} + d_{trop} + \epsilon_P

3. 의사거리 핵심 오차 항(Error Terms) 스케일 분석

위 절대 지배 방정식에서 실제 거리 참값 \rho를 오염 교란시키는 각 항의 물리적 크기 스케일과 제어 상쇄 가능성을 대수학적으로 분해하면 다음과 같다.

1. 위성 시계 오차 (c \cdot dt^s): 물리적 오차 스케일 \sim 2m
   수만 달러의 위성 내부에는 세슘(Cs)이나 루비듐(Rb) 기반의 초정밀 원자 시계가 펌핑 탑재되어 있지만, 아인슈타인의 특수/일반 상대성 역학 이론(초고속 궤도 이동 속도로 인한 시간 지연 및 지구 중력장 왜곡 차이)으로 인해 매일 불가피하게 시간이 앞서가거나 느려진다. 이 편향은 미국 우주군 지상 통제국(Control Segment)에서 매순간 칼만 추정 계산하여 항법 메시지(Navigation Message) 데이터 블록의 다항식 계수(Polynomial Coefficients: a_{f0}, a_{f1}, a_{f2})로 브로드캐스트하므로, 수신기 펌웨어 DSP 단에서 알고리즘적으로 사실상 수 cm 단위까지 상당 부분 완벽히 차감 보상 제거가 가능하다.
2. 수신기 시계 오차 (c \cdot dt_u): 물리적 오차 스케일 수십 \sim 수백 m
   드론의 칩 부피와 단가 발열 제한 문제로 인해 PX4 Pixhawk 메인보드 칩셋 포트들은 일반적인 저가 크리스탈 온도 보상 수정 발진기(TCXO)를 강제 사용한다. 이는 위성 원자 시계 스펙 대비 극도로 부정확하여 가장 거대하고 제어할 수 없는 미지수 오차 폭탄 덩어리를 만들어낸다. 엔지니어링 적인 해결책으로 3차원 궤도 위치(x, y, z) 3개의 변수에 수신기 시계 오차 변수(dt_u) 단 1개를 추가 미지수로 동일 상정하고, 최소 4대 이상의 얽힌 위성 관측 지배 방정식을 일시 연립 행렬망으로 뚫어내어 이 항을 수학적으로 계산 유추 소거해 낸다.
3. 전리층 지연 (d_{ion}): 물리적 오차 스케일 2m \sim 20m 극심한 돌발 변동
   강력한 태양풍(Solar Wind) 이온화 작용으로 인해 고고도 상공(50~1000km) 전리층에 형성된 자유 전자(Free Electrons) 플라즈마 밀접 구역 렌즈이다. RF 주파수가 이를 통과하면 군속도(Group Velocity) 지연 역학 현상에 의해 관측 거리가 늘어나는 착시 지연이 터진다. 주파수 대역의 제곱에 정확히 반비례(1/f^2)하여 지연이 발생하는 수학적 분산 매질(Dispersive Medium) 성향을 확고히 띠므로, L1/L2 이중 주파수 무선기(Dual-frequency Receiver) 하드웨어를 사용할 경우 수학적 역산 조합(Ionosphere-Free Linear Combination)으로 노이즈 파편을 사실상 99% 즉각 소멸시킬 수 있다.
4. 대류권 지연 (d_{trop}): 물리적 오차 스케일 2m \sim 3m
   지표면 근처(0~10km)의 수증기 습윤 밀도와 건조 대기압 공기층을 통과하며 전파가 물리적으로 꺾이고 느려진다. 이 가스 매질은 주파수 대역에 차별 없이 비분산적(Non-dispersive) 굴절 특성을 강하게 띠므로 L1/L2 다중 주파수 보정 스킬로도 수학적 해결이 불가능하다. 대신 하드웨어 내부에서 Saastamoinen과 같은 기후 대기압 역학 모델링 상미분 수식을 통해 머리 위 수직 천정(Zenith) 방향의 지연 상수 스칼라를 예측하고, 각 위성 신호의 수신 앙각(Elevation Angle) 노드에 따른 삼각 함수 매핑 함수(Mapping Function)를 곱해 확률적으로 추정 방어 감소시킨다.
5. 다중경로 및 파동 잔류 노이즈 (\epsilon_P): 물리적 오차 스케일 0.5m \sim 2m
   도심 협곡이나 지표 계곡 등에서 난반사되어 들어온 잔상 시그널 지연파가 원본 직접 신호와 섞이면서 위성 칩 코드 신호 상관 연산 피크(Correlation Peak) 정점을 뭉뚝하게 찌그러뜨리는 지독한 다중경로(Multipath) 에러와 수신 칩셋 하드웨어 노이즈이다. 이는 일정한 모델링으로 파훼하기 가장 까다로운 비선형 랜덤 덩어리이며, RF 노이즈 초크링(Choke Ring) 안테나 접시 설계 등 철저한 하드웨어적인 EMI 디퓨저 필터링 물리 방어막으로 극한 억제해야만 한다.

위와 같이 의사거리 P 관측치에는 수많은 스칼라 오차가 덕지덕지 수렴 및 융합되어 거대하게 층을 이룬다. 아무리 고도화된 PX4 소프트웨어 EKF 추정을 수천 번 루프 수행해도 이 원시 스칼라 의사거리 신호만으로는 결코 1~2m 오차 벽의 한계를 찢고 넘어설 수 없다. 그렇기에 PX4 제어 엔지니어링이 궁극적으로 요구하는 1cm 수준의 완벽에 가까운 자율 제어를 달성하기 위해서는, 필연적으로 이 코드 노이즈 덩어리를 근원 방어 폐기하고 반송파 위상(Carrier Phase) 의 순수 이중 차분 역학인 하드코어 RTK 알고리즘 세계로 도약해야만 하는 당위성이 여기에 존재한다.