数值分析--第02章插值法.pptVIP

数 值 分 析 Numerical Analysis 机械与汽车工程学院 主讲人：孔胜利 kongsl@spu.edu.cn 2011-09-01 第2章 插值法 插值法的概念 拉格朗日插值多项式 Newton插值多项式 等距节点插值 Hermite插值 分段插值和抛物线插值 样条插值 2.1 插值法的概念 举例 已经测得在某处海洋不同深度处的水温如下： 深度（M） 466 741 950 1422 1634 水温（oC） 7.04 4.28 3.40 2.54 2.13 根据这些数据，希望合理地估计出其它深度（如500米，600米，1000 米…）处的水温。 当精确函数 y = f(x) 非常复杂或未知时，在区间[a,b]上一系列节 点 x0 … xm 处测得函数值 y0 = f(x0), …, ym = f(xm)，由此构造一 个简单易算的近似函数 g(x) ? f(x)，满足条件 这个问题称为“插值问题” 这里的 g(x) 称为f(x) 的插值函数。节点 x0 … xm称为插值节点,条 件（*)称为插值条件，区间[a,b]称为插值区间。 2.2 拉格朗日插值 2.3 Newton插值多项式 拉格朗日插值多项式构造简单，形式对称，计算方便，理论分析中 有重要的应用价值。但要想在计算中进一步提高精度，增加节点， 则要重新构造基函数，原来的计算要作废，这对实际计算很不利。 为了克服这个缺点，可把插值多项式表示为如下便于计算的形式 差商（也叫均差） 设 在 上定义，令互异的点 ，相应的函数 值，记两点上的一阶差商为，即 由定义知： 即差商具有对称性。 显然，一阶差商 是一元函数，再考虑它在点 的一阶差商，并记 ，即 称为点 上的二阶差商。 一般地，由m-1阶差商 及 ，再作两点 上的一阶差商，便得到 点上的m阶差商 均差计算表 例题：已知 在 点处的值分别为 计算 解 制差商表 根据问题知插值点x=4.01在 与 之间，故可用前三点 的二次插值多项式计算，即用 计算，代入数据，得 也可以取 作线性插值计算，即 代入数据，得 注：取七位有效数字的真值 2.4 等距节点插值 差分的定义 设函数 在等距节点 上的函数值 为已知，常数 叫做步长，则 分别称为函数 在点 的一阶向前差分，一阶向后差分。 利用一阶差分，可以定义高阶差分。例如： 二阶向前差分 二阶向后差分 一般地， 点的 n 阶向前差分 是 的线性组合。 向后差分 是 的线性组合。 差分表及其应用 Newton前插公式（表初公式） 用插值多项式作近似计算时，当插值点位于表初 附近，可用表初公式构造插值多项式。 令 ，插值点 ，则 表初公式 余项 Newton后插公式（表末公式） 用插值多项式作近似计算时，当插值点位于表末 附近，可用表末公式构造插值多项式。 令 ，插值点 ，则 表末公式 余项 例题： 已知 在7个点处的函数值 试计算 的值。 解 根据函数值作出差分表 由于五阶差分接近于零，可取四次插值多项式计算。插值点0.048位于 附近，故可用表初公式计算。有 ，知 因此 2.5 Hermite插值 Hermite插值 许多实际的插值问题不但要求在节点上函数值相等，而且还要求对应 的导数值也相等，甚至要求更高阶导数也相等，满足这种要求的插值 多项式就是Hermite插值多项式。构造Hermite插值多项式的方法就是 Hermite插值法。 设在节点 上，已知 要求构造满足该条件的插值多项式。 假设待构造的插值多项式H(x)需要满足以下插值条件 这里给出了2n+2个条件，可唯一确定一个次数不超过2n+1的多项式 其形式为 如果根据插值条件来确定系数，显然是非