Quaternion Fast Fourier transform
Y = qfft(X, A, L)
qfft(X, A, L) computes the Quaternion Fast Fourier Transform of the quaternion vector X using transform axis A (direction in 3-space). If X is a matrix, the transforms of the columns are computed.
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 fft 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 is real or complex. If either or both are complex, a complex quaternion FFT is computed.