Loading [MathJax]/jax/output/HTML-CSS/jax.js

物理/多変数解析学

提供: Internet Web School

UNIQ366bd4d6f366c15-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)f(a1,a2,,,ai1,xi,ai+1,,,an)
を考える。
この関数は、一変数なので、任意の点xi での微分係数 
dϕxj=aj,jidxi(xi)limh0,h0ϕxj=aj,ji(xi+h)ϕxj=aj,ji(xi)h
=limh0,h0f(a1,a2,,,ai1,xi+h,ai+1,,,an)f(a1,a2,,,ai1,xi,ai+1,,,an)h
を考えることができる。

定義1(偏微分)
もし、一変数関数 ϕxj=aj,ji(xi)=f(a1,a2,,,ai1,xi,ai+1,,,an) が、ある点xi=aiで微分可能ならば、
関数fは、点a=(a1.a2,,,,an)で,xi について偏微分可能であると言い,
fxi(a)dϕxj=aj,jidxi(ai)
を、f(x) の 点a での変数 xi  についての偏微分係数という。

定義2(偏導関数)
f(x)  がどの点xでも xi に関して偏微分可能であるならば、
任意の点x にその点における xi に関する偏微分係数fxi(x)を対応させると、新しい関数が得られる。
これを、f(x)  の xi に関する偏導関数といい、記号
fxi(x),Dxif(x),fxi(x),f/xi
などで表示する。

以後、簡単のために2変数x,y の関数に限定して議論する。
定理1 合成関数の微分(1)
R2 から R への関数f(x,y) と
R から R への関数g(t) の合成関数 
h(x,y)=g(f(x,y)) 
を考える。
もし、f(x,y) が (x0,y0) で、xに関して偏微分可能で,
g(t) が、t0=f(x0,y0) において微分可能ならば、
h(x,y)=g(f(x,y)) は (x0,y0) で、xに関して偏微分可能であり,
hx(x0,y0)=g(t0)fx(x0,y0)(1)
証明
yを y0 に固定して考えると、一変数関数の合成関数の微分になるので、合成関数の微分公式を適用すればよい。

定理2
f(x,y) を
(x0,y0) を中心とするある半径rの開球体Br(x0,y0)上で、xについて偏微分可能とする。
もし(x,y)Br(x0,y0)の点ならば
x0x の間の ξ が存在して、
f(x,y)f(x,y0)(f(x0,y)f(x0,y0))=(xx0)(fx(ξ,y)fx(ξ,y0))(2)
(注)2次元の開球体Br(x0,y0){(x,y)R2|(x,y)(x0,y0)2<r} は、中心が点(x0,y0) で半径rの円周で囲まれる内部である。
証明
ϕ(x)f(x,y)f(x,y0) とおくと、
式()の左辺=ϕ(x)ϕ(x0)
ϕ(x) は、x0 の近傍で微分可能なので、平均値の定理から、
x0x の間の ξ が存在して、
=(xx0)ϕ(ξ)=(xx0)(fx(ξ,y)fx(ξ,y0))

定理3 
f(x.y) を
(x0,y0) を中心とする開球体Br(x0,y0)上で、xについて偏微分可能とする。
もし(x,y)=(x0+h,y0+k)Br(x0,y0) ならば
f(x,y)=f(x0,y0)+hfx(x0+hθ,y)+kfy(x0,y0+kθ)
を満たす、θ=θ(h,k)(0,1) が存在する。
証明
g(t)f(x0+ht,y)+f(x0,y0+kt) というtの関数を導入する。
すると、
g(1)g(0)=f(x,y)+f(x0,y)(f(x0,y)+f(x0,y0))
=f(x,y)f(x0,y0)
関数 g(t) は、閉区間[0,1] を含む開区間上で微分可能なので、
一変数の微分可能関数の平均値の定理から、
ある数 θ(0,1) が存在して、
g(1)g(0)= g(θ)(10)=g(θ) (a)
故に、f(x,y)f(x0,y0)=g(1)g(0)=g(θ)
関数gの微分は,一変数関数の合成関数の微分公式から
 g(t)=fx(x0+ht,y)h+fy(x0,y0+kt)k(b)
式(a)、(b) から
f(x,y)f(x0,y0)=fx(x0+hθ,y)h+fy(x0,y0+kθ)k(b)
証明終わり

