插值型数值微分的性质及其应用.pdfVIP

1．引言 求导是工程中常遇到的计算，在函数未知或函数表达 30 ， … ， 式复杂时，需借助数值方法，即应用数值微分进行计算。插 = … 嵩 l 值型数值微分是当前数值微分的主要方法之一 ，该方法基 LBdx)J 于插值多项式建立用于求道的数值微分公式，实践中一般 =Ⅵ，B(x) 对插值做等距节点处理。插值型微分有何特点?插值型微分 特别地，取 x=xj，有求低阶数值导数的微分公式 (1)和 公式能否有简单快捷的推导方法? 目前研究主要提出了基 (2)。 于数值精度的概念和基于插值多项式的Lagrange形式推导 fAo(xj)] 出一阶插值型数值微分公式。这两种推导方法都比较复杂， (x) 【wy ．o(，…，嵩 ‘lWA() (1) 且不能反映插值型微分 自身的特性。本文首先给出插值型 【 x【j)j 数值微分的矩阵形式，然后通过分析w(x)的特性，推导出等 B『o(xj)1 距节点条件下插值型数值微分的几个性质，最后说明如何 x) 【、v-yo(，…，嵩 ‘lWB() (2) 基于这些性质快速得到低阶 (一阶和二阶)插值型微分公 B【n)【j)【J 式 。 2．插值型数值微分概述 计算误差分别为 )(j￡)= 导w’()，()= 对区fnq[a，b】上的被插函数 fx【)，取n+1个插值节点a=x0 、v． …x．--b，则被插函数可表示为 f(x)=Ldx)+R(x)。如果记w(x)= 引理：如果Xl—x0==…=xn—xn—l_h，则有 (x-Xo)…(x—x)Ⅱ，则余项R(x)= 导w(x)(∈∈a【，b】)，插值 (1)w’()=(一1)“1(n—j)!lln， (xD) (一1)Ⅱ、v’()。 多项式 x)的拉格朗日形式为 x) [yk ]。 (2)’()2)(j)L【u【=1U一v蓦=1VJJh-1,wx(t,)(一1)x(j)。 (3)w…()= 由于 f (x)=Ln(x)+R (x)，插值型数值微分用 f呻(x) LⅡD【)(x)来计算导数的近似值，相应误差为R (x)。工程中常用 0『