REMI Lab강의 노트

Radar Signal Processing R12

12장 도래각 추정

ESPRIT: 훑지 않고 각 읽기

방향마다 훑지 않고 각을 곧바로 구할 수 있는가? 균일 선형 배열은 안테나 0,…,N−20,\dots,N-2와 1,…,N−11,\dots,N-1의 두 부배열로 나눌 수 있다. 두 번째 부배열은 첫 번째를 안테나 하나만큼 민 것이고, 방향 θ\theta의 스티어링 벡터는 안테나 하나만큼 밀면 ejπsin⁡θe^{j\pi\sin\theta}만큼 돈다.

a2(θ)=ejπsin⁡θ a1(θ)\mathbf a_2(\theta)=e^{j\pi\sin\theta}\,\mathbf a_1(\theta)

신호 부공간의 고유벡터 Es\mathbf E_s도 같은 두 부분 E1\mathbf E_1, E2\mathbf E_2로 나누면, 두 부분은 D×DD\times D 행렬 Ψ\boldsymbol\Psi 하나로 이어진다(E2=E1Ψ\mathbf E_2=\mathbf E_1\boldsymbol\Psi). Ψ\boldsymbol\Psi의 고윳값이 ejπsin⁡θie^{j\pi\sin\theta_i}이다. 최소제곱으로 Ψ=(E1HE1)−1E1HE2\boldsymbol\Psi=(\mathbf E_1^{\mathsf H}\mathbf E_1)^{-1}\mathbf E_1^{\mathsf H}\mathbf E_2를 구하고 그 고윳값의 각에서 θi\theta_i를 구한다. 이 방법을 ESPRIT(estimation of signal parameters via rotational invariance techniques)이라 한다. 각을 훑지 않으므로 계산이 가볍고, 추정값이 각도 격자에 묶이지 않는다.

MUSIC 그림과 같은 R̂(−5°와 10°, 안테나마다 SNR 5 dB, 스냅숏 64개)이고, ①의 초록 상자는 안테나 0…6, 점선 상자는 1…7의 부배열이다. ②의 속 빈 원은 실제 ejπ sin θ, 붉은 점은 R̂에서 구한 Ψ의 고윳값이다.

해 보기

  1. Ψ\boldsymbol\Psi의 고윳값이 실제 ejπsin⁡θe^{j\pi\sin\theta}와 얼마나 가까운지, 각으로 바꾸면 얼마인지 확인하라.
    관찰

    Ψ\boldsymbol\Psi의 고윳값이 실제 자리 가까이에 선다. 각으로 바꾸면 −5.14°-5.14°와 10.13°10.13°이다.

  2. 부배열을 안테나 하나만큼 밀 때 두 신호의 스티어링 벡터가 안테나 하나마다 얼마씩 도는지 확인하라.
    관찰

    −5° 신호는 안테나 하나마다 −15.7°-15.7°, 10° 신호는 31.3°31.3°씩 돈다.