实验四 复合辛普森公式.docxVIP

实验四 复合辛普森公式 一、实验目的 1、理解复合辛普森公式的方法原理； 2、利用复合辛普森公式计算积分。 二、实验步骤 1.算法原理 复合辛普森原理： 将区间 [a, b] 划分为 n 等分，在每个子区间 ?x , x ?上采用辛普森公式，若记 k k ?1 1 x k ?1 2 ? x ? k h, 则得 2 I ? ?b a  f (x)dx ? ?n?1 k ?0  f (x)dx ? h ?n?1 [ f (x 6 k k ?0 记 ) ? 4 f (x  k ?1 2 ) ? f (x k )] ? R n  ( f ). S ? h ?n?1 [ f (x n 6 k k ?0  ) ? 4 f (x  k ?1 2  ) ? f (x )] k ? h [ f (a) ? 4?n?1 6  f (x  k ?1 2 ) ? 2?n?1  f (x k ) ? f (b)], 称为复合辛普森求积公式，其余项为 k ?0 k ?1 R ( f ) ? I ? S n n ? ? h h 4 ?n?1 ( )180 2 ( ) k ?0 f (4) (? k ),? k  ? (x k  , x ). k ?1 于是当 f (x) ? C 4 [a, b] 时，与复合梯形公式相似有 b ? a h R ( f ) ? I ? S n n ? ? ( )4 180 2 f (4) (?),? ? (a, b) 易知误差阶为h 4，收敛性是显然的，实际上，只要 f (x) ? C[a, b] 则可得到收敛性，即 lim S n n?? ? ? b a f (x)dx 此外，由于S n 中求积公系数均为正数，故知辛普森公式计算稳定。2.算法步骤 复合辛普森： 首先将区间 [a, b] 划分为 n 等分，在每个子区间 ?x , x ?上采用辛普森公式，若记 k k ?1 1 x k ?1 2 ? x ? k h, 则得 2 S ? h ?n?1 [ f (x n 6 k ) ? 4 f (x  k ?1 2 ) ? f (x k )] ? h [ f (a) ? 4?n?1 6  f (x  k ?1 2 ) ? 2?n?1  f (x k ) ? f (b)] k ?0 3.程序流程图 k ?0 k ?1 开始 开始 输入数据a,b,n k ? 0, h ? b ? a n ， x ? a ? kh k T ? h ?n?1 [ f (x ) ? f (x )] ? h [ f (a) ? 2?n?1 f (x ) ? f (b)] n 2 k k ?1 2 k k ?0 k ?1 S ? h [ f (a) ? 4?n?1 f (x n 6 k ?1 2 k ?0 ) ? 2?n?1 f (x ) ? f (b)], k k ?1 R ? I ?T 1 n ， R ? I ? S 2 n k ? k ? 1 输出 R , R 1 2 结束 三.程序代码 #includestdio.h #includemath.h #includeconio.h #define eps 0.000001 //计算精度 float func(double x) //% func-- 求积函数（函数文件） { while(x==0) return 1; double y; y=(sin(x)/x); return(y); } double bsimpson(double a,double b) //% a, b -- 求积区间 { int i,n; double h,p,e,s; double t1,t2,s1,s2,x; n=1; h=b-a; t1=h*(func(a)+func(b))/2.0; s1=t1; e=eps+1.0; while(e=eps) {  /*用 t1 代替 s1*/ s=0.0; for(i=0;i=n-1;i++) { x=a+(i+0.5)*h; s=s+func(x); } t2=(t1+h*s)/2.0; s2=(4*t2-t1)/3.0; e=fabs(s2-s1); t1=t2; s1=s2; n=n+n; h=h/2.0; } return(s2); } main() { double a,b,s; printf(input a,b=); scanf(%lf,%lf,a,b); s=bsimpson(a,b); printf(s=%lf\n,s); getch(); return 0; } 四、心得体会 通过本次试验，我理解复合辛普森公式的方法原理，并且能利用复合辛普森公式计算积分，在数值