客服电话
APP下载
迪恩网络APP
随时随地掌握行业动态
官方微信
扫描二维码
关注迪恩网络微信公众号
|
昨天开始看聚类方法，结果才看到第一个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”，可对其进行编号，不编号则是会紧接着新弹出一个绘图界面。
|
2023-10-27
2022-08-15
2022-08-17
2022-09-23
2022-08-13
六六分期app的软件客服如何联系？不知道吗？加qq群【895510560】即可！标题：六六分期
今天小编告诉大家如何处理win10系统火狐flash插件总是崩溃的问题，可能很多用户都不知
今天小编告诉大家如何对win10系统删除桌面回收站图标进行设置，可能很多用户都不知道
今天小编告诉大家如何对win10系统电脑设置节能降温的设置方法，想必大家都遇到过需要
我们在使用xp系统的过程中,经常需要对xp系统无线网络安装向导设置进行设置，可能很多
今天小编告诉大家如何处理win7系统玩cf老是与主机连接不稳定的问题，可能很多用户都不
电脑对日常生活的重要性小编就不多说了，可是一旦碰到win7系统设置cf烟雾头的问题，很
我们在日常使用电脑的时候，有的小伙伴们可能在打开应用的时候会遇见提示应用程序无法
今天小编告诉大家如何对win7系统打开vcf文件进行设置，可能很多用户都不知道怎么对win
今天小编告诉大家如何对win10系统s4开启USB调试模式进行设置，可能很多用户都不知道怎
请发表评论