计算方法及MATLAB2-2分段插值.pptVIP

计算方法及MATLAB实现 * 第2章 插 值 法 第1节 引言 第2节 拉格朗日插值 第3节 均差与牛顿插值多项式 第4节 埃尔米特插值 第5节 分段低次插值 第6节 三次样条插值 * 2.5 分段低次插值 2.5.1 高次插值的龙格振荡现象 这是因为对任意的插值节点，当 时， 不 一定收敛到 . 在次数 增加时逼近 的精度不一定也增加. 根据区间 上给出的节点做出的插值多项式 * 所构造的拉格朗日插值多项式为 以 上的 个等距节点 考虑函数 ，它在 上的各阶导数均 存在. 令 则 (取中节点!) (箕舌线函数!) * 表2-5列出了 时的 的计算结果及 在 上的误差 (误差迅猛增长!) * 可见，随 的增加， 的绝对值几乎成倍增加. 这说明当 时 在 上是不收敛的. 龙格(Runge)证明了，存在一个常数 ，使得当 时， 而当 时 发散. * 取 根据计算画出 及 在 上的图形，见图2-5. 高次插值的龙格振荡现象 (翘尾巴!) * 从图上看到，在 附近, 与 偏离很远， 这说明用高次插值多项式 近似 效 果并不好. 通常不用高次插值，而用分段低次插值. * 下图是用Matlab完成的Lagrange插值： * 2.5.2 分段线性插值 所谓分段线性插值就是通过插值点用折线段连接起来 逼近 由于升高插值多项式的阶数有时并不能达到提高精度 的效果, 所以实际中往往采用分段插值的思想. 分段插值的基本思想是将插值区间划分为若干个小区间, 然后在每个小区间上做满足一定条件的低阶插值. (以直代曲!逐段逼近！) (源自伟大的割圆术！) * 设已知节点 上的函数值 记 求一折线函数 , 满足： 在每个小区间 上是线性函数. 则称 为分段线性插值函数. (步长和最大步长) (闭连续) (插值等式) (分段线性) * 由定义可知 在每个小区间 上可表示为 分段线性插值的误差可利用插值余项得到 其中 线性插值函数 * 由此还可以得到 故 在 上一致收敛到 . 在 上一致成立， * 下图是用Matlab完成的分段线性插值： * * * * 2.6 三次样条插值 * 定义3 若函数 且在每个小区间 上是三次多项式，其中 是给定节点， 2.6.1 三次样条函数 若在节点 上给定函数值 并成立 则称 是节点 上的三次样条函数(Cubic Spline). 计算方法及MATLAB实现