Calculating power spectral density from FFT

  • Context: Undergrad 
  • Thread starter Thread starter GhostLoveScore
  • Start date Start date
  • Tags Tags
    Density Fft Power
Click For Summary

Discussion Overview

The discussion revolves around calculating power spectral density from Fast Fourier Transform (FFT) data, specifically in the context of amateur radio astronomy using an RTL SDR stick. Participants explore the relationship between I/Q data and the FFT output, and how to derive power spectral density from the FFT results.

Discussion Character

  • Technical explanation
  • Exploratory

Main Points Raised

  • One participant expresses interest in understanding the process of calculating power spectral density from FFT, particularly in the context of amateur radio astronomy.
  • Another participant explains that the FFT output is complex and that the power spectrum can be obtained by squaring the magnitude of the FFT output, represented as Z * conj(Z) = abs(Z).^2 = I^2 + Q^2.
  • A question is raised about the distinction between the I and Q in the context of the FFT output versus the raw I/Q data from the SDR.
  • A clarification is provided that in the FFT output, I and Q refer to the real and imaginary components of the FFT output, not the same I and Q from the raw data.

Areas of Agreement / Disagreement

Participants generally agree on the method to calculate power spectral density from the FFT output, but there is some clarification needed regarding the terminology and the distinction between I/Q data and FFT components.

Contextual Notes

There may be limitations in understanding the notation and the relationship between raw I/Q data and FFT output, which could affect the clarity of the discussion.

GhostLoveScore
Messages
149
Reaction score
9
EDIT: Sorry. It's FFT - Fast Fourier Transform, not FTT.
I am interested in doing some amateur radio astronomy. Mainly at 1420MHz, hydrogen line. I have a RTL SDR stick. For those who don't know what that is, it's USB DVB-T receiver that can receive anything between 24 – 1766 MHz.

Now, there is a finished program for what I am asking, but I don't like using something that I don't know how it works. This is the source for finished program https://github.com/mariocannistra/radio-astronomy-fftw

SDR stick outputs I/Q data. I/Q data is just amplitude data but the amplitude is sampled at the same time 90 degrees away from the other sample. Amplitude of the signal is sqrt(I^2+Q^2). (I don't know how to insert equation, sorry). No problem so far.

I can record a sample of the signal and there is a library for Raspberry PI for doing Fast Fourier Transform on that sample. This is where I got stuck. How to find power spectral density from the Fast Fourier Transform?

I got the impression that it's just amplitude squared - I^2+Q^2? Can it be that simple?

I apologize for any nonsense that I may have said here. And of course I wrote FTT instead of FFT in the thread title. If I wrote everything correctly that would be weird. Moderators, please correct that.
 
Physics news on Phys.org
The FFT gives a complex output which is basically Z = I + Qi. So to generate the power spectrum you take Z * conj(Z) = abs(Z).^2 = I^2 + Q^2.
If you're taking the FFT of a real input signal, then the positive and negative frequency parts have equal power, so you can just plot the positive frequency power spectrum and multiply by 2.
 
olivermsun said:
The FFT gives a complex output which is basically Z = I + Qi. So to generate the power spectrum you take Z * conj(Z) = abs(Z).^2 = I^2 + Q^2.
If you're taking the FFT of a real input signal, then the positive and negative frequency parts have equal power, so you can just plot the positive frequency power spectrum and multiply by 2.

I assume that in Z=I + Qi, I and Q are not the same variables as I, Q in raw data received from SDR? In this case Z= I + Qi, I and Q are just real and imaginary components of FFT output?
 
Sorry for the bad notation. Yes, I and Q are real and imaginary components of the FFT output, but they are basically just the coefficients (not the full oscillating signals as in the raw I/Q data outputs).
 

Similar threads

  • · Replies 28 ·
Replies
28
Views
3K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 11 ·
Replies
11
Views
3K
Replies
5
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
Replies
1
Views
2K
  • · Replies 10 ·
Replies
10
Views
3K