Bandstop FIR Filter Design

Design bandstop FIR filter with tunable filter specifications

Since R2023b

• 

Libraries:
DSP System Toolbox / Filtering / Filter Sources

Description

Use the Bandstop FIR Filter Design block to design a bandstop FIR filter using the Hann, Hamming, Blackman, Blackman-Harris, Chebyshev, and Kaiser windowing methods. You can tune the filter order, center frequency, bandwidth, and the sidelobe attenuation during simulation. When the values of these parameters change, the block redesigns the filter and outputs the numerator coefficients in vector form.

Examples

Design and Implement Bandstop Filter with Tunable Frequency Specifications in Simulink

Vary the order and the bandwidth of the bandstop FIR filter during simulation. Design the filter every time the frequency specifications update. Visualize the magnitude response of this varying filter.

Pass a noisy sinusoidal signal through the bandstop FIR filter. Visualize the frequency spectra of the input and the output signals.

Open and Inspect Model

Open the tunable_bandstop_filter model by clicking the Open Model button.

The input signal in the model is a sum of two sine waves with the frequencies of 3 kHz and 15 kHz. The Random Source block adds zero-mean white Gaussian noise with a variance of 0.05 to the sum of sine waves.

The Bandstop FIR Filter Design block designs the filter and the Discrete FIR Filter block implements the filter. The filter order and the filter bandwidth specifications are input through the input ports of the Bandstop FIR Filter Design block.

Run Model

Visualize the spectra of the original signal and the filtered signal in the Spectrum Analyzer. Both the input tones are unaffected since the frequencies fall in the passband region of the filter.

Tune the Filter Order and Filter Bandwidth

Change the filter order to 90 and the filter bandwidth to 0.5 during simulation. To change the parameter values more gradually during simulation, select the Smooth tuned filter parameters check box and specify a smoothing factor in the Bandstop FIR Filter Design block dialog box. The smoothing factor determines the speed at which the parameter values change until they match the desired new value. If you specify a smoothing factor of 0, the block does not smooth the parameter and immediately sets the parameter to the new value. As the smoothing factor approaches 1, the number of smoothing operations, and consequently, the number of filter redesigns increase. This example uses a smoothing factor of 0.6.

When you change the filter order and the filter bandwidth, you can see the magnitude response of the filter change in the Filter Visualizer output. The second tone at 15 kHz is now attenuated as it falls outside the passband region of the filter.

Ports

Input

N — Filter order

even positive integer

Specify the filter order as an even positive integer less than or equal to N_max, where N_max is the value you specify in the Filter maximum order (must be even) parameter. You can change the filter order you input through this port during simulation.

Dependency

To enable this port, select the Specify filter order from input port parameter.

Data Types: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

Fc — Filter center frequency

positive scalar ≤ 1

Specify the center frequency of the filter in the normalized frequency units as a positive scalar less than or equal to 1. You can change the center frequency you input through this port during simulation.

Dependency

To enable this port, select the Specify center frequency from input port parameter.

Data Types: single | double

BW — Filter bandwidth

positive scalar ≤ 1

Specify the bandwidth of the filter in the normalized frequency units as a positive scalar less than or equal to 1. You can change the bandwidth you input through this port during simulation.

Dependency

To enable this port, select the Specify bandwidth from input port parameter.

Data Types: single | double

r — Sidelobe attenuation in dB

positive scalar

Specify the sidelobe attenuation of the Chebyshev type filter as a positive scalar in dB. You can change the sidelobe attenuation you input through this port during simulation.

Dependency

To enable this port:

• Set the Window function parameter to Chebyshev.

• Select the Specify sidelobe attenuation from input port parameter.

Data Types: single | double

beta — Kaiser window parameter

nonnegative scalar

Specify the Kaiser window parameter as a nonnegative scalar. You can change the Kaiser window parameter you input through this port during simulation.

Dependency

To enable this port:

• Set the Window function parameter to Kaiser.

• Select the Specify kaiser window parameter from input port parameter.

