使用Scipy minimize求解约束优化问题时结果错误的原因排查

核心问题拆解 #

1. 约束条件的合理性澄清 #

你问题描述中写的约束是x+2y≥500，但这个约束下目标函数f(x,y)=xy不存在最大值——因为x和y可以无限增大，只要满足x+2y≥500，xy就能无限大。显然你实际要解决的是**x+2y≤500**的约束（这也是你代码注释里明确写的），这个约束下问题才有最优解。

2. 纸笔拉格朗日乘数法的错误 #

你得到的x≈333.4，y≈166.7是错误的，正确的拉格朗日求解过程如下：

• 最大化f(x,y)=xy，约束g(x,y)=x+2y-500=0（最优解必然在约束边界上）
• 构造拉格朗日函数：L = xy - λ(x+2y-500)
• 求偏导并令其为0：
  • ∂L/∂x = y - λ = 0 → y=λ
  • ∂L/∂y = x - 2λ = 0 → x=2λ
  • ∂L/∂λ = -(x+2y-500)=0 → x+2y=500
• 代入求解：将x=2λ、y=λ代入x+2y=500，得2λ+2λ=500 → λ=125，最终x=250，y=125

这个结果和Scipy返回的结果完全一致，说明你的代码是正确的，错误出在纸笔计算的拉格朗日推导过程中（大概率是误将约束的梯度项当成了(1,1)，对应x+y=500的约束，才会得到x≈333.4，y≈166.7的结果）。

3. 代码的正确性验证 #

你的代码逻辑是对的：

• 将最大化xy转化为最小化-xy，符合Scipy的API要求
• 约束条件定义正确：
  • c1函数实现了a+2b ≤500的约束（-(a+2b-500)≥0等价于a+2b≤500）
  • c2、c3实现了a>0、b>0的正性约束
• 初始值选择合理，求解过程收敛到了正确的最优解

结论 #

Scipy返回的x=250，y=125是正确的最优解，你的纸笔计算存在推导错误，代码本身没有问题。

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