Main Content

Filter Analyzer

View, analyze, and compare filters

Since R2024a. Recommended over FVTool.

Description

The Filter Analyzer app is an interactive tool for visualizing, analyzing, and comparing digital filters. Using the app, you can

  • Import filter objects or filter coefficients. For more information, see Import Filter.

  • View, analyze, and compare responses of multiple digital filters. For more information, see Analysis.

  • View a list of filters in the Filters table and the details for each filter in the Filter Information table.

  • Plot a new filter analysis plot in a separate display window.

  • Specify filter sample rate and analysis sample rate separately.

  • Save the state of the current session, including filters, displays, and analysis options, for use in a future session of the app.

You can use the filterAnalyzer command-line interface to add, remove, and update filters, displays, and analysis options.

Filter Analyzer app

Open the Filter Analyzer App

  • MATLAB® Toolstrip: On the Apps tab, under Signal Processing and Communications, click the app icon.

  • MATLAB command prompt: Enter filterAnalyzer.

Examples

expand all

Design sixth-order Chebyshev Type 1 lowpass and highpass filters with a normalized passband edge frequency of 0.4π rad/sample and 12 dB of passband ripple. Express the designs as cascades of fourth-order transfer functions. Use Filter Analyzer to display the magnitude and phase responses of the filters.

[zl,pl,kl] = cheby1(6,12,0.4);
[bl,al] = zp2ctf(zl,pl,kl,SectionOrder=4);
[zh,ph,kh] = cheby1(6,12,0.4,"high");
[bh,ah] = zp2ctf(zh,ph,kh,SectionOrder=4);

filterAnalyzer(bl,al,bh,ah,OverlayAnalysis="phase")

Design a 30th-order FIR bandpass filter to filter signals sampled at 2 kHz. Specify a stopband ranging from 250 Hz to 400 Hz. Use Filter Analyzer to display the poles and zeros of the filter transfer function.

iirbp = designfilt("bandpassfir",FilterOrder=30, ...
    CutoffFrequency1=250,CutoffFrequency2=400, ...
    SampleRate=2000);

filterAnalyzer(iirbp,Analysis="polezero")

Design a tenth-order elliptic bandpass filter with 5 dB of passband ripple and 60 dB of stopband attenuation. Specify passband edge frequencies of 0.2π rad/sample and 0.45π rad/sample. Express the design as a cascade of fourth-order transfer functions.

[z,p,k] = ellip(5,5,60,[0.2 0.45]);
[bb,aa] = zp2ctf(z,p,k,SectionOrder=4);

Design a finite impulse response bandpass filter with 5 dB of passband ripple and asymmetric stopbands for use with signals sampled at 2 kHz.

  • At lower frequencies, the stopband has 80 dB of attenuation and the transition region ranges from 500 Hz to 600 Hz.

  • At higher frequencies, the stopband has 40 dB of attenuation and the transition region ranges from 750 Hz to 900 Hz.

dfir = designfilt("bandpassfir", ...
    SampleRate=2e3,PassbandRipple=5, ...
    StopbandFrequency1=500,PassbandFrequency1=600, ...
    StopbandAttenuation1=80, ...
    PassbandFrequency2=750,StopbandFrequency2=900, ...
    StopbandAttenuation2=40);

Start a Filter Analyzer session to analyze the filters. Import the filters. On the Analyzer tab, click Import Filter.

  • To import the elliptic filter, select Filter Coefficients. Select bb as the Numerator and aa as the Denominator. Specify Ellip for the Filter Name, and leave the Sample Rate as Normalized. Click Import.

  • To import the FIR filter, select Filter Objects, select dfir, and click Import and Close.

Alternatively, open Filter Analyzer by using the command-line interface. By default, the app displays magnitude responses. Only one of the filters has a sample rate, so the app displays the responses using normalized frequencies.

