哈尔滨工程大学工程算法课件04插值方法_课件.pptVIP

* 插值法中的拉格朗日插值多项式是研究数值微积分与微分方程数值解的重要工具。牛顿插值多项式是拉格朗日插值多项式的变形，具有承袭性，比拉格朗日插值多项式节省计算量。分段低次多项式插值由于具有良好的稳定性与收敛性，且算法简单，便于应用。特别是应用广泛的三次样条插值，不但有较好的稳定性和收敛性，而且具有较好的光滑性，从而满足了许多实际问题的要求。 * hermite插值是用一条曲线来逼近，最高次数可能高于三次；三次样条插值是用连续的曲线来逼近，最高次数是三次。 三次样条与分段 Hermite 插值的根本区别在于S(x)自身光滑(考虑了二阶导数)，不需要知道 f 的导数值（除了在2个端点可能需要）；而Hermite插值依赖于f 在所有插值点的导数值。(S(x)为插值基函数,f为你要插值的函数) * * * 若用Pn(x)近似f(x)，则误差为Rn(x)。当x=xi(i=0,1, 2,…,n)时， Rn(xi )=0， p(xi)= yi ；而当x?xi时， Rn(xi ) ? 0， 所以p(xi)? f(x) * 若用Pn(x)近似f(x)，则误差为Rn(x)。当x=xi(i=0,1, 2,…,n)时， Rn(xi )=0， p(xi)= yi ；而当x?xi时， Rn(xi ) ? 0， 所以p(xi)? f(x) * 各阶差商可按下表排列方式逐列进行计算，称下表为差商表 * 实用中常采取等距节点 xi+1 - xi = h，利用此特点，引入差分概念，可简化计算 * f(x)=f(x0)+(x- x0)f[x1,x0]+(x- x0)(x-x1)f[x1,x0,x] f[x1,x0,x] = (x-x2) f[x2,x1,x0,x] +f[x2,x1,x0] f(x)=f(x0)+(x- x0)f[x1,x0] + (x- x0)(x-x1)f[x2,x1,x0] + (x- x0)(x-x1)(x-x2) f[x2,x1,x0,x] * Nn(x) Rn(x) 其中Nn(x)称为牛顿插值多项式 Rn(x)称为牛顿插值余项 * xi f[xi] f[xi,xi+1] f[xi,xi+1,xi+2] f[xi,xi+1,xi+2] x0 f(x0) x1 f(x1) f[x0,x1] x2 f(x2) f[x1,x2] f[x0,x1,x2] x3 f(x3) f[x2,x3] f[x1,x2,x3] f[x0,x1,x2 ,x3] … … … … 4.4.2 牛顿插值公式 * 定义 函数y= f(x)在区间[xi ,xi+1]上的平均变化率 自变量之差和因变量之差之比叫差商 称为f(x)关于xi , xi+1 的一阶差商,并记为f[xi ,xi+1] 二阶差商 m阶差商 * f[xi,xj,xk]是指 f[xi , xj , xk]= f[xj , xk]- f[xi , xj ] xk- xi 一般的,可定义区间[xi, xi+1 ,…, xi+n]上的n阶差商为 差商及其性质 * 差商表 xi f[xi] f[xi,xi+1] f[xi,xi+1,xi+2] f[xi,xi+1,xi+2,xi+3] x0 f(x0) x1 f(x1) f[x0,x1] x2 f(x2) f[x1,x2] f[x0,x1,x2] x3 f(x3) f[x2,x3 ] f[x1,x2,x3] f[x0,x1,x2 ,x3] … … … f[x1,x2]- f[x0,x1] x2 – x0 * xi f[xi] f[xi,xi+1] f[xi,xi+1,xi+2] f[xi,xi+1,xi+2 ,xi+3] 0 0 2 8 3 27 5 125 6 216 例：求 f(xi)= x3在节点 x=0, 2, 3, 5, 6上的各阶差商值 解: 计算得如下表 * 在n+1个节点处各阶差商的计算方法 差商及其性质 * 性质3 若f[x, x0, x1 , …, xk ]是 x 的 m 次多项式, 则 f[x, x0, x1 ,…, xk , xk+1]是 x 的 m-1 次多项式 证：由差商定义 右端分子为 m 次多项式, 且当 x = xk+1 时, 分子为0 ,故分子含有因子 xk+1 – x，与分母相消后，右端为m-1 次多项式。 * 性质4 若 f(x)是n次多项式, 则f [x, x0, x1 , …, xn ]恒为0 证： f (x)是n次多项式，则f [x, x0 ]是 n-1次多项式, f [x, x0, x1 ]是 n-2 次多项式, 依次