实验二怎样计算Pi.docxVIP

数学实验实验报 学院：数学与统计学院 班级：数学与应用数学3班 学号：201370010314 姓名：康萍 时间：2016.04.05 实验二怎样计算兀 一、 实验目的 分别用下列三种方法计算兀的近似值，并比较三种方法的精确度： 数值积分法：通过使用Mathematica7. 0编写梯形公式和辛普森公式的程序语 言计算兀O 泰勒级数法：利用反正切函数泰勒级数计算帀。 蒙特卡罗(Monte Carlo)法：通过使用Mathematica7. 0编写蒙特卡罗公式 的程序语言来计算兀。 二、 实验环境 基于Windows环境下的Mathematica7. 0软件。 三、 实验的基本理论和方法 1、数值积分法 以单位圆的圆心为原点建立直角坐标系，则单位圆在第一象限内的部分G是 一个扇形，由曲线y = J1-兀2(兀丘［0,1］)及两条坐标轴用成，它的面积S = — o算 4 出了S的近似值，它的4倍就是帀的近似值。而扇形面积S实际上就是定积分 ￡^ll-x2dx =彳。 与兀冇关的定积分冇很多,比如 亠 的定积分(*亠必=4就比 肛卩的 1 + F 」。1 +兀 $ 4 定积分更容易计算，更适合于用来计算兀。 一般地，要计算定积分f f(x)dx，也就是计算曲线y = f\x)与直线 y = QyX = a^ = b所围成的曲边梯形G的面积S。为此，用一组平行于y轴的直线 兀二兀(1 in-\,a = x{) v兀］x2 ??? xz?_, xfJ = b)将曲边梯形T分成n个小曲边梯形，总而积S分成这些小曲边梯形的而积z和。如呆取n很大，使每个小曲边 梯形的宽度都很小，可以将它上方的边界/(啟也兀)近似的看作直线段， 将每个小曲边梯形近似的看作梯形来求而积，就得到梯形公式。如果更准确些， 将每个小曲边梯形的上边界近似的看作抛物线段，就得到辛普森公式。具体公式 如下： 梯形公式设分点西,???,￡“将积分区间［Q,b］分成n等份，即 兀=Q + i(b 一 d)/ /2,0 W i 。所有的曲边梯形的宽度都是力=0 — a)/川。记% = f(x{ ) 则第i个曲边梯形的面积S,近似的等于梯形面积i(yt._l+y^0将所有这些梯形 面积加起來就得到 S，— b n 必+力+…+仏+笃儿 这就是梯形公式。 辛普森公式 仍用分点Xz=6Z4-/—(1//7-1)将区间肚方］分成n等份，直 n 线x = x/(lz/?-l)将曲边梯形分成n个小曲边梯形。再做每个小区间［形十為］的 中点兀］/-2似的看作经过三点(x, /(x)、兀=兀1，兀1,齐/——2丿的抛物线段，则可求得S,6nX-i+4X 1其中 中点兀］ /-2 似的看作经过三点(x, /(x) 、 兀=兀1，兀1,齐 /—— 2丿 的抛物线段，则可求得 S, 6n X-i+4X 1 其中 I—— 2 于是得到 晋伉+儿)+2(刃+儿+…儿J+牡 +儿+???+〉 1 n— 2 2 这就是辛普森公式。 2、泰勒级数法 利用反止切函数的泰勒级数 r3 5 2—1 X X / . \jfc-i A arctan x = x H + (-l I + ?… 3 5 v 7 2k-\ 当X的绝对值小于1,最好是远小于1,这样，随着指数的增加，X的幕快速接 近于泰勒级数就会很快收敛,比如,取嗨得到的arctan*就收敛的快,在 arctan-L.1-1 “arctan*罟相差屮取2/2-1 = 63得到的arctan|的近似值的误差就小于准确度度已经非常高 To我们并不知道arctan- “arctan*罟相差 2 多少。IlIa - arctan= --a9 TOC \o 1-5 \h \z 4 / 、Q / 、 Q 71 tan p = tanl —-tz tan tancr I —— 4 2_ ,. 71 7 1 + tan — tan a 1 + lx — 4 2 因此0 = arctan —,即—-arctan — = arctan —，从而得到 4 2 3 (1)71 1 1 (1) ——=arctan — + arctan 一 4 2 3 arctan -比arctan-收敛得更快。利用泰?勒级数计算出arctan -与arctan -的近似 3 4 2 3 值再相加，然后再乘以4,就得到兀的近似值。 还可以考虑用= arctan-来计算兀，它收敛的更快。［ha=-易算出 5 120 tan 4a = 12 119 71 120 1 \ tan 4a 一 tan 1 71 4 一 119 4 , , 71 (120 r / 1 + tan 4a tan 1 + - 4 119 tan la tan 4a - 1 239 TOC \o 1-5 \h