拉格朗日乘数法、条件极值

                     

贡献者: addis; Connection

预备知识 1 导数与函数极值,梯度

   若要求函数 $f(x_1,\dots, x_N)$ 在 $m$ 个约束条件

\begin{equation} g_i(x_1, \dots, x_N) = 0 \qquad (i = 1,\dots, M)~ \end{equation}
下的极值,可用以下的拉格朗日乘数法

   令拉格朗日函数为

\begin{equation} \Lambda(x_1,\dots, x_N) = f + \sum_{i=1}^M \lambda_i g_i~, \end{equation}
其中 $\lambda_i$ 为未知任意常数。同时满足约束条件和
\begin{equation} \frac{\partial \Lambda}{\partial x_i} = 0 \quad (i = 1,\dots, N)~ \end{equation}
的点,就是函数 $f(x_1, \dots, x_N)$ 的极值点稳定点

   注意式 3 有 $N$ 条等式,结合式 1 共有 $M+N$ 条等式组成方程组(注意未必是线性的)。而未知量 $x_1,\dots,x_N$ 和 $\lambda_1,\dots,\lambda_M$ 也有 $M+N$ 个。所以一般情况下存在唯一解或有限个解。

1. 几何理解

   下面举一个二元函数的例子以说明拉格朗日乘数法的几何意义。

图
图 1:例 1

例 1 

   求函数 $f(x,y) = -a(x^2 + y^2)$ 在轨迹 $x+y = c$ 约束下的最值($a,c$ 均为常数)。

   从几何上来理解,函数 $f(x,y)$ 是一个倒置的抛物面,其等高线的俯视图如图 1 蓝色的箭头代表函数值增加最大的方向,即梯度的方向。考察点沿着同一条等高线的位置变化时,函数值不变。曲面上任意一点的梯度为

\begin{equation} \boldsymbol\nabla f = \frac{\partial f}{\partial x} \hat{\boldsymbol{\mathbf{x}}} + \frac{\partial f}{\partial y} \hat{\boldsymbol{\mathbf{y}}} = -2ax \hat{\boldsymbol{\mathbf{x}}} -2ay \hat{\boldsymbol{\mathbf{y}}} = -2a \boldsymbol{\mathbf{r}} ~, \end{equation}
其中 $ \boldsymbol{\mathbf{r}} $ 为位矢。这说明,梯度的方向总是延径向向内,与距离成正比。

   约束条件如图 1 中红线所示,我们要求的就是红线上面函数 $f$ 的极值。在该题的情况下,所求的极值显然在红线离原点最近的地方出现。

   但是如何总结出约束条件下多元函数极值点的一般的判断条件并用梯度表示呢?根据梯度的性质(式 5 ),如果考察点沿红线的切向移动 $ \,\mathrm{d}{ \boldsymbol{\mathbf{l}} } $,函数值的微分为

\begin{equation} \,\mathrm{d}{f} = \boldsymbol\nabla f \boldsymbol\cdot \,\mathrm{d}{ \boldsymbol{\mathbf{l}} } ~. \end{equation}
使用导数和极值的关系,我们要找的极值点满足:函数在极值点处延红线切向的方向导数为零。即
\begin{equation} \boldsymbol\nabla f \boldsymbol\cdot \,\mathrm{d}{ \boldsymbol{\mathbf{l}} } = 0~. \end{equation}
这说明 $ \boldsymbol\nabla f$ 只可能具有极值点处红线的法向量分量而切向分量为 0,即 $ \boldsymbol\nabla f$ 与法向量共线。令约束条件即红线的方程为
\begin{equation} g(x, y) = x+y-c = 0~, \end{equation}
$g(x,y)$ 的梯度就是其法向量。由矢量共线的定义,存在常数 $\lambda$,使
\begin{equation} \boldsymbol\nabla f + \lambda \boldsymbol\nabla g = 0 \qquad \text{即} \qquad \boldsymbol\nabla (f + \lambda g) = 0~, \end{equation}
这就是拉格朗日乘数法。其中 $f + \lambda g$ 就是拉格朗日函数。下面根据该条件,对本题进行求解。本题中拉格朗日函数梯度的两个分量分别为
\begin{equation} \begin{aligned} 0 &= \frac{\partial}{\partial{x}} (f + \lambda g) = -2ax + \lambda~,\\ 0 &= \frac{\partial}{\partial{y}} (f + \lambda g) = -2ay + \lambda~. \end{aligned} \end{equation}
再结合约束条件 $x + y - c = 0$,得到 $x,y,\lambda$ 三个未知量的线性方程组。解得极值点为 $(x,y)=(c/2, c/2)$,$\lambda=ac$。

   以上的推导可以很容易地拓展到多元函数和多个限制条件中去。例如求三元函数 $f(x,y,z)$ 中(可看作空间的标量场),在两个约束条件 $g_1(x,y,z) = 0$,$g_2(x,y,z) = 0$ 下的极值。注意两个约束条件分别代表一个曲面,同时成立代表两个曲面相交的空间曲线。我们要在该曲线上求函数 $f$ 的极值。同样,若曲线上的某点是极值,那必须满足曲线在该点的切线方向 $f$ 的方向导数为零,即 $ \boldsymbol\nabla f$ 与该点法向量共线。然而这是个三维问题,空间曲线的任意一点可以有两条不共线的法向量(不一定互相垂直),该点的任何其他法向量可以由这两个法线的线性组合来表示。即存在常数 $\lambda_1$ 和 $\lambda_2$,使得

