Matlab二维数据插值

Matlab二维数据插值

Matlab可以进行二维数据插值，这里主要介绍interp2和griddata命令。

Matlab具有强大的数值分析功能，这里介绍一下它的二维数据插值。

首先看看interp2的功能，它需要三个二维数组作为初始值，然后在它的数据范围内进行插值。

可以对其范围内的数组和数值进行插值。

[X,Y] = meshgrid(-3:3);V = peaks(X,Y);%给出XYV的值figuresurf(X,Y,V);%作图

可以进行数据插值，下面使用默认方法进行插值：

[xx,yy] = meshgrid(-3:0.1:3);>> vv = interp2(X,Y,V,xx,yy);>> figure;surf(xx,yy,vv)

%可见除了增加数据量，其他没什么变化。

下面使用cubic插值方法：

vv1 = interp2(X,Y,V,xx,yy,'cubic');>> figure;surf(xx,yy,vv1)

%可见数据质量有了明显改善。另外还有nearest，spline等插值方法供选择。

下面介绍griddata，griddata可以在已知XYZ的情况下使用，也就是我们如果只知道三列数值，而想要对这些值范围内的数据进行插值，就可以用这个命令。

架设我们上述X Y V均为一维数据；

aa(1:49,1) = X(1:49);>> aa(1:49,2) = Y(1:49);>> aa(1:49,3) = V(1:49);

我们将XYV全部放在数组aa中，然后使用上述方法对xx和yy进行插值：

vv3 = griddata(aa(:,1),aa(:,2),aa(:,3),xx,yy);

figure;surf(xx,yy,vv3)

可以看到与上述interp2命令得到的结果基本相同，当然也可以进行cubic插值：

vv4 = griddata(aa(:,1),aa(:,2),aa(:,3),xx,yy,'cubic');

figure;surf(xx,yy,vv4)

%效果与interp2的cubic是相同的（如下图），但是这可以用于只有三列数据的情况。