Calculating power spectral density from FFT

  • Thread starter Thread starter GhostLoveScore
  • Start date Start date
  • Tags Tags
    Density Fft Power
Click For Summary
The discussion centers on calculating power spectral density (PSD) from Fast Fourier Transform (FFT) data in the context of amateur radio astronomy at 1420 MHz. The user is utilizing an RTL SDR stick to capture I/Q data, which represents amplitude sampled 90 degrees apart. To derive the PSD from the FFT output, it is clarified that the power spectrum can be obtained by squaring the magnitude of the complex output, specifically using the formula Z * conj(Z) = abs(Z)^2 = I^2 + Q^2. Additionally, when analyzing a real input signal, the positive and negative frequency components have equal power, allowing for the positive frequency spectrum to be plotted and multiplied by 2. The distinction between the I/Q variables in raw data and those in the FFT output is also confirmed, emphasizing that they represent different aspects of the data.
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 4 ·
Replies
4
Views
2K
Replies
5
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
Replies
1
Views
2K
  • · Replies 10 ·
Replies
10
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K