matlab成功的函数拟合方法_更新中.docVIP

一维插值方案 clear 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) %使用一维数据内插值(该题中只能在1900和2010之间进行插值,大于2010和小于1900都%无效)命令 x = 1900:1:2010 y = interp1(year,product,x,spine); plot(year,product,o,x,y) 插值说明: interp1(x,Y,xi,method) %用指定的算法计算插值： ’nearest’：最近邻点插值，直接完成计算； ’linear’：线性插值（缺省方式），直接完成计算； ’spine’：三次样条函数插值。对于该方法，命令interp1 调用函数spline、ppval、mkpp、umkpp。这些命令生成一系列用于分段多项式操作的函 数。命令spline 用它们执行三次样条函数插值； ’pchip’：分段三次Hermite 插值。对于该方法，命令interp1 调用函数pchip，用于对向量x 与y 执行分段三次内插值。该方法保留单调性与 数据的外形； ’cubic’：与’pchip’操作相同； ’v5cubic’：在MATLAB 5.0 中的三次插值。 对于超出x 范围的xi 的分量，使用方法’nearest’、’linear’、’v5cubic’的插值算法，相应地将返回NaN。对其他的方法，interp1 将对超出的分量执行外插值算法。 yi = interp1(x,Y,xi,method,extrap) %对于超出x 范围的xi 中的分量将执行特殊的外插值法extrap。 yi = interp1(x,Y,xi,method,extrapval) %确定超出x 范围的xi 中的分量的外插值extrapval，其值通常取NaN 或0。 例1 clear; x = 0:10; y = x.*sin(x); xx = 0:.25:10; yy = interp1(x,y,xx) plot(x,y,kd,xx,yy) interp2 功能 二维数据内插值（表格查找） [X,Y] = meshgrid(-3:.25:3); Z = peaks(X,Y); [XI,YI] = meshgrid(-3:.125:3); ZZ = interp2(X,Y,Z,XI,YI); surfl(X,Y,Z);hold on; surfl(XI,YI,ZZ+15) axis([-3 3 -3 3 -5 20]); shading flat hold off 功能 三维数据插值interp3（查表） [x,y,z,v] = flow(20); [xx,yy,zz] = meshgrid(.1:.25:10, -3:.25:3, -3:.25:3); vv = interp3(x,y,z,v,xx,yy,zz); slice(xx,yy,zz,vv,[6 9.5],[1 2],[-2 .2]); shading interp;colormap cool clear Z=peaks for w=1:1:100 V=[w/10,0,w/10] contour(Z,V) %C=contour(Z,V) %Clabel(C) Hold on title(等高线及其标注) end end x=0:10 y=0:.1:1 [d,B]=meshgrid(x,y) z=1./(B.*d.^2+1); surf(B,d,z) x=0:0.05:10 y=0:0.05:1 [X,Y]=meshgrid(x,y) Z=( X.^3+ 3.*Y.^2+5*Y); %Z=( X.^2+ 3.*Y.^3+5*Y);% surf(X,Y,Z) %一张普通的三维曲面,有时需要旋转一下才能看到下图的结果; x=0:0.05:1 y=0:0.05:1 [X,Y]=meshgrid(x,y) Z=( X.^2-Y.^2);% Z=( 4*X.^3*Y-4*X.*Y.^3); surf(X,Y,Z) %一张普通的三维曲面,有时需要旋转一下才能看到下图的结果; clear x=-2:0.1:2 y=-2:0.1:2 [X,Y]=meshgrid(x,y) Z=(X.^2+Y.^2).^0.5 for w=1:1:100 V=[w/3,w/pi,w/3] contour(Z,V) hold on end