第4章 多重信号分类算法
4.1 引言
多重信号分类(MUSIC)算法基本思想:将任意阵列输出数据的协方差矩阵进行特征分解,从而得到与信号分量相对应的信号子空间和与信号分量相正交的噪声子空间,利用这两个子空间的正交性估计信号的参数(入射方向、极化信息、信号强度等)。
MUSIC算法在特定的条件下具有很高的分辨力、估计精度及稳定性。
4.2 MUSIC算法
4.2.1 经典MUSIC算法
窄带远场信号的DOA数学模型为
X ( t ) = A ( θ ) S ( t ) + N ( t ) \pmb{X}(t)=\pmb{A}(\theta)\pmb{S} (t)+ \pmb{N} (t) XXX(t)=AAA(θ)SSS(t)+NNN(t)
阵列数据的协方差矩阵为
R = E [ X X H ] = A E [ S S H ] A H + σ 2 I = A R S A H + σ 2 I \begin{aligned} \pmb{R} &=E[\pmb{X}\pmb{X}^H]\\ &=\pmb{A}E[\pmb{S}\pmb{S}^H]\pmb{A}^H+\sigma^2\pmb{I}\\ &= \pmb{A}\pmb{R}_S\pmb{A}^H+ \sigma^2\pmb{I} \end{aligned} RRR=E[XXXXXXH]=AAAE[SSSSSSH]AAAH+σ2III=AAARRRSAAAH+σ2III
由于信号与噪声相互独立,数据协方差矩阵可以分解为与信号、噪声相关的两部分,其中 R S \pmb{R}_S RRRS是信号的协方差矩阵, A R S A H \pmb{A}\pmb{R}_S\pmb{A}^H AAARRRSAAAH是信号部分。
对 R \pmb{R} RRR进行特征分解
R = U S Σ S U S H + U N Σ N U N H \pmb{R}=\pmb{U}_S\Sigma_S\pmb{U}_S^H+\pmb{U}_N\Sigma_N\pmb{U}_N ^H RRR=UUUSΣSUUUSH+UUUNΣNUUUNH
式中, U S \pmb{U}_S UUUS是由大特征值对应的特征矢量张成的子空间,即信号子空间, U N \pmb{U}_N UUUN是由小特征值对应的特征矢量张成的子空间,即噪声子空间。
理想条件下数据空间中的信号子空间与噪声子空间是相互正交的,即信号子空间中的导向矢量也与噪声子空间正交
a H ( θ ) U N = 0 \pmb{a}^H(\theta)\pmb{U}_N=0 aaaH(θ)UUUN=0
经典的MUSIC算法正是基于上述性质提出,考虑到实际接收数据矩阵是有限长的,即数据协方差矩阵的最大似然估计为
R ^ = 1 L ∑ i = 1 L X X H \hat{\pmb{R}}=\frac{1}{L}\sum^{L}_{i=1}\pmb{X}\pmb{X}^H RRR^=L1i=1∑LXXXXXXH
对 R ^ \hat{\pmb{R}} RRR^进行特征分解可以计算得到噪声子空间特征矢量 U ^ N \hat{\pmb{U}}_N UUU^N。由于存在噪声, a ( θ ) \pmb{a}(\theta) aaa(θ)与 U ^ N \hat{\pmb{U}}_N UUU^N并不能完全正交,因此,实际上求DOA是以最小优化搜索实现的,即
θ M U S I C = arg min θ a H ( θ ) U ^ N U ^ N H a ( θ ) \theta_{\rm MUSIC}=\argmin_\theta\pmb{a}^H(\theta)\hat{\pmb{U}}_N\hat{\pmb{U}}_N^H\pmb{a}(\theta) θMUSIC=θargminaaaH(θ)UUU^NUUU^NHaaa(θ)
MUSIC算法的谱估计公式为
P M U S I C = 1 a H ( θ ) U ^ N U ^ N H a ( θ ) P_{\rm MUSIC}=\frac{1}{\pmb{a}^H(\theta)\hat{\pmb{U}}_N\hat{\pmb{U}}_N^H\pmb{a}(\theta)} PMUSIC=aaaH(θ)UUU^NUUU^NHaaa(θ)1
MUSIC算法步骤:
MUSIC算法 |
---|
1.由阵列的接收数据得到数据协方差矩阵 R ^ \hat{\pmb{R}} RRR^; |
2.对 R ^ \hat{\pmb{R}} RRR^进行特征分解; |
3.由 R ^ \hat{\pmb{R}} RRR^的特征值进行信号源数判断; |
4.确定信号子空间 U ^ S \hat{\pmb{U}}_S UUU^S与噪声子空间 U ^ N \hat{\pmb{U}}_N UUU^N; |
5.根据信号参数范围由谱估计公式进行谱峰搜索; |
6.找出极大值点对应的角度就是信号入射方向。 |
注:
(1)非理想情况下得到的协方差矩阵特征值满足:
λ 1 > λ 2 > ⋯ > λ N > λ N + 1 > ⋯ > λ M > \lambda_1>\lambda_2>\cdots>\lambda_N>\lambda_N+1>\cdots>\lambda_M> λ1>λ2>⋯>λN>λN+1>⋯>λM>
但不满足
λ 1 ≥ λ 2 ≥ ⋯ ≥ λ N ≥ λ N + 1 = ⋯ = λ M = σ 2 \lambda_1\geq\lambda_2\geq\cdots\geq\lambda_N\geq\lambda_{N+1}=\cdots=\lambda_M=\sigma^2 λ1≥λ