The Fourier series The Fourier series A large class of phenomena can be described as periodic in nature: waves, sounds, light, radio, water waves etc . It is natural to attempt to describe these phenomena by means of expansions in periodic functions . The Fourier series is an expansion of a function in terms of trigonometric sines and cosines : Suppose f(x) is defined over a finite range –L x L, i.e. f(x) is periodic with period 2L. The trigonometric functions are periodic with period 2 L = 2 , so it is natural to expand these functions in terms of trigonometric functions with an argument [(x / 2L) 2 n], n : ) sin( ) cos( 2 1 ) ( 1 1 0 x L n b x L n a a x f m m m m dx x L n x f L a L L n ) cos( 1 dx x L n x f L b L L n ) sin( 1
22
Embed
The Fourier series A large class of phenomena can be described as periodic in nature: waves, sounds, light, radio, water waves etc. It is natural to attempt.
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
The Fourier seriesThe Fourier series
A large class of phenomena can be described as periodic in nature: waves, sounds, light, radio, water waves etc.
It is natural to attempt to describe these phenomena by means of expansions in periodic functions.
The Fourier series is an expansion of a function in terms of trigonometric sines and cosines:
Suppose f(x) is defined over a finite range –L x L, i.e. f(x) is periodic with period 2L. The trigonometric functions are periodic with period 2L = 2 , so it is natural to expand these functions in terms of trigonometric functions with
an argument [(x / 2L) 2 n], n:
)sin()cos(2
1)(
110 x
L
nbx
L
naaxf
mm
mm
dxxL
nxf
La
L
Ln )cos(1 dxx
L
nxf
Lb
L
Ln )sin(1
Example for Fourier seriesExample for Fourier series
Let’s take the following function:
| x| ; - 8 <= x <= 8 f(x 16n) , n ; otherwise
-f(x) is even in x while sin(x) is odd => bn’s must be zero .
; if n is odd
; if n is even
{f(x)=
22220
4]1)[cos(
2)cos(
2
n
Ln
n
Ldxx
L
nx
La
L
n
0]1)[cos(2
)cos(2
220
nn
Ldxx
L
nx
La
L
n
>=the expansion is: ))cos(18
1(2
)(22
xL
n
n
Lxf
oddn
Example for Fourier series(2)Example for Fourier series(2)
f(x)
x
y
n= 1
n= 3
n= 5 ))cos(18
1(2
)(22
xL
n
n
Lxf
oddn
source: f(x) =|x|, -8 >= x >= 8
Successive approximations of f(x)
Complex version of the Fourier expansionComplex version of the Fourier expansion
The Euler identity :
The inverse equations:
,
sincos iei
)(sin 21 iii ee )(cos 2
1 ii ee
Using the formulas above and some properties of exponential function, the Fourier series can also be written as an expansion in terms of complex
The discrete Fourier transformThe discrete Fourier transform
Motivation: computer applications of the Fourier transform require that all of the definitions and properties of Fourier transforms be translated into analogous statements appropriate to functions represented by a discrete set of sampling points rather than by continuous functions.
Let f(x) be a function.Let {fk = f(xk){ be a set of N function values, k = 0, 1, …, N-1.Let be the separation of the equidistant sampling points.Assumption: N is even.
The discrete Fourier transform is :
The inverse discrete transform is :
xkxk x
1,...,1,0,)(0
/2^
NnfefN
kk
nkNin
1,...,1,0,)(1
1
0
^/2
NkfeN
fN
nn
nkNik
The discrete Fourier transform(2)The discrete Fourier transform(2)
1,...,1,0,)(0
/2^
NnfefN
kk
nkNin
Let’s examine more closely the formula of the discrete Fourier transform :
We know that (it’s called n-th root of unity), so the formula above can be rewritten as:
N
iwN
2
1,...,1,0,0
^
NnfwfN
kk
knn
)(^^
nnn wff
Usage of the Usage of the FFast ast FFourier ourier TTransform - motivationransform - motivation
Let be two polynomials.
1
0
)(n
k
kk xaxA
1
0
)(n
k
kk xbxB,
We want to multiply them: C(x) = A(x)*B(x).
Two ways to do this: 1. where - takes time
2) .a (calculate A(x) and B(x) values in 2n-1 distinct points x0, …, x2n-2; we get two vectors {(x0, y0),…,(x2n-2, y2n-2){, {(x0,z0),…,(x2n-2, z 2n-2){; (b) multiply these vectors: {(x0, y0 z0),…,(x2n-2, y2n-2 z 2n-2){ (c) calculate the polynomial C(x) that passes through the result vector (interpolation) - takes time if use FFT
,)(22
0
n
j
jj xcxC )( 2n
j
kkjkj bac
0
)log( nn
Uniqueness of Uniqueness of C(x)C(x)
Theorem1: for any set {(x0, y0),…,(xn-\, yn-1){ on n distinct points there is a unique polynomial C(x) with degree less than n such that yi = C(xi) for i = 0,1,…, n-1.
Proof: we can write yi = C(xi) = for i = 0,1,…, n-1 as matrices multiplication :
1
0
n
k
kik xc
c0 c1
... cn-1
=
y0 y1
... yn-1
*
The matrix V(x0 ,..., xn-\ ) is called a Vandermonde matrix.Since all x0 ,..., xn-\ are distinct, a discriminate of V isn’t zero, so it is reversible.=> we can calculate ci = V (x0 ,..., xn-\ ) -1 yi .
Discrete FFT and FFTinverseDiscrete FFT and FFTinverse
We will use FFT to execute step (a) and FFTinv to execute step (c).
Let be a polynomial.
We want to execute a step (a): calculate step (a): calculate C(x)C(x) in in nn distinct points distinct points xx00, …, x, …, xn-1n-1..FFT uses special n points – wn
0, wn1 , wn
2 , …, wnn-1.
wn is called the n-th (complex) root of unity.That means that wn
n = 1.
1
0
)(n
k
kk xcxC
Properties of n-th root of unityProperties of n-th root of unity
Let’s examine some properties of wn :
11.. There are exactly n n-th roots of unity: wn0, wn
1 , wn2 , …, wn
n-1.
Each one of them can be presented as . wn is called the principal n-th root of unity.
22. The inverse of wn is wn-1 = wn
n-1: wn * wn-1 = wn
0 = 1; wn * wn
n-1 = wnn = 1.
33 .wdndk = wn
k : wdndk = = = wn
k .
44. wnn/2 = w2 = -1 : wn
n/2 = w(n/2)*2n/2 = (by property 3) w2 = = -1.
nie /2
dkdnie )( /2 knie )( /2
2/2 ie
Properties of n-th root of unity(2)Properties of n-th root of unity(2)
55 .If n > 0 and n is even => (wnk ) 2 = wn/2
k , k = 0, 1, …, n–1:
) .wnk (2 = = = wn/2
k .
) .wnk+n/2 (2 = wn
2k+n = wn2k wn
n = wn2k = (wn
k ) 2 = wn/2k .
66 . : .
knie 2/2 )( knie )( 2//2
geometry series formula
DFFT previewDFFT preview
Let be a polynomial.
We want to execute a step (a): calculate step (a): calculate ppcc(x)(x) in in wn0, wn
1 , wn2 , …, wn
n-1.
Observation: Let be an n-tuple of the coefficients of ppcc(x).Assume that n is power of 2 (if it isn’t, we force it to be by adding 0’s).
Let and
>= , where is the polynomial that is defined by the vector a, and is the polynomial that is defined by the
vector b .
>=
Let t = n/2 => wnt = wn
n/2 = -1 .
>=for i < n/2 :
for i >= n/2 :
1
0
)(n
k
kkc xcxp
),...,,( 110 ncccc
),...,,( 220 nccca ),...,,( 131 ncccb
apbp)()()( 22 xxpxpxp bac
)()()( 22 inb
in
ina
inc wpwwpwp
)()()()()( 2222 inb
jn
ina
inb
tjn
ina
inc wpwwpwpwwpwp
)()()( 22 inb
in
ina
inc wpwwpwp
DFFT algorithmDFFT algorithmDFFT ( , wn) : // n is a power of 2 and wn
n/2 = -1 // array C[0, …, n – 1] // for the answer if n = = 1 then C[0] c[0] else t n / 2 arrays a, b, A, B [0,…,t – 1] // intermediate arrays for i 0 to t – 1 do: a[i] c[2i] b[i] c[2i + 1] // recursive Fourier transform computation A DFFT (a, wn
2 ) B DFFT (b, wn
2 ) // Fourier transform computation of the vector c for i 0 to t – 1 do: temp = wn
We can write yi = C(wni) = for i = 0,1,…, n-1 as matrices multiplication :
1
0
n
k
iknk wc
c0 c1
... cn-1
=
y0 y1
... yn-1
*
The matrix Vij(x0 ,..., xn-\ ) is called a Vandermonde matrix.Since all x0 ,..., xn-\ are distinct, a discriminate of Vij isn’t zero, so it is reversible.=> we can calculate ci = Vij (x0 ,..., xn-\ ) -1 yi .
We want to execute a step step (c) calculate the polynomial C(x) that passes through the result vector (interpolation)