13.1.1 위성 항법 시스템(GNSS)의 측위 오차 모델링 및 지배 방정식

무인항공기(UAV)가 목적지를 향해 자율 비행을 수행할 때, 가장 근원적이고 치명적인 피드백 시스템 변수는 현재 자신의 절대 3차원 공간 좌표이다. 이 공간 좌표는 지구 궤도를 도는 위성 항법 시스템(GNSS: Global Navigation Satellite System)의 RF 신호를 수신하여 삼각 측량 기반으로 계산된다. 그러나 위성에서 송출된 전파 신호가 우주 진공을 가로질러 대기권을 통과하고, 최종적으로 지표면 근처에 위치한 PX4 비행 기체의 GPS 안테나에 도달하기까지는 수많은 물리적 굴절과 지연 법칙의 가혹한 지배를 받는다.

본 절에서는 단일 GNSS 수신기가 겪게 되는 측위 오차의 근본적인 물리/수학적 모델링을 해부하고, 이를 설명하는 핵심 지배 방정식(Governing Equation)을 공학적으로 다각도로 전개한다.

1. 관측 가능량: 의사거리(Pseudorange)와 반송파 위상(Carrier Phase)

수신기가 하늘의 독립된 위성 신호로부터 획득할 수 있는 원시(Raw) 관측 데이터는 역학적으로 크게 두 가지 물리량으로 분류된다.

의사거리(Pseudorange, $P$ ): 위성에서 전송된 디지털 시간 정보(C/A 코드 또는 P 코드)가 수신기에 도달하는 데 걸린 시간차( $\Delta t$ )에 빛의 물리적 속력 상수( $c$ )를 곱하여 얻어낸 대략적인 거리이다. 신호 파장 주기(L1 C/A 코드 기준 약 300m)가 길어 초기 절대 거리 해석이 매우 쉽지만, 미터(m) 단위의 본질적인 측정 열잡음(Noise) 역학을 수반하게 된다.
반송파 위상(Carrier Phase, $\Phi$ ): 데이터를 변조하여 실어 나르는 초고주파 라디오 반송파파(L1, L2 대역 등) 파동 그 자체의 미세한 사이클 단위 위상 변화를 계측하는 방식이다. 파장이 약 19cm 내외로 매우 짧아 극히 초정밀한 관측이 가능해지지만, 최초 신호 획득 시 신호가 날아온 파장의 전체 누적 반복 횟수를 절대 알 수 없는 치명적인 모호정수(Ambiguity) 맹점 문제를 무조건 해결해야 하는 극악의 비선형 난이도를 띤다.

이 두 가지 독립적인 관측량은 근본적으로 동일한 절대 기하학적 거리를 지향하지만, 우주와 대기권의 외부 환경 오차를 받아들이고 왜곡되는 매커니즘이 수식적으로 미세하게 상이하다.

2. 의사거리(Pseudorange)의 지배 방정식 모델링

위성 $s$ 와 수신기 $u$ 사이의 가장 원초적이고 거시적인 의사거리 관측치( $P_u^s$ )를 지배하는 수학적 방정식은 대수학적으로 아래와 같이 전개 증명된다.

$P_u^s = \rho_u^s + c \cdot (dt_u - dt^s) + d_{ion,u}^s + d_{trop,u}^s + \epsilon_{P_u^s}$

수식에 포진한 각 에러 항의 역학적 원천은 다음과 같다.

기하학적 실제 통달 거리 ( $\rho_u^s$ ): 위성의 3차원 절대 궤도 좌표 $(X^s, Y^s, Z^s)$ 와 엣지 에이전트 수신기의 실제 좌표 $(X_u, Y_u, Z_u)$ 간의 완전 무결한 유클리드(Euclidean) 순수 직선거리이다. 수많은 EKF 필터 추정 끝에 우리가 오차를 모두 걷어내고 기체 제어기(Position Controller)에 전달해야 할 유일한 타겟 참값 변수이다.
$\rho_u^s = \sqrt{(X^s - X_u)^2 + (Y^s - Y_u)^2 + (Z^s - Z_u)^2}$
시계 오차 편향 ( $c \cdot (dt_u - dt^s)$ ): 위성 내부의 원자 시계 오차 편향( $dt^s$ )과, 상대적으로 저비용인 쿼츠(Quartz) 진동자 발진기를 사용하는 기체 하드웨어 클록의 오차 편향( $dt_u$ )에 광속( $c$ ) 스칼라가 배수 연산되면서 발생하는 가짜 거리 증폭이다. 방대한 광속 탓에 단 $1\mu s$ (마이크로초)의 수신기 칩셋 시계 오차만 발생해도, 공간 위치 오차는 즉렬적으로 약 300m나 무자비하게 틀어지게 된다.
전리층 지연 ( $d_{ion,u}^s$ ): 지상 50~1000km 고도의 자유 전자(Free Electrons) 플라즈마 밀도 구간을 RF 전파가 통과하며 겪는 굴절 및 물리적 지연(Delay)이다. 태양 흑점 폭발 등 우주 스펙트럼 날씨에 가장 극심하게 요동치며, GNSS 전체 오차의 가장 압도적인 지분을 차지한다.
대류권 지연 ( $d_{trop,u}^s$ ): 지표면 최하단 근처(고도 0~50km)의 건조 공기 굴절률과 습윤(수증기) 밀집도에 의한 물리적 전파 속도 저하 스칼라이다. 기체의 운항 지역 기상 상태(로컬 온도, 기압, 습도) 빈도에 전면적으로 영향을 받게 된다.
잔여 노이즈 집합체 ( $\epsilon_{P_u^s}$ ): 도시 협곡 빌딩이나 산출물 반사각에 의한 다중 경로(Multipath) 중첩 왜곡 스펙트럼과 수신기 복조기 내부 아날로그 열잡음이 유발하는 확률론적 백색 잡음(White Noise) 파편들이다.

