《插值与拟合.ppt

多项式拟合的Matlab命令：Polyfit函数 注意：最邻近插值一般不连续。具有连续性的最简单的插值是分片线性插值。 最邻近插值 x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? y (x1, y1) (x1, y2) (x2, y1) (x2, y2) O 二维或高维情形的最邻近插值，与被插值点最邻近的 节点的函数值即为所求。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 将四个插值点（矩形的四个顶点）处的函数值依次简记为： 分片线性插值 x y ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? (xi, yj) (xi, yj+1) (xi+1, yj) (xi+1, yj+1) O f (xi, yj)=f1，f (xi+1, yj)=f2，f (xi+1, yj+1)=f3，f (xi, yj+1)=f4 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 插值函数为： 第二片(上三角形区域)：(x, y)满足 插值函数为： 注意：(x, y)当然应该是在插值节点所形成的矩形区域内。显然，分片线性插值函数是连续的； 分两片的函数表达式如下： 第一片(下三角形区域)： (x, y)满足 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 双线性插值是一片一片的空间二次曲面构成。 双线性插值函数的形式如下： 其中有四个待定系数，利用该函数在矩形的四个顶点（插值节点）的函数值，得到四个代数方程，正好确定四个系数。 双线性插值 x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? y (x1, y1) (x1, y2) (x2, y1) (x2, y2) O 返回 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 二维插值 对二维插值问题，MATLAB分别给出了针对插值基点为网格节点的插值函数及针对插值基点为散乱节点的插值函数调用格式。 1、 用MATLAB作网格节点数据的插值 对上述问题，MATLAB提供了二维插值函数interp2,其基本格式为： Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. z=interp2(x0,y0,z0,x,y,’method’) 其中x0,y0是自变量。X0,y0的分量值必须是单调递增的。X0和y0分别是m维和n维向量，分别表示已知数据点的横、纵坐标向量，z0是m*n维矩阵，标明相应于所给数据网格点的函数值。向量x,y是待求函数值所给定网格点的的横、纵坐标向量，x,y的值分别不能超出x0,y0的范围。 而method为可选参数，有四种选择： ‘nearest’ 最邻近插值 ‘linear’ 线性插值 ‘spline’ 三次样条插值 ‘cubic’ 三次插值 缺省时, 是线性插值 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 例4：测得平板表面3×5网格点处的温度分别为： 82 81 80 82 84 79 63 61 65 81