Fractional Fourier transform

Fractional Fourier transform

In mathematics, in the area of harmonic analysis, the fractional Fourier transform (FRFT) is a linear transformation 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 an "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 given by Namias (1980), but it was not widely recognized until it was independently reinvented around 1993 by several groups of researchers (Almeida, 1994).

A completely different meaning for "fractional Fourier transform" was introduced by Bailey and Swartztrauber (1991) 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.

See also the chirplet transform for a related generalization of the Fourier transform.

Definition

If the continuous Fourier transform of a function f(t) is denoted by mathcal{F}(f), then mathcal{F}^2(f)=mathcal{F}(mathcal{F}(f)), and in general mathcal{F}^{(n+1)}(f)=mathcal{F}(mathcal{F}^n(f)); similarly, mathcal{F}^{-n}(F) denotes the "n"-th power of the inverse transform mathcal{F}^{-1}(F) of F(omega). The FRFT further extends this definition to handle non-integer powers n=2alpha/pi for any real alpha, denoted by mathcal{F}_alpha(f) and having the properties:

:mathcal{F}_alpha(f) = mathcal{F}^{2alpha/pi}(f)

when n=2alpha/pi is an integer, and:

:mathcal{F}_{alpha+eta}(f) = mathcal{F}_alpha(mathcal{F}_eta(f)) = mathcal{F}_eta(mathcal{F}_alpha(f)).

More specifically, mathcal{F}_alpha(f) is given by the equation:

:mathcal{F}_alpha(f)(omega) = sqrt{frac{1-icot(alpha)}{2pi e^{i cot(alpha) omega^2/2} int_{-infty}^infty e^{-icsc(alpha) omega t + i cot(alpha) t^2/2}f(t), dt.

Note that, for alpha=pi/2, this becomes precisely the definition of the continuous Fourier transform, and for alpha=-pi/2 it is the definition of the inverse continuous Fourier transform.

If alpha is an integer multiple of pi, 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 easily, since mathcal{F}^2(f)=f(-t), mathcal{F}_alpha(f) must be simply f(t) or f(-t) for alpha an even or odd multiple of pi, respectively.

There also exist related fractional generalizations of similar transforms such as the discrete Fourier transform.

Physical Meaning of the Fractional Fourier Transform

The physical meaning of the Fourier transform is as a transformation of a time domain signal into a frequency domain signal. On the other hand, the physical meaning of the inverse Fourier transform is as a transformation of a frequency domain signal into a time domain signal. Apparently, fractional Fourier transforms can transform a signal (either in the time domain or frequency domain) into the domain between time and frequency.

Take the below figure as an example. If the signal in the time domain is rectangular (as below), it will become a sinc function in the frequency domain. But if we apply the fractional Fourier transform to the rectangular signal, the transformation output will be in the domain between time and frequency.


Actually, fractional Fourier transform is a rotation operation on the time frequency distribution. From the definition above, for α=0, there will be no change after applying fractional Fourier transform, and for α=π/2, fractional Fourier transform becomes a Fourier transform, which rotates the time frequency distribution with π/2. For other value of α, 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 α.

Application

Fractional Fourier transform can be used in time frequency analysis and DSP. It is useful to filter noise, but with the condition that it does not overlap with the desired signal in the time frequency domain. Let’s see 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.

ee also

Other time-frequency transforms:

* short-time Fourier transform
* wavelet transform
* chirplet transform

External links

* [http://tfd.sourceforge.net/ DiscreteTFDs -- software for computing the fractional Fourier transform and time-frequency distributions]
*" [http://demonstrations.wolfram.com/FractionalFourierTransform/ Fractional Fourier Transform] " by Enrique Zeleny, The Wolfram Demonstrations Project.

References

* V. Namias, "The fractional order Fourier transform and its application to quantum mechanics," "J. Inst. Appl. Math." 25, 241–265 (1980).
* A. W. Lohmann, "Image rotation, Wigner rotation and the fractional Fourier transform," "J. Opt. Soc. Am." A 10, 2181–2186 (1993).
* Luís B. Almeida, "The fractional Fourier transform and time-frequency representations," "IEEE Trans. Sig. Processing" 42 (11), 3084–3091 (1994).
* Haldun M. Ozaktas, Zeev Zalevsky and M. Alper Kutay. "The Fractional Fourier Transform with Applications in Optics and Signal Processing". John Wiley & Sons (2001). Series in Pure and Applied Optics.
* Soo-Chang Pei and Jian-Jiun Ding, "Relations between fractional operations and time-frequency distributions, and their applications," "IEEE Trans. Sig. Processing" 49 (8), 1638–1655 (2001).
* D. H. Bailey and P. N. Swarztrauber, "The fractional Fourier transform and applications," "SIAM Review" 33, 389-404 (1991). (Note that this article refers to the chirp-z transform variant, not the FRFT.)
*Jian-Jiun Ding, Time frequency analysis and wavelet transform class notes, the Department of Electrical Engineering, National Taiwan University (NTU), Taipei, Taiwan, 2007.


Wikimedia Foundation. 2010.

Игры ⚽ Поможем решить контрольную работу

Look at other dictionaries:

  • Fourier transform — Fourier transforms Continuous Fourier transform Fourier series Discrete Fourier transform Discrete time Fourier transform Related transforms The Fourier transform is a mathematical operation that decomposes a function into its constituent… …   Wikipedia

  • Discrete Fourier transform — Fourier transforms Continuous Fourier transform Fourier series Discrete Fourier transform Discrete time Fourier transform Related transforms In mathematics, the discrete Fourier transform (DFT) is a specific kind of discrete transform, used in… …   Wikipedia

  • Short-time Fourier transform — The short time Fourier transform (STFT), or alternatively short term Fourier transform, is a Fourier related transform used to determine the sinusoidal frequency and phase content of local sections of a signal as it changes over time. STFT… …   Wikipedia

  • Fourier — (pronEng|ˈfʊərieɪ, French pronunciation IPA2|fuʁie) may refer to:*Charles Fourier (1772–1837), a French utopian socialist thinker *Joseph Fourier (1768–1830), a French mathematician and physicist **Mathematics, physics, and engineering terms… …   Wikipedia

  • Discrete-time Fourier transform — In mathematics, the discrete time Fourier transform (DTFT) is one of the specific forms of Fourier analysis. As such, it transforms one function into another, which is called the frequency domain representation, or simply the DTFT , of the… …   Wikipedia

  • Fourier analysis — In mathematics, Fourier analysis is a subject area which grew out of the study of Fourier series. The subject began with trying to understand when it was possible to represent general functions by sums of simpler trigonometric functions. The… …   Wikipedia

  • Fourier operator — The Fourier operator is the kernel of the Fredholm integral of the first kind that defines the continuous Fourier transform.It may be thought of as a limiting case for when the size of the discrete Fourier transform increases without bound while… …   Wikipedia

  • Chirplet transform — Comparison of wave, wavelet, chirp, and chirplet In signal processing, the chirplet transform is an inner product of an input signal with a family of analysis primitives called chirplets. Contents …   Wikipedia

  • List of Fourier-related transforms — This is a list of linear transformations of functions related to Fourier analysis. Such transformations map a function to a set of coefficients of basis functions, where the basis functions are sinusoidal and are therefore strongly localized in… …   Wikipedia

  • Fractional calculus — is a branch of mathematical analysis that studies the possibility of taking real number powers of the differential operator ::D = frac{d}{dx} , and the integration operator J . (Usually J is used in favor of I to avoid confusion with other I like …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”