图的基本概念
利用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,5]; %这里会自动生成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', [] );
邻接矩阵
迪杰斯特拉算法
算法步骤演示
缺点
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之后才有哦
Floyd算法(弗洛伊德算法)
Floyd‐Warshall算法(英语:Floyd‐Warshall algorithm或简写为Floyd algorithm),中文亦称弗洛伊德算法,是解决任意两点间的最短路径的一种算法,可以正确处理无向图或有向图(可以有负权重,但不可存在负权回路)的最短路径问题。
Floyd算法与迪杰斯特拉算法或贝尔曼福特算法相比,能够一次性的求出任意两点之间的最短路径,后两种算法运行一次只能计算出给定的起点和终点之间的最短路径。
当然,Floyd算法计算的时间也要高于后两种算法,其算法核心的步骤由三层循环构成
算法代码:
function [dist,path] = Floyd_algorithm(D) %% 该函数用于求解一个权重邻接矩阵任意两个节点之间的最短路径 % 输入: % D是权重邻接矩阵 % 输出: % dist是最短距离矩阵,其元素dist_ij表示表示i,j两个节点的最短距离 % path是路径矩阵,其元素path_ij表示起点为i,终点为j的两个节点之间的最短路径要经过的节点 n = size(D,1); % 计算节点的个数 % 初始化dist矩阵 dist = D; % 下面我们来初始化path矩阵 path = zeros(n); for j = 1:n path(:,j) = j; % 将第j列的元素变为j end for i = 1:n path(i,i) = -1; % 将主对角线元素变为-1 end % 下面开始三个循环 for k=1:n % 中间节点k从1- n 循环 for i=1:n % 起始节点i从1- n 循环 for j=1:n % 终点节点j从1-n 循环 if dist(i,j)>dist(i,k)+dist(k,j) % 如果i,j两个节点间的最短距离大于i和k的最短距离+k和j的最短距离 dist(i,j)=dist(i,k)+dist(k,j); % 那么我们就令这两个较短的距离之和取代i,j两点之间的最短距离 path(i,j)=path(i,k); % 起点为i,终点为j的两个节点之间的最短路径要经过的节点更新为path(i,k) % 注意,上面一行语句不能写成path(i,j) = k; 这是网上很多地方都容易犯的错误,在PPT11页中会告诉大家为什么不能这么写 end end end end end
核心思想就是,每一次允许从小到大依次从不同的顶点进行中转来更新dist和path.
打印路径
代码:
function [] = print_path(path,dist,i,j) %% 该函数的作用是打印从i到j经过的最短路径 % 输入: % path是使用floyd算法求出来的路径矩阵 % dist是使用floyd算法求出来的最短距离矩阵 % i是起始节点的编号 % j是终点节点的编号 % 输出:无 if i == j warning('起点和终点相同,请检查后重新输入') % 在屏幕中提示警告信息 return; % 不运行下面的语句,直接退出函数 end if path(i,j) == j % 如果path(i,j) = j,则有两种可能: % (1)如果dist(i,j) 为 Inf , 则说明从i到j没有路径可以到达 if dist(i,j) == Inf disp(['从',num2str(i),'到',num2str(j),'没有路径可以到达']) % (2)如果dist(i,j) 不为 Inf , 则说明从i到j可直接到达,且为最短路径 else disp(['从',num2str(i),'到',num2str(j),'的最短路径为']) disp([num2str(i),' ---> ',num2str(j)]) disp(['最短距离为',num2str(dist(i,j))]) end else % 如果path(i,j) ~= j,则说明中间经过了其他节点: k = path(i,j); result = [num2str(i),' ---> ']; % 初始化要打印的这个字符串 while k ~= j % 只要k不等于j, 就一直循环下去 result = [result , num2str(k) , ' ---> ' ]; % i先走到k这个节点处 k = path(k,j); end result = [result , num2str(k)]; disp(['从',num2str(i),'到',num2str(j),'的最短路径为']) disp(result) disp(['最短距离为',num2str(dist(i,j))]) end end
打印所有矩阵函数
function [] = print_all_path(D) %% 该函数的作用是求解一个权重邻接矩阵任意两个节点之间的最短路径,并打印所有的结果出来 % 输入: % D是权重邻接矩阵 % 输出:无 [dist,path] = Floyd_algorithm(D); % 调用之前的Floyd_algorithm函数 n = size(D,1); if n == 1 warning('请输入至少两阶以上的权重邻接矩阵') % 在屏幕中提示警告信息 return; % 不运行下面的语句,直接退出函数 end for i = 1:n for j = 1:n if i ~= j % 不等号用~=表示 print_path(path,dist,i,j); % 调用之前的print_path函数 disp('-------------------------------------------') disp(' ') end end end end
例一:
% PPT第七页的例子 %% 首先将图转换为权重邻接矩阵D n = 5; %一共五个节点 D = ones(n) ./ zeros(n); % 全部元素初始化为Inf for i = 1:n D(i,i) = 0; % 主对角线元素为0 end D(1,2) = 3; D(1,3) = 8; D(1,5) = -4; D(2,5) = 7; D(2,4) = 1; D(3,2) = 4; D(4,3) = -5; D(5,4) = 6; D(4,1) = 2; %% 调用Floyd_algorithm函数求解 [dist,path] = Floyd_algorithm(D) print_path(path,dist,1,5) print_path(path,dist,1,4) print_path(path,dist,3,1) clc disp('下面我们打印任意两点之间的最短距离:') print_all_path(D)
例二:
% 思考题的参考答案 %% 首先将图转换为权重邻接矩阵D n = 9; %一共九个节点 D = zeros(n); % 全部元素初始化为0, 等会你们就知道为什么这样设置啦 % 因为是无向图,所以权重邻接矩阵是一个对称矩阵 D(1,2) = 4; D(1,8) = 8; D(2,8) = 3; D(2,3) = 8; D(8,9) = 1; D(8,7) = 6; D(9,7) = 6; D(9,3) = 2; D(7,6) = 2; D(3,4) = 7; D(3,6) = 4; D(6,4) = 14; D(4,5) = 9; D(6,5) = 10; D = D+D'; % 这个操作可以得到对称矩阵的另一半 for i = 1:n for j = 1:n if (i ~= j) && (D(i,j) == 0) D(i,j) = Inf; % 将非主对角线上的0元素全部变为Inf end end end %% 调用Floyd_algorithm函数求解 [dist,path] = Floyd_algorithm(D)