%ch8.3
clc;
clear;
close all;
%%%%%%%%%%%%%%%%%%%%%%%%%訓練樣本%%%%%%%%%%%%%%%%%%%%%%%%%
P=-1:0.1:1;
T=[-0.9602 -0.5770 -0.0729 0.3771 0.6405 0.6600 0.4609 ...
    0.1336 -0.2013 -0.4344 -0.5000 -0.3930 -0.1647 0.0988 ...
    0.3072 0.3960 0.3449 0.1816 -0.0312 -0.2189 -0.3201];
plot(P,T,'+');           %繪制輸入/輸出矢量點
xlabel('輸入矢量P');
ylabel('目標矢量T');
%%%%%%%%%%%%%%%%%%%%%設計一個徑向基函數網絡%%%%%%%%%%%%%%%
p=-3:.1:3;
a=radbas(p);
plot(p,a)              %繪制輸入輸出二維圖形
xlabel('輸入p');
ylabel('輸出a');
%%%%%%%%%%%%%%%%%%%%%%%徑向基傳遞函數權值之和%%%%%%%%%%%%%
a2=radbas(p-1.5);
a3=radbas(p+2);
a4=a+a2*1+a3*0.5;
plot(p,a,'b-',p,a2,'b-',p,a3,'b-',p,a4,'m--')
xlabel('輸入p');
ylabel('輸出a');
eg=0.02;
sc=1;
net=newrb(P,T,eg,sc);               %創建一個神經網絡
Newrb,neurons=0,sse=3.69051  
%%%%%%%%%%%%%%%%%%%%%%%%網絡測試%%%%%%%%%%%%%%%%%%%%%%%%%%
plot(P,T,'+');
xlabel('輸入');
X=-1:.01:1;
Y=sim(net,X);                      %程序仿真
hold on 
plot(X,Y);
hold off
legend({'目標','輸出'});           %在圖形中插入注釋