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