Skip to main content

How Windowing Affects Your Periodogram


The windowed periodogram is a widely used technique for estimating the Power Spectral Density (PSD) of a signal. It enhances the classical periodogram by mitigating spectral leakage through the application of a windowing function. This technique is essential in signal processing for accurate frequency-domain analysis.

 

Power Spectral Density (PSD)

The PSD characterizes how the power of a signal is distributed across different frequency components. For a discrete-time signal, the PSD is defined as the Fourier Transform of the signal’s autocorrelation function:

Sx(f) = FT{Rx(Ï„)}

Here, Rx(Ï„)}is the autocorrelation function.

FT : Fourier Transform

 

Classical Periodogram

The periodogram is a non-parametric PSD estimation method based on the Discrete Fourier Transform (DFT):

Px(f) = \(\frac{1}{N}\) X(f)2

Here:

  • X(f): DFT of the signal x(n)

  • N: Signal length

However, the classical periodogram suffers from spectral leakage due to abrupt truncation of the signal.

 

Windowing to Mitigate Spectral Leakage

Spectral leakage can be minimized by applying a window function to the signal before computing the DFT. The resulting PSD estimate is called the windowed periodogram:

Pw(f) = \(\frac{1}{NW}\) Xw(f)2

Here:

  • w(n): Window function

  • W: Window normalization factor

Common Window Functions

  • Rectangular Window: Equivalent to the classical periodogram.

w[n]=1, 0≤n≤N−1

w[n]=0, otherwise

Where, N is the window length

  • Hamming Window: Reduces sidelobe amplitudes, improving frequency resolution.

w[n]=0.5(1−cos(\(\frac{\ 2\pi n}{N - 1}\ \))), 0≤n≤N−1

Where, N is the window length

  • Hanning Window: Similar to Hamming but with less sidelobe attenuation.

w[n]=0.54 – 0.46cos(\(\frac{\ 2\pi n}{N - 1}\ \)), 0≤n≤N−1

Where, N is the window length

  • Blackman Window: Offers even greater sidelobe suppression but at the cost of wider main lobes.

