时间:2021-05-19
分析:
高中的时候,我们学习过,可以通过矩形法或者矩形法来求定积分。
思路就是将积分区间划分成n等份,然后将这n等份近似看成矩形(或梯形),然后对所有的矩形(或梯形)的面积进行求和。
简单的例子:
求函数X^2在的定积分
矩形法:
复制代码 代码如下:
#include<iostream>
#include<math.h>
using namespace std;
int main(){
float fun(float x);
float a,b;
cout<<"请输入函数X^2的定积分的下限a和上限b:";
cin>>a>>b;
int n=50;//将区间划分成50份
float h=(b-a)/n;//h是每个区间分大小
float s=0;//s是矩形的面积的和
float i=0;
for(i=a;i<b;i+=h){
s=s+fun(i)*h;
}
cout<<"\n结果是:"<<s<<endl;
cout<<endl;
}
float fun(float x){
return pow(x,2);
}
梯形法:
复制代码 代码如下:
#include<iostream>
#include<math.h>
using namespace std;
int main(){
float fun(float x);
float a,b;
cout<<"请输入函数X^2的定积分的下限a和上限b:";
cin>>a>>b;
int n=50;//将区间划分成50份
float h=(b-a)/n;//h是每个区间分大小
float s=0;//s是矩形的面积的和
float i=0;
for(i=a;i<b;i+=h){
s=s+((fun(i)+fun(i+h))*h)/2;
}
cout<<"\n结果是:"<<s<<endl;
cout<<endl;
}
float fun(float x){
return pow(x,2);
}
一个较复杂的例子
写一个通用函数,用来求sinx 、 cosx 、 e^x 、 x^2 的定积分
分析:fun为用来求定积分的通用函数,调用fun函数的时候,需要将积分的上限,下限,区间划分的份数以及被积函数的指针传递过来。
矩形法:
复制代码 代码如下:
#include<iostream>
#include<math.h>
using namespace std;
int main(){
float fsin( float x);
float fcos( float x);
float fe( float x);
float fpf(float x);
float fun(float a,float b, int n,float (*p)(float x));
float a[4],b[4],r[4];
cout<<"请输入求正弦函数定积分的上限a和下限b:";
cin>>a[0]>>b[0];
r[0]=fun(a[0],b[0],50,fsin);
cout<<"\n结果是:"<<r[0]<<endl;
cout<<"\n请输入求余弦函数定积分的上限a和下限b:";
cin>>a[1]>>b[1];
r[1]=fun(a[1],b[1],50,fcos);
cout<<"\n结果是:"<<r[1]<<endl;
cout<<"\n请输入求以e为底的指数函数定积分的上限a和下限b:";
cin>>a[2]>>b[2];
r[2]=fun(a[2],b[2],50,fe);
cout<<"\n结果是:"<<r[2]<<endl;
cout<<"\n请输入求X^2函数定积分的上限a和下限b:";
cin>>a[3]>>b[3];
r[3]=fun(a[3],b[3],50,fpf);
cout<<"\n结果是:"<<r[3]<<endl;
cout<<endl;
return 0;
}
float fsin(float x){
return sin(x);
}
float fcos(float x){
return cos(x);
}
float fe(float x){
return exp(x);
}
float fpf(float x){
return pow(x,2);
}
float fun(float a,float b,int n,float (*p)(float x)){
float i;
float h=(b-a)/n;
float s=0;
for(i=a;i<b;i+=h){
s=s+p(i)*h;//利用了矩形求面积的公式
}
return s;
}
梯形法:
复制代码 代码如下:
#include<iostream>
#include<math.h>
using namespace std;
int main(){
float fsin( float x);
float fcos( float x);
float fe( float x);
float fpf(float x);
float fun(float a,float b, int n,float (*p)(float x));
float a[4],b[4],r[4];
cout<<"请输入求正弦函数定积分的上限a和下限b:";
cin>>a[0]>>b[0];
r[0]=fun(a[0],b[0],50,fsin);
cout<<"\n结果是:"<<r[0]<<endl;
cout<<"\n请输入求余弦函数定积分的上限a和下限b:";
cin>>a[1]>>b[1];
r[1]=fun(a[1],b[1],50,fcos);
cout<<"\n结果是:"<<r[1]<<endl;
cout<<"\n请输入求以e为底的指数函数定积分的上限a和下限b:";
cin>>a[2]>>b[2];
r[2]=fun(a[2],b[2],50,fe);
cout<<"\n结果是:"<<r[2]<<endl;
cout<<"\n请输入求X^2函数定积分的上限a和下限b:";
cin>>a[3]>>b[3];
r[3]=fun(a[3],b[3],50,fpf);
cout<<"\n结果是:"<<r[3]<<endl;
cout<<endl;
return 0;
}
float fsin(float x){
return sin(x);
}
float fcos(float x){
return cos(x);
}
float fe(float x){
return exp(x);
}
float fpf(float x){
return pow(x,2);
}
float fun(float a,float b,int n,float (*p)(float x)){
float i;
float h=(b-a)/n;
float s=0;
for(i=a;i<b;i+=h){
s=s+((p(i)+p(i+h))*h)/2;//梯形法求面积
}
return s;
}
声明:本页内容来源网络,仅供用户参考;我单位不保证亦不表示资料全面及准确无误,也不保证亦不表示这些资料为最新信息,如因任何原因,本网内容或者用户因倚赖本网内容造成任何损失或损害,我单位将不会负任何法律责任。如涉及版权问题,请提交至online#300.cn邮箱联系删除。
本文为大家分享了C语言使用矩形法求定积分的通用函数,供大家参考,具体内容如下要求:写一个用矩形法求定积分的通用函数,分别求:sin(x),cos(x),e^x。
我就废话不多说了,直接上代码吧!#龙贝格法求积分importmatha=0#积分下限b=1#积分上限eps=10**-5#精度T=[]#复化梯形序列S=[]#S
本文实例为大家分享了C语言程序实现黎曼和求定积分,供大家参考,具体内容如下通过黎曼和解定积分既是把在xy平面中函数曲线与x轴区间区域划分成多个矩形并求它们的面积
题目要求利用指针编写程序,用梯形法计算下列公式中的定积分:参考代码首先说明一下指针的用处:因为所传递的参数均为数字,并不需要使用指针提高效率,故这里使用指针指向
求直方图中的最大矩形面积:例如给定直方图{2,3,1,2,4,2}则直方图中最大矩形面积为x=(3,6),|x|=3,y=2,max面积=6思考:利用枚举法/*