工程数学 第一章 计算方法.ppt

则由前面Hermite分段插值有: 设 , 在结点内部, 3次多项式有2阶连续导数, 故只要在内部结点处具有1, 2阶连续导数即可。 (1) 而对于式(1) , 已有 , 故一阶导数连续。 因此, 只要在内结点 上二阶导也连续, 即 (2) (2) 中包含了未知量 m0, m1, · · ·, mn-1, mn , 边界处 m0 , mn可由边界条件确定。 m1, · · ·, mn-1 可由(2)式的 n-1个方程确定。 下面对(2)进行推导。 m1, · · ·, mn-1 xk xk-1 xk+1 x 二个3次多项式, 先求出2阶导数, 再在 xk 处相等: 令 , 得 , 即 , 其中 , 如果已知边界条件(两端点)为: , , 则得关于 m1, · · ·, mn-1 的方程组 解出 m1, · · ·, mn-1 代入(1)式即得分段3次样条插值函数 。 推完了, 下面举个例子算算看。 [例4.13] 插值到此为止。下面讲拟合。 (本节完) 作业: 练习四 10 (此题可以选做) §4.6 曲线拟合的最小二乘法 巴尔末(J.J.Balmer，1825-1898) 特殊爱好: 数字游戏 职 业: 数学教师，瑞士某女子中学， 兼巴塞尔大学无薪讲师 “我能用公式把任意4个数字有规律地联系起来” 4 — Why? 信息量? 道家: 1生2，2生3，3生万物 4 —万物？全部信息量?! 4 条谱线的波长数据 实验: 氢光谱谱线的波长数据 巴尔末公式: 艰苦努力+天才 引子 1884年 很有规律, 能记住吗? Lagrange插值公式: 下面做些练习: 但必须记住: 满足n+1个插值条件的n次插值多项式是唯一存在的。 不管用待定系数法, 或套用Lagrange插值公式; 还是下一节要讨论的Newton插值公式, 得到的都是同一个多项式!!! [例4.7] , 求 解: [例4.8] 利用100, 121, 144的开方求 的近似值。 解: [例4.9] 解: 已知正弦往函数表 利用抛物插值求 的值, 并估计误差。 [例4.10] 解: 设 , 求其在 , , 插值节点的。 上作为 因过3点的2次插值多项式是唯一的, 因此 从 也可以证明得到。 [例4.11] 解: 设 , 试用Lagrange插值余项公式 写出以-1，0，1，2为插值节点的3 次插值多项式。 得 [例4.11] 解: 插值条件中有导数值， 此时插值余项? 设有函数表 , 求插值多项式， 并写出余项表达式。 设 则 由插值条件得: 得 则 余项： 类似Lagrange插值余项(只是把导数点也用上) 则 , 故可设 设 则 当 t = 0, 1, 2, x 时 1 当 t = 时 即 当 t = 时 得 有 代入有: (本节完) 作业: 练习四 1, 2, 3, 9 上机作业4.1: 编写Lagrange插值公式程序, 并自己做几个数值实例。 下一节讲 对应于点斜式的 插值公式的另一种构造法: Newton 插值公式。 作业第2题提示: 2.1 可设 f(x) = xk , 代入Lagrange插值公式, 此时余项为0. 2.2 可设 f(x) = (t-x)k ,代入Lagrange插值公式, 此时余项为0. 得到的公式中再取 t = x。 §4.3 Newton 插值公式与差商 1. 公式构造: Nn(x ) = ? 2. 余 项: Rn(x) = ? 问题重提: 已知 f (x) x f (x0) x0 x1 f (x1) · · · · · · xn f (xn) 求 n 次多项式 使满足: Nn(x) 。 (1) , i = 0,1, · · ·, n Nn(xi ) = f (xi ) 同样问题, 仅仅是另一种不同表示: 同样结果, Nn(x) Ln(x) ≡ 什么是 ??? 差商 先相减? 后相除? √ 再提思想来源: n = 1: y = P1(x) 2点式 点斜式 f (x) x f (x0) x0 x1 f (x1) f (x) x f (x0) x0 x1 f (x1) xn f (xn) · · · · · · L1(x) N1(x) Lagrange 插值公式 Ln(x) Newton 插值公式 Nn(x) 待定系数法 ? 唯一性 存在性 Rn(x) 余项 看来3种方法都有用。 1. 公式构造: Nn(x ) = ? 推导思想: 递增式构造, 故具有承袭性。 n = 0: 0次多项式 n = 1: 设 则 由 得