方向微分

ei を直交座標系のxi座標軸の正方向の方向・向きを持つ単位長さのベクトルとする(第i直交座標ベクトルと呼ぼう)。
多変数関数y=f(x1,x2,,,xn)の、点x=(x1,x2,,,xn)での偏微分係数 fxi(x) は、
x を、第i座標(座標ベクトルei)に平行に無限に小さい距離移動させるときの、関数fの変化率とみなせる。
式で書くと
fxi(x)=limh0,h0f(x+hei)f(x)h

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

定義 方向微分
関数y=f(x1,x2,,,xn)の、点x=(x1,x2,,,xn)での,a 方向の微分係数とは、
limh0,h0f(x+ha)f(x)h
のことで、
fa(x),fa(x),Daf(x)
などと書く。

命題1
(1) ei 方向の微分は、ei 座標軸(xi座標軸)に関する偏微分である。
ここで、ei はxi座標軸の正方向向きの単位長さのベクトル。
式で書くと、
fei(x)=fxi(x)
(2)α を任意の実数とすると
fαei(x)=αfxi(x)

微分(全微分) 

この§も、記述を簡単にするため、2変数関数で説明する。
一般のn変数の場合への拡張は、記述は複雑になるが、容易である。

多変数実数値関数の微分可能性

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

定義3 微分可能性(全微分可能性)
関数f(x,y)が、或る開集合U上で定義されているとする。
fが 点(x0,y0)U微分可能(あるいは全微分可能)とは、
ある定数c1, c2が存在して、
ノルムが微小な任意のベクトルh=(h1,h2)に対して f(x0+h1,y0+h2)=f(x0,y0)+c1h1+c2h2+δ(h1,h2)(a)(注1参照のこと)
ここで、limh0δ(h1,h2)/h=0(b) (注2参照のこと)
この時、 c(c1,c2) を、fの点(x0,y0)における導値(derivative)または微分係数といい、
f(x0,y0), Df(x0,y0) などと書く。

(注1)δ(h1,h2)は、h1,h2 の関数
(注2)ノルムとしては、どのp-ノルムを用いても良い。
このテキストの「1.4.3  一般のノルムの定義とノルムの同等性」を参照のこと。

定理4
fが 点(x0,y0)U で微分可能ならば、

1)f(x0,y0) で偏微分可能で、
式(a)のc1,c2 はそれぞれ、(x0,y0) でのx、yに関する偏微分係数である。
すなわち、f(x0,y0)=(fx(x0,y0),fy(x0,y0))
2)e を任意のベクトルとすると、
f(x0,y0)e方向に微分可能で、
Def(x0,y0)=Df(x0,y0)e

証明
1)を示そう。
式(a) で、y=y0 と固定すると
f(x,y0)=f(x0,y0)+c1(xx0)+δ(h1,h2)(c)
ここで、
lim(x,y0)(x0,y0)δ(h1,h2)/(x,y0)(x0,y0)=0(d)
(x,y0)(x0,y0) xx0(x,y0)(x0,y0)=|h1| なので
式(d)から、limxx0δ(h1,h2)/|h1|=0(e)
式(c)の両辺を、xx0(0) で割り、整頓すると、
f(x,y0)f(x0,y0)xx0=c1+δ(h1,h2)xx0
この式の両辺の極限xx0をとると、式(e)から
limxx0,xx0f(x,y0)f(x0,y0)xx0=c1
を得る。
この左辺は、xに関する偏微分fx(x0,y0)の定義式である。
式(a) で、x=x0 と固定すると,同様の議論で、
c2=fy(x0,y0) を得る。
1)の証明終わり
2)を証明しよう。
e=0 の時は、Def(x0,y0)=0であることは、方向微分の定義から直ちにわかるので、2)は成り立つ。
e0 の時;
方向微分の定義から
Def(x0,y0)=limt0,t0f((x0,y0)T+te)f((x0,y0)T)t(a)
他方、fが (x0,y0) で全微分可能なので、
f((x0,y0)T+te)f((x0,y0)T)=Df(x0,y0)te+o(te)(b)
式(b)を式(a)の右辺の代入すると、
Def(x0,y0)=limt0,t0(Df(x0,y0)e+o(te)t)=Df(x0,y0)e
これで2)が示せた。
証明終わり

