垂足座標

垂足座標（すいそくざひょう、英: pedal coordinates）は、ユークリッド平面上の点と曲線に対して定義される座標である。具体的には、点（pedal point）Oと曲線Cに対して、C上の点Pを、PとOの距離rと、Pにおける接線とOの距離p（垂直距離）の組(r, p)によって表す座標である。Oとの距離rと、Pにおける法線とOの距離p_cの組で表す座標（contrapedal coordinate）もよく使われる。ただし、${\textstyle p_{c}:={\sqrt {r^{2}-p^{2}}}}$である。曲線が垂足座標によって表現された式をPedal equationという。

いくつかの曲線はPedal equationによって単純に表すことができる。また、Pedal equationによる形式は曲線の曲率などの計算を簡単にする。垂足座標は古典力学や天体力学において、力の問題を解決するのに使うことができる。

方程式

直交座標

Cを直交座標でf(x, y) = 0と与え、Oを原点とする。点(x, y)の垂足座標は次のように計算できる。

$r={\sqrt {x^{2}+y^{2}}}$

$p={\frac {x{\frac {\partial f}{\partial x}}+y{\frac {\partial f}{\partial y}}}{\sqrt {\left({\frac {\partial f}{\partial x}}\right)^{2}+\left({\frac {\partial f}{\partial y}}\right)^{2}}}}.$

CのPedal equationはこれらの式からx,yを消去することによって得られる。

新たな変数zを用いて、曲線を斎次多項式g(x, y, z) = 0で表したとき、pは次のようにより単純な形で書ける。

$p={\frac {\frac {\partial g}{\partial z}}{\sqrt {\left({\frac {\partial g}{\partial x}}\right)^{2}+\left({\frac {\partial g}{\partial y}}\right)^{2}}}}$

z = 1とすればもとの表示を得る。

極座標

Cを極方程式r = f(θ)で与えたとき、偏接角を$\phi$として

$p=r\sin \phi$

$r={\frac {dr}{d\theta }}\tan \phi .$

と計算できる。CのPedal equationはこれらの式からθを消去することによって得られる。

contrapedal coordinatesへの変換は$p_{c}:={\sqrt {r^{2}-p^{2}}}$として、

$\left|{\frac {dr}{d\theta }}\right|={\frac {rp_{c}}{p}},$

が成り立つことより行える。これは、自励微分方程式によって極方程式が、

$f\left(r,\left|{\frac {dr}{d\theta }}\right|\right)=0,$

と表せるならば、pedal equationは次の形となることを意味する。

$f\left(r,{\frac {rp_{c}}{p}}\right)=0.$

例

例として対数螺旋を挙げる。

$r=ae^{{\frac {\cos \alpha }{\sin \alpha }}\theta }.$

$\theta$ における微分によって次の式を得る。

${\frac {dr}{d\theta }}={\frac {\cos \alpha }{\sin \alpha }}ae^{{\frac {\cos \alpha }{\sin \alpha }}\theta }={\frac {\cos \alpha }{\sin \alpha }}r,$

したがって、

$\left|{\frac {dr}{d\theta }}\right|=\left|{\frac {\cos \alpha }{\sin \alpha }}\right|r,$

が成立する。また、垂足座標において

${\frac {r}{p}}p_{c}=\left|{\frac {\cos \alpha }{\sin \alpha }}\right|r,\qquad \Rightarrow \qquad |\sin \alpha |p_{c}=|\cos \alpha |p,$

を得る。または、$p_{c}^{2}=r^{2}-p^{2}$を用いて、

$p=|\sin \alpha |r.$

この方法は、極方程式におけるn階（$n\geq 1$）の自励な微分方程式の解となる曲線に一般化できる。

$f\left(r,|r'_{\theta }|,r''_{\theta },|r'''_{\theta }|\dots ,r_{\theta }^{(2j)},|r_{\theta }^{(2j+1)}|,\dots ,r_{\theta }^{(n)}\right)=0,$

は与えられた曲線の垂足曲線で、垂足座標では次の式で表現される。

