圆周率的解析

1、圆周率的解析数学导报2006年2期上刊登了文1，该文认为“悠久的历史、奇妙的性质使得n蕴含着丰富的数学教育功能，然而在以往的数学教学中，这种功能却常常处于一种闲置状态，被当成一个冰冷的符号在使用。”并从“培养求真精神、激发好奇心、进行研究性学习”三 个方面展开探讨如何开发n中的数学教育功能。笔者读后，很受启发，同时也觉得如果充分利用信息技术(譬如用超级画板)来探索 n ,那就显得更加直观，完全可以达到可视化的效果。超级画板是一款智能化教学软件，功能相当强大，具有动态几何，符号运算，统计图表，编程环境，文稿演示等多种功能。它不仅能够满足中学数学教学和数学实验的几乎所有需求， 而且操作方式简单快捷

2、，易学易用，体现了人性化，智能化，可视化，动态化和程序化。本 文所有课件都是用免费版本的超级画板制作，该软件可在 下载，更多经典案例参见文2。下面以几个实例来说明如何利用超级画板探索n的性质。例1：圆内接正多边形周长与圆周长之间的关系圆周率n导常被定义为圆周长与直径的比，但圆是曲线，周长如何求呢？古代数学家很早就想出了穷竭法，将将圆的周长定义为圆内接正多边形周长的极限，设圆的周长为C,圆内接正多边形周长为P，则有lim巳-C。可若要严格证明，必然涉及“单调有界数列必有极限”这一抽象的极限知识，怎样才能让中学生比较直观而又深刻地理解这一点呢？利用 超级画板，简单的几步操作就能很好地解决这个问题。

3、(1) 点击【作图】|【文本作图】，在文本对话框中输入“Function(rho=3,0,2*pi,500,);Function(rho=3,0,2*pi,floor(n)+1,);Variable(n);”,点击运行则可生成圆和圆内接正多边形以及一把控制多边形边数的变量尺；(2) 拖动参数n使之约等于3，然后双击圆内接多边形，在对话框中勾选“折线段”。此时我们拖动参数 n使之从小到大变化，看到正多边形向圆“逼近”的过程(图1, 2,3),很自然就能体会到刘徽关于“割圆术”的论述 ：“割之弥细，所失弥少。割之又割，以至 于不可割，则与圆合体,而无所失矣。”我们可以进一步作出圆的外切多边形，利用

4、“两边夹”的思想来计算圆的周长。图1D.DD图2图3例2 :圆面积公式与周长公式之间的关系在小学就学过圆的面积公式 S和周长公式C = 2二r，但这两个公式有什么关系呢？求出圆的周长之后，如何来求圆的面积呢？我们可以用超级画板作出圆与矩形互变的动画：首先将圆分成若干等份(图 4),然后将圆弧展开(图 5),接着将上边部分平移一个三 角形的位置(图6)，最后将上边部分插入到下边部分(图 7)。我们直观地看到圆转化成矩C形(近似)的过程，很容易就能得到圆的面积S r = r r = -. r2。同时，我们可以得2出结论：等分份数越多，圆弧就越接近于直线，最后所得图形就越接近于矩形。图4图5图6图7

5、例3圆周率的微积分定义除了用圆周长与直径的比来定义圆周率n，还可以用单位圆的面积来定义 n。这是因为当r = 1时，圆的面积 s二二r2二二，于是又可以用4 j. 1-X2dx表示n。下面是用超级 画板探索j1 Ji-x2dx的过程。(1) 点击【作图】|【文本作图】，在文本对话框中输入“ Function(y=(1-xA2)/2),0,1, n,); Variable( n);”点击运行则可生成一段圆弧和一把变量尺；(2) 双击圆弧，在对话框中勾选“x轴区域”，并填充颜色；(3) 点击【作图】|【文本作图】，在文本对话框中输入“MeasureLowerArea(5); MeasureUppe

6、rArea(5); ”计算得到积分上和，积分下和；此时拖动参数n使之从小到大变化(图 8, 9),我们看到矩形向圆“逼近”，积分上和 与积分下和也越来越接近，导常取积分上和与积分下和的平均值来近似表示一，这也是微4积分中求曲多边形面积一般方法。图9例4利用反正切公式计算圆周率自1673年莱布尼茨用4arctan1计算n以后，数学家们又找到了大量的反正切公式来计算圆周率。在超级画板中，只要在【测量】|【测量表达式】中依次输入“4*arctan(1)”“ 16*arcta n(1/5)-4*arcta n(1/239)”“ 12*arctan(1/4)+4*arctan(1/20)+4*arcta

