二重积分的复化梯形公式

《二重积分的复化梯形公式》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

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

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

3、k2fx,y0+fx,yn+2j=1n-1fx,yj然后将y当做常数，在x方向上计算积分，有abf(x,y0)dx≈h2fx0,y0+fxm,y0+2i=1m-1fxi,y0abf(x,yn)dx≈h2fx0,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,yj5于是abcdf(x,y)dxdy≈hk14fx0,y0

4、+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=0ni=0mcijfxi,yj积分区间的4个角点的系数是14，4个边界的系数是12，内部节点的系数是1。3、算法实例用复化梯形公式计算二重积分s=ab02π01t2+R21-t2(cos2θa2+sin2θb2)dtdθ取R=30m,a=30.6m,b=29.6m,s的真实值为5679.81解：程序具体为：#include

5、eam>#includeusingnamespacestd;constintL(500);constdoubleT(30),A(30.6),B(29.6);voidmain(){floata,b,c,d;floath,k;inti,j,m,n;doubleFx=0,F1=0,F2=0,F3=0,F4=0,F5=0,F6=0;floatx[L]={0};floaty[L]={0};floatR[L][L]={0};5cout<<"请输入二重积分沿x轴方向的上下限"<>a>>b;cout<<"请输

6、入沿x轴方向的等分数"<>m;cout<>c>>d;cout<<"请输入沿y轴方向的等分数"<>n;h=(b-a)/m;k=(d-c)/n;for(i=0;i

7、[i])*(cos(y[j])*cos(y[j])/(A*A)+sin(y[j])*sin(y[j])/(B*B)));}}//4个角点系数F1=R[0][0]+R[m][0]+R[0][n]+R[m][n];//4个边界系数for(i=1;i

8、(1/4*F1+1/2*(F2+F3+F4+F5)+F6);cout<<"二重积分的结果为："<