DIGITAL SIGNAL PROCESSING

(19A04704b)

IV B .TECH I SEM – R19

By

Mr.P.Srinivasulu
Assistant Professor
Department of ECE
VEMU INSTITUTE OF TECHNOLOGY

1
 UNIT-1
Introduction

2
 Signal Processing
• Humans are the most advanced signal processors
– speech and pattern recognition, speech synthesis,…

• We encounter many types of signals in various

applications
– Electrical signals: voltage, current, magnetic and electric fields,…
– Mechanical signals: velocity, force, displacement,…
– Acoustic signals: sound, vibration,…
– Other signals: pressure, temperature,…

• Most real-world signals are analog

– They are continuous in time and amplitude
– Convert to voltage or currents using sensors and transducers

• Analog circuits process these signals using

– Resistors, Capacitors, Inductors, Amplifiers,…

• Analog signal processing examples

– Audio processing in FM radios
– Video processing in traditional TV sets

3
Limitations of Analog Signal Processing
• Accuracy limitations due to
– Component tolerances
– Undesired nonlinearities

• Limited repeatability due to

– Tolerances
– Changes in environmental conditions
• Temperature
• Vibration
• Sensitivity to electrical noise
• Limited dynamic range for voltage and currents
• Inflexibility to changes
• Difficulty of implementing certain operations
– Nonlinear operations
– Time-varying operations

• Difficulty of storing information

4
 Digital Signal Processing
• Represent signals by a sequence of numbers
– Sampling or analog-to-digital conversions
• Perform processing on these numbers with a digital processor
– Digital signal processing
• Reconstruct analog signal from processed numbers
– Reconstruction or digital-to-analog conversion
digital digital
signal signal
analog analog
signal A/D DSP D/A signal

• Analog input – analog output

– Digital recording of music
• Analog input – digital output
– Touch tone phone dialing
• Digital input – analog output
– Text to speech
• Digital input – digital output
– Compression of a file on computer
5
 Pros and Cons of Digital Signal Processing
• Pros
– Accuracy can be controlled by choosing word length
– Repeatable
– Sensitivity to electrical noise is minimal
– Dynamic range can be controlled using floating point numbers
– Flexibility can be achieved with software implementations
– Non-linear and time-varying operations are easier to implement
– Digital storage is cheap
– Digital information can be encrypted for security
– Price/performance and reduced time-to-market

• Cons
– Sampling causes loss of information
– A/D and D/A requires mixed-signal hardware
– Limited speed of processors
– Quantization and round-off errors
6
Analog, digital, mixed signal
processing

7
Digital Signal Processing

8
Sampling and reconstruction

9
Sample and hold (S/H)circuit

10
A/D converter

11
A/D converter

12
Quantization noise

13
D/A convertion

14
D/A convertion

15
Reconstruction

16
Reconstruction

17
Reconstruction

18
Reconstruction

19
 Signals
• Continuous-time signals are functions of a real argument
x(t) where t can take any real value
x(t) may be 0 for a given range of values of t
• Discrete-time signals are functions of an argument that
takes values from a discrete set
x[n] where n  {...-3,-2,-1,0,1,2,3...}
Integer index n instead of time t for discrete-time systems
• x may be an array of values (multi channel signal)
• Values for x may be real or complex

20
Discrete-time Signals and Systems
• Continuous-time signals are defined over a
continuum of times and thus are represented by a
continuous independent variable.
• Discrete-time signals are defined at discrete times
and thus the independent variable has discrete
values.
• Analog signals are those for which both time and
amplitude are continuous.
• Digital signals are those for which both time and
amplitude are discrete.

21
 Analog vs. Digital
• The amplitude of an analog signal can take any real or complex value at each
time/sample

• The amplitude of a digital signal takes values from a

discrete set
1

-1
22
Periodic (Uniform) Sampling
• Sampling is a continuous to discrete-time conversion

-3 -2 -1 0 1 2 3 4

• Most common sampling is periodic

x n   x cnT     n  

• T is the sampling period in second

• fs = 1/T is the sampling frequency in Hz
• Sampling frequency in radian-per-second s=2fs rad/sec
• Use [.] for discrete-time and (.) for continuous time signals
• This is the ideal case not the practical but close enough
– In practice it is implement with an analog-to-digital converters
– We get digital signals that are quantized in amplitude and time

23
 Periodic Sampling
• Sampling is, in general, not reversible
• Given a sampled signal one could fit infinite continuous signals
through the samples
1

0.5

-0.5
-1
0 20 40 60 80 100

• Fundamental issue in digital signal processing

– If we loose information during sampling we cannot recover it
• Under certain conditions an analog signal can be sampled without
loss so that it can be reconstructed perfectly

24
 Representation of Sampling
• Mathematically convenient to represent in two stages
– Impulse train modulator
– Conversion of impulse train to a sequence

s(t)

Convert impulse
xc(t) x train to discrete- x[n]=xc(nT)
time sequence

xc(t) x[n]
s(t)

t n
-3T-2T-T 0 T 2T3T4T -3 -2 -1 0 1 2 3 4

25
 Unit Sample Sequence

[n] = 0, n  0
= 1, n = 0.

… 1
…

0 n

The unit sample sequence plays the same role for discrete-time sequences and
systems that the unit impulse (Dirac delta function) does for continuous-time
signals and systems.

26
 Impulse Function
The impulse function, also known as Dirac’s delta function, is used to
represented quantities that are highly localized in space. Examples include
point optical sources and electrical charges.

The impulse function can be visualized as a narrow spike having infinite

height and zero width, such that its area is equal to unity.

27
Definition of Impulse Function
The impulse function may be defined from its basic properties.

δ( x  x 0 )  0 , x  x0

x2
f ( x ) δ( x  x 0 ) dx  f ( x 0 ), x1  x 0  x 2

x1

Where f(x) is any complex-valued function of x. If f(x) is discontinuous at the

point x0, the value of f(x0) is taken as the average of the limiting values as x
approaches x0 from above and below.

This property is called the sifting property.

28
 Graphical Representation
On graphs we will represent (x-x0) as a spike of unit
height located at the point x0.

δ( x  x 0 )

x0

29
 Sampling Operation
The delta function samples the function f(x).

f(x0)

f ( x ) δ( x  x 0 )

x0

The function f(x) (x-x0) is graphed as a spike of height f(x0) located at the point x0.

30
 Unit Step Sequence
u[n] = 1, n  0
= 0, n < 0.

… 1
…

0 n

u [ n ]   [ n ]   [ n  1]   [ n  2 ] 


u [ n ]    [n  k ]
Conversely, the impulse sequence can be expressed
k 0
as the first backward difference of the unit step
n
sequence:
or u[ n ] 
  [ k ]
k  
 [ n ]  u [ n ]  u [ n  1]
31
 Exponential Sequence
x[n] =An

…
0 n
If we want an exponential sequence that is
zero for n < 0, we can write this as:

x[n ]  A  u[n ]
n

32
 Geometric Series
A one-sided exponential sequence of the form

 n , for n ≥ 0 and  an arbitrary constant

is called a geometric series. The series converges for |a| < 1, and its sum converges
to

1
 n

1  
n0


The sum of a finite number N of terms is
 N 1
1
N

 n

n0 1  

A general form can also be written:


N 2 1
  
N 2 N1

 
n

nN1 1 33
 Sinusoidal Sequence
x [ n ]  A cos(  o n   )

0 n
…

34
Sequence as a sum of scaled, delayed
impulses
a1
a-3

-4 -2 0 4 6 n
a2 a7

p [ n ]  a  3 δ[ n  3 ]  a 1 δ[ n  1]  a 2 δ[ n  2 ]  a 7 δ[ n  7 ]

35
 Sequence Operations
• The product and sum of two sequences are defined as the sample-by-
sample product and sum, respectively.
• Multiplication of a sequence by a number is defined as multiplication
of each sample value by this number.
• A sequence y[n] is said to be a delayed or shifted version of a sequence
x[n] if
y[n] = x[n – nd]
where nd is an integer.
• Combination of Basic Sequences
Ex 1 x[n] = K n, n  0,
= 0, n < 0,
or x[n] = K n u[n].

36
Systems

37
 Systems

T{}
x[n] y[n]

A discrete-time system is a transformation that maps an

input sequence x[n] into an output sequence y[n].

System Characteristics

1. Linear vs. non-linear

2. Causal vs. non-causal
3. Time invariant
38
 System Characteristics

T{}
x[n] y[n]

1. Linear vs. non-linear

2. Time invariant vs. time variant
3. Causal vs. non-causal
4. Stable vs. unstable
5. Memoryless vs. state-sensitive
6. Invertible vs. non-invertible
39
 Discrete-Time Systems
• Discrete-Time Sequence is a mathematical operation that maps a given
input sequence x[n] into an output sequence y[n]

y[n]  T { x[n]} x[n] T{.} y[n]

• Example Discrete-Time Systems

– Moving (Running) Average
y[n]  x[n]  x[n  1]  x[n  2 ]  x[n  3 ]

– Maximum

y[n]  max x[n], x[n  1], x[n  2 ]

– Ideal Delay System

y[n]  x[n  n o ]

40
 Linearity
A linear system is one that obeys the principle of
superposition,
T a1 x [ n ]  a 2 x [ n ]  a 1 y 1 [ n ]  a 2 y 2 [ n ]

where the output of a linear combination of inputs is the

same linear combination applied to the individual outputs.
This result means that a complicated system can be
decomposed into a linear combination of elementary
functions whose transformation is known, and then taking
the same linear combination of the results. Linearity also
implies that the behavior of the system is independent of
the magnitude of the input.
41
 Linear Systems
• Linear System: A system is linear if and only if
T { x 1 [n]  x 2 [n]}  T x 1 [n]  T x 2 [n] (additivit y)

and

T ax [n]   aT x[n] (scaling)

• Examples
– Ideal Delay System
y[n]  x[n  n o ]

T { x 1 [n]  x 2 [n ]}  x 1 [n  n o ]  x 2 [n  n o ]

T { x 2 [n ]}  T x 1 [n]  x 1 [n  n o ]  x 2 [n  n o ]
T ax [n]  ax 1 [n  n o ]
aT x[n]  ax 1 [n  n o ]

42
 Time (Shift) Invariance
A system is said to be shift invariant if the only effect caused by
a shift in the position of the input is an equal shift in the position
of the output, that is

T x [ n  n 0 ]   y [ n  n 0 ]

The magnitude and shape of the output are unchanged, only the
location of the output is changed.

43
 Time-Invariant Systems
• Time-Invariant (shift-invariant) Systems
– A time shift at the input causes corresponding time-shift at output
y[n]  T { x[n]}  y[n  n o ]  T x[n  n o]

• Example
– Square
y n   x[n  n ]
2
Delay the input the output is
y[n]  x[n] 
2 1 o
 
y n - n o   x[n  n  ] 
2
Delay the output gives o

• Counter Example
– Compressor System
Delay the input the output is y n   x[Mn  n ]
1 o
y[n]  x[Mn ]
Delay the output gives y n - n o   x M n  n o 

44
 Impulse Response
When the input to a system is a single impulse, the output is called
the impulse response. Let h[n] be the impulse response, given by

T  [ n ]   h [ n ]

A general sequence f [x] can be represented as a linear combination

of impulses, since


f ( x )  f ( x )   ( x )   
f ( u )  ( x  u ) du

 
f [ n ]  f [ n ]   [ n ]   f [ k ] [ n  k ]
k  

45
 Linear Shift-Invariant Systems
Suppose that T{} is a linear, shift-invariant system with h[n] as its
impulse response.
Then, using the principle
 of superposition,

     
T s [ n ]  T   s [ k ] [ n  k ]    s [ k ]T  [ n  k ]
 k    k  

and finally after invoking shift-invariance

 

T s [ n ]   s [ k ]T  [ n  k ]   s [ k ]h [ n  k]
k   k  

T s [ n ]  s [ n ]  h [ n ]

This very important result says that the output of any linear, shift-
invariant system is given by the convolution of the input with the
impulse response of the system. 46
 Causality
A system is causal if, for every choice of n0, the output
sequence at the index n = n0 depends only on the input
sequence values for n  0.

All physical time-based systems are causal because they are

unable to look into the future and anticipate a signal value
that will occur later.

47
 Causal System
• Causality
– A system is causal it’s output is a function of only the current and
previous samples

• Examples
– Backward Difference
y[n]  x[n]  x[n  1]

• Counter Example
– Forward Difference
y[n]  x[n  1]  x[n]

48
 Stability
A system is stable in the bounded-input, bounded-output
(BIBO) sense if and only it every bounded input produces a
bounded output sequence.
The input x[n] is bounded if there exists a fixed positive finite
value Bx such that

x[n ]  B x   for all n

Stability requires that for any possible input sequence there exist
a fixed positive value By such that

y [ n ]  B y  

49
 Stable System
• Stability (in the sense of bounded-input bounded-output BIBO)
– A system is stable if and only if every bounded input produces a bounded
output

x[n]  B x    y[n]  B y  

• Example
– Square
y[n]  x[n] 
2

if input is bounded by x[n]  B x  



output is bounded by y[n]  B 2x  

• Counter Example
– Log
y[n]  log 10  x[n] 

even if input is bounded by x[n]  B x  

output not bounded for x n   0  y 0   log 1

 x n    

50
 Memory (State)

A system is referred to as memoryless if the output y[n] at every

value of n depends only on the input x[n] at the same value of n.

If the system has no memory, it is called a static system. Otherwise

it is a dynamic system.

51
 Memoryless System
• Memoryless System
– A system is memoryless if the output y[n] at every value of n depends
only on the input x[n] at the same value of n

• Example Memoryless Systems

– Square
x[n] 
2
y[n] 

Sign
– y[n]  sign x[n]

• Counter Example
– Ideal Delay System
y[n]  x[n  n o ]

52
 Invertible System
A system is invertible if for each output sequence we can find a
unique input sequence. If two or more input sequences produce
the same output sequence, the system is not invertible.

53
 Passive and Lossless Systems
A system is said to be passive if, for every finite energy input
sequence x[n], the output sequence has at most the same energy:

   
 x [ n ]
2


2
y[n ]   
n   n  



If the above inequality is satisfied with an equal sign for every

input signal, the system is said to be lossless.

54
 Examples of Systems

