第四章插值与拟合精选文档.ppt

《第四章插值与拟合精选文档.ppt》由会员分享，可在线阅读，更多相关《第四章插值与拟合精选文档.ppt（86页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、第四章插值与拟合本讲稿第一页，共八十六页插值问题的提出插值问题的提出在许多工程实际问题中，常有如下情况：在许多工程实际问题中，常有如下情况：函数关系没有明显的解析表达式，只有实验观函数关系没有明显的解析表达式，只有实验观测来确定与自变量的测来确定与自变量的某些某些相对应的函数值；相对应的函数值；l函数虽然有解析表达式，但是使用很不方便。函数虽然有解析表达式，但是使用很不方便。对上述问题中的函数建立一个简单的便于计算和处对上述问题中的函数建立一个简单的便于计算和处理的近似表达式，即用理的近似表达式，即用一个简单的函数一个简单的函数来近似来近似代替代替这些这些不便处理的函数不便处理的函数插值函数插

2、值函数。4 插值与曲线拟合插值与曲线拟合本讲稿第二页，共八十六页拟合问题的提出拟合问题的提出实际问题中通过测量得到的数据比较多，而且这些数据实际问题中通过测量得到的数据比较多，而且这些数据本身含有一定误差，根据这些数据求取近似函数的方法是本身含有一定误差，根据这些数据求取近似函数的方法是曲线拟合曲线拟合。插值要求找到的近似函数的曲线通过所有的数据点。插值要求找到的近似函数的曲线通过所有的数据点。曲线拟合不要求近似函数的曲线通过所有数据，只要求曲线拟合不要求近似函数的曲线通过所有数据，只要求该曲线能反映数据变化的基本趋势。该曲线能反映数据变化的基本趋势。拟合的主要目的是：拟合的主要目的是：去掉测

3、量数据所含的测量误差。去掉测量数据所含的测量误差。插值与拟合的区别插值与拟合的区别本讲稿第三页，共八十六页插值与拟合插值与拟合插值插值拟合拟合插值插值:过点过点;(适合精确数据适合精确数据)拟合拟合:不过点不过点,整体近似整体近似;(适合经验公式或有误适合经验公式或有误 差的差的数据数据)本讲稿第四页，共八十六页4.1 插值问题的数学提法插值问题的数学提法已知已知a,ba,b上的函数上的函数y=y=f f(x x)在在n+1n+1个互异点处的函数值个互异点处的函数值:fnf2f1f0f(x)xnx2x1x0 x求简单函数求简单函数P Pn n(x x)，使得，使得计算计算f f(x x)可通过

4、计算可通过计算P P(x x)来近似代替。如下图所示。来近似代替。如下图所示。x0 x1x2x3x4xP(x)f(x)本讲稿第六页，共八十六页这就是插值问题这就是插值问题,(*)(*)式为插值条件式为插值条件,称为插值节点称为插值节点由于插值函数由于插值函数 的选择不同，就产生不同类型的插的选择不同，就产生不同类型的插值。若值。若 为代数多项式为代数多项式，就是代数插值，若就是代数插值，若 为三角多项式就称为三角多项式插值，若为三角多项式就称为三角多项式插值，若 为为有理函数就称为有理函数插值。由于代数多项式结有理函数就称为有理函数插值。由于代数多项式结构简单，本章主要介绍构简单，本章主要介绍

5、代数多项式代数多项式插值问题。插值问题。本讲稿第七页，共八十六页2.2.满足插值条件的多项式满足插值条件的多项式 P(x)是否存在且唯一？是否存在且唯一？3.3.用用P(x)P(x)代替代替f(x)f(x)的误差估计，即截断误差的估计；的误差估计，即截断误差的估计；对于多项式插值，我们主要讨论以下几个问题对于多项式插值，我们主要讨论以下几个问题：4.4.当插值节点无限加密时，插值函数是否收敛于当插值节点无限加密时，插值函数是否收敛于f(x)f(x)。1.如何构造出插值多项式如何构造出插值多项式P(x);即插值多项式的常用构造方即插值多项式的常用构造方法有哪些？法有哪些？本讲稿第八页，共八十六页

