用积分算子方法求二阶数值微分.pdfVIP

第30卷 第6期 武 汉 理 工 大 学 学 报 Vo1．30 No．6 2008年6月 JOURNAL OF WLrHAN UNIVERSITY OF TECHN0IDGY Jun．2008 用积分算子方法求二阶数值微分 吕小红，吴传生，高 飞 (武汉理工大学理学院，武汉430070) 摘 要：数值微分就是用离散方法近似地求出函数在某点的导数值。关于数值微分己有许多求解方法，但这些方法都 有各自的局限性，且关于二阶导数近似逼近的方法研究相对较少。基于Groetsch的思想，提出了利用积分算子来构造近 似二阶导数的方法，并将此方法应用于二阶数值微分问题，给出了相应的误差估计。通过数值实验表明，此方法对于二 阶数值微分问题十分有效，而且计算量小。 关键词： 二阶数值微分； 积分算子； 不适定问题； 正则化 中图分类号： O 241．2 文献标识码： A 文章编号：1671—4431(2008)06—0178—03 Solve Second·order Numerical Differentiation by Using Integral Operator LVXiao—hong，WU Chuan—sheng，GAO Fei (School of Sciences，Wuhan University of Technology，Wuhan 430070，China) Abstract： Numerical differentiation was that derivative value of a function at a certain point was approximately solved in dis— crete method，There have been a lot of．solutions to numerical differentiation．However，they have their limitations of their OWII．Moreover，there were relatively few researches on derivative of second order approximation．Based on Groetsch’S thought a method WaS brought forward that could structure approximate second—-order derivative by using integral operator and by apply—- ing it tO second—order numerical problems we provide corresponding error estimate． It WaS made clear by numerical experiments that the method WaS much valid to second—order numerical differential problems；what’S more there were not many calculations． Key words： second—order numerical differentiation； integral operator； ill—posed problem； regularization 求解数值微分问题，就是通过测量函数在一些离散点上的值，求得函数的近似导数。这个问题，产生于 许多数学模型中以及实际问题中。找到一个稳定的求解微分问题的方法，无论是对科学研究，还是实际应 用，都是非常重要的。众所周知，当函数值本身有较小的扰动时，相应的导数的扰动可以是任意大的，或者 说，由函数在有限个点的值求导数的问题是不适定的。对此问题，已有大量的数值算法和收敛性分析，其中 主要方法是利用有限差分并选择合适的步长来近似导数，或者构造一类积分算子来近似导数。如果离散点