数学建模之第五讲多项式插值.pptVIP

附：差商性质 (2) 差商具有对称性,即任意调换节点的次序,差商的值不变 解：构造差商表计算差商 例：已知 求三次牛顿插值多项式近似f(x) 12 15 2 0 7 4 3 1 12 7 -1 15 4 -3.5 13 2 3 -1.25 4 1 0 1 三阶差商 二阶差商 一阶差商 则三次Newton插值多项式为 五、Newton差分型插值多项式 差分定义 依此类推 可以证明 如 差分表 在等距节点的前提下,差商与差分有如下关系 依此类推 由差商与向前差分的关系 如果假设 Newton向前(差分)插值公式 则插值公式 化为 其余项 化为 称 为Newton向前插值公式 插值余项为 插值余项为 根据向前差分和向后差分的关系 如果假设 可得Newton向后插值公式 Newton向后(差分)插值公式 §4 Hermite插值 一、问题的提出 不少实际问题不但要求在节点上函数值相等，而且还要求它的导数值也相等（即要求在节点上具有一阶光滑度），甚至要求高阶导数也相等，满足这种要求的插值多项式就是埃尔米特（Hermite）插值多项式。 下面只讨论函数值与导数值个数相等的情况。 称为二重密切Hermite插值 已知 求一个2n+1次多项式H(x)近似未知函数f，要求H(x)满足 以两个插值结点为例： 已知 求一个3次多项式H3(x)近似未知函数f，要求H3(x)满足 这样的H称为三次Hermite插值多项式。 从问题提出过程可以看出，这里给出了未知函数的四个条件，原本可以利用这4个插值条件确定三次Hermite插值多项式的系数。 二、三次Hermite插值多项式的确定 这里我们采用求拉格朗日插值多项式的基函数方法加以确定。 设 其中， 都是三次多项式 为了保证H3(x)满足插值条件，这四个基函数只需满足下列条件 以 的推导为例说明 设 只需确定四个基函数，则三次Hermite插值多项式也就确定了 可以求出 所以， 类似地可以求出其它三个基函数。 * * Chap4 多项式插值 Polynomial interpolation §1 引言 问题的提出 在实际问题中常遇到这样的函数 ,其在某个区间 上是存在的。但是,通过观察或测量或试验只能得到在 区间上有限个离散点 上的函数值 或者 的函数表达式是已知的,但却很复杂 而不便于计算,希望用一个简单的函数来描述它。 最常用的简单函数是多项式函数,这是因为多项式函数便于计算,同时理论上多项式函数有如下良好性质 Weierstrass逼近定理 0 a b y=f(x) y=P(x) y=f(x)+e y=f(x)-e 美国的人口普查每10年举行一次，下表列出了从1940年到1990年的人口（按千人计） 249633 226542 203302 179323 151326 132165 人口（千人） 1990 1980 1970 1960 1950 1940 年 能否利用这些数据合理地估计人口的数量，比如1965年的人口，甚至2010年的人口。 两个例子 假定某地区某天的气温变化记录如下： 16 17 18 20 22 24 25 27 29 31 32 31 温度 24 23 22 21 20 19 18 17 16 15 14 13 时刻 28 25 23 22 20 18 16 15 14 14 14 14 15 温度 12 11 10 9 8 7 6 5 4 3 2 1 0 时刻 找出这一天气温变化规律 插值问题的数学提法 在n+1个互异结点处的函数值 已知函数 求一个n次多项式 使其满足 从而可以用P(x)近似f(x) 插值条件的意义？ 多项式P(x)称为插值多项式. 多项式插值的几何意义 求经过n+1个点的n次多项式曲线y=P(x), 用来近似未知或复杂函数 x0 x1 x2 x3 x4 x P(x) ? f(x) x0 x1 x2 xn 例如，n=1时求一个一次多项式P1(x)，直线使其经过已知的两点 y x 当n=2时，求一个二次多项式P2(x)（抛物线），使其通过已知的3个点 y x 问题 满足插值条件的插值多项式是否存在?是否唯一? 如何构造?如何估计误差? 插值多项式的存在唯一性 设n次多项式 满足插值条件 则应有下面的方程组 要证明插值多项式存在唯一，只要证明 存在且唯一， 即只要证明该方程组存在唯一解。只需证明其系数矩阵 行列式不为零即可。 该行列式为范德蒙行列式。利用行列式性质可得 故插值多项式存在唯一 §2 Lagrange插值多项式 一、插值基函数 取基函数为 它满足 注意：有几个结点就有几个基函数 二、Lagrange插值多项式 根据基函数的性质，容易证明Ln(x)满足插值条件 特别地，两点一次插值（线性插值）多项式为： 三点二次插值（抛物线