数值计算方法第二章程序.ppt

第二章 插值法 本章主要内容 * 1、插值问题 2、拉格朗日插值 3、牛顿插值 4、埃尔米特插值 5、三次样条插值 第一节 插值问题 问题的提出 * 1 函数表达式过于复杂不便于计算, 而又需要计算许多点处的函数值 2 仅有几个采样点处的函数值, 而又需要知道非采样点处的函数值 □上述问题的一种解决思路：建立复杂函数或者未知函数的一个便于计算的近似表达式. □解决方法－插值法 插值的数学提法 * 实际问题中常常需求函数 的函数值、导数值、零点、值或积分值。但往往不知道其确切表达式，或表达式很复杂，只是知道它在某些点处的函数值与导数值。本章介绍两种方法（插值法与最小二乘法）求其近似式 在众多函数中,多项式最简单、最易计算，已知函数 个互不相同的点处的函数值 ，为求 的近似式，自然应当选 次多项式 使 满足条件 插值的几何意义 * 需要研究的几个问题 满足条件的插值多项式是否存在，是否唯一 如何构造满足条件的插值多项式 用插值多项式 代替 误差如何估计 * 存在唯一性 * 设插值多项式为 带入条件 系数行列式为范德蒙行列式 只有唯一解 对于次数不大于n 的多项式 n次插值多项式就是其本身 误差估计 * 误差估计 * 由公式可见， 如果 在 上有界， 即存在常数 使 ，则必有 拉格朗日多项式插值 * 一、线性插值 当n=1时，已知 设 代入已知条件，得方程组 记 ＋ = 例1.由常用的正弦函数特殊值，利用线性插值求 * 解： 取 则 所以 ※真值为：0.5878 二、抛物插值 * 当n=2时，已知 设 带入3个节点，得到 ＋ ＋ 拉格朗日插值 * 1、插值基函数 定义 如果 为区间 的一组节点，称 次多项式： = = 为插值基函数 * 插值基函数性质 拉格朗日插值多项式 令 * 即 可得 其中 * 例1: 解: * 且 在例1中,如果只给出两个节点169和225,也可以作插值多项式,即1次Lagrange插值多项式,有两个插值基函数, 这种插值方法称为Lagrange线性插值,也可以在n+1个 节点中取相邻的两个节点作线性插值 * 牛顿插值 Lagrange插值多项式的缺点： 插值基函数计算复杂 高次插值的精度不一定高 差商(均差 定义 称 * 依此类推 显然 规定 为 在 处的零阶差商 由微商的定义可知差商是微商的离散形式 * 差商具有如下性质： (2) 差商具有对称性,即任意调换节点的次序,差商的值不变 （3） f（x)为n次多项式，则 为n-1次多项式 * 差商的表格计算方法： * * * * *