最短距离聚类的matlab实现-1(含聚类图-含距离计算).doc

最短距离聚类的matlab实现-1【2013-5-21更新】说明：正文中命令部分可以直接在Matlab中运行，作者(Yangfd09)于2013-5-21 19:15:50在MATLAB R2009a(7.8.0.347)中运行通过%最短距离聚类（含距离计算，含聚类图）%说明：此程序的优点在于每一步都是自己编写的，很少用matlab现成的指令，%所以更适合于初学者，有助于理解各种标准化方法和距离计算方法。%程序包含了极差标准化（两种方法）、中心化、标准差标准化、总和标准化和极大值标准化等标准化方法，%以及绝对值距离、欧氏距离、明科夫斯基距离和切比雪夫距离等距离计算方法。%=导入数据标准化=m(i) m(i)=testdata(i,j); end if testdata(i,j)选择标准化方式，计算距离如果需要计算明科夫斯基距离，把对应程序行前面的“%”删去即可%if(1)% a=input(计算明可夫斯基距离？（Y/N）,s);% if(a=Y|a=y)% d_ming0=zeros(M,M);% p=input(输入参数p：,s);% p=str2double(p);% end%endfor i=1:M for j=1:M for k=1:N d_abs(i,j)=d_abs(i,j)+abs(test(i,k)-test(j,k); d_ou0(i,j)=d_ou0(i,j)+(test(i,k)-test(j,k)2; d_qie0(k)=abs(test(i,k)-test(j,k);% if(a=Y|a=y)% d_ming0(i,j)=d_ming0(i,j)+abs(test(i,k)-test(j,k)p;% end end d_qie(i,j)=max(d_qie0); endendd_ou=sqrt(d_ou0);%if(a=Y|a=y)d_ming=(d_ming0).(1/p);end%=选择距离，进行最短距离聚类test(i,j)&test(i,j)=0) min=test(i,j); x=i;y=j; end endendt=t+1;str=第,num2str(t),步：lambda=,num2str(min), 合并G(,. num2str(a(x),)和G(,num2str(a(y),),记为G(,num2str(t+MM),);disp(str) %提示操作步骤Z(t,:)=a(x),a(y),min; %收集dendrogram()画聚类图时所需的数据 a(x,y)=; %每执行一步，在a中删除被合并的数据号，在末尾顺次新增一个if(isempty(a) a(end+1)=t+MM;end g=zeros(1,M-2); %生成新的距离矩阵ii=0;for i=1:M if(i=x|i=y); else ii=ii+1; g(ii)=(test(x,i)=test(y,i)*test(y,i); end %.两列变一列：两两比较取其小者endtest(x,y,:)=;test(:,x,y)=;test=cat(1,test,g);test=cat(2,test,g;0); M=length(test(1,:);enddendrogram(Z,0); %画树状聚类图title(最短距离聚类);ylabel(lambda);clear M MM i ii j x y g a t min str Z测试数据：(数据较多，显示可能混乱，复制到.mat文件中即可)1249.9910.17875.4299.92207.78677.08545.31691.23927.091313.12537.44616.05538.41429.95583.13128.99424.2557.63702.97615.36740.2582.47685119.85285.8754.383072.341003.56321.5473.39674.5287.59133.9595.38158.920.60.580.230.660.440.630.490.410.460.740.530.360.250.320.270.490.40.410.370.340.590.310.620.790.410.840.210.340.460.390.410.540.520.570.82184.34150.11291.8723.636.53129.54187.97185.32266.39206.998.92141.47142.8262.88215.2333.3468.83148.63238.2367.74121.13114.64160.0729.9772.054.48416.88193.5636.2179.3473.9925.946.5817.167.851999.972264.55688.58273.7881.65582.67842.64596.63418.615452.911307.271200.081062.29251.41655.54575.11230.59628.591149.2528.76750.61309.822334.812036.83114.9771.75858.55589.43224.79360.57366.59294.0971.1366.12184.72279.09112.8135.2320.3310.5856.7970.9235.7148.14431.8566.4344.9650.1723.3646.7541.8816.7746.0765.8438.7360.4732.37176.15184.7914.9711.5289.8956.1219.7952.4240.8916.954.917.4825.492680.661130.19709.59394.31139.66901.7755.68480.376452597.12568.05742.6524.64162.29503.02210.83264.05412.7497.8513.53574.81346.122040.18951.99219.09162.68909.1747.97178.77412.79586.4264.1685.5181.41236.55577.33225.4375.8965.430.93115.2896.5988.44130.92560.54135.79118.0982.436.9668.0665.7547.9675.6796.1769.68131.4859.7304.76189.0337.1819.81160.68120.0741.9784.2388.5255.0921.9317.8651.76说