Inverse fast Mellin transform
X = iffmt(MELLIN, BETA) X = iffmt(MELLIN, BETA, M)
a complex vectorof size N: the Mellin transform to be inverted.
Mellin must have been obtained from FMT with frequency running from FMIN to 0.5 Hz.
a real vector: the Mellin variable issued from FMT (only the two first elements are used).
a positive integer: the number of points of the inverse Mellin transform. (default : N).
a complex column vector of size M: the inverse Mellin transform.
iffmt computes the inverse fast Mellin transform of MELLIN. WARNING : the inverse of the Mellin transform is correct only if the Mellin transform has been computed from FMIN to 0.5 Hz, and if the original signal is analytic.
sig = atoms(128,[64,0.25,32,1]); [MELLIN,BETA] = ffmt(real(sig),0.005,0.5,300); X = iffmt(MELLIN,BETA,128); clf; plot(real(sig),'g'); plot(real(X)); legend(_("initial"),_("reconstructed")) | ![]() | ![]() |