4插值与拟合方法课件-1.docVIP

第4章 插值与拟合?方法 插值与拟合?方法是用有?限个函数值?去推断或表?示函数的方?法，它在理论数?学中提到的?不多。 本章主要介?绍有关解决?这类问题的?理论和方法?，涉及的内容?有多项式插?值，分段插值及?曲线拟合等?。对应的方法?有Lagr?ange插?值，Newto?n插值，Hermi?te插值，分段多项式?插值和线性?最小二乘拟?合。 4.1 实际案例 4.2 问题的描述?与基本概念? 先获得函数?（已知或未知?）在有限个点?上的值 … … 由表中数据?构造一个函?数P(x)作为f (x) 的近似函数?，去参与有关?f (x)的运算。 科学计算中?，解决不易求?出的未知函?数的问题主?要采用插值?和拟合两种?方法。 1）插值问题的?描述 已知函数在?[a,b]上的n+1个互异点?处的函数值?，求f (x) 的一个近似?函数P (x)， （4.1） P (x) 称为f (x)的一个插值?函数; f (x) 称为被插函?数;点为插值节?点; 称为插值条?件; 称为插值余?项。 当插值函数?P (x)是多项式时?称为代数插?值（或多项式插?值）。 一个代数插?值函数P (x)可写为 若它满足插?值条件（4.1），则有线性方?程组 （4.2） 当m=n，它的系数行?列式为范德?蒙行列式 因为插值节?点互异，，故线性方程?组（4.2）有唯一解，于是有定理 当插值节点?互异时，存在一个满?足插值条件?定理 满足插值条?件足插值条件? 令 显然有是次?数≤n的多项式?，且 说明有n+1个零点，由代数基本?定理有H (x) ? 0，由此得。 插值的一个?目的是对函?数作近似计?算。 假设[a, b] 是包含插值?点的最小闭?区间，当用插值函?数P(x)来近似计算?x在[a, b]的函数值时?，称为内插计?算，否则称为外?插或外推计?算。 2）拟合问题的?描述 已知在[a,b]上的n+1个(互异或不互?异)点处的函数?值，求f (x) 的一个近似?函数， 这里是n+1维向量，是某种范数?，，。 求出的称为?拟合函数。 3）插值函数和?拟合函数的?几何解释 1) 插值函数图?示 2）拟合函数图?示 4.3插值法 Lagra?nge插值? Lagra?nge插值?是 n次多项式?插值。 基本思想 将待求的n?次多项式插?值函数改写?成用已知函?数值为系数?的n+1个待定n?次多项式的?线性组合型?式，再利用插值?条件和函数?分解技术确?定n+1个待定n?次多项式形?式求出插值?多项式。 1) 构造原理 已知数表 … … 设n次插值?多项式 （4.3） 式中是与无?关的n次多?项式。 由插值条件?（4.1），有 由于与无关?，可得 （4.4） 为确定，注意到是n?次多项式，由式（4.4）可知 式中a为待?定常数，由确定，于是有 （4.5） 代入式（4.3），有 （4.6） 由n次插值?多项式的唯?一性，可知就是所?求的n次插?值多项式。 式（4.6）称为n次L?agran?ge插值多?项式，而称为 Lagra?nge插值?基函数。 2) 分析 定理4.3. 设函数在[a,b]上有n+1阶导数，是满足插值?条件的n次?插值多项式?，则有对任何?成立 式中。 证明 因为，故有 于是Rn（x）可分解为 （4.8） 为求出k(x)，做辅助函数? （4.9） 则有在时，g(t)=0，即g(t)在[a,b]上有n+2个零点。 显然g(t)在由组成的?n+1个小闭区?间上满足R?olle中?值定理，故g(t) 在[a,b]上有n+1个零点。 类似的有g??(t) 在[a,b]上有n个零?点，反复运用R?olle中?值定理，有在[a,b]上有1个零?点，设为?，则有。 在式（4.9）两边对t求?n+1阶导数，有 将t =? 代入上式，解得 代入式（4.8），即得定理结?果。 定理4.3中若能算?出在[a,b]上的最大值?，则有余项估?计式 （4.10） 此式可以估?计插值函数?在某一点的?误差，若想估计函?数在插值区?间[a,b]上的误差，则要计算出?。在整个插值?区间[a,b]上的误差估?计为 （4.11） 由n次插值?多项式的唯?一性及式（4.7），得到有如下?重要结果 定理4.4 若函数f ( x )在[a,b]上有n+1阶导数，则f ( x )可表示为 对n=1的插值多?项式，称为线性插?值； n=2的插值多?项式称为抛?物线插值或?辛普森插值?. 例4.1 已知的函数?表为 x 3.0 3.1 3.2 3.3 3.4 y=f(x) 1.09861?2 1.13140?2 1.1631