多変数関数の極値とラグランジュの未定乗数法

ヘッセ行列による極値の判別、制約付き最適化を未定乗数法で解く。

この単元の目標
  • 二次偏導関数からヘッセ行列をつくり、極大・極小・鞍点を判別できる
  • 制約つき最大化・最小化をラグランジュの未定乗数法で解ける
  • 未定乗数法の幾何学的意味(勾配の平行条件)を説明できる

極値問題の準備——テイラー近似とヘッセ行列

開領域で二階偏導関数が連続な関数 f(x,y)f(x,y) を考えます。内点の極値候補は、一次の微分がすべて消える停留点 ∇f=0\nabla f = \mathbf{0} に限ります。その点が本当に山の頂上か谷底か、それとも馬の鞍のような鞍点かを見分けるのが二次微分の情報、すなわちヘッセ行列です。1 変数のとき「f′=0f^{\prime} = 0 かつ f′′>0f^{\prime\prime} > 0 なら極小」だったのと同じ役割を、行列が引き継いでいます。

H=(fxxfxyfyxfyy),D=det⁡H=fxxfyy−(fxy)2H = \begin{pmatrix} f_{xx} & f_{xy} \\ f_{yx} & f_{yy} \end{pmatrix}, \qquad D = \det H = f_{xx} f_{yy} - (f_{xy})^2
  • D>0D > 0 かつ fxx>0f_{xx} > 0:極小(ヘッセ行列が正定値=どの方向に歩いても増える)
  • D>0D > 0 かつ fxx<0f_{xx} < 0:極大(負定値=どの方向に歩いても減る)
  • D<0D < 0:鞍点(ある方向では増え、別の方向では減る)
  • D=0D = 0:この判定では決められない(高次の項を調べる必要がある)

🔍 なぜ判別式 D で極値がわかるのか——2 次のテイラー展開

  1. Step 1: 停留点のまわりで 2 次まで展開する
    Δf=f(a+h,b+k)−f(a,b)≈12(fxxh2+2fxyhk+fyyk2)\Delta f = f(a + h, b + k) - f(a, b) \approx \frac{1}{2}\left( f_{xx} h^2 + 2 f_{xy} hk + f_{yy} k^2 \right)

    停留点では一次の項(勾配)が消えるので、二次形式が正定値・負定値・不定符号なら、その符号が十分近くでの増減を決めます。退化する場合は高次項も必要です。

  2. Step 2: 平方完成する
    2Δf≈fxx(h+fxyfxxk)2+fxxfyy−(fxy)2fxxk22\Delta f \approx f_{xx} \left( h + \frac{f_{xy}}{f_{xx}} k \right)^2 + \frac{f_{xx} f_{yy} - (f_{xy})^2}{f_{xx}} k^2

    fxx≠0f_{xx}\ne0 の場合の平方完成です。fxx=0f_{xx}=0, D<0D<0 なら fxy≠0f_{xy}\ne0 なので、混合項から直接二次形式が正負両方の値をとると分かります。

  3. Step 3: 符号を読み取る
    D=fxxfyy−(fxy)2>0,  fxx>0⟹Δf>0 (任意の方向で増える=極小)D = f_{xx} f_{yy} - (f_{xy})^2 > 0, \; f_{xx} > 0 \quad\Longrightarrow\quad \Delta f > 0 \text{ (任意の方向で増える=極小)}

    二つの平方項は非負で、(h,k)≠(0,0)(h,k)\ne(0,0) なら少なくとも一方は正です。二次形式の正定値性により十分小さな変位では高次の剰余より二次項が支配的となり、極小が従います。D<0D<0 なら方向によって二次形式の符号が変わり鞍点です。

ラグランジュの未定乗数法

∇f=λ∇g,g(x,y)=0\nabla f = \lambda \nabla g, \qquad g(x, y) = 0

