×
INTELLIGENT WORK FORUMS
FOR ENGINEERING PROFESSIONALS

Contact US

Log In

Come Join Us!

Are you an
Engineering professional?
Join Eng-Tips Forums!
  • Talk With Other Members
  • Be Notified Of Responses
    To Your Posts
  • Keyword Search
  • One-Click Access To Your
    Favorite Forums
  • Automated Signatures
    On Your Posts
  • Best Of All, It's Free!

*Eng-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

Students Click Here

Lomb Scargle periodogram for large bandwidth signals?

Lomb Scargle periodogram for large bandwidth signals?

Lomb Scargle periodogram for large bandwidth signals?

(OP)
Dear all,

I have to deal with an unevenly sampled time history that is related to a white noise signal.

Now, I can not a priori use FFT methods as the uneven sampling is introducing artifacts that hide the real signal.
I was thinking about using Lomb Scargle processing.

Whereas the preliminary results are pretty good for single frequency component signals, I see that for my case when I turn my input signals into random the things go worse,as the amplitude of the Periodogram seems to go down proportionally with the amount of significant frequency components.

i wounder if somebody already had this issue, and has suggestion on the most adequate chain to be adopted for random signals that are also unevenly sampled.

thanks in advance for your help!
Edo

RE: Lomb Scargle periodogram for large bandwidth signals?

"as the amplitude of the Periodogram seems to go down proportionally with the amount of significant frequency components.
"

is that not what you'd expect?


Here's what i'd do

1 synthesise some white noise with equal time sampling

1a analyse it with fft

1b analyse it with your LSSA code

1c do 1a and 1b agree?

2 extract interpolated signals from signal 1 at your uneven sampling intervals

2a analyse it with your LSSA code


2b do 2a and 1b agree well enough?


3 analyse your actual data with your LSSA code

Cheers

Greg Locock


New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?

RE: Lomb Scargle periodogram for large bandwidth signals?

(OP)
hello

@IRstuff: as of my understanding and of the documentation I have available and is providing examples, a periodogram is able to identify more than 1 frequency component.
@GregLockck: what you indicate is the procedure I followed in order to compare and try to understand the LombScargle "animal".. and I find my self at the point of asking the community :)

What I understood is that
1) LombScargle is a "Power Spectrum" tool (spectrum in "units²"), whereas FFT is an "amplitude" (spectrum in "units") tool.

2) the higher the amount of spectral components (for a given constant time window and sampling), the LombScargle spectrum decreases, whereas I would expect that the amplitudes are also estimated. Therefore I understand that the spectrum restituted by the LombScargle is assessing the product spectral componenent amplitude*probability. But for this I did not find any confirmation or statement.
eg.
1 frequency component with amplitude 1-> 100% probability* amplitude 1
2 frequency components each with amplitude 1-> 50% probability* ampitude 1
3 frequency components each with amplitude 1-> 33% probability* ampitude 1
4 frequency components each with amplitude 1-> 25% probability* amplitude 1

3) What is also unclear to me is if LombScargle is a Power Spectral Density (units²/Hz) or it it "just" a Power Spectrum (units²). But in all documentation I have gone through, I don't find a clear statement ...

thanks in advance for feedbacks,
Edo

RE: Lomb Scargle periodogram for large bandwidth signals?

(OP)
yep. Correct, IRstuff.
and white noise is generated as a train of frequency components by creating the spectrum and then backtransforming it into the time domain and resampling it.
as stated also by GregLocock..

RE: Lomb Scargle periodogram for large bandwidth signals?

I don't read that from Greg's posting. You cannot get a "train" of frequency components, because there isn't any. Random noise is RANDOM; at any given instance, you get random parts of the PSD that randomly sum up to the random time domain signal. Pick a different window in time, and you get a different answer. You will NEVER see the PSD from any single periodogram or FFT, unless it's infinitely long. The PSD will only be as a consequence of collecting gobs of data and superimposing each individual measured spectrum and enveloping that aggregate.

TTFN
FAQ731-376: Eng-Tips.com Forum Policies

Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

Red Flag Submitted

Thank you for helping keep Eng-Tips Forums free from inappropriate posts.
The Eng-Tips staff will check this out and take appropriate action.

Reply To This Thread

Posting in the Eng-Tips forums is a member-only feature.

Click Here to join Eng-Tips and talk with other members! Already a Member? Login



News


Close Box

Join Eng-Tips® Today!

Join your peers on the Internet's largest technical engineering professional community.
It's easy to join and it's free.

Here's Why Members Love Eng-Tips Forums:

Register now while it's still free!

Already a member? Close this window and log in.

Join Us             Close