6、4.2 拉格朗日插值拉格朗日插值可见可见 P(x)是过是过(x0,y0)和和(x1,y1)两点的直线。这种两点的直线。这种插值称为插值称为线性插值线性插值，显然在节点上插值误差为，显然在节点上插值误差为0。)()(001010 xxxxyyyxP +=101xxxx 010 xxxx =y0+y1l0(x)l1(x)拉格朗日1.线性插值（两点插值）线性插值（两点插值）已知函数已知函数 在节点在节点 有函值有函值 ，求作一次多项式求作一次多项式 使得使得1100)(,)(yxPyxP=本讲稿第九页，共八十六页于是线性插值函可以表示为函数值于是线性插值函可以表示为函数值 与基函数的与基函数的线性组

7、合线性组合 与与 称为称为线性插值基函数线性插值基函数。它有如下性质：。它有如下性质：即：即：本讲稿第十页，共八十六页所以所以 例例1 1 已知已知 用线性插值求用线性插值求 近近 似值。似值。基函数分别为基函数分别为:解解插值多项式为插值多项式为本讲稿第十一页，共八十六页2.抛物线插值（三点二次插值）抛物线插值（三点二次插值）已知已知 在节点在节点 上的函数值上的函数值 ，求二次多项式求二次多项式 ，使之满足，使之满足 根据要满足的三个条件，确定三个未知数根据要满足的三个条件，确定三个未知数 ，因此可采用待定系数法。即因此可采用待定系数法。即本讲稿第十二页，共八十六页为避免解线性方程组，下面

8、仿线性插值，用基函为避免解线性方程组，下面仿线性插值，用基函数的方法求解方程组。数的方法求解方程组。设方程组满足条件设方程组满足条件 的方程为的方程为其中基函数应满足其中基函数应满足:xx0 x1x2l0(x)100l1(x)010l2(x)001本讲稿第十三页，共八十六页以以 为例说明基函数的求取方法，当为例说明基函数的求取方法，当 取取 ，时，时，为为0 0，取取 时，时，因此，因此其中其中 可用可用 求出，求出，同理同理所以所以本讲稿第十四页，共八十六页 抛物线插值是三个二次式的线性组合，是抛物线插值是三个二次式的线性组合，是抛物线插值是三个二次式的线性组合，是抛物线插值是三个二次式的线

9、性组合，是 x 的（不高于）的（不高于）的（不高于）的（不高于）二次式，在节点上插值多项式的值和已知函数值相等二次式，在节点上插值多项式的值和已知函数值相等二次式，在节点上插值多项式的值和已知函数值相等二次式，在节点上插值多项式的值和已知函数值相等。本讲稿第十五页，共八十六页设函数设函数 在区间在区间 上有节点上有节点 上的函上的函数值，构造一个次数不超过数值，构造一个次数不超过 次的代数多项式次的代数多项式 ，使，使 。即。即 次代数插值满足在次代数插值满足在 个节点上插值多项式个节点上插值多项式 和被插和被插值函数值函数 相等，而且插值多项式相等，而且插值多项式 的次数不超过的次数不超过

10、次。次。3.n n次代数插值次代数插值次代数插值次代数插值这样的插值多项式能构造出来吗？唯一吗？这样的插值多项式能构造出来吗？唯一吗？本讲稿第十六页，共八十六页插值多项式的存在性和唯一性定理插值多项式的存在性和唯一性定理证证 ：设所求的插值多项式为：设所求的插值多项式为 P(x)=a0+a1x+a2x2+.+anxn则由插值条件式则由插值条件式Pn(xi)=yi (i=0,1,.,n)可得关于系数可得关于系数a0,a1,an的线性代数方程组的线性代数方程组 设节点设节点xi(i=0,1,n)互异互异,则则满足插值满足插值条件条件P(xi)=yi 的的n次多项式次多项式P(x)=a0+a1x+a

11、2x2+.+anxn 存在且唯一存在且唯一.定理本讲稿第十七页，共八十六页 此方程组有此方程组有n n+1+1个方程个方程,n n+1+1个未知数个未知数,其系数行列式其系数行列式是范德蒙行列式，即：是范德蒙行列式，即：本讲稿第十八页，共八十六页由由于于插插值值节节点点 xi 互互不不相相同同,所所有有因因子子 xj-xi 0,所所以以上上述述行行列列式式不不等等于于零零,故故由由克克莱莱姆姆法法则则知知方方程程组组的的解解存存在且唯一在且唯一.即满足条件式即满足条件式 的的n次多项式次多项式存在且唯一。存在且唯一。证毕。证毕。本讲稿第十九页，共八十六页4.4.拉格朗日插值多项式拉格朗日插值多

12、项式拉格朗日插值多项式拉格朗日插值多项式 要求要求n阶插值多项式阶插值多项式，可以通过求方程组的解，可以通过求方程组的解 得到。得到。但由于解线性代数方程组的计算量比较大，构造插值多但由于解线性代数方程组的计算量比较大，构造插值多项式时，仍用项式时，仍用基函数基函数构造。构造。希望能找到满足以下条件的希望能找到满足以下条件的n 次多项式次多项式 l i(x)然后令然后令，则显然有则显然有P(xi)=yi。基函数的构造基函数的构造本讲稿第二十一页，共八十六页其中其中A A为常数为常数,由由li(xi)=1可得可得称之为称之为拉格朗日基函数拉格朗日基函数。li(x)除除xi 点外点外,其余其余xj

13、 都是都是li(x)的根，可设的根，可设与有关，而与无关节点节点f本讲稿第二十二页，共八十六页拉格朗日插值多项式拉格朗日插值多项式 特别地特别地,当当 n=1时又叫线性插值时又叫线性插值,其几何意义为过两点其几何意义为过两点的直线的直线.当当n=2时又叫抛物插值时又叫抛物插值,其几何意义为过三点的其几何意义为过三点的抛物线抛物线.利用拉格朗日基函数式利用拉格朗日基函数式l i(x),构造多项式构造多项式可知其满足可知其满足,称为称为拉格朗日拉格朗日插值多项式插值多项式。Pn(xi)=yi (i=0,1,.,n)应注意应注意,对于插值节点对于插值节点,只要求它们互异只要求它们互异,与与大小次序大

14、小次序无关。无关。本讲稿第二十三页，共八十六页5 插值余项插值余项 截断误差截断误差R(x)=f(x)-P(x)也称为插值多项式的余项。也称为插值多项式的余项。以下为拉格朗日余项定理。以下为拉格朗日余项定理。定理定理 设设f(x)在区间在区间a,b上存在上存在n+1 阶导数阶导数,xi a,b(i=0,1,n)为为n+1个互异节点个互异节点,则对任何则对任何x a,b,有有且与且与x 有关有关)其中其中本讲稿第二十四页，共八十六页注意这里是对t 求导证明：考察截断误差证明：考察截断误差RollesTheorem:若充分光滑，则存在使得。推广：推广：若若使得使得使得使得存在存在使得使得R(x)至

