计算方法 插值方法-2-课堂教学版.pptVIP

* * * * * * 函数f(x)在点x0,x1,…,xn处的k阶差商与点x0,x1,…,xn的排列顺序无关。 * 把x看作是[a,b]上一点，对差商稍做变形可得 依次把后面的式子带入到前面的式子中，可得 相对导数形式的插值余项来讲，差商形式的余项更有一般性，它对f是由离散点给出的情形或f的导数不存在时均适用。 * 函数f(x)在点x0,x1,…,xn处的k阶差商与点x0,x1,…,xn的排列顺序无关。既有： * 在这种意义下，拉格朗日插值可理解为泰勒插值的离散化形式 * 前面讨论了节点任意分布的插值公式，但实际应用时经常遇到等距节点的情形，这时插值公式可以进一步简化，计算也简单得多，为了得到等距节点的插值公式，我们先介绍差分的概念。 分别称为f(x)在xk处以h为步长的向前差分，向后差分及中心差分。 * * 计算差分可列差分表，下表是向前差分表 * 也有的参考书在列查分表的时候使用这种形式，这个时侯选择差分值之需要用到第一行的即可 * 在节点等距的情况下，差商可以用差分来表示。对向前差分： * 在节点等距的情况下，差商可以用差分来表示。对向前差分： * 在节点等距的情况下，差商可以用差分来表示。对向前差分： * 给定f(xi) = yi 的函数表 * 等距节点插值公式有不少实际应用。例如，很多工程设计计算都需要查各种函数表，用计算机计算时就必须解决计算机查表问题。如果把整个函数表存入内存，往往占用存储单元太多；如果用一个解析表达式近似该函数，则可能达不到精度要求。因此，采用存放大间隔函数表，并用插值公式计算函数近似值，也是一种可行的方法。如例7。 * 为节省计算机存储单元，我们采用每500高斯存入一个at值。为了分析使用几阶插值公式合适，应先列出差分表。 * * 这个结果与直接查表得到的值相同，说明用此法在计算机上求at值是可行的。 * * * 八．差商与牛顿插值公式 例6 给出f(x)的函数表，求四次牛顿插值多项式，并由此计算f(0.596)的近似值。 解： p4(x)=0.41075+1.116*(x-0.4)+0.28*(x-0.4)(x-0.55)+ 0.1973*(x-0.4)(x-0.55)(x-0.65)+ 0.03146*(x-0.4)(x-0.55)(x-0.65)(x-0.8) f(0.596)≈ p4(0.596)=0.63192 |R4(0.596)| ≈ |f(x0,x1,x2,x3,x4,x5)*(0.596-0.4)(0.596-0.55) (0.596-0.65)(0.596-0.8)(0.596-0.9) | 0.15 *10-7 5．示例 九．差分与等距节点插值公式 * 1．差分的定义 设节点等距，则可记xk=x0+kh，(k=0,1,…,n) 定义 差分指相邻节点数据的偏差。 Δyk=yk+1-yk ▽yk=yk-yk-1 符号Δ，▽，δ分别称为向前差分算子，向后差分算子及中心差分算子。 步长 向前差分 向后差分 中心差分 九．差分与等距节点插值公式 * 1．差分的定义 已知函数值 f(xi) = yi 一阶差分：Δyi=yi+1-yi 二阶差分：Δ2yi=Δyi+1-Δyi=yi+2-2yi+1+yi m阶差分：Δmyk=Δm-1yk+1-Δm-1yk m阶向后差分：▽myk=▽m-1yk-▽m-1yk-1 中心差分：δyk用到yk+1/2和yk-1/2，这两个值不是函数表上的值。如果用函数表上的值，一阶中心差分应写成： δyk+1/2=yk+1-yk，δyk-1/2=yk-yk-1 二阶中心差分： δ2yk=δyk+1/2-δyk-1/2 九．差分与等距节点插值公式 * 1．差分的定义 yi y0 y1 y2 y3 y4 ┇ Δ Δy0 Δy1 Δy2 Δy3 ┇ Δ2 Δ2y0 Δ2y1 Δ2y2 ┇ Δ3 Δ3y0 Δ3y1 ┇ Δ4 Δ4y0 ┇ 九．差分与等距节点插值公式 * 1．差分的定义 yi y0 y1 y2 y3 y4 ┇ Δ Δy0 Δy1 Δy2 Δy3 ┇ Δ2 Δ2y0 Δ2y1 Δ2y2 ┇ Δ3 Δ3y0 Δ3y1 ┇ Δ4 Δ4y0 ┇ 九．差分与等距节点插值公式 * 2．差分与差商的关系 九．差分与等距节点插值公式 * 3．等距节点插值公式 令x=x0+th，于是 (x-x0) (x-x1)…(x-xk)=t(t-1)(t-2)…(t-k)hk+1 则牛顿插值公式可写作： 称为有限差公式 ，或牛顿前插公式。 余项