0Aliasing happens when a signal is sampled too slowly, causing different signals to look the same. Preventing aliasing helps us keep the true shape of the signal when we digitize it.
1. Choose a sampling rate (Fs) at least twice the highest frequency in the signal (Nyquist rate). 2. Apply a low-pass filter to the signal before sampling to remove frequencies above Fs/2. 3. Sample the filtered signal at the chosen sampling rate.
The Nyquist rate is the minimum sampling rate to avoid aliasing.
Low-pass filtering before sampling is called anti-aliasing filtering.
# Example: Sampling a 1 kHz signal Fs = 3000 # Sampling rate in Hz # Fs is more than 2 * 1000 Hz, so no aliasing expected
# Example: Using a low-pass filter before sampling from scipy.signal import butter, lfilter # Design a low-pass filter with cutoff at 1500 Hz b, a = butter(4, 1500/(Fs/2), btype='low') filtered_signal = lfilter(b, a, original_signal)
This code creates a signal with two frequencies. It applies a low-pass filter before sampling at a rate that would cause aliasing without filtering. The plot shows how filtering helps prevent aliasing.
import numpy as np import matplotlib.pyplot as plt from scipy.signal import butter, lfilter # Create a signal with two frequencies: 500 Hz and 2000 Hz Fs_original = 10000 # Original sampling rate T = 0.01 # seconds # Time vector t = np.linspace(0, T, int(Fs_original * T), endpoint=False) # Signal: sum of 500 Hz and 2000 Hz sine waves signal = np.sin(2 * np.pi * 500 * t) + 0.5 * np.sin(2 * np.pi * 2000 * t) # Sampling rate below Nyquist for 2000 Hz component Fs_sample = 3000 # Sampling rate # Design low-pass filter to prevent aliasing cutoff = Fs_sample / 2 # Nyquist frequency b, a = butter(5, cutoff / (Fs_original / 2), btype='low') filtered_signal = lfilter(b, a, signal) # Downsample by taking every int(Fs_original/Fs_sample) sample step = int(Fs_original / Fs_sample) sampled_signal = filtered_signal[::step] # Time vector for sampled signal t_sampled = t[::step] # Plot original, filtered, and sampled signals plt.figure(figsize=(10, 6)) plt.plot(t, signal, label='Original Signal') plt.plot(t, filtered_signal, label='Filtered Signal') plt.stem(t_sampled, sampled_signal, linefmt='r-', markerfmt='ro', basefmt='k-', label='Sampled Signal', use_line_collection=True) plt.xlabel('Time [s]') plt.ylabel('Amplitude') plt.title('Aliasing and Prevention by Low-pass Filtering') plt.legend() plt.tight_layout() plt.show()
Always sample at least twice the highest frequency to avoid aliasing.
Anti-aliasing filters remove frequencies that cause aliasing before sampling.
Aliasing can cause wrong data interpretation, so prevention is important.
Aliasing happens when sampling too slowly, mixing up signal frequencies.
Use a sampling rate at least twice the highest frequency (Nyquist rate).
Apply a low-pass filter before sampling to remove high frequencies (anti-aliasing).
import numpy as np fs = 800 # sampling frequency f = 400 # signal frequency t = np.arange(0, 1, 1/fs) signal = np.sin(2 * np.pi * f * t) print(signal[:5])What is the main issue with this sampling?
def sample_signal(signal, fs, cutoff):
# Apply low-pass filter
filtered = signal[signal < cutoff]
# Sample at fs
samples = filtered[::int(1000/fs)]
return samples
signal = [100, 200, 300, 400, 500]
fs = 500
cutoff = 250
print(sample_signal(signal, fs, cutoff))
What is the main error?Pass the test and this pattern moves unknown strong while earning +1 Grip toward your certificate.