二维离散数据的数值微分问题的积分算子方法.pdfVIP

第32卷 第4期 武汉理工大学学报 ·信息与管理工程版 2010年8月 JOURNALOFWUT(INFORMATION＆MANAGEMENTENGINEERING) V0l_32No．4 Aug．2010 文章编号：1007—144X(2010)04—0591—04 文献标志码：A 二维离散数据 的数值微分 问题的积分算子方法 朱华平 ，吴传生，吕小红 (武汉理工大学 理学院，湖北 武汉430070) 摘 要：二维离散数据的数值微分问题是典型的不适定问题。基于GROETSCH的思想，提出了利用积分算子 方法来构造近似偏导数和方向导数的方法，并将该方法用于二维离散数据的数值微分问题，并给出了误差估 计。数值实验表明该方法是有效的。 关键词：二维数值微分 ；积分算子 ；不适定问题 中图分类号：0241．2；0241．4 DOI：10．3963／j．issn．1007—144x?2olo．40．020 目前大多数数值微分算法都是基于一维问题 给定步长 h0，对 D上的有界 司积函数 ， 的 ，关于二维数值微分问题的算法非常之少。 而实际上，很多实际问题如图像处理中的边界识 )，)定义积分算子(D (，y) J一 + 别、计算机层析成像(CT)、薄板的腐蚀，以及直管 tCOSOt，Y+tCOS )dt。在一定的条件下，可以用 道的内壁腐蚀等，都需要考虑二维甚至高维的情 (DJ )(，y)作为沿方向l的方向导数 ： 况。与一维数值微分问题一样 ，由有限的二维离散 扰动数据求偏导数及方向导数的问题是不适定的， cl 。 + O 。。 的近似，其中。。 、 y ‘ 因此需用特殊的方法来求解 。已有学者用 Tik— COS 为方向1的方向余弦。则当 =0，卢=订／2 honov正则化方法求解二维数值微分问题，并给出 正则化的构造算法以及相关的误差收敛估计 。 和 ：订／2， ：0时，其分别为偏导数 ， GROETSCH在文献 [4]中提出了用积分算子 ! 方法来构造一阶数值微分的方法，并给出了相应的 ay 。 误差估计。这个方法简单明了，计算简单 。 定理1 若厂(，Y)在D上所有的三阶偏导数 受其启发，笔者用积分算子方法来近似求函数的偏 连续，故所有的三阶偏导数在 D上有界，设它们 导数和方向导数，给出由有限的二维离散扰动数据 的绝对值在D上都不超过M。其给定的扰动数 来构造近似函数的方法以及近似求解函数的偏导 据 (，)，)满足 ilf二厂 I【* ( 。rf(，y)一 数和方向导数算法，并给出误差收敛估计。数值实 f (，Y)J≤ ，则： 验表明该算法是有效的。 y)一 COS of cos卢 1 求方向导数及偏导数的积分算子方法 记 DCR 是一单连通的有界凸区域， ，Y) ≤ 争、。+3，8b 在D上所有的三阶偏导数连续。假定f (，)，)是 证明 由 +tCOSOt，Y+tCOS卢)在 (，Y)的 D上 的可 积 函数，在 意 义 1I厂一厂 ll = Taylor展开式 ：