MATLAB: Has the algorithm for fdesign changed in Signal Processing Toolbox 6.13 (R2010a) as compared to the previous versions

fdesignSignal Processing Toolbox

I am using fdesign.decimator to design an equiripple filter with given specifications. I have recently upgraded from R2007b to R2010a. In R2007b, this code output a filter with 89 coefficients, but the same code in R2010a results in a filter with 127 taps. This is breaking compatibility for my downstream code as I am using the filter and the signal output from the filter as a reference to compare with my filter code written in assembly. Compare the output generated by the following code in R2007a and R2010a:

Fpass = 3400;   % Passband Frequency

Fstop = 4000; % Stopband Frequency Apass = 0.01; % Passband Ripple (dB) Astop = 80; % Stopband Attenuation (dB) Fs = 16000; % Sampling Frequency h = fdesign.decimator(2, 'Lowpass', 'fp,fst,ap,ast', Fpass, Fstop, … Apass, Astop, Fs); Hd = design(h, 'equiripple', … 'FilterStructure', 'firdecim', … 'MinPhase', true, … 'MinOrder', 'any', … 'StopbandShape', 'flat');

Best Answer

This is expected behavior. The filter designed in R2007b does not match the passband specifications very well. In Signal Processing Toolbox 6.11 (R2009a) , there was an improvement in the design of minimum order equiripple lowpass filters which involved measuring the filter response to ensure that the filter order meets the design specifications. This is causing the increment in the filter order in R2009a and later versions. As a workaround, design the filter with non-minimum order specifications thus forcing the same number of coefficients as in R2007b. This however changes the filter specification and thus will return different values for the coefficients.

Related Solutions

MATLAB: Has the algorithm for fdesign changed in Signal Processing Toolbox 6.13 (R2010a) as compared to the previous versions

This is expected behavior. In Signal Processing Toolbox 6.11 (R2009a) , there was an improvement in the design of minimum order equiripple lowpass filters which involved measuring the filter response to ensure that the filter order meets the design specifications. This is causing the increment in the filter order in R2009a and later versions.

As a workaround, design the filter with non-minimum order specifications thus forcing the same number of coefficients as in R2007b. However, this changes the filter specification and will return different values for the coefficients.

MATLAB: Generate C code from a filter which created by fdatool or filterbuilder

Look under File -> Generate Matlab Code -> Filter Design Function. If I choose that option and then make some adjustments in the output file (code generation does not support dfilt objects or fdesign as you discovered)

function b = myfilter %#codegen
%MYFILTER Returns a discrete-time filter object.
%


% MATLAB Code
% Generated by MATLAB(R) 7.14 and the Signal Processing Toolbox 6.17.
%
% Generated on: 01-Jul-2012 05:38:14
%
% Equiripple Lowpass filter designed using the FIRPM function.
% All frequency values are normalized to 1.
Fpass = 0.45;            % Passband Frequency
Fstop = 0.55;            % Stopband Frequency
Dpass = 0.057501127785;  % Passband Ripple
Dstop = 0.001;           % Stopband Attenuation
dens  = 20;              % Density Factor
% Calculate the order from the parameters using FIRPMORD.
[N, Fo, Ao, W] = firpmord([Fpass, Fstop], [1 0], [Dpass, Dstop]);
% Calculate the coefficients using the FIRPM function.
  b  = firpm(N, Fo, Ao, W);

I am able to generate C code

Best Answer

Related Solutions

MATLAB: Has the algorithm for fdesign changed in Signal Processing Toolbox 6.13 (R2010a) as compared to the previous versions

MATLAB: Generate C code from a filter which created by fdatool or filterbuilder

Related Question