lms222.m

应用LMS算法

clc;
clear;
u=0.001;        %  0<u<<1
a1=1;
a2=2;
N=5000;

var=1;
v=0+sqrt(var)*randn(1,N);   %   零均值、方差为var的白噪声

var1=0.01;
v1=0+sqrt(var1)*randn(1,N);   %   零均值、方差为var1的白噪声

w1(2)=0;        %   设定各参数初值
w2(2)=0;
x=v;

for n=2:N
    y(n)=w1(n)*x(n)+w2(n)*x(n-1);
    d(n)=a1*x(n)+a2*x(n-1)+v1(n);
    e(n)=d(n)-y(n);
    w1(n+1)=w1(n)+2*u*x(n)*e(n);
    w2(n+1)=w2(n)+2*u*x(n-1)*e(n);
end

figure(1);
    plot(w1,'g');
    hold on
    plot(w2,'r');
    legend('w1','w2')
    xlabel('n');
    ylabel('w(n)');
    title('w(n)变化的曲线    u=0.001, σ^2=1, a1=1, a2=2');
    hold off

figure(2);
    plot(e.*e);
    xlabel('n');
    ylabel('e^2(n)');
    title('e^2(n)变化的曲线');