「球面座標系」の版間の差分

提供: miniwiki
移動先:案内検索
(1版 をインポートしました)
(内容を「{{テンプレート:20180815sk}}」で置換)
(タグ: Replaced)
 
1行目: 1行目:
[[File:Spherical with grid.svg|thumb|球面座標系]]
+
{{テンプレート:20180815sk}}
'''球面座標系'''(きゅうめんざひょうけい、{{Lang-en|spherical coordinate system}})とは、3次元[[ユークリッド空間]]に定まる[[座標系]]の一つで、一つの[[動径]]座標と二つの[[角度]]座標で表される[[極座標系]]である。第一の角度はある[[座標軸|軸]](通常は {{mvar|z}}-軸を選ぶ)と動径がなす角度で、第二の角度は、その軸に[[垂直]]な[[平面]]にある別の軸(通常は {{mvar|x}}-軸を選ぶ)とこの平面への動径の[[射影]]がなす角度である。通常は動径座標に記号 {{mvar|r}} を用い、第一の角度座標には {{mvar|&theta;}} を、第二の角度座標には {{mvar|&phi;}} を用いて表される。動径座標は {{math|0 &le; ''r'' < &infin;}} の範囲にあり、第一の角度は {{math|0 &le; ''&theta;'' &le; &pi;}} の範囲にある。第二の角度の動く範囲は {{math|&minus;&pi; < ''&phi;'' &le; &pi;}} もしくは {{math|0 &le; ''&phi;'' < 2&pi;}} のどちらかを用いることが多い。
 
 
 
== 座標変換 ==
 
球面座標 {{math|(''r'',''&theta;'',''&phi;'')}} から[[直交座標系|直交直線座標]] {{math|(''x'',''y'',''z'')}} への変換は
 
{{Indent|
 
<math>\begin{cases}
 
x = r\sin\theta\, \cos\phi \\
 
y = r\sin\theta\, \sin\phi \\
 
z = r\cos\theta \\
 
\end{cases}</math>
 
}}
 
で与えられる。第二の角度座標を {{math|&minus;&pi; < ''&phi;'' &le; &pi;}} とする場合は、直交直線座標 {{math|(''x'',''y'',''z'')}} から球面座標 {{math|(''r'',''&theta;'',''&phi;'')}} への変換は
 
{{Indent|
 
<math>\begin{cases}
 
r =\sqrt{x^2+y^2+z^2} \\
 
\theta =\arccos(z/\sqrt{x^2+y^2+z^2}) \\
 
\phi =\sgn(y) \arccos(x/\sqrt{x^2+y^2}) \\
 
\end{cases}</math>
 
}}
 
で与えられる。ここで {{math|sgn}} は[[符号関数]]
 
{{Indent|
 
<math>\sgn(y) =\begin{cases}
 
1 & (y\ge0) \\
 
-1 & (y<0) \\
 
\end{cases} </math>
 
}}
 
である。{{mvar|z}}-軸上 {{math|1=(''x'',''y'') = (0,0)}} において特異性があり、[[分母]]がゼロとなるため {{mvar|&phi;}} が定まらない。さらに[[原点]] {{math|1=(''x'',''y'',''z'') = (0,0,0)}} においては {{mvar|&theta;}} も定まらない。
 
 
 
球面座標 {{math|(''r'',''&theta;'',''&phi;'')}} から直交直線座標 {{math|(''x'',''y'',''z'')}} への変換の式を微分すれば
 
{{Indent|
 
<math>\begin{cases}
 
dx = \sin\theta\, \cos\phi\, dr +r\cos\theta\, \cos\phi\, d\theta -r\sin\theta\, \sin\phi\, d\phi \\
 
dy = \sin\theta\, \sin\phi\, dr +r\cos\theta\, \sin\phi\, d\theta +r\sin\theta\, \cos\phi\, d\phi \\
 
dz = \cos\theta\, dr -r\sin\theta\, d\theta\\
 
\end{cases}</math>
 
}}
 
が得られて、[[ヤコビ行列]]とヤコビ行列式は
 
