function [P,T] = getTestData();
sizeLen = 0.2;
x1 = 0:sizeLen:1;
[line,row] = size(x1);
P = zeros(10000,5);
T = zeros(10000,3);
j = 1;
for x1 = 0:sizeLen:1
for x2 =
0:sizeLen:1
for x3 = 0:sizeLen:1
for x4 = 0:sizeLen:1
for x5 = 0:sizeLen:1
y1=x1*x1+x2*x3+x4*x5;
y2=x1+x2+x3+x4+x5;
y3=sin(x1+x2+x3)+cos(x4+x5);
P(j,:) = [x1,x2,x3,x4,x5];
T(j,:) = [y1,y2,y3];
j = j+1;
end
end
end
end
end
clear;clc;
[P,T] = getTestData();
[line,row] = size(P);
j1=1;
j2=1;
for i=1:line
if
0==mod(i,2)
p1(j1,:) = P(i,:);
t1(j1,:) = T(i,:);
j1 = j1+1;
else
p2(j2,:) = P(i,:);
t2(j2,:) = T(i,:);
j2 = j2+1;
end
end
%
p = p1';
t = t1';
% p = [-1 -1 2 2;0 5 0 5];
% t = [-1 -1 1 1];
net=newff(minmax(P'),[5,3],{'logsig','purelin'},'traingdm');
net.trainParam.show = 50;
net.trainParam.lr = 0.05;
net.trainParam.mc = 0.9;
net.trainParam.epochs = 1000;
net.trainParam.goal = 1e-5;
[net,tr]=train(net, p1', t1');
net.trainParam.show = 50;
net.trainParam.lr = 0.05;
net.trainParam.mc = 0.9;
net.trainParam.epochs = 1000;
net.trainParam.goal = 1e-5;
[net,tr]=train(net, p2', t2');
net.trainParam.show = 50;
net.trainParam.lr = 0.05;
net.trainParam.mc = 0.9;
net.trainParam.epochs = 1000;
net.trainParam.goal = 1e-5;
[net,tr]=train(net, p1', t1');
net.trainParam.show = 50;
net.trainParam.lr = 0.05;
net.trainParam.mc = 0.9;
net.trainParam.epochs = 1000;
net.trainParam.goal = 1e-5;
[net,tr]=train(net, p2', t2');