LMS算法MATLAB代码

clc;
clear all
N=5000;
u=0.0002;
umax=0.01;
umin=0.000002;
p=0.0965;
%p=1;
v=wgn（1Np）;
x（1）=v（1）;
x（2）=0.2*x（1）+v（2）;
for n=3:N;
    x（n）=0.2*x（n-1）-0.9*x（n-2）+v（n）; 
end;
w1（3）=0;
w2（3）=0;
ww1（2）=0;
ww1（3）=0;
ww2（3）=0;
for n=3:N-1;
    e（n）=x（n）-x（n-1）*ww1（n）-x（n-2）*ww2（n）;
    ww1（n+1）=ww1（n）+u*x（n-1）*e（n）;
    ww2（n+1）=ww2（n）+u*x（n-2）*e（n）;
end;
uu=zeros（1N）;
uu（1）=umax;
uu（2）=umax;
uu（3）=umax;
p=zeros（1N）;
A=zeros（1N）;
p（1）=umax;
p（2）=umax;
alpha=0.9;
gama=0.000001;
B=0.9;
e（2）=x（2）-x（1）*ww1（2）;
e（1）=x（1）;
W1=zeros（1N）;
W2=zeros（1N）;
t=1;
L=10;
for i=1:t
for n=3:N-1;
    e（n）=x（n）-x（n-1）*w1（n）-x（n-2）*w2（n）;
    qiuhe=0;
    for