世界上精度最高的椭圆周长初等公式.docx
《世界上精度最高的椭圆周长初等公式.docx》由会员分享,可在线阅读,更多相关《世界上精度最高的椭圆周长初等公式.docx(9页珍藏版)》请在冰点文库上搜索。
世界上精度最高的椭圆周长初等公式
世界上精度最高的椭圆周长初等公式
D
或等于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,表明椭圆周长真值范围是
.
(表中上界值与下界值是笔者用更佳选配系数和次数计算的,因而比此文的上界公式与下界公式精度略高)
敬请各位专家、教授及同仁,提出宝贵意见和建议,在此提前谢过。
再现一次最美丽、精度最高的椭圆周长二级等比公式: