第二章插值法-上海海事大学题材.ppt

第二章 插值法;通常 ;　从几何上看，插值问题就是求过n+1个点;其系数矩阵是n+1阶范德蒙(Vandermonde)行列式; ∵ xi≠xj ,(i≠j),∴此范德蒙行列式的值不为零,方程组有唯一解a0,a1,a2,?,an.; ;如记; 例 已知 ,求;二次插值 已知三点 (xi,yi) (i=0,1,2), 求一个二次多项式 P2(x) , 使 P2 ( xi ) = yi ( i = 0,1,2);用待定系数法可确定 。 ;仍叫拉格朗日二次插值基函数,;∴;;与二次插值类似可得;? n (xi) = yi (i =0,1,2,?,n)。令 Fn(x)=Pn(x) - ?n(x), 则Fn (x)仍为一次数不超过n的多项式。假设Fn(x)不恒等于零, 则方程Fn(x) = 0最多有n个根。但是;定理2 设f(x)在插值区间[a,b]上具有n+1阶连续导数,则对任意x?[a,b],至少存在一点??(a,b),使得;证: 设x为［a,b］内任意一点， 显然,当 x = xi ( i=0,1, … ,n ) 时，因为ωn+1 (xi) = 0， 所以 Rn (xi) = 0,所以公式成立. 当x ≠ xi时，对每一个这样的x均可作一辅助函数:; 这里,n+2个不同零点将划分出n+1个不同的小区间,应用 罗尔定理,每个小区间里有一个? ?(t)的零点,即? ?(t)有n+1个不同零点,同理? ?(t) 有n个不同零点….反复应用罗尔定理有,存在一点ξ,使得;罗尔定理： f(x)在[a,b]上连续,(a,b)上可导,且f (a)=f (b),则至少存在一点?? (a,b),使得f ’(?)=0成立.; 设函数 f (x)在 [a, b] 上有定义 ,xi,xj?[a, b], 且 xi ≠ xj , 则称;为 f (x) 在 xi, xj , xk的二阶差商。;例如一阶差商;;类似地;二 牛顿插值公式 设 x0 x1 , ?, xn 为插值节点 , x 为插值区间内任意一点,以 x, x0, x1 ,?, xn为节点, 依次作一阶、二阶,? , n+1阶差商,可得;*;依次将后式代入前式得;则;视 n 为 k-1, x 为 xk, 即可得差商的性质（4）; 例 已知 ,求 ,?=0.001。;（ 的准确值是0.6428）,三阶牛顿插值已相当准确。;插值节点为等距节点： xk =x0 + kh ，k =0, 1, ……, n ,;如果用函数表上的值，一阶中心差分应写成;而下面定义向后差分，▽表示向后差分算子，;由上面各种算子的定义可得算子间的关系：;验证： n=1 时，;n=3 时，;性质3: 在等距插值的情况下,由定义可得出差分和均差有 如下关系：;验证差分和均差有如下关系???;性质4: 各种差分之间可以互化;把差商与差分关系代入牛顿插值公式，;其余项为：;取节点;向前、向后差分表;例：在微电机设计计算中需要查磁化曲线表，通常给出的表是磁密B ,每间隔100 高斯磁路每厘米长所需安匝数at 的值，下面要解决B 从4000 至11000区间的查表问题。;从差分表中看到三阶差分近似于0 ，计算时只需两阶差分。 当4000≤B≤10500 时用牛顿前插公式； 当10500≤B≤11000 时用牛顿前插公式；;例如，求f（5200）时取 B0= 5000, f0=1.58, Δf0 = 0.11, Δ2 f0 = 0.01 ，;§ 5 埃尔米特(Hermite)插值;先求? i(x),由于?i (xj) =?ij,这一点与拉格朗日插值基函数 li(x)相同.所以?i (x)含有li (x)作为因子,又因为?′i(xj) =0,所以li(x)的根都是Hi(x)的二重根,所以?i (x)含有因子li2(x),显然li2(x)是2n次多项式,因此可令?i (x) =(ax+b) li2 (x);又? i′(x) = a li2(x) + 2(ax+b) li (x) li′(x) ;定理 插值多项式是唯一的;当仅两个节点x0,x1时,埃尔米特插值公式为：;例 求满足 ;其中A 为待定常数，可由条件;显然?(xj)=0 j=0,1,2. 且;Hermite 插值的一般形式:;例：按下表求Hermite 插值多项式;因为;因此所求Her