数值分析第章要点.ppt

* * 二、三弯矩方程 Mi来求S(x)的方法称为三弯矩法。 为参数,这种通过确定 设 在[xi，xi+1]上是一次多项式且可表示为 积分两次并利用 S(xi)=yi S(xi+1)=yi+1 定出积分常数得 对 * 为了使 S(x) 成为所求的三次样条，同时为了 确定参数 Mi ，对S(x) 求导得 * 所以 得 (i=1,2,...,n-1) 由 * 1、I 型样条函数 ，由公式 (6-29)、(6-30) 得 已知 * 即 * 合并式(6-31)、(6-33)有 从中解出Mi(i=0,1,...,n)代入式(6-28) 得 I 型三次样条S(x)。 * 2、II 型样条函数 M0、Mn已知，故方程组(6-31)为 从中解出 Mi(i=1,2,...,n-1) 代入式(6-28) 得 II 型三次样条S(x)。 * 3、III 型样条函数 ，则由公 式(6-29)、(6-30)得 已知 * 整理得 其中 * * 三、三转角方程 设 定 为参数，这种通过确 来求S(x)的方法叫三转角法。 用前面介绍的分段埃尔米特插值，得到在 上S(x)的表达式为 * 所以 同理 * 其中： (6-42) 同三弯矩方程一样，有三种条件： 1、已知 则方程组化为： 由S(x)二阶连续可微，即 * 2、已知 由 可得 由 可得 * 于是有 即矩阵形式为： * 其中 3、已知 则有： * * 则其I 型和II 型三次样条插值函数以及导数的 误差有如下估计式 设 f(x) 在[a，b]上有直到四阶的连续导数， 四、误差估计 * 0 1 2 -1 0 1 0.3679 1.0000 2.7182 0.3679 2.7182 在[－1，1]上已知 三个节点上的函数值和边界点上的导数值，给出三次样条插值多项式，并估计误差 例题 * 例题 * * * * * * * * 基函数满足下列插值条件 H3(x)称为三次埃尔米特插值多项式。 * * 条 件 函 数 函数值 导数值 x0 x1 x0 x1 ?0(x) 1 0 0 0 ?1(x) 0 1 0 0 ?0(x) 0 0 1 0 ?1(x) 0 0 0 1 表2-3 即 * 即 插值点的Lagrange 一次基函数. * 可得满足条件的三次埃尔米特插值多项式为 定理3 满足条件式 的三次埃尔米特插值多项式存在且唯一。 * 二、误差估计 定理4 设f(x)在包含x0、x1的区间[a,b]内存在四阶导数，则当x∈[a,b]时有 且与x有关) * 例1 已知f(x)=x1/2在X=121和144时的函数值及其一阶导数的数据见下表,用埃尔米特插值公式计算1251/2的近似值,并估计其截断误差. * * 得 由 可求得 * 例2 * * 第五节 分段低次插值 * * * * 一、分段线性插值 * 由定义知 在每个小区间 上可以表示为 * 分段线性插值多项式 基函数 * 分段线性插值在区间[xi，xi+1]上的余项估计式为 一致收敛定理：设 ， 在[a,b]上一致收敛到 。 * 例1 已知函数 在区间 上取等 距插值节点（如下表）, 求区间 上分段线性插值函数，并利用它求出 的近似值 1 0.03846 0.05882 0.1 0.2 0.5 1 5 4 3 2 0 解 在每个小区间 上， * 于是 * 二、分段三次Hermite插值 [xi，xi+1]作埃尔米特插值，得 在每个子区间 * 显然，H3(x) 满足插值条件 H3(xi)=yi， 阶导数连续。 在节点 处一 * 分段三次埃尔米特插值在区间[xi，xi+1]上的余 项估计式为 * 一致收敛定理 * -0.16 -0.5 0 0.2 0.5 1 2 1 0 例3 已知函数 ,在区间 上取等距插值节点（如下表),求区间 上分段三次 埃尔米特插值函数,并利用它求出 的近似值 * 解 ,由式(4) ,在每 个小区间 上 于是 * * 定义 给定区间[a，b]的一个划分 a=x0＜x1 ＜...＜xn-1＜xn=b，如果函数S(x)满足： ①在每个小区间[xi，xi+1](i=0,1,...,n-1) 上是三