资源描述：
《基于RLS的自适应滤波器应用仿真.docx》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、基于RLS的自适应滤波器应用仿真南昌大学信息工程学院《现代信号处理》课程作业题目：基于RLS的自适应滤波器应用仿真指导老师：万国金作者：毕圣昭日期：2011-12-058基于RLS的自适应滤波器应用仿真基于RLS的自适应滤波器应用仿真1.自适应滤波原理自适应滤波器是指利用前一时刻的结果，自动调节当前时刻的滤波器参数，以适应信号和噪声未知或随机变化的特性，得到有效的输出，主要由参数可调的数字滤波器和自适应算法两部分组成，如图1所示：图1自适应滤波器原理图x(n)称为输入信号，y(n)称为输出信号，d(n)称
2、为期望信号或者训练信号，e(n)为误差信号，其中，e(n)=d(n)-y(n)。自适应滤波器的系数(权值)根据误差信号e(n)，通过一定的自适应算法不断的进行改变，以达到使输出信号y(n)最接近期望信号。图中参数可调的数字滤波器和自适应算法组成自适应滤波器。自适应滤波算法是滤波器系数权值更新的控制算法，根据输入信号与期望信号以及它们之间的误差信号，自适应滤波算法依据算法准则对滤波器的系数权值进行更新，使其能够使滤波器的输出趋向于期望信号。记数字滤波器脉冲响应为：h(k)=[h0(k)h1(k)…hn-1(
3、k)]T输入采样信号为：x(k)=[x(k)x(k-1)…x(k-n-1)]误差信号为：优化过程就是最小化性能指标J(k),它是误差的平方和：求使J(k)最小的系数向量h(k),即使J(k)对h(k)的导数为零，也就是把J(k)的表达式代入，得：8基于RLS的自适应滤波器应用仿真和由此得出滤波器系数的最优向量：这个表达式由输入信号自相关矩阵和输入信号与参考信号的相关矩阵组成，如下所示，维数都为（n,n）:系数最优向量也可以写成如下形式：自相关和互相关矩阵的递归表达式如下：把的递归表达式代入系数向量表达式，
4、得：即考虑到可以记8基于RLS的自适应滤波器应用仿真用前面得到的表达式求出，并代入上式：或则滤波器系数的递归关系式可以记作其中e(k)表示先验误差。只因为它是由前一个采样时刻的系数算出的，在实际中，很多时候由于h(k)计算的复杂度而不能应用于实时控制。用δ,I代换，其中δ为自适应梯度，I为辨识矩阵（n，n），这时这就是一个最小均方准则问题。2.RLS自适应滤波器的应用仿真：从噪声中提取信号输入信号为：其中b(k)是附加的白噪声。应用于RLS自适应滤波器的算法可描述如下：自适应增益行向量，大小（1，n）;先
5、验误差自适应滤波器系数行向量，大小（1,n）输入信号x(k)的自相关转制矩阵，大小（1,n）自适应滤波输出所研究的滤波器阶数为200，采样周期等于1ms.8基于RLS的自适应滤波器应用仿真程序清单如下：N=1000;n=200;k=12;Ts=1e-1;b=0.8*randn(1,N);fori=1:Nxr(1,i)=sin(k*2*pi*i/N);x(1,i)=xr(1,i)+b(i);endCxx=10000*eye(n);g=zeros(N,n);h=zeros(N,n);e=zeros(1,N);
6、y=zeros(1,N);tr=zeros(1,N);fori=n+1:Ng(i,:)=(Cxx*x(i-n+1:i)'./(1+x(i-n+1:i)*Cxx*x(i-n+1:i)'))';e(1,i)=xr(i)-h(i-1,:)*x(i-n+1:i)';h(i,:)=h(i-1,:)+e(1,i)*g(i,:);Cxx=Cxx-g(i,:)'*x(i-n+1:i)*Cxx;y(1,i)=h(i,:)*x(i-n+1:i)';tr(1,i)=trace(Cxx);endfigure(1)plot(0:N
7、-n,x(1,n:N)),gridtitle('x(k)inputsingnalinV')xlabel('Samples')figure(2)plot(0:N-n,xr(1,n:N),'r'),gridaxis([0800-1.21.2])title('xr(k)referencesingnalinV')xlabel('Samples')figure(3)plot(0:N-n,e(1,n:N)),holdonplot(0:N-n,y(1,n:N),'r'),holdonaxis([0,800,-1.5,1
8、.5])gridtitle('e(k)errorandy(k)outputinV')xlabel('Samples')gtext('e(k)'),gtext('y(k)')8基于RLS的自适应滤波器应用仿真figure(4)plot(0:N-n,h(n:N,1)),holdonplot(0:N-n,h(n:N,2),'r'),holdonaxis([0,800,-0.05,0.05])gridtitle('a(n-1)anda
