计算方法之插值.pptVIP

* 将C和D回代到S(x) （第二次积分式），得 4又由于S(x)在[a,b]上具有二阶连续导数，函数曲线连续光滑，故有： S’(xi-0)=S’(xi+0) 即在点 xi 处左右导数相等。 * 在区间[xi , xi+1]上S(x)的导数函数为： 而在区间[xi-1 , xi]上S(x)的表达式为： * 那么在区间[xi-1 , xi]上S(x)的导数函数为： 由 及 * 由二者相等可得： * 可得 其中 * 5再按不同的边界条件进行讨论： （1）第I类边界条件（例如自然边界条件）： f”(x0)=f”(xn)，即M0=Mn=0，则上式可改写为n-1阶的线性方程组： * （2）第II类边界条件： 给定f’(x0)=S’(x0)=m0，f’(xn)=S’(xn)=mn，将m0和mn依次代入S’(x)在区间[x0 , x1]和[xn-1 , xn]中的表达式，则可得到： 于是得到n+1次线性方程组： * 4、总结 三次样条插值的主要任务是：根据给定的插值条件和边界条件，构造出大M关系式（关于Mi的线性方程组）。而方程组的系数矩阵和常数项矩阵各元素由已知插值条件和边界条件计算而出。 * [例] 给定离散数值表如下，取M0=M3=0构成三次样条插值的M关系式，并计算f(1.25)： xi 1.1 1.2 1.4 1.5 yi 0.4000 0.8000 1.6500 1.8000 解：由题中(xi , yi)的数值可得： h0=0.1，h1=0.2，h2=0.1， 由M0=M3=0的边界条件，得 * 解得：M1=13.125，M2=-31.875。 将M0=0、M1、M2、M3=0代入区间[xi , xi+1]上的S(x)： 特别地：f(1.25)?S(1.25)=1.0436 * 一、问题 埃 尔 米 特 插 值 对于被插值函数 y=f(x) ，已知如下表所示条件： xi x0 x1 …… xn f(xi) f(x0) f(x1) …… f(xn) …… …… 求解插值函数 P(x) ，使满足： （1）P(xi)= f(xi) （2） P’(xi)= f’(xi) * （3） P”(xi)= f”(xi) …… ——此即埃尔米特（Hermite）插值 二、三次埃尔米特插值 H3(x) 对于函数y=f(x)，已知 xi x0 x1 yi=f(xi) y0=f(x0) y1=f(x1) 构造三次多项式H3(x)，使满足插值条件 * 方法——将三次多项式H3(x)预设为 由插值条件 * 由插值条件 * （1）构造 令： 即能满足 再由： * 将A、B回代到 ，其中 同理可得： ，其中 * （2）构造 令： 即能满足 再由 同理可得： 总结 * 例题—用H3(x)近似计算f(x)=lnx在0.60处的函数值f(0.6)，插值节点与对应函数值如下表 xi x0=0.40 x1=0.70 yi=f(xi) y0=-0.916291 y1=-0.356675 mi=f’(xi) m0=2.5 m1=1.42857 解： * 再将x=0.60代入上式计算即可。 代入下式，即得三次埃尔米特插值多项式： * 三、2n+1次埃尔米特插值多项式H2n+1(x) 对于函数y=f(x)，已知 xi x0 x1 …… xn yi=f(xi) y0=f(x0) y1=f(x1) …… yn=f(xn) …… 按下式构造2n+1次埃尔米特插值多项式H2n+1(x) * 习题选讲 习题1.2求绝对误差限、相对误差限及有效数字 定义（有效数字与绝对误差限）： 设有近似数x=?10k?0.a1a2…an或x=?10k?0.a1a2…,若有|e|=|x*-x|?0.5?10k-n，称x具有n位有效数字。 定理（有效数字与相对误差限）：若设有近似数x=?10k?0.a1a2…an，若有n位有效数字，则（相对误差限） * 1） 0.3012= 0.3012?100：k=0，有效数字位数n=4，则绝对误差限|e|=|x*-x|?0.5?10k-n= 0.5?10-4 a=3，相对误差限： 2） 30.12= 0.3012?102：k=2，有效数字位数n=4，则绝对误差限|e|=|x*-x|?0.5?10k-n= 0.5?10-2 a=3，相对误差限： * 3） 30.120= 0.30120?102：k=2，有效数字位数n=5，则绝对误差限|e|=|x*-x|?0.5?10k-n= 0.5?10-3 a=3，相对误差限： 4） 30120= 0.30120?105：k=5，有效数字位数n=5，则绝对误差限|e|=|x*-x|?0.5?10k