实验4数值积分.docxVIP

《数值分析》实验报告 学号  班级  信科  142 班  姓名  韦昕玥 【实验课题】  数值积分 【实验目标】 1、 了解数 分法的意 ； 掌握常用的数 分方法 (特 是梯形法、 Simpson 方法、 Cotes 公式、 Romberg 算法以及 Gauss 求 公式 )的原理。 2、 学会 制复合梯形公式，复合 Simpson 公式， Romberg 算法的程序。 3、 通 体会常用数 分方法的逐步精致化 程。 【理论概述与算法描述】 b n h [ f ( x 1. ) f (x )] 复合梯形公式f (x) dx a k 1 2 k 1 k b h n 1 n 1 2. 4 f (xk 1 ) 2 f ( xk 1) f (b)] 复合 Simpson 公式f ( x)dx [ f (a) a 6 k 1 2 k 0 Cotes 公式 b a 12 f ( x2 ) 32 f (x3) 7 f ( x4 ) b-a C 7 f (x0 ) 32 f ( x1 ) 这里 x=a+kh,h= 90 4 1 f ( x) d x 0 A (f 0 x) 1 A (f1 )x 其中 x0 , x1 高斯点 4. 两点 Gauss 公式 1 Romberg 算法 1. 取 h0 b a 算 T0 h0 f (a) f (b) 2 2. k = 1, 2, (k ) 1 ( k 1) h0 2k 1 f a 2i 1 h0 ? 算下列各步 T0 2 T0 k k 2 i 1 2 3. n = 0, 1, 2, ?k ,= n –1, n –2, ?Tn(k ) 4n Tn(k1 1) Tn( k1) 4n 1 4. 若 Tk (0) Tk (0)1 ， 出 分 Tk(0) ，否 3 【实验问题】 算下列定 分 （1） I 2 sin xdx (精确解： I 1). 0 （2） I 1 sin x (取 f (0) 1 , 参考解： I 0.9460831 ). 0 dx x 【实验过程与结果】 1、 将区间四等分，依次使用复合梯形公式、复合 Simpson 公式、 Cotes 公式、两点 Gauss 公式计算问题（ 1），给出理论误差估计，并比较实际计算误差。 计算值 理论误差 实际计算误差 复合梯形公式 z1 =0.9871 ≤ 0.0202 0.0129 复合 Simpson 公式 z2 =1.0001 ≤ 0.0002 0.0001 Cotes 公式 z3 =0.99999156547 ≤ 0.000012 0.00000843453 两点 Gauss 公式 z4=0.9985 ≤ 0.0022 0.0015 2、 编制复合梯形公式， Simpson 公式的程序，要求对给定精度 tol 10 6 ，终止分割 (见 simpsonT2N.m Simpson.m) 3、 编制 Romberg 算法的程序，要求对给定精度 tol 10 6 ，终止分割。用上述三种 算法分别计算问题（ 1）（2）中的定积分，并比较计算量。 【结果分析、讨论与结论】 经过编程对复合梯形公式、 复合 Simpson 公式、 Cotes 公式、两点算法的实现， 分别算出积分值， 再通过理论的余项公式计算出理论误差，际计算误差是吻合的。  Gauss 公式、 Romberg 发现其情况是与实 【附程序】