{{Indent|
 
<math>\begin{align}
 
\frac{\partial(x,y,z)}{\partial(r,\theta,\phi)}
 
&=\begin{pmatrix}
 
\sin\theta\, \cos\phi & r\cos\theta\, \cos\phi & -r\sin\theta\, \sin\phi \\
 
\sin\theta\, \sin\phi & r\cos\theta\, \sin\phi & r\sin\theta\, \cos\phi \\
 
\cos\theta & -r\sin\theta & 0 \\
 
\end{pmatrix} \\
 
&=\begin{pmatrix}
 
\cos\phi & -\sin\phi & 0 \\
 
\sin\phi & \cos\phi & 0 \\
 
0 & 0 & 1 \\
 
\end{pmatrix}
 
\begin{pmatrix}
 
\sin\theta & \cos\theta & 0 \\
 
0 & 0 & 1 \\
 
\cos\theta & -\sin\theta & 0 \\
 
\end{pmatrix}
 
\begin{pmatrix}
 
1 & 0 & 0 \\
 
0 & r & 0 \\
 
0 & 0 & r\sin\theta \\
 
\end{pmatrix} \\
 
\end{align}</math>
 
}}
 
{{Indent|
 
<math>\left| \frac{\partial(x,y,z)}{\partial(r,\theta,\phi)} \right| =r^2 \sin\theta</math>
 
}}
 
となる。従って球面座標で表した体積素は
 
{{Indent|
 
<math>dV =dx\, dy\, dz =r^2 \sin\theta\, dr\, d\theta\, d\phi</math>
 
}}
 
となる。また、線素の二乗は
 
{{Indent|
 
<math>ds^2 =dx^2 +dy^2 +dz^2 =dr^2 +r^2d\theta^2 +r^2\sin^2\theta\, d\phi^2</math>
 
}}
 
となる。交叉項が現れないため、球座標は各点において動径が増減する方向と二つの角度が増減する方向がそれぞれに直交している[[直交座標系]]である。
 
 
 
== ベクトル解析 ==
 
球面座標 {{math|(''r'',''&theta;'',''&phi;'')}} での[[位置ベクトル]] {{mvar|'''x'''}} の偏微分により
 
{{Indent|
 
<math>\boldsymbol{e}_r =\frac{\partial\boldsymbol{x}}{\partial r},~
 
\boldsymbol{e}_\theta =\frac{1}{r}\frac{\partial\boldsymbol{x}}{\partial\theta},~
 
\boldsymbol{e}_\phi =\frac{1}{r\sin\theta}\frac{\partial\boldsymbol{x}}{\partial\phi}</math>
 
}}
 
を定義する。
 
標準基底 {{mvar|'''e'''{{sub|x}} , '''e'''{{sub|y}} , '''e'''{{sub|z}}}} を用いれば、位置ベクトルの微分は
 
{{Indent|
 
<math>\begin{align}
 
d\boldsymbol{x} &=\boldsymbol{e}_x\, dx +\boldsymbol{e}_y\, dy +\boldsymbol{e}_z\, dz \\
 
&=\begin{pmatrix}
 
\boldsymbol{e}_x & \boldsymbol{e}_y & \boldsymbol{e}_z \\
 
\end{pmatrix}
 
\begin{pmatrix}
 
dx \\ dy \\ dz \\
 
\end{pmatrix} \\
 
&=\begin{pmatrix}
 
\boldsymbol{e}_x & \boldsymbol{e}_y & \boldsymbol{e}_z \\
 
\end{pmatrix}
 
\frac{\partial(x,y,z)}{\partial(r,\theta,\phi)}
 
\begin{pmatrix}
 
dr \\ d\theta \\ d\phi \\
 
\end{pmatrix}
 
\end{align}</math>
 
}}
 
となるので、具体的に
 
{{Indent|
 
<math>\begin{cases}
 
\boldsymbol{e}_r =\boldsymbol{e}_x \sin\theta\, \cos\phi +\boldsymbol{e}_y \sin\theta\, \sin\phi +\boldsymbol{e}_z \cos\theta \\
 
\boldsymbol{e}_\theta =\boldsymbol{e}_x \cos\theta\, \cos\phi +\boldsymbol{e}_y \cos\theta\, \sin\phi -\boldsymbol{e}_z \sin\theta \\
 
\boldsymbol{e}_\phi =-\boldsymbol{e}_x \sin\phi +\boldsymbol{e}_y \cos\phi \\
 
\end{cases}</math>
 
}}
 
