物理/惑星の運動(2)
提供: Internet Web School
(→ 2体の重心は太陽の重心と一致するという近似を行う) |
(→ 微分方程式を解いて、惑星の軌道形を求める ) |
||
(間の15版分が非表示) | |||
13 行: | 13 行: | ||
この系には、系外から力が作用しないため、、<br/> | この系には、系外から力が作用しないため、、<br/> | ||
[[物理/質点の運動と質点系#質点系の重心とその運動|系の重心は等速直線運動をする。]] <br/> | [[物理/質点の運動と質点系#質点系の重心とその運動|系の重心は等速直線運動をする。]] <br/> | ||
- | 従って、この重心を原点とする回転しない座標系も慣性系となる。 | + | 従って、この重心を原点とする回転しない座標系も慣性系となる。<br/> |
- | + | 今後はこの慣性系で観測する。<br/><br/> | |
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
(注1)慣性系の存在はニュートン力学(古典力学)の大前提であり、運動法則の第一法則としてまとめられている。<br/> | (注1)慣性系の存在はニュートン力学(古典力学)の大前提であり、運動法則の第一法則としてまとめられている。<br/> | ||
太陽系の重心に原点を取り、恒星系に対して回転しない座標系は慣性系に近いことが現在分かっている。 | 太陽系の重心に原点を取り、恒星系に対して回転しない座標系は慣性系に近いことが現在分かっている。 | ||
29 行: | 22 行: | ||
(注2)他の惑星の影響を無視し、太陽と当該惑星の2体を質点とみなしても、<br/> | (注2)他の惑星の影響を無視し、太陽と当該惑星の2体を質点とみなしても、<br/> | ||
他の惑星からの影響の少ない地球などの惑星の運動は、かなり正確に予測できることが知られている。<br/> | 他の惑星からの影響の少ない地球などの惑星の運動は、かなり正確に予測できることが知られている。<br/> | ||
- | 正確にするには二体ではなく他の惑星も考慮した多体問題として論じなければならないが、解析は非常に難しくなる。<br/><br/> | + | 正確にするには二体ではなく他の惑星も考慮した多体問題として論じなければならないが、解析は非常に難しくなる。 |
+ | === 2体の重心は太陽の重心と一致するという近似を行う === | ||
+ | 太陽の質量は惑星の質量より圧倒的大きいため二体の重心は太陽の重心とほぼ一致する。 <br/> | ||
+ | *[[wikipedia_ja:大きさ順の太陽系天体の一覧#一覧 |ウィキペディア(大きさ順の太陽系天体の一覧)]] | ||
+ | <br/> | ||
+ | そこで議論を簡単にするため、このテキストでは2体の重心は太陽の重心と一致し、太陽は不動であると仮定する。<br/> | ||
+ | 興味のある読者は、この仮定をしないで、惑星と太陽の運動を解析してみてほしい。<br/> | ||
+ | 議論は殆ど同じように進められます。<br/><br/> | ||
+ | |||
そこで太陽Sを原点Oとする慣性座標系を考える。<br/> | そこで太陽Sを原点Oとする慣性座標系を考える。<br/> | ||
任意の時刻 t の惑星Ptの位置ベクトル→OPtを→r(t)とかく。<br/> | 任意の時刻 t の惑星Ptの位置ベクトル→OPtを→r(t)とかく。<br/> | ||
37 行: | 38 行: | ||
== 運動の第1、第2法則と万有引力の法則からケプラーの3法則を導く == | == 運動の第1、第2法則と万有引力の法則からケプラーの3法則を導く == | ||
- | |||
- | |||
- | |||
- | |||
- | |||
- | |||
=== 惑星は平面上を運動する === | === 惑星は平面上を運動する === | ||
不動の太陽の位置O,時刻tの惑星の位置P(t0)とP(t0)を始点とした初速ベクトル→v0を含む<br/> | 不動の太陽の位置O,時刻tの惑星の位置P(t0)とP(t0)を始点とした初速ベクトル→v0を含む<br/> | ||
158 行: | 153 行: | ||
惑星に働く万有引力を、この極座標で表示すると<br/> | 惑星に働く万有引力を、この極座標で表示すると<br/> | ||
$-\frac{GM}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|} | $-\frac{GM}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|} | ||
- | =-\frac{GM}{ | + | =-\frac{GM}{r^{2}(t)}(\cos \theta (t),\sin \theta (t)\bigr) \qquad (9)$<br/> |
式(8),(9)を惑星の運動方程式<br/> | 式(8),(9)を惑星の運動方程式<br/> | ||
¨→r(t)=−GM‖<br/> | \ddot{\vec r}(t)=-\frac{GM}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|}<br/> | ||
164 行: | 159 行: | ||
\bigl(\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)\bigr)\bigl(\cos \theta (t),\sin \theta (t)\bigr)<br/> | \bigl(\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)\bigr)\bigl(\cos \theta (t),\sin \theta (t)\bigr)<br/> | ||
+\bigl(2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \bigr)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)<br/> | +\bigl(2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \bigr)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)<br/> | ||
- | $= -\frac{GM}{ | + | $= -\frac{GM}{r^{2}(t)}(\cos \theta (t),\sin \theta (t)\bigr)$<br/> |
右辺を左辺に移行すると<br/> | 右辺を左辺に移行すると<br/> | ||
- | $\bigl(\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)+\frac{GM}{ | + | $\bigl(\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)+\frac{GM}{r^{2}(t)}\bigr)\bigl(\cos \theta (t),\sin \theta (t)\bigr)$<br/> |
+\bigl(2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \bigr)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)\equiv 0<br/> | +\bigl(2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \bigr)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)\equiv 0<br/> | ||
\bigl(\cos \theta (t),\sin \theta (t)\bigr)と\bigl(-\sin \theta (t),\cos \theta (t)\bigr)は直交ベクトルなので、この式から<br/><br/> | \bigl(\cos \theta (t),\sin \theta (t)\bigr)と\bigl(-\sin \theta (t),\cos \theta (t)\bigr)は直交ベクトルなので、この式から<br/><br/> | ||
- | $\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)+\frac{GM}{ | + | $\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)+\frac{GM}{r^{2}(t)}\equiv 0 \qquad \qquad (10) $<br/> |
2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \equiv 0 \qquad \qquad (11) <br/> | 2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \equiv 0 \qquad \qquad (11) <br/> | ||
が得られる。<br/><br/> | が得られる。<br/><br/> | ||
175 行: | 170 行: | ||
太陽の位置Oを原点、x軸の正の部分の半直線を角の基準線とする極座標で | 太陽の位置Oを原点、x軸の正の部分の半直線を角の基準線とする極座標で | ||
惑星の運動方程式を表現すると<br/> | 惑星の運動方程式を表現すると<br/> | ||
- | $\ddot {r}-r\dot{\theta}^{2}=-\frac{GM}{ | + | \ddot {r}-r\dot{\theta}^{2}=-\frac{GM}{r^{2}}\qquad \qquad (12) <br/> |
r^2\dot{\theta} = h \quad (任意の定数) \qquad \qquad (13) <br/> | r^2\dot{\theta} = h \quad (任意の定数) \qquad \qquad (13) <br/> | ||
- | となる。<br/> | + | となる。(注2 を参照のこと)<br/> |
証明<br/> | 証明<br/> | ||
式(12)は式(10)から明らか。<br/> | 式(12)は式(10)から明らか。<br/> | ||
188 行: | 183 行: | ||
証明終わり。 \qquad \qquad \qquad \qquad \Box<br/><br/> | 証明終わり。 \qquad \qquad \qquad \qquad \Box<br/><br/> | ||
(注)式(13)は面積速度が一定であることを表現している。<br/> | (注)式(13)は面積速度が一定であることを表現している。<br/> | ||
- | 先に、このことは証明してあるが、その別証明である。 | + | 先に、このことは証明してあるが、その別証明である。<br/> |
+ | (注2)関数の変数tまで書くと、式(12),(13)は、次のように書ける。<br/> | ||
+ | \ddot {r}(t)-r(t)\dot{\theta}^{2}(t)=-\frac{GM}{r^{2}(t)}\qquad \qquad (12') <br/> | ||
+ | r^2(t)\dot{\theta}(t) = h \quad (任意の定数) \qquad \qquad (13') <br/> | ||
- | ==== | + | ==== 微分方程式を解いて、惑星の軌道形を求める ==== |
+ | 連立微分方程式(12)、(13)を解いて、r=r(t),\theta=\theta(t)を求めれば、<br/> | ||
+ | 任意の時刻tにおける惑星の位置がわかる。<br/> | ||
+ | 軌道の形は、この2式からtを消去して、rとθの関数関係<br/> | ||
+ | f(r,\theta)=0, またはr=f(\theta)<br/> | ||
+ | を求めれな良い。<br/> | ||
+ | しかし、残念なことにこれは大変困難である。<br/><br/> | ||
+ | |||
+ | 軌道の形だけを知りたいときには、<br/> | ||
+ | これを回避して簡単な(しかし技巧的な)直接的な解法がある。<br/> | ||
+ | それは、連立微分方程式の解が<br/> | ||
+ | r=f\circ \theta ,\qquad変数まで書けば、r(t)=f\circ \theta (t)\triangleq f(\theta (t)) <br/> | ||
+ | という合成関数になることを仮定して、<br/> | ||
+ | r=f(\theta)を直接求める問題に変換し、解を求める方法である。<br/><br/> | ||
+ | ===== 軌道の形を直接求めるための微分方程式の導出 ===== | ||
+ | \frac{d r(t) }{dt}=\frac{d f\circ \theta (t) }{dt}<br/> | ||
+ | \qquad 合成関数の微分公式を用いて、<br/> | ||
+ | =\frac{d f}{d\theta}(\theta (t))\dot{\theta}(t)<br/> | ||
+ | |||
+ | \qquad 式(13)を\dot{\theta}=\frac{h}{r^2}に変形し、上式に代入すると<br/> | ||
+ | =\frac{df}{d\theta}(\theta (t))\frac{h}{r^2(t)}<br/> | ||
+ | 故に<br/> | ||
+ | \frac{d r(t) }{dt}=\frac{df}{d\theta}(\theta (t))\frac{h}{r^2(t)} \qquad \qquad (b)<br/> | ||
+ | 補題 g(\theta)\triangleq \frac{1}{f(\theta)}とおくと、<br/> | ||
+ | \frac{df(\theta)}{d\theta}= -\frac{1}{g^2(\theta)}\frac{dg(\theta)}{d\theta} | ||
+ | <br/> | ||
+ | 証明は、 | ||
+ | [[wikibooks_ja:高等学校数学III/微分法#商の導関数|ウィキブックス高等学校数学III/微分法#商の導関数]]をご覧ください。<br/><br/> | ||
+ | この補題を用いて式(b)を変形すると<br/> | ||
+ | \frac{d r(t) }{dt}=-\frac{1}{g^2(\theta (t))}\frac{dg}{d\theta}(\theta (t))\frac{h}{r^2(t)} <br/> | ||
+ | \qquad r(t)=f(\theta (t)) なので、g^2(\theta (t))r^2(t)=1となり、<br/> | ||
+ | =-h\frac{dg}{d\theta}(\theta (t)) <br/> | ||
+ | 故に<br/> | ||
+ | \dot{r}(t)=-h\frac{dg}{d\theta}(\theta (t)) \qquad \qquad (c)<br/> | ||
+ | これをさらに時間tで微分すると、合成関数の微分の公式から<br/> | ||
+ | \ddot{r}(t)=-h\frac{d^{2}g}{d\theta ^2}(\theta(t))\dot{\theta}(t)<br/> | ||
+ | \qquad 式(13)から得られる\dot{\theta}=\frac{h}{r^2}を代入して、<br/> | ||
+ | =-h\frac{d^{2}g}{d\theta ^2}(\theta(t)) \frac{h}{r^2(t)} <br/> | ||
+ | 故に<br/> | ||
+ | \ddot{r}(t)=-\frac{h^{2}}{r^2(t)}\frac{d^{2}g}{d\theta ^2}(\theta(t))\qquad \qquad (d)<br/> | ||
+ | が得られる。<br/><br/> | ||
+ | 式(d)と\dot{\theta}=\frac{h}{r^2}を、式(12)に代入して整頓すると次式が得られる。<br/> | ||
+ | \frac{d^{2}g}{d\theta ^2}(\theta(t))+\frac{1}{r(t)}=\frac{GM}{h^2} <br/> | ||
+ | \frac{1}{r(t)}=\frac{1}{f(\theta (t))}=g(\theta (t))なので<br/> | ||
+ | 上式から次式が得られる。<br/> | ||
+ | \frac{d^{2}g}{d\theta ^2}(\theta(t))+g(\theta (t))=\frac{GM}{h^2} <br/> | ||
+ | 従って変数を \theta とする微分方程式<br/> | ||
+ | \frac{d^{2}g}{d\theta ^2}(\theta)+g(\theta )=\frac{GM}{h^2} \qquad \qquad (14)<br/> | ||
+ | の解g=g(\theta)が分かれば<br/> | ||
+ | r=f(\theta)=\frac{1}{g(\theta)} <br/> | ||
+ | が、惑星の軌道の形を表す式になる。 | ||
+ | ===== 微分方程式の解法 ===== | ||
==== 極座標表示の惑星軌道を、直交座標表示に変換 ==== | ==== 極座標表示の惑星軌道を、直交座標表示に変換 ==== |
2018年6月25日 (月) 06:11 時点における最新版
目次[非表示] |
☆☆惑星の運動(2)
2.3節の惑星運動の項では、
惑星の軌道を円と仮定すると、
運動の第1、第2法則のもとでは、
ケプラーの法則と万有引力の法則が同等であることを示した。
この節では軌道について仮定せず、
運動の第1、第2法則のもとで、
万有引力の法則と惑星に関するケプラーの3法則が同等であることを導く。
ニュートンの運動の第一、第二法則による惑星の運動方程式の導出
他の惑星の影響を無視して、太陽と一つの惑星の二体からなる系をある慣性系から観測する。(注1と注2を参照のこと)
太陽の質量をM、惑星の質量をmとする。
この二体を、ともに質点であるとみなし、位置はそれぞれの重心とする。
この系には、系外から力が作用しないため、、
系の重心は等速直線運動をする。
従って、この重心を原点とする回転しない座標系も慣性系となる。
今後はこの慣性系で観測する。
(注1)慣性系の存在はニュートン力学(古典力学)の大前提であり、運動法則の第一法則としてまとめられている。
太陽系の重心に原点を取り、恒星系に対して回転しない座標系は慣性系に近いことが現在分かっている。
慣性系からみた質点系の運動は、テキストの2.3 で、その導入部分を説明している。
(注2)他の惑星の影響を無視し、太陽と当該惑星の2体を質点とみなしても、
他の惑星からの影響の少ない地球などの惑星の運動は、かなり正確に予測できることが知られている。
正確にするには二体ではなく他の惑星も考慮した多体問題として論じなければならないが、解析は非常に難しくなる。
2体の重心は太陽の重心と一致するという近似を行う
太陽の質量は惑星の質量より圧倒的大きいため二体の重心は太陽の重心とほぼ一致する。
そこで議論を簡単にするため、このテキストでは2体の重心は太陽の重心と一致し、太陽は不動であると仮定する。
興味のある読者は、この仮定をしないで、惑星と太陽の運動を解析してみてほしい。
議論は殆ど同じように進められます。
そこで太陽Sを原点Oとする慣性座標系を考える。
任意の時刻\ t\ の惑星P_tの位置ベクトル\vec{OP_t}を\vec{r}(t)とかく。
惑星が太陽からうける力(万有引力の法則を仮定しないと未知)を仮に\vec fと書くと、
惑星の運動は、ニュートンの第2法則から次の微分方程式で形式的には記述されが、
力が未知なので、この方程式は解を持ち得ない。
\ M\frac{d^{2} \vec{r}}{dt^2}(t)=\vec{f}
運動の第1、第2法則と万有引力の法則からケプラーの3法則を導く
惑星は平面上を運動する
不動の太陽の位置O,時刻tの惑星の位置P(t_0)とP(t_0)を始点とした初速ベクトル\vec{v}_0を含む
平面Hを考える。
この平面と直交する一つのベクトル\vec hを考え、
式(1)の両辺の左側からベクトル積として作用させる。
\ \vec h \times \Bigl(M\frac{d \vec{v}}{dt}(t)\Bigr)=\ \vec h \times \Bigl(-\frac{GMm}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|}\Bigr) \qquad \qquad (3)
ベクトル積の性質から
左辺は
\ \vec h \times \Bigl(M\frac{d \vec{v}}{dt}(t)\Bigr)=\ M\vec h \times \Bigl(\frac{d \vec{v}}{dt}(t)\Bigr)=\ M\vec h \times \Bigl(\frac{d^2 \vec{r}}{dt^2}(t)\Bigr)
=\ M\Bigl(\frac{d^2 \vec h \times\vec{r}}{dt^2}(t)\Bigr)(注1をみよ)
右辺は
\ \vec h \times \Bigl(-\frac{GMm}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|}\Bigr) = -\frac{GMm}{\|\vec{r}(t)\|^{3}}\vec h \times \vec{r}(t)
なので、
\ M\Bigl(\frac{d^2 \vec h \times\vec{r}}{dt^2}(t)\Bigr) = -\frac{GMm}{\|\vec{r}(t)\|^{3}}\vec h \times \vec{r}(t)
が得られる。
故に
\ \Bigl(\frac{d^2 \vec h \times\vec{r}}{dt^2}(t)\Bigr) = -\frac{Gm}{\|\vec{r}(t)\|^{3}}\vec h \times \vec{r}(t)\qquad \qquad (4)
\vec hをベクトル積として初期条件に作用させると
\vec h \times \vec{r}(t_0) = \vec h \times \vec{r_0}=0,
\frac{d \vec h \times \vec{r}}{dt}(t_0) =\vec h \times \frac{d\vec{r}}{dt}(t_0) =\vec h \times \vec{v_0} = 0
故に、
\vec h \times \vec{r}(t_0) = 0,
\frac{d \vec h \times \vec{r}}{dt}(t_0) = 0\qquad \qquad (5)
tの関数\vec{x}\triangleq \vec h \times \vec{r}を考えると、
式(4),(5)は
\ \Bigl(\frac{d^2 \vec{x}}{dt^2}(t)\Bigr) = -\frac{Gm}{\|\vec{r}(t)\|^{3}}\vec{x}(t)\qquad \qquad (4')
\vec{x}(t_0) = 0,\quad \frac{d \vec{x}}{dt}(t_0) = 0\qquad \qquad \qquad (5')
常微分方程式(4')を初期条件(5')のもとで解くと、
\vec{x}\triangleq \vec h \times \vec{r}\equiv 0
が得られる。(注2参照)
これよりベクトル\vec{r}(t)は,常にベクトル\vec{h}と直交し、
平面H上にあることが証明された。
(注1)
(注2)証明をみれば明らかなように、
惑星に働く力が、太陽と惑星を結ぶ直線と同じ方向(求心力あるいが遠心力)でありさえすれば、
惑星は平面上を運動する。
惑星の角運動量は保存される
式(1)の両辺に、左側から\vec{r}(t)をベクトル積として作用させると、
\vec{r}(t)\times \ \frac{d^{2} \vec{r}}{dt^2}(t)=-\vec{r}(t)\times \frac{Gm}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|} \qquad \qquad (a)
ここで、ベクトル積の性質から
\quad \frac{d}{dt}\Bigl(\vec{r}\times \frac{d \vec r}{dt}\Bigr)=\frac{d \vec r}{dt}\times \frac{d \vec{r}}{dt}+\vec{r}\times \ \frac{d^{2} \vec{r}}{dt^2}
\quad =\vec{r}\times \ \frac{d^{2} \vec{r}}{dt^2}
\quadまた右辺は
\quad -\vec{r}(t)\times \frac{Gm}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|} = 0(注参照)
なので、式(a)は、
\frac{d}{dt}\Bigl(\vec{r}\times \frac{d \vec r}{dt}\Bigr)=0
となる。
(注)万有引力でなく、惑星に作用する力の向きが原点向きであるか、それと逆向きという条件だけで良い。
\frac{d \vec r}{dt}=\vec v(速度ベクトル)なので、上式は
\frac{d}{dt}\Bigl(\vec{r}\times \vec v \Bigr)=0\qquad \qquad (b)
あるいは、惑星の運動量\vec p=m\vec vを用いて
\frac{d}{dt}\Bigl(\vec{r}\times \vec p \Bigr)=0\qquad \qquad (c)
となる。この2式から、
\vec{r}\times \vec v と\vec{r}\times \vec p
は、ともに時不変(時間が経過しても一定値)であることが分かる。
定義(角運動量)
\vec{r}(t)\times \vec{p}(t) を、(原点からみた)惑星の時刻tにおける角運動量という。
すると上述したことから、次の法則が証明された。
角運動量の保存法則
求心力、ないし遠心力(力の向きが原点向きか、その逆向き)だけを受けて運動する質点の角運動量は保存される。
ケプラーの第2法則の導出
上述したことから惑星は平面上を動くので、
位置ベクトルと運動量ベクトルが、xy平面に載る様な正規直交右手系xyzを選べる。
この座標系のもとで座標表示すると
\vec{r(t)}\times \vec p(t)=\bigl(r_{1}(t),r_{2}(t),0)\bigr)\times \bigl(p_{1}(t),p_{2}(t),0\bigr)
=\bigl(0,0,r_{1}(t)p_{2}(t)-r_{2}(t)p_{1}(t)\bigr)
他方、ベクトル\vec{r}(t) と ベクトル\vec{p}(t)のなす劣角(180度以下のほうの角)を\theta(t)とすると
この2つのベクトルと直交するz方向の\vec{r}(t)\times \vec{p}(t) の成分は
ベクトル積の定義から、
\|\vec{r}(t)\| \times \| \vec{p}(t) \| \sin \theta(t)
なので、
\|\vec{r}(t)\| \times \| \vec{p}(t) \| \sin \theta(t) \qquad \qquad \qquad (6)
が、時不変であることが分かる。
惑星質量は一定なので
式(6)から
\|\vec{r}(t)\| \times \| \vec{v}(t) \| \sin \theta(t)/2 \qquad \qquad \qquad (7)
が、時不変であることが分かる。
式(7)は動径OP(t)が時刻tの瞬間に単位時間あたりに掃く面積を表すので、惑星の面積速度と呼ばれる。
系 ケプラーの第2法則
求心力、ないし遠心力(力の向きが原点向きか、その逆向き)だけを受けて運動する質点の面積速度は一定である。従って、
「惑星の動径(惑星と太陽とを結ぶ線分)が単位時間に描く面積(面積速度)は、一定である」
ケプラーの第1法則の導出
惑星の運動方程式(1)を解いて惑星の軌道を求めよう。
求心力の下での運動を扱うには、極座標が都合がよい。
惑星運動の微分方程式の極座標表示
そこで惑星の位置ベクトルを、x軸の正の部分のつくる半直線を角度の基準とする極座標で
\vec{r}(t)=r(t)\bigl(\cos \theta (t),\sin \theta (t)\bigr)
と表す。
この式を時間tで微分すると
\dot{\vec{r}}(t)={\dot r}(t)\bigl(\cos \theta (t),\sin \theta (t)\bigr)
+r(t){\dot \theta} (t)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)
\qquad 但し、\dot{\vec{r}}(t)\triangleq \frac{d\vec{r}}{dt}(t)
さらにtで微分すると
\ddot{\vec{r}}(t)
=\ddot {r}(t)\bigl(\cos \theta (t),\sin \theta (t)\bigr)
+\dot{r}(t)\dot{\theta} (t)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)
+\dot{r}(t)\dot{\theta} (t)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)
+r(t)\ddot{\theta} (t)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)
-r(t)\dot{\theta}^{2}(t)\bigl(\cos \theta (t),\sin \theta (t) \bigr)
=\bigl(\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)\bigr)\bigl(\cos \theta (t),\sin \theta (t)\bigr)
+\bigl(2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \bigr)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)
故に
\ddot{\vec{r}}(t)
=\bigl(\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)\bigr)\bigl(\cos \theta (t),\sin \theta (t)\bigr)
+\bigl(2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \bigr)\bigl(-\sin \theta (t),\cos \theta (t)\bigr) \qquad \qquad \qquad (8)
が得られた。
惑星に働く万有引力を、この極座標で表示すると
-\frac{GM}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|}
=-\frac{GM}{r^{2}(t)}(\cos \theta (t),\sin \theta (t)\bigr) \qquad (9)
式(8),(9)を惑星の運動方程式
\ddot{\vec r}(t)=-\frac{GM}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|}
に代入すると
\bigl(\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)\bigr)\bigl(\cos \theta (t),\sin \theta (t)\bigr)
+\bigl(2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \bigr)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)
= -\frac{GM}{r^{2}(t)}(\cos \theta (t),\sin \theta (t)\bigr)
右辺を左辺に移行すると
\bigl(\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)+\frac{GM}{r^{2}(t)}\bigr)\bigl(\cos \theta (t),\sin \theta (t)\bigr)
+\bigl(2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \bigr)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)\equiv 0
\bigl(\cos \theta (t),\sin \theta (t)\bigr)と\bigl(-\sin \theta (t),\cos \theta (t)\bigr)は直交ベクトルなので、この式から
\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)+\frac{GM}{r^{2}(t)}\equiv 0 \qquad \qquad (10)
2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \equiv 0 \qquad \qquad (11)
が得られる。
命題;惑星運動の微分方程式の極座標表示
太陽の位置Oを原点、x軸の正の部分の半直線を角の基準線とする極座標で
惑星の運動方程式を表現すると
\ddot {r}-r\dot{\theta}^{2}=-\frac{GM}{r^{2}}\qquad \qquad (12)
r^2\dot{\theta} = h \quad (任意の定数) \qquad \qquad (13)
となる。(注2 を参照のこと)
証明
式(12)は式(10)から明らか。
式(13)の左辺をtで微分すると
\frac{d r^2\dot{\theta}}{dt}=r^2\ddot{\theta} + 2r\dot{r}\dot{\theta}
=r\bigl(r\ddot{\theta} + 2\dot{r}\dot{\theta}\bigr)
\qquad 式(11)から
= 0
左辺の時間微分が零なので、左辺は任意の定数hに等しい。
証明終わり。 \qquad \qquad \qquad \qquad \Box
(注)式(13)は面積速度が一定であることを表現している。
先に、このことは証明してあるが、その別証明である。
(注2)関数の変数tまで書くと、式(12),(13)は、次のように書ける。
\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)=-\frac{GM}{r^{2}(t)}\qquad \qquad (12')
r^2(t)\dot{\theta}(t) = h \quad (任意の定数) \qquad \qquad (13')
微分方程式を解いて、惑星の軌道形を求める
連立微分方程式(12)、(13)を解いて、r=r(t),\theta=\theta(t)を求めれば、
任意の時刻tにおける惑星の位置がわかる。
軌道の形は、この2式からtを消去して、rとθの関数関係
f(r,\theta)=0, またはr=f(\theta)
を求めれな良い。
しかし、残念なことにこれは大変困難である。
軌道の形だけを知りたいときには、
これを回避して簡単な(しかし技巧的な)直接的な解法がある。
それは、連立微分方程式の解が
r=f\circ \theta ,\qquad変数まで書けば、r(t)=f\circ \theta (t)\triangleq f(\theta (t))
という合成関数になることを仮定して、
r=f(\theta)を直接求める問題に変換し、解を求める方法である。
軌道の形を直接求めるための微分方程式の導出
\frac{d r(t) }{dt}=\frac{d f\circ \theta (t) }{dt}
\qquad 合成関数の微分公式を用いて、
=\frac{d f}{d\theta}(\theta (t))\dot{\theta}(t)
\qquad 式(13)を\dot{\theta}=\frac{h}{r^2}に変形し、上式に代入すると
=\frac{df}{d\theta}(\theta (t))\frac{h}{r^2(t)}
故に
\frac{d r(t) }{dt}=\frac{df}{d\theta}(\theta (t))\frac{h}{r^2(t)} \qquad \qquad (b)
補題 g(\theta)\triangleq \frac{1}{f(\theta)}とおくと、
\frac{df(\theta)}{d\theta}= -\frac{1}{g^2(\theta)}\frac{dg(\theta)}{d\theta}
証明は、
ウィキブックス高等学校数学III/微分法#商の導関数をご覧ください。
この補題を用いて式(b)を変形すると
\frac{d r(t) }{dt}=-\frac{1}{g^2(\theta (t))}\frac{dg}{d\theta}(\theta (t))\frac{h}{r^2(t)}
\qquad r(t)=f(\theta (t)) なので、g^2(\theta (t))r^2(t)=1となり、
=-h\frac{dg}{d\theta}(\theta (t))
故に
\dot{r}(t)=-h\frac{dg}{d\theta}(\theta (t)) \qquad \qquad (c)
これをさらに時間tで微分すると、合成関数の微分の公式から
\ddot{r}(t)=-h\frac{d^{2}g}{d\theta ^2}(\theta(t))\dot{\theta}(t)
\qquad 式(13)から得られる\dot{\theta}=\frac{h}{r^2}を代入して、
=-h\frac{d^{2}g}{d\theta ^2}(\theta(t)) \frac{h}{r^2(t)}
故に
\ddot{r}(t)=-\frac{h^{2}}{r^2(t)}\frac{d^{2}g}{d\theta ^2}(\theta(t))\qquad \qquad (d)
が得られる。
式(d)と\dot{\theta}=\frac{h}{r^2}を、式(12)に代入して整頓すると次式が得られる。
\frac{d^{2}g}{d\theta ^2}(\theta(t))+\frac{1}{r(t)}=\frac{GM}{h^2}
\frac{1}{r(t)}=\frac{1}{f(\theta (t))}=g(\theta (t))なので
上式から次式が得られる。
\frac{d^{2}g}{d\theta ^2}(\theta(t))+g(\theta (t))=\frac{GM}{h^2}
従って変数を \theta とする微分方程式
\frac{d^{2}g}{d\theta ^2}(\theta)+g(\theta )=\frac{GM}{h^2} \qquad \qquad (14)
の解g=g(\theta)が分かれば
r=f(\theta)=\frac{1}{g(\theta)}
が、惑星の軌道の形を表す式になる。