• 首页
  • 中国
  • 放大和缩小的计算技巧:以确定圆周率为例

放大和缩小的计算技巧:以确定圆周率为例

2024-06-27 23浏览

手工计算时代需要基于数据精度考虑具体的计算方式,放大、缩小和求平均,就是技巧之一。

(1)Basel问题

自然数平方的倒数之和是收敛的,欧拉得到其*限为π2/6

张江敏.疯狂的绝技——级数加速收敛的艺术http://blog.sciencenet.cn/blog-100379-1075769.html

级数平方衰减,截断到前n项的和Tn,则误差在1/n量级。不过,若想结果精确到小数点后三位,并不需要计算到1000项——对没有计算的余项可进行估算呢。

若圆周率误差为δ,则(π+δ)2/6=π2/6+πδ/3+δ2/6,因而级数的误差和与π的误差δ相当。

由前10项的和T10确定的π为3.049362,修正1/(10+1)之后,则π为3.137525;

修正1/2(10+1)^2=0.004132后,得到π=3.141473

与真值3.141593,仅相差0.000120。

(2)割圆术

半径为1的圆,内接正6边形边长b6=1,

n边形边长bn,b2n=sqrt(2–sqrt(4–bn2));

n=6×2^k边形边长Ak=sqrt(2–Bk),式中

Bk=sqrt(2+sqrt(2+…+sqrt(2+sqrt3)…)),k–1次开方

正2n边形面积S2n=nbn/2,*限是圆周率π。

这就是魏人刘徽公元263年创立的割圆术,其以6×2^4=96边形的边长确定圆周率3.14。

通常认为,南朝祖冲之(429~500年)得到圆内接6×2^11=12288边形的边长,确定圆周率精确至小数点后7位。这是非常困难的:为了保证π的精度,B11=1.9999……需要16~17位有效数字。祖冲之所著《缀术》失传,具体计算方法难以知道;但他肯定知道有效数字随着边数的增加不断减少,进而研究提高精度的途径。

刘徽割圆4次至96边形的结果可以写为

S192=S12+(S24–S12)+(S48–S24)+(S96–S48)+(S192–S96)

S12=3;而**次割圆增加量即S24–S12=0.1058285412;第二次增加量减少,只是**次的0.253240倍,第三次、四次增加量是其前次增加量的0.250805和0.250201倍。

割圆引起的面积增加比值D逐次减小而趋于0.25即1/4。继续割圆的面积增加,若以1/4的比例计算,即弓形面积为(S192–S96)/3是缩小的估算;而以D96的比例计算,即弓形面积为(S192–S96)*D96/(1–D96)是放大的估算。于是有3.14159253<π<3.14159313。以上数据是Excel计算的,所列数字都是有效的。

如果以12位数字计算到192边形(k=5),即可得到祖率3.1415926<π<3.1415927,即比刘徽多割一次即可。我想,当年祖冲之就是这样做的。

又,若注意到割圆引起的面积增加比值与1/4的差异也以1/4的比率减小,因而圆周率会接近于下限,即对上下限求平均并不能提高精度;不过,若依据D96估计D192=0.250050进而计算(S384–S192),即可得到祖率而不必具体计算192边形的边长。

1/8=0.125,1/6=0.166667,而1/7=0.142857…,那么3又1/7即22/7就是圆周率的偏大估计。将1/7修改为P=1/(7+1/n)即n/(7n+1),其在1/8~1/7之间;选择合适的n可以提高精度。计算5次即可得到祖率355/113

n=10,P=0.1408451;n=20,P=0.1418440

n=15,P=0.1415094;n=17,P=0.1416667

*后可确定n=16,P=0.1415929

(3)此情可待成追忆

1994年初搬入两室一厅的新居,有了独立的书房。过年无处可去而翻看旧书,觉得祖冲之似乎不会割圆11次计算圆周率;于是,用十位数字的计算器Casiofx-4500p略作演算,写出“割圆术确定圆周率方法的改进——祖冲之确定圆周率过程之猜测”。敝帚自珍,多次投稿后发表于《安阳师范学院学报》的2003年第2期;也曾在博客介绍http://blog.sciencenet.cn/blog-275648-776036.html。

书上铅笔所写的结果已经不清楚啦。24年也就要过去了。

照片中小册子“夏道行.π和e”1964年第1版,但是1978年重印的,3万6千字,定价0.15元;“张弛.不等式”,5万8千字,定价0.22元。也就是二、三两猪肉的钱——其时含骨统肉一斤0.73元。

1978年暑假后坐船6小时到县城上高中,自己手中也就有了一点零花钱,时常去学校对门的新华书店看看,偶尔也买本书。我觉得读书比做题要好,只是那时候书籍很少啊。

*后给出当年发现的一个结果,只要记得sqrt(5)=2.236就行啊

附录:多次投稿也就是多次退稿啊。下面所展示的投稿 也是被退啊。二十年已经过去了。

从写作到发表用了九年时间。也花了点功夫呢,只是博文里说得简单。现在条件好啊,若做出题目,真不必再三再四投稿,完全可以在博客贴出来。有人看就行。

那时投稿还要到邮局寄纸质文本,有些刊物要审稿费,有些刊物说“6个月未见录用通知可自行处理”,有些刊物因“人力所限来稿概不退还”。文章手工抄写不易,且害怕重抄出错,*初投稿时多说“随信寄上贴足邮票、写好地址的回程信封;若拙稿不能采用,烦请掷回”。当然,