clear
clc
alpha=0.05;
eta=0.2;
ny=2;d=2;nu=3;
n=ny+nu-d+1;
m=6;
w1k1=randn(m,n);%6*4
w1k2=w1k1;
w2k1=randn(1,m);%1*6
w2k2=w2k1;
L=1000;
uk=zeros(nu,1);%2*1
yk=zeros(ny,1);%2*1
fork=1:L
time(k)=k;
u(k)=0.8*sin(0.01*pi*k);
y(k)=uk(2)^3+uk(3)^3+(0.8+yk(1)^3)/(1+yk(1)^2+yk(2)^4);
X=[yk;uk(d:nu)];
o=X;%4*1
net2=w1k1*X;% 6*1
o2=1./(1+exp(-net2));%6*1
ym(k)=w2k1*o2;
e(k)=y(k)-ym(k);
dw2=eta*e(k)*o2';%1*6 e(k)ÊÇÒ»¸öÊý£¬kָϱꡣ
w2=w2k1+dw2+alpha*(w2k1-w2k2);%1*6
dnet2=exp(-net2)./(1+exp(-net2)).^2;%6*1
dw1k1=eta*e(k)*w2k1'.*dnet2*o';
w1=w1k1+dw1k1+alpha*(w1k1-w1k2);
w1k2=w1k1;w1k1=w1;
w2k2=w2k1;w2k1=w2;
for i=nu:-1:2
uk(i)=uk(i-1);
end
uk(1)=u(k);
for i=ny:-1:2
yk(i)=yk(i-1);
end
yk(1)=y(k);
end
subplot(211)
plot(time,y,'r:',time,ym,'k');
xlabel('k'); ylabel('y(k)¡¢y_m(k)');
legend('y(k)','y_m(k)'); %axis([0 L -.4 1.6]);
subplot(212)
plot(time,y-ym,'k');
xlabel('k'); ylabel('e(k)'); axis([0 L -1 1]);
