关于π的研究.docx

关于π的研究 第一部分：π的计算。 1 泰勒级数法 利用反正切函数的泰勒级数的特例麦克劳林级数： arctanx=x-x33+x55…+-1k-1x2k-12k-1 计算π。 将x=1代入上式可以得到： π=4arctan1=4(1-13+15-……) 以上这个无穷级数收敛太慢，不实用，若使其收敛的快一点，可令-1x1.例如arctan15就收敛的较快。通过a=arctan15的值利用倍角公式可得tan2a=512,tan4a=120119.应当注意tan4a约等于1.故4a≈π4。但这还不够准确，应当算出误差b=4a-π4。由tan4a=120119和tanπ4=1得： tanb=tan(4a-π4)=1239.b=arctan1239.故π=16arctan15-arctan1239 2数值积分法 0141+x2dx=4(arctan1-arctan0)=π。可得到π的值。 计算定积分s=abfxdx,也就是计算曲线y=f(x)与直线y=0,x=a,x=b,所围成曲边梯形的面积。用平行于y轴的直线将该曲边梯形平分成n个小曲边梯形这样总面积就等于小曲边梯形的面积之和若n的值取得越大，使每个小曲边梯形的宽度都很小，此时可以将它上方边界看作直线段。将每个小曲边梯形近似地当作梯形求面积，具体如下：设分点为x1……xi将区间n等分。即xi=a+i(b-a)n所有曲边梯形的宽度都是h=b-an,yi=f(xi)则第i个曲边梯形的面积si≈12(yi-1+yi)h.将所有的曲边梯形的面积加一块儿得 s≈i=1n12(yi-1+yi)h 计算得到s≈b-an[y1+y2+…yn-1+y0+yn2]这就是梯形公式 由此可得π的值。 3.BBP计算方法。 π=0∞116n(48n+1-28n+4-18n+5-18n+6) 它打破了传统的圆周率的算法，可以计算圆周率的任意第n位，而不用计算前面的n-1位。这为圆周率的分布式计算提供了可行性。 4.割圆法 设一半径为1的圆，作这个圆的内接正n边形，用此正n边形的周长去近似圆的周长。显然当n→∞时，正n边形的周长就无限趋近于圆周长，求得正n边形周长后除以直径便求出了圆周率。 从几何上观察，可知：正n边形周长随n递增而递增，但始终是个有限值。割法如图1： 设圆半径为1，令半弦长AB=2a，AC=2c，OG和OD分别是等腰△OAB和△OAC的中线。则我们要做的只是求出c关于a的表达式c=c（a）.令GC=b，根据勾股定理有: 进而有 得到此式后，编写计算机程序就很容易了，C语言程序如下： #include stdio.h #include math.h main() { double a,b,c,d,pi; double sqrt(double); int i,j,n; a=0.5; b=0; c=0; d=0.5; scanf(%d,n); for(i=1;i=n;i++) { b=sqrt(1-a*a); c=(1-b)*0.5; d=sqrt(c); a=d; } j=pow(2,n)*3; pi=2*d*j; printf(%d\n,j); printf(%f\n,pi); } 这里有一个问题就是a的初值如何选择？显然越简单直观越好，而已知对于圆内接正六边形的每一条边长等于圆的半径。所以取a=0.5，程序中参数n是对正六边形分割的次数，d的作用是当输入n=0（正六边形）的时候，得到π=3，此所谓的“径圆一三”。将这个文件保存为文本，在linux下用“gcc -lm”命令编译后，打开编译后得到的文件就能执行。 第二部分：π的历史 关于π的历史可分为三个阶段：第一阶段为微积分出现前，第二阶段为微积分出现后计算机出现前，第三阶段为计算机出现后。 第一阶段，大都从几何角度出发，主要是割圆法。关于割圆术, 中西方稍有不同, 但实质一样。中国主要是由圆内接正多边形求π值, 西方则内外夹攻, 用算术平均给出π值。另外在割圆术的应用上,也因人而异, 有的求圆面积, 有的求周长, 有的求面积比等等, 异途同归, 最终都给出π的近似值。中国的割圆术主要以刘徽为先驱,。刘徽给《九章算术》作注时, 首先对“ 周三径一”这句话怀疑, 他认为“ 周三”乃弓之与弦也”。并画图证明“周三”只是圆内接正六形的周长与直径之比, 接着他指出: “ 世传此法, 莫肯精核”。“学者踵古, 习其谬矣他取半径为一尺之圆, 从考虑面积入手, 由内接正6边形开始, 每次边数倍增, 增至1 92 边.得出3.1464625π3他在实际应用上π=3.14. 在古代很多国家的数学家都用过割圆术, 如元前250年在《圆的角量》中, 着眼于周长, 内接外切两侧夹攻, 由正形增至96 边 得π=3.1418.晚