BE2M31CZS / BE2M31DSPA - Exercise - Digital Filtering - FIR Filters
Tasks to do:
Manual design of FIR filter using "Window Method" - low-pass filter
Describe and explain principles of FIR filter design using "Window Method".
Realize step by step a design of low-pass (LP) filter with cutoff frequency
f_c = f_s / 4. Use FIR filter order of M = 30 (or 50,
80), i.e. the length of impulse should be M+1 = 31 (or 51,
1st checked result (1 point):
Display the particular results of a FIR filter design, i.e.:
the part of infinite impulse response of an ideal filter
computed on the basis of Fourier series (use n = -200 : 200),
shortened finite impulse response of the filter with the order
M = 30,
achieved frequency response of designed FIR filter of the
order M (the length of displayed twoside frequency response should
be 1000 samples and display both magnitude response in linear scaling
as well as magnitude response in dBs),
to compare with achieved frequency repsonse of designed filter, add to the same
figure the ideal frequency response with the length of N=1000 samples,
repeat and observe an influence on used weghting window (rectangular vs. Hamming).
Design of FIR filters using MATLAB implementation of window method (fcn fir1)
Design High-Pass FIR filter with the cutoff frequency fc = fs/4,
Use the order M=30 and Hamming window in the first step.
Observe achieved frequency and impulse response of designed filter
(fcn freqz, impz).
2nd checked result (1 point):
Repeat the above mentioned FIR-filter design for various orders
and types of weighting window and display:
magnitude frequency response in dBs for filters of orders M
= 10, 30, 50, 200 and default Hamming window,
magnitude frequency response in dBs for filters of the order M
= 30 and used rectangular, Hamming, and Blackamn weighting window.
3rd checked result (1 point): Compare magnitude frequency response in dBs for designed FIR
filter of the order M = 50 and used Hamming window with
magnitude responses of high-pass IIR filters with same cuttoff
frequency and the order M
= 6 (for all 4 aproximations - ripple in pass-band and
stop-band should be R_p = 1 dB and R_s = 50 dB respectively).
Design of bandpass filter
Design empirically bandpass filter for the band 300
< f < 3400, which can be used for the filtering of an acoustic signal (speech) into the telephone band. Design the filters for sampling frequencies 8kHz, 16kHz and 44,1kHz.
Find the order for which the suppression of the signal in stop band will be minimally 20 dB. Observe the orders required for different sampling frequencies, i.e. fs = 8, 16,
or 44.1 kHz respectively.
4th checked result (1 point):
Work with the speech signals available in
speech_8_16_44.mat (binary MATLAB-format file, it contains 3 signals saved in variables sig8, sig16
and sig44. To load them into MATLAB, use the command "load
observe frequency responses of designed filters for particular
smapling frequencies (compare them with IIR filters for
the same purpose designed at last seminar),
spectrograms of original and filtered signals,
in suitably selected zoomed part observe also a time-shift
between original and filtered signals .