原创力文档- 1、本文档共4页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
function [Wt,Pp]=mintree(n,W)
%求最小生成树，n为顶点个数，W是权值邻接矩阵，不相邻的用inf表示
%Wt是最小生成树的权，Pp（：，1：2）表示最小生成树的两顶点
%Pp(:，4)表示最小生成树的序号
tmpa=find(W~=inf);
[tmpb,tmpc]=find(W~=inf);
w=W(tmpa);
e=[tmpb,tmpc];
[wa,wb]=sort(w);
E=[e(wb,:),wa,wb];
[nE,mE]=size(E);
temp=find(E(:,1)-E(:,2));
E=E(temp,:);
P=E(1,:);
k=length(E(:,1));
while rank(E)0
? ? temp1=max(E(1,2),E(1,1));
? ? temp2=min(E(1,2),E(1,1));
? ? for i=1:k
? ? ? ? if E(i,1)==temp1
? ? ? ? ? ? E(i,1)=temp2;
? ? ? ? end
? ? ? ? if E(i,2)==temp1
? ? ? ? ? ? E(i,2)=temp2;
? ? ? ? end
? ? end
? ? a=find(E(:,1)-E(:,2));
? ? E=E(a,:);
? ? if rank(E)0
? ? ? ? P=[P;E(1,:)];
? ? ? ? k=length(E(:,1));
? ? end
end
Wt=sum(P(:,3));
Pp=[e(P(:,4),:),P(:,3:4)];
for i=1:length(P(:,3))
? ? disp([,e,num2str(P(i,4)),,...
? ? ? ? (v,num2str(P(i,1)),,v,num2str(P(i,2)),)]);
end
axis equal;%画最小生成树 ??
hold on
[x,y]=cylinder(1,n);
xm=min(x(1,:));
ym=min(y(1,:));
xx=max(x(1,:));
yy=max(y(1,:));
axis([xm-abs(xm)*0.15,xx+abs(xx)*0.15,ym-abs(ym)*0.15,yy+abs(yy)*0.15]);
plot(x(1,:),y(1,:),ko);
for i=1:n
? ? temp=[v,int2str(i)];
? ? text(x(1,i),y(1,i),temp);
end
for i=1:nE
? ? plot(x(1,e(i,:)),y(1,e(i,:)),b);
end
for i=1:length(P(:,4))
? ? plot(x(1,Pp(i,1:2)),y(1,Pp(i,1:2)),r);
end
text(-0.35,-1.2,[最小生成树的权为,,num2str(Wt)]);
title(红色连线为最小生成树);
axis off;
hold off;
调用此程序要先编写主程序。
例：
n=6;
W=inf*ones(6);
W(1,[2,3,4])=[6,1,5];
W(2,[3,5])=[5,3];
W(3,[4,5,6])=[5,6,4];
W(4,6)=2;
W(5,6)=6;
[a,b]=mintree(n,W)
命名为example
在窗口中运行此程序得：
e2(v1v3)
e9(v4v6)
e6(v2v5)
e8(v1v4)
e3(v2v1)
a =
? ? 15
b =
? ? ?1 ? ? 3 ? ? 1 ? ? 2
? ? ?4 ? ? 6 ? ? 2 ? ? 9
? ? ?2 ? ? 5 ? ? 3 ? ? 6
? ? ?3 ? ? 6 ? ? 4 ? ? 8
? ? ?2 ? ? 3 ? ? 5???? 3
文档评论（0）