南 京 师 范 大 学 泰 州 学 院
毕 业 论 文（设 计）
（ 一 三 届）
题 目： 漫谈圆周率π值的计算 院（系、部）： 数学科学与应用学院 专 业： 数学与应用数学 姓 名： 学 号 指导教师：
南京师范大学泰州学院教务处 制
南京师范大学泰州学院本科毕业论文
摘要：圆周率?在数学中是一个非常重要的常数，受到广泛关注。古今中外一代代的数学家为计算?献出了自己的智慧和劳动，人类对?值的认识过程，反映了数学和计算技术发展情形的一个侧面。本文首先研究圆周率计算的发展历程，然后给出了利用计算机求解圆周率的三种算法的基本原理。借助Mathematica软件编程给出了每一种算法的计算结果和误差分析，利用Matlab软件对所得数据进行了分析和比较，对三种算法的优缺点进行了讨论，最后阐述了从计算圆周率的过程中得到的启示。
关键词： 圆周率；计算；近似值；数学实验
Abstract: ? is one of the most important constant in mathematics. It has been intensively studied recently. In this paper, we first discuss the development process of the calculating of algorithms to calculate
??,and then we give the basic principles of three
by computer. Using Mathematica and Matlab, we
give the results of the three algorithms and discuss the advantages and disadvantages. Finally, we get some inspiration from the above research.
Keywords: ?；calculation；approximation；mathematical experiment
1
南京师范大学泰州学院本科毕业论文
目 录
1 绪论 ..................................................... 3
1.1研究意义 .......................................................... 3
1.2国内外研究现状 .................................................... 4 1.3本文的研究方法和主要解决问题 ...................................... 4
2 圆周率简介和圆周率计算的四个时期........................... 5
2.1圆周率的简史及其重要性 ............................................ 5 2.2圆周率计算的四个时期 .............................................. 6
2.2.1无算法记录时期 ............................................... 6 2.2.2几何推算时期 ................................................. 7 2.2.3解析计算时期 ................................................. 9 2.2.4计算机运算时期 .............................................. 10
3借助计算机求解圆周率的方法................................ 12
3.1数值积分法 ....................................................... 12 3.1.1 算法原理 ................................................... 12 3.1.2 计算结果及误差分析 ......................................... 14 3.2泰勒级数法 ....................................................... 17
3.2.1 算法原理 ................................................... 17 3.2.2 计算结果及误差分析 ......................................... 18 3.3蒙特卡洛法 ....................................................... 21
3.3.1正方形内投点法 .............................................. 21 3.3.2蒲丰投针法 .................................................. 25 3.3.3 随机整数互素法 ............................................. 27
4从圆周率计算中得到的启示 ................................. 30 谢 辞..................................................... 31 参考文献 .................................................. 32 附 录..................................................... 33
2
南京师范大学泰州学院本科毕业论文
1 绪论
我们知道，平面上圆的周长与直径之比是一个常数，称为圆周率，记作?。在日常生活中，人们经常与?打交道，?的计算伴随着人类的进步而发展，许多数学家在其计算上发费了巨大的精力。有些数学家甚至说：“历史上一个国家所算得的圆周率的准确程度，可以作为衡量一个国家当时数学发展的一面旗帜[1]”。我国伟大的科学家祖冲之（公元429—500年）在前人的基础上深入地研究了圆周率，经过长期坚持不懈的努力，求出了当时世界上最好的近似值，他利用割圆术，求出了精确到小数点后七位数字的圆周率，并且明确指出了圆周率的取值在3.1415926和3.1415927之间，这一举世瞩目的成就在世界上领先了一千多年。在很长一个时期里，计算的?值是数学上一件重要的事情。本课题将从圆周率的简史和重要性开始，了解人们为什么这么执着的计算着圆周率。接着分析圆周率计算的四个时期，更清楚的知道人们为探索圆周率花费的心血，使我们当代大学生受益匪浅，启发我们不仅要学习前人的数学思维方式，更要学习他们孜孜不倦，开拓进取，为科学奋斗终身的精神。最后，利用我们所学的数学知识，运用数学实验方法，结合积分，迭代，随机试验等，在计算机的帮助下，介绍了三种计算圆周率的快捷方法：数值积分法，泰勒级数法和蒙特卡罗法。借助Mathematica软件编程给出了每一种算法的计算结果和误差，利用Matlab软件对所得数据进行了分析，对三种算法的优缺点进行了讨论，
1.1研究意义
作为数学上的一个重要常数, ?不仅用于圆的计算, 而且也在很多的公式中出现, 就我们现在的中学数学教材来说, 数学中的初等几何、高中的立体几何、代数中的三角函数、统计学等等都要用到?, 它是我们最熟悉的无理数；在物理学科中也有很多的公式要用到?, 比如单摆周期T的公式、库仑研究的两个带电质点的相互作用力的公式中也有它的身影, 还有其他的科学分支中也要用到?, 在科学史上有重要的地位。同时从其发展史可以看出在计算圆周率的过程中用到了极限的概念、微积分的思想、概率统计的理论, 我们还要有实数的理论,除了这些以外, 还要靠数学和科学技术的发展, 它展示了数学思想、方法的发展历程, 在寻求圆周率的计算过程中也发现了很多的问题, 从而推动了数学的发展, 促使人们不断为提高计算速度而寻找新的计算方法、改进计算的手段。正如前文中提到的“历史上一个国家所算得的圆周率的准确程度, 可以作为衡量这个国家当时数学发展的一个标志。”因此, 圆周率的发展历史从一个侧面反映了数学的发展历史尤其是算法的发展史, 而且还是计算机科学的发展史, 代表了当时的计算机科学的水平[2]。
几千年来作为数学家们的奋斗目标，古今中外一代代的数学家为此献出了自己的智慧和劳动。在中国有刘徽、祖冲之等，在国外有阿基米德、卡西等。近代科学家如华罗庚、严士健等在其数论论文中也对圆周率问题进行了探讨。古往今来，从未有哪一个数学常数能向圆周率那样吸引众多的学者。圆周率在各个时期的文明中都像一颗闪耀的明珠，它往往能够在一定程度上折射出该文明数学发展的水平[3]。为求得圆周率的值，人
3