现代科学与工程计算-孟大志-第一篇1 2 1.pptVIP

* * * 鲁东大学数 学 与 信 息 学 院 ? 插值法需解决的基本问题 ? 1、插值多项式是否存在？若存在，是否唯一？ 2、若插值多项式存在，应如何构造？ 3、构造的插值多项式与原函数之间的误差有多大？ 1.2 一般多项式插值 1.2.1 一般插值法 定义1.2.1 设某个函数y = f (x) 在区间 [a, b] 上是连续的， 且已知区间 [a, b] 上n +1个互异点处的函数值，即 y i = f ( xi ), i =0, 1, …, n 寻找一个次数不高于n的多项式p( x )使得 p( xi ) = y i , i =0, 1, …, n (1.2.1) 则称p( x ) 为f ( x )的插值多项式，其中xi ( i =0, 1, …, n)称为 插值节点，关系式(1.2.1)称为插值条件。 鲁东大学数 学 与 信 息 学 院 图1.2.1 多项式插值 1．多项式插值问题解的存在唯一性 设n阶插值多项式为 (1.2.2) 由插值条件(1.2.1)得到： 鲁东大学数 学 与 信 息 学 院 (1.2.3) 其系数矩阵的行列式为范德蒙行列式，并可以算出： 这是关于未知量 的线性方程组。由线性代数知识知， 由于xi ≠ xj，所以 。 鲁东大学数 学 与 信 息 学 院 由克莱姆法则知，方程组存在唯一解 即插值多项式p(x)是存在且唯一的。 ， 2．一般插值多项式的构造方法 用于插值的简单函数元素集 + 线性组合结构 → 插值多项式 简单函数元素集是指构成多项式的基函数集合，例如 自然形式（1.2.2）的自然基底 , 、 、 Lagrange型插值多项式的Lagrange型基底 Newton型插值多项式的Newton型基底 Hermite型插值多项式的Hermite型基底 基底不同的插值多项式的表达形式和构造方法也不同。 鲁东大学数 学 与 信 息 学 院 插值多项式的存在唯一性说明，满足插值条件的 多项式存在，并且插值多项式与构造方法无关。 如何构造插值函数才能达到预期的效果呢？ (结构) (集合) 1.2.2 Lagrange插值 已知条件：设某个函数y = f (x) 在区间 [a, b] 上是连续的，且已知 区间 [a, b] 上n +1个互异节点处的函数值，即 yi = f ( xi ), i =0, 1, …, n 目的：构造一个不超过n次的Lagrange插值多项式Ln(x) ,使得 在插值节点处满足条件 1．插值多项式的构造 构造思想： 若求自然形式(1.2.2)的插值多项式问题，只要求 解线性方程组（1.2.3）计算出多项式系数即可。 鲁东大学数 学 与 信 息 学 院 但是求解线性方程组的计算量较大，不适合大型数据的插值构造。 因此希望寻找一种计算量少的确定系数的方法。 则（1.2.2）式的多项式是这组基函数的线性组合。能否 构造另外一组插值基函数，使得构造多项式的线性组合 时不需要计算组合系数，组合系数就取已经给出的函数 值y i，这就是插值Lagrange问题。 如果把自然基底 看作多项式的一组基函数， 因此Lagrange插值方法的核心是构造为了实现以上目的 的插值基函数。 鲁东大学数 学 与 信 息 学 院 插值目标转化: 构造插值多项式→插值基函数 (2) 构造方法： 定义1.2.2 由于Ln(x) 需要满足 Ln(xi) = yi, i=0,1, …, n 可构造函数lk(x)使其在插值节点上 (1.2.4) 称(1.2.4)式为关于节点 xi (i=0,1, …, n) 的Lagrange插值基函数。 成立，即可满足插值条件，其中δkj 为克罗内克δ函数 鲁东大学数 学 与 信 息 学 院 定义1.2.3 使用Lagrange插值的n+1个插值基函数 lk(x) 的线性组合来构造插值多项式，取 yi (i=0,1, …, n) 作为插值多项式的组合系数，使得 则称(1.2.5)式为Lagrange插值多项式。 由(1.2.4)式知 lk(x) 有n个实零点 故可设 由lk(xk)=1知 鲁东大学数 学 与 信 息 学 院 将A代入 得 可见Lagrange插值基函数只与插值节点有关而与函数f (x)无关。 2．线性插值与抛物插值 定义1.2.4 n = 1时的Lagrange插值多项式称为线性插值。 已知函数f (x)在插值节点 处的函数值分别为 由(1.2.6)式知线性插值基函数分别为 鲁东大学数 学 与 信 息