Ideal Delay System y [n ]  x[n  n d ]

 x [ n
2
1
Moving Average System y [ n ]   k]
M 2  M1 1 kM
1

Memoryless non-linear System y [ n ]  x [ n ]2

 n

Accumulator System y [ n ]   x [ k ]
k  

where M is a
Compressor System y [ n ]  x [ Mn ] positive integer


Backward Difference System y [ n ]  x [ n ]  x [ n  1]

55
 Impulse Response of LTI Systems

Find the impulse response by computing the response to [n].

Systems whose impulse responses have only a finite number of

nonzero samples are called finite-duration impulse response
(FIR) systems.

Systems whose impulse responses are of infinite duration are

called infinite-duration impulse response (IIR) systems.

If h[n] = 0 for n < 0, the system is causal.

56
 Impulse Response for Examples
Find the impulse response by computing the response to [n]


Ideal Delay System y [ n ]  δ[ n  n d ] FIR

 1 
 , -M 1  n  M 2

Moving Average System y [ n ]   M  M 1 FIR

2 1
 0, otherwise

n   1, n  0
Accumulator System y [ n ]   δ[ k ]   IIR
k   0, n  0

y[n ]  u[n ]

Backward Difference System y [ n ]  δ[ n ]  δ[ n  1] FIR

57
 Stability Condition for LTI Systems
An LTI system is BIBO stable if and only if its impulse response
is absolutely summable, that is
 
S   h [ k ]  
k  

58
 Stable and Causal LTI Systems
• An LTI system is (BIBO) stable if and only if
– Impulse response is absolute summable


 h k   
k  

– Let’s write the output of the system as

 

y n    hk x n  k    h k  x n  k 
k   k  

– If the input is bounded

x[n]  B x
– Then the output is bounded by  
y n   B x  h k 
k  

– The output is bounded if the absolute sum is finite
• An LTI system is causal if and only if
h k   0 for k  0

59
 Difference Equations
An important subclass of LTI systems are defined by
an Nth-order linear constant-coefficient difference equation:
N M

 ak y[n  k ]   bm x [ n  m ]
0 0

Often the leading coefficient a0 = 1. Then the output y[n] can be

computed recursively from
N M

y [ n ]    a k y [ n  k ]   b m x[n  m ]
k 1 m 0

A causal LTI system of this form can be simulated in

MATLAB using the function filter
y = filter(a,b,x);
60
Accumulator Example
n
Accumulator System y [ n ]   x [ k ]
k  

n 
n 1

y [ n ]   x [ k ]  x [ n ]   x [ k ]  x [ n ]  y [ n  1];
k   k  

b0  1 x[n] + y[n]
a0  1
one sample
a1  1 delay

positive feedback system

61
 Total Solution Calculation
N M

 a k y [ n  k ]   b m x[n  m ]
k 0 m0

The output sequence y[n] consists of a homogeneous solution yh[n] and a

particular solution yp[n].
y [n ]  y h[n ]  y p [n ]

where the homogenous solution yh[n] is obtained from the homogeneous equation:

 a k y h [n  k ]  0
k0

62
 Homogeneous Solution
N
Given the homogeneous equation:
 a k y h [n  k ]  0
k0

Assume that the homogeneous solution is of the form

yh [n ]  
n


then 
N

y h [ n ]   a k

nk
 
n N
a 0 
N
 a 1
N 1
  a N 0
k 0

defines an Nth order characteristic polynomial with roots 1, 2 … N

The general solution is then a sequence yh[n]

N

 Am m
n
yh [n ] 
1

(if the roots are all distinct) The coefficients Am may be found from the
initial conditions. 63
 Particular Solution
The particular solution is required to satisfy the difference equation for a specific
input signal x[n], n ≥ 0.
N M

 a k y [ n  k ]   b m x[n  m ]
k 0 m 0

To find the particular solution we assume for the solution yp[n] a form that depends
on the form of the specific input signal x[n].

y [n ]  y h[n ]  y p [n ]

64
General Form of Particular Solution

Input Signal Particular Solution

x[n] yp[n]
A (constant) K
AMn KMn
AnM K0nM+K1nM-1+…+KM
AnnM An(K0nM+K1nM-1+…+KM)
 A cos(  o n )  K cos(  n )  K sin(  n )
  1 o 2 o

 A sin(  o n ) 

65
 Example (66
Determine the homogeneous solution for

y [ n ]  y [ n  1]  6 y [ n  2 ]  0

yh[n ]  
n
Substitute

  
n n 1
 6
n2
 
n2
 2
   6 0


n2
  3   2  0

Homogeneous solution is then

 3  2 
n n
y [n ]  A   A   A
n n
 A
h 1 1 2 2 1 2

66
 Example (67
Determine the particular solution for

y [ n ]  y [ n  1]  6 y [ n  2 ]  x[ n ]

with x[ n ]  8 u [ n ] and y[-1] = 1 and y[-2] = -1

The particular solution has the form y p [ n ]  




     6  8

which is satisfied by  = -2

67
 Example (68
Determine the total solution for

y [ n ]  y [ n  1]  6 y [ n  2 ]  x[ n ]

with x[ n ]  8 u [ n ] and y[-1] = 1 and y[-2] = -1

The total solution has the form

y [ n ]  y h [ n ]  y p [ n ]  A 1  3   A 2  2   2
n n

  
then 1 1
y [  1]   A 1  A2  2  1 A   1 .8
1
3 2
1 1 A 2  4 .8
y[  3]  A1  A2  2  1
9 4

y [ n ]   1 .8   3   4 .8  2   2
n n

68
 Initial-Rest Conditions

The output for a given input is not uniquely specified.

Auxiliary information or conditions are required.

Linearity, time invariance, and causality of the system will

depend on the auxiliary conditions. If an additional condition is
that the system is initially at rest (called initial-rest conditions),
then the system will be linear, time invariant, and causal.

69
 Zero-input, Zero-state Response
An alternate approach for determining the total solution y[n] of a difference equation
is by computing its zero-input response yzi[n], and zero-state response yzs[n]. Then
the total solution is given by y[n] = yzi[n] + yzs[n].

The zero-input response is obtained by setting the input x[n] = 0 and satisfying the
initial conditions. The zero-state response is obtained by applying the specified input
with all initial conditions set to zero.

70
 Example (1/2)
y [ n ]  y [ n  1]  6 y [ n  2 ]  0 y [  1]  1
y[ 2]   1
y h [ n ]  A 1  1  A 2 2  A1  3   A 2  2 
n n n n

Zero-input response: 
27
A     5 .4
y [0 ]  A  A   y [  1]  6 y [  2 ]   1  6   7 1
5
zi 1 2

y zi [1]  A 1 (  3 )  A 2 ( 2 )   y [ 0 ]  6 y [  1]  7  6  13
8
A 2     1 .6
5

Zero-state response: y [ n ]  y [ n  1]  6 y [ n  2 ]  x[ n ]

with x[ n ]  8 u [ n ]
18
y [ 0 ]  A  A  2  x[ 0 ]  8 A1  5  3 .6
zs 1 2

y zs [1]  A 1 (  3 )  A 2 ( 2 )  2  x[1]  y [ 0 ]  8  8  0 32
71
A2   6 .4
5
 Mitra Example (2/2)
Zero-input response:

yzi [ n ]   5 .4 (  3 ) n  1 .6 ( 2 ) n

Zero-state response:

yzs [ n ]  3 .6 (  3 ) n  6 .4 ( 2 ) n  2

Total solution is

y [ n ]  y z [ n ]  y zs [ n ]

y [ n ]  1 .8 (  3 ) n  4 .8 ( 2 ) n  2

This is the same as before

72
 Impulse Response

The impulse response h[n] of a causal system is the output

observed with input x[n] = [n].
For such a system, x[n] = 0 for n >0, so the particular solution is
zero, yp[n]=0. Thus the impulse response can be generated from the
homogeneous solution by determining the coefficients Am to satisfy
the zero initial conditions (for a causal system).

73
 Example
Determine the impulse response for

y [ n ]  y [ n  1]  6 y [ n  2 ]  x[ n ]

The impulse response is obtained from the homogenous solution:

h [ n ]  A 1  3   A 2  2 
n n

For n=0
y [ 0 ]  y [  1]  6 y [  2 ]  x[ 0 ]

h[ 0 ]  δ[ 0 ]  1 A1  A 2  1
For n=1
y [1]  y [ 0 ]  6 y [  1]  x[1]
 3 A 1  2 A2  A1  A2   0

h[1]  h[ 0 ]  δ[1]  0 3 2
A 1  , A 2 
5 5

3 2
 3  2  n  0
n n
h [ n ]  
5 5 74
 DSP Applications
• Image Processing • Military
– Pattern recognition – Secure communication
– Robotic vision – Radar processing
– Image enhancement – Sonar processing
– Facsimile – Missile guidance
– Satellite weather map • Telecommunications
– Animation – Echo cancellation
• Instrumentation/Control – Adaptive equalization
– Spectrum analysis – ADPCM transcoders
– Position and rate control – Spread spectrum
– Noise reduction – Video conferencing
– Data compression – Data communication
• Speech/audio • Biomedical
– Speech recognition/synthesis – Patient monitoring
– Text to speech – Scanners
– Digital audio – EEG brain mappers
– equalization – ECG analysis
– X-ray storage/enhancement
75
76
Image enhancement

77
 More Examples of Applications
• Sound Recording • Telephone Dialing
Applications Applications
– Compressors and Limiters
• FM Stereo Applications
– Expanders and Noise Gates
– Equalizers and Filters • Electronic Music
– Noise Reduction Systems Synthesis
– Delay and Reverberation – Subtractive Synthesis
Systems – Additive Synthesis
– Special Effect Circuits
• Echo Cancellation in
• Speech Processing
Telephone Networks
– Speech Recognition
– Speech Communication • Interference Cancellation
in Wireless
Telecommunication

78
 Reasons for Using DSP
• Signals and data of many types are • Flexibility in reconfiguring
increasingly stored in digital • Better control of accuracy
computers, and transmitted in requirement
digital form from place to place. In
• Easily transportable and possible
many cases it makes sense to offline processing
process them digitally as well.
• Cheaper hardware in some case
• Digital processing is inherently
stable and reliable. It also offers • In many case DSP is used to
certain technical possibilities not process a number of signals
available with analog methods. simultaneously. This may be done
by interlacing samples (technique
• Rapid advances in IC design and known as TDM) obtained from the
manufacture are producing ever
various signal channels.
more powerful DSP devices at
decreasing cost.
• Limitation in speed & Requirement
in larger bandwidth

79
80
System Analysis

81
Frequency Response

82
Ideal lowpass filter-example

83
Non causal to causal

84
Phase distortion and delay

85
Phase delay

86
Phase delay

87
Group delay

88
Group delay

89
The Z-Transform

90
 Z-Transform Definition
-Counterpart of the Laplace transform for discrete-time signals
-Generalization of the Fourier Transform
-Fourier Transform does not exist for all signals

The z-transform of a sequence x[n] is defined as



Z { x [ n ]}   x [ n ] z  n  X (z)
n  
The inverse z-transform is given by
1
x [ n ]  X ( z ) z n 1 dz
2 j

C

This expression denotes a closed contour integral taken counterclockwise

about the origin and in the region of convergence (ROC). It incorporates the
Cauchy integral theorem from the theory of complex variables. This result is
not given directly in Oppenheim, but may be found in Proakis.
91
 Relationship to Fourier Transform
The z-transform of a sequence x[n] is defined as

 

X ( z )   x [ n ] z  n
n  

The Fourier transform of a sequence x[n] is defined as



X e j    x [ n ]e  j n
n  

j
For z  e the z-transform reduces to the Fourier transform

This domain is a circle of unit radius in the complex plane, that is

|z| = 1.
92
 Convergence of
• DTFT does not always converge
the z-Transform


    x n  e
X e j  j n

n  

– Infinite sum not always finite if x[n] no absolute summable

– Example: x[n] = anu[n] for |a|>1 does not have a DTFT

• Complex variable z can be written as r ej so the z-transform



 j
     j
 n 
   n
 j n

X re  
n  
x n re  
n  
x n r e

• DTFT of x[n] multiplied with exponential sequence r -n

– For certain choices of r the sum maybe made finite



 x n  r - n  
n  

93
 Unit Circle in Complex Z-Plane
-The z-transform is a function of the complex z variable
-Convenient to describe on the complex z-plane
-If we plot z=ej for =0 to 2 we get the unit circle

Im
 j 
ze

Unit Circle


Re
1

z-plane

94
 Region of Convergence (ROC)
For any given sequence, the set of value of z for which the z-transform converges
is called the region of convergence, ROC. The criterion for convergence is that
the z-transform be absolutely summable:


 n

 
x[n ] z  

If some value of z, say, z = z1, is in the ROC, then all values of z on the circle defined
by |z| = |z1| will also be in the ROC. So, the ROC will consist of a ring in the z-plane
centered about the origin. Its outer boundary will be a circle (or the ROC may extend
outward to infinity), and its inner boundary will be a circle (or it may extend inward
to include the origin).

95
 Region of Convergence
Im

r1
r2

Re

z-plane

The region of convergence (ROC) as a ring in the z-plane. For specific cases, the inner boundary
can extend inward to the origin, and the ROC becomes a disk. In other cases, the outer boundary
can extend outward to infinity. 96
 Laurent Series

 




X ( z )   x [ n ] z  n
n  

A power series of this form is a Laurent series. A Laurent series, and therefore
a z-transform, represents an analytic function at every point within the region
of convergence. This means that the z-transform and all its derivatives must be
continuous functions of z within the region of convergence.
P(z)
X ( z ) 
Q (z)

Among the most useful z-transforms are those for which X(z) is a rational
function inside the region of convergence, for example where P(z) and Q(z)
are polynomials. The values of z for which X(z) are zero are the zeros of X(z)
and the values for which X(z) is infinite are the poles of X(z).
97
 Properties of the ROC
1. The ROC is a ring or disk in the z-plane centered at the origin
0  rR  z  rL  

2. The Fourier transform of x[n] converges absolutely if and only if the ROC of the z-
transform of x[n] includes the unit circle.
3. The ROC cannot contain any poles.
4. If x[n] is a finite-duration sequence, then the ROC is the entire z-plane except possibly
z=0 or z=.
5. If x[n] is a right-handed sequence (zero for n < N1 < ), the ROC extends outward from
the outermost (largest magnitude) finite pole in X(z) to (and possibly including infinity).
6. If x[n] is a left-handed sequence (zero for n > N2 > - ), the ROC extends inward from
the innermost (smallest magnitude) nonzero pole in X(z) to (and possibly including) zero.
7. A two-sided sequence is an infinite-duration sequence that is neither right-sided or left-
sided. If x[n] is a two-sided sequence, the ROC will consist of a ring in the z-plane,
bounded on the interior and exterior by a pole and , consistent with property 3, not
containing any poles.
8. The ROC must be a connected region.
98
Properties of ROC Shown Graphically
Finite-Duration Signals

Causal Entire z-plane

Except z = 0

Anticausal Entire z-plane

Except z = infinity

Two-sided Entire z-plane

Except z = 0 and z = 

Infinite-Duration Signals

Causal
|z| > r2

Anticausal |z| < r1

Two-sided
r2 < |z| < r1
99
 Example: Right-Sided Sequence

x[n ]  a nu[n ]


X ( z )   x [ n ] z  n
 n  

 n 1
X (z)   az  1  
1  az 1
n0

ROC az 1
1

or za

100
Example: Left-Sided Sequence

x [ n ]   a n u [  n  1]

nonzero for n  1
 

X ( z )   x [ n ] z  n
n  

1 

 az   1   a 1 z 
n n
X ( z )   1

n   n 0

1 1
 1  
1 1
1 a z 1  az

ROC a 1 z  1

or za
101
Example: Sum of Two Exponential
Sequences
n n
 1   1 
x [ n ]    u [ n ]     u [ n ]
 2   3 


X (z)   x [ n ] z  n
n  

1 1
X ( z )  
1 1
1 1
2
z 1 1
3
z

ROC 1
2
z 1 and  13 z  1

z 1
2

2  1
z 1  1 1 2 zz  
X ( z ) 
1
z

1  z 1 1  z 1 
3 2 12
1 1
 z  1  z  1 
2 3 2 3

Poles at ½ and -1/3, zeros at 0 and 1/12 102

Example: Two-Sided Sequence
 

 1   1 
n

x [ n ]     u [ n ]    u [  n  1]
 3   2 

 1 n 1 1
   u [ n ]  1 z 
 3  1  1
3
z 3
n
1 1 1
   u [  n  1]  z 
2 1  12 z 1 2
2 zz  1

X (z)  12

z  1  z  1

2 3

1 1
ROC  z 
3 2
103
 Example: Finite Length Sequence

a n 0  n  N 1
x [ n ]  
otherwise
 0


X (z)   x [ n ] z  n
n  

1  az

n 1 N
N 1

 az  

1
X (z)  1

n 0 1  az

The sum is finite, so

ROC a   and z  0
Pole-zero plot for N = 16

104
 Z-transforms with the same pole-zero locations
illustrating the different possibilities for the ROC.
Each ROC corresponds to a different sequence.

b) A right-sided sequence c) A left-sided sequence

e) A two-sided sequence
d) A two-sided sequence 105
 Common Z-Transform Pairs
