Периодограмма Lomb-Scargle
[
возвращает оценку степени спектральной плотности (PSD) Lomb-Scargle, pxx
,f
]
= plomb(x
,t
)pxx
, из сигнала, x
, это производится в моменты, заданные в t
T
должен увеличиться монотонно, но нуждаться быть расположенными неравными интервалами. Все элементы t
mustBeNonnegative. pxx
оценен на частотах, возвращенных в f
.
Если x
вектор, он обработан как один канал.
Если x
матрица, затем plomb
вычисляет PSD независимо для каждого столбца и возвращает его в соответствующем столбце pxx
.
x
или t
может содержать NaN
s или NaT
s. Эти значения обработаны как недостающие данные и исключены из расчета спектра.
[
оценивает PSD до максимальной частоты, pxx
,f
]
= plomb(___,fmax
)fmax
, использование любого из входных параметров от предыдущих синтаксисов. Если сигнал производится в N non-NaN
моменты и Δt являются разницей во времени между первым и последним из них, затем pxx
возвращен в round
(fmax
/ min f) точки, где min f = 1 / (4 × N × ts) является наименьшей частотой в который pxx
вычисляется и средним шагом расчета является ts = Δt / (N – 1). fmax
значения по умолчанию к 1 / (2 × ts), который для однородно произведенных сигналов соответствует частоте Найквиста.
[
задает целочисленный фактор сверхдискретизации, pxx
,f
]
= plomb(___,fmax
,ofac
)ofac
. Использование ofac
интерполировать или сглаживать спектр напоминают дополняющий нуль метод для основанных на БПФ методов. pxx
снова возвращен в round
(fmax
/fmin) точки частоты, но минимальная частота, рассмотренная в этом случае, 1 / (ofac
× N × ts). ofac
значения по умолчанию к 4.
[___] = plomb(___,
задает нормализацию периодограммы.spectrumtype
)
Установите spectrumtype
к 'psd'
, или оставьте его незаданным, чтобы получить pxx
как степень спектральная плотность.
Установите spectrumtype
к 'power'
получить спектр мощности входного сигнала.
Установите spectrumtype
к 'normalized'
получить стандартную периодограмму Lomb-Scargle, которая масштабируется два раза отклонением x
.
[___,
возвращает порог уровня мощности, pth
] = plomb(___,'Pd',pdvec
)pth
, таким образом, что пик со значением, больше, чем pth
имеет вероятность pdvec
из того, чтобы быть истинным пиком сигнала а не результатом случайных колебаний. pdvec
может быть вектор. Каждый элемент pdvec
должен быть больше 0 и меньше, чем 1. Каждая строка pth
соответствует элементу pdvec
. pth
имеет то же количество каналов как x
. Эта опция не доступна, если вы задаете выходные частоты в fvec
.
[
возвращает оценку PSD pxx
,w
]
= plomb(x
)x
оцененный в наборе равномерно расположенных с интервалами нормированных частот, w
, охват интервала Найквиста. Используйте NaN
s, чтобы задать недостающие выборки. Все вышеупомянутые опции доступны для нормированных частот. Чтобы получить доступ к ним, задайте пустой массив как второй вход.
plomb(___)
без выходных аргументов строит периодограмму Lomb-Scargle оценка PSD в окне текущей фигуры.
[1] Хорн, Джеймс Х. и Салли Л. Бэлиунас. "Предписание для Анализа Периода Неравномерно Выбранных Временных рядов". Астрофизический Журнал. Издание 302, 1986, стр 757–763.
[2] Lomb, Николас Р. "Анализ Частоты наименьших квадратов Неравноценно Расположенных с интервалами Данных". Астрофизика и Космические исследования. Издание 39, 1976, стр 447–462.
[3] Нажмите, Уильям Х. и Джордж Б. Рибики. "Алгоритм FAST для Спектрального анализа Неравномерно Выборочных данных". Астрофизический Журнал. Издание 338, 1989, стр 277–280.
[4] Scargle, Джеффри Д. "Исследования в Астрономическом Анализе Временных рядов. II. Статистические Аспекты Спектрального анализа Неравномерно Расположенных с интервалами Данных". Астрофизический Журнал. Издание 263, 1982, стр 835–853.
bandpower
| pburg
| pcov
| peig
| periodogram
| pmcov
| pmtm
| pmusic
| pwelch
| pyulear
| spectrogram