• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    迪恩网络公众号

解决MATLAB在K-means提示“KMEANSdoesnotacceptcomplexdata.”错误

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

昨天开始看聚类方法,结果才看到第一个K-means聚类方法就卡壳了。我先大致了解了K-means聚类方法的原理,然后照着老师的代码逐步实现,就在使用kmeans函数进行聚类时,命令窗口报出了“KMEANS does not accept complex data.”的错误,百度也没有查询到相关的解决办法。自己前后经历了反复看原理以及kmeans函数的使用方法,但也没什么进展。直到后面突然想起看报告的错误的源代码,不到五分钟解决了问题。以下为解决方案:

(1)命令窗口报出“错误使用 kmeans (line 159) KMEANS does not accept complex data.”错误,点击“line 159”,显示为以下错误,意思是参数类型不匹配。

(2)我的kmeans聚类的那两行代码是:

bonds = corp(:,{'Coupon','YTM','CurrentYield','RatingNum'});%取出含有'Coupon','YTM','CurrentYield','RatingNum'的列
kidx = kmeans(bonds,numClust,'distance',dist_k);

通过在命令窗口输入“class(bonds)”发现竟然是table类型,肯定不能识别啊,于是我使用table2array()函数将其改成了数组类型,一切问题迎刃而解。

改后的代码是:

bonds = corp(:,{'Coupon','YTM','CurrentYield','RatingNum'});%取出含有'Coupon','YTM','CurrentYield','RatingNum'的列
% bonds = table2array(bonds);
%%%%%!!!!!!!注意:!!!!!!
%%%%%问题就在这,必须将之前table类型的bonds转化为array类型
%%%%%否则将会报“KMEANS does not accept complex data.”这样的错误
%设置类别数量
kidx = kmeans(bonds,numClust,'distance',dist_k);

(3)通过这个问题,才领悟到一种新的改错方式,就是查看错误报告。看来是之前老师教授的知识没有吸收啊,只有遇到问题通过不断尝试才能够真正学会。

下面介绍一下这个K-means聚类MATLAB实例:

【题目背景】

 

下面是MATLAB源代码介绍:

(1)导入数据和预处理数据

 1 load 'BondData.mat';
 2 settle = floor(date);%floor朝无穷大方向取整
 3 %数据预处理
 4 bondData.MaturityN = datenum(bondData.Maturity,'dd-mmm-yyyy');%将时间化为数值
 5 bondData.SettleN = settle*ones(height(bondData),1);
 6 %筛选数据
 7 corp = bondData(bondData.MaturityN>settle&...
 8                 bondData.Type == 'Corp'&...
 9                 bondData.Rating >='CC'&...
10                 bondData.YTM<30&...
11                 bondData.YTM>=0,:);
12 %设置随机数生成方式,保证结果可重现
13 rng('default');

 

(2)探索数据

 1 figure
 2 gscatter(corp.Coupon,corp.YTM,corp.Rating)
 3 set(gca,'lineWidth',2);
 4 xlabel('票面利率')
 5 ylabel('到期收益率')
 6 %选择聚类变量
 7 corp.RatingNum = double(corp.Rating);
 8 bonds = corp(:,{'Coupon','YTM','CurrentYield','RatingNum'});%取出含有'Coupon','YTM','CurrentYield','RatingNum'的列
 9 bonds = table2array(bonds);%!!!!!!!!
10 %%%%%!!!!!!!注意:!!!!!!
11 %%%%%问题就在这,必须将之前table类型的bonds转化为array类型
12 %%%%%否则将会报“KMEANS does not accept complex data.”这样的错误
13 %设置类别数量
14 numClust = 3;
15 %设置用于可视化聚类效果的变量
16 VX = [corp.Coupon,double(corp.Rating),corp.YTM]; 

(3)K-means聚类

 1 dist_k = 'cosine';
 2 kidx = kmeans(bonds,numClust,'distance',dist_k);
 3 %绘制聚类效果图
 4 figure
 5 F1 = plot3(VX(kidx==1,1),VX(kidx==1,2),VX(kidx==1,3),'r*',...
 6            VX(kidx==2,1),VX(kidx==2,2),VX(kidx==2,3),'bo',...
 7            VX(kidx==3,1),VX(kidx==3,3),VX(kidx==3,3),'kd');
 8 set(gca,'lineWidth',2);
 9 grid on;
10 set(F1,'lineWidth',2);
11 xlabel('票面利率','fontsize',12);
12 ylabel('评级得分','fontsize',12);
13 ylabel('到期收益率','fontsize',12);
14 title('K-means聚类')
15 %评估各类别的相关程度
16 dist_metric_k = pdist(bonds,dist_k);
17 dd_k = squareform(dist_metric_k);
18 [~,idx]  = sort(kidx);
19 dd_k = dd_k(idx,idx);
20 figure
21 imagesc(dd_k)
22 set(gca,'linewidth',2);
23 xlabel('数据点','fontsize',12)
24 ylabel('数据点','fontsize',12)
25 title('K-means聚类结果相关程度图','fontsize',12);
26 ylabel(colorbar,['距离矩阵:',dist_k])
27 axis square

 运行结果图如下:

figure1:

figure2:

figure3:

本题中,还遇到了一个小问题,MATLAB无法识别“Figure”,绘图应为“figure”,可对其进行编号,不编号则是会紧接着新弹出一个绘图界面。

 


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
delphi安装一个ActiveX控件发布时间:2022-07-18
下一篇:
Delphi-条件判断那些事发布时间:2022-07-18
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap