数值分析205(免费阅读).ppt

湘潭大学数学与计算科学学院 * 例2 求三次插值样条函数 满足插值条件(见下表) 0 1 2 3 0 0 0 0 及边界条件 解 利用三弯矩方程求解 设 为待定参数， 则关于内剖分节点有 湘潭大学数学与计算科学学院 * (5.14) 其中 再利用边界条件，可得 (5.15) 联立(5.14)和(5.15)，可得如下线性代数方程组 湘潭大学数学与计算科学学院 * (5.16) 由此可以解得 利用(5.6)，可得 的分段表达式为 湘潭大学数学与计算科学学院 * (5.6) 湘潭大学数学与计算科学学院 * 取 为等距剖分， 并对 定义如下两种范数 即 总假设 并记 三次样条插值问题误差估计 湘潭大学数学与计算科学学院 * 定理5.4 设 为Ⅰ型三次插值样条 问题的误差函数， 则有 (5.22) 其中 上一页 下一页 湘潭大学数学与计算科学学院 * §5 分段低次多项式插值 一、分段一次多项式插值 二、分段三次多项式插值 三、三次插值样条函数的基本性质和误差估计 湘潭大学数学与计算科学学院 * 问题的引入： 如果不考虑舍入误差等因素对插值公式的影响 则可以通过增加插值节点的个数 构造相应的高次插值多项式来减少插值误差. 问题：高次插值多项式存在以下两方面的缺陷. 1、收敛性问题 2、数值稳定性问题 湘潭大学数学与计算科学学院 * 1、收敛性问题 Runge现象 例1 设 对于任意正整数n， 的n次Lagrange插值多项式为 其中，等距节点 湘潭大学数学与计算科学学院 * - 5 - 4 - 3 - 2 - 1 0 1 2 3 4 5 - 0.5 0 0.5 1 1.5 2 2.5 n 越大， 端点附近抖动 越大，称为 Runge 现象 Ln(x) ? f (x) ? 湘潭大学数学与计算科学学院 * %marunge.m function marunge() xx=-5:0.05:5; y=1./(1+xx.^2); x1=-5:2:5; y1=1./(1+x1.^2); x2=-5:1:5; y2=1./(1+x2.^2); yy1=malagr(x1,y1,xx); yy2=malagr(x2,y2,xx); plot(xx,yy1,-.); hold on plot(xx,yy2,r-.); plot(xx,y); text(1.5,0.5,L5); text(4.0,1.5,L10); 湘潭大学数学与计算科学学院 * 2、数值稳定性问题 高次插值多项式的计算需要用到高阶差分或差商 差分的误差传播会随阶数的提高越来越严重 利用高次插值多项式来减少插值误差受到了限制 湘潭大学数学与计算科学学院 * 常用的方法： 用分段低次插值多项式来代替高次插值多项式 即将插值函数空间取为分段多项式空间 首先引入区间[a, b]的一个分划(或网格剖分) 称 为剖分节点， 为内剖分节点 为第j 个剖分单元. 小区间 湘潭大学数学与计算科学学院 * 记 和 分别表示 在节点 处的函数值、一阶和二阶 被插函数 导数值. 设关于剖分 的分段多项式插值函数空间为 (5.1) 湘潭大学数学与计算科学学院 * 表示在每个剖分单元为k次多项式 且在整个定义区间 为关于剖分 重数为 的k次多项式 称 [a,b]上l 阶连续可微的函数的全体. 样条空间. (5.1) 特别，当 时， 并称其为(完全)样条函数空间. 为 简记 下面主要讨论 和3的情形. 湘潭大学数学与计算科学学院 * 一、分段一次多项式插值 取插值函数空间为 它表示在每个剖分单元 为线性多项式且在区间[a,b]上连续的函数全体. 注意: (1) 在每个剖分单元上一次多项式空间 的自由度(或维数)为2， 所以其总数为2n； (2) 由于 中的函数在每个内节点上连续, 个约束条件； 所以共有 的自由度为： 湘潭大学数学与计算科学学院 * 分段一次多项式插值问题： 满足插值条件 求 利用线性插值公式得到其分段表达式 这是一折线函数，相应的图象如下图所示： 湘潭大学数学与计算科学学院 * 湘潭大学数学与计算科学学院 * 利用线性插值多项式的误差余项公式 (4.2) 若 在[a,b]上二阶连续可微， 时， 则当 有 从而 湘潭大学数学与计算科学学院 * 由于 于是，在整个区间[a,b]上有 其中 (5.2) 湘潭大学数学与计算科学学院 * 取插值函数空间为 插值函数空间的维数为： 插值问题： 满足插值条件 求 利用两点Hermite插值公式 (3.5) 二、分段三次Hermite多项式插值 湘潭大学数学与计算科学学院 * 的分段表达式 其中 而 湘潭