茉莉花新闻网

中華青年思想與行動的聚合地

椭圆没有周长?

平面椭圆,一个神奇的图形。

小时候的我,觉得椭圆就是一个普普通通的图形。直到我上了高中,接触了圆锥曲线,经历了一番摧残之后,我觉得我似乎认清了椭圆的真面目。而现在,我又碰到了椭圆积分,才发现我真的太天真了......所以,我现在对椭圆充满敬畏之情,不知道何时又会碰到与之有关的更为高深的知识。下面我们就从椭圆的周长开始,慢慢揭开椭圆积分的神秘面纱......


  1. 问题的引入:椭圆的周长

如果我没记错的话求平面曲线的弧长应该是导数的基本应用之一吧。首先,我们来回忆一下计算平面曲线弧长的公式。

设一连续可微的平面曲线

的参数方程为:

我们取这个曲线

上的一段微元并记作

。有勾股定理可得:

而:

带入

的表达式有:

两边同时积分可得:

这就是有关参数方程的弧长公式了。我们先小试牛刀,计算一下圆的周长。

我们知道,圆的标准参数方程是(其中

为圆的半径):

则:

equation?tex=%5Cleft%5C%7B%5Cbegin%7Bmatrix%7Dd+x%3D

代入弧长公式有:

equation?tex=s%3D%5Cint %7B0%7D%5E%7B2%5Cpi%7D%5Csqrt%7B%5Cleft%28+ Rsin%5Cleft+%28+t+%5Cright+%29%5Cright%29%5E%7B2%7D%2B%5Cleft%28+Rcos%5Cleft+%28+t+%5Cright+%29%5Cright%29%5E%7B2%7D%7Ddt

一切过程都十分顺利,那我们再来看看椭圆:

椭圆的标准参数方程大家肯定也不会陌生(其中

为半长轴长,

为半短轴长):

我们仅计算椭圆在第一象限的部分的弧长,之后在乘以

就好了。但是第一象限部分的参数的取值范围会有变化,即在第一象限中

。参数方程的导数为:

equation?tex=%5Cleft%5C%7B%5Cbegin%7Bmatrix%7Dd+x%3D

代入到弧长公式中得到:

equation?tex=s%3D4%5Cint %7B0%7D%5E%7B%5Cfrac%7B%5Cpi%7D%7B2%7D%7D%5Csqrt%7B%5Cleft%28 asin%5Cleft+%28+t+%5Cright+%29%5Cright%29%5E%7B2%7D%2B%5Cleft%28+bcos%5Cleft+%28+t+%5Cright+%29++%5Cright%29%5E%7B2%7D%7Ddt

直到现在,仍一切顺利,我们在化简一下看看:

......嗯?这玩意怎么处理?到这一步会发现根号完全去不掉,原函数也找不到。到此,本文结束。

 

嘿嘿,开个玩笑。聪明的数学家们是不可能就此罢休的,于是他们又开始将上面的式子进一步化简:

equation?tex=%3D4a%5Cint %7B0%7D%5E%7B%5Cfrac%7B%5Cpi%7D%7B2%7D%7D%5Csqrt%7B%5Cleft%28+1 +cos%5E%7B2%7D%5Cleft+%28+t+%5Cright+%29%5Cright%29%2B%5Cfrac%7Bb%5E%7B2%7D%7D%7Ba%5E%7B2%7D%7Dcos%5E%7B2%7D%5Cleft+%28+t+%5Cright+%29%7Ddt

equation?tex=%3D4a%5Cint %7B0%7D%5E%7B%5Cfrac%7B%5Cpi%7D%7B2%7D%7D%5Csqrt%7B1 %5Cleft%28+1 +%5Cfrac%7Bb%5E%7B2%7D%7D%7Ba%5E%7B2%7D%7D%5Cright%29cos%5E%7B2%7D%5Cleft+%28+t+%5Cright+%29%7Ddt

equation?tex=%5Cxrightarrow%5B++1 +%5Cfrac%7Bb%5E%7B2%7D%7D%7Ba%5E%7B2%7D%7D%3A%3D%5Cvarepsilon%5E%7B2%7D+%5D%7B+%5Ctau+%3D%5Cfrac%7B%5Cpi+%7D%7B2%7D t%7D4a%5Cint %7B%5Cfrac%7B%5Cpi%7D%7B2%7D%7D%5E%7B0%7D%5Csqrt%7B1 ++%5Cvarepsilon%5E%7B2%7D++cos%5E%7B2%7D%5Cleft+%28+%5Cfrac%7B%5Cpi+%7D%7B2%7D %5Ctau%5Cright+%29%7Dd%5Cleft%28++%5Cfrac%7B%5Cpi+%7D%7B2%7D %5Ctau%5Cright%29

