差值与拟合方法习题.doc

第五章 差值与拟合方法 的n次插值多项式，则有对任何，存在 使 （5.7） 成立。式中。 在定理5.3中，由于f（n+1) (()一般不知道，插值误差不能算出，但若能计算出|f（n+1) (x)|在[a,b]上的最大值M n+1，则可得余项估计式 （5.10） 利用式（5.10）可以估计插值函数在某一点的误差，若还想估计函数在插值区间[a,b]上的误差，可计算出 这样，可得在整个插值区间[a,b]上的误差估计 （5.11） 由n次插值多项式的唯一性及式（5.7），则有如下结论： 定理5.4 若函数f ( x )在[a,b]上有n+1阶导数，则f ( x )有如下展式： 对n=1的插值多项式，称为线性插值，n=2的插值多项式称为抛物线插值。 【例5.1】 已知的函数表为 x 3.0 3.1 3.2 3.3 3.4 y=f(x) 1.098612 1.131402 1.163151 1.193922 1.223775 试用线性插值和抛物线插值分别计算的近似值，并估计相应的误差。 解 线性插值需要两个节点，内插比外推好，因为，故选，由的Lagrange 插值公式，有 所以有 为保证内插，对抛物线插值，选取三个节点为由n=2的Lagrange 插值公式，有 故有。 为考虑误差，由于当时，有，所以线性插值计算的误差估计为 而当时，，故抛物线插值计算的误差估计为 可见抛物线插值比线性插值精确，这在一般情况下也是成立的。 【例5.2】在上给出的等距节点函数表，若想用二次插值来计算的近似值，并要求截断误差不超过，问此函数表的步长h应为多少？计算时要取多少位有效数字？ 解 设为上的等距节点，要用二次插值应取三个相邻的节点。设是上的任何三个相邻节点，则当时，有 注意到并利用n=2的Lagrange余项定理有函数在的插值余项为 因为 ，且 所以 要，只需，于是得到 取h=0.0057可满足要求，由得，故造表时应在上取1404个等距节点来计算函数值，至少取7位有效数字。 5.3.2 Newton 插值 定义 5.1 已知函数f(x)在n+1个互异节点上的值分别为 记 式中 是中互不相同的k个节点，+1 称为f(x)关于节点的k-1阶差商。 规定零阶差商是函数值本身：，于是一阶差商可写为 利用差商的定义，可以得到有关差商的性质： 性质1 差商可表示为节点处函数值的线性组合：即 式中 证明：用数学归纳法。 当k=1时，，成立。 假设当k=n-1时成立，即有 则当k=n时，由定义 性质3：若在含节点的区间[a，b]上n阶连续可微，则至少存在一点，使得 表5.1 差商表 x f（x） 一阶差商 二阶差商 …… n阶差商 . . . . . . . . ?? 当插值节点是等距节点是，即，则Newton插值公式可以演变成为两个公式：Newton前插公式和Newton后插公式，其用于等距节点插值时更为简单，这两个公式可利用插值节点等距的特点从Newton公式简化得到，具体形式如下： ①Newton前插公式 式中。 称为向前差分。 ②Newton后插公式 式中 ， 称向后差分。 通常对等距插值问题，当要计算的x值靠近时，用Newton前插公式，而当x靠近时，用Newton后插值公式。这两个公式也可用相应的差分表计算相应的差分值。 【例5.3】给定数表： x 1 2 4 5 6 8 f(x) 0 2 8 12 18 28 试分别用二次和四次Newton插值多项式计算的近似值。 解 作差商表： x f(x) 一阶差商 二阶差商 三阶差商 四阶差商 五阶差商 1 0 2 1/3 0 1/30 -1/60 2 2 3 1/3 1/6 -1/12 4 8 4 1 -1/3 5 12 6 -1/3 6 18 5 8 28 用二阶Newton插值近似计算，应选与5.8最近的3个节点，由表中行数据有 所以 要求来计算，取与5.8最近的5个节点 ，由表中行数据有 所以。 Hermite插值 先看如下例子： 【例5.4】设有四阶导数，且， 求函数的一个插值多项式，并用此近似函数来计算的近似值，请尝试给出你所得插值多项式的误差关系式，估计近似计算的误差。 解：仿照Lag