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,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)