三次样条函数计算方法三次样条函数-Read.PPT

华长生制作 第1章 插值 1.6 三次样条函数 1.6 三次样条函数 什么是样条: 是 指飞机或轮船等的制造过程中为描绘 出光滑的外形曲线(放样)所用的工具 样条本质上是一段一段的三次多项式拼合而成的曲线 在拼接处,不仅函数是连续的,且一阶和二阶导数也是连续的 1946年,Schoenberg将样条引入数学,即所谓的样条函数 一、三次样条插值函数 定义1. ------(1) 二、三次样条插值多项式 ------(2) ------(3) ------(4) 少两个条件 并且我们不能只对插值函数在中间节点的状态进行限制 也要对插值多项式在两端点的状态加以要求 也就是所谓的边界条件: 第一类(一阶)边界条件: 第二类(二阶)边界条件 第三类(周期)边界条件 ------(5) ------(6) ------(7) 加上任何一类边界条件(至少两个)后 一般使用第一、二类边界条件, 即 ------(8) 或 常用第二类边界条件 ------(9) 三次样条插值的m关系式 加以整理后可得 ------(10) ------(11) 由条件 由于以上两式相等,得 ------(12) (12)式称为基本方程组 如果问题要求满足第一类(一阶)边界条件: ------(5) ------(5) 基本方程组(12)化为n-1阶方程组 ------(13) 即 将(13)式化为矩阵形式 ------(14) 这是一个三对角方程组 如果问题要求满足第二类(二阶自然)边界条件: ------(6) 由(11)式,可知 ------(15) ----(16) ------(17) ------(18) 与基本方程组(12)联合,并化为矩阵形式,得 -----(19) (19)式与(14)一样,都是三对角方程组,并且都严格对角占优 可以使用追赶法求解,并且解是唯一的 对于问题要求满足第三类(周期)边界条件 请同学们自己思考 现在回到(10)式 设f (x)是定义在 [a, b]区间上的一个二次连续 可微函数， ?为分划： S? (x)在 [xi-1, xi] 上的表达式为： 令 i = 0, 1, 2,…, n 在每一个小区间[xi-1, xi] i = 1,…, n 上都是三次多项式， （6.7） 三次样条插值的M关系式 其中 ，将（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） * * * *