楕円

曖昧さ回避 長円」はこの項目へ転送されています。仏師については「長円 (仏師)」をご覧ください。
円錐切断面の4つのタイプ(放物線(左)、楕円(中央)、円(中央)、双曲線(右))

楕円(だえん、正字: 橢圓: ellipse)とは、平面上のある2定点からの距離の和が一定となるような点の集合から作られる曲線である。

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

概要

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

楕円の内部に2焦点を通る直線を引くとき、これを長軸という。長軸の長さを長径という。長軸と楕円との交点では2焦点からの距離の差が最大となる。また、長軸の垂直二等分線を楕円の内部に引くとき、この線分を短軸という。短軸の長さを短径という。

用語

  • 長軸と短軸の交点は楕円の中心と呼ばれる。
  • 長軸を中心で分けた2つの線分は半長軸と呼ばれ、その長さを長半径という。
  • 短軸を中心で分けた2つの線分は半短軸と呼ばれ、その長さを短半径という。
  • 短径と長径の比は楕円率と呼ばれる。
楕円の長軸(緑線)と短軸(ピンク線)

楕円の方程式

2次元直交座標系で、原点 O が長軸と短軸の交点となる楕円は代数的に次のように書ける。これを標準形という。

x 2 a 2 + y 2 b 2 = 1. {\displaystyle {\frac {x^{2}}{a^{2}}}+{\frac {y^{2}}{b^{2}}}=1.}

a > b > 0 のとき、2a は長軸の長さ(長径)、2b は短軸の長さ(短径)となる。xy 平面上にグラフを書くと横長の楕円となる。また、焦点はx 軸上にあり、その座標は ( a 2 b 2 , 0 ) , ( a 2 b 2 , 0 ) {\displaystyle \left({\sqrt {a^{2}-b^{2}}},0\right),\left(-{\sqrt {a^{2}-b^{2}}},0\right)} となる。

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

頂点の座標は ab のとき ( ± a , 0 ) , ( 0 , ± b ) {\displaystyle (\pm a,0),(0,\pm b)} となる。

同じ楕円は、t媒介変数とする媒介変数表示では、次のように表現できる。

x = a cos t {\displaystyle x=a\,\cos t}
y = b sin t {\displaystyle y=b\,\sin t}
0 t < 2 π {\displaystyle 0\leq t<2\pi }

ただし、t は (x,y) ベクトルのx軸に対する角度ではない(天体力学では離心近点角と呼ばれる)。

「緯度#更成緯度 (reduced latitude)」も参照
媒介変数表示により表された楕円上の点Pと媒介変数tの関係。tは点Pとx軸の角度とは異なる。

また、 u = tan ( t / 2 ) {\displaystyle u=\tan(t/2)} と置くと、

cos ( t ) = 1 u 2 1 + u 2 sin ( t ) = 2 u 1 + u 2 {\displaystyle {\begin{aligned}\cos(t)&={\frac {1-u^{2}}{1+u^{2}}}\\\sin(t)&={\frac {2u}{1+u^{2}}}\end{aligned}}}

となるので、下記の表現でも楕円を表すことができる。この場合uの範囲は[0,1]である。

x = a ( 1 u 2 ) 1 + u 2 y = 2 b u 1 + u 2 {\displaystyle {\begin{aligned}x&={\frac {a(1-u^{2})}{1+u^{2}}}\\y&={\frac {2bu}{1+u^{2}}}\end{aligned}}}

複素平面Cにおいては,Cの二点 a 1 , a 2 {\displaystyle a_{1},a_{2}} からの点 z {\displaystyle z} への 距離 r 1 , r 2 {\displaystyle r_{1},r_{2}} の和が l {\displaystyle l} であるものの軌跡である。

r 1 , 2 = | z a 1 , 2 | r 1 + r 2 = l {\displaystyle {\begin{aligned}&r_{1,2}=|z-a_{1,2}|\\&r_{1}+r_{2}=l\end{aligned}}}

楕円の幾何学的諸量

楕円の形状は離心率 e で表現される。

e = 1 b 2 a 2 {\displaystyle e={\sqrt {1-{\frac {b^{2}}{a^{2}}}}}}

別途、扁平率 f でも表現できる。

f = 1 b a {\displaystyle f=1-{\frac {b}{a}}}

