数值分析-插值.ppt

* * 第四章 插值法 一、问题的提出 在实践中常出现这样的问题，由实验或测量得到一组数据，即 要求出其近似的函数表达式，也就是寻找一个简单函数 ，使 ，这类问题称为插值法。 二、基本概念 设函数y=f(x)在区间[a,b]上有定义，且已知在点 上的值为 ，若存在一个简单的函数 使 成立，则称 为 的插值函数，点 称为插值节点， 区间[a,b]称为插值区间，求 的方法称为插值法，条件（1.1）称为 插值条件。 图1.1 显然插值函数可以很多，其中最简单的是代数多项式， 这种插值函数叫做插值多项式。 于是问题变成： 求一个多项式 使 已知 这样的插值函数存在唯一性 o x0 x1 Xn-1 xn x y y=p(x) y=f(x) ∵ 可写成 由线性代数的知识知道 系数行列式 ≠0 方程有唯一解 又由范德蒙行列式可知 ∴满足条件（1.1）的多项式是存在且唯一的。 拉格朗日插值多项式 1. 点斜式： 对称式： 令 2. 3. 其中 例 解： 已知 例1：已知特殊角 的正弦函数值为 用一次插值， 二次插值多项式近似sinx，并用此近似式求sin500的值。 2. 为节点： ~ ~ 一次插值 为节点： 1. 3. 为节点： ~ ~ ~ ~ 二次插值 为节点： 考虑误差： 已知sin500=0.76604 ~ ~ ~ 由此可见： 1.高次插值比低次插值误差小。 2.内插比外推误差小。 3.节点之间距离越小，误差越小。 讨论误差： 设 在[a,b]上连续， 在（a,b)内存在，节点 是满足插值条件 的插值多项式，则对任意 插值余项 其中 且依赖于x的位置。 结论： 1.n越大，误差越小。 2.节点之间的距离越小，误差越大。 逐步线性插值（埃特金插值） 引进专用符号 表示以 为节点的k次拉格朗日插值公式 过 和 两“点”作线性插值。 两个k-1次插值多项式，把它们看作两点 以此“两点”作线性插值，推出I(X) 如： 埃特金插值公式 具体计算时先列表如下： . . . . . . . . . . . . . . . . . . 其中斜线上是1到4 次的插值多项式。 误差估计： 差商及牛顿插值公式 而 称为f(x)的二阶差商。 定义：给出函数f(x)在n+1个互异节点 上的函数值 称 为f(x)关于点 的一阶差商，简称一阶差商（均差）， 记作 即 称为f(x)的k阶差商。 性质2：在k阶差商 中，任意调换xi和xj的次序，其值不变， 这个性质称为差商的对称性。 性质1：k阶差商 是由函数值 线性组合而成。即： 如k=1 其中 综上所得： . . . 牛顿插值公式 . . . . . . . . . . . . . . . 计算时先列一个差商表如下： 其中斜线上的值就是牛顿插值公式中的系数。 牛顿插值公式是一个递推公式，当需要增加节点时， 只要多计算一行就可以了。 高次插值的误差分析 -1 o 1 0.5- 1.0 1.5- y=f(x) y＝L10(x) x y 例如，在给定区间〔-1，1〕上的函数 ，取等距节点， 如把〔-1，1〕等分，分点为 可以构造10次多项式，用拉格朗日公式写出为 其中 根据计算结果作图， 得出图形如右： （龙格现象） 分段低次插值 设给定n+1个节点 上的函数值 Ⅱ 在每一个小区间 内作二次插值： （从几何上表示用折线近似曲线，也称折线插值） Ⅰ 在每一个小区间 上用一次插值，则有 埃尔米特插值 满足以下条件的插值多项式称为埃尔米特插值多项式 设节点 上的函数值 及导数值 为已知，求插值多项式 满足条件 这里共给出2n+2个条件，可以确定一个次数不到2n+1次的多项式，其形式为： 问题的提出 定义2n+2个基函数： 解决的思路 于是多项式可以写成 *