物理/多変数解析学
提供: Internet Web School
(→微分(全微分) ) |
(→微分(全微分) ) |
||
282 行: | 282 行: | ||
${\bf z} = g(\bf{y})$ を$\bf{R^m}$ の開集合$V$ から${\bf R^n}$への関数とする。<br/> | ${\bf z} = g(\bf{y})$ を$\bf{R^m}$ の開集合$V$ から${\bf R^n}$への関数とする。<br/> | ||
もし関数$f$が点$\bf{x^0}(\in U)$で微分可能で、<br/> | もし関数$f$が点$\bf{x^0}(\in U)$で微分可能で、<br/> | ||
- | $\bf | + | ${\bf y^0}\triangleq f({\bf x^0}) \in V$であり<br/> |
- | 関数$g$が点$\bf | + | 関数$g$が点${\bf y^0}$で微分可能であるならば<br/> |
合成関数<br/> | 合成関数<br/> | ||
- | $\bf | + | ${\bf z} = g\circ f({\bf x})\triangleq g(f({\bf x}))$<br/> |
- | は、点$\bf | + | は、点${\bf x^0}$で微分可能で<br/> |
- | その点の導値 $D(g\circ f)(\bf | + | その点の導値 $D(g\circ f)({\bf x^0})$は<br/> |
- | $D(g\circ f)(\bf | + | $D(g\circ f)({\bf x^0}) = Dg({\bf y^0})Df({\bf x^0}) \quad (注参照のこと)\qquad (a)$<br/> |
である。<br/> | である。<br/> | ||
- | (注)右辺はn×m行列$Dg(\bf | + | (注)右辺はn×m行列$Dg({\bf y^0})$ とm×l行列$Df({\bf x^0})$の行列としての積である。<br/> |
証明<br/> | 証明<br/> | ||
- | 関数$f$が点$\bf | + | 関数$f$が点${\bf x^0}$で微分可能なので、微分可能の定義から<br/> |
- | ノルムの十分小さい任意のl次元ベクトル$\bf | + | ノルムの十分小さい任意のl次元ベクトル${\bf h}$に対して、<br/> |
- | $f(\bf | + | $f({\bf x^0}+{\bf h}) = f({\bf x^0}) + Df({\bf x^0}){\bf h} + \delta(\bf{h})\qquad (a)$<br/> |
- | ここで、$\lim_{\bf | + | ここで、$\lim_{{\bf h}\to {\bf 0}}\frac{\delta({\bf h})}{\|{\bf h} \|}= 0\qquad (b)$<br/> |
- | 同様に、ノルムの十分小さい任意のm次元ベクトル$\bf | + | 同様に、ノルムの十分小さい任意のm次元ベクトル${\bf k}$に対して、<br/> |
- | $g(\bf | + | $g({\bf y^0}+{\bf k}) = g({\bf y^0}) + Dg({\bf y^0}){\bf k} + \delta'(\bf{k})\qquad (c)$<br/> |
- | ここで、$\lim_{\bf | + | ここで、$\lim_{{\bf k}\to {\bf 0}}\frac{\delta'({\bf k})}{\|{\bf k} \|}= 0\qquad (d)$<br/> |
- | + | 式(a)から、<br/> | |
- | $g\circ f(\bf | + | $g\circ f({\bf x^0}+{\bf h}) = g\Bigl(f({\bf x^0}+{\bf h})\Bigr)$<br/> |
- | $ | + | $=g\Bigl(f({\bf x^0}) + Df({\bf x^0}){\bf h} + \delta({\bf h})\Bigr)$<br/> |
- | + | $\qquad$ そこで、${\bf k(h)}=Df({\bf x^0}){\bf h} + \delta({\bf h}) $ とおくと<br/> | |
- | $\bf | + | $=g\Bigl(f({\bf x^0}) + {\bf k(h)}\Bigr)$ |
- | 式(c | + | ${\bf h}$が零ベクトル近づくとき${\bf k}$も零ベクトルに近づくので、<br/> |
- | $g\circ f(\bf | + | 式(c)を適用できて<br/> |
- | $\ | + | $g\circ f({\bf x^0}+{\bf h})$ <br/> |
+ | $=g({\bf y^0}) + Dg({\bf y^0}){\bf k(h)} + \delta'({\bf k(h)})\qquad (c)$<br/> $<br/> | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | $=g\Bigl(f({\bf x^0})\Bigr) + Dg({\bf y^0})\Bigl(Df({\bf x^0}){\bf h} + \delta({\bf h})\Bigr) + \delta'\Bigl(Df(\bf{x^0})\bf{h} + \delta(\bf{h}) \Bigr)\qquad (e)$<br/> | ||
+ | $= g\circ f({\bf x^0})+ Dg({\bf y^0})Df({\bf x^0}){\bf h}+Dg({\bf y^0})\delta({\bf h}) + \delta'({\bf k(h)}) \Bigr) | ||
+ | $\lim_{{\bf h}\to {\bf 0}}\frac{\delta'({\bf k(h)})}{\|{\bf k(h)} \|}= 0$<br/> | ||
=== 高階偏微分 === | === 高階偏微分 === |
2017年12月15日 (金) 17:11時点における版
目次 |
「9.1 多変数解析学」
序
本章の冒頭の偏微分の導入部については下記の本も参考にしてください。
それ以降の内容については、ウィキブックスには殆どないため、
このテクストで今後叙述していく予定です。
多変数の実数値関数の微分
${\bf R^n}=\{(x_1,x_2,,,x_n) \mid x_i\in{\bf R},i=1,2,\cdots n\}$ の開区間
$I^n=\prod_{i=1}^{n}(a_i,b_i)$上で定義された実関数 $y=f(x_1,x_2,,,x_n)$ を考える。
一変数関数の議論から類推するために
以後、$\textbf x:=(x_1,x_2,,,x_n)$とおき、 $y=f(\textbf x)$ と書くこともある。
$I^n \,$上で定義された実数値関数 $\ y=f(\textbf x)=f(x_1,x_2,,,x_n)\,$ の微分について説明する。
一変数の微分から類推すると
微小なベクトル $\textbf h=(h_1,h_2,,,h_n)$ を考え、極限
$\lim_{\vec h \to 0,\vec h\neq 0}\frac{f(\textbf x + \vec h)-f(\textbf x)}{{\bf h} }$
が存在するとき、関数fは微分可能と定義することが考えられる。
しかし残念ながら、
$\vec h$はn次元ベクトルなので、割り算は不可能でありこの定義は無効である。
偏微分
関数$f$ の変数 $\textbf x$ の第i成分 $x_i$ だけを変数とし、
他の変数は任意の実数に固定$\Bigl(x_j = a_j \quad (j\neq i)\Bigr)$して得られる関数
$\phi_{x_j=a_j,j\neq i}(x_i)\triangleq f(a_1,a_2,,,a_{i-1},x_i,a_{i+1},,,a_n) $
を考える。
この関数は、一変数なので、任意の点$x_i $ での微分係数
$\frac{d\phi_{x_j=a_j,j\neq i}}{dx_i}(x_i)\triangleq \lim_{ h \to 0, h\neq 0}\frac{\phi_{x_j=a_j,j\neq i}(x_i+h)-\phi_{x_j=a_j,j\neq i}(x_i)}{h}$
$=\lim_{ h \to 0, h\neq 0}\frac{ f(a_1,a_2,,,a_{i-1},x_{i}+h,a_{i+1},,,a_n)-f(a_1,a_2,,,a_{i-1},x_{i},a_{i+1},,,a_n)}{h}$
を考えることができる。
定義1(偏微分)
もし、一変数関数 $\phi_{x_j=a_j,j\neq i}(x_i)=f(a_1,a_2,,,a_{i-1},x_i,a_{i+1},,,a_n)$ が、ある点$x_i=a_i$で微分可能ならば、
関数$f$は、点$\vec a = (a_1.a_2,,,,a_n)$で,$x_i$ について偏微分可能であると言い,
$\frac{\partial f}{\partial x_i}(\vec a) \triangleq \frac{d\phi_{x_j=a_j,j\neq i}}{dx_i}(a_i)$
を、$f(\textbf x)$ の 点$\vec a$ での変数 $x_i$ についての偏微分係数という。
定義2(偏導関数)
$f(\textbf x)$ がどの点$\textbf x$でも $x_i$ に関して偏微分可能であるならば、
任意の点$\textbf x$ にその点における $x_i$ に関する偏微分係数$\frac{\partial f}{\partial x_i}(\textbf x)$を対応させると、新しい関数が得られる。
これを、$f(\textbf x)$ の $x_i$ に関する偏導関数といい、記号
$f_{x_{i}}(\textbf x),\quad D_{x_i}f(\textbf x),\quad \frac{\partial f}{\partial x_i} (\textbf x),\quad \partial f/\partial x_i$
などで表示する。
以後、簡単のために2変数$x, y$ の関数に限定して議論する。
定理1 合成関数の微分(1)
$R^2$ から $R$ への関数$f(x,y)$ と
$R$ から $R$ への関数$g(t)$ の合成関数
$h(x,y)=g(f(x,y))$
を考える。
もし、$f(x,y)$ が $(x_0,y_0)$ で、xに関して偏微分可能で,
$\qquad g(t)$ が、$t_0=f(x_0,y_0)$ において微分可能ならば、
$h(x,y)=g(f(x,y))$ は $(x_0,y_0)$ で、xに関して偏微分可能であり,
$h_{x}(x_0,y_0)=g'(t_0)f_{x}(x_0,y_0) \qquad \qquad (1)$
証明
yを $ y_0 $ に固定して考えると、一変数関数の合成関数の微分になるので、合成関数の微分公式を適用すればよい。
定理2
$f(x,y)$ を
$(x_0,y_0)$ を中心とするある半径rの開球体$B_{r}(x_0,y_0)$上で、$x$について偏微分可能とする。
もし$(x,y)$ を$B_{r}(x_0,y_0)$の点ならば
$x_0$ と $x$ の間の $\xi$ が存在して、
$f(x,y) - f(x,y_0)-\bigl(f(x_0,y) - f(x_0,y_0) \bigr)= (x-x_0)\bigl(f_{x}(\xi,y)-f_{x}(\xi,y_0)\bigr) \qquad (2)$
(注)2次元の開球体$B_{r}(x_0,y_0)\triangleq \{(x,y)\in R^2 | \|(x,y) - (x_0,y_0) \|_{2} \lt r\} $ は、中心が点$(x_0,y_0)$ で半径rの円周で囲まれる内部である。
証明
$ \phi(x)\triangleq f(x,y) - f(x,y_0)$ とおくと、
式()の左辺$ = \phi(x) - \phi(x_0)$
$\quad \phi(x) $ は、$x_0$ の近傍で微分可能なので、平均値の定理から、
$\quad x_0$ と $x$ の間の $\xi$ が存在して、
$= (x-x_0){\phi}'(\xi) = (x-x_0)\bigl(f_{x}(\xi,y)-f_{x}(\xi,y_0)\bigr) $
定理3
$f(x.y)$ を
$(x_0,y_0)$ を中心とする開球体$B_{r}(x_0,y_0)$上で、xについて偏微分可能とする。
もし$(x,y)=(x_0+h,y_0+k)\in B_{r}(x_0,y_0)$ ならば
$f(x,y) = f(x_0,y_0) + hf_{x}(x_0 + h\theta,y) + kf_{y}(x_0,y_0+ k\theta)$
を満たす、$\theta=\theta(h,k) \in (0,1)$ が存在する。
証明
$ g(t) \triangleq f(x_0+ht,y) + f(x_0,y_0+kt) $ というtの関数を導入する。
すると、
$g(1)-g(0) = f(x,y)+f(x_0,y)-\bigl(f(x_0,y)+f(x_0,y_0) \bigr)$
$\qquad \qquad =f(x,y) - f(x_0,y_0)$
関数 $g(t)$ は、閉区間[0,1] を含む開区間上で微分可能なので、
一変数の微分可能関数の平均値の定理から、
ある数 $\theta \in (0,1)$ が存在して、
$g(1)-g(0) = g'(\theta)(1-0) = g'(\theta) \qquad \qquad (a)$
故に、$ f(x,y) - f(x_0,y_0) = g(1)-g(0) = g'(\theta)$
$\qquad $ 関数gの微分は,一変数関数の合成関数の微分公式から
$\qquad g'(t) = f_{x}(x_0+ht,y)h + f_{y}(x_0,y_0+kt)k \qquad (b)$
式(a)、(b) から
$ f(x,y) - f(x_0,y_0) = f_{x}(x_0+h\theta,y)h + f_{y}(x_0,y_0+k\theta)k \qquad (b)$
証明終わり
方向微分
$\vec{e_i}$ を直交座標系の$x_i$座標軸の正方向の方向・向きを持つ単位長さのベクトルとする(第i直交座標ベクトルと呼ぼう)。
多変数関数$y=f(x_1,x_2,,,x_n)$の、点$\textbf x = (x_1,x_2,,,x_n)$での偏微分係数 $\frac{\partial f}{\partial x_i}(x)$ は、
点$\textbf x $ を、第i座標(座標ベクトル$\vec{e}_i$)に平行に無限に小さい距離移動させるときの、関数fの変化率とみなせる。
式で書くと
$\frac{\partial f}{\partial x_i}(x)
= \lim_{h\to 0,h\neq 0}\frac{f(\textbf x + h\vec{e}_i)-f(\textbf x )}{h}$
このように考えると、点$\textbf x = (x_1,x_2,,,x_n)$を、座標ベクトル$\vec{e}_i$に平行ではなく、
任意に指定するベクトル$\vec a$に平行に微小量動かすときの関数fの変化率を考えることもできることが分かるだろう。
定義 方向微分
関数$y=f(x_1,x_2,,,x_n)$の、点$\textbf x = (x_1,x_2,,,x_n)$での,$\vec a$ 方向の微分係数とは、
$\lim_{h\to 0,h\neq 0}\frac{f(\textbf x + h\vec a)-f(\textbf x )}{h}$
のことで、
$\frac{\partial f}{\partial \vec{a}}(x),\quad f_{\vec a}(x),\quad D_{\vec a}f(x)$
などと書く。
命題1
(1) $\vec{e_i}$ 方向の微分は、$\vec{e_i}$ 座標軸($x_i$座標軸)に関する偏微分である。
ここで、$\vec{e_i}$ は$x_i$座標軸の正方向向きの単位長さのベクトル。
式で書くと、
$\frac{\partial f}{\partial \vec{e_i}}(x) = \frac{\partial f}{\partial x_i}(x) $
(2)$\alpha$ を任意の実数とすると
$\frac{\partial f}{\partial \alpha \vec{e_i}}(x) = \alpha \frac{\partial f}{\partial x_i}(x) $
微分(全微分)
この§も、記述を簡単にするため、2変数関数で説明する。
一般のn変数の場合への拡張は、記述は複雑になるが、容易である。
多変数実数値関数の微分可能性
二変数関数の微分可能性をどう定義したらよいだろうか?
一変数関数の微分の場合、それと同等の条件はいくつか知られているが、
その中で二変数関数に容易に拡張できるものを採用するのが自然である。
1.4.1.1 微分係数の意味 の命題の条件 3)の式(5)が、それに該当する。
定義3 微分可能性(全微分可能性)
関数$f(x,y)$が、或る開集合U(\subset {\bf R^2})上で定義されているとする。
$f$が 点$(x_0,y_0)\in U$ で微分可能(あるいは全微分可能)とは、
ある定数$c_1,\ c_2$が存在して、
ノルムが微小な任意のベクトル$\textbf{h}$ $= (h_1,h_2)$に対して
$f(x_0+h_1,y_0+h_2) = f(x_0,y_0) + c_{1}h_1 + c_{2}h_2 + \delta(h_1,h_2)$(注1参照のこと)$\qquad (a)$
ここで、$\lim_{\bf{h}\to \bf{0}}\delta(h_1,h_2)/\|{\bf h} \| = 0 (注2参照のこと)\qquad (b)$
この時、 $\textbf{c} \triangleq (c_1, c_2)$ を、fの点$(x_0,y_0)$における導値(derivative)または微分係数といい、
$f'(x_0,y_0), \ Df(x_0,y_0)$ などと書く。
(注1)$(x_0,y_0)\in U$ で、Uが開集合なので、
$\|\textbf{h}\|$がある正数より小さければ$(x_0+h_1,y_0+h_2)\in U$となり、
関数$f$は、この点で定義されている。
$\delta(h_1,h_2)$は、$h_1,h_2$ の関数である。
(注2)ノルムとしては、どのp-ノルムを用いても良い。
このテキストの「1.4.3 一般のノルムの定義とノルムの同等性」を参照のこと。
定理4
$f$が 点$(x_0,y_0)\in U$ で微分可能ならば、
1)$f$ は$(x_0,y_0)$ で偏微分可能で、
式(a)の$ c_{1}, c_{2} $ はそれぞれ、点$(x_0,y_0)$ でのx、yに関する偏微分係数である。
すなわち、$f'(x_0,y_0)=(f_{x}(x_0,y_0),f_{y}(x_0,y_0))$
2)${\bf e}=(e_1,e_2)^{T}$ を任意のベクトルとすると、
$f$ は 点$(x_0,y_0)$ で ${\bf e}$方向に微分可能で、
$ D_{\textbf{e}}f(x_0,y_0)=Df(x_0,y_0){\bf e}$
証明
1)を示そう。
式(a) で、$h_2=0$ とすると
$f(x_0+h_1,y_0) = f(x_0,y_0) + c_{1}h_1 + \delta(h_1,0)\qquad (c)$
ここで、
$\lim_{h_1 \to 0,h_1 \neq 0}\frac{\delta(h_1,0)}{|h_1|} = 0\qquad (d) $
式(c)の両辺を、$h_1 (\neq 0)$ で割り、整頓すると、
$ \frac{f(x_0+h_1,y_0) - f(x_0,y_0)}{h_1}= c_{1}+\frac{\delta(h_1,0)}{h_1}$
この式の両辺の極限$h_1\to 0$をとると、式(d)から
$\lim_{h_1 \to 0,h_1 \neq 0}\frac{f(x_0+h_1,y_0) - f(x_0,y_0) }{h_1}= c_{1} $
を得る。
この左辺は、xに関する偏微分$\frac{\partial f}{\partial x}(x_0,y_0)$の定義式である。
式(a) で、$x=x_0$ と固定すると,同様の議論で、
$c_2=\frac{\partial f}{\partial y}(x_0,y_0)$ を得る。
1)の証明終わり
2)を証明しよう。
${\bf e} = {\bf 0}$ の時は、$D_{\textbf{e}}f(x_0,y_0)=0$であることは、方向微分の定義から直ちにわかるので、2)は成り立つ。
${\bf e} \neq \textbf{0}$ の時;
方向微分の定義から
$ D_{\textbf{e}}f(x_0,y_0)=\lim_{t\to 0,t\neq 0}\frac{f(x_0+te_1,y_0+te_2)-f(x_0,y_0)}{t} \qquad (a)$
他方、fが 点$(x_0,y_0)$ で全微分可能なので、
$ f(x_0 + te_1,y_0 + te_2)-f(x_0,y_0)=Df(x_0,y_0)t\textbf{e}+\delta(te_1,te_2))
\qquad (b)$
ここで、$\frac{\delta(te_1,te_2)}{\|t{\bf e}\|} \to 0 \quad (\|t{\bf e}\| \to 0 のとき)$
式(b)を式(a)の右辺の代入すると、
$D_{\bf {e}}f(x_0,y_0)=\lim_{t\to 0,t\neq 0}\Bigl(Df(x_0,y_0){\bf e}+\frac{\delta(te_1,te_2)}{t}\Bigr)=Df(x_0,y_0){\bf e}$
これで2)が示せた。
証明終わり
$f$が微分可能ならば、
$f$の点$(x_0,y_0)$での値と、その近くの点$(x_0+h,y_0+k)$での値の差$f(x_0+h,y_0+k)-f(x_0,y_0)$ は、
$ c_1 h + c_2 k = (c_1,c_2)(h,k)^{T}=\Bigl(f_{x}(x_0,y_0),f_{y}(x_0,y_0)\Bigr)(h,k)^{T}$
で大変精度よく近似できることを意味する。
ここで、ベクトルの右肩についているTという記号は、転置演算を表す記号である。
本テキストの8.1 平面と空間,ベクトルの行列を参照のこと。
定理5
2変数関数関数 $f(x,y)$ を考える。
もし、偏導関数 $f_{x},f_{y}$ の少なくとも一方が $(x_0,y_0)$ で存在し、
他方が、$(x_0,y_0)$ を中心とする半径$\delta$ の開球体 $B_{\delta}(x_0,y_0)$上で存在し、$(x_0,y_0)$ で連続ならば、
$f(x,y)$ は$(x_0,y_0)$ において、微分可能である。
(注)$\delta$はどんなに小さくてもよい。
証明
$f_{x}$が $B_{\delta}(x_0,y_0)$上で存在し、$(x_0,y_0)$ で連続と仮定して、証明すればよい。(他の場合も同様に議論できるから)。
そこで、$f_{x}$が$B_{\delta}(x_0,y_0)$上で存在し、$(x_0,y_0)$ で連続としよう。
$\|\textbf{h}\|_{2}\lt \delta $ を満たす任意の2次元ベクトル${\textbf h}=(h_{1},h_{2})$をとる。
$f(x_0 + h_1,y_0 + h_2)-f(x_0,y_0) $
$= \Bigl(f(x_0+h_1,y_0+h_2)-f(x_0,y_0+h_2)\Bigr)+\Bigl(f(x_0,y_0+h_2)-f(x_0,y_0)\Bigr) \qquad \qquad (a)$
一変数$h_1$の関数
$\phi(h_1)\triangleq f(x_0+h_1,y_0+h_2) \qquad (b)$
を考えると、$\phi(0)=f(x_0,y_0+h_2)$であり、
$f_{x}$が$U_{\delta}(x_0,y_0)$上で存在するので、微分可能な関数である。
一変数の微分可能な関数の平均値の定理から、ある正数$\theta \in (0,1)$ が存在して、
$\phi(h_1)-\phi(0)=h_1{\phi}'(\theta h_1)$
式(b)を用いて、この式を関数$f$を用いて表すと
$f(x_0+h_1,y_0+h_2)-f(x_0,y_0+h_2)=h_{1}D_{x_1}f(x_0+\theta h_1,y_0+h_2) \qquad (c)$
式(a)の右辺の第2項$f(x_0,y_0+h_2)-f(x_0,y_0)$ を考える
関数$f$の$y$についての偏微分$D_{y}f$が$(x_0,y_0)$で存在することから、
$f(x_0,y_0+h_2)-f(x_0,y_0) = h_2 D_{y}f(x_0,y_0)+\delta(h_2)
\qquad (d)$
ここで$\delta(h_2)$は、$\lim_{h_2\to 0,h_2\neq 0}\frac{\delta(h_2)}{|h_2|}=0$をみたす関数。
式(a)の右辺に、式 (c),(d)を代入すると、
$f(x_0+h_1,y_0+h_2)-f(x_0,y_0) $
$=h_{1}D_{x}f(x_0+\theta h_1,y_0+h_2)+h_2D_{y}f(x_0,y_0)+\delta(h_2)$
$=h_{1}D_{x}f(x_0,y_0)+h_2D_{y}f(x_0,y_0)+ h_{1}\Bigl(D_{x}f(x_0+\theta h_1,y_0+h_2))- D_{x}f(x_0,y_0)\Bigr) + \delta(h_2) \qquad (e)$
$\lim_{\bf{h}\to 0,\bf{h}\neq 0}\frac{h_{1}\Bigl(D_{x}f(x_0+\theta h_1,y_0+h_2))- D_{x}f(x_0,y_0)\Bigr) + \delta(h_2) }{\|\bf{h}\|} = 0\qquad (f)$
を示せば、微分可能性の定義から、所要の命題が証明できたことになる。
$\lim_{\bf{h}\to 0,\bf{h}\neq 0}\frac{\delta(h_2)}{\|\bf{h}\|}=0$は明らか。
$\lim_{\bf{h}\to 0,\bf{h}\neq 0}\frac{h_{1}\Bigl(D_{x}f(x_0+\theta h_1,y_0+h_2))- D_{x}f(x_0,y_0)\Bigr)}{\|\bf{h}\|} $
$=\lim_{\bf{h}\to 0,\bf{h}\neq 0}\frac{h_{1}}{\|\bf{h}\|}\Bigl(D_{x}f(x_0+\theta h_1,y_0+h_2))- D_{x}f(x_0,y_0)\Bigr)$
$\qquad \frac{h_{1}}{\|\bf{h}\|}$ は絶対値が1以下の値で
$\qquad D_{x}f$ は、仮定から $(x_0,y_0)$ で連続なので
$\qquad \lim_{\bf{h}\to 0,\bf{h}\neq 0}\Bigl(D_{x}f(x_0+\theta h_1,y_0+h_2))- D_{x}f(x_0,y_0)\Bigr) = 0$が成り立つので
$= 0$
これで式$(f)$ が示せた。定理2の証明終わり。
(注)この定理はn変数関数の場合にも、次のように拡張できる。
定理5d
n変数関数関数 $f(\bf{x})$ を考える$\Bigl(\bf{x}=(x_1,x_2,\cdots x_{n-1},x_n)\Bigr)$。
もし、偏導関数 $\{f_{x_i}\}_{i=1}^{n}$ の少なくとも一つが $\bf{x}=\bf{x^{0}}$ で存在し、
残りの全ての偏導関数が$x^{0}$ を中心とする半径$\delta$ の開球体 $B_{\delta}(\bf{x^{0}})$上で存在し、$\bf{x^{0}}$ で連続ならば、
$f$ は$\bf{x^{0}}$ において、微分可能である。
証明は、同じようにしてできるので省略する。
定義4
n次元空間$\bf{R^{n}}$ の開集合Uで定義される実数値関数
$f(\bf{x})$ が$C^{1}$級 とは、
全ての偏導関数$\{f_{x_i}\}_{i=1}^{n}$がU上で存在し、
かつ、それらがU上の連続関数であること。
U上で定義され実数値をとる$C^{1}$級関数をすべて集めた集合を $C^{1}(U,\bf{R})$ と書く。
(注)n次元空間$\bf{R^{n}}$ の集合Uが開集合であるとは、
Uの任意の要素$\bf{x}$に対して、十分小さな半径rを選ぶと、
$\bf{x}$を中心とし半径rの開球体$B_{r}(\bf{x})$ がUに含まれること。
定理5d の系
$C^{1}$級の関数は微分可能である。
ベクトル値の多変数関数の微分可能性
合成関数の微分を論ずるために、微分可能性をベクトル値関数の場合に拡張する。
本§では行列の初歩的知識が必要である。
$\bf{y} = f(\bf{x})$ をn次元空間$\bf{R^n}$の開集合Uで定義され、m次元空間$\bf{R^m}$に値をとる関数とする。
ベクトル$\bf{y}$ と$\bf{x}$ を座標成分表示した縦ベクトルも同じ記号で表示しておく。
$\bf{y}=\begin{pmatrix} y_1 \\ y_{2} \\ \vdots \\ y_{m} \end{pmatrix}\qquad $
$\bf{x}=\begin{pmatrix} x_1 \\ x_{2} \\ x_{3} \\ \vdots \\ x_{n} \end{pmatrix}\qquad $
関数$\bf{y} = f(\bf{x})$を座標成分表示すると
$y_i = f^{i}({\bf x}) \qquad (i=1,2,3,\cdots m)$
定義5 ベクトル値関数の微分可能性
n変数でm次元空間$\bf{R^m}$に値をとる関数$\bf{y} = f(\bf{x})$が点$\bf{x}$で
微分可能(全微分可能ともいう)とは、
その関数を座標成分表示した、m個のn変数実数値関数
$y_i = f^{i}({\bf x}) \qquad (i=1,2,3,\cdots m)$
が全て、${\bf x}$で微分可能(全微分可能)であること。
定理6
$\bf{y} = f({\bf x})$ をn次元空間$\bf{R^n}$の開集合Uで定義され、m次元空間${\bf R^m}$に値をとる関数とする。
この関数の座標成分表示を $y^{i} = f^{i}({\bf x}) \quad (i=1,2,3,\cdots, n)\qquad \qquad (a)$とする。
1.次の条件1)と 2)は等価である。
1)関数 $f$ が、点 ${\bf x^{0}}(\in U)$ で微分可能である。
2)あるm×n行列Cが存在し、
$B_{r}({\bf x^{0}})\subset U$となるような正数rと、
大きさがrより小さい任意のn次元縦ベクトル ${\bf h}$ に対して
$f({\bf x^{0}}+{\bf h})-f({\bf x^{0}})= C{\bf h} + \delta(h_1,h_2,\cdots,h_n) \qquad (a)$
ここで、$\lim_{{\bf h}\to {bf 0},{\bf h}\neq {\bf 0}}\frac{\delta(h_1,h_2,\cdots,h_n)}{\|\bf{h}\|} = 0 $
3)$C = Df({\bf x^{0}})\triangleq
\begin{pmatrix} Df^{1}({\bf x^{0}}) \\
Df^{2}({\bf x^{0}}) \\
\vdots\\
Df^{m}({\bf x^{0}})\\
\end{pmatrix}$
$=\begin{pmatrix}
\frac{\partial f^{1}}{\partial x_1}(\bf x) & \frac{\partial f^{1}}{\partial x_2}(\bf x) & \frac{\partial f^{1}}{\partial x_3}(\bf x) & \dots & \frac{\partial f^{1}}{\partial x_n}(\bf x)\\
\frac{\partial f^{2}}{\partial x_1}(\bf x) & \frac{\partial f^{1}}{\partial x_2}(\bf x) & \frac{\partial f^{2}}{\partial x_3}(\bf x) & \dots & \frac{\partial f^{2}}{\partial x_n}(\bf x)\\
\frac{\partial f^{3}}{\partial x_1}(\bf x) & \frac{\partial f^{3}}{\partial x_2}(\bf x) & \frac{\partial f^{3}}{\partial x_3}(\bf x) & \dots & \frac{\partial f^{3}}{\partial x_n}(\bf x)\\
\vdots & \vdots & \vdots & \ddots& \vdots\\
\frac{\partial f^{m}}{\partial x_1}(\bf x) & \frac{\partial f^{m}}{\partial x_2}(\bf x) & \frac{\partial f^{m}}{\partial x_3}(\bf x) & \dots & \frac{\partial f^{m}}{\partial x_n}(\bf x)\\
\end{pmatrix}$
証明
容易なので省略する。
定義6
$\bf{y} = f({\bf x})$が点${\bf x}$で微分可能のとき
$C = Df({\bf x^{0}})$ を、関数$f$の${\bf x^{0}}$ での導値(あるいは微分係数)と呼ぶ。
定理7 合成関数の微分
${\bf y} = f({\bf x})$ を${\bf R^l}$ の開集合$U$ から$\bf{R^m}$への関数
${\bf z} = g(\bf{y})$ を$\bf{R^m}$ の開集合$V$ から${\bf R^n}$への関数とする。
もし関数$f$が点$\bf{x^0}(\in U)$で微分可能で、
${\bf y^0}\triangleq f({\bf x^0}) \in V$であり
関数$g$が点${\bf y^0}$で微分可能であるならば
合成関数
${\bf z} = g\circ f({\bf x})\triangleq g(f({\bf x}))$
は、点${\bf x^0}$で微分可能で
その点の導値 $D(g\circ f)({\bf x^0})$は
$D(g\circ f)({\bf x^0}) = Dg({\bf y^0})Df({\bf x^0}) \quad (注参照のこと)\qquad (a)$
である。
(注)右辺はn×m行列$Dg({\bf y^0})$ とm×l行列$Df({\bf x^0})$の行列としての積である。
証明
関数$f$が点${\bf x^0}$で微分可能なので、微分可能の定義から
ノルムの十分小さい任意のl次元ベクトル${\bf h}$に対して、
$f({\bf x^0}+{\bf h}) = f({\bf x^0}) + Df({\bf x^0}){\bf h} + \delta(\bf{h})\qquad (a)$
ここで、$\lim_{{\bf h}\to {\bf 0}}\frac{\delta({\bf h})}{\|{\bf h} \|}= 0\qquad (b)$
同様に、ノルムの十分小さい任意のm次元ベクトル${\bf k}$に対して、
$g({\bf y^0}+{\bf k}) = g({\bf y^0}) + Dg({\bf y^0}){\bf k} + \delta'(\bf{k})\qquad (c)$
ここで、$\lim_{{\bf k}\to {\bf 0}}\frac{\delta'({\bf k})}{\|{\bf k} \|}= 0\qquad (d)$
式(a)から、
$g\circ f({\bf x^0}+{\bf h}) = g\Bigl(f({\bf x^0}+{\bf h})\Bigr)$
$=g\Bigl(f({\bf x^0}) + Df({\bf x^0}){\bf h} + \delta({\bf h})\Bigr)$
$\qquad$ そこで、${\bf k(h)}=Df({\bf x^0}){\bf h} + \delta({\bf h}) $ とおくと
$=g\Bigl(f({\bf x^0}) + {\bf k(h)}\Bigr)$
${\bf h}$が零ベクトル近づくとき${\bf k}$も零ベクトルに近づくので、
式(c)を適用できて
$g\circ f({\bf x^0}+{\bf h})$
$=g({\bf y^0}) + Dg({\bf y^0}){\bf k(h)} + \delta'({\bf k(h)})\qquad (c)$
$
$=g\Bigl(f({\bf x^0})\Bigr) + Dg({\bf y^0})\Bigl(Df({\bf x^0}){\bf h} + \delta({\bf h})\Bigr) + \delta'\Bigl(Df(\bf{x^0})\bf{h} + \delta(\bf{h}) \Bigr)\qquad (e)$
$= g\circ f({\bf x^0})+ Dg({\bf y^0})Df({\bf x^0}){\bf h}+Dg({\bf y^0})\delta({\bf h}) + \delta'({\bf k(h)}) \Bigr)
$\lim_{{\bf h}\to {\bf 0}}\frac{\delta'({\bf k(h)})}{\|{\bf k(h)} \|}= 0$
高階偏微分
(1)二階偏微分
定義 二階偏微分
次は、大変有用な定理である。
定理
${\bf R^n}$の開集合Uで定義された実数値関数fに対し、
点$\textbf{a} \in U$ の近傍W(注参照)で
$ \qquad \qquad f_{x_i,x_j} \ f_{x_j,x_i}$
が共に存在し、$\textbf{a}$において共に連続ならば、
$ \qquad \qquad f_{x_i,x_j}(\textbf{a}) = f_{x_j,x_i}(\textbf{a})$
二階偏微分可能な関数
議論を簡単にするため、この§でも2変数関数$f(\textbf{x})$ で考える。
ここで$\textbf{x}=(x_1,x_2)^{T}$