example2.m

ANN基本知识和应用课件

clear all

x=[-15 -10 -5 0 5 10 15];
y=0.05*x.^3-0.2*x.^2-3*x+20;

x1=[-15:1:15];
y1=0.05*x1.^3-0.2*x1.^2-3*x1+20;

x2=[-20:1:20];
y2=0.05*x2.^3-0.2*x2.^2-3*x2+20;

%--------------------------------------------------

net = newff([-15 15],[3 3 1],{'logsig' 'logsig' 'purelin'});

net.trainParam.epochs = 10000;
net.trainParam.goal = 1e-3;

net = init(net);
net = train(net,x1,y1);

Test=[-20:0.1:20];
Output=sim(net,Test);

%--------------------------------------------------

figure(1);
plot(x,y,'b*',x2,y2,'b')
title('Training Vectors');
xlabel('Input Vector x');
ylabel('Target Vector y for trainning');
legend({'Input Point','Desired Output y=0.05*x^3-0.2*x^2-3*x+20'})

%--------------------------------------------------

figure(2);
plot(x,y,'b*',x2,y2,'b',Test,Output,'r')
title('Testing MLP Network');
xlabel('Input Vector');
ylabel('Output');
legend({'Input point','Desired output','Test output'})