MUSIC算法对信号DOA的应用-Read.DOC

MUSIC算法对信号DOA的应用波达方向（DOA）估计的基本问题就是确定同时处在空间某一区域内多个感兴趣的信号的空间位置（即多个信号到达阵列参考阵元的方向角）。最早的也是最经典的超分辨DOA估计方法是著名的MUSIC方法，MUSIC是多重信号分类（Multiple Signal Classification）的英文缩写。它是由R．O． Schmidt于1979年提出来的，由1986年重新发表的。MUSIC算法利用了信号子空间和噪声子空间的正交性，构造空间谱函数，通过谱峰搜索，检测信号的DOA．它是建立在以下假设基础上的：(1) 阵列形式为线性均匀阵，阵元间距不大于处理最高频率信号波长的二分之一；(2) 处理器的噪声为加性高斯分布，不同阵元间距噪声均为平稳随机过程，独立同分布，空间平稳（各阵元噪声方差相等）；(3) 空间信号为零均值平稳随机过程，它与阵元噪声相互独立；(4) 信号源数小于阵列元数，信号取样数大于阵列元数，信号源为窄带信号，即信号通过天线阵的时间远远小于信号带宽的倒数．5.2.1 MUSIC算法的基本原理 　 图5.1 均匀天线阵列如图5.1，M个天线阵元均匀直线排列，单元间距d为1/2个波长，布置成一个阵列天线。设有P(PM)个互不相关的窄带信号源平面波辐射到线阵上，信源方向分别为 。在第ｎ次采样时刻，得到的数据向量为 　 X（n）=AS(n)+U(n) 　 n=1，2，……N 　　 (5.1)式中X(n)= 为M个阵元输出；　　A= ， 式中　 ，T表示转置， 为载波波长，i=1，2，……，P； 为第i个平面波的复振幅；U(n)= ， 为零均值、方差为 的白噪声，且与信号源不相关；N为采样数。信号和噪声的协方差矩阵分别为　　　S= U= 接收信号的协方差（阵列输出信号协方差）　　 ，以上式中H为共轭转置 　　(5.2)因为 为MXM矩阵，所以能分解为M个特征值和特征向量，把这些特征值和特征向量用 ， (i=l，2，…，M)来表示，则 可表示为　　　 (5.3)这里，V是以 为元素的列矩阵， 是以 为元素的对角矩阵。从这个分析结果，有下面重要性质：[性质1] 各到达波是非相干(信号间相关系数不到l)，设各信号和噪声不相关，在 的特征值里，下面关系成立 　 　　 (5.4)即主要的特征值（信号特征值）个数和到达波束P相等，剩下的特征值（噪声特征值）的大小等于噪声功率。根据这个性质可以估计到达波的个数。进一步，按照特征值分布， 可分为信号功率和噪声功率之和　　 = = 　(5.5)　　 V=[ ]= 由于特征向量相互正交，则由下面第二个重要的性质。[性质2] 对应噪声特征值的特征向量（噪声特征向量）和各到达波的信号向量（信号特征向量）正交 。　 …M， i=1，…P. 　　 (5.6)于是，阵列的空间谱函数可表示为　 　　(5.7)式中分母是信号向量和噪声向量的内积。在性质2成立时的 分母是零， 有一尖峰。MUSIC算法就是通过寻找波峰来估计到达角的。通常把信号特征矢量覆盖的空间称为信号子空间(Signal Subspace)，噪声特征向量覆盖的空间称为噪声子空间（Noise Subspace)。把基于这个原理的估计到达波方向的方法称为部分空间法(Subspace Method)。MUSIC算法就是用信号或噪声子空间进行低秩信息的提取。5.2.2 MUSIC算法的实现　　MUSIC算法的实现步骤:1) 根据N个接收信号矢量得到阵列输出向量的协方差矩阵 　　　 (5.8) 对上面的协方差矩阵进行特征值分解　　　 (5.9)2) 然后按特征值的大小顺序，把与信号个数P相等的特征值和对应的特征向量看作信号子空间，把剩下的（M-P)个特征值和特征向量看作噪声部分空间。　　　 = (5.10)3) 使 变化，按照空间谱 来计算谱函数，通过寻找峰值来得到信号到达方向的估计值。以下给出基于MATLAB的MUSIC算法估计仿真：　　（1）从 入射的三个独立信号源，SNR分别为12dB，10dB，9dB。　　　图5.2 MUSIC算法的谱图 从谱图可以看出：在满足上面的假设前提下，MUSIC算法可精确估计出信号的DOA。尽管MUSIC算法在满足上述假设前提下可以精确估计信号的DOA，但它也有局限性：就是在低SNR和小样本的条件下无法分辨出空间相距比较近的信号。还有就是在现实当中，由于多径效应，接收到的信号一般是高相关信号，甚至是相关信号。当阵列接收到的是相干信号时，MUSIC算法就失去了其有效性，不再能估计出信号的DOA了 。　　（2）如下图，从 入射的三个信号源，SNR分别是20dB、10dB、12dB，其中，后面两个是相干信号。图5