SignalProcessing - Maple Programming Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Science and Engineering : Signal Processing : Visualization : SignalProcessing/LSPeriodogram

SignalProcessing

  

LSPeriodogram

  

plot the Lomb-Scargle periodogram of a signal

 

Calling Sequence

Parameters

Options

Description

Examples

References

Compatibility

Calling Sequence

LSPeriodogram(t, f, options)

Parameters

t

-

rtable ; time

f

-

rtable ; signal

options

-

(optional) one or more of the options listed below

Options

• 

detrend : truefalse : indicates if the data should have any trend removed

• 

frequencyscale : anything 

• 

powerscale : identical("dB", "absolute") 

• 

timeunit : anything 

• 

outputdata : posint : specify the number of the output data points

• 

errors : rtable or realcons 

• 

precise : truefalse : indicates the methods to use

• 

minimumfrequency : float : indicates minimum frequency

• 

maximumfrequency : float : indicates maximum frequency

• 

normalization : identical('standard', 'model', 'log', 'psd') : indicates the normalization methods

• 

oversamplingfactor : numeric : indicates the oversampling factor

• 

nyquistfactor : numeric : indicates the multiple of the average Nyquist frequency

• 

centerdata : truefalse : indicates whether the data needs to subtract the weighted mean

• 

generalized : truefalse : indicates the methods to use for the Lomb Scargle Fast Approximation

Description

• 

The LSPeriodogram(t, f) command plots the power spectrum for one signal. t and f must be a one-dimensional rtable.

• 

The detrend option specifies if the data should have any linear trend removed before being plotted.

• 

The frequencyscale option controls the scaling used on the frequency axis. The default is Hz.

• 

The powerscale option controls the scaling used on the power axis and can be one of "dB" or "absolute". The default is "dB".

• 

The timeunit option specifies the unit of the input time. The default value is s.

• 

The errors option specifies the expected magnitude of the errors, which is used to get the normalized weight for generalized Lomb Scargle periodogram. It can only be a real number or an one-dimensional array. The default value is NULL.

• 

The outputdata option specifies the number of output data points. These output data points correspond to equally spaced frequencies.

• 

The minimumfrequency option specifies the minimum frequency.

• 

The maximumfrequency option specifies the maximum frequency.

• 

If one of maximumfrequency, minimumfrequency and outputdata is specified, then all of them must be provided. If all of them are not provided, then the frequencies used run from 1Toversamplingfactor to nyquistfactorNf, where T=maxtmint and Nf=n2T where n is the number of elements in t. Output frequencies are equally spaced.  

• 

The nyquistfactor option has a default value 1.

• 

The normalization option controls the normalization used on the power and can be one of standard, model, log or psd. The default value is standard. Standard Normalization : Pstandardf=χref2χf2χref2. Model Normalization : Pmodelf=χref2χf2χf2. Logarithmic Normalization:  Plogf=logχref2χf2. PSD Normalization : Ppsdf=χref22χf22. In all normalization methods, let fi be the measurement in time ti with errors σi, then χf2=i=1nfifti2σi2 where n is the number of elements in f, and χref2 is the non-varying reference model which is the χf2 for weighted mean.

• 

The oversamplingfactor option specifies the oversampling factor. The default value is 5.

• 

If the centerdata=false option is provided, then the weighted mean will not be subtracted from signal data. The default value is true.

• 

If the generalized=false option is provided, then the generalized Lomb Scargle Periodogram will not be used. The default value is true.

• 

Additional plotting options as described on the plot/option help page may be included.

Examples

withSignalProcessing:

Generate a time vector

tsortLinearAlgebra:-RandomVector210,generator=0..12Pi,datatype=float8

t0.01969308700552850.04339383144798540.1747061339360980.2531613855537030.2948181219160560.3896813310069070.4486974492913810.5007648600489720.5806959187938450.5838508814800711024 element Vector[column]

(1)

Assume the time unit is second and generate signal using frequencies of 1 Hz and 2 Hz:

f11.0:

f22.0:

sVector210,i→sinf12Piti+1.5sinf22Piti,datatype=float8

s0.4908484649532151.047301634313102.106939846966700.9402276849067050.1597684131121270.8355791393708700.5846321103955350.009611312183662660.7878914338841810.8010748377000021024 element Vector[column]

(2)

LSPeriodogramt,s,frequencyscale=kHz

LSPeriodogramt,s,frequencyscale=Hz

Assume the time unit is day and the frequencies are 1 cycle per day and 2 cycles per day, which are equivalent to 0.000011574 Hz and 0.000023148 Hz

f3124.6060

f30.00001157407407

(3)

f4224.6060

f40.00002314814814

(4)

LSPeriodogramt,s,timeunit=d,frequencyscale=Hz

LSPeriodogramt,s,timeunit=d,frequencyscale=Hz,size=default,square,font=Times,1

References

  

Lomb, Nicholas R. "Least-Squares Frequency Analysis of Unequally Spaced Data." Astrophysics and Space Science. Vol. 39, 1976, pp. 447–462.

  

Press, William H., and Rybicki, George B. "Fast Algorithm for Spectral Analysis of Unevenly Sampled Data." Astrophysical Journal. Vol. 338, 1989, pp. 277–280.

  

Zechmeister, M., and Kürster, M. "The generalised Lomb-Scargle periodogram. A new formalism for the floating-mean and Keplerian periodograms." Astronomy & Astrophysics. Vol. 496, 2009, pp. 577-584

Compatibility

• 

The SignalProcessing[LSPeriodogram] command was introduced in Maple 2019.

• 

For more information on Maple 2019 changes, see Updates in Maple 2019.

See Also

Periodogram

SignalPlot

SignalProcessing

Spectrogram

Statistics:-Detrend