楕円

**楕円**（だえん、、）とは、平面上のある2定点からの距離の和が一定となるような点の集合から作られる曲線である。

基準となる2定点を焦点という。円錐曲線の一種である。

## 概要
2つの焦点が近いほど楕円は円に近づき、2つの焦点が一致したとき楕円はその点を中心とした円になる。そのため円は楕円の特殊な場合であると考えることもできる。

楕円の2焦点を通る直線と楕円の2交点を端点とした線分を長軸という。長軸の長さを長径という。長軸と楕円との交点では2焦点からの距離の差が最大となる。また、長軸の垂直二等分線と楕円の2交点を端点とした線分を短軸という。短軸の長さを短径という。

### 用語
* 長軸と短軸の交点は楕円のと呼ばれる。
* 長軸を中心で分けた2つの線分はと呼ばれ、その長さをという。
* 短軸を中心で分けた2つの線分はと呼ばれ、その長さをという。
* 短径と長径の比はと呼ばれる。

## 楕円の方程式

### 一般形
2次元直交座標系において、楕円の2焦点の座標をそれぞれ (a, b), (c, d)、焦点からの距離の和を k とする。このとき楕円の方程式は、次のように表される。これを**一般形**という。

\sqrt{(x-a)^2 + (y-b)^2} + \sqrt{(x-c)^2 + (y-d)^2} = k 

この方程式は、うまく式変形することにより、必ず

Ax^2 + Bxy + Cy^2 + Dx + Ey + F = 0　　（ただし A, B, C, D, E, F は実数）

という形に表すことができる。証明は以下の通り。

### 標準形
原点 *O* が長軸と短軸の交点となる楕円は、代数的に次のように書ける。これを**標準形**という。
:\frac{x^{2}}{a^{2}} + \frac{y^{2}}{b^{2}} = 1. 
*a* > *b* > 0 のとき、2*a* は長軸の長さ（長径）、2*b* は短軸の長さ（短径）となる。*xy* 平面上にグラフを書くと横長の楕円となる。また、焦点は*x* 軸上にあり、その座標は \left( \sqrt{a^2 - b^2} , 0\right) , \left( -\sqrt{a^2 - b^2} , 0\right) となる。

*b* > *a* > 0 のときは逆に、2*b* が長軸の長さ（長径）、2*a* が短軸の長さ（短径）となる。したがって、*xy* 平面上にグラフを書くと縦長の楕円となる。また、焦点は *y* 軸上にあり、その座標は \left(0, \sqrt{b^2 - a^2} \right) , \left(0, -\sqrt{b^2 - a^2} \right)となる。(*a* = *b*の時は円となる)

頂点の座標は *a* ≠ *b* のとき (\pm a, 0), (0, \pm b)  となる。

### 媒介変数表示
同じ楕円は、*t* を媒介変数とする媒介変数表示では、次のように表現できる。
:x = a\,\cos t
:y = b\,\sin t
:0 \leq t 
ただし、*t* は (x,y) ベクトルのx軸に対する角度ではない（天体力学では離心近点角と呼ばれる）。

また、u = \tan(t/2) と置くと、
:
\begin{align}
\cos(t)&=\frac{1-u^2}{1+u^2} \\
\sin(t)&=\frac{2u}{1+u^2}
\end{align}

となるので、下記の表現でも楕円を表すことができる。この場合uの範囲は[0,1]である。
:
\begin{align}
x&=\frac{a(1-u^2)}{1+u^2} \\
y&=\frac{2bu}{1+u^2}
\end{align}


### 複素数平面
複素平面においては,の二点a_1, a_2からの点zへの
距離r_1, r_2の和がlであるものの軌跡である。
:
\begin{align}
 & r_{1,2} = | z - a_{1,2} | \\
 & r_1 + r_2 = l
\end{align}


## 楕円の幾何学的諸量
楕円の形状は離心率 *e* で表現される。
:e = \sqrt{1 - \frac{b^2}{a^2}}
別途、扁平率 *f* でも表現できる。
:f = 1 - \frac{b}{a}

楕円の面積 *S* は次のように表現できる。
:S = \pi ab\, 

楕円の周長 *C* は *a* > *b* のとき、第二種完全楕円積分を用いて次のように表現できる。
:\begin{align}
C & = 4 \int_{0}^{\frac{\pi}{2}} \sqrt{a^2 \cos^2 t + b^2 \sin^2 t}\,dt \\
 & = 4a \int_{0}^{\frac{\pi}{2}} \sqrt{1-e^2 \sin^2 t}\,dt\\
 & = 2\pi a \sum_{n=0}^\infty {e^{2n}\over 1-2n} \prod_{m=1}^n \left(1-{ 1 \over 2m}\right)^2
