% bp8;例4.8
%
P=[-6 -6.1 -4.1 -4 4 4.1 6 6.1];
T=[0.0 0.0 0.97 0.99 0.01 0.03 1 1];
net=newcf(minmax(P),[1],{'tansig'},'traingd');      
                                   %創建前向BP網絡
net.iw{1,1};net.b{1};
net.trainParam.epochs=300;            %初始化訓練次數
net.trainParam.lr=0.05;
[net tr]=train(net,P,T);                  %訓練網絡
w1=net.iw{1,1},b1=net.b{1},
Wrange=-1:0.1:1;Brange=-2:0.2:2;        %W值的行向量、B值的行向量
ES=errsurf(P,T,Wrange,Brange,'logsig');    %求單神經元的誤差平面（只用在單神經元中）
mesh(ES,[60,30]);                     %作三維網狀面，視角【60，30】 
title('Error Surface Graph')
xlabel('W');
ylabel('B');
zlabel('Error')
figure(2)
[C,h] =contour(Wrange,Brange,ES,6);     %作等高線圖，ES為高
                                   %返回等高線矩陣C，列向量h是線或對象的句柄，
                                   %一條線一個句柄，這些被用作CLABEL的輸入， 
                                   %每個對象包含每個等高線的高度 
clabel(C,h)                          %標上高度值
colormap cool                        %背景的顏色cool
hold on
plot(w1,b1,'r+')
xlabel('W');
ylabel('B');
hold off
figure(3)
plot(tr.perf)                           %畫誤差曲線