.插值与逼近.pptVIP

Ch2 插值与逼近 一、代数插值 内 容 插值要解决的问题 代数插值的概念 Lagrange插值 Newton插值 1.插值要解决的问题 2.代数插值的概念 存在唯一性 常用的代数插值方法 拉格朗日（Lagrange） 牛顿（Newton） 埃特金（Aitken） 埃尔米特（Hermite） …… 3.Lagrange插值 一. 线性插值 已知f(x)两互异节点处的函数值，求f(x)的一次 插值多项式P1(x)。 二. 抛物线插值 已知f(x)三个互异节点处的函数值，求f(x)的二 次插值多项式P2(x)。 三. n次Lagrange插值多项式 4.插值余项 例1 给定函数y=sinx的函数表 总 结 0.64422 0.56464 0.47943 0.38942 sinx 0.7 0.6 0.5 0.4 x 试分别用线性插值和二次插值求sin0.57891的近似值，并估计截断误差。 0.64422 0.56464 0.47943 0.38942 sinx 0.7 0.6 0.5 0.4 x 解： 函数值表 - 0.223144 - 0.356675 - 0.693147 -0.916291 lnx 0.8 0.7 0.5 0.4 x 例2 设 f(x)= lnx , 并假定已给出值表，试近似计算ln(0.6)的值，并指出精度。 解：利用3次Lagrange 插值公式, 简单计算过程如下： *上海理工大学理学院 早在1000多年以前，中国的历法上就已经记载了插值的思想和简单插值的应用。 引例1 男子竞赛奥运会纪录 是否能建立竞赛距离与纪录时间之间的 函数关系，并测算男子1000米纪录。 亚特兰大 挪威 罗达尔 1:42.58 800 1000 悉尼 肯尼亚 恩格尼 3:32.07 1500 亚特兰大 美国 约翰逊 43.49 400 北京 牙买加 博尔特 19.30 200 北京 牙买加 博尔特 9.69 100 散点图 引例2 2.5649 2.4849 2.3979 lnx 13 12 11 x 1.为离散的数据建立连续的数学模型； 2.用简单函数为复杂函数提供好的逼近。 即找到一个既能反应函数的变化规律， 又便于计算的简单函数p(x)来近似代替f(x)。 插值要解决的问题: 反应函数的变化规律： p(xi) = f(xi) = yi， (i=0,1,2,……，n)， p(x)——插值函数，这种问题称为插值问题。 2.便于计算： p(x)的选择： 代数多项式，三角多项式，有理函数…… 若取插值函数p(x)为多项式，则称为代数插值。 假设已知函数 f(x) 在n+1个互异点xi 处的函数值 (或一组对应的离散数据)yi = f(xi) (i=0,1,2,……，n)， 满足条件Pn(xi) = f(xi) = yi ，(i=0,1,2,……,n)， 定义 设函数f(x)在区间[a, b]上有定义，且已知它在 n+1 个互异点 a?x0x1 …… xn ?b处的函数值 yi = f(xi) (i=0,1,2,……,n)，若存在一个次数不超过n次的多 项式 Pn(x)=a0+a1x+……+anxn 则称Pn(x) 为f(x)的n次代数插值多项式。 xi ：插值节点， f(x) ：被插值函数， Pn (xi) = f(xi) = yi ：插值条件，[a, b] ：插值区间 几何意义 问题 1.代数曲线Pn (x)是否存在，有几条? 2.怎样求Pn (x)？ 3.误差 x ：插值点， x在[a, b]内 ：内插， 否则为外插 满足条件Pn (xi) = yi = f(xi) 。 定理1 n次代数插值问题的解是存在且唯一的。 Pn(x)=a0+a1x+……+anxn 已知函数f(x)在n+1个互异节点a?x0x1 …… xn ?b 处的函数值yi = f(xi) (i=0,1,2,……，n)， 则存在唯一一个次数不超过n次的多项式： 证明：设所要构造的插值多项式为： 由插值条件 得到如下线性代数方程组： 此方程组的系数行列式为 范德蒙行列式 因为 所以D ? 0，由Cramer法则，方程组的解存在而且唯一， 因此，Pn(x)由a0, a1,…, an唯一确定。 基本思想——基函数法 把Pn(x)的构造问题转化为n+1个插值基函数（基本插值多项式）lk(x)(k=0,1,…,n)的构造。 y1 y0 y x1 x0 x 线