拉格朗日乘数法深度解析

拉格朗日乘数法的几何本质

拉格朗日乘数的含义

定理1

乘数的灵敏度解释:设 f(c)f^*(c) 是约束 g(x,y)=cg(x,y)=cff 的极值,则

λ=dfdc\lambda = -\frac{df^*}{dc}

即拉格朗日乘数 λ\lambda(取负号后)表示约束右端项 cc 变化一个单位时,目标函数极值的变化率。

在经济学中,λ\lambda 称为影子价格(shadow price),表示资源的边际价值。

几何解释
推论
证明
符号说明

不等式约束与KKT条件

定理2

对于不等式约束优化问题:

minf(x)s.t.gi(x)0(i=1,,m)\min f(x) \quad \text{s.t.} \quad g_i(x) \leq 0 \quad (i=1,\ldots,m)

KKT条件(Karush-Kuhn-Tucker条件)是极值点的必要条件:

  1. 稳定性f(x)+i=1mλigi(x)=0\nabla f(x^*) + \sum_{i=1}^m \lambda_i \nabla g_i(x^*) = 0
  2. 原始可行性gi(x)0g_i(x^*) \leq 0
  3. 对偶可行性λi0\lambda_i \geq 0
  4. 互补松弛λigi(x)=0\lambda_i g_i(x^*) = 0

互补松弛条件意味着:如果约束 gi<0g_i < 0(不起作用),则 λi=0\lambda_i = 0;如果 λi>0\lambda_i > 0,则 gi=0g_i = 0(约束起作用)。

几何解释
推论
证明
符号说明

拉格朗日对偶

拉格朗日对偶函数

原问题:minf(x)\min f(x) s.t. gi(x)0,hj(x)=0g_i(x) \leq 0, h_j(x)=0

拉格朗日函数:L(x,λ,μ)=f(x)+λigi(x)+μjhj(x)L(x, \lambda, \mu) = f(x) + \sum \lambda_i g_i(x) + \sum \mu_j h_j(x)

对偶函数g(λ,μ)=infxL(x,λ,μ)g(\lambda, \mu) = \inf_x L(x, \lambda, \mu)

对偶问题maxλ0,μg(λ,μ)\max_{\lambda \geq 0, \mu} g(\lambda, \mu)

弱对偶定理:对偶问题的最优值 \leq 原问题的最优值。

强对偶定理:在凸优化问题中(满足约束规范),对偶问题最优值 = 原问题最优值。

典型例题

例题1:影子价格的计算

在约束 x+y=cx + y = c 下最大化 f(x,y)=xyf(x,y) = xy,求最优值 f(c)f^*(c) 和拉格朗日乘数 λ\lambda,验证 λ=df/dc\lambda = -df^*/dc

参考答案(3 个标签)
拉格朗日乘数影子价格灵敏度
  1. L=xy+λ(x+yc)L = xy + \lambda(x+y-c)
  2. Lx=y+λ=0,Ly=x+λ=0L_x = y+\lambda=0, L_y = x+\lambda=0,得 x=y=λx=y=-\lambda
  3. 代入约束:2x=c2x=cx=y=c2x=y=\frac{c}{2}λ=c2\lambda = -\frac{c}{2}
  4. 最优值 f(c)=c24f^*(c) = \frac{c^2}{4}
  5. dfdc=c2\frac{df^*}{dc} = \frac{c}{2},而 λ=c2-\lambda = \frac{c}{2},验证了 λ=dfdc\lambda = -\frac{df^*}{dc}

答案f(c)=c24f^*(c)=\frac{c^2}{4}λ=c2\lambda=-\frac{c}{2},满足 λ=dfdc\lambda=-\frac{df^*}{dc}

例题2:KKT条件应用

用KKT条件求解 minf(x,y)=x2+y2\min f(x,y) = x^2 + y^2 s.t. x+y1x + y \geq 1

参考答案(3 个标签)
KKT条件不等式约束二次规划
  1. 标准化约束:g(x,y)=1xy0g(x,y) = 1 - x - y \leq 0
  2. L=x2+y2+λ(1xy)L = x^2+y^2 + \lambda(1-x-y)
  3. KKT条件:
    • L=(2xλ,2yλ)=0\nabla L = (2x-\lambda, 2y-\lambda) = 0x=y=λ2x=y=\frac{\lambda}{2}
    • 1xy01-x-y \leq 0
    • λ0\lambda \geq 0
    • λ(1xy)=0\lambda(1-x-y) = 0
  4. λ=0\lambda=0x=y=0x=y=0,但 100=1>01-0-0=1>0,违反约束可行性。
  5. λ>0\lambda>0,由互补松弛 1xy=01-x-y=0,即 x+y=1x+y=1
  6. 结合 x=y=λ2x=y=\frac{\lambda}{2}2λ2=12\cdot\frac{\lambda}{2}=1λ=1\lambda=1x=y=12x=y=\frac{1}{2}
  7. 最小值 f(12,12)=12f(\frac{1}{2},\frac{1}{2}) = \frac{1}{2}

答案:最小值为 12\frac{1}{2},在 (12,12)(\frac{1}{2},\frac{1}{2}) 处取得。


总结

本文出现的符号

符号类型读音/说明在本文中的含义
λ\lambda拉格朗日乘数lambda等式/不等式约束的乘数
μ\mu乘数mu等式约束的乘数
LL拉格朗日函数Lagrangian拉格朗日函数
g(λ,μ)g(\lambda,\mu)对偶函数dual function拉格朗日对偶函数
ff^*最优值optimal value约束优化的最优值
cc约束右端constraint RHS约束条件的右端项

中英对照

中文术语英文术语音标
拉格朗日乘数法Lagrange multiplier method/ləˈɡrɑːndʒ ˈmʌltɪplaɪər ˈmɛθəd/
KKT条件KKT conditions/ˌkeɪ keɪ ˈtiː kənˈdɪʃənz/
影子价格shadow price/ˈʃædoʊ praɪs/
互补松弛complementary slackness/ˌkɒmplɪˈmɛntəri ˈslæknəs/
对偶理论duality theory/djuːˈæləti ˈθɪəri/
弱对偶weak duality/wiːk djuːˈæləti/
强对偶strong duality/strɒŋ djuːˈæləti/
不等式约束inequality constraint/ˌɪnɪˈkwɒləti kənˈstreɪnt/
起作用约束active constraint/ˈæktɪv kənˈstreɪnt/
灵敏度分析sensitivity analysis/ˌsɛnsəˈtɪvəti əˈnæləsɪs/
非线性规划nonlinear programming/ˌnɒnˈlɪniər ˈproʊɡræmɪŋ/
凸优化convex optimization/ˈkɒnvɛks ˌɒptɪmaɪˈzeɪʃən/