三次样条插值课件.ppt

1、第二章 插值法引例引例2.7.1 三次样条插值函数的概念三次样条插值函数的概念一一 背景背景二、样条函数的定义二、样条函数的定义 例例2.13 定理定理2.8（3 次样条插值函数存在唯一）次样条插值函数存在唯一）2.7.2 三弯矩法三弯矩法边界条件边界条件1（固支边界）（固支边界）边界条件边界条件2（简支边界）（简支边界）边界条件边界条件3（周期边界）（周期边界）例例2.14,2.15 2.7.3 m关系式关系式2.7.4 三次样条插值函数的性质三次样条插值函数的性质2.7 三次样条插值三次样条插值第二章 插值法引例引例:y=sin x 在区间在区间0，上的插值逼近上的插值逼近 1.1.二次插

2、值二次插值 0123400.20.40.60.812.两点埃尔米特插值两点埃尔米特插值 0123400.20.40.60.813.分段埃尔米特插值分段埃尔米特插值x 0 /2 Sin x 010Cos x 101x 0 Sin x 00Cos x 11 第二章 插值法高次插值出现龙格现象高次插值出现龙格现象L-插值（牛顿插值（牛顿插值）插值）Hermite插值插值 分段分段插值插值但分段线性插值在节点处不一定光滑但分段线性插值在节点处不一定光滑分段分段Hermite插值插值但但导数值导数值不容易提取（找到）不容易提取（找到）为得到光滑度更高、应用方便的插值函数为得到光滑度更高、应用方便的插值函

3、数,我们引入我们引入样条样条插值函数插值函数。“样条样条”名词来源于工程中船体、汽车、飞机等名词来源于工程中船体、汽车、飞机等的外形设计：给出外形曲线上的一组离散点的外形设计：给出外形曲线上的一组离散点(样点样点)，如，如(xi,yi)，i=0,1,2,n,将有弹性的将有弹性的细长木条细长木条或或钢条钢条(样条样条)在样点上在样点上固定，使其在其它地方自由弯曲，这样样条所表示的曲线，固定，使其在其它地方自由弯曲，这样样条所表示的曲线，称为称为样条曲线样条曲线(函数函数)。一一 背景背景2.7.1 三次样条插值函数的概念三次样条插值函数的概念第二章 插值法-50500.51x=-5:5;y=1.

4、/(1+x.2);plot(x,y,x,y,o)-50500.51x=-5:5;y=1./(1+x.2);xi=-5:.05:5;yi=spline(x,y,xi);plot(xi,yi,b,x,y,ro)被插值函数被插值函数:211)(xxf -5 x 53/18第二章 插值法x=0,0.0155,0.1485,0.3493,0.6480,1.0547,2.0;y=0,0.1242,0.3654,0.4975,0.5472,0.4781,0;n=length(x);t=0:n-1;tt=0:.25:n-1;xx=spline(t,x,tt);yy=spline(t,y,tt);plot(xx

5、,yy,x,y,o)第二章 插值法相同数据相同数据3 3次样条插值与次样条插值与Lagrange插值效果比较插值效果比较Cubic Spline Interpolation Lagrange Interpolation第二章 插值法 下面介绍应用最广且只有二阶连续导数的三次样条函下面介绍应用最广且只有二阶连续导数的三次样条函数数 在数学上，在数学上，三次样条曲线三次样条曲线表现为近似于一条分段的三次表现为近似于一条分段的三次多项式，它要求在节点处具有一阶和二阶连续导数。多项式，它要求在节点处具有一阶和二阶连续导数。二、样条函数的定义二、样条函数的定义 定义定义 2.8 （三三次样条函数）次样条

6、函数）)(xSb在每一个小区间在每一个小区间1,jjxx上上是次数是次数 1,1,0 nj3 多项式。多项式。baCxSa,)()(2，即具有连续的一阶，二阶导数。，即具有连续的一阶，二阶导数。满足下述条件：满足下述条件：,:10bxxxan )(xS如果函数如果函数 设有对设有对a,b的剖分的剖分的一个的一个3次样条函数。次样条函数。)(xS为关于剖分为关于剖分 则称则称第二章 插值法定义定义2.8*给定区间给定区间a,b上的一个分划上的一个分划:a=x0 x1 xn=b已知已知 f(xj)=yj (j=0,1,n),如果如果 ,),(,),(,),()(1212101nnnxxxxSxxx

7、xSxxxxSxS满足满足:(1)S(x)在在 xj，xj+1上为三次多项式上为三次多项式;(2)S”(x)在区间在区间a，b上连续上连续;(3)S(xj)=yj (j=0,1,n).则称则称 S(x)为三次样条插值函数为三次样条插值函数.第二章 插值法注：注：三次样条与分段三次样条与分段 Hermite 插值的根本区别在于插值的根本区别在于S(x)自自身光滑身光滑，不需要知道，不需要知道 f 的导数值（除了在的导数值（除了在2个端点可能需个端点可能需要）；而要）；而Hermite插值依赖于插值依赖于f 在所有插值点的导数值。在所有插值点的导数值。f(x)H(x)S(x)第二章 插值法 132

8、1,(0,1,1),4iiiiiiiiiiiiiiS xSxix xSxa xb xc xdxx xina b c dnS xn首先指出单靠一张函数表是不能完全确定一个三次样条插值函数的。这是因为：由条件(1)，三次样条插值函数是一个分段三次多项式，若用表示它在第 个区间上的表达式，则这里有四个待定系数,子区间共有 个.要确定需要确定个待定系数。个方程。可得待定系数应满足的可，由条件上连续即个子区间的连接点上连续，只要它们在各间在整个插值区及其导数次多项式另一方面，要求分段三243,2,121nxxxbaxSxSxSn第二章 插值法插值条件插值条件:S(xj)=yj (j=0,1,n)n+1个

9、个连续性条件连续性条件:S(xj+0)=S(xj-0)(j=1,n-1)S(xj+0)=S(xj-0)(j=1,n-1)S (xj+0)=S (xj-0)(j=1,n-1)3(n-1)个个共可建立方程共可建立方程(4n-2)个！个！方程数少于未知数个数方程数少于未知数个数？32000001321111123213211111,(),iiiiiinnnnnna xb xc xdxx xa xb xc xdxx xS xa xb xc xdxx xaxbxcxdxxx第二章 插值法 共有共有24 n个条件个条件,要唯一确定要唯一确定 ,还必须附加还必须附加2 2个条件个条件)(xS这两个条件常在插

10、值区间这两个条件常在插值区间a,ba,b 的边界点的边界点a,ba,b处给出，称处给出，称为为边界条件边界条件。边界条件的类型很多，常见的有：。边界条件的类型很多，常见的有：附加附加2个条件，个条件，有多种给法有多种给法.最常见的给法是最常见的给法是:(a)固支边界固支边界(b)简支边界简支边界 特别地特别地,(自然边界自然边界,三次自然样条三次自然样条););000,nnnSxfxMSxfxM00,nMM000,nnnSxfxm Sxfxm(1)(1)(2)(2)注：注：一般不取一端是一阶导数而另一般不取一端是一阶导数而另一端是二阶导数一端是二阶导数。第二章 插值法)(c第第3种边界条件（周

11、期边界条件）：种边界条件（周期边界条件）：注意：上述注意：上述给出的给出的 个条件是问题本身隐含的，个条件是问题本身隐含的，和和共共 个独立条件须提供，故个独立条件须提供，故 节节点三次样插值点三次样插值问题只有问题只有 个自由度个自由度.(.(请与分段三次请与分段三次HermiteHermite插值比较插值比较!)!)3n 33n 1n 3n 0000000000()nnnnyf xbaSxSxbaSxSxSxSxf xyyS xS xS x 当是周期为的函数时，则要求及其导数均是以为周期的函数。扩充边界条件为：（由周期性知，从而必有故不必再提此要求）-此时称为周期样条函数 第二章 插值法)

