MATLAB拟合关系总结,matlab中函数拟合方法—个人总结

matlab中函数拟合方法—个人总结

目 录 一维插值方案 .2 二维数据内插值(表格查找).3 等高线 .4 三维曲面 .5 等高线 2 6 三维曲面 2 7 matlab 绘制温度场(尚未深入研究).13 二维曲线(非线性)拟合步骤 .18 三维曲线(非线性)拟合步骤 .19 三维曲线的画法 .20 三维曲面的画法 .21 画三维图 3 只有点的数据,没有函数关系式 .23 空间点拟合的基本原理 .27 空间点拟合的最小二乘法 .28 曲面生成后再进行多项式拟合 .37 六点生成曲面 .38 四点生成平面 .39 用三维离散点拟合光滑曲面 1 40 用三维离散点拟合光滑曲面 2 40一维插值方案 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) 1900 1920 1940 1960 1980 2000 2020 50 100 150 200 250 300 插值说明: interp1(x,Y,xi,) %用指定的算法计算插值： ’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,, extrap ) %对于超出 x 范围的 xi 中的分量将执行特殊的外插值法 extrap。 yi = interp1(x,Y,xi,,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) 0 1 2 3 4 5 6 7 8 9 10 -6 -4 -2 0 2 4 6 8 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-4 -2 0 2 4 -4 -2 0 2 4 -10 0 10 20 30 功能 三维数据插值 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]); shadi