Adsp 06 Multirate Processing
Adsp 06 Multirate Processing
Gerhard Schmidt
Christian-Albrechts-Universität zu Kiel
Faculty of Engineering
Institute of Electrical and Information Engineering
Digital Signal Processing and System Theory
Multi-Rate Digital Signal Processing
•Contents
❑ Introduction
❑ Digital processing of continuous-time signals
❑ Efficient FIR structures
❑ DFT and FFT
❑ Digital filters
❑ Multi-rate digital signal processing
❑ Decimation and interpolation
❑ Filters in sampling rate alteration systems
❑ Polyphase decomposition and efficient structures
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 2
Multi-Rate Digital Signal Processing
•Basic Ideas
Why multi-rate systems?
❑ In many practical signal processing applications different sampling rates are present, corresponding to different
bandwidths of the individual signals multi-rate systems.
❑ Often a signal has to be converted from one rate to another. This process is called sampling rate conversion.
❑ Sampling rate conversion can be carried out by analog means, that is D/A conversion followed by A/D conversion
using a different sampling rate
D/A converter introduces signal distortion, and the A/D converter leads to quantization effects.
❑ Sampling rate conversion can also be carried out completely in the digital domain:
Less signal distortions, more elegant and efficient approach.
❑ Topic of this chapter is multi-rate signal processing and sampling rate conversion in the digital domain.
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 3
Multi-Rate Digital Signal Processing
… Inserting the definition of the signal and exploiting that contains a lot of zeros ...
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 5
Multi-Rate Digital Signal Processing
With it follows
Inserting
the result
The z-transform can be obtained as
from above
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 6
Multi-Rate Digital Signal Processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 7
Multi-Rate Digital Signal Processing
Remark: Shifted versions of are weighted with the factor according to the last slide.
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 8
Multi-Rate Digital Signal Processing
Specification for the desired magnitude frequency response of the lowpass anti-aliasing (or decimation) filter:
where denotes the highest frequency that needs to be preserved in the decimated signal.
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 9
Multi-Rate Digital Signal Processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 10
Multi-Rate Digital Signal Processing
•Questions
Questions about sample rate reduction:
Partner work – Please think about the following questions and try to find answers
(first group discussions, afterwards broad discussion in the whole group).
❑ What happens in the spectral domain when you decimate (without filtering) the time-domain signal?
……………………………………………………………………………………………………………………………………………………………………………………..
……………………………………………………………………………………………………………………………………………………………………………………..
❑ Is an anti-aliasing filter always necessary? If not, what are the conditions for applying such a filter?
……………………………………………………………………………………………………………………………………………………………………………………..
……………………………………………………………………………………………………………………………………………………………………………………..
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 11
Multi-Rate Digital Signal Processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 12
Multi-Rate Digital Signal Processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 13
Multi-Rate Digital Signal Processing
We can see that each repeated spectrum is weighted with a complex exponential (rotation) factor (called “twiddle” factor).
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 14
Multi-Rate Digital Signal Processing
Notation: Since the upsampling factor is named with in conformance with the majority of the technical literature in the
following we will denote the length for an FIR filter with .
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 15
Multi-Rate Digital Signal Processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 16
Multi-Rate Digital Signal Processing
The frequency response of does not change by upsampling, however the frequency axis is scaled differently.
The new sampling frequency is now (in terns of for the lower sampling rate) equal to
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 17
Multi-Rate Digital Signal Processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 18
Multi-Rate Digital Signal Processing
where denotes the sampling period. If we instead sample at a much higher rate we have
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 19
Multi-Rate Digital Signal Processing
If is passed through an ideal lowpass filter with cut-off frequency and a gain of , the output of the filter
will be precisely .
Therefore, we can now state our specifications for the lowpass interpolation filter:
Where denotes the highest frequency that needs to be preserved in the interpolated signal
(related to the lower sampling frequency).
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 20
Multi-Rate Digital Signal Processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 21
Multi-Rate Digital Signal Processing
Input-output relation?
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 22
Multi-Rate Digital Signal Processing
With it follows
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 23
Multi-Rate Digital Signal Processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 24
Multi-Rate Digital Signal Processing
•Questions
Motivation of multi-rate structures
Partner work – Please think about the following questions and try to find answers
(first group discussions, afterwards broad discussion in the whole group).
❑ If you would like to convolve a signal at a sample rate of 10 kHz with an impulse response (FIR filter) of 10 seconds
length, how many multiplications and additions do you need per second?
……………………………………………………………………………………………………………………………………………………………………………………..
……………………………………………………………………………………………………………………………………………………………………………………..
❑ Assume that you can split the signal into 10 equally wide bandpass signals (assmuming that you have ideal filters that
are “for free”) and you can use the largest possible subsampling rate, how many multiplications and additions do you
need now (again per second)?
……………………………………………………………………………………………………………………………………………………………………………………..
……………………………………………………………………………………………………………………………………………………………………………………..
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 25
Multi-Rate Digital Signal Processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 26
Multi-Rate Digital Signal Processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 27
Multi-Rate Digital Signal Processing
with
Example for
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 28
Multi-Rate Digital Signal Processing
with
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 29
Multi-Rate Digital Signal Processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 30
Multi-Rate Digital Signal Processing
the input samples appear at the output of the system without any distortion for all .
All in-between samples are determined by interpolation.
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 31
Multi-Rate Digital Signal Processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 32
Multi-Rate Digital Signal Processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 33
Multi-Rate Digital Signal Processing
•Questions
Questions about sample filterbanks:
Partner work – Please think about the following question and try to find answers
(first group discussions, afterwards broad discussion in the whole group).
❑ Please try to derive the equation
by transforming the equation first to the Fourier domain and afterwards to the time domain.
……………………………………………………………………………………………………………………………………………………………………………………..
……………………………………………………………………………………………………………………………………………………………………………………..
……………………………………………………………………………………………………………………………………………………………………………………..
……………………………………………………………………………………………………………………………………………………………………………………..
……………………………………………………………………………………………………………………………………………………………………………………..
……………………………………………………………………………………………………………………………………………………………………………………..
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 34
Multi-Rate Digital Signal Processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 35
Multi-Rate Digital Signal Processing
and the downsampling as . Combining both equations we can write the decimation operation according to
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 36
Multi-Rate Digital Signal Processing
→ Multiplication of with and leads to the result and which are discarded in the decimation process
→ these compositions are not necessary.
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 37
Multi-Rate Digital Signal Processing
The output of the interpolation filter can be obtained as convolution with the length
→ The output sample is obtained by multiplication of with , where a lot of zero multiplications are involved,
which are inserted by upsampling operation.
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 39
Multi-Rate Digital Signal Processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 40
Multi-Rate Digital Signal Processing
denoting the downsampling factor and the z-transform for type-1 polyphase components
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 41
Multi-Rate Digital Signal Processing
Remark: The structure (b) has the same complexity as the direct form structure from the previous section, therefore no
further advantage. However, the polyphase structures are important for digital filter banks.
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 42
Multi-Rate Digital Signal Processing
denoting the upsampling factor, and the type-1 polyphase components of with .
Resulting interpolator structure :
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 44
Multi-Rate Digital Signal Processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 45
Multi-Rate Digital Signal Processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 47
Multi-Rate Digital Signal Processing
→ delay in one branch of the polyphase structure can be replaced with the delay
b. The factor can be shifted before the upsampler, and the factor behind the downsampler:
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 48
Multi-Rate Digital Signal Processing
d. In every branch we now have a decimator (marked with the dashed box), which can again be efficiently realized using
the polyphase structure from the previous section. Thus, each type-1 polyphase component is itself decomposed
again in polyphase components
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 49
Multi-Rate Digital Signal Processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 51
Multi-Rate Digital Signal Processing
Polyphase filters are calculated with 16 kHz sampling rate compared to 48 kHz sampling rate in the original structure.
Rate conversion from 32 kHz to 48 kHz: Exercise!
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 52
Multi-Rate Digital Signal Processing
•Summary – Part 1
❑ Introduction
❑ Digital processing of continuous-time signals
❑ DFT and FFT
❑ Digital filters
❑ Multi-rate digital signal processing
❑ Decimation and interpolation
❑ Filters in sampling rate alteration systems
❑ Polyphase decomposition and efficient structures
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 53
Multi-Rate Digital Signal Processing
•Summary – Part 2
❑ Introduction
❑ Digital processing of continuous-time signals
❑ Efficient FIR structures
❑ DFT and FFT
❑ Digital filters
❑ Multi-rate digital signal processing
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 54
Multi-Rate Digital Signal Processing
•Summary – Part 3
And finally:
Enjoy applying your new knowledge – in the upcoming lectures, during a lab, while working on your thesis and most importantly
during your profession as an engineer.
Digital Signal Processing and System Theory | Advanced Digital Signal Processing| Multi-Rate Digital Signal Processing Slide 55