# 逆関数定理・陰関数定理と条件付き極値：縮小写像で「方程式を解く」

> 微分が可逆なら写像は局所的に可逆である、という逆関数定理を縮小写像の原理で証明し、そこから陰関数定理とラグランジュの乗数法を導く。反例と図で各仮定の役割を確かめます。
> https://rikai.mugen-giken.com/mathematics/calculus/inverse-implicit-functions

## 0. この記事の要点

連立方程式を「解けるか、解けるなら一通りか」と問われたとき、式が一次式なら行列を調べれば答えが出ます。式が曲がっていても、一点のすぐ近くだけなら、その点での一次近似（接線や接平面の考え方）で同じことが言えるのではないか、というのがこの記事の話です。読み終えると、方程式を実際に解かなくても、どの点の近くで「ある変数を別の変数の関数として表せるか」を微分で判定でき、制約付きの最大・最小を求める乗数法の理由も分かります。

- **逆関数定理**: $C^1$ 級写像 $F:\mathbb{R}^n \supset U \to \mathbb{R}^n$ の点 $a$ でのヤコビ行列 $DF(a)$ が正則なら、$F$ は $a$ の近くで $C^1$ 級の逆写像を持ちます。「一次近似が可逆なら、元の写像も局所的に可逆」という主張です。
- 証明は**縮小写像の原理**で行います。方程式 $F(x)=y$ を「$x=\Phi_y(x)$ という不動点問題」に書き換えると、$y$ が $F(a)$ に近い限り $\Phi_y$ が縮小写像になり、解が存在して一意に決まります。
- **陰関数定理**は逆関数定理の言い換えです。$F(x,y)=0$ で $\partial F/\partial y$ が可逆なら、その点の近くで $y$ を $x$ の $C^1$ 級関数として解け、$Dg=-(F_y)^{-1}F_x$ が成り立ちます。
- **ラグランジュの乗数法**は陰関数定理の応用です。制約 $g=0$ の下での極値点では $\nabla f=\sum \lambda_i \nabla g_i$ が成り立ちます。ただし $\nabla g_i$ が一次独立という仮定を落とすと反例があります。
- 導関数が連続でない例 $x+2x^2\sin(1/x)$ は、$f'(0)=1\neq 0$ なのに $0$ の近傍で単射になりません。$C^1$ 級という仮定は本質的です。

<Aside type="tip" title="この記事で使う事前知識">

