计算方法B 第1章 插值.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * Home Work 推导三次样条函数插值在周期边界条件下的m和M关系式的边界条件 P40:12,15 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 误差 性质3 差商性质总结 证明作为作业 有关上机作业 提交计算机作业的ftp是00， 可以直接在“我的电脑”里输入00登录， 也可以用flashfxp等客户端软件登录，登录不需要输入用户名密码，匿名登录即可。 ftp的权限只能提交作业，而不能下载或删除作业， 如果他们发现自己提交的作业有问题，直接重新提交一份即可。 有关上机作业2 为了避免病毒起见，ftp我设置只能提交.rar, .zip, .7z三种压缩包格式，其余的文档类型不能提交，所以他们每个人的作业必须打包压缩后才能提交。 每个人作业的文件名最好包含姓名和学号。 有关上机作业3 我在ftp里设置了“第一次作业”、“第二次作业”等目录，每次作业放在相应的文件夹内。 最好让所有同学不论是否是否已通过邮件发送前两次作业，都往ftp提交作业，邮件每天太多，比较容易漏看，ftp则一目了然，比较清楚。 1.4 Hermite插值 有时候，构造插值函数除了函数值的条件以外，还需要一定的 连续性条件，如一阶导数值等，这种插值称为Hermite插值。 称为二重密切Hermite插值 例：设 x0 ? x1 ? x2, 已知 f(x0)、 f(x1)、 f(x2) 和 f ’(x1), 求多项式 P(x) 满足 P(xi) = f (xi)，i = 0, 1, 2，且 P’(x1) = f ’(x1), 并估计误差。 模仿 Lagrange 多项式的思想，设 解：首先，P 的阶数 = 3 ? + = 2 1 3 ) ( ) ( ) ( ) ( ) ( = 0 i i i x h x1 f ’ x h x f x P ? h0(x) 有根 x1, x2，且 h0’(x1) = 0 ? x1 是重根。 ) ( ) ( ) ( 2 2 1 0 0 x x x x C x h - - = 又: h0(x0) = 1 ? C0 h2(x) 与h0(x) 完全类似。 其中 hi(xj) = ?ij , hi’(x1) = 0, (xi) = 0, ’(x1) = 1 ? h1 ? h1 h1(x) 有根 x0, x2 ? ) )( )( ( ) ( 2 0 1 x x x x B Ax x h - - + = 由余下条件 h1(x1) = 1 和 h1’(x1) = 0 可解。 (x) ? h1 有根 x0, x1, x2 ? ? h1 ) )( )( ( ) ( 2 1 0 1 x x x x x x C x - - - = ? h1 又: ’(x1) = 1 ? C1 可解。 与 Lagrange 分析完全类似 仿照Lagrange插值的做法，首先确定多项式插值空间的维数， 注意到，我们的条件共有2(n+1)个条件，所以，最高次数为2n+1 对二重密切Hermite插值 整个构造步骤如下： 1、确定多项式的最高项次数，就是函数空间的维数 2、假设一组基函数，列出插值多项式 3、列出基函数满足的公式（画表），求基函数 称为 构造基函数方法 误差分析 类似Lagrange插值的分析方法 二重密切Hermite插值误差 例：在[?5, 5]上考察 的Ln(x)。取 - 5 - 4 - 3 - 2 - 1 0 1 2 3 4 5 - 0.5 0 0.5 1 1.5 2 2.5 n 越大， 端点附近抖动 越大 Ln(x) ? f (x) ? 是否次数越高越好呢？ Lab02 Lagrange插值 对函数 构造插值，并求 插值节点取为： (1) (2) 对N=5，10，20，40比较以上两组节点的结果。 Chebyshev点 Home Work P40 4, 6 分段低阶插值 Runge现象 例： 等距节点构造10次Lagrange插值多项式 -0.90 0.04706 1.57872 -0.70 0.07547 -0.22620 -0.50 0.13793 0.25376 -0.30 0.30769 0.23535 1901年，Runge 等距高次插值，数值稳定性差，本身是病态的。 分段线性插值 每个小区间上，作线性插值 (1) (2) 在每个小区间上为一个不高于1次