| 5 | 1/1 | 返回列表 |
| 查看: 1124 | 回復: 4 | ||
[求助]
跪求rbf程序調試
|
|
為什么我運行下面的程序,老提示我 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; %訓練樣本 TestSamNum=7; % Overlap=[0.25 0.1 0.2]; %輸入輸出矩陣 %訓練數(shù)據(jù) SamIn=train(:,1:3)'; SamOut=train(:,4:6)'; %預測數(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ù)中心聚類,并得到新的類內均值 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

木蟲 (正式寫手)
新蟲 (初入文壇)
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 083200學碩321分一志愿暨南大學求調劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 279分求調劑 一志愿211 +11 | chaojifeixia 2026-03-19 | 12/600 |
|
|
[考研] 311求調劑 +5 | 冬十三 2026-03-18 | 5/250 |
|
|
[考研] 22408 344分 求調劑 一志愿 華電計算機技術 +4 | solanXXX 2026-03-20 | 4/200 |
|
|
[考研] 330求調劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 316求調劑 +5 | 梁茜雯 2026-03-19 | 5/250 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
[考研] 招收調劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
[考研] 286分人工智能專業(yè)請求調劑愿意跨考! +3 | lemonzzn 2026-03-17 | 4/200 |
|
|
[考研] 一志愿中國海洋大學,生物學,301分,求調劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
[考研] 生物學調劑招人。! +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考研] 288求調劑,一志愿華南理工大學071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 266求調劑 +5 | 陽陽哇塞 2026-03-14 | 10/500 |
|
|
[考研] 286求調劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 334求調劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
[考研] 材料專碩326求調劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 一志愿,福州大學材料專碩339分求調劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
[考研] 288求調劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|
|
[考研] 本科南京大學一志愿川大藥學327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 一志愿哈工大材料324分求調劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|