%?RBF網(wǎng)絡(luò)逼近，多維輸入情況
%?作者：nuaazdh
%?時(shí)間：2012年4月11日?20:11:25
function?[test_out]=RBF_NN_MIC（train_itrain_otest_i）
%?train_i:訓(xùn)練樣本輸入，行數(shù)表示其維數(shù)，列數(shù)表示樣本個(gè)數(shù)
%?train_o:訓(xùn)練樣本輸出
%?test_i:測(cè)試樣本輸入
%?test_o:測(cè)試樣本輸出
%%?參數(shù)設(shè)定和初始化
size_tr=size（train_i2）;?%訓(xùn)練樣本個(gè)數(shù)
size_te=size（test_i2）;??%測(cè)試樣本個(gè)數(shù)
%train_o=cos（train_i）;%因變量y
in_num=size（train_i1）;%輸入層神經(jīng)單元個(gè)數(shù)
hid_num=7;?????????%隱含層神經(jīng)單元個(gè)數(shù)
w=rand（1hid_num）;%隱含層-輸出層的權(quán)值
b=ones（1hid_num）;%節(jié)點(diǎn)基寬
c=4.0*rand（in_numhid_num）;%基函數(shù)中心
h=zeros（1hid_num）;%RBF網(wǎng)絡(luò)的徑向基向量
eta_w=0.080;%學(xué)習(xí)率η
eta_b=0.010;
eta_c=0.010;
alpha=0.1;%動(dòng)量因子α
%beta=0.1;%變步長(zhǎng)調(diào)整系數(shù)
maxgen=20000;%最大迭代次數(shù)
error_goal=1e-4;%誤差精度要求
mse=zeros（1maxgen）;%均方誤差
wk_1=w;wk_2=w;%前1?2?時(shí)刻的權(quán)值
bk_1=b;bk_2=b;%前1?2?時(shí)刻的基寬值
ck_1=c;ck_2=c;%前1?2時(shí)刻的中心值
delta_c=zeros（in_numhid_num）;
%[~ce]=kmeans（train_ihid_num）;
%c=ce‘;
W=cell（maxgen1）;
B=cell（maxgen1）;
C=cell（maxgen1）;
%%?開始訓(xùn)練網(wǎng)絡(luò)
gen=0;
for?i=1:maxgen
????gen=gen+1;%當(dāng)前迭代次數(shù)
????%%?訓(xùn)練樣本
????for?j=1:size_tr
????????for?k=1:hid_num
????????????h（k）=exp（-norm（train_i（:j）-c（:k））^2/（2*b（k）^2））;
????????end
????????ye=w*h‘;
????????%修正權(quán)值
????????e=train_o（j）-ye;%當(dāng)前輸出誤差????????
????????for?k=1:hid_num????????????????
????????????w（k）=wk_1（k）+eta_w*e*h（k）+alpha