二重积分的复化梯形公式

1、摘要求函数在给定区间上的定积分，在微积分学中已给出了许多计算方法，但是，在实际问题计算中，往往仅给出函数在一些离散点的值，它的解析表达式没有明显的给出，或者，虽然给出解析表达式，但却很难求得其原函数。这时我们可以通过数值方法求出函数积分的近似值。前面已经对一重积分做了相应的学习和实验，但是正如我们所知的，在处理具体问题时，或在数学中有二重积分，更甚至n重积分，当遇到这类问题时，该如何处理？在微积分内容的学习中计算二重积分是用化为累次积分的方法进行的。计算二重积分也是计算累次数值积分的过程。为了简化计算，本实验仅对矩形域上的二重积分进行讨论。本实验将通过复化梯形公式来实现二重积分的求解，编程实现

2、二重积分的复化梯形公式。关键词：二重积分 复化梯形公式1、实验目的1) 通过本次实验体会并学习复化梯形公式在处理二重积分过程中的优点2) 通过对二重积分的复化梯形公式进行编程实现，提高自己的编程能力。3) 用实验报告的形式展现，提高自己在写论文方面的能力2、算法流程对矩形区域a,bc,d进行分割，沿x轴方向将a,bm等分，沿y轴方向将c,dn等分，步长分别为h=b-am, k=d-cn先用复化梯形公式计算cdf(x,y)dy,计算中将x当作常数，有cdf(x,y)dyk2fx,y0+fx,yn+2j=1n-1fx,yj然后将y当做常数，在x方向上计算积分，有abf(x,y0)dxh2fx0,y

3、0+fxm,y0+2i=1m-1fxi,y0abf(x,yn)dxh2fx0,yn+fxm,yn+2i=1m-1fxi,ynabj=1n-1fx,yjdx=j=1n-1abf(x,yj)dx h2j=1n-1fx0,yj+fxm,yj+2i=1m-1fxi,yj=h2j=1n-1fx0,yj+fxm,yj+4j=1n-1i=1m-1fxi,yj于是 abcdf(x,y)dxdyhk14fx0,y0+fxm,y0+fx0,yn+fxm,yn+12i=1m-1fxi,y0+i=1m-1fxi,yn+j=1n-1fx0,yj+j=1n-1fxm,yi+j=1n-1i=1m-1fxi,yj=hkj=0

4、ni=0mcijfxi,yj积分区间的4个角点的系数是14，4个边界的系数是12，内部节点的系数是1。3、算法实例用复化梯形公式计算二重积分 s=ab0201t2+R21-t2(cos2a2+sin2b2)dtd取R=30m,a=30.6m,b=29.6m,s的真实值为5679.81解：程序具体为：#include #include using namespace std;const int L(500);const double T(30),A(30.6),B(29.6);void main()float a,b,c,d;float h,k;int i,j,m,n;double Fx=0,F

5、1=0,F2=0,F3=0,F4=0,F5=0,F6=0;float xL=0;float yL=0;float RLL=0;cout请输入二重积分沿x轴方向的上下限ab;cout请输入沿x轴方向的等分数m;coutendl;cout请输入二重积分沿y轴方向的上下限cd;cout请输入沿y轴方向的等分数n;h=(b-a)/m;k=(d-c)/n;for(i=0;im+1;i+)xi=a+h*i;for(j=0;jn+1;j+)yj=c+k*j;for(i=0;im+1;i+)for(j=0;jn+1;j+)Rij=sqrt(xi*xi+T*T*(1-xi*xi)*(cos(yj)*cos(yj)/(A*A)+sin(yj)*sin(yj)/(B*B); /4个角点系数F1=R00+Rm0+R0n+Rmn;/4个边界系数for(i=1;im;i+)F2=F2+Ri0;F3=F3+Rin;for(j=1;jn;j+)F4=F4+R0j;F5=F5+Rmj;for(i=1;im;i+)for(j=1;jn;j+)F6=F6+Rij;Fx=A*B*h*k*(1/4*F1+1/2*(F2+F3+F4+F5)+F6);cout二重积分的结果为：Fxendl;运行结果：4、对结果进行分析通过用编程实现对上例的求解，可以看出结果较为准确，但是由于二重积分的原理问题，所以本身存在误差，而且