第七节样条插值.pptVIP

* §7 三次样条插值 上面讨论的分段低次插值函数都有一致收敛性，但光滑性较差，对于像高速飞机的机翼形线，船体放样等型值线往往要求有二阶光滑度，即有二阶连续导数，早期工程师制图时，把富有弹性的细长木条（所谓样条）用压铁固定在样点上，在其他地方让它自由弯曲，然后画下长条的曲线，称为样条曲线。它实际上是由分段三次曲线并接而成，在连接点即样点上要求二阶导数连续，从数学上加以概括就得到数学样条这一概念。下面我们讨论最常用的三次样条函数。 * 设在区间[a,b]上取n+1个节点a= x0 x1 … xn =b,其对应的函数值为f(xi)= yi , (i=0,1,2,…,n) , 现求一定义在[a,b]上的函数S(x),使其满足： 1) S(x)在每一个小区间[xi-1 , xi], (i=1,2,…,n)上为三次多项式, 2) S(x)在[a , b]上二阶连续可微, 即S(x) ∈C2[a, b], 3) S(xi) = yi , (i=0,1,2,…,n), 则称 S(x)为 f (x)的三次样条插值函数. 根据S (x ) 在 [a ,b ] 上二阶导数连续，在节点x j ( j = 1 , , 2 , n ? 1) L处应满足连续性条件: * 通常可在区间[a ,b ] 端点x 0 = a , x n = b 上各加一个条件 （称为边界条件），可根据实际问题的要求给定。常见的有以下三种： 1° 已知两端的一阶导数值，即 2° 两端的二阶导数已知，即 其特殊情况 ， 称为自然边界条件。 * 为确定S(x)在各区间上的表达式, 自然会想到前节介绍的分段埃尔米特插值.在每一小区间上,分段埃尔米特插值也恰为三次多项式,不过在那里, f ′(xi)=y′ i 为已知的. 如果能利用S(x)在[a, b]上二阶连续可微将y′ i确定,那么由(6.6), S(x) 便可唯一地表示出来了. 为此先记mi =f ′(xi), (i=0,1,2,…,n), 在区间[xi-1 , xi]上, S(x)的表达式为： 3° 当f (x ) 是以xn ? x0为周期的周期函数时，则要求 S(x ) 也是周期函数。这时边界条件应满足: 而此时 y n = y0 ,这样确定的样条函数S(x ) ，称为周期样条函数 三转角方程 * 样条插值函数的建立 记hi-1= xi - xi-1 , 求 S (x)在 [xi-1 , xi ]上的二阶导函数 S”(x) * 再利用S(x)在[xi , xi+1]上的表达式可计算出 由S(x)二阶连续可微,即S (xi-) = S (xi+),得 记： 将上式整理得方程组 * 这是关于n+1个未知量mi,(i=0,1,…,n)的n-1个线性方程组,该方程组有无穷多组解,在实际问题中,往往根据具体情况补充两个附加条件—通常称为端点条件, 便可唯一确定一组解. 常见的端点条件有： 1) 曲线在两端点x0 , xn 处的导数值为已知的, 即f′(x0)= m0, f′(xn)=mn,方程组为n-1个未知数, n-1个方程,有唯一解. 2) 函数f (x0)在两端点x0 , xn 处的二阶导数为0, 即 f “ (x0)= f ” (xn) =0, 即 S (x0)= S (xn) =0 由S (x0+ ) =0可得 * 得方程组 也有唯一解. 由S (xn- ) =0可得 3) 周期端点条件( 当f (x0)=f (xn)时 )：  S(x0)=S(xn), S′(x0)=S′(xn), S (x0)=S (xn)于是有： * 与前面n-1个方程联立, 也可唯一地解出m0,m1, …,mn. 对于前两种端点条件所对应的方程组可化为下列称之为三对角方程组的形式： 该方程组的系数矩阵具有严格对角优势(对角线元素的绝对值大于该行其它元素绝对值之和),因而是非奇异的,所以方程组有唯一的解,其解法可用追赶法进行. * 最后,把计算三次样条函数的步骤归纳如下： (1) 按所给的边界条件及所给节点(xi ,yi )计算?i ,μi , di (2) 在给定的边界条件下解方程组 (6.27) 计算出mi (i=0,1,2,…,n) 将所求得的mi 代入分段插值公式, 求出各小区间 [xi, xi+1 ]上的样条函数 Si(x), (4) 计算插值区间[ a, b ]上的样条插值函数 S(x) 的值. * 三弯矩方程: 三次样条插值函数 S (x )可以有多种表达方法，有时用二阶导数值 S “ (x j )