15、少有至少有 个根个根n+1=niixxxKxR0)()()(任意固定任意固定 x xi (i=0,n),构造构造=niixtxKtRt0)()()()(t)有有 n+2 个不同的根个不同的根 x0 xn x!)1()()()1(+-+nxKRxnx x=+!)1)()()()1()1(nxKPfxnxnx xx x!)1()()()1(+=+nfxKxnx x0=证毕证毕证毕证毕常用余项定理研究常用余项定理研究常用余项定理研究常用余项定理研究插值的误差估计插值的误差估计插值的误差估计插值的误差估计。本讲稿第二十五页，共八十六页其中其中LagrangeLagrange插值公式的标准型公式插值公式

16、的标准型公式:插值余项插值余项：本讲稿第二十六页，共八十六页 线性插值线性插值线性插值线性插值 ，余项，余项，余项，余项 ，其中：其中：其中：其中：对对对对 求极值：求极值：求极值：求极值：得得得得 为极小值。为极小值。为极小值。为极小值。即即即即取取取取 ，则，则，则，则 。取绝对值：取绝对值：取绝对值：取绝对值：，则：，则：，则：，则：所以所以所以所以 其中：其中：其中：其中：用余项定理求用余项定理求线性插值线性插值线性插值线性插值余项及其估计式。余项及其估计式。余项及其估计式。余项及其估计式。n次插值次插值余项及其估计式。余项及其估计式。其中其中本讲稿第二十七页，共八十六页的抛物插值多项

17、式的抛物插值多项式,且计算且计算f(3)的近似值并估计误差的近似值并估计误差。例例 设设解解 插值为多项式插值为多项式本讲稿第二十八页，共八十六页于是于是因为因为可得可得误差估计误差估计本讲稿第二十九页，共八十六页例：例：已知已知分别利用分别利用 sin x 的的1次、次、2次次 Lagrange 插值计算插值计算 sin 50 并估计误差。并估计误差。解：解：n=1分别利用分别利用x0,x1 以及以及 x1,x2 计算计算利用利用这里这里而而sin 50 =0.7660444)185(50sin10 p pL0.77614外推外推 的实际误差的实际误差 0.010010.01001利用利用s

18、in 50 0.76008,内插内插 的实际误差的实际误差 0.00596 0.00596内插通常优于外推。选择要计算的x所在的区间的端点，插值效果较好。本讲稿第三十页，共八十六页1 Lagrange Polynomialn=2)185(50sin20 p pL0.76543sin 50 =0.76604442次插值的实际误差次插值的实际误差 0.00061 0.00061高次插值通常优于低次插值本讲稿第三十一页，共八十六页如果如果 是次数不超过是次数不超过 次的多项式，取次的多项式，取 个个 节点插值时，插值多项式就是其自身。节点插值时，插值多项式就是其自身。插值多项式插值多项式 只与数据只

19、与数据 ，有关，与节点排有关，与节点排列顺序无关。列顺序无关。个节点的插值多项式不超过个节点的插值多项式不超过 次次。拉格朗日插值小结拉格朗日插值小结：内插比外插精度高。当求某点的函数值时，插值内插比外插精度高。当求某点的函数值时，插值节点应尽可能靠近该节点应尽可能靠近该 点，此时余项小。点，此时余项小。当节点数变化时，需重新计算全部基函数当节点数变化时，需重新计算全部基函数。本讲稿第三十二页，共八十六页插值精度提高的条件插值精度提高的条件插值点与节点靠近插值点与节点靠近内插精度一般比外推高内插精度一般比外推高插值点适当多插值点适当多本讲稿第三十三页，共八十六页4.4 牛顿插值牛顿插值拉格朗日

20、插值的优点是插值多项式特别容易建立，拉格朗日插值的优点是插值多项式特别容易建立，缺点是增加节点时原有多项式不能利缺点是增加节点时原有多项式不能利 用，必须重新建用，必须重新建立，即所有基函数都要重新计算，这就造成计算量的立，即所有基函数都要重新计算，这就造成计算量的浪费；浪费；能否将能否将 P(x)改写成改写成的形式，希望每加一个节点时，的形式，希望每加一个节点时，只附加一项只附加一项上去即可。上去即可。?牛顿插值牛顿插值本节介绍这种插值公式的建立、特点和应用。这本节介绍这种插值公式的建立、特点和应用。这 要用要用到到差商差商的概念。的概念。本讲稿第三十四页，共八十六页一、差商及其基本性质一、

21、差商及其基本性质定义定义1称为称为 f(x)在在xi、xj点的点的一阶差商一阶差商.记为记为函数函数 在区间在区间 上的平均变化率上的平均变化率为函数为函数f(x)在点在点xi、xj、xk的的二阶差商二阶差商.记为记为同样地，称一阶差商的平均变化率同样地，称一阶差商的平均变化率本讲稿第三十五页，共八十六页差商的计算步骤与结果可列成差商表，如下差商的计算步骤与结果可列成差商表，如下一般地，一般地，n-1阶差商的差商阶差商的差商 称为称为f(x)在在x0,x1,xn点的点的 n 阶差商。阶差商。由此可知高阶差商是由比它低一阶的两个差商组成。由此可知高阶差商是由比它低一阶的两个差商组成。本讲稿第三十

22、六页，共八十六页xk函数值一阶差商二阶差商三阶差商.x0f(x0)x1 f(x1)f x0,x1x2 f(x2)f x1,x2f x0,x1,x2x3 f(x3)f x2,x3f x1,x2,x3f x0,x1,x2,x3.利用差商的递推定义，可以用递推来计算差商利用差商的递推定义，可以用递推来计算差商本讲稿第三十七页，共八十六页121520f(x)7431x例例1:已知：已知：计算三阶差商计算三阶差商f1,3,4,7-1.25-3.5-112741315412301三阶差商三阶差商二阶差商二阶差商一阶差商一阶差商f(x)x解：解：作差商表作差商表所以所以f1,3,4,7=-1.25本讲稿第三

23、十八页，共八十六页 fx0,x1,x2,.,xn =fx1,x2,.,xn,x0 性质性质1 差商可以表示为函数值的线性组合，即差商可以表示为函数值的线性组合，即=fx1,x0,x2,.,xn=这一性质称之为差商这一性质称之为差商的对称性的对称性。性质性质2 (对称性对称性)差商与节点的顺序无关，如差商与节点的顺序无关，如这一性质可以用数学归纳法证明。（这一性质可以用数学归纳法证明。（P124）本讲稿第三十九页，共八十六页二、牛顿插值多项式二、牛顿插值多项式如此继续下去，可得一系列等式如此继续下去，可得一系列等式设设x是是a，b上一点，由一阶差商定义上一点，由一阶差商定义同理，由二阶差商定义同

24、理，由二阶差商定义得得得得本讲稿第四十一页，共八十六页依次把后式代入前式，依次把后式代入前式，本讲稿第四十二页，共八十六页P(x)R(x)令：令：最后得最后得本讲稿第四十三页，共八十六页牛顿插值公式特点牛顿插值公式特点则：则：是：是x的的n次多项式，称为次多项式，称为牛顿插值多项式牛顿插值多项式：是：是 最后一项，称为最后一项，称为牛顿插值余项。牛顿插值余项。增加一个节点时，只需增加一个节点时，只需 再增加一项，再增加一项，其余各项都不变。其余各项都不变。在节点上，多项式在节点上，多项式 等于被插函数等于被插函数 的的值，即值，即 所以此时的余项所以此时的余项 ，满足插值条件。，满足插值条件。

25、本讲稿第四十四页，共八十六页牛顿插值公式特点牛顿插值公式特点(续续)取取n1个节点时，牛顿插值多项式次数不超过个节点时，牛顿插值多项式次数不超过n次，各次，各项系数是各阶差商。项系数是各阶差商。由插值多项式的唯一性知由插值多项式的唯一性知,拉格朗日插值多项式和拉格朗日插值多项式和牛顿插值多项式是相等牛顿插值多项式是相等 的的。只是算法不同。只是算法不同。拉格朗日插值余项和牛顿插值余项是相等拉格朗日插值余项和牛顿插值余项是相等 的，即：的，即：本讲稿第四十五页，共八十六页121520f(x)7431x例例2:已知：已知：求满足以上插值条件的牛顿型插值多项式。求满足以上插值条件的牛顿型插值多项式。

26、解：在解：在例例1 1中，我们已计算出中，我们已计算出则牛顿三次插值多项式为则牛顿三次插值多项式为 本讲稿第四十六页，共八十六页4.7 分段低次插值分段低次插值分段低次插值分段低次插值问题的提出：问题的提出：在拉格朗日插值法中，为了使插值多项式更好的在拉格朗日插值法中，为了使插值多项式更好的逼近被插函数，往往要增加插值节点，提高插值多项式逼近被插函数，往往要增加插值节点，提高插值多项式的次数。当插值区间较大时，对于高次插值的次数。当插值区间较大时，对于高次插值在插值区间在插值区间两端会发生剧烈振荡两端会发生剧烈振荡。高次代数插值所发生的这种现象。高次代数插值所发生的这种现象称为称为Runge现

27、象现象.在上个世纪初由在上个世纪初由Runge发现发现.本讲稿第四十七页，共八十六页一、一、高阶插值的龙格现象高阶插值的龙格现象例：例：在在 5,5上考察上考察 的的Pn(x)。取。取-5-4-3-2-1 0 1 2 3 4 5-0.5 0 0.5 1 1.5 2 2.5 n越大，端点附近抖动越大，称为Runge现象Pn(x)f(x)分段低次插值分段低次插值P2（x）P5（x）P10（x）既要增加插值节点，减小插值区间，又不增加插值多项式的次既要增加插值节点，减小插值区间，又不增加插值多项式的次数以减少误差，我们可数以减少误差，我们可 以以采用分段插值采用分段插值的办法。的办法。本讲稿第四十八

28、页，共八十六页二、二、分段线性插值分段线性插值原理原理：将插值区间将插值区间分成若干个小的子段分成若干个小的子段，在每个子段上进，在每个子段上进行低次插值，然后相互连接，用连接相邻节点的行低次插值，然后相互连接，用连接相邻节点的折线逼近折线逼近被插函数。被插函数。定义定义已知函数已知函数 ，给定节点，给定节点的函数值的函数值求一个分段函数求一个分段函数S(x)S(x)，使其满足，使其满足：S(x)在在a,b上连续；上连续；S(xi)=yi，i0,1,nS(x)在每个子段在每个子段 是是线性函数线性函数。称函数称函数S(x)S(x)为为分段线性插值函数分段线性插值函数。本讲稿第四十九页，共八十六

29、页S(x)在子段在子段 上有（线性插值）：上有（线性插值）：在区间在区间 上有：上有：其中其中显然显然本讲稿第五十页，共八十六页！失去了原函数的光滑性。失去了原函数的光滑性。分段线性插值特点分段线性插值特点算法简单，能保证收敛。算法简单，能保证收敛。能获得任意精度的插值。能获得任意精度的插值。局部性质。局部性质。本讲稿第五十一页，共八十六页1例例3 已知函数已知函数 在区间在区间 上取等上取等 距插值节点（如下表）距插值节点（如下表）,求区间求区间 上分段线性上分段线性插值函数，并利用它求出插值函数，并利用它求出 的近似值的近似值 0.038460.058820.10.20.51 54320

30、解解 在每个小区间在每个小区间 上上，本讲稿第五十二页，共八十六页于是于是本讲稿第五十三页，共八十六页4.3 4.3 逐次线性插值逐次线性插值 （本节为了解内容本节为了解内容）逐次线性插值解决拉格朗日插值为提高精度增加插值逐次线性插值解决拉格朗日插值为提高精度增加插值节点时，要重新计算全部基函数，整个插值多项式的结构节点时，要重新计算全部基函数，整个插值多项式的结构都会改变的问题。为使计算有都会改变的问题。为使计算有“承袭性承袭性”，可用逐次线性，可用逐次线性插值或称迭代插值的办法解决。逐次线性插值是重复地进行插值或称迭代插值的办法解决。逐次线性插值是重复地进行线性插值产生从低次到高次的拉格朗

31、日插值多项式序列，直线性插值产生从低次到高次的拉格朗日插值多项式序列，直到获得合适的计算结果，避免增加节点从头开始计算，常用到获得合适的计算结果，避免增加节点从头开始计算，常用的的埃特金埃特金(Aitken)算法算法和和内维尔内维尔(Neville)算法算法。本讲稿第五十四页，共八十六页4.3.1 三个节点的情形三个节点的情形已知已知 f(x)在三个互异节点在三个互异节点x0,x1,x2的函数值的函数值y0,y1,y2用（x0,y0),(x1,y1)做插值做插值xx0 x1 x2yy0 y1 y2用（x0,y0),(x2,y2)做插值做插值用（x1,P01),(x2,P02)做插值做插值本讲稿

