Floyd算法_计算最短距离矩阵和路由矩阵_查询最短距离和路由_matlab实验报告

一、实验目的
利用MATLAB实现Floyd算法，可对输入的邻接距离矩阵计算图中任意两点间的最短距离矩阵和路由矩阵，且能查询任意两点间的最短距离和路由。

二、实验原理
Floyd 算法适用于求解网络中的任意两点间的最短路径：通过图的权值矩阵求出任意两点间的最短距离矩阵和路由矩阵。

优点是容易理解，可以算出任意两个
节点之间最短距离的算法，且程序容易实现，缺点是复杂度达到，不适合计算大量数据。

Floyd 算法可描述如下：
给定图G及其边(i , j ) 的权w, j (1 < i < n ,1 <j < n)
FO：初始化距离矩阵W°)和路由矩阵R0)。

其中：
F1：已求得VF和R k-1)，依据下面的迭代求W和R k)
F2:若k<n,重复F1 ;若k>n,终止。

？？
三、实验内容
1、用MATLAB仿真工具实现Floyd算法：给定图G及其边(i , j ) 的权
w, j (1 < i < n ,1 < j < n)，求出其各个端点之间的最小距离以及路由。

(1)尽可能用 M 函数分别实现算法的关键部分，用 M 脚本来进行算法结果验证;
(2)分别用以下两个初始距离矩阵表示的图进行算法验证：
分别求出WT和R7)。

2、根据最短路由矩阵查询任意两点间的最短距离和路由
(1)最短距离可以从最短距离矩阵的3 (i,j)中直接得出；
(2)相应的路由则可以通过在路由矩阵中查找得出。

由于该程序中使用的是前向矩阵，因此在查找的过程中，路由矩阵中r(i,j) 对应的值为Vi到Vj路由上的下一个端点，这样再代入r(r(i,j),j) ，可得到下下个端点，由此不断循环下去，即可找到最终的路由。

(3)对图1，分别以端点对V4 和V6, V3 和V4 为例，求其最短距离和路由；对图2，分别以端点对V1和V7, V3和V5, V1和V6为例，求其最短距离和路由。

3、输入一邻接权值矩阵，求解最短距离和路由矩阵，及某些点间的最短路径。

四、采用的语言
MatLab 源代码：【】 function [w r] = func1(w) n=length(w);
x = w;
r = zeros(n,1);% 路由矩阵的初始化 for i=1:1:n
for j=1:1:n if x(i,j)==inf r(i,j)=0;
else r(i,j)=j;
end, end end; %迭代求出 k 次 w 值 for k=1:n a=w; s = w;
for i=1:n for j=1:n w(i,j)=min(s(i,j),s(i,k)+s(k,j));
end end %根据k-1次值和k次w值求出k次r值 for i=1:n
for
j=1:n if
i==j
r(i,j)=0
;
elseif w(i,j)<a(i,j) r(i,j)=r(i,k);
else r(i,j)=r(i,j);
end end end end;
【】 function [P u]=func2(w,k1,k2) n = length(w);
U = w;
m = 1;
while m <= n
for i = 1:n;
for j = 1:n;
if U(i,j)>U(i,m) + U(m,j) U(i,j) = U(i,m) + U(m,j);
end
end
end
m = m + 1;
end
u = U(k1,k2);
P1=zeros(1,n);
k = 1;
P1(k) = k2;
V = ones(1,n) * 100; kk = k2;
while kk~=k1
for i = 1:n
V(1,i) = U(k1,kk) - w(i,kk); if V(1,i) == U(k1,i) P1(k+1)=i;
kk=i; k=k+1;
end
end
end k=1;
wrow = find(P1~=0); for j=length(wrow):(-1):1
P(k) = P1(wrow(j)); k=k+1;
P;
【】
w1=[0 100 100 100 ;
100 0 100 5 100 2;
100 100 0 100 100 4 ; 5 100 0 100 100;
100 100 0 100;
100 4 100 0 100;
2 100 100 100 0]; w2=[0 2 100 100 100;
0 100 100 100;
2 100 0 100 5 100 ;
100 100 0 100 100 4;
100 5 100 0 100;
100 100 100 0 ;
100 100 4 100 0]; [W1 R1] = func1(w1) [W2 R2] = func1(w2)
【】 w=input(' 输入权值矩阵 w=');
k1=input(' 输入端点 1:k1=');
k2=input(' 输入端点 2:k2=');
w
[W R] = func1(w)
[P u]=func2(w,k1,k2);
disp(['k1 、k2 间最短路： ',num2str(P)]);
disp(['k1 、k2 间最短距离： ',num2str(u)]);
五、数据结构
1. 主要函数最短距离、路由函数：function [w r] = func1(w) n=length(w); x = w;
r = zeros(n,1);% 路由矩阵的初始化
for i=1:1:n
for j=1:1:n if x(i,j)==100 r(i,j)=0;
else
r(i,j)=j;
end,
end
end;
%迭代求岀k次W直
for k=1:n
a=w;
s = w;
for i=1:n
for j=1:n w(i,j)=min(s(i,j),s(i,k)+s(k,j));
end
% 根据k-1次值和k次w值求岀k次r值for i=1:n
for j=1:n
if i==j
r(i,j)=0;
elseif w(i,j)<a(i,j) r(i,j)=r(i,k);
else
r(i,
j)=r
(i,j
);
end end end end; 最短路径函数：function [P u]=func2(w,k1,k2) n = length(w); U = w; m = 1; while m <= n for i = 1:n;
for j = 1:n;
if U(i,j)>U(i,m) + U(m,j) U(i,j) = U(i,m) + U(m,j);
end end end m = m + 1; end u = U(k1,k2); P1=zeros(1,n); k = 1; P1(k) = k2; V = ones(1,n) * 100; kk = k2; while kk~=k1 for i = 1:n
V(1,i) = U(k1,kk) - w(i,kk); if V(1,i) == U(k1,i)
P1(k+1)=i; kk=i; k=k+1;
end end end
k=1;
wrow = fin d(P1 〜=0);
for j=le ngth(wrow):(-1):1 P(k) = P1(wrow(j)); k=k+1;
end
P;
2.算法的流程图
Floyd算法：
六、实验结论与分析通过上图可知，V4和V6之间最短距离是，最短路由是
V4->V1—>V7->V2— >V6, 3和V4之间最短距离是，最短路由是 V3->V7— >V1->V4
通过上图可知，，点对V1和V7之间最短距离是，最短路由是 V1—>V3->V7 端点对
V3和V5之间最短距离是，最短路由是 V — >V1 — >V2 — >V5 端点对V1和V6之
间最短距离是，最短路由是 V1 — >V — >V5 — >V6
七、遇到的问题及解决方法
(1)图的等价表示方法；
(2)两点间的最短路径查询算法。

八、实验心得
通过本次实验实现了用计算机语言编写 Floys 本掌握了算法的实现方法，对MatLab 编程语言更加熟悉，培养了算法设计与优化能力。

此次实验我受益匪浅。