MATLAB實驗五聚類方法與聚類有效性

《MATLAB實驗五聚類方法與聚類有效性》由會員分享，可在線閱讀，更多相關《MATLAB實驗五聚類方法與聚類有效性（27頁珍藏版）》請在裝配圖網上搜索。

1、廣西財經學院學生實驗報告實驗課程名稱統計軟件及應用課程編碼1022201630實驗課程類型基礎實驗（）專業實驗（V ）綜合實訓（）實驗項目名稱實驗5聚類方法與聚類有效性實驗時間實驗地點學生班級指導老師學生姓名學號學生成績（分數）實驗目的和要求：1. 熟練掌握應用MATLAB軟件計算譜系聚類與K均值聚類的命令。2. 熟練掌握模糊C均值聚類與模糊減法聚類的MATLAB實現。3. 掌握最優聚類數的理論及其實現。實驗方法和步驟(包括實驗所需設備及軟件名稱)：所需設備：筆記本電腦一臺。操作系統：Windows8.1中文版軟件名稱:MATLAB (R2016a)1. 計算各樣品間的歐氏距離、馬氏距離和加權

2、平方距離2. 運用譜系聚類法進行聚類，包括確定最優聚類數，選擇合適的類間距離，同時作 出譜系圖。3. 運用K均值聚類法進行聚類。4. 運用模糊C均值和模糊減法聚類進行聚類。5. 綜合分析以上不同的聚類法所得的聚類結果，能得到什么樣的結論。實驗內容和數據記錄：輸入數據 x=18738.96 778.36 1976.49 運行如下的matlab程序d1=pdist(x) %歐氏距離 d2=pdist(x,seuclidean); % 加權平方距離 d3=pdist(x,mahalanobis); %馬氏距離zl=linkage(dl)%類間距離為最短距離figure(l)%譜系聚類圖H=dendr

3、ogram(zl)z2=linkage(dl,complete)%類間距離為最長距離figure(2)H2=dendrogram(z2)z3=linkage(dl,average)%類間距離為類平均距離figure(3)H3=dendrogram(z3)z4=linkage(x,centroid)%重心距離figure(4)H4=dendrogram(z4) figure(5)z5=linkage(x,ward)%離差平方和H5=dendrogram(z5)R=cophenet(zl,dl),cophenet(z2,dl),cophenet(z3,dl),cophenet(z4,dl),cop