楕円の面積 S は次のように表現できる。

S = π a b {\displaystyle S=\pi ab\,}

楕円の周長 Ca > b のとき、第二種完全楕円積分を用いて次のように表現できる。

C = 4 0 π 2 a 2 cos 2 t + b 2 sin 2 t d t = 4 a 0 π 2 1 e 2 sin 2 t d t = 2 π a n = 0 e 2 n 1 2 n m = 1 n ( 1 1 2 m ) 2 {\displaystyle {\begin{aligned}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{aligned}}}

また n = f / ( 2 f ) {\displaystyle n=f/(2-f)} とおき、二項係数を使って、次のようにも表現できる(Gauss-Kummer級数)[1]

C = 2 π a 1 + n i = 0 ( 1 / 2 i ) 2 n 2 i . {\displaystyle {\begin{aligned}C={\frac {2\pi a}{1+n}}\sum _{i=0}^{\infty }{\binom {1/2}{i}}^{2}n^{2i}.\end{aligned}}}

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

C = 8 π Q 5 / 4 n = 0 ( 1 12 ) n ( 5 12 ) n ( v 1 + n v 2 ) r n ( n ! ) 2 {\displaystyle 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 = 432 ( a 2 b 2 ) 2 ( a b ) 6 b a Q 3 {\displaystyle r={\tfrac {432\left(a^{2}-b^{2}\right)^{2}(a-b)^{6}ba}{Q^{3}}}\,}
Q = b 4 + 60 a b 3 + 134 a 2 b 2 + 60 a 3 b + a 4 {\displaystyle Q=b^{4}+60ab^{3}+134a^{2}b^{2}+60a^{3}b+a^{4}\,}
v 1 = b a ( 15 b 4 + 68 a b 3 + 90 a 2 b 2 + 68 a 3 b + 15 a 4 ) {\displaystyle v_{1}=ba\left(15b^{4}+68ab^{3}+90a^{2}b^{2}+68a^{3}b+15a^{4}\right)\,}
v 2 = a 6 b 6 + 126 a b 5 + 1041 a 2 b 4 + 1764 a 3 b 3 + 1041 a 4 b 2 + 126 a 5 b {\displaystyle v_{2}=-a^{6}-b^{6}+126ab^{5}+1041a^{2}b^{4}+1764a^{3}b^{3}+1041a^{4}b^{2}+126a^{5}b\,}

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

C π [ 3 ( a + b ) ( 3 a + b ) ( a + 3 b )   ] = π [ 3 ( a + b ) 10 a b + 3 ( a 2 + b 2 )   ] {\displaystyle 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 π ( a + b ) [ 1 + 3 ( a b a + b ) 2 10 + 4 3 ( a b a + b ) 2 ] {\displaystyle 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 {\displaystyle x=a\,\cos t,\,y=b\,\sin t}

で表した時、 t = t 1 {\displaystyle t=t_{1}} から t = t 2 {\displaystyle t=t_{2}} までの弧長 L {\displaystyle L}

L = t 1 t 2 ( d x d t ) 2 + ( d y d t ) 2 d t = t 1 t 2 a 2 sin 2 t + b 2 cos 2 t d t {\displaystyle {\begin{aligned}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{aligned}}}

で求められる。これは、 a , b {\displaystyle a,b} の大小関係に関係なく成立する。

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

その前に、媒介変数表示について、補足しておく。 楕円の媒介変数表示には、通常

x = a cos t , y = b sin t {\displaystyle x=a\,\cos t,\,y=b\,\sin t}

が用いられる。この場合、t = 0 では、点 ( a , 0 ) {\displaystyle (a,0)} をとり、 t = π / 2 {\displaystyle \pi /2} では点 ( 0 , b ) {\displaystyle (0,b)} をとるので、tはx軸の正の部分を基準線とする反時計方向の角度になっている。

一方、媒介変数表示は

x = a sin t , y = b cos t {\displaystyle x=a\,\sin t,\,y=b\,\cos t}

とする事もでき、この場合、t = 0 では、点 ( 0 , b ) {\displaystyle (0,b)} をとり、 t = π / 2 {\displaystyle \pi /2} では点 ( a , 0 ) {\displaystyle (a,0)} をとるので、tはy軸の正の部分を基準線とする時計方向の角度になっている。

第二種不完全楕円積分を

E ( φ , k ) = 0 φ 1 k 2 sin 2 ϕ   d ϕ {\displaystyle {\begin{aligned}E(\varphi ,k)&=\int _{0}^{\varphi }{\sqrt {1-k^{2}\sin ^{2}\phi }}~d\phi \end{aligned}}}
0 φ π / 2 {\displaystyle 0\leq \varphi \leq \pi /2}

と表記する。 さらに、楕円上の点を指定する指標として、 ( x , y ) {\displaystyle (x,y)} ベクトルのx軸に対する角度 θ {\displaystyle \theta } も導入する。

( tan θ = y / x , π / 2 θ π / 2 {\displaystyle \tan \theta =y/x,-\pi /2\leq \theta \leq \pi /2} )

A) 0 < b a {\displaystyle 0<b\leq a} の時

e = 1 ( b / a ) 2 {\displaystyle e={\sqrt {1-(b/a)^{2}}}}

楕円(の右半分)を

x = a sin u , y = b cos u , 0 u π {\displaystyle x=a\,\sin u,\,y=b\,\cos u,\,0\leq u\leq \pi }

で表す。 a E ( u , e ) {\displaystyle a\,E(u,e)} は点 ( 0 , b ) {\displaystyle (0,b)} から u {\displaystyle u} が与える点までの弧長となっている。

この時

y / x = tan θ = ( b / a ) / tan u {\displaystyle y/x=\tan \theta =(b/a)/\tan u}
tan u = ( b / a ) / tan θ {\displaystyle \tan u=(b/a)/\tan \theta }
L = u 1 u 2 ( d x d u ) 2 + ( d y d u ) 2 d u = u 1 u 2 a 2 cos 2 u + b 2 sin 2 u d u = a u 1 u 2 1 e 2 sin 2 u d u {\displaystyle {\begin{aligned}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{aligned}}}

となる。 E ( u , e ) {\displaystyle E(u,e)} が点 ( a , 0 ) {\displaystyle (a,0)} を最大の終点とする積分になる事を考慮し、 場合分けをし積分範囲を決めると、次のようになる。

i) 0 θ 1 θ 2 π / 2 {\displaystyle 0\leq \theta _{1}\leq \theta _{2}\leq \pi /2}
L = a ( E ( u 2 , e ) E ( u 1 , e ) ) {\displaystyle L=a\,(E(u_{2},e)-E(u_{1},e))}
u 1 = t h 2 u ( θ 2 ) , u 2 = t h 2 u ( θ 1 ) {\displaystyle u_{1}=th2u(\theta _{2}),\,u_{2}=th2u(\theta _{1})}
ii) π / 2 θ 1 < 0 θ 2 π / 2 {\displaystyle -\pi /2\leq \theta _{1}<0\leq \theta _{2}\leq \pi /2}
L = a ( 2 E ( π / 2 , e ) E ( u 2 , e ) E ( u 1 , e ) ) {\displaystyle L=a\,(2E(\pi /2,e)-E(u_{2},e)-E(u_{1},e))}
u 1 = t h 2 u ( θ 2 ) , u 2 = t h 2 u ( θ 1 ) {\displaystyle u_{1}=th2u(\theta _{2}),\,u_{2}=th2u(-\theta _{1})}
iii) π / 2 θ 1 θ 2 < 0 {\displaystyle -\pi /2\leq \theta _{1}\leq \theta _{2}<0}
L = a ( E ( u 2 , e ) E ( u 1 , e ) ) {\displaystyle L=a\,(E(u_{2},e)-E(u_{1},e))}
u 1 = t h 2 u ( θ 1 ) , u 2 = t h 2 u ( θ 2 ) {\displaystyle u_{1}=th2u(-\theta _{1}),\,u_{2}=th2u(-\theta _{2})}

