二维带边界偏导数值的复化数值积分公式.pdfVIP

第 3O卷第 3期 大 学 数 学 Vo1．30。№ ．3 2014年 6月 CoLLEGEM ATHEM ATICS Jun．2014 二维带边界偏导数值的复化数值积分公式 董 灏， 聂玉峰 (西北工业大学 理学院，陕西 西安 710129) [摘 要]构建了一类二维带边界偏导数值的复化数值积分公式，给出了所建立的两种数值积分公式的 稳定性分析、误差分析和代数精度．与二维复化四点高斯数值积分公式相对 比，所建立的带边界偏导数值的复 化梯形、复化辛普森求积公式在达到相 同精度时所需积分节点大大减少 ，积分的时间复杂度也随之大大减少 ， 实例验证结果 良好 ． [关键词]边界偏导数值；复化梯形；复化辛普森；积分精度；积分节点 [中图分类号]O2414 [文献标识码]A [文章编号]1672—1454(2014)03—0031—07 1 引 言 众所周知 ，在科学工程计算中需要进行大量的定 (重)积分的计算 ，但常常很难得到积分 的精确值 ， 只能求其近似值，因而研究和构造高精度有效的数值积分公式非常有意义．本文通过在原有的复化梯 形、复化辛普森公式的基础上通过添加被积 函数未被利用的偏导数值 ，充分利用被积 函数的 自身信息， 并利用同一内结点上相邻区间公式中偏导数项的系数互相抵消这一特点，得到 了比复化高斯数值积分 公式精度更高，所需积分节点大大减少的数值积分公式． 2 二维情形下公式的建立和分析 在文献[2]中已建立了一重积分 If(x)dx的带端点导数值的Hermite数值积分公式，具体公式 如下 ： If(x)dx—I-f(--1)+厂(1)]+寺 ·If(-1)一f(1)]+ ·f“()， ∈[一1，1]．(1) J — l u -V 1 厂()d 一 1 一 15·[7厂(一1)+16f(0)+7，(1)]+击·If(一1)一 (1)]+ ·，’(， 刁∈[一1，1]． (2) 利用积分变换和累次积分就可以建立二维情形下的带边界偏导数值的复化求积公式 ，下面给出公式的 建立过程． 2．1 二维情形下数值积分公式的建立 考虑到二维情况下二：重积分的积分域相比一维情况复杂得多．文章 以矩形域为研究对象 ，因为对于 复杂的积分域也可 以通过积分区间变换或者区间划分来产生矩形积分域，如可将 圆形积分域通过参数 变换 (，)一 (r，)变为矩形积分域；曲线边界积分域可通过矩形网格逼近，剩余的非矩形区域用其它 积分公式如三角形上的数值积分公式计算． [收稿 日期]2o13—09—16 [基金项 目]陕西高等教育教学改革研究项 目(13BY12) 32 大 学 数 学 第 30卷 二维 情彤 卜矩 彤域 上 的带边 界偏 导数 值 的复化 梯形 、复化 平晋 森 数值 积 分 公 式 只需要 在 二 重积 分 累次积分时使用两次一维情形下的带端点导数值的Hermite数值积分公式(1)，(2)就可以得到 ，即 -r：Ir：，，dd—j：(-r：厂z，)．，d)d． 下面给出积分区间为[n，6]×[c，]上的二重积分 ．r，(z，)dd的带边界偏导数值的复化梯 形 、复化辛普森数值积分公式． 2．1．1 二维带边界偏导数值的复化梯形公式 记四点点集 G一 {(一1，一 1)，(一1，1)，(1，一1)，(1，1)}，记点 P为集合G 内的点 ，定义算符 ， 8 ，8， 一 _{1’； 萎i，， 一 ； 嚣 ，’，