数值积分方法-Read.pptVIP

* * 数学实验 怎样计算? 的值 ? 圆周率是人类获得的最古老的数学概念之一，早在大约3700年前（即公元前1700年左右）的古埃及人就已经在 用256/81（约3.1605）作为π的近似值了。几千年来，人们一直没有停止过求π的努力。 你也许能写出 ? = 3.1415926535 实际计算问题 8628034825342117068 但是你会计算?的值吗？你又能用几种 用Mathematica容易求出?到几百位 In[1] =N[Pi,100] Out[2] 3.1415926535897932384626433832795028841971 6939937510582097494459230781640628620899 方法计算?的值？ 古 典 方 法 分 析 方 法 数值积分方法 概 率 方 法 古典方法 用什么方法来计 算π的近似值呢？显然，不可能仅根据圆周率的定义，用圆的周长去除以直径。起先，人们采用的都是用圆内接正多边形和圆外切正多边形来逼近的古典方法。 6边形 12边形 24边形 圆 ? 阿基米德曾用圆内接 96边形和圆外切96边形夹逼的方法证明了 由 和 导出 ? 公元5世纪，祖冲之指出 比西方得到同样结果几乎早了1000年 ? 十五世纪中叶，阿尔·卡西给出π的16位小数，打破了祖冲之的纪录 ? 1579年,韦达证明 ? 1630年,最后一位用古典方法求π的人格林伯格也只求到了π的第39位小数 从十七世纪中叶起，人们开始用更先进的分析方法来求π的近似值，其中应用的主要工具是收敛的无穷乘积和无穷级数，在本节中我们将介绍一些用此类方法求π近似值的实例。 分析方法 利用幂级数表达式 取 取 ? 1656年，沃里斯(Wallis)证明 ? 在微积分中我们学过泰勒级数，其中有 当x = 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 左边三个正方形组成的矩形中， 由 和 可得 A C B D 问题: 能不能算得更快一点、更精确一点？ Machin公式 简单公式 麦琴(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=50; 将[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; 交充要条件为 从Buffon落针实验谈起： 平行线距离为a，针长度为l(a); 设针中心到较近平行线距离为x， 针与平行线夹角? , 针与平行线相 概率方法 Monte Carlo 法 设计方案1 次数很大,落针应均匀分布，落针中心在A区个数与 (左图正弦曲线下方面积D) 总数之比为A面积与矩形总面积之比=2