IIR wave digital filter decimator
mfilt.iirwdfdecim will be removed in a
future release. Use
hm = mfilt.iirwdfdecim(c1,c2,...)
hm = mfilt.iirwdfdecim(c1,c2,...) constructs
an IIR wave digital decimator given the coefficients specified in
the cell arrays c1, c2, and so on. The IIR decimator
a polyphase IIR filter where each phase is a cascade wave digital
allpass IIR filter.
Each cell array
ci contains a set of vectors
representing a cascade of allpass sections. Each element in one cell
array is one section. For more information about the contents of each
cell array, refer to
The contents of the cell arrays are the same for both filter constructors
mfilt.iirwdfdecim interprets them same way
The following exception applies to interpreting the contents
of a cell array — if one of the cell arrays
only one vector, and that vector comprises a series of 0s and one
element equal to 1, that cell array represents a
dfilt.delay section with latency equal
to the number of zeros, rather than a
This exception occurs with quasi-linear phase IIR decimators.
Usually you do not construct IIR wave digital filter decimators explicitly. Instead, you obtain an IIR wave digital filter decimator as a result of designing a halfband decimator. The first example in the following section illustrates this case.
Design an elliptic halfband decimator with a decimation factor
equal to 2. Both examples use the
structure (an input argument to the
to design the final decimator.
The first portion of this example generates a filter specification
d that stores the specifications for the
tw = 100; % Transition width of filter to design, 100 Hz. ast = 80; % Stopband attenuation of filter 80 dB. fs = 2000; % Sampling frequency of the input signal. m = 2; % Decimation factor. d = fdesign.decimator(m,'halfband','tw,ast',tw,ast,fs);
Now perform the actual design using
hm is an
Hm = design(d,'ellip','FilterStructure','iirwdfdecim'); % Note that realizemdl requires Simulink realizemdl(hm) % Build model of the filter.
Design a linear phase halfband decimator for decimating a signal by a factor of 2.
tw = 100; % Transition width of filter, 100 Hz. ast = 60; % Filter stopband attenuation = 80 dB fs = 2000; % Input signal sampling frequency. m = 2; % Decimation factor. d = fdesign.decimator(m,'halfband','tw,ast',tw,ast,fs);
d to design the final filter
hm = design(d,'iirlinphase','filterstructure',... 'iirwdfdecim'); % Note that realizemdl requires Simulink realizemdl(hm) % Build model of the filter.
The models that
for each example appear below. At this level, the realizations of
the filters are identical. The differences appear in the subsystem
blocks Phase1 and Phase2.
This is the Phase1 subsystem from the halfband model.
Phase1 subsystem from the linear phase model is less revealing—an allpass filter.