MATLAB实验五聚类方法与聚类有效性.docx
《MATLAB实验五聚类方法与聚类有效性.docx》由会员分享，可在线阅读，更多相关《MATLAB实验五聚类方法与聚类有效性.docx（18页珍藏版）》请在冰豆网上搜索。
MATLAB实验五聚类方法与聚类有效性
广西财经学院学生实验报告
实验课程名称
统计软件及应用
课程编码
1022201630
实验课程类型
基础实验（）专业实验（√）综合实训（）
实验项目名称
实验5聚类方法与聚类有效性
实验时间
实验地点
学生班级
指导老师
学生姓名
学号
学生成绩
（分数）
实验目的和要求：
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,'complete'）;
c=cluster（z2,j）;
k=1;
ifk<=j
b=find（c==k）;l=length（b）-1;
ifb>0
a=x（b,:
）;
c（j）=sum（l*var（a））+c（j）;
end
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,:
）==max（u））
index1=find（u（3,:
）==max（u））
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.4219205860884071.3477039465516901.091949050370483
171至175列
1.2360089327577701.3931673639749821.3496103650457791.3329448535697191.308012777795003
176至180列
0.3619474469449950.3737362994010620.183********94350.4206773715568740.268221126735013
181至185列
0.3933256356633770.3709075412646660.3804747923778920.3637044180306310.307659057414210
186至190列
0.4309383276479360.2923147669311290.4049227598814370.5032078308909750.436159152421453
191至195列
0.4959536228751640.3650133397713020.5125071603860770.4453637340208560.419443729539494
196至200列
0.4723524199038260.6828453427416780.3238338767192220.7113897960401180.389141099020651
201至205列
0.7038473387958500.6688651110104340.6867661844856660.188********85790.583979714304358
206至210列
0.7162930941688600.5543206175544980.7051662081367770.8276560899564510.738249128890106
211至215列
0.6040590153163850.6616530847166060.8066117528594530.7604533608019630.716200330664542
216至220列
0.7254663529695920.3652729912695430.0665413550959100.368504*********0.036125213937083
221至225列
.0354********
226至230列
0.1699846591136980.0395853165959300.151********60750.0926948598736740.410884559740324
231至235列
0.0704184559259860.1428437380181570.0786998375284220.1085609229695470.163882274584532
236至240列
0.3889550258641740.1443053398388290.3864606720780260.3484259637369180.371371989686891
241至245列
0.2587777633453080.2616904434135870.3938557225672870.2411148725130820.384445310876593
246至250列
0.5132330559960070.4215405843640680.3641133989130310.3445061586241960.483744175487416
251至255列
0.4427894278593380.4008911554536970.4046479836351590.3837395360697150.057648999835210
256至260列
.0512********
261至265列
0.1900449343655340.0287059511251590.152********27800.0934914235264390.391714238933180
266至270列
0.0970180890298300.0981832671232730.0805060798511520.0870405570179790.107824153792181
271至275列
0.3931153415131490.3479024700602170.3926830540957940.2595500016528610.263568578760823
276至280列
0.3946868988261460.2038788291633050.3832770539974970.5190171475558400.439475726744948
281至285列
0.2441873866070890.3231293647070780.4650554722632560.4394049473754250.420855311294749
286至290列
0.3698136051932110.0487897156683660.0449190930117690.6293299939673940.133********7578
291至295列
.0715********
296至300列
0.4299322157189430.0990997493286440.1277225474025630.0615434570608440.078356027324004
301至305列
0.1629921221133100.0674844245881970.5873686646255480.0892001746074530.070018094097169
306至310列
.0373********
311至315列
.0661********
316至320列
0.6187926263216130.131********81430.0778115766001430.2076638730954430.064298249066363
321至325列
0.151********84710.0829504986965120.4425758742679050.1296465827316710.156********2594
326至330列
.0990********
331至335列
0.4556132966178660.6237699796808110.7583659864570670.6707919458803900.438091151887367
336至340列
0.5730508128028440.7127180542486630.6824125099923070.6438038042051320.613768478690459
341至345列
0.138********27010.0860075421402100.1248891169998410.2616557506209260.178********4348
346至350列
.0977********
351至355列
0.1683696963648750.2069838768261920.0495471253051070.155********24010.087561318166186
356至360列
0.4008568275893530.1224184592657500.1028055650293310.0952334755587550.073619402591980
361至365列
0.1103830356486000.1858647545932260.3176741579291590.2487442339050300.277842430541125
366至370列
0.1192896945632770.2692045608510380.2375235009993750.2413761311024770.201163372737186
371至375列
0.1408442328815770.0829652488515520.4051893687906930.0852313985395050.109984209262057
376至380列
0.0656867028476850.0882407416899930.133********11980.1163559226038790.532430895066205
381至385列
0.2025046367691370.1218936193531060.0820789524969710.152********45510.236004846439644
386至390列
0.468048*********0.152********34200.132********52290.0896320357126850.100460646369611
391至395列
0.189********65200.3533579794174740.4330760741324790.4543365498086190.453273975819261
396至400列
0.3153847937171350.1627089182620300.1227610252889740.1653009194076060.148312946036413
401至405列
.0890********
406列
0.171888599182726
R=
0.9436312853835230.9169606699394090.9398331718148770.9376606808098590.747075679381588
optimaln=
12
a=3
7
5
10
8
2
9
6
11
5
6
5
6
5
5
5
4
10
5
1
5
2
7
12
5
7
7
7
5
b=
1.0e+04*
1.0957620000000000.0788260000000000.0205940000000000.326592000000000
0.7602275000000000.1006115000000000.0116865000000000.359959000000000
1.8738960000000000.0778360000000000.0452750000000000.770787000000000
1.5186390000000000.2405920000000000.0701250000000000.338295000000000
0
升级会员 