- **ヤコビ行列と $C^1$ 級**: 多変数写像の各成分を各変数で偏微分して並べた行列が、その点での一次近似です。偏導関数がすべて連続なら $C^1$ 級です。→ [多変数関数の微分と偏微分](/mathematics/calculus/multivariable-differentiation)
- **行列式と正則行列**: 正方行列の行列式が $0$ でないとき、逆行列が存在します（正則）。→ [行列式](/mathematics/linear-algebra/determinants#def-determinant)
- **平均値の定理**: 区間上の微分可能な関数について、増え方の平均と等しい傾きを持つ点が間にある、という定理です。→ [平均値の定理とテイラーの定理](/mathematics/calculus/mean-value-and-taylor)
- **コーシー列**: 項どうしの距離がいくらでも小さくなる点列のことで、実数ではこれが必ず収束します。→ [コーシー列](/mathematics/real-analysis/completeness-and-cauchy#def-cauchy)
- **コンパクト**: 有界閉集合上の連続関数が最大値・最小値をとる、という性質の土台になる概念です。→ [コンパクト](/mathematics/real-analysis/metric-spaces#def-compactness-metric)

</Aside>

## 1. 動機：「方程式は解けるか」を微分で判定する

連立方程式 $F(x)=y$ を考えます。$F$ が線形、$F(x)=Ax$ なら、答えは[行列式](/mathematics/linear-algebra/determinants#def-determinant)で決まります。$\det A\ne 0$ なら任意の $y$ に対して $x=A^{-1}y$ と一意に解け、$\det A=0$ なら解けない $y$ や解が一意でない $y$ が出ます。

非線形の $F$ ではどうでしょう。円 $x^2+y^2=1$ を「$y$ について解く」ことを考えます。点 $(0,1)$ の近くでは $y=\sqrt{1-x^2}$ と解けますが、円全体を 1 つの関数 $y=g(x)$ で表すことはできません。また点 $(1,0)$ の近くでは、$x$ を決めても $y=\pm\sqrt{1-x^2}$ の 2 つの値があり、そもそも $y$ を $x$ の関数として書けません。**解ける場所と解けない場所**を、方程式を実際に解かずに見分けたい、というのがこの記事の出発点です。

答えは 1 次近似にあります。円の式 $F(x,y)=x^2+y^2-1$ では $\partial F/\partial y=2y$ で、これが $0$ でない点（$y\ne0$）では局所的に $y$ について解け、$0$ になる点（$(\pm1,0)$）では解けません。**1 次近似（線形写像）が可逆なら、元の非線形方程式も局所的には同様に解ける**。これが逆関数定理・陰関数定理の中身です。

用途は広く、曲面や多様体を方程式 $g=0$ で定義したときに、それが局所的に「なめらかな図形」であること（[多様体入門](/mathematics/manifolds/smooth-manifolds) の出発点）や、制約付き最適化のラグランジュの乗数法の根拠になります。この記事では、これら 3 つの定理を、縮小写像の原理という 1 つの道具から順に証明します。

## 2. 準備：ノルム、ヤコビ行列、平均値不等式

$\mathbb{R}^n$ の元を列ベクトルとし、$\|x\|=\sqrt{x_1^2+\cdots+x_n^2}$ をユークリッドノルムとします。$m\times n$ 行列 $A=(a_{ij})$ に対して

$$
\|A\| = \Bigl(\sum_{i,j} a_{ij}^2\Bigr)^{1/2}
$$
とおきます（フロベニウスノルム）。コーシー・シュワルツの不等式（<Ref to="mathematics/linear-algebra/inner-product-spaces#thm-cauchy-schwarz" />）を各行に適用すると $(Ax)_i$ の絶対値が第 $i$ 行のノルムと $\|x\|$ の積以下になることが分かります。これを 2 乗して $i$ について足せば

$$
\|Ax\|\le \|A\|\,\|x\|, \qquad \|AB\|\le \|A\|\,\|B\|
$$
が得られます（2 つ目は $AB$ の各列に 1 つ目を使う）。以後、行列のノルムとはこれを指します。

<Definition id="def-jacobian-map" title="ヤコビ行列と $C^1$ 級写像">
$U\subset\mathbb{R}^n$ を開集合、$F=(F_1,\dots,F_m):U\to\mathbb{R}^m$ とします。ある $m\times n$ 行列 $A$ が存在して

$$
\lim_{h\to 0}\frac{\|F(a+h)-F(a)-Ah\|}{\|h\|}=0
$$
が成り立つとき、$F$ は $a$ で**微分可能**といい、$A$ を $DF(a)$ と書きます。このとき $DF(a)$ の $(i,j)$ 成分は偏微分係数 $\partial F_i/\partial x_j(a)$ です。すべての偏導関数が $U$ 上で存在して連続なとき、$F$ は $C^1$ 級であるといいます。
</Definition>

2 変数の場合の議論（<Ref to="mathematics/calculus/multivariable-differentiation#def-differentiable-n" /> と <Ref to="mathematics/calculus/multivariable-differentiation#thm-c1-implies-differentiable-n" />）は、変数の個数にも成分の個数にも依らずそのまま通ります。以下では、次の 3 つを $n$ 変数・ベクトル値の場合にも成り立つ事実として使います。

1. $C^1$ 級ならば各点で微分可能で、$x\mapsto DF(x)$ は（行列成分の連続関数として）連続である。
2. 連鎖律（<Ref to="mathematics/calculus/multivariable-differentiation#thm-chain-rule-n" />）: $F$ が $a$ で、$G$ が $F(a)$ で微分可能なら $D(G\circ F)(a)=DG(F(a))\,DF(a)$。
3. 微分可能な $F$ について、偏微分係数の行列 $DF(a)$ は一意に定まる。

次の補題は、1 変数の平均値の定理 <Ref to="mathematics/calculus/mean-value-and-taylor#thm-mvt" /> をベクトル値に拡張したものです。ベクトル値関数では「$F(y)-F(x)=DF(c)(y-x)$ となる $c$」は一般に存在しませんが、**不等式の形**なら成り立ちます。

<Lemma id="lem-mean-value-ineq" title="ベクトル値関数の平均値不等式">
$U\subset\mathbb{R}^n$ を開集合、$F:U\to\mathbb{R}^m$ を $C^1$ 級とし、$x,y\in U$ を結ぶ線分 $\{x+t(y-x):0\le t\le1\}$ が $U$ に含まれるとします。この線分上のすべての点 $z$ で $\|DF(z)\|\le M$ ならば

$$
\|F(y)-F(x)\|\le M\|y-x\|
$$
です。
</Lemma>

<Proof of="lem-mean-value-ineq">
$d=y-x$、$v=F(y)-F(x)$ とおき、$v\neq0$ としてよい（$v=0$ なら自明）。スカラー関数 $\varphi(t)=\langle v, F(x+td)\rangle$（$0\le t\le1$）を考えます。連鎖律より $\varphi$ は微分可能で $\varphi'(t)=\langle v, DF(x+td)\,d\rangle$、そして $\varphi'$ は連続です。平均値の定理（<Ref to="mathematics/calculus/mean-value-and-taylor#thm-mvt" />）から、ある $\theta\in(0,1)$ で

$$
\|v\|^2=\langle v,F(y)-F(x)\rangle=\varphi(1)-\varphi(0)=\varphi'(\theta)=\langle v,DF(x+\theta d)\,d\rangle
$$
です。コーシー・シュワルツの不等式と上のノルム不等式から右辺は $\|v\|\,\|DF(x+\theta d)\|\,\|d\|\le \|v\|\,M\,\|d\|$ 以下です。$\|v\|>0$ で割れば $\|v\|\le M\|d\|$ を得ます。
</Proof>

## 3. 縮小写像の原理

方程式を解く道具として、次の原理を使います。1 次近似で解いた近似解を「補正して」何度も当てはめれば、真の解に収束する、という発想を定理にしたものです。

<Definition id="def-contraction" title="縮小写像">
$X\subset\mathbb{R}^n$ とします。写像 $T:X\to X$ が、ある定数 $q\in[0,1)$ について

$$
\|T(x)-T(x')\|\le q\,\|x-x'\| \qquad (x,x'\in X)
$$
を満たすとき、$T$ を**縮小写像**（縮小率 $q$）といいます。
</Definition>

<Theorem id="thm-contraction" title="縮小写像の原理（バナッハの不動点定理）">
$X\subset\mathbb{R}^n$ を空でない閉集合、$T:X\to X$ を縮小率 $q<1$ の縮小写像とします。このとき $T(x^*)=x^*$ となる点 $x^*\in X$ がただ 1 つ存在します。さらに任意の $x_0\in X$ から $x_{k+1}=T(x_k)$ で定めた列は $x^*$ に収束し、$\|x_k-x^*\|\le \dfrac{q^k}{1-q}\|x_1-x_0\|$ が成り立ちます。
</Theorem>

<Proof of="thm-contraction">
**列がコーシー列になること。** $\|x_{k+1}-x_k\|=\|T(x_k)-T(x_{k-1})\|\le q\|x_k-x_{k-1}\|$ を繰り返すと $\|x_{k+1}-x_k\|\le q^k\|x_1-x_0\|$ です。$l>k$ に対して三角不等式と等比級数の和から

$$
\|x_l-x_k\|\le\sum_{j=k}^{l-1}q^j\|x_1-x_0\|\le \frac{q^k}{1-q}\|x_1-x_0\|\xrightarrow[k\to\infty]{}0
$$
となります。したがって $(x_k)$ は[コーシー列](/mathematics/real-analysis/completeness-and-cauchy#def-cauchy)で、各成分は実数のコーシー列です。実数の完備性（<Ref to="mathematics/real-analysis/completeness-and-cauchy#thm-completeness" />）により各成分が収束し、$x_k\to x^*\in\mathbb{R}^n$ です。$X$ は閉集合なので $x^*\in X$ です。

**不動点であること。** $T$ は $\|T(x)-T(x')\|\le q\|x-x'\|$ より連続です。$x_{k+1}=T(x_k)$ の両辺で $k\to\infty$ とすれば $x^*=T(x^*)$ です。上の不等式で $l\to\infty$ とすれば誤差評価も得られます。

**一意性。** $x^*,x^{**}$ がともに不動点なら $\|x^*-x^{**}\|=\|T(x^*)-T(x^{**})\|\le q\|x^*-x^{**}\|$ で、$q<1$ より $\|x^*-x^{**}\|=0$ です。
</Proof>

注意すべきは、$q<1$ と「$T$ が $X$ を $X$ に写す」の 2 つが両方必要なことです。$T(x)=x/2$ は $X=(0,1]$ 上の縮小写像ですが不動点がなく、$X$ が閉でないために定理が使えません。逆関数定理の証明では、この 2 つの確認が議論の大半を占めます。

## 4. 逆関数定理

<Theorem id="thm-inverse-function" title="逆関数定理">
$U\subset\mathbb{R}^n$ を開集合、$F:U\to\mathbb{R}^n$ を $C^1$ 級写像とし、$a\in U$ で $DF(a)$ が正則（$\det DF(a)\ne0$）であるとします。このとき $a$ の開近傍 $V\subset U$ と $b=F(a)$ の開近傍 $W$ が存在して、$F|_V:V\to W$ は全単射であり、その逆写像 $G:W\to V$ も $C^1$ 級です。さらに

$$
DG(y)=\bigl(DF(G(y))\bigr)^{-1}\qquad (y\in W)
$$
が成り立ちます。
</Theorem>

方程式 $F(x)=y$ を $x=x-A^{-1}(F(x)-y)$ と書き換えるのが証明の鍵です。$A=DF(a)$ とおくと、右辺の $x$ の写像 $\Phi_y$ は $F=y$ の解のところで不動点を持ち、その微分は $I-A^{-1}DF(x)$ です。$x=a$ で $0$ ですから、$a$ の近くで小さく、縮小写像になります。

証明は 7 つのステップに分かれますが、役割は 3 つです。ステップ 1〜4 は「$a$ の近くで、$y$ ごとに解 $x$ がただ 1 つ見つかる」ことを示します。ステップ 1 は、微分が $A$ に近い範囲（閉球の半径 $r$）を決める作業です。$DF$ が連続なので、$a$ の近くでは $DF(x)$ が $A$ にいくらでも近く、条件を満たす $r$ がとれます。ステップ 3 の $s$ は、$\Phi_y(a)$ が $a$ から $r/2$ より動かないように $y$ の範囲を決めたもので、縮小の分の $r/2$ と合わせて、ちょうど $r$ 以内に収まります。ステップ 5 は、$y$ が少し動いても解 $x$ が大きくは動かないこと、ステップ 6・7 は、その逆写像が微分可能で、導関数が連続であることを示します。

<Proof of="thm-inverse-function">
$A=DF(a)$、$b=F(a)$ とおきます。

**ステップ 1: $r$ の選択。** $x\mapsto DF(x)$ は連続で $\det$ は成分の多項式ですから、十分小さな $r>0$ をとって閉球 $\bar B=\{x:\|x-a\|\le r\}\subset U$ とし、次の 2 条件を満たすようにできます。

- すべての $x\in\bar B$ で $\|A^{-1}\|\,\|DF(x)-A\|\le\tfrac12$。
- すべての $x\in\bar B$ で $\det DF(x)\ne0$。

**ステップ 2: 縮小写像の構成。** $y\in\mathbb{R}^n$ に対し $\Phi_y(x)=x-A^{-1}(F(x)-y)$ とおきます。$\Phi_y(x)=x$ は $F(x)=y$ と同値です。微分は $D\Phi_y(x)=I-A^{-1}DF(x)=A^{-1}(A-DF(x))$ なので、上のノルム不等式から $\|D\Phi_y(x)\|\le\frac12$（$x\in\bar B$）です。$\bar B$ は凸ですから、<Ref to="lem-mean-value-ineq" /> より

$$
\|\Phi_y(x)-\Phi_y(x')\|\le\tfrac12\|x-x'\|\qquad(x,x'\in\bar B)
$$
です。

**ステップ 3: $\Phi_y$ が $\bar B$ を $\bar B$ に写す。** $s=r/(2\|A^{-1}\|)$ とし、$\|y-b\|<s$ とします。$\Phi_y(a)=a-A^{-1}(b-y)$ ですから $\|\Phi_y(a)-a\|\le\|A^{-1}\|\|y-b\|<r/2$ です。$x\in\bar B$ に対し

$$
\|\Phi_y(x)-a\|\le\|\Phi_y(x)-\Phi_y(a)\|+\|\Phi_y(a)-a\|<\tfrac12 r+\tfrac12 r=r
$$
となります。ゆえに $\Phi_y(x)$ は開球 $B(a,r)$ に入り、とくに $\bar B$ に入ります。

**ステップ 4: 存在と一意性。** $\bar B$ は空でない閉集合ですから、縮小写像の原理（<Ref to="thm-contraction" />）により、$\|y-b\|<s$ を満たす各 $y$ に対して $F(x)=y$ の解 $x\in\bar B$ がただ 1 つ存在し、しかもステップ 3 から $x\in B(a,r)$ です。そこで $W=B(b,s)$、$V=B(a,r)\cap F^{-1}(W)$ とおきます。$F$ は連続ですから $V$ は開集合です。各 $y\in W$ に対し $\bar B$ 内の解 $x$ は $B(a,r)$ にあり $F(x)=y\in W$ ですから $x\in V$ です。よって $F|_V:V\to W$ は全単射で、逆写像 $G:W\to V$ が定まります。

**ステップ 5: $G$ はリプシッツ連続。** $x_1,x_2\in V$ に対し $y_i=F(x_i)$ とおくと $x_i=\Phi_{y_i}(x_i)$ ですが、同じ $y$ で比べるため $\Phi_0$ を使います（$\Phi_0(x)=x-A^{-1}F(x)$ で、ステップ 2 の評価は $y$ に依らない）。$x_1-x_2=\Phi_0(x_1)-\Phi_0(x_2)+A^{-1}(F(x_1)-F(x_2))$ なので

$$
\|x_1-x_2\|\le\tfrac12\|x_1-x_2\|+\|A^{-1}\|\,\|F(x_1)-F(x_2)\|,
$$
すなわち $\|x_1-x_2\|\le2\|A^{-1}\|\,\|F(x_1)-F(x_2)\|$ です。したがって $\|G(y_1)-G(y_2)\|\le 2\|A^{-1}\|\,\|y_1-y_2\|$ で、とくに $G$ は連続です。

**ステップ 6: $G$ は微分可能。** $y_0\in W$、$x_0=G(y_0)\in V$ とし、$B_0=DF(x_0)$ とおきます（ステップ 1 より正則）。$F$ の微分可能性から $F(x)=F(x_0)+B_0(x-x_0)+\varepsilon(x)$、$\|\varepsilon(x)\|/\|x-x_0\|\to0$（$x\to x_0$）です。$y\in W$、$x=G(y)$ に対し

$$
G(y)-G(y_0)-B_0^{-1}(y-y_0)=(x-x_0)-B_0^{-1}\bigl(F(x)-F(x_0)\bigr)=-B_0^{-1}\varepsilon(x)
$$
です。ステップ 5 より $\|x-x_0\|\le2\|A^{-1}\|\|y-y_0\|$、また $y\to y_0$ のとき $x\to x_0$（$G$ の連続性）ですから

$$
\frac{\|G(y)-G(y_0)-B_0^{-1}(y-y_0)\|}{\|y-y_0\|}\le\|B_0^{-1}\|\cdot2\|A^{-1}\|\cdot\frac{\|\varepsilon(x)\|}{\|x-x_0\|}\to0
$$
となります（$y\ne y_0$ なら $x\ne x_0$）。したがって $G$ は $y_0$ で微分可能で $DG(y_0)=B_0^{-1}=DF(G(y_0))^{-1}$ です。

**ステップ 7: $C^1$ 級。** 行列の逆行列は[余因子行列](/mathematics/linear-algebra/determinants#def-cofactor)と行列式で書けます（<Ref to="mathematics/linear-algebra/determinants#prop-adjugate" />）。すなわち $B^{-1}=(\det B)^{-1}\mathrm{adj}\,B$ で、成分は $B$ の成分の多項式を $\det B\ne0$ で割った形ですから、正則行列の集合上で連続です。$DG=(DF\circ G)^{-1}$ は連続関数 $G$、$DF$ と逆行列の合成ですから連続です。よって $G$ は $C^1$ 級です。
</Proof>

<Example id="ex-polar-inverse" title="極座標写像">
$F(r,\theta)=(r\cos\theta,\ r\sin\theta)$ は $\mathbb{R}^2$ 上の $C^\infty$ 級写像で

$$
DF(r,\theta)=\begin{pmatrix}\cos\theta&-r\sin\theta\\ \sin\theta&r\cos\theta\end{pmatrix},\qquad \det DF=r
$$
です。$r\ne0$ の各点で <Ref to="thm-inverse-function" /> が使え、$F$ は局所的に $C^1$ 級の逆写像を持ちます。しかし $F$ は大域的には単射でなく（$\theta$ に $2\pi$ を足しても同じ点）、$r=0$ では $\det DF=0$ で局所的にも単射ではありません（$\{0\}\times\mathbb{R}$ 全体が原点に潰れる）。**逆関数定理は局所的な主張であり、大域的な可逆性を保証しない**ことを、この例が示しています。なお、重積分の変数変換公式（<Ref to="mathematics/calculus/multiple-integrals#thm-change-of-variables" />）が仮定する「全単射で $\det J_\Phi\neq0$」を極座標が満たすのは、$r>0$、$0<\theta<2\pi$ のように定義域を制限したときです。
</Example>

<Remark id="rem-c1-necessary" title="$C^1$ 級という仮定を外すと何が壊れるか">
$f(x)=x+2x^2\sin(1/x)$（$x\ne0$）、$f(0)=0$ とします。$f$ は $\mathbb{R}$ 上で微分可能で、$f'(0)=\lim_{h\to0}(1+2h\sin(1/h))=1$ です。しかし $x\ne0$ では

$$
f'(x)=1+4x\sin\frac1x-2\cos\frac1x
$$
で、$x_k=1/(2\pi k)$ では $\sin(1/x_k)=0$、$\cos(1/x_k)=1$ より $f'(x_k)=-1<0$ です。$f'$ は $x\ne0$ で連続ですから $x_k$ の近くで $f$ は狭義減少です。$x_k\to0$ なので、$0$ のどんな近傍でも $f$ は単調でなく、単射でもありません。$f'(0)\ne0$ で 1 次近似は可逆なのに、局所逆写像が存在しない例です。ステップ 1 で $DF(x)$ が $a$ の近くで $A$ に近いことを使ったように、証明は $DF$ の連続性に依存しています。
</Remark>

## 5. 陰関数定理

<Theorem id="thm-implicit-function" title="陰関数定理">
$\Omega\subset\mathbb{R}^n\times\mathbb{R}^m$ を開集合とし、点を $(x,y)$（$x\in\mathbb{R}^n,\ y\in\mathbb{R}^m$）と書きます。$F:\Omega\to\mathbb{R}^m$ を $C^1$ 級、$(a,b)\in\Omega$ で $F(a,b)=0$ とし、$y$ に関する偏微分のなす $m\times m$ 行列 $D_yF(a,b)=\bigl(\partial F_i/\partial y_j(a,b)\bigr)$ が正則であるとします。このとき $a$ の開近傍 $X$、$b$ の開近傍 $Y$ と $C^1$ 級写像 $g:X\to Y$ が存在して、$X\times Y\subset\Omega$ であり

$$
\{(x,y)\in X\times Y:F(x,y)=0\}=\{(x,g(x)):x\in X\}
$$
が成り立ちます。さらに $g(a)=b$ であり、$D_xF$ を $x$ に関する偏微分のなす $m\times n$ 行列として

$$
Dg(x)=-\bigl(D_yF(x,g(x))\bigr)^{-1}D_xF(x,g(x))
$$
です。
</Theorem>

証明の方針は、$x$ をそのまま残して $F$ の値だけを新しい座標にする写像 $\Psi(x,y)=(x,F(x,y))$ を作ることです。$\Psi$ は入力と出力の次元が等しいので逆関数定理が使え、その逆写像を使って「$F=0$ となる $y$」を $x$ の関数として読み取れます。

<Proof of="thm-implicit-function">
$\Psi:\Omega\to\mathbb{R}^n\times\mathbb{R}^m$ を $\Psi(x,y)=(x,F(x,y))$ と定めます。$\Psi$ は $C^1$ 級で、ヤコビ行列はブロック行列

$$
D\Psi(x,y)=\begin{pmatrix}I_n&0\\ D_xF&D_yF\end{pmatrix}
$$
です。ブロック三角行列の行列式は対角ブロックの行列式の積ですから $\det D\Psi=\det D_yF$ で、$(a,b)$ で $\ne0$ です。逆関数定理（<Ref to="thm-inverse-function" />）により、$(a,b)$ の開近傍 $V\subset\Omega$ と $\Psi(a,b)=(a,0)$ の開近傍 $W$ が存在して、$\Psi:V\to W$ は全単射で、逆写像 $\Psi^{-1}$ は $C^1$ 級です。

$\Psi$ は第 1 成分を変えないので、$\Psi^{-1}(x,z)=(x,H(x,z))$ の形で、$H$ は $C^1$ 級です。$W$ は $(a,0)$ の近傍ですから、$a$ の開近傍 $X$ と $0$ の開近傍 $Z$ で $X\times Z\subset W$ となるものがあります。$g(x)=H(x,0)$（$x\in X$）とおくと、$\Psi(x,g(x))=(x,0)$ すなわち $F(x,g(x))=0$ で、$g$ は $C^1$ 級、$g(a)=b$ です。

一意性を示します。$V$ の点 $(x,y)$ で $x\in X$、$F(x,y)=0$ なら $\Psi(x,y)=(x,0)=\Psi(x,g(x))$ ですから、$\Psi$ の単射性から $y=g(x)$ です。そこで $Y$ を $b$ の開近傍で、$X$ を必要なら縮めて $X\times Y\subset V$ かつ $g(X)\subset Y$ となるように選べば（$g$ の連続性と $V$ の開性より可能）、主張の集合の等式が成り立ちます。

最後に導関数です。恒等式 $F(x,g(x))=0$ の両辺を連鎖律（<Ref to="mathematics/calculus/multivariable-differentiation#thm-chain-rule-n" />）で微分して $D_xF+D_yF\,Dg=0$ を得ます。$D_yF(x,g(x))$ が正則であることは、$X$ を縮めて $\det D_yF\neq0$ となるようにとってあれば成り立つので、両辺に逆行列を掛けて公式を得ます。
</Proof>

定理の内容を円で図にすると次のようになります。点 $(a,b)$ のまわりを長方形 $X\times Y$ で切り取ると、その中では曲線 $F=0$ が $x$ の関数 $y=g(x)$ のグラフだけになります。

<Figure caption="円 x²+y²=1 の上の点 (a,b)=(0,1) のまわりを長方形 X×Y で切り取る。長方形の中では曲線が y=g(x) のグラフ（1 本の曲線）になる">

<svg viewBox="0 0 400 260" width="100%" role="img" aria-label="円の上端の点のまわりを長方形で囲んだ図。長方形の左右の辺を横切る円弧が、長方形の中で x の関数のグラフになっている">
  <circle cx="160" cy="190" r="100" fill="none" stroke="currentColor" stroke-width="1.2" opacity="0.5" />
  <rect x="110" y="70" width="100" height="50" fill="none" stroke="currentColor" stroke-width="1.2" stroke-dasharray="5 4" />
  <path d="M110,103.4 A100,100 0 0 1 210,103.4" fill="none" stroke="var(--sl-color-accent)" stroke-width="3" />
  <circle cx="160" cy="90" r="4.5" fill="var(--sl-color-accent)" />
  <g fill="currentColor" font-size="13">
    <text x="146" y="62">(a, b)</text>
    <text x="156" y="138">X</text>
    <text x="218" y="100">Y</text>
    <text x="216" y="118">y = g(x)</text>
  </g>
</svg>

</Figure>

<Example id="ex-implicit-circle" title="円と陰関数定理">
$F(x,y)=x^2+y^2-1$ では $D_xF=2x$、$D_yF=2y$ です。$(a,b)=(0,1)$ では $D_yF=2\ne0$ なので、<Ref to="thm-implicit-function" /> により $y=g(x)$ と解けます。実際 $g(x)=\sqrt{1-x^2}$ で、公式は $g'(x)=-2x/(2g(x))=-x/\sqrt{1-x^2}$ を与え、直接の微分と一致します。$(a,b)=(1,0)$ では $D_yF=0$ で定理は使えず、実際 $x=1$ の近くで $y$ は一意に決まりません。ただし $D_xF=2\neq 0$ ですから、役割を入れ替えて $x=\sqrt{1-y^2}$ と**$x$ を $y$ の関数として**解くことはできます。条件は「どの変数について解くか」ごとに判定します。
</Example>

<Example id="ex-implicit-folium" title="デカルトの葉線">
$F(x,y)=x^3+y^3-3xy$ とすると $D_xF=3x^2-3y$、$D_yF=3y^2-3x$ です。点 $(3/2,3/2)$ は $F=0$ を満たし（$27/8\cdot 2-27/4=0$）、$D_yF=9/4\ne0$ ですから、その近くで $y=g(x)$ と解けます。$g'(3/2)=-D_xF/D_yF=-(9/4)/(9/4)=-1$ です。原点では $D_xF=D_yF=0$ で、どちらの変数についても解けません。実際、葉線は原点で自分自身と交わる点（結節点）を持ちます。
</Example>

## 6. 条件付き極値とラグランジュの乗数法

関数 $f$ の極値を、制約 $g(x)=0$ を満たす点の中でだけ調べる問題を考えます。たとえば「周の長さが一定の長方形のうち面積が最大のもの」です。制約が陰関数定理で解けるなら、$m$ 個の変数を消去して、変数が $n-m$ 個の通常の極値問題に帰着できます。ただし解くのが面倒なとき、解を陽に求めずに同じ結論を出せるのがラグランジュの乗数法です。

<Definition id="def-constrained-extremum" title="条件付き極値">
$U\subset\mathbb{R}^n$ を開集合、$f:U\to\mathbb{R}$、$g=(g_1,\dots,g_m):U\to\mathbb{R}^m$（$m<n$）とし、$M=\{x\in U:g(x)=0\}$ とおく。点 $a\in M$ が $f|_M$ の**極大点**であるとは、$a$ のある近傍 $N$ が存在して $f(x)\le f(a)$（$x\in M\cap N$）が成り立つことをいう。極小点も同様で、あわせて**条件付き極値点**（制約 $g=0$ の下での極値点）と呼ぶ。
</Definition>

<Theorem id="thm-lagrange" title="ラグランジュの乗数法">
$U\subset\mathbb{R}^n$ を開集合、$f:U\to\mathbb{R}$ と $g=(g_1,\dots,g_m):U\to\mathbb{R}^m$（$m<n$）を $C^1$ 級とします。$a\in M=\{g=0\}$ が $f|_M$ の条件付き極値点であり、かつ $Dg(a)$ の階数が $m$（すなわち $\nabla g_1(a),\dots,\nabla g_m(a)$ が一次独立）ならば、実数 $\lambda_1,\dots,\lambda_m$ が存在して

$$
\nabla f(a)=\sum_{i=1}^m\lambda_i\nabla g_i(a)
$$
が成り立ちます。この $\lambda_i$ を**ラグランジュ乗数**といいます。
</Theorem>

<Proof of="thm-lagrange">
**座標の入れ替え。** $Dg(a)$ は階数 $m$ の $m\times n$ 行列ですから、一次独立な $m$ 個の列があり、それを取り出した $m\times m$ 行列は正則です（<Ref to="mathematics/linear-algebra/determinants#thm-regular" />）。座標の番号を付け替えて、これが最後の $m$ 個の列であるとしてよい。$x=(x',y)$（$x'\in\mathbb{R}^{n-m}$、$y\in\mathbb{R}^m$）、$a=(a',b)$ と書くと、$D_yg(a)$ が正則です。

**制約を解く。** 陰関数定理（<Ref to="thm-implicit-function" />）により、$a'$ の開近傍 $X$、$b$ の開近傍 $Y$ と $C^1$ 級写像 $\varphi:X\to Y$ が存在し、$X\times Y$ 内では $M=\{(x',\varphi(x'))\}$ で、$\varphi(a')=b$、$D\varphi=-(D_yg)^{-1}D_{x'}g$ です。

**自由な $n-m$ 変数の通常の極値問題に帰着。** $h(x')=f(x',\varphi(x'))$（$x'\in X$）とおきます。$(x',\varphi(x'))\in M$ ですから、$a$ が $f|_M$ の極値点であることから、$a'$ は $h$ の（通常の）極値点です。$h$ は $C^1$ 級で、各 $j$ について $t\mapsto h(a'+te_j)$ は $t=0$ で極値をとる微分可能な関数ですから、フェルマーの補題（<Ref to="mathematics/calculus/mean-value-and-taylor#lem-fermat" />）により $\partial h/\partial x'_j(a')=0$ です。連鎖律で書けば

$$
0=\nabla_{x'}h(a')=\nabla_{x'}f(a)+D\varphi(a')^{\mathsf T}\nabla_yf(a)=\nabla_{x'}f(a)-D_{x'}g(a)^{\mathsf T}\bigl(D_yg(a)^{\mathsf T}\bigr)^{-1}\nabla_yf(a).
$$

**乗数の定義。** $\lambda=(D_yg(a)^{\mathsf T})^{-1}\nabla_yf(a)\in\mathbb{R}^m$ とおきます。定義より $\nabla_yf(a)=D_yg(a)^{\mathsf T}\lambda$ で、上の式から $\nabla_{x'}f(a)=D_{x'}g(a)^{\mathsf T}\lambda$ です。2 つをまとめると $\nabla f(a)=Dg(a)^{\mathsf T}\lambda=\sum_i\lambda_i\nabla g_i(a)$ です。
</Proof>

幾何的な意味を確認します。$m=1$ なら、$\nabla f(a)$ と $\nabla g(a)$ が平行、すなわち $f$ の等位線が制約曲線 $g=0$ に接する、という条件です。制約曲線に沿って動いて $f$ が増減するなら、その点は極値ではありません。曲線に沿った増減が 0 になるのは、$\nabla f$ が曲線の接線方向と直交する（＝法線方向 $\nabla g$ と平行になる）ときです。

<Figure caption="単位円 x²+y²=1 上で f=x+y を最大化する。等位線 x+y=c が円に接する点で f は極値をとり、そこで ∇f と ∇g が平行になる">

<svg viewBox="0 0 400 300" width="100%" role="img" aria-label="単位円と、f=x+y の等位線である斜めの直線。直線が円に接する右上の点で、f の勾配ベクトルと円の法線ベクトルが一致している図">
  <g stroke="currentColor" fill="none" stroke-width="1" opacity="0.5">
    <path d="M30,150 H320" />
    <path d="M170,285 V15" />
  </g>
  <circle cx="170" cy="150" r="100" fill="none" stroke="var(--sl-color-accent)" stroke-width="2" />
  <g stroke="currentColor" fill="none" stroke-width="1.2" stroke-dasharray="5 4">
    <path d="M190,28.6 L290,128.6" />
    <path d="M90,70 L250,230" />
    <path d="M50,171.4 L150,271.4" />
  </g>
  <circle cx="240.7" cy="79.3" r="4.5" fill="var(--sl-color-accent)" />
  <circle cx="99.3" cy="220.7" r="4.5" fill="var(--sl-color-accent)" />
  <g stroke="currentColor" fill="none" stroke-width="2">
    <path d="M240.7,79.3 L270,50" />
  </g>
  <path d="M276,44 L262,48 L272,58 Z" fill="currentColor" />
  <g fill="currentColor" font-size="13">
    <text x="276" y="40">∇f ∥ ∇g</text>
    <text x="294" y="140">f = √2</text>
    <text x="256" y="236">f = 0</text>
    <text x="20" y="190">f = −√2</text>
    <text x="316" y="168">x</text>
    <text x="176" y="20">y</text>
  </g>
</svg>

</Figure>

<Example id="ex-lagrange-circle" title="円周上の $x+y$ の最大・最小">
$f=x+y$、$g=x^2+y^2-1$ とします。$M$ は円周で[コンパクト](/mathematics/real-analysis/metric-spaces#def-compactness-metric)ですから、連続関数 $f$ は $M$ 上で最大値と最小値をとります（<Ref to="mathematics/topology/compactness#cor-extreme-value" />）。それらの点では $\nabla g=(2x,2y)\ne0$ ですから <Ref to="thm-lagrange" /> が使え、$(1,1)=\lambda(2x,2y)$、つまり $x=y=1/(2\lambda)$ です。$g=0$ に代入して $x=y=\pm1/\sqrt2$。値は $f=\pm\sqrt2$ で、最大値は $\sqrt2$、最小値は $-\sqrt2$ です。
</Example>

<Example id="ex-lagrange-amgm" title="3 変数の相加相乗平均の不等式">
$x,y,z\ge0$、$x+y+z=3$ のもとで $f=xyz$ の最大値を求めます。この集合は有界閉集合（三角形）ですから最大値が存在します。境界（$x$、$y$、$z$ のどれかが $0$）では $f=0$ ですから、最大点は $x,y,z>0$ の内部にあります。内部で制約 $g=x+y+z-3$ について $\nabla g=(1,1,1)\neq0$ ですから <Ref to="thm-lagrange" /> が使え、$(yz,xz,xy)=\lambda(1,1,1)$ です。$x,y,z>0$ より $xyz=\lambda x=\lambda y=\lambda z$ で、$\lambda\neq 0$ なので $x=y=z=1$、最大値は $1$ です。よって $x+y+z=3$ のとき $xyz\le1$ で、任意の $u,v,w>0$ に対し $x=3u/s$ などと正規化すれば $uvw\le((u+v+w)/3)^3$ を得ます。
</Example>

<Remark id="rem-lagrange-rank" title="階数の仮定が必要な理由">
$f(x,y)=x$、$g(x,y)=y^2-x^3$ とします。$M=\{y^2=x^3\}$ は尖点をもつ曲線で、$M$ 上では $x=y^{2/3}\ge0$ ですから、$f|_M$ は原点で最小値 $0$ をとります。しかし $\nabla g(0,0)=(0,0)$ で、$\nabla f=(1,0)$ はそのスカラー倍になりません。階数条件 $\mathrm{rank}\,Dg(a)=m$ が落ちる点では、<Ref to="thm-implicit-function" /> が使えず、制約集合 $M$ が曲線として解けないため、結論が成り立たなくなります。ラグランジュの乗数法を使うときは、候補点に加えて**$\nabla g=0$ となる点**も別に調べる必要があります。
</Remark>

<Aside type="caution">
乗数法は極値点の**必要条件**にすぎません。$\nabla f=\lambda\nabla g$ を満たす点がすべて極値点とは限りません。極値の存在がコンパクト性などで保証されている場合に、候補を絞って比較するのが典型的な使い方です。
</Aside>

## 7. 演習

<Exercise id="exr-inverse-jacobian" difficulty="標準">
$F(x,y)=(e^x\cos y,\ e^x\sin y)$ について、すべての点で $DF$ が正則であることを示し、点 $(0,0)$ の近傍での逆写像 $G=F^{-1}$ のヤコビ行列 $DG(1,0)$ を求めよ。また、$F$ が $\mathbb{R}^2$ 全体では単射でないことを示せ。
<Solution>
$DF=\begin{pmatrix}e^x\cos y&-e^x\sin y\\ e^x\sin y&e^x\cos y\end{pmatrix}$、$\det DF=e^{2x}\ne0$ です。$F(0,0)=(1,0)$ で $DF(0,0)=I$ ですから、<Ref to="thm-inverse-function" /> により $(0,0)$ の近傍で逆写像 $G$ が存在し、$DG(1,0)=DF(0,0)^{-1}=I$ です。また $F(x,y+2\pi)=F(x,y)$ ですから $F$ は単射ではありません。逆関数定理は局所的な主張だと再確認できます。
</Solution>
</Exercise>

<Exercise id="exr-implicit-derivative" difficulty="標準">
方程式 $x^2+xy+y^3=3$ は点 $(1,1)$ を通る。この点の近くで $y$ を $x$ の関数 $y=g(x)$ として解けることを示し、$g'(1)$ と $g''(1)$ を求めよ。
<Solution>
$F=x^2+xy+y^3-3$ とすると $F(1,1)=0$、$F_y=x+3y^2=4\ne0$ ですから <Ref to="thm-implicit-function" /> により解けます。$F_x=2x+y=3$ で $g'(1)=-3/4$ です。恒等式 $x^2+xg+g^3=3$ を $x$ で 2 回微分すると、1 回目は $2x+g+xg'+3g^2g'=0$、2 回目は $2+2g'+xg''+6g(g')^2+3g^2g''=0$ です。$x=1$、$g=1$、$g'=-3/4$ を代入して $2-3/2+g''+6\cdot9/16+3g''=0$、すなわち $4g''=-(1/2+27/8)=-31/8$ より $g''(1)=-31/32$ です。
</Solution>
</Exercise>

<Exercise id="exr-lagrange-distance" difficulty="標準">
平面 $x+2y+3z=14$ 上の点のうち、原点に最も近い点とその距離を、ラグランジュの乗数法で求めよ。
<Solution>
$f=x^2+y^2+z^2$、$g=x+2y+3z-14$ とおきます。$\nabla g=(1,2,3)\ne0$ で、$\nabla f=\lambda\nabla g$ から $(2x,2y,2z)=\lambda(1,2,3)$、つまり $(x,y,z)=\tfrac{\lambda}{2}(1,2,3)$ です。$g=0$ に代入して $\tfrac\lambda2\cdot14=14$ より $\lambda=2$、点は $(1,2,3)$ で、$f=14$、距離は $\sqrt{14}$ です。最小値の存在は、$f$ が平面上で $\|x\|\to\infty$ のとき $\infty$ に発散するため、十分大きな球と平面の共通部分（コンパクト）で最小をとることから分かります。したがって候補が 1 つだけなので、それが最小点です。
</Solution>
</Exercise>

<Exercise id="exr-contraction-fail" difficulty="易">
次の写像が縮小写像の原理（<Ref to="thm-contraction" />）の仮定のうちどれを満たさないか述べ、不動点の有無を答えよ。(1) $X=\mathbb{R}$、$T(x)=x+1$。(2) $X=(0,1]$、$T(x)=x/2$。(3) $X=[0,\infty)$、$T(x)=\sqrt{x^2+1}$。
<Solution>
(1) $|T(x)-T(x')|=|x-x'|$ で縮小率が $1$（$q<1$ を満たさない）。不動点なし。(2) $X$ が閉集合でない（$0\notin X$）。不動点なし（候補 $0$ は $X$ の外）。(3) $T'(x)=x/\sqrt{x^2+1}<1$ で、平均値の定理から $|T(x)-T(x')|<|x-x'|$ ですが、$\sup T'=1$ で一様な $q<1$ は取れません。不動点の方程式 $x=\sqrt{x^2+1}$ は解を持ちません。つまり「各点で縮小」では不十分で、定数 $q<1$ が要ります。
</Solution>
</Exercise>

## 参考文献

- 杉浦光夫『解析入門 I』東京大学出版会、1980。
- 杉浦光夫『解析入門 II』東京大学出版会、1985。
- 松坂和夫『解析入門 中』岩波書店、1968。
- W. Rudin, *Principles of Mathematical Analysis*, 3rd ed., McGraw-Hill, 1976 — Chapter 9（逆関数定理と陰関数定理）。
