第1章插值方法.pptVIP

1.6分段三次埃尔米特插值第29页，共74页，星期日，2025年，2月5日分段三次埃尔米特插值的算法：输入插值节点xi(i=0,…n)计算插值节点上的原函数f(xi)和其输入插值点的个数mfori=0ton-1h=(xi+1-xi)/(m-1)fork=0tomx=xi+h*kll(k)=h1f(xi)+h2f(xi+1)+h3f’(xi)+h4f’(xi+1)yy(k)=1/(1+x2)nextk绘制曲线nexti第30页，共74页，星期日，2025年，2月5日1.7三次样条插值“样条”是指在工业设计过程中为了描绘出光滑的外形曲线所用的一种工具，即一个具有弹性的细长木条。绘图员首先把一串有序点标在平面上，让样条依次通过这串点，并在每个点上压一块重的压铁，然后用铅笔沿着这样一根样条绘出一条光滑的曲线。第31页，共74页，星期日，2025年，2月5日1.7三次样条插值人们由此得到启发：如果建立起样条的数学模型，并编成程序存放在计算机中，那么只要输入一串有序点列的坐标，再经过一些简单的计算，绘图机便能划出一条光滑的曲线。这样的数学模型经简化后为分段三次多项式，在节点处左右两端的曲线的切线和曲率是连续的，这种模型称为三次样条函数。第32页，共74页，星期日，2025年，2月5日1.7三次样条插值定义：给定[a,b]的划分：a=x0x1…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页，星期日，2025年，2月5日其中第34页，共74页，星期日，2025年，2月5日第35页，共74页，星期日，2025年，2月5日其中第36页，共74页，星期日，2025年，2月5日第37页，共74页，星期日，2025年，2月5日第38页，共74页，星期日，2025年，2月5日第39页，共74页，星期日，2025年，2月5日第40页，共74页，星期日，2025年，2月5日第41页，共74页，星期日，2025年，2月5日第42页，共74页，星期日，2025年，2月5日第43页，共74页，星期日，2025年，2月5日第44页，共74页，星期日，2025年，2月5日第45页，共74页，星期日，2025年，2月5日第46页，共74页，星期日，2025年，2月5日第47页，共74页，星期日，2025年，2月5日第48页，共74页，星期日，2025年，2月5日第49页，共74页，星期日，2025年，2月5日第50页，共74页，星期日，2025年，2月5日第51页，共74页，星期日，2025年，2月5日第52页，共74页，星期日，2025年，2月5日第1页，共74页，星期日，2025年，2月5日1.1拉格朗日插值公式拉格朗日插值公式的基本思想是:把Pn(x)的构造问题转化为n+1个插值基函数Li(x)（i=0,1,…,n)的构造。插值函数Pn(x)是这n＋1插值基函数的线性组合，其组合系数就是对应点上的函数值。第2页，共74页，星期日，2025年，2月5日拉格朗日插值的算法：InputX,(Xi,Yi),I=1,1,…,n0?S,1?kFORK=1TOnt=1FORj=1TOnjk.F..T.t*(X-Xj)/(Xk-Xj)?tS+t*Yk?SOutputS第3页，共74页，星期日，2025年，2月5日绘制拉格朗日插值基函数图形输入插值节点xi(0,…n)输入插值点的个数mh=(xm-x0)/(m-1)计算m个插值点上的插值基函数l(i,k)和ll(k)=∑l(i,k)fork=0tomxx(k)=x0+k*hll(k)=0fori=0tonl(i,k)=1