equation?tex=%3D4a%5Cint %7B%5Cfrac%7B%5Cpi%7D%7B2%7D%7D%5E%7B0%7D%5Csqrt%7B1 ++%5Cvarepsilon%5E%7B2%7D++sin%5E%7B2%7D%5Cleft+%28+%5Ctau%5Cright+%29%7D%5Cleft%28++ d+%5Ctau%5Cright%29%3D4a%5Cint %7B0%7D%5E%7B%5Cfrac%7B%5Cpi%7D%7B2%7D%7D%5Csqrt%7B1 ++%5Cvarepsilon%5E%7B2%7D++sin%5E%7B2%7D%5Cleft+%28+%5Ctau%5Cright+%29%7Dd+%5Ctau%5Cleft%28+%5Cbigstar+%5Cright%29

其中: equation?tex=%5Csqrt%7B1 +%5Cfrac%7Bb%5E%7B2%7D%7D%7Ba%5E%7B2%7D%7D%7D%3A%3D%5Cvarepsilon 叫做椭圆的离心率

式还可做变量代换:

则:

equation?tex=dx%3A%3D+cos%5Cleft%28+%5Ctau+%5Cright%29d%5Ctau%3D%5Csqrt%7B1 sin%5Cleft%28+%5Ctau+%5Cright%29%5E%7B2%7D%7Dd%5Ctau%5Cxrightarrow%5B++%5D%7B+sin%5Cleft%28+%5Ctau+%5Cright%29%3Dx%7D%5Csqrt%7B1

equation?tex=%5CRightarrow+d%5Ctau%3D%5Cfrac%7Bdx%7D%7B%5Csqrt%7B1

则有变量代换后的积分:

equation?tex=4a%5Cint %7B0%7D%5E%7B1%7D%5Csqrt%7B1 %5Cvarepsilon%5E%7B2%7D+x%5E%7B2%7D%7Dd+%5Ctau%3D4a%5Cint %7B0%7D%5E%7B1%7D%5Cfrac%7B%5Csqrt%7B1 +%5Cvarepsilon%5E%7B2%7Dx%5E%7B2%7D%7D%7D%7B%5Csqrt%7B1

还可以写的更有强迫症一点:

equation?tex=4a%5Cint %7B0%7D%5E%7B1%7D%5Cfrac%7B%5Csqrt%7B1 +%5Cvarepsilon%5E%7B2%7Dx%5E%7B2%7D%7D%7D%7B%5Csqrt%7B1 x%5E%7B2%7D%7D%7Ddx%3D4a%5Cint %7B0%7D%5E%7B1%7D%5Cfrac%7B1 +%5Cvarepsilon%5E%7B2%7Dx%5E%7B2%7D%7D%7B%5Csqrt%7B1 +%5Cvarepsilon%5E%7B2%7Dx%5E%7B2%7D%7D%5Csqrt%7B1

到现在,椭圆积分的雏形已经出现了。

2. 椭圆积分的诞生

经过 等数学家的研究,椭圆积分的知识体系渐渐完善,直到 的出现彻底彻底完善了椭圆积分的知识体系。

我们先观察 式,这个式子是椭圆周长的积分公式,而它可以被拆成两部分:

equation?tex=4a%5Cint %7B0%7D%5E%7B1%7D%5Cfrac%7B1 +%5Cvarepsilon%5E%7B2%7Dx%5E%7B2%7D%7D%7B%5Csqrt%7B1 +%5Cvarepsilon%5E%7B2%7Dx%5E%7B2%7D%7D%5Csqrt%7B1 x%5E%7B2%7D%7D%7Ddx%3D4a%5Cint %7B0%7D%5E%7B1%7D%5Cfrac%7Bdx%7D%7B%5Csqrt%7B1 +%5Cvarepsilon%5E%7B2%7Dx%5E%7B2%7D%7D%5Csqrt%7B1 x%5E%7B2%7D%7D%7D 4a%5Cint %7B0%7D%5E%7B1%7D%5Cfrac%7B%5Cvarepsilon%5E%7B2%7Dx%5E%7B2%7D%7D%7B%5Csqrt%7B1 %5Cvarepsilon%5E%7B2%7Dx%5E%7B2%7D%7D%5Csqrt%7B1

我们将拆开后的第一部分拿出来,并去掉积分上下限和系数得到不定积分:

equation?tex=I %7B1%7D%3A%3D%5Cint%5Cfrac%7Bdx%7D%7B%5Csqrt%7B1 +%5Cvarepsilon%5E%7B2%7Dx%5E%7B2%7D%7D%5Csqrt%7B1

再将第二部分拿出来,去掉系数和积分上下限得到另一个不定积分:

