EIE3510Assignment 3 - Final
EIE3510Assignment 3 - Final
Assignment 3
1. (25 points) You are asked to go through the process of deriving a digital IIR filter from a
Butterworth low-pass filter, and examine the frequency responses of both filters.
(a) Derive the transfer function Hc (s) of a 4th-order Butterworth low-pass filter with cut-off
frequency of 50 Hz.
(b) Find out all poles of Hc (s).
(c) Assuming sampling frequency of 400 Hz, determine the corresponding poles of IIR filter H(z)
using the impulse invariance approach.
(d) Derive H(z) with pass-band ripple being equal to 0 dB.
(e) Plot the magnitude and phase responses of Hc (s) and H(z).
2. (20 points) The system function of a discrete-time system is
2 1
H(z) = − .
1− e−0.2 z −1 1− e−0.4 z −1
(a) Assume that this discrete-time filter was designed by the impulse invariance method with
Td = 2; i.e., h[n] = 2hc (2n), where hc (t) is real. Find the system function Hc (s) of a continuous-time
filter that could have been the basis for the design. Is your answer unique? If not, find another
system function Hc (s)
(b) Assume that H(z) was obtained by the bilinear transform method with Td = 2. Find the
system function Hc (s) that could have been the basis for the design. Is your answer unique? If not,
find another Hc (s).
3. (20 points, 2020Spring final exam) The figure below shows the magnitude response of a discrete-
time Butterworth filter. It is observed that the cutoff frequency ωc = 0.2π.
(a) What is the order of this filter?
2
(b) Suppose that the pole-zero of the magnitude-squared response function |Hc (jΩ)| of the analog
Butterworth filter is given in the following diagram (see Reference on Page 2). We use such an analog
prototype to design a discrete-time Butterworth filter H(z) with cutoff frequency ωc = 0.2π. Suppose
impulse invariance mapping is used for the CT to DT conversion. Write down the poles of H(z).
(c) If we want to apply this filter to process some audio signal such that the continuous-time cutoff
frequency is 4kHz. What is the corresponding sampling rate?
(d) If we further require the filter has 20dB attenuation at 5kHz. Specify the order N .
4. (15 points, please also submit the code and plots)
You are required to use MATLAB functions to obtain IIR and FIR filters that satisfy the following
conditions:(Hint: you can use function fir1(), cheby1(), ellip())
- Passband: 0.3π to 0.6π, Stop band 0.26π, 0.64π
- Pasband ripple: < 0.5 dB
1
- Stopband attenuation: > 25 dB
- Delay: as low as possible
(a) Design a Chebyshev IIR filter and an Elliptic IIR filter. Plot their magnitude, phase, and
group delay responses;
(b) Design a linear-phase FIR filter using the windowing or the optimisation approach. Plot the
filter’s magnitude, phase, and group delay responses;
(c) Compare the three filters (Chebyshev IIR filter, Elliptic IIR filter,linear-phase FIR filter )in
terms of frequency response and computation complexity.
5. (20 points, please also submit the code and plots)
Use Matlab function to design a low-pass filter:
Passband: 0 − 0.4π, ripple < 0.01 Stopband: 0.6π − π, ripple < 0.001
You should design the above continuous time filter in forms of Butterworth, Chebyshev I, Cheby-
shev II, and Elliptic filters and plot the whole band and passband seperately.