数值分析-三次样条插值.ppt

华长生制作 三次样条插值 3.6 三次样条插值 什么是样条: 是 指飞机或轮船等的制造过程中为描绘 出光滑的外形曲线(放样)所用的工具 样条本质上是一段一段的三次多项式拼合而成的曲线 在拼接处,不仅函数是连续的,且一阶和二阶导数也是连续的 1946年,Schoenberg将样条引入数学,即所谓的样条函数 一、三次样条插值函数 定义1. ------(1) 二、三次样条插值多项式 ------(2) ------(3) ------(4) 少两个条件 并且我们不能只对插值函数在中间节点的状态进行限制 也要对插值多项式在两端点的状态加以要求 也就是所谓的边界条件: 第一类(一阶)边界条件: 第二类(二阶)边界条件 第三类(周期)边界条件 ------(5) ------(6) ------(7) 加上任何一类边界条件(至少两个)后 一般使用第一、二类边界条件, 即 ------(8) 或 常用第二类边界条件 三次样条插值的M方法 设f (x)是定义在 [a, b]区间上的一个二次连续可微函数， ?为分划： S? (x)在 [xi-1, xi] 上的表达式为： 令 i = 0, 1, 2,…, n 在每一个小区间[xi-1, xi] i = 1,…, n 上都是三次多项式， （6.7） 其中 ，将（6. 7）两次积分得： Ai 和Bi 为积分常数。 因为 所以它满足方程： （6.8） 求 Mi，确定S (x)的表达式。微分（6.8）式 于是 由 得 各项除以hi + hi+1，并记 则（6.9）可以写为 （6.9） 端点条件 最后一个方程。若取 M0 = Mn=0，称为三次自然样条。 （1）给定 补充（6.9）的第一个和 有 （2）给定两端点导数值 分别补充为方程组（6. 9）的第一个和最后一个方程组。 解方程组 经补充后的方程组（6. 9）为 其中，对端点条件（1），有 （6.9） 对端点条件（2），有 （6. 10 ）是一个三对角方程组，可用追赶法解之。 此方程组系数严格对角占优 ！从而存在唯一解。 求出了Mi (i = 0, 1,…, n)，也就求得了S (x)在各个 小区间的表达式Si (x)（i = 0, 1, 2,…, n） 若取等距节点 hi = h, i = 1,…, n –1 例1. 对于给定的节点及函数值 解： 由M关系式 将上述结果代入 定理 . 最后，介绍一个有用的结果 ------(9) 三次样条插值的m方法 加以整理后可得 ------(10) ------(11) 由条件 由于以上两式相等,得 * * * *