华中科技大学ADSP仿真作业(LMS-LSL--PHD-MUSIC)Word文档格式.doc

1、2.1题目62.2数据模型及参数选取62.3算法模型72.3.1PHD算法简介72.3.2MUSIC算法简介82.4实验仿真92.4.1仿真结果92.4.2结果分析102.5实验小结11附录：源代码12附录A：实验一12附录B：实验二14实验一1.1 题目Implement LSL algorithm and LMS algorithm based on figure 3.30（P92） and figure 3.31（P93） .Model and parameters see page 91.课本91页模型及参数为：序列由2阶自回归模型产生：，是均值为0、方差为1的高斯白噪声序列。用LMS

2、算法和LSL算法来估计模型参数、，对两种算法的性能进行比较，并分析LSL算法中初始预测误差剩余对收敛性能的影响。1.2 数据模型及参数选取1）高斯白噪声可以用中的函数产生均值为0、方差为1的标准正态分布随机矩阵来实现，用模型产生信号，激励源是产生的高斯白噪声；2）信号点数N=500，用500个信号来估计滤波器系数；3）依据题意，LMS算法中选取，并取=0.1,1, 10来观察初始预测误差剩余对LSL收敛性能的影响。1.3 算法模型1.3.1 自适应滤波器自适应滤波器由参数可调的数字滤波器和自适应算法两部分组成，如Error! Reference source not found.所示。输入信号

3、通过参数可调的数字滤波器后产生输出信号，将其与参考信号进行比较，形成误差信号。通过某种自适应算法对滤波器参数进行调整，最终使的均方值最小。此时得到的是对的最佳估计。1.31自适应滤波器原理1.3.2 LMS算法简介LMS算法采用平方误差最小的原则代替最小均方误差最小的原则，使用LMS算法来得到权系数的公式如下：（1.31）其中：（1.32）式中，为时刻自适应滤波器的权矢量，为自适应滤波器的阶数；为时刻自适应滤波器的参考输入矢量，由最近个信号采样值构成，；是期望的输出值；为自适应滤波器的参考响应与输出相应之差，即误差信号；是控制自适应速度与稳定性的增益常数，又叫收敛因子或步长因子。1.3.3 L

4、SL算法简介LSL算法步骤如下：1） 初始化这里是前向和后向预测误差剩余的初始值。如果这个值没有给出，可以任意选择。2） 迭代计算（按时间）；3）迭代计算（按阶）同阶的嵌套着按时间进行迭代计算。各阶前向和后向反射系数可分别由以下二式计算：M是给定的滤波器的阶。采用LSL自适应算法，利用给出的数据可计算出格型滤波器的1阶和2阶前向和后向反射系数、和，然后由这些反射系数计算出估计、。经推导得 （1.33） （1.34）1.4 实验仿真1.4.1 仿真结果图1.41 AR（2）所产生信号图1.42自适应权系数、的收敛轨迹 图1.43初始预测误差剩余对LSL收敛性能的影响1.4.2 结果分析1）由图1

5、.42可看出，一是LSL算法和LMS算法计算得到的都收敛于1.558，都收敛于-0.81，但是LSL算法很明显地比LMS算法收敛得更快，二是LMS算法得到的、的收敛轨迹有所波动，而LSL算法所得到收敛轨迹更加平稳；2）由图1.43，在=0.1和1两种情况下，很快就收敛到，而当选为10时，收敛稍许慢了一些，但三种情况的收敛速度都要比LMS算法快的多。1.5 实验小结通过本次实验，对LMS算法和LSL算法的原理有了更深刻的理解，对影响LMS和LSL算法性能的一些参数也有了一定的认识，懂得了如何通过调节参数来获得想要的性能。17实验二22.1 题目给定，是方差为的复白噪声。信噪比定义为1）用PHD方

6、法估计两个频率，使用M+1=3阶自相关矩阵；2）MUSIC方法估计频率，选取N=12，M=2。2.2 数据模型及参数选取1）依题意，可通过提供的函数可对复正弦信号加高斯白噪声得到。2）PHD算法取样自相关点数正弦信号个数自相关矩阵阶数复白噪声功率 = 0.01信噪比待估计频率3）MUSIC算法复白噪声功率 = 0.04信噪比SNR =13.9794dB2.3 算法模型2.3.1 PHD算法简介PHD（Pisarenko Harmonic Decomposition）即谐波分解方法。在的特殊情况下，是矩阵，是M阶对角阵，信号的自相关矩阵可表示为 （2.31） 设是中的任一矢量，表示由映射 （2.

