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

1、精选优质文档-倾情为你奉上精选优质文档-倾情为你奉上专心-专注-专业专心-专注-专业精选优质文档-倾情为你奉上专心-专注-专业 中国矿业大学（北京）理学院数值分析实验报告实验名称不动点迭代法求方程的近似根实验时间5月8日组长签名龙纯鹏班级信息与计算科学（1）班学号成绩组员签名一、实验目的,内容 二、相关背景知识介绍 三、代码四、数值结果 五、计算结果的分析 六、计算中出现的问题，解决方法及体会一、实验目的、内容实验目的：取不同的步长h，分别用复合梯形及复合辛普森求积公式计算积分，通过这个实验清楚地认识到在同样的等分份数下，复合辛普森公式的近似程度明显优于复合梯形公式。 内容：分别用复化梯形公式

2、和复化Simpson公式计算下列积分 (1) ，等分数n 自己取定；并与准确值比较精度并估计误差； (2) ，等分数n 自己取定。相关背景知识介绍（1）算法原理或计算公式 : 设将区间a,b划分为n等份，步长,选取等距节点构造出的插值型求积公式，称为牛顿柯特斯公式。由于牛顿柯特斯公式在n8时具有不稳定性，故不可能通过提高阶的方法来提高求积精度。为了提高精度通常可把积分区间分成若干个子区间，再在子区间上用低阶求积公式。当n=1时，就是我们熟悉的梯形公式 ，在每个子区间(k=0,1,n-1)采用梯形公式，则得即复合梯形公式为当n=2时，就是辛普森公式如下：在每个子区间上采用辛普森公式就得：即复合辛

3、普森公式为.代码（Matlab）syms x f = x*sqar(1+x2) int(f, 0, pi*2) -2*picleara=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);endt=h/2*(g(a)+g(b)+2*t);t =0.7468cleara=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);endt=h/6*(p(a)+p(b)+

4、t);tt = -6.2832cleara=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);endt=h/6*(l(a)+l(b)+t);t t =0.8444四、步骤及数值结果x1=a+k*hf1=f1+f(x1) k+x1=a+k*hf1=f1+f(x1) k+kn?kn?YYt1=h/2*(f(a)+2*f1+f(b)r1=fabs(-4.0/9-t1)Nt1=h/2*(f(a)+2*f1+f(b)r1=fabs(-4.0/

5、9-t1)NJ=1J=1x2=a+j*(1.0/2)*hx2=a+j*(1.0/2)*hJ%2=0?J%2=0?NYNYf2=f2+f(x2)f2=f2+f(x2)f3=f3+f(x2)f3=f3+f(x2)t2=h/6*(f(a)+4*f3+2*f2+f(b)r2=fabs(-4.0/9-t2)t2=h/6*(f(a)+4*f3+2*f2+f(b)r2=fabs(-4.0/9-t2)j2n?j2n?YY输出t1,r1输出t1,r1输出 t2,r2NN -6.2832 0.8444五、计算结果的分析利用梯形公式的余项公式，可得复合梯形公式的截断误差为：,梯形公式的精确度比较低，收敛也比较慢，因此，梯形公式并不直接用来计算积分，而是为其它的积分法(如龙贝格积分法)提供初始数据，在那里，由梯形公式得出的这些不够准确的近似值，将被一些简单的运算加工后变得非常准确。六、 计算中出现的问题，解决方法及体由于在实际计算时，不宜使用高阶的牛顿柯特斯公式，但若积分区间较大，单独用一个低阶的牛顿柯特斯公式来计算积分的近似值，显然精度不好，为了提高数值求积的精确度，可利用积分对区间的可加性来解决这个问题，这就是通常采用的复合求积法。所谓复合求积法，其指导思想就是先将积分区间分成若干