fa = filterAnalyzer(bb,aa,dfir,FilterNames=["ellip" "dfir"]);

Add a display and use it to plot the magnitude responses and the phase responses of the filters.

  • On the Analyzer tab, click New Display.

  • Expand the Analysis gallery so the Overlay Analysis section is visible, and click Phase.

  • Add the filters by clicking the eye icons on the Filters table.

Alternatively, use the filterAnalysisOptions object and the addDisplays and showFilters functions.

opts = filterAnalysisOptions(OverlayAnalysis="phase");
addDisplays(fa,AnalysisOptions=opts)
showFilters(fa,true,FilterNames=["ellip" "dfir"])

Add another display and show the cumulative magnitude responses of the cascaded transfer functions that specify the elliptic filter.

  • Click New Display to add the display and click the eye icon for the elliptic filter.

  • On the Display Options tab, click CTF View ▼ and select Cumulative.

Alternatively, use the command-line interface.

addDisplays(fa,CTFAnalysisMode="cumulative")
showFilters(fa,true,FilterNames="ellip")

Add one more display and show the filter specification mask for the FIR filter.

  • On the Analyzer tab, click New Display and click the eye icon for the FIR filter. The app displays the frequencies in units of Hz.

  • For a digitalFilter object, the app shows the specification mask by default. To remove the mask, on the Display Options tab, click Mask ▼ and clear Specification.

Alternatively, use the command-line interface.

addDisplays(fa)
showFilters(fa,true,FilterNames="dfir")

Related Examples

Parameters

expand all

Use the Filters table to edit information about the filters under analysis.

  • Name — Double-click an entry in this column to edit the name of the filter. You can also right-click the filter in the Filters table and select Rename.

  • Line — Click an entry in this column to edit the color used to display that filter response.

  • Eye — Click an entry in this column to add the corresponding filter to the current display or to remove it. You can also (since R2024b) drag filters from the Filters table onto any display for visualization.

  • Cursor — Click an entry in this column to show or remove data cursors for the corresponding filter. Pause on an entry in this column to select between showing one cursor and showing two.

  • Sample Rate — Click an entry in this column to edit the sample rate for the corresponding filter.

To duplicate a filter, right-click it in the Filters table and select Duplicate.

To delete a filter from the app, right-click it in the Filters table and select Delete.

Click Import Filter to import the filter you want to analyze. Filter Analyzer supports these filter types.

Filter Coefficients

You can use Filter Analyzer to analyze filters specified as numerator and denominator coefficients. If you specify the coefficients as the L-row matrices

B=[b11b12b1,m+1b21b22b2,m+1bL1bL2bL,m+1],A=[a11a12a1,n+1a21a22a2,n+1aL1aL2aL,n+1],

Filter Analyzer assumes you have specified the filter as a sequence of L cascaded transfer functions (CTF), such that the full transfer function of the filter is

H(z)=b11+b12z1++b1,m+1zma11+a12z1++a1,n+1zn×b21+b22z1++b2,m+1zma21+a22z1++a2,n+1zn××bL1+bL2z1++bL,m+1zmaL1+aL2z1++aL,n+1zn,

where m ≥ 0 is the numerator order of the filter and n ≥ 0 is the denominator order.

  • If L = 1, then B and A are row vectors that specify the transfer function of an IIR filter.

  • If you specify both B and A as column vectors, Filter Analyzer assumes they represent the transfer function of an IIR filter.

  • If B is a scalar, Filter Analyzer assumes you specified the filter as a cascade of all-pole IIR filters with each section having a scaling gain equal to B.

  • If A is a scalar, Filter Analyzer assumes you specified the filter as a cascade of FIR filters with each section having a scaling gain equal to 1/A.

Note

  • To convert second-order section matrices to cascaded transfer functions, use the sos2ctf function.

  • To convert a zero-pole-gain filter representation to cascaded transfer functions, use the zp2ctf function.

Coefficients and Gain

