物理を解説 ♪
ツイッター用のシェアボタン フェイスブック用のシェアボタン はてなブックマーク用のシェアボタン ライン用のシェアボタン
ツイッター用のシェアボタン フェイスブック用のシェアボタン はてなブックマーク用のシェアボタン ライン用のシェアボタン

偏微分の座標変換

計算は面倒だが理屈は簡単。
作成:2002/7/21
更新:2008/6/30

偏微分の変換

偏微分を含んだ式の座標変換というのは物理でよく使う.この計算は微分演算子の変換の方法さえ分かっていればまるで問題ない.

例えばデカルト座標から極座標へ変換するときの偏微分の変換式は, 数式 となるのであるが,なぜそうなるのかというところまで理解できぬまま,そういうものなのだとごまかしながら公式集を頼りにしている人が結構いたりする.学生時分の私がそうであったし,最近,読者の方からもこれについての質問を受けたので今回の説明には需要があるに違いないと判断する.

以下ではこのような変換の導き方と,なぜそのように書けるのかという考え方を説明する.式だけ示されても困る人もいるだろうから,ついでに使い方も説明しておこう.


考え方

関数f(x,y,z)xで偏微分した量∂f/∂xがあるとする.これと全く同じ量を極座標だけを使って表したい.

そのためにまずは,関数f(x,y,z)に含まれる変数x,y,zのそれぞれに次の変換式を代入してやろう. 数式 そうすることで,fの変数は(r,θ,φ)へと変わる.これによって関数の形は変わってしまうので,別の記号を使ったり,f'などと表した方がいいのかも知れないが,ここでは引き続き,変換後の関数をもfで表すことにしよう.

こういうことは特に物理では良くある.関数の記号はその形を区別するためではなく,その関数が表す物理的な意味を表すために付けられていたりすることが多いからだ.「力F」とか「ポテンシャルU」だとか「電場E」だとか,たとえ座標変換によってその関数の形が変わっても,それが表すものの内容は変わらないから,記号を変えないで使うことが多いのである.

話を進めよう.今やf(r,θ,φ)となったこの関数は,もはやxで偏微分することは出来ない.ではどうしたらいいだろう.

∂f/∂xというのは,変数のうちのxだけが変化したときのfの変化率を表していたのだった.今は,xが微小変化したらr,θ,φのいずれもが変化する可能性がある.そのことによるf(r,θ,φ)の微小変化は次のように表されるだろう.分からなければ前回の「全微分」の記事を参照してほしい. 数式 xが微小変化したことによるfの変化率を求めたいのだから,この両辺をdxで割ってやればいい. 数式 しかしこの式の表記はあまり正しくない.今は変数x,y,zのうちのxだけを変化させたという想定なので,両辺にある常微分は,この場合,すべて偏微分で書き表されるべき量なのだ. 数式 これで,xによる偏微分をr,θ,φによる偏微分の組み合わせによって表す関係が導かれたことになる.ちょっと分かりにくいだろうか.すっきりさせてみよう.ここまで関数fを使って説明してきたが,この話は別にfでなくともどんな関数でもいいわけで,この際,書くのを省いてしまうことにしよう. 数式 ただfを省いただけではないことに気が付かれただろうか.各項内で順序を入れ替えてある.fを省いただけだと∂/∂rなどは「微分演算子」になり,そのすぐ後に来るものを微分しなさいという意味になってしまうので都合が悪いからである.例えば第 1 項のfを省いてそのままの順序にしておくと,この後に来る関数に∂r/∂xを掛けてからその全体をrで微分しなさいという,意図しない意味にとられてしまう.それで式の意味を誤解されないように各項内での順序を変えておいたわけだ.

ここまではxによる偏微分を考えてきたが,他の変数についても全く同じことである. 数式


使い方

単なる繰り返しになるかも知れないが,念のためにまとめとして書いておこう.

例えば,デカルト座標で表された関数f(x,y,z)xで偏微分したものがあり,これを極座標で表された形に変換したいとする.

∂f/∂xというのは,∂/∂x fという具合に分けて書ける.この∂/∂xの部分に先ほど求めた式を代わりに入れてやればいいのだ.つまり, 数式 という具合に計算できるということである.関数fが各項に入って 3 つに増えてしまう事については全く気にしなくていい.この計算で正しいのだから.

あとは計算しやすいように,関数fを極座標を使って表してやればいい.これは簡単であろう.関数の中に含まれているx,y,z,(2) 式を代入してやれば,この関数は極座標r,θ,φだけで表された関数になる.