equation?tex=I %7B2%7D%3A%3D%5Cint%5Cfrac%7Bx%5E%7B2%7D%7D%7B%5Csqrt%7B1 %5Cvarepsilon%5E%7B2%7Dx%5E%7B2%7D%7D%5Csqrt%7B1

另外还有个一个不定积分:

equation?tex=I %7B3%7D%3A%3D%5Cint%5Cfrac%7Bdx%7D%7B%5Cleft%28+x a+%5Cright%29%5Csqrt%7B1 +%5Cvarepsilon%5E%7B2%7Dx%5E%7B2%7D%7D%5Csqrt%7B1

这三个不定积分便是 所总结得到的。若将上面的三个不定积分做变量代换: 则:

equation?tex=%5Ctilde%7BI%7D %7B1%7D%3A%3DF%5Cleft%28+k%2C%5Cphi+%5Cright%29%3D%5Cint%5Cfrac%7Bdx%7D%7B%5Csqrt%7B1 +k%5E%7B2%7Dx%5E%7B2%7D%7D%5Csqrt%7B1 x%5E%7B2%7D%7D%7D%3D%5Cint %7B0%7D%5E%7B%5Cphi%7D%5Cfrac%7Bcos%5Cleft%28+%5Cphi%5Cright%29d%5Cphi%7D%7B%5Csqrt%7B1 +k%5E%7B2%7Dsin%5E%7B2%7D%5Cleft%28+%5Cphi%5Cright%29%7D%5Csqrt%7B1 sin%5E%7B2%7D%5Cleft%28+%5Cphi%5Cright%29%7D%7D%3D%5Cint %7B0%7D%5E%7B%5Cphi%7D%5Cfrac%7Bd%5Cphi%7D%7B%5Csqrt%7B1 +k%5E%7B2%7Dsin%5E%7B2%7D%5Cleft%28+%5Cphi%5Cright%29%7D%7D

equation?tex=%7BI%7B%7D%27%7D %7B2%7D%3A%3D%5Cint%5Cfrac%7Bx%5E%7B2%7D%7D%7B%5Csqrt%7B1 k%5E%7B2%7Dx%5E%7B2%7D%7D%5Csqrt%7B1 x%5E%7B2%7D%7D%7Ddx%3D%5Cint %7B0%7D%5E%7B%5Cphi%7D%5Cfrac%7Bsin%5E%7B2%7D%5Cleft%28+%5Cphi%5Cright%29cos%5Cleft%28+%5Cphi%5Cright%29%7D%7B%5Csqrt%7B1 k%5E%7B2%7Dsin%5E%7B2%7D%5Cleft%28+%5Cphi%5Cright%29%7D%5Csqrt%7B1 sin%5E%7B2%7D%5Cleft%28+%5Cphi%5Cright%29%7D%7Dd%5Cphi%3D%5Cint %7B0%7D%5E%7B%5Cphi%7D%5Cfrac%7Bsin%5E%7B2%7D%5Cleft%28+%5Cphi%5Cright%29%7D%7B%5Csqrt%7B1 k%5E%7B2%7Dsin%5E%7B2%7D%5Cleft%28+%5Cphi%5Cright%29%7D%7Dd%5Cphi

equation?tex=%5Ctilde%7BI%7D %7B2%7D%3D%5Ctilde%7BI%7D %7B1%7D k%5E%7B2%7D%7BI%7B%7D%27%7D %7B2%7D%3DE%5Cleft%28+k%2C%5Cphi+%5Cright%29%3D%5Cint %7B0%7D%5E%7B%5Cphi%7D%5Cfrac%7Bd%5Cphi%7D%7B%5Csqrt%7B1 k%5E%7B2%7Dsin%5E%7B2%7D%5Cleft%28+%5Cphi%5Cright%29%7D%7D %5Cint %7B0%7D%5E%7B%5Cphi%7D%5Cfrac%7Bk%5E%7B2%7Dsin%5E%7B2%7D%5Cleft%28+%5Cphi%5Cright%29%7D%7B%5Csqrt%7B1 k%5E%7B2%7Dsin%5E%7B2%7D%5Cleft%28+%5Cphi%5Cright%29%7D%7Dd%5Cphi%3D%5Cint %7B0%7D%5E%7B%5Cphi%7D%5Csqrt%7B1 k%5E%7B2%7Dsin%5E%7B2%7D%5Cleft%28+%5Cphi%5Cright%29%7Dd%5Cphi

