数学建模社区-数学中国

标题: 图中最短路径解决方法 Dijkstra [打印本页]

作者: 2744557306    时间: 2023-12-22 11:02
标题: 图中最短路径解决方法 Dijkstra
这段代码是Dijkstra算法的一个实现,用于在带权图中找到从一个起始点到其他所有点的最短路径。下面是对代码的逐步解释:
1 M" `7 p  ]* Dclear all8 o/ M7 J5 Y/ P/ U. h) [
%图论最短路问题的Dijkstra算法
  j( d$ c, u" q( i, ~) G%邻接矩阵(点与点的关系)' d! M) e+ R  q
w=[0,2,4,inf,inf,inf,inf;  3 X; k' q' G6 F/ {0 K
   2,0,inf,3,3,1,inf;
6 X) R8 `- G/ h3 m6 B7 \7 u/ z   4,inf,0,2,3,1,inf;                     
4 }; X9 `1 |8 X- ?   inf,3,2,0,inf,inf,1;                    6 w4 r- ?7 m( h) O) z
   inf,3,3,inf,0,inf,3;4 F) D3 t1 E' P  r1 U
   inf,1,1,inf,inf,0,4;1 W( U1 `* t  ^( {
   inf,inf,inf,1,3,4,0];
* E$ ]) e1 c) h" K' U. y- Un=size(w,1);%记录图中点数2 K* y) u( v/ z. q: K8 G/ v6 R

, G5 \8 k* a9 _8 H" yfor i=1:n# O8 ^% m7 w* n3 ^* `
    l(i)=w(1,i); %为l(v)赋初值
9 P5 _* Y0 J% y& M    z(i)=1;      %为z(v)赋初值1
& U9 P# H) k6 n) T. D6 R2 Gend. p. @* s! U9 W, J. X
4 w; ?& o4 V- F' l/ B
s=[];            %s集合
" h2 H& G2 }1 V$ ^s(1)=1;          %s集合的第1个元素为起点
1 J+ f  b- U& y6 D; au=s(1);0 s8 v2 F$ S) O7 s% j3 h4 f; N% S4 s3 y
k=1;             %k记录集合s中点的数量& K4 C% J- g  |
& W0 {+ z; S0 I4 Z
while k<n       %当集合s未包含所有元素的时候执行循环
* B& _) |' T" i0 o    for i=1:n    %更新一遍l(v),z(v)9 q& y9 P' l  y2 D8 w7 n* X% K
        if l(i)>l(u)+w(u,i), o" F# [8 w4 H' m
            l(i)=l(u)+w(u,i);
8 {+ l( T3 p$ N            z(i)=u;4 v0 U# n7 ?+ n. P. l
        end( u) Z8 `7 [9 q, p
    end" \0 R1 s5 o( q1 i- V
( a" P* Q6 n! A+ R; b
    %找l(i)中最小的v加入s集合
$ D" n, b0 @# ^8 a  g+ ^% i    ll=l;
( G  I( K# l. \$ ]4 A& o6 u    for i=1:n) J. S/ E7 C3 L/ W
        for j=1:k
/ @% U/ V1 s0 E* Q            if i==s(j), O" v6 d9 Z$ r2 O: \
                ll(i)=inf; %去除掉已经在s集合中的点* N( ]& x( R, q. Q, j
            end 3 A, J* s( S9 r* f  h/ _
        end' t% t1 I  S4 Q2 q7 x8 _
    end; B, Y  _$ i) w5 _
    [lv,v]=min(ll); %求最小的l(v)
" J" J& y3 w6 z/ K5 _# @6 N    s(k+1)=v;       %加入集合s
) U& x5 `) M* [& S    u=v;" f0 S0 z% e. F) ?
    k=k+1;
( c& s" p# F# B$ {+ uend
) J1 n/ `0 j& b6 P
  C# h" M! G. d# R# m" C# S+ S7 R& ifprintf('最短路为:%1d->%1d->%1d->%1d\n',z(z(z(7))),z(z(7)),z(7),7)  C$ j; O3 Y2 I* R

8 D0 \. N% B' [解释:
0 b' m6 w/ E9 q& i8 Z, L+ }+ g. B
  |$ r2 o. H2 h6 D3 B1.清除变量: clear all 语句清除所有之前定义的变量,以确保从干净的状态开始执行程序。
+ e2 J+ p( K4 f, y2.邻接矩阵: 图被表示为邻接矩阵 w,其中 w(i, j) 表示从顶点 i 到顶点 j 的边的权重。 inf 表示没有直接的边。
, y* Z+ k# x7 X! q3.初始化: l 数组存储从起始点到每个点的当前最短路径长度,z 数组存储路径中的前一个顶点。初始时,将起始点到每个点的距离初始化,并将前一个顶点初始化为起始点。) b5 i9 C1 m  V/ g4 [0 |
4.主循环: 使用 while 循环,每次选择一个新的顶点加入集合 s,直到 s 包含所有顶点。在每次循环中,通过更新 l 和 z 数组来逐步找到最短路径。
4 E# e9 q4 W6 X1 A$ d5.输出: 最后,使用 fprintf 打印从起始点到指定点的最短路径。在这里,路径是通过回溯 z 数组得到的。1 f' Y$ M1 z$ S' q9 @

# N9 s3 t+ t  f$ I1 a! l注意:这段代码的输出是针对特定的终点(顶点7)进行的,你可能需要根据你的需求更改这个值。
/ \! s+ }7 u$ j9 S% c: v( _' l! s4 ^$ c, C1 F

$ w) o! o7 e0 l1 m8 @( F4 L0 y# B

Dijkstra.m

913 Bytes, 下载次数: 0, 下载积分: 体力 -2 点

售价: 1 点体力  [记录]  [购买]






欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5