祖冲之与圆周率的计算.docVIP

祖冲之和圆周率的计算 所谓“圆周率”是指一个圆的周长与其直径的比值。古今中外，许多人致力于圆周率的研究与计算。为了计算出圆周率的越来越好的近似值，一代代的数学家为这个神秘的数贡献了无数的时间与心血。……，一直到正三○七二边形，算出圆周率等于三点一四一六。 祖冲之（公元429－500年），中国南北朝时著名的数学家、天文学家22∕7、密率：355∕113。 以前人计算圆周率，是要探究圆周率是否循环小数。自从1761年Lambert证明了圆周率是无理数，1882年Lindemann证明了圆周率是超越数后，圆周率的神秘面纱就被揭开了。现在人计算圆周率多是为了验证计算机的计算能力。 古人计算圆周率，一般是用割圆法。这种基于几何的算法计算量大，速度慢，吃力不讨好。随着数学的发展，数学家们在进行数学研究时有意无意地发现了许多计算圆周率的公式。英国天文学教授John Machin于1706年发现Machin公式，他利用这个公式计算到了100位的圆周率。还有很多类似于Machin公式的反正切公式。Machin公式： 分析其中的arctgx公式可以知道，这是一个级数公式，而在程序设计中则可以用一个“累加器”算法来实现。 用流程图来表现，则在流程图中，必定有判别框，并根据判别条件成立与否分别设置了重复部分操作内容的分支流程。 二、算法的程序实现 为了实现这个算法，则需要编制相应的程序，在程序中除了需要用到赋值语句、输入输出语句、其它计算语句外，还必须用到循环语句。 范例：我使用VB来编写程序实现这个算法。 算法中用到了一条输入语句、两个循环语句、一个输出语句以及多个赋值语句。 （1）建立窗体和输入、输出、命令按钮组件对象。 （2）编写“Command1”触发的程序代码。 在“Private Sub command1_click ”和“End Sub”之间输入以下的程序代码。 Dim i As Integer, n As Integer, pi As Double, arc1 As Double, arc2 As Double, x As Single i Text1.Text x 1 / 5 arc1 0 For n 1 To I arc1 arc1 + -1 ^ n - 1 * x ^ 2 * n - 1 / 2 * n - 1 Next x 1 / 239 arc2 0 For n 1 To I arc2 arc2 + -1 ^ n - 1 * x ^ 2 * n - 1 / 2 * n - 1 Next pi 16 * arc1 - 4 * arc2 Label1.Caption pi 第一行，定义了两个整数类型的数值变量I和n，一个单精度浮点数变量x，以及三个双精度浮点数变量pi、arc1和arc2。其中pi用于表示圆周率的值。 第二行，将text1文本框中的数据转换为整型数值并赋值给整型变量i.。 第三行，将x赋值为“1/5”。 第四行，将arc1赋值为“0”。因为arc1是一个乘加器，所以其初值应该是0。 第五行，表示开始一个循环，循环变量n从1开始，步长为1，依次取值到I，一共循环I次。 第六行，arc1 arc1 + -1 ^ n - 1 * x ^ 2 * n - 1 / 2 * n - 1 ， 这是一个累加器的算法，它将变量arc1的原值加上表达式的值，然后将加法运算的结果重新赋值给变量arc1作为arc1的新值。 Machin公式中的级数代数式，转换成表达式则为： -1 ^ n - 1 * x ^ 2 * n - 1 / 2 * n - 1 第七行。NEXT，表示循环变量n增加一个步长的值1，然后判断“n I”是否成立，如果成立则继续循环，否则不再循环直接执行下一个语句。 第八行，将x赋值为“1/239”。“Loop While I n”表示当I n成立时继续循环，从第五行“Do”的下面一行继续执行。如果I n不成立，即I比n大时，则不再循环，直接执行下一行即第九行的语句，从而结束循环。 第九行，将arc2赋值为“0”。Arc2也是一个乘加器。 第十~十二行，通过循环计算arc2的值。 第十三行，pi 16 * arc1 - 4 * arc2。通过Machin公式计算圆周率的值。 第十四行，将表示圆周率的变量pi的值赋值给“Label1”组件对象的“Caption”属性，输出圆周率。 （3）运行程序。 将第一个文本框中的“Text1”删除，重新输入“9”；然后单击“Command1”，就能在原来“Label1”的位置上输出pi的值“3．1座机电话号码8682”。输入“10”时，输出pi的值“3.1座机电话