[2018年最新整理]Ch2插值法.docVIP

  • 5
  • 0
  • 约2.93千字
  • 约 12页
  • 2018-02-20 发布于浙江
  • 举报
[2018年最新整理]Ch2插值法

Ch2. 插值法 §1. 插值问题 引例 矿井中某处的瓦斯浓度与该处距地面的距离有关,现用仪器测得从地面到井下500米每隔50米的瓦斯浓度数据,根据这些数据完成下列工作:(1)寻找一个函数,要求从此函数中可近似求得从地面到井下500米之间任意一点处的瓦斯浓度;(2)估计井下600米处的瓦斯浓度。 第一个问题可归结为“已知函数在处的值,求函数在区间内其它点处的值”,这种问题适宜用插值方法解决。 但对第二个问题不宜用插值方法,因为600米已超出所给数据范围,用插值函数外推插值区间外的数据会产生较大的误差。 解决第二个问题的常用方法是,根据地面到井下500处的数据求出瓦斯浓度与地面到井下距离之间的函数关系,由求井下600米处的瓦斯浓度。 定义 设在中个点处的值为已知,现根据上述数据构造一个简单函数,使,这种问题称为插值问题。 ,分别称为被插值函数、插值函数、插值节点和插值条件。若为多项式,则此问题称为多项式插值或代数插值。 定理1 在插值节点处,取给定值,且次数不高于的插值多项式是存在且唯一的。 证 令,则根据插值条件有下列等式: (关于的阶线性方程组), 其系数行列式是范德蒙(Vandermonde)行列式 。 根据克莱姆法则,此方程组存在唯一解,即存在且唯一。 §2. Lagrange插值 1. 线性插值与抛物插值 (1) 线性插值 , 其中称为线性插值的基函数。 (2) 抛。 分别令,即得 , 故 ,其中称为抛物插值的基函数。 2. Lagrange插值多项式 定义 对个插值节点,令 , 则显然。此时,满足,称之为Langrange插值多项式,称为Lagrange插值的基函数。 编程时宜用。 3. 插值公式的余项 定理2 设上连续,在内可导,则以插值多项式逼近的截断误差(即余项) 。 例1 已知函数的数据如下,分别用线性插值和二次插值求的近似值。 0.5 0.6 0.7 -0.693147 -0.510826 -0.356675 解 , 。 。 §3. 逐次线性插值法 对插值节点及对应的函数值,用表示一个非负整数序列,将个节点所确定的不高于次的插值多项式记为,则 , 即次插值多项式可以用两个次插值多项式通过线性插值获得—逐次线性插值。 Aitken算法: 例2 根据下表近似计算在处的值。 , , , , 。 同理可得。 ,类似可得。 ,故。 §4. 均差与牛顿插值公式 1. 均差(差商)及其性质 定义 称为关于的一阶均差; 称为关于的二阶均差; 类似地可定义阶均差 。 定理 ,即阶均差可表示为的线性组合,从而阶均差与节点的排列次序无关。 证 当时,右左 不妨假设成立,则 。 2. 牛顿插值公式 设插值节点上的插值多项式为 。 分别令,则有 ; ; 从而,,,, ,故 —— 牛顿插值公式。 例3 P34 0.40 0.55 0.65 0.80 0.90 1.05 K=0 0.41075 0.57815 0.69675 0.88811 1.02652 1.25382 K=1 1.11600 1.18600 1.27573 1.38410 1.51533 K=2 0.28000 0.35893 0.43348 0.52493 K=3 0.19733 0.21300 0.22863 K=4 0.03134 0.03126 K=5 -0.00012 , 。 。 §5. 差分与等距节点插值公式 1. 差分及其性质 定义 对等距节点,记。 向前差分:; 向后差分:; 中心差分:; 二阶差分:; 阶差分:。 不变算子与移位算子,即。 由,得。 (1) 差分类似于微分的性质 (2) 函数值与差分可相互线性表示 , 。 (3) 差分与均差的关系 。 证 时,。 假设时成立,则 。 2. 等距节点插值公式 令,则,, 。 代入牛顿插值公式得, 。 §6. Runge现象与高次插值的讨论 1. Runge现象 例4 ,节点,求插值多项式。 解 ——10次Langrange插值多项式。 结果如下: 2. 讨论 (1) 节点的增多固然能使插值函数在更多的地方与相等,但在两个节点之间不一定能很好地逼近,有时差异很大,所以在实际中,高次插值(7次以上)很少使用; (2) 可将分成若干小区间,在小区间内用低次(二次,三次)插值,即分段低§7. 三次样条插值 1. 三次样条插值函数 能够逐段表示成三次多项式且二阶导数连续(具有二阶光滑度)的函数。 定义 设,且在上为三次多项式,其中,则称为上的三次样条函数。 若对给定的,满足,则称为三次样条插值函数。 边界条件: 第一种边界条件: 固支梁条件; 第二种边界条件: 简

文档评论(0)

1亿VIP精品文档

相关文档