复化梯形

将积分区间[a,b]划分n等分，步长，求积节点，在每个小区间上应用梯形公式

然后将它们累加求和，作为所求积分I的近似值.

记

式为复化梯形求积公式，下标n表示将区间n等分。

算法流程

算法代码

[cpp] view plaincopy

double f(double x){
if(x==0) return 1;
else return (sin(x)/x);
}
double FuhuaTixing(int n,double a,double b){
double h = (b-a)/n;
double x = a;
double s = 0;
for(int k=0; k< n-1; k++){
x += h;
s += f(x);
}
double T = (f(a)+s*2+f(b))*h/2;
return T;
}
int main(){
char ans='n';
do{
cout<<"请输入积分区间(a,b)："<<endl;
double a;
double b;
cin>>a>>b;
cout<<"请输入等分份数n： "<<endl;
int n; cin>>n;
cout<<"由复化梯形公式球的结果："<<FuhuaTixing(n,a,b)<<endl;
cout<<"是否要继续？(y/n)";
cin>>ans;
}while(ans == 'y');
return 0;
}

复化辛复生

将积分区间[a,b]划分n等分，记子区间的中点为在每个小区间上应用辛普森公式，则有

其中

记

式为复化辛普森求积公式。

算法流程

算法代码

[cpp] view plaincopy

double f(double x){
if (x==0)
return 1;
else return (sin(x)/x);
}
double Xinfusheng(double a,double b,int n){
double h = (b-a)/n;
double x = a+1/2*h;
double s = 4*f(x);
for(int k=1;k<n;k++){
x += 1/2*h;
s += 4*f(x);
x += 1/2*h;
s += 2*f(x);
}
double T=(f(a)+s+f(b))*h/6;
return T;
}
int main(){
char ans='n';
do{
cout<<"请输入积分区间(a,b)："<<endl;
double a;
double b;
cin>>a>>b;
cout<<"请输入等分份数n： "<<endl;
int n;
cin>>n;
cout<<"由复化梯形公式球的结果："<<Xinfusheng(a,b,n)<<endl;
cout<<"是否要继续？(y/n)";
cin>>ans;
}while(ans == 'y');
return 0;
}

实验过程原始记录

分别用复化梯形公式和复化辛浦生公式计算定积分

取n=2，4，8，16，精确解为0.9460831

实验结果及分析

1、用复化梯形公式和复化辛甫生公式都能得到较为准确的结果，且等分份数越多，结果的精度越高，梯形公式虽然在等分16份时得到精度与等分4份时相同，但已经越来越接近精确解。辛甫生公式由于C++运算中双精度数值（double）只有7位有效数字的限制，增加等分份数时不容易看出其精度的增加。
2、比较两种方法运算的结果，复化辛甫生公式等分2份时实际要计算5个点的函数值，与复化梯形公式等分4份时计算量基本相同，但得到精度明显复化辛甫生公式要精确很多。

3、复化梯形公式和复化辛甫生公式对于光滑性较差的被积函数都能得到较为精确的结果，而且公式简单，十分利于编译简单的程序由计算机运算，因而得到广泛的应用。
4、实验中的主要误差来自于计算机浮点运算中的截余。

（转载请注明作者和出处：http://blog.csdn.net/xiaowei_cqu 未经允许请勿用于商业用途）

微分求积：复化梯形、复化辛浦生相关推荐

【数值分析】微分求积：复化梯形、复化辛浦生
本科课程参见:<软件学院那些课> 复化梯形将积分区间[a,b]划分n等分,步长,求积节点,在每个小区间上应用梯形公式然后将它们累加求和,作为所求积分I的近似值. 记
复化辛甫生求积法matlab,数值分析习题.doc
数值分析习题第一章绪论姓名学号班级习题主要考察点:有效数字的计算.计算方法的比较选择.误差和误差限的计算. 1 若误差限为,那么近似数0.003400有几位有效数字?(有效数字的计算) 2 ...
python用复化梯形公式积分_复化梯形求积分——用Python进行数值计算
用程序来求积分的方法有很多,这篇文章主要是有关牛顿-科特斯公式. 学过插值算法的同学最容易想到的就是用插值函数代替被积分函数来求积分,但实际上在大部分场景下这是行不通的. 插值函数一般是一个不超过n次 ...
计算方法之数值积分方法——复化梯形法，复化辛普森法，龙贝格法，三点高斯公式附matlap程序下载
数值积分复化求积法就是将求积区间[a,b]划分为n等份,步长h=(b-a)/n,等分点为xi=a+ih,i=0,1,2,-,n.然后用低阶求积公式求的每个字段[xi,xi+1]上的积分值I,然后再将 ...
c语言simpson积分计算方法,数值分析复化Simpson积分公式和复化梯形积分公式计算积分的通用程序...
数值分析复化Simpson积分公式和复化梯形积分公式计算积分的通用程序数值分析第五次程序作业 PB09001057 孙琪 [问题] 分别编写用复化Simpson积分公式和复化梯形积分公式计算积分的通 ...
乘幂法、复化梯形以及二分法求特征值【相关例题python实现】
目录写在前面的话乘幂法例题代码结果二分法求对称三对角矩阵特征值代码结果复化梯形例题代码结果写在前面的话考试结束,<数值计算方法>编程作业发出来给大家看看~ 乘幂 ...
复化梯形c语言计算定积分,第六章函数和宏定义实验(2)
1.利用复化梯形公式计算定积分 2.计算Ackerman函数 3.编写计算x的y次幂的递归函数getpower(int x,int y),并在主程序中实现输入输出 4.编写计算学生年龄的递归函数 5. ...
用c语言写复化梯形积分公式,第六章函数和宏定义实验（2）
C程序设计实验报告姓名:童飞实验地点:教学楼514教室实验时间:5月16日实验项目: 1.利用复化梯形公式计算定积分 2.计算Ackerman函数 3.编写计算x的y次幂的递归函数 ...
[计算机数值分析]复化梯形的递推公式的变步长算法求积分
如下图变步长的算法框图,其中T1和T2分别代表二分前后的积分值. 问题描述: 运行示例: 源码: #include<iostream> #include<cmath> doub ...

微分求积：复化梯形、复化辛浦生

复化梯形

算法流程

算法代码

复化辛复生

算法流程

算法代码

实验过程原始记录

实验结果及分析

（转载请注明作者和出处：http://blog.csdn.net/xiaowei_cqu 未经允许请勿用于商业用途）

微分求积：复化梯形、复化辛浦生相关推荐

最新文章

热门文章