CVXPY优化报错DCPError：目标函数不符合DCP规则求助

问题根源 #

出现cvxpy.error.DCPError: Problem does not follow DCP rules错误，核心原因有两个：

1. 语法错误：代码中Flow(Pressure*Pressure_Conversion)误用了函数调用语法，变量间乘法应使用*运算符，而非括号。
2. 违反DCP规则：修正语法后的目标函数为(Flow * Pressure) / Efficiency形式，属于变量乘积除以变量的非凸表达式，而CVXPY默认要求最小化问题的目标函数必须是凸函数，不符合DCP规范。

修复步骤 #

1. 修正基础语法错误 #

将目标函数中的函数调用改为变量乘法：

obj = Costofelectricity * (((Flow * (Pressure * Pressure_Conversion)) / (BlowerFactorConstant * Efficiency)) * Power_Conversion)

2. 处理非凸问题 #

针对非凸的目标函数，有两种可行解决路径：

路径一：使用非凸求解器

CVXPY支持部分非凸求解器（如Scipy的SLSQP、Gurobi等），只需在求解时指定即可：

# 使用Scipy的SLSQP求解器（需提前安装scipy）
results = prob.solve(solver=cp.SCIPY, method='SLSQP')

路径二：重新建模为凸问题（推荐）

通过变量替换将非凸问题转化为符合DCP规则的凸问题：

• 令Efficiency_inv = 1/Efficiency，原Efficiency ∈ [0.5, 0.7]对应Efficiency_inv ∈ [1/0.7, 1/0.5] ≈ [1.428, 2]
• 替换后目标函数变为线性凸表达式，所有CVXPY默认凸求解器均可处理

完整修正代码（路径二示例） #

import cvxpy as cp

# 定义优化变量
Efficiency_inv = cp.Variable()
Flow = cp.Variable()
Pressure = cp.Variable()

# 常量与转换系数
BlowerFactorConstant = 229
Pressure_Conversion = 1
Power_Conversion = 0.7457
Costofelectricity = 12.28

# 凸目标函数
obj = Costofelectricity * (Flow * Pressure * Pressure_Conversion * Power_Conversion) / BlowerFactorConstant * Efficiency_inv

# 约束条件
constraints = [
    Pressure >= 2.53, Pressure <= 7.6,
    Flow >= 550, Flow <= 1650,
    Efficiency_inv >= 1/0.7, Efficiency_inv <= 1/0.5,
    Flow >= 0, Pressure >= 0
]

# 构建并求解问题
objective = cp.Minimize(obj)
prob = cp.Problem(objective, constraints)
results = prob.solve()

print("最小化成本结果：", results)
print("最优流量：", Flow.value)
print("最优压力：", Pressure.value)
print("最优效率：", 1/Efficiency_inv.value)

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