User:WakingLili/Sandbox/FT for Ava

The Fourier transform, named after Joseph Fourier, is a mathematical transform with many applications in physics and engineering. Very commonly it transforms a mathematical function of time, f(t), into a new function, sometimes denoted by or F, whose argument is frequency with units of cycles/s (hertz) or radians per second. The new function is then known as the Fourier transform and/or the frequency spectrum of the function f. The Fourier transform is also a reversible operation. Thus, given the function one can determine the original function, f. (See Fourier inversion theorem.) f and are also respectively known as time domain and frequency domain representations of the same "event". Most often perhaps, f is a real-valued function, and is complex valued, where a complex number describes both the amplitude and phase of a corresponding frequency component. In general, f is also complex, such as the analytic representation of a real-valued function. The term "Fourier transform" refers to both the transform operation and to the complex-valued function it produces.

In the case of a periodic function (for example, a continuous but not necessarily sinusoidal musical sound), the Fourier transform can be simplified to the calculation of a discrete set of complex amplitudes, called Fourier series coefficients. Also, when a time-domain function is sampled to facilitate storage or computer-processing, it is still possible to recreate a version of the original Fourier transform according to the Poisson summation formula, also known as discrete-time Fourier transform. These topics are addressed in separate articles. For an overview of those and other related operations, refer to Fourier analysis or List of Fourier-related transforms.

Definition

edit

There are several common conventions for defining the Fourier transform ƒ̂ of an integrable function f : RC (Kaiser 1994, p. 29), (Rahman 2011, p. 11). This article will use the definition:

 ,   for every real number ξ.

When the independent variable x represents time (with SI unit of seconds), the transform variable ξ represents frequency (in hertz). Under suitable conditions, f is determined by ƒ̂ via the inverse transform:

    for every real number x.

The statement that f can be reconstructed from ƒ̂ is known as the Fourier inversion theorem, and was first introduced in Fourier's Analytical Theory of Heat (Fourier 1822, p. 525), (Fourier & Freeman 1878, p. 408), although what would be considered a proof by modern standards was not given until much later (Titchmarsh 1948, p. 1). The functions f and ƒ̂ often are referred to as a Fourier integral pair or Fourier transform pair (Rahman 2011, p. 10).

For other common conventions and notations, including using the angular frequency ω instead of the frequency ξ, see Other conventions and Other notations below. The Fourier transform on Euclidean space is treated separately, in which the variable x often represents position and ξ momentum.

Introduction

edit
 
The Fourier transform relates the function's time domain, shown in red, to the function's frequency domain, shown in blue. The component frequencies, spread across the frequency spectrum, are represented as peaks in the frequency domain.

The motivation for the Fourier transform comes from the study of Fourier series. In the study of Fourier series, complicated but periodic functions are written as the sum of simple waves mathematically represented by sines and cosines. The Fourier transform is an extension of the Fourier series that results when the period of the represented function is lengthened and allowed to approach infinity.(Taneja 2008, p. 192)

Due to the properties of sine and cosine, it is possible to recover the amplitude of each wave in a Fourier series using an integral. In many cases it is desirable to use Euler's formula, which states that eiθ= cos(2πθ) + isin(2πθ), to write Fourier series in terms of the basic waves e2πiθ. This has the advantage of simplifying many of the formulas involved, and provides a formulation for Fourier series that more closely resembles the definition followed in this article. Re-writing sines and cosines as complex exponentials makes it necessary for the Fourier coefficients to be complex valued. The usual interpretation of this complex number is that it gives both the amplitude (or size) of the wave present in the function and the phase (or the initial angle) of the wave. These complex exponentials sometimes contain negative "frequencies". If θ is measured in seconds, then the waves e2πiθ and e−2πiθ both complete one cycle per second, but they represent different frequencies in the Fourier transform. Hence, frequency no longer measures the number of cycles per unit time, but is still closely related.

There is a close connection between the definition of Fourier series and the Fourier transform for functions f which are zero outside of an interval. For such a function, we can calculate its Fourier series on any interval that includes the points where f is not identically zero. The Fourier transform is also defined for such a function. As we increase the length of the interval on which we calculate the Fourier series, then the Fourier series coefficients begin to look like the Fourier transform and the sum of the Fourier series of f begins to look like the inverse Fourier transform. To explain this more precisely, suppose that T is large enough so that the interval [−T/2,T/2] contains the interval on which f is not identically zero. Then the n-th series coefficient cn is given by:

 

Comparing this to the definition of the Fourier transform, it follows that cn = ƒ̂(n/T) since f(x) is zero outside [−T/2,T/2]. Thus the Fourier coefficients are just the values of the Fourier transform sampled on a grid of width 1/T. As T increases the Fourier coefficients more closely represent the Fourier transform of the function.

Under appropriate conditions, the sum of the Fourier series of f will equal the function f. In other words, f can be written:

 

where the last sum is simply the first sum rewritten using the definitions ξn = n/T, and Δξ = (n + 1)/Tn/T = 1/T.

This second sum is a Riemann sum, and so by letting T → ∞ it will converge to the integral for the inverse Fourier transform given in the definition section. Under suitable conditions this argument may be made precise (Stein & Shakarchi 2003).

In the study of Fourier series the numbers cn could be thought of as the "amount" of the wave present in the Fourier series of f. Similarly, as seen above, the Fourier transform can be thought of as a function that measures how much of each individual frequency is present in our function f, and we can recombine these waves by using an integral (or "continuous sum") to reproduce the original function.

Example

edit

The following images provide a visual illustration of how the Fourier transform measures whether a frequency is present in a particular function. The function depicted f(t) = cos(6πt) et2 oscillates at 3 hertz (if t measures seconds) and tends quickly to 0. (The second factor in this equation is an envelope function that shapes the continuous sinusoid into a short pulse. Its general form is a Gaussian function). This function was specially chosen to have a real Fourier transform which can easily be plotted. The first image contains its graph. In order to calculate ƒ̂(3) we must integrate e−2πi(3t)f(t). The second image shows the plot of the real and imaginary parts of this function. The real part of the integrand is almost always positive, because when f(t) is negative, the real part of e−2πi(3t) is negative as well. Because they oscillate at the same rate, when f(t) is positive, so is the real part of e−2πi(3t). The result is that when you integrate the real part of the integrand you get a relatively large number (in this case 0.5). On the other hand, when you try to measure a frequency that is not present, as in the case when we look at ƒ̂(5), the integrand oscillates enough so that the integral is very small. The general situation may be a bit more complicated than this, but this in spirit is how the Fourier transform measures how much of an individual frequency is present in a function f(t).

Properties of the Fourier transform

edit

Here we assume f(x), g(x) and h(x) are integrable functions, are Lebesgue-measurable on the real line, and satisfy:

 

We denote the Fourier transforms of these functions by   ,    and    respectively.

Basic properties

edit

The Fourier transform has the following basic properties: (Pinsky 2002).

Linearity
For any complex numbers a and b, if h(x) = af(x) + bg(x), then   
Translation
For any real number x0, if     then   
Modulation
For any real number ξ0 if   then   
Scaling
For a non-zero real number a, if h(x) = f(ax), then        The case a = −1 leads to the time-reversal property, which states: if h(x) = f(−x), then  
Conjugation
If     then   
In particular, if f is real, then one has the reality condition , that is,   is a hermitian function.
And if f is purely imaginary, then   
Integration
Substituting   in the definition, we obtain
 

That is, the evaluation of the Fourier transform in the origin ( ) equals the integral of f all over its domain.

Convolution theorem

edit

The Fourier transform translates between convolution and multiplication of functions. If f(x) and g(x) are integrable functions with Fourier transforms   and   respectively, then the Fourier transform of the convolution is given by the product of the Fourier transforms   and   (under other conventions for the definition of the Fourier transform a constant factor may appear).

This means that if:

 

where ∗ denotes the convolution operation, then:

 

In linear time invariant (LTI) system theory, it is common to interpret g(x) as the impulse response of an LTI system with input f(x) and output h(x), since substituting the unit impulse for f(x) yields h(x) = g(x). In this case,   represents the frequency response of the system.

Conversely, if f(x) can be decomposed as the product of two square integrable functions p(x) and q(x), then the Fourier transform of f(x) is given by the convolution of the respective Fourier transforms   and  .

Cross-correlation theorem

edit

In an analogous manner, it can be shown that if h(x) is the cross-correlation of f(x) and g(x):

 

then the Fourier transform of h(x) is:

 

As a special case, the autocorrelation of function f(x) is:

 

for which

 

Other notations

edit

Other common notations for ƒ̂(ξ) include:

 

Denoting the Fourier transform by a capital letter corresponding to the letter of function being transformed (such as f(x) and F(ξ)) is especially common in the sciences and engineering. In electronics, the omega (ω) is often used instead of ξ due to its interpretation as angular frequency, sometimes it is written as F(jω), where j is the imaginary unit, to indicate its relationship with the Laplace transform, and sometimes it is written informally as F(2πf) in order to use ordinary frequency.

The interpretation of the complex function ƒ̂(ξ) may be aided by expressing it in polar coordinate form

 

