世界上精度最高的椭圆周长初等公式.doc

世界上精度最高的椭圆周长初等公式

世界上精度最高的椭圆周长初等公式
成都七中高中远程教学 周钰承 15983417081

根据微积分基本原理，可以写出椭圆周长的定积分公式，但由于被积函数的原函数不是初等函数，所以椭圆周长没有标准的初等公式。但数学家们推导、证明了下面这个椭圆周长标准公式：
（1）
公式（2）中，。
这个公式表明，椭圆周长的主要部分为，我们可以把（1）中括号里从第二项起称为椭圆率多项式：
（2）
通常，我们要计算椭圆周长，必须先给出一个精确度。
假如要求我们误差率低于，我们设需要计算到椭圆率多项式第n项，不妨设，则椭圆率多项式（2）中，第n+1项及其以后无穷多项之和必须满足下列不等式：

因为（注意）：

所以只须：

（3）
取满足不等式（4）的最小整数。为此，我们需要一个带有函数的学生计算器，根据精确度要求，首先计算出我们应该计算到第几项。计算所得的椭圆周长值在给定误差率的情况下是精确的。注意：计算到椭圆率多项式第n项，就是标准公式（2）括号中算到2n次方项；若n为负数或者小于2，就算到椭圆率多项式（3）第2项，即公式（2）中括号里的4次方项。例如：n>-1.86745.则标准公式（2）中，中括号里应该算到4次方项。因为误差公式证明中n大于或等于2是前提条件。
需要知道的是，多数情况下求椭圆周长，只须计算到前两三项，因而往往可以笔算。
但是，当时，算得:,即用到椭圆率多项式第58项即116次方项，才能保证误差率小于万分之一。为此，我们可以构建一个新的函数模型，用以解决很小即很大时的计算问题。
我们把椭圆率多项式中的系数简化得：
（4）
观察（5），由于，所以随着增大而减小；各项系数逐渐变小，但与等比数列相比，“小得越来越慢”。根据（5）式的这些特点，我们构造一个多项式函数，使它与前三项相同，同时为了方便运算，我们从第二项起各项系数为等比数列：
（5）
变形为无穷等比数列求和（其系数从项开始为等比数列）：

用近似代替，代入标准公式（1）得

从而得到一个椭圆周长的近似公式：
（6）
公式（6）中，。这是我们在下一课时计算椭圆周长要用到的近似公式。为了突出这个公式，我们称（6）为椭圆周长一级等比公式。

近似公式如果没有误差估计是没有实际意义的。这个初等公式的精度如何呢？为此，我们介绍一下椭圆周长误差率定理：椭圆周长真值C满足下列不等式：
是椭圆周长的一个下界公式；
是椭圆周长的一个上界公式。
是椭圆周长二级等比公式。
限于篇幅，此处不给出详细证明过程，只对证明思路作简要介绍。
证明思路：上界公式与下界公式的级数展开式与椭圆标准公式相比，和它以前的所有系数完全相同，称第一部分。而从到的系数，称第二部分。以的系数为首项，公比为（100：144）的等比数列作为下界公式的第二部分，公比为（103：144）为等比数列作为上界公式的第二部分，可用完全归纳法证明：下界公式每项系数小于标准公式中相同次数的系数（仅有一项例外，但可用前一项系数补足），上界公式每项系数大于标准公式中相同次数的系数。从项开始称为第三部分。上界公式与下界公式中的项主要作用是保持与二级等比公式形式上的统一性，所以它们并不是最佳选配的系数和次数。可用数学归纳法证明上界公式和下界公式第三部分的正确性。
注意二级等比公式中，第二个分式的分母中为，不是40次方，而是4次方，它的出现是从到渐变过程中产生的。二级等比公式中系数和次数是最佳选配的，不可更易。
上界公式与下界公式均超过了目前所有的椭圆周长初等公式（包括中国椭圆周长公式）的精度；用上界公式与下界公式及两边夹定理，可以求出椭圆周长的精确值，这是上界公式与下界公式的主要优点，它能够让我们判断，我们用程序计算标准椭圆周长公式时，累加一百万项后精确度如何，特别是当接近于1时。
二级等比公式是比上界公式和下界公式精度还高得多的椭圆周长初等公式。笔者预言，这将是地球上精度最高的初等公式，永远不会再出现比这个公式更简洁、更美丽、更实用、精度更高的椭圆周长初等公式。我们可以怀疑用程序累加项名达公式五百万项的结果，但不可怀疑这个仅用学生计算器就可计算椭圆周长的二级等比公式。

一级等比C1
真值下界W1
两边夹得真值C
真值上界W2
二级等比C2
1
0.00
3.9924406639
3.9997223959
4.000（舍入值）
4.0000073192
4.0000258422
1
0.05
4.0171654582
4.0194052969
4.0194（不足值）
4.0194730943
4.0194213666
1
0.10
4.0631510073
4.0639698166
4.0640（过剩值）
4.0639850957
4.0639741647
1
0.15
4.1257797162
4.1260991688
4.12610（不足）
4.1261027235
4.1261002185
1
0.20
4.2018801742
4.2020085931
4.202009（舍入）
4.2020094578
4.2020088474
1
0.25
4.2891586239
4.2892107981
4.289211（舍入）
4.2892110135
4.2892108607
1
0.30
4.3858889139
4.3859100446
4.3859100（不）
4.3859100983
4.3859100601
1
0.35
4.4907311329
4.4907395661
4.4907396（过）
4.4907395793
4.4907395699
1
0.40
4.6026192371
4.6026225174
4.60262252（舍）
4.6026225205
4.6026225183
1
0.45
4.7206880143
4.7206892437
4.720689244舍
4.7206892444
4.7206892439
1
0.50
4.8442236721
4.8442241102
4.8442241102舍
4.8442241103
4.8442241102
1
0.55
4.9726292543
4.9726294004
4.9726294004舍
4.9726294005
4.9726294004
1
0.60
5.105399728
5.1053997727
5.1053997727
5.1053997727
5.1053997727
1
0.65
5.2421035874
5.2421035996
5.2421035996
5.2421035996
5.2421035996
1
0.70
5.3823689786
5.3823689815
5.3823689815
5.3823689815
5.3823689815
1
0.75
5.5258730396
5.5258730402
5.5258730402
5.5258730402
5.5258730402
1
0.80
5.6723335777
5.6723335778
5.6723335778
5.6723335778
5.6723335778
1
0.85
5.8215024802
5.8215024803
5.8215024803
5.8215024803
5.8215024803
1
0.90
5.9731604325
5.9731604325
5.9731604325
5.9731604325
5.9731604325
1
0.95
6.1271126366
6.1271126366
6.1271126366
6.1271126366
6.1271126366
1
1.00
6.2831853072
6.2831853072
6.2831853072
6.2831853072
6.2831853072

一级等比C1
真值下界W1
两边夹得真值
真值上界W2
二级等比C2

备注：舍入值若为4.202009，表明椭圆周长真值范围是；
不足值若为4.12610，表明椭圆周长真值范围是；
过剩值若为4.0640，表明椭圆周长真值范围是.
（表中上界值与下界值是笔者用更佳选配系数和次数计算的，因而比此文的上界公式与下界公式精度略高）
敬请各位专家、教授及同仁，提出宝贵意见和建议，在此提前谢过。再现一次最美丽、精度最高的椭圆周长二级等比公式：