第2章 插值法part-1解析.ppt

■2. 已知y=f(x)在若干点的函数值, 求经过各点(xi, yi)的简单函数P(x). ■ 3. 函数f(x)难以计算,希望在区间[a,b]找到简单函数P(x)代替f(x),以便计算. 例子: f(x)=arctan(x), 在[-1,1]区间上, 可用 简单函数多项式P(x)=0.987733x-0.202335x3 代替 f(x) 插值法的研究内容 存在唯一性定理： 其系数行列式为Vandermonde行列式 现在求 及 ， 令 其中 从而有： 由此得： ， 故： ， * 由 的表达式可得： 于是得到： 同理可得 * 例：已知 分别利用 sin x 的1次、2次 Lagrange 插值计算 sin 50? 并估计误差。 解： n = 1 分别利用x0, x1 以及 x1, x2 计算 ?利用 这里 而 sin 50? = 0.7660444… ) 18 5 ( 50 sin 1 0 ? ? p L 0.77614 外推 /* extrapolation */ 的实际误差 ? ?0.01001 利用 sin 50? ? 0.76008, 内插 /* interpolation */ 的实际误差 ? 0.00596 内插通常优于外推。选择要计算的 x 所在的区间的端点，插值效果较好。 * n = 2 ) 18 5 ( 50 sin 2 0 ? ? p L 0.76543 sin 50? = 0.7660444… 2次插值的实际误差 ? 0.00061 高次插值通常优于低次插值 Return * 第二章 习题 P. 42 第2题 第2题改用L二次 插值计算, 估计误 差 P.42 第3题（线性插值用x=0.5, 0.6两点；二次插 值用x=0.5, 0.6, 0.7三点）,估计误差 P.42 第4题 第6题 第8题 p.43 第19题 * * * * 性质4：设f(x)在[a,b] 存在 n 阶导数，且 ，则  使得: 例. f(x)=3x4-5x3+4x-8, 则 f[0,1,2,3,4] f[1,2,3,4,5,6] =f(4)/4! =3, = f(5)/5! =0 差商的计算 一阶差商 二阶差商 三阶差商 四阶差商 * 差 商 表 已知 计算三阶差商 解：列表计算 算例 * § 2.3.2 牛顿插值公式 根据差商的定义，把 x 看成[a,b]上的一点， 可得： * 牛顿插值公式 把后一式 代入前一式 根据差商的定义，把 x 看成[a,b]上的一点， 显然Nn(x)满足插值条件，且次数不超过n,它就 是插值多项式，第i 次项系数为： 其中 我们称 为牛顿插值多项式. 已知 的函数表，求4 次牛顿插值多项式, 并求 算例 * 从表中可以看到4 阶差商几乎为0，故取4次插值多项式即可， 于是： 0.40 0.41075 0.55 0.57815 1.11600 0.65 0.69675 1.18600 0.28000 0.80 0.88811 1.27573 0.35893 0.19733 0.90 1.02652 1.38410 0.43348 0.21300 0.03134 1.05 1.25382 1.51533 0.52493 0.22863 0.03126 -0.00012 解：列表计算 已知f (x)的函数表，求4 次牛顿插值 多项式,并求f (0.596). 算例: 0.40 0.41075 0.55 0.57815 1.11600 0.65 0.69675 1.18600 0.28000 0.80 0.88811 1.27573 0.35893 0.19733 0.90 1.02652 1.38410 0.43348 0.21300 0.03134 1.05 1.25382 1.51533 0.52493 0.22863 0.03126 -0.00012 解：列表计算 已知 的函数表，求4 次牛顿插值多项式, 并求 算例 截断误差为： §2.4 差分与等距节点插值公式 设函数 f (x)在等距节点上的值 fi= f (xi)为已知， 这里h为常数，称为步长。 下面来讨论