第三讲三次样条函数

1、计算方法计算方法第第3讲讲 样条函数样条函数本讲主要问题本讲主要问题一、一、一、一、一、一、样条函数样条函数样条函数样条函数样条函数样条函数二、三次样条插值二、三次样条插值二、三次样条插值二、三次样条插值二、三次样条插值二、三次样条插值三、三次样条函数的构造三、三次样条函数的构造三、三次样条函数的构造三、三次样条函数的构造三、三次样条函数的构造三、三次样条函数的构造 分段插值存在着一个缺点分段插值存在着一个缺点, 就是会导致插就是会导致插值函数在子区间的端点值函数在子区间的端点(衔接处衔接处)不光滑不光滑, 即导即导数不连续数不连续, 对于一些实际问题对于一些实际问题, 不但要求一阶导数不但要

2、求一阶导数连续连续, 而且要求二阶导数连续而且要求二阶导数连续. 为了满足这些要求为了满足这些要求, 人们人们引入了样条插值的概念引入了样条插值的概念. 所谓所谓所谓所谓所谓所谓 “样条样条样条样条样条样条” ( ( (splinesplinespline) ) )是工程绘图中的一种工具是工程绘图中的一种工具是工程绘图中的一种工具是工程绘图中的一种工具是工程绘图中的一种工具是工程绘图中的一种工具, , , 它是有它是有它是有它是有它是有它是有弹弹弹弹弹弹性的细长木条性的细长木条性的细长木条性的细长木条性的细长木条性的细长木条. . . 绘图时绘图时绘图时绘图时绘图时绘图时, , , 用细木条连

3、接相近的几个结点用细木条连接相近的几个结点用细木条连接相近的几个结点用细木条连接相近的几个结点用细木条连接相近的几个结点用细木条连接相近的几个结点, , , 然后然后然后然后然后然后再进行拼接再进行拼接再进行拼接再进行拼接再进行拼接再进行拼接, , ,连接全部结点连接全部结点连接全部结点连接全部结点连接全部结点连接全部结点, , , 使之成为一条光滑曲线使之成为一条光滑曲线使之成为一条光滑曲线使之成为一条光滑曲线使之成为一条光滑曲线使之成为一条光滑曲线, , , 且在结且在结且在结且在结且在结且在结点处具有连续的曲率点处具有连续的曲率点处具有连续的曲率点处具有连续的曲率点处具有连续的曲率点处具

4、有连续的曲率. . . 样条函数就是对这样的曲线进行数学模拟得到的样条函数就是对这样的曲线进行数学模拟得到的样条函数就是对这样的曲线进行数学模拟得到的样条函数就是对这样的曲线进行数学模拟得到的样条函数就是对这样的曲线进行数学模拟得到的样条函数就是对这样的曲线进行数学模拟得到的. . . 它除它除它除它除它除它除了要求给出各个结点处的函数值外了要求给出各个结点处的函数值外了要求给出各个结点处的函数值外了要求给出各个结点处的函数值外了要求给出各个结点处的函数值外了要求给出各个结点处的函数值外, , , 只需提供两个边界点处只需提供两个边界点处只需提供两个边界点处只需提供两个边界点处只需提供两个边界

5、点处只需提供两个边界点处导数信息导数信息导数信息导数信息导数信息导数信息, , , 便可满足对光滑性的不同要求便可满足对光滑性的不同要求便可满足对光滑性的不同要求便可满足对光滑性的不同要求便可满足对光滑性的不同要求便可满足对光滑性的不同要求. . . 定义定义定义定义定义定义 设设设设设设f f f( ( (x x x) ) )是区间是区间是区间是区间是区间是区间 a a a, , , b b b 上的一个连上的一个连上的一个连上的一个连上的一个连上的一个连续可微函数续可微函数续可微函数续可微函数续可微函数续可微函数, , , 在区间在区间在区间在区间在区间在区间 a a a, , , b b

6、 b 上给定一组节点上给定一组节点上给定一组节点上给定一组节点上给定一组节点上给定一组节点: : : a a a= = =x x x0 0 0 x x x1 1 1 x x x2 2 2 x x xn n n= = =b b b设函数设函数设函数设函数设函数设函数S S S( ( (x x x) ) )满足条件满足条件满足条件满足条件满足条件满足条件: : :一、样条函数一、样条函数 (1) (1) (1) S S S( ( (x x x) ) )在每个子区间在每个子区间在每个子区间在每个子区间在每个子区间在每个子区间 x x xi i i , , , x x xi i i+1+1+1( (

7、(i i i=0, 1, 2, =0, 1, 2, =0, 1, 2, , , , n n n 1)1)1)上是次数不超过上是次数不超过上是次数不超过上是次数不超过上是次数不超过上是次数不超过mm m的多项式的多项式的多项式的多项式的多项式的多项式; ; ; (2) (2) (2) S S S( ( (x x x) ) )在区间在区间在区间在区间在区间在区间 a a a, , , b b b 上有上有上有上有上有上有mm m 1 1 1阶连续导数阶连续导数阶连续导数阶连续导数阶连续导数阶连续导数. . .则称则称则称则称则称则称S S S( ( (x x x) ) )是定义在是定义在是定义在是

8、定义在是定义在是定义在 a a a, , , b b b 上的上的上的上的上的上的mm m次样条函数次样条函数次样条函数次样条函数次样条函数次样条函数, , , x x x0 0 0, , , x x x1 1 1, , , x x x2 2 2, , , , , , x x xn n n称为称为称为称为称为称为样条节点样条节点样条节点样条节点样条节点样条节点, , , 其中其中其中其中其中其中x x x1 1 1, , , , , , x x xn n n 1 1 1称为称为称为称为称为称为内结点内结点内结点内结点内结点内结点, , , x x x0 0 0, , ,x x xn n n 称

9、为称为称为称为称为称为边界节点边界节点边界节点边界节点边界节点边界节点。 当当当当当当mm m=3=3=3时时时时时时, , , 便成为最常用的便成为最常用的便成为最常用的便成为最常用的便成为最常用的便成为最常用的三次样条函数三次样条函数三次样条函数三次样条函数三次样条函数三次样条函数. . . 样条插值的样条插值的样条插值的样条插值的样条插值的样条插值的思想思想思想思想思想思想: : : 逐段选取适当的逐段选取适当的逐段选取适当的逐段选取适当的逐段选取适当的逐段选取适当的低次多项式低次多项式低次多项式低次多项式低次多项式低次多项式, , , 按一定的光滑性要求连接起来按一定的光滑性要求连接起

10、来按一定的光滑性要求连接起来按一定的光滑性要求连接起来按一定的光滑性要求连接起来按一定的光滑性要求连接起来构成插值函数构成插值函数构成插值函数构成插值函数构成插值函数构成插值函数. . .二、三次样条插值二、三次样条插值 定义定义定义定义定义定义 设给定区间设给定区间设给定区间设给定区间设给定区间设给定区间 a a a, , , b b b 上上上上上上n n n+1+1+1个点个点个点个点个点个点 a a a= = =x x x0 0 0 x x x1 1 1 x x x2 2 2 x x xn n n= = =b b b, , , 以及相应的函数值以及相应的函数值以及相应的函数值以及相应的

11、函数值以及相应的函数值以及相应的函数值 y y yi i i=f=f=f( ( (x x xi i i), ), ), i=i=i=0, 1, , 0, 1, , 0, 1, , n n n. . . 如果如果如果如果如果如果函数函数函数函数函数函数S S S( ( (x x x) ) )满足满足满足满足满足满足: : : (1) (1) (1)在每个子区间在每个子区间在每个子区间在每个子区间在每个子区间在每个子区间 x x xk k k , , , x x xk k k+1+1+1( ( (k k k=0,1,=0,1,=0,1,n n n 1)1)1)上上上上上上, , , S S S(

12、( (x x x) ) )是不超过三次的多项式是不超过三次的多项式是不超过三次的多项式是不超过三次的多项式是不超过三次的多项式是不超过三次的多项式, , , 且且且且且且S S S( ( (x x xi i i ) ) )= = =y y yi i i, , , i=i=i=0, 1, , 0, 1, , 0, 1, , n n n; ; ; (2) (2) (2) S S S( ( (x x x) ) )、 S S S ( ( (x x x) ) )、 S S S ( ( (x x x) ) )在在在在在在 a a a, , , b b b 上连续上连续上连续上连续上连续上连续. . .则称

13、则称则称则称则称则称S S S( ( (x x x) ) )是是是是是是f f f( ( (x x x) ) )在节点在节点在节点在节点在节点在节点x x x0 0 0, , , x x x1 1 1, , , x x x2 2 2, , , , , , x x xn n n上的上的上的上的上的上的三次样条插三次样条插三次样条插三次样条插三次样条插三次样条插 值函数值函数值函数值函数值函数值函数. . . 例例例例例例1 1 1 给定区间给定区间给定区间给定区间给定区间给定区间0, 30, 30, 3上上上上上上 3 3 3 个点的函数个点的函数个点的函数个点的函数个点的函数个点的函数值值值值

14、值值 f f f(0)=0, (0)=0, (0)=0, f f f(1)=2, (1)=2, (1)=2, f f f(3)=4, (3)=4, (3)=4, 试求数试求数试求数试求数试求数试求数 a a a, , , b b b, , , c c c, , , d d d, , , 使函数使函数使函数使函数使函数使函数 S S S( ( (x x x) ) )为给定点上的三次样条插值为给定点上的三次样条插值为给定点上的三次样条插值为给定点上的三次样条插值为给定点上的三次样条插值为给定点上的三次样条插值函数函数函数函数函数函数. . . 其中其中其中其中其中其中232,01( ).1, 13

15、xxdxS xaxbxcxx答案答案答案答案答案答案: : :1,4,2,0.abcd 给定给定给定给定给定给定n n n+1+1+1个样点个样点个样点个样点个样点个样点( ( (x x xi i i, , , y y yi i i )( )( )(i=i=i=0, 1, , 0, 1, , 0, 1, , n n n), ), ),确定一个三次样条插值函数需要确定一个三次样条插值函数需要确定一个三次样条插值函数需要确定一个三次样条插值函数需要确定一个三次样条插值函数需要确定一个三次样条插值函数需要4 4 4n n n个独个独个独个独个独个独立条件立条件立条件立条件立条件立条件. . . 在定

16、义中在定义中在定义中在定义中在定义中在定义中, , , 已指定了已指定了已指定了已指定了已指定了已指定了4 4 4n n n 2 2 2个条件个条件个条件个条件个条件个条件, , , 即即即即即即00(), ()()(),(1,2,.1)()(),()(),nniiiiiiiS xyS xyS xSxyinSxSxSxSx所以所以所以所以所以所以, , , 一般需补充指定一般需补充指定一般需补充指定一般需补充指定一般需补充指定一般需补充指定2 2 2个个个个个个边界条件边界条件边界条件边界条件边界条件边界条件. . .三、三次样条函数的构造三、三次样条函数的构造 三弯矩插值法三弯矩插值法 记记

17、记记记记 MMMi i i = = = S S S ( ( (x x xi i i), ), ), f f f( ( (x x xi i i)=)=)= f f fi i i= = = y y yi i i , , , 考虑它在任考虑它在任考虑它在任考虑它在任考虑它在任考虑它在任一一一一一一区间区间区间区间区间区间 x x xi i i, , , x x xi i i+1+1+1 上的形式上的形式上的形式上的形式上的形式上的形式. . . 根据三次样条的定义可知根据三次样条的定义可知根据三次样条的定义可知根据三次样条的定义可知根据三次样条的定义可知根据三次样条的定义可知, , , S S S(

18、 ( (x x x) ) )的二阶导数的二阶导数的二阶导数的二阶导数的二阶导数的二阶导数 S S S ( ( (x x x) ) )在每一个子区间在每一个子区间在每一个子区间在每一个子区间在每一个子区间在每一个子区间 x x xi i i, , , x x xi i i+1+1+1( ( ( i i i=0, 1, =0, 1, =0, 1, 2, 2, 2, , , , n n n 1)1)1)上都是线性函数上都是线性函数上都是线性函数上都是线性函数上都是线性函数上都是线性函数. . . 于是在于是在于是在于是在于是在于是在 x x xi i i, , , x x xi i i+1+1+1

19、上上上上上上 S S S( ( (x x x)=)=)=S S Si i i( ( (x x x) ) )的二阶导数表示成的二阶导数表示成的二阶导数表示成的二阶导数表示成的二阶导数表示成的二阶导数表示成111( ) , ,iiiiiiiixxxxSxMMxxxhh其中其中其中其中其中其中 h h hi i i= = = x x xi i i+1+1+1 x x xi i i . . . 对对对对对对S S S ( ( (x x x) ) )连续积分两次连续积分两次连续积分两次连续积分两次连续积分两次连续积分两次, , , 并利用插值并利用插值并利用插值并利用插值并利用插值并利用插值条件条件条件

20、条件条件条件S S S( ( (x x xi i i)= )= )= y y yi i i , , , 得到得到得到得到得到得到3311111()()( )66 ()()66iiiiiiiiiiiiiiiixxxxS xMMhhyMyMhxxhxxhh1,iixx x只要能求出所有的只要能求出所有的只要能求出所有的只要能求出所有的只要能求出所有的只要能求出所有的 MMMi i i, , , 就能求出三次就能求出三次就能求出三次就能求出三次就能求出三次就能求出三次样条插值函数样条插值函数样条插值函数样条插值函数样条插值函数样条插值函数S S S( ( (x x x). ). ).221111()

21、()( )22 6iiiiiiiiiiiixxxxS xMMhhyyMMhh 下面考虑下面考虑下面考虑下面考虑下面考虑下面考虑 MMMi i i 的求法的求法的求法的求法的求法的求法. . . 由连续性由连续性由连续性由连续性由连续性由连续性 S S S ( ( (x x xi i i )= )= )= S S S ( ( (x x xi i i+ + +) ) ), (, (, (i i i=1, 2, , =1, 2, , =1, 2, , n n n 1) 1) 1) 得得得得得得 i i iMMMi i i 1 1 1+2+2+2MMMi i i+ + + i i iMMMi i i+

22、1+1+1= = = d d di i i 1111111,6()iiiiiiiiiiiiiiihhhyyyydhhhh 其中其中其中其中其中其中该方程组有该方程组有该方程组有该方程组有该方程组有该方程组有n n n 1 1 1个方程个方程个方程个方程个方程个方程, , , 但有但有但有但有但有但有n n n+1+1+1个变量个变量个变量个变量个变量个变量MMMi i i. . . 下面介绍几种常用的边界条件下面介绍几种常用的边界条件下面介绍几种常用的边界条件下面介绍几种常用的边界条件下面介绍几种常用的边界条件下面介绍几种常用的边界条件 第第第第第第1 1 1型边界条件型边界条件型边界条件型边

23、界条件型边界条件型边界条件: : : 已知已知已知已知已知已知f f f( ( (x x x) ) )在两端点的导数在两端点的导数在两端点的导数在两端点的导数在两端点的导数在两端点的导数f f f ( ( (a a a) ) )和和和和和和f f f ( ( (b b b) ) ), , , 要求要求要求要求要求要求 S S S ( ( (a a a) = ) = ) = f f f ( ( (a a a), ), ), S S S ( ( (b b b) = ) = ) = f f f ( ( (b b b) ) ) 第第第第第第2 2 2型边界条件型边界条件型边界条件型边界条件型边界条件型

24、边界条件: : : 已知已知已知已知已知已知f f f( ( (x x x) ) )在两端点的二阶导数在两端点的二阶导数在两端点的二阶导数在两端点的二阶导数在两端点的二阶导数在两端点的二阶导数f f f ( ( (a a a) ) )和和和和和和f f f ( ( (b b b) ) ) , , ,要求要求要求要求要求要求 S S S ( ( (a a a)=)=)=MMM0 0 0 = = = f f f ( ( (a a a), ), ), S S S ( ( (b b b)=)=)=MMMn n n= = = f f f ( ( (b b b) ) )特别当特别当特别当特别当特别当特别当

25、 S S S ( ( (a a a)= )= )= S S S ( ( (b b b) =0) =0) =0时时时时时时, , , S S S( ( (x x x) ) )称为称为称为称为称为称为自然三次样条自然三次样条自然三次样条自然三次样条自然三次样条自然三次样条. . . 第第第第第第3 3 3型边界条件型边界条件型边界条件型边界条件型边界条件型边界条件： 已知已知已知已知已知已知f f f( ( (x x x) ) )是以是以是以是以是以是以b b b a a a为周期的周期函数为周期的周期函数为周期的周期函数为周期的周期函数为周期的周期函数为周期的周期函数, , , 要求要求要求要求

26、要求要求S S S( ( (x x x) ) )满满满满满满 足周期条件足周期条件足周期条件足周期条件足周期条件足周期条件 S S S( ( (a a a) = ) = ) = S S S( ( (b b b), ), ), S S S ( ( (a a a+ + +)= )= )= S S S ( ( (b b b ), ), ), S S S ( ( (a a a+ + +)= )= )= S S S ( ( (b b b ) ) ) 三次样条插值问题加上第三次样条插值问题加上第三次样条插值问题加上第三次样条插值问题加上第三次样条插值问题加上第三次样条插值问题加上第 i i i 型边界型边

27、界型边界型边界型边界型边界条件称为第条件称为第条件称为第条件称为第条件称为第条件称为第 i i i型插值问题型插值问题型插值问题型插值问题型插值问题型插值问题( ( (i i i1, 2, 3). 1, 2, 3). 1, 2, 3). 可可可可可可以证明第以证明第以证明第以证明第以证明第以证明第 i i i 型插值问题的解是存在且唯一的型插值问题的解是存在且唯一的型插值问题的解是存在且唯一的型插值问题的解是存在且唯一的型插值问题的解是存在且唯一的型插值问题的解是存在且唯一的. . .他们对应如下的方程组他们对应如下的方程组他们对应如下的方程组他们对应如下的方程组他们对应如下的方程组他们对应如