equation?tex=%5Ctilde%7BI%7D %7B3%7D%3A%3D%5Cpi%5Cleft+%28n%2C+k%2C+%5Cphi+%5Cright+%29%3D%5Cint %7B0%7D%5E%7B%5Cphi%7D%5Cfrac%7Bd%5Cphi%7D%7B%5Cleft%28+1%2Bn%5E%7B2%7Dsin%5E%7B2%7D%5Cleft%28+%5Cphi+%5Cright%29+%5Cright%29%5Csqrt%7B1 k%5E%7B2%7Dsin%5E%7B2%7D%5Cleft%28+%5Cphi%5Cright%29%7D%7D (这个我不知道怎么来的...)

上面的 分别叫做 第一类,第二类,第三类椭圆积分。

之后 又定义了三类 椭圆积分,是将 椭圆积分里面的 换回 得到的,即:

equation?tex=F%5Cleft%28+k%2Cx+%5Cright%29%3D%5Cint %7B0%7D%5E%7Bx%7D%5Cfrac%7Bdx%7D%7B%5Csqrt%7B1 +k%5E%7B2%7Dx%5E%7B2%7D%7D%5Csqrt%7B1

equation?tex=E%5Cleft%28+k%2Cx%5Cright%29%3D%5Cint %7B0%7D%5E%7Bx%7D%5Cfrac%7Bx%5E%7B2%7D%7D%7B%5Csqrt%7B1 k%5E%7B2%7Dx%5E%7B2%7D%7D%5Csqrt%7B1

equation?tex=%5Cpi%5Cleft+%28n%2C+k%2C+x%5Cright+%29%3D%5Cint %7B0%7D%5E%7Bx%7D%5Cfrac%7Bdx%7D%7B%5Cleft%28+1%2Bn%5E%7B2%7Dx%5E%7B2%7D%5Cright%29%5Csqrt%7B1 k%5E%7B2%7Dx%5E%7B2%7D%7D%5Csqrt%7B1

参数 叫做椭圆积分的模。

特别的,当 时,这三类椭圆积分都称为完全椭圆积分,否则称为不完全椭圆积分,即:

完全 椭圆积分:

完全 椭圆积分:

3. 椭圆的周长公式

椭圆并非没有周长,只不过没有精确值罢了。对于其周长公式,是一个无穷级数的形式:

equation?tex=C%3D2%5Cpi+a%5Cleft+%28+1 %5Csum %7B+i%3D1%7D%5E%7B%5Cinfty%7D%5Cleft+%28+%5Cprod %7Bj%3D1%7D%5E%7Bi%7D%5Cfrac%7B2j 1%7D%7B2j%7D+%5Cright+%29%5E%7B2%7D%5Cfrac%7B%5Cvarepsilon%5E%7B2i%7D%7D%7B2i 1%7D+%5Cright+%29

其中: 为半长轴长, 为椭圆的离心率。这个级数是由第二类椭圆积分展开所得到的。(可惜我不会展开)。可见,当离心率为零时,级数退化为圆的周长公式。

当然,椭圆的周长公式有几个近似公式:

  • 利用算数平均值近似:(精度较低)

v2 88e0961fcb67f1484f62f8b73e34ac70 b
Fehler(误差)。误差与离心率和半短轴与半长轴之比的关系。图片来源:维基百科。
  • 利用平方均值近似:(精度一般)

v2 3f6d01d06c757546af43f205a0d1671a b
Fehler(误差)。误差与离心率和半短轴与半长轴之比的关系。图片来源:维基百科。
  • 近似公式:(精度很高)

equation?tex=C%5Capprox+%5Cpi+%5Cleft+%28+a%2Bb+%5Cright+%29%5Cleft+%28+1%2B+%5Cfrac%7B3%5Clambda+%5E%7B2%7D%7D%7B10%2B%5Csqrt%7B4 3%5Clambda+%5E%7B2%7D%7D%7D%5Cright+%29%2C%5Clambda+%3D%5Cfrac%7Ba b%7D%7Ba%2Bb%7D

v2 d19992289e5d187cc810433571172559 b
Fehler(误差),Bereich(区间)。误差与离心率之间的关系。图片来源:维基百科。
v2 d353fc881f5573ba3c20b25badceecf4 b

还有没关注专栏《数学及自然科学}》的朋友们请赶快关注吧!您的支持是我最大的动力

同类信息

查看全部

茉莉花论坛作为一个开放社区,允许您发表任何符合社区规定的文章和评论。

茉莉花新闻网

        中国茉莉花革命网始创于2011年2月20日,受阿拉伯之春的感召,大家共同组织、发起了中国茉莉花革命。后由数名义工无偿坚持至今,并发展成为广受翻墙网民欢迎的新闻聚合网站并提供论坛服务。

新闻汇总

邮件订阅

输入您的邮件地址:

linkedin facebook pinterest youtube rss twitter instagram facebook-blank rss-blank linkedin-blank pinterest youtube twitter instagram