[精华]复合辛普森公式[精华]复合辛普森公式
实验5 复合辛普森公式
李涛 201226100108 计自1201
一、实验目的
482-4, 用复合辛普森公式计算积分，使误差不超过(注意所给积分特点，1，cosxdx10,0
做出相应的处理后再计算) 二、实验步骤
1.算法原理
复合辛普森原理:
将区间划分为n等分，在每个子区间上采用辛普森公式，若记,,x,x[a,b]kk，1
1x,x，h,则得 k，12k2
n,1b, I,f(x)dx,f(x)dx,,ak,0
n,1h, ,[f(x)，4f(x)，f(x)]，R(...
![[精华]复合辛普森公式](https://swf.iask.com/kpfZ4ovPzcA.svg?ssig=rKn24R6Q77&Expires=1719909824&KID=sina,ishare&range=0-6347)
[精华]复合辛普森公式
实验5 复合辛普森公式
李涛 201226100108 计自1201
一、实验目的
482-4, 用复合辛普森公式计算积分，使误差不超过(注意所给积分特点，1，cosxdx10,0
做出相应的处理后再计算) 二、实验步骤
1.算法原理
复合辛普森原理:
将区间划分为n等分，在每个子区间上采用辛普森公式，若记,,x,x[a,b]kk，1
1x,x，h,则得 k，12k2
n,1b, I,f(x)dx,f(x)dx,,ak,0
n,1h, ,[f(x)，4f(x)，f(x)]，R(f).,kkkn，126k,0
记
n,1h, S,[f(x)，4f(x)，f(x)],nkk，12k6k,0
n,1n,1h, ,[f(a)，4f(x)，2f(x)，f(b)],,,k，12k6k,0k,1称为复合辛普森求积公式，其余项为
n,1hh4(4), R(f),I,S,,()f(,),,,(x,x).,nnkkkk，11802k,0
4于是当时，与复合梯形公式相似有 f(x),C[a,b]
b,ah4(4), R(f),I,S,,()f(,),,,(a,b)nn1802
4易知误差阶为，收敛性是显然的，实际上，只要则可得到收敛性，即hf(x),C[a,b]
b, limS,f(x)dxn,an,,
此外，由于中求积公系数均为正数，故知辛普森公式计算稳定。 Sn
2.算法步骤
复合辛普森:
首先将区间划分为n等分，在每个子区间上采用辛普森公式，若记,,x,x[a,b]kk，1
1x,x，h,则得 k，12k2
n,1n,1n,1hh S,[f(x)，4f(x)，f(x)],[f(a)，4f(x)，2f(x)，f(b)],,,nkk，12kk，12k66k,0k,0k,1
算法过程:
这里将辛普森公式写为Sn()函数，然后在Solve()函数里依次计算S1,S2,S4,S6.......当相
邻的精度小于eps时退出循环，则S2n保存结果。 三.程序代码
#include
#include
#define eps 1e-6
using namespace std;
double f(double x){
return sqrt(1+cos(x)*cos(x)); }//被积函数
double Sn(double a,double b,double n){
double h=(a+b)/(2*n);
double sum=0;
for(int k=1;k<=n;k++){
sum+=2*f(a+(2*k-1)*h);
sum+=f(a+2*k*h);
}
sum=(sum*2+f(a)-f(b))*h/3;
return sum;
}//辛普森公式
double Solve(double a,double b){
int i=1;
double S1n=Sn(a,b,i);
double S2n=Sn(a,b,2*i);
while(fabs(S1n-S2n)>eps){
cout<<"n = "<
本文档为【[精华]复合辛普森公式】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。
该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。
[版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
浙公网安备 33021202002483