28、下的方程组: : :000111122221111220222 nnnnnnnMdMdMdMdMd对于第对于第对于第对于第对于第对于第1 1 1型插值问题型插值问题型插值问题型插值问题型插值问题型插值问题: : : 00100001111, 6 (),1, 6(). nnnnnnndyyhyhdyyyhh对于第对于第对于第对于第对于第对于第2 2 2型插值问题型插值问题型插值问题型插值问题型插值问题型插值问题: : : 0000, 2, 0, 2. nnndydy对于第对于第对于第对于第对于第对于第3 3 3型插值问题型插值问题型插值问题型插值问题型插值问题型插值问题: : : 0112 nn

29、nnnnMMMMMd1110111(), 6()() 其其中中nnnnnnnnnnhhhdyyhyyhhh 以上各组条件与前述方程组联立以上各组条件与前述方程组联立以上各组条件与前述方程组联立以上各组条件与前述方程组联立以上各组条件与前述方程组联立以上各组条件与前述方程组联立, , ,可以解出未知参数可以解出未知参数可以解出未知参数可以解出未知参数可以解出未知参数可以解出未知参数 MMM0 0 0, , , MMM1 1 1, , , , , , MMMn n n, , , 然然然然然然后代入后代入后代入后代入后代入后代入S S S( ( (x x x) ) ) 表达式表达式表达式表达式表达式

30、表达式, , , 即可求得样条函数即可求得样条函数即可求得样条函数即可求得样条函数即可求得样条函数即可求得样条函数. . . 上面构造方法中上面构造方法中上面构造方法中上面构造方法中上面构造方法中上面构造方法中 MMMi i i 相应于力学中细梁在相应于力学中细梁在相应于力学中细梁在相应于力学中细梁在相应于力学中细梁在相应于力学中细梁在 x x xi i i 处处处处处处截面的弯矩截面的弯矩截面的弯矩截面的弯矩截面的弯矩截面的弯矩, , , 每一个方程中又至多出现相邻的三个每一个方程中又至多出现相邻的三个每一个方程中又至多出现相邻的三个每一个方程中又至多出现相邻的三个每一个方程中又至多出现相邻

31、的三个每一个方程中又至多出现相邻的三个MMMi i i, , , 通常称为通常称为通常称为通常称为通常称为通常称为三弯矩法三弯矩法三弯矩法三弯矩法三弯矩法三弯矩法. . . 求三次样条插值函数的步骤归纳为求三次样条插值函数的步骤归纳为求三次样条插值函数的步骤归纳为求三次样条插值函数的步骤归纳为求三次样条插值函数的步骤归纳为求三次样条插值函数的步骤归纳为: : : (1) (1) (1)确定边界条件确定边界条件确定边界条件确定边界条件确定边界条件确定边界条件, , , 判定是第几型插判定是第几型插判定是第几型插判定是第几型插判定是第几型插判定是第几型插值问题值问题值问题值问题值问题值问题. .

32、. (2) (2) (2)根据所确定的条件计算各值根据所确定的条件计算各值根据所确定的条件计算各值根据所确定的条件计算各值根据所确定的条件计算各值根据所确定的条件计算各值, , , 形成方程组形成方程组形成方程组形成方程组形成方程组形成方程组. . . (3) (3) (3)解方程组解方程组解方程组解方程组解方程组解方程组, , , 求得求得求得求得求得求得MMM0 0 0, , , MMM1 1 1 , , , MMM2 2 2, , , MMMn n n. . . (4) (4) (4)将求得的将求得的将求得的将求得的将求得的将求得的 MMMi i i 值代回值代回值代回值代回值代回值代回 S S S( ( (x x x) ) )的表达式中的表达式中的表达式中的表达式中的表达式中的表达式中, , , 从而可从而可从而可从而可从而可从而可求得函数求得函数求得函数求得函数求得函数求得函数 y y y= = =f f f( ( (x x x) ) )在任一点的近似值在任一点的近似值在任一点的近似值在任一点的近似值在任一点的近似值在任一点的近似值S S S( ( (x x x) ) ). . . 例例例例例例2 2 2 给定函数表给定函数表给定函数表给定函数表给定函数表给定函数