1、圆周率圆周率,一般以来表示,是一个在数学及物理学普遍存在的数学常数。它定义为圆形之周长与直径之比。它也等于圆形之面积与半径平方之比。是精确计算圆周长、圆面积、球体积等几何形状的关键。分析学上,可定义为是最小的x0 使得 sin(x)=0。数学简史2常用的 近以值包括疏率:22/7 及密率:355/113。这两项均由祖冲之给出。约等于(精确到小数点后第100位)3.14159 26535 89793 23846 26433 83279 50288 41971 69399 37510 58209 74944 59230 78164 06286 20899 86280 34825 34211 706
2、80数学简史3古希腊欧几里得的几何原本(约公元前3世纪初)中提到圆周率是常数,中国古算书周髀算经(约公元前2世纪)中有径一而周三的记载,也认为圆周率是常数。历史上曾采用过圆周率的多种近似值,早期大都是通过实验而得到的结果,如古埃及纸草书(约公元前1700)中取=(4/3)43.1604。数学简史4第一个用科学方法寻求圆周率数值的人是阿基米德,他在圆的度量(公元前3世纪)中用圆内接和外切正多边形的周长确定圆周长的上下界,从正六边形开始,逐次加倍计算到正96边形,得到(3+(10/71)(3+(1/7),开创了圆周率计算的几何方法(亦称古典方法,或阿基米德方法),得出精确到小数点后两位的值。数学简
3、史5中国数学家刘徽在注释九章算术时(公元263年)只用圆内接正多边形就求得的近似值,也得出精确到两位小数的值,他的方法被后人称为割圆术,其中有求极限的思想。数学简史6南北朝时代的数学家祖冲之利用割圆术进一步得出精确到小数点后7位的值(公元466年),给出不足近似值3.1415926和过剩近似值3.1415927,还得到两个近似分数值,密率355/113和约率22/7,这一纪录在世界上保持了一千年之久。为纪念祖冲之对中国圆周率发展的贡献,将这一推算值用他的名字被命名为“祖冲之圆周率”,简称“祖率”。数学简史7其中的密率在西方直到1573才由德国人奥托得到,1625年发表于荷兰工程师安托尼斯的著作
4、中,欧洲称之为安托尼斯率。阿拉伯数学家卡西在15世纪初求得圆周率17位精确小数值,打破祖冲之保持近千年的纪录。德国数学家柯伦于1596年将值算到20位小数值,后投入毕生精力,于1610年算到小数后35位数,该数值被用他的名字称为鲁道夫数。数学简史8除的数值计算外,它的性质探讨也吸引了众多数学家。1761年瑞士数学家兰伯特第一个证明是无理数。1794年法国数学家勒让德又证明了2也是无理数。到1882年德国数学家林德曼首次证明了是超越数,由此否定了困惑人们两千多年的化圆为方尺规作图问题。还有人对的特征及与其它数字的联系进行研究。如1929年苏联数学家格尔丰德证明了e 是超越数等等。数学简史9在历史
5、上在历史上,有不少数学家都对有不少数学家都对圆周率作出过研究圆周率作出过研究,当中著名当中著名的有的有阿基米德、托勒密、张衡、阿基米德、托勒密、张衡、祖冲之祖冲之等。他们在自己的国家等。他们在自己的国家用各自的方法用各自的方法,辛辛苦苦地去辛辛苦苦地去计算圆周率的值。下面计算圆周率的值。下面,就是就是世上各个地方对圆周率的研究世上各个地方对圆周率的研究成果。成果。数学简史10研究圆周率历史的几个阶段起承转接数学简史11起起【起】即为圆周率的起源,那究竟是谁先发现它?古巴比伦人从计算周界发现:一块出土于 1936 年的黏土块上记载,在古巴比伦时期(约公元前 1900-1600 年),巴比伦人相信
6、六边形的周界为0;57,36(以底数 60 计,亦即=96/100=24/25)乘以它的外接圆的周界:六边形周界=24/25 其外接圆周界=24/25 直径由此,得出相信是最古老的圆周率的近似值:巴比伦巴比伦=25/8=3.125=25/8=3.125数学简史12承承【承】是承继安提丰和布赖森的穷举法而发展的一个时期:以多边形找寻圆周率的值数学简史13古希腊西那库斯的阿基米德(Archimedes of Syracuse,公元前 287-212 年),是第一个有系统地找出圆周率的近似值和圆周率的上下限的数学家。他采用了安提丰和布赖森的穷举法,但他的研究重点则在多边形的周界。阿基米德在圆的度量(
7、The Measurement of the Circle)中,提出三个有关圆的定理。即:3.14084.3.14285.数学简史14刘徽是独立开创以多边形面积迫近圆面积的穷举法割圆术来找出圆周率的值的。最后,刘徽更求得正 3072 边形的面积,从而得出:=3927/1250=3.1416 即 的值准确至小数后三个位,后人称为徽率。数学简史15祖冲之运用了刘徽的割圆术及他无比的耐性与坚持(当时并没有算盘等计算工具,只能靠小竹子帮助计算,但他实质的计算方法则无从确定),算到:3.1415926 3.1415927他还发现了约率:祖冲之更取 =22/7(=3.14.)作为约率密率:=355/113
8、(=3.1415929)作为密率,以表示圆周率的近似值。祖率:是圆周率的值准确至小数后 7 个位,后称3.1415926。数学简史16转转【转】是寻求圆周率的一个转折点。圆周率的计算有了新的突破以解析表达式表示及求出圆周率的值。数学简史17接接接是紧接着以上发现的很多计算圆周率值的公式所延伸的一个时期:随着科技的突飞猛进,计算机的发明,令圆周率的计算速度有了新的突破。数学简史18圆周率的研究方法古人计算圆周率,一般是用割圆法。即用圆的内接或外切正多边形来逼近圆的周长。Archimedes用正96边形得到圆周率小数点后3位的精度;刘徽用正3072边形得到5位精度;Ludolph Van Ceul
9、en用正262边形得到了35位精度。这种基于几何的算法计算量大,速度慢,吃力不讨好。随着数学的发展,数学家们在进行数学研究时有意无意地发现了许多计算圆周率的公式。下面挑选一些经典的常用公式加以介绍。除了这些经典公式外,还有很多其它公式和由这些经典公式衍生出来的公式,就不一一列举了。数学简史191、Machin公式:这个公式由英国天文学教授John Machin于1706年发现。他利用这个公式计算到了100位的圆周率。Machin公式每计算一项可以得到1.4位的十进制精度。因为它的计算过程中被乘数和被除数都不大于长整数,所以可以很容易地在计算机上编程实现。数学简史202、Ramanujan公式:
10、1914年,印度数学家Srinivasa Ramanujan在他的论文里发表了一系列共14条圆周率的计算公式,这是其中之一。这个公式每计算一项可以得到8位的十进制精度。1985年Gosper用这个公式计算到了圆周率的17,500,000位。数学简史213、AGM(Arithmetic-Geometric Mean)算法:Gauss-Legendre公式:初值:重复计算:最后计算:数学简史224、Borwein四次迭代式:初值:重复计算:最后计算:这个公式由Jonathan Borwein和Peter Borwein于1985年发表,它四次收敛于圆周率。数学简史235、Bailey-Borwei
11、n-Plouffe算法:这个公式简称BBP公式,由David Bailey,Peter Borwein和Simon Plouffe于1995年共同发表。它打破了传统的圆周率的算法,可以计算圆周率的任意第n位,而不用计算前面的n-1位。这为圆周率的分布式计算提供了可行性。1997年,Fabrice Bellard找到了一个比BBP快40的公式:数学简史24圆周率的新纪录圆周率的最新计算纪录由两位日本人Daisuke Takahashi和Yasumasa Kanada所创造。他们在日本东京大学的IT中心,以Gauss-Legendre算法编写程序,利用一台每秒可执行一万亿次浮点运算的超级计算机,从
12、日本时间1999年9月18日19:00:52起,计算了37小时21分04秒,得到了圆周率的206,158,430,208(3*236)位十进制精度,之后和他们于1999年6月27日以Borwein四次迭代式计算了46小时得到的结果相比,发现最后45位小数有差异,因此他们取小数点后206,158,430,000位的值为本次计算结果。这一结果打破了他们于1999年4月创造的68,719,470,000位的世界纪录。数学简史25人们追寻圆周率 的历史至今已有四千年,由发现圆周和直径的比为一常数,进而以多边形迫近圆的方法求 值,转成发现更多计算及表示 的公式、级数再随着电脑的发明与科技的发展,圆周率值的位数得以突飞猛进总结总结数学简史262022年10月1日数学简史27