%利用神經網絡進行函數逼近
clear all
x=0:0.1*pi:4*pi;
y=sin(x);
%設定迭代次數
net.trainparam.epochs=10000;
%網絡初始化
net=newff([0,4*pi],[8,8,8,8,1],{'tansig','logsig','logsig','tansig','tansig'});
%訓練網絡
[net,tr,y1,e]=train(net,x,y);

X=0:0.01*pi:4*pi;
%網絡泛化
y2=sim(net,X);

subplot(2,1,2);
plot(X,y2);
title('網絡產生')
grid on
subplot(2,1,1);
plot(x,y,'o');
title('原始數據')
grid on