w[n]=0.42 – 0.5(cos(\(\frac{\ 2\pi n}{N - 1}\ \)) + 0.08(cos(\(\frac{\ 4\pi n}{N - 1}\ \)), 0≤n≤N−1

Where, N is the window length

 

Implementation Steps

  1. Segment the Signal: Divide the signal into overlapping or non-overlapping segments of length N.

  2. Apply a Window Function: Multiply each segment by a window function w(n).

  3. Compute the DFT: Calculate the DFT of the windowed segments.

  4. Average the Periodograms: For overlapping segments, average the periodograms to reduce variance.

     

Properties of the Windowed Periodogram

  • Bias: Windowing introduces bias in the PSD estimate as the window modifies the signal spectrum.

  • Variance: Averaging periodograms (Welch method) reduces variance but decreases frequency resolution.

  • Trade-Off: The choice of window affects the trade-off between spectral resolution and leakage suppression.

     

    MATLAB Code

    clc;
    clear;
    close all;

    fs = 48000;
    t = 0:1/fs:0.02;
    f_ping = 12000;

    % Base sine wave
    sine_wave = sin(2*pi*f_ping*t)';

    % Apply windows
    w_rect = ones(size(sine_wave));
    w_hann = hann(length(sine_wave));
    w_hamming = hamming(length(sine_wave));
    w_blackman = blackman(length(sine_wave));

    % Windowed signals
    s_rect = sine_wave .* w_rect;
    s_hann = sine_wave .* w_hann;
    s_hamming = sine_wave .* w_hamming;
    s_blackman = sine_wave .* w_blackman;

    % FFT
    Nfft = 4096;
    f = fs*(0:Nfft/2-1)/Nfft;

    % Function to compute and normalize spectrum
    get_norm_fft = @(sig) abs(fft(sig, Nfft))/max(abs(fft(sig, Nfft)));

    S_rect = get_norm_fft(s_rect);
    S_hann = get_norm_fft(s_hann);
    S_hamming = get_norm_fft(s_hamming);
    S_blackman = get_norm_fft(s_blackman);

    % Mainlobe power (±2 bins around peak)
    mainlobe_bins = 2;

    % Function to compute power ratio
    compute_power_ratio = @(S) ...
    deal( ...
    sum(S.^2), ... % Total power
    max(1, find(S == max(S), 1)), ... % Peak bin
    @(peak_bin) sum(S(max(1,peak_bin-mainlobe_bins):min(Nfft,peak_bin+mainlobe_bins)).^2), ...
    @(total, main) 10*log10((total-main)/main) ... % dB sidelobe/mainlobe ratio
    );

    % Calculate ratios
    [total_r, peak_r, get_main_r, get_slr_r] = compute_power_ratio(S_rect);
    main_r = get_main_r(peak_r); slr_r = get_slr_r(total_r, main_r);

    [total_h, peak_h, get_main_h, get_slr_h] = compute_power_ratio(S_hann);
    main_h = get_main_h(peak_h); slr_h = get_slr_h(total_h, main_h);

    [total_ham, peak_ham, get_main_ham, get_slr_ham] = compute_power_ratio(S_hamming);
    main_ham = get_main_ham(peak_ham); slr_ham = get_slr_ham(total_ham, main_ham);

    [total_b, peak_b, get_main_b, get_slr_b] = compute_power_ratio(S_blackman);
    main_b = get_main_b(peak_b); slr_b = get_slr_b(total_b, main_b);

    % Display Results
    fprintf('Window | Mainlobe Power | Sidelobe Power | Sidelobe/Main (dB)\n');
    fprintf('------------|----------------|----------------|--------------------\n');
    fprintf('Rectangular | %14.4f | %14.4f | %18.2f\n', main_r, total_r - main_r, slr_r);
    fprintf('Hann | %14.4f | %14.4f | %18.2f\n', main_h, total_h - main_h, slr_h);
    fprintf('Hamming | %14.4f | %14.4f | %18.2f\n', main_ham, total_ham - main_ham, slr_ham);
    fprintf('Blackman | %14.4f | %14.4f | %18.2f\n', main_b, total_b - main_b, slr_b);

    % Plot
    figure;
    plot(f, 20*log10(S_rect(1:Nfft/2)), 'k'); hold on;
    plot(f, 20*log10(S_hann(1:Nfft/2)), 'r');
    plot(f, 20*log10(S_hamming(1:Nfft/2)), 'g');
    plot(f, 20*log10(S_blackman(1:Nfft/2)), 'b');
    legend('Rectangular','Hann','Hamming','Blackman');
    xlim([f_ping-3000 f_ping+3000]); ylim([-100 5]);
    xlabel('Frequency (Hz)'); ylabel('Magnitude (dB)');
    title('Windowing Effects on Spectrum');
    grid on;

    Output 

    Window      | Mainlobe Power | Sidelobe Power | Sidelobe/Main (dB)
    ------------|----------------|----------------|--------------------
    Rectangular |         3.5771 |         4.9562 |               1.42
    Hann        |         4.3630 |         8.4370 |               2.86
    Hamming     |         4.2367 |         7.3928 |               2.42
    Blackman    |         4.4940 |        10.2410 |               3.58

     

     








Applications

  • Signal Processing: Analyzing frequency content of time-varying signals.

  • Communications: Evaluating spectrum occupancy in wireless systems.

  • Bioinformatics: Investigating periodicities in biological signals (e.g., EEG, ECG).

  • Seismology: Characterizing seismic wave frequencies.

     

    Further Reading

    1. Periodogram in MATLAB

Try Interactive Online Simulator to see the Effect

📊 MATLAB & ONLINE SIMULATIONS (Web-Based)
🎓 GATE & UGC - NET SOLUTIONS
📡 WIRELESS & SIGNAL PROCESSING


Contact Us

Name

Email *

Message *

Popular Posts

Q-function in BER vs SNR Calculation (with Simulation)

Q-function in BER vs. SNR Calculation In digital communications and signal processing, the Q-function plays a significant role in predicting system reliability. It allows engineers to quantify the probability that Gaussian noise will exceed a specific threshold, causing a bit error. What is the Q-function? The Q-function is a mathematical function representing the tail probability of the standard normal (Gaussian) distribution. It is the complementary cumulative distribution function (CCDF) of a standard Gaussian distribution. Q(x) = (1 / √(2Ï€)) ∫â‚“∞ e^(-t² / 2) dt Q-Function Interactive Simulator Move the slider to see how the "Tail Probability" (the area in red) changes. This area represents the Probability of Error (BER) . Threshold Distance ( x ) — (Simulates Increasing SNR) x = 1.0 Q(x) = 0.1587 ...
Read more

BER vs SNR for M-ary QAM, M-ary PSK, QPSK, BPSK, ...(MATLAB Code + Simulator)

Bit Error Rate (BER) & SNR Guide Analyze communication system performance with our interactive simulators and MATLAB tools. 📘 Theory 🧮 Simulators 💻 MATLAB Code 📚 Resources BER Definition SNR Formula BER Calculator MATLAB Comparison 📂 Explore M-ary QAM, PSK, and QPSK Topics ▼ 🧮 Constellation Simulator: M-ary QAM 🧮 Constellation Simulator: M-ary PSK 🧮 BER calculation for ASK, FSK, and PSK 🧮 Approaches to BER vs SNR What is Bit Error Rate (BER)? The BER indicates how many corrupted bits are received compared to the total number of bits sent. It is the primary figure of merit f...
Read more

Frequency Shift Keying (FSK) Modulation & Demodulation (with Simulation)

Frequency Shift Keying (FSK) Theoretical Foundations: Frequency Shift Keying (FSK) is a discrete frequency modulation scheme wherein the digital information is encoded via instantaneous shifts in the carrier signal's frequency. The fundamental implementation is Binary FSK (BFSK), which maps binary data onto two distinct, discrete spectral states. A binary '1' (the "mark" state) is represented by a carrier frequency \( f_1 \), while a binary '0' (the "space" state) corresponds to frequency \( f_2 \). Each symbol is sustained for a bit interval denoted by \( T_b \). FSK Transmitter Characterization: The mathematical model for the modulated BFSK output \( s(t) \) is defined as: \[ s(t) = \begin{cases} A_c \cos(2\pi f_1 t), & \text{for } m = 1 \\ A_c \cos(2\pi f_2 t), & \text{for } m = 0 \end{cases} \] ...
Read more

RMS Delay Spread, Excess Delay Spread and Multi-path ...(with MATLAB + Simulator)

📘 Overview of Delay Spread and Multi-path 🧮 Excess Delay spread 🧮 Power delay Profile 🧮 RMS Delay Spread 📚 Further Reading 📂 Other Topics on RMS Delay Spread, Excess Delay ... 🧮 Multipath Components or MPCs 🧮 Online Simulator for Calculating RMS Delay Spread 🧮 Why is there significant multipath in the case of very high frequencies? 🧮 Why RMS Delay Spread is essential for wireless communication? 🧮 Why the Power Delay Profile is essential? 🧮 MATLAB Codes for Calculating Different Types of delay Spreads Delay Spread, Excess Delay Spread, and Multipath (MPCs) The fundamental distinction between wireless and wired connections is that in wireless connections signal reaches at receiver thru multipath signal propagation rather than directed transmission like co-axial cable. Wireless Communication has no set communication path between the transmitter and the receiver. The line...
Read more

OFDM Symbols and Subcarriers Explained

This article explains how OFDM (Orthogonal Frequency Division Multiplexing) symbols and subcarriers work. It covers modulation, mapping symbols to subcarriers, subcarrier frequency spacing, IFFT synthesis, cyclic prefix, and transmission. Step 1: Modulation First, modulate the input bitstream. For example, with 16-QAM , each group of 4 bits maps to one QAM symbol. Suppose we generate a sequence of QAM symbols: s0, s1, s2, s3, s4, s5, …, s63 Step 2: Mapping Symbols to Subcarriers Assume N sub = 8 subcarriers. Each OFDM symbol in the frequency domain contains 8 QAM symbols (one per subcarrier): Mapping (example) OFDM symbol 1 → s0, s1, s2, s3, s4, s5, s6, s7 OFDM symbol 2 → s8, s9, s10, s11, s12, s13, s14, s15 … OFDM sym...
Read more

FFT Butterfly Method Explained (with Example of 4-point DFT)

  FFT Using Butterfly Method Given: x[n] = {0, 1, 2, 3} Step 1: Split into Even & Odd Even indices: x e = {0, 2} Odd indices: x o = {1, 3} Step 2: 2-point DFT For any {a, b}: DFT = {a + b, a - b} Even Part: E = {0+2, 0-2} = {2, -2} Odd Part: O = {1+3, 1-3} = {4, -2} Step 3: Combine Using Butterfly X[k] = E[k] + W k O[k] X[k + N/2] = E[k] - W k O[k] For N = 4: W 0 = 1 W 1 = -j Final Calculations X[0] = 2 + 4 = 6 X[2] = 2 - 4 = -2 X[1] = -2 + (-j)(-2) = -2 + 2j X[3] = -2 - (-j)(-2) = -2 - 2j Final Answer: X[k] = {6, -2 + 2j, -2, -2 - 2j} Try Interactive Online Simulations Interactive FFT Online Simulator (For understanding Fundamentals)  Interactive FFT Online Simulator (Analyze .CSV, .MP3, .MP4, etc. Further Reading Fourier Transform OFDM Return to Fourier Transform Main Page →
Read more

Orthogonal Time Frequency Space (OTFS) (with MATLAB)

In OTFS (Orthogonal Time Frequency Space) modulation — a scheme designed for high-Doppler and time-varying wireless channels — the terms ISFFT and SFFT are key mathematical transformations used to move between different representation domains. Figure: OTFS block diagram 1. ISFFT — Inverse Symplectic Finite Fourier Transform Purpose: Transforms data symbols from the delay-Doppler domain to the time-frequency domain . \[ X[n, m] = \frac{1}{\sqrt{NM}} \sum_{k=0}^{N-1} \sum_{l=0}^{M-1} x[k, l] \, e^{j2\pi \left( \frac{nk}{N} - \frac{ml}{M} \right)} \] Here, \( N \) is the number of Doppler bins (time slots), and \( M \) is the number of delay bins (subcarriers). The ISFFT maps each data symbol from the delay-Doppler grid (where the channel is sparse and easier to equalize) to the time-frequency grid (where standard multicarrier modulation like OFDM can be applied). 2. SFFT — Symplectic Finite Fourier Transform Purpose: Performs the reverse operation ...
Read more