最短路问题及其求解方法

最短路问题是最重要的优化问题之一。可以解决很多生产生活中的具体问题，如：管道的铺设、线路的安排、厂区的布置、设备的更新等。类似最短距离（最小费用、最小流量，最短路线）相关的问题都可以抽象为最短路问题。

常见的最短路问题可以分为两类：网络图上的最短路问题，以及网格地图上的最短路问题。其中网络图上的最短路问题主要包括：固定起始点最短路问题（解决该问题的基础算法为Djkstra算法），以及任意两点间的最短路问题（可使用Floyd算法解决）。

一、固定起点最短路（Djkstra算法）

最短路指的是两个节点之间权重最小的路径。 例如下面的网络图，u0到其余各节点的最短路形状为一颗树。这主要是因为最短路是一条路径（路上的节点不重复，如果重复了，可以通过删除相同节点间的路使通路进一步变短），且最短路的任一段也是最短路（如图2：u0到u5的最短路由u0->u1->u4加上edge（u4，u5）。则此时，这段路上u0->u1->u4也是u0到u4的最短路。否则，如果u0->u1->u4的这段路不是u0到u4最短路的话，则可以找到更短的路从u0到u4，再加上edge（u4，u5),构造更短的u0到u5的路，这与u0->u1->u4加上edge（u4，u5）是u0到u5的最短路矛盾。

图1

图2

求固定起点u0到任意点v的最短路的算法（Dijkstra算法）思想

算法的过程就是从始点出发，按照树生长的方式，逐步延长路径，并在每一步都保持这种延长最短。

具体算法步骤（略）

编程实现步骤：

第一步：编写Dijkstra算法的函数

function [min,path]=dijkstra(w,start,terminal)

%输入变量w为所求图的带权邻接矩阵，start、terminal分别为路径的起点和终点的编号。

%返回path为从start到termial的最短路径以及长度min

n=size(w,1); label(start)=0; f(start)=start;

%n为所求图的顶点个数，label存放到各点的最短路径，f(v)表示v的父顶点用来还原路径

%初始化将除了start以外的顶点label均设置为无穷大

for i=1:n

if i~=start

label(i)=inf;

end

end

%s数组存放已经搜好的顶点集，初始化只有start

s(1)=start; u=start;

while length(s)<n

%遍历一遍顶点，将不在顶点集中的顶点选出来进行下面的if判定

for i=1:n

ins=0;

for j=1:length(s)

if i==s(j)

ins=1;

end

end

%判断是否有中继顶点使得它们之间的距离更短，如果有的话更新距离并更新前驱结点

if ins==0

v=i;

if label(v)>(label(u)+w(u,v))

label(v)=(label(u)+w(u,v)); f(v)=u;

end

end

end

v1=0;

k=inf;

%同上再次进行遍历，找到目前最短的路径顶点v1，放入顶点集并改变u的值

for i=1:n

ins=0;

for j=1:length(s)

if i==s(j)

ins=1;

end

end

if ins==0

v=i;

if k>label(v)

k=label(v); v1=v;

end

end

end

s(length(s)+1)=v1;

u=v1;

end

min=label(terminal); path(1)=terminal;

i=1;

%按倒序结果推出最短路径

while path(i)~=start

path(i+1)=f(path(i));

i=i+1 ;

end

path(i)=start;

L=length(path);

%翻转得到最短路径

path=path(L:-1:1);

第二步：建立一个包含上述网络连接关系的txt文件:net1.txt（可以先存在excle文件中，再复制到txt）其中，第一列为连边起始点，第二列为连边起终点，第三列为权重。由于无向图的连边有对称性，因此只存一个方向即可。

122

131

148

246

251

347

379

455

461

472

563

589

674

686

783

之所以建立这样的文件，主要是为了让我们的算法更符合实际问题需求。实际问题相关的网络一般比较

第三步：编写主程序

%生成无向图的邻接矩阵，运行dijikstra函数

clear all

%%%%%%%%%读取存储在net1文件中的网络

dataname = strvcat('net1');

datapath = 'E:\数学建模\数学建模微课陈星\2023年最短路讲座\把实际问题处理成邻接矩阵示例\'; %文件地址

thisdatapath = strcat(datapath,dataname,'.txt'); % 把文件名与文件地址连起来

linklist=load(thisdatapath); % 加载包含网络信息的txt数据�

edge=sortrows(linklist(:,1:3),1); %按第一列从小到达为点排序�

edge=unique(edge,'rows');%删除重复的边�

%%%%%%%%%生成矩阵

Ne =length(edge(:,1)); %网络节点数

vs=max(max(edge)); %�网络边数

a=zeros(vs,vs); %定义初始邻接矩阵 ：元素全部为零

ind=find(a==0);%定义初始邻接矩阵 ：把所有元素变为inf

a(ind)=inf;

for i=1:Ne

i;

a(edge(i,1),edge(i,2))=edge(i,3);

a(edge(i,2),edge(i,1))=edge(i,3);

end

for i=1:vs

a(i,i)=0;

end

%%%%调用dijkstra函数

start=1;terminal=6;

[min,path]=dijkstra(a,start,terminal);

min,path