手把手玩转Benders分解：两阶段鲁棒优化的暴力美学-编程阁

基于benders分解算法的两阶段鲁棒问题求解关键词：两阶段鲁棒 benders分解法鲁棒优化参考文档：《Solving two-stage robust optimization problems using a column-and-constraint generation method》（问题背景是这个文献，benders分解过程见CSDN）仿真平台：MATLAB YALMIP+CPLEX 优势：代码注释详实，适合参考学习，非目前烂大街的微网两阶段规划版本，请仔细辨识！主要内容：代码构建了两阶段鲁棒优化模型，并用文档中的相对简单的算例，进行benders分解算法的验证，此篇文献是benders分解算法的入门级文献，其经典程度不言而喻，几乎每个搞benders分解算法的两阶段鲁棒的人都绕不过此篇文献，所以萌新们或者新手们赶紧冲起来学习吧！注意这个编程语言是MATLAB

最近被实验室师兄按头安利了这篇Benders分解的经典文献，说是搞鲁棒优化的没人能绕开它。亲自撸完代码后发现——真香！今天咱们就用最粗暴的方式，把这个优雅的数学算法扒个精光。

先搞点前戏：两阶段问题长啥样？

想象你是个厂长，先要决定生产线布局（第一阶段），然后根据实际市场需求调整生产计划（第二阶段）。但坑爹的是市场需求充满不确定性，这时候就得用鲁棒优化让方案在最差情况下也能扛得住。

Benders分解的核心就一句话：主问题挖坑，子问题填土

主问题负责生成"看起来很美"的初始方案，子问题专门找茬——验证这个方案在不确定参数的最坏情况下会不会翻车。要是翻车了，就生成一个Benders割平面扔回主问题，逼它重新做人。

上代码！先看主问题怎么造作

%% 主问题建模 y = sdpvar(3,1); % 第一阶段决策变量 theta = sdpvar(1); % 辅助变量 Constraints = [y >= 0, sum(y) <= 100]; % 基础约束 Objective = 2*y(1) + 3*y(2) + theta; % 目标函数

这里的阴险之处在于theta变量，它像个骗子一样假装自己代表了第二阶段的最优值。子问题的工作就是不断戳穿这个骗局，逼theta现出原形。

子问题才是真·魔鬼

function [cut, status] = solve_sub(y) u = sdpvar(2,1); % 不确定参数 x = sdpvar(2,1); % 第二阶段决策 % 不确定性集合：||u||_1 <= 3 Uncertainty = [u >= -1, u <= 2, sum(abs(u)) <= 3]; % 第二阶段约束 SubConstraints = [x >= 0, x <= 10, 3*x(1) + 2*x(2) <= 15 + y'*[1;2;3]]; SubObjective = -(4*x(1) + 5*x(2) + [0.5;1.5]'*u); % 注意符号取反 options = sdpsettings('verbose',0,'solver','cplex'); optimize([SubConstraints,Uncertainty], SubObjective, options); % 提取对偶变量生成cut dual_values = getdual(SubConstraints); cut = (theta >= ... ); % 具体公式参考文献 end

子问题求解时有个骚操作：把原问题的max-min转换成单层优化。这里用到了强对偶定理，相当于给不确定性参数u套上紧箍咒，让它再怎么作妖也翻不出五指山。

主从调教的迭代过程

while abs(UB - LB) > 1e-4 % 解主问题 optimize(MasterConstraints, MasterObjective); y_val = value(y); LB = value(MasterObjective); % 解子问题 [new_cut, status] = solve_sub(y_val); UB = min(UB, 2*y_val(1)+3*y_val(2)+value(theta)); % 添加新约束 MasterConstraints = [MasterConstraints, new_cut]; end

这个循环就像在玩"大家来找茬"：主问题每次给出一个方案，子问题就变身杠精，专门找最恶劣的工况来打脸。每次迭代都会在方案里埋下针对这种恶劣工况的防御措施。

新手必踩的三大坑：