最小费用最大流问题

合集 - 数学建模算法(27)

1.整数规划2024-07-072.非线性规划2024-10-113.非线性规划之飞行管理问题2024-10-114.遗传算法与直接搜索2024-10-115.图论基础：最小路径问题2024-07-076.网络最大流量问题2024-07-07

7.最小费用最大流问题2024-07-05

8.旅行商（TSP）问题2024-07-059.微分方程模型2024-07-0910.数据预处理2024-10-1111.插值方法2024-07-0612.拟合算法2024-07-0813.数据初始化2024-10-1114.层次分析法2024-07-1315.模糊综合评价2024-07-1416.CRITIC指标客观赋权方法2024-10-1117.Topsis评价法2024-10-1118.熵权法2024-10-1119.灰色关联度分析2024-10-1120.聚类分析2024-10-1121.主成分分析2024-10-1122.因子分析2024-10-1123.判别分析2024-10-1124.典型相关分析2024-10-1125.对应分析2024-10-1126.多维标度法2024-10-1127.亮度变化与空间滤波2024-10-11

收起

什么是最小费用最大流

最小费用最大流问题，是指在图中找到一条费用最小的途径，同时它的流量最大。
用数学表示为，
$x=(x_{ij})$表示一个可行流
$w=(w_{ij})$表示单位通过改弧的费用
那么一个可行流的总费用为$\sum _{(v_i,v_j)\in A}{w}_{ij}\cdot x_{ij}$

最小费最大流问题相比单纯的最大流问题还有一个限制条件，总费用要最小，这就要求在调整流量的同时要让费用尽可能小。

线性规划语言描述

用数学规划的语言来解最小费用最大流为：

所以我们可以先求最大流，然后在发点、收点流量已知的情况下用线性规划求解最小费用。

matlab代码实现

以这张网络图为例

（u,v）分别表示网络的流量和费用。

%% 通过最大流算法求出最大流量% 定义网络容量矩阵
networkCapacity = [
    0, 9, 5, 0, 0, 0;
    0, 0, 2, 0, 0, 5;
    0, 0, 0, 9, 0, 0;
    0, 6, 0, 0, 0, 10;
    8, 0, 7, 0, 0, 0;
    0, 0, 0, 0, 0, 0
];

% 定义源点和汇点
sourceNode = 5;
sinkNode = 6;

% 调用最大流算法
[maxFlowValue, flowMatrix] = maxflow(networkCapacity, sourceNode, sinkNode);%求解结果为14，后面用得到.

%% 线性规划求最小费用
% 定义目标函数的系数向量
costCoefficients = [2; 8; 5; 2; 3; 1; 4; 6; 7];

% 定义等式约束的系数矩阵和常数向量，中间节点出入流量和为0
equalityConstraintsMatrix = [
    1, 1, 0, 0, 0, 0, 0, 0, 0;
    1, 0, -1, -1, 0, 0, 0, 0, 0;
    0, 1, 1, 0, -1, 1, 0, 0, 0;
    0, 0, 0, 1, 0, -1, 1, -1, 0;
    0, 0, 0, 0, 1, 0, -1, 0, -1;
    0, 0, 0, 0, 0, 0, 0, 1, 1
];
equalityConstraintsVector = [14; 0; 0; 0; 0; 14];%s，t点的初始流量最大

% 定义变量的下界和上界
variableLowerBounds = [0; 0; 0; 0; 0; 0; 0; 0; 0];
variableUpperBounds = [8; 7; 5; 9; 9; 2; 6; 5; 10];

% 调用 linprog 函数求解线性规划问题
[optimalSolution, optimalObjectiveValue, exitFlag, outputStructure] = ...
    linprog(costCoefficients, [], [], equalityConstraintsMatrix, ...
    equalityConstraintsVector, variableLowerBounds, variableUpperBounds);

% 显示结果
fprintf('Optimal solution:\n'); disp(optimalSolution);
fprintf('Optimal objective value:\n'); disp(optimalObjectiveValue);
fprintf('Exit flag:\n'); disp(exitFlag);
fprintf('Output structure:\n'); disp(outputStructure);