数值分析课件第六章_数值插值方法.ppt

例1. 解: 作为多项式插值，三次已是较高的次数，次数再高就有可能发生Runge现象，因此，对有n+1个节点的插值问题，我们需要使用分段插值等其他方法。 二. 添加“修正项”的算法 先求出不带导数的牛顿插值 ，然后每加一个导数就添加一项 “修正项 ”，用这种方法可方便地求出更一般的埃尔米特插值多项式 P137 习题五：5 ，11，13，16 本章作业 2. 分段二次插值的误差估计 由于 那么分段二次插值L2(x)的余项为： 例: 解: (1) 分段线性Lagrange插值的公式为 同理 (2) 分段二次Lagrange插值的公式为 样条：是 指飞机或轮船等的制造过程中为描绘出光滑的外形曲线(放样)所用的工具。 样条本质上是一段一段的三次多项式拼合而成的曲线，在拼接处，不仅函数是连续的，且一阶和二阶导数也是连续的。 1946年，Schoenberg将样条引入数学，即所谓的样条函数。 因分段线性插值导数不连续，为克服在节点处不光滑的缺陷，故引入样条插值概念。 三、三次样条插值 一、三次样条插值函数的定义 定义： 给定区间[a,b]上的一个划分： a = x0 x1…xn=b， 已知函数f(x)在点xj上的函数值为 f (xj) = yj, ( j= 0,1,2,···,n) 如果存在分段函数 (1)S(x)在每一个子区间[xj-1 , xj ] ( j= 0,1,2,···,n)上是一个三次多项式； (2) S(x)在每一个内接点xj ( j= 1,2,···,n-1)上具有直到二阶的连续导数； 则称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，bj，cj，dj} (j=1,2,…,n)。 由(2)知，S(x)，S`(x)，S``(x)在内节点x1,x2,…,xn-1上连续，则 j=1,2,…,n-1 可得3n-3个方程，又由条件(3) j=0,1,…,n 得n+1个方程，共可得4n-2个方程。 要确定4n个未知数，还差两个方程。 通常在端点x0 =a, xn =b处各附加一个条件，称边界条件，常见有两种： (1)二阶边界条件： (2)一阶边界条件： 共4n个方程，可唯一地确定4n个未知数。 例 已知f(x)：f(-1)=1,f(0)=0,f(1)=1，求f(x)在[-1,1]上的三次自然样条插值函数， f(x) 在-1,1处的二阶导数为0。 解 设 由插值条件和函数连续条件得： 由一阶及二阶导数连续得： 由自然边界条件得： 联立上面8个方程，求解得 故 用二阶导数值构造三次样条插值函数 （又称三弯矩算法） 记 hj=xj-xj－1， (j=1,2,…,n) 设S``(xj)=Mj，(j=0,1,2,…,n) 因为S(x)是三次多项式，故S``(x)是线性函数。按线性插值公式可得 积分两次得 二、三次样条插值函数的建立 其中c1，c2为积分常数。 将S(xj-1)=yj-1， S(xj)=yj代入上式，可确定c1，c2。 故 上式称M表达式，只需确定Mi，即可确定三次样条插值函数。 将M表达式两端对x求导，得 令x=xj，得左导数 令x=xj-1，得右导数 故 因S(x)一阶导数连续，即 整理得 记 则得Mj的n-1个方程： (j=1,2,…,n-1) 对固定边界条件： 得 或记为 其中 则求Mi的方程的矩阵形式为 上述方程组是关于Mj(j=0,1,…,n)的三对角方程组。 该方程组是严格对角占优的三对角方程组，有唯一解，可用追赶法求解。 上述方程组是关于Mj(j=1,…,n-1)的三对角方程组。 Mj在力学上解释为细梁在xj截面处的弯矩，称为S(x)的矩，故称三弯矩方程组。 该方程组是严格对角占优的三对角方程组，有唯一解，可用追赶法求解。 对于二阶边界条件： 6.5 埃尔米特插值（Hermite） Newton插值和Lagrange插值虽然构造比较简单，但都存在插值曲线在节点处有尖点，不光滑，插值多项式在节点处不可导等缺点。为使插值曲线更好地逼近原曲线, 除在节点处的函数值相等 外，还要求节点处的一阶、 二阶导数等也相等 （从而有公切线、 相同的凹凸性等) 埃尔米特插值的基本思想为 设a ≤ x0 x1……xn≤b上， (j=0,1,2…,n)