1

我有两个 3D 数组： 

 A=[
12751   4616    15915
15864   4622    15667
12877   4683    21050
15816   4668    21253
11374   5006    18495
16995   5466    18493
11638   4880    20023
17078   4938    20006
11576   4886    17011
];

和

 B=[
12402   2138    15743
10285   3175    15851
10237   3084    21052
12130   2129    21299
8074    3802    18505
14515   1623    18497
8415    3713    19856
14462   1120    20061
8340    3711    17145
14483   1157    16990];

我想在 Matlab 中使用它们之间的样条曲线或 3D 插值。我Interp3应该如何定义？VVI = interp3(X,Y,Z,V,XI,YI,ZI)

或Y在yy = spline(x,Y,xx).

4

1 回答 1

0

我不太了解您的数据，是（或反之亦然）B的函数吗？A此外，这些阵列似乎是 1D，而不是 3D。你能澄清一下吗？

在你的函数调用中

yy = spline(x, Y, xx)

Y是您要插值的因变量：Y是的函数，x上述函数调用的结果是返回Yat的值xx。作为一维的示例，请尝试

x = linspace(0., 2.*pi, 100);
Y = sin(x);

% What is the value of Y (or sin(x)) at pi/2?
xx = pi/2.;
yy = spline(x, Y, xx); % This should result in yy = 1.000

查看spline文档以获取更多信息和使用此功能的示例。

现在这个函数只适用于 1D 拟合，并且（我认为）等价于yy = interp1(x, Y, xx, 'spline'). 如果要进行三维查找，则必须使用interp3，它将上述示例推广到 3D。因此，除了一个独立的坐标 之外x，我们还有两个y、z和 三个坐标用于我们要执行查找的点xx：yy和zz。您正在插值的函数必须是坐标的 3D 函数(x, y, z)。尝试，例如：

x = linspace(-1., 1., 100); y = x; z = x;
[X, Y, Z] = meshgrid(x, y, z);
s = exp(-sqrt(X.^2 + Y.^2 + Z.^2));

sinterp = interp3(x, y, z, s, 0., 0., 0.) % Should give sinterp = 0.9827
于 2012-04-27T15:24:54.397 回答