第1章插值方法.pptVIP

1.6 分段三次埃尔米特插值 第29页，共74页，编辑于2022年，星期五 分段三次埃尔米特插值的算法： 输入插值节点xi(i=0,…n) 计算插值节点上的原函数f(xi)和其 输入插值点的个数m for i=0 to n-1 h=(xi+1-xi)/(m-1) for k=0 to m x=xi+h*k ll(k)=h1f(xi)+h2f(xi+1)+h3f’(xi)+h4f’(xi+1) yy(k)=1/(1+x2) next k 绘制曲线 next i 第30页，共74页，编辑于2022年，星期五 1.7 三次样条插值 “样条”是指在工业设计过程中为了描绘出光滑的外形曲线所用的一种工具，即一个具有弹性的细长木条。 绘图员首先把一串有序点标在平面上，让样条依次通过这串点，并在每个点上压一块重的压铁，然后用铅笔沿着这样一根样条绘出一条光滑的曲线。 第31页，共74页，编辑于2022年，星期五 1.7 三次样条插值 人们由此得到启发：如果建立起样条的数学模型，并编成程序存放在计算机中，那么只要输入一串有序点列的坐标，再经过一些简单的计算，绘图机便能划出一条光滑的曲线。 这样的数学模型经简化后为分段三次多项式，在节点处左右两端的曲线的切线和曲率是连续的，这种模型称为三次样条函数。 第32页，共74页，编辑于2022年，星期五 1.7 三次样条插值 定义：给定[a,b]的划分：a=x0<x1<…<xn=b,如果函数S(x)在区间［a,b]上满足以下条件： 1.在每一个子区间（Xi,Xi+1）(i=1,2, …,n-1)上S(x)是三次多项式； 2. S(x)在区间［a,b]上具有连续的二阶导数； 3. S(x)=Yi (i=0,1,2, …,n), S′(x0)=y′0, S′(xn)=y′n 我们称S(x)为三次样条函数。 第33页，共74页，编辑于2022年，星期五 其中 第34页，共74页，编辑于2022年，星期五 第35页，共74页，编辑于2022年，星期五 其中 第36页，共74页，编辑于2022年，星期五 第37页，共74页，编辑于2022年，星期五 第38页，共74页，编辑于2022年，星期五 第39页，共74页，编辑于2022年，星期五 第40页，共74页，编辑于2022年，星期五 第41页，共74页，编辑于2022年，星期五 第42页，共74页，编辑于2022年，星期五 第43页，共74页，编辑于2022年，星期五 第44页，共74页，编辑于2022年，星期五 第45页，共74页，编辑于2022年，星期五 第46页，共74页，编辑于2022年，星期五 第47页，共74页，编辑于2022年，星期五 第48页，共74页，编辑于2022年，星期五 第49页，共74页，编辑于2022年，星期五 第50页，共74页，编辑于2022年，星期五 第51页，共74页，编辑于2022年，星期五 第52页，共74页，编辑于2022年，星期五 第1页，共74页，编辑于2022年，星期五 1.1 拉格朗日插值公式 拉格朗日插值公式的基本思想是:把Pn(x)的构造问题转化为n+1个插值基函数Li(x)（i=0,1,…,n)的构造。插值函数Pn(x)是这n＋1插值基函数的线性组合，其组合系数就是对应点上的函数值。 第2页，共74页，编辑于2022年，星期五 拉格朗日插值的算法： Input X, (Xi,Yi), I=1,1,…,n 0?S, 1?k FOR K=1 TO n t=1 FOR j=1 TO n j<>k .F. .T. t*(X-Xj)/(Xk-Xj) ?t S+t*Yk ?S Output S 第3页，共74页，编辑于2022年，星期五 绘制拉格朗日插值基函数图形 输入插值节点xi(0,…n) 输入插值点的个数 m h=(xm-x0)/(m-1) 计算m个插值点上的插值基函数l(i,k)和 ll(k)=∑l(i,k) for k=0 to m xx(k)=x0+k*h ll(k)=0 for i=0 to n l(i,k)=1 for j=0 to n if j<>i