in terms of the two real functions A(ξ) and φ(ξ) where:

 

is the amplitude and

   

is the phase (see arg function).

Then the inverse transform can be written:

 

which is a recombination of all the frequency components of f(x). Each component is a complex sinusoid of the form e2πixξ  whose amplitude is A(ξ) and whose initial phase angle (at x = 0) is φ(ξ).

The Fourier transform may be thought of as a mapping on function spaces. This mapping is here denoted   and   is used to denote the Fourier transform of the function f. This mapping is linear, which means that   can also be seen as a linear transformation on the function space and implies that the standard notation in linear algebra of applying a linear transformation to a vector (here the function f) can be used to write   instead of  . Since the result of applying the Fourier transform is again a function, we can be interested in the value of this function evaluated at the value ξ for its variable, and this is denoted either as   or as  . Notice that in the former case, it is implicitly understood that   is applied first to f and then the resulting function is evaluated at ξ, not the other way around.

In mathematics and various applied sciences it is often necessary to distinguish between a function f and the value of f when its variable equals x, denoted f(x). This means that a notation like   formally can be interpreted as the Fourier transform of the values of f at x. Despite this flaw, the previous notation appears frequently, often when a particular function or a function of a particular variable is to be transformed.

For example,   is sometimes used to express that the Fourier transform of a rectangular function is a sinc function,

or   is used to express the shift property of the Fourier transform.

Notice, that the last example is only correct under the assumption that the transformed function is a function of x, not of x0.

Other conventions

edit

The Fourier transform can also be written in terms of angular frequency: ω = 2πξ whose units are radians per second.

The substitution ξ = ω/(2π) into the formulas above produces this convention:

 

Under this convention, the inverse transform becomes:

 

Unlike the convention followed in this article, when the Fourier transform is defined this way, it is no longer a unitary transformation on L2(Rn). There is also less symmetry between the formulas for the Fourier transform and its inverse.

Another convention is to split the factor of (2π)n evenly between the Fourier transform and its inverse, which leads to definitions:

 
 

Under this convention, the Fourier transform is again a unitary transformation on L2(Rn). It also restores the symmetry between the Fourier transform and its inverse.

Variations of all three conventions can be created by conjugating the complex-exponential kernel of both the forward and the reverse transform. The signs must be opposites. Other than that, the choice is (again) a matter of convention.

Summary of popular forms of the Fourier transform
ordinary frequency ξ (hertz) unitary  

 

angular frequency ω (rad/s) non-unitary  

 

unitary  

 

As discussed above, the characteristic function of a random variable is the same as the Fourier–Stieltjes transform of its distribution measure, but in this context it is typical to take a different convention for the constants. Typically characteristic function is defined  .

As in the case of the "non-unitary angular frequency" convention above, there is no factor of 2π appearing in either of the integral, or in the exponential. Unlike any of the conventions appearing above, this convention takes the opposite sign in the exponential.

Tables of important Fourier transforms

edit

The following tables record some closed form Fourier transforms. For functions f(x), g(x) and h(x) denote their Fourier transforms by ƒ̂,  , and   respectively. Only the three most common conventions are included. It may be useful to notice that entry 105 gives a relationship between the Fourier transform of a function and the original function, which can be seen as relating the Fourier transform and its inverse.

Functional relationships

edit

The Fourier transforms in this table may be found in Erdélyi (1954) or Kammler (2000, appendix).

Function Fourier transform
unitary, ordinary frequency
Fourier transform
unitary, angular frequency
Fourier transform
non-unitary, angular frequency
Remarks
   

 

 

 

 

 

Definition
101         Linearity
102         Shift in time domain
103         Shift in frequency domain, dual of 102
104         Scaling in the time domain. If   is large, then   is concentrated around 0 and   spreads out and flattens.
105         Duality. Here   needs to be calculated using the same method as Fourier transform column. Results from swapping "dummy" variables of x and ξ or ω or μ.
106        
107         This is the dual of 106
108         The notation   denotes the convolution of f and g — this rule is the convolution theorem
109         This is the dual of 108
110 For   purely real       Hermitian symmetry.   indicates the complex conjugate.
111 For   a purely real even function  ,   and   are purely real even functions.
112 For   a purely real odd function  ,   and   are purely imaginary odd functions.

Square-integrable functions

edit

The Fourier transforms in this table may be found in (Campbell & Foster 1948), (Erdélyi 1954), or the appendix of (Kammler 2000).

Function Fourier transform
unitary, ordinary frequency
Fourier transform
unitary, angular frequency
Fourier transform
non-unitary, angular frequency
Remarks
   

 

 

 

 

 

