ソースを表示

=「9.1　多変数解析学」　=
== 序 ==
本章の冒頭の偏微分の導入部については下記の本も参考にしてください。<br/>
*[[wikibooks_ja:解析学基礎/多変数関数の微積分|多変数関数の微積分]]
それ以降の内容については、ウィキブックスには殆どないため、<br/>
このテクストで今後叙述する予定です。<br/>
==多変数の実数値関数の微分   ==
${\bf R^n}=\{(x_1,x_2,,,x_n)  \mid x_i\in{\bf R},i=1,2,\cdots n\}$ の開区間 <br/>
$I^n=\prod_{i=1}^{n}(a_i,b_i)$上で定義された実関数 $y=f(x_1,x_2,,,x_n)$ を考える。<br/>
一変数関数の議論から類推するために<br/>
以後、$\vec{x}:=(x_1,x_2,,,x_n)$とおき、　$y=f(\vec{x})$　と書くこともある。<br/>
$I^n \,$上で定義された実数値関数　$\ y=f(\vec{x})=f(x_1,x_2,,,x_n)\,$　の微分について説明する。<br/>
一変数の微分から類推すると<br/>
微小なベクトル　$\vec h=(h_1,h_2,,,h_n)$　を考え、極限<br/>
$\lim_{\vec h \to 0,\vec h\neq 0}\frac{f(\vec x + \vec h)-f(\vec x)}{{\bf h} }$<br/>
が存在するとき、関数ｆは微分可能と定義することが考えられる。<br/>
しかし残念ながら、<br/>
$\vec h$はｎ次元ベクトルなので、割り算は不可能でありこの定義は無効である。
===偏微分===
関数$f$　の変数　$\vec{x}$　の第i成分　$x_i$　だけを変数とし、<br/>
他の変数は任意の実数に固定$\Bigl(x_j = a_j \quad (j\neq i)\Bigr)$して得られる関数<br/>
$\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) $<br/>
を考える。<br/>
この関数は、一変数なので、任意の点$x_i $　での微分係数　<br/>
$\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}$<br/>
$=\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}$<br/>
を考えることができる。<br/><br/>
定義（偏微分）<br/>
もし、一変数関数　$\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$で微分可能ならば、<br/>
関数ｆは、点$\vec a = (a_1.a_2,,,,a_n)$で,$x_i$　について'''偏微分可能'''であると言い,<br/>
$\frac{\partial f}{\partial x_i}(\vec a) \triangleq \frac{d\phi_{x_j=a_j,j\neq i}}{dx_i}(a_i)$<br/>
を、$f(\vec{x})$　の 点$\vec a$　での変数　$x_i$ 　についての'''偏微分係数'''という。<br/><br/>
'''定義（偏導関数）'''<br/>
$f(\vec{x})$ 　がどの点$\vec{x}$でも　$x_i$　に関して偏微分可能であるならば、<br/>
任意の点$x_i$　にその点の偏微分係数$\frac{d\phi^i}{dx_i}(x_i)$を対応させると、新しい関数が得られる。<br/>
これを、$f(\vec{x})$ 　の　$x_i$　に関する偏導関数といい、記号<br/>
$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$<br/>
などで表示する。<br/><br/>
*[[wikipedia_ja:偏微分 |ウィキペディア(偏微分)]]
以後、簡単のために2変数　x.y　の関数に限定して議論する。
定理（合成関数の微分）<br/>
$R^2$　から　$R$　への関数$f(x,y)$　と<br/>
$R$　から　$R$　への関数$g(t)$　の合成関数　<br/>
$h(x,y)=g(f(x,y))$　<br/>
を考える。<br/>
もし、$f(x,y)$　が　$(x_0,y_0)$　で、ｘに関して偏微分可能で,<br/>
$\qquad g(t)$　が、$t_0=f(x_0,y_0)$　において微分可能ならば、<br/>
$h(x,y)=g(f(x,y)$　は　$(x_0,y_0)$　で、ｘに関して偏微分可能であり,<br/>
$h_{x}(x_0,y_0)=g'(t_0)f_{x}(x_0,y_0)  \qquad \qquad ()$<br/>
証明<br/>
ｙを $ y_0 $ に固定して考えると、一変数関数の合成関数の微分になるので、合成関数の微分公式を適用すればよい。<br/><br/>
====　平均値の定理　 ====
定理<br/>
f(x,y)　を <br/>
$(x_0,y_0)$ の[[wikipedia_ja:開|開]]で[[wikipedia_ja:凸集合|凸]]な近傍 U 上で、ｘについて偏微分可能とする。<br/>
もし$(x,y)$ を近傍Ｕの点ならば<br/>
$x_0$ と $x$ の間の $\xi$　が存在して、<br/>
$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  ()$ <br/>
(注）例えば、中心$(x_0,y_0)$、半径ｒの小さな開球体$S_{r}(x_0,y_0)\triangleq \{(x,y)\in R^2 | \|(x,y) - (x_0,y_0) \| \lt r\} $　など。 <br/> 
証明 <br/>
$ \phi(x)\triangleq f(x,y) - f(x,y_0)$ とおくと、<br/>
式（）の左辺$ = \phi(x) - \phi(x_0)$<br/>
$\quad \phi(x) $　は、$x_0$　の近傍で微分可能なので、[[wikipedia_ja:平均値の定理|平均値の定理]]から、<br/>
$\quad x_0$ と $x$ の間の $\xi$　が存在して、<br/>
$= (x-x_0){\phi}'(\xi) = (x-x_0)\bigl(f_{x}(\xi,y)-f_{x}(\xi,y_0)\bigr) $<br/><br/>
定理　<br/>
f(x,y)　を <br/>
$(x_0,y_0)$ の[[wikipedia_ja:開|開]]で[[wikipedia_ja:凸集合|凸]]な近傍 U 上で、ｘについて偏微分可能とする。<br/>
もし$(x,y)=(x_0+h,y_0+k)$ を近傍Ｕの点ならば<br/>
$f(x,y) = f(x_0,y_0) + hf_{x}(x_0 + h\theta,y) + kf_{y}(x_0,y_0+ k\theta)$<br/>を満たす、$\theta=\theta(h,k) \in (0,1)$　が存在する。<br/>
証明<br/>
$ g(t) \triangleq f(x_0+ht,y) + f(x_0,y_0+kt) $ というｔの関数を導入する。<br/>
すると、<br/>
$g(1)-g(0) = f(x,y)+f(x_0,y)-\bigl(f(x_0,y)+f(x_0,y_0) \bigr)$ <br/>
$\qquad \qquad =f(x,y) - f(x_0,y_0)$ <br/>
関数　$g(t)$　は、閉区間[0,1]　を含む開区間上で微分可能なので、<br/>
一変数の微分可能関数の平均値の定理から、<br/>
ある数　$\theta \in (0,1)$　が存在して、<br/>
$g(1)-g(0) =　g'(\theta)(1-0) = g'(\theta)　\qquad \qquad (a)$<br/>
故に、$ f(x,y) - f(x_0,y_0) = g(1)-g(0) = g'(\theta)$<br/>
$\qquad $ 関数ｇの微分は,一変数関数の合成関数の微分公式から<br/>
$\qquad　g'(t) = f_{x}(x_0+ht,y)h + f_{y}(x_0,y_0+kt)k  \qquad (b)$<br/>
式(a)、(b) から<br/>
$ 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)$<br/>
証明終わり<br/><br/>
====　高階偏微分 ====
（１）二階偏微分<br/>
定義　二階偏微分<br/>

次は、大変有用な定理である。<br/>
定理<br/>
${\bf R^n}$の開集合Ｕで定義された実数値関数ｆに対し、<br/>
点$\textbf{a} \in U$ の近傍Ｗ（注参照）で<br/>
$ \qquad \qquad f_{x_i,x_j} \ f_{x_j,x_i}$<br/>
が共に存在し、$\textbf{a}$において共に連続ならば、<br/>
$ \qquad \qquad f_{x_i,x_j}(\textbf{a}) = f_{x_j,x_i}(\textbf{a})$<br/>

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

このように考えると、点$\vec x = (x_1,x_2,,,x_n)$を、座標ベクトル$\vec{e}_i$に平行ではなく、<br/>
任意に指定するベクトル$\vec a$に平行に微小量動かすときの関数ｆの変化率を考えることもできることが分かるだろう。<br/><br/>
'''定義　方向微分'''<br/>
関数$y=f(x_1,x_2,,,x_n)$の、点$\vec x = (x_1,x_2,,,x_n)$での,$\vec a$　方向の微分係数とは、<br/>
$\lim_{h\to 0,h\neq 0}\frac{f(\vec x + h\vec a)-f(\vec x )}{h}$<br/>
のことで、<br/>
$\frac{\partial f}{\partial \vec{a}}(x),\quad f_{\vec a}(x),\quad D_{\vec a}(x)$<br/>
などと書く。<br/><br/>
命題<br/>
(1)　$\vec{e_i}$　方向の微分は、$\vec{e_i}$　座標軸($x_i$座標軸）に関する偏微分である。<br/>
ここで、$\vec{e_i}$　は$x_i$座標軸の正方向向きの単位長さのベクトル。<br/>
式で書くと、<br/>
$\frac{\partial f}{\partial \vec{e_i}}(x) = \frac{\partial f}{\partial x_i}(x) $<br/>
（２）$\alpha$　を任意の実数とすると<br/>
$\frac{\partial f}{\partial \alpha \vec{e_i}}(x) = \alpha \frac{\partial f}{\partial x_i}(x) $<br/>

===微分(全微分）　===
この§も、2変数関数で説明する。<br/>
二変数関数の微分可能性をどう定義したらよいだろうか？<br/>
一変数関数の微分の場合、それと同等の条件はいくつか知られているが、<br/>
その中で二変数関数に容易に拡張できるものを採用するのが自然である。<br/>
[[物理/解析入門（１）実数の性質、連続関数、導関数と微分#微分係数の意味|1.4.1.1 微分係数の意味]]　の命題の条件　３）が、それに該当する。<br/><br/>
'''定義１；微分可能性'''(全微分可能性）<br/>
関数f(x,y)が、或る開集合Ｕ上で定義されているとする。<br/>
ｆが 点$(x_0,y_0)\in U$ で'''微分可能'''とは、<br/>
ある定数$c_1,\ c_2$が存在して、<br/>
$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)$<br/>
ここで、$\lim_{(x,y)\to (x_0,y_0)}\delta(x,y;x_o,y_0)/\|(x,y)-(x_0,y_0) \| = 0 \qquad (b)$ <br/>
(注）[[物理/平面と空間,ベクトルの性質#一般のノルムの定義とノルムの同等性|「一般のノルムの定義とノルムの同等性」]]のノルムの同等性から、<br/>
ノルムとしては、どのｐ－ノルムを用いても良い。<br/><br/>
定理１<br/>
ｆが 点$(x_0,y_0)\in U$ で微分可能ならば、偏微分可能で、<br/>
式(a)の$ c_{1}, c_{2} $　はそれぞれ、$(x_0,y_0)$ でのｘ、ｙに関する偏微分係数である。<br/>
証明<br/>
式(a)　で、$y=y_0$　と固定すると<br/>
$f(x,y_0) = f(x_0,y_0) + c_{1}(x-x_0) + \delta(x,y_0;x_o,y_0)\qquad (c)$<br/>
$\lim_{(x,y_0) \to (x_0,y_0)}\delta(x,y_0;x_o,y_0)/\|(x,y_0)-(x_0,y_0) \|  $<br/>
$= \lim_{x \to x_0}\delta(x,y_0;x_o,y_0)/|x - x_0| = 0　\qquad (d)$<br/>
式(c)の両辺を、$x-x_0(\neq 0)$ で割り、整頓すると、<br/>
$ \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}$<br/>
この式の両辺の極限$x\to x_0$をとると、<br/>
$\lim_{x\to x_0x\neq x_0}\frac{f(x,y_0) - f(x_0,y_0) }{x-x_0} = c_{1} $<br/>
を得る。<br/>
この左辺は、xに関する偏微分$\frac{\partial f}{\partial x}(x_0,y_0)$である。<br/>式(a)　で、$x=x_0$　と固定すると,同様の議論で、<br/>
$c_2=\frac{\partial f}{\partial y}(x_0,y_0)$　を得る。<br/>
証明終わり<br/><br/>
ｆが微分可能ならば、<br/>
ｆの点$(x_0,y_0)$での値と、その近くの点$(x_0+h,y_0+k)$での値の差$f(x_0+h,y_0+k)-f(x_0,y_0)$　は、<br/>
$ 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}$<br/>
で大変精度よく近似できることを意味する。<br/>
定理２；<br/>
関数 $f(x,y)$ が点$(x_0,y_0)$ の開近傍上で偏微分可能で、<br/>
$f_{x},f_{y}$ の一方、あるいは両方が、$(x_0,y_0)$で連続ならば、<br/>
$f(x,y)$ は$(x_0,y_0)$ において、微分可能である。<br/>
証明<br/>

系；$C^{1}$級の関数は微分可能<br/>
定理３　（合成関数の微分）<br/>
２つの２変数の実関数$x_{1}=x_{1}(\xi_{1},\xi_{2}),\ x_{2}=x_{2}(\xi_{1},\xi_{2})$を、<br/>
共に、点$(\xi_{1}^{0},\xi_{2}^{0})$ において微分可能、<br/>
2変数の実関数 $f(x_1,x_2)$ が、<br/>
点 $(x_{1}^{0},x_{2}^0)=\bigl( x_{1}(\xi_{1}^{0},\xi_{2}^{0}),x_{2}(\xi_{1}^{0},\xi_{2}^{0})\bigr) $ において微分可能とする。<br/>
すると、合成関数<br/>
$\qquad \qquad g(\xi_{1},\xi_{2})\triangleq f\bigl( x_{1}(\xi_{1},\xi_{2}),x_{2}(\xi_{1},\xi_{2})\bigr) $<br/>
は、$(\xi_{1}^{0},\xi_{2}^{0})$ で微分可能であり、<br/>
$g_{\xi_{1}}(\xi_{1}^{0},\xi_{2}^{0})$<br/>
$=\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}) $<br/>
$g_{\xi_{2}}(\xi_{1}^{0},\xi_{2}^{0})$<br/>
$=\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}) $<br/>
<br/>