32、第五十五页，共八十六页 上式即是拉格朗日二次插值多项式。上式即是拉格朗日二次插值多项式。两个线性两个线性插值的结果再进行插值的结果再进行线性线性插值，得到抛物插值，得到抛物线性线性插值。插值。本讲稿第五十六页，共八十六页三个节点的情形写成表格的形式三个节点的情形写成表格的形式函数值一阶插值二阶插值本讲稿第五十七页，共八十六页4.3.2 埃特金算法埃特金算法 两个线性插值的结果，再进行线性插值可得抛物两个线性插值的结果，再进行线性插值可得抛物线插值，这个过程可以继续下去，一般地，利用两线插值，这个过程可以继续下去，一般地，利用两个个 次插值次插值 和和 进行线性插值，可得进行线性插值，可得 次次

33、插值插值 ，用基函数的形式表示，用基函数的形式表示 当当 ，时，得到时，得到 ，当，当 ，时，得到时，得到 ，当，当 ，时，得到时，得到 。反复执行这一算式，可以逐步。反复执行这一算式，可以逐步构造出如下的插值顺序。按这个顺序可求出某点的构造出如下的插值顺序。按这个顺序可求出某点的插值。插值。本讲稿第五十八页，共八十六页函数值一阶插值二阶插值三阶插值埃特金插值表埃特金插值表本讲稿第五十九页，共八十六页埃特金插值特点埃特金插值特点 埃特金插值埃特金插值是将一个高次插值过程归结为线性插是将一个高次插值过程归结为线性插是将一个高次插值过程归结为线性插是将一个高次插值过程归结为线性插值的多次重复。埃特

