自适应滤波实验报告

《自适应滤波实验报告》由会员分享，可在线阅读，更多相关《自适应滤波实验报告（11页珍藏版）》请在装配图网上搜索。

1、LMS自适应滤波实验报告姓名： 学号：_ 日期：实验内容：利用自适应滤波法研究从宽带信号中提取单频信号的方法。设 x t st Acos2 f1tBcos2 f2t，st 是宽带信号，A，B，f1，f2，任选（1）要求提取两个单频信号；（2） 设f2 t f，要求提取单频信号 cos 2 f?t，研究f的大小对提取单频信 号的影响。1. 自适应滤波器原理自适应滤波器理论是现代信号处理技术的重要组成部分，它对复杂信号的处理具有独特的功能。自适应滤波器在信号处理中属于随机信号处理的范畴。在一些信号和噪声特性无法预知或他们是随时间变化的情况下，自适应滤波器通过自适应滤波算法调整滤波器系数，使得滤波器

2、的特性随信号和噪声的变化，以达到最优滤波的效果，解决了固定全系数的维纳滤器和卡尔曼滤波器的不足。（1）自适应横向滤波器所谓自适应滤波，就是利用前一时刻已获得的滤波器参数等结果，自动调节现时刻的滤波器参数，以适应信号和噪声未知或随时间变化的统计特性，从而实现最优滤波。 自适应滤波器由两个部分组成：滤波器结构和调节滤波器系数的自适应算法。自适应滤波器的特点是自动调节自身的冲激响应， 达到最优滤波，此算法适用于平稳和非平稳随机信号，并且不要求知道信号和噪声的统计特性。一个单输入的横向自适应滤波器的原理框图如图所示:实际上这种单输入系统就是一个fir网络结构，其输出 y n用滤波器单位脉冲响应表示成下

3、式：N 1y nw m x n mm 0这里w n称为滤波器单位脉冲响应，令：i m i,Wj w i 1 ,x xn i 1 , n用j表示,上式可以写成NyjWiXji 1这里Wi也称为滤波器加权系数。用上面公式表示其输出，适用于自适应线性组合器，也适用于FIR滤波器。将上式表示成矩阵形式：式中TTyjXjW w XjTTWw1,w?,.,wn , X jX1j,X2j ,.,XNj误差信号表示为Tejdj yj dj W Xj(2) 最小均方(LMS算法Widrow等人提出的最小均方算法，是用梯度的估计值代替梯度的精确值，这种算法简 单易行，因此获得了广泛的应用。LMS算法的梯度估计值用

4、一条样本曲线进行计算，公式如下:2jej2ej2ej2ejW1w2Wn因为ejdj WTXj所以2ejW12ejw22 TejWnXjj 2eXjWji Wj 2 e XjFIR滤波器中的第i个权系数的计算公式为Wj i,i Wj,i2 ejXj,i i 1,2,3,，N2. 编程思想在本题目中，要求在受宽带信号干扰的输入信号中分别提取出两个单频信号，因此，可采用KLS自适应滤波方法提取。在仿真中设置如下：（1）单频信号幅度A=2, B=3,相位3（2）调整步长分别设为 mu仁0.0001, mu2=0.00001 ；（3）滤波器阶数M=50;（4） 宽带信号设为均值为0，方差为1的白噪声；（

5、5） 提取单频信号 F1 Acos2 f1t，设置期望信号为理想信号，即为F1 ；（6） 提取单频信号F2Bcos2 f?t，设置期望信号为理想信号，即为F2 ；然后，根据LMS算法，更新滤波器系数 W3. 源代码%自适应滤波LMS算法%Delte 为 f2-f1fun ction Adaptive_Filter_LMS（f1,Delte）mu1=0.0001;mu2=0.00001; %步长Ts=0.0001; %采样间隔N=3000;%言号长度M=50;%阶数A=2;B=3;f2=f1+Delte;Phi=pi/3;%频率幅度相位设置s=1*ra ndn (1,N);%宽带信号k=1:N;

6、%两单频信号F1=A*cos(2*pi*f1*k*Ts+Phi);F2=B*cos(2*pi*f2*k*Ts);F=F1+F2;X=s+F1+F2;Y1=zeros(1,N);Y2=zeros(1,N);W1= zeros(1,M);W2=zeros(1,M);e1=zeros(1,N);e2=zeros(1,N);%提取信号F1for n=M:Nx=X( n:-1: n-M+1);Y1( n)=W1*x;e1( n)=F1( n)-Y1( n);W1=W1+2*mu1.*e1( n).*x;end%提取信号F2for n=M:Nx=X( n:-1: n-M+1);Y2( n)=W2*x;e2

7、( n)=F2( n)-Y2( n);W2=W2+2*mu2.*e2( n) .*x;endfigure(1)subplot(4,1,1);plot(F1(1:N);ylim(-2 2);title(%滤波器输入值%滤波器输出值%误差值%系数调整%滤波器输入值%滤波器输出值%误差值%系数调整单频信号波形F1);subplot(4,1,2);单频信号波形F2);plot(F2(1:N);ylim(-3 3);title(subplot(4,1,3);plot(F(1:N); title(无噪声两单频信号波形);subplot(4,1,4);plot(X(1:N);title(受宽带干扰的输入信号

8、波形);figure (2)subplot(4,1,1);plot(Y1(1:N); title(提取的单频信号波形F1);subplot(4,1,2);plot(e1(1:N); ti tle(提取单频信号波形F1误差);subplot(4,1,3);plot(Y2(1:N); title(提取的单频信号波形F2);subplot(4,1,4);plot(e2(1:N); ti tle(提取单频信号波形F2误差);4. 实验结果及分析(1) f1 100Hz, f2 200Hz，提取单频信号单频信号波形F1单频信号波形F2V I-ttti-20-211 J-J押y05001000150020

9、0025003000提取的单频信号波形F2提取单频信号波形F2误差分析：从图中可以看出，该仿真提取出了两单频信号，随着对滤波器权系数的更新，误差慢慢减小。通过比较可以看出，提取结果存在误差，且F1提取效果较差，F2提取效果较好。(2)改变频率差值，提取单频信号fi 100Hz, f 200Hz单频信号波形F1单频信号波形F2050010001500200025003000无噪声两单频信号波形受宽带干扰的输入信号波形提取的单频信号波形F10-5500250030002050025003000100015002000提取的单频信号波形F20-2-5 L0500100015002000提取单频信号波

10、形F2误差25003000-5 L0Ir50010001500200025003000100015002000提取单频信号波形F1误差图4f1100Hz, f 50Hz-2 L0单频信号波形F1Il / J5002500L100015002000单频信号波形F2300020-2-50100015002000无噪声两单频信号波形250030005000500A.io100015002000受宽带干扰的输入信号波形25003000提取的单频信号波形F15 -c匚1t0 - sW/eTWJwWz也沁卫打人叮丫 -5 ELLI:050010001500200025003000提取单频信号波形F1误差5 LILIII:0 -:*rArAr 讨了 .、，严、；吩.L-5 IIr050010001500200025003000提取的单频信号波形F2提取单频信号波形F2误差505001000150020002500-53000分析：当两信号频差较大时，如本次仿真中f 200Hz，滤波器较好的提取出两信号，但频差较小时，滤波器的效果很差，不能很好的提取出两单频信号。