第二章样条.pptVIP

* * 样条插值 分段插值存在着一个缺点,就是会导致插值函数在子区间的端点(衔接处)不光滑,即导数不连续,对于一些实际问题,不但要求一阶导数连续,而且要求二阶导数连续。为了满足这些要求,人们引入了样条插值的概念。 所谓“样条”(SPLINE)是工程绘图中的一种工具,它是有弹性的细长木条,绘图时,用细木条连接相近的几个结点,然后再进行拼接,连接全部结点,使之成为一条光滑曲线,且在结点处具有连续的曲率。样条函数就是对这样的曲线进行数学模拟得到的。它除了要求给出各个结点处的函数值外,只需提供两个边界点处导数信息,便可满足对光滑性的不同要求。 一、样条函数的定义 设f(x)是区间[a,b]上的一个连续可微函数,在区间[a,b]上给定一组基点: a=x0x1x2???xn=b 设函数s(x)满足条件 (1) s(x)在每个子区间[xi , xi+1](i=0,1,2,???,n-1)上是次数不超过m的多项式; (2) s(x)在区间[a , b]上有m-1阶连续导数; 则称s(x)是定义在[a ,b]上的m次样条函数。x0，x1，x2, ???称为样条结点,其中x1,???,xn-1称为内结点, x0 , xn 称为边界结点。当m=3时,便成为最常用的三次样条函数。 二、三次样条插值函数 设y = f(x)在点 x0，x1，x2, ??? xn的值为y0，y1，y2, ???yn，若函数S(x)满足下列条件 S(xi)=f(xi) =yi , i=0,1,2,???,n (1.1) 则称S(x)为函数f(x)的三次样条插值函数, 简称 三次样条。 构造三次样条插值函数的方法有很多，这里介绍一个常用的方法：三弯矩插值法 记Mi = S″(xi), f(xi)= fi= yi ,考虑它在任一区间[xi,xi+1]上的形式.根据三次样条的定义可知 ,S(x)的二阶导数S″(x)在每一个子区间[xi,xi+1] ( i=0,1,2,???,n-1)上都是线性函数. 于是在[xi,xi+1] 上S(x)=Si(x)的二阶导数表示成 (1.2) 其中 hi= xi+1–xi . 对S″(x)连续积分两次,并利用插值条件S(xi)= yi ,得到 三次样条函数的构造 x x i , x i+1 S”(x) M i , M i+1 因此，只要能求出所有的{M i}，就能求出样条插值函数S(x). 下面考虑Mi的求法 则由连续性 S(xi-)= S(xi+) ,(i=1,2,……,n-1) 得 μiMi-1+2Mi+λiMi+1= di 其中 上面的方程组有n-1个方程，但有n+1个变量Mi，故需两个方程才能求唯一解，为此引入下列边界条件 下面介绍几种常用的边界条件 第一型边界条件： 已知f(x)在两端点的导数f(a)和f(b) ，要求 S′(a) = f′(a) , S′(b) = f′(b) 第二型边界条件： 已知f(x)在两端点的二阶导数f(a)和f(b) ，要求 S″(a)=M0 = f″(a) , S″(b)=Mn= f″(b) 特别当 S″(a)= S″(b) =0时，S(x)称为自然三次样条 第三型边界条件： 已知f(x)是以b -a为周期的周期函数 ，要求S(x)满 足周期条件 S (a) = S(b) , S′(a+)= S′(b-) , S″(a+)= S″(b-) 三次样条插值问题