通过 GPML 实现 4-D 高斯过程回归的问题

Question

我参考链接https://stats.stackexchange.com/questions/105516/how-to-implement-a-2-d-gaussian-processes-regression-through-gpml-matlab并创建一个二维高斯过程回归。 我想创建一个 4-d 高斯过程回归，但是 'meshgrid' 只允许 3 个inputs([X,Y,Z] = meshgrid(x,y,z)) ; 如何在meshgrid中添加另一个输入？

3-d 代码是这样的：

X1train = linspace(-4.5,4.5,10);
X2train = linspace(-4.5,4.5,10);
X3train = linspace(-4.5,4.5,10);
X = [X1train' X2train' X3train'];
Y = [X1train + X2train + X3train]';

%Testdata
[Xtest1, Xtest2, Xtest3] = meshgrid(-4.5:0.1:4.5, -4.5:0.1:4.5, -4.5:0.1:4.5);
 Xtest = [Xtest1(:) Xtest2(:) Xtest3(:)];

% implement regression 
[ymu ys2 fmu fs2] = gp(hyp, @infExact, [], covfunc, likfunc, X, Y, Xtest);

如果我创建了一个 X4train，那就意味着我需要一个 Xtest4，我该如何将 Xtest4 添加到 meshgrid 中？

GPML 代码来自http://www.gaussianprocess.org/gpml/code/matlab/doc/

Answer 1

您可以使用ndgrid创建 n 维网格，但请记住，它不会直接创建与 meshgrid 相同的输出，您必须先对其进行转换。 （文档中也解释了如何做到这一点）