插值法培训教程文件.pptVIP

一般地 (5-29) (5-30) 插值法 其中 为二项式展开式的系数. 32 (3)各种差分之间可以互化,这由差分表即可看出. 向后差分与中心差分化成向前差分的公式如下: (5-31) (5-32) 例如, (4)可用差分表示差商.例如 插值法 33 一般地有 (5-33) 插值法 34 结合式(5-21),可得差分与导数的关系 (5-34) 5.2.4 等距节点插值公式 如果插值节点是等距的,则插值公式可用差分表示. 设已知节点 将式(5-33) 代入插值公式(5-19), 得 35 p41 二阶(向前)差分为: 一般地, m阶向前差分为 (5-23) 插值法 37 一阶向前差分为: 设有等距节点 步长h 为常数, (5-24) 称为 在 处以 为步长的向后差分. 插值法 称为 在 处以h 为步长的二阶向后差分. 一般地, m阶向后差分为 (5-26) 38 p31 差分具有以下性质: (1)各阶差分均具有线性性,即若 则对任意正整数m,都有 (2)各阶差分均可表示成函数值的线性组合, 插值法 39 (5-29) (5-30) 其中 (3)各种差分之间可以互化,这由差分表即可看出. (4)可用差分表示差商. 插值法 40 (5-33) 结合式(5-21),可得差分与导数的关系 (5-34) 5.2.4 等距节点插值公式 如果插值节点是等距的,则插值公式可用差分表示. 设已知节点 将式(5-33) 代入插值公式(5-19), 得 (5-19) 41 (5-35) 若令 则上式又可变形为 (5-36) 式(5-36)称为Newton向前插值公式,其余项为 42 p45 (5-37) 完全类似地,也可以用向后差分表示Newton插值公式.令 则有 (5-38) 插值法 式(5-38)称为Newton向后插值公式,其余项为 43 p45 (5-39) Newton向前,向后插值公式均是Newton插值公式在 等距节点时的变形. 实际计算时,也可列表进行. 将表5-6 中对角线上的差分值与对应右端因子乘积求和即得 Newton向前插值公式,而Newton向后插值公式则为最后 的节点所在行的各阶差分值与对应列下端因子乘积之 和. 插值法 44 插值法 45 p42 p43 53 例4 已知 的函数表如下, 分别用Newton向前, 向后插值公式,求 的近似值. [解] 取 按表5-6计算得 插值法 46 插值法 47 p50 p51 Newton向前插值公式为 故 代入上式得 插值法 又 48 由式(5-37),误差为 插值法 49 Newton向后插值公式为 将 代入上式得 插值法 50 p47 查表可得 如果取 用二阶Newton 向后插值公式,则得 将 代入上式得 插值法 51 p47 其误差为 为便于递推计算,当插值点x接近数据表头时, 用向前插值公式,当插值点x接近数据表尾时,则采用向 一般 后插值公式. 插值法 52 例 已知等距节点及相应点上的函数值如下表 0 1 2 3 0.4 0.6 0.8 1.0 1.5 1.8 2.2 2.8 分别Newton向前向后插值公式计算 及 的近似值. [解] 取 先构造向前差分表 53 插值法 54 Newton向前插值公式为 p56 故 代入上式得 插值法 又 55 由式(5-37),误差为 插值法 56 Newton向后插值公式为 p54 将 代入上式得 插值法 57 其误差为 插值法 58 若 则 * §5.2 牛顿(Newton)插值 如果将图5-1中的直线用点斜式方程表示,即把线性 插值公式改写成以下形式 (5-13) 由此导出插值多项式的又一种形式--牛顿插值公式. 在导出Newton公式前,先介绍公式表示中所需要用到的 差商,差分的概念及性质. 插值法 1 5.2.1 差商 定义5.1 设有函数 为一系列互不 相等的点, 称为 关于点 的一阶差商(也称均差),记为 即 类似于高阶导数的定义,称一阶差商的差商 插值法 (或 ) 2 为 关于点 的二阶差商,记为 一般地,称 为 关于点 的 k阶差商, (5-14) 插值法 (或 ) (或 ) 时称 当 关于点 的零阶差商 记为 为 3 (1) 各阶差商均有线性性,即若 则对任意正整数 k 都有 (5-15) (2) k 阶商差 可表示成 的线性组合 (5-16) 其中 插值法 容易证明,差商具有下述性质: 4 事实上, 是显然的.当 时 一般情况可用归纳法证明. (3) 各阶差商均具有对称性,即改变节点的位置, 差商值不变.如 插值法 5 (4)若 是n次多项式,则一阶差商 也是 是n次多项式,则 于是 可分解为 其中 为 n-1次多项式.所以 为 次多项式. 计算各阶差商,可按表5-1