第二章 插值法-数值分析.pptVIP

1. 插值概念 2-1 插值多项式的存在唯一性 证明 插值条件(2.2)等价于线性方程组 例 给定数据如下 ： 后面为链接 ? 分段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 */ ? 构造三次样条插值函数的三转角方程 /* method of bending moment */ §4 Newton’s Interpolation 注：? 由唯一性可知 Nn(x) ? Ln(x)， 只是算法不同，表达形式不同，故其余项也相同，即 ? 实际计算过程为 f (x0) f (x1) f (x2) … f (xn?1) f (xn) f [x0, x1] f [x1, x2] … … … … f [xn?1, xn] f [x0, x1 , x2] … … … … f [xn?2, xn?1, xn] f [x0, …, xn] f (xn+1) f [xn, xn+1] f [xn?1, xn, xn+1] f [x1, …, xn+1] f [x0, …, xn+1] 增加 如果精度不够，增加节点xn+1,同时表中增加一行，三角形斜边上即为所要求的各次项系数。 §1 Lagrange Polynomial 练习： 设f(x)=3x2+5 ,求f[x0,x1,x2] , f[x0,x1,x2, x3] 练习： 练习： §4 Ne