资源描述:
《MATLAB实验五聚类方法与聚类有效性.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、广西财经学院学生实验报告实验课程名称统计软件及应用课程编码实验课程类型基础实验()专业实验(√)综合实训()实验项目名称实验5聚类方法与聚类有效性实验时间实验地点学生班级指导老师学生姓名学号学生成绩(分数)实验目的和要求:1.熟练掌握应用MATLAB软件计算谱系聚类与K均值聚类的命令。2.熟练掌握模糊C均值聚类与模糊减法聚类的MATLAB实现。3.掌握最优聚类数的理论及其实现。实验方法和步骤(包括实验所需设备及软件名称):所需设备:笔记本电脑一台。操作系统:Windows8.1中文版软件名称:MATLAB(R2016a)1.计算各样品间的欧氏距
2、离、马氏距离和加权平方距离2.运用谱系聚类法进行聚类,包括确定最优聚类数,选择合适的类间距离,同时作出谱系图。3.运用K均值聚类法进行聚类。4.运用模糊C均值和模糊减法聚类进行聚类。5.综合分析以上不同的聚类法所得的聚类结果,能得到什么样的结论。实验内容和数据记录:输入数据x=[18738.96778.36...1976.49]运行如下的matlab程序d1=pdist(x)%欧氏距离d2=pdist(x,'seuclidean');%加权平方距离d3=pdist(x,'mahalanobis');%马氏距离z1=linkage(d1)%类间距
3、离为最短距离figure(1)%谱系聚类图H=dendrogram(z1)z2=linkage(d1,'complete')%类间距离为最长距离figure(2)H2=dendrogram(z2)z3=linkage(d1,'average')%类间距离为类平均距离figure(3)H3=dendrogram(z3)z4=linkage(x,'centroid')%重心距离figure(4)H4=dendrogram(z4)z5=linkage(x,'ward')%离差平方和figure(5)H5=dendrogram(z5)R=[cophen
4、et(z1,d1),cophenet(z2,d1),cophenet(z3,d1),cophenet(z4,d1),cophenet(z5,d1)][n,p]=size(x);n1=n-1;formatlongc=zeros(n1-1,1);forj=2:n1d1=pdist(x);z2=linkage(d1,'complete');c=cluster(z2,j);k=1;ifk<=jb=find(c==k);l=length(b)-1;ifb>0a=x(b,:);c(j)=sum(l*var(a))+c(j);endendendR2=1-c/
5、sum(n1*var(x));optimaln=find(R2==max(R2))[ab]=kmeans(x,12)n1=length(find(a==1))n2=length(find(a==2))n3=length(find(a==3))n4=length(find(a==4))n5=length(find(a==5))n6=length(find(a==6))n7=length(find(a==8))n9=length(find(a==9))n10=length(find(a==10))n11=length(find(a==11))n12
6、=length(find(a==12))[centeru]=fcm(x,3);index1=find(u(1,:)==max(u))index1=find(u(2,:)==max(u))index1=find(u(3,:)==max(u))C=subclust(x,0.6)运行程序,可以得出结果如下d1=1.0e+04*1至5列1.38041.65910.64160.65601.97886至10列0.66660.49400.49011.41140.690911至15列1.02990.41311.42471.30451.287416至20列0.6
7、5920.33211.51850.43561.103921至25列1.54171.88620.38330.98361.357226至30列1.65461.81021.34040.82940.370231至35列0.11080.13471.34830.11940.929036至40列0.72830.65600.16510.93580.920041至45列0.10380.10970.26240.58350.882546至50列0.10250.50320.28280.06430.388551至55列0.33510.55820.60280.08720.
8、500056至60列0.13730.55020.68521.22710.395761至65列0.66690.65530.03250.43840.715