201         The rectangular pulse and the normalized sinc function, here defined as sinc(x) = sin(πx)/(πx)
202         Dual of rule 201. The rectangular function is an ideal low-pass filter, and the sinc function is the non-causal impulse response of such a filter.
203         The function tri(x) is the triangular function
204         Dual of rule 203.
205         The function u(x) is the Heaviside unit step function and a>0.
206         This shows that, for the unitary Fourier transforms, the Gaussian function exp(−αx2) is its own Fourier transform for some choice of α. For this to be integrable we must have Re(α)>0.
207         For a>0. That is, the Fourier transform of a decaying exponential function is a Lorentzian function.
208         Hyperbolic secant is its own Fourier transform
209    

   

 

   

 

   

  is the Hermite's polynomial. If a = 1 then the Gauss-Hermite functions are eigenfunctions of the Fourier transform operator. For a derivation, see Hermite polynomial. The formula reduces to 206 for n = 0.

Distributions

edit

The Fourier transforms in this table may be found in (Erdélyi 1954) or the appendix of (Kammler 2000).

Function Fourier transform
unitary, ordinary frequency
Fourier transform
unitary, angular frequency
Fourier transform
non-unitary, angular frequency
Remarks
   

 

 

 

 

 

301         The distribution δ(ξ) denotes the Dirac delta function.
302         Dual of rule 301.
303         This follows from 103 and 301.
304         This follows from rules 101 and 303 using Euler's formula:  

 

305         This follows from 101 and 303 using  

 

306        
307        
308         Here, n is a natural number and   is the n-th distribution derivative of the Dirac delta function. This rule follows from rules 107 and 301. Combining this rule with 101, we can transform all polynomials.
309         Here sgn(ξ) is the sign function. Note that 1/x is not a distribution. It is necessary to use the Cauchy principal value when testing against Schwartz functions. This rule is useful in studying the Hilbert transform.
310  

 

      1/xn is the homogeneous distribution defined by the distributional derivative  
311         This formula is valid for 0 > α > −1. For α > 0 some singular terms arise at the origin that can be found by differentiating 318. If Re α > −1, then   is a locally integrable function, and so a tempered distribution. The function   is a holomorphic function from the right half-plane to the space of tempered distributions. It admits a unique meromorphic extension to a tempered distribution, also denoted   for α ≠ −2, −4, ... (See homogeneous distribution.)
312         The dual of rule 309. This time the Fourier transforms need to be considered as Cauchy principal value.
313         The function u(x) is the Heaviside unit step function; this follows from rules 101, 301, and 312.
314         This function is known as the Dirac comb function. This result can be derived from 302 and 102, together with the fact that  

  as distributions.

315         The function J0(x) is the zeroth order Bessel function of first kind.
316         This is a generalization of 315. The function Jn(x) is the n-th order Bessel function of first kind. The function Tn(x) is the Chebyshev polynomial of the first kind.
317           is the Euler–Mascheroni constant.
318         This formula is valid for 1 > α > 0. Use differentiation to derive formula for higher exponents. u is the Heaviside function.

Two-dimensional functions

edit
Function Fourier transform
unitary, ordinary frequency
Fourier transform
unitary, angular frequency
Fourier transform
non-unitary, angular frequency
400    

 

 

 

 

 

401        
402        
Remarks

To 400: The variables ξx, ξy, ωx, ωy, νx and νy are real numbers. The integrals are taken over the entire plane.

To 401: Both functions are Gaussians, which may not have unit volume.

To 402: The function is defined by circ(r)=1 0≤r≤1, and is 0 otherwise. This is the Airy distribution, and is expressed using J1 (the order 1 Bessel function of the first kind). (Stein & Weiss 1971, Thm. IV.3.3)

Formulas for general n-dimensional functions

edit
Function Fourier transform
unitary, ordinary frequency
Fourier transform
unitary, angular frequency
Fourier transform
non-unitary, angular frequency
500    

 

 

 

 

 

501    
 
 
 
 
 
502    
503    
Remarks

To 501: The function χ[0,1] is the indicator function of the interval [0, 1]. The function Γ(x) is the gamma function. The function Jn/2 + δ is a Bessel function of the first kind, with order n/2 + δ. Taking n = 2 and δ = 0 produces 402. (Stein & Weiss 1971, Thm. 4.15)

To 502: See Riesz potential. The formula also holds for all α ≠ −n, −n − 1, ... by analytic continuation, but then the function and its Fourier transforms need to be understood as suitably regularized tempered distributions. See homogeneous distribution.

