Properties of the Fourier Transform
Total Page:16
File Type:pdf, Size:1020Kb
FOURIER BOOKLET -1 2 The Fourier Transform The definition of a one dimensional continuous function, denoted by f (x), the Fourier transform is defined by: ∞ F(u) = f (x) exp(−ı2pux)dx (1) Z−∞ with the inverse Fourier transform defined by; ∞ f (x) = F(u) exp(ı2pux)du (2) Z−∞ where it should be noted that the factors of 2p are incorporated into the transform kernel1. Some insight to the Fourier transform can be gained by considering the case of the Fourier transform of a real signal f (x). In this case the Fourier transform can be separated to give, F(u) = Fr(u) + ıFı(u) (3) where we have, ∞ Fr(u) = f (x) cos(2pux)dx Z−∞ ∞ Fı(u) = − f (x) sin(2pux)dx Z−∞ So the real part of the Fourier transform is the decomposition of f (x) in terms of cosine func- tions, and the imaginary part a decomposition in terms of sine functions. The u variable in the Fourier transform is interpreted as a frequency, for example if f (x) is a sound signal with x measured in seconds then F(u) is its frequency spectrum with u measured in Hertz (s−1). NOTE: Clearly (ux) must be dimensionless, so if x has dimensions of time then u must have dimensions of time−1. This is one of the most common applications for Fourier Transforms where f (x) is a detected signal (for example a sound made by a musical instrument), and the Fourier Transform is used to give the spectral response. 2.1 Properties of the Fourier Transform The Fourier transform has a range of useful properties, some of which are listed below. In most cases the proof of these properties is simple and can be formulated by use of equation 1 and equation 2.. The proofs of many of these properties are given in the questions and solutions at the back of this booklet. Linearity: The Fourier transform is a linear operation so that the Fourier transform of the sum of two functions is given by the sum of the individual Fourier transforms. Therefore, F fa f (x) + bg(x)g = aF(u) + bG(u) (4) 1There are various definitions of the Fourier transform that puts the 2p either inside the kernel or as external scaling factors. The difference between them whether the variable in Fourier space is a “frequency” or “angular frequency”. The difference between the definitions are clearly just a scaling factor. The optics and digital Fourier applications the 2p is usually defined to be inside the kernel but in solid state physics and differential equation solution the 2p constant is usually an external scaling factor. School of Physics Fourier Transform Revised: 10 September 2007 FOURIER BOOKLET -2 where F(u) and G(u) are the Fourier transforms of f (x) and and g(x) and a and b are constants. This property is central to the use of Fourier transforms when describing linear systems. Complex Conjugate: The Fourier transform of the Complex Conjugate of a function is given by F f f ∗(x)g = F∗(−u) (5) where F(u) is the Fourier transform of f (x). Forward and Inverse: We have that F fF(u)g = f (−x) (6) so that if we apply the Fourier transform twice to a function, we get a spatially reversed version of the function. Similarly with the inverse Fourier transform we have that, F −1 f f (x)g = F(−u) (7) so that the Fourier and inverse Fourier transforms differ only by a sign. Differentials: The Fourier transform of the derivative of a functions is given by d f (x) = ı2puF(u) (8) F dx and the second derivative is given by d2 f (x) = −(2pu)2 F(u) (9) F dx2 This property will be used in the DIGITAL IMAGE ANALYSIS and THEORY OF IMAGE PRO- CESSING course to form the derivative of an image. Power Spectrum: The Power Spectrum of a signal is defined by the modulus square of the Fourier transform, being jF(u)j2. This can be interpreted as the power of the frequency com- ponents. Any function and its Fourier transform obey the condition that ∞ ∞ j f (x)j2 dx = jF(u)j2 du (10) Z−∞ Z−∞ which is frequently known as Parseval's Theorem2. If f (x) is interpreted at a voltage, then this theorem states that the power is the same whether measured in real (time), or Fourier (frequency) space. 2.2 Two Dimensional Fourier Transform Since the three courses covered by this booklet use two-dimensional scalar potentials or images we will be dealing with two dimensional function. We will define the two dimensional Fourier transform of a continuous function f (x;y) by, F(u;v) = f (x;y) exp(−ı2p(ux + vy)) dxdy (11) Z Z 2Strictly speaking Parseval's Theorem applies to the case of Fourier series, and the equivalent theorem for Fourier transforms is correctly, but less commonly, known as Rayleigh's theorem School of Physics Fourier Transform Revised: 10 September 2007 FOURIER BOOKLET -3 with the inverse Fourier transform defined by; f (x;y) = F(u;v) exp(ı2p(ux + vy)) dudv (12) Z Z where the limits of integration are taken from −¥ ! ¥3 Again for a real two dimensional function f (x;y), the Fourier transform can be considered as the decomposition of a function into its sinusoidal components. If f (x;y) is considered to be an image with the “brightness” of the image at point (x0;y0) given by f (x0;y0), then variables x;y have the dimensions of length. In Fourier space the variables u;v have therefore the dimensions of inverse length, which is interpreted as Spatial Frequency. NOTE: Typically x and y are measured in mm so that u and v have are in units of mm−1 also referred to at lines per mm. The Fourier transform can then be taken as being the decomposition of the image into two di- mensional sinusoidal spatial frequency components. This property will be examined in greater detail the relevant courses. The properties of one the dimensional Fourier transforms covered in the previous section con- vert into two dimensions. Clearly the derivatives then become ¶ f (x;y) F = ı2puF(u;v) (13) ¶x and with ¶ f (x;y) F = ı2pvF(u;v) (14) ¶y yielding the important result that, F ∇2 f (x;y) = −(2pw)2 F(u;v) (15) where we have that w2 = u2 + v2. So that taking the Laplacian of a function in real space is equivalent to multiplying its Fourier transform by a circularly symmetric quadratic of −4p2w2. The two dimensional Fourier Transform F(u;v), of a function f (x;y) is a separable operation, and can be written as, F(u;v) = P(u;y)exp(−ı2pvy)dy (16) Z where P(u;y) = f (x;y) exp(−ı2pux)dx (17) Z where P(u;y) is the Fourier Transform of f (x;y) with respect to x only. This property of separability will be considered in greater depth with regards to digital images and will lead to an implementation of two dimensional discrete Fourier Transforms in terms of one dimensional Fourier Transforms. 2.3 The Three-Dimensional Fourier Transform In the three dimensional case we have a function f (~r) where ~r = (x;y;z), then the three- dimensional Fourier Transform F(~s) = f (~r) exp(−ı2p~r :~s) d~r Z Z Z 3Unless otherwise specified all integral limits will be assumed to be from −¥ ! ¥ School of Physics Fourier Transform Revised: 10 September 2007 FOURIER BOOKLET -4 where ~s = (u;v;w) being the three reciprocal variables each with units length−1. Similarly the inverse Fourier Transform is given by f (~r) = F(~s) exp(ı2p~r:~s) d~s Z Z Z This is used extensively in solid state physics where the three-dimensional Fourier Transform of a crystal structures is usually called Reciprocal Space4. The three-dimensional Fourier Transform is again separable into one-dimensional Fourier Trans- form. This property is independent of the dimensionality and multi-dimensional Fourier Trans- form can be formulated as a series of one dimensional Fourier Transforms. 4This is also referred to as~k-space where~k = 2p~s School of Physics Fourier Transform Revised: 10 September 2007.