plomb
Lomb-Scargle 주기도
구문
설명
[
는 pxx
,f
] = plomb(x
,t
)t
에 지정된 시점에 샘플링된 신호 x
의 Lomb-Scargle 파워 스펙트럼 밀도(PSD) 추정값 pxx
를 반환합니다. t
는 단조 증가(Monotonically Increasing)해야 하지만, 간격이 균일할 필요는 없습니다. t
의 모든 요소는 음수가 아니어야 합니다. pxx
는 f
로 반환된 주파수에서 계산됩니다.
x
가 벡터이면 단일채널로 처리됩니다.x
가 행렬이면plomb
는 각 열의 PSD를 개별적으로 계산하고 그 값을pxx
의 해당 열로 반환합니다.
x
또는 t
에는 NaN
또는 NaT
가 포함될 수 있습니다. 단, 이러한 값은 누락된 데이터로 처리되므로 스펙트럼 계산에서 제외됩니다.
[
는 위 구문의 입력 인수를 사용하여 최대 주파수 pxx
,f
] = plomb(___,fmax
)fmax
까지 PSD를 추정합니다. 신호가 NaN
이 아닌 N개의 시점에서 샘플링되고 Δt가 첫 번째 시점과 마지막 시점의 시간 차분인 경우, pxx
는 round
(fmax
/fmin)개의 지점에서 반환됩니다. 여기서 fmin = 1/(4 × N × ts)는 pxx
가 계산된 가장 작은 주파수이고, 평균 샘플 시간은 ts = Δt/(N – 1)입니다. fmax
는 디폴트 값이 1/(2 × ts)이며, 균일하게 샘플링된 신호인 경우 나이퀴스트 주파수에 해당합니다.
plomb(___)
(출력 인수 없음)는 현재 Figure 창에 Lomb-Scargle 주기도 PSD 추정값을 플로팅합니다.
예제
입력 인수
출력 인수
세부 정보
참고 문헌
[1] 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.
[2] Lomb, Nicholas R. "Least-Squares Frequency Analysis of Unequally Spaced Data." Astrophysics and Space Science. Vol. 39, 1976, pp. 447–462.
[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] 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.