用Music算法估计非相干信号的仿真.docxVIP

用Music算法估计非相干信号的仿真Music算法的原理1.1Music算法的基本步骤考虑均匀线阵由M个无方向性阵元构成 (如图 1 所示), 阵元间距为d(d不大于λ/ 2),设有N(NM)个窄带信号源平面波辐射到线阵上,波长为 λ , 信源方向分别为 。则第 k 次快拍得到的数据向量为: (1)图 1 M元均匀线阵列模型式中，X(k)=为M个阵元输出,K为快拍次数。,其中 T 表示转置,为第 i 个平面波的复振幅,为零均值、方差为的高斯白噪声, 假设信号和噪声以及各个阵元上的噪声互不相关, 则阵列输出向量的协方差矩阵为 （2）式中 是为信号相关矩阵,为 M阶单位矩阵。对式(2)所定义的协方差矩阵进行特征分解后, 可以得到 M个特征值: 及相应的特征矢量本原理可知, 方向矩阵张成的空间为。从而 MUSIC空间谱函数可定义为 1.2 Music算法的DOA估计流程图图2 Music算法测向估计流程图1.3计算机仿真采样点数为1024，SNR=15dB，均匀线阵，阵元数M=8，阵元间距d=/2，一维独立信源且信源数N=2，信号频率为：。Matlab代码如下：%%%%%%%% MUSIC算法 %%%%%%%%%%[U,V]=eig(R);UU=U(:,1:M-p);theta=-90:0.5:90;%%谱峰搜索for ii=1:length(theta) AA=zeros(1,le