fが微分可能ならば、
fの点(x0,y0)での値と、その近くの点(x0+h,y0+k)での値の差f(x0+h,y0+k)f(x0,y0) は、
c1h+c2k=(c1,c2)(h,k)T=(fx(x0,y0),fy(x0,y0))(h,k)T
で大変精度よく近似できることを意味する。
ここで、ベクトルの右肩についているTという記号は、転置演算を表す記号である。
本テキストの8.1 平面と空間,ベクトルの行列を参照のこと。

定理5;
2変数関数関数 f(x,y) を考える。
もし、偏導関数 fx,fy の少なくとも一方が (x0,y0) で存在し、 他方が、(x0,y0) を中心とする半径δ の開球体 Bδ(x0,y0)上で存在し、(x0,y0) で連続ならば、
f(x,y)(x0,y0) において、微分可能である。
(注)δはどんなに小さくてもよい。
証明
fxが Bδ(x0,y0)上で存在し、(x0,y0) で連続と仮定して、証明すればよい。(他の場合も同様に議論できるから)。
そこで、fxBδ(x0,y0)上で存在し、(x0,y0) で連続としよう。
h2<δ を満たす任意の2次元ベクトルh=(h1,h2)をとる。
f((x0,y0)+h)f(x0,y0)
=(f(x0+h1,y0+h2)f(x0,y0+h2))+(f(x0,y0+h2)f(x0,y0))(a)
一変数h1の関数
ϕ(h1)f(x0+h1,y0+h2)(b)
を考えると、ϕ(0)=f(x0,y0+h2)であり、
fxUδ(x0,y0)上で存在するので、微分可能な関数である。
一変数の微分可能な関数の平均値の定理から、ある正数θ(0,1) が存在して、
ϕ(h1)ϕ(0)=h1ϕ(θh1)
式(b)を用いて、この式を関数を用いて表すと
f(x0+h1,y0+h2)f(x0,y0+h2)=h1Dx1f(x0+θh1,y0+h2)(c)
式(a)の右辺の第2項f(x0,y0+h2)f(x0,y0) を考える
関数fyについての偏微分Dyf(x0,y0)で存在することから、
f(x0,y0+h2)f(x0,y0)=h2Dyf(x0,y0)+δ(h2)(d)
ここで0(h2)は、limh20,h20δ(h2)|h2|=0をみたす関数
式(a)の右辺に、式 (c),(d)を代入すると、
f(x0+h1,y0+h2)f(x0,y0)
=h1Dxf(x0+θh1,y0+h2)+h2Dyf(x0,y0)+δ(h2)
=h1Dxf(x0,y0)+h2Dyf(x0,y0)+h1(Dxf(x0+θh1,y0+h2))Dxf(x0,y0))+δ(h2)(e)

limh0,h0h1(Dxf(x0+θh1,y0+h2))Dxf(x0,y0))+δ(h2)h=0(f)
を示せば、微分可能性の定義から、所要の命題が証明できたことになる。
limh0,h0δ(h2)h=0は明らか。
limh0,h0h1(Dxf(x0+θh1,y0+h2))Dxf(x0,y0))h
=limh0,h0h1h(Dxf(x0+θh1,y0+h2))Dxf(x0,y0))
h1h は絶対値が1以下の値で
Dxf は、仮定から (x0,y0) で連続なので
limh0,h0(Dxf(x0+θh1,y0+h2))Dxf(x0,y0))=0が成り立つので
=0
これで式(f) が示せた。定理2の証明終わり。
(注)この定理はn変数関数の場合にも、次のように拡張できる。
定理5’
n変数関数関数 f(x) を考える(x=(x1,x2,xn1,xn))
もし、偏導関数 {fxi}ni=1 の少なくとも一つが x=x0 で存在し、
残りの全ての偏導関数がx0 を中心とする半径δ の開球体 Bδ(x0)上で存在し、x0 で連続ならば、
fx0 において、微分可能である。
証明は、同じようにしてできるので省略する。

定義4
n次元空間Rn の開集合Uで定義される実数値関数
f(x) がC1とは、
全ての偏導関数{fxi}ni=1がU上で存在し、
かつ、それらがU上の連続関数であること。
U上で定義され実数値をとるC1級関数をすべて集めた集合を C1(U,R) と書く。 
(注)n次元空間Rn の集合Uが開集合であるとは、
Uの任意の要素xに対して、十分小さな半径rを選ぶと、
xを中心とし半径rの開球体Br(x) がUに含まれること。

