计算方法 第二章 插值与曲线拟合.pptVIP

Confidential, for review only 引言 许多实际问题都用函数 来表示某种内在规律的数量关系,其中相当一部分函数是通过实验或观测得到的.虽然 在某个区间[a，b]上是存在的，有的还是连续的，但却只能给出[a,b]上一系列点 这只是一张函数表；有的函数虽然有解析表达式,但由于计算复杂,使用不方便,通常也构造一个函数表。如三角函数表、对数表、平方根表、立方根表等等。 2.1 问题的提出 函数解析式未知,通过实验观测得到的一组数据, 即在某个区间[a, b]上给出一系列点的函数值 yi= f(xi) 或者给出函数表 例2.3 已知x=1, 4, 9 的平方根值, 用抛物插值公式, 求 由Rolle定理知： 的相邻两个零点之间至少存在一个零点，即 在（a,b）内至少有n+1个互异零点。 同理对 应用Rolle定理知： 在（a,b）内至少有n个互异零点，如此反复应用Rolle定理n+1次知： 在（a,b）内至少有一个零点 。 2.4 .1 差商及其性质 定义 函数y= f(x)在区间[xi ,xi+1]上的平均变化率 2.4 .1 差商及其性质 f[xi,xj,xk]是指 差商表 例2.11 求 f(xi)= x3在节点 x=0, 2, 3, 5, 6上的各阶差商值 解: 计算得如下表 f(x)=f(x0)+(x- x0)f[x1,x0]+(x- x0)(x-x1)f[x1,x0,x] 如当n=1时， f(x) = f(x0) + (x- x0)f[x1,x0] + (x- x0)(x-x1) f[x1,x0,x] Nn(x)= f(x0) + (x- x0)f[x1,x0] 例 2.12 已知 x = 1, 4, 9 的平方根值，求 解： 例2.15 求 并估计其误差 在区间[ 4 , 9 ]上， 特例（n=1) 例题 增加新节点x，并且f(x)为(n+1)阶可导时，有 (??[x0,x1,…,xn]) (??[x0,x1,…,xn,x]) |f(x)(n+1)|?Mn+1 2.4.2 牛顿插值余项 4.4 .1 差商及其性质 例2.13 已知 x=0, 2, 3, 5 对应的函数值为 y=1, 3, 2, 5 , 作三次Newton插值多项式。 xi f(xi) 一阶差商 二阶差商 三阶差商 0 1 2 3 1 3 2 -1 -2/3 5 5 3/2 5/6 3/10 ∴ 所求的三次Newton插值多项式为 如已知x = 0, 2, 3, 5, 6时，对应的函数值为y = 1, 3, 2, 5, 6,作差商表 -11/120 -1/4 -1/6 1 6 6 3/10 5/6 3/2 5 5 -2/3 -1 2 3 1 3 2 1 0 四阶差商 三阶差商 二阶差商 一阶差商 y x 四次牛顿插值多项式为 4.4 .1 差商及其性质 例2.14 已知 f(x) = x7+ x4+ 3x+ 1 求 f [20, 21, … 27 ] 及 f [20, 21, … 27, 28 ] 分析：本题 f(x)是一个多项式, 故应利用差商的性质 解: 由差商与导数之间的关系 解：作函数 f(x) = 取 x0=4, x1=9, x2=6.25 , 建立差商表 2.5 6.25 3 9 2 4 f[xi,xi+1,xi+2] f [xi,xi+1,] f(x) x N2(7)= 2 + (7-4)*0.2 + (7-4)*(7-9)*(-0.00808) = 2.64848 f 3(x) = Rn (x) 余式近似 0.5 *10 -2, N2(7) = 2.64848 可舍入为2.65 | f(x)(n+1) | ? Mn+1 由 许多实际问题不