Processing math: 100%
スポンサーリンク

「偏微分?」「高階微分?」「Ck級?」【解析学の基礎シリーズ】偏微分編 その1

微分法

本記事の内容

本記事は偏微分、高階偏微分、Ck級について解説する記事です。

本記事を読むにあたり、1変数関数の微分と多変数関数の極限について知っている必要があるため、以下の記事も合わせてご覧ください。

※シリーズ化しているので、シリーズの一部の記事のリンクを掲載しています。

↓多変数関数の極限の記事

多変数関数の微分のちょっとしたイントロ①

1変数実数値関数と同様にして、多変数関数にも微分法を考えます。
むしろ、多変数関数の微分法を考えることのほうが多い印象があります。
というのも、1変数実数値関数f:RRというのは、何か1つの実数xに対して、実数f(x)を対応させる写像のことなので、値がxにしか依存しません。
つまり、何か1つの対象が変数xにしか縛られない、ということでxさえ分かってしまえば対象が何であるか確定する、ということになります。

しかし、現実はそう単純ではありません。
例えば、インクを水に垂らしたときに、そのインクが水の中でどんどん広がっていく現象を記述したりする方程式を拡散方程式といいますが、拡散方程式は変数が複数存在する微分方程式(微分が含まれる方程式)です。
ちなみに、拡散方程式は花粉の広がり方などに使われたりもします。

他には、流体中の粒子の運動を記述したNavier-Stokes方程式も変数が複数存在する微分方程式(これを偏微分方程式といいます)です。

(※ちなみに、Navier-Stokes方程式はNewtonの古典力学の運動方程式を流体力学版へと改変したものです。つまり、Navier-Stokes方程式も運動方程式です。)

このように、実は変数がただ1つだけの方程式よりも複数存在する方程式の方がより多く出現します(と断言してしまうのはマズいかもしれないけれど、筆者はそういう印象です)。

余談 ちなみに…Navier-Stokes方程式はNewtonの古典力学の運動方程式を流体力学版へと改変したものです。
つまり、Navier-Stokes方程式も運動方程式です。

少々筆者のパーソナルな話になってしまいますが、高校生のときの筆者が進学する大学を選ぶときに、「世の中の現象を数式で表せるなんて凄い!面白そう!これがやりたい!」と思って大学を選びました。
とはいえ、実際やってみたらより純粋な数学のほうが面白くなってしまって専門にはしませんでしてけれどね(笑)

多変数関数の微分のちょっとしたイントロ②(微分から偏微分へ)

1変数の微分から多変数の微分についての橋渡しをしたいと思います(うまくできるかな…)。

1変数の微分法のチャラい復習

微分可能、微分係数

I=(a,b)Rの開区間とし、f:IRとする。 このとき、tIに対して c=limh0f(t+h)f(t)h というcRが存在するとき、ft微分可能であるといい、cf導値、または微分係数という。このとき、 c=f(t)=dfdt(t)=ddtf(t)=df(t)dt=(f(t))=(Df)(t) と書く。

でした。
詳しくは【解析学の基礎シリーズ】1変数実数値関数の微分編 その1を御覧ください。

偏微分のイントロ+全微分との関係

多変数の場合だと、変化量hは実数ではなくベクトルになります。
従って、このままでは使えません。

実は、多変数関数には微分が2種類存在します。

  • 全微分f(a)
  • 偏微分fxj(a)

実は、1変数の微分f(a)とよく対応するのは、全微分の方です(それ故同じ記号を使っています)。

例えば、1変数実数値関数fのグラフy=f(x)上の点(a,f(a))における接線の方程式は
y=f(a)(xa)+f(a)
ですが、多変数実数値関数fのグラフy=f(x)上の点(a,f(a))においては接線ではなく接平面になります。
その方程式は
y=f(a)(xa)+f(a)
です。
形式的には違いはなくて、覚えるという意味では苦労がほぼ無いと思います。

また、fがベクトル値f=(f1f2fm)の場合、
fxj(a)=(f1xj(a)f2xj(a)fmxj(a))
となります。

全微分の偏微分の関係性はある種、簡単です。

f(a)=(fx1(a) fx2(a)  fxn(a) )=(f1x1(a)f1x2(a)f1xn(a)f2x1(a)f2x2(a)f2xn(a)fmx1(a)fmx2(a)fmxn(a))=(fixj(a))
です。
つまり、