34、金插值的每个数据均为插值结果，值的多次重复。埃特金插值的每个数据均为插值结果，值的多次重复。埃特金插值的每个数据均为插值结果，值的多次重复。埃特金插值的每个数据均为插值结果，从这些数据的一致程度即可判断插值结果的精度。这从这些数据的一致程度即可判断插值结果的精度。这从这些数据的一致程度即可判断插值结果的精度。这从这些数据的一致程度即可判断插值结果的精度。这样可以逐步生成插值结果，每做一步检查一下计算结样可以逐步生成插值结果，每做一步检查一下计算结样可以逐步生成插值结果，每做一步检查一下计算结样可以逐步生成插值结果，每做一步检查一下计算结果的精度，如不满足要求，则增加一个节点再算，直果的精度，如

35、不满足要求，则增加一个节点再算，直果的精度，如不满足要求，则增加一个节点再算，直果的精度，如不满足要求，则增加一个节点再算，直到满足要求为止。在节点较多时，用这种方法可降低到满足要求为止。在节点较多时，用这种方法可降低到满足要求为止。在节点较多时，用这种方法可降低到满足要求为止。在节点较多时，用这种方法可降低插值次数插值次数插值次数插值次数本讲稿第六十页，共八十六页 的近似值为的近似值为1.5。本讲稿第六十一页，共八十六页已知已知f(x)在三个互异点在三个互异点 0,1,2的函数值的函数值1,3,9用（0,1),(1,3)作插值作插值用（0,1),(2,9)作插值作插值用（1,P01),(2,

36、P02)作插值作插值本讲稿第六十二页，共八十六页一般得到的观测数据本身不一定完全可靠，个别数据的误差甚一般得到的观测数据本身不一定完全可靠，个别数据的误差甚至可能很大，且数据很多。曲线拟合是从已知的一大堆数据中找出至可能很大，且数据很多。曲线拟合是从已知的一大堆数据中找出规律，即设法构造一条曲线（拟合曲线）反映数据点总的趋势。规律，即设法构造一条曲线（拟合曲线）反映数据点总的趋势。4.8 4.8 曲线拟合的最小二乘法曲线拟合的最小二乘法 曲线拟合问题曲线拟合问题曲线拟合同插值法的区别曲线拟合同插值法的区别:曲线拟合考虑实验曲线拟合考虑实验数据带有测量误数据带有测量误差差，同时，同时测量数据又多

