求最大可靠路的算法

%% ****************求最大可靠路的算法****************************
% 说明：                                                                                                          //
% 在给定的通信或运输网络中。往往要利用已知的各弧（或边）的可靠性概率，             //
% 求出指定顶点S到T的一条可靠性概率达到最大的路，此路称为由S到T的最大              //
% 可靠路。                                                                                                       //
% 算法思想：（可将概率最大转换为熵最小，即最短路）                                           //
% 已知网络G上各顶点间有向路（或路）的完好概率为0<=pij<=1。设一条由               //
% 顶点S至顶点T的有向路（或路）所经过的顶点序号依次为{S,S1，S2，...,SK,T}       //
% 则这条路所经过的各弧（或边）的完好概率分别是pss1,ps1s2,...,pskT。这条         //
% 路的总完好概率pST是它经过的所有弧（或边）完好概率的乘积。因此，目标可         //
% 以转换为寻找一条使pST取极大值的有向路（或路）的问题。                                 //
% 参数说明：                                                                                                   //
% A 表示网络的完好概率矩阵                                                                              //
% P 表示找到的路径                                                                                          //
% p 表示路径P对应的概率                                                                                  //
% f=0 表示找到路，否则f=1                                                                               //

%%
function [P p f] = p_pathf(A)

[m n] = size(A);
f = 0;
B = zeros(m,n);

%对原矩阵进行转换
for i =1:m
    for j =1:n
        if A(i,j)>0 && A(i,j)<1
            B(i,j) = - log(A(i,j));
        elseif A(i,j) == 0
            B(i,j) = inf;
        end
    end
end
% 利用Warshall-Floyd算法求最短路
[P d] = f_path(B);
if d < inf
    p = 1;
    for i =1:(length(P)-1)
        p = p*A(P(i),P(i+1)); %计算最短路的完好概率
    end
else
    P = 0;p = 0;f=1;
end