使用CVXPY求解[1,10]区间最小化1/x触发DCPError的原因咨询

我尝试用CVXPY求解以下优化问题：

# Define the variable
x = cp.Variable()

# Define the objective function
objective = cp.Minimize(1/x)

# Define the constraints
constraints = [x >= 1, x <= 10]

# Formulate the problem
problem = cp.Problem(objective, constraints)

# Solve the problem
problem.solve()

# Print the results
print("Optimal value of x:", x.value)
print("Optimal minimum value:", 1/x.value)

运行后抛出如下错误：

DCPError: Problem does not follow DCP rules. Specifically:
The objective is not DCP. Its following subexpressions are not:
1.0 / var8856

我知道1/x并非全局凸函数，但它在区间[1,10]内是凸的，为何CVXPY仍会抛出DCPError？

解答 #

• CVXPY的DCP（Disciplined Convex Programming）规则是基于表达式符号与凸性属性的静态句法检查，不会结合变量的约束区间分析局部凸性，只判断表达式本身的全局凸性/凹性属性。
• 虽然1/x在x>0区间内是凸函数，但直接写1/x不符合DCP规则要求——CVXPY需要使用它预定义的、符合凸性规范的原子函数来表达目标或约束，而非普通算术除法。
• 正确写法是用cp.inv_pos(x)替代1/x，这个原子函数专门被定义为x>0时的凸函数，且你的约束x>=1满足它的定义域要求。修改后的目标函数应为：

  objective = cp.Minimize(cp.inv_pos(x))
  

• 本质上，DCP规则是一套严格的句法规范，目的是确保问题能被可靠求解，它不做“结合约束推导局部凸性”这类复杂动态分析，必须通过合法原子函数明确表达凸性。

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