%% ****************求最大可靠路的算法**************************** % 说明: // % 在给定的通信或运输网络中。往往要利用已知的各弧(或边)的可靠性概率, // % 求出指定顶点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