で表される。
 
 
 
標準内積を考えれば
 
{{Indent|
 
<math>| \boldsymbol{e}_r |^2 =| \boldsymbol{e}_\theta |^2
 
=| \boldsymbol{e}_\phi |^2 =1</math>
 
}}
 
{{Indent|
 
<math>\boldsymbol{e}_r\cdot \boldsymbol{e}_\theta
 
=\boldsymbol{e}_r\cdot \boldsymbol{e}_\phi
 
=\boldsymbol{e}_\theta\cdot \boldsymbol{e}_\phi =0</math>
 
}}
 
となり、これらは正規直交基底である。
 
任意のベクトル場 {{mvar|'''A'''}} は
 
{{Indent|
 
<math>A_r =\boldsymbol{e}_r \cdot \boldsymbol{A},~
 
A_\theta =\boldsymbol{e}_\theta \cdot \boldsymbol{A},~
 
A_\phi =\boldsymbol{e}_\phi \cdot \boldsymbol{A}</math>
 
}}
 
{{Indent|
 
<math>\boldsymbol{A} =A_r \boldsymbol{e}_r +A_\theta \boldsymbol{e}_\theta +A_\phi \boldsymbol{e}_\phi</math>
 
}}
 
によって成分表示される。
 
ベクトル場の球面座標による微分は
 
{{Indent|
 
<math>\frac{\partial\boldsymbol{A}}{\partial r}
 
=\frac{\partial A_r}{\partial r} \boldsymbol{e}_r
 
+\frac{\partial A_\theta}{\partial r} \boldsymbol{e}_\theta
 
+\frac{\partial A_\phi}{\partial r} \boldsymbol{e}_\phi</math>
 
}}
 
{{Indent|
 
<math>\frac{\partial\boldsymbol{A}}{\partial\theta}
 
=\left( \frac{\partial A_r}{\partial\theta} -A_\theta \right) \boldsymbol{e}_r
 
+\left( \frac{\partial A_\theta}{\partial\theta} +A_r\right) \boldsymbol{e}_\theta
 
+\frac{\partial A_\phi}{\partial\theta} \boldsymbol{e}_\phi</math>
 
}}
 
{{Indent|
 
<math>\frac{\partial\boldsymbol{A}}{\partial\phi}
 
=\left( \frac{\partial A_r}{\partial\phi}
 
-A_\phi\sin\theta \right) \boldsymbol{e}_r
 
+\left( \frac{\partial A_\theta}{\partial\phi}
 
-A_\phi\cos\theta \right) \boldsymbol{e}_\theta
 
+\left( \frac{\partial A_\phi}{\partial\phi}
 
+A_r\sin\theta +A_\theta\cos\theta \right) \boldsymbol{e}_\phi
 
</math>
 
}}
 
で与えられる。
 
 
 
=== スカラー場の勾配 ===
 
スカラー場 {{math|''f''('''''x''''')}} の[[勾配 (ベクトル解析)|勾配]]は
 
{{Indent|
 
<math>df =(\mathrm{grad}\, f)\cdot d\boldsymbol{x}</math>
 
}}
 
で定義されるベクトル場である。球面座標で表した位置ベクトルの微分が
 
{{Indent|
 
<math>d\boldsymbol{x} =\boldsymbol{e}_r\, dr +r\boldsymbol{e}_\theta\, d\theta +r\sin\theta\, \boldsymbol{e}_\phi\, d\phi </math>
 
}}
 
であることから、球面座標系でのスカラー場 {{mvar|f}} の勾配は
 
{{Indent|
 
<math>\mathrm{grad}\, f =\boldsymbol{e}_r \frac{\partial f}{\partial r}
 
+\frac{\boldsymbol{e}_\theta}{r} \frac{\partial f}{\partial\theta}
 
+\frac{\boldsymbol{e}_\phi}{r\sin\theta} \frac{\partial f}{\partial\phi}</math>
 
}}
 
