N=500;M=20;n=1;a1=-0.8;h=zeros(M,n+1,3);e=zeros(M,n,3);ford=1:3ifd==1delta=0.01;elsedelta=0.05*(d-1);end;fork=1:Mb=0.2*randn(1,N);y(1)=1;fori=2:Ny(i)=-a1*y(i-1)+b(i);endfori=n+1:Ne(k,i,d)=y(i)-h(k,i,d)*y(i-1);h(k,i+1,d)=h(k,i,d)+delta*y(i-1)*e(k,i,d);endendendford=1:3fori=1:Nem(i,d)=0;hm(i,d)=0;f