As of R2020b, there is no available method in MATLAB that allows the interpolation a multiinput-multioutput function based on scattered data. 'scatteredInterpolant' can only be used for up to 3 dimensions.
This workflow can be achieved by training a deep neural network on the provided pairs of tuples exploiting the 'featureInputLayer' for the input layer of the corresponding network. An example is provided in the following, that should serve as a demonstration and should be adjusted for a particular training task:
%% Number of input and output variables
numInput = n1;
numOutput = n2;
%% Definition of the training data
x = A(:, 1:numInput);
y = A(:, numInput + 1:end);
%% Definition of a sample layer architecture
layers = [ ...
featureInputLayer(numInput, "Name", "myFeatureInputLayer")
fullyConnectedLayer(8, "Name", "myFullyConnectedLayer1")
tanhLayer("Name", "myTanhLayer")
fullyConnectedLayer(numOutput, "Name", "myFullyConnectedLayer2")
regressionLayer("Name", "myRegressionLayer")
];
%% Definition of sample training options
opts = trainingOptions('adam', ...
'MaxEpochs', 1000, ...
'InitialLearnRate', 0.01,...
'Shuffle', 'every-epoch', ...
'Plots', 'training-progress', ...
'MiniBatchSize', 128, ...
'Verbose', false);
%% Training of the network
[trainedNet, info] = trainNetwork(x, y, layers, opts);
%% Compute a predicted output tuple on a new input tuple
% Please replace 'xHat' with the input vector of the actual application case
% in order to be able to run the following commands
% xHat = [xHat1 xHat2 ... xHatn1];
% yHat = predict(trainedNet, xHat)
