第二章 插值new.ppt

从Aitken插值公式向算法转化要考虑的问题是： (1) 插值公式右端n-1次多项式应如何处理； (2) 插值表中的元素应设置多少个存储单元； (3) 插值表中第k列第i行元素的计算公式。 ?  Aitken插值算法如下：  input x,(xi,yi),i=0,1,…,n 1 k L： for i=k,k+1,…,n do  end  ifk≠nthen k+1 k, go to L  ifk=n,output yn  　　Aitken插值算法为二重循环。外循环为k循环，用于计算Aitken插值表中的第k列；内循环为i循环，用于计算Aitken插值表中的第k列中的第i个元素。 Newton插值算法如下： inputx,(xi,yi),i=0,1,…,n。 y=y0,t=1。 forj=1,…,n do t=t*(x-xj-1)  for i=0,…,n-j do  end y=y+y0*t  end output (x,y),(xi,yi),i=0,1,…,n。 Newton插值算法中的j循环由三部分组成：计算(x-xj)的累积，存入t单元；内套一个i循环用来依次计算差商表中的各阶差商，存入yi单元；y单元用于存放Newton公式中各项累加之和。 (2) 差商具有对称性,即任意调换节点的次序, 差商的值不变。如 或者说： 该性质证明用数学归纳法 该性质由性质(1)得到 注意： 由性质(2)知k阶插商 又可以写成： 因为: 差商的计算方法(表格法): 规定函数值为零阶差商 差商表 Chashang.m 二、Newton插值公式 设插值多项式 满足插值条件 则待定系数为 称 定义3. 由插值多项式的唯一性,Newton插值多项式 与Lagrange插值多项式是相同的 为k次多项式 Newton插值法的优点是计算较简单,尤其是增加节点时, 计算只要增加一项,这是Lagrange插值无法比的. 但是Newton插值仍然没有改变Lagrange插值的插值曲线 在节点处有尖点，不光滑，插值多项式在节点处不可导 等缺点 ? 等距节点公式 /* Formulae with Equal Spacing */ 向前差分 /* forward difference */ i i i f f f - = ? + 1 i k i k i k i k f f f f 1 1 1 1 ) ( - + - - ? - ? = ? ? = ? 向后差分 /* backward difference */ 1 1 1 - - - ? - ? = ? i k i k i k f f f i?1 i i f f f - = ? 中心差分 /* centered difference */ 其中 当节点等距分布时: ? 差分的重要性质： ? 线性：例如 ? 若 f (x)是 m 次多项式，则 是 次多项式，而 ? 差分值可由函数值算出： ? = - + - = D n j j k n j k n f j n f 0 ) 1 ( ? = - + - - = ? n j n j k j n k n f j n f 0 ) 1 ( 其中 /* binomial coefficients */ ? 函数值可由差分值算出： k j n j k n f j n f D ? = + = 0 ? k k k h k f x x f ! ] , ... , [ 0 0 D = k n k k n n n h k f x x x f ! ] , ... , , [ 1 ? = - - k k k h f f 0 ) ( ) ( D = x 由 Rn 表达式 * 插值法是一种古老的数学方法，它来自生产实践.早在一千多年前，我国科学家在研究历法时就应用了线性插值与二次插值，但它的基本理论却是在微积分产生以后才逐步完善的，其应用也日益广泛.特别是由于计算机的使用和航空、造船、精密机械加工等实际问题的需要，使插值法在理论上和实践上得到进一步发展.尤其是近几十年发展起来的样条(Spline)插值，获得了极为广泛的应用，并成为计算机图形学的基础.　　本章主要讨论如何求插值多项式、分段插值函数、三次样条插值、插值多项式的存在唯一性及误差估计 实际问题中碰到的函数是各种各样的。有的表达很复杂，有的甚至给不出数学式子，而只是给出了一些离散数据——譬