Loading [MathJax]/jax/element/mml/optable/MathOperators.js

物理/多変数解析学

提供: Internet Web School

UNIQ2dfce60676f55679-MathJax-2-QINU2 による版

目次

[非表示]

「9.1 多変数解析学」 

本章の冒頭の偏微分の導入部については下記の本も参考にしてください。

それ以降の内容については、ウィキブックスには殆どないため、
このテクストで今後叙述する予定です。

多変数の実数値関数の微分

Rn={(x1,x2,,,xn)xiR,i=1,2,n} の開区間
In=ni=1(ai,bi)上で定義された実関数 y=f(x1,x2,,,xn) を考える。
一変数関数の議論から類推するために
以後、x:=(x1,x2,,,xn)とおき、 y=f(x) と書くこともある。
In上で定義された実数値関数  y=f(x)=f(x1,x2,,,xn) の微分について説明する。
一変数の微分から類推すると
微小なベクトル h=(h1,h2,,,hn) を考え、極限
limh0,h0f(x+h)f(x)h
が存在するとき、関数fは微分可能と定義することが考えられる。
しかし残念ながら、
hはn次元ベクトルなので、割り算は不可能でありこの定義は無効である。

偏微分

関数f の変数 x の第i成分 xi だけを変数とし、
他の変数は任意の実数に固定(xj=aj(ji))して得られる関数
ϕxj=aj,ji(xi)
を考える。
この関数は、一変数なので、任意の点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)}{\bf 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)}{\bf h}
を考えることができる。

定義(偏微分)
もし、一変数関数 \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(\vec{x}) の 点\vec a での変数 x_i  についての偏微分係数という。

定義(偏導関数)
f(\vec{x})  がどの点\vec{x}でも x_i に関して偏微分可能であるならば、
任意の点x_i にその点の偏微分係数\frac{d\phi^i}{dx_i}(x_i)を対応させると、新しい関数が得られる。
これを、f(\vec{x})  の x_i に関する偏導関数といい、記号
f_{x_{i}}(\vec{x}),\quad D_{x_i}f(\vec{x}),\quad \frac{\partial f}{\partial x_i} (\vec{x}),\quad \partial f/\partial x_i
などで表示する。

以後、簡単のために2変数 x.y の関数に限定して議論する。 定理(合成関数の微分)
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 ()
証明
yを y_0 に固定して考えると、一変数関数の合成関数の微分になるので、合成関数の微分公式を適用すればよい。

 平均値の定理 

定理
f(x,y) を
(x_0,y_0)な近傍 U 上で、xについて偏微分可能とする。
もし(x,y) を近傍Uの点ならば
x_0x の間の \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 ()
(注)例えば、中心(x_0,y_0)、半径rの小さな開球体S_{r}(x_0,y_0)\triangleq \{(x,y)\in R^2 | \|(x,y) - (x_0,y_0) \| \lt r\}  など。
証明
\phi(x)\triangleq f(x,y) - f(x,y_0) とおくと、
式()の左辺 = \phi(x) - \phi(x_0)
\quad \phi(x)  は、x_0 の近傍で微分可能なので、平均値の定理から、
\quad x_0x の間の \xi が存在して、
= (x-x_0){\phi}'(\xi) = (x-x_0)\bigl(f_{x}(\xi,y)-f_{x}(\xi,y_0)\bigr)

定理 
f(x,y) を
(x_0,y_0)な近傍 U 上で、xについて偏微分可能とする。
もし(x,y)=(x_0+h,y_0+k) を近傍Uの点ならば
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)
証明終わり

 高階偏微分

(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})

方向微分

\vec{e_i} を直交座標系のx_i座標軸の正方向の方向・向きを持つ単位長さのベクトルとする(第i直交座標ベクトルと呼ぼう)。
多変数関数y=f(x_1,x_2,,,x_n)の、点\vec x = (x_1,x_2,,,x_n)での偏微分係数 \frac{\partial f}{\partial x_i}(x) は、
\vec x  を、第i座標(座標ベクトル\vec{e}_i)に平行に無限に小さい距離移動させるときの、関数fの変化率とみなせる。
式で書くと
\frac{\partial f}{\partial x_i}(x) = \lim_{h\to 0,h\neq 0}\frac{f(\vec x + h\vec{e}_i)-f(\vec x )}{h}

