数值计算Chapter4 插值与拟合.pptVIP

内插和外推 §3 分段低次插值 三种插值的比较 拉格朗日插值（高次多项式插值）： 曲线光滑；误差估计有表达式 收敛性不能保证（振荡现象） 用于理论分析，实际意义不大 分段线性插值： 收敛性良好 只用两个节点，且线性，简单实用 曲线不光滑 三次样条插值:(*) 曲线2阶光滑，收敛性有保证 实际中应用广泛 误差估计较难 则方程组可表示为 例4 给定数据表 用一次、二次前插公式计算 的近似值。 解 介于 与 之间，又用前插公式， 故取 又 所以 构造向前差分表： 1.6487 0.5 0.1569 0.0150 1.4918 0.4 0.0016 0.1419 0.0134 1.3499 0.3 0.1285 1.2214 0.2 用一次前插公式计算: 用二次前插公式计算: 例：在[?5, 5]上考察 的Ln(x)。取 Ln(x) ? f (x) ? - 5 - 4 - 3 - 2 - 1 0 1 2 3 4 5 - 0.5 0 0.5 1 1.5 2 2.5 n 越大， 端点附近抖动 越大，称为 Runge 现象 故实际计算时，常常把插值区间分成若干个小区间，在每个小区 间上使用低次插值，这就是分段低次插值。 给定节点 处的函数值 后， 要计算 在 处的近似值， 可先选取两个节点 和 ， 使 然后在 上做线性插值（即n=1）， 这种方法称为分段线性插值（折线插值）。 若取x附近的三个点进行二次插值，就叫做分段二次插值，又 叫做分段抛物插值。 即 一、分段线性插值与分段二次插值 (分段Lagrange插值) 例子 用分段线性插值法求插值,并观察插值误差. 1.在[-6,6]中平均选取5个点作插值(xch11) 4.在[-6,6]中平均选取41个点作插值(xch14) 2.在[-6,6]中平均选取11个点作插值(xch12) 3.在[-6,6]中平均选取21个点作插值(xch13) 二、三次样条插值 1、定义 对于给定数表 其中 若函数 满足： ⑴ 在每个子区间 上都是不高于三次 的多项式； ⑵ 在 上都连续； ⑶ 则称 为函数 关于节点 的三次样条插值函数. 由定义中的条件⑴， 记 为 在第i 个小区间 上的表达式， 设 从而，有 其中 是待定常数。 由条件⑵， 在各节点上连续， 在每一 个节点处的左、右极限相等，即 这样，给出了 个等式， 再由条件⑶，给出了 个等式， 这样，定义中共给出了 个等式， 还需两个条件， 就可以确定出4n个待定常数，从而得到 由于小区间共有n个，所以共有4n个待 定常数。 ⑵ 端点处的二阶导数值 特别地， 称为自然边界条件。 相应的，满足自然边界条件的三次样条插值函数称为自然样 条插值函数。 ⑶ 若 是以 为周期的函数，则可要求 都是以 为周期的函数， ⑴ 端点处的一阶导数值 相应的，满足这一边界条件的三次样条插值函数又称为周期 样条插值函数。 这两个条件通常在区间[a，b] 的端点上给出的已知，称为边界条件或端点条件. 边界条件类型很多，常用的有： 设 在各节点处的二阶导数为： 因为在每个子区间 上 高于三次的多项式， 都是不 所以 必是一阶多项式或常数， 记 有 在 上用线性插值，有 连续积分两次，利用 确定出积分中的 常数后，得 2、求三次样条插值函数 的方法（三弯距法） 由此，只要确定出n+1个待定参数 就可得到各 小区间上的 进而得到 为了求Mi（i=0，1，2…n），利用条件（2）可得关于Mi的方程组 其中 （i=1,2,…，n-1） 这是含n+1个未知数、n-1个方程的线性方程组。要求解，还需两个用两个边界条件。 记 解出 后，代入三对角方程组 解方程组，得