7、32）得到的一组矢量，令 （2.33）于是 （2.34）由此看出，的秩最多是M，若可任意选择，则秩恰为M。由于是矩阵，这就是要求他在中。利用式（2.33），只需证明是满秩的或=是线性独立的。但 （2.35）它不能是零，除非所有是零（假设所有是非零且所有的频率不同）。这就证明了的秩是且有唯一的零特征值。由于 （2.36）与有相同的特征矢量，而的特征值是，所以的特征值是。现考虑对于最小特征值的特征矢量，有 （2.37）由式（2.36），得到 （2.38）左乘，得到 （2.39）由于P是正定的，于是得出 （2.310）, （2.311）上式可写成展开形式, （2.312）这就是说，如果给定，就能够求

8、出多项式, （2.313）位于单位圆上的零点，这些零点的角度就是正弦波的频率。2.3.2 MUSIC算法简介MUSIC（Multiple Signal Classification）即多信号分类算法，这种算法是以下列谱估计 （2.314）的峰作为正弦波频率的估计的。式中，表示频率，表示信号矢量。理论上，当时，。由于存在估计误差，故在正弦波频率上或者附近将有一个峰，与式（2.314）的M个最大峰对应的频率便是正弦波的频率估计。式（2.314）的最大值也可以利用信号子空间的特征矢量来进行计算。为此，注意到该式分母最小时最大，但分母可表示成 （2.315）该式利用了。式（2.314）的最大化等效于下

9、式 （2.316）的最大化。可见，无论是原先提出的噪声子空间特征矢量，还是式（2.315）的信号子空间矢量，都可以用来确定正弦波的频率估计。2.4 实验仿真2.4.1 仿真结果PHD算法中解方程结果（运行五次程序）：0.45170.45910.45030.45200.45490.52550.53690.52470.52570.5298表2.41解方程求得的频率估计值图2.41作图法估计频率值0.43660.45330.44590.45050.46200.51500.52650.52600.5403表2.42 PHD作图法捕捉到的频率值图2.42 MUSIC求估计频率值0.46940.47080

10、.47350.47250.47120.51180.51230.51040.5134表2.43 MUSIC方法捕捉到的频率值2.4.2 结果分析1）PHD算法中，解方程的方法，求解方程（2.412）位于单位圆上的零点，并计算的角度，然后把变换到，再将角频率除以，就可以得到估计的频率值。由表 2.1可以看出，估计出来的两个频率值与实际值0.5和0.52存在一定的误差。PHD的另一种估计频率的方法是，计算式（2.410）的值，在两个频率附近，值应该为0，但实际上由于估计误差，不会精确地等于0，如果把计算的结果倒数，那么在估计的频率附近，将会出现尖峰，捕捉尖峰所对应的横轴值，就可以得到估计的频率值。由

11、图2.41、表2.44与表2.42可以看出，两个峰值频点与实际都存在一定的误差，但都在可以接受的范围之内。2）图2.42为MUSIC 方法估计频率，其方法与PHD的第二种方法类似，由表2.45可看出，该方法估计的频率值与标准值0.5和0.52之间也存在误差。2.5 实验小结通过本次实验，对频率估计方法有了一定的了解，掌握了PHD和MUSIC算法估计频率方法。源代码clear;close all;% 由AR（2）模型生成输入信号% 初始化参数a1 = 1.558;a2 = -0.81; N = 500; %信号点数 %信号及白噪声信号序列的初始化x = zeros（1,N）; %信号的初始化Noise = randn（1,N）; %白噪声的初始化，均值为0，方差为1x（1） = Noise（1）;x（2） = a1*x（1）+Noise（2）;for n = 3:N x（n） = a1*x（n-1）+a2*x（n-2）+Noise（n）;end% LMS算法估计a1、a2的收敛曲线 u = 0.005;