Data Types: single | double

Output

Num — Numerator coefficients

row vector

Numerator coefficients, returned as a row vector of length N_max + 1, where N_max is the maximum filter order that you specify using the Filter maximum order (must be even) parameter.

This equation represents the vector of numerator coefficients b.

$$b=\left[b_0,b_1,b_2,\mathrm{...},b_{N\max }\right]$$

This equation represents the FIR filter in the transfer function form.

$$H(z)=b_0+b_1{z}^{-1}+b_2{z}^{-2}+\mathrm{...}+b_{N\max }{z}^{-N\max }$$

When the actual filter order N is less than the maximum filter order N_max , the last N_max − N elements in the vector are trivial with a value of 0.

The data type of this port depends on the value of the Output data type parameter.

Data Types: single | double

Parameters

Smooth tuned filter parameters — Smooth tuned filter parameters

off (default) | on

Since R2024a

When you select this check box, the block updates the tunable filter parameters such as the filter order, center frequency, bandwidth, sidelobe attenuation, and Kaiser window parameter gradually by repeatedly designing the filter with intermediate values of the parameter until the parameter reaches the desired value. This process is called parameter smoothing. For more information, see Parameter Smoothing.

When you clear this check box, the block redesigns the filter with the new parameter value immediately.

Smoothing factor — Smoothing factor

0.6 (default) | nonnegative scalar less than 1

Since R2024a

Specify the smoothing factor α as a nonnegative scalar less than 1. When α = 0, no smoothing occurs. As α approaches 1, the number of smoothing operations, and consequently, the number of filter redesigns increase.

Tunable: Yes

Dependencies

To enable this parameter, select the Smooth tuned filter parameters parameter.

Specify filter order from input port — Flag to specify filter order from input port

off (default) | on

Select this parameter to specify the filter order from the input port N. When you clear this parameter, you can specify the filter order in the block dialog box using the Filter order (must be even) parameter.

Filter order (must be even) — Filter order

100 (default) | even positive integer

Specify the filter order as an even positive integer that is less than or equal to the value of the Filter maximum order (must be even) parameter.

Tunable: Yes

Dependency

To enable this parameter, clear the Specify filter order from input port parameter.

Data Types: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

Filter maximum order (must be even) — Maximum order of filter

100 (default) | even positive integer

Specify the maximum order of the filter as an even positive integer. The value you specify in the Filter order (must be even) parameter must be less than or equal to the value you specify in the Filter maximum order (must be even) parameter.

Data Types: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

Specify center frequency from input port — Flag to specify center frequency from input port

off (default) | on

Select this parameter to specify the center frequency from the input port Fc. When you clear this parameter, you can specify the center frequency in the block dialog box using the Center frequency parameter.

Center frequency — Center frequency

0.5 (default) | positive scalar ≤ 1

Specify the center frequency of the filter in the normalized frequency units as a positive scalar less than or equal to 1.

Tunable: Yes

Dependency

To enable this parameter, clear the Specify center frequency from input port parameter.

Data Types: single | double

Specify bandwidth from input port — Flag to specify bandwidth from input port

off (default) | on

Select this parameter to specify the bandwidth from the input port BW. When you clear this parameter, you can specify the bandwidth in the block dialog box using the Bandwidth parameter.

Bandwidth — Filter bandwidth

0.1 (default) | positive scalar ≤ 1

Specify the bandwidth of the filter in the normalized frequency units as a positive scalar less than or equal to 1.

Tunable: Yes

Dependency

To enable this parameter, clear the Specify bandwidth from input port parameter.

Data Types: single | double

Window function — Window function

Hann (default) | Hamming | Blackman | Blackman-Harris | Chebyshev | Kaiser

Specify the window function to use to design the filter:

• Hann

• Hamming

• Blackman

• Blackman-Harris

• Chebyshev

• Kaiser

Specify sidelobe attenuation from input port — Flag to specify sidelobe attenuation from input port

off (default) | on

