图论MATLAB算法汇编.docVIP

第一章：Dijkstra算法 求下面赋权图（左图）中顶点u0到其余顶点的最短路。 其邻接矩阵W为： function dijkstra %注：此程序仅作参考，欢迎批评指正。 clc clear %Dijkstra算法: %%%%给邻接矩阵赋值%%%%%%%%%%%% a [0,1,2,inf,7,inf,4,8; 1,0,2,3,inf,inf,inf,7; 0,0,0,1,5,inf,inf,inf; 0,0,0,0,3,6,inf,inf; 0,0,0,0,0,4,3,inf; 0,0,0,0,0,0,6,4; 0,0,0,0,0,0,0,2;]; for i 2:8 for j 1:i-1 a i,j a j,i ; end end dot size a,1 ;%节点数 fprintf \t邻接矩阵的标准形式： ; a fuquantu a;%在赋权图中用到 fprintf \t其中，inf代表无穷大∞,a i,j 代表第i个节点到第j个节点的权。 1≤i≤%d;1≤j≤%d ,dot,dot ; %%%%%%%%%%%%%%%%%%%%%%%%% U 1 1;%U代表已确定节点的集合 B [1:dot],[1:dot] inf; for i 1:dot [l,n i ] min min a :,U ,[],2 ; L n i l; a :,n i a :,n i +l; B :,U a :,U ; [t,m i ] min min B ; a n i ,: inf; % a if i 2 U i-1 n i ; end end % n,L near [m;n;L] fprintf \t其中，near 3,j 表示第near 2,j 个节点到第1个节点的最短路，\n该最短路径通过边v near 2,j ,near 1,j 。 1≤j≤%d \n,dot ; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %为了直观表示，做出赋权图，为了方便起见节点位置由随机数产生（各个节点的横纵坐标不会相等） a fuquantu; Y 3*randperm dot ;%将1~dot打乱顺序，保证各个节点的横纵坐标不会相等 X 3*randperm dot ; hold on axis equal for i 1:dot-1 for j i+1:dot if a i,j ~ inf % plot [X i ,Y i ],[X j ,Y j ] plot [X i ,X j ],[Y i ,Y j ],linewidth,3 ; hold on text 2*X i +X j /3, 2*Y i +Y j /3,[num2str a i,j ],Rotation,180/pi*atan Y j -Y i / X j -X i ,FontSize,18,VerticalAlignment,cap %将权值标在1/3等分点位置，方向与倾斜角对应，便于区分 end end end hold on for i 2:dot plot [X near 1,i ,X near 2,i ],[Y near 1,i ,Y near 2,i ],linewidth,3,color,r ; % text 2*X near 1,i +X near 2,i /3, 2*Y near 1,i +Y near 2,i /3,[num2str a near 1,i ,near 2,i ],Rotation,180/pi*atan Y near 2,i -Y near 1,i / X near 2,i -X near 1,i ,FontSize,18 end plot X,Y,o,markersize,20,MarkerEdgeColor,g,MarkerFaceColor,g ; for i 1:dot text X i ,Y i ,[n,num2str i-1 ],HorizontalAlignment,center ;%为了美观，保证 x,y 在文本text中心 end 注:标 部分是通过网上查找的资料再借助doc命令写出来的 运行结果： 邻接矩阵的标准形式： a 0 1 2 Inf 7 Inf 4 8 1 0 2 3 Inf Inf Inf 7 2 2 0 1 5 Inf Inf Inf Inf 3 1 0 3 6 Inf Inf 7 Inf 5 3 0 4 3 Inf Inf Inf Inf 6 4 0 6 4 4 Inf Inf Inf 3 6 0 2 8 7