如何解决MATLAB优化报错"无法从OptimizationConstraint转逻辑值"

报错原因 #

MATLAB基于问题的优化框架中，Ppve、Pelet都属于未赋值的优化变量，建模阶段Ppve(i) >= Pelet(i)返回的是优化约束对象，不是可判断真假的布尔值，因此不能直接用于普通if逻辑判断。

解决方法 #

使用整数规划的大M约束法，引入二进制辅助变量实现分支逻辑，不需要用if判断，直接把逻辑转化为线性约束即可。

具体实现步骤 #

1. 新增二进制优化变量，以及足够大的常数M（取值大于所有变量的最大可能差值即可）
2. 将你的条件逻辑转化为4组线性约束：
   • 当二进制变量z(i)=1时，满足Ppve(i) >= Pelet(i)且Pbe(i)==0
   • 当二进制变量z(i)=0时，满足Ppve(i) < Pelet(i)且Pbe(i)>0

修改后代码示例 #

变量声明部分新增代码

% 新增二进制辅助变量，z取值只能是0或1
z = optimvar('z',1,n,'Type','integer','LowerBound',0,'UpperBound',1);
% 定义大M，取足够大的正数即可，这里取光伏和负载上限的2倍
M = max([daily_Ppv, Pemax]) * 2;

约束部分替换原BAT约束的循环代码

% 初始化4组逻辑约束
PowerFlow.Constraints.BATLogic1 = optimconstr(n);
PowerFlow.Constraints.BATLogic2 = optimconstr(n);
PowerFlow.Constraints.BATLogic3 = optimconstr(n);
PowerFlow.Constraints.BATLogic4 = optimconstr(n);

for i = 1:n
    % z=1时生效：Ppve(i) >= Pelet(i)
    PowerFlow.Constraints.BATLogic1(i) = Ppve(i) - Pelet(i) >= -M*(1 - z(i));
    % z=1时生效：Pbe(i) <=0，结合你原本Pbe的下界0，等价于Pbe(i)==0
    PowerFlow.Constraints.BATLogic2(i) = Pbe(i) <= M*(1 - z(i));
    % z=0时生效：Pelet(i) > Ppve(i)
    PowerFlow.Constraints.BATLogic3(i) = Pelet(i) - Ppve(i) >= -M*z(i);
    % z=0时生效：Pbe(i) >=极小值，等价于你要的Pbe(i)>0，阈值可根据精度调整
    PowerFlow.Constraints.BATLogic4(i) = Pbe(i) >= 1e-3*(1 - z(i));
end

注意事项 #

• 引入二进制变量后，需要选择支持整数规划的求解器（如intlinprog、GA等），不能使用纯连续变量求解器
• 1e-3是Pbe的最小取值阈值，可根据你的工程精度要求调整，只要是大于0的极小值即可

内容的提问来源于stack exchange，提问作者Lorenzo Vallerini