| 5 | 1/1 | 返回列表 |
| 查看: 1125 | 回復(fù): 4 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
跪求rbf程序調(diào)試
|
||
|
為什么我運(yùn)行下面的程序,老提示我 In mean at 219 ,Divide by zero? 就是這條語(yǔ)句 Centers(:,i)=mean(SamIn(:,Index)')'; 求高人指點(diǎn) clc clear data=xlsread('date.xls'); train=data(1:70, ;test=data(71:77, ;SamNum=70; %訓(xùn)練樣本 TestSamNum=7; % Overlap=[0.25 0.1 0.2]; %輸入輸出矩陣 %訓(xùn)練數(shù)據(jù) SamIn=train(:,1:3)'; SamOut=train(:,4:6)'; %預(yù)測(cè)數(shù)據(jù) TestSamIn=test(:,1:3)'; TestSamOut=test(:,4:6)'; %節(jié)點(diǎn)數(shù) InDim=3; %輸入節(jié)點(diǎn) ClusterNum=38; %隱含節(jié)點(diǎn) %初始化中心 Centers=SamIn(:,1:ClusterNum); NumberInClusters=zeros(ClusterNum,3); IndexInClusters=zeros(ClusterNum,SamNum); %迭代求解均值中心 while 1 NumberInClusters=zeros(ClusterNum,3); IndexInClusters=zeros(ClusterNum,SamNum); %計(jì)算樣本點(diǎn)到中心的距離 for i=1:SamNum AllDistance=dist(Centers',SamIn(:,i)); [MinDist,Pos]=min(AllDistance); NumberInClusters(Pos)=NumberInClusters(Pos)+1; IndexInClusters(Pos,NumberInClusters(Pos))=i; end OldCenters=Centers; %根據(jù)中心聚類,并得到新的類內(nèi)均值 for i=1:ClusterNum Index=IndexInClusters(i,1:NumberInClusters(i)); Centers(:,i)=mean(SamIn(:,Index)')'; end %判斷中心有沒有變化 EqualNum=sum(sum(Centers==OldCenters)); if EqualNum==InDim*ClusterNum break end end AllDistances=dist(Centers',Centers); Maximun=max(max(AllDistance)); for i=1:ClusterNum AllDistances(i,i)=Maximun+1; end Spreads=Overlap*min(AllDistances)'; Distance=dist(Centers',Samln); SpreadsMat_l=repmat(Spreads(:,1),1,SamNum); SpreadsMat_2=repmat(Spreads(:,2),1,SamNum); SpreadsMat_3=repmat(Spreads(:,3),1,SamNum); HiddenUnitOut_1=radbas(Distance./SpreadsMat_1); HiddenUnitOut_2=radbas(Distance./SpreadsMat_2); HiddenUnitOut_3=radbas(Distance./SpreadsMat_3); HiddenUnitOutEx_1=[HiddenUnitOut_l' ones(SamNum,1)]'; HiddenUnitOutEx_2=[HiddenUnitOut_2' ones(SamNum,1)]'; HiddenUnitOutEx_3=[HiddenUnitOut_3' ones(SamNum,1)]'; W2Ex_1=SamOut(1, *pinv(HiddenUnitOutEx_1);W2Ex_2=SamOut(2, *pinv(HiddenUnitOutEx_2);W2Ex_3=SamOut(3, *pinv(HiddenUnitOutEx_3);W2_1=W2Ex_l(:,1:ClusterNum); B2_1=W2Ex_l(:,ClusterNum+1); W2_2=W2Ex_2(:,1:ClusterNum); B2_2=W2Ex_2(:,ClusterNum+1); W2_3=W2Ex_3(:,1:ClusterNum); B2_3=W2Ex_3(:,ClusterNum+1); TestDistance=dist(Centers',TestSamln); TestSpreadsMat_l=repmat(Spreads(:,1),1,TestSamNum); TestSpreadsMat_2=repmat(Spreads(:,2),1,TestSamNum); TestSpreadsMat_3=repmat(Spreads(:,3),1,TestSamNum); TestHiddenUnitOut_1=radbas(TestDistance./TestSpreadsMat_1); TestHiddenUnitout_2=radbas(TestDistance./TestSpreadsMat_2); TestHiddenUnitOut_3=radbas(TestDistance./TestSpreadsMat_3); TestNNOut_1=W2_1*TestHiddenUnitOut_1+repmat(B2_1,[1,TestSamNum]); TestNNOut_2=W2_2*TestHiddenUnitOut_2+repmat(B2_2,[1,TestSamNum]); TestNNOut_3=W2_3*TestHiddenUnitOut_3+repmat(B2_3,[l,TestSamNurn]); TestNNOut(1, =TestNNOut_l;TestNNOut(2, =TestNNOut_2;TestNNOut(3, =TestNNOut_3;plot(TestSamIn,TestNNOut(1, ,'k-');hold on |
新蟲 (初入文壇)
至尊木蟲 (著名寫手)
weibo.com/138147022

木蟲 (正式寫手)
新蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +8 | 申子申申 2026-03-19 | 13/650 |
|
|
[考研] 一志愿華中農(nóng)業(yè)071010,總分320求調(diào)劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級(jí)已過,有科研經(jīng)歷 +13 | 曦熙兮 2026-03-15 | 13/650 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生(接收計(jì)算機(jī)專業(yè)調(diào)劑) +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院,一志愿華科070200學(xué)碩,346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|