第4章插值逼近研究报告.ppt

本章要求 1．熟练掌握拉格朗日插值、牛顿插值的基本思想与误差估计； 2．掌握分段低次插值的概念算法； 3．掌握三次样条插值及其算法的计算机实现； 重点：多项式插值与误差估计、三次样条算法的计算机实现。 解之得， 显然， l0(x), l1(x) 满足：线性插值函数可以写成节点上函数值的线性组合，即 L1(x) = l0(x) y0 + l1(x) y1 ．称l0(x), l1(x) 分别为x0, x1的插值基函数。 二次插值的误差 例4.1.2 给定sin11°=0.190809,sin12°=0.207912,  sin13°=0.224951, 构造二次插值，并计算  sin11°30′。 例4.1.3 要制作三角函数sin x的值表，已知表值有四位小数，要求用线性插值引起的截断误差不超过表值的舍入误差，试确定其最大允许的步长。 解 f (x)= sin x, 设 xi-1, xi为任意两个插值节点，最大允许步长 记为 h = hi = xi －xi-1， 先构造 l0(x)。由上表知， x1 , x2, …, xn 为 l0(x) 的零点，设 于是可得，基函数为 所以我们得到 n 次Lagrange插值多项式： 例4.1.4 已知插值点 (-2.00,17.00), (0.00,1.00), (1.00,2.00), (2.00,17.00), 求三次插值，并计算 f (0.6)。 解 先计算4个节点上的基函数： 三次Lagrange插值多项式为： 一般地，k 阶差商 f [x0, x1,…, xk] 定义为： 一阶差商：f (x)关于点x0，x1的一阶差商记为 f [x0, x1]， 性质 2 各阶差商具有对称性, 即改变差商中节点的次序不会  改变差商的值。设i0, i1, …, ik为0, 1, …, k的任一排列, 则 性质1 k 阶差商 f [x0, x1,…, xk]可表成节点上函数值 f (x0), f (x1), …, f (xk) 的线性组合，即 性质 4 若 f (x)在 [a, b] 存在 n +1阶导数，xi∈ [a, b] , i = 0,1,…,n, 固定 x∈[a, b], 则 n+1 阶差商与导数 存在如下关系：  由性质1知，任意改变节点的次序，只改变公式右端求和的 次序，故其值不变。例如,由定义知， 由差商的定义 一阶差商是由节点上函数值定义的，二阶差商是由一阶差商定义的，…，依此构造差商表： 例 计算 (?2, 17), (0, 1), (1, 2), (2, 19)的一至三阶差商。 解 由表易知，   利用差商与导数的关系（性质4）！ 解 二次Newton插值 给定三个互异插值点(xi , f (xi)), i =0,1,2. 设 n 次Newton插值公式 给定n+1个插值点(xi, f (xi)), i = 0, 1, 2,…, n, xi 互异， 插值误差为： 2. 牛顿插值的误差不要求函数的高阶导数存在，所以更具有一般性。它对 f (x)是由离散点给出的函数情形或 f (x)的导数不存在的情形均适用。 3. 引入记号: f [x0] = f (x0) , t0(x) = 1, t1(x) = x ? x0, t2(x) = (x ? x0)(x ? x1), …,tn(x) = (x ? x0)(x ? x1) … (x ? xn?1)， 于是 n 次Newton插值公式可表为 ti(x) = ti?1(x)(x ?xi?1), i =1,2,…, n； ti(xj) = 0， j i，  ti(xj) ≠0，j ≥i. 4. Newton插值具有承袭性质，即 例 给定四个插值点(?2,17), (0,1), (1,2), (2,19), 计算 N2(0.9), N3(0.9)。 解 首先列差商表 Newton插值的算法 用牛顿插值公式，首先要计算各阶差商值，然后再计算插值。 牛顿插值由承