在细化网格的 3D 网格上插入标量场

Question

我有函数 f=f(x,y,z) 在大小为 MxNxP 的点网格上进行数值计算。

我需要在（显然）不同维度的更精细网格上插入这些值。就像我对 interp1 所做的那样，当我在一个区间内只有一些函数值并且我想对其进行插值以找到更多点的值。

我尝试执行此代码：

XX=linspace(0,45,50);
YY=linspace(0,0.5,50);
ZZ=linspace(0,0.5,50);

U_1 = interp3(X,Y,Z,f,XX,YY,ZZ,'linear')

其中 X Y Z 是我具有 f 值的点的向量。 XX,YY,ZZ 是我想要解决的点。

我尝试使用 griddata 和 interp3，但收到以下错误消息：

使用 griddedInterpolant 时出错
网格向量不定义与给定值匹配的点网格。
interp3 中的错误（第 133 行）
F = griddedInterpolant({X, Y, Z}, V, method,extrap);

你能帮帮我吗？

Answer 1

您似乎有一个不规则的网格数据集。在这种情况下，您需要使用griddata() 函数进行插值。

请注意您定义的点数组 X、Y、Z 具有相同的维度。插值点数组 XX,YY,ZZ 也必须具有相同的维度。 这意味着以下表达式必须是true（或1）：

all(size(X)==size(Y)) && all(size(X)==size(Z)))

为了创建插值，您可以使用 meshgrid。将插值网格的范围调整到实际数据范围是明智的：

[XX,YY,ZZ]=meshgrid(linspace(min(X(:)),max(X(:)),50), ...
                    linspace(min(Y(:)),max(Y(:)),50), ...
                    linspace(min(Z(:)),max(Z(:)),50));

这可能会得到您期望的结果。