Mathematica计算Pi的值.doc

 PAGE 5 Mathematica计算的值 姓名: 学号: 班级: 实验目的 学习使用Mathematica软件的一些基本功能计算的值，以下通过三种不同的方法求解： 1.数值积分法 2.泰勒级数法 3.蒙特卡洛（Monte Carlo）方法 实验的基本原理和方法 1.Mathematica中常用绘图函数Plot在绘制高次函数时的方法； 2.计算圆周率的数值积分法、泰勒级数法、蒙特卡洛（Monte Carlo）方法，并且利用特定的公式来计算圆周率 。 实验的内容和步骤 （1）数值积分法计算 半径为1的圆称为单位圆，它的面积等于。只要计算出单位圆的面积，就算出了。在坐标轴上画出以圆点为圆心，以1为半径的单位圆（如下图），则这个单位圆在第一象限的部分是一个扇形，而且面积是单位圆的1/4，于是，我们只要算出此扇形的面积，便可以计算出。 Mathematica输入如下： Plot[{4(1-x*x)},{x,0,1}] 图1 在计算扇形面积时，很容易想到使用数学分析中积分的方法，第一象限中的扇形由曲线及两条坐标轴围成，实际操作中，我们不能准确地计算它的面积，于是就通过分割的方法，将其划分为许多小的梯形，通过利用梯形的面积近似于扇形面积来计算。 利用Mathematica编程计算上式: 运行结果如下： 图2 从而得到 的近似值为3.14159265358979323846264338328，可以看出，用这种方法计算所得到的值是相当精确的。n越大，计算出来的扇形面积的近似值就越接近的准确值。 2.泰勒级数法计算 反正切函数的泰勒级数 计算，实验运行如下： 从实验过程可以看出，这种方法花费的时间很长。原因是当x=1时得到的的展开式收敛太慢。要使泰勒级数收敛得快，容易想到，应当使x的绝对值小于1，最好是远比1小。例如，因为 ，所以我们可以计算出的值，从而得到的值。这样，就使得收敛速度加快。改进后可以看出，泰勒级数法得到的结果比数值分析法精确到小数点后更多位。 （3）蒙特卡洛法计算 在数值分析法中，我们利用求单位圆的1/4面积来得到 ，从而得到 。单位圆的1/4是一个扇形，它是边长为1的??位正方形的一部分，单位正方形的面积。只要能够求出扇形的面积S在正方形的面积中所占的比例，就能立即得到S，从而得到的值。下面的问题归结为如何求k的值，这就用到了一种利用随机数来解决此种问题的蒙特卡洛方法，其原理就是在正方形中随机的投入很多点，是所投的每个点落在正方形中每一个位置的机会均等，看其中有多少个点落在扇形内。降落在扇形内的点的个数 与所投店的总数 的比可以近似的作为 的近似值。利用Mathematica编程如下： 从运行结果来看，蒙特卡洛法的计算结果为3.136，虽然精确度不太高，但运行时间短，在很多场合下，特别是在对精确度要求不高的情况下很有用的。 （4）利用麦琴给出，推出=4()。对比以上方法，这种简单的直接用公式求的的方法要简单得多，所以用处更广。 实验结果和结果分析： 实验过程中通过对于不同方法计算的值，对比比较，可以从中发现不同算法对计算结果及其运算精度，运行效率的影响，并不断改进方法，选择最优解法。实验达到了预期的结果。