clc,clear;
a=zeros(6);
a(1,2)=50;a(1,4)=40;a(1,5)=25;a(1,6)=10;
a(2,3)=15;a(2,4)=20;a(2,6)=25;
a(3,4)=10;a(3,5)=20;
a(4,5)=10;a(4,6)=25;
a(5,6)=55;
a=a+a';
ll=find(a==0);
a(ll)=inf;
pb(1:length(a))=0;pb(1)=1;index1=1;index2=ones(1,length(a));
d(1:length(a))=inf;d(1)=0;temp=1; %pb是用来标记是否走过（每次取最小路径所到达的地点），index1是指遍历地点的路径，index2是指到底地点前一个地点，然后从index1

                                  %中找到它的最短路径。d是记录从1到每个地点的最短路径。
while sum(pb)<length(a)tb=find(pb==0);                  %取出没走过的点d(tb)=min(d(tb),d(temp)+a(temp,tb));            %用上一次的最小值，和前者与这次的步径相加相比较，得出当前的最短路径tmpb=find(d(tb)==min(d(tb)));             temp=tb(tmpb(1));                %取出目前最短的路径，然后选定为下一个起点pb(temp)=1;index1=[index1,temp];temp2=find(d(index1)==d(temp)-a(temp,index1));      %确定到目前点之前一点，确定的方法是用到此点的最短路径减去这次走的路径与之前的最短路径作比较。（按

                                                        %index1顺序），这样就可以得到上次来到此地的点。index2(temp)=index1(temp2(1));
end
d,index1,index2

clc;clear;
a=zeros(7);a(1,2)=50;a(1,3)=60;
a(2,4)=65;a(2,5)=40;
a(3,4)=52;a(3,7)=45;
a(4,5)=50;a(4,6)=30;a(4,7)=42;
a(5,6)=70;a=a+a';
a(find(a==0))=inf;
result=[];p=1;tb=2:length(a);
while length(result)~=length(a)-1temp=a(p,tb);                  %取出p向量里储存的顶点到其他未选中的点的边temp=temp(:);d=min(temp);                   %从中取出最小的，加入到result，将这个边的终点加入向量p[jb,kb]=find(a(p,tb)==d);j=p(jb(1));k=tb(kb(1));result=[result,[j;k;d]];p=[p,k];tb(find(tb==k))=[];
endresult

clc;clear;a(1,2)=50;a(1,3)=60;
a(2,4)=65;a(2,5)=40;
a(3,4)=52;a(3,7)=45;
a(4,5)=50;a(4,6)=30;a(4,7)=42;
a(5,6)=70;[i,j,b]=find(a);
data=[i';j';b'];                                %data存储每条边的端点和权值，index记录每条边的端点
index=data(1:2,:);
loop=max(size(a))-1;
result=[];while length(result)<looptemp=min(data(3,:));flag=find(data(3,:)==temp);flag=flag(1);v1=data(1,flag);v2=data(2,flag);              %每次找到之后，将数据存入data,而index就开始将线变成点，主要就是用来检验已经取过的树有没有已经走过这两点，

                                                 %从而避免圈的出现。if index(1,flag)~=index(2,flag)result=[result,data(:,flag)];endindex(find(index==v2))=v1;data(:,flag)=[];index(:,flag)=[];
endresult