본문 바로가기

ECEF6

[INS] 관성항법 방정식 (Kinematics of Inertial Navigation) 관성항법시스템 (INS, inertial navigation system)은 관성센서인 가속도계와 자이로스코프에서 측정된 가속도(specific force, 단위 질량당 힘)와 각속도를 이용하여 운동체의 위치, 속도, 자세를 추정하는 항법시스템의 한 종류다. 관성항법시스템은 단독으로 위치, 속도, 자세를 추정할 수 있기 때문에 외부 센서나 신호를 사용할 수 없거나 또는 신뢰할 수 없는 상황, 예를 들면 수중이나 지하에서 특히 유용하다. 관성항법의 구현 방식에는 운동체의 회전 운동과 기계적으로 분리된 짐벌(Gimbal)에 관성센서를 장착하는 짐벌 (Gimbaled INS) 방식과, 운동체에 관성센서를 직접 부착하는 스트랩다운 (Strapdown INS) 방식으로 나눌 수 있다. 최신 관성항법시스템은 대부분.. 2024. 3. 3.
NED 좌표계와 LLH 좌표계간의 속도 변환식 미사일의 속도 \(\vec{V}\) 는 ECEF 기준, 즉 지면 기준의 상대적인 속도이므로 다음과 같다. \[ \vec{V} = \frac{^ed\vec{r}}{dt} \tag{1} \] 여기서 \(\vec{r}\) 는 지구 중심에서 미사일까지의 위치벡터다. 한편 미사일의 질량 중심점을 원점으로 하고 속도 방향을 x축, 지표면과 수평인 평면에서 속도 방향의 오른쪽 방향을 y축으로 하는 좌표계를 미사일 운동 좌표계 \(\{d\}\) 로 정의하므로 속도벡터를 좌표계 \(\{d\}\) 로 표시하면 다음과 같다. \[ V^d= \begin{bmatrix} V \\ 0 \\ 0 \end{bmatrix} \tag{2} \] 여기서 \(V\) 는 \(\vec{V}\) 의 크기다. NED(North-East-Down).. 2022. 1. 9.
ECEF-LLH 좌표계 상호 변환 매트랩 코드 LLH 좌표계에서 ECEF좌표계로 좌표변환하는 문제를 알고리즘 형태로 정리하면 다음과 같다. 입력: 위도 (\(\lambda_{lat}\)), 경도 (\(\lambda_{lon}\)), 높이 (\(h\)) 1. 접선반경 (\(R_{tr}\)) 계산: \(R_{tr}=\frac{ R_{eq}}{ \sqrt{1-e_{er}^2 \sin^2 \lambda_{lat}}}\) 2. 벡터 \(r^e\) 계산: \(r^e= \begin{bmatrix} (R_{tr}+h) \cos \lambda_{lat} \cos \lambda_{lon} \\ (R_{tr}+h) \cos \lambda_{lat} \sin \lambda_{lon} \\ \left( R_{tr} (1-e_{er}^2 )+h \right) \sin \la.. 2022. 1. 1.
ECEF 좌표계와 LLH 좌표계 지구중심지구고정 좌표계(ECEF, earth-centered earth-fixed frame)는 지구의 중심에 원점이 위치하며 지구에 고정되어 있어서 지구와 함께 자전하는 좌표계이다. 지구와 함께 자전한다는 점에서 ECI 좌표계와는 다르다. 기호로는 {e}로 표시한다. 좌표계의 \(\hat{e}_1-\hat{e}_2\) 평면은 지구의 적도면에 위치한다. \(\hat{e}_3\) 축은 ECI 좌표계의 \(\hat{i}_3\) 와 같은 방향으로 지구의 자전축 방향이며 \(\hat{e}_1\) 축은 지구 적도와 그리니치(Greenwich) 자오선이 만나는 점을 향한다. \(\hat{e}_2\) 축은 오른손 법칙에 의해 정해진다. ECI 좌표계에 대한 ECEF 좌표계의 각속도 벡터는 \(^i \vec{\omeg.. 2021. 12. 30.
ECEF 좌표계에서 미사일 운동 방정식 유도 지구 중심에서 미사일의 위치까지의 위치 벡터를 \(\vec{r}\) 이라고 하고 미사일을 질량 \(m\) 인 질점이라고 가정하면, 뉴턴의 운동법칙에 의해서 미사일 운동 방정식은 다음과 같이 주어진다. \[ \frac{^id}{dt} \left( m \frac{^id\vec{r}}{dt} \right) = \vec{L}+\vec{D}+m \vec{g} \tag{1} \] 여기서 \(\vec{L}\) 은 양력, \(\vec{D}\) 는 항력, \(\vec{g}\) 는 중력가속도다. 식 (1)에서 중요한 점은 질량 \(m\) 이 상수가 아니라 시간의 함수라는 것이다. 그럼에도 불구하고 식 (1)을 아래 식과 같이 미분하면 안된다. \[ \frac{^id}{dt} \left( m \frac{^id \vec{r}.. 2021. 12. 21.
미사일 좌표계의 정의 미사일 운동 방정식을 세우기 위해서는 상황에 따라 다음과 같이 여러 개의 좌표계가 필요하다. (1) ECI (earth-centered inertial)와 ECEF (earth-centered earth-fixed) 좌표계: ECI 좌표계 \(\{i\}\) 와 ECEF 좌표계 \(\{e\}\) 좌표계는 다음 그림과 같이 정의한다. ECI 좌표계에 대한 ECEF 좌표계의 각속도 벡터는 \[ ^i \vec{\omega}^e = \omega_{ie} \hat{e}_3 \tag{1} \] 이며 지구자전 각속도 \(\omega_{ie}\) 는 약 \(360^0/day\) 로서 WGS-84(World Geodetic System 1984)의 국제 표준값은 \(\omega_{ie} = 7.291151467 \time.. 2021. 12. 20.