ここで

t h 2 u ( θ ) = tan 1 ( ( b / a ) / tan θ ) {\displaystyle th2u(\theta )=\tan ^{-1}((b/a)/\tan \theta )}
(ただし、 t h 2 u ( 0 ) = π / 2 , t h 2 u ( ± π / 2 ) = 0 {\displaystyle th2u(0)=\pi /2,\,th2u(\pm \pi /2)=0} とする)

である。

B) 0 < a b {\displaystyle 0<a\leq b} の時

e = 1 ( a / b ) 2 {\displaystyle e={\sqrt {1-(a/b)^{2}}}}

楕円(の右半分)を

x = a cos v , y = b sin v , π / 2 v π / 2 {\displaystyle x=a\,\cos v,\,y=b\,\sin v,\,-\pi /2\leq v\leq \pi /2}

で表す。 b E ( v , e ) {\displaystyle b\,E(v,e)} は点 ( a , 0 ) {\displaystyle (a,0)} から v {\displaystyle v} が与える点までの弧長となっている。

この時

y / x = tan θ = ( b / a ) tan v {\displaystyle y/x=\tan \theta =(b/a)\tan v}
tan v = ( a / b ) tan θ {\displaystyle \tan v=(a/b)\tan \theta }
L = v 1 v 2 ( d x d v ) 2 + ( d y d v ) 2 d v = v 1 v 2 a 2 sin 2 v + b 2 cos 2 v d v = b v 1 v 2 1 e 2 sin 2 v d v {\displaystyle {\begin{aligned}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{aligned}}}

