暴强Dijkstra算法求任意两点间最短路径.docxVIP

精品文档 精品文档 PAGE 精品文档 效果展示： 开头输入的是点的序列号（表示第几个点），显示的是最短路径的走法（同样以点的序列号显示，表示途径的第几个点）。 %编写m文件 function[distance,path]=dijkstra(A,s,e) [DISTANCE,PATH]=DIJKSTRA(A,S,E) returnsthedistanceandpathbetweenthestartnodeandtheendnode. A:adjcentmatrix s:startnode e:endnode initialize辞谑療縵疡茑济瀉鬮荪祕纰煩緩捫。 n=size(A,1); D=A(s,:); path=[]; visit=ones(1,n); visit(s)=0; parent=zeros(1,n);廈郦藥笔穑談缌銼謎亘级煩箧议钭。  %nodenumber %distancevector %pathvector %nodevisibility %sourcenodeisunvisible %parentnode 钌钡辄癢贵窍蒔齟綜热钆灤蹺铙餉。 %theshortestdistance fori=1:n-1 %BlueSethasn-1nodes temp=zeros(1,n); count=0; forj=1:n ifvisit(j) temp=[temp(1:count)D(j)]; else temp=[temp(1:count)inf]; end count=count+1; end [value,index]=min(temp); j=index;visit(j)=0; fork=1:n ifD(k)D(j)+A(j,k) D(k)=D(j)+A(j,k); parent(k)=j; end瓚阄瀠計呐陧懸纾莢辫释币铈铼緊。 end end distance=D(e); theshortestdistancepathifparent(e)==0 return;end貝貓泞縲屆礦偽阳庫褳鹫苈叁锇录。 path=zeros(1,2*n); t=e;path(1)=t;count=1; whilet~=st0 p=parent(t); path=[ppath(1:count)]; t=p; count=count+1; end ifcount=2*n鵠擋颊貰繆给灭續滎户璽堅鐠绉俩。  %pathpreallocation error([Thepathpreallocationlengthistooshort.,... Pleaseredefinepathpreallocationparameter.]); end path(1)=s; path=path(1:count); 骂覓順愠铵蟄儷憑慍肠芻齲谬邏陣。 %算法实现 clc;clear;closeall; %%载入设置数据 lines=load(); %点与点之间的距离矩阵 A=lines; A(find(A10))=inf; %对步长的限制，根据自己的要求决定！我们在此选择 10. A就是连接矩阵，其中对角线为0，表示本身 有连接关系的就对应线的长度 没有连接关系的就对应inf 铍晉诚崂调辔缛鳐钡誊诹馄韓慳躜。 %%下面的是dijstra 算法，有两种方式可以调用 s=input( 输入起点); %起点（点的序号） e=input( 输入终点); %终点（点的序号） [distance,path0]=dijkstra(A,s,e); fprintf(\nUseDijkstratheMinDistanceis:%.5f\n,distance); fprintf(\nUseDijkstratheMinDistancepathis:\n); disp(path0); A1=A; A1(isinf(A1))=0; [d,p,pred]=graphshortestpath(sparse(A1),s,e);藹鮪繆蕭蒞緹钲網矶闸毀韦赔蓽鴟。 fprintf(\nUsegraphshortestpaththeMinDistanceis:%.5f\n,d); fprintf(\nUsegraphshortestpaththeMinDistancepathis:\n);disp(p); fori=1:length(path0) ifi==length(path0) temp=[path0(1)path0(i)]; else temp=[path0(i)path0(i+1)]; end end慚蝎鹧囪厣魇鍆龆穑锢诳录梦辽雙。