enbw

R2026b

Equivalent noise bandwidth

Syntax

bw = enbw(window)

bw = enbw(window,fs)

Description

bw = enbw(window) returns the two-sided equivalent noise bandwidth, bw, for a uniformly sampled window, window. The equivalent noise bandwidth is normalized by the noise power per frequency bin.

example

bw = enbw(window,fs) returns the two-sided equivalent noise bandwidth, bw, in Hz.

example

Examples

collapse all

Equivalent Noise Bandwidth of Hamming Window

Open Live Script

Determine the equivalent noise bandwidth of a Hamming window 1000 samples in length.

bw = enbw(hamming(1000))

bw = 
1.3638

Equivalent Noise Bandwidth of Flat Top Window

Open Live Script

Determine the equivalent noise bandwidth in Hz of a flat top window 10000 samples in length. The sample rate is 44.1 kHz.

bw = enbw(flattopwin(10000),44.1e3)

bw = 
16.6285

Equivalent Rectangular Noise Bandwidth

Open Live Script

Obtain the equivalent rectangular noise bandwidth of a 128-sample Hann window.

Generate the window and compute its discrete Fourier transform over 2048 frequencies. Shift the transform so it is centered at zero frequency and compute its squared magnitude.

lw = 128;
win = hann(lw);
lt = 2048;
windft = fftshift(fft(win,lt));

ad = abs(windft).^2;
mg = max(ad);

Specify a sample rate of 1 kHz. Use enbw to compute the equivalent noise bandwidth of the window and verify that the value coincides with the definition.

fs = 1000;

bw = enbw(win,fs)

bw = 
11.8110

bdef = sum((win).^2)/sum(win)^2*fs

bdef = 
11.8110

Plot the squared magnitude of the window. Overlay a rectangle with height equal to the peak of the squared magnitude and width equal to the equivalent noise bandwidth.

freq = -fs/2:fs/lt:fs/2-fs/lt;

plot(freq,ad, ...
    bw/2*[-1 -1 1 1],mg*[0 1 1 0],'--')
xlim(bw*[-1 1])

Figure contains an axes object. The axes object contains 2 objects of type line.

Verify that the area of the rectangle contains the same total power as the window.

Adiff = trapz(freq,ad)-bw*mg

Adiff = 
1.4552e-11

Repeat the computation using normalized frequencies. Find the equivalent noise bandwidth of the window. Verify that enbw gives the same value as the definition.

bw = enbw(win)

bw = 
1.5118

bdef = sum((win).^2)/sum(win)^2*lw

bdef = 
1.5118

Plot the squared magnitude of the window. Overlay a rectangle with height equal to the peak of the squared magnitude and width equal to the equivalent noise bandwidth.

freqn = -1/2:1/lt:1/2-1/lt;

plot(freqn,ad, ...
    bw/2*[-1 -1 1 1]/lw,mg*[0 1 1 0],'--')
xlim(bw*[-1 1]/lw)

Figure contains an axes object. The axes object contains 2 objects of type line.

Verify that the area of the rectangle contains the same total power as the window.

Adiff = trapz(freqn,ad)-bw*mg/lw

Adiff = 
0

Input Arguments

collapse all

`window` — Window vector
real-valued row or column vector

Uniformly sampled window vector, specified as a row or column vector with real-valued elements.

Example: hamming(1000)

Data Types: double | single

`fs` — Sample rate
positive scalar

Sample rate, specified as a positive scalar.

Output Arguments

collapse all

`bw` — Equivalent noise bandwidth
positive scalar

Equivalent noise bandwidth, specified as a positive scalar.

Data Types: double | single

More About

collapse all

Equivalent Noise Bandwidth

The total power contained in a window is the area under the curve produced by the squared magnitude of the Fourier transform of the window. The equivalent noise bandwidth of the window is the width of a rectangle with height equal to the peak of the squared magnitude of the Fourier transform and area equal to the total power.