Sequence Transform ROC
 [n ] 1 All z
Sequence Transform ROC
1 1  [cos  ] z  1
u[n ] z1 [cos  n ]u [ n ]
1  z 1 0 0
z1
1  [ 2 cos  ] z 1  z  2
0
1
 u [  n  1] z1 [sin  ] z  1
1 z 1 [sin  n ]u [ n ] 0
z1
1  [ 2 cos  ] z 1  z  2
0

m 0
 [n  m ] z * 1  [ r cos  ] z  1
[cos  n ]u [ n ] 0
z r
1
1  [ 2 r cos  ] z 1
 r 2 z 2
0
n
a u[n ] za 0
1  az 1
[ r sin  ] z  1
1 [sin  n ]u [ n ] 0
z r
1  [ 2 r cos  ] z 1  r 2 z  2
0
za
 a n u [  n  1] 1  az 1 0

az
1 an, 0  n  N  1 1  a N zN
na n u [ n ] z 0
za 1  az 1
1  az 1  0, otherwise
2

1
az
 na n u [  n  1] za
1  az 
2
1

106
*All z except 0 (if m > 0) or  (if m<0)
 Z-Transform Properties (1/2)
Linearity

ax 1 [ n ]  bx 2 [ n ] Z  aX 1 ( z )  bX 2 (z) ROC contains Rx  Rx

1 2

Time Shifting

x [ n  n 0 ] Z  z
 n0
X (z) ROC = Rx (except for possible addition
or deletion of 0 or )

Multiplication by an Exponential Sequence

n Z  z 
z x [ n ]   X   ROC  z 0 R x
0
 z 0 

Differentiation of X(z)