\begin{equation} \boldsymbol\nabla f + \lambda_1 \boldsymbol\nabla g_1 + \lambda_2 \boldsymbol\nabla g_2 = 0~. \end{equation}
\begin{equation} \boldsymbol\nabla (f + \lambda_1 g_1 + \lambda_2 g_2) = 0~, \end{equation}
其中 $f + \lambda_1 g_1 + \lambda_2 g_2$ 就是拉格朗日函数。当然,限制条件的数目不一定是 2,若只有一个限制条件 $g_1 = 0$,就是求曲面 $g_1 = 0$ 上函数 $f$ 取得极值的点,条件是函数的梯度 $ \boldsymbol\nabla f$ 与曲面的法向量 $ \boldsymbol\nabla g_1$ 共线。

例 2 

   椭圆方程为 $x^2/a^2 + y^2/b^2 = 1$,求其内接长方型的最大面积(长方形的边平行于 $x$ 轴和 $y$ 轴)。

   令长方型的右上角与椭圆的切点为 $(x,y)$,则长方型的面积为 $S = 4xy$。约束函数为 $g = x^2/a^2 + y^2/b^2 - 1$。拉格朗日函数为

\begin{equation} \Lambda(x,y) = S + \lambda g = 4xy + \lambda \left(\frac{x^2}{a^2} + \frac{y^2}{b^2} - 1 \right) ~. \end{equation}
令 $x, y$ 偏导都为 0,得
\begin{equation} \begin{cases} 4y + 2\lambda x/a^2 = 0\\ 4x + 2\lambda y/b^2 = 0 \end{cases}~, \end{equation}
两式消去 $\lambda$,得 $x/y = a/b$。再联立椭圆方程,解得 $(x,y) = (a/\sqrt{2}, b/\sqrt{2})$,面积最大值为 $2ab$。

2. 物理意义

预备知识 2 力场、保守场、势能

   如果从力学角度审视约束下质点运动的稳定位置,拉格朗日乘数法将会显得十分自然。1

   考虑一个特例。平面 $R^2$ 内一条曲线 $C:f(x,y)=0$ 上的轨道,轨道上有一个质点。

   假设 $f$ 可微,力场为势场(例如,引力场),忽略摩擦、空气阻力等非理想条件。

   显然,轨道对质点的约束力 $R$ 必沿法线方向。由于曲线 $C$ 的法向量有形式 $( \partial f/\partial x , \partial f/\partial y )$,由向量共线的定义,我们有理由设约束力为

\begin{equation} \boldsymbol{\mathbf{R}} = \left(\lambda \frac{\partial f}{\partial x} ,\lambda \frac{\partial f}{\partial y} \right) ~. \end{equation}

   既然力场是势场,这里假定一个引力场,那么我们可以将引力表示为引力势的梯度,即

\begin{equation} \boldsymbol{\mathbf{F}} = - \boldsymbol\nabla p = \left(- \frac{\partial p}{\partial x} ,- \frac{\partial p}{\partial y} \right) ~. \end{equation}
$p=p(x,y)$ 是 $(x,y)$ 处的势能。当势能取极值时,质点处于平衡位置(合外力 $ \boldsymbol{\mathbf{R}} + \boldsymbol{\mathbf{F}} =0$ 的位置),换言之,
\begin{equation} \boldsymbol{\mathbf{R}} + \boldsymbol{\mathbf{F}} = 0 ~. \end{equation}
写成关于分量的方程组
\begin{equation} \lambda \frac{\partial f}{\partial x} - \frac{\partial p}{\partial x} = 0~, \quad \lambda \frac{\partial f}{\partial y} - \frac{\partial p}{\partial y} = 0 ~. \end{equation}
这正是拉格朗日乘数法的形式,也即条件极值的必要条件。


1. ^ 参见卓里奇《数学分析》第八章第七节第 12 题


致读者: 小时百科一直以来坚持所有内容免费无广告,这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。 因此,我们请求广大读者热心打赏 ,使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 20 元,我们一周就能脱离亏损, 并在接下来的一年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款, 他们的付出帮助了 99% 的读者免费获取知识, 我们在此表示感谢。

                     

友情链接: 超理论坛 | ©小时科技 保留一切权利