Note: This page has been translated by MathWorks. Please click here

To view all translated materals including this page, select Japan from the country navigator on the bottom of this page.

To view all translated materals including this page, select Japan from the country navigator on the bottom of this page.

Lomb-Scargle periodogram

```
[pxx,f]
= plomb(x,t)
```

```
[pxx,f]
= plomb(x,fs)
```

```
[pxx,f]
= plomb(___,fmax)
```

```
[pxx,f]
= plomb(___,fmax,ofac)
```

```
[pxx,fvec]
= plomb(___,fvec)
```

`[___] = plomb(___,spectrumtype)`

`[___,pth] = plomb(___,'Pd',pdvec)`

```
[pxx,w]
= plomb(x)
```

`plomb(___)`

`[`

returns the Lomb-Scargle power spectral density (PSD) estimate,
`pxx`

,`f`

]
= plomb(`x`

,`t`

)`pxx`

, of a signal, `x`

, that is
sampled at the instants specified in `t`

.
`t`

must increase monotonically but need not be uniformly
spaced. All elements of `t`

must be nonnegative.
`pxx`

is evaluated at the frequencies returned in
`f`

.

If

`x`

is a vector, it is treated as a single channel.If

`x`

is a matrix, then`plomb`

computes the PSD independently for each column and returns it in the corresponding column of`pxx`

.

`x`

or `t`

can contain `NaN`

s
or `NaT`

s. These values are treated as missing data
and excluded from the spectrum computation.

`[`

estimates
the PSD up to a maximum frequency, `pxx`

,`f`

]
= plomb(___,`fmax`

)`fmax`

, using
any of the input arguments from previous syntaxes. If the signal is
sampled at *N* non-`NaN`

instants,
and Δ*t* is
the time difference between the first and the last of them, then `pxx`

is
returned at `round`

(`fmax`

/ *f*_{min}) points,
where *f*_{min} = 1/(4 × *N* × *t _{s}*) is
the smallest frequency at which

`pxx`

is computed
and the average sample time is `fmax`

defaults
to 1/(2 × `[`

specifies
an integer oversampling factor, `pxx`

,`f`

]
= plomb(___,`fmax`

,`ofac`

)`ofac`

. The use
of `ofac`

to interpolate or smooth a spectrum resembles
the zero-padding technique for FFT-based methods. `pxx`

is
again returned at `round`

(`fmax`

/*f*_{min})
frequency points, but the minimum frequency considered in this case
is 1/(`ofac`

× *N* × *t _{s}*).

`ofac`

defaults
to 4.`[___,`

returns
the power-level threshold, `pth`

] = plomb(___,'Pd',`pdvec`

)`pth`

, such that a peak
with a value larger than `pth`

has a probability `pdvec`

of
being a true signal peak and not the result of random fluctuations. `pdvec`

can
be a vector. Every element of `pdvec`

must be greater
than 0 and smaller than 1. Each row of `pth`

corresponds
to an element of `pdvec`

. `pth`

has
the same number of channels as `x`

. This option
is not available if you specify the output frequencies in `fvec`

.

`[`

returns the PSD estimate
of `pxx`

,`w`

]
= plomb(`x`

)`x`

evaluated at a set of evenly spaced normalized
frequencies, `w`

, spanning the Nyquist interval.
Use `NaN`

s to specify missing samples. All of the
above options are available for normalized frequencies. To access
them, specify an empty array as the second input.

`plomb(___)`

with no output arguments
plots the Lomb-Scargle periodogram PSD estimate in the current figure
window.

[1] Lomb, Nicholas R. “Least-Squares Frequency
Analysis of Unequally Spaced Data.” *Astrophysics
and Space Science.* Vol. 39, 1976, pp. 447–462.

[2] Scargle, Jeffrey D. “Studies in Astronomical
Time Series Analysis. II. Statistical Aspects of Spectral Analysis
of Unevenly Spaced Data.” *Astrophysical Journal.* Vol.
263, 1982, pp. 835–853.

[3] Press, William H., and George B. Rybicki.
“Fast Algorithm for Spectral Analysis of Unevenly Sampled Data.” *Astrophysical
Journal.* Vol. 338, 1989, pp. 277–280.

[4] Horne, James H., and Sallie L. Baliunas. “A Prescription
for Period Analysis of Unevenly Sampled Time Series.” *Astrophysical
Journal.* Vol. 302, 1986, pp. 757–763.

`bandpower`

| `pburg`

| `pcov`

| `peig`

| `periodogram`

| `pmcov`

| `pmtm`

| `pmusic`

| `pwelch`

| `pyulear`

| `spectrogram`

Was this topic helpful?