Z dX ( z )
nx [ n ]    z ROC  R x
dz 107
 Z-Transform Properties (2/2
Conjugation of a Complex Sequence

x  [ n ] Z X *
( z* ) ROC Rx

Time Reversal 

1  * 1
x [  n ]   X  * 
* Z ROC 
 z  Rx

Convolution of Sequences

x 1 [ n ]  x 2 [ n ] Z  X 1 ( z ) X 2 ( z ) ROC contains Rx  R

1 x2

Initial-Value Theorem

x [ 0 ]  lim X ( z ) provided that x[n] is zero for n < 0, i.e. that x[n] is causal.
z   108
 Inverse z-Transform
Method of inspection – recognize certain transform pairs.

Partial Fraction Expansion

M
k  1  c
M
z
1 
 bk z b0 k 1
k

X (z)  k 0 Factor to X (z) 

N
a0  1  d
N

z  1 
 a k z
k
 k
k 1
k 0





A k  1  d k z  1  X ( z )
N
X (z)   Ak
1
where
1 d z zdk
k 0 k

109
 Example: Second-Order Z-Transform
1
X ( z )  1
1  1
4
z 1
1  1
2
z 1
 z 
 2




A1 A2
X (z)   1
1
1 z
1 1 12
4
z

1
A   1
 14   1
1
1 1
2

1
A  2
 12   1
2
1 1
4

1 2
X (z)  
1 1
1 1
4
z 1 1
2
z
110
 Partial Fraction Expansion
M N N
Ak
If M > N X (z)   B z r
r
 
1 d z
1

r0 k 1 k

Br can be obtained by long division of numerator by denominator, stopping

when the remainder is of lower degree than the denominator.

If M > N and X(z) has multiple-order poles, specifically a pole of order s at z=di

M N
X ( z )   B z  r  
N s
Ak Cm
 
1  d 
r 1 m
r0 k 1,k  i 1  dkz m 1 z 1
i

1  d sm  
s
 1

C m  sm  sm 1  d iw X w 
( s  m )! (  d i )  dw  wd 1

111
 Example
X (z)  1  1 2
1
 1  2 z 
z
2
z 1
z 
1  1
z 1 1  z 1  1 3
z 1  1
2
z
2
2 2

2
A1 A2 B  z   2
X ( z )  B 0  1
 0
1 z 2 
1 1
2 z 1 z 1 2

1  ( 1 ) 1 
2
3 
2
A A
X (z)  B  1
 2 A  2
  9
1
 12 
1
1 1  z 1
1
0 1
2
z 1 1
 9 8 (1  1) 2
X ( z )  2  1
 A   8
1 2 z 1
1 z 1
1 
1
2
1 1
2

 12 
n
x [ n ]  2 [ n ]  u[n ]  8u[n ]
112
 Power Series Expansion
 

X ( z )   x [ n ] z  n
n  

Note that  [ n  m ] Z  z

m

Example :

X ( z )  z 2 1  z  1 1  z  2   z 2 
1
1
2
1
2
z 1 1
2
z

x[n ]  2 [n  2 ]  1
2
 [ n  1]   [ n ]  1
 [ n  1]
2

113
 Example
X ( z )  log 1  az 1
 z a

Expand in power series:

 1
 
(  1)
n 1 n
a z
n

log 1  az  
n 1 n

 (  1) n 1 a n , 
n 1
x[n ]   n
 0, n0

114
 Contour Integration
Cauchy integral theorem

 1, k 1
z
k
1
2 j
dz  
C
0. k 1

C is a counterclockwise contour that encircles the origin.

Then one can show that 

1
x [ n ]  X ( z ) z n 1 dz
2 j

C

x[n ]   residues of X(z)z n 1

at the poles inside C 

115

Residue Calculation
If X(z) is a rational function of z, we can write

X  z z n 1   (z)
z  d0
s

Then one can show that

X 
1
Res (z)z n 1
at z  d 0   d 
s  1 !
0

116
 Quick Review of LTI Systems
• LTI Systems are uniquely d etermined by their impulse response
y n    x k  h n  k   x k   h k 
k  
• We can write the input-output relation also in the z-domain

Y z   H z X z 

• Or we can define an LTI system with its frequency response

   H e X e 
Y e j j j

• H(ej) defines magnitude and phase change at each frequency

• We can define a magnitude response

Y e j     
 H e j X e j
• And a phase response

Y e  j
 
H e j
 X e  
j

117
 Phase Distortion and Delay
• Remember the ideal delay system
hi n    n  n d   
DTFT
 H i e  
j e
 j n d

• In terms of magnitude and phase response


H id e  j
 1


 H id e j     n d
  



• Delay distortion is generally acceptable form of distortion

– Translates into a simple delay in time
• Also called a linear phase response
– Generally used as target phase response in system design
• Ideal lowpass or highpass filters have zero phase response
– Not implementable in practice
118
 System Functions for Difference Equations
• Ideal systems are conceptually useful but not implementable
• Constant-coefficient difference equations are
– general to represent most useful systems
– Implementable
– LTI and causal with zero initial conditions
N M

 a y nk
 k    b x n k
 k 
k 0 k 0

• The z-transform is useful in analyzing difference equations

• Let’s take the z-transform of both sides

N M

 a z k
Y z    b z  k X z 
 k k
k0 k0

 
   
 
N k M k

  ak z Y z    b k z X z
 k 0   k 0 

119
 System Function
• Systems described as difference equations have system functions of the
form
 1  
M
  
M
Y z  bk
z
k b c z 1
k
H z  
X z    0 
kN
0 k 1
 N
a
 a k z
k  0 
 1  d k
z 1 
k0 k 1

• Example
H z   1  z
1 2
 
1  2 z 1  z  2

Y z 
 3  1  1  1 z  1  3 z  2 X z 


 1 1 
z 1  z
1   
  




 2  4  4 8


 1 
1 1  3 2
z  z

Y z   1  
2 z  1  z  2 X z 
 4 8 


y n   y n  1   y n  2   x n   2 x n  1   x n  2 
1 3
4 8

120
 Stability and Causality
• A system function does not uniquely specify a system
– Need to know the ROC
• Properties of system gives clues about the ROC
• Causal systems must be right sided
– ROC is outside the outermost pole
• Stable system requires absolute summable impulse response


 h n   
k  

– Absolute summability implies existence of DTFT

– DTFT exists if unit circle is in the ROC
– Therefore, stability implies that the ROC includes the unit circle
• Causal AND stable systems have all poles inside unit circle
– Causal hence the ROC is outside outermost pole
– Stable hence unit circle included in ROC
– This means outermost pole is inside unit circle
– Hence all poles are inside unit circle
121
 Example
• Let’s consider the following LTI system

y n   y n  1   y n  2   x n 
5
2
• System function can be written as





H z  
1 1 
1


1  z 
1  2 z  1 
 2 


• Three possibilities for ROC

– If causal ROC1 but not stable ROC : z  2
1
– If stable ROC2 but not causal
1
– If not causal neither stable ROC3 ROC :  z  2
2
2
1
ROC 3
: z 
2
122
 Structures for Discrete-Time Systems

• Block Diagram Representation of Linear Constant-Coefficient

Difference Equations
• Signal Flow Graph Representation of Linear Constant-Coefficient
Difference Equations
• Basic Structures for IIR Systems
• Transposed Forms
• Basic Network Structures for FIR Systems
• Lattice Structures

123
 Introduction

• Example: The system function of a discrete-time system is

b  b z 1
H ( z )  0 1
, |z| > |a|
1
1  az

• Its impulse response will be

h[n] = b0anu[n] + b1an-1u[n-1]
• Its difference equation will be
y[n] –ay[n-1] = b0x[n] + b1x[n-1]
Since this system has an infinite-duration impulse response, it is not possible
to implement the system by discrete convolution. However, it can be rewritten
in a form that provides the basis for an algorithm for recursive computation.
y[n] = ay[n-1] + b0x[n] + b1x[n-1]

124
Block Diagram Representation of Linear Constant-coefficient
Difference Equations

x[n] a ax[n] Multiplication of a

sequence by a constant

z-1 Unit delay

x[n] x[n-1]

+
Addition of two sequences
x1[n] x1[n] + x2[n]

x2[n]

125
 Block Diagram Representation of Linear Constant-coefficient
Difference Equations


M
bk z  k
Y (z)
H (z)  k 0
 

N X (z)
z k
1  a
k 1
k

 
1 
 1 
  M k 
H (z)  H 2 ( z)H 1( z)   b z 
   k 0 k
N

  1   a k
z k 
  k 1 
V ( z)  H 1( z) X (z)
Y (z)  H 2 ( z )V ( z )  
2 
M k   1 

H ( z )  H 1 ( z ) H 2 ( z )     bk z  
 k 0   1  N


k
 akz 

   k 1

W ( z )  H 2 ( z ) X ( z )

Y ( z )  H 1 ( z )W ( z )

126
 Block diagram representation for a general Nth-order
difference equation:
Direct Form 127



x[n] b0 v[n] y[n]
+ +
z-1 z-1
b1 a1
x[n-1]
+ + y[n-1]

z-1 z-1
x[n-2] y[n-2]
bM-1 aN-1
z-1 + +
x[n-M] z-1
bM aN
y[n-N]

127
 Block diagram representation for a general Nth-order
difference equation:
Direct Form 128

x[n] w[n] b0 y[n]

+ +
z-1
a1 b1
+ w[n-1] w[n-1] +
z-1 z-1
w[n-2] w[n-2]
aN-1 bM-1
+ +
z-1 z-1
aN bM
w[n-N] w[n-M]

128
Combination of delay units (in case N = M)

x[n] w[n] b0 y[n]

+ +
z-1
a1 b1
+ +
z-1

aN-1 bN-1
+ +
z-1
aN bN

129
 Block Diagram Representation of Linear Constant-
coefficient Difference Equations 2
• An implementation with the minimum number of delay elements is
commonly referred to as a canonic form implementation.
• The direct form I is a direct realization of the difference equation
satisfied by the input x[n] and the output y[n], which in turn can be
written directly from the system function by inspection.
• The direct form II or canonic direct form is an rearrangement of the
direct form I in order to combine the delay units together.

130
Signal Flow Graph Representation of Linear Constant-
coefficient Difference Equations

a
Attenuator
x[n] d e y[n]
a z-1
z-1 c
Delay Unit b
z-1

Node: Adder, Separator, Source, or Sink

131
 Basic Structures for IIR Systems

• Direct Forms
• Cascade Form
• Parallel Form
• Feedback in IIR Systems

132
 Basic Structures for IIR Systems

• Direct Forms

N M

y [ n ]   a k y [ n  k ]   b k x[n  k ]
k 1 k 0


k
bk z
k 0
H (z)  N

 a
k
1 k
z
k 1

133
Direct Form I (M = N)

134
Direct Form II (M = N)

135
Direct Form II

136
x[n] y[n]
+ +
z-1 z-1
2 0.75
+ +
z-1 z-1
-0.125

x[n] y[n]
+ +
z-1
0.75 2
+ +
1  2 z 1  z 2
z-1 H ( z ) 
-0.125 1  0 .75 z 1  0 .125 z  2

137
x[n] y[n]

z-1 2 0.75 z-1

z-1 z-1
-0.125

x[n] y[n]

z-1
0.75 2

z-1
-0.125

138
 Basic Structures for IIR Systems 2

• Cascade Form
(1  g z 1 )
M 1 M 2

(1  h k z 1 )( 1  h kz  1 )
 k 
H (z)  A k 1 k 1


N1
(1  c z 1 ) N 2

(1  d z 1 )( 1  d  z  1 )
k 1
k  k 1
k k

where M = M1+2M2 and N = N1+2N2 .

• A modular structure that is advantageous for many types of
implementations is obtained by combining pairs of real factors and
complex conjugate pairs into second-order factors.
N  b z 1  b 2
s
b z
H (z)   0k 1k
2
2k

k 1 1  a 1 k z 1  a 2 k z
where Ns is the largest integer contained in (N+1)/2.

139
 w1[n] y1[n] w2[n] y2[n] w3[n] y3[n]
x[n] b01 b02 b03 y[n]

z-1 z-1 z-1

a11 b11 a12 b12 a13 b13

zb-1 zb-1 zb-1

a 21 a 22 a 23

1  2 z 1  z 2 (1  z  1 )( 1  z  1 )
H (z)  
1  0 .75 z  1  0 .125 z  2 (1  0 . 5 z 1 )( 1  0 .25 z  1 )

x[n] y[n]

z-1
z-1 0.5 z-1 z-1 0.25

x[n] y[n]

z-1 z-1
0.5 0.25

140
 Basic Structures for IIR Systems 3

• Parallel Form

B (1  e z )
NP N1 N2 1
A
H (z)  C k
z
k
 1  c kz 1


k k
(1  d z 1 )( 1  d z )
 1

k0 k 1 k k 1 k k

where N = N1+2N2 . If M ณ N, then NP = M - N; otherwise, the first summation in

right hand side of equation above is not included.
• Alternatively, the real poles of H(z) can be grouped in pairs :
 e z 1

e
H (z) 
N NS
P
C z k 

0k 1k
k 1 2

k 0 k 1 1  a 1k z  a 2k z

where NS is the largest integer contained in (N+1)/2, and if NP = M - N is

negative, the first sum is not present.

141
 C0

w1[n] b y1[n]
01

a11 z-1 b
11

a21 z-1

x[n] w2[n] b y2[n] y[n]

02

a12 z-1 b
12

a22 z-1

w3[n] b y3[n]
03

a13 z-1 b
13

a23 z-1 Parallel form structure for

sixth order system (M=N=6).

142
 1 2
H (z)  1  8  7  8 z 1
z2z  2
1  0 .75 z 1  0 .125 z  2 1  0 .75 z 1  0 .125 z
18 25
 8  
1  0 .5 z 1
1  0 .25 z 1

x[n] y[n]

-7 8
z-1
0.75 8 x[n] 18 y[n]

z-1 z-1
-0.125 0.5

-25

0.25 z-1

143
 Transposed Forms

• Transposition (or flow graph reversal) of a flow graph is accomplished by

reversing the directions of all branches in the network while keeping the
branch transmittances as they were and reversing the roles of the input and
output so that source nodes become sink nodes and vice versa.

x[n] y[n] y[n] x[n]

z-1 z-1
a a
x[n] y[n]

z-1 a

144
x[n] b0 y[n]

z-1 b1 a1 z-1

z-1
z-1 b2 a2

bN-1 aN-1

z-1
z-1 bN aN

x[n] b0 y[n]

z-1 a1 b1 z-1

z-1
z-1 a2 b2

aN-1 bN-1

z-1
z-1 aN bN

145
x[n] b0 y[n]
z-1
a1 b1
z-1
a2 b2

aN-1 bN-1

z-1 x[n] b0 y[n]

aN bN
z-1
b1 a1
z-1
b2 a2

bN-1 aN-1

z-1
bN aN

146
 Basic Network Structures for FIR Systems

• Direct Form
– It is also referred to as a tapped delay line structure or a transversal filter
structure.
• Transposed Form
• Cascade Form
M M S
1 2
H ( z )   h [ n ] z  n   ( b0 k  b1 k z  b2k z

)
n0 k 1

where MS is the largest integer contained in (M + 1)/2. If M is odd, one of

coefficients b2k will be zero.

147
 Direct Form

• For causal FIR system, the system function has only zeros (except for
poles at z = 0) with the difference equation:
y[n] = SMk=0 bkx[n-k]
• It can be interpreted as the discrete convolution of x[n] with the
impulse response
h[n] = bn , n = 0, 1, …, M,
0 , otherwise.

148
Direct Form (Tapped Delay Line or Transversal Filter)

x[n] b0 y[n]

z-1 b1

z-1 b2

bN-1

z-1 bN

x[n] z-1 z-1 z-1

h[0] h[1] h[2] h[M-1] h[M] y[n]

149
 Transposed Form of FIR Network

z-1 z-1 z-1 z-1 y[n]

h[M] h[M-1] h[M-2] h[2] h[1] h[0]

x[n]

150
 Cascade Form Structure of a FIR System

x[n] b01 b02 b0Ms y[n]

z-1 b11 z-1 b12 z-1 b1Ms

z-1 b21 z-1 b22 z-1 b2Ms

151
 Structures for Linear-Phase FIR Systems

• Structures for Linear Phase FIR Systems :

– h[M-n] = h[n] for n = 0, 1, ..., M
• For M is an even integer : Type I
y[n] = Sk=0(M/2)-1 h[k](x[n-k] + x[n-M+k]) + h[M/2]x[n-M/2]
• For M is an odd integer : Type II
y[n] = Sk=0(M-1)/2 h[k](x[n-k] + x[n-M+k])
– h[M-n] = -h[n] for n = 0, 1, ..., M
• For M is an even integer : Type III
y[n] = Sk=0(M/2)-1 h[k](x[n-k] - x[n-M+k])
• For M is an odd integer : Type IV
y[n] = Sk=0(M-1)/2 h[k](x[n-k] - x[n-M+k])

152
 Direct form structure for an FIR linear-phase when M is even.

x[n] z-1 z-1 z-1

z-1 z-1

y[n] h[0] h[1] h[2] h[(M/2)-1] h[M/2]

Direct form structure for an FIR linear-phase when M is odd.

x[n] z-1 z-1 z-1

z-1

z-1 z-1 z-1

y[n] h[0] h[1] h[2] h[(M-3)/2] h[(M-1)/2]

153
 Lattice Structures

• Theory of autoregressive signal modeling :

– Lattice Structure
• Development of digital filter structures that are analogous to analog
filter structures :
– Wave Digital Filters
• Another structure development approach is based on state-variable
representations and linear transformations.

154
 Lattice Structures 2

  
• FIR Lattice 
  
Y (z)  N m 
H ( z )   A ( z )  1   a m z
 
X (z)  m 1 


The coefficients {ki} in the lattice structures are referred

to as the k-parameters, which are called reflection
coefficients or PARCOR coefficients.
– When used in signal modeling, the k-parameters
are estimated from a data signal .
– Given a set of k-parameters, the system function
can be found and therefore the impulse response :
by recurrence formula for computing A(z) in terms
of the intermediate system functions

155
Reflection coefficients or PARCOR coefficients structure

x[n] e0[n] e1[n] e2[n] eN-1[n] eN[n] y[n]

-k1 -k2 -kN

-k1 -k2 -kN
z-1 z-1 z-1
e~0[n] e~1[n] e~2[n] e~N-1[n] e~ [n]

Signal flow graph of an FIR lattice system

e0[n] = e~ [n] = x[n]

0
ei~[n] = ei-1[n] – kie~i-~1[n-1], i = 1, 2, …, N,
e [n] = -k e [n] + e [n-1]
i i i-1 i-1
y[n] = eN[n]

156
 A recurrence formula for computing A(z) = H(z) = Y(z)/X(z) can be
obtained in terms of intermediate system functions:

 
 a
Ei (z) i
A (z)   1 (i)
m m 
m 1
z 
i
O 

E (z) 


By recursive technique:
a (i) = k ,
i i
am(i) = am(i-1) - ki ai-m(i-1) ,
m = 1, 2, ..., (i-1)

Or by reverse recursive technique:

ki = a i(i)
am(i-1) = [am(i) + k ia i-m ]/[1 – ki ], m = 1, 2, …, (i – 1)
(i) 2

157
 Example:

A(z) = (1-0.8jz-1)(1+0.8jz-1)(1-0.9z-1) = 1 – 0.9z-1 + 0.64z-2 – 0.576z-3.

Then, a1(3) = 0.9, a 2(3) = -0.64, and a (3) = 0.576
The k-parameter can be computed as follow:
k3 = a (3) = 0.576
3
a (2) = [a (3) + k a (3)]/[1 – k 2] = 0.79518245
1 1 3 2 3
a (2) = [a (3) + k3a1
(3)]/[1 – k 2] = - 0.18197491
2 2 3
k2 = a2(2) = - 0.18197491
+ k2a1 ]/[1 – k2 ] = 0.67275747
a (1) = [a (2) (2) 2
1 1
k1 = a1(1) = 0.67275747

158
x[n] y[n]
-0.6728 +0.182 -0.576
-0.6728 +0.182 -0.576

z-1 z-1 z-1

z–1 z–1 z-1

x[n]
-0.9 0.64 -0.576
y[n]

159
 All-Pole Lattice

• A lattice system with an all-pole system function H(z) = 1/A(z) can be

developed from the FIR lattice .
– all roots of A(z) must be inside the unit circle: |ki| < 1, i = 1, ..., N
ei[n] ei-1[n]
eN[n] = x[n]
Ei-1[n] = ei[n] + kie~ [n-1], i = N, (N-1), …1, ki
i-1
-ki
e~ [n] = -k e [n] + e~ [n-1] z-1 z-1
i i i-1 i-1
~ e~ [n] e~ [n]
y[n] = e0[n] = e 0[n] i i-1

x[n] e0[n] e1[n] e2[n] eN-1[n] eN[n] y[n]

kN-2
kN k1
-kN -k -k1
z-1 N-2
z-1 z-1
e~N[n] e~N-1[n] e~1[n] e~ [n] 160
 Basic all-pole lattice structures

• Three-multiplier form
• Four-multiplier, normalized form
N

 cos i
i 1
H (z) 
A( z)

• Four-multiplier, Kelly-Lochbaum form : was first derived as an acoustic tube

model for speech synthesis.

 (1  k i )
H (z)  i 1

A( z)

161
ei[n] ei-1[n]

-ki ki
Three-multiplier form
e’i[n] (1 - ki2) e’i-1[n]

ei[n] cos qi ei-1[n]

-sin qi sin qi
Four-multiplier, normalized
form
e’i[n] cos qi e’i-1[n]

ei[n] (1 + ki) ei-1[n]

-ki ki Four-multiplier, Kelly-

Lochbaum form
e’i[n] (1 - ki) e’i-1[n]

162
 Lattice Systems with Poles and Zeros

N i 1
Y (z) c z A (z ) B(z)
H (z)  
 
i i

X (z) i0 A( z) A( z )
N
m
B ( z )   bm z
m 0


(i)
bm  c m  ci a i m
i m 1

x[n] = eN[n] eN-1[n] eN-2[n] e1[n] e0[n]

Section Section Section

e’N[n] N e’N-1[n] N - 1e’N-2[n] e’1[n] e’01[n]

cN cN-1 cN-2 c1 c0

y[n]
163
Example of lattice IIR filter with poles and zeros

x[n] y[n]

0.9
z-1 3

-0.64
z-1 3

0.576
z-1

x[n]

0.576 -0.182 0.6728

-0.576 0.182 -0.6728

z-1 z-1 5.4612 -1

z4.5404
3.9

y[n]

164
 UNIT-2
DFS , DFT & FFT

165
Fourier representation of signals

166
Fourier representation of signals

167
Fourier representation of signals

168
Discrete complex exponentials

169
 Discrete Fourier Series ~
• Given a periodic sequence x [n] with period N so that
~ ~
x [n]  x [n  rN ]

• The Fourier series~ repre1sentation can be written as

 X k e
~
x [n]  j 2  / N kn

N k

• The Fourier series representation of continuous-time

periodic signals require infinite many complex exponentials
• Not that for discrete-time periodic signals we have
j2  / N k  mN n j2  / N kn j2  mn  j2  / N kn
e  e e  e

• Due to the periodicity of the complex exponential we only

need N exponentials for discrete time Fourier series
~
N 1 j 2  / N kn
~ 1
x [n] 
N  X k  e
k 0

170
 Discrete Fourier Series Pair
• A periodic sequence in terms of Fourier series coefficients
1 ~  
N 1 j 2  / N kn
~
x [n] 
 X k e
N k 0

• The Fourier series coefficients can be obtained via

~
 N 1
~
 j 2  / N kn

X k  
n0
x [n]e

• For convenience we sometimes use

 j 2  / N 
• Analysis equation WN  e

~
 

N 1 kn
~
X k   x [n]W N
• Synthesis equation n0

~ 1 N 1
~  kn
x [n] 
 X k  W N
N k 0

171
Fourier series for discrete-time periodic
signals

172
Discrete-time Fourier series
(DTFS)

173
Fourier representation of aperiodic
signals

174
Discrete-time Fourier transform
(DTFT)

175
Discrete Fourier Transform

176
Discrete Fourier Transform

177
Discrete Fourier Transform

178
Discrete Fourier Transform

179
Discrete Fourier Transform

180
Summary of properties

181
 DFT Pair & Properties
• The Discrete Fourier Transform pair
N 1 N 1
 1
X k  
 j 2  / N kn
X k  e
j2 / N  kn
 x[n]e x[ n ] 
N 
n0 k 0

182
Circular convolution

183
Modulo Indices and Periodic
Repetition

184
Overlap During Periodic
Repetition

185
Periodic repetition: N=4

186
Periodic repetition: N=4

187
Modulo Indices and the Periodic
Repetition

188
Modulo Indices and the Periodic
Repetition

189
Modulo Indices and the Periodic
Repetition

190
Modulo Indices and the Periodic
Repetition

191
Circular convolution

192
Circular convolution

193
Circular convolution-another
interpretation

194
Using DFT for Linear
Convolution

195
Using DFT for Linear Convolution

196
Using DFT for Linear Convolution

197
Using DFT for Linear
Convolution

198
Using DFT for Linear
Convolution

199
Using DFT for cicular
Convolution

200
Using DFT for cicular
Convolution

201
Using DFT for cicular
Convolution

202
Filtering of Long Data Sequences

203
Filtering of Long Data Sequences

204
Over-lap Add

205
Over-lap Add

206
Over-lap Add

207
Over-lap Add

208
Over-lap Add

209
Over-lap Add

210
Over-lap Add

211
Over-lap save

212
Over-lap save

213
Over-lap save

214
Over-lap save input segment stage

215
Over-lap save input segment stage

216
Over-lap save input segment stage

217
Over-lap save filtering stage

218
Over-lap save output blocks

219
Over-lap save output blocks

220
Over-lap save output blocks

221
Over-lap save

222
Over-lap save

223
Relationships between CTFT,
DTFT, & DFT

224
Relationships between CTFT,
DTFT, & DFT

225
Relationships between CTFT, DTFT,
& DFT

226
Fast Fourier Transform

227
 Discrete Fourier Transform
• The DFT pair was given as
N 1 N 1 
1  
X k    x[n]e
 j 2  / N kn x[n] 
N
 X k e j 2  / N kn

n0 k 0

• Baseline for computational complexity:

– Each DFT coefficient requires
• N complex multiplications
• N-1 complex additions
– All N DFT coefficients require
• N2 complex multiplications
• N(N-1) complex additions

• Complexity in terms of real operations

• 4N2 real multiplications
• 2N(N-1) real additions

• Most fast methods are based on symmetry properties

– Conjugate symmetry e
 j2  / N k N  n 
 e
 j2  / N kN
e
 j2  / N k   n 
 e
j2  / N kn

– Periodicity in n and k e
 j2  / N kn
 e
 j2  / N k  n  N 
 e
j2  / N k  N n

228
Direct computation of DFT

229
Direct computation of DFT

230
231
FFT

232
 Decimation-In-Time FFT Algorithms
• Makes use of both symmetry and periodicity
• Consider special case of N an integer power of 2
• Separate x[n] into two sequence of length N/2
– Even indexed samples in the first sequence
– Odd indexed samples in the other sequence
N 1 N 1 N 1

X k    x[n]e  x[n]e
 j 2  / N kn  j 2  / N kn  j 2  / N kn
 x[n]e  
n0 n even n odd

• Substitute variables n=2r for n even and n=2r+1 for odd

N / 2 1 N / 2 1

X k    x[2 r ]W 2 rk  x[2 r
2 r  1 k
N
  1]W N
r0 r0

N / 2 1 N / 2 1

 x[2 r ]W rk  x[2 r
k rk
 N /2
 W N  1]W N / 2
r0 r0

 G k   W k H
N
k 
• G[k] and H[k] are the N/2-point DFT’s of each subsequence
233
 Decimation In Time

• 8-point DFT example using

decimation-in-time
• Two N/2-point DFTs
– 2(N/2)2 complex multiplications
– 2(N/2)2 complex additions
• Combining the DFT outputs
– N complex multiplications
– N complex additions
• Total complexity
– N2/2+N complex multiplications
– N2/2+N complex additions
– More efficient than direct DFT
• Repeat same process
– Divide N/2-point DFTs into
– Two N/4-point DFTs
– Combine outputs

234
 Decimation In Time Cont’d
• After two steps of decimation in time

• Repeat until we’re left with two-point DFT’s

235
 Decimation-In-Time FFT Algorithm
• Final flow graph for 8-point decimation in time

• Complexity:
– Nlog2N complex multiplications and additions

236
 Butterfly Computation
• Flow graph constitutes of butterflies

• We can implement each butterfly with one multiplication

• Final complexity for decimation-in-time FFT

– (N/2)log2N complex multiplications and additions

237
 In-Place Computation
• Decimation-in-time flow graphs require two sets of registers
– Input and output for each stage
• Note the arrangement of the input indices
– Bit reversed indexing

X 0   x 0  
0 X 0 000   x 000 
X 0 1   x 4   X 0 001   x 100 

X 0 2   x 2   X 0 010   x 010 

X 0 3   x 6   X 0 011   x 110 

X 0 4   x 1   X 0 100   x 001 

X 0 5   x 5   X 0 101   x 101 

X 0 6   x 3   X 0 110   x 011 

X 0 7   x 7   X 0 111   x 111 

238
 Decimation-In-Frequency FFT Algorithm
• The DFT equation 


N 1

X k    x[n]W nk N
n0

• Split the DFT equation into even and odd frequency indexes
N 1 N / 2 1 N 1

X 2 r    x[n]W n2r
N   x[n]W n2r
N   x[n]W n2r
N
n0 n0 nN / 2

• Substitute variables to get

N / 2 1 N / 2 1 N / 2 1

X 2 r    x[n]W n 2 r  x[n   x[n]  x[n  N / 2 ] W

n  N / 2 2 r
N
  N / 2 ]W N
 nr
N/2
n0 n0 n0

• Similarly for odd-numbered frequencies N / 2 1

X 2 r  1     x[n]  x[n  N / 2 ] W n 2 r  1 
N /2

n0

239
 Decimation-In-Frequency FFT Algorithm
• Final flow graph for 8-point decimation in frequency

240
 UNIT-3
IIR filters

241
 Filter Design Techniques
• Any discrete-time system that modifies certain frequencies
• Frequency-selective filters pass only certain frequencies
• Filter Design Steps
– Specification
• Problem or application specific
– Approximation of specification with a discrete-time system
• Our focus is to go from spec to discrete-time system
– Implementation
• Realization of discrete-time systems depends on target technology
• We already studied the use of discrete-time systems to implement a
continuous-time system
– If our specifications are given in continuous time we can use

x n  y n 
xc(t) C/D H(ej) D/C yr(t)


H e j   H c j / T    

242
 Digital Filter Specifications


• Only the magnitude approximation problem

• Four basic types of ideal filters with magnitude responses
as shown below (Piecewise flat)

243
 Digital Filter Specifications
• These filters are unealisable because (one of the
following is sufficient)
– their impulse responses infinitely long non-
causal
– Their amplitude responses cannot be equal to a
constant over a band of frequencies
Another perspective that provides some
understanding can be obtained by looking at the
ideal amplitude squared.

244
 Digital Filter Specifications
• The realisable squared amplitude response transfer
function (and its differential) is continuous in
• Such functions 
– if IIR can be infinite at point but around that
point cannot be zero.
– if FIR cannot be infinite anywhere.
• Hence previous differential of ideal response is
unrealisable

245
 Digital Filter Specifications
• For example the magnitude response of a digital
lowpass filter may be given as indicated below

246
 Digital Filter Specifications

• In the passband 0   we require

p

that G ( e j
) 1 with a deviation  p

1 p  G (e j
)  1  p,   p

• In the stopband  s     we re qs uire

that G ( e  )  0 with a deviation
j

G (e j
)   s,  s    
247


Digital Filter Specifications 




Filter specification parameters

•  p
- passband edge frequency
• s - stopband edge frequency
•  p - peak ripple value in the passband
• s - peak ripple value in the stopband

248
 Digital Filter Specifications

• Practical specifications are often given in

terms of loss function (in dB)
• G ( )   20 log 10 G ( e )
j

• Peak passband ripple

 p   20 log 10 (1   p )
dB
• Minimum stopband attenuation
 s   20 log 10 ( s ) dB

249
 Digital Filter Specifications

• In practice, passband edge frequency F p and

stopband edge frequency arFes specified in
Hz
• For digital filter design, normalized bandedge
frequencies need to be computed from
specifications in Hz using  2 F
 p 
p

p
 2 Fp T
FT FT
s 2 Fs
s    2 F s T
FT FT 250
 Digital Filter Specifications
• Example - Let F p
 7 kHz, F s
 3
kHz, and F  25
T
kHz
• Then
2  ( 7  10 3 )
 p  3  0 .56 
25  10 
2  ( 3  10 3 ) 
 s  3  0 .24 
25  10

251
 IIR Digital Filter Design






Standard approach
(1) Convert the digital filter specifications into
an analogue prototype lowpass filter
specifications
(2) Determine the analogue lowpass filter
transfer function H a ( s )
(3) Transform H a ( s ) by replacing the complex
variable to the digital transfer function
G (z)
252
 IIR Digital Filter Design
• This approach has been widely used for the
following reasons:
(1) Analogue approximation techniques are
highly advanced
(2) They usually yield closed-form
solutions
(3) Extensive tables are available for
analogue filter design
(4) Very often applications require digital
simulation of analogue systems
253
 IIR Digital Filter Design

• Let an analogue transfer function be

Pa ( s )
H (s) 
a
Da (s)

where the subscript “a” indicates the

analogue domain
• A digital transfer function derived from this
is denoted as
P (z)
G ( z ) 
D (z)
254
 IIR Digital Filter Design
• Basic idea behind the conversion of H a ( s ) intoG ( z )
is to apply a mapping from the s-domain to the z-
domain so that essential properties of the analogue
frequency response are preserved
• Thus mapping function should be such that
– Imaginary (j  ) axis in the s-plane be
mapped onto the unit circle of the z-plane
– A stable analogue transfer function be mapped
into a stable digital transfer function

255
Specification for effective frequency response of a continuous-time lowpass
filter and its corresponding specifications for discrete-time system.

dp or d1 passband ripple
ds or d2 stopband ripple
Wp, wp passband edge frequency
Ws, ws stopband edge frequency
e2 passband ripple parameter

1 – dp = 1/1 + e2

BW bandwidth = wu – wl
wc 3-dB cutoff frequency
wu, wl upper and lower 3-dB cutoff
frequensies
Dw transition band = |wp – ws|
Ap passband ripple in dB
=  20log10(1  dp)
As stopband attenuation in dB
= -20log10(ds)

256
 Design of Discrete-Time IIR Filters

• From Analog (Continuous-Time) Filters

– Approximation of Derivatives
– Impulse Invariance
– the Bilinear Transformation

257
 Reasons of Design of Discrete-Time IIR Filters from
Continuous-Time Filters

• The art of continuous-time IIR filter design is highly advanced and,

since useful results can be achieved, it is advantageous to use the
design procedures already developed for continuous-time filters.

• Many useful continuous-time IIR design methods have relatively

simple closed-form design formulas. Therefore, discrete-time IIR
filter design methods based on such standard continuous-time design
formulas are rather simple to carry out.

• The standard approximation methods that work well for continuous-

time IIR filters do not lead to simple closed-form design formulas
when these methods are applied directly to the discrete-time IIR case.

258
 Characteristics of Commonly Used Analog Filters

• Butterworth Filter
• Chebyshev Filter
– Chebyshev Type I
– Chebyshev Type II of Inverse Chebyshev Filter

259
 Butterworth Filter
• Lowpass Butterworth filters are all-pole filters characterized by the magnitude-squared
frequency response

|H(W)|2 = 1/[1 + (W/Wc)2N] = 1/[1 + e2(W/Wp)2N]

where N is the order of the filter, Wc is its – 3-dB frequency (cutoff frequency), Wp is
the bandpass edge frequency, and 1/(1 + e2) is the band-edge value of |H(W)|2.

• At W = Ws (where Ws is the stopband edge frequency) we have

1/[1 + e2(Ws/Wp)2N] = d22
and
N = (1/2)log10[(1/d22) – 1]/log10(Ws/Wc) = log10(d/e)/log10(Ws/Wp)
where d2= 1/1 + d22.

• Thus the Butterworth filter is completely characterized by the parameters N, d2, e, and
the ratio Ws/Wp.

260
 Butterworth Lowpass Filters
• Passband is designed to be maximally flat
• The magnitude-squared function is of the form
 
2 1  1
 j 
2
H  H s 
1   j / j c 
c
1  s / j c 2 N
c 2N

s k   1   j c 
1 / 2N  j / 2 N  2 k  N  1 
  ce for k  0,1,...,2N -1

261
Frequency response of lowpass Butterworth filters

262
 Chebyshev Filters
• The magnitude squared response of the analog lowpass Type I Chebyshev
filter of Nth order is given by:

|H(W)|2 = 1/[1 + e2TN2(W/W )].

where TN(W) is the Chebyshev polynomial of order N:
TN(W) = cos(Ncos-1 W), |W|  1,
= cosh(Ncosh-1 W), |W| > 1.

The polynomial can be derived via a recurrence relation given by

Tr(W) = 2WTr-1(W) – Tr-2(W), r  2,
with T0(W) = 1 and T1(W) = W.

• The magnitude squared response of the analog lowpass Type II or inverse

Chebyshev filter of Nth order is given by:
|H(W)|2 = 1/[1 + e2{TN(Ws/Wp)/ TN(Ws/W)}2].

263
 Chebyshev Filters
• Equiripple in the passband and monotonic in the stopband
• Or equiripple in the stopband
2 1
and monotonic in the passband
Hc  j   V N x   cos N cos 1
x 
1   V N  / 
2 2
c 




264






Frequency response of
lowpass Type I Chebyshev filter

|H(W)|2 = 1/[1 + e2TN2(W/W )]

Frequency response of
lowpass Type II Chebyshev filter

|H(W)|2 = 1/[1 + e2{TN2(Ws/Wp)/TN2(Ws/W)}]

265
N = log10[( 1 - d 2 +  1 – d 2(1 + e2))/ed ]/log [(W /W ) +  (W /W )2 – 1 ]
2 2 2 10 s p s p
= [cosh-1(d/e)]/[cosh-1(Ws/Wp)]

for both Type I and II Chebyshev filters, and where

d2 = 1/  1 + d2.

• The poles of a Type I Chebyshev filter lie on an ellipse in the s-plane with major
axis r1 = Wp{(b2 + 1)/2b] and minor axis r1 = Wp{(b2 - 1)/2b] where b is related to
e according to
b = {[ 1 + e2 + 1]/e}1/N
• The zeros of a Type II Chebyshev filter are located on the imaginary axis.

266
 Type I: pole positions are
xk = r2cosfk
yk = r1sinfk
fk = [p/2] + [(2k + 1)p/2N]
r1 = Wp[b2 + 1]/2b
r2 = Wp[b2 – 1]/2b

b = {[ 1 + e2 + 1]/e}1/N

Type II: zero positions are

sk = jWs/sinfk
and pole positions are

vk = Wsxk/ xk 2 + yk 2

wk = Wsyk/ xk2 + y 2k

b = {[1 +  1 – d 22 ]/d }1/N

Determination of the pole locations k = 0, 1, …, N-1.

for a Chebyshev filter.
267
 Approximation of Derivative Method

• Approximation of derivative method is the simplest one for converting an

analog filter into a digital filter by approximating the differential equation by
an equivalent difference equation.
– For the derivative dy(t)/dt at time t = nT, we substitute the backward difference
[y(nT) – y(nT – T)]/T. Thus
dy ( t ) y ( nT )  y ( nT T ) y [ n ]  y [ n  1]
 
dt t  nT
T T

where T represents the sampling period. Then, s = (1 – z-1)/T

– The second derivative d2y(t)/dt2 is derived into second difference as follow:

dy ( t ) y [ n ]  2 y [ n  1 ]  y [ n  2 ]
 2
dt t  nT
T
which s2 = [(1 – z-1)/T]2. So, for the kth derivative of y(t), sk = [(1 – z-1)/T]k.

268
 Approximation of Derivative Method
• Hence, the system function for the digital IIR filter obtained as a result of the
approximation of the derivatives by finite difference is
H(z) = Ha(s)|s=(z-1)/Tz

• It is clear that points in the LHP of the s-plane are mapped into the
corresponding points inside the unit circle in the z-plane and points in the
RHP of the s-plane are mapped into points outside this circle.
– Consequently, a stable analog filter is transformed into a stable digital filter due
to this mapping property.

jW
Unit circle

s-plane
z-plane
269
 Example: Approximation of derivative method

Convert the analog bandpass filter with system function

Ha(s) = 1/[(s + 0.1)2 + 9]
Into a digital IIR filter by use of the backward difference for the derivative.

Substitute for s = (1 – z-1)/T into Ha(s) yields

H(z) = 1/[((1 – z-1)/T) + 0.1)2 + 9]

 T 2

