计算最小费用最大流MATLAB源代码,文件名为mp_mc.m
function[Mm,mc,Mmr]=mp_mc(a,c)
A=a; %各路径最大承载流量矩阵
C=c; %各路径花费矩阵
Mm=0; %初始可行流设为零
mc=0; %最小花费变量
mcr=0;
mrd=0;
n=0;
while mrd~=inf %一直叠代到以花费为权值找不到最短路径
for i=1:(size(mcr',1)-1)
if a(mcr(i),mcr(i+1))==inf
ta=A(mcr(i+1),mcr(i))-a(mcr(i+1),mcr(i));
else
ta=a(mcr(i),mcr(i+1));
end
n=min(ta,n); %将最短路径上的最小允许流量提取出来
end
for i=1:(size(mcr',1)-1)
if a(mcr(i),mcr(i+1))==inf
a(mcr(i+1),mcr(i))=a(mcr(i+1),mcr(i))+n;
else
a(mcr(i),mcr(i+1))=a(mcr(i),mcr(i+1))-n;
end
end
Mm=Mm+n; %将每次叠代后增加的流量累加,叠代完成时就得到最大流量
for i=1:size(a,1)
for j=1:size(a',1)
if i