lsl lms.m

此程序为自适应信号处理中

clear 
close all

clear;
close all
N=1024;
M=3;
a1=1.558;
a2=-0.81;
kf=zeros(M,N);
kb=zeros(M,N);
w1=zeros(1,N);
w2=zeros(1,N);
ef=zeros(M,N);
eb=zeros(M,N);
delta=zeros(M,N);
dd=1.00;
epf=dd*ones(M,N);
epb=dd*ones(M,N);
gma=ones(M,N);
u=0.002;
x= randn(1,N);
xu(1)=x(1);
xu(2)=a1*xu(1)+x(2);
for n=3:N
    xu(n)=x(n)+a1*xu(n-1)+a2*xu(n-2);
end

for n=2:N
   eb(1,n)=xu(n-1);
   ef(1,n)=xu(n-1);
   epb(1,n)=epf(1,n-1)+xu(n-1)^2;
   epf(1,n)=epf(1,n-1)+xu(n-1)^2;
   gma(1,n)=1;
   
   for m=2:3
       delta(m,n)=delta(m,n-1)+eb(m-1,n-1)*ef(m-1,n)/gma(m-1,n);
       ef(m,n)=ef(m-1,n)-delta(m,n)*eb(m-1,n-1)/epb(m-1,n-1);
       eb(m,n)=eb(m-1,n-1)-delta(m,n)*ef(m-1,n)/epf(m-1,n);
       epf(m,n)=epf(m-1,n)-delta(m,n)^2/epb(m-1,n-1);
       epb(m,n)=epb(m-1,n-1)-delta(m,n)^2/epf(m-1,n);
       gma(m,n-1)=gma(m-1,n-1)-eb(m-1,n-1)^2/epb(m-1,n-1);
       kf(m,n)=delta(m,n)/epf(m-1,n);
       kb(m,n)=delta(m,n)/epb(m-1,n-1);
   end
     w1(n)=kb(2,n)-kf(2,n)*kb(3,n);
     w2(n)=kb(3,n);
end

        

wt1=zeros(1,N);
wt2=zeros(1,N);
e=zeros(1,N);
for n=3:N-1
    y(n)=wt1(n)*xu(n-1)+wt2(n)*xu(n-2);
    e(n)=xu(n)-y(n);
    wt1(n+1)=wt1(n)+2*u*e(n)*xu(n-1);
    wt2(n+1)=wt2(n)+2*u*e(n)*xu(n-2);
end


plot(w1)
hold on
plot(w2)
plot(wt1,'r')
plot(wt2,'r')
hold off 
axis([0,1024,-2,2])
title('LSL和LMS权系数变化曲线比较')
xlabel('n');
ylabel('w(n)');