end end % Reshape the Pseudo-Inverse to match "D" so "C(:,:,1)" provides a convolution % kernel which performs Savitzky-Golay filtering. % "C" will be a 3-dimensional matrix. C = reshape(C ...