12、,1,0(),(,(nixfxii 且且;10bxxxan (1)如果如果是定义在是定义在上函数且已知上函数且已知)(xfy 函数表函数表)(xf,ba 定理定理2.8（3 次样条插值函数存在唯一次样条插值函数存在唯一）唯一唯一3 3次样条插值函数次样条插值函数)(xS,且满足且满足）。）或或（或或（cba)()(xf,ba (2)给定边界条件给定边界条件）或或（或或（cba)(，则，则于于存在存在第二章 插值法例例 2.13 已知 f(1)=1,f(0)=0,f(1)=1.求1，1 上的三次自然样条(满足自然边界条件).解解 设,)(1001222232112131xdxcxbxaxdxcx

13、bxaxS则有:S(-1)=a1+b1c1+d1=f(-1)=1,S(0)=d1=f(0)=0,S(1)=a2+b2+c2+d2=f(1)=1,S(0-0)=d1=S(0+0)=d2,S-(0)=c1=S+(0)=c2,S-(0)=b1=S+(0)=b2 由自然边界条件:S(0)=6a1+2b1=0,S(1)=6a2+2b2=0 解方程组,得 a1=-a2=1/2,b1=b2=3/2,c1=c2=d1=d2=0第二章 插值法,)(1023210123212323xxxxxxxS问题的解-1-0.8-0.6-0.4-0.200.20.40.60.8100.20.40.60.81x=-1,0,1;

14、y=1,0,1;f1=inline(0.5*x.3+1.5*x.2);f2=inline(-0.5*x.3+1.5*x.2);t1=-1:.1:0;t2=0:.1:1;p1=f1(t1);p2=f2(t2);plot(x,y,o,t1,t2,p1,p2,r)Hold on,plot(t1,t2,t1,t2.2)y=x2第二章 插值法三次样条插值函数三次样条插值函数 可以有多种表达式，可以有多种表达式，有时用二阶导数值有时用二阶导数值表示时，使用更方便。表示时，使用更方便。在力学上解释为细梁在力学上解释为细梁在在 处的弯矩，并且得到的弯矩与相邻两个弯矩处的弯矩，并且得到的弯矩与相邻两个弯矩有关，

15、故称用有关，故称用 表示表示 的算法为的算法为三弯矩算法三弯矩算法。)(xS()(0,1,)iiSxMinMixiMi)(xS2.7.2 构造三次样条插值函数的构造三次样条插值函数的三弯矩法三弯矩法 -三次样条插值函数的二阶导数表示三次样条插值函数的二阶导数表示第二章 插值法是是三三次次样样条条因因为为)(xSj11,()(,),jjjjjjxx xhS xSxxx令1(,)jjjxSxx所以是在上一次函数,插插值值函函数数，),2,0,1j)(nMxSjj ，（，（令令由两点拉格朗日插值由两点拉格朗日插值可表示为可表示为,)(11 jjjjjjMhxxMhxxxS参数参数(2.46)对对上上

16、式积分式积分,得得22111()()(),22jjjjjjxxxxS xMMchh(2.47)(2.48)再积分再积分,得得331112()()(),66jjjjjjxxxxS xMMc xchh,1 jjxxxjjjjjxxhxxx 11,第二章 插值法 由条件由条件11)(,)(jjjjyxSyxS，确定积分常数，确定积分常数12,c c(2.49)(2.47)(2.48)21221111211(),61(),6jjjjjjjjjjS xh Mc xcyS xh Mc xcy22111()()(),22jjjjjjxxxxS xMMchh 331112()()(),66jjjjjjxxxx

17、S xMMc xchh,1 jjxxx111112111(),61().6jjjjjjjjjjjjjjjjyych MMhy xy xch x Mx Mh第二章 插值法 将将上式上式代入代入(2.48)得到得到三三次样条插值函数的表达式次样条插值函数的表达式331122111()()()66()(),66jjjjjjjjjjjjjjjjxxxxS xMMhhM hxxMhxxyyhh(2.50)由上讨论可知由上讨论可知,只要确定只要确定Mj(j=0,1,n)这这n+1个值个值,就就可定出三次样条插值函数可定出三次样条插值函数S(x)。为了确定为了确定Mj(j=0,1,n),对对S(x)求导得求

18、导得221111()()()226jjjjjjjjjjjjxxxxyyMMS xMMhhhh 1,jjxxx(2.51)1,jjxxx第二章 插值法13311112211111112211111111()()()()66()(),66()()()22 6,jjjjjjjjjjjjjjjjjjjjjjjjjjjjjjS xxxxxS xMMhhMhxxM hxxyyhhxxxxS xMMhxxhyyMMhh 类似地可求出在区间上的表达式,从而得1,jjxxx(2.52)1,jjxxx第二章 插值法1111111()(0)(0),636 1,1,jjjjjjjjjjjjjjjS xS xS xhh

19、hhyyyyMMMhhjn利用在内接点的连续性,即可得 0,:(0)(0)njjMMS xS x为了求要用导数连续条件1111111(2.51):(0),36(2.52):(0),63jjjjjjjjjjjjjjjjhhyyS xMMhhhyyS xMMh由得由得(2.53)(2.54)(2.55)第二章 插值法1111111,636jjjjjjjjjjjjjhhhhyyyyMMMhh（2.55）(1,2,1)jn上式两边同乘以上式两边同乘以 ,即得方程即得方程 16jjhh11111111162jjjjjjjjjjjjjjjjjhhyyyyMMMhhhhhhhh11111111166,.jj

20、jjjjjjjjjjjjjjjjjjhhhhhhyyyydf xxxhhhh若记若记 （2.56）第二章 插值法1111111 2,1,1,6,.jjjjjjjjjjjjjjjjjjMMMdjnhhdf xx xhhhh其中所得方程可简写成所得方程可简写成10112121223212111222nnnnnnMMMdMMMdMMMd（2.58）即即 （2.57）个个方方程程1 n 三弯矩方程三弯矩方程第二章 插值法 这是一个含有这是一个含有n+1+1个未知数、个未知数、n-1-1个方程的线性方个方程的线性方程组程组.要完全确定要完全确定Mi(i=0,1,n)的值还需要补充两个的值还需要补充两个条

21、件条件,这两个条件通常根据实际问题的需要，根据插这两个条件通常根据实际问题的需要，根据插值区间值区间 a,b 的两个端点处的边界条件来补充。的两个端点处的边界条件来补充。第二章 插值法由由(2.53),得得0001100063fhyyMhMh 由由(2.54),得得1111136 nnnnnnnnhyyfMhMh(1)若若已知，已知，,)(,)(000nnnmfxSmfxS 11()36jjjjjjjjyyhhSxMMh11111()36jjjjjjjjyyhhS xMMh )(621111 nnnnnnnhyyfhMM)(620001010fhyyhMM 0d nd 则令则令j=0，令令j=

22、n，边界条件边界条件1（固支边界）（固支边界）-第二章 插值法001,1,nnjjM令得满足的方程组0011111111212212nnnnnnMdMdMdMd(2.59)对角占优的三对角带状矩阵第二章 插值法(2)若若nnnfMxSfMxS )(,)(000已知，已知，代入方程代入方程(2.58)，只只需解需解n-1个方程个方程 nnnnnnnnnfdddfdMMMM112201112211222212222 (2.60)边界条件边界条件2（简支边界）（简支边界）-对角占优的三对角带状矩阵第二章 插值法 (3)对第三类边界条件：对第三类边界条件：),0()0(0 nxSxS)0()0(0 n

23、xSxS 0M,nM 11001101110)(3166 nnnnnnnhyyhyyMhhMhMh两边同除以两边同除以得得,610 nhh)(62110011011011100 nnnnnnnnnhyyhyyhhMMhhhMhhh(j=n)1111111 (0),(2.53)36 (0),(2.54)63jjjjjjjjjjjjjjjjhhyyS xMMhhhyyS xMMh由和可得(j=n)(j=0)边界条件边界条件3（周期边界）（周期边界）-第二章 插值法)(62110011011011100 nnnnnnnnnhyyhyyhhMMhhhMhhh)(61100110 nnnnnhyyhyy

24、hhd,1101 nnnnhhh ,100 nnhhh 令令nnnnndMMM 211 得得又由又由 nMM 0,三弯矩方程可写为三弯矩方程可写为 nnnnnnnnddddMMMM1211211122112222 (2.61)第二章 插值法 nnnnnnnnddddMMMM1211211122112222 (2.61)，nnnnnnnddddMMMM110110111102222 (2.59)nnnnnnnnnfdddfdMMMM11220112211222212222 (2.60)小结：在三个边界条件下的三弯矩方程小结：在三个边界条件下的三弯矩方程第二章 插值法说明：说明：（1）方程组方程组

25、(2.59)(2.61)系数矩阵都是严格对角占优矩系数矩阵都是严格对角占优矩 阵，因此方程组阵，因此方程组(2.59)(2.61)有唯一解有唯一解 （2）Mj 在力学上为细梁在在力学上为细梁在xj处处截面截面处的处的弯矩弯矩,且弯矩与且弯矩与相邻相邻的两个弯矩有关的两个弯矩有关,故方程组故方程组(2.59)(2.61)称为称为三弯矩三弯矩方程。方程。Mj 在数学上称为在数学上称为曲率曲率。实际上实际上,方程组方程组(2.59)(2.61)的系数矩阵是一类特殊的系数矩阵是一类特殊的矩阵，在后面线性方程组的解法中，将专门介绍这类的矩阵，在后面线性方程组的解法中，将专门介绍这类方程组的解法和性质。方

26、程组的解法和性质。第二章 插值法.yfx例2.14：已知函数的函数值如下：142.75.05.125.1SSxS，使它满足边界条件上求三次样条插值函数，在区间 x-1.5 0 1 2 y0.125 -1 1 9,10,1,2,32.59iiiiidMMi 解：先根据给定条件和边界条件算出，写出确定的线性方程组，在本例中给出的是第种边界条件，确定的线性方程组形如。第二章 插值法8,.2,.75.0,1,1,5.1,3221103211xxfxxfxxfhhhxxhiii由所给定的函数表知：121212,1,1(2.56)0.6.0.5.0.4.0.5.6.6.18iiig indd 于是由和的算

27、式知：0001013323316,66,36ndddfxxyhdyfxxh 由第边界条件下与的计算公式知：第二章 插值法 16.4.4.53210MMMMMxxSii：上的值在各个节点解方程组得 的表达式为：在上知：由区间上的表达式在各个子式即得代入最后将10,3,1,3xxxSnixSxSMii36186.6621005.025.0004.026.00012,32103210MMMMMMMM的方程组为：故确定第二章 插值法 1021111121003101311016666hxxhMyhxxhMyhxxMhxxMxS 2,1,3642,1,0.1223322xxxxxSxxxS 213642

28、101205.11223223xxxxxxxxxxS在本例中，将代入整理后可得：故所求三次样条插值函数为：0,5.1,124,5.1.125.0.0.5.1231101010 xxxxSMMyyxx第二章 插值法例2.15 已知的函数值如下：x 1 2 4 5 f(x)1 3 4 2在区间1,5上求三次样条插值函数S(x),使它满足边界条件 0)5(,0)1(SS解:这是在第二种边界条件下的插值问题，故确定 的方程组形如（2.60）所示，由已知边界条件,有 则得求解 的方程组为 3210,MMMM0)(,0)(333000 MyxSMyxS21,MM11122222MdMd nnnnnnnnn

29、fdddfdMMMM11220112211222212222 第二章 插值法根据给定数据和边界条件算出 与 21,12,dd121321hhh2,21,2,322110 xxfxxfxxf32,3232222121hhhhhh112010166 1(,)(2)33 2df x xf x xhh 2231212661(,)(2)532df x xf x xhh 第二章 插值法523233222121MMMM则得方程组 解得 49,4321MM又 030 MM即得S(x)在各子区间上的表达式,由式（2.51）知,S(x)在 上的表达式为代入式(2.50)3,2,1)(ixSi10,xx130113

30、1016)(6)()(hxxMhxxMxS102111112100)(6)(6hxxhMyhxxhMy将 代入上式化简后得 43,0,1,3,1,2,11011010MMhyyxx第二章 插值法1478381)(231xxxxS同理S(x)在 上的表达式为 21,xx1478381)(232xxxxSS(x)在 上的表达式为 32,xx1949184583)(233xxxxS第二章 插值法故所求的三次样条插值函数S(x)在区间 上的表达式为 5,1)54(1949184583)42(1478381)21(1478381)(232323xxxxxxxxxxxxxS第二章 插值法练习练习 设在节点

31、设在节点 上，函数上，函数 的值为的值为 ，。试求三试求三次样条插值函数次样条插值函数 ，满足条件，满足条件)3,2,1,0(iixi)(xf5.0)(,0)(10 xxff5.1)(,2)(32 xxff)(xS.3.3)(,3.0)()2(,1)(,2.0)()1(3030 xxxxSSSS11121211211110.5(1,2)0.53,66jjjjjjjjjjjjjjjjjhhhhjhhddyyyydhhhh 解解 （1）是固支边界，先求，再求解，可知）是固支边界，先求，再求解，可知,jjjd jM第二章 插值法对第一类边界条件对第一类边界条件03()0.2,()1,S xS x 1

32、000003233226()1.86()3yydfhhyydfhh(2.59),代入有 3638.1215.025.05.025.0123210MMMM0123:0.36,2.52,3.72,0.36MMMM 解得代入三次样条插值函数的表达式（代入三次样条插值函数的表达式（2.50），经化简有），经化简有0011122233dddd2222第二章 插值法3232320.480.180.2,()1.04(1)1.26(1)1.28(1)0.5,0.68(2)1.86(2)0.68(2)2,xxxS xxxxxxx 0,1 1,22,3xxx（）是简支条件，不过要注意（）是简支条件，不过要注意 的

33、不同。由于的不同。由于 和和 已已知，故可以化简得知，故可以化简得dd3030,M3M0331122111()()()66()(),66jjjjjjjjjjjjjjjjxxxxS xMMhhM hxxMhxxyyhh代入三次样条插值函数的表达式（代入三次样条插值函数的表达式（2.50），经化简有），经化简有第二章 插值法.3.153.6411421 MM由此解得由此解得 。122.7,4.5MM 将将 代入三次样条插值函数的表达式（代入三次样条插值函数的表达式（2.50），经），经化简有化简有MMMM3210,2)2(45.0225.213.1,5.0)1(35.1135.112.1,15.0

34、15.05.0)(232323xxxxxxxxxxS 3,2 2,1 1,0 xxx331122111()()()66()(),66jjjjjjjjjjjjjjjjxxxxS xMMhhM hxxMhxxyyhh第二章 插值法练习练习 已知离散点:(1.1,0.4000),(1.2,0.8000),(1.4,1.6500),(1.5,1.8000),取自然边界条件 M0=Mn=0,构造三次样条插值函数，并计算 f(1.25).解解 n=3.h0=x1-x0=0.1,h1=0.2,h2=0.1,因此，分段的三次样条插值函数为,由（2.50）由(2.56)计算得1112220.6667,0.333

35、3,5,0.3333,0.6667,55.dd 0311220,from(2.60)yields13.125,20.66675,solve31.875.0.6667255MMMMMM 320321322()21.87572.187583.187532.875,1.1,1.2()()37.5141.625173.312569.725,1.2,1.4()53.125239.0625359.5625178.95,1.4,1.5SxxxxxS xS xxxxxSxxxxx 1(1.25)(1.25)(1.25)1.0336,1.251.2,1.4.fSS第二章 插值法 上述三次样条插值的上述三次样条插

36、值的基本思想和特点是基本思想和特点是：先利用一阶导数先利用一阶导数 在内节点在内节点 上的连续性上的连续性以及边界条件以及边界条件,列出确定二阶导数列出确定二阶导数 的线性方的线性方程组（力学上称为程组（力学上称为三弯矩方组三弯矩方组），由此解出），由此解出 ，再用，再用 来表达来表达S(x)。,0,1,iiMSxin1,2,1ix in SxiMiM实际上，还可以通过别的途径来求取三次样条插值函数。实际上，还可以通过别的途径来求取三次样条插值函数。iimm如：可以先利用二阶导数在内节点上的连续性及边界条件，列如：可以先利用二阶导数在内节点上的连续性及边界条件，列出确定一阶导数出确定一阶导数

37、的线性方程组（力的线性方程组（力学上称为三转角方程组），由此解出学上称为三转角方程组），由此解出 ，再用，再用 表达表达 S(x),在在某些情况下，这种方法比前者更简单适用。某些情况下，这种方法比前者更简单适用。,0,1,iimSxin第二章 插值法2.7.3 mm关系式关系式 用一阶导数表示的样条插值函数 给定插值点(xi,yi),设S(xi)=mi,i=0,1,2,n,则 xi,xi+1上的三次Hermite插值为 2211111112211111()1212()()(2.62)iiiiiiiiiiiiiiiiiiiiiiiixxxxxxxxS xyyxxxxxxxxxxxxxxmxxmx

38、xxx 令 hi=xi+1-xi,S(x)C2a,b,对(2.62)求二阶导数1123231122612612()()()2626()().iiiiiiiiiiiiiiiiSxxxyxxyhhhhxxmxxmhhhh第二章 插值法 令 xi+=xi+0,在 xi,xi+1上得到 xi 点的右导数，同理，在 xi-1,xi 上构造三次样条插值 S(x),在 xi-1,xi上得点 xi 的左导数，11226642()(2.63)iiiiiiiiiSxyymmhhhh 112211116624()(2.64)iiiiiiiiiSxyymmhhhh2(),from(0)(+0),we obtain b

39、y(2.63),(2.64),iiS xC a bS xS x第二章 插值法11+2+=,1,2,1(2.65)iiiiiimmmein 三种边界条件：三种边界条件：1111where ,1,(2.66)3()()iiiiiiiiiiiiiiihhheyyyyhh 00(1)Given(),().nnS xm S xm000(2)Given(),().Especially,0.nnnSxMSxMMM00000(3)If ()is periodic function,i.e.()(),then()is also periodic.Thus we have()(),()(),that is,.nn

40、nnnf xf xf xS xS xS xSxSxmmMM第二章 插值法111012222223311112222nnnnnnnnnfemmememfe由此可解得由此可解得m1,m2,mn-1，从而得，从而得 S(x)的表达式的表达式.（2.66）对于边界条件对于边界条件(1),两个方程两个方程则则m1,m2,mn-1满足方程组满足方程组 00,nnmf mf第二章 插值法 对于边界条件对于边界条件(2),可导出两个方程可导出两个方程:0001102000112111001010111426(0)()246(0)()3,23,.222nnnnnnnnnnnnnnnMSxmmyyhhhMSxmm

41、yyhhhffhymmx xhymmxx(2.67)第二章 插值法若令若令fhxxefhxxennnnnff 2,3,2,31100100(0,1,)imin则（则（2.65）和（）和（2.67）可合并成矩阵形式）可合并成矩阵形式 eeeemmmmnnnnnn1101101111212212 （2.68）可解出可解出从而得从而得 S(x)的表达式的表达式.第二章 插值法由由(2.65)和和(2.6)可解出可解出 ，方程组的矩阵形式为，方程组的矩阵形式为),1,0(nimi 对于边界条件（对于边界条件（3），可得），可得emmmmmnnnnnn2110 （2.69）).,(3,110010101

42、xxxxehhhhhhnnnnnnnnnnff 其中其中eeeemmmmnnnnnnnn1211212122112222（2.70）第二章 插值法在实际应用中，如果不需要规定内节点处的一阶导数值，那么使在实际应用中，如果不需要规定内节点处的一阶导数值，那么使用三次样条插值函数会得到很好的效果。三次样条插值函数用三次样条插值函数会得到很好的效果。三次样条插值函数 不不仅在内节点处的二阶导数是连续的，而且仅在内节点处的二阶导数是连续的，而且 逼近逼近 具有很好的具有很好的收敛性，也是数值稳定的。下面给出三次样条插值函数的一些重要性收敛性，也是数值稳定的。下面给出三次样条插值函数的一些重要性质。质。

43、)(xS)(xf)(xS2.7.4 三次样条插值函数的性质三次样条插值函数的性质第二章 插值法值函数值函数 ,则有估计式则有估计式 定理定理2.9 设函数设函数,)(4baxfC 记记)(max,)(max110)4(4xxfMiinibxahx 则对任意则对任意,bax 满足边界条件（满足边界条件（2.44）或（）或（2.45）的三次样条插）的三次样条插)(xS（2.69）2,1,0,4)()(4)()(kkkxxMhCSfkk其中其中.81,241,3845210CCC由于误差估计与收敛性定理的证明比较复杂，下面只给出误差估计的结论。由于误差估计与收敛性定理的证明比较复杂，下面只给出误差估

44、计的结论。第二章 插值法 误差估计式（误差估计式（2.69）除可以用于误差估计外）除可以用于误差估计外,它进一步表明，它进一步表明，当当 时，在插值区间时，在插值区间 上，对于满足边界上，对于满足边界条件（条件（2.44）或（）或（2.45）的插值函数）的插值函数 ，不仅，不仅 一致收一致收敛于敛于 ，而且，而且 一致收敛于一致收敛于 ，一致收敛于一致收敛于 。,)(4baxfC)(xS)(xS)(xf )(xS )(xf )(xf)(xS,ba 用三次样条绘制的曲线不仅有很好的光滑度，而且当节用三次样条绘制的曲线不仅有很好的光滑度，而且当节点逐渐加密时，其函数值在整体上能很好地逼近被插函数，

45、点逐渐加密时，其函数值在整体上能很好地逼近被插函数，相应的导数值也收敛于被插函数的导数，不会发生龙格现象相应的导数值也收敛于被插函数的导数，不会发生龙格现象。因此三次样条在计算机辅助设计中有广泛的应用。因此三次样条在计算机辅助设计中有广泛的应用。第二章 插值法 插值插值 名称名称 插值插值 条件条件 构造构造 方法方法优点优点简便计简便计算方法算方法 缺点缺点L-插值插值先求基函数先求基函数再求插值函数再求插值函数 对称对称秦九韶秦九韶算法算法(1)计算量大计算量大(2)数值不稳定、不数值不稳定、不 收敛收敛(Runge现象现象)nixfxLiin,1,0),()(本章内容：本章内容：H-插值

46、插值一、一、同上同上二、利用牛二、利用牛 顿插值顿插值 收敛收敛 同上同上高次有高次有Runge现象现象 分段分段 插值插值 基函数法基函数法收敛收敛 同上同上 光滑性差光滑性差三次样三次样条插值条插值 基函数法基函数法收敛收敛 同上同上 光滑性光滑性 有改进有改进,2,1;,1,0),()(),()()()(22 lnixfxHxfxHililniin改进方法：列维尔算法、埃特金算法、牛顿法改进方法：列维尔算法、埃特金算法、牛顿法nimxIxfxIiiiih,1,0,)(),()(nixfxSii,1,0),()(注：注：分段插值中乍看上去没有构造插值基函数，但实际上线性插分段插值中乍看上去没有构造插值基函数，但实际上线性插32)(,)(,1jjjxxxxxx 值用了线性值用了线性L-L-插值基函数，三次样条插值的基函数可为：插值基函数，三次样条插值的基函数可为：