\documentclass{jsarticle}
\title{ITPASS　数値計算実習課題 その1}
\author{理学部・地球惑星科学科　  高木 瑠美　　joho07}
\begin{document}
\maketitle{}
\def\diff{\mathrm d}
\def\ddt#1{\frac{\diff #1}{\diff t}}
\section*{解答}
	万有引力の法則 
		\[
		  F = -\frac{GMm}{r^{2}}
		\]
	質量が$m_{1}$である中心星（位置ベクトル$\bf r_{1}$）
	、質量が$m_{2}$である惑星(位置ベクトル$\bf r_{2}$)のみで構成される惑星系において、
	惑星の軌道を計算する。
	考える系において、支配的な力は万有引力のみであるとする。
	$G$は万有引力定数。
\section*{1.}
	万有引力の法則より、\\
		中心星に働く力
			\[
			 {\bf F_{12}} = -\frac{Gm_{1}m_{2}}{|{\bf r_{1}-r_{2}}|^{3}}({\bf r_{1}-r_{2}})
			\]
		惑星に働く力
			\[
			{\bf F_{21}} = -\frac{Gm_{1}m_{2}}{|{\bf r_{2}-r_{1}}|^{3}}({\bf r_{2}-r_{1}})
			\]
		となっている。\\
	ここで、相対ベクトル$\bf r = r_{2} - r_{1}$を定義する。\\
	この系において、運動方程式は以下のようになる。\\
		中心星の運動方程式
			\begin{equation}
			m_{1}\frac{d^{2}{\bf r_{1}}}{d{t^{2}}} = {\bf F_{12}} 
			= -\frac{Gm_{1}m_{2}}{|{\bf r}|^{3}}{\bf r} \label{1}
			\end{equation}
		惑星の運動方程式
			\begin{equation}
			m_{2}\frac{d^{2}{\bf r_{2}}}{d{t^{2}}} = {\bf F_{21}} 
			= \frac{Gm_{1}m_{2}}{|{\bf r}|^{3}}{\bf r} \label{2}
			\end{equation}
	$(\ref{2})式 \times \frac{1}{m_{2}} - (\ref{1})式\times \frac{1}{m_{1}}$より、\\
		\[
		\frac{d^{2}{\bf r_{2}}}{d{t^{2}}} - \frac{d^{2}{\bf r_{1}}}{d{t^{2}}} =
		 -G\frac{m_{1}}{|{\bf r}|^{3}}{\bf r}-G\frac{m_{2}}{|{\bf r}|^{3}}{\bf r}
		\]
		\begin{equation}
		\raisebox{.2ex}{.}\raisebox{1.2ex}{.}\raisebox{.2ex}{.}
		\hspace{2em}
		\frac{d^{2}{\bf r}}{d{t^{2}}} = -G\frac{m_{1} + m_{2}}{|{\bf r}|^{3}}{\bf r} \label{3}
		\end{equation}
	換算質量$μ = \frac{m_{1}m_{2}}{m_{1} + m_{2}}$を定義し、惑星に働く力を$\bf F(r)$と置き直す。
	（中心星に働く力は$-\bf F(r)$となる。）\\
	すると、($\ref{3}$)式は、以下のように書ける。
		\begin{equation}
		μ\frac{d^{2}\bf r}{dt^{2}} = \bf F(r) \label{4}
		\end{equation}
	($\ref{4}$)式は、質量$μ$の1つの質点がその位置によって定まる力を受けたときの運動方程式に等しい。\\
	したがって、2体問題が、換算質量$μ$と相対位置$\bf r$を用いれば、1体問題に帰着されることとなる。\\
	また、$\bf F(r)$を使って、惑星の運動方程式と中心星の運動方程式を書き直すと以下のようになる。
		\begin{equation}
		m_{2}\frac{d^{2}{\bf r}}{d{t^{2}}} = \bf F(r) \label{12}
		\end{equation}
		\begin{equation}
		m_{1}\frac{d^{2}{\bf r}}{d{t^{2}}} = -\bf F(r) \label{11}
		\end{equation}
	上式($\ref{12}$)式と(\ref{11})式を足し合わせると、
		\begin{equation}
		m_{1}\frac{d^{2}{\bf r}}{d{t^{2}}} + m_{2}\frac{d^{2}{\bf r}}{d{t^{2}}} = \bf 0
		\end{equation}
	これより、$\frac{d^{2}{\bf r_{G}}}{d{t^{2}}} = 
		\frac{m_{1}\frac{d^{2}{\bf r}}{d{t^{2}}} + m_{2}\frac{d^{2}{\bf r}}{d{t^{2}}}}{m_{1} + m_{2}} = \bf 0$
	となるから、質量中心（この系での重心）は、等速直線運動することが分かる。
\section*{2.}
	1.の運動方程式を成分に分ける。$(\ref{3})$式について\\
	相対ベクトル${\bf r} = (x , y)$に対して、速度は以下のように定義する。
		\[
		{\bf v} \equiv (v_{x} , v_{y}) = \left(\ddt{x} , \ddt{y}\right)
		\]
	これより、相対ベクトルに対する加速度は、
		\[
		{\bf a} \equiv \left(\ddt{v_{x}} , \ddt{v_{y}}\right)
		\]
	したがって、$(\ref{3})$式は、
		\[
		\left(\begin{array}{cc}
		\ddt{v_{x}} \\
		\ddt{v_{y}}
		\end{array} 
		\right)
		= -G\frac{(m_{1} + m_{2})}{(x^{2} + y^{2})^{\frac{3}{2}}}
		\left(\begin{array}{cc} 
		x \\
		y 
		\end{array} 
		\right)
		\]
		\[
		\left \{
		\begin{array}{cc}
		\ddt{v_{x}} = -G\frac{(m_{1} + m_{2})}{(x^{2} + y^{2})^{\frac{3}{2}}}x \\
		\ddt{v_{y}} = -G\frac{(m_{1} + m_{2})}{(x^{2} + y^{2})^{\frac{3}{2}}}y
		\end{array}
		\right.
		\]
	となる。\\
	\\
	\\
以上、ITPASS 数値計算実習課題 その1 の解答である。
\end{document}