2.插值法-1程序.ppt

第2章 插值法 离散变量连续化 许多实际问题都用函数 y=f (x) 表示某种内在规律的数量关系, 数学上的一条曲线。 部分函数是通过实验或观测得到的，只能给出函数在定义区间上部分离散点的数值。 有的函数表达式计算复杂，使用不方便，需要根据已有的函数表构造一个能反映原函数特性，又便于计算的函数。通常选择一类简单函数如多项式函数、分段函数等 插值函数 拉格朗日插值法 拉格朗日插值多项式 插值余项与误差估计 均差与牛顿插值公式 牛顿插值公式 差分与等距节点插值 等距节点插值公式 分段低次插值 并非插值多项式次数越高，精度也越高 一般采用分段低次插值 如：一次(线性)插值 二次(抛物线)插值 三次插值 插值结点的选择 已知结点及其函数值(xi,yi),i=0,1,2,…,n，求f(x) x应在插值结点之间，如: xjxxk 插值区间应尽量小，如： xjxxi+1 分段低次插值 2.6.2 分段线性插值 分段二次插值 由相邻三个插值节点构造二次插值多项式 构成分段二次插值 表2-5列出了 时的 的计算结果及 在 上的误差 可见，随 的增加， 的绝对值几乎成倍增加. 这说明当 时 在 上是不收敛的. Runge证明了，存在一个常数 ，使得当 时， 而当 时 发散. 均差有如下的基本性质： （3.3） 这个性质可用归纳法证明. 1° k 阶均差可表为函数值 的线 性组合， 这性质也表明均差与节点的排列次序无关，称为均差 的对称性. 即 均差的性质 3° 若 在 上存在 n 阶导数，且节点 （3.5） 这公式可直接用微分中值定理证明. （3.4） 2° 由性质1°及（3.2）可得 则 阶均差与导数关系如下： 均差性质（续） 均差计算可列均差表如下（表2-1）. 均差表 根据均差定义，把 看成 上一点， 可得 只要把后一式代入前一式，就得到 其中 （3.6） 返回 （3.7） 是由（2.10）定义的. 显然，由（3.6）确定的多项式 满足插值条件， 且次数不超过 ， 称 为牛顿（Newton）均差插值多项式. 系数 就是均差表2-1中加横线的各阶均差，它比拉格朗日插值计算量省，且便于程序设计. 其系数为 它就是形如（3.1）的多项式， 但（3.7）更有一般性，它在 是由离散点给出的情形或 导数不存在时也是适用的. （3.7）为插值余项，由插值多项式唯一性知，它与 拉格朗日插值多项式的余项应该是等价的. 事实上，利用均差与导数关系式就可以证明这一点. 牛顿插值多项式的优点还在于它的递进性，当增加 插值节点时，只要在原来插值多项式的基础上增加一项 即可. 首先根据给定函数表造出均差表. 给出 的函数表（见表2-2），求4次牛顿插 值多项式，并由此计算 的近似值. 例2 从均差表看到4阶均差近似常数，5阶均差近似为0. 故取4次插值多项式 做近似即可. 于是 按牛顿插值公式，将数据代入 截断误差 这说明截断误差很小，可忽略不计. 实际应用时经常遇到等距节点的情形，这时插值公式 可以进一步简化，计算也简单得多. 差分及其性质 设函数 在等距节点 上 的值 为已知，这里 为常数，称为步长. 为了得到等距节点的插值公式，先介绍差分的概念. 记号 （4.1） （4.2） 分别称为 在 处以 为步长的向前差分，向后差分 符号 ， 分别称为向前差分算子，向后差分算子 定义3 利用一阶差分可定义二阶差分为 一般地可定义 阶差分为 引入两个常用算子符号： 不变算子 及移位算子 于是，由 定义如下： 可得： 不变算子 移位算子 向前差分 向后差分 差分基本性质. 性质1 其中 为二项式展开系数. 例如 各阶差分均可用函数值表示. （4.4） （4.5） 性质2 例如，可用向前差分表示 ， 所以 （4.6） 可用各阶差