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

巡山小妖精
517次浏览
2020年12月13日 03:50
最佳经验
本文由作者推荐

红豆简谱-浪漫情书

2020年12月13日发(作者:袁滋)


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

根据微积分基本原理,可以写出椭圆周长的定积分公式,但由于被积函数的原函数不
是初等函数,所以椭圆周长没有标准的初等公式。但数学家们推导、证明了下面这个椭圆周
长标 准公式:


1

2
2

1

2
4

13

2
6

5!!

2
8

C

(ab)

1











(1)

24

2 46

8!!



2



公式(2)中,


ab

ab
这个公式表明,椭圆周长的主要部分为

(ab)
,我们可以把(1)中括号里从 第二
项起称为
椭圆率多项式


1

1

4

13

6

5!!

8
f(

)


2










(2)
224246

8!!

通常,我们要计算椭圆周长,必 须先给出一个精确度。
假如要求我们误差率低于
d
,我们设需要计算到椭圆率多项式 第n项,不妨设
n2

则椭圆率多项式(2)中,第n+1项及其以后无穷多项之和 必须满足下列不等式:
2222

(2n1)!!

2n2< br>
(2n1)!!

2n4

(2n3)!!

2n6





d

(2n2)!!




(2n4)!!




(2n6)!!



 
因为(注意
n2
):
222

(2n1)! !

2n2

(2n1)!!

2n4
< br>(2n3)!!

2n6




< br>

(2n2)!!




(2n 4)!!




(2n6)!!





(2n1)!!

2n2

(2n1)!!

2n4

(2n1)!!

2n 6








(2n 2)!!




(2n2)!!




(2n2)!!




< br>
1

2n2

1

2n4

1

2n6













256

256

256

1

2n2
 
256
1

2
1

2n2
d< br> 所以只须:
256
1

2
222
222
222

2n2
256(1

2
)d

ln256(1

2
)d
n1
(3)
2l n




n
取满足不等式(4)的最小整数。为此,我们需要 一个带有函数的学生计算器,根
据精确度要求,首先计算出我们应该计算到第几项。计算所得的椭圆周长 值在给定误差率
d
的情况下是精确的。注意:计算到椭圆率多项式
f(
)
第n项,就是标准公式(2)括号中
算到2n次方项;若n为负数或者小于2,就算到椭 圆率多项式(3)第2项,即公式(2)
中括号里的4次方项。例如:n>-1.86745.则标准公 式(2)中,中括号里应该算到4次方项。
因为误差公式证明中n大于或等于2是前提条件。
需要知道的是,多数情况下求椭圆周长,只须计算到
f(

)
前两三项,因而 往往可
以笔算。
ln256(1

2
)

1 57.42
,但是,当

0.95,d0.0001
时,算得:
n
2ln

即用到椭圆率多项式第58项即116次方项,才能保证误差率小于万 分之一。为此,我们可
以构建一个新的函数模型,用以解决

b
很小即
很大时的计算问题。
a
我们把椭圆率多项式
f(

)
中的系数简化得:
1
2
1
4
1
6
25
8
49
10
441
12


3


4

7


8


10


(4)
4
44444
观察(5),由于
0

 1
,所以

2n
随着
n
增大而减小;各项系数逐渐变小,但
f(

)
与等比数列相比,“小得越来越慢”。根据(5)式的这些特点, 我们构造一个多项式函数
g(

)

使它与
f(

)
前三项相同,同时为了方便运算,我们从第二项起各项系数为等比数列:
g(
)
1
2
1
4
1
6
1
8< br>1
10
1
12


3


4


5


6


7


(5)
4
44444
4
变形为无穷等比数列求和( 其系数从

项开始为等比数列):
1
2
1
4
11 11


3

(1

2

2< br>
4

3

6

4

8< br>)
44
4444
111111


2

3

4
[1(

)
2
(

)
4
(

)
6
(

)
8
]
42222
4
111


2

3

4
(当n时
,
n
0)
1< br>4
4
1(

)
2
2

1
2

4



4
6416

2
(16

2
4

4
)

4< br>
6416

2
16

2
3

4

6416

2
g(

)


16

2
3

4

g(

)
近似代替
f(

)
,代入标准公式(1)得: 6416

2

16

2
3
< br>4

C

(ab)


1
6 416

2





从而得到一个椭圆周长的近似公式:

643

4
< br>C
1


(ab)


6416

2


(6)

公式(6)中,


ab
。这是我们在下一课时计算椭圆周长要用到的近似公式。
ab
为 了突出这个公式,我们称(6)为椭圆周长
一级等比公式


近似公式如果 没有误差估计是没有实际意义的。这个初等公式的精度如何呢?为此,
我们介绍一下椭圆周长
误 差率定理
:椭圆周长真值C满足下列不等式:

643

4
643

4
9

8
(14429

2


40
)

9

8(14432

2


40
)


(ab)


14
C

(ab)


142240

2240

6416
2(144100



)6416

2(144 103



)


643

4
9

8
(14432

2


40
)

是椭圆周长的一个
下界公式

W1


(ab)


142240

2(144100



)

6416


643

4
9

8
(14429< br>
2


40
)

是椭圆周长的一个
上界公式

W
2


(ab)


142240