Parseval's theorem states that the total energy for the N-sample window w(n) equals the sum of the squared magnitudes of the time-domain window samples or the integral over frequency of the squared magnitude of the Fourier transform of the window:

$\int_{- 1 / 2}^{1 / 2} {| W (f) |}^{2} d f = \sum_{n} {| w (n) |}^{2} .$

The power spectrum of the window has peak magnitude at f = 0:

${| W (0) |}^{2} = {| \sum_{n} w (n) |}^{2} .$

To find the width of the equivalent rectangle, and thus the equivalent noise bandwidth, divide the area by the height:

$\frac{\int_{- 1 / 2}^{1 / 2} {| W (f) |}^{2} d f}{{| W (0) |}^{2}} = \frac{\sum_{n} {| w (n) |}^{2}}{{| \sum_{n} w (n) |}^{2}} .$

When a sample rate f_s is specified, enbw returns the previous expression multiplied by f_s:

$ENBW = f_{s} \frac{\sum_{n} {| w (n) |}^{2}}{{| \sum_{n} w (n) |}^{2}} .$

For normalized frequencies, the function divides by the bin width, f_s/N:

$ENBW = N \frac{\sum_{n} {| w (n) |}^{2}}{{| \sum_{n} w (n) |}^{2}} .$

See Equivalent Rectangular Noise Bandwidth for an example that, for a given window:

Compares the result of enbw to the definition for a specified sample rate and for normalized frequencies.
Plots the equivalent rectangular bandwidth over the magnitude spectrum of the window.
Verifies that the total power contained in the window is the same as the power contained in a rectangle with height equal to the peak squared magnitude of the Fourier transform of the window and width equal to the equivalent noise bandwidth.

References

[1] harris, fredric j. “On the Use of Windows for Harmonic Analysis with the Discrete Fourier Transform.” Proceedings of the IEEE^®. Vol. 66, January 1978, pp. 51–83.

Extended Capabilities

expand all

C/C++ Code Generation
Generate C and C++ code using MATLAB® Coder™.

GPU Code Generation
Generate CUDA® code for NVIDIA® GPUs using GPU Coder™.

Thread-Based Environment
Run code in the background using MATLAB® `backgroundPool` or accelerate code with Parallel Computing Toolbox™ `ThreadPool`.

The enbw function fully supports thread-based environments. For more information, see Run MATLAB Functions in Thread-Based Environment.

Version History

Introduced in R2013a

expand all

R2026b: Use thread-based environments

You can use enbw in thread-based environments. For more information, see Run MATLAB Functions in Thread-Based Environment.

enbw

Syntax

Description

Examples

Equivalent Noise Bandwidth of Hamming Window

Equivalent Noise Bandwidth of Flat Top Window

Equivalent Rectangular Noise Bandwidth

Input Arguments

window — Window vector real-valued row or column vector

fs — Sample rate positive scalar

Output Arguments

bw — Equivalent noise bandwidth positive scalar

More About

Equivalent Noise Bandwidth

References

Extended Capabilities

C/C++ Code Generation Generate C and C++ code using MATLAB® Coder™.

GPU Code Generation Generate CUDA® code for NVIDIA® GPUs using GPU Coder™.

Thread-Based Environment Run code in the background using MATLAB® backgroundPool or accelerate code with Parallel Computing Toolbox™ ThreadPool.

Version History

R2026b: Use thread-based environments

See Also

`window` — Window vector
real-valued row or column vector

`fs` — Sample rate
positive scalar

`bw` — Equivalent noise bandwidth
positive scalar

C/C++ Code Generation
Generate C and C++ code using MATLAB® Coder™.

GPU Code Generation
Generate CUDA® code for NVIDIA® GPUs using GPU Coder™.

Thread-Based Environment
Run code in the background using MATLAB® `backgroundPool` or accelerate code with Parallel Computing Toolbox™ `ThreadPool`.