$f(p,p_{c},p_{c}p_{c}',p_{c}(p_{c}p_{c}')',\dots ,(p_{c}\partial _{p})^{n}p)=0,$

ただし、微分の変数は$p$。

力の問題

古典力学における力の問題は、垂足座標を用いれば簡単に解決することができる。

力学系を次のように考える。

${\ddot {x}}=F^{\prime }(|x|^{2})x+2G^{\prime }(|x|^{2}){\dot {x}}^{\perp },$

ただし、中心が$F$であるローレンツ力$G$が発生している平面において、位置と速度がそれぞれ$x,{\dot {x}}$の試験粒子の発展を述べている。

次の量はこの系において保存されている。

$L=x\cdot {\dot {x}}^{\perp }+G(|x|^{2}),\qquad c=|{\dot {x}}|^{2}-F(|x|^{2}),$

このとき$x$の軌跡は、原点をpedal pointとする垂足座標において、次のように与えられる。

${\frac {\left(L-G(r^{2})\right)^{2}}{p^{2}}}=F(r^{2})+c,$

この形式は2017年、ペトル・ブラシュケ（Petr Blaschke）によって発見された。

例

ケプラー問題を考える。

${\ddot {x}}=-{\frac {M}{|x|^{3}}}x,$

垂足座標によって、即座に解にたどり着くことができる。

${\frac {L^{2}}{2p^{2}}}={\frac {M}{r}}+c,$

ここで$L$は粒子の角運動量、$c$はそのエネルギー。また、等式は垂足座標において円錐曲線を表すことが分かった。

逆に試験粒子が曲線C上を動くために与える必要のある力を演繹することもできる。

特別な曲線

正弦波螺旋

正弦波螺旋は次の式で書かれる曲線である。

$r^{n}=a^{n}\sin(n\theta )$

偏接角は

$\psi =n\theta$

である。したがってpedal equationは、

$pa^{n}=r^{n+1}.$

となる。正弦波曲線の仲間のpedal equationを下の表に示した。

n	曲線	垂足点	Pedal eq.
任意	半径aの円	中心	$pa^{n}=r^{n+1}$
1	直径aの円	円周上の点	pa = r2
−1	直線	直線からa離れた点	p = a
1⁄2	カージオイド	尖点	p2a = r3
−1⁄2	放物線	焦点	p2 = ar
2	ベルヌーイのレムニスケート	中心	pa2 = r3
−2	直角双曲線	中心	rp = a2

渦巻

$r=c\theta ^{\alpha },$

の形で与えられる渦巻は次の式を満足する。

${\frac {dr}{d\theta }}=\alpha r^{\frac {\alpha -1}{\alpha }},$

したがって、その垂足座標における表示は次のようになる。

${\frac {1}{p^{2}}}={\frac {\alpha ^{2}c^{\frac {2}{\alpha }}}{r^{2+{\frac {2}{\alpha }}}}}+{\frac {1}{r^{2}}}.$

下の表は特別の場合。

$\alpha$	曲線	Pedal point	Pedal eq.
1	アルキメデスの螺旋	原点	${\frac {1}{p^{2}}}={\frac {1}{r^{2}}}+{\frac {c^{2}}{r^{4}}}$
−1	双曲螺旋	原点	${\frac {1}{p^{2}}}={\frac {1}{r^{2}}}+{\frac {1}{c^{2}}}$
1⁄2	放物螺旋	原点	${\frac {1}{p^{2}}}={\frac {1}{r^{2}}}+{\frac {c^{4}}{4r^{6}}}$
−1⁄2	リチュース	原点	${\frac {1}{p^{2}}}={\frac {1}{r^{2}}}+{\frac {r^{2}}{4c^{4}}}$

サイクロイド

エピサイクロイドまたはハイポサイクロイドは次のパラメトリック方程式で表せる。

$x(\theta )=(a+b)\cos \theta -b\cos \left({\frac {a+b}{b}}\theta \right)$

$y(\theta )=(a+b)\sin \theta -b\sin \left({\frac {a+b}{b}}\theta \right),$

原点におけるpedal equationは、

$r^{2}=a^{2}+{\frac {4(a+b)b}{(a+2b)^{2}}}p^{2}$