If you have a scaling gain separate from the coefficient values, you can enter it in Filter Analyzer using the Import Filters dialog box. At the command line, you can specify the coefficients and gain as a cell array of the form {B,A,g}, where B and A are as defined in the Filter Coefficients section.

The gain can be a scalar overall gain or a vector of section gains.

  • If the gain is a scalar, Filter Analyzer applies the value uniformly to all the cascade filter sections.

  • If the gain is a vector, it must have one more element than the number of filter sections in the cascade. Filter Analyzer applies a scale value to each of the filter sections and applies the last value uniformly to all the filter sections.

If you specify the coefficient matrices and gain vector as

B=[b11b12b1,m+1b21b22b2,m+1bL1bL2bL,m+1],A=[a11a12a1,n+1a21a22a2,n+1aL1aL2aL,n+1],g=[g1g2gLgS],

Filter Analyzer uses the transfer function

H(z)=gS(g1b11+b12z1++b1,m+1zma11+a12z1++a1,n+1zn×g2b21+b22z1++b2,m+1zma21+a22z1++a2,n+1zn××gLbL1+bL2z1++bL,m+1zmaL1+aL2z1++aL,n+1zn).

digitalFilter Objects

You can use Filter Analyzer to analyze digitalFilter objects. Use designfilt to generate or edit digital filters based on frequency-response specifications.

Filter System Objects

If you have DSP System Toolbox™, you can use Filter Analyzer to analyze these filter System objects.

System object™
dsp.AllpassFilter (DSP System Toolbox)
dsp.AllpoleFilter (DSP System Toolbox)
dsp.CICCompensationDecimator (DSP System Toolbox)
dsp.CICCompensationInterpolator (DSP System Toolbox)
dsp.CICDecimator (DSP System Toolbox)
dsp.CICInterpolator (DSP System Toolbox)

dsp.ComplexBandpassDecimator (DSP System Toolbox) (since R2024b)

dsp.CoupledAllpassFilter (DSP System Toolbox)

dsp.DCBlocker (DSP System Toolbox) (since R2024b)

dsp.Delay (DSP System Toolbox)
dsp.Differentiator (DSP System Toolbox)
dsp.DigitalDownConverter (DSP System Toolbox)
dsp.DigitalUpConverter (DSP System Toolbox)
dsp.FIRDecimator (DSP System Toolbox)
dsp.FIRFilter (DSP System Toolbox)
dsp.FIRHalfbandDecimator (DSP System Toolbox)
dsp.FIRHalfbandInterpolator (DSP System Toolbox)
dsp.FIRInterpolator (DSP System Toolbox)
dsp.FIRRateConverter (DSP System Toolbox)
dsp.FarrowRateConverter (DSP System Toolbox)
dsp.FilterCascade (DSP System Toolbox)
dsp.FourthOrderSectionFilter (DSP System Toolbox)Filter Analyzer does not support fixed-point arithmetic for this object
dsp.HighpassFilter (DSP System Toolbox)
dsp.IIRFilter (DSP System Toolbox)
dsp.IIRHalfbandDecimator (DSP System Toolbox)
dsp.IIRHalfbandInterpolator (DSP System Toolbox)
dsp.LowpassFilter (DSP System Toolbox)
dsp.NotchPeakFilter (DSP System Toolbox)Filter Analyzer does not support fixed-point arithmetic for this System object
dsp.ParallelFilter (DSP System Toolbox)
dsp.SOSFilter (DSP System Toolbox)
dsp.VariableBandwidthFIRFilter (DSP System Toolbox)
dsp.VariableBandwidthIIRFilter (DSP System Toolbox)

If you also have DSP HDL Toolbox™, you can use Filter Analyzer to analyze these filter System objects.

System object

dsphdl.FIRFilter (DSP HDL Toolbox) (since R2024b)

dsphdl.BiquadFilter (DSP HDL Toolbox) (since R2024b)

