부록: 기초 연습문제 및 해설
문제 1: 위성 위치 벡터 표현
위성항법시스템에서 위성의 3차원 위치를 정확히 기술하기 위해서는 지구중심좌표계(ECEF, Earth-Centered Earth-Fixed)에서의 위성 위치를 나타내야 한다. 일반적으로 위성의 위치 벡터를 $\mathbf{r}_{\text{sat}}$라고 할 때, 이는 다음과 같이 정의할 수 있다.
문제 요구 사항
위성의 궤도 요소(예: 궤도 경사각, 승교점 적경, 근지점 편각 등)를 알고 있을 때, 특정 시각 $t$에서의 $\mathbf{r}_{\text{sat}}$를 유도하는 공식을 간단히 유도하라.
지구 자전 효과를 고려했을 때, 지구중심좌표계에서 위성 궤도 위치를 어떻게 보정해야 하는지 설명하라.
해설 개요
$\mathbf{r}_{\text{sat}}$ 벡터 구성 요소의 해석
지구 회전에 따른 보정 계수 설명
(자세한 중간 계산 과정 및 수식 유도는 다음 섹션에서 단계별로 안내)
문제 2: 측위 방정식 (Pseudo-Range 식)
$N$개의 위성을 사용하여 지상 수신기의 위치를 구할 때, 각 위성으로부터의 관측치(유사거리, pseudo-range)를 $\rho_i$라 하고, 수신기의 위치 벡터를 $\mathbf{x} = (x, y, z)$라고 하자. 이때, 각 위성 위치 $\mathbf{r}_{i} = (x_i, y_i, z_i)$와의 거리 관계식은 다음과 같이 표현된다.
여기서 $c$는 빛의 속도, $d t$는 수신기 시계 오차이다.
문제 요구 사항
위의 측위 방정식을 벡터 형태로 간단히 표현하라.
측위 방정식에서 시계 오차($d t$)가 미치는 영향과 그 추정 방법에 대해 서술하라.
$N$개의 측위 방정식을 행렬 형태로 나타낼 때, 어떻게 선형화하여 최소자승법(LSQ)으로 풀 수 있는지 대략 설명하라.
해설 개요
$\rho_i$에 대한 오차 요소 분해
오차 항 추정을 위한 행렬 식 구성
최소자승법을 활용한 해 구하기 전개
(해설 전개는 뒤에서 단계별로 연습문제 풀이와 함께 상세 제시)
문제 3: 측위 정확도 분석
지상 수신기의 추정 위치를 $\hat{\mathbf{x}}$라고 할 때, 위성기반 측위에서는 위성 배열(위성 기하도)과 측정 잡음 특성에 따라 추정 위치의 오차 공분산이 달라진다. 일반적으로 다음과 같은 공분산 행렬을 사용하여 표현한다.
문제 요구 사항
위성 기하배열에 따른 공분산 행렬의 변화 요인을 설명하라.
$\mathbf{P}$ 행렬을 대각화(diagonalization)하거나 주성분분석(PCA) 관점에서 해석하면 어떤 정보를 얻을 수 있는지 기술하라.
수신기 안테나 높이 오차, 위성 궤도 오차 등이 $\mathbf{P}$에 미치는 정성적 영향(커짐 또는 작아짐 등)을 설명하라.
해설 개요
공분산 행렬이 가지는 기하학적 의미
위성 기하도(DOP)와 공분산의 연관성
안테나 높이 및 위성 궤도 변동 요소 처리 아이디어
문제 4: 동적 모델 기반 항법 해석
위성항법시스템에서 지상 수신기의 이동 궤적(예: 차량, 항공기 등)을 추정할 때, 단순하게 시각별 독립 측위 결과를 연결하는 방법과 달리 동적 모델(운동학적/동역학적 관계식)을 고려할 수 있다. 이를 통해 시간에 따른 수신기 궤적을 더욱 정확하고 안정적으로 추정하는 방안이 제시된다.
문제 요구 사항
운동학적 모델(예: 등가속도, 등의속도 등)을 가정했을 때, 측위 결과에 어떻게 반영되는지 수식으로 표현하라.
동역학적 모델(예: 차량에 대한 타이어 회전 모델, 항공기에 대한 6자유도 방정식 등)을 고려할 때, 측위 결과 오차가 어떻게 감소되는지 정성적으로 설명하라.
$t_{k}$ 시점에서 $\mathbf{x}(t_{k})$의 추정값을 가지고 $t_{k+1}$ 시점의 $\mathbf{x}(t_{k+1})$를 예측(Prediction)하고, 이후 GNSS 관측으로 보정(Update)하는 과정을 선형 칼만필터(Kalman Filter) 관점에서 간략하게 설명하라.
해설 개요
동적 모델에 따른 상태(위치, 속도, 가속도 등) 전이 방정식 기술
GNSS 관측식(측위 방정식)과의 결합을 통한 추정 오차 감소 과정
선형 칼만필터 필수 변수가 되는 예측 공분산, 잔차 공분산 등 행렬적 표현
문제 5: GNSS 기반 시간 동기화
위성항법시스템의 또 다른 중요한 활용 분야는 시각 동기화이다. 위성 시각을 기준으로 지상국의 시계를 맞추는 것이 필요한 상황(예: 통신 시스템 동기화, 전력망 동기화 등)에서, GNSS 위성 신호가 제공하는 정밀 시각 정보를 활용할 수 있다.
문제 요구 사항
시각 동기화를 위한 기본 식을 유도하라. 위성에서 전송되는 시각 정보와 수신기 로컬 시계의 관계식을 간단히 수식으로 나타내라.
지상 수신기의 시계 보정량($d t$)이 $\alpha(t)$라는 함수적 형태를 가진다고 할 때, 이를 연속적으로 추적하여 추정하기 위한 방안(예: 2차 이상 다항 보정, 또는 PLL 기반 보정 등)을 개략적으로 설명하라.
위성 시각 오차와 지상 시각 오차(로컬 시계 드리프트 등)가 결합될 경우, 동기화가 얼마나 민감해지는지 정성적으로 논하라.
해설 개요
위성항법시스템의 시각 정보: GPS Time, GLONASS Time, Galileo Time 등의 차이점
지상 시계 드리프트가 시스템 전체에 미치는 영향
보정 알고리즘의 안정성(짧은 구간 vs. 긴 구간 관점)
문제 6: 항공기용 SBAS 및 GBAS 개념
항공기의 항법 안정성을 높이기 위해 보정신호를 제공하는 SBAS(Satellite Based Augmentation System)와 GBAS(Ground Based Augmentation System)가 있다. SBAS는 광역 보정(WAAS, EGNOS 등)을 제공하고, GBAS는 공항 근방에서 높은 정확도를 확보하기 위한 국지 보정(Local Differential GNSS)을 제공한다.
문제 요구 사항
SBAS가 제공하는 광역 보정 정보(위성 궤도 보정, 위성 시계 보정, 전리층 보정 등)를 간단한 수식 형태로 나타내고, 수신기 측위 식에 어떻게 반영되는지 설명하라.
GBAS 방식에서 사용되는 국지 기준국의 오차 모델을 행렬 형태로 기술하고, 측정 보정이 단일 기준국에서 다수 항공기에 동시에 적용될 때 발생할 수 있는 잠재적 오차 항목을 언급하라.
SBAS와 GBAS를 결합 운용할 경우, 보정 정보가 상호 보완적으로 작용하는 메커니즘을 해석하라.
해설 개요
SBAS 보정 메시지 형식(Range Correction, Ionosphere Correction 등)
GBAS 기준국 관측 방정식과 항공기 측정치 결합
결합 운용 시 오차 상쇄 효과
문제 7: 다중경로(Multipath) 오차 모델링
위성항법시스템에서 수신 신호는 지면, 건물, 기타 구조물 등에 반사되어 직접 경로 이외의 반사 경로(다중경로)로도 수신기에 도달할 수 있다. 이는 측위 오차에 큰 영향을 미치므로, 다중경로 오차를 모델링하고 보정하는 과정이 중요하다.
문제 요구 사항
위성–수신기 간 직접 경로와 다중경로 신호가 섞여 들어올 때, 코드측정치(유사거리 측정치)에 어떻게 영향을 주는지 간단한 수식으로 표현하라.
다중경로 모델링에서 고려되어야 하는 파라미터(반사 표면의 특성, 입사각, 편파 등)를 나열하고, 각각이 측정 오차에 미치는 정성적 영향을 설명하라.
실내/도심 협곡 환경에서 발생하는 다중경로 오차가 SBAS, GBAS, 또는 RTK(실시간 이동측위)에 어떤 문제를 야기할 수 있는지 언급하라.
해설 개요
다중경로 신호의 위상 및 지연 특성
반사 계수, 회절 효과 등을 포함한 간단한 반사 모델
보정 및 차분기법(Differential GNSS, RTK)과 다중경로의 관계
문제 8: 캐리어 위상 측정(Carrier Phase Measurement)
GNSS 측정에는 코드측정치뿐만 아니라 캐리어(반송파) 위상 측정이 활용된다. 캐리어 위상은 높은 주파수를 가지므로 더 정밀한 측위를 가능하게 하지만, 정수 모호성(ambiguity) 문제를 해결해야 한다.
문제 요구 사항
캐리어 위상 측정치가 코드측정치보다 더 정밀한 이유를 간단히 수식으로 설명하라.
위상 측정 식에 나타나는 정수 모호성( $N$ )의 물리적 의미를 기술하고, 이를 해결하기 위한 기본 아이디어(예: LAMBDA 방법, TCAR 등)를 개략적으로 설명하라.
정수 모호성을 잘못 추정했을 때, 측정치에 나타나는 단절(jump)의 크기를 파장(파장 길이, 예: L1 파장 약 19 cm)과 연관지어 서술하라.
해설 개요
캐리어 위상 측정식: $\rho_\phi = \lambda \left( N + \phi \right)$
초기화 단계에서 정수 모호성 해를 찾는 전형적 알고리즘
오차가 발생했을 경우 전체 측위에 미치는 영향
문제 9: 이중 차분(Differential) 측정 기법
GNSS 측정에서 오차를 줄이기 위해 사용하는 대표적인 방법이 이중 차분(Double Difference, DD) 기법이다. 이중 차분은 기준국과 이동국, 그리고 서로 다른 위성 간의 측정치 차분을 통해 위성 시계 오차, 수신기 시계 오차, 전리층·대류권 지연 등을 제거하거나 감소시킨다.
문제 요구 사항
이중 차분 기법에서 “Single Difference”와 “Double Difference”를 정의하는 식을 각각 제시하라.
이중 차분을 적용했을 때 제거되는/감소되는 오차 항목과 그대로 남아있는(또는 새로 생성되는) 오차 항목을 구분하라.
이중 차분 기법을 행렬 형태(예: DD 측정 벡터 $\mathbf{y}{DD}$, DD 설계행렬 $\mathbf{H}{DD}$)로 나타내어, 모호성 추정에 어떻게 활용되는지 간단하게 서술하라.
해설 개요
Single Difference: 위성 간 차분, 또는 수신기 간 차분
Double Difference: 위성 간 차분 + 수신기 간 차분 결합
남는 오차 항목(예: 멀티패스, 측정 노이즈 등), 제거 가능한 항목(위성 시계 오차 등)
문제 10: 정밀 단독 측위(Precise Point Positioning, PPP)
GNSS 측정 자료를 활용하여 단일 수신기(별도의 기준국 없이)만으로도 센티미터~데시미터급 정밀도를 달성하는 기법이 정밀 단독 측위(Precise Point Positioning)이다. 이를 위해서는 정밀 궤도력(Precise Orbit)과 정밀 시계(Precise Clock) 보정 정보가 필요하며, GNSS 측정 식에 다양한 오차 항을 추가로 포함한 모델링이 요구된다.
문제 요구 사항
PPP의 기본 식을 코드 측정치와 캐리어 위상 측정치 각각에 대해 간단히 제시하라.
코드 측정치(유사거리) 식에 정밀 궤도력, 정밀 시계, 전리층, 대류권 보정 항을 추가한 형태를 나타내라.
캐리어 위상 식에는 정수 모호성($N$)과 추가 보정 항을 포함하라.
위성 시계와 수신기 시계 보정 이외에, 지구 자전 보정, 상대론적 보정, 전리층·대류권 지연 보정 등이 PPP 식에 어떻게 반영되는지 서술하라.
PPP 해를 수렴시키는 데 걸리는 시간(수십 분~수 시간 규모)에 영향을 주는 주요 요인(위성 기하, 잡음 레벨, 모호성 해결 방법 등)을 정성적으로 언급하라.
해설 개요
PPP 정밀도 달성을 위해 필요한 정밀 궤도/시계 보정값
전리층 단일 주파수/이중 주파수 모델, 대류권 지연 모델링
정수 모호성 부호화(오차 모델)와 수렴 시간 단축 기법
문제 11: 전리층 및 대류권 오차 보정
GNSS 신호는 전리층과 대류권을 통과하면서 굴절/지연이 발생한다. 전리층은 전하 입자의 분포에 따라 신호 주파수별로 다른 지연 양상을 보이고, 대류권은 온도, 압력, 습도에 의해 신호 경로 길이가 변화한다. 정확한 보정을 위해 여러 모델(예: Klobuchar 모델, NeQuick 모델, Saastamoinen 모델 등)이 제안되어 있다.
문제 요구 사항
전리층 지연이 이중 주파수 관측치에서 어떻게 제거 또는 감소되는지, 간단한 식으로 나타내라.
예: $\mathbf{L}\mathrm{ion-free} = \alpha \mathbf{L}\mathrm{L1} + \beta \mathbf{L}_\mathrm{L2}$ 형태
대류권 지연(건조성분, 습윤성분)을 분리하여 모델링하는 식을 블록 수식으로 나타내라.
보정 모델(전리층, 대류권)이 정확하지 않을 경우, 잔여 오차가 측위 결과에 미치는 정성적 영향을 서술하라.
해설 개요
전리층 지연의 주파수 의존적 특성: 코드/위상 측정치 구분
대류권 지연 분해(건조/습윤 성분) 및 천정 관측 각도 변화 고려
모델 잔차 오차가 집단적으로 큰 편차를 유발할 가능성
문제 12: GNSS/INS 융합 측위
단독 GNSS 기반 측위는 신호 가시성이 떨어지거나 (도심 협곡, 터널 등) 다중경로가 심한 환경에서 오차가 커질 수 있다. 이를 보완하기 위해 관성항법장치(INS)와 GNSS를 융합하여 일관된 항법 해결책을 제공하는 방식이 있다.
문제 요구 사항
GNSS/INS 융합 시스템에서 상태벡터 $\mathbf{x}$가 위치, 속도, 자세 등으로 구성될 때, 이를 동적 모델 형태로 나타내라.
예: $\mathbf{x}_{k+1} = \mathbf{F}(\mathbf{x}_k, \mathbf{u}_k) + \mathbf{w}_k$ (여기서 $\mathbf{u}_k$는 가속도/각속도 측정치, $\mathbf{w}_k$는 시스템 노이즈)
GNSS 관측 모델(측위 방정식)과 INS 동적 모델을 칼만필터(또는 확장 칼만필터, EKF)로 결합하는 과정을 간략히 식으로 제시하라.
예: 예측 단계, 업데이트 단계에서의 상태 및 공분산 갱신 식
INS 단독 항법 시 누적되는 바이어스, 스케일 팩터 등의 오차가 GNSS로부터 어떻게 보정(추정)되는지 정성적으로 기술하라.
해설 개요
INS 측정방정식: 가속도/각속도 적분 -> 위치, 자세 추정
GNSS 관측방정식과의 결합 -> 실제 위치 오차 추정 및 보정
하이브리드 항법 시스템의 신뢰성/연속성 확보
Last updated