このように考えると、点\vec x = (x_1,x_2,,,x_n)を、座標ベクトル\vec{e}_iに平行ではなく、
任意に指定するベクトル\vec aに平行に微小量動かすときの関数fの変化率を考えることもできることが分かるだろう。

定義 方向微分
関数y=f(x_1,x_2,,,x_n)の、点\vec x = (x_1,x_2,,,x_n)での,\vec a 方向の微分係数とは、
\lim_{h\to 0,h\neq 0}\frac{f(\vec x + h\vec a)-f(\vec x )}{h}
のことで、
\frac{\partial f}{\partial \vec{a}}(x),\quad f_{\vec a}(x),\quad D_{\vec a}f(x)
などと書く。

命題
(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変数関数で説明する。
二変数関数の微分可能性をどう定義したらよいだろうか?
一変数関数の微分の場合、それと同等の条件はいくつか知られているが、
その中で二変数関数に容易に拡張できるものを採用するのが自然である。
1.4.1.1 微分係数の意味 の命題の条件 3)の式(5)が、それに該当する。

定義1;微分可能性(全微分可能性)
関数f(x,y)が、或る開集合U上で定義されているとする。
fが 点(x_0,y_0)\in U微分可能(あるいは全微分可能)とは、
ある定数c_1,\ c_2が存在して、
f(x,y) = f(x_0,y_0) + c_{1}(x-x_0) + c_{2}(y-y_0) + \delta(x,y;x_o,y_0)\qquad (a)
ここで、\lim_{(x,y)\to (x_0,y_0)}\delta(x,y;x_o,y_0)/\|(x,y)-(x_0,y_0) \| = 0 \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.4.3  一般のノルムの定義とノルムの同等性」から、
ノルムとしては、どのp-ノルムを用いても同等である。

