二重积分的数值解.pdf

二重积分的数值解 摘 要 在定积分数值计算公式的基础上，推广得到矩形区域上二重积分的中点公 式、梯形公式和Simpson公式，并对其误差进行了估计。为了减小误差，将上述 三个公式运用了复化求积的思想，得到复化中点公式、复化梯形公式和复化 Simpson 公式，并对其误差进行了估计，由此得到二重积分的自适应复化求积法。 最后对复化梯形公式进行了Richardson外推。 关键字 二重积分；自适应复化求积法；外推算法；误差分析 1 引言 积分的计算是自然科学的中的一个基本问题。如果被积函数存在原函数，则 可由牛顿-莱布尼茨公式来实现。但是，在许多实际问题中，被积函数找不到用 初等函数表示的原函数，或者被积函数仅由实验测量或数值计算给出的若干离散 [1-2] 点上的量值，于是引入了数值积分。一元的数值积分的近似计算已解决 ,二重 [3] [4] [5] 积分的中点公式 、梯形公式 和Simpson公式 也有简单的介绍。 本文研究了二重积分的求积公式及其复化公式和外推算法。本文对三重积分 也是适用的,最后给出算例验证理论结果。 2 二重积分的求积公式 考虑矩形区域D: a≤ x≤b,c≤ y≤ d 上的连续函数f(x,y)的二重积分 { } d b I= ∫∫f(x,y)dxdy= ∫ ∫f(x,y)dxdy. (2.1) D c a 2.1 中点公式及其截断误差 [1] 将(2.1)式应用中矩形公式 两次,可得： d b I= ∫∫ f(x,y)dxdy c a (2.2) a+b c+d ≈ (d−c)(b−a)f( , ). 2 2 称上式为二重积分的中点公式。 下面给出误差的分析，如果f(x,y)在区域D内具有两阶连续偏导数，已知定 [2] 积分的中点公式的截断误差 为： -1 - b a+b (b−a)3 ′′ ∫f(x)dx−(b−a)f( 2 )= 24 f (ξ),ξ∈[a,b]. a 故可以得到: b 3 2 a+b (b−a) ∂ f(ξ,y) f(x,y)dx−(b−a)f( ,y)= 1 , ξ ∈[a,b], (2.3) ∫ 2 24 ∂x2 1 a d a+b