浅谈圆周率的由来和应用

  本文通过对圆周率各个时期由来的认识，深刻的理解到圆周率的历史价值，包括通过π找出各种表达式，通过 计算圆的面积和周长，一些函数的定义，积分的计算，指数的构成等都要用到 ；最后还介绍了圆周率在测试计算机性能上的应用和圆周率在C语言上的应用．

  很早以前,人们看出,圆的周长和直经的比是个与圆的大小无关的常数,并称之为圆周率.1600年,英国威廉.奥托兰特首先使用 表示圆周率,因为 是希腊之“圆周〞的第一个字母,而δ是“直径〞的第一个字母,当δ=1时,圆周率为 .1706年英国的琼斯首先使用 .1737年欧拉在其著作中使用 .后来被数学家广泛承受,一直没用至今.

  公元前150年左右,另一位古希腊数学家托勒密用弦表法(以1 的圆心角所对弦长乘以360再除以圆的直径)给出了 的近似值3.1416.

  公元200年间,我国数学家徽提供了求圆周率的科学方法----割圆术 〔如图2所示〕,表达了极限观点.徽与阿基米德的方法不一样,他只取“接〞不取“外切〞.利用圆面积不等式推出结果,起到了事半功倍的效果.而后,祖冲之在圆周率的计算上取得了世界领头羊,求得“约率〞和“密率〞(又称祖率)得到3.1415926 3.1415927.可惜,祖冲之的计算方式后来失传了.人们推测他用了徽的割圆术,但终究用什么方法,还是一个谜.

  例如，1777年，法国数学家蒲丰研究投针问题，将一根长为l的的针任意投到画有间距为 的平行线的平面上 ，他得到得结论是：该针与任一平行线相交的概率是 ，圆周率与随机现象产生了密切联系即 在概率中也有作用．在数学中还有一个重要公式 将圆周率与虚数单位 联系起来．

  它现在可以被人们用来测试或检验超级计算机的各项性能，特别是运算速度与计算过程的稳定性．这对计算机本身的改良至关重要．

  就在几年前，当 公司推出奔腾 时，发现它有一点小问题，这问题正是通过运行 的计算而找到的．这正是超高精度的 计算直到今天仍然有重要意义的原因之一．

  众所周知，圆周率一般以 来表示，是一个在数学及物理学都会存在的数学常数．它定义为圆形之周长与直径之比．它也等于圆形之面积与半径平方之比．是准确计算圆周长、圆面积、球体积等几何形状的关键值．圆周率是一个常数〔约等于3.1415926〕，是代表圆周长和直径的比例．它是一个无理数，即是一个无限不循环小数．圆周率在生产实践中应用十分普遍，在科学不很兴旺的古代，计算圆周率是一件相当复杂和困难的工作．俗话说得好，“有理走遍天下，无理寸步难行〞圆周率 好比这个“理〞．有了圆周率 不仅解决了困惑众多数学家的三大著名几何问题之一的化圆为方的不可能性更为后续的数学研究奠定了根底．因此，圆周率的理论和计算在某些特定的程度上反映了一个国家的数学水平．

  由于没初始化随机种子，所以该寒暑在执行过程中并不具有“随机性〞，由于在TC和在VC6.0中初始化随机种子的函数并不一样，本程序为保持兼容性，未参加那些容．

  在网上流传很久的“精简〞代码，如果你进入谷歌的主页，并且搜索“外星人计算 的程序〞这个词汇，很容易看到这样一个程序：

  人们试图从统计上得悉 的各位数字是否有某种规律.竞争还在继续,正如有人所说,数学家探索中的进程也像 这个数一样:永不循环,无止无休……

  1579年法国韦达发现了关系式，首次摆脱了几何学的旧方法，寻求到了 的解析表达式．1650年瓦里斯把 表示成无穷乘积，无穷连分数，无穷级数等各种值表达式纷纷出现，值计算精度也迅速增加．稍后，莱布尼茨发现接着欧拉证明了这些公式的计算量都很大．尽管形式格外的简单， 值的计算方式的最大突破是找到了它的反正切函数表达式．1706年英国数学家麦欣首先发现了其计算速度远超于方典算法．

  该程序当然违背了C语言的代码规，不过如果你运行它，你会惊奇的发现它得到了圆周率小数点后面800位数字．这段代码来自“国际C语言混乱代码大赛〞〔The International Obfuscated C Code Contest 〕．

  〔1〕圆周率现在可以被人们用来测试或检验超级计算机的各项性能，特别是运算速度与计算过程的稳定性；

  蒙特卡洛算法简单描述：在树枝积分法中，利用求单位圆的 的面积来求得 ，从而得到 ．单位圆的 面积是一个扇形，它是边长为1单位正方形的一局部．只要求出扇形面积 在正方形面积S中占的比例 就立即能得到 ，从而得到 的值．方法是在正方形中随机投入很多个点落在扇形．将落在扇形的点数m与所投点的总数n的比 作为k的近似值． 落在扇形的充要条件是 ．

  随机模拟方法就是这样一种计算圆周率的方法．虽然这种方法并不能“非常好〞的逼近圆周率，但是在其他领域，随机模拟有很大的作用．这种算法也称为蒙特卡洛算法．

  为了探索圆周率在各种方面上的广泛应用，通过网上以及书籍上的查找知道了圆周率各种各样的由来方法，还知道了可以由圆周率找到各种级数和函数等多种表达式，还包括通过 来计算圆的周长和面积，以及一些函数的定义、积分的计算、指数的构要用到 ，最终得到了圆周率在测试计算机性能上的应用以及圆周率在C语言方面中的应用等等．

  某个古代文牍员以不同长度的半径画了一些圆，他取了每个圆的直径〔将半径加倍〕仅仅是为了好玩．他决定以每个圆的直径为单位长度在圆周上丈量．令人惊奇的是，不管圆的大小如何，圆周总是直径的3倍多一点．由于 与圆的特殊关系，故数学家设计用来计算出圆的面积和周长的新方法．

  例 一个圆形花坛的直径是4米，沿它的外侧铺一条1米宽的小路，求这条小路的面积。〔准确到0.1平方米〕

  15世纪,的数学家阿尔.卡西通过分别计算圆接和外接正3 2 边形周长,把 值推到小数点后16位,打破了祖冲之保持了上千年的记录.

  本世纪50年代以后,圆周率 的计算开场借助于电子计算机,从而出现了新的突破.目前有人宣称已经把 计算到了亿位甚至十亿位以上的有效数字.

  这个简单的公式就可以用循环实现来实现圆周率的近似计算〔虽然它并不真正适用于计算圆周率〕，下面来分析一下．

  我们用C语言中的 类型的变量来存放最终的结果，而上面公式中的n作循环变量再适合不过了，在实际编码中我们取 从1变化到30000.由于 类型本身的精度限制， 其实也没有必要取太大的值．在循环中， 每次增加2，而要注意类型转换和在循环累加中项的符号变化，这里为了清楚起见，我们增加了变量s，由它产生了正负符号的变化．下面来看代码：

  C语言课程盗了循环章节，它的魅Baidu Nhomakorabea就逐渐显现了出来，许多小程序和算法都可以让学生去尝试时限．而我们在高等数学泰勒公式章节学习的圆周率公式这样一个时间段就派上了用场．

  中学时我们就已经知道： ，从而 ，如果我们应用泰勒公式将 展开，就可以得到

  是一个很重要的常数.一位德国数学家评论道：“历史上一个国家所算得的圆周率的准确程度,可以做为衡量这个这家当时数学开展水平的重要标志.〞古今中外很多数学家都孜孜不倦地寻求过 值〔如图1所示〕的计算方法.

  公元前200年间古希腊数学家阿基米德首先从理论上给出 值的正确求法.他用圆外切与接多边形的周长从大、小两个方向上同时逐步逼近圆的周长,巧妙地求得 ．