全微分係数は偏微分係数を成分とする行列のこと。

です。

要するに、今までの1変数実数値関数の微分係数は接線の傾きだったのに対して、多変数の場合は接平面の傾きに対応する、ということです。
そして、多変数の微分も結局の所は1変数の微分と形式的にはさして変わらないぜということです。

直感的な偏微分のイメージ

本質的に1変数の場合と何が違うか、というと、1変数の場合は微分に方向が1つしかありませんでした。
というのも変数が1つだったからです。
しかしながら、多変数の場合は方向が複数あります。
例えば2変数だったとすると、x方向(横の方向)とy方向(奥行きの方向)の2つの方向があります。
それ故、xで微分する、yで微分する、ということになるのです。

先程の記号fxfyはそれぞれx方向のfの変化量、y方向のfの変化量を表しています。

で、偏微分って何スか?

「偏微分とはこれです!」と言いたいところですが、まずは例を挙げます。

偏微分

例1.(2変数2次関数) 実定数a,,b,c,d,p,q,rRに対して、
f(x,y)=ax2+bxy+cy2px+qy+r((x,y)R2)
として、f:R2Rを定めるとき、
fx=2ax+by+p
です。

「おいおい。偏微分についてまだ何も言ってないのにいきなりヘンチクリンな記号を使って書かれたってわかんねんよ。」となると思いますが、大丈夫です。
後でしっかり解説しますが、ここでは1変数関数のxでの微分をdfdxと書いたのに対して、多変数関数のxでの微分をfxと書いた、と思ってくれてOKです。

多変数関数に対して、ある1つの変数xで微分することをxで偏微分する、といいますが、偏微分する際は他の変数(今回はy)を定数とみなして微分します。
同様にして
fy=bx+2cy+q2fx2=x(fx)=2a2fyx=y(fx)=b2fxy=x(fy)=b2fy2=y(fy)=2c
です。

ちなみに、fの全微分f(x,y)
f(x,y)=(fx fy)=(2ax+by+p bc+2cy+q)
となります。
しかし、これは(2ax+by+pbc+2cy+q)
ではありません。
これはf(x,y)という記号で表されます。

で、結局偏微分って何なんですか?

そんな怒らないでくださいよ。ちゃんと説明しますから。

偏微分可能、偏微分係数

ΩRnの開集合、f:ΩRma=(a1an)Ωj={1,,n}とする。fが点aで変数xjについて偏微分可能であるとは、極限 limh0f(a+hej)f(a)h が存在することを言う。ただし、ejは第j成分が1で、それ以外の成分すべて0であるようなRnのベクトルとする。すなわち、 ej=(00100)j である。
 このとき、極限値Rmfの点aでの変数xjについての偏微分係数(the partial derivative of f in the direction xj at a)と呼び、 fxj(a),xjf(a),fxj(a) などの記号で表す。

さて、上記の
f(a+hej)f(a)h
をベクトルの記法を使わずに成分でもって表すと、
f(a+hej)f(a)h=f(a1,,aj1,aj+h,aj+1,.an)f(a1,,aj1,aj,aj+1,.an)h
となります。

偏導関数、高階微分、Ck

1変数実数値関数と同様にして、高階微分、Ck級というコンセプトが定まります。

偏導関数、高階微分、Ck

ΩRnの開集合、f:ΩRmとする。
  1. 偏微分可能
  2. j{1,,n}とする。fΩxjについて偏微分可能であるとは、xΩに対して、fxで、変数xjについて偏微分可能であることをいう。このとき、写像 Ωxfxj(x)Rmfの変数xjに関する偏導関数と呼び、 fxj,xjf,fxj などの記号で表す。

  3. 1階偏導関数
  4. fx1,,fxnf1階偏導関数と呼ぶ。

  5. 記号
  6. i,j{1,,n}とする。fΩで変数xjについて偏微分可能で、偏導関数fxjΩで変数xjについて偏微分可能であるとき、xj(fxj)2fxixj,fxjxi などの記号で表す。i=jの場合、つまり2fxjxj2fx2jとも書く。

  7. 2階偏導関数
  8. 2fxixj (i,j=1,,n)f2階偏導関数と呼ぶ。

  9. k階偏導関数
  10. 4.と同様にして任意のkN (k3)に対して、fk階偏導関数が定まる。

  11. Ck
  12. kNとする。fΩCk級であるとは、fΩk階のすべての偏導関数を持ち、かつすれら全てとf自身がΩで連続であることを言う。

  13. C
  14. fΩC級であるとは、kNに対して、fΩCk級であることをいう。

  15. C0
  16. fΩC0級であるとは、fΩで連続であることをいう。f自身をf0階偏導関数ともいう。