Filter Analyzer supports these analysis types. To access the options available for each analysis, use the Analysis Options button in the Analysis Options toolstrip section.

Frequency-Domain Analyses

AnalysisIconOptionsDescription
Magnitude response

Magnitude Mode: Select dB, Linear, Squared, or Zero Phase

Normalize Magnitude: Toggle on or off

  • The frequency response is the filter transfer function evaluated at the unit circle, H(z) = H(e).

  • The frequency response is generally complex-valued. You can choose to plot the magnitude of the frequency response, either in linear units or in decibels, or its magnitude squared.

  • The zero-phase response Hr(ω) is related to the frequency response by H(e) = Hr(ω)e(ω), where ϕ(ω) is the continuous phase.

  • For more information, see freqz and zerophase.

Phase response

Phase Units: Select Radians or Degrees

Phase Display: Select Phase or Continuous

  • The phase response is the frequency-dependent phase angle ϕ(ω) of the filter transfer function.

  • For more information, see phasez, angle, and unwrap.

Group delay response

Group Delay Units: Select Samples or Time
  • The group delay is the derivative of the phase response with respect to frequency.

  • The group delay is the average delay of the filter as a function of frequency.

  • For more information, see grpdelay.

Phase delay response

Phase Units: Select Radians or Degrees
  • The phase delay is the phase response divided by the frequency.

  • The phase delay is the time delay the filter imposes on each frequency component of the input signal.

  • For more information, see phasedelay.

Magnitude estimate

Magnitude Mode: Select dB, Linear, or Squared

Normalize Magnitude: Toggle on or off

Number of Trials: Specify a numeric scalar

  • Filter Analyzer estimates the magnitude response by filtering sinusoids with uniformly distributed random frequencies and computing the ratio of output to input.

  • For more information, see freqrespest (DSP System Toolbox) and freqrespopts (DSP System Toolbox).

Round-off noise power spectrum

Number of Trials: Specify a numeric scalar
  • Filter Analyzer can compute the power spectral density due to the round-off noise produced by quantization errors within the filter. The result is an average over the specified number of trials.

  • For more information, see noisepsd (DSP System Toolbox) and noisepsdopts (DSP System Toolbox).

Time-Domain Analyses

AnalysisIconOptionsDescription
Impulse response

Specify Length: Select Auto or User-defined
  • The impulse response is the response of the filter to a unit impulse.

  • The transfer function is the Fourier transform of the impulse response.

  • For more information, see impz.

Step response

Specify Length: Select Auto or User-defined
  • The step response is the response of the filter to a step input.

  • The step response is the cumulative sum of the impulse response.

  • For more information, see stepz.

Other Analyses

AnalysisIconOptionsDescription
Pole-zero plot

N/A
  • The pole-zero plot shows the poles and zeros of the filter transfer function on the z-plane.

  • Filter Analyzer uses the symbol o to denote a zero and the symbol x to denote a pole.

  • For more information, see zplane.

Filter coefficients

Coefficients Format: Select Decimal, Hexadecimal, or Binary
  • The coefficients depend on the filter structure (cascaded transfer functions, direct-form, or lattice).

  • For filters specified as cascaded transfer functions, Filter Analyzer displays the numerator coefficients, denominator coefficients, and gain of each section of the cascade separately.

  • You must have a Fixed-Point Designer™ license to use the Binary option.

Filter information

N/A
  • The filter information includes the structure, phase, stability, and implementation cost of the filter.

Frequency Normalization

In Filter Analyzer, you can display filter responses using normalized frequencies or in terms of a sample rate of your choice. To select how to display frequencies, in the Analysis Options section of the Analyzer tab, set Frequency Normalization to one of these values.

  • Auto — If some filters do not have a sample rate, the app analyzes filter responses using normalized frequencies measured in rad/sample. If all filters have a sample rate, the app analyzes filter responses using cyclical frequencies measured in Hz.

  • Normalized — The app analyzes filter responses using normalized frequencies measured in rad/sample.

  • Unnormalized — The app analyzes filter responses using cyclical frequencies measured in Hz.

