数学建模学习笔记--图论中最短路算法与程序实践

1.Floyd算法

1）根据已知的部分节点之间的连接信息，建立初始距离矩阵B(i,j),其中没有给出距离的赋予一个充分大数值，以便于更新。（i,j=1,2,…n）
2）进行迭代计算。对任意点（i，j），若存在k，使B（i，k)+B(k,j)<B(i,j),则更新B（i，j）=B(i,k)+B(k,j)
3)直到所有点距离不再更新停止计算，得到最短路径矩阵B（i，j）

for k=1:n
for i=1 :n
for j=1:n
t=B(i,k)+B(k,j);
if t<B(i,j) B(i,j)=t; end
end
end
end

实例： 已知50个点之间相互连接信息见表1及续表，求最短距离矩阵

n=50; %Matlab实现的Floyd算法
A=zeros(n,n);
for i=1:n
for j=1:n
if(i==j) A(i,j)=0;
else A(i,j)=100000;
end
end
end %赋直接距离信息
A(1,2)=400;A(1,3)=450; A(2,4)=300;A(2,21)=230; A(2,47)=140;A(3,4)=600;
A(4,5)=210;A(4,19)=310;A(5,6)=230;A(5,7)=200; A(6,7)=320; A(6,8)=340;
A(7,8)=170;A(7,18)=160;A(8,9)=200;A(8,15)=285; A(9,10)=180; A(10,11)=150;
A(10,15)=160; A(11,12)=140; A(11,14)=130; A(12,13)=200; A(13,34)=400;
A(14,15)=190;A(14,26)=190; A(15,16)=170; A(15,17)=250; A(16,17)=140;
A(16,18)=130; A(17,27)=240; A(18,19)=204; A(18,25)=180; A(19,20)=140;
A(19,24)=175; A(20,21)=180; A(20,24)=190; A(21,22)=300; A(21,23)=270;
A(21,47)=350;A(22,44)=160;A(22,45)=270;A(22,48)=180;A(23,24)=240;
A(23,29)=210;A(23,30)=290;A(23,44)=150;A(24,25)=170;A(24,28)=130;
A(26,27)=140;A(26,34)=320;A(27,28)=190;A(28,29)=260;A(29,31)=190;
A(30,31)=240;A(30,42)=130;A(30,43)=210;A(31,32)=230;A(31,36)=260;
A(31,50)=210;A(32,33)=190;A(32,35)=140;A(32,36)=240;A(33,34)=210;
A(35,37)=160;A(36,39)=180;A(36,40)=190;A(37,38)=135;A(38,39)=130;
A(39,41)=310;A(40,41)=140;A(40,50)=190;A(42,50)=200;A(43,44)=260;
A(43,45)=210;A(45,46)=240;A(46,48)=280;A(48,49)=200;
for j=1:n
for i=1:j-1
A(j,i)=A(i,j); %使矩阵对称
end
end
B=A;
%利用Floyd算法计算最短距离矩阵
for k=1:n
for i=1 :n
for j=1:n
t=B(i,k)+B(k,j);
if t<B(i,j) B(i,j)=t; end
end
end
end

本文来源：中国大学mooc西北工业大学数学建模课程