数值分析-插值方法.ppt

第一章 插值方法 高 云 插值方法的意义 插值方法的应用 对于早期的插值问题来说，f(x)通常是已 知函数，比如对数函数，指数函数，三角 函数等，并且已经有了这些函数值列表， 插值法可以用来计算那些不在表中的点处 的函数值。对于这一类问题来说，现在已 经不需要用插值方法来计算。 插值方法的应用 对于现在的许多实际问题来说，我们并不 知道f(x)的具体形式，所对应的函数值可 能是由测量仪器或其他物理设备中直接读 出来的，f(x)只是一个数学概念意义下的 函数。 （比如：图像的方法处理，天气预报，机床加工等方面） 泰勒插值 泰勒插值余项 拉格朗日插值 问题2 求作 n 次多项式 pn（x），使满足条件 pn（xi）= yi，i = 0，1，…，n （2） 这就是所谓拉格朗日（ Lagrange）插值. 线性插值 抛物线插值 插值问题的可解性 插值问题的可解性 多项式插值定理 运用基函数法求拉格朗日问题 基函数表 构造基函数 基函数的一般形式 基函数插值的一般表达式 插值余项 如何推导插值余项 插值余项 插值误差举例 插值误差举例 拉格朗日插值的几点问题 特殊情况 误差的事后估计 埃特金算法的迭代原理 埃特金算法 1.5 牛顿插值公式 提出的原因： 1 拉格朗日插值每增加一个新点都要重新计算插值公式。 2 埃特金算法虽具有承袭性，但算法是递推型的，不便于进行理论上的分析 3 牛顿公式具有承袭性并且理论推导严密 插商（均差）及其性质 插商及其性质 插商及其性质 插商及其性质 插商（均差）的性质 插商（均差）的性质 插商（均差）的性质 插商（均差）的性质 插商（均差）的性质 例题1 插商形式的插值公式 插商形式的插值公式 插商形式的插值公式 插商形式的插值公式 牛顿插值公式 牛顿插值公式 牛顿插值公式 例题2 例题3 例题3 例题3 拉格朗日插值与牛顿插值的比较 差分形式的插值公式 差分的重要性质 差分的性质 称 为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 假设 f（x）在含有点 x0的区间［a，b］内有直到 n +1阶导数，则当 x∈［a，b］时，对于由式（1）给出的 pn（x），成立 式中 ξ界于 x0与 x之间，因而 ξ∈［a，b］. (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 满足条件 , 存在 使得 。 且 推广：若 使得 使得 罗尔定理 : 若 在［ ］连续，在 充分光滑， ? 问题： 对于相同的插值公式，内插与外推哪一个的精度高。 插值点越多得到插值公式的精度越高？ 拉格朗日插值对于不同的初始函数，在相同点上的插值公式也不同。 多项式插值是唯一的插值方式？ 基函数的形式只和插值点的x坐标相关，和y值无关。 由n个点插值得到的基函数的次数必定是n-1次的多项式 基本假设与依据是 假设 f″（x）在［a，b］内改变不大 事后估计法 1．差商（均差）的定义 定义1：设有函数f (x)以