となる。ベクトル微分演算子を
 
{{Indent|
 
<math>\nabla =\boldsymbol{e}_r \frac{\partial}{\partial r}
 
+\frac{\boldsymbol{e}_\theta}{r} \frac{\partial}{\partial\theta}
 
+\frac{\boldsymbol{e}_\phi}{r\sin\theta} \frac{\partial}{\partial\phi}</math>
 
}}
 
で定めれば
 
{{Indent|
 
<math>\mathrm{grad}\, f =\nabla f</math>
 
}}
 
と書ける。
 
 
 
=== ベクトル場の発散 ===
 
球面座標系でのベクトル場 {{mvar|'''A'''}} の[[発散 (ベクトル解析)|発散]]は
 
{{Indent|
 
<math>\begin{align}
 
\mathrm{div}\, \boldsymbol{A} &= \nabla\cdot \boldsymbol{A}
 
=\boldsymbol{e}_r\cdot\frac{\partial\boldsymbol{A}}{\partial r}
 
+\frac{\boldsymbol{e}_\theta}{r} \cdot\frac{\partial\boldsymbol{A}}{\partial\theta}
 
+\frac{\boldsymbol{e}_\phi}{r\sin\theta} \cdot\frac{\partial\boldsymbol{A}}{\partial\phi} \\
 
&= \left( \frac{\partial A_r}{\partial r} +\frac{2}{r} A_r\right)
 
+\frac{1}{r} \left( \frac{\partial A_\theta}{\partial\theta} +A_\theta\cot\theta \right)
 
+\frac{1}{r\sin\theta} \frac{\partial A_\phi}{\partial\phi} \\
 
&=\frac{1}{r^2} \frac{\partial(r^2A_r)}{\partial r}
 
+\frac{1}{r\sin\theta} \frac{\partial(A_\theta\sin\theta)}{\partial\theta}
 
+\frac{1}{r\sin\theta} \frac{\partial A_\phi}{\partial\phi} \\
 
\end{align}</math>
 
}}
 
となる。
 
 
 
=== ベクトル場の回転 ===
 
球面座標系でのベクトル場 {{mvar|'''A'''}} の[[回転 (ベクトル解析)|回転]]は
 
{{Indent|
 
<math>\begin{align}
 
\mathrm{rot}\, \boldsymbol{A} &=\nabla\times\boldsymbol{A}
 
=\boldsymbol{e}_r\times\frac{\partial\boldsymbol{A}}{\partial r}
 
+\frac{\boldsymbol{e}_\theta}{r} \times\frac{\partial\boldsymbol{A}}{\partial\theta}
 
+\frac{\boldsymbol{e}_\phi}{r\sin\theta} \times\frac{\partial\boldsymbol{A}}{\partial\phi} \\
 
&=\frac{\boldsymbol{e}_r}{r\sin\theta} \left[
 
\frac{\partial(A_\phi\sin\theta)}{\partial\theta} -\frac{\partial A_\theta}{\partial\phi} \right]
 
+\frac{\boldsymbol{e}_\theta}{r} \left[
 
\frac{1}{\sin\theta}\frac{\partial A_r}{\partial\phi}
 
-\frac{\partial(rA_\phi)}{\partial r} \right]
 
+\frac{\boldsymbol{e}_\phi}{r} \left[
 
\frac{\partial(rA_\theta)}{\partial r}
 
-\frac{\partial A_r}{\partial\theta} \right] \\
 
\end{align}</math>
 
}}
 
となる。
 
 
 
 
 
== 関連項目 ==
 
* [[天球座標系]]
 
* [[極座標系]]
 
* [[球面]]
 
 
 
{{DEFAULTSORT:きゆうめんさひようけい}}
 
[[Category:座標]]
 
[[Category:数学に関する記事]]
 
{{Math-stub}}
 
 
 
[[fi:Koordinaatisto#Pallokoordinaatisto]]
 
[[it:Sistema di riferimento#Il sistema sferico]]
 

2018/10/29/ (月) 00:07時点における最新版



楽天市場検索: