[法学]第2章插值法.ppt

(5.5) 故 (5.6) 将（5.5）中的 改为 ，即得 在子区间 上的表达式 ，并由此得： 将(5.6),(5.7)代入(5.4)整理后得 (5.7) 两边同乘以 ，即得方程组 若记 (5.8) 则所得方程组可简写成 即 (5.9) 这是一个含有 n+1个未知数、n-1 个方程的线性方程组。要确定 的值，还需用到边界条件。在第 (1) 种边界条件下，由于 和 已知，可以得到包含 另外两个线性方程。由(5.5)知， 在子区间 上的导数为 故由条件 立即可得 即 (5.10) 同理，由条件 可得 （5.11） 将(5.9)、(5.10)、(5.11)合在一起，即得确定 的线性方程组： (5.12) 其中 (5.13) 已知，在方程组(5.13)中实际上只包含有n-1个未知数 ，并且可以改写成 在第(2) 种边界条件下，由 (5.14) 在第(3)种边界条件下，由 直接可得 (5.15) 由条件 可得 注意到 和 ，上式整理后得 若记 则所得方程可简写成 (5.16) 将 (5.9) 、(5.15) 、(5.16) 合在一起，即得确定 的线形方程组： (5.17) 利用线性代数知识，可以证明方程组 (5.12) 、(5.14) 及(5.17)的系数矩阵都是非奇异的，从而都有唯一确定的解。 针对不同的边界条件，解相应的方程组(5.12) 、 (5.14)或(5.17) ，求出 的值，将它们代入（5.3），就可以得到 在各子区间上的表达式。综上分析，有 解 因为0.12介于0.1与0.2之间，故取 ，此时 。为求 , 构造差分表2—2。表中长方形框中各数依次为 在 处的函数值和各阶差分。若用线性插值求sin﹙0.12﹚的近似值，则由前插公式﹙3.2﹚立即可得 用二次插值得 用三次插值得： 因 很接近，且由差分表2—2可以看出，三阶差分接近于常数（即 接近于零），故取 作为 的近似值，此时由余项公式（3.3）可知其截断误差 3.2 向后差分与牛顿向后插值公式 在等距节点 下， 除了向前差分外，还可引入向后差分和中心差分，其定义和记号分别如下： 在点 处以h为步长的一阶向后差分和m阶向后差分分别为 在 点处以为步长的一阶中心差分和m阶中心差分分别为 其中 　 各阶向后差分与中心差分的计算，可通过构造向后差分表与中心差分表来完成?参见表2－２?。 利用向后差分，可简化牛顿插值多项式（３.1），导出与牛顿前插公式?3.2?类似的公式，即，若将节点的排列次序看作 ，那么?３.1)可写成 ? 根据插值条件 , 可得到一个用向后差分表示的插值多项式 其中t0，插枝多项式(3.4)称为牛顿向后插值公式，简称后插公式。它适用于计