定理1
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)\textbf{e} を任意のベクトルとすると、
fは(x_0,y_0)\textbf{e}方向に微分可能で、
D_{\textbf{e}}f(x_0,y_0)=Df(x_0,y_0)\textbf{e}
証明
1)を示そう。
式(a) で、y=y_0 と固定すると
f(x,y_0) = f(x_0,y_0) + c_{1}(x-x_0) + \delta(x,y_0;x_o,y_0)\qquad (c)
\lim_{(x,y_0) \to (x_0,y_0)}\delta(x,y_0;x_o,y_0)/\|(x,y_0)-(x_0,y_0) \|
= \lim_{x \to x_0}\delta(x,y_0;x_o,y_0)/|x - x_0| = 0 \qquad  \qquad \qquad (d)
式(c)の両辺を、x-x_0(\neq 0) で割り、整頓すると、
\frac{f(x,y_0) - f(x_0,y_0) - c_{1}(x-x_0)}{x-x_0)}=\frac{\delta(x,y_0;x_o,y_0)}{x-x_0}
この式の両辺の極限x\to x_0をとると、
\lim_{x\to x_0x\neq x_0}\frac{f(x,y_0) - f(x_0,y_0) }{x-x_0} = 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)を証明しよう。
\textbf{e} = \textbf{0} の時は、D_{\textbf{e}}f(x_0,y_0)=0であることは、方向微分の定義から直ちにわかるので、2)は成り立つ。
\textbf{e} \neq \textbf{0} の時;
方向微分の定義から
D_{\textbf{e}}f(x_0,y_0)=\lim_{t\to 0,t\neq 0}\frac{f\bigl((x_0,y_0)^{T} + t\textbf{e}\bigr)-f\bigl( (x_0,y_0)^{T} \bigr)}{t} \qquad (a)
他方、fが (x_0,y_0) で全微分可能なので、
f\bigl( (x_0,y_0)^{T} + t\textbf{e}\bigr)-f( (x_0,y_0)^{T} )=Df(x_0,y_0)t\textbf{e}+o(\|t\textbf{e}\|) \qquad \qquad (b)
式(b)を式(a)の右辺の代入すると、
D_{\textbf{e}}f(x_0,y_0)=\lim_{t\to 0,t\neq 0}\bigl(Df(x_0,y_0)\textbf{e}+\frac{o(\|t\textbf{e}\|}{t}\bigr)=Df(x_0,y_0)\textbf{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}
で大変精度よく近似できることを意味する。
定理2;
関数 f(x,y) が点(x_0,y_0) の開近傍上で偏微分可能で、
f_{x},f_{y} の一方、あるいは両方が、(x_0,y_0)で連続ならば、
f(x,y)(x_0,y_0) において、微分可能である。
証明
この証明には、次の補題を用いる。
補題
一松15ページ
f_{x}(x_0,y_0)で連続とする。
系;C^{1}級の関数は微分可能
定理3 (合成関数の微分)
2つの2変数の実関数x_{1}=x_{1}(\xi_{1},\xi_{2}),\ x_{2}=x_{2}(\xi_{1},\xi_{2})を、
共に、点(\xi_{1}^{0},\xi_{2}^{0}) において微分可能、
2変数の実関数 f(x_1,x_2) が、
(x_{1}^{0},x_{2}^0)=\bigl( x_{1}(\xi_{1}^{0},\xi_{2}^{0}),x_{2}(\xi_{1}^{0},\xi_{2}^{0})\bigr) において微分可能とする。
すると、合成関数
\qquad \qquad g(\xi_{1},\xi_{2})\triangleq f\bigl( x_{1}(\xi_{1},\xi_{2}),x_{2}(\xi_{1},\xi_{2})\bigr)
は、(\xi_{1}^{0},\xi_{2}^{0}) で微分可能であり、
g_{\xi_{1}}(\xi_{1}^{0},\xi_{2}^{0})
=\frac{\partial f}{\partial x_1}(x_{1}^{0},x_{2}^0)\frac{\partial x_{1}}{\partial \xi_{1}}(\xi_{1}^{0},\xi_{2}^{0}) + \frac{\partial f}{\partial x_2}(x_{1}^{0},x_{2}^0)\frac{\partial x_{2}}{\partial \xi_{1}}(\xi_{1}^{0},\xi_{2}^{0})
g_{\xi_{2}}(\xi_{1}^{0},\xi_{2}^{0})
=\frac{\partial f}{\partial x_1}(x_{1}^{0},x_{2}^0)\frac{\partial x_{1}}{\partial \xi_{2}}(\xi_{1}^{0},\xi_{2}^{0}) + \frac{\partial f}{\partial x_2}(x_{1}^{0},x_{2}^0)\frac{\partial x_{2}}{\partial \xi_{2}}(\xi_{1}^{0},\xi_{2}^{0})

\vec{\xi}^{0}=(\xi_{1}^{0},\xi_{2}^{0}),\ \vec{x}^{0}=(x_{1}^{0},x_{2}^{0}) とおけば、上式は
g_{\xi_{j}}(\vec{\xi}^{0}) = \sum_{i=1}^{2}f_{x_i}(\vec{x}^{0})(x_i)_{\xi_{j}}(\vec{\xi}^{0}) \quad (j=1,2)
あるいは、
g_{\xi_{j}}(\vec{\xi}^{0}) = \sum_{i=1}^{2}\frac{\partial f}{\partial x_i}(\vec{x}^{0})\frac{\partial x_i}{\partial \xi_{j}}(\vec{\xi}^{0}) \quad (j=1,2)
と書ける。
これは、さらに容易にわかる、どの点の関数値かを省略すれば
g_{\xi_{j}} = \sum_{i=1}^{2}f_{x_i}\cdot (x_i)_{\xi_{j}} \quad (j=1,2)
あるいは
g_{\xi_{j}} = \sum_{i=1}^{2}\frac{\partial f}{\partial x_i}\frac{\partial x_i}{\partial \xi_{j}} \quad (j=1,2)
と略記できる。
証明

個人用ツール