Skip to main content

Aliasing in Signal Processing: Nyquist-Shannon Theorem & Best Practices for Downsampling and Upsampling


Aliasing in Signal Processing

1. Nyquist-Shannon Sampling Theorem (Standard)

This is the foundational principle for preventing aliasing in both downsampling and upsampling:

  • The signal must be sampled at least at twice the highest frequency component in the signal (i.e., the sampling rate should be greater than or equal to twice the maximum frequency, known as the Nyquist rate).

Mathematically:

fs ≥ 2 × fmax

Where:

  • fs is the sampling rate.
  • fmax is the maximum frequency present in the signal.

For a signal with frequency content up to fmax, the Nyquist rate is 2 × fmax. If you violate this and sample below this threshold, aliasing will occur.

 

Aliasing in Signal Processing

Aliasing is a phenomenon that occurs when a continuous signal (or sampled signal) is undersampled during digital signal processing. This results in different signals becoming indistinguishable from each other when sampled at a lower rate. Aliasing is a major issue in both downsampling (decimation) and upsampling (interpolation), although it manifests differently in each case.

 

1. Aliasing during Downsampling (Decimation)

Downsampling involves reducing the sampling rate of a signal. When you downsample a signal, you take fewer samples per unit time. If the signal contains frequencies higher than half the new sampling rate (the Nyquist frequency), those high-frequency components get folded back into the lower frequency range, causing aliasing.

Conditions for Aliasing in Downsampling:

  • The original signal contains frequency components above half the new sampling rate (i.e., above the Nyquist frequency).
  • The Nyquist-Shannon sampling theorem states that to avoid aliasing, the original signal should be band-limited to half of the sampling rate (i.e., the Nyquist frequency).
  • If you downsample without low-pass filtering (anti-aliasing filter) the signal to remove these high-frequency components, aliasing will occur.

Example:

Suppose your original signal has a sampling rate of 1000 Hz. This means that the Nyquist frequency is 500 Hz.

If you want to downsample to a new sampling rate of 300 Hz, the Nyquist frequency for this rate is 150 Hz.

Any frequency component in the signal above 150 Hz will alias and appear as lower-frequency components.

To prevent aliasing during downsampling:

  • Pre-Filtering (Anti-Aliasing Filter): Apply a low-pass filter to the signal before downsampling. This filter should remove any frequency components above the Nyquist frequency of the new sampling rate.

For example, if you downsample from a sampling rate of fs = 1000 Hz to fs' = 300 Hz, the filter should remove all frequencies above 150 Hz (half of the new sampling rate).

  • Downsample the Signal: After the low-pass filtering step, downsample the signal by selecting every M-th sample, where M is the downsampling factor.

Mathematical Standard for Downsampling:

fnew = fold / M

Where:

  • fnew is the new, lower sampling rate.
  • fold is the original, higher sampling rate.
  • M is the downsampling factor.

 

2. Aliasing during Upsampling (Interpolation)

Upsampling, or interpolation, involves increasing the sampling rate of a signal. During upsampling, additional samples are inserted between the existing samples to increase the sample rate, but the frequency content of the signal does not change. Upsampling itself doesn't cause aliasing directly, but it can introduce artifacts if not handled properly.

Aliasing in Upsampling occurs due to:

  • The introduction of high-frequency components (above the Nyquist frequency of the new sampling rate) into the signal after the interpolation process. This can happen because interpolation can generate frequencies higher than the Nyquist frequency of the target sampling rate.

Example:

Suppose your original signal has a sampling rate of 200 Hz, and you upsample to 800 Hz.

The Nyquist frequency of the original signal is 100 Hz (half of 200 Hz), so the signal may have frequencies between 0 Hz and 100 Hz.

After upsampling to 800 Hz, the Nyquist frequency is 400 Hz.

If the signal is upsampled by inserting samples, the interpolated signal may introduce high-frequency components (between 100 Hz and 400 Hz) that could cause aliasing.

How to Prevent Aliasing during Upsampling:

  • Zero Padding: Insert M-1 zeros between each original sample. For example, if you want to upsample by a factor of M = 3, insert two zeros between each original sample.
  • Apply Interpolation Filter: After zero padding, apply a low-pass filter (such as an FIR filter) to remove any high-frequency components that might have been introduced during the interpolation process. This filter should have a cutoff at the Nyquist frequency of the new sampling rate.

Mathematical Standard for Upsampling:

fnew = M × fold

Where:

  • fnew is the new, higher sampling rate.
  • fold is the original, lower sampling rate.
  • M is the upsampling factor.

 

