円柱座標系 (えんちゅうざひょうけい、英 : cylindrical coordinate system ; 円筒座標系 とも)は三次元の座標系 であって、点 の位置を
特別に選ばれた基準軸からの距離、
特別に選ばれた基準方向に対する軸から測った方向、
基準軸に直交する特別に選ばれた基準平面からの距離
A cylindrical coordinate system with origin O , polar axis A , and longitudinal axis L . The dot is the point with radial distance ρ = 4 , angular coordinate φ = 130° , and height z = 4 .
の三者によって決定する。ただし基準平面からの「距離」はその点が基準平面の(表または裏の)どちら側に面するかによって正または負の値を持つものとする。
円筒座標系の原点 (origin ) とは、上記三つの座標成分がすべて 0 として与えることができるような点を言う。これは上記の基準平面と基準軸の交点になる。
この基準軸は円筒軸 (cylindrical axis ; 円柱軸) や緯線 (longitudinal axis ) などと呼ばれ、極軸 (polar axis ) とは区別される(極軸あるいは始線とは、基準平面上に載っている半直線 で、原点から出てその平面上の基準方向を指し示すものを表すために用いられる)。緯線に垂直な任意の方向を射線 (radial lines ; 放射状の直線) と総称する。
円筒軸からの距離は動径距離 (radial distance ; 放射距離) や動径 (radius ) と呼ばれ、円筒軸周りの偏角座標は角度位置 (angular position ) や方位角 (azimuth ) などと言う。動径成分と方位角成分を併せて極座標成分 (polar coordinates ) と言い、これはその点を通り基準平面に平行な平面上の二次元の極座標系 に対応する。第三座標は基準平面を水平面と見るとき高さ (height ) や高度 (altitude ) と呼んだり、緯度 (longitudinal position [1] ) や軸位置 (axial position [2] ) などとも言う。
円筒座標系は緯線周りの何らかの回転対称性 を持つ物体や現象(例えば、丸い断面を持つ直線パイプを流れる水流や、金属円柱 の熱分布、長い真っ直ぐなワイヤー内の電荷 から出る電場 、天文学における降着円盤 など)との関連で有意である。
円筒座標系は「円筒極座標系」("cylindrical polar coordinates"[3] ) とか「極円筒座標系」("polar cylindrical coordinates"[4] ) などと呼ばれたりもする。銀河の星の位置を指定するために用いられることもある(「銀河中心的円柱極座標系」("galactocentric cylindrical polar coordinates"[5] ))。
与えられた点 P の三つの座標成分 (ρ , φ , z ) は以下のように与えられる:
軸距離 (axial distance ) または動径距離 (radial distance ) ρ は z -軸から点 P までのユークリッド距離 ;
方位角 (azimuth ) φ は基準平面上の基準方向と原点から点 P の基準平面への正射影へ結んだ直線との間の角;
軸座標 (axial coordinate ) または高さ (height ) z は基準平面から点 P までの符号付き距離。
一意な円筒座標成分
極座標系 の場合と同様に、円筒座標 (ρ , φ , z ) の点と同じ点を表す同値な座標成分は無数に存在する— 具体的に書けば、n を整数 として (ρ , φ + n ×360°, z ) および (−ρ , φ + (2n + 1)×180°, z ) は同じ点を表す— (さらに言えば、動径成分 ρ が 0 ならば、方位角は任意にとれる)。
各点に対して座標が一意に定まることが必要となるような場面では、例えば、動径成分は非負 (ρ ≥ 0 ) で、方位角 φ は360°間隔の特定区間 (例えば [ −180°, +180°] や [ 0, 360°] など)に入るなどと約束すればよい。
表記の慣習
円筒座標系の記法は一通りではない。ISO 標準規格31-11 では (ρ , φ , z ) は ρ が動径座標、φ が方位角、z が高さとなるように求めている。しかし動径成分はしばしば r や s と書かれるし、方位角が θ や t と書くこともよくある。もう一つの成分も、高さを h で表したり、円筒軸が水平線と見なされるときには x と書くなど、文脈に依存した特定の文字が任意に用いられる。
The coordinate surfaces of the cylindrical coordinates (ρ , φ , z ) . The red cylinder shows the points with ρ = 2 , the blue plane shows the points with z = 1 , and the yellow half-plane shows the points with φ = −60° . The z -axis is vertical and the x -axis is highlighted in green. The three surfaces intersect at the point P with those coordinates (shown as a black sphere); the Cartesian coordinates of P are roughly (1.0, −1.732, 1.0).
Cylindrical coordinate surfaces. The three orthogonal components, ρ (green), φ (red), and z (blue), each increasing at a constant rate. The point is at the intersection between the three colored surfaces.
具体的な状況において、そして多くの数学的図式化において、偏角成分の正の方向は正の高さを持つ任意の点から見て反時計回り に測るものとする。
円筒座標系は数ある三次元座標系の一つであり、ほかの三次元座標系と以下に述べるような変換公式のもとで読み替えることができる。
直交座標系
円筒座標系と直交座標系との間の変換は、円筒座標系の基準面を xy -直交座標平面(方程式 z = 0 の表す平面)と見なし、かつ円筒軸を z -直交座標軸と見なすと容易に扱える。こうすれば、両座標系で z -座標は共通になり、円筒座標 (ρ ,φ ,z ) と直交座標 (x ,y ,z ) の間の変換則は二次元の極座標系の場合と同じく
{
x
=
ρ
cos
φ
y
=
ρ
sin
φ
z
=
z
{\displaystyle {\begin{cases}x=\rho \cos \varphi \\y=\rho \sin \varphi \\z=z\end{cases}}}
および
ρ
=
x
2
+
y
2
φ
=
{
0
if
x
=
0
and
y
=
0
arcsin
y
ρ
if
x
≥
0
arctan
y
x
if
x
>
0
−
arcsin
y
ρ
+
π
if
x
<
0
{\displaystyle {\begin{aligned}\rho &={\sqrt {x^{2}+y^{2}}}\\\varphi &={\begin{cases}0&{\text{if }}x=0{\text{ and }}y=0\\\arcsin {\dfrac {y}{\rho }}&{\text{if }}x\geq 0\\\arctan {\dfrac {y}{x}}&{\text{if }}x>0\\-\arcsin {\dfrac {y}{\rho }}+\pi &{\text{if }}x<0\end{cases}}\end{aligned}}}
で相互に読み替え可能である。ここで逆正弦函数 arcsin は [ −π/2, +π/2] = [ −90°, +90°] の範囲の角を返す枝をとるものとする。これらの公式から、方位角 φ は [ −90°, +270°] の範囲になることがわかる。(他の変換公式は極座標系 の項を参照のこと。)
多くのモダンなプログラミング言語には、与えられた x, y から(上記のような場合分けをすることなく)(−π, π) の範囲での正しい方位角 φ を計算する函数が提供されている。[注釈 1]
球面座標系
球面座標系 (半径 r , 仰角 (elevation) または傾斜角 (inclination) θ , 方位角 φ ) の円筒座標系への変換は
で、円筒座標系から球面座標系への変換は
で得られる。
円筒極座標系を含む多くの問題において線素および体積要素がこの座標系でどのように表されるかを知っていることは有意である(それにより曲線経路や体積を含む問題を積分によって解く方法が考えられるようになる)。
円筒座標系における線素 (英語版 ) は
d
r
:=
d
ρ
ρ
^
+
ρ
d
φ
φ
^
+
d
z
z
^
{\textstyle {\mathit {dr}}:={\mathit {d\rho }}\,{\hat {\boldsymbol {\rho }}}+\rho \,{\mathit {d\varphi }}\,{\hat {\boldsymbol {\varphi }}}+{\mathit {dz}}\,{\hat {\boldsymbol {z}}}}
で与えられる。
円筒座標系における体積要素 は
d
V
:=
ρ
d
ρ
d
φ
d
z
{\textstyle {\mathit {dV}}:=\rho \,{\mathit {d\rho }}\,{\mathit {d\varphi }}\,{\mathit {dz}}}
となる。
動径成分 ρ が一定となる曲面(垂直円筒)上の面素(面積要素)は
d
S
ρ
:=
ρ
d
φ
d
z
{\textstyle {\mathit {dS}}_{\rho }:=\rho \,{\mathit {d\varphi }}\,{\mathit {dz}}}
になる。
方位角 φ が一定の曲面(垂直半平面)上の面素は
d
S
φ
:=
d
ρ
d
z
{\textstyle {\mathit {dS}}_{\varphi }:={\mathit {d\rho }}\,{\mathit {dz}}}
である。
高さ z が一定の曲面(水平面)上の面素は
d
S
z
:=
ρ
d
ρ
d
φ
{\textstyle {\mathit {dS}}_{z}:=\rho \,{\mathit {d\rho }}\,{\mathit {d\varphi }}}
で得られる。
微分作用素のナブラ を円筒座標系で表現すれば、以下のように勾配 ・発散 ・回転 およびラプラス作用素 はそれぞれ、
∇
f
=
∂
f
∂
ρ
ρ
^
+
1
ρ
∂
f
∂
φ
φ
^
+
∂
f
∂
z
z
^
∇
⋅
A
=
1
ρ
∂
∂
ρ
(
ρ
A
ρ
)
+
1
ρ
∂
A
φ
∂
φ
+
∂
A
z
∂
z
∇
×
A
=
(
1
ρ
∂
A
z
∂
φ
−
∂
A
φ
∂
z
)
ρ
^
+
(
∂
A
ρ
∂
z
−
∂
A
z
∂
ρ
)
φ
^
+
1
ρ
(
∂
∂
ρ
(
ρ
A
φ
)
−
∂
A
ρ
∂
φ
)
z
^
∇
2
f
=
1
ρ
∂
∂
ρ
(
ρ
∂
f
∂
ρ
)
+
1
ρ
2
∂
2
f
∂
φ
2
+
∂
2
f
∂
z
2
{\displaystyle {\begin{aligned}\nabla f&={\frac {\partial f}{\partial \rho }}{\hat {\boldsymbol {\rho }}}+{\frac {1}{\rho }}{\frac {\partial f}{\partial \varphi }}{\hat {\boldsymbol {\varphi }}}+{\frac {\partial f}{\partial z}}{\hat {\boldsymbol {z}}}\\[8px]\nabla \cdot {\boldsymbol {A}}&={\frac {1}{\rho }}{\frac {\partial }{\partial \rho }}\left(\rho A_{\rho }\right)+{\frac {1}{\rho }}{\frac {\partial A_{\varphi }}{\partial \varphi }}+{\frac {\partial A_{z}}{\partial z}}\\[8px]\nabla \times {\boldsymbol {A}}&=\left({\frac {1}{\rho }}{\frac {\partial A_{z}}{\partial \varphi }}-{\frac {\partial A_{\varphi }}{\partial z}}\right){\hat {\boldsymbol {\rho }}}+\left({\frac {\partial A_{\rho }}{\partial z}}-{\frac {\partial A_{z}}{\partial \rho }}\right){\hat {\boldsymbol {\varphi }}}+{\frac {1}{\rho }}\left({\frac {\partial }{\partial \rho }}\left(\rho A_{\varphi }\right)-{\frac {\partial A_{\rho }}{\partial \varphi }}\right){\hat {\boldsymbol {z}}}\\[8px]\nabla ^{2}f&={\frac {1}{\rho }}{\frac {\partial }{\partial \rho }}\left(\rho {\frac {\partial f}{\partial \rho }}\right)+{\frac {1}{\rho ^{2}}}{\frac {\partial ^{2}f}{\partial \varphi ^{2}}}+{\frac {\partial ^{2}f}{\partial z^{2}}}\end{aligned}}}
と書くことができる。
円筒対称性を持つ系におけるラプラス方程式 の解は円筒調和函数 (英語版 ) と言う。
注釈
For example, this function is called by atan2 (y ,x ) in the C programming language, and atan(y ,x ) in Common Lisp .
Morse, Philip M. ; Feshbach, Herman (1953). Methods of Theoretical Physics, Part I . New York City : McGraw-Hill . pp. 656–657. ISBN 0-07-043316-X . LCCN 52-11515
Margenau, Henry ; Murphy, George M. (1956). The Mathematics of Physics and Chemistry . New York City: D. van Nostrand. p. 178. ISBN 9780882754239 . LCCN 55-10911 . OCLC 3017486
Korn, Granino A.; Korn, Theresa M. (1961). Mathematical Handbook for Scientists and Engineers . New York City: McGraw-Hill. pp. 174–175. LCCN 59-14456 . ASIN B0000CKZX7
Sauer, Robert; Szabó, István (1967). Mathematische Hilfsmittel des Ingenieurs . New York City: Springer-Verlag . p. 95. LCCN 67-25285
Zwillinger, Daniel (1992). Handbook of Integration . Boston : Jones and Bartlett Publishers . p. 113. ISBN 0-86720-293-9 . OCLC 25710023
Moon, P.; Spencer, D. E. (1988). “Circular-Cylinder Coordinates (r, ψ, z)”. Field Theory Handbook, Including Coordinate Systems, Differential Equations, and Their Solutions (corrected 2nd ed.). New York City: Springer-Verlag. pp. 12–17, Table 1.02. ISBN 978-0-387-18430-2
List of canonical coordinate transformations
Vector fields in cylindrical and spherical coordinates
Del in cylindrical and spherical coordinates