H ( z )  1  0 .2 T  9 .01 T
1
2

1
2 ( 1  0 .1 T )
z  1
z
2
1  0 .2 T  9 .01 T 2
1  0 .2 T  9 .01 T 2

T can be selected to satisfied specification of designed filter. For example, if T = 0.1,

the poles are located at
p1,2 = 0.91  j0.27 = 0.949exp[ j16.5o]

270
 Filter Design by Impulse Invariance
• Remember impulse invariance
– Mapping a continuous-time impulse response to discrete-time
– Mapping a continuous-time frequency response to discrete-time
h n   T d h c nT d


   2 

H e j  H c
j  j k 
k    Td Td 


• If the continuous-time filter is bandlimited to


H  j 
c
 0    / Td
j   
H e  H c
j 

  
 T d 


• If we start from discrete-time specifications Td cancels out

– Start with discrete-time spec in terms of 
– Go to continuous-time =  /T and design continuous-time filter
– Use impulse invariance to map it back to discrete-time = T
• Works best for bandlimited filters due to possible aliasing

271
 Impulse Invariance of System Functions
• Develop impulse invariance relation between system functions
• Partial fraction expansion of transfer function
Ak
s 
N
H
c
 
k 1 s  sk
• Corresponding impulse response
N s t
 A k e k t  0
