flch3eg8.m

刚刚接触神经网络的朋友们看了这本书会有一定得帮助的！！

clear
w0=[0.05 0.2 0.1], w1=[0.2 -0.02 0.01]
x=[1 1 1]'; M=80; N=50;%x=[0.1 0.2 0.1]'
a1=0.5; a2=0.5; uk0=1; q=0.024;
y1=0.1; y0=0.2;
d=randn(M,1);
for  k=1:M;
        d1=d(k);
    yp=1-exp(-k/2);
    for  i=1:N
        I=w1*x-0.1; 
        uk=1/(1+exp(-I));
 y2=0.33*y1+0.132*y0+0.5*uk+0.038*uk0;
        y8=0.33*y1+0.132*y0+0.51*uk+0.031*uk0+q*d1;
  y0=y1; y1=y2; uk0=uk;
  e1=yp-y8;
  w2=w1+0.682*e1*x'+0.113*(w1-w0);
  w0=w1; w1=w2;
 % e3=yp-y2
  %e1=[(yp-y8)^.2]/2;
  if e1<=0.001,break
  else
  end
  
    end
  y3(k)=y2;
  y4(k)=y8;
  yp1(k)=yp;
  w(k)=q*d1;
  e(k)=e1;
 % e2(k)=e3;
end
    w2
    %grapher
    k=1:M;
    subplot(4,1,1)
    plot(k,yp1,k,y3,k,y3,'r*')
    ylabel('yp1 , y3')
    legend('yp1=input','y3=output')
    subplot(4,1,2)
    plot(k,yp1,k,y4,k,y4,'r+')
    ylabel('yp1 , y4')
    legend('yp1=input','y4=output')
    subplot(4,1,3)
    plot(k,w)
    ylabel('w') 
    subplot(4,1,4)
    plot(k,e)
    xlabel('k'),ylabel('e')