第四章 函数的插值与拟合法优秀课件.ppt

《第四章 函数的插值与拟合法优秀课件.ppt》由会员分享，可在线阅读，更多相关《第四章 函数的插值与拟合法优秀课件.ppt（46页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第四章 函数的插值与拟合法第1页，本讲稿共46页定义定义 4.1 设设 y=f(x)在区间在区间a,b上连续，在上连续，在a,b内内n+1个互不个互不 相同的点相同的点 上取值上取值 .求一个性态较好的简单函数求一个性态较好的简单函数P(x)，使，使 得得 则称则称P(x)为为f(x)的的插值函数插值函数a,b-插值区间插值区间-插值节（结）点插值节（结）点f(x)-被插函数被插函数（4-1）-插值条件插值条件求插值函数求插值函数P(x)的方法的方法-插值法插值法一、插值函数一、插值函数4.1 引言引言 第2页，本讲稿共46页yxo插值插值第3页，本讲稿共46页(1)当当P(x)为次数不超过为

2、次数不超过n次的代数多项式时，相应的插值法称为次的代数多项式时，相应的插值法称为 多项式插值多项式插值；(2)当当P(x)为三角多项式时，相应的插值法称为为三角多项式时，相应的插值法称为三角插值三角插值；(3)当当P(x)为分段解析函数时，相应的插值法称为为分段解析函数时，相应的插值法称为分段插值分段插值。其中三角插值主要用于处理周期函数。其中三角插值主要用于处理周期函数。本章仅介绍本章仅介绍最基本的多项式插值最基本的多项式插值。第4页，本讲稿共46页定理定理 4.1 在在 n+1 个互异点个互异点 上满足插值条上满足插值条 件件(4-1)的次数不超过的次数不超过n次的插值多项式次的插值多项式

3、 存在且惟一。存在且惟一。所以，解存在且惟一，这说明由式所以，解存在且惟一，这说明由式(4-2)表示的表示的 存在存在且惟一，证毕。且惟一，证毕。证证二、多项式插值的唯一性二、多项式插值的唯一性设设 有有 第5页，本讲稿共46页 4.2 插值多项式的构造插值多项式的构造 一、一、基本插值多项式基本插值多项式 定义定义 下列表函数下列表函数 的插值多项式的插值多项式 叫做以叫做以 为节点的基本插为节点的基本插值多项式。值多项式。由定义可知由定义可知xx0 x1-xi-1xixi+1-xny00-010-0 4.2.1 拉格朗日插值多项式拉格朗日插值多项式(4-4)第6页，本讲稿共46页解解-n次

4、次Lagrange插值基函数插值基函数第7页，本讲稿共46页注：注：1.n+1个节点个节点n+1个基本插值多项式。个基本插值多项式。2.仅与节点有关，与仅与节点有关，与f(x)无关。无关。二二、Lagrange插值多项式插值多项式求下列列表函数的多项式求下列列表函数的多项式Ln(x)xx0 x1-xi-1xixi+1-xnyy0y1-yi-1yiyi+1-yn-n次拉格朗日插值多项式次拉格朗日插值多项式 解解 第8页，本讲稿共46页线性插值线性插值(n=1)，抛物插值抛物插值 (n=2)注：注：1.是是 的线性组合。的线性组合。2.与节点的排列顺序无关。与节点的排列顺序无关。第9页，本讲稿共4

5、6页例：已知列表函数例：已知列表函数,并计算并计算f(0.5)的计算值。的计算值。解：解：x-1012y111-5第10页，本讲稿共46页第11页，本讲稿共46页三、三、Lagrange插值多项式的余项插值多项式的余项 定理定理 4.2（误差估计定理）（误差估计定理）第12页，本讲稿共46页注注（1）余项公式主要用于理论分析。实际使用时，代）余项公式主要用于理论分析。实际使用时，代 之以误差估计式之以误差估计式 （2）插值节点的选取应尽量靠近插值点，以使）插值节点的选取应尽量靠近插值点，以使 尽可能小，以减小误差。尽可能小，以减小误差。推论推论 第13页，本讲稿共46页例例 4.1 给定函数表

6、给定函数表试分别用线性插值和抛物插值求试分别用线性插值和抛物插值求ln 1.46的近似值并估计误差。的近似值并估计误差。x1.21.31.41.51.61.7lnx0.1823220.2623640.3364720.4054650.4700040.530628解解作线性插值作线性插值 得得第14页，本讲稿共46页 作抛物插值作抛物插值 第15页，本讲稿共46页 4.2.2 牛顿均差插值多项式牛顿均差插值多项式 一、均差，均差表一、均差，均差表 定义定义 设设 f(x)在在a,b上连续，及自变量上连续，及自变量 第16页，本讲稿共46页第17页，本讲稿共46页节点节点 一阶均差一阶均差 二阶均差

