智能控制第二版刘金琨思考题与习题7-1

网络训练程序

clear all;

close all;

xite=0.50;

alfa=0.05;

w2=rands(6,1);

w2_1=w2;w2_2=w2;

w1=rands(2,6);

w1_1=w1;w1_2=w1;

dw1=0*w1;

I=[0,0,0,0,0,0]';

Iout=[0,0,0,0,0,0]';

FI=[0,0,0,0,0,0]';

k=0;

E=1.0;

NS=3;

while E>=1e-020

k=k+1;

times(k)=k;

for s=1:1:NS

xs=[1,0;

0,0;

0,1];

ys=[1,0,-1]';

x=xs(s,:);

for j=1:1:6

I(j)=x*w1(:,j);

Iout(j)=1/(1+exp(-I(j))); end

y1=w2'*Iout;

el=0;

y=ys(s,:);

el=el+0.5*(y(1)-y1(1))^2;

es(s)=el;

E=0;

if s==NS

for s=1:1:NS

E=E+es(s);

end

end

ey=y-y1;

w2=w2_1+xite*Iout*ey+alfa*(w2_1-w2_2);

for j=1:1:6

S=1/(1+exp(-I(j)));

FI(j)=S*(1-S);

end

for i=1:1:2

for j=1:1:6

dw1(i,j)=xite*FI(j)*x(i)*ey(1)*w2(j,1) end

end

w1=w1_1+dw1+alfa*(w1_1-w1_2);

w1_2=w1_1;w1_1=w1;

w2_2=w2_1;w2_1=w2;

end %End of for

Ek(k)=E;

end %End of while

figure(1);

plot(times,Ek,'r');

xlabel('k');ylabel('E');

save wfile w1 w2;

网络测试程序

clear all;

load wfile w1 w2;

%N Samples

x=[1,0.1;

0.5,0.5;

0.1,1];

for i=1:1:3

for j=1:1:6

I(i,j)=x(i,:)*w1(:,j);

Iout(i,j)=1/(1+exp(-I(i,j)));

end

end

y=w2'*Iout';

y=y'