または、

$p^{2}=A(r^{2}-a^{2})$

である。ただし、

$A={\frac {(a+2b)^{2}}{4(a+b)b}}.$

特別の場合を以下の表に示した。b=a⁄nとする。

n	曲線	Pedal eq.
1, −1⁄2	カージオイド	$p^{2}={\frac {9}{8}}(r^{2}-a^{2})$
2, −2⁄3	ネフロイド	$p^{2}={\frac {4}{3}}(r^{2}-a^{2})$
−3, −3⁄2	デルトイド	$p^{2}=-{\frac {1}{8}}(r^{2}-a^{2})$
−4, −4⁄3	アステロイド	$p^{2}=-{\frac {1}{3}}(r^{2}-a^{2})$

他の曲線

他の有名な曲線においては下の表にまとめる。

曲線	方程式	Pedal point	Pedal eq.
直線	$ax+by+c=0$	原点	$p={\frac {|c|}{\sqrt {a^{2}+b^{2}}}}$
点	$(x_{0},y_{0})$	原点	$r={\sqrt {x_{0}^{2}+y_{0}^{2}}}$
円	$|x-a|=R$	原点	$2pR=r^{2}+R^{2}-|a|^{2}$
伸開線	$r={\frac {a}{\cos \alpha }},\ \theta =\tan \alpha -\alpha$	原点	$p_{c}=|a|$
楕円	${\frac {x^{2}}{a^{2}}}+{\frac {y^{2}}{b^{2}}}=1$	中心	${\frac {a^{2}b^{2}}{p^{2}}}+r^{2}=a^{2}+b^{2}$
双曲線	${\frac {x^{2}}{a^{2}}}-{\frac {y^{2}}{b^{2}}}=1$	中心	$-{\frac {a^{2}b^{2}}{p^{2}}}+r^{2}=a^{2}-b^{2}$
楕円	${\frac {x^{2}}{a^{2}}}+{\frac {y^{2}}{b^{2}}}=1$	焦点	${\frac {b^{2}}{p^{2}}}={\frac {2a}{r}}-1$
双曲線	${\frac {x^{2}}{a^{2}}}-{\frac {y^{2}}{b^{2}}}=1$	焦点	${\frac {b^{2}}{p^{2}}}={\frac {2a}{r}}+1$
対数螺旋	$r=ae^{\theta \cot \alpha }$	極	$p=r\sin \alpha$
デカルトの卵形	$|x|+\alpha |x-a|=C,$	焦点	$b:=C^{2}-\alpha ^{2}|a|^{2}$として、 ${\frac {(b-(1-\alpha ^{2})r^{2})^{2}}{4p^{2}}}={\frac {Cb}{r}}+(1-\alpha ^{2})Cr-((1-\alpha ^{2})C^{2}+b)$
カッシーニの卵形線	$|x||x-a|=C,$	焦点	${\frac {(3C^{2}+r^{4}-|a|^{2}r^{2})^{2}}{p^{2}}}=4C^{2}\left({\frac {2C^{2}}{r^{2}}}+2r^{2}-|a|^{2}\right).$
カッシーニの卵形線	$|x-a||x+a|=C,$	中心	$2Rpr=r^{4}+R^{2}-|a|^{2}.$

関連項目

• 垂足曲線

出典

• R.C. Yates (1952). “Pedal Equations”. A Handbook on Curves and Their Properties. Ann Arbor, MI: J. W. Edwards. pp. 166 ff. 
• J. Edwards (1892). Differential Calculus. London: MacMillan and Co.. pp. 161 ff.. https://archive.org/details/in.ernet.dli.2015.109607 

• P. Blaschke (2017). “Pedal coordinates, dark Kepler and other force problems”. Journal of Mathematical Physics 58/6. arXiv:1704.00897. doi:10.1063/1.4984905. https://zenodo.org/record/897629/files/article.pdf. 
• Petr, Blaschke. "Pedal coordinates, solar sail orbits, Dipole drive and other force problems". arXiv:2005.05767。

外部リンク

• Weisstein, Eric W. "Pedal coordinates". mathworld.wolfram.com (英語).