Multirate Signal Processing PDF

Module 5 - Multirate Signal Processing
Patrick A. Naylor
Digital Signal Processing p.1/25
Contents
Applications of multirate signal processing Fundamentals decimation interpolation Resampling by rational fractions Multirate identities Polyphase representations Maximally decimated lter banks aliasing amplitude and phase distortion perfect reconstruction conditions
Introduction
In single-rate DSP systems, all data is sampled at the same rate no change of rate within the system. In multirate DSP systems, sample rates are changed (or are different) within the system Multirate can offer several advantages reduced computational complexity reduced transmission data rate.
Example: Audio sample rate conversion

recording studios use 192 kHz CD uses 44.1 kHz wideband speech coding using 16 kHz master from studio must be rate-converted by a factor 44.1 192
Example: Oversampling ADC

Consider a Nyquist rate ADC in which the signal is sampled at the desired precision and at a rate such that Nyquists sampling criterion is just satised. Bandwidth for audio is 20 Hz < f < 20 kHz Antialiasing lter required has very demanding specication |H (j )| = 0 dB, f < 20 kHz |H (j )| < 96 dB, f 44.1 kHz 2
Requires high order analogue lter such as elliptic lters that have very nonlinear phase characteristics hard to design, expensive and bad for audio quality.
Nyquist Rate Conversion Anti-aliasing Filter.
Consider oversampling the signal at, say, 64 times the Nyquist rate but with lower precision. Then use multirate techniques to convert sample rate back to 44.1 kHz with full precision. New (over-sampled) sampling rate is 44.1 64 kHz. Requires simple antialiasing lter |H (j )| = 0 dB, f < 20 kHz |H (j )| < 96 dB, f (44.1 64) 44.1 kHz 2
Could be implemented by simple lter (eg. RC network) Recover desired sampling rate by downsampling process.
Oversampled Conversion Antialiasing Filter
Overall System
This is a simplied version In these lectures we will study blocks like G(z ) and 64
Example: Subband Coding

Consider quantizing the samples of a speech signal. How many bits are required? In general, 16 bits precision per sample is normally used for audio. This gives an adequate dynamic range. In practice, certain frequency bands are less important perceptually because they contain less signicant information bands with less information or lower perceptual importance may be quantized with lower precision - fewer bits. Divide the spectrum of the signal into several subbands then allocate bits to each band appropriately.
16 bits per sample, 10 kHz sampling frequency gives 160 kbits/s Divide into 2 bands: high frequency and low frequency subbands. High frequencies of speech are less important to intelligibility. Therefore use only 8 bits per sample The sampling frequency can be reduced by a factor of 2 since bandwidth is halved, still satisfying Nyquist criterion. 5 16 + 5 8 = 120 kbits/s 4:3 compression Reconstructed signal has no noticeable reduction is signal quality.
Fundamental Multirate Operations

Downsampling by a factor M lter and M-fold decimator
Upsampling by a factor L L-fold expander and lter
M-fold Decimator
For an input sequence x(n), select only the samples which occur at integer multiples of M . The other samples are thrown away.
yD (n) = x(M n) Aliasing will occur in yD (n) unless x(n) is sufciently bandlimited loss of information.
Eg. M = 2
L-fold Expander
For an input sequence x(n), insert L 1 zeros between each sample.
yE (n) = x(M n) x(n) can always be recovered from yE (n) no loss of information, no aliasing.
Eg. L = 2
Frequency Domain View of the Expander

From the denition of the z-transform
YE (z ) =
n=
yE (n)z n yE (kL)z kL
k=
=
k=
x(k )z kL = X (z L )
For frequency response write z = ej giving YE (ej ) = X (ejL ) YE is a compressed version of X Multiple images of X (ej ) are created in YE (ej ) between = 0 and = 2
To use the expander for interpolation, a lowpass lter is applied after the expander to remove the images (shaded).
Frequency Domain View of the Decimator

From the denition of the z-transform

YD (z ) =
n=
yD (n)z n =
n=
x(M n)z n
Let x1 (n) = Then

x(n) 0
if n is an integer multiple of M otherwise
YD (z ) =
n=
x1 (M n)z n =
k=
x1 (k )z k/M
since x1 (k ) = 0 unless k is a multiple of M .
Therefore 1 1/M YD (z ) = X1 (z )= M
M 1 k X z 1/M WM k=0
as will be shown on the next slide and using k WM = ej 2k/M
For frequency response write z = ej to give 1 YD (ej ) = M

M 1
X ej (2k)/M
k=0
We arrive at the previous expression for YD (z ) by considering a new sequence cM (n) = and then writing x1 (n) = cM (n)x(n) Further consideration of cM (n) tells us that cM (n) is the inverse Fourier transform of unity and can be written 1 cM (n) = M
M 1 kn WM k=0
1 if n is an integer multiple of M 0 otherwise
Then 1 X 1 (z ) = M 1 = M 1 = M
M 1 kn n x(n)WM z k=0 n= M 1
x(n)
k=0 n= M 1 k X zWM k=0
n k WM z
from the denition of the z-transform. So nally 1 YD (z ) = M

M 1 k X z 1/M WM k=0
What does YD (z ) =
1 M
M 1 k=0
k X z 1/M WM represent?
stretching of X (ej ) to X (ej/M ) creating M 1 copies of the stretched versions shifting each copy by successive multiples of 2 and superimposing (adding) all the shifted copies dividing the result by M
To use a decimation process we must rst bandlimit the signal to | | < M .

Summary
Downsampling
Upsampling

Multirate Signal Processing PDF

Uploaded by

Copyright:

Available Formats

Multirate Signal Processing PDF

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Multirate Signal Processing PDF

Uploaded by

Copyright:

Available Formats

What are some applications of multirate signal processing?

What are some applications of multirate signal processing?

What are the advantages of using multirate signal processing?

What are the advantages of using multirate signal processing?

Module 5 - Multirate Signal Processing

Digital Signal Processing p.1/25

Digital Signal Processing p.2/25

Digital Signal Processing p.3/25

Example: Audio sample rate conversion

Digital Signal Processing p.4/25

Example: Oversampling ADC

Digital Signal Processing p.5/25

Nyquist Rate Conversion Anti-aliasing Filter.

Digital Signal Processing p.6/25

Digital Signal Processing p.7/25

Oversampled Conversion Antialiasing Filter

Digital Signal Processing p.8/25

Example: Subband Coding

Digital Signal Processing p.10/25

Digital Signal Processing p.11/25

Fundamental Multirate Operations

Upsampling by a factor L L-fold expander and lter

Digital Signal Processing p.12/25

Digital Signal Processing p.13/25

Digital Signal Processing p.14/25

Digital Signal Processing p.15/25

Digital Signal Processing p.16/25

Frequency Domain View of the Expander

Frequency Domain View of the Decimator

Let x1 (n) = Then

if n is an integer multiple of M otherwise

since x1 (k ) = 0 unless k is a multiple of M .

Digital Signal Processing p.19/25

as will be shown on the next slide and using k WM = ej 2k/M

For frequency response write z = ej to give 1 YD (ej ) = M

Digital Signal Processing p.20/25

1 if n is an integer multiple of M 0 otherwise

Digital Signal Processing p.21/25

from the denition of the z-transform. So nally 1 YD (z ) = M

Digital Signal Processing p.22/25

Digital Signal Processing p.23/25

To use a decimation process we must rst bandlimit the signal to | | < M .

Digital Signal Processing p.25/25

You might also like