いやあ。長い。長かった。
しかしながら、読んでみると、「ああ、1変数のときと同じじゃんね」となったかと思います。
強いて言えば、多変数ベクトル値関数の極限を知っている必要がありますが、これについてはすでに学んでいるので大丈夫かと思います。

ちなみに、チャラく復習すると、以下でした。

多変数ベクトル値関数の収束

ΩRnRnの領域、f:ΩRmを写像(関数)、aˉΩARmとする。xaのときf(x)Aに収束する(このことをf(x)A (xa)と表す)とは、 (ϵ>0)(δ>0) s.t. (xˉΩ:0<|xa|<δ|f(x)A|<ϵ) が成り立つことをいう。 Aのことをxaのときのf(x)極限と呼び、 limxaf(x)=A と表す。 また、xaのときにf(x)が極限を持つ(あるいは「f(x)が収束する」)とは、あるARmが存在して、xaのときにf(x)Aに収束することをいう。

また、多変数ベクトル値関数の連続は以下でした。

多変数ベクトル値関数の連続

ΩRnの領域、f:ΩRmとする。
  • aΩとする。faΩで連続(continuous at a)であるとは、 limxaf(x)=f(a) が成り立つことをいう。
  • すなわち、 (ϵ>0)(δ>0) s.t. (xΩ; 0<|xa|<δ|f(x)f(a)|<ϵ) が成り立つことをいう。
  • fΩで連続である(continuous on Ω)とは、任意のaΩに対して、faで連続であることをいう。
  • すなわち、 (aΩ)(ϵ>0)(δ>0) s.t. (xΩ: 0<|xa|<δ|f(x)f(a)|<ϵ) が成り立つことをいう。

いっちょ偏微分可能であることを証明してみっか。

してみましょう。

例2. f(x,y)={xyx2+y2((x,y)(0,0))0((x,y)=(0,0))とします。

このとき、fが原点O=(0,0)で偏微分可能であるかを調べてみます(まあ、偏微分可能なんですけどね)。

まず、fが原点でxで偏微分可能であるかを調べてみます。

limh0f((0,0)+h(1,0))f(0,0)h=limh0f(h,0)f(0,0)h=limh01h(h0h2+02f(0,0))=limh01h(0h2+y20)=0
従って、極限が存在するので、fは原点でxで偏微分可能です。

同様にして、
limh0f((0,0)+h(0,1))f(0,0)h=limh0f(0,h)f(0,0)h=limh01h(h002+h2f(0,0))=limh01h(002+h20)=0
となるので、fは原点でyで偏微分可能です。

例3.f(x,y)={xx2+y2((x,y)(0,0))0((x,y)=(0,0))とします。

このとき、f(0,0)で偏微分可能かを考えてみます。

まず、fが原点でxで偏微分可能であるかを調べてみます。

limh0f((0,0)+h(1,0))f(0,0)h=limh0f(h,0)f(0,0)h=limh01h(hh2+02f(0,0))=limh01h(hh20)=limh01hhh2=limh01h2=
従って、極限が存在しないので、fは原点でxで偏微分可能ではありません。

同様にして、

limh0f((0,0)+h(0,1))f(0,0)h=limh0f(0,h)f(0,0)h=limh01h(002+h2f(0,0))=limh01h(00)=0
従って、極限が存在するので、fは原点でyで偏微分可能です。

今回は、偏微分、高階微分、Ck級について解説して、実際に偏微分可能かどうかを実例をもとに計算してみました。

偏微分は結局の所、形式的にはある変数に注目したとき、それ以外を定数と見て微分するということで、これはその変数方向の変化量を表しています。

次回は、偏微分の交換可能性について解説します。

乞うご期待!質問、コメントなどお待ちしております!

コメントをする

タイトルとURLをコピーしました