となる。 E ( v , e ) {\displaystyle E(v,e)} が点 ( a , 0 ) {\displaystyle (a,0)} を始点とする積分になる事を考慮し、 場合分けをし積分範囲を決めると、次のようになる。

i) 0 θ 1 θ 2 π / 2 {\displaystyle 0\leq \theta _{1}\leq \theta _{2}\leq \pi /2}
L = b ( E ( v 2 , e ) E ( v 1 , e ) ) {\displaystyle L=b\,(E(v_{2},e)-E(v_{1},e))}
v 1 = t h 2 v ( θ 1 ) , v 2 = t h 2 v ( θ 2 ) {\displaystyle v_{1}=th2v(\theta _{1}),\,v_{2}=th2v(\theta _{2})}
ii) π / 2 < θ 1 < 0 θ 2 π / 2 {\displaystyle -\pi /2<\theta _{1}<0\leq \theta _{2}\leq \pi /2}
L = b ( E ( v 2 , e ) + E ( v 1 , e ) ) {\displaystyle L=b\,(E(v_{2},e)+E(v_{1},e))}
v 1 = t h 2 v ( θ 1 ) , v 2 = t h 2 v ( θ 2 ) {\displaystyle v_{1}=th2v(-\theta _{1}),\,v_{2}=th2v(\theta _{2})}
iii) π / 2 θ 1 θ 2 < 0 {\displaystyle -\pi /2\leq \theta _{1}\leq \theta _{2}<0}
L = b ( E ( v 2 , e ) E ( v 1 , e ) ) {\displaystyle L=b\,(E(v_{2},e)-E(v_{1},e))}
v 1 = t h 2 v ( θ 2 ) , v 2 = t h 2 v ( θ 1 ) {\displaystyle v_{1}=th2v(-\theta _{2}),\,v_{2}=th2v(-\theta _{1})}

ここで

t h 2 v ( θ ) = tan 1 ( ( a / b ) tan θ ) {\displaystyle th2v(\theta )=\tan ^{-1}((a/b)\tan \theta )}
(ただし、 t h 2 v ( ± π / 2 ) = ± π / 2 {\displaystyle th2v(\pm \pi /2)=\pm \pi /2} とする)

である。

作図法

糸を使った作図例
アルキメデスの楕円コンパスen:Trammel of Archimedesを使った作図例
楕円は内トロコイドの特殊な場合として表される。図は rc = 10, rm = 5, rd = 1 の場合。

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

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

歴史

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

脚注

  1. ^ Weisstein, Eric W. "Gauss-Kummer Series". mathworld.wolfram.com (英語).
  2. ^ Cetin Hakimoglu-Brown iamned.com math page

参考文献

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

関連項目

外部リンク

ウィキメディア・コモンズには、楕円に関連するメディアおよびカテゴリがあります。
  • 楕円歯車: オランダ数学者en:Frans van Schootenが著書 "De organica conicarum sectionum in plano descriptione tractatus" において紹介している機構
典拠管理データベース: 国立図書館 ウィキデータを編集
  • フランス
  • BnF data
  • ドイツ
  • イスラエル
  • アメリカ