图论模型——图论最短路径求解:迪杰斯特拉算法和贝尔曼‐福特算法。

个人学习笔记,课程为数学建模清风付费课程

目录

一、图论基础

1.1图的基本概念

1.2在线作图

1.2.1网站

1.2.2MATLAB

1.3无向图的权重邻接矩阵 

1.4有向图的权重邻接矩阵 

二、迪杰斯特拉算法 

2.1概念

2.2步骤 

2.3缺点

三、贝尔曼‐福特算法

3.1概念

3.2负权回路

3.3代码 

3.3.1计算最短路径

3.3.2返回任意两点的距离矩阵

3.3.3找给定范围内所有的点

3.3.4示例代码 


一、图论基础

1.1图的基本概念

       图论中的图( Graph )是由若干给定的点及连接两点的线所构成的图形,这种图形通常用来描述某些事物之间的某种特定关系,用点代表事物,用连接两点的线表示相应两个事
物间具有这种关系。
       一个图可以用数学语言描述为 G(V(G),E(G)) V(vertex) 指的是图的顶点集,E(edge) 指的是图的边集。根据边是否有方向,可将图分为有向图和无向图。
       另外,有些图的边上还可能有权值,这样的图称为有权图。

 1.2在线作图

1.2.1网站

https://csacademy.com/app/graph_editor/

1.2.2MATLAB

注:
1 Matlab 做出来的图不是很漂亮,要是节点比较少,还是推荐大家使用在线作图。
2 )该函数在 2015b 之后的版本才支持,如果运行出错请下载新版本 Matlab
%% Matlab作无向图
% (1)无权重(每条边的权重默认为1)
% 函数graph(s,t):可在 s 和 t 中的对应节点之间创建边,并生成一个图
% s 和 t 都必须具有相同的元素数;这些节点必须都是从1开始的正整数,或都是字符串元胞数组。
s1 = [1,2,3,4];
t1 = [2,3,1,1];
G1 = graph(s1, t1);
plot(G1)
% 注意哦,编号最好是从1开始连续编号,不要自己随便定义编号
s1 = [1,2,3,4];
t1 = [2,3,1,1];
G1 = graph(s1, t1);
plot(G1)

% 注意字符串元胞数组是用大括号包起来的哦
s2 = {'学校','电影院','网吧','酒店'};
t2 = {'电影院','酒店','酒店','KTV'};
G2 = graph(s2, t2);
plot(G2, 'linewidth', 2)  % 设置线的宽度
% 下面的命令是在画图后不显示坐标
set( gca, 'XTick', [], 'YTick', [] );  

%% Matlab作有向图
% 无权图 digraph(s,t)
s = [1,2,3,4,1];
t = [2,3,1,1,4];
G = digraph(s, t);
plot(G)
set( gca, 'XTick', [], 'YTick', [] );  

1.3无向图的权重邻接矩阵 

结论:
1 )无向图对应的权重邻接矩阵 D 是一个对称矩阵;
2 )其主对角线上元素为 0.
3 𝐷 ij 表示第 i 个节点到第 j 个节点的权重。
%% Matlab作无向图
% (2)有权重
% 函数graph(s,t,w):可在 s 和 t 中的对应节点之间以w的权重创建边,并生成一个图
s = [1,2,3,4];
t = [2,3,1,1];
w = [3,8,9,2];
G = graph(s, t, w);
plot(G, 'EdgeLabel', G.Edges.Weight, 'linewidth', 2) 
set( gca, 'XTick', [], 'YTick', [] );  

1.4有向图的权重邻接矩阵 

结论:
1 )有向图对应的权重邻接矩阵 D 是一般不再是对称矩阵;
2 )其主对角线上元素为 0.
3 𝐷 ij 表示第 i 个节点到第 j 个节点的权重。
%% Matlab作有向图
% 有权图 digraph(s,t,w)
s = [1,2,3,4];
t = [2,3,1,1];
w = [3,8,9,2];
G = digraph(s, t, w);
plot(G, 'EdgeLabel', G.Edges.Weight, 'linewidth', 2) 
set( gca, 'XTick', [], 'YTick', [] );  

二、迪杰斯特拉算法 

2.1概念

https://www.bilibili.com/video/av54668527

2.2步骤 

2.3缺点

可以用于有向图  但不能处理负权重

在本例中,按照该算法,从1到2的最短路径为2,路径为1->2,但是实际上最短路径为1->3->2,距离为1

 

若要修复该缺点,则需要使用Bellman‐Ford(贝尔曼福特)算法

三、贝尔曼福特算法

3.1概念

用于求解单源、有负权边的最短路问题

实现通过m次迭代求出从起点到终点不超过m条边构成的最短路径

有兴趣的同学可以参考下面两份资料弄懂其实现原理:

3.2负权回路

       在一个图里每条边都有一个权值(有正有负 ) 如果存在一个环(从某个点出发又回到自己的路径),而且这个环上所有权值之和是负数,那这就是一个负权环,也叫负权回路。
       存在负权回路的图是不能求两点间最短路的,因为只要在负权回路上不断兜圈子,所得的最短路长度可以任意小。
含有负权重的无向图都是负权回路。

庆幸的是,含有负权重的图特别少见,且一旦出现负权重,也往往是在有向图中。因此大家不用担心算法求解不出来的问题。

3.3代码 

3.3.1计算最短路径

可选的算法:

3.3.2返回任意两点的距离矩阵

3.3.3找给定范围内所有的点

3.3.4示例代码 

%% 注意:以下代码需要较新版本的matlab才能运行(最好是2016版本及以上哦)
% 如果运行出错请下载新版的matlab代码再运行

% 注意哦,Matlab中的图节点要从1开始编号,所以这里把0全部改为了9
% 编号最好是从1开始连续编号,不要自己随便定义编号
s = [9 9 1 1 2 2 2 7 7 6 6  5  5 4];
t = [1 7 7 2 8 3 5 8 6 8 5  3  4 3];
w = [4 8 3 8 2 7 4 1 6 6 2 14 10 9];
G = graph(s,t,w);
plot(G, 'EdgeLabel', G.Edges.Weight, 'linewidth', 2) 
set( gca, 'XTick', [], 'YTick', [] );  
[P,d] = shortestpath(G, 9, 4)  %注意:该函数matlab2015b之后才有哦

% 在图中高亮我们的最短路径
myplot = plot(G, 'EdgeLabel', G.Edges.Weight, 'linewidth', 2);  %首先将图赋给一个变量
highlight(myplot, P, 'EdgeColor', 'r')   %对这个变量即我们刚刚绘制的图形进行高亮处理(给边加上r红色)

% 求出任意两点的最短路径矩阵
D = distances(G)   %注意:该函数matlab2015b之后才有哦
D(1,2)  % 1 -> 2的最短路径
D(9,4)  % 9 -> 4的最短路径

% 找出给定范围内的所有点  nearest(G,s,d)
% 返回图形 G 中与节点 s 的距离在 d 之内的所有节点
[nodeIDs,dist] = nearest(G, 2, 10)   %注意:该函数matlab2016a之后才有哦

  • 19
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值