C=dsft(F);
dsft(F) computes the discrete symplectic Fourier transform of F. F must be a matrix or a 3D array. If F is a 3D array, the transformation is applied along the first two dimensions.
Let F be a \(K \times L\) matrix. Then the DSFT of F is given by
for \(m=0,\ldots,L-1\) and \(n=0,\ldots,K-1\).
The dsft is its own inverse.
H. G. Feichtinger, M. Hazewinkel, N. Kaiblinger, E. Matusiak, and M. Neuhauser. Metaplectic operators on cn. The Quarterly Journal of Mathematics, 59(1):15--28, 2008.