ヒルの方程式のソースを表示

{{Otheruses|宇宙工学における方程式|数学における微分方程式の形|ヒル微分方程式|生化学における方程式|ヒルの式}}

'''ヒルの方程式'''（Hill's equations あるいは Clohessy-Wiltshire equations）とは、[[人工衛星の軌道|衛星軌道]]を動く[[質点]]のそばを動く別の質点との相対運動を表す式であり、[[宇宙ステーション]]のような[[人工衛星]]との[[ランデブー (宇宙開発)|ランデブー]]のための[[基礎方程式]]となる<ref name="Clohessy-Wiltshire">{{Cite journal|last1=Clohessy |first1=W. H. |last2=Wiltshire |first2=R. S. |date=1960 |title=Terminal Guidance System for Satellite Rendezvous |url=https://doi.org/10.2514/8.8704 |journal=Journal of the Aerospace Sciences |volume=27 |issue=9 |pages=653–658 |doi=10.2514/8.8704}}</ref>。

== 歴史 ==
1878年に[[ジョージ・ウィリアム・ヒル]]によって地球に対する月の相対的な軌道に関する初期の研究が行われ、その過程でこの方程式は導出された<ref>{{cite journal | vauthors=((Hill, G. W.)) | journal=American Journal of Mathematics | title=Researches in the Lunar Theory | volume=1 | issue=1 | pages=5–26 | publisher=Johns Hopkins University Press | date= 1878 | issn=0002-9327 | doi=10.2307/2369430| jstor=2369430 }} </ref>。 この周期的な要素を含む方程式に関する研究により、後に数学で[[ヒル微分方程式]]として一般化された。

1960年に W. H. Clohessy と R. S. Wiltshire によって軌道上でのランデブーを実行するための一般的な運動方程式として再発見された<ref name="Clohessy-Wiltshire"/>。

== 概説 ==
単純化のため、[[円軌道]]に限定する。軌道上の質点（ターゲット）に、ターゲットとともに動く座標原点をおき、軌道中心から軌道の直径の外向きに x軸をとる。（すなわち、地球を回る人工衛星ならば、 つねに地表側が -x方向となるような[[座標系]]とする。）また軌道運動の進行方向に y軸をとり、 xy面に垂直に z軸をとり、軌道運動の[[角速度]]を<math>\omega</math>とする。ここで[[試験粒子|テスト粒子]]（チェイサー）と中心天体間の[[重力]]のみを考えた[[ニュートンの運動方程式]]を立てて、ターゲットとともに軌道を動く[[回転座標系]]からみた運動方程式を線形化する。このようにして求めたチェイサーの動きをあらわす運動方程式であるヒルの方程式は、以下のようになる。
{{Indent|<math> \begin{align}
F_x &= \ddot x - 2\omega\dot y  - 3 \omega^2 x \\
F_y &= \ddot y + 2\omega\dot x \\
F_z &= \ddot z  + \omega^2 z 
\end{align}</math>}}
ただし<math> F </math>は単位質量あたりの外力、すなわち外力による加速度とする。

式の形から以下のことがわかる。

#<math>\dot{y} > 0</math> と <math>F_x > 0</math> が同じ効果がある。つまり、進行方向へ動くと軌道上方への力が加わったように感じられる（[[コリオリの力]]）。
#<math>\dot{x} > 0</math> と <math>F_y < 0 </math>が同じ効果がある。つまり、軌道の上方に動くと進行方向逆向きの力が加わったように感じられる。
#z軸の動きは他と独立で、[[単振動]]運動となる。

外力 F が無いとき、チェイサーの xy面での動きは、 y軸方向に次第に移動していく楕円の運動となり、'''C-W解'''（Clohessy-Wiltshire解）と呼ばれる。
:<math>
\begin{bmatrix} x \\ y \\ z \end{bmatrix} =
\begin{bmatrix}
 4-3\cos \omega t & 0 & 0 \\
 6(\sin \omega t - \omega t) & 1 & 0 \\
 0 & 0 & \cos \omega t
\end{bmatrix}
\begin{bmatrix} x_0 \\ y_0 \\ z_0 \end{bmatrix} +
\frac{1}{\omega}
\begin{bmatrix}
 \sin \omega t & 2(1- \cos \omega t) & 0 \\
 -2(1- \cos \omega t) & 4 \sin \omega t -3 \omega t & 0 \\
 0 & 0 & \sin \omega t
\end{bmatrix}
\begin{bmatrix} \dot x_0 \\ \dot y_0 \\ \dot z_0 \end{bmatrix} 
</math>.
xy面で移動する楕円の進行方向の動きを使って、チェイサーがターゲットの進行方向前方から接近することを +V bar 接近（後方からの場合は −V bar接近）、下方から接近することを +R bar 接近（上方からの場合は −R bar 接近）という（ R は radius から）。

==参考文献==
<references />
* 茂原正道・木田隆、『宇宙工学入門II —宇宙ステーションと惑星間飛行のための誘導・制御』、培風館、1998、ISBN 4-563-03530-0
* 木田隆・小松敬治・川口淳一郎、『人工衛星と宇宙探査機』（宇宙工学シリーズ３）、コロナ社、2001、ISBN 4-339-01223-8

{{DEFAULTSORT:ひるのほうていしき}}
[[Category:天体力学]]
[[Category:軌道]]
[[Category:人工衛星の軌道]]
[[Category:微分方程式]]
[[Category:物理学の方程式]]