7、n(1/1985)” 就可很容易地算出结果(图 10 )。4arctan(l) = 3J-I159265 16arrtHii( -)-4 aid a n(吕亍)=3 1415926512arctan( -)+4arctan(丽)+4irctan( 抚= 3,141592&5图10例5利用无穷级数表达式计算圆周率如果你觉得按照例4的作法计算圆周率太直接，没什么看头的话，我们可以作出分步的动画。譬如利用欧拉公式二 1 二二二622 32 4212n来作动画。（1）在【测量】|【测量表达式】中测量 x，n和piA2/6的值；（2 ）单击右键，作出 n的动画，参数范围为 1到1,动画类型为一次

8、运动；再作出 x的 动画，参数范围为 0到0,动画类型为一次运动；（3）作出n的动画，参数范围为 n到n +1，动画类型为一次运动；再作出 x的动画，参 数范围为x到x+1/nA2，动画类型为一次运动；作完上述准备工作之后，先点击左边两个动画，使得n= 1, x = 0;然后在左边的程序区中，输入“ Move（10,11）; ”，接着反复按“CTRL+ENTER执行该语句。每执行一次，n就增加1, x的值就增加1/n2，图11是执行100次后的结果，可以预料继续执行下去，x的值2会越来越接近一。61> a HoueaOJH ; t AT兀类似地，我们可以利用: =2 22j44-61 3

9、密宙宙西1111 1 二- - -357911来作动412画，甚至利用连分数4 =1 一-2+匚 52 22。导过作这些动画，除了感慨n的形式各种各样之外，还将发现不同的无穷级数表达式的收敛速度是有很大差别的。例6利用概率实验计算圆周率随着新课程的改革，概率统计已经进入高中数学教材，而很多的实验如果按照原始的方法来做，需要花费很多宝贵的时间。借助于超级画板，就能使原本费时费力的实验过程变得 轻而易举。与其它动态几何软件相比，超级画板有着更多的数学函数，譬如其中的随机函数就使得此处的随机动画变成可能，统计表格功能更是为模拟概率统计方面的实验提供了方 便。根据几何概型的原理，我们可以用超级画板模拟

10、投豆实验（图12）。在正方形中随机取4k的结果就是圆周率的近似值。类似的方法还可以模拟蒲丰投针实验(图13）。蒲丰投针实验涮坦畑十71皿模鯉豆实验投i熬断-剜配图12饲内肚数k <S5lMt| jucantk- Ml K iHt-aii图13点，记录投豆的次数n和落在圆内次数k，当投豆次数足够多时，计算例7编程探讨圆周率的近似分数我国在圆周率的研究方面作出了杰出贡献,最为大家熟悉的是祖冲之将圆周率的值精确小的“住楼下”355就行了，另外“ 355113729215272827、& 9是连续的自然数,7+8= 15到小数点后7位，但对于祖冲之提出的密率有什么优越性，可能大家认识得还

11、不是很清楚。所谓密率，就是用355来作为n的近似值。密率极具创造性，简单易记，准确有趣。简单113易记是显然的，将1, 3, 5三个奇数各重复一遍后“平均”斩为两段，再让大的“住楼上”，也是密率有趣之处。那么密率的准确表现在哪里？很多数学家都对此问题作了论述，参见文3、4。其中文4用一个初等不等式非常巧妙地证明了若既约分数q比355更接近p 113圆周率二，则分母p 一定要比16586要大”，证明过程虽不复杂，但恐怕一般的数学爱好者 难以想到。下面我们就借助计算机的高速运算能力，用笨办法来解决。355“要找比355更接近于n的分数，而且分母尽可能小”转化成数学语言就是“找一对113正整数p,

12、q,使满足(1)n，也就是<|叫p|113P比栄更接近于；(2)要求p尽可能小”，进一步转化成超级画板的程序语言如下：d=z-pi; a=z;for(p=3;p<113;p=p+1)for(q=floor(p*(pi-d);q<floor(p*(pi+d)+1;q=q+1)if(q/p-pi)A2<(z-pi)A2)a=q/p; elsea;我们在程序区按“ CTRL+ENTER执行此程序，计算机返回“ >> (355)/(113) ”，这说明 在分母小于113的情况下，找不到符合要求的分数。假如我们将“ for(p=3;p<113;p=p+1)”一行

13、中的p=3改为p=113,p<113改为p<1000,运行之后，计算机仍然返回“>> (355)/(113) ”。52163355耐心地做下去，最终发现符合要求的分数是52163。进一步计算，我们就会感受到型 的16604113神奇所在：虽然n的另外一个渐近分数 366的简单程度与它差不多，但与n的绝对误差却106约是它的312倍；而绝对误差仅比它约小0.2 %的n另一渐近分数52163，却比它复杂得16604多。只有深刻感受密率的奥妙，我们才会发自内心敬重祖冲之这一伟大的数学家。利用超级画板，我们还可以作很多关于 n的探索，远不止本文所列的几个例子，感兴趣的朋友可以自行探索。笔者曾向齐民友先生演示超级画板，齐先生看后大发感慨，说现在很多数学教学