logo

타원의 둘레 📂기하학

타원의 둘레

대부분의 자료에서 제2 종 타원적분이 어떻게 유도되는지 그 과정이 자세하게 나와있지 않다. 있더라도 틀린 경우가 많아1 ‘정확하고’ ‘자세한’ 내용을 직접 작성했다. 참고로 보아스 수리물리학 3판의 내용도 틀렸다.

공식

장반경이 aa, 단반경이 bb, 이심률이 k2k^{2}타원의 둘레는 다음과 같이 계산된다.

E=4a0π21k2cos2θdθ,k2=a2b2a2 E=4a\int _{0} ^{{\textstyle \frac{\pi}{2}}} \sqrt{ 1-k^{2}\cos^{2} \theta } d\theta,\quad k^{2}=\frac{a^{2}-b^{2}}{a^{2}}

유도

우선 타원위의 점 P=(x,y)P=(x,y)를 각도를 통해 표현해보자. 다만 이 때의 각도는 PP와 원점 OOxx축이 이루는 각도가 아니다. PP에서 그은 수직선(수평선)이 원 x2+y2=a2(x2+y2=b2)x^{2}+y^{2}=a^{2}(x^{2}+y^{2}=b^{2})와 만나는 점이 OOxx축과 이루는 각도이다. 아래의 그림을 보자. 타원 x2a2+y2b2=1(a>b>0){\textstyle \frac{x^{2}}{a^{2}}}+\frac{y^{2}}{b^{2}}=1 (a>b>0)과 두 원 x2+y2=a2x^{2}+y^{2}=a^{2}, x2+y2=b2x^{2}+y^{2}=b^{2}이 그려져있다.

5F20FF7A0.png

위 그림을 통해 x=acosθx=a\cos \theta임은 자명하다. 이를 타원의 방정식에 대입해보면 다음과 같다.

a2cos2θa2+y2b2=1    cos2θ+y2b2=1 \begin{align*} &&\frac{a^{2}\cos^{2}\theta}{a^{2}}+\frac{y^{2}}{b^{2}}&=1 \\ \implies &&\cos ^{2}\theta+\frac{y^{2}}{b^{2}}&=1 \end{align*}

따라서 y=bsinθy=b\sin \theta. 그러면 타원의 둘레 EE를 아래와 같이 계산할 수 있다.

E=02π(dxdθ)2+(dydθ)2dθ=02π(asinθ)2+(bcosθ)2dθ=02πa2sin2θ+b2cos2θdθ=02πa2a2cos2θ+b2cos2θ=02πa1a2b2a2cos2θdθ \begin{align*} E &= \int _{0} ^{2\pi} \sqrt{ \left( \frac{dx}{d\theta} \right)^{2} + \left( \frac{dy}{d\theta} \right)^{2}}d\theta \\ &= \int _{0} ^{2\pi} \sqrt{ \left( -a\sin \theta \right)^{2} + \left( b\cos \theta \right)^{2}}d\theta \\ &=\int _{0} ^{2\pi} \sqrt{ a^{2}\sin^{2}\theta+ b^{2}\cos^{2} \theta } d\theta \\ &=\int _{0} ^{2\pi} \sqrt{ a^{2}-a^{2}\cos^{2}\theta+ b^{2}\cos^{2} \theta } \\ &=\int _{0} ^{2\pi} a\sqrt{ 1-\frac{a^{2}-b^{2}}{a^{2}}\cos^{2} \theta } d\theta \end{align*}

이때 a2b2a2(b<a)\sqrt{\frac{a^{2}-b^{2}}{a^{2}}}(b<a)는 타원의 이심률이며 ϵ\epsilon이나 kk로 표기한다. 1사분면의 둘레만 계산해서 4배하면 둘레가 나오므로, 최종적으로 타원의 둘레는 다음과 같다.

E=4a0π21k2cos2θdθ,k2=a2b2a2 E=4a\int _{0} ^{{\textstyle \frac{\pi}{2}}} \sqrt{ 1-k^{2}\cos^{2} \theta } d\theta,\quad k^{2}=\frac{a^{2}-b^{2}}{a^{2}}

이때 만약 타원이 yy축 쪽으로 더 길어서 b>a>0b>a>0인 조건이라면 아래 그림과 같다.

2.png

그러면 타원의 둘레는 다음과 같다.

E=4b0π21k2sin2θdθ,k2=b2a2b2 E=4b\int _{0} ^{{\textstyle \frac{\pi}{2}}} \sqrt{ 1-k^{2}\sin^{2} \theta } d\theta ,\quad k^{2}=\frac{b^{2}-a^{2} }{b^{2}}

제2 종 타원 적분

이때 위의 적분을 특별히 제2 종 타원 적분, 혹은 완전 제2 종 타원적분이라 부르고 아래와 같이 표기한다.

E(k)=0π21k2sin2θdθ E(k) = \int_{0}^{{\textstyle \frac{\pi}{2}}}\sqrt{ 1-k^{2}\sin^{2} \theta } d\theta

계산

위의 타원 적분은 초등함수로 나타낼 수 없고 수치적인 계산으로만 구할 수 있다. kk의 값에 따른 적분값을 구하면 아래와 같다.

untitled.png b=1.1547b=1.1547이고 a=1a=1이면 k=0.5k=0.5이다. E(0.5)=1.351E(0.5)=1.351이므로 타원 x2+y21.15472=1 x^{2}+\frac{ y^{2}}{1.1547^{2}}=1의 둘레는 다음과 같다.

4bE(k)=4×1.1547×1.351=6.239 4bE(k)=4\times 1.1547 \times 1.351=6.239


  1. 위키에는 유도 과정이 없고 나무위키는 틀렸다. ↩︎