插值与拟合科学与工程计算(数值分析)精要.ppt

科学和工程计算 第2章 插值法 插值法 插值法是一种古老的数学方法，早在一千多年前的隋唐时期定制历法时就广泛应用了二次插值。刘焯将等距节点的二次插值应用于天文计算。 插值理论却是在17世纪微积分产生后才逐步发展起来的，Newton插值公式理论是当时的重要成果。 由于计算机的使用以及航空、造船、精密仪器的加工，插值法在理论和实践上都得到进一步发展，获得了广泛的应用。 插值法 引言 拉格朗日插值 均差与牛顿插值公式 差分与等距节点插值 埃尔米特插值 分段低次插值 三次样条插值 科学和工程计算 第2章 插值和拟合 Lagrange插值多项式的缺点 解决 小结 曲线拟合 早期工程师制图时，把富有弹性的细长木条（所谓样条）用压铁固定在样点上，在其它地方让它自由弯曲，然后画下长条的曲线，称为样条曲线。它实际上是由分段三次曲线并接而成，在连接点即样点上要求二阶导数连续，从数学上加以概括就得到数学样条这一概念。 三次样条插值 样条本质上是一段一段的三次多项式拼合而成的曲线 在拼接处,不仅函数是连续的,且一阶和二阶导数也是连续的 一、三次样条插值函数 定义1. ------(1) 注：三次样条与分段 Hermite 插值的根本区别在于S(x)自身光滑，不需要知道 f 的导数值（除了在2个端点可能需要）；而Hermite插值依赖于f 在所有插值点的导数值。 f(x) H(x) S(x) 要求出S(x)，则在每个小区间上 要确定4个待定系数，共有n个小区间，所以应确定4n个参数。 共(n+1)+(3n-3)=4n-2个条件，因此还需要两个条件才能确定S(x) 可在区间端点a,b上各加一个条件（边界条件），具体要根据实际问题要求给定； 已知两端的一阶导数值 2. 两端的二阶导数已知 其特殊情况为 3. 当f(x)是 为周期的周期函数时，则要求S(x)也是周期函数，这时边界条件应满足： 这样确定的样条函数S(x)称为周期样条函数； 加上任何一类边界条件(至少两个)后 一般使用第一、二类边界条件, 常用第二类边界条件 样条插值函数的建立 即 或 可直接利用分段三次Hermit插值，只要假定 可得 加以整理后可得 ------(10) 由条件 ------(11) 由于以上两式相等,得 基本方程组 如果问题要求满足第一类(一阶)边界条件: 即 ------(12) 基本方程组化为n-1阶方程组 将上式化为矩阵形式 ------(13) ------(14) 这是一个三对角方程组 如果问题要求满足第二类(二阶自然)边界条件: 由(11)式,可知 ------(15) ----(16) ------(17) ------(18) 与基本方程组(12)联合,并化为矩阵形式,得 -----(19) (19)式与(14)一样,都是三对角方程组,解是唯一的； 例1. 对于给定的节点及函数值 解: 由(12)式可得 由(19)式得基本方程组 将上述结果代入(10)式 定理 . 最后，介绍一个有用的结果 当函数只在有限点集上给定函数值，要在包含该点击的区间上用公式给出函数的简单表达式，这些都涉及到在区间[a,b]上用简单函数逼近已知复杂函数的问题，这就是函数逼近问题。 插值法就是函数逼近问题的一种 拟解决的问题： 计算复杂的函数值 已知有限点集上的函数值，给出在包含该点集的区间上函数的简单表达式 函数逼近——对函数类A中给定的函数f(x)，记作 要求在另一类简单的便于计算的函数类B中求函数 使p(x)与f(x)的误差在某种度量意义下最小。 逼近问题 函数逼近 曲线拟合 实例：考察某种纤维的强度与其拉伸倍数的关系,下表 是实际测定的24个纤维样品的强度与相应的拉伸倍数 是记录: 纤维强度随拉伸 倍数增加而增加 并且24个点大致分 布在一条直线附近 必须找到一种度量标准来衡量什么曲线最接近所有数据点 (1) 仍然是已知 x1 … xm ; y1 … ym, 求一个简单易算的近似函数 P(x) ? f(x)。 但是 ① m 很大； ② yi 本身是测量值，不准确，即 yi ? f (xi) 这时没必要取 P(xi) = yi , 而要使 P(xi) ? yi 总体上尽可