数学建模：8 图论

在线作图网页：

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

       Dijkstra 局限性是不可以处理带有负权值的图

        Floyd、Bellman-Ford算法可以解决带有负权值的图（有条件）：

• 有负值边的无向图 —— 不可解决（在负值边无限走，则无限小）
• 有环，环中有负值边，环边的权值和为负 —— 不可解决（在环上兜圈子，无限小）

       Dijkstra、Bellman-Ford：单源最短路径

       Bellman-Ford：任意两点间最短路

Matlab作图

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%无向图graph
% 点的编号必须从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', [] ); 

% 有权重
w = [3,8,9,2];
G = graph(s, t, w);
% 在边上标明权重
plot(G, 'EdgeLabel', G.Edges.Weight, 'linewidth', 2) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%有向图digraph
% 有权图 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) 
% G.Edges.Weight 显示边的权重
set( gca, 'XTick', [], 'YTick', [] );

Matlab求最短路径

% Matlab中的图节点要从1开始编号
% 编号从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) % 作图G
set( gca, 'XTick', [], 'YTick', [] ); 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 两点间最短路径
[P,d] = shortestpath(G, 9, 4)  % 9到4的最短路径
% P是最短路径向量 d是最短路径值

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 高亮最短路径
myplot = plot(G, 'EdgeLabel', G.Edges.Weight, 'linewidth', 2);  %首先将图赋给一个变量
highlight(myplot, P, 'EdgeColor', 'r')   %对这个变量进行高亮处理（给边加上r红色）

%% 求出任意两点的最短路径矩阵
D = distances(G) 
D(1,2)  % 1 -> 2的最短路径
D(9,4)  % 9 -> 4的最短路径

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

Dijkstra

三个数据结构：Visited - 记录该点是否被访问；Distance - 到起点的距离；Parent - 上一步的点

两个点集：被选择过的点集O、未被选的点集N

• 每次在 N 中选择离点集 O 距离最小的点加入 O

Bellman-Ford

       贝尔曼‐福特算法不再将节点区分为是否已访问的状态，因为贝尔曼‐福特模型是利用循环来进

行更新权重的，每循环一次，贝尔曼福特算法都会更新所有的节点的信息

原理：

https://blog.csdn.net/a8082649/article/details/81812000

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