计算方法—插值法解释.ppt

第二章 插 值 法 （ Interpolation） 2.1 引言 2.2 拉格朗日插值 2.3 均差与牛顿插值公式 2.5 分段低次插值 2.6 三次样条插值 2.4 埃尔米特插值 Chapter2 插值法 Chapter2 插值法 表示两个变量x,y内在关系一般由函数式y=f(x)表达。但在实际问题中的函数是多种多样的,有下面两种情况： （1）由实验观测而得的一组离散数据(函数表)， 显然这种函数关系式y=f(x)存在且连续，但未知。 （2）函数解析表达式已知，但计算复杂，不便使用。通常也列函数表，如y=sin(x),y=lg(x) 2.1 引言 办法是：根据所给的y=f(x)的函数表，构造一个简单的连续函数P(X)近似替代f(x)。 Chapter2 插值法 近似代替即逼近的方法有很多种：插值方法、最佳一致逼近、最佳平方逼近、曲线拟合。 由于问题的复杂性,直接研究函数f(x)可能 很困难,但为了研究函数的变化规律,有时要求 不在表上的函数值，怎么办？ 简单连续函数P(x)指可用四则运算计算的函数: 如有理函数(分式函数) ,多项式或分段多项式。 Chapter2 插值法 插值问题的数学提法： 已知函数y=f(x)在n+1个点x0,x1, …,xn上的函数值yi=f(xi)(i=0,1, …,n)，求一个简单函数y=P(x)，使其满足P(xi)=yi,（i=0,1, …,n）。即要求该简单函数曲线要经过y=f(x)上已知的这n+1个点(x0,y0),(x1,y1), …,(xn,yn),同时在其它x∈[a,b]上要估计误差R(x)=f(x)-P(x)。 Chapter2 插值法 重要术语 对于n+1个基点的插值问题，我们称： f(x) 为被插值函数； P(x)为插值函数； x0,x1,…,xn为插值基点或插值节点； P(xk)=f(xk),k=0,1,…,n为插值条件； [a,b]为插值区间。 注释：对于早期的插值问题来说，f(x)通常是已知的，比如对数函数，指数函数，三角函数等，这些问题现在已经不用插值法来计算了；对于现在的许多实际问题来说，我们并不知道f(x)的具体形式，所对应的函数值可能是由测量仪器或其他物理设备中直接读出来的，f(x)只是一个概念中的函数 。 Chapter2 插值法 多项式插值 对于n+1个基点的插值问题，如果要求插值函数是次数不超过n的多项式，记为Pn(x)，则相应的问题就是多项式插值，并且把Pn(x)称为插值多项式。 实际上，我们所考虑的插值函数通常都是多项式函数或分段多项式函数。由于次数不超过n的多项式的一般形式为： Pn(x)=a0+a1x+a2x2+…+anxn 所以只要确定了n+1个系数a0，a1， …，a2，an，我们便确定了一个插值多项式。 Chapter2 插值法 y x x0 x1 x2 … xn-1 xn y0 y1 y2 多项式插值的几何意义： 多项式Pn(x) ，其几何曲线过给定的y=f(x)的n+1个点(xi,yi) i=0,1,2,…,n。 yn yn-1 Chapter2 插值法 Chapter2 插值法 2-1 插值多项式的唯一性 已知y=f(x)的函数表， 且xi(i=0,1,…,n)两两互异，xi∈ [a,b]。 求次数不超过n的多项式 使得Pn(xi)=yi,i=0,1,2…,n 此问题中Pn(x)是否存在？存在是否唯一？如何求？ 显然关键是确定多项式Pn(x)的系数a0,a1, …,an。 Chapter2 插值法 2-1 插值多项式的唯一性 定理：在n+1个互异的插值节点x0,x1,…,xn 上满足插值条件Pn(xi)=yi,i=0,1,2…,n的次数不超过n的代数多项式Pn(x)存在且唯一。 分析：为求 主要考虑插值条件 Chapter2 插值法 2-1 插值多项式的唯一性 证明：由插值条件，有 其系数矩阵的行列式为 关于未知量a0,a1, …,an的 非齐次线性方程组 Chapter2 插值法 例 给定f(x)的函数表，求f(x)的次数不超过3的插值多项式。 解：设 则， 解方程组得a0=10,a1=5,a2=-10,a3=2 即P3(x)=10+5x-10x2+2x3 当n=20,在109次/秒的计算机上计算需几万年！ Chapter2 插值法 2-2 线性插值与抛物插值 问题的提法： 已知函数y=f(x)的函数表 求次数不超过1的多项式L1(x)=a0+a1x 满足插值条件L1(xk)=yk, L1(xk+1)=yk+1。 分析