资源描述：
《music、esprit、mvdr算法的谱估计.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、课程（论文）题目：MUSIC、ESPRIT、MVDR算法的谱估计内容：1算法原理1.1MUSIC算法MUSIC算法利用信号子空间和噪声子空间的正交性，构造空间谱函数，通过谱峰搜索，估计信号频率。由，且矩阵可逆得，。又由于矩阵为正定的对角矩阵，方程两边可再同时左乘，推出，，。这就表明，信号频率向量与噪声子空间的特征向量正交。信号角频率的估计可以由扫描函数的K个峰值位置确定。1.2ESPRIT算法ESPRIT算法即基于旋转不变技术的信号参数估计。连续M个时刻的观测值可表示为向量形式。定义随机过程，且向量和矩阵分别为，，则。 向量的自相关矩阵为，向量和的互
2、相关矩阵为。对进行特征分解，找到的最小特征值。定义矩阵：，可以通过求解方程式来求得到矩阵的广义特征值。当时，矩阵是奇异的；而时，是满秩的。矩阵对的广义特征值恰为，这些根的相位即为信号的频率估计。1.3MVDR算法MVDR算法即最小方差无失真响应算法，是有别于经典功率谱估计和参数模型估计的另一类信号频率估计方法。定义向量，假定信号通过一个M抽头的FIR滤波器，则输出信号为：的功率可以表示为。 为求得滤波器的系数，需要满足在对给定的某一频率处，无失真地通过，且最小。此时， 则。并不是真正意义上的功率谱，但它描述了信号真正谱的相对强度，可
3、以由此估计信号频率。2算法实现采用空间谱估计的典型代表MUSIC算法，ESPRIT算法和MVDR算法，对含有高斯白噪声的复正弦信号进行频率估计。选取的信号数，阵元数，采样数0。待检测信号的归一化频率为，仿真的待检测信号为：假设都是零均值，方差为1的白噪声，采样数为N，且彼此之间相互独立，则阵列响应矩阵为假定噪声为零均值，方差为1的高斯白噪声，采样数为N，则待检测信号引入的噪声为于是仿真信号为：。2.1MUSIC算法利用MUSIC算法进行频率估计时，首先求出仿真信号的自相关矩阵，然后对自相关矩阵进行特征值分解得到特征值空间和特征向量空间，然后对特征值空间进行升序排序
4、，取较小的个特征值所对应的特征向量组成噪声子空间G，然后求得,MUSIC算法程序段如下所示：R=X*X'/N;[UD]=eig(R);[d,index]=sort(diag(D));U=U(:,index);fori=-N:Na=exp(-j*2*pi*[0:M-1]'*(i*0.5/N));Pmusic(i+1+N)=abs(1/(a'*U(:,1:end-p)*U(:,1:end-p)'*a));endplot(omg,10*log10(Pmusic/max(Pmusic)));%画出MUSIC算法图形2.1ESPRIT算法ESPRIT算法要先构造相关矩阵和，
5、然后对进行特征值分解得到最小特征值即为噪声的方差，通过对矩阵对进行广义特征值分解，最接近单位圆的K个特征值相位即为信号的频率估计。ESPRIT算法程序段如下所示：S1=U(1:end-1,8:10);S2=U(2:end,8:10);S12=[S1S2];[UuDd]=eig((S12'*S12));[dd,ind]=sort(diag(Dd));Uu=Uu(:,ind);dt=-Uu(1:3,1:3)*inv(Uu(4:6,1:3));dd=eig(dt);t_esprit=asin(-angle(dd)/(2*pi));plot(t_esprit,[000],
6、'*','color','red')%画出ESPRIT算法图形，红色表示2.3MVDR算法MVDR算法要先求自相关矩阵，然后将峰搜索矩阵带入最小方差谱估计公式通过观察谱峰得到信号的频率估计值。MVDR算法程序段如下所示：fori=-N:Na=exp(-j*2*pi*[0:M-1]'*(i*0.5/N));Pmvdr(i+N+1)=abs(1/(a'*inv(R)*a));endplot(omg,10*log10(Pmvdr/max(Pmvdr)),'--','color','black')2.4程序流程图各算法程序流程图如下所示：图2.4.1MUSIC算法图2.4
7、.2ESPRIT算法图2.4.3MVDR算法2仿真结果使用MATLAB软件进行仿真，得出MUSIC算法、ESPRIT算法、MVDR算法的仿真图形分别如图3.1、图3.2、图3.3所示。为了将这三种算法的谱估计效果更加直观地比较，将三种算法的仿真结果图形放在一张图中进行比较，如图3.4所示。另外，MUSIC算法得出的频率估计为，运算时间time=1.415；ESPRIT算法的频率估计为，运算时间time=0.0251；MVDR算法的频率估计为，运算时间time=0.0910。图3.1MUSIC仿真结果图3.2ESPRIT仿真结果图3.3MVDR仿真结果图3.4各种算
8、法仿真比较
