Ⅰ. 서 론
여러 임무를 수행해야 하는 다기능 레이다에서 단일 펄스의 송수신만으로 표적의 거리와 속도 탐지가 가능하다면 짧은 드웰(dwell) 시간으로 표적의 거리와 속도를 탐지할 수 있기 때문에 효율적인 자원관리 측면에서 이점이 많다. 여러 개의 가간섭성(coherent) 펄스들을 이용한 속도 탐지는 도플러 분해능 및 클러터 제거 측면에서 유리하지만 짧은 드웰을 형성하기 힘들다. PCM(pulse code modulation) 파형을 이용해 단일 펄스로 거리와 속도를 탐지하는 방법은 있다[1]. 그러나 PCM 파형은 펄스 압축 후에 신호의 부엽 레벨 특성이 좋지 않고 정합 필터와 표적 신호의 도플러와의 차이에 대한 민감성이 높아 신호의 도플러 값이 커질수록 펄스 압축 후에 주/부엽 특성이 나빠지고, SNR(signal to noise ratio) 손실이 커진다[2]. LFM (linear frequency modulation) 파형은 펄스 압축 후에 신호의 부엽 레벨이 비교적 낮고, 표적 신호의 도플러에 의해 신호의 주/부엽 특성이 많이 변하지 않는 좋은 도플러 내성 때문에 높은 도플러 속도를 가지는 표적 탐지에도 유리하다[2]. 반면에 좋은 도플러 내성은 단일 펄스 LFM 파형으로 표적의 도플러 속도를 추정하는 것을 어렵게 한다. LFM 파형에서 속도를 추정하는 기존 연구들은 존재한다[3]~[5]. 참고문헌 [3]은 LFM 파형의 펄스 압축에서 나타나는 거리-도플러 결합(range-Doppler coupling) 현상을 이용한 것이며, 참고문헌 [4]는 정합 필터(matched filter)와 신호의 도플러 주파수와의 불일치에 의한 SNR 손실 현상을 이용한 것이다. 참고문헌 [3]과 [4]의 방법들은 레이다 반사 신호의 시간 축 영역에서 표적 신호가 있는 부분을 최적으로 추출할 수 있어야 최적으로 동작한다. 그러나 참고문헌 [3]과 [4]에는 최적의 알고리즘 동작을 위한 표적 신호 추출 방법이 기술되지 않아 실제 적용에 제한적이다. 참고문헌 [5]의 방법은 부분적으로 상승 스윕(sweep) 주파수와 하강 스윕 주파수를 가지는 LFM 파형(triangular FM 파형)으로 표적의 도플러를 추출한다. 참고문헌 [5]에서 사용하는 LFM 파형에 의해 펄스 압축 후 표적 신호의 피크(peak)가 표적의 도플러 크기에 따라 동일 크기의 두 개의 피크로 분리되어 나타난다. 시간 축에서 두 개의 피크로 분리된 폭의 크기로부터 표적의 도플러를 구할 수 있다. 그러나 참고문헌 [5]의 방법은 신호의 피크가 두 개로 분리됨으로 신호 피크의 크기에 손실이 발생하며, 여러 표적이 모여 있을 때 한 표적에 대한 한 쌍의 피크를 찾는 과정 등이 필요한 단점이 있다.
본 논문에서는 단위 펄스 LFM 파형으로 표적의 도플러를 추정하는 방법을 제안한다. 제안한 방법은 표적 신호의 도플러를 추정하기 위하여 LFM 파형의 거리-도플러 결합 현상과 정합 필터와 신호의 도플러 주파수와의 불일치에 의한 SNR 손실현상을 이용한다. 그리고 펄스 압축 외에 추가적인 두 번의 정합 필터링 처리로 표적의 도플러를 추정한다. 제안한 방법의 적용을 위한 표적 탐지 및 표적 신호 추출 방법 등의 과정을 기술하였으며, 시뮬레이션을 통해 제안 방법의 기능과 성능을 확인하였다.
Ⅱ. 도플러를 가진 LFM 반향 신호의 펄스 압축
레이다에서 송신하는 펄스 LFM 파형은 다음과 같다.
여기에서 xt(t)는 레이다 송신 신호, Tp는 펄스 폭, β는 스윕 주파수 대역폭, fc는 송신 주파수, At는 송신 신호 크기, 그리고 0≤t<Tp일 때에 Π(t/Tp)=1로 이외에는 Π(t/Tp)=0로 정의한다. R0 거리에 위치한 표적에 대한 레이다 반향 신호를 수신하여 송신 주파수 성분을 제거한(기저대역) 신호에서 LFM 반향 신호 부분을 추출하여 표현하면 다음과 같다.
여기에서 B(θr)=Arexp(jθr), Ar는 수신 신호 크기, θr는 신호 위상, Π((t−τ)/Tp)의 값은 0≤t-τ<Tp일 때에 1로 이외에는 0으로 정의, τ=2R0/C, C는 광속, fd는 수신 신호의 도플러 성분, 그리고 fd=2vd/λc로 표현되며, vd는 도플러 속도, λc는 송신 파장이다. 수신 신호 yR(t)를 기준신호 exp(jπ(β/Tp)(t-Tp/2)2)로 정합 필터링 처리하여 펄스 압축을 수행하는데, 그림 1은 β=40 MHz, Tp=200 μs인 LFM 파형을 이용한 경우에vd=6,400 m/s인 표적 반향 신호를 펄스 압축한 결과를 보여준다.
그림 1에서 보듯이 표적 속도가 0 m/s인 경우를 기준으로 했을 때 표적 속도 6,400 m/s인 경우에는 표적 전력 피크 위치가 거리 쪽으로 이용한 것을 볼 수 있다. 표적 피크의 크기도 좀 더 줄어든 것으로 볼 수 있다. 즉, LFM 파형의 펄스 압축에서 거리-도플러 결합에 의한 표적 피크의 거리 이동 현상과 정합 필터와 신호의 도플러 주파수와의 불일치에 의한 SNR 손실현상을 관찰할 수 있다. 본 논문은 LFM 펄스 압축에서 이런 거리-도플러 결합과 표 적 SNR 손실 현상을 이용하여 표적의 도플러 속도를 추출한다. 그림 2는 거리-도플러 결합 현상에 의한 표적 도플러와 표적 피크 이동 거리 관계를 나타낸다. 표적 도플러 fd에 의한 표적 피크 이동 거리 Δr은 Δr=−fdCTp/(2β)으로 표현할 수 있다[2]. 표적 피크의 이동 거리의 양으로 표적의 도플러를 탐지한다면 도플러 속도 탐지 해상도 Δv는 파형의 거리 해상도 및 레이다 신호의 샘플링 주파수 Fs에 의해 결정되며, 식 (3)과 같다(거리 해상도와 샘플링 주기를 동일하게 설정한 것으로 가정).
만약 Fs=40 MHz, Tp=200 μs, β=40 MHz, λc=0.03이면 Δv=74.95 m/s이며, 이런 속도 탐지 해상도는 너무 낮은 값이다(Δv는 표적 전력 피크의 한 거리 셀 이동에 해당하는 표적 도플러 속도를 나타낸다.). 식 (3)에서 속도 해상도는 β에 비해 Fs와 Tp의 값이 커진다면 좋아진다는 것을 알 수 있다. 본 논문에서는 LFM 펄스 압축 후에 표적 속도가 0 m/s인 표적 피크를 기준으로 표적 피크의 상대적 이동 거리로부터 표적의 속도를 추출하며, 보간(interpolation)을 통해 추출하는 속도의 해상도를 높인다.
Ⅲ. 단일 LFM 펄스로 표적 속도 탐지 처리 과정 설계
그림 3은 설계한 단일 LFM 펄스로 표적의 도플러 속도를 탐지하는 과정을 나타내었다. 그림 3의 첫 번째 과정은 식 (3)으로 결정되는 속도 탐지 해상도 간격으로 관심 속도 값의 최소치부터 최대치까지에 해당하는 도플러 속도를 포함하는 LFM 파형 신호들을 만들어 일렬로 연결하여 만든 새로운 원형 신호 C(n)을 만드는 과정이며 처음 한 번만 만들어 계속 사용하면 된다. 신호 C(n)은 관심 있는 모든 도플러 속도를 포함하는 신호이고, 이 신호를 이용하여 정합 필터링 처리를 하면 시간 축에서 연속한 LFM 신호들 중에서 표적의 속도에 해당하는 LFM 신호 지점에서 최대의 전력 피크가 나타날 것이다. 신호 C(n)은 시간 축에서 식 (4)와 같이 표현할 수 있다.
여기에서 K=ceil(Vmax/Δv), Vmax는 최대 관심 속도, ceil(x)는 x보다 크고 가까운 정수, f(k)=(Δv/3)2k/λc, α(k)=(k+3K)(Tp+3K/Fs)+(k+3K+1)(K/Fs)이다. 그림 3의 단일 LFM 파형 운영과 레이다 반사 신호 회득 과정은 식 (1)의 신호를 전송하여 식 (2)의 반향 신호를 포함하는 기저대역의 레이다 반사 신호를 획득하는 과정이다. 그 다음으로 일반적인 펄스 압축 등의 레이다 신호처리를 통해 표적 탐지하여 표적 거리를 추출한다[2]. 그리고 획득한 한 PRI(pulse repetition interval) 내의 레이다 반사 신호에서 표적 거리에 위치해 있는 표적 신호의 디지털 데이터 Y1(n)을 추출한다(1차 추출). 표적 디지털 신호 데이터 Y1(n)을 추출할 때는 표적의 도플러 속도에 의한 표적 전력 피크의 거리 셀(cell) 이동 현상(거리-도플러 결합 현상)을 고려하여 식 (5)와 같이 한 PRI의 데이터 R(n) 내의 l의 색인(index) 위치에서 데이터를 취한다.
여기에서 Nt는 한 PRI 샘플 내에서 탐지된 표적의 전력 피크 위치의 색인 값(한 PRI 내에서 탐지된 표적 위치), Np=TpFs는 펄스폭에 해당하는 데이터 샘플 수이다. 즉, 표적의 최대 속도 Vmax에 의한 −K 만큼의 표적 전력 피크의 셀 이동과 최저 속도 −Vmax에 의한 K 만큼의 셀 이동을 고려하여 표적의 신호 데이터를 모두 가져올 수 있도록 여유롭게 취하였다. 그림 3에서 다음 과정은 Y1(n)을 정합 필터로 해서 C(n)을 정합 필터링하여 C(n)을 펄스 압축 처리하는 과정이다. 그림 4는 C(n)을 펄스 압축 처리한 예이다. 그림 4에 대한 시뮬레이션 변수 값은 Fs=40 MHz, Tp=200 μs, β=40 MHz, λc=0.03, Vmax=6,400 m/s, 표적 속도 3,000 m/s, 그리고 표적 SNR은 43.9 dB이다. 그림 4에는 최대 피크 위치를 검은 점으로 찍어서 거리 색인 위치 X 가 표현되어 있다. 그러면 X 값을 이용하여 식 (6)과 같이 최대 피크 위치에 해당하는 LFM 파형의 도플러 속도를 확인(C(n)에 대한 식 (4)에서 도플러 속도의 색인 값 k를 이용)하여 1차로 표적 속도 v1을 추출할 수 있다.
여기서 round(x)는 x를 반올림한 정수값이다. 그림 4 에서 X=3,154,000, K=86, 그리고 Δv=74.95 m/s임을 이용하면, v1=2,998 m/s로 1차 표적 속도 값이 구해진다. 이 값은 실제 표적 속도 3,000 m/s에 비해 2 m/s의 오차가 있다. 1차 표적 속도 값 추출에서 어떤 경우에는 최대 37.48 m/s (Δv의 반)까지 측정오차가 날 수 있다. 이런 표적 속도 측정 오차를 줄이기 위해서 보간 기법을 이용할 것이다. 그림 3에서 1차 표적 속도 추출 값을 이용한 표적 디지털 신호의 2차 추출 과정은 한 PRI 내의 표적 피크 셀 위치와 표적 속도에 의한 거리 셀 이동을 고려하여 레이다 반사 신호에서 표적의 신호 데이터가 있는 영역 부분만을 정확하게 취하는 과정이다. 식 (7)은 식 (6)에서 구한 −(3K-Imax)의 값을 이용한 Y2(n) 획득을 표현했다.
이렇게 2차로 취한 표적 신호 데이터 Y2(n)을 보간 기법으로 M배 샘플 수를 증가시킴으로써 탐지 속도 해상도를 증가시킨다. 마찬가지로 M배 샘플 수를 증가시킨 도플러 속도 0인 LFM 파형 신호 데이터를 정합 필터로 사용하여 M배 샘플 수를 증가시킨 Y2(n)을 정합 필터링 처리를 하여 최종적으로 탐지 속도 해상도가 개선된 표적 속도를 추출하게 된다. 그림 5는 Y2(n) 추출과 보간법을 통해 개선된 탐지 속도 해상도로 속도를 추출하는 중간 과정의 예이다. 그림 5는 그림 4와 같은 시뮬레이션 변수 상황에서 M=80배로 보간 후에 펄스 압축하고, 그 펄스 압축 전력 값 결과의 데이터 순서를 L/2만큼 순환 이동(circularly shifting)하여 최대 전력 피크 위치를 표시한 그림이다. 식 (3)을 이용하면 개선된 속도 탐지 해상도 ΔvM=0.9369 m/s이다. 여기서 L은 그림 5와 같은 펄스 압축 결과의 총 데이터 길이를 나타낸다. 그림 5의 경우의 예는 L=1,048,576이며, 만약 표적 속도가 0이면 표적 전력 피크는 그림 5의 결과, 상황에서 L/2=524,288의 색인 위치에 존재하게 된다. 즉, 표적 전력 피크의 색인 값이 L/2에서 얼마나 이동하였는지의 값으로 표적의 속도를 측정할 수 있다. 그림 5의 경우는 L/2-521,086=3,202 셀만큼 전력 피크가 이동하였으며, 속도 v2=3,202×ΔυM인 2,999.8 m/s로 최종적인 표적 속도 값이 추출된다. 실제 표적 속도가 3,000 m/s이므로 1차로 탐지했던 속도 2,998 m/s에 비해 2차의 최종적인 속도 v2의 측정 오차가 상당히 줄어든 것을 확인할 수 있다.
Ⅳ. 시뮬레이션 결과 및 분석
이 장에서는 시뮬레이션을 통해 제안 방법의 기능과 성능을 확인한다. 먼저 그림 6은 다양한 표적 속도에서 제안한 방법으로 표적을 탐지한 결과를 보여준다. 시뮬레이션 조건은 Fs=40 MHz, 200 μs, β=40 MHz, λc= 0.03, Vmax=6,400 m/s, M=80이다. 그리고 표적 SNR 47.3 dB 및 0 m/s에서 6,400 m/s까지 37 m/s 간격으로 총 173가지의 표적 속도에서 표적 속도를 탐지해 보았다. 그림 6의 x축은 실제 표적 속도 vt를 나타내며, y축은 제안 방법으로 추정한 속도 ve와 실제 표적 속도 vt와의 차이 값의 절대치로서 속도 값 추정 오차를 나타낸다. 그림 6에서 보듯이 1차 속도 값 추정에서보다 2차 속도 값 추정에서 추정 오차가 상당히 줄어든 것을 볼 수 있으며, 2차의 최종적인 속도 값 추정에서 추정 오차의 절대치는 1.14 m/s 이하를 유지하는 것을 볼 수 있다(그림 6에서 ΔvM=0.9369 m/s임).
그림 6에서 1차 속도 추정치의 평균제곱근오차는 21.81 m/s이며, 최종적인(2차) 속도 추정치의 평균제곱근오차는 0.51 m/s이다.
그림 7은 SNR에 따른 표적 속도 탐지 성능을 나타낸다. 그림 7의 시뮬레이션 조건은 표적의 속도와 SNR을 제외하면 그림 6과 같다. 그림 7에서 실제 표적 속도는 3,000 m/s이며, x축은 표적의 SNR, y축은 제안 방법으로 추정한 속도 ve와 실제 표적 속도 vt와의 차이 값의 절대치로서 속도 값 추정 오차를 나타낸다. 그림 7(a)는 표적 SNR 10 dB에서 55 dB까지 일 때의 표적 속도 추정 오차를 나타내며, 추정 오차가 적으면서 계속 안정적으로 추출되기 시작하는 지점을 검은 점으로 표시하고, 그 지점의 SNR과 추정 오차를 표시하였다. 그림 7(a)의 점선 사 각형 영역을 확대하여 그림 7(b)에 나타내었다. 그림 7에 서 보면 표적 SNR이 25 dB 이하인 지점은 표적 속도 추정 오차가 크며, 표적 속도 추정이 안정적이지 않다. 표적 SNR이 25 dB 이상이 되면 표적 속도 추정 오차가 (어떤 한 지점에서 오차가 1,800 m/s인 것을 제외하면) 450 m/s이하로 계속 떨어진다. 표적 SNR이 31 dB 이상이 되면 표적 속도 추정 오차가 (어떤 두 지점에서 오차가 400 m/s 정도인 것을 제외하면) 73 m/s 이하로 계속 떨어진다. 표적 SNR이 34 dB 이상이 되면 표적 속도 추정 오차가 73 m/s 이하를 계속 유지한다. 그리고 표적 SNR이 39 dB 이상이 되면 표적 속도 추정 오차가 1.14 m/s 이하를 계속 유지한다. 제안한 방법의 특성상 일단 1차 속도 추정이 잘 되어 1차 속도 추정 오차가 Δv=74.95 m/s보다 작으면 2차 속도 추정 오차는 1.14 m/s 이하가 됨을 확인하였다. 그림 7에서 보듯이 안정적인 표적 속도 탐지 성능을 보장하기 위해서는 상당히 높은 SNR을 요구한다. 그러나 높은 SNR (39 dB 이상)이 보장되면 표적 속도 탐지 해상도는 아주 좋다.
그림 8은 표적이 밀집되어 표적 신호들이 서로 영향을 미칠 때의 상황을 시뮬레이션하여 표적 속도 탐지 성능을 확인한 결과를 나타낸다. 그림 8(a)의 시뮬레이션 조건은 표적 SNR, 표적 속도 및 표적의 개수가 3개이며, 서로 밀집되어 있다는 것 외에는 그림 6과 같다. 그림 8(a)는 표적 3개가 파형 펄스폭의 거리 영역 안에 밀집해 존재하여 표적 신호의 한 펄스폭 안에 3개의 표적 신호가 섞여 들어 있는 경우이다. 3개 표적의 속도는 각각 2,250 m/s, 3.000m/s, 3,750 m/s이며, 표적의 SNR은 각각 43.9 dB, 42.0 dB, 36.0 dB이다. 그림 8(a)의 결과는 속도가 2,250 m/s이고, SNR이 43.9 dB인 표적을 탐지하고, 표적 속도를 추출한 결과를 그림 5와 같은 형식으로 나타낸 것이다. 탐지된 표적 속도는 (L/2-5,218,206)ΔυM으로서 2,325.4 m/s이다. 표적의 SNR이 충분하지만 밀집된 다른 표적 신호에 영향을 받아서 표적 속도 추정 오차가 커졌다(75.4 m/s). 표적 속도를 추출하는데 쓰이는 펄스폭 영역의 표적 신호에 SNR이 1.9 dB 밖에 차이가 나지 않는 다른 표적 신호가 같이 섞여 있다. 그림 8(b)의 시뮬레이션 조건은 표적 SNR을 제외하고는 그림 8(a)와 같다. 그림 8(b)에서 표적 SNR은 각각 50.0 dB, 42.0 dB, 36.0 dB이다. 그림 8(b)의 결과는 속도가 2,250 m/s이고, SNR이 50.0 dB인 표적을 탐지 및 속도를 추출한 결과를 그림 8(a)처럼 나타낸 것으로 탐지된 표적 속도는 2,249.5 m/s이다. 펄스폭 영역 안에 다른 표적 신호가 섞여 있지만 약 8 dB 이상의 차이로 탐지하고자 하는 표적 신호의 전력이 다른 신호의 전력보다 높으니 제안한 방법의 정상적인 성능으로 표적 속도가 탐지되었다.
Ⅴ. 결 론
레이다에서 단일 펄스 LFM 파형을 운영하여 획득한 단일 펄스의 LFM 레이다 반사 신호에서 표적의 속도를 탐지하는 방법을 제시하였다. 제안한 방법으로 표적 속도를 탐지하기 위한 처리 과정의 흐름을 보였으며, 처리 흐름에서 표적 거리 탐지를 위한 정합 필터링 처리 외에 추가적인 2번의 정합 필터링 처리를 하여 표적 속도를 추출하였다. 제안 방법은 표적이 저속 영역에 있든지 고속 영역에 있든지 잘 탐지했으며, 다만 속도 탐지 오차는 본 논문에서 정의한 속도 탐지 해상도에 의해 좌우되었다. 본 논문에서 정의한 속도 탐지 해상도는 펄스 폭 등의 파형 설계 및 샘플링 주파수에 의해 좌우되며, 파형 설계를 통해 속도 탐지 해상도를 높이는 것은 제한이 따르기 때문에 본 논문에서는 LFM 파형의 거리-도플러 결합 현상을 기반으로 보간 기법을 이용해 표적 속도 탐지 해상도를 높였다. 시뮬레이션을 통해 제안한 방법의 표적 속도 탐지 기능과 속도 추정 성능을 확인하였다. 제안한 방법이 최상의 성능을 보이려면 다소 높은 SNR(39 dB 이상)을 요구하였다. 시뮬레이션을 통해 제안한 방법의 1차 표적 속도 추출에서 파형에 의해 결정되는 속도 탐지 해상도 이하의 오차로 표적 속도를 안정적으로 탐지한다면 최종적인 2차 표적 속도 추출에서 적은 오차로 표적 속도를 탐지함을 확인하였다. 제안한 알고리즘의 개선을 위해 낮은 SNR에서도 1차 표적 속도 추출이 안정적으로 이루어지는 방법이 연구되어야 할 것이다. 제안한 방법은 짧은 드웰 시간으로 표적의 거리와 속도를 탐지할 수 있는 파형의 설계를 가능하게 하기 때문에 제한된 시간에 많은 기능을 수행하는 다기능 레이다에 표적 SNR에 따라 알고리즘을 단계적으로 적용하면 이용될 수 있다. 물론 응용 분야마다 단계적 알고리즘 적용 방법이 연구되어야 한다. 그리고 제안한 알고리즘의 현 상황에서 적용하기 적합한 대상 표적은 장거리에서 움직이는 RCS(radar cross sections)가 높은 표적(예를 들면, 장거리 탄도탄)이다. 장거리 표적은 넓은 송신 펄스폭을 원하고, 넓은 송신 펄스폭은 높은 SNR과 좋은 속도 탐지 해상도를 얻게 한다.