lms.m

本代码是我自己写的

k=[1:5]
N=5;
xk=cos(2*pi*k/N);
dk=2*sin(2*pi*k/N);
xk1=cos(2*pi*(k-1)/N);
R=[mean(xk*xk'),mean(xk*xk1');mean(xk1*xk'),mean(xk1*xk1')];
P=[mean(dk*xk');mean(dk*xk1')];
Xk1=[xk(1);xk1(1)];
Xk2=[xk(2);xk1(2)];
Xk3=[xk(3);xk1(3)];
Xk4=[xk(4);xk1(4)];
Xk5=[xk(5);xk1(5)];
W=[0;0];
ek1=dk(1)-Xk1'*W;
ek2=dk(2)-Xk2'*W;
ek3=dk(3)-Xk3'*W;
ek4=dk(4)-Xk4'*W;
ek5=dk(5)-Xk5'*W;
i=0;
for x=1
    r=rem(i,5);
    if r=0
    W=W+2*0.05*ek1*Xk1;
    elseif r=1
       W=W+2*0.05*ek2*Xk2; 
    elseif r=2
        W=W+2*0.05*ek3*Xk3;
    elseif r=3
        W=W+2*0.05*ek4*Xk4;
    else r=4
        W=W+2*0.05*ek5*Xk5;      
    end
    i=i+1;
    x=le(i,10000);
end
plot(i,W)