复合求积公式-上机实践报告.docVIP

中国矿业大学（北京）理学院 数值分析实验报告 实验名称 不动点迭代法求方程的近似根 实验时间 5月8日 组长签名 龙纯鹏 班级 信息与计算科学（1）班 学号 1110座机电话号码 成绩 组员签名 1110座机电话号码 1110座机电话号码 1110座机电话号码 1110座机电话号码 1110座机电话号码 一、实验目的,内容 二、相关背景知识介绍 三、代码 四、数值结果 五、计算结果的分析 六、计算中出现的问题，解决方法及体会 一、实验目的、内容 实验目的：取不同的步长h，分别用复合梯形及复合辛普森求积公式计算积分，通过这个实验清楚地认识到在同样的等分份数下，复合辛普森公式的近似程度明显优于复合梯形公式。 内容：分别用复化梯形公式和复化Simpson公式计算下列积分 1 ，等分数n 自己取定；并与准确值比较精度并估计误差； 2 ，等分数n 自己取定。 相关背景知识介绍 （1）算法原理或计算公式 : 设将区间[a,b]划分为n等份，步长,选取等距节点构造出的插值型求积公式，称为牛顿—柯特斯公式。由于牛顿—柯特斯公式在n8时具有不稳定性，故不可能通过提高阶的方法来提高求积精度。为了提高精度通常可把积分区间分成若干个子区间，再在子区间上用低阶求积公式。 当n 1时，就是我们熟悉的梯形公式 ，在每个子区间 k 0,1,```,n-1 采用梯形公式，则得 即复合梯形公式为 当n 2时，就是辛普森公式如下： 在每个子区间上采用辛普森公式就得： 即复合辛普森公式为 . 代码（Matlab） syms x f x*sqar 1+x^2 int f, 0, pi*2 -2*pi clear a 0; b 1; n 100; h b-a /n; t 0; for k 1:1:n-1; x k a+k*h; t t+g x k ; end t h/2* g a +g b +2*t ; t 0.7468 clear a 0; b 2*pi; n 100; h b-a /n; t 0; for k 1:1:n-1; x k*2 a+k*h; x k*2+1 x k*2 +h/2; t t+4*p x k*2+1 +2*p x k*2 ; end t h/6* p a +p b +t ; t t -6.2832 clear a 0; b 2*pi; n 100; h b-a /n; t 0; for k 1:1:n-1; x k*2 a+k*h; x k*2+1 x k*2 +h/2; t t+4*l x k*2+1 +2*l x k*2 ; end t h/6* l a +l b +t ; t t 0.8444 四、步骤及数值结果 -6.2832 0.8444 五、计算结果的分析 利用梯形公式的余项公式，可得复合梯形公式的截断误差为：,梯形公式的精确度比较低，收敛也比较慢，因此，梯形公式并不直接用来计算积分，而是为其它的积分法 如龙贝格积分法 提供初始数据，在那里，由梯形公式得出的这些不够准确的近似值，将被一些简单的运算加工后变得非常准确。 六、 计算中出现的问题，解决方法及体 由于在实际计算时，不宜使用高阶的牛顿——柯特斯公式，但若积分区间较大，单独用一个低阶的牛顿——柯特斯公式来计算积分的近似值，显然精度不好，为了提高数值求积的精确度，可利用积分对区间的可加性来解决这个问题，这就是通常采用的复合求积法。所谓复合求积法，其指导思想就是先将积分区间分成若干个小区间，在每个小区间上采用低阶求积公式进行计算，然后把所有子区间的计算结果相加得出新的求积公式，这种公式就叫做复合求积公式。 教 师 评 语 指导教师： 2013 年 5 月 8 日 4 x1 a+k*h f1 f1+f x1 k++ k n? Y t1 h/2* f a +2*f1+f b r1 fabs -4.0/9-t1 N J 1 x2 a+j* 1.0/2 *h J%2 0? N Y f2 f2+f x2 f3 f3+f x2 t2 h/6* f a +4*f3+2*f2+f b r2 fabs -4.0/9-t2 j 2n? Y 输出t1,r1 输出 t2,r2 N