Томсоновский спектр мультизаострения
Hs = spectrum.mtm
Hs = spectrum.mtm(TimeBW)
Hs = spectrum.mtm(DPSS,Concentrations)
Hs = spectrum.mtm(...,CombineMethod)
Примечание
Использование spectrum.mtm не рекомендуется. Использование pmtm вместо этого.
Hs = spectrum.mtm возвращает Томсоновский объект спектра мультизаострения по умолчанию, Hs это задает параметры для Томсоновского мультизаострения спектральный алгоритм оценки, который использует линейную или нелинейную комбинацию модифицированных периодограмм. Периодограммы вычисляются с помощью последовательности ортогональных заострений (окна в частотном диапазоне) заданный от дискретных вытянутых сфероидальных последовательностей (dpss). Этот объект использует следующие значения по умолчанию:
| PropertyName | Значение по умолчанию | Описание |
|---|---|---|
| 4 | Продукт времени и полосы пропускания для дискретных вытянутых сфероидальных последовательностей (или последовательностей Slepian) используемый в качестве окон данных |
|
| Алгоритм для объединения отдельных спектральных оценок. Допустимыми значениями является |
Hs = spectrum.mtm(TimeBW) возвращает объект спектра, Hs с продуктом полосы пропускания требуемого времени.
Hs = spectrum.mtm(DPSS,Concentrations) возвращает объект спектра, Hs с заданными dpss данными заостряется и их концентрации.
Примечание
Можно или задать продукт полосы пропускания времени (TimeBW) или DPSS данные заостряются и их Concentrations. Смотрите dpss и pmtm для получения дополнительной информации.
Hs = spectrum.mtm(...,CombineMethod) возвращает объект спектра, Hs, с заданным методом для объединения спектральных оценок. Обратитесь к приведенной выше таблице для допустимого CombineMethod значения.
Задайте косинус 200 Гц, добавьте шум и просмотрите его спектральную плотность мощности с помощью Томсоновского алгоритма мультизаострения с продуктом полосы пропускания времени 3,5.
Fs=1000; t=0:1/Fs:.3; x=cos(2*pi*t*200)+randn(size(t)); Hs=spectrum.mtm(3.5); psd(Hs,x,'Fs',Fs)
Вышеупомянутый пример мог быть сделан путем определения заострений данных и концентраций вместо продукта полосы пропускания времени.
Fs=1000; t=0:1/Fs:.3; x=cos(2*pi*t*200)+randn(size(t)); [e,v]=dpss(length(x),3.5); Hs=spectrum.mtm(e,v); psd(Hs,x,'Fs',Fs)
periodogram | pmtm | pwelch