Assuming that you want to interpolate between the corresponding elements of those 4x4 matrices, then you could do something like this:
C = {randi(9,4,4), randi(9,4,4), randi(9,4,4)}; % 4x4 matrices in a cell array
C{:}
L = [1;3;2]; % load
W = [7;8;9]; % omega
M = cat(3,C{:});
[X,Y] = meshgrid(1:numel(C),1:16);
F = scatteredInterpolant(Y(:),L(X(:)),W(X(:)),M(:),'linear','linear');
Interpolate the entire matrix:
Lq = 1.5;
Wq = 7.5;
out = reshape(F({1:16,Lq,Wq}),4,4)
Lq = 2;
Wq = 8;
out = reshape(F({1:16,Lq,Wq}),4,4)