37、测量数据又多，若用插值法时得到的插值多项式次数将很，若用插值法时得到的插值多项式次数将很高，不实用高，不实用。曲线拟合得到的多项式不必通过每一点。曲线拟合得到的多项式不必通过每一点。本讲稿第六十三页，共八十六页设已知某物理过程（函数关系）设已知某物理过程（函数关系）的一组的一组观测（实验）数据观测（实验）数据 ，要求在某，要求在某特特定函数类定函数类（例如（例如多项式多项式）中找出一个函数）中找出一个函数 ，作，作为为 的近似函数，使得在的近似函数，使得在 上的误差（或称上的误差（或称残差残差）按某种按某种度量标准度量标准为最小，这就是拟合问题，也称为最小，这就是拟合问题，也称曲线曲线拟合。拟

38、合。曲线拟合的定义曲线拟合的定义本讲稿第六十四页，共八十六页记向量记向量 ，要求残差，要求残差 按某种度量按某种度量标准为最小，即要求向量标准为最小，即要求向量 的某种范数的某种范数 最小。例最小。例如，可要求如，可要求 的的1-1-范数范数 或或 -范数范数 即：即：或或 最小。但由于最小。但由于不便于微分运算不便于微分运算，因此，通常用，因此，通常用2-2-范数。范数。1 1 为最小，这种要求误差平方和最小的拟合称为为最小，这种要求误差平方和最小的拟合称为曲线拟曲线拟合的最小二乘法合的最小二乘法。一、最小二乘拟合原理一、最小二乘拟合原理本讲稿第六十五页，共八十六页最小二乘法的最小二乘法的几

