![Korean Institute of Information Technology](/images/basic/pubreader_back.gif)
기동특성을 사용한 우주표적들 간의 변별 성능 분석
![](/_common/images/orcid.gif)
![](/_common/images/orcid.gif)
![](/_common/images/orcid.gif)
![](/_common/images/orcid.gif)
![](/_common/images/orcid.gif)
![](/_common/images/orcid.gif)
초록
탄도미사일의 탐지 및 추적은 우주표적들 간의 변별 문제를 해결해야 가능하다. 대표적인 우주표적으로는 탄두와 기만체가 있으며, 이들은 형상정보가 매우 유사하기 때문에 레이더에서 획득되는 RCS(Radar Coss Sction) 및 HRRP(High Resolution Range Profile) 정보만으로는 정확한 변별이 어렵다. 본 연구에서는 효율적인 우주표적 간의 변별을 위하여 기동특성으로부터 탄도계수를 추출하는 기법을 제안한다. 제안된 기법은 탄두와 기만체의 기동 궤적을 칼만필터로 각각 추적한 후, 고도 100km이하의 내기권에서 각 표적들의 탄도계수를 실시간 추정하여 가장 낮은 값을 가지는 표적을 탄두로 변별한다. 물리광학(Physical Optics) 기법을 이용하여 실제 미세운동에 따른 표적의 RCS 변화를 이용한 시뮬레이션 결과, 내기권에서 추출된 탄도계수로부터 효율적인 변별이 수행되는 것을 확인할 수 있었다.
Abstract
Detection and tracking of ballistic missile can be possible when the discrimination between space-targets is possible. Since representative space-targets such as warhead and decoy have similar shapes, radar crosssection (RCS) and high resolution range profile are not suitable for discrimination. In this paper, an effective discrimination method based on the ballistic factor is proposed. The proposed method uses Kalman-filter to track a real-time maneuvering trajectory and extracts the ballistic factor from an estimated trajectory at altitudes less than 100km. Since the ballistic factor of warhead is generally larger than decoys, effective discrimination is conducted. In simulations using the RCS values predicted by the physical optics in which targets are engaged in high-maneuver and micro-motion, the proposed method efficiently discriminated targets.
Keywords:
micro-doppler, ballistic missile, ballistic factor, radar, ISARⅠ. 서 론
탄도미사일은 일반적으로 추진, 중간 및 종착의 3단계로 구성된 기동궤적을 가진다[1]. 먼저 추진단계에서 모든 연료를 소모하여 고도 100km 이상의 외기권으로 진입한 후, 소모된 연료탱크는 분리되고 표적은 매우 낮은 대기밀도의 중간단계에 놓이게 된다. 이때 탄도미사일은 지상레이더 및 조기경보레이더에 의한 탐지 및 추적을 피하기 위해 탄두와 기만체로 구성된 우주표적들로 분리된다[2-4]. 기만체의 경우, 탄두와 유사한 형상을 가지는 우주표적이기 때문에 대표적인 레이더 반사신호인 레이더 단면적(RCS, Radar Cross Section), 고해상도 거리측면도(HRRP, High Resolution Range Profile) 및 역합성개구면레이더(ISAR, Inverse Synthetic Aperture Radar) 영상[5]으로는 탄두와의 정확한 변별이 불가능하다. 따라서 이들을 보다 효율적으로 변별하기 위해서는 형상학적 정보가 아닌 다른 특징이 추출되어야 한다.
변별을 위한 대표적인 특징으로는 미세운동 및 기동 특성이 있다. 먼저 미세운동은 형상학적 특징 및 자세제어장치 유무로 인하여 발생하며, 이를 바탕으로 한 변별 연구는 최근까지 다양하게 수행되어 왔다. 하지만 미세운동 특징은 시간-주파수 영상을 기반으로 많은 계산량을 요구하기 때문에 실용성 측면에서 문제점이 제기되어 왔다[6]. 반면에 기동특성의 경우에는 국내외적으로 변별 연구 결과가 많지는 않지만, 실시간 추적을 통한 특징 추출이 가능하다는 이점이 있다. 따라서 기동특성으로부터 효율적인 특징을 추출 할 수 있을 경우, 보다 실용적인 변별이 가능할 것이다.
우주표적들의 기동특성은 기동궤적의 중간단계에서 낮은 대기밀도로 인하여 중력과 항력에 의해서만 정의된다[7]. 이때 중력가속도의 영향을 제외하면 탄두와 기만체 간의 기동특성 차이점은 오직 항력 가속도에 의해서만 결정된다고 볼 수 있다. 하지만, 중간단계에서는 매우 낮은 대기밀도로 인하여 이와 비례하는 항력 가속도가 매우 작으며, 사실상 중력가속도에 의한 영향이 지배적으로 나타난다. 그러나 고도 100km이하의 종착단계에 진입할 경우, 대기밀도는 기하급수적으로 증가하면서 탄두와 기만체간의 항력가속도 차이도 크게 나타나기 시작한다. 이때 항력가속도의 차이점은 주로 표적의 물리적 특징인 질량, 표면적 및 항력계수로 인해 발생되며, 이 3가지 변수들로 구성된 수치 값을 일반적으로 탄도계수라고 정의한다. 그리고 보통 탄두는 핵과 같은 무거운 물질들로 구성되기 때문에 상대적으로 기만체보다 질량이 큰 특징을 가지며, 이 특징은 질량과 반비례 관계인 탄도계수에서도 쉽게 확인될 수 있다. 따라서 종착단계에서 나타난 우주표적들 중 가장 낮은 탄도계수를 가지는 표적을 탄두로 선택 할 경우, 효율적인 변별을 수행할 수 있게 된다.
본 논문에서는 종착단계에서 탄두와 기만체간의 변별을 위하여 탄도계수를 특징으로 추출하는 알고리즘을 바탕으로 변별 성능을 분석하였다. 먼저 제안된 기법은 표적의 기동궤적을 추정하기 위하여 칼만필터(KF, Kalman Filter)를 사용한 추적단계로 시작한다. 이때 우주표적들로부터 획득되는 RCS는 표적 형상 및 관측각도에 따라 매우 민감하게 변화기 때문에 KF 성능에도 큰 영향을 준다. 따라서 탄두와 기만체의 형상에 대한 고주파 수치해석기법인 물리광학(PO, Physical Optics)를 사용하여 RCS를 분석한 후, 기동 및 미세운동을 가지는 우주표적과 레이더 간의 관측각도에 따른 RCS 변화를 고려하여 레이더 신호모델링을 수행하였다. 이후 KF로 추정한 기동궤적 정보를 바탕으로 실시간 탄도계수 추출을 수행한 뒤, 최종적으로 우주표적들로부터 추출된 탄도계수들 중 가장 작은 값을 가지는 표적을 탄두로 변별한다. 본 논문에서는 제안된 기법의 효용성 검증을 위하여 원뿔형상 탄두와 원뿔, 원통형상의 2가지 기만체들을 사용한 변별 시뮬레이션을 수행하였다.
Ⅱ. 수학적 모델링
2.1 기동궤적
중간 및 종착단계에 있는 탄두의 기동궤적 모델링을 위해 먼저 공기 저항력 D 을 다음과 같이 고려한다[7].
(1) |
여기서 ρ(t)는 대기밀도, H 는 고도, S 는 표면적, CD는 항력계수 그리고 V 는 속도이다. 이때 공기 저항력에 대한 항력가속도 aD는 뉴턴 제2법칙에 따라 다음과 같이 정의된다.
(2) |
여기서 m는 질량이다. 만약 초기발사각도 ϕ으로 추진될 경우, 수평 및 수직방향의 가속도 ax ,az는 다음과 같다.
(3) |
여기서 g0는 중력가속도, re는 지구반경 그리고 r 는 지구중심에서 표적까지 거리이다. 따라서 최종적으로 수평 및 수직방향의 궤적 X (t) ,Z (t)은 두 번의 적분연산을 통하여 다음과 같이 정의된다.
(4) |
여기서 V0는 초기발사속도이다. 기만체의 기동궤적은 표적 형상 정보에 맞는 3가지 변수 m, S 및 CD를 선택하여 탄두 궤적 모델링과 동일한 방식으로 구성된다.
2.2 미세운동
탄두의 대표적인 미세운동은 원추(Coning) 및 장동(Nutation)이 있으며, 장동은 미세한 진동운동이기 때문에 본 연구에서는 관측각도 변화에 큰 영향을 주는 원추운동만을 고려하였다. 원추운동은 기동 방향을 기준으로 원추각도 θc만큼 표적이 기울어진 상태에서 기동방향을 회전축 으로 회전하는 운동이며, 로드리게스 회전공식을 적용할 경우, 시간 t 에 대하여 회전 행렬 Rc을 다음과 같이 간단히 표현할 수 있다[8][9].
(5) |
여기서 I 는 단위행렬, Ωc는 원추운동 각속도 그리고 는 다음과 같다.
(6) |
따라서 실시간 탄두의 위치벡터 는 지역좌표계 내에서 다음과 같이 정의된다(그림 1).
(7) |
여기서 는 탄두의 초기위치벡터이다.
기만체의 경우, 대표적인 미세운동으로 떨림(Wobble)을 고려하였다(그림 2).
이는 임의의 회전축 벡터를 법선으로 가지는 평면(떨림운동 평면)에서 정현파적으로 ± θw0 범위에서 회전하는 것으로 가정할 수 있으며, 이러한 떨림 운동 평면을 정의하기 위해서는 새로운 좌표계를 사용해야 한다. 이를 위하여 지역좌표계에서 회전축 벡터를 z′축으로 가정할 경우, 떨림운동 평면은 x′y′평면으로 정의되며 각 축의 단위벡터를 으로 계산한 후, 변환행렬 을 다음과 같이 구성할 수 있다.
(8) |
이를 기만체 위치벡터와 곱할 경우, 새로운 좌표계로 변환된다. 이때 변환된 좌표계에서 떨림운동의 회전특성을 적용시키기 위한 떨림운동 회전행렬 B 는 다음과 같이 정의된다.
(9) |
여기서 θw0는 떨림운동 크기, Ωw는 떨림운동의 각 속도이다. 변환된 새로운 좌표계에서 회전 행렬 B 를 이용하여 정현파적으로 회전시킨 후, 원래의 좌표계로 역 변환하는 행렬인 AT를 곱하여 최종적으로 떨림운동으로 인한 실시간 위치 좌표 를 다음과 같이 모델링하였다[7].
(10) |
여기서 는 기만체의 초기위치벡터이다.
2.3 RCS
탄두와 기만체는 대표적으로 원뿔형상을 가지기 때문에 X밴드(10GHz)를 사용하는 레이더 입장에서 PO를 사용한 고주파 근사를 통해 RCS 분석이 가능하다. 먼저 원뿔 형상의 경우 RCS σf (θ)는 다음과 같이 정의된다[10].
(11) |
여기서 λ는 파장, h는 원뿔 높이, df는 반지름 그리고 θ는 관측각도이다. 그리고 추가적으로 원통형상 기만체의 RCS σc(θ)는 다음과 같다[10].
(12) |
여기서 dc는 원통의 반지름이다. 이때 미세운동 및 기동하는 표적의 위치벡터와 레이더 시선방향벡터 간의 관측각도 θ를 계산할 수 있을 경우, 실시간 RCS 변화를 계산할 수 있다. 먼저 국부좌표계에서 미세운동 및 기동하는 탄두의 위치벡터는 식 (4)와 식 (7)을 사용하여 다음과 같이 정의할 수 있다.
(13) |
만약 레이더 시선방향벡터가 일 경우, 탄두에 대한 관측각도 θ는 제 2 코사인법칙에 따라 다음과 같이 정의될 수 있다.
(14) |
여기서 • 는 내적 연산자이다. 기만체의 경우, 식 (13)에서 기만체의 기동궤적 및 식 (10)을 사용하여 수행된다.
Ⅲ. 제안된 기법
제안된 기법은 탄두와 기만체들에 대한 KF 기반의 기동궤적 추적단계와 탄도계수 특징 추출단계로 나뉜다.
3.1 기동궤적 추적단계
본 연구에서는 표적의 RCS 변화에 따른 탄도계수 특징추출 성능을 집중적으로 분석하기 위하여 확장 칼만필터(Extended Kalman-Filter)가 아닌 다소 간단한 KF를 추적기법으로 사용하였다. 먼저 초기 6차원 위치 및 속도벡터 와 초기 6×6 측정오차 공분산 행렬 P0에 대한 예측과정을 시작으로 다음과 같이 수행된다.
(15) |
(16) |
여기서 와 는 k번째 예측 값들이며, Q 는 시스템 내에서 발생하는 잡음오차 공분산 행렬, 그리고 A 는 상태벡터를 예측하기 위한 시스템 행렬로 다음과 같이 정의된다.
(17) |
여기서 dt 는 추적 주기이다. 이때 예측된 상태벡터로 위치 및 속도벡터를 추정하기 위해 k번째 칼만 이득 Kk를 아래와 같이 정의한다.
(18) |
위 식에서 R 는 레이더에 들어오는 잡음오차 공분산 행렬, H = I6 × 6 는 관측 행렬로 정규직교 단위 행렬이다. 따라서 k번째 위치 및 속도벡터는 다음과 같이 추정된다.
(19) |
위 식에서 는 현재 레이더로 측정된 표적의 위치벡터이며, 표적이 탄두일 경우 다음과 같이 정의된다.
(20) |
식 (20)은 표적의 RCS 크기 정보에 따라 탐지 성공 유무가 결정된다. 본 연구에서는 RCS 탐지 임계값에 따른 기동 및 미세운동을 가지는 표적의 탐지 및 추적 가능성을 함께 분석한다.
3.2 탄도계수 추출단계
탄도계수 γ 는 식 (1)에서 정의한 질량, 표면적 및 항력계수로 구성되며, 다음과 같이 식 (2)와 연관성을 가진다.
(21) |
이때 항력가속도 aD , 속도 V 및 고도 H 를 추정할 경우, 표적 고유의 물리적 특징인 탄도계수를 추출 할 수 있다. 만약 KF를 통해 표적의 실시간 위치 및 속도벡터 를 획득할 경우, 3가지 변수들은 다음과 같이 추정된다.
(22) |
(23) |
(24) |
는 속도벡터의 미분연산을 통해 추정된 가속도벡터, 그리고 는 실시간 위치벡터로 계산된 중력가속도이다. 제안된 알고리즘 절차는 그림 3과 같다.
Ⅳ. 시뮬레이션 결과
시뮬레이션을 위하여 먼저 표 1, 2의 탄두와 기만체들에 대한 형상학적 특징 및 기동변수들을 바탕으로 2차원 기동 궤적을 구성하였다(그림 4). 이 때 2대의 기만체들은 외기권 내 중간단계 지점인 224.45초에서 탄두와 함께 분리된 후, 330초에서 종착단계에 진입하게 된다. 이때 각각의 미세운동도 기동궤적과 함께 모델링하였다(표 3).
본 연구에서는 관측각도 변화에 따른 탄도계수 추출 성능 분석을 위하여 지상레이더 위치에 대한 2가지 시나리오를 고려하였다. 첫 번째 시나리오에서 레이더는 미사일 사거리 중앙 [250,0,0] km에 위치되고, 두 번째 시나리오는 미사일 종착지 근처 [450,0,0] km에 위치된다.
4.1 [250, 0, 0] km의 지상 레이더
첫 번째 시나리오에서 종착단계 진입 시점인 330~331초 동안 나타나는 우주표적들의 RCS 변화 특성은 그림 5와 같다. 고도 100km인 330초에서 탄두의 초기관측각도는 157°이고, 최대 RCS는 약 5dB로 나타났다(그림 5(a)). 이와 동시에 미세운동에 의하여 관측각도는 90°~157°까지 변화하였으며, 그 결과 최대 53dB의 RCS변화폭을 가지면서 주기적으로 나타남을 알 수 있었다. 반면에 기만체들의 경우, 미세운동 역학특성으로 인해 상대적으로 극심하게 RCS가 변화함을 알 수 있었다(그림 5). 이때 RCS 탐지 임계값을 -30dB로 설정할 경우, 탄두의 RCS는 330.06초에서 최대 값이 나타나고 0.2초 후에 다시 최대 RCS가 나타남을 확인하였다(그림 5(a)). 이는 원추운동 각속도 Ωc =31.4(=2π(5))이기 때문에 RCS가 주기적으로 나타나는 것이다. 따라서 KF를 사용한 최적의 탄두 탐지 및 추적 방법은 330.06초부터 0.2초 간격으로 수행되어야 한다.
만약 최적화된 KF(표 4)를 사용 할 경우, 종착단계 구간인 330~400초에서 그림 6과 같이 탄도계수가 추출되었다. 탄도계수의 이론값은 식 (21)에 표 1의 질량, 표면적 및 항력계수 값들을 대입 할 경우 추정될 수 있으며, 그 결과 이론값 1.82×10-4은 제안된 기법으로 추정된 결과 1.84×10-4와 유사함을 확인할 수 있었다. 마찬가지로 2가지 기만체들의 이론값 6.38×10-4, 6.83×10-4도 추정값 0.64×10-4 , 0.69×10-4과 유사하였으며, 탄두의 탄도계수가 가장 작은 값을 가짐을 확인할 수 있었다.
4.2 [450, 0, 0] km의 지상 레이더
두 번째 시나리오에서는 지상 레이더가 미사일 종착지 부근에 위치하기 때문에, RCS 변화특성이 다르게 나타났다(그림 7). 먼저 330초에서 탄두의 초기관측각도는 70°이고, 미세운동에 의하여 12°~128°의 관측각도 변화가 나타났다. 이때 그림 7(a)에서 RCS 변화폭은 첫 번째 시나리오와 유사하지만, 변화주기는 두 배로 증가함을 볼 수 있었다. 이러한 점은 원뿔 형상의 RCS가 관측각도 75°에서 최대값을 가지기 때문이다. 반면에 기만체의 경우에는 원통 기만체의 최대 RCS가 첫 번째 시나리오보다 50dB만큼 증가하였다(그림 7(c)). 이는 원통 RCS가 관측각도 90°에서 매우 크게 나타나기 때문이다.
앞서 첫 번째 시나리오에서 최적화 된 KF를 그대로 사용할 경우, 두 번째 시나리오에서 우주표적들의 탐지 및 추적이 불가능하였다. 이는 매 탐지순간마다 나타나는 RCS가 RCS 탐지 임계값 –30dB 보다 낮게 나타나기 때문이다. 따라서 RCS를 단순히 통계적 모델링으로 수행하거나 미세운동영향을 무시할 경우, 추적 및 탐지가 실패할 확률이 매우 높음을 의미한다.
Ⅴ. 결 론
본 연구에서는 기동특성에서 추출된 탄도계수를 기반으로 우주표적들 간의 변별 기법 성능을 분석하였다. 일반적으로 기존의 추적기법들은 단순히 표적들을 통계적 모델링으로 가정하지만, 본 연구에서는 실제 우주표적들의 형상들에 PO 기반의 수치해석을 통한 RCS 분석과 기동 및 미세운동에 의한 관측각도 변화를 함께 고려하였다. 시뮬레이션 결과, 탐지 및 추적의 성공여부는 관측각도 변화에 따른 RCS에 매우 의존적으로 나타났으며, RCS 변화특성에 최적화된 KF를 사용할 경우 탄도계수가 성공적으로 추출됨을 확인할 수 있었다.
향후 본 연구에서는 우주표적들의 특징들을 바탕으로 최적의 탐지 및 추적기법을 개발할 예정이다. 특히, 실제 추적 환경 특성을 고려하기 위하여 모노펄스 레이더 기반의 거리, 방위각 및 고각 추정과 최적화된 비선형 확장칼만필터를 통한 탄도계수 추출을 수행할 것이다.
Acknowledgments
방위사업청과 국방과학연구소가 지원하는 레이다/IR 표적식별 특화연구실 사업의 일환으로 수행되었음
References
- N. George, and A. Theodore, "Future challenges to ballistic missile defense", IEEE spectrum, Sep), (1997.
- Y. X. Liu, X. Li, and Z. W. Zhuang, "The presession signature of aerospace target and it’s application in radar target recognition", J. Develop. Natural Sci., 14(11), p1329-1332, (2004).
-
H. Rathore, M. Katukuri, D. Rao, and B. N. Murthy, "Real-time embedded software design for onboard computer of anti-ballistic missile system", ICCCT. 2014 International Conference on, p1-5, Dec), (2014.
[https://doi.org/10.1109/iccct2.2014.7066741]
- T. T. H. Le, J. H. Kim, J. H. Kim, and H. W. Kim, "An Approach to Build an Efficient Intrusion Detection Classifier", Journal of Platform Technology, 3(4), p43-52, Dec), (2015.
-
I. O Choi, J. H. Jung, K. T. Kim, and S. H. Park, "An improved method to recognize range profiles using statistical hypothesis test and principal component analysis", Journal of KIIT, 15(1), p115-122, Jun), (2013.
[https://doi.org/10.14801/jkiit.2017.15.1.115]
-
I. O Choi, S. H. Kim, J. H. Jung, K. T. Kim, and S. H. Park, "An efficient method to extract the micro-motion parameter of the missile using the time-frequency image", Journal of KIEES, 27(6), p557-565, Jun), (2016.
[https://doi.org/10.5515/kjkiees.2016.27.6.557]
-
H. S. Kim, "A Requirement Assessment Algorithm for Anti-Ballistic Missile Considering Ballistic Missile's Flight Characteristics", KIMST, (2011).
[https://doi.org/10.9766/kimst.2011.14.6.1009]
- V. C. Chen, F. Li, S. Ho, and H. Wechsler, "Micro-Doppler effect in radar", IEEE Trans. Aero. Electron Syst., 42(1), p2-21, Oct), (2006.
-
I. O. Choi, S. H. Kim, J. H. Jung, K. T. Kim, and S. H. Park, "An Efficient Method to Extract the Micro-Motion Parameter of the Missile Using the Time-Frequency Image", JEES, 27(6), p557-565, Jun), (2016.
[https://doi.org/10.5515/kjkiees.2016.27.6.557]
- B. Mahafza, Radar Systems and Design using MATLAB, Champman & Hall/CRC Press LLC, Jan), (2000.
2012년 8월 : 부경대학교 전자정보통신공학과(공학사)
2014년 8월 : 부경대학교 전자공학과(공학석사)
2015년 9월 ~ 현재 : 포항공과 대학교 전자전기공학과 박사과정
관심분야 : 미세도플러 분석, 레이다 표적인식, 레이다 신호처리, 패턴인식, 스펙트럼 추정
2015년 2월 : 부경대학교 전자공학과(공학사)
2017년 2월 : 부경대학교 전자공학과(공학석사)
2017년 2월 ~ 현재 : 포항공과 대학교 전자전기공학과 연구원
관심분야 : 레이다 표적식별, 레이다 신호처리, 미세도플러 분석 등
1998년 2월 : 경북대학교 전자공학과(공학사)
2001년 2월 : 경북대학교 전자공학과(공학석사)
2006년 2월 : 경북대학교 전자공학과(공학박사)
2011년 1월 ~ 현재 : 국방과학연구소 연구원
관심분야 : 레이다 신호처리, 레이다 표적 인식
1991년 2월 : 공군사관학교 전자공학과(공학사)
1995년 2월 : 서울대학교 전자공학과(공학사)
1998년 2월 : 포항공과대학교 전자전기공학과(공학석사)
2007년 2월 : 포항공과대학교 전자전기공학과(공학박사)
2013년 2월 : 고려대학교 경영학과(MBA)
2013년 3월 ~ 현재 : 한택에너지 대표이사
2016년 1월 ~ 현재 : KAIST 국방무인화기술 특화연구센터 연구 교수
관심분야 : 풍력단지 EMI/EMC 분석, 포탄 RCS 분석, 레이다 표적식별, 방산경제 및 R/D 전략 등
1994년 2월 : 포항공과대학교 전자전기공학과(공학사)
1996년 2월 : 포항공과대학교 전자전기공학과(공학석사)
1999년 2월 : 포항공과대학교 전자전기공학과(공학박사)
2002년 3월 ~ 2011년 2월 : 영남대학교 전자공학과 교수
2011년 3월 ~ 현재 : 포항공과대학교 전자전기공학과 교수
2012년 9월 ~ 현재 : 레이다/IR 표적식별 특화연구실 실장
관심분야 : 레이다 신호 처리 및 영상, 레이다 표적인식 및 패턴인식, 전자기수치해석 및 RCS 측정
2004년 2월 : 포항공과대학교 전자전기공학과(공학사)
2007년 2월 : 포항공과대학교 전자전기공학과(공학석사)
2010년 2월 : 포항공과대학교 전자전기공학과(공학박사)
2010년 9월 ~ 현재 : 부경대학교 전자공학과 부교수
관심분야 : 미세도플러 분석, 풍력단지 EMI/EMC 분석, 레이다표적인식, 레이다영상, 레이다 신호처리 등