《数学建模公选课课件》第七讲 多项式插值与曲线拟合.docVIP

插值与拟合 §1多项式插值问题 已知函数在n+1个互异结点处的函数值，如下表所示： xi x0 x1 …… xn yi=f(xi) y0 y1 …… yn 求一个n次多项式Pn(x)，使得Pn(xi)=yi，i=1,2,……,n。并利用Pn(x)近似未知函数f(x)。 从几何上看就是寻找一条n次多项式曲线y= Pn(x)，使其通过平面上已知的n+1个点： 一、Lagrange插值 其中， 二、Newton插值 其中， 随着插值结点的增多，插值多项式的次数也增加。然而多项式次数越高，近似效果未必越好，反而容易出现高次插值的Runge现象，为此需要考虑下面的分段插值问题。 三、分段插值 1、分段线性插值 在相邻两个结点[xk,xk+1]内，求一条线段近似函数f(x)，x∈[xk,xk+1]。 2、分段抛物插值 在相邻三个结点之间用抛物线近似未知函数。 四、样条插值 分段线性插值虽然避免了高次多项式插值的Runge现象，然而在插值结点处又产生了新问题：不光滑。为了克服这一现象，引入三次样条插值：在相邻两个结点之间用三次多项式函数si(x)近似未知未知函数，并保证在插值结点处满足衔接条件： 五、Matlab插值命令 yi=interp1(x，y，xi，method) (x，y)：插值节点；xi：被插值点；yi：xi处的插值结果； method：插值方法；‘nearest’ 最邻近插值；‘linear’ 线性插值； ‘spline’ 三次样条插值；‘cubic’ 立方插值。缺省时：分段线性插值。 注意：所有的插值方法都要求x是单调的，并且xi不能够超过x的范围。 例1：在1-12的11小时内，每隔1小时测量一次温度，测得的温度依次为：5，8，9，15，25，29，31，30，22，25，27，24。试估计每隔1/10小时的温度值。 Matlab命令 hours=1:12; temps=[5 8 9 15 25 29 31 30 22 25 27 24]; h=1:0.1:12; t=interp1(hours,temps,h,spline); plot(hours,temps,+,h,t,hours,temps,r:) %作图 xlabel(Hour),ylabel(Degrees Celsius’) 练习：下列数据表示从1790年到2000年的美国人口数据，利用这些数据给出1790年—2000年每隔5年的美国人口数据，并预测2005年、2010年人口数，与实际值比较。（人口单位：千人） 年份 1790 1800 1810 1820 1830 1840 人口 3929 5308 7240 9638 12866 17069 年份 1850 1860 1870 1880 1890 1900 人口 23192 31443 38558 50156 62948 75995 年份 1910 1920 1930 1940 1950 1960 人口 91972 105711 122755 131669 150697 179323 年份 1970 1980 1990 2000 2005 2010 人口 203212 226505 248710 281416 ？ ？ §2二维插值问题 已知函数在若干结点处的函数值， ， 找一个简单二元函数使其通过已知的m×n个结点，并用此函数近似f(x, y)。 一、结点类型 1、网格型结点 2、散乱结点 二、插值方法 1、最邻近插值 二维或高维情形的最邻近插值，与被插值点最邻近的节点的函数值即为所求。 注意：最邻近插值一般不连续。具有连续性的最简单的插值是分片线性插值。 2、分片线性插值 将四个插值点（矩形的四个顶点）处的函数值依次简记为： f (xi, yj)=f1，f (xi+1, yj)=f2，f (xi+1, yj+1)=f3，f (xi, yj+1)=f4 分两片的函数表达式如下： 第一片(下三角形区域)： (x, y)满足 插值函数为： 第二片(上三角形区域)：(x, y)满足 插值函数为： 注意：(x, y)当然应该是在插值节点所形成的矩形区域内。显然，分片线性插值函数是连续的； 3、双线性插值 双线性插值是一片一片的空间二次曲面构成。 双线性插值函数的形式如下： 其中有四个待定系数，利用该函数在矩形的四个顶点（插值节点）的函数值，得到四个代数方程，正好确定四个系数。 三、用MATLAB作网格节点数据的插值 z=interp2(x0,y0,z0,x,y,’method’