13.1.1.1.1 의사거리 기본 방정식: P = \rho + c(dt_u - dt_s) + d_{ion} + d_{trop} + \epsilon_P

PX4-Autopilot의 로컬 GPS 드라이버(px4_gps 모듈 데몬)가 외부 위성으로부터 수신된 직렬 텔레메트리 스트림データを 추적 해석할 때, 가장 최초로 도출해 내는 수학적 대수학의 뼈대가 바로 수신기 엣지 안테나와 단일 위성 간의 거리를 산출하는 의사거리(Pseudorange) 기본 지배 방정식이다. 이 방정식은 단순히 목적지까지의 기하학적 유클리드 공간 거리를 넘어서, 우주와 대기권을 무자비하게 관통하며 누적 적층되는 모든 물리적 오차 스칼라들의 거대한 집합체를 의미한다.

본 절에서는 의사거리 지배 방정식 자체의 수학적 유도 원리와 생성 궤적을 학술적으로 해부하고, 펌웨어가 내부 EKF2 필터를 통해 처리해 내야 할 각 독립 오차 항들의 공학적 의미를 명확히 정립 분석한다.

1. 관측 시간차 지연(Time Delay) 모델링과 기하학적 참거리(\rho)

일반적인 위성 항법 측위(Positioning)의 대원칙은 ’시간이 곧 물리적 거리이다’라는 광학적 공리에서 출발한다. 상공의 위성 s가 자신의 신호 발생 시점(Time of Emission)과 궤도 데이터를 디지털 구형파로 싣고 지상을 향해 전파를 송출하면, 로버 수신기 u는 도달된 신호 패턴 파형을 자신의 내부 RF 클록(Time of Reception)과 교차 상관 연산(Correlation)하여 총 도달 시간차 \Delta t를 고도의 정밀도로 계측하게 된다.

만약 진공 상태 우주에서 완벽한 신의 시계, 즉 참된 위성 항법 시스템 시간(True GPS Time, t_{true})을 위성과 수신기가 단 1나노초의 오차도 없이 완전히 동기화 공유하고 있다고 가정한다면, 두 안테나 사이를 가로지르는 순수 기하학적 3차원 직선거리(True Geometric Range, \rho)는 광속 상수 c를 스칼라배 곱하여 다음과 같이 완벽하게 산출 증명된다.

\rho = c \cdot (t_{true, Rx} - t_{true, Tx})

이 기하학적 거리 \rho는 위성의 3차원 좌표 성분 (X^s, Y^s, Z^s)과 에이전트 수신기의 미지 좌표 (X_u, Y_u, Z_u)를 연결하는 공간 유클리드(Euclidean) 거리 공식의 핵심으로, 수신기 상태 벡터 연쇄 방정식이 찾아내야 할 최종 타겟 참값 우물이다.
\rho = \sqrt{(X^s - X_u)^2 + (Y^s - Y_u)^2 + (Z^s - Z_u)^2}

2. 하드웨어 시계 편향 역학(Clock Bias): (dt_u - dt_s)

그러나 가혹한 물리 엔지니어링 환경의 실제 하드웨어 회로 상에서 완벽한 시스템 시간 동기화란 절대 성립하지 않는다.

• 위성의 공간 오차 (dt_s): 수천만 달러의 위성 인프라는 약 10^{-13} 초 단위 분해능 안정성을 지닌 세슘(Cs) 원자 시계를 탑재함에도 불구하고, 아인슈타인 일반 상대성 체계(지구 중력장 장력 왜곡 차이)와 특수 상대성 역학(고속 궤조 비행체)의 지배에 의해 매일 절대 시간 t_{true}와 미세하게 틀어지는 오차 궤적 편향인 dt_s를 무조건적으로 잉태한다.
• 수신기의 하드웨어 오차 (dt_u): PX4 펌웨어가 얹어지는 무인기 메인보드(Pixhawk, Cube 등) 칩셋 내부의 값싼 석영 크리스탈 수정 진동자(TCXO)는 절대 우주 시간과 비교 불가한 엄청난 시간 괴리를 지속적으로 뿜어내며 오차 편향 dt_u를 발생 증폭시킨다.

• 실제 송신 시간: t_{Tx} = t_{true, Tx} + dt^s
• 실제 수신 시간: t_{Rx} = t_{true, Rx} + dt_u

결과적으로 수신기 디지털 보드가 기계적으로 측정한 시간차 \Delta t를 절대 거리 역산에 대입하면 양측의 가짜 시계 오차가 곱해진 기형적인 상태가 수학적으로 도출된다.
P_{raw} = c \cdot (t_{Rx} - t_{Tx}) = c \cdot (t_{true, Rx} - t_{true, Tx}) + c \cdot dt_u - c \cdot dt_s
\Rightarrow P_{raw} = \rho + c(dt_u - dt_s)