Analysis Sample Rate

You can choose a reference sample rate, to use to compare the filters plotted in a display to each other, by entering a value for Analysis Sample Rate. You can also choose the highest sample rate among all the filters in the display by selecting Max. To edit the sample rate units, click the pencil icon.

Analysis Options

Select Analysis Options to open the Analysis Options pane, which contains these options for all frequency-domain analyses:

  • Frequency Scale — To display responses in linear frequency scale, select Linear. To use a logarithmic frequency scale, select Log.

  • Frequency Range — To display responses over positive frequencies only, select One Sided. To display responses over the whole frequency range, select Two Sided. To display responses over the whole frequency range centered at zero, select Centered. To display responses over a custom range of frequencies, select User-defined.

  • Number of Points — Select the number of discrete Fourier transform points you want to use to display frequency-domain responses.

For a list of options available with each analysis type, see Analysis.

Filter Analyzer supports these display options:

  • Legend — Click to show or hide legends in the active display or in all displays.

  • Grid — Click to show or hide the grid in the active display or in all displays.

  • Hide Cursors — Click to hide the cursors in the active display or in all displays. To show the cursor corresponding to a given filter, click the cursor column of a signal in the Filters table.

  • Mask — Click to show or hide the spectral mask in the active display. You can use a standard filter specification mask or you can define your own.

    Specify a custom mask as a set of frequencies and a corresponding set of values. To specify your selection, click User Settings. You can use normalized frequencies or express frequencies in units of Hz, kHz, MHz, or GHz. You can specify values for the magnitude of interest, the square of the magnitude, or the magnitude expressed in dB. The frequencies and the values must be finite.

    Note

    To use standard specification masks, filters must have design metadata. To include metadata in your design, use a digitalFilter object or a filter System object.

If you specify a filter as a cascade of transfer functions, you can choose the way that Filter Analyzer displays section-by-section responses. This option applies only for displays with one filter. In the Display Options tab, click CTF View to toggle between displaying the overall response and displaying the responses section by section. Click CTF View ▼ to select one of these options:

  • Individual — Show filter responses of individual sections.

  • Cumulative — Show filter responses of cumulative sections.

  • User-defined — Show filter responses of selected sections or combinations of sections. To specify your selection, click User Settings and specify Sections a cell array.

    For example, {[1 2 3] [4 5 6]} directs the app to display responses for a cascade of sections 1, 2, and 3 and a cascade of sections 4, 5, and 6.

Programmatic Use

expand all

filterAnalyzer(filt1,...,filtn) plots the responses of the specified filters in the Filter Analyzer app.

  • If Filter Analyzer is not open, this syntax opens the app and plots the responses.

  • If Filter Analyzer is open, this syntax plots the responses in a new display in the app.

Specify input filters as coefficient matrices, cell arrays, digitalFilter objects, or System objects. For more information, see Import Filter.

filterAnalyzer(filt1,...,filtn,Name=Value) specifies nondefault options using one or more name-value arguments.

filterAnalyzer(filename) opens a Filter Analyzer session stored in the specified MAT file filename. If Filter Analyzer is already open, this syntax replaces the current app session with the new session.

filterAnalyzer(filename,"append") appends the filters stored in the specified MAT file filename to the current Filter Analyzer session. If Filter Analyzer is not open, this syntax is equivalent to the previous syntax.

[fa,dispnums] = filterAnalyzer(___) returns a handle object for the Filter Analyzer and the numbers corresponding to the newly added displays, using any combination of input arguments from previous syntaxes.. You can also obtain the handle fa by typing fa = getFilterAnalyzerHandle at the command line.

Version History

Introduced in R2024a

expand all