附录1：
function [T c]=Primf(a) l=length(a); a(a==0)=inf; k=1:l;
listV(k)=0; listV(1)=1; e=1;
while (e if listV(i)==1; for j=1:l if listV(j)==0 && min>a(i,j) min = a(i,j);b=a(i,j); s=i;d=j; end end end end listV(d)=1; distance(e)=b; source(e)=s; destination(e)=d; e=e+1; end T=[source;destination]; for g=1:e-1 c(g)=a(T(1,g),T(2,g)); end c; 15 附录2： sample=[ 0 25 10 100 14.41176471 55.88235294 16.47058824 35.29411765 17.28395062 27.16049383 17.77777778 22.22222222 18.42105263 15.78947368 20 0 21.53846154 71.15384615 22.09302326 13.95348837 23.20610687 21.3740458 24.20382166 28.02547771 26.20689655 41.37931034 28.18181818 54.54545455 29.06779661 87.28813559 30 10 32 45 32.265625 94.140625 32.72727273 84.84848485 34.05063291 93.67088608 35 100 36.02739726 93.15068493 37.77777778 81.48148148 38.0952381 29.76190476 38.91891892 18.91891892 39.28571429 5.357142857 41.42857143 21.42857143 42.28571429 51.42857143 42.66666667 18.33333333 45.39877301 30.67484663 46.08695652 26.08695652 47 7.5 47.25490196 90.19607843 47.36 17.6 48.8 8 50 0 51.42857143 31.42857143 57.24137931 10.34482759 60.84507042 15.49295775 63.22580645 64.51612903 70.4 14 16 72.28070175 70.1754386 73.97260274 34.24657534 76 56 82.22222222 62.22222222 85 50 85.10638298 11.70212766 87.40740741 79.62962963 89.48717949 56.41025641 92.24489796 82.65306122 97.08333333 77.08333333 100.2325581 80.23255814 102.8888889 75.55555556 103.1578947 37.89473684 105.125 78.75 111.3513514 38.91891892 118.8235294 27.45098039 120 100 123.3333333 24.44444444 123.9175258 28.86597938 127.6470588 25.88235294 131.7241379 70.68965517 132.9411765 67.64705882 142.8571429 42.85714286 146 35 147.0588235 32.35294118 160 0 200 50]; 17 附录3： %%%%%Kmeans 聚类函数 function [D,L,cluster_center]=kmeans(k) load dot dist=[];sample=result;dingdian=[]; %%%%%%%%%%%%%%顶点确定 %%画出所有点 figure;hold on for i=1:size(sample,1) plot(sample(i,1),sample(i,2),'k*') end %%%%中心点 x_mean=[mean(result(:,1)),mean(result(:,2))] d=abs2(sample,x_mean); dingdian=find(d==max(d)); n=1; s=zeros(size(d)); while n for j=1:length(dingdian) s=s+abs2(sample,sample(dingdian(j),:)) end [a,b]=sort(s,'descend') for i=1:length(s) if length(find(dingdian==b(i)))==0 n=n+1; dingdian(n)=b(i); break; end end end for i=1:length(dingdian) plot(sample(dingdian(i),1),sample(dingdian(i),2),'ro') end hold off %%%%初始顶点dingdian t=[];q=dingdian;c=[]; for j=1:length(q) c(j,:)=sample(q(j),:); end for i=1:2000 dist=[]; for j=1:length(q) dist=[dist,abs2(sample,c(j,:))]; 18 end [A,B]=min(dist,[],2); cluster_idx=cell(k,1) for j=1:k cluster_idx{j}=find(B==j); end temp=c; for j=1:k c(j,:)= sum(sample(cluster_idx{j},:))/length(cluster_idx{j}); end d1=[]; for j=1:k d1(j)=abs2(c(j,:),temp(j,:)); end if length(find(d1>0.001))==0 cluster_center=c; break; end end cluster=cell(k,1) x=[]; for i=1:k x=sample(cluster_idx{i},:); cluster{i}=x; end hf=figure('NumberTitle','off','name','cluster by Kmeans'); co={'b*','r*','g*','k*','bo','ro','go','ko'}; h=[]; for i=1:k xx=cluster{i}; hold on; h(i) = scatter(xx(:,1),xx(:,2),50,co{i}); scatter(cluster_center(i,1),cluster_center(i,2),'ko'); end switch n case 1 legend(h,'cluster 1') case 2 legend(h,'cluster 1','cluster 2') case 3 legend(h,'cluster 1','cluster 2','cluster 3') case 4 legend(h,'cluster 1','cluster 2','cluster 3','cluster 4') 19