分数阶光滑函数线性和二次插值公式余项估计.pdfVIP

2014年 11月 计 算 数 学 第 36卷第 4期 NOV．，2014 MATHEMATICA NUMERICA SINICA Vo1．36，No．4 分数阶光滑函数线性和二次插值公式 余项估计冰1) 王同科 佘海艳 刘志方 (天津师范大学数学科学学院，天津300387) 摘 要 本文在局部分数阶导数定义的基础上给出了高阶局部分数阶导数定义，并据此得到了一般形 式的分数阶Taylor公式．用该公式给出了分数阶光滑函数线性和二次插值公式余项的表达式，并 进一步导出了分段线性插值的收敛阶估计．针对分数阶导数临界阶计算困难的问题，本文利用线 性插值余项设计了一种外推算法，能够比较准确地求出函数在某点的局部分数阶导数的临界阶． 最后通过编写算法的Mathematica程序，验证了理论分析的正确性，并用实例说明了算法的有效 性． 关键词：局部分数阶导数；分数阶Taylor公式；线性和二次插值余项；临界阶估计 MR ~2ooo)主题分类：65D05，26A33 1．引 言 插值法是利用函数 f(x)在某区间中若干点的函数值，构造适当的特定函数，在这些点上 取已知值，在其他点上用这些特定函数的值作为f(x)的近似值的方法．设f(x)在某区间上有 定义， ，i=0，l，… ，rt为该区间上互不相同的 佗+1个点，多项式插值就是求一个不超过 rt 次的多项式p()，满足p(xi)=f(xi)，i=0，1，… ，n．通常有多种途径构造插值多项式，常用 的有Lagrange插值多项式和Newton插值多项式，其余项主要有三种形式，分别为微分型余 项、差商型余项和积分型余项．当f(x)在区间 [a，b]上n+1次可导 (记为f(x)∈C叶 [a，b]) 时，Lagrange插值多项式的微分型余项为 ，(n+1)r、 R()=，()一p()= 了 (—xo)(x-X1)… (— n)， ∈(n，6)· 当f(x)∈Via，b】时，Newton插值多项式的差商型余项为 R(x)=f(x)一p(x)=fix，Xo，Xl，… ，X】(—xo)(x—X1)·-·(X—X)． 当f(x)∈C [a，b](1 n+1)时，插值的积分型余项为 【】 1 n R()= )一p()= ∑ (t)k-1 (t)dt， 其中li(x)为对应节点Xi的Lagrange插值基函数． 插值余项的不同表示有着不同的应用，当被插值函数f(x)仅连续时，只能用差商型余项， 当f(x)∈C a，b]时用积分型余项，当f(x)∈Cn+a【，b]时还可以使用微分型余项． 2013年 1o月 10日收到． )基金项 目：国家 自然科学基金 资助项 目 394 计 算 数 学 对于更复杂的插值问题，许多学者从不同侧面做了深入的研究．文献 f21对于低度光滑函 数给出了插值余项及其导数的误差估计，文献 [3]总结了一些插值的收敛性问题，文献 [4]研 究了连续函数基于非等距节点线性插值的误差估计问题，文献 [5】给出了分段光滑函数的误 差估计式． 从插值的实际计算效果来看，我们发现，对于一类非充分光滑的函数，最简单的情形如 f(x)=X。，X∈[0，1】，0 1，其分段线性插值函数在 =0附近误差较大．对于这类在插 值区间上包含不可导点的函数，虽然可以使用差商型余项表示式，但无法给出插值余项的精细 刻画形式，如何准确地刻画其插值余项需要进一步研究．目前，对于这类函数，通常使用分数 阶导数描述其光滑程度． 分数阶导数的研究已有长期的历史，其定义也有多种，主要有 Riem