Chap4Sec1 插值及数据拟合1n.pptVIP

杨东武 ydw_1978@126.com 数据插值 多项式插值问题的一般提法 插值多项式的唯一性 拉格朗日插值的基函数构造法 插值余项 分段插值 数据插值 多项式插值问题的一般提法 插值多项式的唯一性 拉格朗日插值的基函数构造法 插值余项 分段插值 数据插值 多项式插值问题的一般提法 插值多项式的唯一性 拉格朗日插值的基函数构造法 插值余项 分段插值 数据插值 多项式插值问题的一般提法 插值多项式的唯一性 拉格朗日插值的基函数构造法 插值余项 分段插值 数据插值 多项式插值问题的一般提法 插值多项式的唯一性 拉格朗日插值的基函数构造法 插值余项 分段插值 数据插值 多项式插值问题的一般提法 插值多项式的唯一性 拉格朗日插值的基函数构造法 插值余项 分段插值 * 机电工程学院 上节课内容回顾 2. 试写出埃特金法、牛顿迭代法以及弦截法的迭代公式？ 3. 应用牛顿迭代法的前提条件有哪些？ 1. 求解非线性方程时迭代公式 的收敛条件是什么？ 4. 试写出计算 的牛顿法迭代公式？ 以近似计算函数值为例来说明 散点上的函数值，即已知函数表 例1：设在实际问题中，某些变量之间的函数关系是存在的，但通常不能用式子表示，只能由实验、观测得到 在一系列离 如何计算 ？ 我们希望 寻求一个简单且易于计算的函数 来近 似 ，即 ，一般 可选为多 多项式、三角多项式、有理函数或样条函数等。 例2：有些函数虽有表达式，但较复杂，计算函数值不经济，这时也希望用简单的函数来逼近。 当精确函数 y = f(x) 非常复杂或未知时，在一系列节点 x0 … xn 处测得函数值 y0 = f(x0), … yn = f(xn)，由此构造一个简单易算的近似函数 g(x) ? f(x)，满足条件g(xi) = f(xi) (i = 0, … n)。这里的 g(x) 称为f(x) 的插值函数。最常用的插值函数是 …? 多项式。 x0 x1 x2 x3 x4 x g(x) ? f(x) §4.1多项式插值问题的一般提法 §4.2 拉格朗日(Lagrange)插值 n i y x P i i n , ... , 0 , ) ( = = 求 n 次多项式 使得 条件：无重合节点，即 思考题： 这样的n次多项式能构造出来吗？如何去做？ n i y x P i i n , ... , 0 , ) ( = = 求 n 次多项式 使得 注:一次多项式插值 --- 过两点直线。 二次多项式插值 --- 过三点抛物线。 若不将多项式次数限制为 n ，则插值多项式不唯一。 思考：通过解线性方程组求取插值多项式系数的方法是否是最好的？ 每改动一个节点的函数值（节点位置）是否就得重新求解一次线性方程组吗？ n = 1 已知 x0 , x1 ; y0 , y1 ，求 使得 1 1 1 0 0 1 ) ( , ) ( y x P y x P = = 可见 P1(x) 是过 ( x0 , y0 ) 和 ( x1, y1 ) 两点的直线。 ) ( ) ( 0 0 1 0 1 0 1 x x x x y y y x P - - - + = 1 0 1 x x x x - - 0 1 0 x x x x - - = y0 + y1 l0(x) l1(x) ? = = 1 0 ) ( i i i y x l 称为拉格朗日插值基函数 ， 满足条件 li(xj)=?ij /* Kronecker Delta */ 二. 拉格朗日插值的基函数构造法 n ? 1 希望找到li(x)，i = 0, …, n 使得 li(xj)=?ij ；然后令 ? = = n i i i n y x l x P 0 ) ( ) ( ，则显然有Pn(xi) = yi 。 ? = - = - - - = n j j ? i j i n i i i x x C x x x x x x C x l 0 0 ) ( ) )...( )...( ( ) ( 拉格朗日插值多项式，常记为Ln(x) 拉格朗日插值基函数 与 有关，而与 无关 节点 f 是n次多项式。 除外 除外 li(x) 每个 li 有 n 个根 x0 … xi … xn， ? - = = j=0 j ? i j i i i i x x C x l ) ( 1 1 ) ( n 插值余项 设节点 在[a , b]内存在, 考察截断误差