🔍 なぜ「勾配が平行」でよいのか——等高線の幾何学

  1. Step 1: 制約は曲線、目的関数は等高線
    制約 g(x,y)=0 は曲線、f(x,y)=c は等高線族\text{制約 } g(x, y) = 0 \text{ は曲線、} f(x, y) = c \text{ は等高線族}

    曲線上を動きながら f の値をできるだけ大きくしたい状況を考える。

  2. Step 2: 勾配は等高線に垂直
    ∇f⊥f(x,y)=c,∇g⊥g(x,y)=0\nabla f \perp f(x,y)=c, \qquad \nabla g \perp g(x,y)=0

    勾配ベクトルの定義から、関数が最も急増加する方向は等高線の法線方向。

  3. Step 3: 制約曲線に沿った方向微分が消える
    ddtf(r(t))=∇f⋅r′(t)=0,r′(t)∥接線\frac{d}{dt} f(\mathbf{r}(t)) = \nabla f \cdot \mathbf{r}^{\prime}(t) = 0, \qquad \mathbf{r}^{\prime}(t) \parallel \text{接線}

    制約曲線に沿って動いている間、極值では f の値が一時的に変化しなくなる。つまり ∇f は曲線の接線と直交。

  4. Step 4: ∇g が非零なら、法線方向はそのスカラー倍
    ∇f∥∇g⟺∇f=λ∇g\nabla f \parallel \nabla g \quad\Longleftrightarrow\quad \nabla f = \lambda \nabla g

    f,gf,g が連続微分可能で ∇g≠0\nabla g\ne0 の点では、∇g は制約曲線の法線だから、∇f はそのスカラー倍。比例定数 λ が「未定乗数」。これと制約式を連立すれば未知数 3 個の連立方程式から候補を探せます。必要条件なので、候補の値の比較や特異点・境界の確認を別に行います。

例題(制約つき最大値)

周囲の長さが 20 m の長方形のうち、面積を最大にするものを未定乗数法で求めよ。

解答・解説

f(x,y)=xyf(x,y) = xy, 制約 g=x+y−10=0g = x + y - 10 = 0。∇f=λ∇g\nabla f = \lambda \nabla g より y=λy = \lambda, x=λx = \lambda、ゆえに x=y=5x = y = 5。制約を代入すると xy=x(10−x)=25−(x−5)2≤25xy=x(10-x)=25-(x-5)^2\le25。よって面積は **最大 25 平方メートル**(正方形が最適)。

練習問題
  1. f(x,y)=x3+y3−3xyf(x, y) = x^3 + y^3 - 3xy の極値を求めよ。

    ヒント

    fx=3x2−3y=0f_x = 3x^2 - 3y = 0, fy=3y2−3x=0f_y = 3y^2 - 3x = 0。実数解は (0,0)(0,0) と (1,1)(1,1)。

    解答・解説

    D=fxxfyy−fxy2=6x⋅6y−(−3)2D = f_{xx} f_{yy} - f_{xy}^2 = 6x \cdot 6y - (-3)^2。(0,0)(0,0) では D=−9<0D = -9 < 0 で**鞍点**。(1,1)(1,1) では D=36−9=27>0D = 36 - 9 = 27 > 0, fxx=6>0f_{xx} = 6 > 0 で**極小値 −1-1**

  2. f(x,y)=x+yf(x, y) = x + y を円周 x2+y2=4x^2 + y^2 = 4 上で最大・最小にせよ。

    ヒント

    1=λ2x1 = \lambda 2x, 1=λ2y1 = \lambda 2y より x=yx = y。

    解答・解説

    x=y=2x = y = \sqrt{2} で**最大 222\sqrt{2}**、x=y=−2x = y = -\sqrt{2} で**最小 −22-2\sqrt{2}**。(コーシー・シュワルツの等号条件と一致)

  3. 体積 V が一定の直方体のうち表面積を最小にするものの形状を述べよ。

    ヒント

    xyz=Vxyz = V の下で 2(xy+yz+zx)2(xy + yz + zx) を最小化。

    解答・解説

    **立方体(すべての辺が等しい)**。対称性から x=y=zx = y = z が解となり、これは不等式 (xy+yz+zx)/3≥(x2y2z2)1/3(xy + yz + zx)/3 \ge (x^2y^2z^2)^{1/3} の等号ケースにも対応する。

確認クイズ(極値と未定乗数法)

問1. 停留点で D < 0 のとき、その点は何と呼ばれるか。

問2. ラグランジュの未定乗数法の条件 ∇f = λ∇g は何を表すか。

問3. 二階偏導関数が連続な関数の停留点で、ヘッセ行列が正定値ならどうなるか。

📝 メモこの端末のみに保存
メモ欄はページ読み込み後に使えるようになります。