To 503: This is the formula for a multivariate normal distribution normalized to 1 with a mean of 0. Bold variables are vectors or matrices. Following the notation of the aforementioned page,   and  

See also

edit

References

edit
  • Boashash, B., ed. (2003), Time-Frequency Signal Analysis and Processing: A Comprehensive Reference, Oxford: Elsevier Science, ISBN 0-08-044335-4
  • Bochner S., Chandrasekharan K. (1949), Fourier Transforms, Princeton University Press
  • Bracewell, R. N. (2000), The Fourier Transform and Its Applications (3rd ed.), Boston: McGraw-Hill, ISBN 0-07-116043-4.
  • Campbell, George; Foster, Ronald (1948), Fourier Integrals for Practical Applications, New York: D. Van Nostrand Company, Inc..
  • Condon, E. U. (1937), "Immersion of the Fourier transform in a continuous group of functional transformations", Proc. Nat. Acad. Sci. USA, 23: 158–164.
  • Duoandikoetxea, Javier (2001), Fourier Analysis, American Mathematical Society, ISBN 0-8218-2172-5.
  • Dym, H; McKean, H (1985), Fourier Series and Integrals, Academic Press, ISBN 978-0-12-226451-1.
  • Erdélyi, Arthur, ed. (1954), Tables of Integral Transforms, vol. 1, New Your: McGraw-Hill
  • Fourier, J. B. Joseph (1822), Théorie Analytique de la Chaleur, Paris: Chez Firmin Didot, père et fils
  • Fourier, J. B. Joseph; Freeman, Alexander, translator (1878), The Analytical Theory of Heat, The University Press {{citation}}: |first2= has generic name (help)CS1 maint: multiple names: authors list (link)
  • Grafakos, Loukas (2004), Classical and Modern Fourier Analysis, Prentice-Hall, ISBN 0-13-035399-X.
  • Hewitt, Edwin; Ross, Kenneth A. (1970), Abstract harmonic analysis. Vol. II: Structure and analysis for compact groups. Analysis on locally compact Abelian groups, Die Grundlehren der mathematischen Wissenschaften, Band 152, Berlin, New York: Springer-Verlag, MR 0262773.
  • Hörmander, L. (1976), Linear Partial Differential Operators, Volume 1, Springer-Verlag, ISBN 978-3-540-00662-6.
  • James, J.F. (2011), A Student's Guide to Fourier Transforms (3rd ed.), New York: Cambridge University Press, ISBN 978-0-521-17683-5.
  • Kaiser, Gerald (1994), A Friendly Guide to Wavelets, Birkhäuser, ISBN 0-8176-3711-7
  • Kammler, David (2000), A First Course in Fourier Analysis, Prentice Hall, ISBN 0-13-578782-3
  • Katznelson, Yitzhak (1976), An introduction to Harmonic Analysis, Dover, ISBN 0-486-63331-4
  • Knapp, Anthony W. (2001), Representation Theory of Semisimple Groups: An Overview Based on Examples, Princeton University Press, ISBN 978-0-691-09089-4
  • Pinsky, Mark (2002), Introduction to Fourier Analysis and Wavelets, Brooks/Cole, ISBN 0-534-37660-6
  • Polyanin, A. D.; Manzhirov, A. V. (1998), Handbook of Integral Equations, Boca Raton: CRC Press, ISBN 0-8493-2876-4.
  • Rudin, Walter (1987), Real and Complex Analysis (Third ed.), Singapore: McGraw Hill, ISBN 0-07-100276-6.
  • Rahman, Matiur (2011), Applications of Fourier Transforms to Generalized Functions, WIT Press, ISBN 1845645642.
  • Stein, Elias; Shakarchi, Rami (2003), Fourier Analysis: An introduction, Princeton University Press, ISBN 0-691-11384-X.
  • Stein, Elias; Weiss, Guido (1971), Introduction to Fourier Analysis on Euclidean Spaces, Princeton, N.J.: Princeton University Press, ISBN 978-0-691-08078-9.
  • Taneja, HC (2008), "Chapter 18: Fourier integrals and Fourier transforms", Advanced Engineering Mathematics:, Volume 2, New Delhi, India: I. K. International Pvt Ltd, ISBN 8189866567.
  • Titchmarsh, E (1948), Introduction to the theory of Fourier integrals (2nd ed.), Oxford University: Clarendon Press (published 1986), ISBN 978-0-8284-0324-5.
  • Wilson, R. G. (1995), Fourier Series and Optical Transform Techniques in Contemporary Optics, New York: Wiley, ISBN 0-471-30357-7.
  • Yosida, K. (1968), Functional Analysis, Springer-Verlag, ISBN 3-540-58654-7.
edit