基于LMS和RLS的自适应滤波器的应用仿真的设计与开发.docx

《基于LMS和RLS的自适应滤波器的应用仿真的设计与开发.docx》由会员分享，可在线阅读，更多相关《基于LMS和RLS的自适应滤波器的应用仿真的设计与开发.docx（15页珍藏版）》请在第一文库网上搜索。

1、湖南大学计算机与通信学院课程作业2题 目：基于LMS和RLS的自适应滤波器的应用仿真基于LMS和RLS的自适应滤波器应用仿真1. 自适应滤波原理自适应滤波器是指利用前一时刻的结果，自动调节当前时刻的滤波器参数，以适应信号和噪声未知或随机变化的特性，得到有效的输出，主要由参数可调的数字滤波器和自适应算法两部分组成，如图1所示自适应算法图1自适应滤波器原理图x(n)称为输入信号，y(n)称为输出信号，d(n)称为期望信号或者训练信号，e(n)为误差僖号，其中，e(n)=d(n)-y(n).自适应滤波器的系数(权值)根据误差信号e(n),通过一定的自适应算法不断的进行改变，以达到使输出信号y(n)最

2、接近期望信号图中参数可调的数字滤波器和自适应算法组成自适应滤波器。自适应滤波算法是滤波器系数权值更新的控制算法，根据输入信号与期望信号以及它们之间的误差信号，自适应滤波算法依据算法准则对滤波器的系数权值进行更新，使其能够使滤波器的输出趋向于期望信号。原理记数字滤波器脉冲响应为：h(k) = h0(k) hi (k) *(!,(Z:)x(/:) - x(k)xT(k)h(k -1)其中e(k) = y(k) - xT (k)h(k -1)e(k)表示先验误差。只因为它是由前一个采样时刻的系数算出的，在实际中，很多时候由于h(k)计算的复杂度而不能应用于实时控制。用8 ,1代换c;(Z),其中:5

3、为自适应梯度，I为辨识矩阵(n, n)这时h(k) = h(k 1) + 3x(k)e(k)这时就是一个最小均方准则问题。2. LMS自适应滤波器举例自回归过程的自适应预估器自回归过程是用来描述伴随一些可能性规律出现的统计现象的瞬时估计的随机过程。一阶自回归模型的公式如下：y(Z) = qy(Z l) + (Z)%是模型的唯一参数，b(k)是零均值白噪声。用一个自适应滤波器生成一个可以对参数出进行一步预测的一阶自适应预估器。LMS算法可由如下方程表示：yd)Ae(k) = y(k)CL(左 + 1)= i)e(Z)Cl d取N个点估计参数,为获取平均值重复M次。而且分别对6二0.01, 8=0

4、. 05, 6二0. 1进行计算。参数为固定在-0. 6。程序清单如下：N=500;M=20;n=l:al=-O.8;h=zeros (M, n+1, 3);c二zeros (M, n, 3);for d=l:3if d=l delta=0. 01;else delta=0. 05*(dT);end;for k=l:Mb=0. 2*randn (1, N);y =1;for i=2:Ny(i)=-al*y(i-l)+b(i);endfor i=n+l:Ne(k, i, d)=y(i)-h(k, i, d)*y(i-l);h(k, i + 1, d) =h(k, i, d) +delta*y (

5、i-l)*e(k, i, d);endendendfor d=l:3for i=l:Nem(i, d)=0;hm(i, d)=0;for j=l:Mem(i, d) =em(i, d)+e(j, -2;hm(i, d)=hm(i, d)+h(j, i, d);endendendfigure(l)semilogy(1:150, em(l:150, 1), hold onsemilogy(1:150, cm(1:150, 2), r), hold onsemilogy (1:150, em(l: 150, 3), g), hold offaxis(0 150 0. 01 1), gridtitle

6、 C Mean square error )xlabel ( Samples*)gtext ( leftarrowd=0. Of );gtext C leftarrowd=0.05,);gtext(5leftarrowd=0. T);figure (2), plot (1 :N, hm(l :N, 1), hold onplot (1 :N, hm(l :N, 2), * r1), hold onplot (1 :N, hm(l :N, 3), g), hold off, gridtitle( Filter coeffcient evalution，)xlabelC Samples),gtex

