基于MATLAB的Kruskal避圈算法求最小生成树.docVIP

基于MATLAB的Kruskal避圈算法求最小生成树 function E=Kruskal(w)%图论最小生成树Kruskal避圈算法(使用时根据题目修改w和n)%w为邻接矩阵[m,n]=size(w);k=1;for i=1:n-1for j=i+1:n??? if w(i,j)~=0?????? x(1,k)=w(i,j);%记录边?????? x(2,k)=i;%记录起点?????? x(3,k)=j;%记录终点?????? k=k+1;??? endendendk=k-1;%统计边数 k为边数 %步骤一%冒泡法给边的大小排序for i=1:kfor j=i+1:k??? if x(1,i)x(1,j)??? a=x(1,i);x(1,i)=x(1,j);x(1,j)=a;??? a=x(2,i);x(2,i)=x(2,j);x(2,j)=a;??? a=x(3,i);x(3,i)=x(3,j);x(3,j)=a;??? endendend%给各点标号赋初值for i=1:nl(i)=i;end%初始时选e1加入集合EE(1,1)=x(1,1); %E矩阵的第一行记录最小生成树的边长E(2,1)=x(2,1); %E矩阵的第二行记录边的起点E(3,1)=x(3,1); %E矩阵的第三行记录边的终点a=min([l(E(2,1)),l(E(3,1))]);l(E(2,1))=a;l(E(3,1))=a;b=1;%记录E中边数 for i=2:k%步骤四?? if b==n-1 %如果树中边数达到n-1????? break???? %算法终止?? end%步骤二if l(x(2,i))~=l(x(3,i)) %如果两顶点标号不同??? b=b+1;?????????????????? %将这条边加入E??? E(1,b)=x(1,i);??? E(2,b)=x(2,i);??? E(3,b)=x(3,i); %步骤三for j=1:n??????????????? %对于所有顶点??? if l(j)==max([l(E(2,b)),l(E(3,b))])%如果该顶点的标号,等于=,新加入边中的顶点标号较大的值??????? l(j)=min([l(E(2,b)),l(E(3,b))]);%将其改为较小的那一个以避圈??? endendendendE