3. 반송파 위상(Carrier Phase)의 지배 방정식 모델링

센티미터(cm) 단위의 RTK 자율 비행을 제어하기 위해 펌웨어가 궁극적으로 최우선 추적 해석해야 할 반송파 위상 관측치( $\Phi_u^s$ ) 방정식은 파장 단위(관측 거리 단위 $m$ )로 정규화 스케일링할 경우 다음과 같은 지배 방정식을 갖는다.

$\Phi_u^s = \rho_u^s + c \cdot (dt_u - dt^s) - d_{ion,u}^s + d_{trop,u}^s + \lambda N_u^s + \epsilon_{\Phi_u^s}$

의사거리 역학 모델링과 매우 흡사한 뼈대 구조를 공유하나, 파동 역학적으로 완전히 거스를 수 없는 결정적 변곡점 차이가 존재한다.

전리층 위상 진전 ( $-d_{ion,u}^s$ ): 군속도(Group Velocity)의 제어를 받는 앞선 디지털 코드 신호 펄스는 대기 마찰처럼 체류 지연을 겪지만, 파동 자체의 위상속도(Phase Velocity) 지배를 받는 전자기 반송파 물결은 전리층 플라즈마 분산 단일 매질을 통과할 때 수학적으로 오히려 진공 광속을 역으로 초과하는 ‘가속 진전(Phase Advance)’ 착시 측정 현상이 도출된다. 따라서 전리층 오차 항의 측정 부호가 기이하게도 음수( $-$ )로 반전되어 기록된다는 점이 지배 방정식의 가장 큰 수리적 특징이다.
모호정수 폭탄의 등장 ( $\lambda N_u^s$ ): 반송파 고유의 파장 길이( $\lambda$ , L1의 경우 약 0.19m)와, 수신기 전원이 켜지고 최초로 해당 위성을 락킹(Tracking)할 당시부터 우주를 가로지른 파장의 전체 누적 정수 폭 궤도 개수인 사이클 미지수( $N_u^s$ )의 곱이다. 이 수치는 단단한 바위처럼, 일단 추적이 시작되어 시간이 지나 파동 락킹이 풀리지(Cycle Slip) 않는 이상 그 항체 통신 세션 수명 안에서는 절대 불변하는 고정된 영구 상수(Constant Integer) 로 존재한다.

4. EKF2 필터 통합 통찰과 오차 소거 철학

위에서 전개된 지배 방정식의 물리적 오차 항들( $d_{ion}$ , $d_{trop}$ , $dt_s$ )은 아무리 하이엔드 CPU를 달고 연산하여도, 하늘을 나는 ‘단일 기체 수신기’ 1대 하나만으로는 완벽히 모델링 예측하거나 수치적으로 씻어내는 것이 수학적으로 구조상 불가능하다. Klobuchar 대기 모델이나 SBAS(위성 기반 오차 광역 보정 시스템)를 통해 통계학적 일부분의 추정 연산 감쇄는 가능하지만, 이는 여전히 미터 급의 운명적 한계를 지닌다.

비행 기체가 치명적인 1~2cm 급의 RTK 체계를 오점 없이 완성하기 위해서는, 결국 저 복잡한 오차 모델 지배 방정식을 무조건 두 개의 물리적 수신기(지상 베이스와 비행 로버)에서 무동기화로 동시에 차감(Difference) 연산시켜, 끈적한 공간 대기 오차와 위성/하드웨어 시계 오차 편향을 수식 양변에서 한 번에 완전히 궤멸 소거(Cancellation)해 버려야만 한다. 이렇게 극한으로 정화된 상대 이중 차분 방정식을 쥐고서야, 비로소 PX4 내부 RTK 데몬과 EKF2 아키텍처는 노이즈가 탈색 증발된 깨끗한 상태 벡터 융합(State Fusion) 업데이트를 수행하며, 무인지상관제시스템(GCS)이 열망하는 오차 제로의 자율 비행 위치 궤도 통제권을 완전하게 획득하게 되는 것이다.