数学建模 π的计算 .pdfVIP

实验报告 实验目的 1. 掌握用MATLAB 软件求π值的方法，并对结果作初步分析； 2. 掌握用MATLAB 软件实现蒙特卡洛法的过程。 实验内容 运用数值积分法、泰勒级数法、模拟蒲丰实验和随机整数互素法 计算 π 值，取 n = （k=3, 4, …,7），记录所得结果和程序 运行时间，并作简要分析。 实验过程 一、数值积分法 单位圆的面积等于以单位圆的圆心为原点建立直角坐标系， 单位圆在第一象限内的部分G 是一个扇形，其面积 S = π/4 只需计算出 S 的近似值，它的4 倍就是 π 的近似值 1 、 梯形公式法 将扇形 G 分为 n 个宽度相等的部分 Gk （1 ≤ k ≤ n ） Gk ：曲边梯形 左右边界平行，上方边界为曲线 n→ ∞ Gk→梯 G G 形 Gk 梯形面积： 等价于 π 取 n =10k （k=3, 4, …,7 ），按上述方法，通过计算扇形面积计 算 π 的近似值。 程序 Clear tic n=1e3(3, 4, …,7); x=linspace(0,1,n+1); y=(1-x.^2).^0.5; pai=(0.5*y(1)+0.5*y(n+1)+sum(y(2:n)))*4/n toc 2、辛普森（Simpson ）公式法 仍用分点 xi = a + i(b-a)/n (1≤ i ≤ n-1) 将区间[a, b]分成 n 等份， 直线 x = xi (1≤ i ≤ n-1) 将曲边梯形分成 n 个小曲边梯形，再作 每个小区间[xi-1, xi]的中点 将第 i 个小曲边梯形的上边界 y=f(x) (xi-1≤x≤xi) 近似地看作经 过这三点的抛物线段，则可求得： 其中 于 是 得 到 S ; 即 程序： clear tic n=1e3(3,4,5,6,7); x=linspace(0,1,n+1); x2=x(1:n)+1/(2*n); y=4./(1+x.^2); y2=4./(1+x2.^2); pai=((y(1)+y(n+1))+2*sum(y(2:n))+4*sum(y2))/(6*n) toc 二、泰勒级数法 考虑反正切函数的泰勒级数 取 x =1 ，n =1e3(3,4,5,6,7) 程序： clear tic k=1e3(3,4,5,6,7); x=1; n=1:k; pai=sum((-1).^(n-1).*x.^(2.*n-1)./(2.*n-1)) toc 结果显示：花费时间长！准确度差！！ 原因：x =1 时得到的arctan1 的展开式收敛太慢！ Maqin 公式法 提高收敛速度： x 1 随着指数的增加，x 的幂快速接近于0 ， 泰勒级数就会快速收敛 取 2k-1 = 63 ，其误差已经非常小，表明收敛速度非常快 有 再考虑收敛更快的： /4=arctan 1/2+arctan 1/3 程序： clear tic k=1e3(3,4,5,6,7); n=1:k; x=1/2; at_1=sum((-1).^(n-1).*x.^(2.*n-1)./(2.*n-1)); x=1/3; at_2=sum((-1).^(n-1).*x.^(2.*n-1)./(2.*n-1)); pai_1=4*at_1+4*at_2 toc /4=4arctan 1/5-arctan 1/239 程序： clear tic k=1e3(3,4,5,6,7); n=1:k; x=1/5; at_5=sum((-1).^(n-1).