MATLAB插值_拟合全.pptVIP

Matlab差值拟合 plot(X1,Y1,LineSpec1,X2,Y2,LineSpec2…) 将按顺序分别画出由三参数定义Xi,Yi,LineSpeci的线条。其中参数LineSpeci指明了线条的类型，标记符号，和画线用的颜色。 图形窗口的控制命令 命令1 subplot 用法 ：subplot(m,n,p) 将一图形窗口分成m*n个小窗口，在第p个小窗口中创建一坐标轴。 命令2 hold 用法 ：hold on 保留当前图形与当前坐标轴 ，后面的图形命令只能在当前存在的坐标轴中增加图形 。不覆盖前面图形 。 hold off 覆盖前面图形 。 举例 t = 0:pi/20:2*pi; plot(t,t.*cos(t),-.r*) hold on plot(exp(t/100).*sin(t-pi/2),--mo) plot(sin(t-pi),:bs) hold off 插值 自变量x与因变量y的函数y = f(x)的关系式有时不能直接写出表达式，而只能得到函数在若干个点的函数值。当要求知道观测点之外的函数值时，需要估计函数值在该点的值。 插值命令 命令1 interp1 功能 一维数据插值。该命令对数据点之间计算内插值。它找出一元函数f(x)在中间点的数值。其中函数f(x)由所给数据决定。 格式 yi = interp1(x,Y,xi)%返回插值向量yi，每一元素对应于参量xi，同时由向量x与Y的内插值决定。 yi = interp1(x,Y,xi,method) %用指定的算法计算插值 ’nearest’：最近邻点插值，直接完成计算； ’linear’：线性插值（缺省方式）。 ’spline’：三次样条函数插值。 ’pchip’：分段三次Hermite插值。 举例 x = 0:10; y = x.*sin(x); xx = 0:.25:10; yy = interp1(x,y,xx); plot(x,y,kd,xx,yy) year = 1900:10:2010; product = [75.995 91.972 105.711 123.203 131.669 150.697 179.323 203.212 226.505 249.633 256.344 267.893 ]; p1995 = interp1(year,product,1995) x = 1900:1:2010; y = interp1(year,product,x,pchip); plot(year,product,o,x,y) 结果 p1995 = 252.9885 曲线拟合： 根据n个离散的点，拟合出一条曲线y=F(x)，每个点到F(x）的距离两两相乘的积最小。 多项式拟合 命令：polyfit 用法：polyfit(x, y, n)描述n阶多项式y=f(x)的数据进行最小二乘曲线拟合 x=[0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1]; y=[-.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2]; n=2; p=polyfit(x, y, n) p = -9.8108 20.1293 -0.0317 其解是y = －9.8108x2 ＋20.1293x－0.0317。 为了将曲线拟合解与数据点比较，让我们把二者都绘成图。 xi=linspace(0, 1, 100); z=polyval(p, xi); plot(x, y, o , x, y, xi, z, : ) 非线性数据（曲线）拟合 命令：lsqcurvefit 用法：x = lsqcurvefit(fun,x0,xdata,ydata) 举例： x = [3.6 7.7 9.3 4.1 8.6 2.8 1.3 7.9 10.0 5.4]; y = [16.5 150.6 263.1 24.7 208.5 9.9 2.7 163.9 325.0 54.3]; 求系数 先建立拟合函数文件，并保存为myfun.m function F = myfun(x,xdata) F = x(1)*xdata.^2 + x(2)*sin(xdata) + x(3)*xdata.^3; 输入数据 xdata = [3.6 7.7 9.3 4.1 8.6 2.8 1.3 7.9 10.0 5.4]; ydata = [16.5 150.6 263.1 24.7 208.5 9