Exponenta Banner

wsst

Вейвлет synchrosqueezed преобразовывает

свернуть все на странице

Синтаксис

sst = wsst(x)
[sst,f] = wsst(x)
[___] = wsst(x,fs)
[___] = wsst(x,ts)
[___] = wsst(___,wav)
wsst(___)
[___] = wsst(___,Name,Value)

Описание

пример

sst = wsst(x) возвращается вейвлет synchrosqueezed преобразовывают, sst, который вы используете, чтобы исследовать данные в плоскости частоты времени. synchrosqueezed преобразовывают, уменьшал энергетическое смазывание, когда по сравнению с непрерывным вейвлетом преобразовывают. Вход, x, должен быть 1D вещественный сигнал по крайней мере с четырьмя выборками. wsst вычисляет synchrosqueezed, преобразовывают использование аналитического вейвлета Morlet.

пример

[sst,f] = wsst(x) возвращает вектор из частот, f, в циклах на выборку. Частоты соответствуют строкам sst.

пример

[___] = wsst(x,fs) вычисляет synchrosqueezed, преобразовывают использование заданной частоты дискретизации, fs, в Гц, чтобы вычислить synchrosqueezed преобразовывают. Если вы задаете f вывод , wsst возвращает частоты в Гц. Можно использовать любую предыдущую комбинацию выходных значений.

пример

[___] = wsst(x,ts) использование a duration ts с положительным, скалярным входом, как интервал выборки. Длительность может быть в годах, днях, часах, минутах или секундах. Если вы задаете ts и f вывод , wsst возвращает частоты в f в циклах в единицу времени, где единица измерения времени выведена из заданной длительности.

пример

[___] = wsst(___,wav) использует аналитический вейвлет, заданный wav чтобы вычислить synchrosqueezed преобразовывают. Допустимыми значениями является 'amor' и 'bump', которые задают аналитический Morlet и вейвлет удара, соответственно.

wsst(___) без графиков выходных аргументов synchrosqueezed преобразовывают в зависимости от времени и частоты. Если вы не задаете частоту дискретизации, fs, или интервал, ts, synchrosqueezed преобразовывают, построен в циклах на выборку. Если вы задаете частоту дискретизации, synchrosqueezed преобразовывают, построен в Гц. Если вы задаете интервал выборки с помощью длительности, график находится в циклах в единицу времени. Единицы измерения времени выведены из длительности.

[___] = wsst(___,Name,Value) возвращается synchrosqueezed преобразовывают с дополнительными опциями, заданными одним или несколькими Name,Value парные аргументы.

Примеры

свернуть все

Скрипт Open Live Script

Получите вейвлет synchrosqueezed преобразование речевой выборки использование значений по умолчанию.

load mtlb;
sst = wsst(mtlb);
Скрипт Open Live Script

Получите вейвлет 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)
Скрипт Open Live Script

Получите и постройте вейвлет 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');

Figure contains an axes object. The axes object with title Wavelet Synchrosqueezed Transform contains an object of type surface.

Скрипт Open Live Script

Получите вейвлет synchrosqueezed преобразование данных о солнечном пятне с помощью вейвлета Morlet по умолчанию. Задайте интервал выборки, чтобы быть одним годом.

load sunspot
wsst(sunspot(:,2),years(1))

Figure contains an axes object. The axes object with title Wavelet Synchrosqueezed Transform contains an object of type surface.

Скрипт Open Live Script

Получите и постройте вейвлет synchrosqueezed преобразование данных о солнечном пятне с помощью вейвлета удара. Задайте интервал выборки, чтобы быть 1 для одной выборки в год.

load sunspot
wsst(sunspot(:,2),years(1),'bump')

Figure contains an axes object. The axes object with title Wavelet Synchrosqueezed Transform contains an object of type surface.

Входные параметры

свернуть все

Входной сигнал в виде строки или вектор-столбца. x должно быть 1D, вещественный сигнал по крайней мере с четырьмя выборками.

Частота дискретизации в виде положительной скалярной величины.

Интервал выборки, также известный как период выборки в виде a duration с входом положительной скалярной величины. Допустимая длительность years, days, hours, seconds, и minutes. Вы не можете использовать календарную длительность (caldays, calweeks, calmonths, calquarters, или calyears). Вы не можете задать оба ts и fs.

Пример: sst = wsst(x,hours(12))

Аналитический вейвлет использовался для расчета synchrosqueezed, преобразовывают в виде одного из следующего:

  • 'amor' — Аналитический вейвлет Morlet

  • 'bump' — Ударьте вейвлет

Аргументы name-value

Задайте дополнительные разделенные запятой пары Name,Value аргументы. Name имя аргумента и Value соответствующее значение. Name должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.

Пример: 'VoicesPerOctave',26

Количество речи на октаву, чтобы использовать в synchrosqueezed преобразовывает в виде разделенной запятой пары, состоящей из 'VoicesPerOctave' и ровное целое число от 10 до 48. Продуктом количества речи на октаву и количества октав является количество шкал. Количество октав зависит от размера входа x и floor(log2(numel(x)))-1.

Опция, чтобы расширить входной сигнал симметрично в виде разделенной запятой пары, состоящей из 'ExtendSignal' и любой false или true. Расширение сигнала симметрично может смягчить граничные последствия. Если вы задаете false, затем сигнал не расширен. Если вы задаете true, затем сигнал расширен.

Выходные аргументы

свернуть все

Synchrosqueezed преобразовывают, возвращенный как матрица. По умолчанию synchrosqueezed преобразовывают, использует floor(log2(numel(x)))-1 октавы, 32 речи на октаву и аналитический вейвлет Morlet. sst Na-by-N матрица, где Na является количеством шкал, и N является количеством отсчетов в x. Количеством по умолчанию шкал является 32*(floor(log2(numel(x)))-1).

Частоты synchrosqueezed преобразовывают, возвращенный как вектор. Частоты соответствуют строкам sst. Если вы не задаете fs или ts, частоты находятся в циклах на выборку. Если вы задаете fs, частоты находятся в Гц. Если вы задаете ts, частоты находятся в циклах в единицу времени. Длина вектора частоты совпадает с количеством sst 'Строки' . Если вы задаете ts как интервал выборки, ts используется для расчета преобразование шкалы к частоте в f.

Ссылки

[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.

Смотрите также

| | | | | | |

Темы

Введенный в R2016a

Документация Wavelet Toolbox

Поддержка