2.插值法-2概述.ppt

2.5 埃尔米特插值 两点三次插值多项式 分段三次埃尔米特插值 2.7 三次样条插值 2.7.2 样条插值函数的建立 2.7.3 误差界与收敛性 小结 插值的概念 拉格朗日插值 均差与牛顿插值公式 等距节点插值 分段低次插值 埃尔米特插值 三次样条插值函数 类似地可求出 在区间 上的表达式，从而得 利用 可得 （7.10） 其中 （7.11） 对第一种边界条件（7.3），可导出两个方程 （7.12） 如果令 （7.13） 那么（7.10）及（7.12）可写成矩阵形式 对第二种边界条件（7.4），直接得端点方程 （7.14） 如果令 , 也可以写成（7.13）的矩阵形式. 则（7.10）和（7.14） 对于第三种边界条件（7.5），可得 其中 （7.15） （7.10）和（7.15）可以写成矩阵形式 （7.16） （7.13）和（7.16）是关于 的三对角 方程组， 在力学上解释为细梁在 截面处的弯矩，称 为 的矩， （7.13）和（7.16）的系数矩阵为严格对角占优阵，有 唯一解, 求解方法可见第5章第4节追赶法，将解得结果代入 （7.8）的表达式即可. 方程组（7.13）和（7.16）称为三弯矩方程. 设 为定义在 上的函数，在节点 试求三次样条函数 ，使它满足边界条件 例5 上的值如下： 由此得矩阵形式的方程组（7.13）为 解 由（7.11）及（7.12） 求解得 代入（7.8）得 （曲线见图2-6） 图2-6 给定函数 节点 用三次样条插值求 取 直接上机计算可求出 在表2-6所列各点的值. 例6 下图是用Matlab完成的样条插值（附程序）： 附：样条插值程序 n=11; m=61; x=-5:10/(m-1):5; y=1./(1+x.^2); z=0*x; x0=-5:10/(n-1):5; y0=1./(1+x0.^2); y1=interp1(x0, y0, x, ’spline’); plot(x, z, ’r’, x, y, ’k:’, x, y1, ’r’) gtext(‘Spline’), gtext(‘y=1/(1+x^2)’) title(‘Spline’) 注：interp1(x0, y0, x, ’spline’)为Matlab中现成的样条插值程序. 也可以将三种插值结果画在一起： 定理4 设 为满足第一种或第二 则有估计式 （7.17） 其中 种边界条件（7.3）或（7.4）的三次样条函数， 令 这个定理不但给出了三次样条插值函数 的误差 估计. 还说明了当 时， 及其一阶导数 和 二阶导数 均分别一致收敛于 ， 及 * 有些实际的插值问题不但要求在节点上函数值相等， 下面只讨论函数值与导数值个数相等的情况. 满足这种要求的插值多项式就是埃尔米特插值多项式. 而且还要求对应的导数值也相等，甚至要求高阶导数也相等. (5.1) 这里共有 个插值条件，可唯一确定一个次数不超过 的多项式 ， 问题是求插值多项式 ， 设在节点 上， 现在仍采用求拉格朗日插值多项式的基函数方法. 满足条件 其形式为 问题描述 将满足条件（5.1）的插值多项式 写成用插 值基函数表示的形式 (5.3) 先求出 个插值基函数 及 ， 每一个基函数都是 次多项式， 且满足条件 (5.2) 令 由条件（5.2），有 下面的问题就是如何求出这些基函数 及 利用拉格朗日插值基函数 解出 由于 整理得 于是 (5.4) 两端取对数再求导，得 同理，可得 (5.5) 可以证明满足条件（5.1）的插值多项式是惟一的. 用反证法，假设 及 均满足条件（5.1）， 这样， 有 重根，但 是不高于 次的多 项式， 于是 在每个节点 上的值及导数值均为零，即