6讲-ch2-5样条插值[录课].ppt

第五节 样条插值法;拉格朗日(Lagrange) 插值多项式，是插值基函数的 线性组合，形式简单，结构对称，便于理论分析。;牛顿(Newton)插值引入了差商的概念，在增加新的 节点时，只是增加一项，前面结果可再利用。;称为龙格Runge现象。;?;带导数的插值; 插值函数在子区间的端点(衔接处)不光滑, 从而导数不连续。; (1) S(x)在每个子区间[xi , xi+1](i=0,1,2,???,n-1)上是次数不超过m的多项式; (2) S(x)在区间[a , b]上有m-1阶连续导数;;Spline; 所谓“样条”(Spline)是工程绘图中的一种工具,它是有弹性的细长木条,绘图时,用细木条连接相近的几个结点,然后再进行拼接,连接全部结点,使之成为一条光滑曲线,且在结点处具有连续的曲率。样条函数就是对这样的曲线进行数学模拟得到的。 它除了要求给出各个结点处的函数值外,只需提供两个边界点处导数信息,便可满足对光滑性的不同要求。; 设y = f(x)在点 x0，x1，x2, ??? xn的值为y0，y1，y2, ???yn，若函数S(x)满足下列条件 （1）S(xi)=f(xi) =yi , i=0,1,2,???,n (1.1) （2）在每个子区间[xi , xi+1](i=0,1,2,???,n-1)上S(x)是三次多项式，记为 （3）S(x)在[a,b]上二阶连续可微。 则称S(x)为函数f(x)的三次样条插值函数, 简称 三次样条。; ;再次积分得：;化简得：;由上式可解出：;; 因此，只要能求出所有的{M i}，就能求出样 条插值函数S(x).下面考虑Mi的求法;化简得;;三次样条插值函数的构造; 上面的方程组有n-1个方程，但有n+1个变量Mi，故还需两个方程才能求唯一解，为此引入下列边界条件; 第三型边界条件： 已知f(x)是以b -a为周期的周期函数 ，要求S(x)满足 周期条件;化简得;联立方程组：;;说明; 已知函数f(x)的数值表如下： x 2 4 6 f(x) 3 7 13 f′(x) 1 -1 试求f(x) 在[2,6]上的三次样条插值函数 ; 得方程组 2 M0 + M1 = 3 0.5 M0 + 2M1 +0.5 M2 = 1.5 M1 +2 M2 = -12 解得：M0 =0.25 , M1 =2.5 M2 = -7.25;故所求的三次样条插值函数;三次样条插值函数的构造;Maths程序如下：;不能单独用c;g1=Plot[%[[1]],{x,2,4}] g2=Plot[%%[[2]],{x,4,6}] g3=ListPlot[B,Prolog-AbsolutePointSize[15]] Show[g1,g2,g3,Prolog-AbsolutePointSize[15]];三次样条插值函数的构造;样条插值图形;2、若满足S″(a)=M0 = f″(a) , S″(b)=Mn= f″(b），;三次样条插值函数的构造; 总结以上论述，可得求三次样条的步骤为：;;对条件(1)有;方程组为：;对条件(2)有;三次样条插值函数的构造;内容小结;二 元 样 条 曲 面