Fourier Analysis

Total Page:16

File Type:pdf, Size:1020Kb

Fourier Analysis FOURIER ANALYSIS Lucas Illing 2008 Contents 1 Fourier Series 2 1.1 General Introduction . 2 1.2 Discontinuous Functions . 5 1.3 Complex Fourier Series . 7 2 Fourier Transform 8 2.1 Definition . 8 2.2 The issue of convention . 11 2.3 Convolution Theorem . 12 2.4 Spectral Leakage . 13 3 Discrete Time 17 3.1 Discrete Time Fourier Transform . 17 3.2 Discrete Fourier Transform (and FFT) . 19 4 Executive Summary 20 1 1. Fourier Series 1 Fourier Series 1.1 General Introduction Consider a function f(τ) that is periodic with period T . f(τ + T ) = f(τ) (1) We may always rescale τ to make the function 2π periodic. To do so, define 2π a new independent variable t = T τ, so that f(t + 2π) = f(t) (2) So let us consider the set of all sufficiently nice functions f(t) of a real variable t that are periodic, with period 2π. Since the function is periodic we only need to consider its behavior on one interval of length 2π, e.g. on the interval (−π; π). The idea is to decompose any such function f(t) into an infinite sum, or series, of simpler functions. Following Joseph Fourier (1768-1830) consider the infinite sum of sine and cosine functions 1 a0 X f(t) = + [a cos(nt) + b sin(nt)] (3) 2 n n n=1 where the constant coefficients an and bn are called the Fourier coefficients of f. The first question one would like to answer is how to find those coefficients. To do so we utilize the orthogonality of sine and cosine functions: Z π Z π 1 cos(nt) cos(mt) dt = [cos((m − n)t) + cos((m + n)t)] dt −π −π 2 8 2π; m = n = 0 <> = π; m = n 6= 0 :>0; m 6= n ( 2π; m = n = 0 = (4) πδmn; m 6= 0 2 General Introduction Similarly, Z π Z π 1 sin(nt) sin(mt) dt = [cos((m − n)t) − cos((m + n)t)] dt −π −π 2 ( 0 m = 0 = (5) πδmn m 6= 0 and Z π Z π 1 sin(nt) cos(mt) dt = [sin((m − n)t) + sin((m + n)t)] dt −π −π 2 = 0 (6) Using the orthogonality and the assumed expression for the infinite series given in Eq. (3), it follows that the Fourier coefficients are 1 Z π an = f(t) cos(nt) dt (7) π −π 1 Z π bn = f(t) sin(nt) dt (8) π −π This initial insight by Fourier was followed by centuries of a work on the second obvious question: Are the RHS and LHS in Eq. (3) actually the same? Clearly one needs to determine for which class of functions f the infinite series on the right hand side of Eq. (3) will converge. That is, what is a sufficiently nice function f? The precise answer is not of concern here, it suffices to know that the Fourier series exists and converges for periodic functions of the type you are used to, e.g. functions for which first and second order derivatives exists almost everywhere, that are finite and have at most a finite number of discontinuities and zero crossings in the interval (−π; π). When determining a the Fourier series of a periodic function f(t) with period T , any interval (t0; t0 + T ) can be used, with the choice being one of conve- nience or personal preference. For example, in the rescaled time coordinates considering the interval (0; 2π) works just as well as considering (−π; π) as we have done. If a function is even so that f(t) = f(−t), then f(t) sin(nt) is odd. (This follows since sin(nt) is odd and an even function times an odd function is an odd function.) Therefore, bn = 0 for all n. Similarly, if a function is odd so that f(t) = −f(−t), then f(t) cos(nt) is odd. (This follows since cos(nt) is even and an even function times an odd function is an odd function.) Therefore, an = 0 for all n. 3 1. Fourier Series ' $%& $%! $%" 123.4/51 $%# $ !! !" !# $ # " ! ()*+,-./0/ Figure 1: A full-wave-rectifier converts a sinusoidal input, sin(!t), to j sin(!t)j. Example - Rectified sine wave: A first step in converting AC-power from the power-grid to the DC-power that most devices need is to utilize a full-wave rectifier, such as the diode bridge shown in Fig. 1, which converts a sinusoidal input to an output that is the absolute value of the input sine-wave. One notes immediately that for a sinusoidal input, the output of the rectifier is periodic with half of the period of the input. The fundamental frequency of the output is twice the input frequency. How can that be? The reason is that the circuit is not a linear circuit. The presence of diodes makes this circuit nonlinear and allows the circuit to shift power from the fundamental frequency to twice its value. One might wonder whether that is all that is happening. Does the output have contributions (power) at other frequencies? To answer this we look at the Fourier series of the output. Since the output f = j sin(!t)j is even, i.e. f(t) = f(−t), no terms of the form sin(n!t) will appear in the answer. It suffices to determine the an coefficients. For a0 one obtains 1 Z 0 1 Z π a0 = − sin(!t) d(!t) + sin(!t) d(!t) π −π π 0 2 Z π 4 = sin(!t) d(!t) = (9) π 0 π 4 Discontinuous Functions and for the remaining an one gets 2 Z π an = sin(!t) cos(n!t) d(!t) π 0 2 Z π 1 = [− sin((n − 1)!t) + sin((n + 1)!t)] d(!t) π 0 2 1 1 −1 = fcos(nπ − π) − 1g + fcos(nπ + π) − 1g π n − 1 n + 1 ( − 4 1 ; n even = π n2−1 (10) 0; n odd: Note, that the sine and cosine functions are orthogonal on the interval (−π; π). They are not orthogonal on the interval (0; π) and we do get a nonzero contribution for even n. To summarize the result, 1 2 4 X cos(n!t) j sin(!t)j = − : (11) π π n2 − 1 n=2;4;6;::: For an input with frequency f0, the output has a DC-offset, the part that we really care about when building a DC-voltage supply. It has no contribution at f = f0. It does have contributions at frequencies 2f0; 4f0; 6f0;:::. 1.2 Discontinuous Functions In the above example, Eq. (11), the nth coefficient decreases as 1=n2. This decay of the coefficients is in contrast to the Fourier series of a square wave 1 4 X 1 f (t) = sin(n!t) (12) sw π n n=1;3;5;::: where the nth coefficient falls off as 1=n. This is true in general 1 1. If f(t) has discontinuities, the nth coefficient decreases as 1=n. The convergences is slow and many terms need to be kept to approximate such a function well. 1G. Raisbeck, Order of magnitude of Fourier coefficients. Am. Math. Mon. 62, 149-155 (1955). 5 1. Fourier Series Figure 2: The Gibbs phenomenon is an overshoot (or "ringing") of Fourier series and other eigenfunction series occurring at simple discontinuities. Shown in color are the first few partial sums of the square-wave Fourier series. (Source math- world.wolfram.com/GibbsPhenomenon.html) This means that a function generator that generates square waves through the addition of sinusoidal waveforms needs to have a bandwidth (max. freq. it can generate) that is large compared to the frequency of the square-wave that is generated. 2. If f(t) is continuous (although possibly with discontinuous derivatives) the nth coefficient decreases as 1=n2. There is another consequence of a discontinuity in f(t) that can cause trouble in practical applications, where one necessarily only adds a finite number of sinusoidal terms. The nth partial sum of the Fourier series of a piecewise continuously differentiable periodic function f behaves at a jump discontinuity in a peculiar manner. It has large oscillations near the jump, which might increase the maximum of the partial sum above that of the function itself. It turns out that the Fourier series exceeds the height of a square wave by about 9 percent. This is the so-called Gibbs phenomenon, shown in Fig. 2. Increasing the number of terms in the partial sum does not decrease the magnitude of the overshoot but moves the overshoot extremum point closer and closer to the jump discontinuity. You will have the opportunity to explore this in the lab. 6 Complex Fourier Series 1.3 Complex Fourier Series At this stage in your physics career you are all well acquainted with complex numbers and functions. Let us then generalize the Fourier series to complex functions. To motivate this, return to the Fourier series, Eq. (3): 1 a0 X f(t) = + [a cos(nt) + b sin(nt)] 2 n n n=1 1 int −int int −int a0 X e + e e − e = + a + b 2 n 2 n 2i n=1 1 −∞ a0 X an − ibn X a−m + ib−m = + eint + eimt (13) 2 2 2 n=1 m=−1 where we substituted m = −n in the last term on the last line. Equation (13) clearly suggests the much simpler complex form of the Fourier series +1 X in(2πf0)t x(t) = Xn e : (14) n=−∞ with the coefficients given by 1 Z T=2 −in(2πf0)t Xn = x(t) e dt (15) T −T=2 Here, the Fourier series is written for a complex periodic function x(t) with arbitrary period T = 1=f0.
Recommended publications
  • The Fast Fourier Transform
    The Fast Fourier Transform Derek L. Smith SIAM Seminar on Algorithms - Fall 2014 University of California, Santa Barbara October 15, 2014 Table of Contents History of the FFT The Discrete Fourier Transform The Fast Fourier Transform MP3 Compression via the DFT The Fourier Transform in Mathematics Table of Contents History of the FFT The Discrete Fourier Transform The Fast Fourier Transform MP3 Compression via the DFT The Fourier Transform in Mathematics Navigating the Origins of the FFT The Royal Observatory, Greenwich, in London has a stainless steel strip on the ground marking the original location of the prime meridian. There's also a plaque stating that the GPS reference meridian is now 100m to the east. This photo is the culmination of hundreds of years of mathematical tricks which answer the question: How to construct a more accurate clock? Or map? Or star chart? Time, Location and the Stars The answer involves a naturally occurring reference system. Throughout history, humans have measured their location on earth, in order to more accurately describe the position of astronomical bodies, in order to build better time-keeping devices, to more successfully navigate the earth, to more accurately record the stars... and so on... and so on... Time, Location and the Stars Transoceanic exploration previously required a vessel stocked with maps, star charts and a highly accurate clock. Institutions such as the Royal Observatory primarily existed to improve a nations' navigation capabilities. The current state-of-the-art includes atomic clocks, GPS and computerized maps, as well as a whole constellation of government organizations.
    [Show full text]
  • B1. Fourier Analysis of Discrete Time Signals
    B1. Fourier Analysis of Discrete Time Signals Objectives • Introduce discrete time periodic signals • Define the Discrete Fourier Series (DFS) expansion of periodic signals • Define the Discrete Fourier Transform (DFT) of signals with finite length • Determine the Discrete Fourier Transform of a complex exponential 1. Introduction In the previous chapter we defined the concept of a signal both in continuous time (analog) and discrete time (digital). Although the time domain is the most natural, since everything (including our own lives) evolves in time, it is not the only possible representation. In this chapter we introduce the concept of expanding a generic signal in terms of elementary signals, such as complex exponentials and sinusoids. This leads to the frequency domain representation of a signal in terms of its Fourier Transform and the concept of frequency spectrum so that we characterize a signal in terms of its frequency components. First we begin with the introduction of periodic signals, which keep repeating in time. For these signals it is fairly easy to determine an expansion in terms of sinusoids and complex exponentials, since these are just particular cases of periodic signals. This is extended to signals of a finite duration which becomes the Discrete Fourier Transform (DFT), one of the most widely used algorithms in Signal Processing. The concepts introduced in this chapter are at the basis of spectral estimation of signals addressed in the next two chapters. 2. Periodic Signals VIDEO: Periodic Signals (19:45) http://faculty.nps.edu/rcristi/eo3404/b-discrete-fourier-transform/videos/chapter1-seg1_media/chapter1-seg1- 0.wmv http://faculty.nps.edu/rcristi/eo3404/b-discrete-fourier-transform/videos/b1_02_periodicSignals.mp4 In this section we define a class of discrete time signals called Periodic Signals.
    [Show full text]
  • GIBBS PHENOMENON SUPPRESSION and OPTIMAL WINDOWING for ATTENUATION and Q MEASUREMENTS* Cheh Pan Department of Geophysics, Stanfo
    SLAC-PUB-6222 September 1993 (Ml GIBBS PHENOMENON SUPPRESSION AND OPTIMAL WINDOWING FOR ATTENUATION AND Q MEASUREMENTS* Cheh Pan Department of Geophysics, Stanford University Stanford, CA 94305-2215 and Stanford Linear Accelerator Center Stanford, CA 94309 INTRODUCTION There are basically four known techniques, spectral ratios (Hauge 1981, Johnston and Toksoz 1981, Moos 1984, Goldberg et al 1985, Patten 1988, Sams and Goldberg 1990), forward modeling (Chuen and Toksoz 1981), inversion (Cheng et al 1986), and first pulse rise time (Gladwin and Stacey 1974, Moos 1984), that have been used in the past to measure the attenuation coefficient or quality factor from seismic and acoustic data. Problems have been encountered in using the spectral techniques, which included: (1) the correction for geometrical divergence of the acoustic wave front; (2) the suppression of the , Gibbs phenomenon or the ringing effect in the spectra; and (3) the elimination of contamination from interfering wave modes. Geometrical corrections have been presented by Patten (1988) and Sams and Goldberg (1990) for the borehole acoustics case. The remaining difficulties in the application of the spectral ratios technique come mainly from the suppression of the Gibbs phenomenon and optimal windowing of wave modes. This report deals with these two problems. e - Submitted to Geophysics * Work supP;orted by Department of Energy contract DE-ACO3-76SFOO5 15 : I FUNDAMENTALS The amplitudes R&j) of a seismic signal at frequencyfrecorded by a receiver at a distance or offset x from the source can be represented as, R(x,f)= A(f)G( (1) where A(/) is the source term, G(x) is the geometrical divergence which is assumed to be independent of frequency, and a is the attenuation coefficient.
    [Show full text]
  • Fourier Series
    Academic Press Encyclopedia of Physical Science and Technology Fourier Series James S. Walker Department of Mathematics University of Wisconsin–Eau Claire Eau Claire, WI 54702–4004 Phone: 715–836–3301 Fax: 715–836–2924 e-mail: [email protected] 1 2 Encyclopedia of Physical Science and Technology I. Introduction II. Historical background III. Definition of Fourier series IV. Convergence of Fourier series V. Convergence in norm VI. Summability of Fourier series VII. Generalized Fourier series VIII. Discrete Fourier series IX. Conclusion GLOSSARY ¢¤£¦¥¨§ Bounded variation: A function has bounded variation on a closed interval ¡ ¢ if there exists a positive constant © such that, for all finite sets of points "! "! $#&% (' #*) © ¥ , the inequality is satisfied. Jordan proved that a function has bounded variation if and only if it can be expressed as the difference of two non-decreasing functions. Countably infinite set: A set is countably infinite if it can be put into one-to-one £0/"£ correspondence with the set of natural numbers ( +,£¦-.£ ). Examples: The integers and the rational numbers are countably infinite sets. "! "!;: # # 123547698 Continuous function: If , then the function is continuous at the point : . Such a point is called a continuity point for . A function which is continuous at all points is simply referred to as continuous. Lebesgue measure zero: A set < of real numbers is said to have Lebesgue measure ! $#¨CED B ¢ £¦¥ zero if, for each =?>A@ , there exists a collection of open intervals such ! ! D D J# K% $#L) ¢ £¦¥ ¥ ¢ = that <GFIH and . Examples: All finite sets, and all countably infinite sets, have Lebesgue measure zero. "! "! % # % # Odd and even functions: A function is odd if for all in its "! "! % # # domain.
    [Show full text]
  • FOURIER TRANSFORM Very Broadly Speaking, the Fourier Transform Is a Systematic Way to Decompose “Generic” Functions Into
    FOURIER TRANSFORM TERENCE TAO Very broadly speaking, the Fourier transform is a systematic way to decompose “generic” functions into a superposition of “symmetric” functions. These symmetric functions are usually quite explicit (such as a trigonometric function sin(nx) or cos(nx)), and are often associated with physical concepts such as frequency or energy. What “symmetric” means here will be left vague, but it will usually be associated with some sort of group G, which is usually (though not always) abelian. Indeed, the Fourier transform is a fundamental tool in the study of groups (and more precisely in the representation theory of groups, which roughly speaking describes how a group can define a notion of symmetry). The Fourier transform is also related to topics in linear algebra, such as the representation of a vector as linear combinations of an orthonormal basis, or as linear combinations of eigenvectors of a matrix (or a linear operator). To give a very simple prototype of the Fourier transform, consider a real-valued function f : R → R. Recall that such a function f(x) is even if f(−x) = f(x) for all x ∈ R, and is odd if f(−x) = −f(x) for all x ∈ R. A typical function f, such as f(x) = x3 + 3x2 + 3x + 1, will be neither even nor odd. However, one can always write f as the superposition f = fe + fo of an even function fe and an odd function fo by the formulae f(x) + f(−x) f(x) − f(−x) f (x) := ; f (x) := . e 2 o 2 3 2 2 3 For instance, if f(x) = x + 3x + 3x + 1, then fe(x) = 3x + 1 and fo(x) = x + 3x.
    [Show full text]
  • Configuring Spectrogram Views
    Chapter 5: Spectrographic Analysis Configuring spectrographic views To create a new spectrogram, spectrogram slice, or selection spectrum view, click on the appropriate New View button in the view toolbar (Figure 5.3) or choose a view type from the View > New menu. Figure 5.3. New View buttons Figure 5.3. The New View buttons, in Raven’s view toolbar. 112 Raven Pro 1.4 User’s Manual Chapter 5: Spectrographic Analysis A dialog box appears, containing parameters for configuring the requested type of spectrographic view (Figure 5.4). The dialog boxes for configuring spectrogram and spectrogram slice views are identical, except for their titles. The dialog boxes are identical because both view types calculate a spectrogram of the entire sound; the only difference between spectrogram and spectrogram slice views is in how the data are displayed (see “How the spectrographic views are related” on page 110). The dialog box for configuring a selection spectrum view is the same, except that it lacks the Averaging parameter. The remainder of this section explains each of the parameters in the configuration dialog box. Figure 5.4 Configure Spectrogram dialog Figure 5.4. The Configure New Spectrogram dialog box. Window type Each data record is multiplied by a window function before its spectrum is calculated. Window functions are used to reduce the magnitude of spurious “sidelobe” energy that appears at frequencies flanking each analysis frequency in a spectrum. These sidelobes appear as a result of Raven Pro 1.4 User’s Manual 113 Chapter 5: Spectrographic Analysis analyzing a finite (truncated) portion of a signal.
    [Show full text]
  • An Introduction to Fourier Analysis Fourier Series, Partial Differential Equations and Fourier Transforms
    An Introduction to Fourier Analysis Fourier Series, Partial Differential Equations and Fourier Transforms Notes prepared for MA3139 Arthur L. Schoenstadt Department of Applied Mathematics Naval Postgraduate School Code MA/Zh Monterey, California 93943 August 18, 2005 c 1992 - Professor Arthur L. Schoenstadt 1 Contents 1 Infinite Sequences, Infinite Series and Improper Integrals 1 1.1Introduction.................................... 1 1.2FunctionsandSequences............................. 2 1.3Limits....................................... 5 1.4TheOrderNotation................................ 8 1.5 Infinite Series . ................................ 11 1.6ConvergenceTests................................ 13 1.7ErrorEstimates.................................. 15 1.8SequencesofFunctions.............................. 18 2 Fourier Series 25 2.1Introduction.................................... 25 2.2DerivationoftheFourierSeriesCoefficients.................. 26 2.3OddandEvenFunctions............................. 35 2.4ConvergencePropertiesofFourierSeries.................... 40 2.5InterpretationoftheFourierCoefficients.................... 48 2.6TheComplexFormoftheFourierSeries.................... 53 2.7FourierSeriesandOrdinaryDifferentialEquations............... 56 2.8FourierSeriesandDigitalDataTransmission.................. 60 3 The One-Dimensional Wave Equation 70 3.1Introduction.................................... 70 3.2TheOne-DimensionalWaveEquation...................... 70 3.3 Boundary Conditions ............................... 76 3.4InitialConditions................................
    [Show full text]
  • Lecture 11 : Discrete Cosine Transform Moving Into the Frequency Domain
    Lecture 11 : Discrete Cosine Transform Moving into the Frequency Domain Frequency domains can be obtained through the transformation from one (time or spatial) domain to the other (frequency) via Fourier Transform (FT) (see Lecture 3) — MPEG Audio. Discrete Cosine Transform (DCT) (new ) — Heart of JPEG and MPEG Video, MPEG Audio. Note : We mention some image (and video) examples in this section with DCT (in particular) but also the FT is commonly applied to filter multimedia data. External Link: MIT OCW 8.03 Lecture 11 Fourier Analysis Video Recap: Fourier Transform The tool which converts a spatial (real space) description of audio/image data into one in terms of its frequency components is called the Fourier transform. The new version is usually referred to as the Fourier space description of the data. We then essentially process the data: E.g . for filtering basically this means attenuating or setting certain frequencies to zero We then need to convert data back to real audio/imagery to use in our applications. The corresponding inverse transformation which turns a Fourier space description back into a real space one is called the inverse Fourier transform. What do Frequencies Mean in an Image? Large values at high frequency components mean the data is changing rapidly on a short distance scale. E.g .: a page of small font text, brick wall, vegetation. Large low frequency components then the large scale features of the picture are more important. E.g . a single fairly simple object which occupies most of the image. The Road to Compression How do we achieve compression? Low pass filter — ignore high frequency noise components Only store lower frequency components High pass filter — spot gradual changes If changes are too low/slow — eye does not respond so ignore? Low Pass Image Compression Example MATLAB demo, dctdemo.m, (uses DCT) to Load an image Low pass filter in frequency (DCT) space Tune compression via a single slider value n to select coefficients Inverse DCT, subtract input and filtered image to see compression artefacts.
    [Show full text]
  • Parallel Fast Fourier Transform Transforms
    Parallel Fast Fourier Parallel Fast Fourier Transform Transforms Massimiliano Guarrasi – [email protected] MassimilianoSuper Guarrasi Computing Applications and Innovation Department [email protected] Fourier Transforms ∞ H ()f = h(t)e2πift dt ∫−∞ ∞ h(t) = H ( f )e−2πift df ∫−∞ Frequency Domain Time Domain Real Space Reciprocal Space 2 of 49 Discrete Fourier Transform (DFT) In many application contexts the Fourier transform is approximated with a Discrete Fourier Transform (DFT): N −1 N −1 ∞ π π π H ()f = h(t)e2 if nt dt ≈ h e2 if ntk ∆ = ∆ h e2 if ntk n ∫−∞ ∑ k ∑ k k =0 k=0 f = n / ∆ = ∆ n tk k / N = ∆ fn n / N −1 () = ∆ 2πikn / N H fn ∑ hk e k =0 The last expression is periodic, with period N. It define a ∆∆∆ between 2 sets of numbers , Hn & hk ( H(f n) = Hn ) 3 of 49 Discrete Fourier Transforms (DFT) N −1 N −1 = 2πikn / N = 1 −2πikn / N H n ∑ hk e hk ∑ H ne k =0 N n=0 frequencies from 0 to fc (maximum frequency) are mapped in the values with index from 0 to N/2-1, while negative ones are up to -fc mapped with index values of N / 2 to N Scale like N*N 4 of 49 Fast Fourier Transform (FFT) The DFT can be calculated very efficiently using the algorithm known as the FFT, which uses symmetry properties of the DFT s um. 5 of 49 Fast Fourier Transform (FFT) exp(2 πi/N) DFT of even terms DFT of odd terms 6 of 49 Fast Fourier Transform (FFT) Now Iterate: Fe = F ee + Wk/2 Feo Fo = F oe + Wk/2 Foo You obtain a series for each value of f n oeoeooeo..oe F = f n Scale like N*logN (binary tree) 7 of 49 How to compute a FFT on a distributed memory system 8 of 49 Introduction • On a 1D array: – Algorithm limits: • All the tasks must know the whole initial array • No advantages in using distributed memory systems – Solutions: • Using OpenMP it is possible to increase the performance on shared memory systems • On a Multi-Dimensional array: – It is possible to use distributed memory systems 9 of 49 Multi-dimensional FFT( an example ) 1) For each value of j and k z Apply FFT to h( 1..
    [Show full text]
  • Evaluating Fourier Transforms with MATLAB
    ECE 460 – Introduction to Communication Systems MATLAB Tutorial #2 Evaluating Fourier Transforms with MATLAB In class we study the analytic approach for determining the Fourier transform of a continuous time signal. In this tutorial numerical methods are used for finding the Fourier transform of continuous time signals with MATLAB are presented. Using MATLAB to Plot the Fourier Transform of a Time Function The aperiodic pulse shown below: x(t) 1 t -2 2 has a Fourier transform: X ( jf ) = 4sinc(4π f ) This can be found using the Table of Fourier Transforms. We can use MATLAB to plot this transform. MATLAB has a built-in sinc function. However, the definition of the MATLAB sinc function is slightly different than the one used in class and on the Fourier transform table. In MATLAB: sin(π x) sinc(x) = π x Thus, in MATLAB we write the transform, X, using sinc(4f), since the π factor is built in to the function. The following MATLAB commands will plot this Fourier Transform: >> f=-5:.01:5; >> X=4*sinc(4*f); >> plot(f,X) In this case, the Fourier transform is a purely real function. Thus, we can plot it as shown above. In general, Fourier transforms are complex functions and we need to plot the amplitude and phase spectrum separately. This can be done using the following commands: >> plot(f,abs(X)) >> plot(f,angle(X)) Note that the angle is either zero or π. This reflects the positive and negative values of the transform function. Performing the Fourier Integral Numerically For the pulse presented above, the Fourier transform can be found easily using the table.
    [Show full text]
  • Fourier Transforms & the Convolution Theorem
    Convolution, Correlation, & Fourier Transforms James R. Graham 11/25/2009 Introduction • A large class of signal processing techniques fall under the category of Fourier transform methods – These methods fall into two broad categories • Efficient method for accomplishing common data manipulations • Problems related to the Fourier transform or the power spectrum Time & Frequency Domains • A physical process can be described in two ways – In the time domain, by h as a function of time t, that is h(t), -∞ < t < ∞ – In the frequency domain, by H that gives its amplitude and phase as a function of frequency f, that is H(f), with -∞ < f < ∞ • In general h and H are complex numbers • It is useful to think of h(t) and H(f) as two different representations of the same function – One goes back and forth between these two representations by Fourier transforms Fourier Transforms ∞ H( f )= ∫ h(t)e−2πift dt −∞ ∞ h(t)= ∫ H ( f )e2πift df −∞ • If t is measured in seconds, then f is in cycles per second or Hz • Other units – E.g, if h=h(x) and x is in meters, then H is a function of spatial frequency measured in cycles per meter Fourier Transforms • The Fourier transform is a linear operator – The transform of the sum of two functions is the sum of the transforms h12 = h1 + h2 ∞ H ( f ) h e−2πift dt 12 = ∫ 12 −∞ ∞ ∞ ∞ h h e−2πift dt h e−2πift dt h e−2πift dt = ∫ ( 1 + 2 ) = ∫ 1 + ∫ 2 −∞ −∞ −∞ = H1 + H 2 Fourier Transforms • h(t) may have some special properties – Real, imaginary – Even: h(t) = h(-t) – Odd: h(t) = -h(-t) • In the frequency domain these
    [Show full text]
  • Fourier Analysis
    Chapter 1 Fourier analysis In this chapter we review some basic results from signal analysis and processing. We shall not go into detail and assume the reader has some basic background in signal analysis and processing. As basis for signal analysis, we use the Fourier transform. We start with the continuous Fourier transformation. But in applications on the computer we deal with a discrete Fourier transformation, which introduces the special effect known as aliasing. We use the Fourier transformation for processes such as convolution, correlation and filtering. Some special attention is given to deconvolution, the inverse process of convolution, since it is needed in later chapters of these lecture notes. 1.1 Continuous Fourier Transform. The Fourier transformation is a special case of an integral transformation: the transforma- tion decomposes the signal in weigthed basis functions. In our case these basis functions are the cosine and sine (remember exp(iφ) = cos(φ) + i sin(φ)). The result will be the weight functions of each basis function. When we have a function which is a function of the independent variable t, then we can transform this independent variable to the independent variable frequency f via: +1 A(f) = a(t) exp( 2πift)dt (1.1) −∞ − Z In order to go back to the independent variable t, we define the inverse transform as: +1 a(t) = A(f) exp(2πift)df (1.2) Z−∞ Notice that for the function in the time domain, we use lower-case letters, while for the frequency-domain expression the corresponding uppercase letters are used. A(f) is called the spectrum of a(t).
    [Show full text]