具体的な式を求める方法

ここまでは理論だ.あとは,∂r/∂xなどの部分を具体的に計算して求めてやれば,(1) 式のようなものが得られるはずである.そのためには,(x,y,z)(r,θ,φ)の間の関係式を使ってやればいいだろう.しかし次の関係を使って微分を計算するのは少々面倒なのだ. 数式 これで計算できないこともない.面倒だが逆関数の微分を使ってやればいいだけの話だ.しかし別の方法もある.(2) 式のようなすっきりした関係式を使う方法だ.それをこれから説明しよう.どちらの方法が簡単かは場合によって異なる.

ここで注意しなければならないことだが,例えば∂r/∂xを計算したいというので,xrで偏微分して・・・つまり∂x/∂rを計算してからその逆数を取ってやるなどという方法は使えない.微分というのは微小量どうしの割り算に過ぎないとは言ってきたが,偏微分の場合には多少意味合いが異なる.∂r/∂xyzを固定したときのrの微小変化であるが,∂x/∂rを計算する場合にrを微小変化させるとyzも変化してしまっているからである.私は以前,恥ずかしながらこのやり方で間違った結果を導いて悩み込んでしまった.

ではどうすればよいか.立場を逆転させて考えてやればいい.ここまでデカルト座標から極座標への変換を考えてきたが,極座標からデカルト座標への変換を考えれば次のようになるはずである. 数式 この計算は非常に楽であって結果はこうなる. 数式 これを連立方程式と見て逆に解いてやれば求めるものが得られる.・・・と簡単には言うものの,これは大変な作業になりそうである.それで線形代数の知識を使うのである.もともと線形代数というのは連立 1 次方程式を楽に解くために発展した学問なのだ.この式を行列形式で書いてやれば, 数式 であり,ここで出てくる 3 × 3 行列の逆行列さえ求めてやれば,それを両辺にかけることで望む形式に持っていける.では 3 × 3 行列の逆行列はどうやって求めたらいいのかそれはここでは説明しないが「クラメルの公式」「余因子行列」などという言葉を頼りにして教科書を調べてやればすぐに見つかるだろう.簡単に書いておけば,余因子行列を転置したものを元の行列の行列式で割ってやればいいだけの話だ.計算の結果は 数式 のようになり,これは初めに掲げた (1) の変換式と同じものになっている.この考えで極座標や円筒座標に限らず,どんな座標系についても計算できる.がんばれ.


2 階微分はちょっと厄介だ

もう少し説明しておかないと私は安心して眠れない.2 階微分を計算するときに間違う人がいるのではないかと心配だからだ.関数fxで 2 階微分したもの∂^2f/∂x^2,次のように分けて書くことが出来る. 数式 微分演算子が 2 つ重なるということは,fxで微分したもの全体をさらにxで微分しなさいということであり,ちゃんと意味が通っている.そう言えば高校生のときに数学の先生が,「微分の記号って言うのは実にうまく定義されているなぁ」と一人で感動していたのは,多分これのことだったのだろう.

2 階微分の座標変換を計算するときにはこの意味を崩さないように気を付けなくてはならない.どういう事かと言えば・・・そうだなぁ.(1) 式の中で∂/∂zの変換式 数式 が一番簡単そうなので例としてこれを使うことにしよう.つまり,∂^2/∂z^2というのが∂/∂zを二つ重ねたものだからといって,次のように普通に掛け算をしたのでは間違いだということである. 数式 2 行目までは問題ない.3 行目で間違いを犯している.演算子の変形は,後に必ず何かの関数が入ることを意識して行わなくてはならないのである.

例えば,∂/∂x gという形の演算子があったとする.この関数gも演算子の一部であって,これはこの後に来る関数にまずgを掛けてからその全体をxで偏微分するという意味である.分かり易いように関数fを入れて試してみよう.これは, 数式 のように計算することであろう.だからここから関数fを省いて演算子のみで表したものは 数式 という具合に変形しなければならないことが分かる.

このことを頭において先ほどの式を正しく計算してみよう.演算子の後に積の形がある時には積の微分公式を使って変形する.掛ける順番によって結果が変わることにも気を付けなくてはならない. 数式 これだけ分かっていれば,もう大抵の座標変換は問題ないだろう.しかし・・・面倒くさいな.やっぱり公式集に頼るか.(笑



趣味の物理学書店の案内バナー

EMAN物理note出張所の案内バナー