数值分析san(研究生).pptVIP

2.3 差商与牛顿插值 牛顿插值多项式的构成 2.6 分段插值法 给定 (x∈［-5,5］)。取等距节点 xi=-5+i(i=0,1,…,10), 试建立插值多项式L10(x), 并作图形, 观察L10(x)对f(x)的逼近效果。 分段三次埃尔米特插值 为了避免Runge现象的发生, 很自然地会想到把区间［-5, 5］等分为10个小区间, 在每一个小区间内应用低次插值。但由于每个小区间只有两个端点（插值节点）, 按照已知的方法, 得到的将是一个分段线性插值函数。 已知xi,f(xi),f(xi)(i=0,1,…,n),求分段三次插值函数H(x)满足 H(xi)=f(xi),H(xi)=f(xi) i=0,1,…,n 为了得到插值函数,考虑任意子区间xi,xi+1］,i∈(0,1,…,n-1), 采用Lagrange插值函数结构, 在第i个子区间上H(x)=f(xi)h1(x)+f(xi+1)h2(x)+f(xi)h3(x)+f(xi+1)h4(x) 这样，就把H(x)的构造问题转化为四个插值基函数hk(x)(k=1,2,3,4)的构造问题。 “样条”一词本来是指在飞机或轮船设计过程中为了 描绘出光滑的外形曲线所用的一种工具，即一个具 有弹性的细长木条。事实上，在作了某些近似简化 后，样条的数学模型并不复杂，它只是分段的三次 多项式曲线：在相邻两块压铁之间是三次多项式曲 线；在压铁处，左右两段曲线的切线和曲率是连续的。 （3） 差商和导数 若 f(x) 在[a,b] 上存在n 阶导数，由余项表达式可得，n 阶差商与导数有如下关系 按差商的定义 4.3.2牛顿插值公式 1。牛顿插值公式的建立 牛顿插值多项式 (f(x)的前n+1项） 牛顿插值余项 (f(x)的最后一项） ……. 1 乘积因子 二阶差商 一阶差商 零阶差商 2。牛顿插值的特点 (1) P(x)次数不超过n次，项数不超过n+1项。各项系数是差商表上对角线的各阶差商值。 (2) P(x)满足插值条件，在节点上f(xi )=P(xi ) . (3) 增加一个节点，只需增加一项。 n次牛顿插值多项式 计算牛顿插值多项式的步骤 (1) 作差商表. (2) 写出牛顿插值多项式（表中对角线上各差商值就是P(x)的各项系数）。 (3) 计算插值点的近似值。 余项公式 n次牛顿插值多项式 又证:牛顿插值多项式 根据线性插值的点斜式 可得牛顿差商型线性插值多项式： 设 由 可得 牛顿差商型二次插值多项式： n次牛顿插值多项式 解 先作差商表 xi f(xi ) 1阶 2阶 3阶 4阶 0.40 0.55 0.65 0.80 0.90 1.1160 1.1860 1.2757 1.3841 0.41075 0.57815 0.69675 0.88811 1.02652 0.2800 0.3588 0.4336 0.191 0.214 0.034 由Newton公式得四次插值多项式 2.4 差分与等距插值 前面介绍了节点任意分布的插值公式，但实际应用时经常遇到等距节点的情形，这时插值公式可进一步简化，计算也简单很多。为此，我们先介绍差分的概念。 1. 差分的定义 2.5 埃尔米特插值 * 第2章 插值法 2.1 问题的提法 2.2 拉格朗日插值 2.3 差商与牛顿插值 2.4 差分与等距节点的牛顿插值 2.5 埃尔米特（Hermite）插值 2.6 分段插值法 2.7 三次样条插值 定义 函数y=f(x)在区间[a,b]上有函数值 即 x0 x1 x2 …… xn y0 y1 y2 …… yn 其中x0 ,x1,x2 ,…,xn是区间[a,b]上的互异点，要构造一个 简单的函数 作为f(x)的近似表达式，使满足 (插值条件 ) 这类问题称为插值问题。 -----f(x)的插值函数, f(x) -----被插值函数 x0 ,x1,x2 ,…,xn -----插值节点， [a,b]称为插值区间 求插值函数的方法称为插值法。 若x∈[a,b]，可计算f(x) 的近似值φ(x), 则 x