Select this parameter to specify the sidelobe attenuation from the input port r. When you clear this parameter, you can specify the sidelobe attenuation in the block dialog box using the Sidelobe attenuation (dB) parameter.

Dependency

To enable this parameter, set the Window function parameter to Chebyshev.

Sidelobe attenuation (dB) — Sidelobe attenuation in dB

60 (default) | positive scalar

Specify the sidelobe attenuation of the Chebyshev type FIR filter as a positive scalar in dB.

Tunable: Yes

Dependency

To enable this parameter:

• Set the Window function parameter to Chebyshev.

• Clear the Specify sidelobe attenuation from input port parameter.

Data Types: single | double

Specify kaiser window parameter from input port — Flag to specify Kaiser window parameter from input port

off (default) | on

Select this parameter to specify the Kaiser window parameter from the input port beta. When you clear this parameter, you can specify the Kaiser window parameter in the block dialog box using Kaiser window parameter.

Dependency

To enable this parameter, set the Window function parameter to Kaiser.

Kaiser window parameter — Kaiser window parameter

0.5 (default) | nonnegative scalar

Specify the Kaiser window parameter as a nonnegative scalar.

Tunable: Yes

Dependency

To enable this parameter:

• Set the Window function parameter to Kaiser.

• Clear the Specify kaiser window parameter from input port parameter.

Data Types: single | double

Sample time (−1 for inherited) — Sample time in seconds

-1 (default) | positive scalar | Inf

Specify the sample time as -1 (inherited), Inf (constant sample time), or a positive scalar. This parameter determines when the block produces the output and updates its internal state during simulation. For more details, see What Is Sample Time? (Simulink).

Data Types: single | double

Output data type — Data type of filter coefficients

double (default) | single

Specify the data type of the filter coefficients that the block outputs through the Num port. You can set the data type to double or single.

Simulate using — Type of simulation to run

Interpreted execution (default) | Code generation

Specify the type of simulation to run. You can set this parameter to:

• Interpreted execution –– Simulate model using the MATLAB^® interpreter. This option shortens startup time.

• Code generation –– Simulate model using generated C code. The first time you run a simulation, Simulink^® generates C code for the block. The C code is reused for subsequent simulations as long as the model does not change. This option requires additional startup time but provides faster subsequent simulations.

Block Characteristics

Data Types	double | integer | single
Direct Feedthrough	no
Multidimensional Signals	no
Variable-Size Signals	no
Zero-Crossing Detection	no

More About

Parameter Smoothing

Parameter smoothing gradually changes the values of filter parameters until the desired value is reached.

If you select the Smooth tuned filter parameters check box and update the value of tunable filter parameters during simulation, the block repeatedly redesigns the filter with intermediate values of the parameter until the parameter reaches the desired value. The block uses a first-order IIR filter to smooth the parameter.

$$g[n]=\alpha g[n-1]+\left(1-\alpha \right)t[n],$$

where:

• g[n] is the instantaneous smoothed parameter

• t[n] is the discontinuous target parameter

• α is the smoothing factor. When α = 0, no smoothing occurs. As α approaches 1, the number of smoothing operations, and consequently, the number of filter redesigns increase.

Extended Capabilities

C/C++ Code Generation

Generate C and C++ code using Simulink® Coder™.

When you generate code from this block, you can tune the filter design specifications in the generated code.

Version History

Introduced in R2023b

R2024a: Change in the default value of Simulate using parameter

The default value of the Simulate using parameter is now Interpreted execution. With this change, the block uses the MATLAB interpreter for simulation by default.

R2024a: Tunable Bandstop FIR Filter block renamed

The Tunable Bandstop FIR Filter has been renamed as Bandstop FIR Filter Design block.

R2024a: Support for smoothing tunable filter parameters

When tunable filter parameters such as the filter order, center frequency, bandwidth, sidelobe attenuation, and Kaiser window parameter change during simulation, the block can smooth these parameters until the parameters reach the desired value. For more information, see Parameter Smoothing.