DFRFT: A Classified Review of Recent Methods with Its Application

In the literature, there are various algorithms available for computing the discrete fractional Fourier transform (DFRFT). In this paper, all the existing methods are reviewed, classified into four categories, and subsequently compared to find out the best alternative from the view point of minimal computational error, computational complexity, transform features, and additional features like security. Subsequently, the correlation theorem of FRFT has been utilized to remove significantly the Doppler shift caused due to motion of receiver in the DSB-SC AM signal. Finally, the role of DFRFT has been investigated in the area of steganography.


Introduction
Due to the inadequate performance of Fourier transform for the analysis of nonstationary signal, normally encountered in communication systems when either transmitter or receiver is moving, fractional Fourier transform (FRFT) can be utilized.The concept of FRFT was first introduced by N. Wiener in 1929.However, the FRFT was recognized as a "transform method" by mathematical bodies after the work of Victor Namias in 1980 in which the concept of FRFT had been introduced by considering fractional power of eigenfunctions of the ordinary FT.The mathematical description of FRFT was given by McBride and Keer in 1987.And finally, a general definition of FRFT for all classes of signals was given by Cariolaro et al.To compute the FRFT of any signal its discrete version was needed which initiates the work for defining discrete FRFT (DFRFT).However, due to availability of different approaches for evaluating DFRFT, a detailed study is needed in the context of suitable definition of DFRFT.This study along with comparison between the existing methods has been included in Section 3.This paper includes five sections.The basic concept of fractional Fourier transform is described in Section 2. The detailed study of discrete FRFT (DFRFT) along with the comparative analysis between the different methods of evaluating DFRFT is included in Section 3. Subsequently, the utility of FRFT in DSB-SC signal detection under noisy condition and the application of DFRFT in steganography are being dealt in Section 4. Finally, conclusions are made in Section 5.

