三次样条函数的构造方法.docVIP

三次样条函数的构造方法 第15卷第2期 2000年6月 江南学院学报 JO URNA L OF JIANGNAN COLL EGE V ol. 15N o. 2June ,2000 文章编号:1008-8547(2000) 02-0090-04 三次样条函数的构造方法 谢文龙 (江南学院数理部, 江苏无锡214063) 摘要:三次样条函数是一类在工程技术上应用十分广泛的插值函数, 而且它的构造也具有特色。若利用幂级数的泰勒展开形式, 则可以直接构造出在小区间上的三次样条函数s (x ) , 利用插值条件及在插值节点处的一阶与二阶导数的连续性可确定出s (x ) 中系数的关系式, 再加上两个边界条件通过解线性方程组即得三次样条函数的分段表达式, 最后估计了它的余项, 给出了误差限。关键词:三次样条函数, 线性方程组, 误差限中图分类号:O241. 4; T B112 文献标识码:A 在插值函数中有一类非常重要的函数———, 是分段低次多项式, 且光滑性较好, ; ) 的函数值及提供(n -1) () 。其中最常见的有三次样条插值函数[1], 如三弯矩方法, 三转角方法[2]等, 本文采用形, 直接利用函数的幂级数展开形式, 同样可以得到所要求的三次样条插值函数。 若对区间[a ,b]作任一分划π:aΦx 0 23 s (x ) =y j +a 1(x -x j ) +a 2(x -x j ) +a 3(x -x j ) (1) 且s (x ) 满足条件s (x j ) =y j ,s (x j +1) =y j +1, 以下文中均记 h j =x j +1-x j 及f [xj , x j +1]=f [xj -1 , = x j +1-x j h j h j +h j-1 , x j , x j +1]= 。 1利用一阶导数值求样条函数 (x j ) =m j , s ′(x j +1) =m j +1, 其中m j , m j +1未知(欲求) , 若s ′ 则s (x ) 可写成 23 s (x ) =y j +m j (x -x j ) +a 2(x -x j ) +a 3(x -x j ) (2) 收稿日期:1999-06-12; 修订日期:2000-02-26作者简介:谢文龙(1964-) , 男, 江苏江阴人, 讲师. ?90? 由s (x j +1) =y j +1可得 m j +a 2h j +a 3h j 2=f [xj ,x j +1] (3) 对(2) 式求导得 (x ) =m j +2a 2(x -x j ) +3a 3(x -x j ) 2s ′ (x j +1) =m j +1可得由s ′ m j +2a 2h j +3a 3h j 2=m j +1 (4) 由(3) 、(4) 式联立方程可解得a 2, a 3并代入(2) 式, 得 (x ) =y j +m j (x -x j ) +h -j 1(3f [xj , x j +1]-2m j -m j +1) (x -x j ) 2s ′ 23 +h j -(-2f [xj , x j +1]+m j +m j +1) (x -x j ) (5) (x j +0) 的表达式若对(5) 式求导二次, 得到s ″ (x j +0) =s ″(x j ) =2h -j 1(3f [xj , x j +1]-2m j -m j +1) (6) s ″ (x j -0) 的表达式同理可得在区间[xj -1, x j ]上的三次样条函数s (x ) 及二阶导数值s ″ (x j -0) =2h -j -11(-3f [xj-1, x j ]+m j -1+2m j ) (7) s ″ (x j +0) =s ″(x j -0) , 化简并整理可得由样条函数的二阶导数连续性知s ″ m j -1+2m j +m j+1=3f [xj-1, x j ]+f [xj , x j ) h j +h j -1h j +h j -1h j +h j -1h j +h j -1 若记 μj =h j +h j -1, λj =h j +h j -1, d j =3j 1, j j f j , x +1则上式化为 (j =1,2, …, n -1) j j -1+2j j +1j , 常见的有以下两种: (x 0) =f 0″,s ″(x n ) =f n ″(1) 给出两个端点的二阶导数值, 即s ″ (x 0) =2h -01(3f [x0, x 1]-2m 0-m 1) =f 0″则由(6) 式知s ″ (8) 化简得同理由(7) 式知化简得 2m 0+m 1=d 0=3f [x0, x 1]- h 0f 0″2 1 (9) (x n ) =2h n --11(-3f [