Fourier Series and Transforms Orthogonal functions Fourier Series
- Slides: 30
Fourier Series and Transforms • • • Orthogonal functions Fourier Series Discrete Fourier Series Fourier Transform Chebyshev polynomials Scope: we are trying to approximate an arbitrary function and obtain basis functions with appropriate coefficients. Fourier Theory Modern Seismology – Data processing and inversion 1
Fourier Series The Problem we are trying to approximate a function f(x) by another function gn(x) which consists of a sum over N orthogonal functions F(x) weighted by some coefficients an. Fourier Theory Modern Seismology – Data processing and inversion 2
The Problem. . . and we are looking for optimal functions in a least squares (l 2) sense. . . a good choice for the basis functions F(x) are orthogonal functions. What are orthogonal functions? Two functions f and g are said to be orthogonal in the interval [a, b] if How is this related to the more conceivable concept of orthogonal vectors? Let us look at the original definition of integrals: Fourier Theory Modern Seismology – Data processing and inversion 3
Orthogonal Functions . . . where x 0=a and x. N=b, and xi-xi-1= x. . . If we interpret f(xi) and g(xi) as the ith components of an N component vector, then this sum corresponds directly to a scalar product of vectors. The vanishing of the scalar product is the condition for orthogonality of vectors (or functions). fi Fourier Theory gi Modern Seismology – Data processing and inversion 4
Periodic functions Let us assume we have a piecewise continuous function of the form . . . we want to approximate this function with a linear combination of 2 periodic functions: Fourier Theory Modern Seismology – Data processing and inversion 5
Orthogonality. . . are these functions orthogonal ? . . . YES, and these relations are valid for any interval of length 2. Now we know that this is an orthogonal basis, but how can we obtain the coefficients for the basis functions? from minimising f(x)-g(x) Fourier Theory Modern Seismology – Data processing and inversion 6
Fourier coefficients optimal functions g(x) are given if . . . with the definition of g(x) we get. . . leading to Fourier Theory Modern Seismology – Data processing and inversion 7
Fourier approximation of |x|. . . Example. . . leads to the Fourier Serie . . and for n<4 g(x) looks like Fourier Theory Modern Seismology – Data processing and inversion 8
Fourier approximation of x 2. . . another Example. . . leads to the Fourier Serie . . and for N<11, g(x) looks like Fourier Theory Modern Seismology – Data processing and inversion 9
Fourier - discrete functions. . . what happens if we know our function f(x) only at the points it turns out that in this particular case the coefficients are given by . . the so-defined Fourier polynomial is the unique interpolating function to the function f(xj) with N=2 m Fourier Theory Modern Seismology – Data processing and inversion 10
Fourier - collocation points. . . with the important property that. . . in our previous examples. . . f(x)=|x| => f(x) - blue ; g(x) - red; xi - ‘+’ Fourier Theory Modern Seismology – Data processing and inversion 11
Fourier series - convergence f(x)=x 2 => f(x) - blue ; g(x) - red; xi - ‘+’ Fourier Theory Modern Seismology – Data processing and inversion 12
Fourier series - convergence f(x)=x 2 => f(x) - blue ; g(x) - red; xi - ‘+’ Fourier Theory Modern Seismology – Data processing and inversion 13
Gibb’s phenomenon f(x)=x 2 => f(x) - blue ; g(x) - red; xi - ‘+’ The overshoot for equispaced Fourier interpolations is 14% of the step height. Fourier Theory Modern Seismology – Data processing and inversion 14
Chebyshev polynomials We have seen that Fourier series are excellent for interpolating (and differentiating) periodic functions defined on a regularly spaced grid. In many circumstances physical phenomena which are not periodic (in space) and occur in a limited area. This quest leads to the use of Chebyshev polynomials. We depart by observing that cos(n ) can be expressed by a polynomial in cos( ): . . . which leads us to the definition: Fourier Theory Modern Seismology – Data processing and inversion 15
Chebyshev polynomials - definition. . . for the Chebyshev polynomials Tn(x). Note that because of x=cos( ) they are defined in the interval [-1, 1] (which - however can be extended to ). The first polynomials are Fourier Theory Modern Seismology – Data processing and inversion 16
Chebyshev polynomials - Graphical The first ten polynomials look like [0, -1] The n-th polynomial has extrema with values 1 or -1 at Fourier Theory Modern Seismology – Data processing and inversion 17
Chebyshev collocation points These extrema are not equidistant (like the Fourier extrema) k x(k) Fourier Theory Modern Seismology – Data processing and inversion 18
Chebyshev polynomials - orthogonality. . . are the Chebyshev polynomials orthogonal? Chebyshev polynomials are an orthogonal set of functions in the interval [-1, 1] with respect to the weight function such that . . . this can be easily verified noting that Fourier Theory Modern Seismology – Data processing and inversion 19
Chebyshev polynomials - interpolation. . . we are now faced with the same problem as with the Fourier series. We want to approximate a function f(x), this time not a periodical function but a function which is defined between [-1, 1]. We are looking for gn(x) . . . and we are faced with the problem, how we can determine the coefficients ck. Again we obtain this by finding the extremum (minimum) Fourier Theory Modern Seismology – Data processing and inversion 20
Chebyshev polynomials - interpolation. . . to obtain. . . surprisingly these coefficients can be calculated with FFT techniques, noting that . . . and the fact that f(cos ) is a 2 -periodic function. . . which means that the coefficients ck are the Fourier coefficients ak of the periodic function F( )=f(cos )! Fourier Theory Modern Seismology – Data processing and inversion 21
Chebyshev - discrete functions. . . what happens if we know our function f(x) only at the points in this particular case the coefficients are given by . . . leading to the polynomial. . . with the property Fourier Theory Modern Seismology – Data processing and inversion 22
Chebyshev - collocation points - |x| f(x)=|x| => f(x) - blue ; gn(x) - red; xi - ‘+’ 8 points 16 points Fourier Theory Modern Seismology – Data processing and inversion 23
Chebyshev - collocation points - |x| f(x)=|x| => f(x) - blue ; gn(x) - red; xi - ‘+’ 32 points 128 points Fourier Theory Modern Seismology – Data processing and inversion 24
Chebyshev - collocation points - x 2 f(x)=x 2 => f(x) - blue ; gn(x) - red; xi - ‘+’ 8 points The interpolating function gn(x) was shifted by a small amount to be visible at all! 64 points Fourier Theory Modern Seismology – Data processing and inversion 25
Chebyshev vs. Fourier - numerical Chebyshev Fourier f(x)=x 2 => f(x) - blue ; g. N(x) - red; xi - ‘+’ This graph speaks for itself ! Gibb’s phenomenon with Chebyshev? Fourier Theory Modern Seismology – Data processing and inversion 26
Chebyshev vs. Fourier - Gibb’s Chebyshev Fourier f(x)=sign(x- ) => f(x) - blue ; g. N(x) - red; xi - ‘+’ Gibb’s phenomenon with Chebyshev? YES! Fourier Theory Modern Seismology – Data processing and inversion 27
Chebyshev vs. Fourier - Gibb’s Chebyshev Fourier f(x)=sign(x- ) => f(x) - blue ; g. N(x) - red; xi - ‘+’ Fourier Theory Modern Seismology – Data processing and inversion 28
Fourier vs. Chebyshev Fourier collocation points periodic functions domain limited area [-1, 1] basis functions interpolating function Fourier Theory Modern Seismology – Data processing and inversion 29
Fourier vs. Chebyshev (cont’d) Chebyshev Fourier coefficients • Gibb’s phenomenon for discontinuous functions • limited area calculations some properties • Efficient calculation via FFT • infinite domain through periodicity • grid densification at boundaries • coefficients via FFT • excellent convergence at boundaries • Gibb’s phenomenon Fourier Theory Modern Seismology – Data processing and inversion 30
- Orthogonality of trigonometric functions
- Fourier series and orthogonal functions
- Orthogonal functions in fourier series
- Lie2d
- Fourier series and integrals
- Discontinuous forcing functions
- Fourier sine series
- Orthogonal series expansion
- Fourier series trigonometric form
- Fourier sine and cosine series examples
- In half range cosine series we assume the function to be
- What advice does agamemnon give odysseus?
- Tabel transformasi laplace invers
- Transforming data into information
- This transforms a bare stage into the world of the play
- What is drama
- Z transform time shifting property
- Transforms eroded parts of earth's surface into lakes
- Energy transformation and conservation
- Which phase transforms srs document
- Transforms of derivatives
- Friction transforms mechanical energy to
- Z transform of delta function
- Image transforms
- Digital image processing
- Photosynthesis transforms light energy into chemical energy
- Series de fourier
- Formula series de fourier
- Series de fourier complejas
- Frequency domain to time domain
- Inverse fourier transform table