Fractional Fourier Transform (FRFT)
The FRFT is a mathematical tool which maps a signal from one domain to other in time-frequency plane (e.g., FRFT becomes FT when it maps a signal in time domain to frequency domain), shown in Figure 1.The FRFT of a signal () with angle parameter , represented by   () is defined as [1]: The angle parameter associated with the FRFT is given as: When  = 0, the FRFT becomes identity operator, and for  = /2 ( = 1), the FRFT becomes Fourier transform (FT).
In signal processing, cross-correlation is a measure of similarity of two waveforms as a function of a time lag applied to one of them.The cross-correlation theorem for the FRFT is given as follows [2].Let   () be defined as weighted crosscorrelation of two functions () and (), and   (),   (), and  , (), be defined as FRFT of (), (), and   (), respectively.Then, ()  ( + )  (+)cot   ←→  , () = √ 2 1 −  cot   −(/2) 2 cot    (−)   () . ( This cross-correlation theorem of FRFT is utilized for detection of DSB-SC signal when the signal encounters a Doppler shift due to motion of transmitter or receiver.

Discrete Fractional Fourier Transform (DFRFT)
Many researchers have tried to come up with some method of evaluating DFRFT.On this journey, many methods and algorithms had come into existence.But none got acclamation due to deficiency of either not satisfying some of the prime properties that its continuous type possessses or nonexistence of a closed-form expression.Previously, many classifications have already been documented [3].However, in this section, all the approaches available in the literature are included and divided into four classes on the basis of their methodology of evaluation.These classes are (i) sampling-based method (ii) linear combination method (iii) eigenvector-based method (iv) weighted summation-based method.
3.1.Sampling-Based Method.In this approach, by considering the samples of the kernel matrix under some predefined constraint, two different algorithms are available in the literature and they are named as sampling type DFRFT and closed-form type DFRFT.These two methods are described below: 3.1.1.Sampling Type DFRFT.The sampling theorem for the FRFT of band-limited and time-limited signals is followed from Shannon's sampling theorem [4][5][6][7][8][9].With this approach, the simplest way to derive the DFRFT is sampling the continuous FRFT and computing it directly from the samples, but by sampling the continuous FRFT directly, the resultant discrete transform obtained will lose many important properties.
The most serious problem with this type DFRFT will be noncompliance of unitary and reversibility properties.In addition, the DFRFT obtained by direct sampling of the FRFT lacks closed-form expressions and is nonadditive, which affects its application domain.In order to maintain some of the FRFT properties, a type of DFRFT was derived as a special case of the continuous FRFT as reported in [10].Specifically, the input function was assumed as a periodic, equally spaced impulse train.Since this type of DFRFT is a special case of continuous FRFT, many properties of the FRFT exist and have the fast algorithm.However, this type of DFRFT cannot be defined for all values of  due to various imposed constraints.At the same time, Ozaktas et al. [11] had proposed two innovative approaches for obtaining the DFRFT through the sampling of the FRFT.Both methods were based on the idea of manipulating the expression of the FRFT by sampling appropriately.
In the first method, the DFRFT of a signal () can be obtained by multiplying the signal with a chirp function as followes: For this chirp multiplication, the bandwidth and timebandwidth product of   () can be as large as twice of ().Thus, the sampling of   () will be done at intervals of (1/2Δ).If the samples of () are spaced at (1/Δ), then this needs to be interpolated before multiplying it by the samples of the chirp function to obtain the desired samples of   ().
The next step was to convolve the signal   () with another chirp function.To perform this convolution, the chirp signal was replaced by its band-limited version because   () is assumed as a band-limited signal.
The band-limited version of chirp function can be expressed in the form of Fresnel integral.Finally, again performing a chirp multiplication, the DFRFT of signal () can be obtained as: The convolution operation in the above equation can be achieved by sampling () and performing the convolution using the fast Fourier transform (FFT).Overall, the procedure starts with  samples spaced at (1/Δ), which uniquely characterizes the function () and returns the same for   .For  and   , denote column vectors with  elements containing the samples of () and its DFRFT, respectively.The above procedure can be given in a matrix notation as: Here,  and  are matrices representing the decimation and interpolation operations, Λ is a diagonal matrix that corresponds to chirp multiplication, and  corresponds to the convolution operation.The above method had allowed obtaining the samples of the th transform in terms of the samples of the original function.
First method includes the analysis of one step in terms of Fresnel integral, which requires large number of computation time particularly in inverse DFRFT.To remove this discrepancy, a second method was introduced by Ozaktas et al. [11].In this method, one assumption was taken that the Wigner distribution of function () could be zero outside a circle of diameter Δ centered at the origin and by limiting the order "" to the interval 0.5 ≤ || ≤ 1.5.Simultaneously, the amount of vertical shear in Wigner space resulting from the chirp modulation is bounded by Δ/2.The FRFT can be written as: Thus,   2 cot  () can be written by using Shannon's interpolation formula as: where  = (Δ) 2 ; putting (9) in (8) and interchanging the order of integration and summation, )) . ( And after some algebraic manipulation, (10) can be restructured as: ) . ( This summation can be interpreted as the convolution of  (/2Δ) 2 csc  and the chirp modulated function  (cot −csc )(/2Δ) 2 (/2Δ).This convolution can be computed by fast Fourier transform (FFT).And the overall methodology can be written in matrix form as: where the   matrix was given as: In the same article, Ozaktas et al. [11] gave alternate methods to calculate the DFRFT by introducing scaling process.All these methods resulted in the same time complexity.However, in practice, these alternate methods are not preferable because they require coordinate scaling.Scaling is not an advantageous process if the original data has already been sampled.This is because scaling requires additional interpolations and so forth, which will require additional computation.
However, in both cases, the Wigner distribution of the signal was confined to a circle of diameter Δ around the origin.Thus, there might be several discrete fractional Fourier transform matrices, some of which were more elegantly expressible than the others, which gave the same result within the accuracy of this approximation.All these matrices would yield results that were in increasingly better agreement as the signal energy contained in the circle was increasingly closer to the total energy.[12] had derived a different type of DFRFT which was neat in concept and reversible but lacks the additivity property.Due to the orientation of practical usage, two types of DFRFT were derived which are different in parameterization.The parameters of the first type have directly linked to the continuous FRFT and suit the applications of computing the continuous FRFT.On the other hand, the second type has the simpler parameters set and allows more elegant expression for the operator kernels and more is suitable for other applications of DFRFT, such as the filter design, pattern recognition, and phase retrieval.Considering the expression of continuous FRFT:

Closed-Form Type DFRFT. Pei and Ding
And to derive the DFRFT, sampling of () and   () with the interval Δ and Δ needs to be performed as: where  = −, . . .,  and  = −, . . ., .By using this sampling, FRFT can be converted [13] as: The above can be written as: where And   () can be reversible only if the inverse transform of   (, ), for all  ≥ , would be Hermitian (conjugate and transpose), given [12] as: From ( 17) and (19), After some manipulation, The summation for  in (21) should be equal to ( − ), then where || is some integer prime to (2 + 1); in this case, (18) becomes After normalization, the transform matrix [12] can be given as: The definition of DFRFT is further divided into two parts [12].The first part was defined for "sin  > 0" as: And the second part was defined for "sin  < 0" as: This definition of DFRFT possesses reversibility property and the periodicity property: The DFRFT of type 1 has a very important advantage in terms of its efficiency in calculating and implementing the DFRFT [12].Due to two chirp multiplications and one FFT required for the implementation of this type of DFRFT, the total number of the multiplication operations required was "2 + (/2)log 2 , " for length of the output as  = 2 + 1.
Subsequently, another type of DFRFT, named as type 2, was introduced by modifying the definitions obtained for discrete affine Fourier transform (DAFT) and presented as: For the DFRFT of type 2, the parameter "" is used to control the variation of the chirp in frequency domain [13].The DFRFT of type 2 also needed "2 + (/2)log 2 " multiplication operations.It does not follow the additivity property but observes the convertibility operation.
These definitions of DFRFT are better to previous ones, that is, sampling type, linear combination type, and eigenfunction type, as it require less number of multiplications, having less complexity.Also it has a closed form expression being rich in properties.

Linear Combination Type DFRFT.
In [14][15][16][17], the discrete fractional Fourier transform was derived by using the linear combination of identity operation ( 0 ), discrete Fourier transform ( /2 ), time inverse operation (  ), and inverse discrete Fourier transform ( 3/2 ).In [15], the concept of DFRFT was given by linear combination of Lagrange interpolation polynomial of degree 3.For any function (V) where And   (V) is the Lagrange interpolation polynomial of degree 3 taking value "one" at V =   and "zero" otherwise.It follows that Finally, by taking the principle th-root of eigenvalue   , the fractional transform can be given as: where Subsequently, in [14,15], it had also been shown that the operator defined by ( 17) is unitary; satisfies angle additivity and angle multiplicity properties, and the operator is periodic with a fundamental period of four.One analogous approach of this method was presented in [17], in which the DFRFT was calculated by first chirping the signal in time domain, then taking its Fourier transform, and finally chirping again in frequency domain.Depending on the character of the chirped signal, the embedded FT can be one of the two known classes in discrete case-discrete-time and periodic discretetime.
For discrete-time DFRFT, the expression for DFRFT of signal () was given as: And its inverse transform as: where   = 2/ csc  and   (, ) is the kernel of FRFT.
Similarly, for periodic discrete-time signal () with period   , the expression for DFRFT was given as: And its inverse transform as: where   =  and   =  = 2/ csc  In this type of DFRFT, the transform matrix is orthogonal and satisfies the additivity property along with the reversibility property.However, the main problem with this method is that the transform result is not matched to the continuous FRFT.Besides, it performed very similarly to the original Fourier transform or the identity operation but lost the important characteristic of fractionalization.

Eigenvector Type DFRFT.
In this classification, the basic approach of eigenvalues and eigenvectors of the kernel matrix of the FRFT is considered to define the discrete FRFT.Here, the eigenvalues of the discrete Fourier transform (DFT) are considered first and then by obtaining the powers of these eigenvectors of DFT matrix, the corresponding eigenvalues and eigenvectors of DFRFT matrix are derived.This classification includes three different approaches, named as, eigenvector decomposition method, random type method, and  method DFRFT.

3.3.1.
Eigenvector Decomposition Type DFRFT.The authors had derived another type of discrete fractional Fourier transform as reported in [13,[18][19][20][21][22][23][24][25] by searching the eigenvectors and eigenvalues of the DFT matrix followed by computing the fractional power of the DFT matrix based on Hermite function.This type of DFRFT worked very similarly to the continuous FRFT, and it fulfills the properties of orthogonality, additivity, and reversibility.The fractional power of matrix could be calculated from its eigendecomposition and the power of eigenvalues.Unfortunately, there exist two types of ambiguity in deciding the fractional power of the DFT kernel matrix.

(i) Ambiguity in Deciding the Fractional Powers of
Eigenvalues.The square roots of unity are "1" and "-1" from elementary mathematics.This indicates that there exists root ambiguity in deciding the fractional power of eigenvalues.(ii) Ambiguity in Deciding the Eigenvectors of the DFT Kernel Matrix.The DFT eigenvectors constitute four major eigensubspaces; therefore, the choices for the DFT eigenvectors to construct the DFRFT kernel are multiple and not unique.
Because of the above mentioned ambiguities, several DFRFT kernel matrices were possible which could obey the rotational properties.The idea for developing the DFRFT is to find the discrete form of the following: where  indicates the rotation angle associated with FRFT of signal and   () represents the th-order normalized Hermite function with unit variance.The th-order normalized Hermite function with variance  is defined as: where, ℎ  (⋅) represents the th-order Hermite polynomial.However, the normalized Hermite functions with unitary variance   (⋅) indicate the eigenfunction of the FRFT.With the help of Hermite functions, the FRFT of a signal can be computed [24] as: The FRFT can be interpreted as a weighting summation of Hermite functions, as given in (40).The weighting coefficients are obtained by multiplying the phase term  − and the inner product of the input signal with the corresponding Hermite functions [24].The development of this DFRFT was based on the eigendecomposition of the DFT kernel.The eigenvalues of DFT kernel "" can be given as {1, −, −1, }, and corresponding to each eigenvalue one eigenvector can be evaluated.The continuous FRFT has a Hermite function with unitary variance as its eigenfunction.The corresponding eigenfunction property for the DFT can be given as: where û represents the eigenvector of DFT corresponding to the th-order discrete Hermite functions.In order to retain the eigenfunction properties, the unit variance Hermite function is sampled with a period   = √2/.In the case of continuous FRFT, the terms of the Hermite functions are summed up from order zero to infinity.However, for the discrete case, only eigenvectors for the DFT Hermite eigenvectors can be added.
The selection of the DFT Hermite eigenvectors is usually made from low to high orders, due to small approximation error of the low DFT Hermite eigenvectors [24].
Finally, the transform kernel of DFRFT can be defined as: where for N is even, û represents the normalized eigenvector corresponding to the kth-order discrete Hermite function, and  can be defined as follows.
For  odd, And for  even, The convergence for the eigenvectors obtained from these matrixes is not so fast for the high-order Hermite functions.To ensure orthogonality of the DFT Hermite eigenvectors û , the Gram-Schmidt algorithm (GSA), or the orthogonal procrustes algorithm (OPA) can be used [24].The GSA minimizes the errors between the samples of the Hermite functions and the orthogonal DFT Hermite eigenvectors.On the other hand, the OPA minimizes the total errors between those samples.The main difference between the approach proposed in [24] and similar approaches proposed in [21][22][23] is found in the eigenvectors obtained.The eigenvectors obtained by method reported in [24] were discrete Mathieu function, although the Mathieu functions can converge to Hermite functions as obtained by the methods reported in [21][22][23].Given the number of points  and the rotation angle , the DFT Hermite eigenvectors can be computed followed by determining the eigenvalues of DFRFT.The computation of the DFRFT can be implemented only by a transform kernel matrix multiplication.The complexity of computing the DFRFT is ( 2 ) same as in the DFT case.By adjusting the rotational angles, the method for implementing the DFRFT of a signal  can be given as: The coefficients   's are the inner products of signal and eigenvectors, and they can be computed in advance.If rotation angle changed then only the diagonal matrix  2/ needs to be recomputed.Additionally, the authors [24] investigated the relationship between the FRFT and the DFRFT and found that for a sampling period equal to √2/, the DFRFT performs a circular rotation of the signal in the timefrequency plane.However, the DFRFT becomes an elliptical rotation in the continuous time-frequency plane for sampling periods different from √2/ [24].Therefore, for these elliptical rotations an angle modification and a postphase compensation in the DFRFT have been required to obtain results similar to the continuous FRFT [24].This approach has been extended to the so-called multiple-parameter discrete fractional Fourier transform (MPDFRFT) [20,25].In fact, the MPDFRFT maintains all desired properties and reduces to the DFRFT when all of its order parameters are the same.
A similar approach to [24] has been proposed in [19].However, authors in [19] believe that the discrete time counterparts of the continuous time Hermite-Gaussians maintained the same properties because these discrete time counterparts exhibited better approximations than the other proposed approaches [18].In order to resolve this issue about the approximation of Hermite-Gaussian functions, a nearly tridiagonal commuting matrix of the DFT and a corresponding version of the DFRFT were proposed in [13].Most of the eigenvectors of this proposed nearly tridiagonal matrix result in a good approximation of the continuous Hermite-Gaussian functions by providing a smaller approximation error in comparison to the previous approaches.[26] had introduced a modified method to calculate the DFRFT based on the random DFT eigenvectors and eigenvalues.In this method, a new commuting matrix with random DFT eigenvectors was constructed.Then, a random discrete fractional Fourier transform (RDFRFT) kernel matrix with random DFT eigenvectors and eigenvalues was developed.The magnitude and phase of the transformed output are found to be random in RDFRFT.This can be considered as a special feature of RDFRFT.Previously, the DFT was randomized to define the discrete random Fourier transform (DRFT) in [27] by taking random powers of eigenvalues of the DFT matrix.However, the eigenvectors of the DRFT were not random and were computed using the same method as proposed in [23].Although, the phase of the DRFT was random, its magnitude was not random, as shown in [27].The eigen-decomposition of DFT matrix "" was given as:

Random Type DFRFT. Pei and Hsue
The DFRFT with one order parameter "" was defined by [25]: The definition of RDFRFT was based on the definition of multiple-parameter DFRFT (MPDFRFT) [25], given as: where the 1 ×  parameter vector â was defined as: This order parameter of the MPDFRFT can be taken as independent random numbers to define an MPDFRFT with the following random eigenvalues [26],     ,  = 0, 1, . . .,  − 1, and random eigenvalues were given as: The generation of RDFRFT is given in the following section.In this generation of RDFRFT, the DFT-commuting matrix with orthonormal random DFT eigenvectors is used.
First defining a  ×  real random Matrix "" being "" symmetric and generated by  = , where  represents circular reversal matrix, which is given by: where,  −1 is the ( − 1) × ( − 1) reversal matrix, having nonzero entries "ones" placed on the antidiagonal locations.
The RDFRFT kernel matrix  â  satisfied various properties of the FRFT like-additivity, FT convertibility, and angle additivity.The computational complexity of the -point RDFRFT was ( 2 ) same as that of the MPDFRFT, but RDFRFT has an advantage over MPDFRFT in terms of one more free parameter associated with it as parameter vector â.The magnitude and phase responses of the transformed output, with the RDFRFT, are random because the eigenvectors and eigenvalues of the RDFRFT matrix are random [26], as shown in (52).This aspect establishes its advantageous role in security areas like-cryptography, encryption, and watermarking.

3.3.3.
The  Method DFRFT.In this work [28], authors focus on a different approach to derive eigenvectors of the DFT matrix in order to build the DFRFT matrix that mimics the properties of its continuous counterpart.A straightforward and refined derivation of the eigenvectors of the DFT matrix is proposed without using any commuting matrices, but using only the factored form of a basic property of the DFT matrix stating that four consecutive Fourier transforms is the identity transform where  4  and   are the centered DFT (CDFT) and the identity matrices of order , respectively.We utilize factorization of the CDFT matrix to obtain new matrices, whose columns are eigenvectors of the CDFT matrix.Thereafter, GSA was used to find the orthonormal eigenvectors obtained by factorization of the CDFT matrix.By employing a DFT-shift matrix , eigenvectors of the ordinary-DFT matrix are obtained by using the CDFT matrix.Later on, the  and  matrices together are also used to boost the performance.
In [24], Pei et al. compute the samples of the Hermite-Gauss functions and employ Gram-Schmidt algorithm (GSA) to orthogonalize these samples by projecting on the  matrix, which is equivalent to orthogonalizing the samples of the Hermite-Gauss functions without using the  matrix.Hanna et al. [29] claim to find the eigenvectors of the DFT matrix using a similar method.However, they use the ordinary-DFT matrix as a basis to their work, but Serbes and Durak-Ata [28] employs the CDFT matrix.Additionally, the eigenvectors are nonorthogonal and dissimilar to the samples of Hermite-Gauss functions.Serbes and Durak-Ata [28] work is different from these methods in the sense that the non-orthogonal eigenvectors of the CDFT matrix is determined without utilizing any  matrix or samples of the Hermite-Gauss functions, and then orthogonalization of eigenvectors is obtained by employing the GSA.
The DFT maps the signal from [0,  − 1] discrete input space to [0, 2] discrete frequency space, whereas the CDFT maps [−( − 1)/2, ( − 1)/2] to [−, ].This allows CDFT to define even and odd functions such as Hermite-Gauss functions, that is contrary to the ordinary DFT definition.In order to approximate the samples of the continuous FrFT and to imply the rotation property, Hermite-Gauss-like eigenvectors of the DFT matrix have to be obtained.Serbes and Durak-Ata [28] use the CDFT matrix instead of the ordinary DFT matrix, because the Hermite-Gauss-like vectors are eigenvectors of only the CDFT matrix, not the DFT matrix, that is, when a discrete Hermite-Gauss-like eigenvector is transformed using the ordinary DFT matrix, the output is the shifted version of the Hermite-Gauss vector multiplied by a complex sinusoidal.
The centered-DFRFT matrix can be obtained by calculating the Hermite-Gaussian-like eigenvectors of DFRFT as followes: where The above mentioned method satisfies unitary and angle additivity properties.Also, it reduces to the ordinary DFT when  = 1 and approximates to the samples of the continuous FRFT for fractional values of .
Performance of the Serbes and Durak-Ata [28] algorithm can be improved by using it in combination with  and  matrices.Pei et al. [13] uses the  method together with the Grunbaum's  matrix [30] as a linear combination, in which it has been observed that Hermite-Gauss-like eigenvectors of ( + 15) are more accurate than both  and  alone.Serbes and Durak-Ata [28] proposed that eigenvectors of linear combinations of , , and   matrices as  1  +  1  +  1   produce more accurate eigenvectors.Pei and Hsue [25] showed that using the linear combination of  and  as ( + 15) produces better results than the  or  alone.Whereas, Serbes and Durak-Ata [28] showed that using ( + 30 − 7  ) produces a more accurate approximation to the continuous FRFT.Although, better commuting matrix generation method is available in the literature [31][32][33], but all these methods are explained in the context of DFT and some other transform only, excluding DFRFT.Therefore, the betterment offered by these methods can be utilized in devising the new algorithm of DFRFT based on these commuting matrixes.[34] had developed a new method for DFRFT computation.The idea of the developed method was to compute the DFRFT at any angle by a weighted summation of the DFRFT's with the special angles.In this method DFRFT computation with odd point of length can be realized by the weighted summation of the DFRFTs in special angles.The special angles were multiples of 2/ for the odd case and multiples of 2/( + 1) for even length.Regardless of even-or oddlength cases, the weighting coefficients were obtained from an IDFT operation.The -point IDFT was needed for odd length, and the ( + 1) -point was computed for even length.For discrete signal "" having odd length "", the DFRFT of "" for rotation angle "" can be computed as:

Weighted Summation Type DFRFT. Yeh and Pei
where the weighting coefficients  , was given as: The weighted summation  , has a closed-form solution: Similarly, for discrete signal "" having even length , the DFRFT of "" for rotation angle "" can be computed by the following expression: The weighting coefficients,  , can be computed as: (61) Similar to the odd length, the weighted summation  , for even case has also a closed-form solution: The eigenvector decomposition methods require  2 storages to store the transform kernel.If the computing angle is changed, the transform kernel is also changed and needs to be recomputed.In this algorithm, the DFRFT of the specified angle was stored.The DFRFT at any angle can be obtained by weighted summation of these specified DFRFTs.The weighted computation in this algorithm still takes ( 2 ) multiplications; therefore, its computation load is still ( 2 ).In [34], two implementation methods for the DFRFT computation algorithm were introduced.One was called the parallel method and the other cascade method.Similar to the parallel method, the special angle and numbers of terms of cascade forms are also different for the even and odd cases.Both of the implementation methods can have the advantages over each other for different application.The parallel method was suitable for the signal, whose DFRFTs with special angles were already known.Hence, the computation of the DFRFT will become only a linear combination of the DFRFT's in special angles.And, the cascade method means that the computation of the DFRFT can be realized from the DFRFT with only one specified angle.If the DFRFT with the specified angle can be computed efficiently, the computation of the DFRFT will become efficient.Therefore, it has been found suitable for VLSI implementation [34].
Although, the existence of one well-accepted closedform definition of DFRFT is not available; however many definitions with their pros and cons are accessible.A comparison of all the proposed classes of DFRFT based on some properties and complexity is shown in Table 1.In this Table, yes is included against the method which possesses given property/parameter, otherwise NO is mentioned.
From the above mentioned comparative analysis, it can be recapitulated that RDFRFT provides a better solution for computation of the DFRFT.The computation of RDFRFT requires more time than sampling type, linear combination type, and closed form type; however, it satisfies more properties than any other method with very less number of imposed constraints.Previously, IEEE has conferred to use the Eigenvectors decomposition type DFRFT as the most appropriate definition of DFRFT.In the year 2006, Pei et al. [25] had established the superiority of their method (MPDFRFT) over the earlier reported methods.Subsequently, in the year 2009, the MPDFRFT was modified by Pei et al. [26] and given as RDFRFT.It was also established that RDFRFT has one more degree of freedom, while satisfying the same set of properties as satisfied by MPDFRFT with equal computational complexity.Although, Serbes and Durak-Ata [28] proposed a modified version of eigenvalue decomposition technique with slight betterment but this does not redeem the advantage offered by RDFRFT, that is, security of the data.This makes the RDFRFT as a better choice in security applications.

Application of DFRFT
In this section, the application of DFRFT has been investigated in two areas.First, the detection of Doppler shifted DSB-SC signal will be performed than DFRFT has been applied in the image processing application, particularly in steganography.

Removal of Doppler Shift in DSB-SC AM Signal.
The weighted autocorrelation theorem defined for FRFT [2] is utilized to remove the unwanted frequency components introduced in double sideband suppressed-carrier amplitude modulated (DSB-SC-AM) signal.It is assumed that a Doppler component is introduced in the carrier frequency due to motion of transmitter and/or receiver.The transmitted DSB-SC amplitude-modulated signal represented by "()" with modulating signal frequency as   and carrier frequency as   is given as: Similarly, the received signal, which is corrupted by noise (Doppler frequency component-  ), represented as "()" is given by:  () = cos {2  } cos {2 (  +   ) } . ( The values of "  ,   , and   " are taken as 1 Hz, 10 Hz and as 10% of carrier frequency, respectively.The simulation is performed on the platform of MATLAB software (version 7.1) on a system having configuration: processor as Pentium-4, Intel (R) CPU-1.8GHz processor having 1 GB RAM.The effect of Doppler shift in the carrier is illustrated in Figure 2(a), where Fourier transform of original signal and Doppler corrupted signal is plotted.
It is evident from Figure 2(a) that for a received Doppler corrupted DSB-SC signal, the magnitude spectrum of the signal in Fourier domain is not giving the exact information about original frequency component present in the transmitted signal due to only one peak present in the spectrum in Fourier domain.Whereas, when the same corrupted DSB-SC signal is analyzed in fractional Fourier domain, the magnitude spectrum of FRFT of weighted auto-correlation of corrupted DSB-SC signal gives two clear peak at an angle corresponding to the value of  = 0.16, as shown in Figure 2(b).Subsequently, the actual frequency components present in transmitted DSB-SC can be estimated by observing the location of two peaks in this FRFT domain.From the observed value of FRFT domain variable " peak " at which the peak is obtained,   2. From Table 2, it can be clearly observed that the calculated value of higher and lower frequency components present in DSB-SC AM signal is matched nicely with the estimated value of higher and lower frequency components even when the signal is corrupted with the presence of a Doppler frequency component in the signal itself.

Data Hiding by Steganography in Fractional Fourier
Domain.Steganography is the technique of hiding confidential information within any other media.The objective of steganography is to hide a secret message within a covermedia in such a way that others cannot determine the presence of the hidden message.Hiding information into a media requires the following elements: (i) the cover media (C) that will hold the hidden data, (ii) the secret message (M), may be text, image, or any type of data, (iii) the stego function and its inverse, (iv) stego-key (K) may be used to hide and unhide the message.
In this exercise, the cover and secret messages are assumed as Lena image and Mandrill image, respectively.The stego function is defined as Here, first, the DFRFT (F  ) of cover image and secret image is obtained, then the steganographed image can be achieved by using the stego function defined in (66), where F −1  is representing the inverse DFRFT and  is used as stego weight.Similarly, the inverse stego function for obtaining the secret information can be written as: In this process, the values of stego weight "" and DFRFT order "" have been assumed as 0.05 and 0.2, respectively.For retrieval of secret image, the DFRFT order " 1 " in (67) should be equal to the negative of the DFRFT order "" in (66).The results obtained are shown in Figure 3.In this simulation exercise, the detection of hidden image is also performed for two other DFRFT angle parameters, both are in other side of actual, that is, at  1 = −0.4 and  1 = −0.22 in Figures 3(e) and 3(f), respectively.It can be easily verified from these results that for obtaining the hidden image correctly, the correct value of angular parameter of DFRFT should be known; hence, this correct value can also be used as a key other than stego weight.This extra key is only possible with the DFRFT and not with DFT, as there no such variation of angle is possible with DFT.This again confirms the superiority of DFRFT over DFT in the security applications.

Conclusion
In this paper, the different approaches of obtaining discrete fractional Fourier transform (DFRFT) have been analyzed.These different approaches are classified in four categories based on their methodology of evaluation.Further, with the help of a comparative analysis, it is established that the random type DFRFT emerged out as a better choice for calculating DFRFT.Thereafter, the random type DFRFT along with weighted auto-correlation theorem defined for the FRFT has been employed for the successful removal of the Doppler frequency component introduced in the DSB-SC AM signal due to time-variant motion of transmitter or receiver.Subsequently, the role of DFRFT in steganography has also been successfully performed with the availability of one extra key, in the form of DFRFT angular parameter.The results obtained are in the conformity of the superiority of the fractional Fourier transform over Fourier transform.

Figure 2 :
Figure 2: Removal of Doppler component introduced in DSB-SC signal.(a) FRFT of weighted autocorrelation of DSB-SC signal and corrupted DSB-SC signal at the FRFT order of  = 1 (i.e., Fourier transform domain) and (b) FRFT of weighted auto-correlation of corrupted DSB-SC signal at different FRFT order.

Table 1 :
Comparison of Different Types of DFRFT Algorithms.