数学建模--（5）图论模型

最新推荐文章于 2024-01-16 09:39:55 发布

西柚小萌新

最新推荐文章于 2024-01-16 09:39:55 发布

阅读量1.5k

点赞数 1

分类专栏：数学建模文章标签：图论算法数据结构

本文链接：https://blog.csdn.net/qq_58602552/article/details/130160545

版权

数学建模专栏收录该内容

17 篇文章 22 订阅

订阅专栏

本讲将简要介绍图论中的基本概念，并主要讲解图论中的最短路径问题。根据图的不同，我们将学习两种不同的算法：迪杰斯特拉Dijkstra算法和Bellman‐Ford（贝尔曼‐福特）算法。

图的基本概念

图论中的图（Graph）是由若干给定的点及连接两点的线所构成的图形，这种图形通常用来描述某些事物之间的某种特定关系，用点代表事物，用连接两点的线表示相应两个事物间具有这种关系。
一个图可以用数学语言描述为：

G(V(G),E(G))。
%V(vertex)指的是图的顶点集，
%E(edge)指的是图的边集。

根据边是否有方向，可将图分为有向图和无向图。另外，有些图的边上还可能有权值，这样的图称为有权图。

在线作图

CS AcademyCSAcademy is a next generation educational platform. Discover computer science with interactive lessons and a seamless online code editor.https://csacademy.com/app/graph_editor/matlab作图

G1=graph(s,t,w)
plot(G1)
%在s和t中对应的节点之间创建边，并生成一个图,w是权重

digraph(s,t,w)
%有向图

无向图的权重邻接矩阵

有向图的权重邻接矩阵

迪杰斯特拉算法

迪杰斯特拉(Dijkstra)算法是由荷兰计算机科学家狄克斯特拉于1959年提出的。是寻找从一个顶点到其余各顶点的最短路径算法，可用来解决最短路径问题。

具体的过程参照：数据结构–Dijkstra算法最清楚的讲解

更正：第2、3步中，B(23)应为B(13)。

迪杰斯特拉算法采用贪心算法的策略，将所有顶点分为已标记点和未标记点两个集合，从起始点开始，不断在未标记点中寻找距离起始点路径最短的顶点，并将其标记，直到所有顶点都被标记为止。

用法

计算赋权图中各对顶点之间最短路径，显然可以调用Dijkstra算法。

具体方法是：每次以不同的顶点作为起点，用Dijkstra算法求出从该起点到其余顶点的最短路径，反复执行n-1次这样的操作，就可得到从每一个顶点到其他顶点的最短路径。这种算法的时间复杂度 O(n^3)

缺点

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

贝尔曼-福特（Floyd)算法

第二种解决这一问题的方法是由R. W. Floyd提出的算法，称为Floyd算法。其时间复杂度也O(n^3)
) ，但形式上要简单些。

对于赋权图 G = ( V , E , A 0 ) 其中，顶点集 V= {v 1 , . . . , v n }，邻接矩阵：

用法

Floyd算法的基本思想是递推产生一个矩阵序列 A 1 , . . . A k , . . . A n ，其中矩阵A的第 i 行第 j 列元素 A k ( i , j ) 表示从顶点 v i 到vj的路径上所经过的顶点序号不大于 k的最短路径长度。

计算时利用迭代公式：
举一个简单的例子帮助大家理解：

信息初始化：

第一次迭代时，以 a 为起点判断矩阵是否需要更新：

第二次迭代，加入顶点 b，即k=1时判断矩阵是否需要更新：

第三次迭代，加入顶点 c，即k=2时判断矩阵是否需要更新：

matlab计算最短路径

[P,d] = shortestpath(G,start,end [,'Method',algorithm] )
Dists = graphallshortestpaths(DG)  %求所有节点最短路径

找给定范围内所有的点

[nodeIDs,dist] = nearest(G,s,d [,'Method',algorithm])

例题1

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

W = [10,5,1,2,4,7,6,3,9,2];
% sparse生成稀疏矩阵，除了标注的元素外，其余都是0
DG = sparse([1,1,2,2,3,4,4,5,5,5],...  % sparse第一个矩阵，三个.代表没写完，下一行接着
            [2,5,3,5,4,1,3,2,3,4],W);  % sparse第二个矩阵
% sparse里第一个和第二个矩阵相同位置的元素值就是非零元素的值，W是每条边的权值
% dist是最短路径的值，path是最短路径的节点顺序，pred是每一个节点的最短路径的终点前一个节点
[dist,path,pred] = graphshortestpath(DG,1,3);  % 节点1到3的最短路径
point_name = ["city1","city2","city3","city4","city5"];
p = biograph(DG,point_name,'ShowWeights','on');
h = view(p);  % biograph生成图对象，view显示该图
disp(DG);

例题2

用Floyd算法求解和之前同样的问题。

clear;clc;
n=6; a=zeros(n);
a(1,2)=50;a(1,4)=40;a(1,5)=25;a(1,6)=10;
a(2,3)=15;a(2,4)=20;a(2,6)=25; a(3,4)=10;a(3,5)=20;
a(4,5)=10;a(4,6)=25; a(5,6)=55;
a=a+a'; 
a(a==0)=inf;  % 把所有零元素替换成无穷
a([1:n+1:n^2])=0;  % 对角线元素替换成零，Matlab中数据是逐列存储的★
path=zeros(n);
for k=1:n
   for i=1:n
      for j=1:n
         if a(i,j)>a(i,k)+a(k,j)
            a(i,j)=a(i,k)+a(k,j);
            path(i,j)=k;
         end 
      end
   end
end
a, path

例题3

%% 注意：以下代码需要较新版本的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之后才有哦

% 注意：以下代码需要较新版本的matlab才能运行（最好是2016版本及以上哦）
% 如果运行出错请下载新版的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', [] );  

% （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', [] );  

%% 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', [] );  

% 有权图 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', [] );  

% % 注意：代码文件仅供参考，一定不要直接用于自己的数模论文中
% % 国赛对于论文的查重要求非常严格，代码雷同也算作抄袭