4、henet(z5,dl) n,p=size(x);nl=n-l;format longc=zeros(nl-l,l);for j=2:nldl=pdist(x);z2=linkage(dl,complete);c=cluster(z2,j);k=l;if k0a=x(b,:); c(j)=sum(l*var(a)+c(j);endendendR2=l-c/sum(nl*var(x); optimaln=find(R2=max(R2)a b=kmeans(x,12)nl=length(find(a=l)n2=length(find(a=2)n3=length(find(a=3)n4=length

5、(find(a=4)n5=length(find(a=5)n6=length(find(a=6)n7=length(find(a=8)n9=length(find(a=9)nl0=length(find(a=10)nll=length(find(a=ll)nl2=length(find(a=12)center u=fcm(x,3);indexl=find(u(l,:)=max(u)indexl=find(u(2,:)=max(u)index1=find(u(3,:)=max(u)C=subclust(x,0.6)運行程序，可以得出結果如下d1 =1.0e+04 *1至5列1.054814297

6、3538041.0536055827765910.9777118148365601.2111848001997886至10列0.3461552688866660.6883249145149401.0292842654641140.73303333003690911至15列1.0734700415502990.7439237227841311.0249571659230451.05298184271287416至20列0.9693904595564160.4925463408249011.0562140917342471.0878677862451850.5834071512665920.950

7、9511780533210.8963324449343561.06030772292103921至25列1.1647116408854171.0992146647188620.9953628265798361.15752300829357226至30列1.1004932145565461.0940982673581020.0333757898782940.17364455184370231至35列0.0984267922011080.1573129298913470.3895668683011940.71242172157929036至40列0.0552813312972830.4003368

8、497465600.4146914371493580.04417580167920041至45列0.0841926099310380.0630093894510970.1594750478758350.1100077109388250.9384023489138331.0944067560834041.3099663031734830.0985303144316510.64754648281262446至50列0.2194672460710250.0719716903650320.0823835893306430.46540024553388551至55列0.1206152562033510.

9、1513772847655820.1026581297608720.20466613962500056至60列0.1526726632013730.1179988077355021.3030696899022710.39404612203395761至65列0.7085418551066690.0303080753265530.0687788751143840.39732235452715266至70列0.0215304520156920.0569854592418100.6354443718760910.14055482251064971至75列0.0853751838533890.1986

10、353631255020.1157249761028280.0675024158160280.1654842649468520.3921396353303250.0606266861291300.0409643256993200.1230715472723080.08078875539949976至80列0.4354086639181630.0940855129709140.0538039409801920.09750183090588681至85列0.1726965841729360.1768896254334891.1936114887185870.32679722239792786至90

11、列0.5906985425045500.1283552173345520.1263305566005310.26809313403740991至95列0.1409927551011040.0975691222364950.5001841452485270.03929525330113096至100列0.1317136776230930.1053946908292820.2686241203894390.187087534517936101至105列0.3209704927092210.1182625393182470.1291570459285900.2953545592368600.2682

12、577714624500.1346713267997310.1298793841030980.2192761594930010.1933482546443080.154875073640015106至110列0.1506598186445210.2389081801571471.2399697523927750.3022076553216350.637409200317190111至115列0.1151763955374540.3405810919531500.1712308319987960.3788144870962570.120340580375034116至120列0.13222515

13、45584270.1087998640991800.5912668714328580.1556672164779730.180389627675762121至125列0.2105402848031700.1465916191635800.2049961258780270.1577459629974730.481534817258316126至130列0.1597791375680820.2458654030440230.1640936451542230.1546132503345040.269272797965186131至135列1.4650382328099840.538981045296

14、5860.8577890025600700.1940149139344700.543588674688868136至140列0.1800663830174860.5569315302943440.2095369383044430.172477754235148141至145列0.7901699504524330.2952092071904260.3601919009458710.175762002898806146至150列0.0656301490399040.1407853692753620.2505274391279330.144672404300198151至155列0.12966728

15、02290540.1546532849020670.9512767696080880.668288245898430156至160列1.2763056565141440.9463582478041810.9390821052203051.3024813538941740.1652851090177210.1743475489561010.5651359600750960.2559164102553801.308961889668298161至165列1.2647622419099171.2981450904028410.7258148351756121.1827290327268541.321

16、279812330454166至170列1.1241520515010421.3014828111219911.4219205860884071.3477039465516901.091949050370483171至175列1.2360089327577701.3931673639749821.3496103650457791.3329448535697191.308012777795003176至180列0.3619474469449950.3737362994010620.1833081278094350.4206773715568740.268221126735013181至185列0

17、.3933256356633770.3709075412646660.3804747923778920.3637044180306310.307659057414210186至190列0.4309383276479360.2923147669311290.4049227598814370.5032078308909750.436159152421453191至195列0.4959536228751640.3650133397713020.5125071603860770.4453637340208560.419443729539494196至200列0.4723524199038260.682

18、8453427416780.3238338767192220.7113897960401180.389141099020651201至205列0.7038473387958500.6688651110104340.6867661844856660.1885853909585790.583979714304358206至210列0.7162930941688600.5543206175544980.7051662081367770.8276560899564510.738249128890106211至215列0.6040590153163850.6616530847166060.8066117

19、528594530.7604533608019630.716200330664542216至220列0.7254663529695920.3652729912695430.0665413550959100.3685043901271730.036125213937083221至225列0.6076451984135150.0354068569489020.0658342142126720.1125738703740790.085724986404198226至230列0.1517251394660750.1699846591136980.0395853165959300.09269485987

20、36740.410884559740324231至235列0.0786998375284220.0704184559259860.1428437380181570.1085609229695470.163882274584532236至240列0.3864606720780260.3889550258641740.1443053398388290.3484259637369180.371371989686891241至245列0.3938557225672870.2587777633453080.2616904434135870.2411148725130820.384445310876593

21、246至250列0.3641133989130310.5132330559960070.4215405843640680.4046479836351590.6245656491282560.1520647879227800.0805060798511520.3926830540957940.3445061586241960.483744175487416251至255列0.4427894278593380.4008911554536970.3837395360697150.057648999835210256至260列0.0512226184317050.0786732596439220.12

22、88500048932870.028235542495231261至265列0.1900449343655340.0287059511251590.0934914235264390.391714238933180266至270列0.0970180890298300.0981832671232730.0870405570179790.107824153792181271至275列0.3931153415131490.3479024700602170.2595500016528610.263568578760823276至280列0.3946868988261460.203878829163305

23、0.3832770539974970.5190171475558400.439475726744948281至285列0.2441873866070890.3231293647070780.4650554722632560.4394049473754250.420855311294749286至290列0.3698136051932110.0487897156683660.0449190930117690.6293299939673940.133767935937578291至295列0.0715131513918940.1970848846055930.0332291342649790.13

24、02039968203740.082543809222739296至300列0.4299322157189430.0990997493286440.1277225474025630.0615434570608440.078356027324004301至305列0.1629921221133100.0674844245881970.5873686646255480.0892001746074530.070018094097169306至310列0.1520630974793030.0373698232133900.1080304463473150.383838054463077311至315列

25、0.0661008495709400.1430465281088640.1010890042338930.138034152259504316至320列0.6187926263216130.1317404499081430.2076638730954430.064298249066363321至325列0.1510705776284710.0829504986965120.1296465827316710.156881851002594326至330列0.0990324355148350.0531717164571540.4993352104208150.630932674244249331至

26、335列0.4556132966178660.6237699796808110.1740009468451250.0973853481279400.0778115766001430.4425758742679050.1800958084687150.7583659864570670.6707919458803900.438091151887367336至340列0.5730508128028440.7127180542486630.6824125099923070.6438038042051320.613768478690459341至345列0.1389082088827010.086007

27、5421402100.1248891169998410.2616557506209260.178825537214348346至350列0.3311920610310580.0977826839783000.2227851798437230.1858805434627310.164190880992216351至355列0.1683696963648750.2069838768261920.0495471253051070.1558865495224010.087561318166186356至360列0.4008568275893530.1224184592657500.1028055650

28、293310.0952334755587550.073619402591980361至365列0.1103830356486000.1858647545932260.3176741579291590.2487442339050300.277842430541125366至370列0.1192896945632770.2692045608510380.2413761311024770.201163372737186371至375列0.1408442328815770.0829652488515520.0852313985395050.109984209262057376至380列0.065686

29、7028476850.0882407416899930.1163559226038790.532430895066205381至385列0.2025046367691370.1218936193531060.1521713227845510.236004846439644386至390列0.4680485103042210.1526844051434200.0896320357126850.1004606463696110.2375235009993750.4051893687906930.1333660817411980.0820789524969710.132390588555229391

30、至395列0.1890113320465200.3533579794174740.4330760741324790.4543365498086190.453273975819261396至400列0.3153847937171350.1627089182620300.1227610252889740.1653009194076060.148312946036413401至405列0.0890409706034250.1464867294501450.1279582342407080.1132834078804130.170858945039468406列0.171888599182726R =

31、0.9436312853835230.9169606699394090.9398331718148770.9376606808098590.747075679381588optimaln =12a =3751081.0e+04 *1.0957620000000000.3265920000000000.7602275000000000.3599590000000001.8738960000000000.7707870000000001.5186390000000000.3382950000000000.9235123000000000.3223403000000001.2643100000000

32、000.4165040000000000.8646406000000000.3361372000000000.9494590000000000.4610320000000002.1791110000000000.6199770000000001.0302815000000000.3173745000000000.0788260000000000.1006115000000000.0778360000000000.2405920000000000.1073434000000000.1793046666666670.1100680000000000.1483300000000000.1399140

33、000000000.1434855000000000.0205940000000000.0116865000000000.0452750000000000.0701250000000000.0243162000000000.0535706666666670.0274496000000000.0248040000000000.0369120000000000.038289500000000nil =1n12 =1Iteration count = 1, obj. fen = 155508663.694128Iteration count = 2, obj. fcn = 116346760.288

34、898Iteration count = 3, obj. fcn = 95702808.742790Iteration count = 4, obj. fcn = 78436271.654292Iteration count = 5, obj. fcn = 68565426.886294Iteration count = 6, obj. fcn = 62591744.851415Iteration count = 7, obj. fcn = 58038028.433253Iteration count = & obj. fcn = 53529919.373342Iteration count

35、= 9, obj. fcn = 49298407.252189Iteration count = 10, obj. fcn = 46451107.925572Iteration count = 11, obj. fcn = 45148844.219200Iteration count = 12, obj. fcn = 44685340.250062Iteration count = 13, obj. fcn = 44538617.475096Iteration count = 14, obj. fcn = 44494374.468347Iteration count = 15, obj. fc

36、n = 44481304.303276Iteration count = 16, obj. fcn = 44477478.128347Iteration count = 17, obj. fcn = 44476362.782714Iteration count = 18, obj. fcn = 44476038.325805Iteration count = 19, obj. fcn = 44475944.039470Iteration count = 20, obj. fen = 44475916.655147Iteration count = 21, obj. fcn = 44475908

37、.704029Iteration count = 22, obj. fen = 44475906.395761Iteration count = 23, obj. fcn = 44475905.725710Iteration count = 24, obj. fcn = 44475905.531215Iteration count = 25, obj. fcn = 44475905.474761Iteration count = 26, obj. fcn = 44475905.458374Iteration count = 27, obj. fcn = 44475905.453618Itera

38、tion count = 28, obj. fcn = 44475905.452238Iteration count = 29, obj. fcn = 44475905.451837Iteration count = 30, obj. fcn = 44475905.451721Iteration count = 31, obj. fcn = 44475905.451687Iteration count = 32, obj. fcn = 44475905.451677 index1 =1至19列23456101214151618192021222325262720至21列2829實驗結論和體會：

39、1. 計算各樣品間的歐氏距離、馬氏距離和加權平方距離如上結果所示，篇幅有限， 馬氏距離和加權平方距離得出的結果就不再寫在實驗報告里。2. 根據譜系聚類的有效性，作出檢驗，好確定選擇哪種類間距離最好。得出的R值 可以知道，選擇類間距離為最短距離時效果最好。3. 確定最優的分類數，可以知道最優的分類數是optimaln=12，代入然后進行K均 值聚類法進行聚類，模糊C均值和模糊減法聚類法進行聚類。體會：數據分析工作實在是不容易，要發揚我國的傳統美德，艱苦奮斗才能把統計 搞好。學生簽名： 年 月 日 實驗指導教師評語：教師簽字： 年 月 日注：1.指導教師評語和學生成績一欄由指導教師填寫，其它欄目內容均由學生填寫。2.“實驗項目名稱”要與該實驗課程教學大綱中的“實驗項目”相對應。