数字分析模型ppt.ppt

数值分析模型 第一节 插值法 在生产和科学研究中,经常出现这样的问题:由实验或测量得到的某一函数 在一系列点 处的值 ,需要构造一个简单函数 作为函数 的近似表达式: ,使得 这类问题称为插值问题. 称为被插值函数, 称为插值函数 称为插值节点;式(6-1)称为插值条件. 常用的插值函数是多项式与样条函数. 拉格朗日(lagrange)插值 取n次多项式pn(x)作为插值函数 pn(x)=a0+a1x+a2x2+…+anxn 利用插值条件有： 其系数行列式为n+1阶范德蒙行列式，因插值节点 互不相同，所以方程组的解存在且唯一。 其系数行列式为范德蒙(Vandermonde)行列式: 由于插值节点 互不相同,所以上述行列式不等于0,故由克莱姆(Cramer)法则知,方程组(6-3)的解存在而且是唯一的. 实际上比较简单的方法不是解方程组(6-3),而是构造一组插值基函数.为此,首先求满足条件 的n次多项式 因为式（6-4）表明，除 点以外，其他所有的节点都是n次多项式 的零点，故设 其中A为待定常数。由 可得 所以 称之为拉格朗日插值基函数。 利用插值基函数（6-5），可以构造多项式 就是满足插值条件 的拉格郎日插值问题的解，称式（6-6）为拉格朗日插值多项式。 特别地，当n=1时称为线性插值，其插值多项式为： 满足 从几何上看, 为过两点 的直线。 当n=2时，称为抛物线插值，其插值多项式为： 满足 从几何上看 为过点 和 的一条抛物线。 埃尔米特插值 许多插值问题不但要求在节点上函数值相等，而且还要求对应的导数值也相等，甚至要求高阶导数也相等，满足这种要求的插值多项式被称为埃尔米特（Hermite）插值多项式设在节点 上， 要求插值多项式 ，满足条件 由于（6-7）式给出了2n+2个条件，因此可以唯一确定一个次数不超过2n+1的多项式 ，其形式为 根据（6-7）式来确定 显然非常复杂。 仿照拉格朗日插值多项式的基函数方法，可先求插值基函数 ，及 共用2n+2个，每一个基函数都是2n-1次多项式，且满足条件 于是满足条件（6-7）的插值多项式 可写成 由条件（6-8）式显然有 利用拉格朗日插值基函数 令 其中 为（6-5）式所表示的基函数。 由条件（6-8）式可得 整理得: 解出 对 两边取对数求导可得 于是 同理 仿照拉格朗日插值余项的证明方法，若 在 内的2n+2阶导数存在，则其插值余项为 其中 三次样条插值 分段线性插值，具有良好的稳定性和收敛性，但光滑性较差。在数学上若函数（曲线）的k阶导数存在且连续，则称该曲线具有k阶光滑性。易见，分段线性插值不光滑，这影响了它在某些工程技术实际问题中的应用。例如：在船体、飞机等外形曲线的设计中，不仅要求曲线连续而且还要求曲线的曲率连续，这就要求插值函数具有连续的二阶导数。为解决这一类问题，就产生了三次样条插值。 所谓样条（Spline），本来是指一种绘图工具，它是一种富有弹性的细长木条，在飞机或轮船制造过程中，被用于描绘光滑的外形曲线。使用时，用压铁将其固定在一些给定的节点上，在其他地方任其自然弯曲，然后依样画下的光滑曲线，就称为样条曲线。它实际上是由分段三次曲线拼接而成，在连续点即节点上，不仅函数自身是连续的，而且它的一阶和二阶导数也是连续的。从数学上加以概括，可得到样条函数的定义如下： 三次样条函数记作S（x），a≤x≤b，满足： ①在每个小区间[xi，xi+