LMS滤波器matlab实现

1、矩阵在LMS自适应滤波器的应用1.自适应滤波原理自适应滤波器是指利用前一时刻的结果，自动调节当前时刻的滤波器参数， 以适应信号和噪声未知或随机变化的特性，得到有效的输出，主要由参数可调的 数字滤波器和自适应算法两部分组成，如图1所示x(n)称为输入信号，y(n)称为输出信号，d(n)称为期望信号或者训练信号, e(n)为误差信号，其中，e(n)=d(n)-y (n).自适应滤波器的系数(权值)w(i)根据 误差信号e(n)，通过一定的自适应算法不断的进行改变，以达到使输出信号y(n) 最接近期望信号1。图中参数可调的数字滤波器和自适应算法组成自适应滤波器。自适应滤波算 法是滤波器系数权值更新的

2、控制算法，根据输入信号与期望信号以及它们之间的 误差信号，自适应滤波算法依据算法准则对滤波器的系数权值进行更新，使其能 够使滤波器的输出趋向于期望信号。将上述表示成矩阵形式，输出人二片二旷勺，式中，w = D 口2 ,违,X = 7/, 口述、误差信号被作为加权系数的控制信号，rjL=rjL-均方误差（性能函数）E（D7）= E（口口一 U/j）2 =- 2LQ口 + rP,（1）其中口口口是期望信号d和输入信号X的互相关矩阵。上式表明，当输入信号时平稳随机信号时，均方误差函数E（述）是权系数W 的二次函数，它是一个中间上凹的超拋物形曲面，具有唯一最小值的函数用梯度向量，使性能函数到达它的最小

3、点，V =呂欝，包欝，,为求最佳权系数，令Vo=0,即由信号处理理论可知，当滤波器的单位脉冲响应取最佳值时，其误差信号和 输入信号是正交的。或者对（1）求导，得到2口 ，口一 2匚7口口令它为0,可得到最佳权矢量的表达式口 =口口（3）此时，均方误差取最小值，将（1）中均方误差为2,将最优权值系数带入，可得2 =+ 口口 口口口-2LfjU+ UP 口口 = 口口 + (7-口丫 式- 口）记 V = - , 则 得 到 4 = +口口（4）将口口口进行分解，得到口口口 = Ha口,其中，Q为正交矩阵或者特征矩阵， 八是由特征值组成的对角矩阵，用下式表示A = Diag（7 口，口口），并带入

4、（4）式有，4 = + EJP口 A tzP （5）最.徒？降法的递推公式为口庆1 = LJ +（一V ）,（6）-V是性能函数下降最快的方向。口为步长因子，用来控制迭代的速度，口 的值越大，则速度越快，但是精度下降，口的值太小，则算法运算量增大2。因此需要根据要求作出适当的选择。7丁 2口口口口一2口带入(6),+ (2口口口口+2口口龙两边都减去口，做变换，得到+1= 一 2口口口己口, 递推有=、- 2口口口己口口0(7)得到权系数的递推公式，口=七 ujj- 2口小z7)(8)2. LMS自适应滤波器举例本实验通过一个二阶自回归过程来研究实时数据集平均对LMS算法的影响，AR模型的差分

5、方程为：y(n) + 口口(口) + G口_!)=其中0=1.558； 4=-0. 81；口(口)是零均值方差为的白噪声；图3为AR模型及其二阶自适应线性预测模型，S3: AR模型将权初始值设为0,根据LMS算法的基本步骤可以写出该算法的matlab程序如下:clearclose allclcal=1.588;a2 二-0.81；u=0. 001；N=1024；G=100；e=zeros(l,N)；wl=zeros(l,N+l)；w2=zeros(l,N+l)；y=zeros(l,N)；ee=zeros(l,N);%每个点的误差平方ep=zeros (1,N) ；%每个点的误差平方累积eq=z

6、eros (1,N) ；%每个点的100次误差平方均值wll=zeros(1,N+l) ；%wl权值的累积w22=zeros (1, N+l); %w2权值的累积 for g=l：G v=randn(l,N)； xg=v ；x(2)=x *al+v (2)；for n=3：Nx (n) =al*x (n-l)+a2*x (n-2) +v (n);endfigure (1)plot (x)titleC输入信号x)for n=3：Ny (n) =wl (n) *x (n-1) +w2 (n) *x (n-2)；e (n) =x (n) -y (n)；wl (n+l) =wl (n) +2*u*e

7、(n) *x (n-1)；w2 (n+l) =w2 (n) +2*u*e (n) *x (n-2);ee(n) =e (n厂 2；endwll=wl+wll;w22=w2+w22；ep=ep+ee；endeq=ep/G；Wl=wll/G；W2=w22/G；figure(2)subplot (2,1,1)plot (wl)hold onsubplot(2,1,2)plot (Wl)hold onsubplot (2,1,1)plot (w2,frf)titleCwl与w2的收敛曲线,u二0. 001)hold onsubplot(2,1,2)plot (W2,r)；title C100次平均后w

8、l与w2的收敛曲线,u=0. 001) figure(3) subplot(2,1,1)plot (e)titled误差曲线(学习曲线)u=0. 001f) subplot(2,1,2)plot (eq)titleClOO次平均误差曲线(学习曲线)u=0. oor)下面对结果进行分析:图5为初始权值为0, u二0.001时的误差曲线，即所谓的学习曲线，收敛速度 较快的为100次集平均学习曲线，而起伏较大的为单次实现的学习曲线.图6平方误差曲线(学习曲线)u=0. 001图6为误差的收敛曲线，可以看出100次平均误差的收敛效果较快。3.结论信号通过系统的响应，在数学上可以等效为矩阵的运算，这种等效对于求解 输出信号是很方便的。应用矩阵的变换和向量的性质求解自适应滤波器的权系数，通过观测矩阵中 值的变化能比较直观的反映出滤波器在迭代过程中性能