实验 怎样计算圆周率.pptVIP

实验 怎样计算圆周率 实验2 怎样计算圆周率 一、实验目的 本实验主要涉及数值积分法、Taylor级数法和Monte Carlo法，要求尝试利用所知道的数学知识来进行近似计算。 二、实际问题 众所周知，圆周率 是平面上圆的周长与直径之比，它等于3.141 592 6…。古代人把3作为它的近似值。古希腊阿基米德（Archimedes）曾得到 。我国宋代的祖冲之得到 的近似值 （约率）和 （密率），后者化为小数后等于3.141 592…，与 的准确值的误差在 以下。 但是，你是否知道怎样计算 的近似值？是否尝试过利用所知道的数学知识来计算 的近似值，自己当一回阿基米德或祖冲之？我们建议了下面几种方法。这些方法很简单，一学就会，你不妨在计算机上试一试，体验一下自己得到 值时的喜悦。相信你学会这些方法之后不会满足于只用来计算 ，而会利用它们来做别的事情。你也可以自己想出一些别的方法来计算 。当然，现在有更快更好的方法计算 ，可以计算出 的十万位以上的近似值。但这些方法就不那么简单了，不可能在这里介绍。如果有兴趣，可以从别的资料去了解。 但是，你是否知道怎样计算 的近似值？是否尝试过利用所知道的数学知识来计算 的近似值，自己当一回阿基米德或祖冲之？我们建议了下面几种方法。这些方法很简单，一学就会，你不妨在计算机上试一试，体验一下自己得到 值时的喜悦。相信你学会这些方法之后不会满足于只用来计算 ，而会利用它们来做别的事情。你也可以自己想出一些别的方法来计算 。当然，现在有更快更好的方法计算 ，可以计算出 的十万位以上的近似值。但这些方法就不那么简单了，不可能在这里介绍。如果有兴趣，可以从别的资料去了解。 三、预备知识 1. 数值积分法 定积分 计算出这个积分的数值，也就得到了 的值。 一般地，要计算定积分 ，也就是计算曲线 与直线 所围成的曲边梯形T的面积。 为此，用一组平行于y轴的直线  将曲边梯形T分成 n个小曲边梯形，总面积S分成这些小曲边梯形的面积之和。如果取n很大，使每个小曲边梯形的宽度都很小，可以将它上方的边界 近似看成直线段，将每个小曲边梯形近似地当作梯形来求面积，就得到梯形公式。如果更准确些，将每个小曲边梯形的上边界近似地看作抛物线段，就得到辛普森公式。具体公式如下： 梯形公式：设分点 将积分区间 n等分，即 ，所有的曲边梯形的宽度都是 ，记 ，则第i个曲边梯形的面积 近似地等于梯形面积 ，将所有这些梯形面积加起来就得到 这就是梯形公式。 Simpson公式：仍用分点 将积分区间 分成n等分，再作每个小区间  的中点 ，将第i个小曲边梯形的上边界近似地看作经过三点 的抛物线段，则可求得 其中 ，于是得到 这就是Simpson公式。 2. Taylor级数法 利用反正切函数的Taylor级数 计算 。 将 代入上面的级数可以得到 这似乎可以用来计算 。但是，这个无穷级数 应当比1小，最好是远比1小。 不实用。要使Taylor级数收敛快， 收敛太慢， 比如， 得较快。但算出它的值对计算 有何帮助呢？ 就收敛 也就是说： 与 有什么关系呢？ 我们有 ，由正切的倍角公式可得 ， ， 。可以用 作 这还不够准确。我们把误差 求出来对 进行修正。由 及 可求得 的近似值。但 ， 于是得到 这称为Maqin公式，利用 的Taylor展开式 ， 的近似值，就可以由Maqin 由于是交错级数，当| x | 1时，可以简单地用 来估计截断误差。 求出 公式求出 的近似值。 Taylor级数是无穷级数，实际计算时必然 只能取它的前n项，导致截断误差 3. Mo