- 在线时间
- 0 小时
- 最后登录
- 2010-9-11
- 注册时间
- 2010-4-19
- 听众数
- 3
- 收听数
- 0
- 能力
- 0 分
- 体力
- 128 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 42
- 相册
- 0
- 日志
- 1
- 记录
- 2
- 帖子
- 6
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 8
升级   38.95% 该用户从未签到
- 自我介绍
- 我是一个数模的喜爱者
 群组: matlab共享与进阶 群组: MATLAB讨论交流群 群组: MATLAB 群组: Matlab讨论组 |
%最短路问题(迪克斯特拉算法). y6 E, ]2 u/ W
clear;
' K6 S1 X8 T( g1 ^# [clc;
/ ?% A2 }- w' r; y' e& ~M=10000;
( ^1 D: b; T5 n& s% O1 Va(1,:)=[0,50,M,40,25,10];% I& ^# Q4 W8 b6 I+ M
a(2,:)=[zeros(1,2),15,20,M,25];
4 H7 P# J$ f- S0 A. [4 n' Ma(3,:)=[zeros(1,3),10,20,M];
1 I0 |( Q+ g O. ]6 c, }# ca(4,:)=[zeros(1,4),10,25];
6 q8 l7 c) t1 R9 b/ _9 [7 Q) {a(5,:)=[zeros(1,5),55];& v) Y6 `8 U/ o- ^& Q$ |
a(6,:)=zeros(1,6);6 G! \0 `. J* e
a=a+a';
5 p P3 _: U; C- H%length函数返回的是矩阵中行数和列数中较大的长度) @! Y. L d b6 G0 ~$ I4 @" b
pb(1:length(a))=0;%生成一个一行多列的0矩阵
, U- e+ U( @7 D8 ipb(1)=1; %第一个数赋值为1( q( d1 O9 Y$ r0 k& V, n
index1=1;
6 |. @9 M0 A: i( |* n5 Lindex2=ones(1,length(a));%生成一个一行多列的1矩阵3 j, ^' p x5 p; D( {0 V X
d(1:length(a))=M; d(1)=0; temp=1;
1 n. q/ U5 u2 i! D8 X%sum函数为各行之和% Y9 W) I* X$ p' e# m( ]' X
while sum(pb)<length(a)
* g" n: ~' b9 Q* i tb=find(pb==0);%find函数返回的是pb矩阵中数值为0的坐标(2,3,4,5,6....)' [0 G* v& x% }7 F0 L2 ^
d(tb)=min(d(tb),d(temp)+a(temp,tb));& L1 _1 f& h) I x$ y$ F
tmpb=find(d(tb)==min(d(tb)));
* t8 m& R# X T1 \ temp=tb(tmpb(1));: G4 ?) e6 B/ y- h3 }) ^
pb(temp)=1;8 l% t5 b. k5 f. k% m+ k$ `: T
index1=[index1,temp];- }1 _4 t7 s# V
index=index1(find(d(index1)==d(temp)-a(temp,index1)));
2 I; \! A8 ^5 w& @ if length(index)>=28 n7 K9 m) A6 X3 ~ Z6 O
index=index(1);
- c; h. l1 }" V& C* d; ? end: e% G+ E+ C t+ `
index2(temp)=index;9 x0 J4 B2 B% g. e
end
2 |' ?" Q: u& a5 L3 I/ {- Hd,index1,index2
5 y- z* L0 v d& c7 x迪杰斯特拉算法用于求解一个有向图(也可以是无向图,无向图是有向图的一种特例)的一个点(称之为原点)到其余各点(称之为周边点)的最短路径问题。算法构思很是巧妙(我这么认为),简直达到了“无心插柳柳成荫”的境界。算法本身并不是按照我们的思维习惯——求解从原点到第一个点的最短路径,再到第二个点的最短路径,直至最后求解完成到第n个点的最短路径,而是求解从原点出发的各有向路径的从小到大的排列(如果这个有向图中有环1-2-3-1算法岂不是永无终结之日了??!!),但是算法最终确实得到了从原点到图中其余各点的最短路径,可以说这是个副产品,对于算法的终结条件也应该以求得了原点到图中其余各点的最短路径为宜。清楚了算法的这种巧妙构思后,理解算法本身就不是难题了。! L$ _% O( n$ C1 y( P3 r
算法把一个图(G)中的点划分成了若干部分:
_7 Z) q/ G( Y; P7 l& U 1):原点(v);
% ~, M1 q* n( {2 L9 x2 S3 N, I 2):所有周边点(C);, X" H" U4 V4 U& C) A: {
另外有一个辅助集合S,从v到S中的点的最短路径已经求得。S的最初状态是空集。
7 H4 ]* X/ K1 r0 a. R$ d- G 这样就可以进一步划分图(G):
& |- ]6 S& X4 F5 t 1):原点(v);# r$ a) ~ T6 i! m1 N9 x) J1 l# U- R
2):已求出v至其最短路径的周边点(S);
1 c8 e& W0 d2 `7 \ 3):尚未求出v至其最短路径的周边点(Other=C-S);" A6 M+ I9 l1 ?6 }) m6 \
算法的主体思想:. ?6 ^* n5 j3 o( W, t
A、找到v——Other所有路径中的的最短路径vd=v——d(Other的一个元素);: |8 ?' y# Y( ~1 L4 X* C
B、找到v——S——Other所有路径中的的最短路径vi=v——i(Other的一个元素);
# p0 o4 ?+ ^8 O3 b: zC、比较vd和vi如果vd<=vi则将d加入S且从Other中删除,否则将i加入S且从Other中删除。6 |# F# t p9 ?: Y, q
重复以上步骤直至Other为空集。
c! L8 n: i9 c+ k( M7 a s/ ^3 u9 `- ^0 E/ Y) }8 G% T n* B9 [
我们求得的最短路径是升序排列的,那为什么下一条最短路径就存在于v—— ( n) q$ ~) e& U+ ^0 {% T7 p
L; F e7 v" V |
zan
|