数值分析3-插值方法教学课件.ppt

插商及其性质 称 为f (x)在点x0, x1,…, xn处的n阶差商。 插商及其性质 f (x0) f (x1) f (x2) … f (xn?1) f (xn) f [x0, x1] f [x1, x2] … … … … f [xn?1, xn] f [x0, x1 , x2] … … … … f [xn?2, xn?1, xn] f [x0, …, xn] xn+1 f (xn+1) f [xn, xn+1] f [xn?1, xn, xn+1] f [x1, …, xn+1] f [x0, …, xn+1] xi yi 一阶差商 二阶差商 n 阶差商 … … 由差商定义可知：高阶差商是两个低一阶差商的差商。 x0 x1 x2 xn-1 xn 插商（均差）的性质 插商（均差）的性质 插商（均差）的性质 插商（均差）的性质 插商（均差）的性质 例题1 已知 插商形式的插值公式 插商形式的插值公式 插商形式的插值公式 插商形式的插值公式 牛顿插值公式 牛顿插值公式 牛顿插值公式 例题2 例题3 例题3 例题3 拉格朗日插值与牛顿插值的比较 差分形式的插值公式 向前差分 i i i f f f - = ? + 1 i k i k i k i k f f f f 1 1 1 1 ) ( - + - - ? - ? = ? ? = ? 向后差分 1 1 1 - - - ? - ? = ? i k i k i k f f f i?1 i i f f f - = ? 中心差分 其中 当节点等距分布时: 差分的重要性质 性质3：若 f (x)是 m 次多项式，则 是 性质1：常数的差分等于零 性质2：差分算子为线性算子 次多项式，且 性质4： 这个性质类比于 差分的性质 性质5： （类比于分部积分法则 ） 性质6：当节点xk是等距时，差分差商存在着关系： 第一章 插值方法 插值方法的意义 插值方法的应用 对于早期的插值问题来说，f(x)通常是已 知函数，比如对数函数，指数函数，三角 函数等，并且已经有了这些函数值列表， 插值法可以用来计算那些不在表中的点处 的函数值。对于这一类问题来说，现在已 经不需要用插值方法来计算。 插值方法的应用 对于现在的许多实际问题来说，我们并不 知道f(x)的具体形式，所对应的函数值可 能是由测量仪器或其他物理设备中直接读 出来的，f(x)只是一个数学概念意义下的 函数。 （比如：图像的方法处理，天气预报，机床加工等方面） 泰勒插值 条件： 泰勒插值余项 定理 1 假设 f（x）在含有点 x0的区间［a，b］内有直到 n +1阶导数，则当 x∈［a，b］时，对于由式（1）给出的 pn（x），成立 式中 ξ界于 x0与 x之间，因而 ξ∈［a，b］. 拉格朗日插值 问题2 求作 n 次多项式 pn（x），使满足条件 pn（xi）= yi，i = 0，1，…，n （2） 这就是所谓拉格朗日（ Lagrange）插值. 线性插值 (x0 ,y0) (x1 ,y1) P1(x) f(x) x0 x1 可见 P1(x) 是过 ( x0 , y0 ) 和 ( x1, y1 ) 两点的直线。 抛物线插值 x0 x1 x2 p2(x) ? f(x) f(x) 因过三点的二次曲线为抛物线，故称为抛物插值。 插值问题的可解性 设所求的插值多项式为 可建立关于系数 a0，a1，…，an的线性方程组 待定系数法 插值问题的可解性 克莱姆法则 范德蒙行列式 多项式插值定理 定理 (唯一性) 满足 的 n 阶插值 多项式是唯一存在的。 运用基函数法求拉格朗日问题 基函数的一般形式 要使得 则要求 依此类推要满足初始条件，所有基函数必须满足下列条件。 基函数表 构造基函数 由已知条件，假设 又因为 则 基函数的一般形式 即 基函数插值的一般表达式 插值余项 在[a , b]内存在, 考察截断误差 设节点 ，且 f 满足条件 , 存在 使得 。 且 推广：若 使得 使得 罗尔定理 : 若 在［ ］连续，在 充分光滑， 如何推导插值余项 插值余项 插值误差举例 插值误差举例 拉格朗日插值的几点问题 ? 问题： 对于相同的插值公式，内插与外推哪一个的精度高。 插值点越多得到插值公式的精度