7、二阶均差 三阶均差三阶均差 四阶均差四阶均差 五阶均差五阶均差第18页，本讲稿共46页例例 4.3 试用列表法对下例表格函数求试用列表法对下例表格函数求 f1,3,5,7x012345678f(x)107-3-19-39-59-71-599列表计算得列表计算得 xif(xi)一阶均差一阶均差二阶均差二阶均差三阶均差三阶均差13577-19-59-59-13-200-1.7551.125所以所以 f 1,3,5,7=1.125解解第19页，本讲稿共46页二、均差的性质二、均差的性质这性质又称为均差这性质又称为均差关于自变量对称关于自变量对称第20页，本讲稿共46页第21页，本讲稿共46页三、牛顿

8、均差插值多项式三、牛顿均差插值多项式第22页，本讲稿共46页第23页，本讲稿共46页Nn(x)称为称为牛顿均差插值多项式牛顿均差插值多项式。证证第24页，本讲稿共46页第25页，本讲稿共46页例例1：已知：已知x0.400.550.650.800.901.05y0.410750.578150.696150.888111.026521.25386第26页，本讲稿共46页0 0.40 0.410751 0.55 0.57815 1.11602 0.65 0.69615 1.1860 0.28003 0.80 0.88811 1.2757 0.3583 0.1974 0.90 1.02652 1.3

9、848 0.4336 0.214 0.0345 1.05 1.25386 1.5156 0.5260 0.231 0.034 0 k xk yk 一阶一阶 二阶二阶 三阶三阶 四阶四阶 五阶五阶 解解:(1)(2)与与0.596最接近的三个节点最接近的三个节点第27页，本讲稿共46页例例 4.4 给定表格函数给定表格函数 x12345f(x)0.50.1751.31-1.49510.36（1）试用二次牛顿均差插值法求）试用二次牛顿均差插值法求 f(2.8)的近似值；的近似值；（2）设）设 f(x)=-1.166 已知，试用（已知，试用（1）中构造的插值多项）中构造的插值多项 式求式求 x 的近

10、似值。的近似值。解解(1)选取节点选取节点x=2,3,4kxkf(xk)一阶均差一阶均差二阶均差二阶均差三阶均差三阶均差020.1751.135-1.97-0.9131.31-2.0805-1.0724-1.495-0.665310.5第28页，本讲稿共46页 4.3 分段低次插值（略）分段低次插值（略）第29页，本讲稿共46页 4.4 最小二乘法最小二乘法 4.4.1 最小二乘法的提出最小二乘法的提出 xx0 x1-xi-1xixi+1-xnyy0y1-yi-1yiyi+1-yn已知立表函数：已知立表函数：函函数数插插值值问问题题现在的问题：现在的问题：对上面立表函数求一近似函数（不要求插值

11、）对上面立表函数求一近似函数（不要求插值）y=P(x)满足：满足：（1）反映立表函数的一般趋势）反映立表函数的一般趋势（2）要求偏离立表函数很小）要求偏离立表函数很小-函数的拟合函数的拟合第30页，本讲稿共46页yxo插值插值拟合拟合方差方差第31页，本讲稿共46页 4.4.2 数据的多项式最小二乘拟合数据的多项式最小二乘拟合 xx0 x1-xi-1xixi+1-xnyy0y1-yi-1yiyi+1-yn已知一组数据：已知一组数据：解解-这个多项式称为这组数据的这个多项式称为这组数据的 最小二乘拟合多项式最小二乘拟合多项式第32页，本讲稿共46页第33页，本讲稿共46页它称为它称为法方程组法方

12、程组（或（或正规方程组正规方程组）。第34页，本讲稿共46页解得解得第37页，本讲稿共46页例例4.5 试对以下数据进行多项式拟合试对以下数据进行多项式拟合 xi12345678yi1.13.88.715.624.6 37.4 49.664.2解解第38页，本讲稿共46页第39页，本讲稿共46页注注第40页，本讲稿共46页第41页，本讲稿共46页x-2012y1.905.113.8例例 4.6 用最小二乘法求形如用最小二乘法求形如 y=ax+bx2 的多项式，使与下的多项式，使与下列数据拟合（得数保留三位小数）列数据拟合（得数保留三位小数）解解第42页，本讲稿共46页例例 4.7 给定数据给定

13、数据试求形如试求形如 y=a+bx2 的拟合多项式（得数保留三位小数）。的拟合多项式（得数保留三位小数）。解解x-5012y52.63.45.510.5第43页，本讲稿共46页 4.4.3 最小二乘法的应用例最小二乘法的应用例 例例 4.8 试对下表数据进行试对下表数据进行指数指数拟合拟合 xi123456yi15.420.426.837.248.863.3 4.4.3.1 数据的指数拟合数据的指数拟合 解解xi123456Zi=ln yi2.73443.01553.28843.61633.88774.1479第44页，本讲稿共46页xi123456yi15.420.426.837.248.863.3xi123456Zi=ln yi2.73443.01553.28843.61633.88774.1479第45页，本讲稿共46页 4.4.3.2 超定方程组的最小二乘解超定方程组的最小二乘解 解解例例 4.9 试求以下超定方程组的最小二乘解。试求以下超定方程组的最小二乘解。第46页，本讲稿共46页