39、何意义几何意义 求在给定点求在给定点 x0,x1,xm 处与点处与点(x0,y0),(x1,y1),(xm,ym)的距离平方和最小的曲线：的距离平方和最小的曲线：y=F(x),本讲稿第六十六页，共八十六页二、多项式拟合二、多项式拟合这样的曲线拟合叫这样的曲线拟合叫多项式拟合多项式拟合。满足上式的满足上式的y y 叫叫最小二最小二乘拟合多项式乘拟合多项式.特别地特别地,当当m m=1=1时时,一次多项式拟合又叫一次多项式拟合又叫直直线拟合线拟合。当当N=m+1时，所得的拟合多项式就是拉格朗日或牛顿时，所得的拟合多项式就是拉格朗日或牛顿插值多项式。插值多项式。对给定的一组数据对给定的一组数据(x

40、xi i,y yi i)()(i i=1,=1,N N),),寻求寻求m m次次多项式（多项式（N N m）使总误差满足使总误差满足目标函数本讲稿第六十七页，共八十六页由于目标函数由于目标函数J 可看作是关于可看作是关于的多元函数，故拟合多项式的问题变为求极值问题。的多元函数，故拟合多项式的问题变为求极值问题。由求极值的由求极值的必要条件必要条件得得即即对a0偏导对a1偏导例本讲稿第六十八页，共八十六页上述方程组是关于上述方程组是关于a aj j的线性方程组，通常称为的线性方程组，通常称为正则方程组。正则方程组。它有唯一解。由方程组可求出系数它有唯一解。由方程组可求出系数a aj j。多项式拟

