MATLAB: Are settings saved incorrectly when I import a filter from FDATool to SPTool using the Filter Design Toolbox 3.1 (R14SP1) and the Signal Processing Toolbox 6.2.1 (R14SP1)

butterworthfdatoolimportSignal Processing Toolboxsptool

I am creating a lowpass Butterworth filter using FDATool. I export the filter to the workspace, open SPTool, and import the filter. When I view the filter, the frequency axis for the magnitude response is incorrect. This seems to be due to the fact that some of the filter specifications have not been imported.

Best Answer

This bug has been fixed in Release 14 Service Pack 3 (R14SP3). For previous product releases, read below for any possible workarounds:

This is expected behavior in the Signal Processing Toolbox 6.2.1 (R14SP1) in the way that filters are imported into SPTool from FDATool. This is due to the fact that the Fs value is set to 2 by default when a filter is exported from FDATool to SPTool.

Related Solutions

MATLAB: Am I unable to view the coefficients of the arbitrary group delay filter in the FDATOOL GUI using Signal Processing Toolbox 6.2.1 (R14SP1)

This bug has been fixed in Release 14 Service Pack 2 (R14SP2). For previous releases, please read below for any possible workarounds:

We have verified that there is a bug in Signal Processing Toolbox 6.2.1 (R14SP1) in the way that FVTOOL handles viewing the coefficients of Arbitrary Group Delay filters. FDATOOL uses an integrated version of FVTOOL to allow you to analyze the filters you have designed. The "[b,a]" item in the FDATOOL toolbar is one of these integrated tools.

To work around this issue, export the filter coefficients to the workspace to view their values.

MATLAB: Do I observe a shift in mean when using the DECIMATE function in MATLAB 6.5 (R13)

The shift in the mean is caused by the order of the filter used in DECIMATE. By default, an 8th order Chebyshev Type I Low Pass Filter is used as an anti-aliasing filter before removing samples. However, there is a cost function that evaluates the magnitude response of the filter at the cut-off frequency to determine the filter order. Depending on the machine used (Intel Pentium, AMD Athlon, etc) and the operating system (Windows XP, Linux, etc), different numerical libraries produce different results in the cost function and as a result, different filter orders are selected for the same DECIMATE command.

Even order filters generate a 0.05 dB attenuation at f = 0, which then produces the shift in the mean. Odd order filters do not exhibit this behavior at f = 0 and hence do not produce the shift in the mean.

The workaround is to use DECIMATE with the FIR filter option. The FIR filter does not have a 0.05 dB attenuation at f = 0 and hence the shift in mean does not appear in the downsampled signal.

The reason for evaluating the cost function that computes the magnitude response of the filter at the cut-off frequency is that as the passband-edge frequency decreases, the filter that is designed may degenerate for the default 8th order filter. You can see the distortion by following these steps:

1. Start MATLAB

2. Open decimate.m

3. Put a breakpoint on line 126 of decimate.m

4. In the Command Window, enter the following command:

decimate([1:10],61)

5. When MATLAB breaks at line 126, view the 8th order Chebyshev Type I IIR filter by entering the following command:

fvtool(b,a)

6. Zoom into the cut-off frequency of the filter and you will notice that the passband has significant distortion.

The distortion is caused by forming the transfer function from poles and zeros and the roundoff error that accumulates from the convolutions involved. By lowering the filter order, we can reduce this distortion significantly. You can verify this by resuming the MATLAB debugger and viewing the frequency response of the 7th, 6th, etc order filters at the cutoff frequency.

Best Answer

Related Solutions

MATLAB: Am I unable to view the coefficients of the arbitrary group delay filter in the FDATOOL GUI using Signal Processing Toolbox 6.2.1 (R14SP1)

MATLAB: Do I observe a shift in mean when using the DECIMATE function in MATLAB 6.5 (R13)

Related Question