hc t   k  1
 0 t  0
• Impulse response of discrete-time filter
   
n
   
sk Td
N

s k nT N
d

h n  T dh c nT d  
k 1
Td A k e u n   Td A k e u n
• System function k 1

 N TdA k

H z   sk Td 1
k 1 1  e z
s Td
• Pole s=sk in s-domain transform into pole at e
k

272
 Impulse Invariant Algorithm

• Step 1: define specifications of filter

– Ripple in frequency bands
– Critical frequencies: passband edge, stopband edge, and/or cutoff frequencies.
– Filter band type: lowpass, highpass, bandpass, bandstop.
• Step 2: linear transform critical frequencies as follow
W = w/Td
• Step 3: select filter structure type and its order: Bessel, Butterworth, Chebyshev
type I, Chebyshev type II or inverse Chebyshev, elliptic.
• Step 4: convert Ha(s) to H(z) using linear transform in step 2.
• Step 5: verify the result. If it does not meet requirement, return to step 3.

273
 Example: Impulse invariant method

Convert the analog filter with system function

Ha(s) = [s + 0.1]/[(s + 0.1)2 + 9]
into a digital IIR filter by means of the impulse invariance method.

The analog filter has a zero at s = - 0.1 and a pair of complex conjugate poles at pk = - 0.1  j3.
Thus,
H s  
1 1

a
2
 2

s  0 .1  j 3 s  0 .1  j 3

H  z  
1 1
Then
2
 0 .1 T 1  2
 0 .1 T  j 3T 1
1 e 1 e
j 3T
e z e z

274
Frequency response
of digital filter.

Frequency response
of analog filter.

275
 Disadvantage of previous
techniques: frequency
warping  aliasing effect
and error in specifications
of designed filter (frequencies)
So, prewarping of frequency
is considered.












276
 Example
• Impulse invariance applied to Butterworth
0 . 89125 
 H e j   1 0    0 . 2 


H e j
  0.17783 0 . 3     

• Since sampling rate Td cancels out we can assume Td=1

• Map spec to continuous time

0 . 89125  H  j   1 0    0 . 2 

H  j   0.17783 0 . 3     

• Butterworth filter is monotonic so spec will be satisfied if

H c j0 .2    0 .89125 and H c j0 .3    0.17783

H
c
 j 
2

1   j
1

/ j

c
2 N
• Determine N and c to satisfy these conditions

277
 Example Cont’d
• Satisfy both constrains
 0 .2  2 N  1 2  0 .3  2 N  1 2
1       and 1      
   
 c   0 . 89125   c   0 . 17783 

• Solve these equations to get

N  5 . 8858  6 and  c
 0 . 70474
• N must be an integer so we round it up to meet the spec
• Poles of transfer function
1 / 12  j   j 2 k  11 
 1 
/ 12
sk   ce for k  0,1,...,11
• c

The transfer function

H s   0 . 12093
s 2
 0 . 364 s  0 . 4945 s 2
 0 . 9945 s  0 . 4945 s 2
 1 . 3585 s  0 . 4945 
• Mapping to z-domain
0 . 2871  0 . 4466 z  1  2 . 1428  1 . 1455 z  1
H z  

1 2
1  1 . 2971 z  0 . 6949 z 1  1 . 0691 z  1  0 . 3699 z  2
1 . 8557  0 . 6303 z  1

1  0 . 9972 z  1  0 . 257 z  2
278
Example Cont’d

279
 Filter Design by Bilinear Transformation
• Get around the aliasing problem of impulse invariance
• Map the entire s-plane onto the unit-circle in the z-plane
– Nonlinear transformation
– Frequency response subject to warping
• Bilinear transformation
s 
2  1  z  1
 

  z  1 
 T d  1 
• Transformed system function
H z   H  2  1   1 
 z
c  
 1  
T d 1  z 
 
• Again Td cancels out so we can ignore it
• We can solve the transformation for z as
1  T / 2 s 1  T / 2  j T /2
z  d

 d d s    j
1  T d
/ 2 s 1   T d / 2  j T d / 2

• Maps the left-half s-plane into the inside of the unit-circle in z

– Stable in one domain would stay in the other

280
 Bilinear Transformation
• On the unit circle the transform becomes
1  j T d / 2 j
z  e
1  j T d / 2



• To derive the relation between  and 

2  1  e  j  2  2 e  j / 2 j sin  / 2   2j   
s       j    tan  
 j / 2 



cos  / 2  


T d   j  Td  2 e
1  e  Td  2 


• Which yields 
2      Td 
  tan   or   2 arctan  
Td 2   2 

281
Bilinear Transformation

282
 Example
• Bilinear transform applied to Butterworth
0 . 89125  H e j   1 0    0 . 2 

H e j   0.17783 0 . 3     
• Apply bilinear transformation to specifications
2  0 .2  
0 . 89125  H  j   1 0    tan  
 Td  2 
2  0 .3  
H  j   0.17783 tan      
Td  2 
• We can assume Td=1 and apply the specifications to
2 1
H  j  
1   /  
2N
• To get c
c

 2 tan 0 . 1  2 N  1 2  2 tan 0 . 15  2 N  1 2

1   
  and 1      
  
  0 . 89125    0 . 17783 
 c   c 

283
 • Example Cont’d
Solve N and  c
  1
2
  
2

log    1   1  

  1    0 . 89125 
  0 . 17783 
 
 
    0 . 766
N   5 . 305  6 c
2 log  tan 0 . 15   tan 0 . 1  

• The resulting transfer function has the following poles

 j 2 k  11 
sk    1   j c  
1 / 12 / 12
 ce for k  0,1,...,11

• Resulting in
0 . 20238
H s  
c
s 2
 0 . 3996 s  0 . 5871 s 2
 1 . 0836 s  0 . 5871 s 2
 1 . 4802 s  0 . 5871 


• Applying the bilinear transform yields

0 . 0007378 1  z 
1 6

H z    
1  1 . 2686 z
1 
 0 . 7051 z  2 1  1 . 0106 z  1  0 . 3583 z  2 
  1

 
1  0 . 9044 z
1
 0 . 2155 z  2

284
Example Cont’d

285
 IIR Digital Filter: The bilinear
transformation 



• To obtain G(z) replace s by f(z) in H(s)

• Start with requirements on G(z)
G(z) Available H(s)
Stable Stable
Real and Rational in z Real and Rational
in s
Order n Order n
L.P. (lowpass) cutoff L.P. cutoff  T
c c

286
 Bilinear Transformation
• Mapping of s-plane into the z-plane

287
 Bilinear Transformation
• For z  e j
with unity scalar we have
1  e  j 
j    j tan(  / 2 )
 j
1  e

or   tan(  / 2 )

288
 Bilinear Transformation
• Mapping is highly nonlinear
• Complete negative imaginary axis in the s-
plane from    to   0 is mapped into
the lower half of the unit circle in the z-plane
from z   1 to z  1
• Complete positive imaginary axis in the s-
plane from   0 to   is mapped into the
upper half of the unit circle in the z-plane
from z  1 toz   1
289
 Bilinear Transformation
• Nonlinear mapping introduces a distortion
in the frequency axis called frequency
warping
• Effect of warping shown below

290
 Spectral Transformations
• To transform G L ( z ) a given lowpass transfer
function to another transfer function G D ( zˆ)
that may be a lowpass, highpass, bandpass or
bandstop filter (solutions given by
Constantinides)
• 1 has been used to denote the unit delay in
z 1
the prototype lowpass filter G L ( z ) and zˆ
to denote the unit delay in the transformed
filter G D ( zˆ) to avoid confusion
291
 Spectral Transformations
• Unit circles in z- and zˆ -planes defined by
z  e j , zˆ  e jˆ
• Transformation from z-domain to
zˆ-domain given by

• Then
z  F ( zˆ)

G D ( zˆ)  G L { F ( zˆ)}

292
 Spectral Transformations
• From z  F ( zˆ) , thusz  F ( zˆ) , hence
  1, if z  1

F ( zˆ)   1, if z  1
 1, if z  1

• Therefore1 / F ( zˆ) must be a stable allpass function
1 L  1   * zˆ
𝑙
    , 𝑙  1
F ( zˆ) 𝑙 1  zˆ   𝑙 

293
 Lowpass-to-Lowpass
Spectral Transformation
• To transform a lowpass filterG L ( z ) with a cutoff
frequency  c to another lowpass filter G D ( zˆ)
with a cutoff frequency ˆc, the transformation is
1 1 1   zˆ
z  
F ( zˆ) zˆ  

• On the unit circle we havjê

 j e 
e 
1   e  jˆ
which yields 1   
tan(  / 2 )   tan( ˆ / 2 )
 
 1    294





Lowpass-to-Lowpass


Spectral Transformation
• Solving we get sin (   ˆ ) / 2 
  c c

sin (  c  ˆc ) / 2 
• Example - Consider the lowpass digital filter
0 .0662 (1  z 1 ) 3
G L ( z )  1 1 2
(1  0 .2593 z )( 1  0 .6763 z  0 .3917 z )
which has a passband from dc to 0 . 25  with
a 0.5 dB ripple
• Redesign the above filter to move the
passband edge to
0 . 35  295
 Lowpass-to-Lowpass
Spectral Transformation
• Here sin( 0 .05  )
     0 .1934
sin( 0 .3  )
•
Hence, the desired lowpass transfer function is
( zˆ)  G L ( z )
1
GD  1  zˆ  0 .1934
z  1
1  0 .1934 zˆ

-10
Gain, dB

G (z) G (z)
L D
-20

-30

-40
0 0.2 0.4 0.6 0.8 1

Lowpass-to-Lowpass



Spectral Transformation

