nnliu1seperate0.m

一个bp实现回归的nn程序

[e-mail]
%%jiangnan university
%%2007.7.15
%%rbf
clc;
clear;
 
 
f1=[a1 b1 c1   ]';
f2=[d1 g1]';
for i=1:3
    f1(i,:)=(f1(i,:)-min(f1(i,:)))/(max(f1(i,:))-min(f1(i,:)));
end
f2(1,:)=(f2(1,:)-min(f2(1,:)))/(max(f2(1,:))-min(f2(1,:)));
f1_train=[f1(:,2) f1(:,4) f1(:,5) f1(:,7) f1(:,8) f1(:,10) f1(:,11) f1(:,13) f1(:,14) f1(:,16) f1(:,17) f1(:,19) f1(:,20) f1(:,22) f1(:,23)];
f2_train=[f2(:,2) f2(:,4) f2(:,5) f2(:,7) f2(:,8) f2(:,10) f2(:,11) f2(:,13) f2(:,14) f2(:,16) f2(:,17) f2(:,19) f2(:,20) f2(:,22) f2(:,23)];
f1_test=[f1(:,1) f1(:,3) f1(:,6) f1(:,9) f1(:,12) f1(:,15) f1(:,18) f1(:,21) f1(:,24) f1(:,27)];
f2_test=[f2(:,1) f2(:,3) f2(:,6) f2(:,9) f2(:,12) f2(:,15) f2(:,18) f2(:,21) f2(:,24) f2(:,27)];
no=[3 9 12 25 47];
for i=1:5
%net=newff(minmax(f1),[no(i),2],{'tansig','logsig'});%build a netPURELIN
net=newff(minmax(f1),[no(i),2],{'PURELIN','logsig'},'TRAINBFG','LEARNGD');
%net=newff(minmax(f1),[no(i),2],{'tansig','PURELIN'});%aaa=minmax(f1);
net.trainparam.epochs=1500;
 net.trainParam.goal = 0.00000;
net=init(net);
 
net=train(net,f1_train,f2_train);
temp=sim(net,f1_test);
y(2*i-1,:)=temp(1,:);
y(2*i,:)=temp(2,:);
end
y1=[y(1,:);y(2,:)];
y2=[y(3,:);y(4,:)];
y3=[y(5,:);y(6,:)];
y4=[y(7,:);y(8,:)];
y5=[y(9,:);y(10,:)];
for i=1:10
    error1(i)=norm(y1(1,i)-f2_test(1,i));%just 1line
    error2(i)=norm(y2(1,i)-f2_test(1,i));
    error3(i)=norm(y3(1,i)-f2_test(1,i));
    error4(i)=norm(y4(1,i)-f2_test(1,i));%good effect
    error5(i)=norm(y5(1,i)-f2_test(1,i));
end
figure(1)
plot(1:10,error1,'rx');
%figure(2)
hold on;
plot(1:10,error2,'--');
hold on;
%figure(3)
plot(1:10,error3,'+');
hold on;
%figure(4)
plot(1:10,error4,'y+');
hold on;
%figure(5)
plot(1:10,error5,'g:');
hold on;

for i=1:5
figure(i+5)
%for i=1:5
plot(1:10,y(2*i-1,:),'g:',1:10,f2_test(1,:),'k-.')
hold on
end