Frequency spectrum interpretation in Signal Processing - Step-by-Step Execution

Concept Flow - Frequency spectrum interpretation

Input Signal

↓

Apply Fourier Transform

↓

Obtain Frequency Spectrum

↓

Analyze Magnitude & Phase

↓

Identify Dominant Frequencies

↓

Interpret Signal Characteristics

The process starts with an input signal, applies Fourier Transform to get the frequency spectrum, then analyzes magnitude and phase to identify dominant frequencies and interpret the signal.

Execution Sample

Signal Processing

import numpy as np
import matplotlib.pyplot as plt
t = np.linspace(0, 1, 1000, endpoint=False)
signal = np.sin(2*np.pi*5*t) + 0.5*np.sin(2*np.pi*10*t)
frequency_spectrum = np.fft.fft(signal)
magnitude = np.abs(frequency_spectrum)

This code creates a signal with two sine waves at 5 Hz and 10 Hz, computes its frequency spectrum using FFT, and calculates the magnitude of each frequency component.

Execution Table

Step	Action	Input/Variable	Output/Result	Explanation
1	Create time vector t	t = np.linspace(0,1,1000, endpoint=False)	t array from 0 to 1 sec	Defines time points for signal sampling
2	Generate signal	signal = sin(2π5t) + 0.5sin(2π10t)	Signal array with combined sine waves	Signal contains 5 Hz and 10 Hz components
3	Apply FFT	frequency_spectrum = fft(signal)	Complex frequency spectrum array	Transforms signal to frequency domain
4	Calculate magnitude	magnitude = abs(frequency_spectrum)	Magnitude array of frequencies	Shows strength of each frequency
5	Identify peaks	Look at magnitude array	Peaks at indices ~5 and ~10	Dominant frequencies at 5 Hz and 10 Hz
6	Stop	N/A	N/A	All steps completed

💡 All steps executed to transform signal and identify dominant frequencies

Variable Tracker

Variable	Start	After Step 2	After Step 3	After Step 4	Final
t	undefined	array 0 to 1 sec	array 0 to 1 sec	array 0 to 1 sec	array 0 to 1 sec
signal	undefined	signal with 5Hz & 10Hz	signal with 5Hz & 10Hz	signal with 5Hz & 10Hz	signal with 5Hz & 10Hz
frequency_spectrum	undefined	undefined	complex array	complex array	complex array
magnitude	undefined	undefined	undefined	magnitude array	magnitude array

Key Moments - 3 Insights

Why does the frequency spectrum have complex numbers?

The FFT output is complex because it encodes both amplitude and phase of each frequency component, as shown in step 3 of the execution_table.

How do we know which frequencies are dominant from the magnitude array?

Dominant frequencies correspond to peaks in the magnitude array, as seen in step 5 where peaks near indices 5 and 10 indicate 5 Hz and 10 Hz components.

Why do we use the absolute value of the FFT output?

Taking the absolute value gives the magnitude (strength) of each frequency component, ignoring phase, which helps identify dominant frequencies (step 4).

Visual Quiz - 3 Questions

Look at the execution_table, what is the output after applying FFT (step 3)?

A A real-valued array representing signal amplitude

B A complex array representing frequency components

C A time vector from 0 to 1 second

D A magnitude array of frequencies

💡 Hint

Refer to step 3 in execution_table where FFT output is described as a complex array.

At which step do we identify the dominant frequencies in the signal?

A Step 5

B Step 2

C Step 3

D Step 4

💡 Hint

Step 5 in execution_table mentions identifying peaks in the magnitude array.

If the signal had a stronger 10 Hz component, how would the magnitude array change?

A The magnitude array would become complex

B The peak near index 5 would be higher

C The peak near index 10 would be higher

D The time vector t would change

💡 Hint

Refer to variable_tracker and execution_table steps 4 and 5 about magnitude peaks corresponding to frequency strength.

Concept Snapshot

Frequency Spectrum Interpretation:
- Input signal is transformed using Fourier Transform (FFT).
- FFT output is complex, encoding amplitude and phase.
- Magnitude (absolute value) shows strength of frequencies.
- Peaks in magnitude identify dominant frequencies.
- Used to analyze signal characteristics in frequency domain.

Full Transcript

Frequency spectrum interpretation involves converting a time-based signal into its frequency components using Fourier Transform. The process starts by defining a time vector and creating a signal composed of sine waves at different frequencies. Applying FFT transforms the signal into a complex frequency spectrum, which contains both amplitude and phase information. Taking the absolute value of this spectrum gives the magnitude, showing how strong each frequency is. By examining peaks in the magnitude array, we identify dominant frequencies in the signal. This helps us understand the signal's characteristics in the frequency domain.

Practice

1. What does a peak in a frequency spectrum represent in signal processing?

easy

A. A strong presence of a specific frequency in the signal

B. The time duration of the signal

C. The amplitude of the original time signal

D. The noise level in the signal

Solution

1. Step 1: Understand frequency spectrum meaning

   The frequency spectrum shows how much of each frequency is present in the signal.

2. Step 2: Interpret peaks in the spectrum

   A peak means that frequency has a strong presence or high amplitude in the signal.

3. Final Answer:

   A strong presence of a specific frequency in the signal -> Option A

4. Quick Check:

   Peak = strong frequency presence [OK]

Hint: Peaks show strong frequencies in the signal [OK]

Common Mistakes:

• Confusing time duration with frequency
• Thinking amplitude in time domain equals peak in spectrum
• Assuming noise causes peaks

2. Which of the following is the correct way to compute the frequency spectrum of a signal in Python using numpy?

easy

A. np.fft.freq(signal)

B. np.fft.ifft(signal)

C. np.fft.fft(signal)

D. np.fft.spectrum(signal)

Solution

1. Step 1: Recall numpy FFT functions

   np.fft.fft() computes the Fourier Transform to get frequency spectrum.

2. Step 2: Check other options

   np.fft.ifft() computes inverse FFT, others do not exist.

3. Final Answer:

   np.fft.fft(signal) -> Option C

4. Quick Check:

   FFT = np.fft.fft() [OK]

Hint: Use np.fft.fft() for frequency spectrum [OK]

Common Mistakes:

• Using inverse FFT instead of FFT
• Using non-existent numpy functions
• Confusing frequency and time domain functions

3. Given the Python code:

import numpy as np
signal = np.sin(2 * np.pi * 5 * np.linspace(0, 1, 100))
spectrum = np.fft.fft(signal)
peak_freq_index = np.argmax(np.abs(spectrum))
print(peak_freq_index)

What is the printed output?

medium

A. 50

B. 95

C. 0

D. 5

Solution

1. Step 1: Understand the signal frequency

   The signal is a sine wave at 5 Hz sampled 100 times over 1 second.

2. Step 2: Find peak frequency index in FFT

   FFT peak index corresponds to frequency 5, so index 5 is max magnitude.

3. Final Answer:

   5 -> Option D

4. Quick Check:

   Peak index = frequency = 5 [OK]

Hint: Peak index matches sine wave frequency [OK]

Common Mistakes:

• Confusing index with frequency value
• Ignoring sampling rate and length
• Mistaking negative frequency indices

4. The following code attempts to find the main frequency in a signal but gives wrong results:

import numpy as np
signal = np.sin(2 * np.pi * 10 * np.linspace(0, 1, 100))
spectrum = np.fft.fft(signal)
peak = np.argmax(spectrum)
print(peak)

What is the error?

medium

A. Using np.fft.fft instead of np.fft.ifft

B. Using np.argmax on complex spectrum instead of magnitude

C. Signal frequency is too high for sampling rate

D. Not normalizing the spectrum before argmax

Solution

1. Step 1: Check how peak is found

   np.argmax(spectrum) finds max of complex numbers, which is incorrect.

2. Step 2: Correct method to find peak

   Use np.argmax(np.abs(spectrum)) to find peak magnitude index.

3. Final Answer:

   Using np.argmax on complex spectrum instead of magnitude -> Option B

4. Quick Check:

   Peak must be on magnitude, not complex values [OK]

Hint: Use magnitude (abs) before argmax [OK]

Common Mistakes:

• Argmax on complex array directly
• Confusing FFT and IFFT
• Ignoring magnitude for peak detection

5. You have a signal composed of two sine waves at 3 Hz and 7 Hz sampled at 50 Hz for 2 seconds. Which code snippet correctly identifies both frequencies from the frequency spectrum?

hard

A. import numpy as np
t = np.linspace(0, 2, 100)
signal = np.sin(2*np.pi*3*t) + np.sin(2*np.pi*7*t)
spectrum = np.fft.fft(signal)
freqs = np.fft.fftfreq(len(t), d=1/50)
peaks = freqs[np.argsort(np.abs(spectrum))[-2:]]
print(np.sort(np.abs(peaks)))

B. import numpy as np
t = np.linspace(0, 2, 100)
signal = np.sin(2*np.pi*3*t) + np.sin(2*np.pi*7*t)
spectrum = np.fft.ifft(signal)
freqs = np.fft.fftfreq(len(t), d=1/50)
peaks = freqs[np.argmax(np.abs(spectrum))]
print(peaks)

C. import numpy as np
t = np.linspace(0, 2, 50)
signal = np.sin(2*np.pi*3*t) + np.sin(2*np.pi*7*t)
spectrum = np.fft.fft(signal)
freqs = np.fft.fftfreq(len(t), d=1/50)
peaks = freqs[np.argsort(np.abs(spectrum))[-2:]]
print(np.sort(np.abs(peaks)))

D. import numpy as np
t = np.linspace(0, 2, 100)
signal = np.sin(2*np.pi*3*t) + np.sin(2*np.pi*7*t)
spectrum = np.fft.fft(signal)
freqs = np.fft.fftfreq(len(t), d=1/100)
peaks = freqs[np.argsort(np.abs(spectrum))[-2:]]
print(np.sort(np.abs(peaks)))

Solution

1. Step 1: Check sampling and signal length

   Sampling at 50 Hz for 2 seconds means 100 samples, so t has 100 points.

2. Step 2: Verify FFT and frequency calculation

   Use np.fft.fft(signal) and np.fft.fftfreq with d=1/50 for correct frequency bins.

3. Step 3: Identify two largest peaks

   Sort magnitudes and pick last two indices for peaks, then print sorted absolute frequencies.

4. Final Answer:

   Code snippet D correctly identifies both frequencies -> Option A

5. Quick Check:

   Correct sampling, FFT, freq bins, and peak selection [OK]

Hint: Match sample count and rate, use fftfreq with correct d [OK]

Common Mistakes:

• Using inverse FFT instead of FFT
• Wrong sample count or sampling interval
• Incorrect frequency bin calculation