插值与拟合总结.pdfVIP

插值与拟合 在生产和科学实验中，自变量 与因变量 间的函数关系 有时不能写出解析表达式，而只能得到函数在若干点的函数 值或导数值，或者表达式过于复杂而需要较大的计算量。当 要求知道其它点的函数值时，需要估计函数值在该点的值。 为了完成这样的任务，需要构造一个比较简单的函数 ，使 函数在观测点的值等于已知的值，或使函数在该点的导数值 等于已知的值，寻找这样的函数 有很多方法。根据测量数 据的类型有以下两类处理观测数据的方法。 （1）测量值是准确的，没有误差，一般用插值。 （2 ）测量值与真实值有误差，一般用曲线拟合。 在 MATLAB 中，无论是插值还是拟合，都有相应的函数来 处理。 一、插 值 1、一维插值： 已知离散点上的数据集 ，即已知在点集 X= 上的函数值 Y= ，构造一个解析函数（其图形为一曲线）通过这些点， 并能够求出这些点之间的值，这一过程称为一维插值。 MATLAB 命令： yi=interp1(X, Y, xi, method) 该命令用指定的算法找出一个一元函数 ，然后以 给出 处的 值。 xi 可以是一个标量，也可以是一个向量，是向量时，必 须单调， method 可以下列方法之一： ‘nearest：最近邻点插值，直接完成计算；’ ‘spline：三次样条函数插值；’ ‘linear：线性插值（缺省方式）’ ，直接完成计算； ‘cubic：三次函数插值；’ 例 1：已知某产品从 1900 年到 2010 年每隔 10 年的产量 为： 75.995, 91.972, 105.711, 123.203, 131.699, 150.697, 179.323, 203.212, 226.505, 249.633, 256.344, 267.893 ， 计算出 1995 年的产量，用三次样条插值的方法，画出每隔 一年的插值曲线图形，同时将原始的数据画在同一图上。 解：程序如下 year=1900:10:2010; product=[75.995, 91.972, 105.711, 123.203, 131.699, 150.697, 179.323, 203.212, 226.505, 249.633, 256.344, 267.893] p1995=interp1(year,product,1995) x=1900:2010; y=interp1(year,product,x,cubic); plot(year,product,o,x,y); 计算结果为： p1995=252.9885 。 2 、二维插值 已知离散点上的数据集 ，即已知在点集 上的函数值 ， 构造一个解析函数（其图形为一曲面）通过这些点，并能够 求出这些已知点以外的点的函数值，这一过程称为二维插 值。 MATLAB 函数： Zi=interp2(X,Y ,Z,Xi,Yi,method) 该命令用指定的算法找出一个二元函数 ，然后以 给出 处 的值。 X 是 m 维向量，指明所给数据网格点的横坐标， y 是 n 维向量， 指明所给数据网格点的纵坐标， z 是 mxn 维矩阵， 标明相应于所给数据网格点的函数值。向量 xi ，yi 是给定的