Lecture
Lecture
2
Course
Outcomes
S.NO Description Blooms’ Taxonomy
Level
• Overview
• Activities
Lectures - Course Notes/Literature
Homeworks/Exercise sessions
Paper Review Assignment
Exam
Why• study DSP ?
Analog Systems vs. Digital Systems
IN OUT IN 2 OUT
A/D +2 D/A
=4
• Main message:
Consumer electronics products have become
supercomputers (Mops/sec…Gops/sec), packed with DSP
functionalities…
DSP in applications : GSM
Cellular mobile telephony (e.g. GSM)
decision
x + x
1,0,1,1,0,…
a 1/a
noise
• While in practice…
DSP in applications : GSM
• DSP for digital communications (`physical layer’ ) :
!!
• In practice… .59,.41,.76,.05,.37,…
noise
?? 1,0,1,1,0,…
a
`Multipath’
Channel b
Δ
c
+
≈ Δ
Δ
Δ
Δ Δ
d
+
e
Δ Δ Δ Δ
DSP• GSM
in Channel
applications : GSM
Estimation/Compensation (continued)
• Multi-path channel is modeled with short (3…5 taps) FIR filter
a
b
Δ
IN[k] c OUT[k]
OUT[1] IN [1] 0 0 0 0
OUT[2] IN [ 2] IN [1] 0 0 0 a Δ Δ +
d
OUT[3] IN [3] IN [2] IN [1] 0 0 b Δ Δ Δ
OUT[4] IN [ 4] IN [3] IN [2] IN [1] 0 . c e
OUT[5] IN [5] Δ Δ Δ Δ
IN [4] IN [3] IN [2] IN [1] d
e
OUT[ N ] 0 0 0 0 IN [ N 4]
=convolution
DSP• GSM
in Channel
applications : GSM
Estimation/Compensation (continued)
2
OUT [1] TS [1] 0 0 0 0
OUT [ 2 ]
TS [ 2 ] TS [1] 0 0 0
a
OUT [ 4 ]
OUT [ 3]
TS
TS [ 3] TS [ 2 ] TS [1] 0 0
b
min OUT [ 5]
TS [[ 54]] TS [ 3] TS [ 2 ] TS [1] 0
d
. c
a, b, c, d , e TS [ 4 ] TS [ 3] TS [ 2 ] TS [1]
e
OUT [ N ] 0 0 0 0 TS [ N 4 ]
2
See PART-III on optimal filtering
DSP• GSM
in Channel
applications : GSM
Estimation/Compensation (continued)
= SPECTACULAR !!
DSP• GSM
in Channel
applications : GSM
Estimation/Compensation
= SPECTACULAR !!
DSP• GSM
in Channel
applications : GSM
Estimation/Compensation
X 1000
(0.3…1.5km)
c
se
FFT-rate = 4.3215 kHz
s/
op
(i.e. >4000 times 512-point FFTs per second !!!!)
0M
>2
• basic sampling rate is 2.21 MHz (=512*4.3215k)
8.84 MHz A/D or D/A (multi-rate structure)
c
se
s/
• fixed HP/LP/BP front-end filtering for frequency duplex
op
0M
See PART-I on filter banks/transmultiplexers
>6
• adjustable time-domain equalization filter (TEQ)
e.g. 32 taps @ 2.21 MHz
filter initialization via least-squares/eigenvalue procedure
• adaptive frequency-domain equalization filters (FEQ)
See PART-III on adaptive filters
= BOX FULL OF DSP/MATHEMATICS !!
DSP in applications : VDSL
VDSL specs
• 4096-point (I)FFT’s (or `similar’) for DMT-modulation
c
se
s/
FFT-rate = 4.3215 kHz
op
M
!!!!)
60
>1
(i.e. >4000 times 4096-point FFTs per second
• etc..
0
0 0.5 1 1.5 2 2.5 3
• FIR filters
• Linear-phase FIR filters
• FIR design by optimization
Weighted least-squares design, Minimax design
• FIR design in practice
`Windows’, Equiripple design, Software (Matlab,…)
• IIR filters
• Poles and Zeros
• IIR design by optimization
Weighted least-squares design, Minimax design
• IIR design in practice
Analog IIR design : Butterworth/Chebyshev/elliptic
Analog->digital : impulse invariant, bilinear transform,…
Software (Matlab)
FIR Filters
FIR filter = finite impulse response filter
• Impulse response
Linear Phase FIR Filters
e jx e jx 2. cos x
frequency response is
zero at zero
at zeroat0, 0
LP/HP/BP LP/BP BP HP
PS: IIR filters can NEVER have linear-phase property ! (proof see literature)
FIR Filter Design by Optimization
(I) Weighted Least Squares Design :
• Select one of the basic forms that yield linear phase
e.g. Type-1
• Optimization criterion is
W ( ) 0
where is a weighting function
FIR Filter Design by Optimization
• …This is equivalent to
Passband Ripple
1
0.8
Ad ( ) 1, P (pass - band) 0.6 Passband Cutoff -> <- Stopband Cutoff
0
0 0.5 1 1.5 2 2.5 3
i.e.
W ( ) 1, P W ( ) , S
1
xOPT Q . p Compare to p.10
+++ : simple
- - - : unpredictable behavior in between sample freqs.
FIR Filter Design by Optimization
• This is often supplemented with additional constraints, e.g. for pass-band and stop-band ripple
control : A( P ,i ) 1 P , for pass - band freqs. P ,1 , P , 2 ,... ( P is pass - band ripple)
A( S ,i ) S , for stop - band freqs. S ,1 , S , 2 ,... ( S is stop - band ripple)
AP .x bP
subject to (=pass-band constraints)
A .x b
S S
(=stop-band constraints)
= `Quadratic Programming’ problem
FIR Filter Design by Optimization
(II) `Minimax’ Design :
• Select one of the basic forms that yield linear phase
e.g. Type-1
• Optimization criterion is
W ( ) 0
where is a weighting function
• Leads to `Semi-Definite Programming’ (SDP) problem, for which efficient interior-
point algo’s & software are available.
FIR Filter Design by Optimization
• Conclusion:
(I) weighted least squares design
(II) minimax design
provide general `framework’, procedures to translate filter design problems into
standard optimization problems
• In practice (and in textbooks):
emphasis on specific (ad-hoc) procedures :
- filter design based on `windows’
- equiripple design
FIR Example
Filter: Low-pass
Design using `Windows’
filter design
• Ideal low-pass filter is
1 C
H d ( )
0 C
• Hence ideal time-domain impulse response is
• Time-domain multiplication with a window function w[k] corresponds to frequency domain convolution
with W(z) : h[k ] hd [k ].w[k ]
H ( z ) H d ( z ) *W ( z )
• Candidate windows : Han, Hamming, Blackman, Kaiser,…. (see textbooks, see DSP-I)
• Window choice/design = trade-off between side-lobe levels (define peak pass-/stop-band ripple) and
width main-lobe (defines transition bandwidth)
FIR Equiripple Design
• Starting point is minimax criterion, e.g.
• Based on theory of Chebyshev approximation and the `alternation theorem’, which (roughly) states that
the optimal d’s are such that the `max’ (maximum weighted approximation error) is obtained at L o+2
extremal frequencies…
46
REFERENCE VIDEOS
https://nptel.ac.in/courses/104/106/104106093/
https://www.coursera.org/lecture/materials-science/crystallograp
hy-and-the-electron-microscope-YhPAB
www.youtube.com/watch?v=mwlWIfNvt1w
www.youtube.com/watch?v=zdmEaXnB-5Q
https://www.youtube.com/watch?v=A56UtkYYExg
https://www.youtube.com/watch?v=x3LOTaEVLak
https://www.youtube.com/watch?v=1DHCV3LgITY
47
THANK YOU
For queries
Email: [email protected]