여기서 곱해지는 광속 상수 c의 무자비한 거대함(약 2.997 \times 10^8 m/s) 때문에, 임베디드 수신기 칩 단에서 단 1\mu s(마이크로초) 단위의 비동기 처리 지연 버퍼만 발생해도 기체의 관측 위치 좌표는 순식간에 약 300m라는 어마어마한 수치로 터져서 붕괴되어 버린다.

3. 대기 공간의 굴절 오차 패널티: d_{ion} + d_{trop}

위성 궤도로부터 쏟아진 전파가 우주의 고진공 레이어를 통과해 지구의 두꺼운 대기 분자 질량 속으로 돌입 강하하는 순간, 매질의 굴절률(Refractive Index) 변화에 따른 전파속도(Wave Velocity) 지연 오차가 치명적으로 터져 나온다.

• 전리층 굴절 지연 스칼라 (d_{ion}): 지상 50 \sim 1,000 km 고도 상공에 걸쳐서 태양의 극한 복사 에너지로 인해 파생된 플라즈마(가스 이온 및 자유 전자)가 군속도 전파 펄스의 파동 발목을 잡아 강력하게 지연시킨다. 총 전자 밀도(TEC: Total Electron Content)와 주파수의 대역 밴드 제곱에 철저하게 비례하는 분산 매질(Dispersive Medium) 역학 특성으로, 태양풍 요동 시 최대 20m \sim 30m 이상의 변동성 오류 덩어리를 여과 없이 선사한다.
• 대류권 굴절 지연 스칼라 (d_{trop}): 지상 0 \sim 50 km 권역 대기권 내의 순수 건조 대기(질소/산소 기체 분자) 압력 밀도와 스펙트럼 기재 기상의 수증기(습윤 밀도)에 의한 물리적 저항 인자이다. 주파수 파장의 대역에 일절 관계없이 모든 종류의 위성 관측 데이터에 똑같이 획일적인 2 \sim 3m 안팎의 비분산적(Non-dispersive) 굴절 패널티를 중첩으로 더한다.

4. 다중 경로 에코 및 잔여 스펙트럼 노이즈: \epsilon_P

위와 같은 절대 우주 대수학적 거대한 방정식 항들 외에도, 의사거리 최후단에는 수리적인 공식 모델링으로 잡아내기 극히 까다로운 종단 노이즈 집합체 \epsilon_P가 영원한 꼬리표처럼 뒤틀려 합산된다.
도심 협곡 비행 중 콘크리트 글래스 빌딩 벽과 산악 계곡에 부딪혀 난반사된 지연 시그널 다발들(Multipath Error)이 직접파 신호 동조 피크와 섞여 들어가 뭉개지거나, 안테나 모듈 수신기 LNA(Low Noise Amplifier) 회로 보드 내부 소실 소자의 아날로그 열화학적 열잡음(Receiver White Noise) 파편들이 이 항에 귀속되어 통상 0.5m \sim 2m 가량의 변동 오차율을 지속적으로 토해낸다.

5. 지배 방정식의 최종 결합과 독립 제어의 운명적 한계

결론적으로 이상의 기하학적 참값 스칼라, 시계 왜곡 편향 스칼라, 대기권 굴절 패널티 집합 변수, 그리고 통계적 확률 잔류 에러 뭉치를 단일 수학 선형식으로 모두 누적 합산하면, 마침내 수신기가 해석해 내는 완전한 형태의 의사거리 지배 방정식 정체가 명확한 수식으로 도출된다.

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

그러나 아무리 PX4 펌웨어의 수학적 두뇌인 EKF2 필터 다이내믹스가 최첨단 칼만 상태 연산 행렬을 초당 수백 번 가혹하게 가동 증폭하여 텐서 오차를 역산해 낸다 하여도, 이 C/A 코드 기반의 단방향 원거리 의사거리 P 방정식 하나에만 매몰된 단일 수신기 상태로는 결코 무거운 저 노이즈 굴절항(d_{ion}, d_{trop})의 잔재 뿌리를 100% 온전히 증발 소거할 수 없다. 이를 기반으로 센티미터 수준의 비행 정밀 성능을 뽑아내는 것은 물리 공학적으로 명백한 어불성설이다.

따라서 비행 궤도의 모든 오차 단상을 근원적으로 궤멸시키고 PX4-Autopilot의 완벽 자율화를 성취하기 위해, 무인항공 시스템은 필연적이고 운명적으로 지상 베이스 스테이션과의 반송파 위상(Carrier Phase) 이중 차분 행렬 구조를 투입 융합하는 차세대 RTK 측위 네트워크 기술 플랫폼으로 그 설계적 도약을 강제 받게 되는 것이다.