7、t ( d=0. Of ), gtext C d=O. 05), gtext ( d=0. f )得到的如下结果图：图2.平均方差误差图3.滤波器系数曲线系数以时间常数的指数曲线收敛，6越大，时间常数越小3. RLS自适应滤波器的应用仿真从噪声中提取信号输入信号为：x(k) = & cos(w0r + 裔)+ b(k)其中b(k)是附加的白噪声。xr(k) = AQ cos(“/ + 4)应用于RLS自适应滤波器的算法可描述如下：e(k) = y(k)- h7 (k V)x(k)h(k) = h(k 1) + g(k)e(k)c： *) = c(k -l)-g(k)xT (幻 c二(k -1)

8、y(k) = hT(k)x(k)所研究的滤波器阶数为200,程序清单如下：N=1000;n=200;k=12;Ts=le-1b=0.8*randn(1, N);for i=l:Nxr(1A i)=sin(k*2*pi*i/N);x (1, i) =xr (1, i) +b (i);自适应增益行向量，大小(1, n);先验误差自适应滤波器系数行向量，大小(l,n)输入信号x(k)的自相关转制矩阵，大小(1, n)自适应滤波输出采样周期等于ims.endCxx=10000*eye(n);g=zeros (N,n);h=zeros (N,n);e=zeros(1,N);y=zeros(1,N);tr

9、=zeros(1,N);for i=n+l:Ng(i, :) = (Cxx*x (i-n+1:i) 1./(1+x(i-n+1:i)*Cxx*x (i-n+1:i) f) 1 ;e (1z i)=xr (i)-h (i-lz :)*x (i-n+1:i) 1;h(i, :)=h(i-lz :)+e(1,i)*g (i,:);Cxx=Cxx-g(i, :) 1*x (i-n+1:i)*Cxx;y(lz i)=h(iz:)*x (i-n+1:i) *;tr (1,i)=trace (Cxx);endfigure(1)plot (0:N-nz x (lz n:N),gridtitle(1x(k) i

10、nput singnal in V1)xlabel(1 Samples 1)figure(2)plot(0:N-n,xr(1,n:N),r),gridaxis ( 0 800 -1.2 1.2)title(!xr(k) reference singnal in V1)xlabel(1 Samples 1)figure(3)plot(0:N-nz e(1,n:N)z hold onplot(0:N-n,y(1,n:N)Jr，)，hold ongridtitle (1e(k) error and y(k) output in V1)xlabel(1 Samples 1)gtext(* e(k) )

11、,gtext(y(k)figure(4)plot(0:N-n,h(n:NZ1),hold onplot(0:N-n,h(n:N,2)z r1),hold offgridtitle (1 a (n-1) and a (n-2) coeffcients evolution 1)xlabel(1 Samples 1)figure(5)numl=fliplr(h(N,:);sysl=tf(numl,1, Ts);bode(sysl),hold offtitle(1 Synthesized filter 1)xlabel(1 Frequency in rad/s1)ylabel ( 1 Phase in

12、 degree;Module in dB1)figure (6)semilogy(0:N-n,tr(n:N)rgridtitle(1Cxx matrix trace 1)xlabel(r Samples 1)实验结果图如下：图4,输入信号x(k)图5参考信号xr(k)图6误差e(k)和输出信号y(k)图7.滤波器系数a(n-l)和a(n-l)变化曲线系数的变化曲线在200步时有一个超调，这是由于h(k)向量为零，所以200步以后仅代表x值。获得的滤波器的传递函数也类似于LMS滤波器的传递函数，相应的预测也类似。它的中心频率调整为正弦信号频率，即75rad/s,如下图所示Eile Edi t Jix ew Inse