第7章_插值法.ppt

结束 对于分段线性插值的余项估计有下列结果 定理7.3 设给定节点为a=x0＜x1＜…＜xn=b, f(xi)=yi, f ”(x)在［a,b］上存在,则对任意的x∈［a,b］有: * （证明） 结束 * 7.5.2 分段埃尔米特插值 分段线性插值,从整体上看,当M2+∞时,较小的h可使插值函数?(x)较好地逼近被插函数f(x),但在节点xi处失去了光滑性,实际上?(x)为连接(xi-1, yi-1)与(xi , yi ), (i=1,2,…,n)的折线.如果不仅要求?(x)在节点满足?(xi )= yi,而且要求? ′(xi)=f′(xi ),则可采用分段Hermite插值： 设函数f(x)在[a,b]上有定义,且对给定的n+1个节点a= x0 x1 … xn =b上 f(xi )= yi , f′(xi )=y′i, 求作一个插值函数H(x),使H(x)满足： 1) H(xi )= xi , H ′(xi)=y′i , (i=1,2,…,n) 2) H(x)在[xi-1 , xi], (i=1,2,…,n)上为三次多项式 由前面对两节点x0,x1上埃尔米特插值的讨论,用于区间[xi-1, xi]上,不难构造出插值基函数： 结束 * 结束 * 结束 * 则分段三次埃尔米特插值函数为 分段三次埃尔米特插值的余项 定理7.4 设给定节点为a= x0 x1 … xn =b, f(xi)= yi , f′(xi)=y′i ,f(4)(x)在[a,b]上连续,则对任意x?[a,b],有 其中 (证明类似定理7.3,见教材.) 结束 * §7.6 样 条 插 值 样条(Spline)插值,是一种既能克服高次多项式插值的缺陷,又能保证一定光滑性的插值,即具有一定的光滑性的分段插值. 样条一词来源于工程中的样条曲线.绘图员在描绘过已知点(样点)的光滑曲线时,往往用细长的木条或金属条(样条)把邻近的几个点连接成一条光滑的曲线.这一过程,就力学而言,相当于细梁(样条)在若干点处受到集中载荷时的挠曲.这样描绘出的曲线,实际上是一段段的三次曲线拼成的曲线,在拼接处,不仅连续,而且它的一阶及二阶导函数也是连续的,通常把这样得到的函数叫做样条函数. 结束 * 7.6.1 样条插值的基本概念 设在区间[a,b]上取n+1个节点a= x0 x1 … xn =b,其对应的函数值为f(xi)= yi , (i=1,2,…,n) ,现求一定义在[a,b]上的函数S(x),使其满足： 1) S(x)在每一个小区间[xi-1 , xi], (i=1,2,…,n)上为三次多项式, 2) S(x)在[a,b]上二阶连续可微,即S(x) ∈C2[a,b], 3) S(xi) = yi , (i=1,2,…,n), 则称S(x)为f(x)的三次样条函数. 为确定S(x)在各区间上的表达式,自然会想到前节介绍的分段埃尔米特插值.在每一小区间上,分段埃尔米特插值也恰为三次多项式,不过在那里,f ′(xi)=y′ i为已知的.如果能利用S(x)在[a,b]上二阶连续可微将y′ i确定,那么由(7.19), S(x)便可唯一地表示出来了.为此先记mi =f ′(xi), (i=1,2,…,n),在区间[xi-1 , xi]上, S(x)的表达式为： 结束 * 7.6.2 样条插值公式 记hi-1= xi - xi-1 ,求S”(x)的二阶导函数 结束 * 再利用S(x)在[xi , xi+1]上的表达式可计算出 由S(x)二阶连续可微,即S”(xi-) = S”(xi+),得 记： 将上式整理得方程组 结束 * 这是关于n+1个未知量mi,(i=0,1,…,n)的n-1个线性方程组,该方程组有无穷多组解,在实际问题中,往往根据具体情况补充两个附加条件—通常称为端点条件,便可唯一确定一组解.常见的端点条件有： 1)曲线在两端点x0 , xn处的导数值为已知的,即f′(x0)= m0, f′(xn)=mn,方程组为n-1个未知数,n-1个方程,有唯一解. 2)函数f (x0)在两端点x0 , xn处的二阶导数为0, 即f”(x0)= f”(xn) =0, 即S”(x0)= S”(xn) =0 由S”(x0+ ) =0可得 由S”(xn- ) =0可得 得方程组 也有唯一解. 结束 * 3) 周期端点条件(当f (x0)=f (xn)时)：S(x0)=S(xn),S′(x0)=S′(xn),S”(x0)=S”(xn) 于是有： 与前面n-1个方程联立,也可唯一地解出m0,m1, …,mn. 对于前两种端点条件所对应的方程组可化为下列称之为三对角方程组