Summary of Aliasing:

  • Downsampling (decimation): Aliasing occurs when the signal contains frequency components above the Nyquist frequency of the new sampling rate, which will fold back into the lower frequencies. Low-pass filtering (anti-aliasing filter) before downsampling is necessary to avoid aliasing.
  • Upsampling (interpolation): Aliasing occurs if the interpolation introduces high-frequency components that exceed the Nyquist frequency of the new sampling rate. This can be avoided by using a low-pass filter to remove any high-frequency components before or after interpolation.

 

Further Reading

  1.  

People are good at skipping over material they already know!

View Related Topics to







Contact Us

Name

Email *

Message *

Popular Posts

BER vs SNR for M-ary QAM, M-ary PSK, QPSK, BPSK, ...

📘 Overview of BER and SNR 🧮 Online Simulator for BER calculation of m-ary QAM and m-ary PSK 🧮 MATLAB Code for BER calculation of M-ary QAM, M-ary PSK, QPSK, BPSK, ... 📚 Further Reading 📂 View Other Topics on M-ary QAM, M-ary PSK, QPSK ... 🧮 Online Simulator for Constellation Diagram of m-ary QAM 🧮 Online Simulator for Constellation Diagram of m-ary PSK 🧮 MATLAB Code for BER calculation of ASK, FSK, and PSK 🧮 MATLAB Code for BER calculation of Alamouti Scheme 🧮 Different approaches to calculate BER vs SNR What is Bit Error Rate (BER)? The abbreviation BER stands for Bit Error Rate, which indicates how many corrupted bits are received (after the demodulation process) compared to the total number of bits sent in a communication process. BER = (number of bits received in error) / (total number of tran...

MATLAB Code for ASK, FSK, and PSK

📘 Overview & Theory 🧮 MATLAB Code for ASK 🧮 MATLAB Code for FSK 🧮 MATLAB Code for PSK 🧮 Simulator for binary ASK, FSK, and PSK Modulations 📚 Further Reading ASK, FSK & PSK HomePage MATLAB Code MATLAB Code for ASK Modulation and Demodulation % The code is written by SalimWireless.Com % Clear previous data and plots clc; clear all; close all; % Parameters Tb = 1; % Bit duration (s) fc = 10; % Carrier frequency (Hz) N_bits = 10; % Number of bits Fs = 100 * fc; % Sampling frequency (ensure at least 2*fc, more for better representation) Ts = 1/Fs; % Sampling interval samples_per_bit = Fs * Tb; % Number of samples per bit duration % Generate random binary data rng(10); % Set random seed for reproducibility binary_data = randi([0, 1], 1, N_bits); % Generate random binary data (0 or 1) % Initialize arrays for continuous signals t_overall = 0:Ts:(N_bits...

Antenna Gain-Combining Methods - EGC, MRC, SC, and RMSGC

📘 Overview 🧮 Equal gain combining (EGC) 🧮 Maximum ratio combining (MRC) 🧮 Selective combining (SC) 🧮 Root mean square gain combining (RMSGC) 🧮 Zero-Forcing (ZF) Combining 🧮 MATLAB Code 📚 Further Reading  There are different antenna gain-combining methods. They are as follows. 1. Equal gain combining (EGC) 2. Maximum ratio combining (MRC) 3. Selective combining (SC) 4. Root mean square gain combining (RMSGC) 5. Zero-Forcing (ZF) Combining  1. Equal gain combining method Equal Gain Combining (EGC) is a diversity combining technique in which the receiver aligns the phase of the received signals from multiple antennas (or channels) but gives them equal amplitude weight before summing. This means each received signal is phase-corrected to be coherent with others, but no scaling is applied based on signal strength or channel quality (unlike MRC). Mathematically, for received signa...

Constellation Diagrams of ASK, PSK, and FSK

📘 Overview of Energy per Bit (Eb / N0) 🧮 Online Simulator for constellation diagrams of ASK, FSK, and PSK 🧮 Theory behind Constellation Diagrams of ASK, FSK, and PSK 🧮 MATLAB Codes for Constellation Diagrams of ASK, FSK, and PSK 📚 Further Reading 📂 Other Topics on Constellation Diagrams of ASK, PSK, and FSK ... 🧮 Simulator for constellation diagrams of m-ary PSK 🧮 Simulator for constellation diagrams of m-ary QAM BASK (Binary ASK) Modulation: Transmits one of two signals: 0 or -√Eb, where Eb​ is the energy per bit. These signals represent binary 0 and 1.    BFSK (Binary FSK) Modulation: Transmits one of two signals: +√Eb​ ( On the y-axis, the phase shift of 90 degrees with respect to the x-axis, which is also termed phase offset ) or √Eb (on x-axis), where Eb​ is the energy per bit. These signals represent binary 0 and 1.  BPSK (Binary PSK) Modulation: Transmits one of two signals...

Comparisons among ASK, PSK, and FSK | And the definitions of each

📘 Comparisons among ASK, FSK, and PSK 🧮 Online Simulator for calculating Bandwidth of ASK, FSK, and PSK 🧮 MATLAB Code for BER vs. SNR Analysis of ASK, FSK, and PSK 📚 Further Reading 📂 View Other Topics on Comparisons among ASK, PSK, and FSK ... 🧮 Comparisons of Noise Sensitivity, Bandwidth, Complexity, etc. 🧮 MATLAB Code for Constellation Diagrams of ASK, FSK, and PSK 🧮 Online Simulator for ASK, FSK, and PSK Generation 🧮 Online Simulator for ASK, FSK, and PSK Constellation 🧮 Some Questions and Answers Modulation ASK, FSK & PSK Constellation MATLAB Simulink MATLAB Code Comparisons among ASK, PSK, and FSK    Comparisons among ASK, PSK, and FSK Comparison among ASK, FSK, and PSK Parameters ASK FSK PSK Variable Characteristics Amplitude Frequency ...

BER performance of QPSK with BPSK, 4-QAM, 16-QAM, 64-QAM, 256-QAM, etc

📘 Overview 📚 QPSK vs BPSK and QAM: A Comparison of Modulation Schemes in Wireless Communication 📚 Real-World Example 🧮 MATLAB Code 📚 Further Reading   QPSK provides twice the data rate compared to BPSK. However, the bit error rate (BER) is approximately the same as BPSK at low SNR values when gray coding is used. On the other hand, QPSK exhibits similar spectral efficiency to 4-QAM and 16-QAM under low SNR conditions. In very noisy channels, QPSK can sometimes achieve better spectral efficiency than 4-QAM or 16-QAM. In practical wireless communication scenarios, QPSK is commonly used along with QAM techniques, especially where adaptive modulation is applied. Modulation Bits/Symbol Points in Constellation Usage Notes BPSK 1 2 Very robust, used in weak signals QPSK 2 4 Balanced speed & reliability 4-QAM ...

MATLAB code for BER vs SNR for M-QAM, M-PSK, QPSk, BPSK, ...

🧮 MATLAB Code for BPSK, M-ary PSK, and M-ary QAM Together 🧮 MATLAB Code for M-ary QAM 🧮 MATLAB Code for M-ary PSK 📚 Further Reading MATLAB Script for BER vs. SNR for M-QAM, M-PSK, QPSK, BPSK % Written by Salim Wireless clc; clear; close all; num_symbols = 1e5; snr_db = -20:2:20; psk_orders = [2, 4, 8, 16, 32]; qam_orders = [4, 16, 64, 256]; ber_psk_results = zeros(length(psk_orders), length(snr_db)); ber_qam_results = zeros(length(qam_orders), length(snr_db)); for i = 1:length(psk_orders) psk_order = psk_orders(i); for j = 1:length(snr_db) data_symbols = randi([0, psk_order-1], 1, num_symbols); modulated_signal = pskmod(data_symbols, psk_order, pi/psk_order); received_signal = awgn(modulated_signal, snr_db(j), 'measured'); demodulated_symbols = pskdemod(received_signal, psk_order, pi/psk_order); ber_psk_results(i, j) = sum(data_symbols ~= demodulated_symbols) / num_symbols; end end for i...

Theoretical vs. simulated BER vs. SNR for ASK, FSK, and PSK

📘 Overview 🧮 Simulator for calculating BER 🧮 MATLAB Codes for calculating theoretical BER 🧮 MATLAB Codes for calculating simulated BER 📚 Further Reading BER vs. SNR denotes how many bits in error are received for a given signal-to-noise ratio, typically measured in dB. Common noise types in wireless systems: 1. Additive White Gaussian Noise (AWGN) 2. Rayleigh Fading AWGN adds random noise; Rayleigh fading attenuates the signal variably. A good SNR helps reduce these effects. Simulator for calculating BER vs SNR for binary ASK, FSK, and PSK Calculate BER for Binary ASK Modulation Enter SNR (dB): Calculate BER Calculate BER for Binary FSK Modulation Enter SNR (dB): Calculate BER Calculate BER for Binary PSK Modulation Enter SNR (dB): Calculate BER BER vs. SNR Curves MATLAB Code for Theoretical BER % The code is written by SalimWireless.Com clc; clear; close all; % SNR va...