Introduction to Digital Signal Processing
Introduction to Digital Signal Processing
Chapter 1-NA
xử lí tín hiệu số (Trường Đại học Sư phạm Kỹ Thuật Thành phố Hồ Chí Minh)
CHAPTER 1:
INTRODUCTION
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Duration: 1 hr
Outline:
1. Signals
3. Why DSP?
* * * * * * *
Learning Digital Signal Processing is not something you accomplish; it’s a journey
you take.
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Signals ?
data, etc.
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
1-D signals
EEG
ECG
Speech signal
Color image
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
(indexed image)
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
ICA
ECA
Imagin
Plaque
g
location CCA
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
2.5-
2.5-D video signals
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
scanners, etc.
quality.
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
APPLICATIONS
Analyze Filter
Signal estimation Signal synthesis
Signal recognition
Signal modeling
Measures Processed
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Radar
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Biomedical
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Speech compression
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Speech
recognition
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Communication
mobile phone
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Image processing
David Lindsay
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Image processing
Image compression: reducing the redundancy in the
image data
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Image processing
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Music
Playback
Manipulation/mixing
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Noise removal
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Duration: 1 hr
Outline:
1. Signals
3. Why DSP?
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
original signal
"digital system“.
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Discrete-
Discrete-time signal vs.
continuous
continuous--time signal
Continuous-time signal:
- sound, voice…
Discrete-time signal:
second, …)
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
00 10 00 10 11
CT signal DT signal
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Performed by:
circuits (ASIC)
accelerators
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
{ 1 2 4 -9 5 3 }
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Duration: 1 hr
Outline:
1. Signals
3. Why DSP?
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
More reliable
assistant tools)
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
quantization error
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Duration: 2 hr
Outline:
1. A/D conversion
2. D/A conversion
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
ADC
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Sampling
Analog Digital
Sampling
world world
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Sampling
aliasing
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
antialiasing filter
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Anti-
Anti-aliasing filter
magnitude
frequency
0 W 2W =f 3W 4W
s
magnitude
spectrum
frequency
0 W 2W =f 3W 4W
s
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Audio CDs:
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
0 1 2 3 4
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Quantization
Quantization step
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
1.5V
1.1V
1.25V
1.0V
0.82V
0.5V
0.0V
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
voltages between
001 0.375 0.1875-0.5625
0 and 3V. The
to digital values.
111 2.625 2.4375-3.0
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Quantization parameters
Number of bits: N
N
Resolution: the gap between levels Q = R/2
N
Dynamic range = 20log(R/Q) = 20log(2 ) = 6.02N dB
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Error
Noise
Q/2
Q
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Non
Non--uniform quantization
variable
Output Non-
Q value is directly uniform
proportional to signal
amplitude SNR is
constant
Uniform
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
A s1 ( t ) 1
, 0 ≤ s1 ( t ) ≤
s 2 ( t ) = 1 + ln A
A
1 + ln( A s ( t ) )
s2(t)
1
1
, < s1 ( t ) ≤ 1
1 + ln A A 1.0
A=87.6
A=5
- 1.0 A=1
0 1.0 s1(t)
- 1.0
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
word
30-31 1111 15 31
62-63 1111 31 63
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
ITU G.711 A-
A-law curve
1.0
1
7/8
2
6/8
3
5/8
4
4/8
5
3/8
6
2/8 7
1/8
8
0
1/16 1/8 1/4 1/2 1.0
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Sign bit: 0
Part 1: 010
Part 2: 1110
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Duration: 2 hr
Outline:
1. A/D conversion
2. D/A conversion
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
DAC
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Anti-
Anti-imaging filter
magnitude
Anti-imaging
Images
filter
0 W 2W =f 4W = 2f
s s
frequency
Original two-sided
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
HW
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
HW
T ≤ Ts /(log2 L)
(b) When is the equality sign valid?
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
HW
n 0 1 2 3 4 5 6 7 8
Sample(V) 0.5715 4.9575 0.6250 3.6125 4.0500 0.9555 2.8755 1.5625 2.7500
sample.
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
HW
(c) How many bits per sample are needed for the number of levels
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
HW
(a) -3.45198 V
(b) 1.01119 V
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
HW
for the code 000 and a 5 V output for the code 111,
Draw the zero order hold output from the converter for
111 101 011 101 000 001 011 010 100 110
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Lecture #2
The concept of frequency in
CT & DT signals
Duration: 2 hrs
Outline:
1. CT sinusoidal signals
2. DT sinusoidal signals
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Mathematical description of CT
sinusoidal signals
Functions:
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
themselves different
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Properties of CT sinusoidal
signals (cont)
For f = 0 T = ∞
p
For f = ∞ T = 0
p
A j ( Ωt +θ ) A − j ( Ωt +θ )
xa (t ) = A cos(Ωt + θ ) = e + e
2 2
The frequency range for CT signal:
-∞ < f < +∞
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Mathematical description of DT
sinusoidal signals
1 .5
Plot:
0 .5
Amplitude
-0 . 5
-1
-1 . 5
-2
0 5 10 15 20 25 30 35 40
T i m e in d e x n
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
x(n + N) = x(n) ∀n
A cos[ 2πF0 ( n + N ) + θ)] = A cos( 2πF0 n + θ) ∀n
2πF0 N = 2πk
k
F0 =
N
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
obtained when:
Ω = π (or Ω = −π)
or, equivalently, 1 1
F= (or F = − )
2 2
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
0.5 0.5
0 0
-0.5 -0.5
-1 -1
0 5 10 15 20 25 30 0 5 10 15 20 25 30
F = 1/2
F = 3/6 0
FF =03/4= 3/4
1 1
0.5 0.5
0 0
-0.5 -0.5
-1 -1
0 5 10 15 20 25 30 0 5 10 15 20 25 30
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
x (t) x (nT)
a a
A cos(2πf nT + θ)
A cos( 2πf t + θ) 2 πf n
= A cos
+ θ
fS
f
F= Normalized
fs frequency
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
CT signals DT signals
ω = 2πf Ω = 2πF
− ∞ < ω < +∞
− ∞ < f < +∞
f − π ≤ Ω ≤ +π
F=
fs − 1 / 2 ≤ F ≤ +1 / 2
− π / T ≤ ω ≤ +π / T
− f s / 2 ≤ f ≤ +f s / 2 1
T=
fs
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
Exercise
x(t ) = 3 cos100π t , t[ s ]
a) Determine the minimum sampling rate required to avoid aliasing
b) Suppose that the signal is sampled at the rate fs = 200 Hz. What is
the DT signal obtained after sampling?
c) Suppose that the signal is sampled at the rate fs = 75 Hz. What is
the DT signal obtained after sampling?
d) What is the frequency 0 < f < fs/2 of a sinusoidal signal that yields
samples identical to those obtained in part (c)?
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
HW
Downloaded by T L (lethiloi040915@[Link])
lOMoARcPSD|13340603
HW
Downloaded by T L (lethiloi040915@[Link])