| 查看: 1278 | 回复: 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]; %输入输出矩阵 %训练数据 SamIn=train(:,1:3)'; SamOut=train(:,4:6)'; %预测数据 TestSamIn=test(:,1:3)'; TestSamOut=test(:,4:6)'; %节点数 InDim=3; %输入节点 ClusterNum=38; %隐含节点 %初始化中心 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; %根据中心聚类，并得到新的类内均值 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 |
回复此楼» 猜你喜欢
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
青基如何破局
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
基金不中，共勉
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有7人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助matlab程序调试
已经有5人回复
- matlab程序调试
已经有7人回复
- MATLAB程序调试
已经有9人回复
- RBF神经网络，训练完成后，预测结果误差太大了，怎么回事呢？？？
已经有11人回复
- MATLAB RBF神经网络
已经有14人回复
- matlab程序出错调试
已经有7人回复
- 软件调试的艺术
已经有26人回复
- 程序调试 vc 6.0多文件工程里面为什么总是提示没有定义呢
已经有28人回复
- Matalab 程序调试
已经有6人回复
- 无法找到调试信息
已经有3人回复
- 求助一下程序调试过程出现的错误
已经有4人回复
- 螺旋波波头 程序调试 除零错误
已经有12人回复
- C语言调试时出现Unhandled exception 0xc0000005:Access Violation怎么回事
已经有4人回复
- matlab调试问题
已经有3人回复
- GDB单步调试问题
已经有8人回复
- 【求助】用C#编写的Windows应用程序中，用到了数组，不能成功调试
已经有4人回复
- 【求助】VS2005数值计算程序调试【急求！！！！加金币】
已经有9人回复
- 【求助】VC++数值算法集应用到MFC中报错，程序调试求助
已经有8人回复
- 【求助】VC++程序算法集调试报错
已经有8人回复
youth0826
至尊木虫 (著名写手)
weibo.com/138147022
- 应助: 172 (高中生)
- 金币: 10458.4
- 散金: 145
- 红花: 4
- 帖子: 1726
- 在线: 149.4小时
- 虫号: 707170
- 注册: 2009-02-23
- 专业: 模式识别
2楼2013-03-06 11:38:14
想飞的猪
木虫 (正式写手)
- 应助: 40 (小学生)
- 金币: 12225.4
- 散金: 1831
- 红花: 9
- 帖子: 584
- 在线: 156.7小时
- 虫号: 584434
- 注册: 2008-08-01
- 性别: GG
- 专业: 机器人学及机器人技术
3楼2013-03-06 17:47:01
4楼2013-05-08 10:31:25
5楼2015-03-27 21:32:45
投票:
