Fractional Fourier Transform
Fractional Fourier Transform
In mathematics, in the area of harmonic analysis, the fractional Fourier transform (FRFT) is a family of
linear transformations generalizing the Fourier transform. It can be thought of as the Fourier transform to
the n-th power, where n need not be an integer — thus, it can transform a function to any intermediate
domain between time and frequency. Its applications range from filter design and signal analysis to phase
retrieval and pattern recognition.
The FRFT can be used to define fractional convolution, correlation, and other operations, and can also be
further generalized into the linear canonical transformation (LCT). An early definition of the FRFT was
introduced by Condon,[1] by solving for the Green's function for phase-space rotations, and also by
Namias,[2] generalizing work of Wiener[3] on Hermite polynomials.
However, it was not widely recognized in signal processing until it was independently reintroduced around
1993 by several groups.[4] Since then, there has been a surge of interest in extending Shannon's sampling
theorem[5][6] for signals which are band-limited in the Fractional Fourier domain.
A completely different meaning for "fractional Fourier transform" was introduced by Bailey and
Swartztrauber[7] as essentially another name for a z-transform, and in particular for the case that
corresponds to a discrete Fourier transform shifted by a fractional amount in frequency space (multiplying
the input by a linear chirp) and evaluating at a fractional set of frequency points (e.g. considering only a
small portion of the spectrum). (Such transforms can be evaluated efficiently by Bluestein's FFT algorithm.)
This terminology has fallen out of use in most of the technical literature, however, in preference to the
FRFT. The remainder of this article describes the FRFT.
Introduction
The continuous Fourier transform of a function is a unitary operator of space that maps
the function to its frequential version (all expressions are taken in the sense, rather than pointwise):
More precisely, let us introduce the parity operator that inverts , . Then the
following properties hold:
The FRFT provides a family of linear transforms that further extends this definition to handle non-integer
powers of the FT.
Definition
Note: some authors write the transform in terms of the "order a " instead of the "angle α", in which case the
α is usually a times π/2. Although these two forms are equivalent, one must be careful about which
definition the author uses.
For any real α, the α-angle fractional Fourier transform of a function ƒ is denoted by and defined
by
Formally, this formula is only valid when the input function is in a sufficiently nice space (such as L1 or
Schwartz space), and is defined via a density argument, in a way similar to that of the ordinary Fourier
transform (see article), in the general case.[8]
If α is an integer multiple of π, then the cotangent and cosecant functions above diverge. However, this can
be handled by taking the limit, and leads to a Dirac delta function in the integrand. More directly, since
must be simply f(t) or f(−t) for α an even or odd multiple of π respectively.
For α = π/2 , this becomes precisely the definition of the continuous Fourier transform, and for α = −π/2
it is the definition of the inverse continuous Fourier transform.
The FRFT argument u is neither a spatial one x nor a frequency ξ. We will see why it can be interpreted as
linear combination of both coordinates (x,ξ). When we want to distinguish the α-angular fractional
domain, we will let denote the argument of .
Remark: with the angular frequency ω convention instead of the frequency one, the FRFT formula is the
Mehler kernel,
Properties
The α-th order fractional Fourier transform operator, , has the properties:
Additivity
Linearity
Integer Orders
Inverse
Commutativity
Associativity
Unitarity
Time Reversal
Then
that is,
Transform of a scaled function
Then,
Notice that the fractional Fourier transform of cannot be expressed as a scaled version of .
Rather, the fractional Fourier transform of turns out to be a scaled and chirp modulated version of
where is a different order.
Fractional kernel
Here again the special cases are consistent with the limit behavior when α approaches a multiple of π .
symmetry:
inverse:
additivity:
Related transforms
There also exist related fractional generalizations of similar transforms such as the discrete Fourier
transform.
Generalizations
The Fourier transform is essentially bosonic; it works because it is consistent with the superposition
principle and related interference patterns. There is also a fermionic Fourier transform.[13] These have been
generalized into a supersymmetric FRFT, and a supersymmetric Radon transform.[13] There is also a
fractional Radon transform, a symplectic FRFT, and a symplectic wavelet transform.[14] Because quantum
circuits are based on unitary operations, they are useful for computing integral transforms as the latter are
unitary operators on a function space. A quantum circuit has been designed which implements the
FRFT.[15]
Interpretation
The usual interpretation of the Fourier transform is as a
transformation of a time domain signal into a frequency domain
signal. On the other hand, the interpretation of the inverse Fourier
transform is as a transformation of a frequency domain signal into a
time domain signal. Fractional Fourier transforms transform a signal
(either in the time domain or frequency domain) into the domain
between time and frequency: it is a rotation in the time–frequency A rect function turns into a sinc
domain. This perspective is generalized by the linear canonical function as the order of the fractional
transformation, which generalizes the fractional Fourier transform Fourier transform becomes 1
and allows linear transforms of the time–frequency domain other
than rotation.
Take the figure below as an example. If the signal in the time domain is rectangular (as below), it becomes
a sinc function in the frequency domain. But if one applies the fractional Fourier transform to the
rectangular signal, the transformation output will be in the domain between time and frequency.
Fractional Fourier transform
The fractional Fourier transform is a rotation operation on a time–frequency distribution. From the
definition above, for α = 0, there will be no change after applying the fractional Fourier transform, while for
α = π/2, the fractional Fourier transform becomes a plain Fourier transform, which rotates the time–
frequency distribution with π/2. For other value of α, the fractional Fourier transform rotates the time–
frequency distribution according to α. The following figure shows the results of the fractional Fourier
transform with different values of α.
Time/frequency distribution of fractional Fourier transform
Application
Fractional Fourier transform can be used in time frequency analysis and DSP.[16] It is useful to filter noise,
but with the condition that it does not overlap with the desired signal in the time–frequency domain.
Consider the following example. We cannot apply a filter directly to eliminate the noise, but with the help
of the fractional Fourier transform, we can rotate the signal (including the desired signal and noise) first. We
then apply a specific filter, which will allow only the desired signal to pass. Thus the noise will be removed
completely. Then we use the fractional Fourier transform again to rotate the signal back and we can get the
desired signal.
Fractional Fourier transform in DSP
Thus, using just truncation in the time domain, or equivalently low-pass filters in the frequency domain, one
can cut out any convex set in time–frequency space. In contrast, using time domain or frequency domain
tools without a fractional Fourier transform would only allow cutting out rectangles parallel to the axes.
Fractional Fourier transforms also have applications in quantum physics. For example, they are used to
formulate entropic uncertainty relations,[17] in high-dimensional quantum key distribution schemes with
single photons,[18] and in observing spatial entanglement of photon pairs.[19]
They are also useful in the design of optical systems and for optimizing holographic storage efficiency.[20]
See also
Least-squares spectral analysis
Fractional calculus
Mehler kernel
References
1. Condon, Edward U. (1937). "Immersion of the Fourier transform in a continuous group of
functional transformations" (https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1076889). Proc.
Natl. Acad. Sci. USA. 23 (3): 158–164. Bibcode:1937PNAS...23..158C (https://ui.adsabs.har
vard.edu/abs/1937PNAS...23..158C). doi:10.1073/pnas.23.3.158 (https://doi.org/10.1073%2
Fpnas.23.3.158). PMC 1076889 (https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1076889).
PMID 16588141 (https://pubmed.ncbi.nlm.nih.gov/16588141).
2. Namias, V. (1980). "The fractional order Fourier transform and its application to quantum
mechanics". IMA Journal of Applied Mathematics. 25 (3): 241–265.
doi:10.1093/imamat/25.3.241 (https://doi.org/10.1093%2Fimamat%2F25.3.241).
3. Wiener, N. (April 1929). "Hermitian Polynomials and Fourier Analysis". Journal of
Mathematics and Physics. 8 (1–4): 70–73. doi:10.1002/sapm19298170 (https://doi.org/10.10
02%2Fsapm19298170).
4. Almeida, Luís B. (1994). "The fractional Fourier transform and time–frequency
representations". IEEE Trans. Signal Process. 42 (11): 3084–3091.
Bibcode:1994ITSP...42.3084A (https://ui.adsabs.harvard.edu/abs/1994ITSP...42.3084A).
doi:10.1109/78.330368 (https://doi.org/10.1109%2F78.330368).
5. Tao, Ran; Deng, Bing; Zhang, Wei-Qiang; Wang, Yue (2008). "Sampling and sampling rate
conversion of band limited signals in the fractional Fourier transform domain". IEEE
Transactions on Signal Processing. 56 (1): 158–171. Bibcode:2008ITSP...56..158T (https://u
i.adsabs.harvard.edu/abs/2008ITSP...56..158T). doi:10.1109/TSP.2007.901666 (https://doi.o
rg/10.1109%2FTSP.2007.901666). S2CID 7001222 (https://api.semanticscholar.org/CorpusI
D:7001222).
6. Bhandari, A.; Marziliano, P. (2010). "Sampling and reconstruction of sparse signals in
fractional Fourier domain". IEEE Signal Processing Letters. 17 (3): 221–224.
Bibcode:2010ISPL...17..221B (https://ui.adsabs.harvard.edu/abs/2010ISPL...17..221B).
doi:10.1109/LSP.2009.2035242 (https://doi.org/10.1109%2FLSP.2009.2035242).
hdl:10356/92280 (https://hdl.handle.net/10356%2F92280). S2CID 11959415 (https://api.sem
anticscholar.org/CorpusID:11959415).
7. Bailey, D. H.; Swarztrauber, P. N. (1991). "The fractional Fourier transform and applications".
SIAM Review. 33 (3): 389–404. doi:10.1137/1033097 (https://doi.org/10.1137%2F1033097).
(Note that this article refers to the chirp-z transform variant, not the FRFT.)
8. Missbauer, Andreas (2012). Gabor Frames and the Fractional Fourier Transform (https://web.
archive.org/web/20181103131426/https://www.univie.ac.at/nuhag-php/bibtex/open_files/135
86_Missbauer_Diplom_final.pdf) (PDF) (MSc). University of Vienna. Archived from the
original (https://www.univie.ac.at/nuhag-php/bibtex/open_files/13586_Missbauer_Diplom_fi
nal.pdf) (PDF) on 3 November 2018. Retrieved 3 November 2018.
9. Candan, Kutay & Ozaktas 2000.
10. Ozaktas, Zalevsky & Kutay 2001, Chapter 6.
11. Somma, Rolando D. (2016). "Quantum simulations of one dimensional quantum systems".
Quantum Information and Computation. 16: 1125–1168. arXiv:1503.06319v2 (https://arxiv.or
g/abs/1503.06319v2).
12. Shi, Jun; Zhang, NaiTong; Liu, Xiaoping (June 2012). "A novel fractional wavelet transform
and its applications". Sci. China Inf. Sci. 55 (6): 1270–1279. doi:10.1007/s11432-011-4320-x
(https://doi.org/10.1007%2Fs11432-011-4320-x). S2CID 3772011 (https://api.semanticschol
ar.org/CorpusID:3772011).
13. De Bie, Hendrik (1 September 2008). "Fourier transform and related integral transforms in
superspace". Journal of Mathematical Analysis and Applications. 345 (1): 147–164.
arXiv:0805.1918 (https://arxiv.org/abs/0805.1918). Bibcode:2008JMAA..345..147D (https://ui.
adsabs.harvard.edu/abs/2008JMAA..345..147D). doi:10.1016/j.jmaa.2008.03.047 (https://do
i.org/10.1016%2Fj.jmaa.2008.03.047). S2CID 17066592 (https://api.semanticscholar.org/Co
rpusID:17066592).
14. Fan, Hong-yi; Hu, Li-yun (2009). "Optical transformation from chirplet to fractional Fourier
transformation kernel". Journal of Modern Optics. 56 (11): 1227–1229. arXiv:0902.1800 (http
s://arxiv.org/abs/0902.1800). Bibcode:2009JMOp...56.1227F (https://ui.adsabs.harvard.edu/a
bs/2009JMOp...56.1227F). doi:10.1080/09500340903033690 (https://doi.org/10.1080%2F09
500340903033690). S2CID 118463188 (https://api.semanticscholar.org/CorpusID:1184631
88).
15. Klappenecker, Andreas; Roetteler, Martin (January 2002). "Engineering Functional Quantum
Algorithms". Physical Review A. 67 (1): 010302. arXiv:quant-ph/0208130 (https://arxiv.org/ab
s/quant-ph/0208130). doi:10.1103/PhysRevA.67.010302 (https://doi.org/10.1103%2FPhysR
evA.67.010302). S2CID 14501861 (https://api.semanticscholar.org/CorpusID:14501861).
16. Sejdić, Ervin; Djurović, Igor; Stanković, LJubiša (June 2011). "Fractional Fourier transform
as a signal processing tool: An overview of recent developments". Signal Processing. 91 (6):
1351–1369. doi:10.1016/j.sigpro.2010.10.008 (https://doi.org/10.1016%2Fj.sigpro.2010.10.0
08). S2CID 14203403 (https://api.semanticscholar.org/CorpusID:14203403).
17. Huang, Yichen (24 May 2011). "Entropic uncertainty relations in multidimensional position
and momentum spaces". Physical Review A. 83 (5): 052124. arXiv:1101.2944 (https://arxiv.o
rg/abs/1101.2944). Bibcode:2011PhRvA..83e2124H (https://ui.adsabs.harvard.edu/abs/201
1PhRvA..83e2124H). doi:10.1103/PhysRevA.83.052124 (https://doi.org/10.1103%2FPhysR
evA.83.052124). S2CID 119243096 (https://api.semanticscholar.org/CorpusID:119243096).
18. Walborn, SP; Lemelle, DS; Tasca, DS; Souto Ribeiro, PH (13 June 2008). "Schemes for
quantum key distribution with higher-order alphabets using single-photon fractional Fourier
optics". Physical Review A. 77 (6): 062323. doi:10.1103/PhysRevA.77.062323 (https://doi.or
g/10.1103%2FPhysRevA.77.062323).
19. Tasca, DS; Walborn, SP; Souto Ribeiro, PH; Toscano, F (8 July 2008). "Detection of
transverse entanglement in phase space". Physical Review A. 78 (1): 010304(R).
arXiv:0806.3044 (https://arxiv.org/abs/0806.3044). doi:10.1103/PhysRevA.78.010304 (http
s://doi.org/10.1103%2FPhysRevA.78.010304).
20. Pégard, Nicolas C.; Fleischer, Jason W. (2011). "Optimizing holographic data storage using
a fractional Fourier transform" (http://www.opticsinfobase.org/abstract.cfm?URI=ol-36-13-255
1). Optics Letters. 36 (13): 2551–2553. Bibcode:2011OptL...36.2551P (https://ui.adsabs.harv
ard.edu/abs/2011OptL...36.2551P). doi:10.1364/OL.36.002551 (https://doi.org/10.1364%2F
OL.36.002551). PMID 21725476 (https://pubmed.ncbi.nlm.nih.gov/21725476).
Bibliography
Candan, C.; Kutay, M. A.; Ozaktas, H. M. (May 2000). "The discrete fractional Fourier
transform" (http://repository.bilkent.edu.tr/bitstream/11693/11130/1/10.1109-78.839980.pdf)
(PDF). IEEE Transactions on Signal Processing. 48 (5): 1329–1337.
Bibcode:2000ITSP...48.1329C (https://ui.adsabs.harvard.edu/abs/2000ITSP...48.1329C).
doi:10.1109/78.839980 (https://doi.org/10.1109%2F78.839980). hdl:11693/11130 (https://hd
l.handle.net/11693%2F11130).
Ding, Jian-Jiun (2007). Time frequency analysis and wavelet transform (Class notes). Taipei,
Taiwan: Department of Electrical Engineering, National Taiwan University (NTU).
Lohmann, A. W. (1993). "Image rotation, Wigner rotation and the fractional Fourier
transform". J. Opt. Soc. Am. A (10): 2181–2186. Bibcode:1993JOSAA..10.2181L (https://ui.a
dsabs.harvard.edu/abs/1993JOSAA..10.2181L). doi:10.1364/JOSAA.10.002181 (https://doi.
org/10.1364%2FJOSAA.10.002181).
Ozaktas, Haldun M.; Zalevsky, Zeev; Kutay, M. Alper (2001). The Fractional Fourier
Transform with Applications in Optics and Signal Processing (http://www.ee.bilkent.edu.tr/~h
aldun/wileybook.html). Series in Pure and Applied Optics. John Wiley & Sons. ISBN 978-0-
471-96346-2.
Pei, Soo-Chang; Ding, Jian-Jiun (2001). "Relations between fractional operations and time–
frequency distributions, and their applications". IEEE Trans. Signal Process. 49 (8): 1638–
1655. Bibcode:2001ITSP...49.1638P (https://ui.adsabs.harvard.edu/abs/2001ITSP...49.1638
P). doi:10.1109/78.934134 (https://doi.org/10.1109%2F78.934134).
Saxena, Rajiv; Singh, Kulbir (January–February 2005). "Fractional Fourier transform: A
novel tool for signal processing" (https://web.archive.org/web/20110716112239/http://journa
l.library.iisc.ernet.in/vol200501/paper2/11.pdf) (PDF). J. Indian Inst. Sci. 85: 11–26. Archived
from the original (http://journal.library.iisc.ernet.in/vol200501/paper2/11.pdf) (PDF) on 16
July 2011.
External links
DiscreteTFDs -- software for computing the fractional Fourier transform and time–frequency
distributions (https://tfd.sourceforge.net/)
"Fractional Fourier Transform (http://demonstrations.wolfram.com/FractionalFourierTransfor
m/)" by Enrique Zeleny, The Wolfram Demonstrations Project.
Dr YangQuan Chen's FRFT (Fractional Fourier Transform) Webpages (https://web.archive.o
rg/web/20090215081459/http://mechatronics.ece.usu.edu/foc/FRFT/)
LTFAT - A free (GPL) Matlab / Octave toolbox (https://ltfat.sourceforge.net/) Contains several
version of the fractional Fourier transform (https://ltfat.sourceforge.net/doc/fourier/ffracft.php).