296
 Lowpass-to-Lowpass
Spectral Transformation
• The lowpass-to-lowpass transformation
1 1 1   zˆ
z  
F ( zˆ) zˆ  

can also be used as highpass-to-highpass,

bandpass-to-bandpass and bandstop-to-
bandstop transformations

297
 Lowpass-to-Highpass
Spectral Transformation
• Desired transformation
1
z 1   zˆ  
1
1   zˆ
• The transformation parameter  is given by
cos (  c  ˆc ) / 2 
 
cos (  c  ˆc ) / 2 
where  c is the cutoff frequency of the lowpass
filter and ˆc is the cutoff frequency of the desired
highpass filter

298
 Lowpass-to-Highpass
Spectral Transformation
• Example - Transform the lowpass filter
0 .0662 (1  z 1 ) 3
G L ( z )  1 1 2
(1  0 .2593 z )( 1  0 .6763 z  0 .3917 z )

• with a passband edge at 0 . 25  to a highpass

filter with a passband edge at 0 .55 
• Here    cos( 0 .4 ) / cos( 0 .15  )   0 .3468
• The desired transformation is
1
z 1   zˆ  0 .3468
1
1  0 .3468 zˆ
299
 Lowpass-to-Highpass
Spectral Transformation

• The desired highpass filter is

1
G D ( zˆ)  G ( z ) z 1

zˆ  0 .3468
1
1 0 .3468 zˆ




Gain, dB







     
Normalized frequency
300
 Lowpass-to-Highpass
Spectral Transformation
• The lowpass-to-highpass transformation can
also be used to transform a highpass filter with
a cutoff at  c to a lowpass filter with a cutoff
at ˆc
• and transform a bandpass filter with a center
frequency at  o to a bandstop filter with a
center frequency at ˆo

301
 Lowpass-to-Bandpass
Spectral Transformation
• Desired transformation
2  2   1
1
zˆ  zˆ 
z 1   1
 1   1
2
2  1
zˆ  1
zˆ 
1 1

302
 Lowpass-to-Bandpass
Spectral Transformation
• The parameters  and  are given by

cos (ˆ c 2  ˆc 1 ) / 2 

 
cos (ˆ c 2  ˆc 1 ) / 2 


  cot (ˆ c 2  ˆc 1 ) / 2  tan(  c / 2 )

where  c is the cutoff frequency of the lowpass
filter, and ˆc 1 and ˆc 2 are the desired upper and
lower cutoff frequencies of the bandpass filter

303
 Lowpass-to-Bandpass
Spectral Transformation
• Special Case - The transformation can be
simplified if  c  ˆc 2  ˆc 1
• Then the transformation reduces to
1
 zˆ  
z 1   zˆ1
1
1   zˆ
where   cos ˆo with ˆo
denoting the
desired center frequency of the bandpass filter

304
 Lowpass-to-Bandstop
Spectral Transformation
• Desired transformation

2  2 1 1  
zˆ  zˆ 
1 1  1  
z 
1  2 2 1
zˆ  1
zˆ 
1  1  

305
 Lowpass-to-Bandstop
Spectral Transformation
• The parameters  and  are given by
cos (ˆ c 2  ˆc 1 ) / 2 
 
cos (ˆ c 2  ˆc 1 ) / 2 

  tan (ˆ c 2  ˆc 1 ) / 2  tan(  c / 2 )

where  c is the cutoff frequency of the
lowpass filter, and ˆc 1 and ˆc 2 are the desired
upper and lower cutoff frequencies of the
bandstop filter

306
 UNIT-4
FIR Filters

307
 Selection of Filter Type

• The transfer function H(z) meeting the

specifications must be a causal transfer
function
• For IIR real digital filter the transfer
1
function is a real rational function of z
1 2 M
H ( z )  p 0  p1 z  p 2 z   p M z
1 2 N
d 0  d1z  d 2 z   d N z
• H(z) must be stable and of lowest order N or
M for reduced computational complexity 308
 Selection of Filter Type

• FIR real digital filter transfer function is a

polynomial in z  1 (order N) with real
coefficients N
n
H ( z )   h [ n ] z
n0
• For reduced computational complexity, degree N
of H(z) must be as small as possible
• If a linear phase is desired then we must have:

h[ n ]   h[ N  n ]
309
 Selection of Filter Type
• Advantages in using an FIR filter -
(1) Can be designed with exact linear phase
(2) Filter structure always stable with quantised
coefficients
• Disadvantages in using an FIR filter - Order of an
FIR filter is considerably higher than that of an
equivalent IIR filter meeting the same
specifications; this leads to higher computational
complexity for FIR

310
 FIR Filter Design
Digital filters with finite-duration impulse response (all-zero, or FIR filters)
have both advantages and disadvantages compared to infinite-duration
impulse response (IIR) filters.
FIR filters have the following primary advantages:

• They can have exactly linear phase.

• They are always stable.
• The design methods are generally linear.
• They can be realized efficiently in hardware.
• The filter startup transients have finite duration.

The primary disadvantage of FIR filters is that they often require a much
higher filter order than IIR filters to achieve a given level of performance.
Correspondingly, the delay of these filters is often much greater than for an
equal performance IIR filter.
 FIR Design
FIR Digital Filter Design
Three commonly used approaches to FIR
filter design -
(1) Windowed Fourier series approach
(2) Frequency sampling approach
(3) Computer-based optimization methods

312
 Finite Impulse Response Filters
• The transfer function is given by
N 1
n
H ( z )   h ( n ). z
n0

• The length of Impulse Response is N

• All poles are at z  0.
• Zeros can be placed anywhere on the z-
plane

313
 FIR: Linear phase

For phase linearity the FIR transfer

function must have zeros outside the
unit circle

314
 Linear Phase
• What is linear phase?
• Ans: The phase is a straight line in the passband of
the system.
• Example: linear phase (all pass system)
• I Group delay is given by the negative of the slope
of the line

315
 Linear phase
• linear phase (low pass system)
• Linear characteristics only need to pertain to
the passband frequencies only.

316
 FIR: Linear phase
• For Linear Phase t.f. (order N-1)
• h (n)   h( N  1  n )

• so that for N even:

N
 2
1 N 1 

H ( z )   h ( n ). z  n   h ( n ). z n

n0 n N
2

N 1 N 1
2 2
  h ( n ). z n
  h ( N  1  n ). z  ( N 1  n )

n0 n0
N 1

 
2
  h ( n ) z n
 z m m  N 1 n
n0 317
 FIR: Linear phase
• for N odd:
N 2 1 
1


 N  1    N 2 1 
 
 

n  
H ( z )   h ( n ). z  z m  h   z
n0  2 


• I) On C : z 1 we have for N even, and

+ve sign 
j T  j T  N  1  N 2  1  
H (e )  e 2
2 h ( n ). cos  T n  N  1  
 
.     


n0   2  
318
 FIR: Linear phase
• II) While for –ve sign
j T  j T  N  1  N 2  1  
H (e )  e 2
j 2 h ( n ). sin T n  N  1  
 
.     


n0   2  

• [Note: antisymmetric case adds  / 2 rads to

phase, with discontinuity at  0 ]
• III) For N odd with

+ve signN  1
 j T  N 1 

H (e j T

 2 
   
)  e h
  
  2 
N 3
 
2   N  1   
  2 h ( n ). cos  T  n    

n0   2    319

 FIR: Linear phase
• IV) While with a –ve sign
j T  j T  N 1 
 2  N 2 3 


H (e )   j .h ( n ). sin   N  1   
  2  T  n 

e
   
 n0   2   
 


• [Notice that for the antisymmetric case to have

linear phase we require
 N  1 
h   0.
 2 

The phase discontinuity is as for N even]
320
 FIR: Linear phase
• The cases most commonly used in filter
design are (I) and (III), for which the
amplitude characteristic can be written as a
polynomial in
T
cos
2

321
 Summary of Properties
K

H   F   a k cos k  
j 0  j N / 2
e e
k0





Type I II III IV
Order N even odd even odd
Symmetry symmetric symmetric anti-symmetric anti-symmetric
Period 2 4 2 4
0 0 0 /2 /2
F() 1 cos(/2) sin() sin(/2)
K N/2 (N-1)/2 (N-2)/2 (N-1)/2
H(0) arbitrary arbitrary 0 0
H() arbitrary 0 0 arbitrary


Design of FIR filters: Windows 





(i) Start with ideal infinite duration h ( n )

(ii) Truncate to finite length. (This produces
unwanted ripples increasing in height near
discontinuity.)
~
(iii) Modify to h ( n )  h ( n ). w ( n )
Weight w(n) is the window

323
 Design of FIR filters: Windows
• Simplest way of designing FIR filters
• Method is all discrete-time no continuous-time involved
• Start with ideal frequency response
  j
 

 
 j n  j j n

  1

Hd e  
n  
hd n e

2   
H e e h d n 
d
d 

• Choose ideal frequency response as desired response

• Most ideal impulse responses are of infinite length
• The easiest way to obtain a causal FIR filter from ideal is
h n  0 n  M
h n    d

 0 else
• More generally
1 0 n  M
h n   h n w n  where w n  

0 else
d

324
 Properties of Windows
• Prefer windows that concentrate around DC in frequency
– Less smearing, closer approximation
• Prefer window that has minimal span in time
– Less coefficient in designed filter, computationally efficient
• So we want concentration in time and in frequency
– Contradictory requirements
• Example: Rectangular window
sin  M  1  / 2 
  j

M
1  e
  j n  j M  1

 e
 j M / 2

W e 

n0
e 1 e
 j
sin  / 2 

325
 Windowing distortion
• increasing window length generally reduces the
width of the main lobe
• peak of sidelobes is generally independent of M

326
 Windows
Commonly used windows

• Rectangular 1 n  N  1
1  2 n
• Bartlett N
 2  n 
2

• Hanning 1  cos  

• Hamming  N  2  n 
0 .54  0 .46 cos  
•  N 
 2 n   4 n 
• Blackman 0 .42  0 .5 cos    0 .08 cos  
 N   N 
•  2 
 2 n 
• Kaiser J 0   1     J0 ( )

  N  1 



327
 Rectangular Window
• Narrowest main lob
– 4/(M+1)
– Sharpest transitions at
discontinuities in frequency

• Large side lobs

– -13 dB
– Large oscillation around
discontinuities

• Simplest window possible

1 0 n M
w n   
0 else

328
 Bartlett (Triangular) Window
• Medium main lob
– 8/M

• Side lobs
– -25 dB

• Hamming window
performs better
• Simple equation
 2n / M 0 n  M /2

w n   2  2 n / M M /2  n  M
 else
 0

329
 Hanning Window
• Medium main lob
– 8/M

• Side lobs
– -31 dB

• Hamming window performs

better

• Same complexity as
Hamming
1   2 n 
w n    1  cos    0 n  M
2   M  

 0 else

330
 Hamming Window
• Medium main lob
– 8/M

• Good side lobs

– -41 dB

• Simpler than Blackman 


  2 n  
w n   0 . 54  0 . 46 cos   0 n M
  M 

 0 else

331
 Blackman Window
• Large main lob
– 12/M

• Very good side lobs

– -57 dB

• Complex equation
  2 n   4  n 
 0 . 42  0 . 5 cos    0 . 08 cos   0 n M
 
w n    M   M 

 0 else

332
 Kaiser Window Filter Design Method
• Parameterized equation
forming a set of windows
– Parameter to change main-lob
width and side-lob area trade-off
  2 
n M / 2  
I  1    
w n   0
 M /2  
   0 n  M
I  

0

 0 else

– I0(.) represents zeroth-order

modified Bessel function of 1st
kind

33
Comparison of windows

334
 Kaiser window
• Kaiser window
β Transition Min. stop
width (Hz) attn dB
2.12 1.5/N 30
4.54 2.9/N 50
6.76 4.3/N 70
8.96 5.7/N 90
335
 Example
• Lowpass filter of length 51 and c   /2
Lowpass Filter Designed Using Hann window Lowpass Filter Designed Using Hamming window
0 0

-50 -50

-100 -100

0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1

/ /
Lowpass Filter Designed Using Blackman window
0

-50

-100

0 0.2 0.4 0.6 0.8 1 336

/
 Frequency Sampling Method
• In this approach we are given H ( k ) and
need to find H ( z )
• This is an interpolation problem and the
solution is given in the DFT part of the
course N
1 N 1 1 z
H ( z )   H ( k ).

2 
N k 0 j k
1  e N
. z 1

• It has similar problems to the windowing

approach 337
 FIR Digital Filter Order Estimation

Kaiser’s Formula:
 20 log 10
(  p s )  13
N  1
14 .6 (s   p ) / 2 


• ie N is inversely proportional to transition

band width and not on transition band
location

338
 UNIT-5
Multirate signal processing &
Finite Word length Effects

339
Single vs Multirate Processing

340
Basic Multirate operations: Decimation
and Interpolation

341
M-fold Decimator

342
 Sampling Rate Reduction by an Integer Factor:
Downsampling
• We reduce the sampling rate of a sequence by “sampling” it
x n   x nM   x nMT 
d c


• This is accomplished with a sampling rate compressor

• We obtain xd[n] that is identical to what we would get by

reconstructing the signal and resampling it with T’=MT
• There will be no aliasing if
 
  N
T' MT

343
 Frequency Domain Representation of Downsampling
• Recall the DTFT of x[n]=xc(nT)   2  k  
1 
X e j   X  j   
T k  
c
  T T  
• The DTFT of the downsampled signal can similarly written as
1    
X e  
j
2 r  
 1     2  r  
 
 
d


X c  j     X c  j   
T' r     T' T'  MT r     MT MT  

• Let’s represent the summation index as

r  i  kM where -   k   and 0  i M
   2 k 2  i   
1 M 1  1 
X e 
j
    X  j
c
    
d


M i  0  T r     MT T MT   
• And finally


1 M 1
   2 i 

e  
j 
X j
 Xe 
M M 



d M i 0

 

344

Frequency Domain Representation of Downsampling

345
Aliasing

346
Frequency Domain Representation of Downsampling w/ Prefilter

347
Decimation filter

348
L-fold Interpolator

349
 Increasing the Sampling Rate by an Integer Factor:
Upsampling
• We increase the sampling rate of a sequence interpolating it
x i n   x n / L   x nT
c
/ L 

