《插值方法基本思想》PPT课件.ppt_第1页
《插值方法基本思想》PPT课件.ppt_第2页
《插值方法基本思想》PPT课件.ppt_第3页
《插值方法基本思想》PPT课件.ppt_第4页
《插值方法基本思想》PPT课件.ppt_第5页
已阅读5页,还剩48页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

插值法基本思路 张兴元2011年8月 一元多项式插值 教学内容插值问题插值问题求解方法 重点 线性插值二次插值n次插值分段线性插值Hermite插值分段三次Hermite插值样条插值函数 难点 要求掌握以上方法的原理及其在MATLAB中的实现方法 插值问题 1 提法已知n 1个节点 xj yj j 0 1 n 其中xj互不相同 不妨设a x0 x1 xn b 求任一插值点x xj 处的插值y xj yj 可以看成是由某个函数y g x 产生的 g的解析表达式可能十分复杂 或不存在封闭形式 也可以未知 2 求解的基本思路构造一个相对简单的函数y f x 使f x 通过全部节点 即f xj yj j 0 1 n 再用f x 计算插值 即y f x f x 称为插值函数 如果f x 为k次多项式 f x 就是插值多项式 此时插值为代数插值 如果f x 为有理函数 就是有理插值 如果f x 为三角函数 则为三角插值 多项式插值 线性插值 y f x 函数表 线性插值 两点式方程 Lagrange插值 是l0 x 和l1 x 的线性组合 基函数 线性插值 点斜式方程 均差 Newton插值 一阶均差的一般定义 线性插值 余项 两种不同的构造方式 Lagrange和Newton 效果一样吗 此处一样 两种不同的构造方式 Lagrange和Newton 可以推广到多个点吗 可以 多项式插值 二次插值 y f x 函数表 二次插值 Lagrange基函数方法 Lagrange插值 二次插值 Newton均差法 二阶均差 Newton插值 二次插值 余项 例1 已知 试利用插值法近似计算 解 有几位有效数字 多项式插值 n次插值 y f x 函数表 xi互不相同 存在吗 唯一吗 如何构造 n次插值 存在性 唯一性 存在且唯一 n次插值 插值多项式的构造 方法一 Lagrange型插值多项式 基函数 基函数的特点 Lagrange插值多项式 n次插值 插值多项式的构造 方法二 Newton型插值多项式 均差表或差商表 Newton插值多项式 n次插值 插值余项与事后误差估计 插值余项 其中 事后误差估计方法 误差 n次插值 示例 例2 基于5个点 k cos k k 0 1 2 3 4 1 构造f x cos x 的差商表 2 并用差商表找出牛顿插值多项式的系数 3 写出四次牛顿插值多项式N4 x 4 计算N4 2 5 解 第一步 明确插值点 xk yk 第二步 构造差商表 第三步 写出相应的牛顿插值多项式N4 x 第四步 计算近似值N4 2 5 多项式插值的震荡性质 用Lagrange插值多项式LN x 近似f x a x b 虽然随着节点个数的增加 LN x 的次数N更大 多数情况下误差 RN x 会变小 但是N增加时 LN x 的光滑性变坏 有时会出现很大的震荡 理论上 当N 时 在 a b 内并不能保证LN x 处处收敛于f x Runge给出了一个有名的例子 多项式插值的震荡性质 Runge给出了一个有名的例子 多项式插值的震荡性质 高次插值多项式的这些缺陷 促使人们转而寻求简单的低次多项式插值 多项式插值 分段线性插值 是线性函数 多项式插值 分段线性插值 分段线性插值函数为 余项估计为 多项式插值 分段线性插值 分段线性插值多项式L1 x 的图像上是连接各插值点的一条折线 如右图 y sin x 的插值逼近图形变化 特点 曲线的光滑性较差 在节点处有尖点 增加节点 减小步长 会改善效果 若f x 在 a b 上连续 则 多项式插值 Hermite插值 考虑只有两个节点的插值问题 如何选择基函数 多项式插值 Hermite插值 希望插值系数与Lagrange插值一样简单 假设 其中 多项式插值 Hermite插值 可知 由 可得 Lagrange插值基函数 类似可得 即 将以上结果代入 多项式插值 Hermite插值 多项式插值 Hermite插值 得两个节点的三次Hermite插值公式 多项式插值 Hermite插值的插值余项 两点三次Hermite插值的余项为 例3 多项式插值 Hermite插值的插值余项 解 作为多项式插值 三次已是较高的次数 次数再高就有可能发生Runge现象 因此 对有n 1节点的插值问题 我们可以使用分段两点三次Hermite插值 多项式插值 Hermite插值的插值余项 多项式插值 分段三次Hermite插值 可构造两点三次Hermite插值多项式 多项式插值 分段三次Hermite插值 其中 分段三次Hermite插值多项式 余项为 多项式插值 样条函数插值 分段插值的思想及优缺点1 思想 将图形分段 每段为一个低阶多项式Sk x 并在相邻点之间进行多项式插值 组成一个分段的多项式曲线 2 分类 1 分段线性插值优点 简单 缺点 连续但不光滑 曲率不连续变化 2 分段二次多项式插值优点 简单 缺点 偶数点x2k处曲率变化很大 曲率不连续变化 3 改进方法 利用分段三次样条插值 分段三次多项式 连续 光滑 曲率连续变化 多项式的次数较低 多项式插值 样条函数插值 什么是样条 是指飞机或轮船等的制造过程中为描绘出光滑的外形曲线 放样 所用的工具 样条本质上是一段一段的三次多项式拼合而成的曲线 在拼接处 不仅函数是连续的 且一阶和二阶导数也是连续的 1946年 Schoenberg将样条引入数学 即所谓的样条函数 多项式插值 样条函数插值 1 三次样条插值函数的定义 多项式插值 样条函数插值 2 确定三次样条插值函数的条件 多项式插值 样条函数插值 多项式插值 样条函数插值 3 三次样条插值函数的构造方法 3 1用节点处一阶导数表示的三次样条插值函数 多项式插值 样条函数插值 多项式插值 样条函数插值 多项式插值 样条函数插值 整理后 得到 引入记号 则有方程组 多项式插值 样条函数插值 该方程组为三对角方程组 可以利用追赶法求解 多项式插值 样条函数插值 多项式插值 样条函数插值 多项式插值 样条函数插值 3 2用节点处二阶导数表示的三次

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论