%要逼近的函數樣本點
P = -1:.1:1;
T=sin(pi*P);
plot(P,T,'+');
title('待逼近的函數樣本點');
xlabel('輸入值');
ylabel('目標值');
%徑向基傳遞函數
p = -3:.1:3;
a = radbas(p);
figure;
plot(p,a)
title('徑向基傳遞函數');
xlabel('輸入');
ylabel('輸出'); 
%徑向基傳遞函數及其加權和
a2 = radbas(p-1.5);
a3 = radbas(p+2);
a4 = a + a2*1 + a3*0.5;
figure;
plot(p,a,'b--',p,a2,'b--',p,a3,'b--',p,a4,'m-')
title('徑向基函數的加權和');
xlabel('輸入');
ylabel('輸出');
%建立網絡
eg = 0.02; %設置平方和誤差參數
sc = 1;    %設置展開常量
net=newrb(P,T,eg,sc);
%仿真網絡并繪出結果
figure;
X=-1:.01:1;
Y=sim(net,X);
hold on;
plot(X,Y);
plot(P,T,'+');
title('仿真結果');
xlabel('輸入');
ylabel('網絡輸出及目標輸出');
hold off;