第五章数值插值方法.pptVIP

第五章 插值方法 第五章 插值方法 插值的基本概念 Lagrange插值 分段低次插值 均差与Newton插值 Hermite插值 三次样条插值 5.1 插值的基本概念 例. 某地区某年夏季时节间隔 30 天的日出日落时间为 日照时间的变化设为 y x a0+ a1x + a2x2, 设 Pn x a0 + a1x + a2x2 + ……+ anxn 是y f x 在[a,b]上的n+1个互异节点x0,x1,…,xn的插值多项式，则求Pn x 问题归结为求系数a0,a1,…,an。 故Pn x 存在且唯一。 5.2 Lagrange插值 给定插值节点 x0, x1， y0 f x0 ，y1 f x1 . 求线性插值多项式L1 x a0+ a1x，使满足: L1 x0 y0 , L1 x1 y1. 显然 L x l0 x y0+l1 x y1+l2 x y2 满足条件 L2 xj yj j 0,1,2 设有n+1个互异节点x0 x1 … xn，且 yi f xi i 0,1,2…,n 构造Ln x ，使 Ln xj yj j 0,1,2,…,n 由n 1,2时的讨论可得 取x0 4, y0 2，x1 9, y1 3 ，x2 16, y2 4. 定理 设 f x 在[a,b]上具有n阶连续导数, 且 f n+1 x 存在,节点a ≤ x0 x1 … xn≤b， Ln x 是满足条件Ln xj yj j 0,1,2,…,n 的插值多项式，则对任何x ∈ [a,b]，插值余项 当 f x 是n次多项式时, Ln x f x 。即n次多项式的n次插值函数即为该n次多项式本身。 这说明截断误差很小。 设节点x0 x1 … xn，分段插值函数Hn x 在两个相邻节点构成的小区间[xj , xj+1 ] j 0,1,…n-1 上满足条件： 其中 1 S x 在每一个子区间[xj-1 , xj ] j 0,1,2,···,n 上是一个三次多项式； 2 S x 在每一个内接点xj j 1,2,···,n-1 上具有直到二阶的连续导数； 例 已知f x ：f -1 1，f 0 0，f 1 1，求f x 在[-1,1]上的三次自然样条插值函数。 P137 习题五：1 ，2，3，4，7，8 Lagrange 插值基函数 同理可得 将以上结果代入 得两个节点的三次Hermite插值公式： 二、三次Hermite插值的余项 定理： 设f x 在区间[a,b]上有定义，f x 在 a,b 内有4阶导数，H3 x 是满足插值条件 j 0,1 的三次Hermite插值函数，则对任意的x∈[a,b]，H x 的插值余项为 证明： 由 i 0,1 可知，x0，x1均为R3 x 的二重零点，因此可设 其中K x 待定 构造辅助函数 i 0,1 因此φ t 至少有5个零点。 连续4次使用Rolle定理可得，至少存在一点ξ∈[x0,x1]，使得 即 所以，两点三次Hermite插值的余项为 例1. 解: 作为多项式插值，三次已是较高的次数，次数再高就有可能发生Runge现象，因此，对有n+1个节点的插值问题，我们可以使用分段两点三次Hermite插值。 三、分段三次Hermite插值 用三次Hermite插值，当x? [xj , xj+1 ]时，有 5.6 三次样条插值 样条：是 指飞机或轮船等的制造过程中为描绘出光滑的外形曲线 放样 所用的工具。 样条本质上是一段一段的三次多项式拼合而成的曲线，在拼接处，不仅函数是连续的，且一阶和二阶导数也是连续的。 1946年，Schoenberg将样条引入数学，即所谓的样条函数。 因分段线性插值导数不连续，埃尔米特插值导数连续但需要已知，故引入样条插值概念。 （一）三次样条插值函数的定义： 定义： 给定区间[a,b]上的一个划分： a x0 x1 … xn b， 已知函数f x 在点xj上的函数值为 f xj yj, j 0,1,2,···,n 如果存在分段函数 则称S x 为节点x0,x1,…, xn 上的三次样条函数。 若S x 在节点x0,x1,…, xn 上还满足插值条件： 3 S xj yj j 0,1,2,···,n 则称S x 为三次样条插值函数。（即全部通过样点的二阶连续可微的分段三次多项式函数） 满足下述条件： （二）三次样条插值多项式的确定： 由 1 知， S x 在每一个小区间[xj-1 , xj ] 上是一三次多项式，若记为Sj x ，则可设 要确定函数S x 的表达式，须确定4n个未知系数 aj，