%取训练样本组和测试样本组
xx1=[xx(1,:);xx(2,:);xx(4,:);xx(5,:);xx(6,:);xx(7,:);xx(8,:);xx(10,:);xx(11,:);xx(12,:);xx(13,:);xx(15,:);xx(16,:);xx(17,:);xx(18,:);xx(19,:);xx(20,:);xx(22,:);xx(23,:);xx(24,:);xx(25,:);xx(26,:);xx(27,:);xx(28,:);xx(30,:);xx(31,:);xx(32,:);xx(34,:);xx(35,:);xx(36,:);];
yy1=[yy(1,:);yy(2,:);yy(4,:);yy(5,:);yy(6,:);yy(7,:);yy(8,:);yy(10,:);yy(11,:);yy(12,:);yy(13,:);yy(15,:);yy(16,:);yy(17,:);yy(18,:);yy(19,:);yy(20,:);yy(22,:);yy(23,:);yy(24,:);yy(25,:);yy(26,:);yy(27,:);yy(28,:);yy(30,:);yy(31,:);yy(32,:);yy(34,:);yy(35,:);yy(36,:);];
xx2=[xx(3,:);xx(9,:);xx(14,:);xx(21,:); xx(29,:);xx(33,:)];
yy2=[yy(3,:);yy(9,:);yy(14,:);yy(21,:); yy(29,:);yy(33,:)];
%TRAINBPX训练程序式
p=xx1'; t=yy1';
[R,Q]=size(p);[S2,Q]=size(t);S1=8;[w1,b1]=rands(S1,R);[w2,b2]=rands(S2,S1);
disp_fqre=20;max_epoch=6000;err_goal=0.001;lr=0.55;im=1.05;dm=0.7;mc=0.85;er=1.04;
TP=[disp_fqre max_epoch err_goal lr im dm mc er];
[w1,b1,w2,b2, te,tr]=trainbpx(w1,b1,'logsig',w2,b2,'logsig', p,t,TP);
%仿真数据处理
pp=xx2';
tt=[0.97 0.84 0.53 0.40 1.42 1.88;11.46 14.65 5.15 3.09 11.18 13.86;8.48 9.04 7.16 8.34 14.11 17.23;1.51 1.79 1.04 0.8 2.03 2.61];
aa=simuff(pp,w1,b1,'logsig',w2,b2,'logsig');m=aa(1,:);n=aa(2,:);l=aa(3,:);g=aa(4,:);y5=tt(1,:);y6=tt(2,:); y7=tt(3,:); y8=tt(4,:);
kk=[];
y1(i)=0.4+(m(i)-0.1)*(1.88-0.4)/0.8;
y2(i)=3.09+(n(i)-0.1)*(14.65-3.09)/0.8;
for i=1:6
y3(i)=6.26+(l(i)-0.1)*(17.23-6.26)/0.8;
y4(i)=0.8+(g(i)-0.1)*(2.61-0.8)/0.8;
kk=[kk;y1;y2;y3;y4];
end
for j=1:6
ll=[];
y9(j)=(abs((y5(j)-kk(1,j))/y5(j)))*100;
y10(j)=(abs((y6(j)-kk(2,j))/y6(j)))*100;
y11(j)=(abs((y7(j)-kk(3,j))/y7(j)))*100;
y12(j)=(abs((y8(j)-kk(4,j))/y8(j)))*100;
ll=[ll,y9;y10;y11;y12];
end