\end{align}
また n=f/(2-f) とおき、二項係数を使って、次のようにも表現できる（Gauss-Kummer級数）。
:\begin{align}
C =  \frac{2\pi a}{1+n} \sum_{i=0}^\infty \binom{1/2}{i}^2 n^{2i} .
\end{align}

計算機で計算する場合に有用な式としては、分母が \tfrac{27}{1024} \left (\tfrac{a-b}{a+b} \right )^{8}  の率で消える式が次のように導出されている。

:
C = \frac{8\pi}{Q^{5/4}}\sum_{n=0}^\infty  \frac{ \left(\tfrac{1}{12} \right)_{n} \left(\tfrac{5}{12} \right)_{n} \left(v_{1}+nv_{2} \right)r^{n}}{ \left(n! \right)^{2}} 

:: r = \tfrac{432 \left(a^{2}-b^{2} \right)^{2}(a-b)^{6}ba}{Q^3}\, 
:: Q = b^{4}+60ab^{3}+134a^{2}b^{2}+60a^{3}b+a^{4}\, 
:: v_{1} = ba \left(15b^{4}+68ab^{3}+90a^{2}b^{2}+68a^{3}b+15a^{4} \right)\, 
:: v_{2} = -a^{6}-b^{6}+126ab^{5}+1041a^{2}b^{4}+1764a^{3}b^{3}+1041a^{4}b^{2}+126a^{5}b\,

近似式としては、シュリニヴァーサ・ラマヌジャンによる次の二式がある。簡便なものとしては、

:C \approx \pi \left[3(a+b) - \sqrt{(3a+b)(a+3b)}~\right] = \pi \left[3(a+b) - \sqrt{10ab+3 \left(a^2+b^2 \right)} ~ \right]

があり、さらに良い近似として、次式がある。

:C \approx \pi (a + b) \left[ 1 + \frac{3 \left( \dfrac{a-b}{a+b} \right)^2}{10 + \sqrt{4 - 3 \left( \dfrac{a-b}{a+b} \right)^2}} \right]

より一般的には、対応する角度の関数としての、周長の一部である楕円弧長は、第二種不完全楕円積分で表される。

## 楕円弧長と第二種不完全楕円積分の関係の詳細
楕円を媒介変数表示
:x = a\,\cos t,\,y = b\,\sin t
で表した時、t = t_1からt = t_2までの弧長Lは
:\begin{align}
L & =  \int_{t_1}^{t_2} \sqrt{\left({dx \over dt}\right)^2 + \left({dy \over dt}\right)^2}\,dt \\
 & = \int_{t_1}^{t_2} \sqrt{a^2 \sin^2 t + b^2 \cos^2 t}\,dt
\end{align}
で求められる。これは、a,bの大小関係に関係なく成立する。

この式は第二種不完全楕円積分で表す事ができるが、
a,bの大小関係やt_1,t_2の範囲により場合分けが必要になる為、以下に詳述する。

その前に、媒介変数表示について、補足しておく。
楕円の媒介変数表示には、通常
:x = a\,\cos t,\,y = b\,\sin t
が用いられる。この場合、t = 0 では、点(a, 0)をとり、
t = \pi/2では点(0, b)
をとるので、tはx軸の正の部分を基準線とする反時計方向の角度になっている。

一方、媒介変数表示は
:x = a\,\sin t,\,y = b\,\cos t
とする事もでき、この場合、t = 0 では、点(0, b)をとり、
t = \pi/2では点(a, 0)
をとるので、tはy軸の正の部分を基準線とする時計方向の角度になっている。

第二種不完全楕円積分を
:\begin{align}
E(\varphi,k) &= \int_0^\varphi \sqrt{1-k^2\sin^2\phi} ~ d\phi\end{align}

:0 \leq \varphi \leq \pi/2
と表記する。
さらに、楕円上の点を指定する指標として、(x,y)ベクトルのx軸に対する角度\thetaも導入する。
:(\tan\theta = y/x, -\pi/2 \leq \theta \leq \pi/2)

A) 0  の時
:e = \sqrt{1 - (b/a)^2}
楕円(の右半分)を
:x = a\,\sin u,\, y = b\,\cos u,\,0 \leq u \leq \pi
で表す。a\,E(u,e)は点(0,b)からuが与える点までの弧長となっている。

この時
:y/x = \tan\theta = (b/a)/\tan u
:\tan u = (b/a)/\tan\theta
:\begin{align}
L & =  \int_{u_1}^{u_2} \sqrt{\left({dx \over du}\right)^2 + \left({dy \over du}\right)^2}\,du \\
 & = \int_{u_1}^{u_2} \sqrt{a^2 \cos^2 u + b^2 \sin^2 u}\,du \\
 & = a\,\int_{u_1}^{u_2} \sqrt{1-e^2 \sin^2 u}\,du
