Fourier series for periodic signals • Example 1: Determine the spectra of the signal... Fourier series for periodic signals • Solution of example 1: 6... Fourier transform of aperiodic s
Trang 1Xử lý tín hiệu số Fourier Transform
Trang 2• Frequency analysis of discrete time signal
• Properties of Fourier transform
• Frequency domain characteristics of LTI systems
• Discrete Fourier Transform
• Fast Fourier Transform
2
Trang 31.Frequency analysis of discrete time signal
1.1 Fourier series for periodic signals
– Given a periodic signal x(n) with period N
(DTFS)
Trang 41.1 Fourier series for periodic signals
• The spectrum of a signal x(n) which is periodic with period N,
is a periodic sequence with period N
4
Trang 51.Frequency analysis of discrete time signal
1.1 Fourier series for periodic signals
• Example 1: Determine the spectra of the signal
Trang 61.1 Fourier series for periodic signals
• Solution of example 1:
6
Trang 81.2 Fourier transform of aperiodic signals
• The Fourier transform of a finite energy signal x(n) is defined
as
• X(w) is periodic with period 2𝜋:
8
Trang 91.Frequency analysis of discrete time signal
1.2 Fourier transform of aperiodic signals
• In summary, the Fourier transform pair of a discrete time is as follows
• Uniform convergence is guaranteed if x(n) is absolutely
summable
Trang 101.2 Fourier transform of aperiodic signals
• The spectrum X(w) is, in general, a complex valued function
of frequency
• The energy density spectrum of x(n) is
10
Trang 111.Frequency analysis of discrete time signal
1.2 Fourier transform of aperiodic signals
• Example 2:
Trang 121.2 Fourier transform of aperiodic signals
• Solution of example 2:
12
Trang 131.Frequency analysis of discrete time signal
1.2 Fourier transform of aperiodic signals
• Solution of example 2 (cont’d):
Trang 15
2 Properties of Fourier transform
Trang 1616
Trang 172 Properties of Fourier transform
• Example 3:
Trang 18• Solution of Example 3:
18
Trang 192 Properties of Fourier transform
• Solution of Example 3:
Trang 20• Solution of Example 3 (cont’d):
20
a=0.8
Trang 212 Properties of Fourier transform
• Example 4:
Trang 22• Solution of Example 4:
22
Trang 233 Frequency domain characteristics of LTI
systems
• The response of any relaxed-system to arbitrary input signal is:
• Excite the system with the complex exponential
• Obtain the response
Trang 24• The Fourier transform of the unit sample response h(k) of the system
• The function H(𝜔) exists if the system is BIBO stable
• The response of the system to the complex exponential is
24
Trang 253 Frequency domain characteristics of LTI
systems
• Example
Trang 273 Frequency domain characteristics of LTI
systems
• In general, H(𝜔) is a complex function of 𝜔, hence
• Expressed in terms of real and image components
Trang 28• Example
28
Trang 293 Frequency domain characteristics of LTI
systems
• The impulse response
• It follows that
• Hence
Trang 3030
Trang 313 Frequency domain characteristics of LTI
systems
• Example
• Determine the response of the system (with given impulse
response h(n)) to the input signal x(n)
Trang 32• Solution
• The frequency response
• With each term in the input signal
• The response to the input signal
32
Trang 334 Discrete Fourier Transform
4.1 Frequency domain sampling
• Aperiodic finite- energy signals have continuous spectra
• Sample X(𝜔) periodically at a spacing of 𝛿𝜔 radians, take N equidistant samples in the interval 0 ≤ 𝜔 ≤ 2𝜋 with spacing
𝜔 = 2𝜋/𝑁
Trang 344.1 Frequency domain sampling
34
Trang 354 Discrete Fourier Transform
4.2 Discrete Fourier Transform
• A finite-duration sequence x(n) of length L has the DFT
where N ≥ L
• Recover x(n) from its DFT (inverse DFT)
Trang 36• 4.2 Discrete Fourier Transform
• Exercise
– Compute the 8-point DFT of
36
Trang 374 Discrete Fourier Transform
• 4.2 Discrete Fourier Transform
• Exercise
– Compute the 6-point DFT of
Trang 384.2 Properties of the DFT
Denote
• Periodicity, Linearity, Circular symmetries
• Time reversal
• Circular time shift
• Circular frequency shift
• Complex conjugate properties
• Circular correlation
• Parseval ’s theorem
38
Trang 394 Discrete Fourier Transform
Trang 414 Discrete Fourier Transform
4.2 Properties of the DFT
• Circular symmetries
– x’(n) is related to x(n) by a circular shift
Trang 454 Discrete Fourier Transform
Trang 464.3 Multiplication of two DFTs and Circular convolution
– Suppose we have two finite-duration sequences of length
N, x1(n) and x2(n) Their respective N-point DFTS are:
– Multiply two sequences, the result is a DFT, say X3(k), of a sequence x3(n)
46
Trang 474 Discrete Fourier Transform
4.3 Multiplication of two DFTs and Circular convolution
– The inverse of X3(k) is
– Circular convolution
Trang 484.3 Multiplication of two DFTs and Circular convolution
• Example
48
Trang 494 Discrete Fourier Transform
4.3 Multiplication of two DFTs and Circular convolution
• Solution
Trang 504.3 Multiplication of two DFTs and Circular convolution
50
Trang 514 Discrete Fourier Transform
4.3 Multiplication of two DFTs and Circular convolution
Trang 524.3 Multiplication of two DFTs and Circular convolution
52
Trang 534 Discrete Fourier Transform
4.3 Multiplication of two DFTs and Circular convolution
Trang 544.3 Multiplication of two DFTs and Circular convolution
• The circular of the two sequences yields the sequence
54
Trang 554 Discrete Fourier Transform
4.3 Multiplication of two DFTs and Circular convolution
• Example
• By means of DFT and IDFT, determine the sequence x3(n)
corresponding to the circular convolution of x1(n) and x2(n)
Trang 564.3 Multiplication of two DFTs and Circular convolution
• Solution
• Compute the DFTs
56
Trang 574 Discrete Fourier Transform
4.3 Multiplication of two DFTs and Circular convolution
Trang 584.3 Multiplication of two DFTs and Circular convolution
• Solution
• Thus,
58
Trang 594 Discrete Fourier Transform
4.3 Multiplication of two DFTs and Circular convolution
Trang 60• Exercise
60
Trang 615 Fast Fourier Transform
• The DFT can be written in the formula
Trang 62• Direct computation of DFT does not exploit the symmetry and periodicity properties
62
Trang 635 Fast Fourier Transform
5.1 Radix-2 FFT algorithms
• Split N-point data sequence x(n) into two N/2-point sequence
• The N-point DFT of x(n) can be expressed as
Trang 655 Fast Fourier Transform
5.1 Radix-2 FFT algorithms
• We can repeat the process for each of the sequences f1(n)
and f2(n)
Trang 665.1 Radix-2 FFT algorithms
• Thus
• The total number of multiplication is reduced again
66
Trang 675 Fast Fourier Transform
5.1 Radix-2 FFT algorithms
• Comparison of computational complexity
Trang 685.1 Radix-2 FFT algorithms
• Perform 8-point DFT in 3 stages
68
Trang 705.1 Radix-2 FFT algorithms
• Basic butterfly computation
70
Trang 715 Fast Fourier Transform
5.1 Radix-2 FFT algorithms
• Compute 8-point DFT (using FFT algorithm) of
x(n)={1, 2, 3, 0, 0, 0, 0, 0}
Trang 725.2 Radix-4 FFT algorithms
72