Quaternion Fast Fourier transform
Y = qfft2(X, A, L)
qfft2(X, A, L) computes the Quaternion Fast Fourier Transform of the quaternion matrix X using transform axis A (direction in 3-space).
L specifies the handedness of the transform ('L' or 'R') - determined by the position of the complex exponential relative to X. ('L' gives a transform with the exponential on the left of the signal.) See the related function fft2 which supplies a default axis and handedness.
The transform axis, A must be a pure quaternion (real or complex) but it need not have unit modulus (although the transform itself is computed using a unit-modulus axis, so that the axis is a root of -1).
This function uses the MATLAB® fft function to compute two or four complex FFTs depending on whether X and A are real or complex. If either or both are complex, a complex quaternion FFT is computed.