Ⅰ. 서 론
최근 스텔스 기능을 보유한 저피탐 항공기가 지속적으로 개발됨에 따라 종래의 레이다 시스템으로는 이러한 항공기에 대한 추적이 어려운 실정이다. 이에 대한 대안으로 저주파 대역 방송신호를 이용한 수동형 레이다에 관한 연구가 활발히 진행되었다[1]～[7]. 저주파 대역에서는 특정 주파수 대역에 적합하도록 설계된 전파흡수도료(radar absorbing material: RAM)나 저피탐 형상 기술(RAS) 등이 효과적이지 못하다고 알려져 있다. 또한 동시에 다수의 방송신호를 기반으로 멀티스태틱(multistatic) 레이다를 구현할 수 있으므로 모노스태틱(monostatic) 레이다에 비해 표적의 위치를 보다 정밀히 추적할 수 있다. 이와 더불어 저주파 레이다 시스템은 생산비용을 크게 감소시킬 수 있는 장점이 있다. 국내의 여러 무선 방송서비스 중 FM 라디오 방송은 오랜 기간동안 서비스가 되어 넓은 커버리지와 다수의 송출 방송국을 보유하고 있다. 따라서 본 논문에서는 FM 대역 방송신호 기반의 멀티스태틱 수동형 레이다, 즉 passive coherent location(PCL)의 개발을 위한 산란신호 모의 방법에 대해 제안하였다.
본 논문의 2장에서는 산란신호 모의 방법에 대해 살펴본 후, 3장에서는 이를 구현하는 방법에 대해 설명한다. 이후 4장에서는 구현된 산란신호를 검증하고, 본 논문의 결론을 맺는다.
Ⅱ. 산란신호 모의
멀티스태틱 PCL 장비를 통해 수신되는 FM 방송신호는 송신소로부터 최단거리 경로를 통해 PCL 수신기로 전달되는 기준신호, 주변 지형물들의 반사를 통해 수신되는 클러터 신호, 그리고 움직이는 표적에 반사되어 수신되는 표적반사신호로 구분할 수 있다.
국내의 FM 방송신호는 다음과 같이 변조된다[8]. 우선 FM 변조 대상이 되는 기저대역 신호는 스테레오포닉(stereophonic) 변조 과정을 거친다.
이때 기저대역 신호는 인간이 들을 수 있는 가청 주파수 대역의 신호로 좌/우 채널에 따라 L(t)와 R(t)로 표현되며, C1과 C2는 각각 0.45, 0.1의 상대적 크기를 가진다. 또한 fp는 파일롯 톤의 주파수로 19 kHz이며, 복조시에 사용된다. 따라서 좌/우 합채널은 DC 부근으로, 좌/우 차채널은 파일롯 톤의 2배인 38 kHz로 AM 변조된다. 이후 스테레오포닉 변조된 신호 m(t)는 규약에 따라 아래와 같이 FM 변조가 진행된다.
fΔ는 변조주파수로 국내 FM 방송 규약에 따라 75 kHz가 된다. 이때 IQ 데이터는 x(t)에 반송파 exp(−j2πfct)를 곱한 후 기저대역 신호만 추출하면 된다. 따라서 FM 변조된 IQ 신호는 다음과 같다.
송신소에서 수신기로 직접 수신되는 신호는 기준신호로 활용된다. 따라서 기준신호는 xIQ신호를 지연시킴으로써 쉽게 계산할 수 있다. 송신소와 수신기 사이의 거리를 Lref라고 하면 지연시간 Δtdirect는 c0Lref가 된다. 여기에서 c0는 빛의 속도를 나타내며, 신호세기의 감쇠는 Friis 방정식에 따라 나타낼 수 있다[9]. 따라서 기준신호의 IQ yref는 다음과 같다.
여기에서 Pt는 FM 방송국의 송출전력, Gt, Gr은 각각 송/수신 안테나 이득, 그리고 λ는 파장이다.
PCL 장비의 수신신호는 수신기로 직접 들어오는 신호 외에 주변의 산이나 건물 등 여러 장애물들로부터 반사되어 들어오는 신호도 존재한다. 이러한 신호들은 기준신호나 표적반사신호와 구분하기 위해 클러터 신호로 통칭한다. 클러터들의 레이다 반사 면적(radar cross section: RCS), 지연시간, 송신소와 클터터 사이의 거리, 클러터와 수신기 사이의 거리를 각각 σclutter, Δtclutter, Lt, Lr로 나타내면 레이다 방정식을 이용하여 다음과 같이 나타낼 수 있다[10].
그림 1과 같이 바이스태틱(bistatic) 레이다의 경우, 송신기와 수신기가 서로 떨어져 있으므로 표적의 이동속도 vr을 수신기에서 관찰된 표적의 속도로 치환하여 도플러 주파수를 계산해야 한다. 하지만 이를 등가적으로 전파의 이동경로(Dt + Dr)를 시간에 따른 변화량으로도 나타낼 수 있다[11].
이렇게 계산된 도플러 주파수는 레이다 방정식에 적용하여 아래와 같이 계산할 수 있다.
도플러 주파수는 시간에 따라 변하는 함수이기 때문에 이를 적분한 후 초기 반송파의 위상값 ϕfc(t0)를 보상해야 정확한 도플러 효과를 반영할 수 있다. 만약 표적 이동거리가 fs/c0보다 클 경우, IQ 데이터 xIQ 를 1/fs 만큼 천이시킨 후 도플러 주파수 적분항과 위상을 적절히 초기화해야 한다. 여기에서 fs는 샘플링 주파수이다.
Ⅲ. 개발 지원장비 구현
그림 2는 멀티스태틱 PCL 개발지원장비의 블록다이어 그램을 나타낸다. 개발지원장비는 운용 S/W, IQ 생성기, 모의신호 발생기로 구성된다.
운용자는 송신소와 수신기 위치, 주변 클러터 환경 설정, 표적의 비행 시나리오 등을 결정하면, 운용 S/W는 이에 따라 제어 메시지를 생성한다. IQ 생성기는 제어 메시지에 따라 FM 라디오 모의신호를 생성한 후, 이를 바탕으로 송신소로부터 직접 수신하는 기준신호, 주변 클러터에 산란되어 수신하는 클러터신호, 그리고 시나리오에 따라 움직이는 표적에 반사된 표적반사신호를 생성한다. 이후 각 신호들은 입사방향에 따라 배열 안테나의 어레이 매니폴드가 곱해진 후 합쳐진다. 배열 안테나 수신특성이 합성된 신호는 모의신호 발생기 내 digital analog convertor (DAC)의 동적범위를 고려해 높은 레벨과 낮은 레벨신호로 분리되어 모의신호 발생기로 전달된다. 모의신호 발생기에서 생성된 다중채널 FM 대역 신호는 실시간으로 생성되어 향후 멀티스태틱 PCL 주장비 개발에 이용된다. 따라서 주어진 자원내에서 가능한 정확히 산란신호를 모의해야 한다.
본 연구에서는 ray-tracing 기반의 wireless insite 상용 프로그램[12]을 이용하여 송신소들과 수신기 사이의 신호지연, 감쇠량 및 수신방향을 계산하였다. 그림 3은 식장산 FM 송신소로부터 수신한 신호의 시뮬레이션 결과(실선)와 측정결과(파선)를 나타낸다. 측정은 실측한 IQ데이터를 자기자신과 교차상관(cross correlation)한 결과이며, 클러터 반사신호가 수신된 영역에서는 다른 곳보다 상관계 수가 높게 계산된다. 측정은 총 5회 수행되었으며, 수신한 음성신호의 종류에 따라 다소 차이가 발생한다. 하지만 시뮬레이터를 통해 예측한 클러터 위치 및 반사신호 세기와 유사하게 상관계수가 측정된 것을 확인할 수 있다. 따라서 ray-tracing 기법을 이용한 클러터 분석이 FM대역에서도 효과적으로 사용될 수 있음을 확인하고, 이를 DB로 구축하여 활용하였다.
앞절에서 서술한 식 (6)과 식 (7)에 의해 계산되는 도플러 효과는 정확하지만, 많은 연산량을 필요로 한다. 따라서 실시간으로 IQ 데이터를 생성해야 하는 개발지원장비에서 사용하기에는 다소 부적합하다. 따라서 각 IQ샘플을 전파 도래시간 Δttarget만큼 지연을 시켜 도플러 효과를 다음과 같이 계산하였다.
이때 신호의 지연은 스테레오포닉 변조부터 FM 변조까지 모든 과정에서 Δttarget만큼 지연이 되어야 한다. 우선 스테레오포닉 변조는 음성신호의 좌/우 채널 신호를 Δttarget만큼 지연시킨다. 이때 음성신호 대역폭은 20 kHz 내외이지만 샘플링 대역폭은 그보다 수 배 이상 높으므로 FIR 필터를 기반으로 한 매트랩의 fractional delay 함수 등을 이용하여 효과적으로 신호지연을 시킬 수 있다[13]. 이때 발생되는 오차가 식 (7)과의 주요한 차이가 된다. 그 외에는 삼각함수의 지연이 되므로 해석적으로 정확한 지연이 가능하다.
이후 스테레오포닉 변조된 신호 mtarget(t)는 FM 변조가 된다. 이때 반드시 반송파의 위상 지연 exp(−j2πfcΔttarget)를 함께 계산해 주어야 한다.
마지막으로 표적반사신호는 표적의 바이스태틱 RCS를 고려하여 다음과 같이 계산한다.
식 (7)과 (10)의 검증을 위해 시변 도플러 주파수 편이 현상을 삼각함수를 이용해 해석적으로 비교하였다. 사용된 삼각함수는 각각 100 MHz, 100 MHz+125 kHz, 100 MHz+250 kHz 의 주파수 성분을 가지며, 진폭은 1, 0.5, 0.2 이다. 이때 신호원과 표적 사이의 왕복거리와 이에 따른 도플러 주파수를 그림 4(a)와 그림 4(b)에 각각 도시하였다. 시간에 따른 왕복거리가 포물선의 형태를 가지므로 이에 대한 미분인 도플러 주파수는 선형특성을 갖는다. 각각 왕복거리로부터 계산된 지연시간을 바탕으로 계산된 신호와 도플러 주파수를 곱하여 계산된 신호의 차이를 그림 4(c)에 도시하였다. 각 신호들의 차이가 수치해석 오차인 0.4 % 이하로 서로 잘 일치하는 것을 확인하였다.
표적의 RCS 분석에는 WIPL-D 상용 EM해석 툴을 이용하였다[14]. WIPL-D는 고차항 기저함수 (high order basis function)을 기반으로 유전체 등이 포함된 해석을 비교적 빠른시간에 정확히 해석할 수 있다. 그림 4는 98 MHz에 서 전파의 입사방향(화살표)에 따른 표적의 바이스태틱 RCS를 나타낸다. 표적은 PEC 구조위에 두께 1 mm인 RAM 도료가 있을 경우와 없을 경우에 대해 각각 분석을 하였다. FM 대역에서는 두 RCS가 2 dB 이내로 거의 동일하므로 FM 방송을 이용하는 PCL장비가 저피탐 스텔스 기능의 표적도 탐지가 가능함을 확인하였다. 본 연구에서는 각기 다른 표적의 바이스태틱 RCS를 분석하여 DB로 활용하였으며, 사용자가 임의로 RCS를 지정하거나 추가할 수 있도록 구성하였다.
주장비는 배열 안테나를 이용하여 FM 방송신호와 표적반사신호를 수신한다. 따라서 본 연구에서도 안테나의 어레이 매니폴드를 이용하여 각각의 안테나로 수신한 전파를 모델링하였다. 안테나는 FM 방송대역에서 반파장 길이를 갖는 다이폴을 적용한 원형 배열안테나이다. 그림 6은 각각 89 MHz에서 측정된 어레이 매니폴드(실선)와 시뮬레이션(점선) 결과를 나타내며, 두 결과가 서로 잘 일치한다. 또한 안테나 사이의 커플링 효과에 의해 다소 찌그러지는 구간이 발생한다. 이는 안테나로 입사되는 전파를 정확히 모델링하기 위해서는 안테나 간 커플링 효과를 고려해야 함을 의미한다. 안테나를 통해 수신된 신호는 다음과 같다.
여기에서 ych1은 1번 안테나를 통해 수신하는 IQ데이터를 의미한다. 따라서 기준신호와 N개의 클러터신호, M개의 표적신호가 각각 도래각 ϕ에 따른 어레이매니폴드 AF(ϕ)를 곱한 후 이를 더해서 구해진다.
Ⅳ. 모의신호 분석
최종적으로 IQ생성기에서 발생된 IQ신호가 적절히 PCL 신호를 모의하는지 확인하기 위해 ARD (Amplitude Range & Doppler) 지도를 도시하였다. 음성신호는 다양한 주파수 성분이 포함되어 있는 음악신호를 이용하였으며, 사용된 샘플은 약 0.15초 가량이다. ARD 지도는 다음과 같이 계산된다[10].
이때 도플러 주파수의 편이현상을 간단히 검증하기 위해 클러터 반사신호는 제외하였다. 표적은 ‘x’로 나타내었으며, ARD 지도에서 표적이 적절히 모의된 것을 확인 할 수 있다. 다만 표적이 위치하지 않는 부분에서도 값이 크게 표시되어 모호성이 증가된 부분들이 있다. 이는 이러한 모호성이 적도록 송신신호를 설계할 수 있는 능동형 레이다와 달리 임의의 FM 방송신호를 이용해야 하는 수동형 레이다이기 때문인 것으로 사료된다.
Ⅴ. 결 론
본 논문에서는 FM 방송신호를 이용한 PCL 장비의 산란신호 모의를 위한 방법을 서술하였다. Ray-tracing기법을 이용하여 송신소-수신기 사이의 주변환경에 의한 클러터 신호를 모의하고, 측정값과 유사함을 보였다. 또한 표적반사신호의 경우, 시간에 따라 변하는 도플러 효과를 정확히 고려하기 위해 IQ데이터를 시간영역에서 지연시키는 방법을 제안하였다. 또한 EM 시뮬레이션 툴을 이용하여 대표적인 표적의 RCS를 DB로 구축하고, 안테나의 어레이 매니폴드를 고려한 산란신호 모의를 수행하였다. 최종적으로 ARD 지도를 도시하여 표적의 반사신호가 적절히 모의된 것을 확인하였다. 추후 본 논문에서 설명된 산란신호 모의기법은 운용 S/W 및 모의신호 발생기와 통합되어 PCL 개발지원장비를 구성하게 된다. 그리고 이를 기반으로 PCL 주장비의 RF 성능 테스트, 클러터 및 기준신호 제거 테스트, 안테나 지향 테스트 등 성능 시험에 사용될 예정이다.
