matlab与聚类分析

《matlab与聚类分析》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第五章聚类分析安庆师范学院胡云峰习题5.3.下表给出我国历年职工人数（单位：万人），请用有序样品的fisher法聚类。年份全民所有制集体所有制195215802319541881121195624235541958453266219605044925196233031012196434651136196639391264196841701334197047921424197256101524197460071644197668601813197874512048198080192425解第一步数据标准化后计算直径D程序：X=[

2、158023;1881121;2423554;4532662;5044925;33031012;34651136;...39391264;41701334;47921424;56101524;60071644;68601813;...74512048;80192425];stdr=std(X);[n,m]=size(X);X=X./stdr(ones(n,1),:);[np]=size(X);D=zeros(n,n);fori=1:1:n;forj=1:1:n;ifi

3、*sum(X(t,:));y=zeros(1,j-i+1);fors=i:1:jy(s)=(X(s,:)-xgjunzhi)*(X(s,:)-xgjunzhi)';ends=i:1:j;D(i,j)=sum(y);elseD(i,j)=0;endendendD=D'输出结果矩阵太大，所以用excel处理了一下D=0000000000000000.022567000000000000000.448980.2457800000000000002.06321.39810.600240000000000003.92562.6511.1

4、8020.11098000000000004.50223.00911.42380.569530.4086200000000005.1793.43531.66480.825760.538310.020440000000006.08234.0211.9761.0230.633430.127810.047757000000007.03114.65022.32551.23130.7550.263410.112750.01245600000008.33225.57622.90941.60451.05310.606190.338810.1

5、31220.06003200000010.3127.10344.01172.41261.77721.37930.923140.526640.315410.0994010000012.6968.99725.44223.51142.75482.35531.6691.04570.654960.256320.03671000016.29111.9987.86885.50384.56864.11933.10322.14681.47070.771220.308580.1276200021.11716.12811.3218.42987.23

6、166.64875.21163.83122.77931.68770.88810.460160.10709002822.16716.52812.97811.38610.5468.55966.6275.07163.45392.17481.34430.598320.199510我们只看下三角所有元素，其它元素理解为空第二步我们计算损失函数矩阵L程序：%设计一个把样品分为两类的程序,以及对应最后一类分割点D=D';L=zeros(n-1,n-1);alp=zeros(n-1,n-1);form=2:n;s=zeros(1,m-1);f

7、orj=2:ms(1,j-1)=D(1,j-1)+D(j,m);endL(m-1,1)=min(s(1,1:m-1));forj=1:m-1ifL(m-1,1)==s(1,j);alp(m-1,1)=j+1;endendend%分为k类fork=3:n;form=k:ns=zeros(1,m-k+1);forj=k:m;s(1,j-k+1)=L(j-2,k-2)+D(j,m);endL(m-1,k-1)=min(s(1,1:m-k+1));forj=1:m-k+1ifL(m-1,k-1)==s(1,j);alp(m-1,k-1

8、)=j+k-1;endendendend输出结果这里由于表太大，用excel处理一下L=000000000000000.02256700000000000000.448980.0225670000000000000.559960.133550.0225670000000000