.怎样计算圆周率的值.pptVIP

* * 数学实验 怎样计算 ? 的值 ? 哪里有数,哪里就有美. - Proclus 知其然,更知其所以然. -中国先哲 本实验内容参考了中国科技大学 说明： 李尚志等编写的数学实验教材 你也许能写出 ? = 3.1415926535 实际问题 8628034825342117068 但是你会计算?的值吗？你又能用几种 ?―圆周率, 我们十分熟悉的常数. 用Mathematica容易求出?到几百位 In[1] N[Pi,100] Out[2] 3.1415926535897932384626433832795028841971 6939937510582097494459230781640628620899 方法计算?的值？ 方法1 利用幂级数表达式 用Mathematica计算 [Out8] 3.14154265358982449 In[1] k=1000; S1=N[4*Sum[(-1)^(n-1)/(2n-1),{n,1,k}],18] [Out2] 3.14059265383979293 In[3] k=10000; [Out4] 3.14149265359004324 In[5] k=15000; [Out6] 3.14152598692320065 In[7] k=20000 问题: 能不能算得更快一点、更精确 Machin公式 简单公式 一点？ 再用Mathematica Out[8] 3.141592653589793238462643383279 Clear[k,n,S] In[1] k=10; S2=N[4*Sum[(-1)^(n-1)*(1/2)^(2n-1)/(2n-1) +(-1)^(n-1)*(1/3)^(2n-1)/(2n-1),{n,1,k}],20] Out[2] 3.14159257960635121097 In[3] k=20; Out[4] 3.1415926535897574098 In[5] k=50; Out[6] 3.14159265358979323846 In[7] k=500; 方法2 将[0,1]区间n等分，取xk=k/n, 利用数值积分 yk= 1/ (1+xk2) Out[8] 3.141592646923126571795976843597 Mathematica In[1] y[x_]:=4/(1+x^2); n=100; S3=N[1/(2*n)*(Sum[2*y[k/n],{k,1,n-1}]+y[0]+y[1]),30] Out[2] 3.1415759869231285559229513739 In[3] n=500; Out[4] 3.141591986923126571922960843596 In[5] n=1000; Out[6] 3.141592486923126571797960843597 In[7] n=5000; 方法3 与平行线相交充要条件为 ? 线下方面积D) 次数很大,落针应均匀分布，落针中心在D与 Monte Carlo 法 从Buffon落针实验谈起： 平行线距离为1，针长度为1; 设针中心到较近平行线距离 为y ，针与平行线夹角? , 针 (左图正弦曲 总数之比为D面积与总面积之比为2/? o D 设计方案 ?＝4 m/n 计算机模拟：产生区间[0,1]上数目为n的一组 在正方形 0 x 1, 0 y1 上随机的投大量的点，那么 落在四分之一园内的点数 数m与在正方形内的点数n 之比m/n应为这两部分图形 面积之比＝?/4,故 随机数(x,y)，计算满足x2+y21的点数m Mathemetica Out[2] 3.14736 In[1] n=10000; S4= Block[{i,m=0}, For[i=n,i0,i--, m=m+If[Random[]^2+Random[]^2=1,1,0]]; N[4*m/n,10]] Out[2] 3.1352 In[1] n=50000; Out[2] 3.15336 In[1] n=100000;