用Matlab插值

Question

我有一个具有不同值的向量。 一些值是零，有时它们甚至一个接一个地出现。 我需要针对另一个具有相同大小的向量绘制该向量，但是不能包含零。 我可以对向量进行某种插值的最佳方法是什么，怎么做？ 我尝试阅读有关mat-lab的插值的知识，但我对实现插值的理解不够好。 如果有可能逐步向我解释它，我将不胜感激，因为我是这个程序的新手。

谢谢

Answer 1

从包含两个等长向量x ， y的数据集开始，其中将排除y等于零的y值，首先选择不包括零的子集：

incld = y~=0;

然后，您可以对该子集进行插值：

yn = interp1(x(incld),y(incld),x);

实施例的结果，绘制x针对y （绿色）和x对yn （红色）：

编辑

请注意，根据插值的定义，如果终端点为零，则必须单独进行处理，例如，通过在上一行之前运行以下命令：

if y(1)==0, y(1) = y(find(y~=0,1,'first'))/2; end
if y(end)==0, y(end) = y(find(y~=0,1,'last'))/2; end

编辑＃2

这是上面的2D版本，其中数组X和Y是与2D数组Z的条目相对应的坐标：

[nr nc]=size(Z);
[X Y] = meshgrid([1:nc],[1:nr]);
X2 = X;
Y2 = Y;
Z2 = Z;
excld = Z==0;
X2(excld) = [];
Y2(excld) = [];
Z2(excld) = [];
ZN = griddata(X2,Y2,Z2,X,Y);

ZN包含插值点。

在下图中，深蓝色补丁显示零。 左边在插值之前，右边在插值之后：