微分学
常微分
1つだけ変数をもつ関数f(x)を考えよう。この関数の常微分(ordinary derivative)df/dxから何が得られるだろうか?答えは、変数xを微小量dxだけ動かしたとき、関数f(x)がどれだけ急激に変化するかということ、である。これは次の式からも分かる。
df=(dfdx)dx
もしもxを微小量dxだけ増加させれば、fも微小量dfだけ変化する。微分はこの比例係数をあらわしている。例えば、図1.17(a)では関数f(x)はxについて緩やかに変化しているから、微分もこれに対応して小さい。一方で、図1.17(b)では関数f(x)はxについて急激に変化しているから、x=0から離れるにつれて微分も大きくなっている。
幾何学的には、導関数df/dxはxに対するfのグラフの傾きをあらわしている。
図1.17(a)
図1.17(b)
勾配
次に3つの変数を持つ関数を考えよう。ここでは室内の温度の関数T(x,y,z)を考えてみる。部屋のある角を原点にして座標軸を設ければ、室内の各点の温度が分かるだろう。先ほどの導関数df/dxを一般化して変数が3つの関数T(x,y,z)でも使えるようにしたい。
微分は、位置が微小量変化したときに関数がどれくらい急激に変化するかをあらわしていると考えられる。しかし、今の場合は変数が1つの場合と違って、位置がどの方向に微小変化するのかという問題があるため、より複雑になっている。
もし位置が高くなれば温度は急激に上がり、位置が水平に変化すれば温度は大して変わらないと思われる。実際、温度Tがどのように変化するのか?という問題は、どの方向に変化させるのかで答えが変わってくるのだ。
幸運なことに、実はこれは見た目ほど難しい問題ではない。偏微分の基本定理に依れば、
dT=(∂T∂x)dx+(∂T∂y)dy+(∂T∂z)dz
が成り立つ。このことは、3つの変数全てがそれぞれdx、dy、dzだけ微小変化したときに温度Tがどのように変化するかをあらわしている。無限個の微分が必要ではないということは注目に値する。つまり、3次元座標の方向それぞれに沿った偏微分(partial derivative)がありさえすれば良いのだ。
(34)は以下の式のように、内積を連想させる。
dT=(∂T∂xˆx+∂T∂yˆy+∂T∂zˆz)⋅(dxˆx+dyˆy+dzˆz)=(∇T)⋅(dl)
但し、
∇T:=∂T∂xˆx+∂T∂yˆy+∂T∂zˆz
であり、これはTの勾配(gradient)と呼ばれる。この∇Tは3つの成分を持ったベクトル量であり、(35)は(33)を3次元の場合に一般化した微分であったのだということに注目せよ。
勾配は、幾何学的には、任意のベクトルのように大きさと向きを持つものである。その幾何学的意味を明らかにするために、まず内積の式(35)を(1)を用いて以下のように書き直す。
dT=∇T⋅dl=|∇T||dl|cosθ
但し、θは∇Tとdlの間の角である。もし、長さ|dl|を固定してθを変化させたとき、cosθ=1となる角度θ=0でdTは最大となる。つまり、距離|dl|を固定したなら、dTは∇Tと同じ方向へ動かしたときに最大となる。従って、以下の2つが言える。
- 勾配∇Tは関数Tが最も増加する方向を向く。
- 大きさ|∇T|はこの最も増加する方向に沿った傾き(変化の割合)を与えてくれる。
丘の中腹に立っている自分を想像してみよう。
周りを見回して最も急な上り坂を見出せ。それは勾配の向きである。そして、その方向の傾斜、つまり平均変化率を測ってみよ。それはこの勾配の大きさである(先ほどと違って、ここでの関数は丘の高さであるから、緯度と経度のみに依存している。
つまりこの関数は2変数ではなく3変数であるが、勾配の幾何学的な意味はこちらの場合の方がより理解しやすい。)。(37)から分かるように、最も急な下り坂の方向は最も急な上り坂の方向の正反対の方向である。
これに対して、直角の方向、θ=90∘では勾配は0となる。つまり、勾配は等高線に対して直交している。外面はこれらの性質を持たないと思いつくかもしれないが、外面にはねじれがあり、そこが微分不可能な関数に対応している。
勾配が0というのは、どのような意味があるのだろうか?
もし、(x,y,z)で∇T=0であるならば、(x,y,z)から微小変位した点においてdT=0である。これはつまり、関数T(x,y,z)の停留点(stationary point)である。
これは最大値(頂点)、最小点(谷底)、鞍点(峠)、あるいは山の肩のいずれかである。これは1変数の関数の場合に導関数が0の点が最大点、最小点、あるいは変曲点を示していたことと類似している。
特に、3変数関数の極値を見つけたい場合にはその勾配を0に設定する必要がある。
問題
位置ベクトルの大きさr=√x2+y2+z2の勾配を求めよ。
解説
以下のように計算する。
∇r=∂r∂xˆx+∂r∂yˆy+∂r∂zˆz=122x√x2+y2+z2ˆx+122y√x2+y2+z2ˆy+122z√x2+y2+z2ˆz=xˆx+yˆy+zˆz√x2+y2+z2=rr=ˆr
この結果は、原点からの距離は動径方向で最も急激に増加し、増加の比率は動径方向で1になることを示している。これは理にかなっている結果である。