MATLAB实验五聚类方法与聚类有效性Word文档格式.docx
- 文档编号:18234972
- 上传时间:2022-12-14
- 格式:DOCX
- 页数:18
- 大小:71.31KB
MATLAB实验五聚类方法与聚类有效性Word文档格式.docx
《MATLAB实验五聚类方法与聚类有效性Word文档格式.docx》由会员分享,可在线阅读,更多相关《MATLAB实验五聚类方法与聚类有效性Word文档格式.docx(18页珍藏版)》请在冰豆网上搜索。
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.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
4
1
b=
1.0957620000000000.0788260000000000.0205940000000000.326592000000000
0.7602275000000000.1006115000000000.0116865000000000.359959000000000
1.8738960000000000.0778360000000000.0452750000000000.770787000000000
1.5186390000000000.2405920000000000.0701250000000000.338295000000000
0
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB 实验 五聚类 方法 有效性