BP神经网络函数模拟

clear;

%%定义样本数据
x1=20*rand（1100）-10;
x2=20*rand（1100）-10;
x3=20*rand（1100）-10;
x4=20*rand（1100）-10;

%%划分训练数据与测试数据
p = [x1;x2];
t = [x1+x2-4];
intest = [x3;x4];
outest = [x3+x4-4];

%%创建并训练网络
net = newff（pt[16]）;
net = train（netpt）;
yy = sim（netintest）;

%%输出图像
l1 = linspace（-101050）;
l2 = linspace（-101050）;
[X1X2] = meshgrid（l1l2）;
P = [X1（:）‘;X2（:）‘];
Z = sim（netP）;
Y = reshape（Z5050）;
figure（1）
surf（X1X2Y）;hold on;
axis（[-1010-1010min（Y（:））max（Y（:））]）;
youtest = sim（netintest）;
plot3（intest（1:）intest（2:）youtest‘b*‘）
xlabel（‘X1‘）;
ylabel（‘X2‘）;
zlabel（‘Y‘）;
legend（‘拟合结果‘‘样本点‘）;
title（‘Y=X1+X2－4‘）
hold off 

%%计算、显示相对误差等
e = （outest-yy）./outest;
figure（2）;
plot（1:10