MATLAB实验五聚类方法与聚类有效性Word文档下载推荐.docx
《MATLAB实验五聚类方法与聚类有效性Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《MATLAB实验五聚类方法与聚类有效性Word文档下载推荐.docx(19页珍藏版)》请在冰豆网上搜索。
学生成绩
(分数)
实验目的和要求:
1.熟练掌握应用MATLAB软件计算谱系聚类与K均值聚类的命令。
2.熟练掌握模糊C均值聚类与模糊减法聚类的MATLAB实现。
3.掌握最优聚类数的理论及其实现。
实验方法和步骤(包括实验所需设备及软件名称):
所需设备:
笔记本电脑一台。
操作系统:
Windows8.1中文版软件名称:
MATLAB(R2016a)
1.计算各样品间的欧氏距离、马氏距离和加权平方距离
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)%类间距离为最短距离
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=[cophenet(z1,d1),cophenet(z2,d1),cophenet(z3,d1),cophenet(z4,d1),cophenet(z5,d1)]
[n,p]=size(x);
n1=n-1;
formatlong
c=zeros(n1-1,1);
forj=2:
n1
d1=pdist(x);
z2=linkage(d1,'
c=cluster(z2,j);
k=1;
ifk<
=j
b=find(c==k);
l=length(b)-1;
ifb>
a=x(b,:
c(j)=sum(l*var(a))+c(j);
end
end
R2=1-c/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=length(find(a==12))
[centeru]=fcm(x,3);
index1=find(u(1,:
)==max(u))
index1=find(u(2,:
index1=find(u(3,:
C=subclust(x,0.6)
运行程序,可以得出结果如下
d1=
1.0e+04*
1至5列
1.0548142973538041.0536055827765910.9693904595564160.9777118148365601.211184800199788
6至10列
0.3461552688866660.6883249145149400.4925463408249011.0292842654641140.733033330036909
11至15列
.0734********
16至20列
0.5834071512665920.9509511780533211.0878677862451850.8963324449343561.060307722921039
21至25列
1.1647116408854171.0992146647188620.9384023489138330.9953628265798361.157********3572
26至30列
1.1004932145565461.0940982673581021.0944067560834040.0333757898782940.173********3702
31至35列
0.0984267922011080.157********13471.3099663031734830.3895668683011940.712421721579290
36至40列
.0552********
41至45列
0.0841926099310380.0630093894510970.6475464828126240.159********58350.110007710938825
46至50列
0.2194672460710250.0719716903650320.1157249761028280.0823835893306430.465400245533885
51至55列
0.1206152562033510.151********55820.0675024158160280.1026581297608720.204666139625000
56至60列
0.152********13730.1179988077355020.1654842649468521.3030696899022710.394046122033957
61至65列
0.7085418551066690.0303080753265530.3921396353303250.0687788751143840.397322354527152
66至70列
.021*********
71至75列
.0853********
76至80列
0.4354086639181630.0940855129709140.1291570459285900.0538039409801920.097501830905886
81至85列
0.1726965841729360.176********34890.2953545592368601.1936114887185870.326797222397927
86至90列
0.5906985425045500.1283552173345520.2682577714624500.1263305566005310.268093134037409
91至95列
0.1409927551011040.0975691222364950.134********97310.5001841452485270.039295253301130
96至100列
0.131********30930.1053946908292820.1298793841030980.2686241203894390.187********7936
101至105列
0.3209704927092210.1182625393182470.2192761594930010.1933482546443080.154875073640015
106至110列
0.150********45210.2389081801571471.2399697523927750.3022076553216350.637409*********
111至115列
0.1151763955374540.3405810919531500.1712308319987960.3788144870962570.120340580375034
116至120列
0.132********84270.1087998640991800.5912668714328580.155********79730.180********5762
121至125列
0.2105402848031700.1465916191635800.2049961258780270.157********74730.481534817258316
126至130列
0.159********80820.2458654030440230.1640936451542230.1546132503345040.269272797965186
131至135列
1.4650382328099840.5389810452965860.8577890025600700.1940149139344700.543588674688868
136至140列
0.180********74860.5569315302943440.1652851090177210.2095369383044430.172477754235148
141至145列
0.7901699504524330.2952092071904260.1743475489561010.3601919009458710.175********8806
146至150列
0.0656301490399040.1407853692753620.5651359600750960.2505274391279330.144672404300198
151至155列
0.129667*********0.1546532849020670.2559164102553800.9512767696080880.668288245898430
156至160列
1.2763056565141440.9463582478041811.3089618896682980.9390821052203051.302481353894174
161至165列
1.2647622419099171.2981450904028410.7258148351756121.182********68541.321279812330454
166至170列
1.1241520515010421.3014828111219911