定理5'の系
C1級の関数は微分可能である。

ベクトル値の多変数関数の微分可能性

合成関数の微分を論ずるために、微分可能性をベクトル値関数の場合に拡張する。
本§では行列の初歩的知識が必要である。

y=f(x) をn次元空間Rnの開集合Uで定義され、m次元空間Rmに値をとる関数とする。
ベクトルy とx を座標成分表示した縦ベクトルも同じ記号で表示しておく。
y=(y1y2ym) x=(x1x2x3xn)
関数y=f(x)を座標成分表示すると
yi=fi(x)(i=1,2,3,m)

定義5 ベクトル値関数の微分可能性
n変数でm次元空間Rmに値をとる関数y=f(x)が点xで 微分可能(全微分可能ともいう)とは、
その関数を座標成分表示した、m個のn変数実数値関数
yi=fi(x)(i=1,2,3,m)
が全て、xで微分可能(全微分可能)であること。
定義6
y=f(x)が点xで微分可能のとき
fixでの導値(微分係数)Dfi(x)は、n次元横ベクトルであり、
このm個の横ベクトルを縦に並べたm×n行列をDf(x) と書き、
y=f(x)の点xでの導値(微分係数)と呼ぶ。
Df(x)= (f1x1(x)f1x2(x)f1x3(x)f1xn(x)f2x1(x)f1x2(x)f2x3(x)f2xn(x)f3x1(x)f3x2(x)f3x3(x)f3xn(x)fmx1(x)fmx2(x)fmx3(x)fmxn(x))

定理6
y=f(x) をn次元空間Rnの開集合Uで定義され、m次元空間Rmに値をとる関数とする。
この関数の座標成分表示を yi=fi(x)(i=1,2,3,,n)(a)とする。
この時次の2条件は等価である。
1)関数 f が、点 x0(U) で微分可能で
各 2)あるm×n行列Cが存在し、
Br(x0)Uとなるようなrと、
大きさがrより小さい任意のn次元縦ベクトル h に対して
f(x0+h)f(x0)=Ch+δ(h)(a)
ここで、limhbf0,h0δ(h)h=0

定理 (合成関数の微分)
2つの2変数の実関数x1=x1(ξ1,ξ2), x2=x2(ξ1,ξ2)を、
共に、点(ξ01,ξ02) において微分可能、
2変数の実関数 f(x1,x2) が、
(x01,x02)=(x1(ξ01,ξ02),x2(ξ01,ξ02)) において微分可能とする。
すると、合成関数
g(ξ1,ξ2)f(x1(ξ1,ξ2),x2(ξ1,ξ2))
は、(ξ01,ξ02) で微分可能であり、
gξ1(ξ01,ξ02)
=fx1(x01,x02)x1ξ1(ξ01,ξ02)+fx2(x01,x02)x2ξ1(ξ01,ξ02)
gξ2(ξ01,ξ02)
=fx1(x01,x02)x1ξ2(ξ01,ξ02)+fx2(x01,x02)x2ξ2(ξ01,ξ02)

ξ0=(ξ01,ξ02), x0=(x01,x02) とおけば、上式は
gξj(ξ0)=2i=1fxi(x0)(xi)ξj(ξ0)(j=1,2)
あるいは、
gξj(ξ0)=2i=1fxi(x0)xiξj(ξ0)(j=1,2)
と書ける。
これは、さらに容易にわかる、どの点の関数値かを省略すれば
gξj=2i=1fxi(xi)ξj(j=1,2)
あるいは
gξj=2i=1fxixiξj(j=1,2)
と略記できる。
証明

 高階偏微分

(1)二階偏微分
定義 二階偏微分

次は、大変有用な定理である。
定理
Rnの開集合Uで定義された実数値関数fに対し、
aU の近傍W(注参照)で
fxi,xj fxj,xi
が共に存在し、aにおいて共に連続ならば、
fxi,xj(a)=fxj,xi(a)

二階偏微分可能な関数

議論を簡単にするため、この§でも2変数関数f(x) で考える。
ここでx=(x1,x2)T

C2級の関数

個人用ツール