Exponenta Banner
exponenta event banner

wsst

Вейвлет, синхронизированное преобразование

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

Синтаксис

sst = wsst (x)
[sst, f] = wsst (x)
[___] = wsst (x, fs)
[___] = wsst (x, ts)
[___] = wsst (___, wav)
wsst (___)
[___] = wsst (___, имя, значение)

Описание

пример

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

пример

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

пример

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

пример

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

пример

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

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

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

Примеры

свернуть все

Открыть сценарий в реальном времени

Получение вейвлет-синхронного преобразования образца речи с использованием значений по умолчанию.

load mtlb;
sst = wsst(mtlb);
Открыть сценарий в реальном времени

Получить вейвлет-синхронизированное преобразование речевого сигнала и сравнить исходные и восстановленные сигналы.

Загрузите речевой сигнал и получите его синхронизированное преобразование.

load mtlb
soundsc(mtlb,Fs)
dt = 1/Fs;
t = 0:dt:numel(mtlb)*dt-dt;
[sst,f] = wsst(mtlb,Fs);

Постройте график синхронизированного преобразования.

pcolor(t,f,abs(sst))
shading interp
xlabel('Seconds')
ylabel('Frequency (Hz)')
title('Synchrosqueezed Transform')

Получить обратное синхронизированное преобразование и воспроизвести восстановленный речевой сигнал.

xrec = iwsst(sst);
soundsc(xrec,Fs)
Открыть сценарий в реальном времени

Получить и построить график импульсного синхроскопического преобразования квадратичной чирп. Чирп отбирают при частоте 1000 Гц.

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. The axes with title Wavelet Synchrosqueezed Transform contains an object of type surface.

Открыть сценарий в реальном времени

Получить импульсное синхронное преобразование данных солнечных точек с использованием используемого по умолчанию вейвлета Морле. Укажите интервал выборки, равный одному году.

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

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

Открыть сценарий в реальном времени

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

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

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

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

свернуть все

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

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

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

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

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

  • 'amor' - Аналитический вейвлет Морле

  • 'bump' - Ударный импульс

Аргументы пары «имя-значение»

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

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

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

Опция для симметричного расширения входного сигнала, определяемого как разделенная запятыми пара, состоящая из 'ExtendSignal' и либо false или true. Симметричное расширение сигнала может ослабить граничные эффекты. При указании false, то сигнал не выдвигается. При указании true, то сигнал удлиняется.

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

свернуть все

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

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

Ссылки

[1] Daubechies, I., J. Lu и H.-T. У. «Синхронизированные вейвлет-преобразования: эмпирический инструмент, похожий на декомпозицию». Прикладной и вычислительный гармонический анализ. Том 30, номер 2, 2011, стр. 243-261.

[2] Такур, Г., Э. Бревдо, Н. С. Фучкар и Х.-Т. У. «Алгоритм Synchrosqueezing для изменяющегося во времени спектрального анализа: свойства надежности и новые приложения палеоклимата». Обработка сигналов. Том 93, номер 5, 2013, стр. 1079-1094.

См. также

| | | | | | |

Темы

Представлен в R2016a

Документация инструментария вейвлет

Поддержка