Степень мультизаострения спектральная оценка плотности
pxx = pmtm(x)
pxx = pmtm(x,nw)
pxx = pmtm(x,nw,nfft)
[pxx,w] = pmtm(___)
[pxx,f] = pmtm(___,fs)
[pxx,w] = pmtm(x,nw,w)
[pxx,f] = pmtm(x,nw,f,fs)
[___] = pmtm(___,method)
[___] = pmtm(x,e,v)
[___] = pmtm(x,dpss_params)
[___] = pmtm(___,'DropLastTaper',dropflag)
[___] = pmtm(___,freqrange)
[___,pxxc] = pmtm(___,'ConfidenceLevel',probability)
pmtm(___)
возвращает оценку степени спектральной плотности (PSD) мультизаострения Томсона, pxx
= pmtm(x
)pxx
, входного сигнала, x
. Когда x
является вектором, он обработан как один канал. Когда x
является матрицей, PSD вычисляется независимо для каждого столбца и хранится в соответствующем столбце pxx
. Заострения являются дискретным вытянутым сфероидальным (DPSS), или Slepian, последовательностями. Полупропускная способность времени, nw
, продукт равняется 4. По умолчанию pmtm
использует первые 2 × nw
– 1 последовательность DPSS. Если x
с действительным знаком, pxx
является односторонней оценкой PSD. Если x
с комплексным знаком, pxx
является двухсторонней оценкой PSD. Число точек, nfft
, в дискретном преобразовании Фурье (DFT) являются максимумом 256 или следующая степень двойки, больше, чем длина сигнала.
[
возвращает вектор частоты, pxx
,f
] = pmtm(___,fs
)f
, в циклах в единицу времени. Частота дискретизации, fs
, является количеством выборок в единицу времени. Если модуль времени является секундами, то f
находится в циклах/секунда (Гц). Для сигналов с действительным знаком f
охватывает интервал [0, fs
/2], когда nfft
даже и [0, fs
/2), когда nfft
нечетен. Для сигналов с комплексным знаком f
охватывает интервал [0, fs
). fs
должен быть четвертым входом к pmtm
. Чтобы ввести частоту дискретизации и все еще использовать значения по умолчанию предыдущих дополнительных аргументов, задайте эти аргументы как пустые, []
.
[
возвращает двухстороннее мультизаострение оценки PSD на частотах, заданных в векторе, pxx
,f
] = pmtm(x
,nw
,f
,fs
)f
. Векторный f
должен содержать по крайней мере два элемента, потому что в противном случае функция интерпретирует его как nfft
. Частоты в f
находятся в циклах в единицу времени. Частота дискретизации, fs
, является количеством выборок в единицу времени. Если модуль времени является секундами, то f
находится в циклах/секунда (Гц).
[___] = pmtm(
использует массив ячеек, x
,dpss_params
)dpss_params
, чтобы передать входные параметры dpss
кроме числа элементов в последовательностях. Число элементов в последовательностях является первым входным параметром к dpss
и не включено в dpss_params
. Примером этого использования является pxx = pmtm(randn(1000,1),{2.5,3})
.
[___] = pmtm(___,'DropLastTaper',
задает, роняет ли dropflag
)pmtm
последнее заострение в вычислении мультизаострения оценка PSD. dropflag
является логическим. Значением по умолчанию dropflag
является true
, и последнее заострение не используется в оценке PSD.
[___] = pmtm(___,
возвращает мультизаострение оценка PSD по частотному диапазону, заданному freqrange
)freqrange
. Допустимыми опциями для freqrange
является 'onesided'
, 'twosided'
и 'centered'
.
[___,
возвращает pxxc
] = pmtm(___,'ConfidenceLevel',probability
)probability
× 100% доверительных интервалов для оценки PSD в pxxc
.
pmtm(___)
без выходных аргументов строит мультизаострение оценка PSD в окне текущей фигуры.
[1] Персиваль, D. B., и А. Т. Уолден, спектральный анализ для физических приложений: мультизаостритесь и обычные одномерные методы. Кембридж, Великобритания: Издательство Кембриджского университета, 1993.
[2] Томсон, D. J. “Оценка спектра и гармонический анализ”. Продолжения IEEE®. Издание 70, 1982, стр 1055–1096.
dpss
| periodogram
| pwelch