Dijstra算法求一个点到另一个点之间的最短路径和距离

最新推荐文章于 2020-12-11 15:00:00 发布

为啥不能修改昵称啊

最新推荐文章于 2020-12-11 15:00:00 发布

阅读量2.6k

点赞数

分类专栏：数学建模

本文链接：https://blog.csdn.net/weixin_43845922/article/details/108174312

版权

数学建模专栏收录该内容

11 篇文章 1 订阅

订阅专栏

Dijstra

输入邻接矩阵
- 结果截图
调用的函数

输入邻接矩阵

% weight=    [0     2     8     1   Inf   Inf   Inf   Inf   Inf   Inf   Inf;
%             2     0     6   Inf     1   Inf   Inf   Inf   Inf   Inf   Inf;
%             8     6     0     7     5     1     2   Inf   Inf   Inf   Inf;
%             1   Inf     7     0   Inf   Inf     9   Inf   Inf   Inf   Inf;
%           Inf     1     5   Inf     0     3   Inf     2     9   Inf   Inf;
%           Inf   Inf     1   Inf     3     0     4   Inf     6   Inf   Inf;
%           Inf   Inf     2     9   Inf     4     0   Inf     3     1   Inf;
%           Inf   Inf   Inf   Inf     2   Inf   Inf     0     7   Inf     9;
%           Inf   Inf   Inf   Inf     9     6     3     7     0     1     2;
%           Inf   Inf   Inf   Inf   Inf   Inf     1   Inf     1     0     4;
%           Inf   Inf   Inf   Inf   Inf   Inf   Inf     9     2     4     0;];
weight = [0 50 Inf 40 25 10;
    50 0 15 20 Inf 25;
    Inf 15 0 10 20 Inf;
    40 20 10 0 10 25;
    25 Inf 20 10 0 55;
    10 25 Inf 25 55 0];
for i = 2:6
    fprintf('c1到c%d的最近距离和路径如下:',i);
      [dis, path]=dijkstra(weight,1,i)
      fprintf('\n');
end
% weight为邻接矩阵
% 1 为起始点， i为终止点

结果截图

在这里插入图片描述

调用的函数

function [min,path]=dijkstra(w,start,terminal)
n=size(w,1); label(start)=0; f(start)=start;
for i=1:n
   if i~=start
       label(i)=inf;
end, end
s(1)=start; u=start;
while length(s)<n
   for i=1:n
      ins=0;
      for j=1:length(s)
         if i==s(j)
            ins=1;
         end,  
      end
      if ins==0
         v=i;
         if label(v)>(label(u)+w(u,v))
            label(v)=(label(u)+w(u,v)); 
         f(v)=u;
         end, 
      end, 
   end   
v1=0;
   k=inf;
   for i=1:n
         ins=0;
         for j=1:length(s)
            if i==s(j)
               ins=1;
            end, 
         end
         if ins==0
            v=i;
            if k>label(v)
               k=label(v);  v1=v;
            end,  
         end,  
   end
   s(length(s)+1)=v1;  
   u=v1;
end
min=label(terminal); path(1)=terminal;
i=1; 
while path(i)~=start
      path(i+1)=f(path(i));
      i=i+1 ;
end
path(i)=start;
L=length(path);
path=path(L:-1:1);