Вейвлет synchrosqueezed преобразовывает
sst = wsst(x)sst, который вы используете, чтобы исследовать данные в плоскости частоты времени. synchrosqueezed преобразовывают, уменьшал энергетическое смазывание, когда по сравнению с непрерывным вейвлетом преобразовывают. Вход, x, должен быть 1D сигнал с действительным знаком по крайней мере с четырьмя выборками. wsst вычисляет synchrosqueezed, преобразовывают использование аналитического вейвлета Morlet.
[___] = wsst( использование a x,ts)duration ts с положительным, скалярным входом, как интервал выборки. Длительность может быть в годах, днях, часах, минутах или секундах. Если вы задаете ts и f вывод , wsst возвращает частоты в f в циклах в единицу времени, где единица измерения времени выведена из заданной длительности.
[___] = wsst(___, использует аналитический вейвлет, заданный wav)wav чтобы вычислить synchrosqueezed преобразовывают. Допустимыми значениями является 'amor' и 'bump', которые задают аналитический Morlet и вейвлет удара, соответственно.
wsst(___) без графиков выходных аргументов synchrosqueezed преобразовывают как функция времени и частоты. Если вы не задаете частоту дискретизации, fs, или интервал, ts, synchrosqueezed преобразовывают, построен в циклах на выборку. Если вы задаете частоту дискретизации, synchrosqueezed преобразовывают, построен в Гц. Если вы задаете интервал выборки с помощью длительности, график находится в циклах в единицу времени. Единицы измерения времени выведены из длительности.
[___] = wsst(___, возвращается synchrosqueezed преобразовывают с дополнительными опциями, заданными одним или несколькими Name,Value)Name,Value парные аргументы.
Получите вейвлет synchrosqueezed преобразование речевой выборки использование значений по умолчанию.
load mtlb;
sst = wsst(mtlb);Получите вейвлет synchrosqueezed, преобразование речи сигнализируют и сравнивают исходные и восстановленные сигналы.
Загрузите речь, сигнализируют и получают ее synchrosqueezed, преобразовывают.
load mtlb
soundsc(mtlb,Fs)
dt = 1/Fs;
t = 0:dt:numel(mtlb)*dt-dt;
[sst,f] = wsst(mtlb,Fs);График synchrosqueezed преобразовывает.
pcolor(t,f,abs(sst)) shading interp xlabel('Seconds') ylabel('Frequency (Hz)') title('Synchrosqueezed Transform')
Получите инверсию synchrosqueezed, преобразовывают и проигрывают восстановленный речевой сигнал.
xrec = iwsst(sst); soundsc(xrec,Fs)
Получите и постройте вейвлет synchrosqueezed преобразование квадратичного щебета. Щебет производится на уровне 1 000 Гц.
load quadchirp; [sst,f] = wsst(quadchirp,1000); hp = pcolor(tquad,f,abs(sst)); hp.EdgeColor = 'none'; title('Wavelet Synchrosqueezed Transform'); xlabel('Time'); ylabel('Hz');
Получите вейвлет synchrosqueezed преобразование данных о солнечном пятне с помощью вейвлета Morlet по умолчанию. Задайте интервал выборки, чтобы быть одним годом.
load sunspot
wsst(sunspot(:,2),years(1))
Получите и постройте вейвлет synchrosqueezed преобразование данных о солнечном пятне с помощью вейвлета удара. Задайте интервал выборки, чтобы быть 1 для одной выборки в год.
load sunspot wsst(sunspot(:,2),years(1),'bump')
[1] Daubechies, я., Дж. Лу и H.-T. Ву. "Вейвлет Synchrosqueezed преобразовывает: подобный эмпирическому разложению моды инструмент". Примененный и Вычислительный Гармонический Анализ. Издание 30, Номер 2, 2011, стр 243–261.
[2] Thakur, G., Э. Бревдо, N. S. Fučkar и H.-T. Ву. "Алгоритм Synchrosqueezing для изменяющегося во времени спектрального анализа: свойства робастности и новые приложения палеоклимата". Обработка сигналов. Издание 93, Номер 5, 2013, стр 1079–1094.