数值分析 Ch2-插值法.pptVIP

1. 插值概念 2-1 插值多项式的存在唯一性 证明 插值条件(2.2)等价于线性方程组 例 给定数据如下 ： 后面为链接 §4 Piecewise Polynomial Approximation 7-2 分段线性插值 /* piecewise linear interpolation */ 在每个区间 上，用1阶多项式 (直线) 逼近 f (x): 记 ，易证：当 时， 一致 失去了原函数的光滑性。 ? 分段Hermite插值 /* Hermite piecewise polynomials */ 给定 在 上利用两点的 y 及 y’ 构造3次Hermite函数 导数一般不易得到。 How can we make a smooth interpolation without asking too much from f ? Headache … 7-3 分段三次Hermite插值 /* Hermite piecewise polynomials */ §8 三次样条 /* Cubic Spline */ 分段插值法具有一致的收敛性, 但它只保证插值函数整体的连续性, 但在连接处不一定光滑，不能够满足精密机械设计（如船体、飞机、汽车等的外形曲线设计）对函数光滑性的要求。 早期的工程技术人员在绘制给定点的曲线时，使用一种具有弹性的细长木条（或金属条），称之为样条（Spline），强迫它弯曲通过已知点。弹性力学理论指出样条的挠度曲线具有二阶连续的导函数，并且在相邻给定点之间为三次多项式，即为数学上的三次样条插值曲线。 定义 　　　设 。三次样条函数 , 且在每个 上为三次多项式 /* cubic polynomial */。若它同时还满足 ，则称为 f 的三次样条插值函数 /* cubic spline interpolant */. 注：三次样条与分段 Hermite 插值的根本区别在于S(x)自身光滑，不需要知道 f 的导数值（除了在2个端点可能需要）；而Hermite插值依赖于f 在所有插值点的导数值。 f(x) H(x) S(x) 三次样条插值函数 在每一个小区间上是3次的多项式, 在整个插值区间上有4n个系数. 且有4n-2个约束: 内节点 边界节点 要确定4n个系数，还需附加2个约束条件. 常用的约束条件有以下三类： 此时一般有 成立. ③.周期性边界条件 , ②.弯矩边界条件 特别的称 为自然边界条件. ①.转角边界条件 ? 构造三次样条插值函数的三转角方程 /* method of bending moment */ ? 构造三次样条插值函数的三转角方程 /* method of bending moment */ ? 构造三次样条插值函数的三转角方程 /* method of bending moment */ ? 构造三次样条插值函数的三转角方程 /* method of bending moment */ 1.线性： 2.差商可以表示为函数值的线性组合： 3. 对称性：由2知，差商的值与节点的顺序无关！ 4. 差商的另一种定义：由2，3及均差定义可得 §4 Newton’s Interpolation §4 Newton’s Interpolation ? 牛顿插值 /* Newton’s Interpolation */ … … … … Nn(x)—n次多项式，满足： Nn(xi)= f(xi) Rn(x)—插值余项，满足Rn(xi)=0，i=0,…,n ai = f [ x0, …, xi ] (1) (2) (n) (n)?(n-1) ?… ?(2) ?(1) §4 Newton’s Interpolation