• This is accomplished with a sampling rate expander

• We obtain xi[n] that is identical to what we would get by

reconstructing the signal and resampling it with T’=T/L
• Upsampling consists of two steps
– Expanding
 x n / L  n  0 ,∓ L ,∓ 2L ,...   
n  


xe    x k  n  kL
 0 else k  

– Interpolating

350
 Frequency Domain Representation of Expander
• The DTFT of xe[n]can be written as
  j     j n

   j Lk
 j L

Xe e     x k  n  kL e   x k e X e
n    k    k  

• The output of the expander is frequency-scaled

351
Input-output relation on the Spectrum

352
Periodicity and spectrum images

353
 Frequency Domain Representation of Interpolator
• The DTFT of the desired interpolated signals is

• The extrapolator output is given as

• To get interpolated signal we apply the following LPF

354
Interpolation filters

355
Fractional sampling rate convertor

356
Fractional sampling rate convertor

357
 Changing the Sampling Rate by Non-Integer Factor

• Combine decimation and interpolation for non-integer factors

Interpolator Decimator

x[n] xe[n] Lowpass filter xi[n] Lowpass filter xo[n] xd[n]

L Gain = L Gain = 1 M
Cutoff = p/L Cutoff = p/M

T T/L T/L T/L TM/L

• The two low-pass filters can be combined into a single one

Lowpass filter
x[n] xe[n] Gain = L xo[n] xd[n]
L Cutoff = M
min(p/L, p/M)

T T/L T/L TM/L 358

 Time Domain
• xi[n] in a low-pass filtered version of x[n]
• The low-pass filter impulse response is
h n   sin  n / L 
i
n / L
• Hence the interpolated signal is written as
     sin  n  kL  / L 
xin 

k  
x k
 n  kL /L
• Note that h i0   1

h in   0 n  ∓ L, ∓ 2L,...

• Therefore the filter output can be written as

x n   x n / L   x nT /L   x nT ' for n  0, ∓ L, ∓ 2L,...
i c c

359
Sampling of bandpass signals

360
Sampling of bandpass signals

361
Over sampling -ADC

362
363
364
365
Sub band coding

366
Sub band coding

367
Digital filter banks

368
Finite Word length Effects

369
 Finite Wordlength Effects
• Finite register lengths and A/D converters
cause errors in:-
(i) Input quantisation.
(ii) Coefficient (or multiplier)
quantisation
(iii) Products of multiplication truncated
or rounded due to machine length

370
 Finite Wordlength Effects
• Quantisation
Output
eo ( k )

Q
ei ( k )
Input

Q Q
  ei , o ( k ) 
2 2

371
 Finite Wordlength Effects
• The pdf for e using rounding
1
Q

Q Q
 2 2 p ( e ). de  E { e 2 }

Q 2
• Noise power 
2
  e
2

or Q 2
2
Q
 2

12

372
 Finite Wordlength Effects
• Let input signal be sinusoidal of unity
amplitude. Then total signal power P  1
2

• If b bits used for binary then Q  2 2b

2b
so that  
2 2 3

• Hence  2  3  2b
P   .2
2

or SNR  1 .8  6 b dB

373
 Finite Wordlength Effects
• Consider a simple example of finite
precision on the coefficients a,b of second
order system with poles  e  j



1
H ( z ) 
1 2
1  az  bz
 1
H ( z ) 
1  2  cos  . z 1   2 . z  2

• where a  2  cos  b   2

374
 Finite Wordlength Effects
bit pattern 2  cos  ,  2 

000 0 0
• 001 0.125 0.354
010 0.25 0.5
011 0.375 0.611
100 0.5 0.707
101 0.625 0.791
110 0.75 0.866
111 0.875 0.935
1.0 1.0 1.0
375
 Finite Wordlength Effects
• Finite wordlength computations

INPUT OUTPU
T
+

376
 Limit-cycles; "Effective Pole"
Model; Deadband

• Observe that for H (z)  1

1 2
(1  b1 z  b2 z )

• instability occurs when b 2  1

• i.e. poles are
• (i) either on unit circle when complex
• (ii) or one real pole is outside unit
circle.
• Instability under the "effective pole" model
is considered as follows
377
 Finite Wordlength Effects
• In the time domain with H (z) 
Y (z)
X (z)
•
y ( n )  x ( n )  b1 y ( n  1)  b 2 y ( n  2 )

• With b 2  1 for instability we have

Q b 2 y ( n  2 )  indistinguishable from y ( n  2 )
• Where Q   is quantisation

378
 Finite Wordlength Effects
• With rounding, therefore we have
b 2 y ( n  2 )  0 .5 y (n  2)
are indistinguishable (for integers)
or b 2 y ( n  2 )  0 .5  y ( n  2 )
• Hence  0 .5
y ( n  2 ) 
1  b2

• With both positive and negative numbers

 0 .5
y (n  2) 
1  b2
379
 Finite Wordlength Effects
 0 .5
• The range of integers
1  b2

constitutes a set of integers that cannot be

individually distinguished as separate or from the
asymptotic system behaviour.
• The band of integers  0 .5 0 .5 
 ,  
 1  b2 1  b2 
is known as the "deadband".
• In the second order system, under rounding, the
output assumes a cyclic set of values of the
deadband. This is a limit-cycle.
380
 Finite Wordlength Effects
• Consider the transfer function
G (z)  1
1 2
(1  b1 z  b2 z )

y k  x k  b1 y k  1  b 2 y k  2

• if poles are complex then impulse response

is given by h k

k

hk  . sin ( k  1) 
sin 

381
 Finite Wordlength Effects
  cos 1   
• Where   b  b1 
2  2 b 2 
• If b 2  1 then the response is sinusiodal
with frequency

1  b
cos 1
 1 
 
T  2 


• Thus product quantisation causes instability

implying an "effective “ b 2  1 .

382
 Finite Wordlength Effects

• Notice that with infinite precision the

response converges to the origin

• With finite precision the reponse does not

converge to the origin but assumes
cyclically a set of values –the Limit Cycle

383
 Finite Wordlength Effects
• Assume e1 ( k ) ,e 2 ( k ) ….. are not
correlated, random processes etc.

2
2
   e
2
 hi 2 ( k ) Q
e
2
0i
k 0 
12
Hence total output noise power
2b
 2   2
  2
 2. 2

2k sin 2
( k 1) 
0 01 02   .
12 k 0 sin 2



b
• Where Q  2 and
k sin ( k  1) 
h1 ( k )  h 2 ( k )   . ; k  0
sin 
384
 Finite Wordlength Effects

• ie
2b
 2  2 1   2
. 1 
0  2 
6 1   1  4
 2 2 cos 2 

385
 Finite Wordlength Effects






B(n+1)
• For FFT A(n)

B(n) -
B(n+1)
W(n)

A ( n  1)  A ( n )  W ( n ). B ( n )
B ( n  1)  A ( n )  W ( n ). B ( n )

A(n)
A(n+1)
B(n+1)
B(n) B(n)W(n)

386
 Finite Wordlength Effects
• FFT
2 2
A ( n  1)  B ( n  1)  2

2 2
A ( n  1)  2 A(n)

A(n )  2 A(n)

• AVERAGE GROWTH: 1/2 BIT/PASS

387
 Finite Wordlength Effects
IMAG 1.0
• FFT
-1.0 1.0
REAL

-1.0
A x ( n  1)  A x ( n )  B x ( n ) C ( n )  B y (n )S (n )

A x ( n  1)  A x ( n )  B x ( n ) C ( n )  B y (n) S (n )

A x ( n  1)
 1 . 0  C ( n )  S ( n )  2 . 414 ....
Ax (n )

• PEAK GROWTH: 1.21.. BITS/PASS

388
 Finite Wordlength Effects
• Linear modelling of product quantisation
x(n) ~
x (n)
Q 

• Modelled as
x(n) + ~
x (n)  x(n)  q (n)

q(n)

389
 Finite Wordlength Effects
• For rounding operations q(n) is uniform
distributed between Q
 2 , and where Q is
Q
2

the quantisation step (i.e. in a wordlength of

bits with sign magnitude representation or
mod 2, Q  2 ).b

• A discrete-time system with quantisation at

the output of each multiplier may be
considered as a multi-input linear system

390
 Finite Wordlength Effects

q1 ( n )... q 2 ( n )... q p ( n )

 x ( n ) h(n)  y ( n )
• Then 
 p
  r) 
y ( n )   x ( r ). h ( n  r )    q  ( r ). h  ( n  
r0   1  r  0 

• where h ( n ) is the impulse response of the

system from  the output of the multiplier
to y(n).

391
 Finite Wordlength Effects
• For zero input i.e. x ( n )  0 ,  n we can write
p 
y ( n )   qˆ .  h ( n  r )
 1 r0

• where qˆ is the maximum of q ( r ) ,   , r

which is not more than Q
 2
Q p
 
 .  r) 

• ie y (n )  
2   1 n  0 h ( n 





392
 Finite Wordlength Effects 


• However
 
 h ( n )   h ( n )
n0 n0

• And hence
pQ 
y ( n )  .  h ( n )
2 n0

• ie we can estimate the maximum swing at

the output from the system parameters and
quantisation level
393
Finite Precision Numerical
Effects

394
 Quantization in Implementing Systems
• Consider the following system

• A more realistic model would be

• In order to analyze it we would prefer

395
 Effects of Coefficient Quantization in IIR Systems

• When the parameters of a rational system are quantized

– The poles and zeros of the system function move
• If the system structure of the system is sensitive to
perturbation of coefficients
– The resulting system may no longer be stable
– The resulting system may no longer meet the original specs
• We need to do a detailed sensitivity analysis
– Quantize the coefficients and analyze frequency response
– Compare frequency response to original response
• We would like to have a general sense of the effect of
quantization

396
 M
Effects on Roots M

 bkz k
Quantiza  b̂ k z  k

H z   k 0
N Ĥ  z   k 0
N

 a tion 1   â
k
1  k z k
z k
k 1 k 1

• Each root is affected by quantization errors in ALL coefficient

• Tightly clustered roots can be significantly effected
– Narrow-bandwidth lowpass or bandpass filters can be very
sensitive to quantization noise
• The larger the number of roots in a cluster the more sensitive it
becomes
• This is the reason why second order cascade structures are less
sensitive to quantization error than higher order system
– Each second order system is independent from each other

397
 Poles of Quantized Second-Order Sections
• Consider a 2nd order system with complex-conjugate pole pair

• The pole locations after quantization will be on the grid point

 3-bits

7-bits 

398
 Coupled-Form Implementation of Complex-Conjugate Pair

• Equivalent implementation of the

second order system

• But the quantization grid this time is

399
 Effects of Coefficient Quantization in FIR Systems
• No poles to worry about only zeros
• Direct form is commonly used for FIR systems
M

H z    h n z  n
n0

• Suppose the coefficients are quantized M

M

Ĥ z    ĥ n z  n  H  z    H z   H z     h  n  z  n
n0
n0

• Quantized system is linearly related to the quantization error

• Again quantization noise is higher for clustered zeros

• However, most FIR filters have spread zeros

400
 Round-Off Noise in Digital Filters
• Difference equations
implemented with finite-
precision arithmetic are
non-linear systems
• Second order direct form I
system
• Model with quantization
effect
• Density function error
terms for rounding

401
 Analysis of Quantization Error
• Combine all error terms to single location to get

e n   e 0 n   e 1 n 

 e 2 n   e 3n   e 4n 



2B
2
• The variance of e[n] in the general case is  2
e  M  1  N 
12
N

• The contribution of e[n] to the output is f n    a f n

k
 k   e n 
k 1

• The variance of the output error term f[n] is

2B
n  H ef z   1 / A z 
2
M  1  N 2 
 2
f
  h
ef
12 n  

402
 Round-Off Noise in a First-Order System
• Suppose we want to implement the following stable system
b
H z   a  1
1
1  az
• The quantization error noise variance is
2B  2 2B  2B  
2  M  1  N 2
 h n   2 2  2 2  1 

2n
a




f ef
12 12 2

12 n   n0
1  a 
• Noise variance increases as |a| gets closer to the unit circle
• As |a| gets closer to 1 we have to use more bits to compensate for the
increasing error

403
 Zero-Input Limit Cycles in Fixed-Point Realization of IIR Filters

• For stable IIR systems the output will decay to zero when the input
becomes zero
• A finite-precision implementation, however, may continue to oscillate
indefinitely
• Nonlinear behaviour very difficult to analyze so we sill study by example
• Example: Limit Cycle Behavior in First-Order Systems
y n   ay n  1   x n  a  1

• Assume x[n] and y[n-1]

are implemented by 4 bit register

404
 Example Cont’d
y n   ay n  1   x n  a  1

• Assume that a=1/2=0.100b and the input is

x n  
7
 n   0 . 111 b  n 
8



• If we calculate the output for values of n

n y[n] Q(y[n])
0 7/8=0.111b 7/8=0.111b
1 7/16=0.011100b 1/2=0.100b
2 1/4=0.010000b 1/4=0.010b
3 1/8=0.001000b 1/8=0.001b
4 1/16=0.00010b 1/8=0.001b

• A finite input caused an oscillation with period 1

405
 Example: Limit Cycles due to Overflow
• Consider a second-order system realized by
ŷ n   x n   Q a 1 ŷ n  1   Q a 2 ŷ n  2 
– Where Q() represents two’s complement rounding
– Word length is chosen to be 4 bits
• Assume a1=3/4=0.110b and a2=-3/4=1.010b
• Also assume
ŷ  1   3 / 4  0 . 110 b and ŷ  2    3 / 4  1 . 010 b

• The output at sample n=0 is

ŷ 0   0 . 110 b  0.110b  1 . 010 b  1.010b
 0.100100b  0.100100b
• After rounding up we get
ŷ 0   0.101b  0.101b  1.010b  -3/4

• Binary carry overflows into the sign bit changing the sign
• When repeated for n=1
ŷ 0   1.010b  1.010b  0 . 110  3 / 4

406
 Avoiding Limit Cycles
• Desirable to get zero output for zero input: Avoid limit-cycles
• Generally adding more bits would avoid overflow
• Using double-length accumulators at addition points would
decrease likelihood of limit cycles
• Trade-off between limit-cycle avoidance and complexity
• FIR systems cannot support zero-input limit cycles

407