\end{align}
となる。E(u,e)が点(a,0)を最大の終点とする積分になる事を考慮し、
場合分けをし積分範囲を決めると、次のようになる。
:i) 0 \leq \theta_1 \leq \theta_2 \leq \pi/2
::L = a\,(E(u_2,e)-E(u_1,e))
::u_1 = th2u(\theta_2),\,u_2 = th2u(\theta_1)
:ii) -\pi/2 \leq \theta_1 
::L = a\,(2E(\pi/2,e)-E(u_2,e)-E(u_1,e))
::u_1 = th2u(\theta_2),\,u_2 = th2u(-\theta_1)
:iii) -\pi/2 \leq \theta_1 \leq \theta_2 
::L = a\,(E(u_2,e)-E(u_1,e))
::u_1 = th2u(-\theta_1),\,u_2 = th2u(-\theta_2)
ここで
:th2u(\theta) = \tan^{-1}((b/a)/\tan\theta)
:（ただし、th2u(0) = \pi/2,\,th2u(\pm\pi/2) = 0とする）
である。

B) 0  の時
:e = \sqrt{1 - (a/b)^2}
楕円(の右半分)を
:x = a\,\cos v,\,y = b\,\sin v,\,-\pi/2 \leq v \leq \pi/2
で表す。b\,E(v,e)は点(a,0)からvが与える点までの弧長となっている。

この時
:y/x = \tan\theta = (b/a)\tan v
:\tan v = (a/b)\tan\theta
:\begin{align}
L & =  \int_{v_1}^{v_2} \sqrt{\left({dx \over dv}\right)^2 + \left({dy \over dv}\right)^2}\,dv \\
 & = \int_{v_1}^{v_2} \sqrt{a^2 \sin^2 v + b^2 \cos^2 v}\,dv \\
 & = b\,\int_{v_1}^{v_2} \sqrt{1-e^2 \sin^2 v}\,dv
\end{align}
となる。E(v,e)が点(a,0)を始点とする積分になる事を考慮し、
場合分けをし積分範囲を決めると、次のようになる。
:i) 0 \leq \theta_1 \leq \theta_2 \leq \pi/2
::L = b\,(E(v_2,e)-E(v_1,e))
::v_1 = th2v(\theta_1),\,v_2 = th2v(\theta_2)
:ii) -\pi/2 
::L = b\,(E(v_2,e)+E(v_1,e))
::v_1 = th2v(-\theta_1),\,v_2 = th2v(\theta_2)
:iii) -\pi/2 \leq \theta_1 \leq \theta_2 
::L = b\,(E(v_2,e)-E(v_1,e))
::v_1 = th2v(-\theta_2),\,v_2 = th2v(-\theta_1)
ここで
:th2v(\theta) = \tan^{-1}((a/b)\tan\theta)
:（ただし、th2v(\pm\pi/2) = \pm\pi/2とする）
である。
## 作図法

2つの焦点に、焦点間距離よりも長い1本の糸の両端をそれぞれ固定し、糸が張る状態で節に取り付けた筆記具を動かす。この他、楕円コンパス、楕円テンプレートなどを使って作図はできる。

また、内トロコイドの特殊な場合に楕円が描画される。

## 歴史
中国語で楕円の楕は「木の切り株」の意味で「木の切り口」の
形から名付けられたと考えられている。
日本では田畑の実際の形から「飯櫃」「平卵形」などと呼ばれていたが、関孝和は「側円」と呼んだ。江戸時代には側円と呼ばれ明治になって楕円と呼ばれるようになった。

## 脚注


## 参考文献
*『曲線の事典 性質・歴史・作図法』礒田正美、Maria G. Bartolini Bussi編、田端毅、讃岐勝、礒田正美著:共立出版、2009年 ISBN 9784320019072

## 関連項目
* 円錐曲線: 楕円は円錐曲線のひとつに分類される。円錐曲線には放物線、双曲線も含まれる。
* 楕円軌道: 惑星、衛星、人工衛星の軌道は楕円軌道を描く。
** ケプラーの法則
* 楕円曲線: 楕円形をした曲線のことを指している用語ではないので、注意が必要である。
** 楕円曲線暗号
* オーバル
**オーバルオフィス：アメリカ大統領の執務室（ホワイトハウス）。
* 楕円積分: 元々、楕円やレムニスケートの周長から研究されてきた。
* 楕円関数: 第一種楕円積分の逆関数を発端として研究されてきた。
* 楕円体
** 回転楕円体
** 地球楕円体
* 子午線弧
* ラグビーボール: 長径280 - 300[mm], 長径方向の外周740 - 770[mm], 短径方向の外周580 - 620[mm]と規定されている。(, 短径そのものは規定されていない)

## 外部リンク

*楕円歯車: オランダの数学者:en:Frans van Schootenが著書 "''De organica conicarum sectionum in plano descriptione tractatus''" において紹介している機構