はじめに
大学に入るとベクトル演算をする機会が増えるかと思います。
ベクトルの計算は往々にして煩雑ですが、特に球座標などの曲がった座標系の場合のベクトル演算はさらに煩雑なものになります。
ここでは、3次元球座標における、ベクトルの勾配(gradient)、発散(divergence)、回転(rotation)、ラプラシアンの表式と簡単な導出方法を、備忘録を兼ねてまとめておきたいと思います。
まず先に式を書くと以下の通りになります。
3次元球座標系におけるスカラー関数の勾配は
![](gradient.png)
![](divergence.png)
![](rotation.png)
![](laplacian.png)
ここで、
![](f.png)
![](A.png)
- 導出のための下準備
- 座標変換
- 単位ベクトルの変換
- 偏微分演算子の変換
- 単位ベクトルの微分
- 単位ベクトルの外積
- 球座標での勾配(gradient)
- 球座標での発散(divergence)
- 球座標での回転(rotation)
- 球座標でのラプラシアン
- PDFファイル
導出のための下準備
座標変換
直交座標から球座標
への座標変換を行います。
両者の関係は、
![](xyzrtp.png)
![](rtpxyz.png)
![](thetaphi.png)
単位ベクトルの変換
直交座標における単位ベクトルと、球座標における単位ベクトル
の関係を導出しておきます。
球座標での単位ベクトルはそれぞれ、
![](unitvec.png)
この式に、
![](rvec.png)
![](rtpunit.png)
![](xyzunit.png)
![](rtpunit2.png)
![](rtpunit3.png)
![](rtpunit4.png)
![](matA.png)
偏微分演算子の変換
続いて、直交座標での偏微分演算子と球座標での偏微分演算子の対応関係を求めます。 偏微分演算子には以下の対応関係があります。
![](ddx.png)
![](xyz.png)
![](rtp.png)
![](ddx2.png)
![](ddxtan-1.png)
![](ddxddyddz.png)
![](ddxddyddz2.png)
単位ベクトルの微分
単位ベクトルを微分した際の関係式を導出しておきます。
この先の導出において、単位ベクトルの微分は非常に重要です。
直交座標では、単位ベクトルを
で偏微分しても0になりますが、球座標ではそうなりません。
「単位ベクトルの変換」の部分で、球座標での単位ベクトルを
を使って表したので、その式を
で偏微分します。
偏微分を行ったあとに
を
で書き直してやると、以下の関係式を得ます。
![](erbibun.png)
![](etbibun.png)
![](epbibun.png)
単位ベクトルの外積
ベクトル関数の回転を求める際には単位ベクトル同士の外積も必要なので先に示しておきます。 球座標の単位ベクトルの外積は
![](unitgaiseki.png)
![](rtp.png)
球座標での勾配(gradient)
まず、任意のスカラー関数の勾配(gradient)を導出します。
直交座標での勾配は、直交座標での偏微分演算子と単位ベクトルを用いて
![](gradf.png)
![](gradf2.png)
![](gradf3.png)
![](gradient.png)
また、この式からスカラー関数を省いて微分ベクトル演算子の形で書くと、
![](nabla.png)
球座標での発散(divergence)
次に、任意のベクトル関数の発散(divergence)を導出します。
ベクトル関数は次のように成分表記することが出来ます。
![](vecA.png)
![](divA.png)
![](divA2.png)
![](divA3.png)
![](divergence.png)
球座標での回転(rotation)
続いて、任意のベクトル関数の回転(rotation)を導出します。
先ほど同様に考えると、
![](rotA.png)
![](rotA2.png)
![](rotA3.png)
![](rotation.png)
球座標でのラプラシアン
最後に、球座標でのラプラシアンの表式を導出します。
任意のスカラー関数にラプラシアンを作用させたものは、
![](laplace.png)
![](laplace2.png)
![](laplace3.png)
![](laplacian.png)
また、演算子の部分だけ取り出すと、
![](laplace4.png)
以上の計算はPDF文書にまとめました。
ここに載せたものより少しだけ詳細な計算も書いてあります。
PDFファイル
2017年04月24日 更新