Getting the power spectral density from a plot

  • Context: Graduate 
  • Thread starter Thread starter Malamala
  • Start date Start date
Click For Summary

Discussion Overview

The discussion revolves around obtaining the power spectral density (PSD) from time series data collected from a self-heterodyne measurement involving a laser. Participants explore various methods and transformations that could be applied to the data, including the Welch method and FFT, while addressing concerns about windowing and segment length.

Discussion Character

  • Exploratory, Technical explanation, Debate/contested, Mathematical reasoning

Main Points Raised

  • One participant describes the process of obtaining the PSD from a laser measurement and expresses uncertainty about the necessary transformations, such as windowing and detrending.
  • Another participant suggests using the Welch method to compute the PSD and emphasizes the importance of considering window length.
  • A later reply reiterates the suggestion of the Welch method and highlights that the visual output of the PSD can vary significantly based on the segment length parameter.
  • One participant proposes trying the FFT with a typical 256-point window and squaring the frequency values as an alternative method for finding the PSD.

Areas of Agreement / Disagreement

Participants generally agree on the use of the Welch method and FFT for computing the PSD, but there is no consensus on the optimal parameters for windowing and segment length, leading to differing opinions on the best approach.

Contextual Notes

The discussion includes uncertainty regarding the effects of different transformations on the PSD results, and participants express varying levels of familiarity with the computational methods involved.

Who May Find This Useful

Individuals interested in signal processing, particularly in the context of laser measurements and power spectral density analysis, may find this discussion relevant.

Malamala
Messages
348
Reaction score
28
Hello! I have some time series from a self-heterodyne measurement. Basically i have a laser and I combine it with a delayed version of itself, shifted by ##\Omega = 2\pi \times 80## MHz. The delay is ##\tau = 100##ns. I remove the DC so the signal from the beat note is:

$$V(t) = V_0\cos(\Omega t + \phi(t+\tau)-\phi(t))$$
where ##V## is the oscilloscope signal (in V), and ##\phi(t)## is the phase noise of the laser. What I am interested in, is getting the power spectral density (PSD) of this ##\phi(t+\tau)-\phi(t)## term. I have the data on my laptop (it was recorded for ##2## ms, with a sampling rate of ##500## mega sample per second). I understand that, ideally, I just need to multiply this by ##sin(\Omega)## and ##cos(\Omega)## and get the in-phase (I) and quadrature (Q) signal, then compute ##arctan(Q/I)##, and then what I need is just the FFT of this. But after reading online more there seems to be other things that I might need to do, like windowings, applying a hilbert transform, detrending, welch transforms and doing some and not the others change the result. And I am not sure which transformations I need to include (or how can I figure out which I should include). I've never done this in practice before and it's a bit overwhelming. Can someone guide me a bit on how to proceed?
 
Physics news on Phys.org
Hi. The Welch method computes the PSD. You will need to consider your window length. Depending on the software you are using, it can be a rather quick and easy command. What computational software do you have?
 
Last edited:
osilmag said:
Hi. The Welch method computes the PSD. You will need to consider your window length. Depending on the software you are using, it can be a rather quick and easy command. What computational software do you have?
I am using python. Welch method is built it, but it has this parameter which defines how many points I have in each segment and the obtained PSD changes visually based on that. And it only looks like what I would expect, based on other papers, only when I include all the points in my data set in the segment (i.e. I have only one segment). So I am not sure if that is correct.
 
You could try the FFT with a typical 256 pt window and then square the value of each frequency. That is kind of the original way of finding it. Compare how that looks with the welch method.
 

Similar threads

  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 0 ·
Replies
0
Views
2K
  • · Replies 6 ·
Replies
6
Views
3K
  • · Replies 1 ·
Replies
1
Views
4K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 8 ·
Replies
8
Views
3K
Replies
9
Views
2K
  • · Replies 7 ·
Replies
7
Views
3K
Replies
1
Views
2K
Replies
1
Views
2K