41、合的一般方法可归纳为以下几步：多项式拟合的一般方法可归纳为以下几步：写出正规方程组，求出系数，得到写出正规方程组，求出系数，得到拟合多项式拟合多项式注注：当：当m m较大时，方程组一般病态。较大时，方程组一般病态。由已知数据画出函数粗略的图形由已知数据画出函数粗略的图形散点图，确散点图，确定拟合多项式的次数定拟合多项式的次数m；列表计算列表计算本讲稿第六十九页，共八十六页例例1 测得铜导线在温度测得铜导线在温度 时的电阻时的电阻 如下表求电阻如下表求电阻R与温度与温度T的近似函数关系。的近似函数关系。85.1083.9082.3580.8079.2577.8076.30 50.045.140.

42、036.030.125.019.1 6543210i解：把表中数据画在坐标纸上，会看到数据点的分布解：把表中数据画在坐标纸上，会看到数据点的分布可用一条直线来近似描述。因此用线性函数拟合。可用一条直线来近似描述。因此用线性函数拟合。本讲稿第七十页，共八十六页列表如下列表如下0 20029.4459325.83565.5245.34255.0002500.0085.1050.063783.8902034.0183.9045.153294.0001600.0082.3540.042908.8001296.0080.8036.032385.425906.0179.2530.121945.000625

43、.0077.8025.011457.330364.8176.3019.1 Ri i本讲稿第七十一页，共八十六页故得故得R与与T的拟合直线为的拟合直线为 R=70.572+0.291T正规方程组为正规方程组为解方程组得解方程组得本讲稿第七十二页，共八十六页列表如下列表如下解解 设拟合曲线方程为设拟合曲线方程为 例例2 已知实验数据如下表已知实验数据如下表试用最小二乘法求它的二次拟合多项式试用最小二乘法求它的二次拟合多项式 i4321124510 yi1098765431 xi876543210本讲稿第七十三页，共八十六页 10251472531730173813253 4004010000100

44、0100410824327656172981397128164096512642864972401343491753661296216361645010625125252536416256641644245158127953110101111010 i本讲稿第七十四页，共八十六页得正规方程组得正规方程组解得解得故拟合多项式为故拟合多项式为本讲稿第七十五页，共八十六页三、超定方程组的最小二乘解三、超定方程组的最小二乘解当方程组中方程的个数多于未知量的个数时，称此方当方程组中方程的个数多于未知量的个数时，称此方程组为超定方程组。一般来说程组为超定方程组。一般来说 ，超定方程组无解（此时为，超定方程

45、组无解（此时为矛盾方程组），这时需要寻找方程组的一个矛盾方程组），这时需要寻找方程组的一个“最近似最近似”的的解。解。设线性方程组设线性方程组(mn)本讲稿第七十六页，共八十六页把方程组写成矩阵形式把方程组写成矩阵形式如果有向量如果有向量x x使得使得达到最小，则称达到最小，则称x x为为超定方程组的最小二乘解超定方程组的最小二乘解。定义残差向量定义残差向量取目标函数取目标函数本讲稿第七十七页，共八十六页曲线拟合的最小二乘法可以看成求超定方程组的最小二乘曲线拟合的最小二乘法可以看成求超定方程组的最小二乘解的问题解的问题。利用矩阵运算可得利用矩阵运算可得从而有正则方程组从而有正则方程组解以上方程

46、组解以上方程组可得超定方程组的解。可得超定方程组的解。将给定的数据将给定的数据(x xi i,y yi i)()(i i=1,=1,n n),),代入代入m m次多项式次多项式本讲稿第七十八页，共八十六页得到矛盾方程组得到矛盾方程组写成矩阵形式写成矩阵形式本讲稿第七十九页，共八十六页对应正则方程组对应正则方程组求正则方程组求正则方程组可得到所给数据的拟合多项式。可得到所给数据的拟合多项式。本讲稿第八十页，共八十六页要求熟练掌握的内容：要求熟练掌握的内容：线性插值、抛物线插值的算法。线性插值、抛物线插值的算法。拉格朗日插值多项式的构成，基函数及其特点；拉格朗日插值多项式的构成，基函数及其特点；差商的定义；牛顿插值公式。差商的定义；牛顿插值公式。最小二乘法拟合原理及多项式拟合的算法。最小二乘法拟合原理及多项式拟合的算法。要求掌握的内容要求掌握的内容：拉格朗日插值多项式的余项及结果的误差估计；拉格朗日插值多项式的余项及结果的误差估计；了解高次插值的龙贝格现象和分段插值。了解高次插值的龙贝格现象和分段插值。本章要求本章要求本讲稿第八十一页，共八十六页作作 业业P186 第第3题题 第第4题题 本讲稿第八十二页，共八十六页函数类（即函数的形式）函数类（即函数的形式）l代数多项式代数多项式 l 三角函数三角函数l其它正交多项式其它正交多项式 本讲稿第八十三页，共八十六页