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