数模学习日记——最短路径

1.例题

在下图中，用点表示城市，现有共五个城市。点与点之间的连线表示城市间有道路相连。连线旁的数字表示道路的长度。先计划从城市到城市铺设一条天然气管道，请设计出最小长度管道铺设方案。

---

2.定义

从图中的某个顶点出发，到达另一个顶点的所经过的边的权重之和最小的一条路径

• 图：边和节点组成的结构，在数学建模中例如本题目中道路和城市
• 边带有方向的是有向图，否则为无向图
• 权重：每条边都有之对应的值，本题中边道路，边的权重就是道路长度，当然是越小越好

---

3.解决方法

求解最短路径：Dijkstra算法，或MATLAB的graphshortestpath函数（目前matlab有shortestpath与shortestpathtree函数）

---

4.语法

①graphshortestpath函数

[dist,path,pred] = graphshortestpath(G,S)
[ ___ ] = graphshortestpath(G,S,D)
[ ___ ] = graphshortestpath(___,Name,Value)

②shortestpath函数

P = shortestpath(G,s,t)
P = shortestpath(G,s,t,'Method',algorithm)
[P,d] = shortestpath(___)
[P,d,edgepath] = shortestpath(___)

•  P = shortestpath(G,s,t) 计算从源节点 s 处开始到目标节点 t 处结束的最短路径。如果图进行了加权（即 G.Edges 包含变量 Weight），则这些权重用作沿图中各边的距离。否则，所有边距离都视为 1。
• P = shortestpath(G,s,t,'Method',algorithm) 可选择性地指定在计算最短路径时使用的算法。例如，如果 G 是加权图，则 shortestpath(G,s,t,'Method','unweighted') 将忽略 G 中的边权重，而将所有边权重视为 1。
• [P,d] = shortestpath(___) 还使用上述语法中的任何输入参数返回最短路径的长度 d。
• [P,d,edgepath] = shortestpath(___) 还返回从 s 到 t 的最短路径上所有边的边索引 edgepath。

③shortestpathtree函数

TR = shortestpathtree(G,s)
TR = shortestpathtree(G,s,t)
TR = shortestpathtree(___,Name,Value)
[TR,D] = shortestpathtree(___)
[TR,D,E] = shortestpathtree(___)

•  TR = shortestpathtree(G,s) 返回有向图 TR，其中包含从源节点 s 到图中所有其他节点的最短路径树。如果图进行了加权（即 G.Edges 包含变量 Weight），则这些权重用作沿图中各边的距离。否则，所有边距离都视为 1。

• TR = shortestpathtree(G,s,t) 计算多个源或目标节点之间的最短路径树：s可以是单个源节点，t 可以指定多个目标节点。 s可以指定多个源节点，t 可以指定单个目标节点。

• TR = shortestpathtree(___,Name,Value) 使用由一个或多个名称-值对组参数指定的其他选项，这些选项使用上述语法中的任意输入参数组合。例如，shortestpathtree(G,s,'OutputForm','vector') 返回用于描述最短路径树的数值向量

• [TR,D] = shortestpathtree(___) 还返回树中各节点之间的最短路径距离

• [TR,D,E] = shortestpathtree(___) 还返回逻辑向量 E，指示图中的每条边是否在 TR 中。

---

5.题目解决

使用graphshortestpath函数

%sparse生成稀疏矩阵
W = [10,5,2,1,4,6,7,3,9,2]
DG = sparse([1,1,2,2,3,4,4,5,5,5], ...
    [2,5,5,3,4,3,1,2,3,4],W)
%biograph生成图对象；view显示该图
point_name = ["城市1","城市2","城市3","城市4","城市5"]
h = view(biograph(DG,point_name,'ShowWeights','on'))
%dist是最短路径的值
%path是最短路径的节点顺序
%pred是到每一个节点的最短路径的终点前一个节点
[dist,path,pred] = graphshortestpath(DG,1,3)
%将最短路径的节点和边缘标记为红色并增加线宽
%getedgesbynodeid得到图h的指定边的句柄
%第一个参数是图，第二个是边的出发点，第三个是边的入点
%句柄确保能找到对应的东西
%get查询图的属性，h.Nodes(path),'ID'得到图h中最短路径的边
%set函数设置图形属性
edges = getedgesbynodeid(h,get(h.Nodes(path),'ID'));
set(edges,'LineColor',[1,0,0])  %RGB数值，红绿蓝
set(edges,'LineWidth',2)