3D DFT matrix

We know that dft can be achieved by the dft matrix, in 2D situation, fft(x) equals to F x F_T, where F is given by dftmtx(size(x)), so I want to know if it comes to 3D situation, how can we get such a matrix to achieve 3D dft?