最短路径Dijkstra算法原理及Matlab实现「建议收藏」

最短路径Dijkstra算法原理及Matlab实现「建议收藏」图论的基础知识不再阐述。最短路径算法主要有二Dijkstra算法Floyd算法Dijkstra算法研究的是从初始点到其他每一结点的最短路径而Floyd算法研究的是任意两结点之间的最短路径以下图为例,首先介绍Dijstra的原理红字为各结点的编号,蓝字为各结点之间的距离首先定义几个变量结点个数n;二维矩阵M(nxn),距离矩阵,连通的结点间即为距离,不…

大家好,又见面了,我是你们的朋友全栈君。

图论的基础知识不再阐述。
最短路径算法主要有二

  1. Dijkstra算法
  2. Floyd算法

Dijkstra算法研究的是从初始点到其他每一结点的最短路径
而Floyd算法研究的是任意两结点之间的最短路径

以下图为例,首先介绍Dijstra的原理
最短路径Dijkstra算法原理及Matlab实现「建议收藏」
红字为各结点的编号,蓝字为各结点之间的距离

  1. 首先定义几个变量

    结点个数n;
    二维矩阵M(nxn),距离矩阵,连通的结点间即为距离,不连通的结点间为正无穷,和自己的距离为0;
    一维矩阵pb(1xn),若第i点已找到最短路径,则fp(i)=1,否则等于0,对于初始结点,fp=1;
    距离矩阵d(1xn),若第i点已找到最短路径,则的d(i)=最短距离,否则为0,初始结点d=0;
    上一结点矩阵path(1xn),若第i点找到了最短路径,则path存放这一条最短路径的前一个结点,通过对每一点的回溯,可以找到最短路径。

  2. 根据距离写出以下距离矩阵

    这里写图片描述

  3. 确定初始点为v1,则pb(1)=1;

    在图中,结点上,我们将已找到最短路径的点标为它的最短距离,(可以理解为v1点已找到最短路径,距离为0),未找到的其余点表为正无穷(即表示不连通)。
    这里写图片描述
    在与v1连通的点中,即在矩阵m的第1行,寻找最小值,最小值所在列即确定的最短路径的结点,如同v2最短,pb(2)=1,d(2)=1,对于已找到最短路径的v2上一结点为v1,path(2)=1;
    这里写图片描述
    接着,在

    • 与v1连通的,且未找到最短距离的节点的距离
    • 与v2连通的,且未找到最短距离节点的距离+v2的最短距离

    以上两种中寻找最短距离,最短为v6,pb(6)=1;d(6)=2;path(6)=1;
    这里写图片描述
    重复以上步骤,在

    • 与v1连通的,且未找到最短距离的节点的距离
    • 与v2连通的,且未找到最短距离节点的距离+v2的最短距离
    • 与v6连通的,且未找到最短距离节点的距离+v2的最短距离
      以上三种中寻找最短路径,最短为v3,pb(3)=1;d(3)=3);path(3)=6;
      这里写图片描述
      我们可以发现,所要寻找的最短路径即为
      对于已找到最短路径的点(包括初始结点),在与其连通的,未找到最短路径的结点中,将之间距离与圆圈中的距离(即上一结点已找到的最短路径)相加,求得的最小值。
      如果有多个相同的最短距离,任取其中一个。
      最终最短路径即距离如下图
      这里写图片描述

附上代码:

clc;clear;
n=6;   %设置矩阵大小
temp=1;  %设置起始点
m=zeros(6);%定义n阶零矩阵
m(1,2)=1;m(1,6)=2;%设置矩阵中非零非无穷的值
m(2,1)=1;m(2,3)=4;m(2,6)=4;
m(3,2)=4;m(3,4)=2;m(3,6)=1;
m(4,3)=2;m(4,5)=3;m(4,6)=3;
m(5,4)=3,m(5,6)=5;
m(6,1)=2;m(6,2)=4;m(6,3)=1;m(6,4)=3;m(6,5)=5;

for i=1:n
    for j=1:n
       if(m(i,j)==0)
           m(i,j)=inf;
       end
    end
end
for i=1:n
    m(i,i)=0;
end
pb(1:length(m))=0;pb(temp)=1;%求出最短路径的点为1,未求出的为0
d(1:length(m))=0;%存放各点的最短距离
path(1:length(m))=0;%存放各点最短路径的上一点标号
while sum(pb)<n %判断每一点是否都已找到最短路径
 tb=find(pb==0);%找到还未找到最短路径的点
 fb=find(pb);%找出已找到最短路径的点
 min=inf;
 for i=1:length(fb)
     for j=1:length(tb)
         plus=d(fb(i))+m(fb(i),tb(j));  %比较已确定的点与其相邻未确定点的距离
         if((d(fb(i))+m(fb(i),tb(j)))<min)
             min=d(fb(i))+m(fb(i),tb(j));
             lastpoint=fb(i);
             newpoint=tb(j);
         end
     end
 end
 d(newpoint)=min;
 pb(newpoint)=1;
 path(newpoint)=lastpoint; %最小值时的与之连接点
end
d
path

路径只需向上回溯就可得到。

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/137532.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • SQL Server中QUOTENAME函数的使用

    SQL Server中QUOTENAME函数的使用–函数QUOTENAME–功能:返回带有分隔符的Unicode 字符串,分隔符的加入可使输入的字符串成为有效的Microsoft SQL Server 2005 分隔标识符。–语法QUOTENAME ( character_string [ , quote_character ] ) –举例说明:–比如你有一个表,名字叫index–你有一个动

    2022年7月25日
    7
  • laravel 授权使用gate门类

    laravel 授权使用gate门类

    2021年10月24日
    36
  • redis集群实战

    redis集群实战

    2021年9月15日
    45
  • error: failed to push some refs to ‘git@gitee.com:xxx/xxx.git’ 问题 解决办法

    error: failed to push some refs to ‘git@gitee.com:xxx/xxx.git’ 问题 解决办法原因分析:本地和远程的文件应该合并后才能上传本地的新文件解决办法:1.gitpullgiteemaster先拉下来,自动合并(类似于SVN的update)2.gitpushgit

    2022年7月2日
    17
  • Stata Kendall 相关系数作图

    Stata Kendall 相关系数作图StataKendall相关系数作图回答Superficial.的问题,测试CSDN的markdown发帖功能如何插入一段漂亮的代码片生成一个适合你的列表创建一个表格设定内容居中、居左、居右SmartyPants创建一个自定义列表如何创建一个注脚注释也是必不可少的KaTeX数学公式新的甘特图功能,丰富你的文章UML图表FLowchart流程图导出与导入导出导入此帖目的有二:回答Superficial.的问题,测试CSDN的markdown发帖功能如何插入一段漂亮的代码片去博客设置页面,选择

    2022年6月17日
    24
  • 【stm32f407】SPI实验 驱动W25Q128「建议收藏」

    一.SPI介绍SPI是英语SerialPeripheralinterface的缩写,顾名思义就是串行外围设备接口。是Motorola首先在其MC68HCXX系列处理器上定义的。SPI接口主要应用在EEPROM,FLASH,实时时钟,AD转换器,还有数字信号处理器和数字信号解码器之间。SPI,是一种高速的,全双工,同步的通信总线,并且在芯片的管脚上只占用四根线,节约了芯片的管脚…

    2022年4月7日
    40

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号