| 5 | 1/1 | 返回列表 |
| 查看: 1123 | 回復(fù): 4 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
[求助]
跪求rbf程序調(diào)試
|
||
|
為什么我運行下面的程序,老提示我 In mean at 219 ,Divide by zero? 就是這條語句 Centers(:,i)=mean(SamIn(:,Index)')'; 求高人指點 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ù)測數(shù)據(jù) TestSamIn=test(:,1:3)'; TestSamOut=test(:,4:6)'; %節(jié)點數(shù) InDim=3; %輸入節(jié)點 ClusterNum=38; %隱含節(jié)點 %初始化中心 Centers=SamIn(:,1:ClusterNum); NumberInClusters=zeros(ClusterNum,3); IndexInClusters=zeros(ClusterNum,SamNum); %迭代求解均值中心 while 1 NumberInClusters=zeros(ClusterNum,3); IndexInClusters=zeros(ClusterNum,SamNum); %計算樣本點到中心的距離 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ā)表 | |
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 278求調(diào)劑 +6 | 煙火先于春 2026-03-17 | 6/300 |
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩,總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊 2026-03-18 | 6/300 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研] 316求調(diào)劑 +5 | 梁茜雯 2026-03-19 | 5/250 |
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生(接收計算機(jī)專業(yè)調(diào)劑) +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
[考研] 281求調(diào)劑(0805) +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 0703化學(xué)調(diào)劑 +10 | 妮妮ninicgb 2026-03-15 | 14/700 |
|
|
[考研] 0854可跨調(diào)劑,一作一項核心論文五項專利,省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|