2(144103



)
6416


643

4
9

8
(14432

2
13

40
)

是椭圆周长
二级等比公式

C
2


(ab)


14

224
2(1441 00



)

6416

限于篇幅 ,此处不给出详细证明过程,只对证明思路作简要介绍。
证明思路:上界公式与下界公式的级数展开式 与椭圆标准公式相比,

和它以前的所
有系数完全相同,称第一部分。而从



的系数,称第二部分。以

的系数为首项,
公比为(1 00:144)的等比数列作为下界公式的第二部分,公比为(103:144)为等比数列
作为上界公 式的第二部分,可用完全归纳法证明:下界公式每项系数小于标准公式中

相同
次数的 系数(仅有一项例外,但可用前一项系数补足),上界公式每项系数大于标准公式中

相同次数 的系数。从

项开始称为第三部分。上界公式与下界公式中的

项主要作用是
保持与二级等比公式形式上的统一性,所以它们并不是最佳选配的系数和次数。可用数学归
纳法 证明上界公式和下界公式第三部分的正确性。
注意二级等比公式中,第二个分式的分母中为


,不是40次方,而是4次方,它的
4
4040
12381010



渐变过程中产生的。出现是从
100


103
二级等比公式中系数和次数是最佳选配的,
不可更易。
上界公式 与下界公式均超过了目前所有的椭圆周长初等公式(包括中国椭圆周长公式)
的精度;用上界公式与下界 公式及两边夹定理,可以求出椭圆周长的精确值,这是上界公式
与下界公式的主要优点,它能够让我们判 断,我们用程序计算标准椭圆周长公式时,累加一
百万项后精确度如何,特别是当

接 近于1时。
二级等比公式是比上界公式和下界公式精度还高得多的椭圆周长初等公式。笔者预言,这将是地球上精度最高的初等公式,永远不会再出现比这个公式更简洁、更美丽、更实用、
精度更高 的椭圆周长初等公式。我们可以怀疑用程序累加项名达公式五百万项的结果,但不
可怀疑这个仅用学生计 算器就可计算椭圆周长的二级等比公式。
a

b

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

0.00 3.9924406639
3.9997223959
4.000(舍入值) 4. 4.
0.05 4.
4.
4.0194(不足值) 4. 4.
4.
4.0640(过剩值)0.10 4. 4. 4.
0.15 4.1257797162
4.1260991688
4.12610(不足) 4.1261027235 4.1261002185
0.20 4.2018801742
4.2020085931
4.202009(舍入) 4.2020094578 4.2020088474
0.25 4.2891586239
4.2892107981
4.289211(舍入) 4.2892110135 4.2892108607
0.30 4.3858889139
4.3859100446
4.3859100(不) 4.3859100983 4.3859100601
0.35 4.4907311329
4.4907395661
4.4907396(过) 4.4907395793 4.4907395699
0.40 4.6026192371
4.6026225174
4.60262252(舍) 4.6026225205 4.6026225183
0.45 4.7206880143
4.7206892437
4.720689244舍 4.7206892444 4.7206892439
0.50 4.8442236721
4.8442241102
4.8442241102舍 4.8442241103 4.8442241102
0.55 4.9726292543
4.9726294004
4.9726294004舍 4.9726294005 4.9726294004
0.60 5.105399728
5.1053997727
5.1053997727 5.1053997727 5.1053997727
5.2421035996 5.2421035996
5.3823689815 5.3823689815
5.5258730402 5.5258730402
5.6723335778 5.6723335778
5.8215024803 5.8215024803
5.9731604325 5.9731604325
6.1271126366 6.1271126366
6.2831853072 6.2831853072
真值上界W2 二级等比C2
0.65 5.2421035874
5.2421035996
5.2421035996
0.70 5.3823689786
5.3823689815
5.3823689815
0.75 5.5258730396
5.5258730402
5.5258730402
0.80 5.6723335777
5.6723335778
5.6723335778
0.85 5.8215024802
5.8215024803
5.8215024803
0.90 5.9731604325
5.9731604325
5.9731604325
0.95 6.1271126366
6.1271126366 6.1271126366
1.00 6.2831853072
6.2831853072 6.2831853072
2
2
b

一级等比C1 真值下界W1 两边夹得真值

备注:舍入值若为4.202009,表明椭圆周长真值范围是
4.2020085x4. 2020095

不足值若为4.12610,表明椭圆周长真值范围是
4.126100x4.12611

过剩值若为4.0640,表明椭圆周长真值范围是
4.06390x4.06400
.
(表中上界值与下界值是笔者用更佳选配系数和次数计算的,因而比此文的上界公式与下界
公式 精度略高)
敬请各位专家、教授及同仁,提出宝贵意见和建议,在此提前谢过。再现一次最美丽、精度最高的椭圆周长二级等比公式:



643

4< br>9

8
(14432

2
13

40
)


C
2


(ab)


14

224
2(144100


)

6416

湖北郧阳医学院-莫旗教育局


啤酒洗头的好处-金达莱花歌词


网店产品-有趣的实验


多情总为无情伤-我爱mm


parker钢笔-淘宝学


大宋皇帝列表-工资级别


美白乳液排行榜-黄土高原成因


最经典网名-上海牌照拍卖网