计算方法0215.pptx

计算方法 2月15日 1 第2章 插值法 引言 拉格朗日插值 均差与牛顿插值多项式 埃尔米特插值 分段低次插值 三次样条插值 2 什么是插值问题？ 简单地说，给定 (x0,y0) ，(x1,y1) ， …… (xn,yn)， 给定 x，确定y =? 类似地，对于多变元函数. 给定 (x00, x01, …x0n, y0) … (xk0, xk1, …, xkn, yk) 以及给定 x0, x1, …, xk, 如何确定y? 3 严格定义 4 首选多项式插值 方便性 易于计算 易于微分 易于积分 唯一性 不管用什么办法获得，满足条件的n阶多项式插值函数（给定n+1个插值条件）只有一个 5 多项式插值函数的唯一性 6 而 恰为范德蒙(Vandermonde) 行列式。由线性代数知： 7 如何确定多项式插值函数 解析方法 解方程求解多项式系数ak 数值方法： 拉格朗日、牛顿…… 8 线性插值（一次多项式） n=1 2个插值节点 一次多项式插值函数 满足2个插值节点约束 9 线性插值基函数 一次多项式插值函数的另一种写法 10 线性插值基函数的性质 11 抛物线插值（二次多项式） n=2 3个插值节点 二次多项式插值函数 满足2个插值节点约束 12 二次插值基函数特性 仿照一次（线性）插值基函数构造一次插值多项式的原理，如果也有满足条件 的3个二次函数，则也可类似构造出二次插值函数，形如： 13 二次插值基函数 14 二次插值多项式 15 n次多项式插值 16 n次插值基函数 17 拉格朗日插值多项式 基于上述的n次插值基函数，可得到相应的n次多项式插值函数，即拉格朗日插值多项式： 18 拉格朗日插值法的tips 特殊情况下，n次拉格朗日插值多项式的次数可能小于n． 唯一性。定理1 在次数不超过n的多项式集合Hn中，满足n+1个插值节点约束条件的插值多项式Ln(x) ∈Hn是存在唯一的。（证明自学） 特例： 19 一个重要的结论 在上述特例中，最特的一个： 20 拉格朗日插值余项 21 证明: 只给出思路,详见课本 按定义, Rn(x)有根x0, x1, …, xn. 所以,可以假设 Rn(x) = K(x)?n+1(x). 其中K(x)待定． 把x也看作一个确定的值，构造一个函数 ?(t)=f(t) -Ln(t) - K(t)?n+1(t). 则?(t)有(n+2)个根x0, x1, …, xn, x.由Roll定理， ?’(t)在?(t)的两个根之间有一个根，所以在[a, b]上?’(t)有至少n+1个根. …… 反复应用Roll定理，到最后可以解得K(x).最终定理得证明． 22 拉格朗日插值法应用 注意： ??(a, b)一般是不好确定的．但是若我们能求出　　　在(a, b)内的界，则可得到截断误差的界！ 例2　见课本28页，演示见example201.m. 23 拉格朗日插值法的缺陷 拉格朗日插值法的优点是公式结构紧凑．不足在于当插值节点增减时全部插值基函数都要重新计算，这在实际计算中很不方便． 有没有其它方法，它可以克服该缺点，即增加节点时，在现有公式的基础上求解，从而降低计算量？ 有！再来? 24 习题 P48 2, 4, 5. 25