Exponenta Banner
exponenta event banner

cwtfreqbounds

Максимальная и минимальная частота или период CWT

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

Синтаксис

[minfreq, maxfreq] = cwtfreqbounds (N)
[minfreq, maxfreq] = cwtfreqbounds (N, Fs)
[maxperiod, minperiod] = cwtfreqbounds (N, Ts)
[___] = cwtfreqbounds (___, имя, значение)

Описание

пример

[minfreq,maxfreq] = cwtfreqbounds(N) возвращает минимальную и максимальную частоты полосы пропускания вейвлетов в циклах/выборке для сигнала длиной N. Минимальная и максимальная частоты определяются для импульса Морса по умолчанию (3,60). Минимальная частота определяется так, чтобы два стандартных временных отклонения вейвлета по умолчанию охватывали N- точечный сигнал на самой грубой шкале. Максимальная частота такова, что самый высокий частотный полосовой фильтр падает до ½ его пиковой величины на частоте Найквиста.

[minfreq,maxfreq] = cwtfreqbounds(N,Fs) возвращает полосовые частоты в герцах для частоты дискретизации Fs.

[maxperiod,minperiod] = cwtfreqbounds(N,Ts) возвращает периоды полосы пропускания для периода выборки Ts. maxperiod и minperiod скалярные длительности с тем же форматом, что и Ts. Если количество стандартных отклонений установлено таким образом, что log2(maxperiod/minperiod) < 1/NV где NV - количество голосов на октаву, maxperiod настраивается на minperiod*2^(1/NV).

пример

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

Примеры

свернуть все

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

Получение минимальной и максимальной частот полосы пропускания для сигнала с 1000 отсчетами с использованием значений по умолчанию.

[minfreq,maxfreq] = cwtfreqbounds(1000)
minfreq = 0.0033

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

Получают минимальную и максимальную частоты полосы пропускания импульса для импульса Морса по умолчанию для сигнала длиной 10000 и частоты дискретизации 1kHz. Установите отсечку на 100% так, чтобы наибольший частотный полосовой фильтр в Nyquist достиг пика.

sigLength = 10000;
Fs = 1e3;
[minfreq,maxfreq] = cwtfreqbounds(sigLength,Fs,'cutoff',100);

Создание банка фильтров с использованием значений, возвращаемых cwtfreqbounds. Постройте график частотной характеристики. Следует отметить, что наибольший частотный полосовой фильтр имеет пики на частоте Найквиста 500 Гц.

fb = cwtfilterbank('SignalLength',sigLength,'SamplingFrequency',Fs,'FrequencyLimits',[minfreq maxfreq]);
freqz(fb)

Figure contains an axes. The axes with title CWT Filter Bank contains 106 objects of type line.

Создайте второй банк частот, идентичный первому, но используйте ограничения частоты по умолчанию. Постройте график частотной характеристики. В случае вейвлета Морса банк фильтров CWT использует пороговое значение по умолчанию 50% в Nyquist.

fb2 = cwtfilterbank('SignalLength',sigLength,'SamplingFrequency',Fs);
figure
freqz(fb2)

Figure contains an axes. The axes with title CWT Filter Bank contains 104 objects of type line.

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

свернуть все

Длина сигнала, заданная как положительное целое число, большее или равное 4.

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

Типы данных: double

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

Типы данных: duration

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

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

Пример: 'Wavelet',"amor",'CUTOFF',75

Вейвлет анализа, используемый для определения минимальной и максимальной частот или периодов, указанных как 'Morse', 'amor', или 'bump'. Эти строки определяют аналитический вейвлет Морса, Морлета и ударного импульса соответственно. Вейвлет по умолчанию - это аналитический вейвлет Морса (3,60).

Для вейвлетов Морса можно также параметризовать вейвлет с помощью TimeBandwidth или WaveletParameters пары имя-значение.

Пример: 'Wavelet','bump'

Процент пиковой величины при Nyquist, указанный как скаляр между 0 и 100. Настройка CUTOFF на 0 указывает, что частотно-вейвлет-отклик затухает до 0 в Nyquist. Настройка CUTOFF до 100 указывает на то, что значение вейвлет-полосовых фильтров достигает пика в Nyquist.

Типы данных: double

Количество стандартных временных отклонений, используемых для определения минимальной частоты (самая длинная шкала), определяемое как положительное целое число, большее или равное 2. Для морсовского, аналитического морлетовского и ударного вейвлетов четыре стандартных отклонения обычно гарантируют, что вейвлет затухает до нуля на концах поддержки сигнала. Приращение 'StandardDeviations'на кратные 4, например 4 * M, гарантирует, что M целых вейвлетов умещаются в пределах длины сигнала. Если не указано, 'StandardDeviations'по умолчанию - 2. Если количество стандартных отклонений установлено так, что log2 (minfreq/maxfreq) > -1/NV где NV - количество голосов на октаву, minfreq настраивается на maxfreq* 2 ^ (-1/NV).

Типы данных: double

Это свойство доступно только для чтения.

Временная полоса пропускания для вейвлетов Морса, заданная как положительный скаляр. Это свойство допустимо только в том случае, если Wavelet свойство - 'morse'. Это свойство определяет параметр временной полосы волны Морса с параметром симметрии (гамма), зафиксированным на 3. TimeBandwidth является положительным числом, строго большим 3 и меньшим или равным 120. Чем больше параметр временной полосы пропускания, тем более распределенным является вейвлет во времени и более узким является вейвлет по частоте. Стандартное отклонение импульса Морса во времени составляет приблизительно sqrt(TimeBandwidth/2). Стандартное отклонение частоты составляет приблизительно 1/2*sqrt(2/TimeBandwidth).

TimeBandwidth и WaveletParameters не могут быть указаны оба свойства.

Пример: 'TimeBandwidth',4

Типы данных: double

Это свойство доступно только для чтения.

Параметры импульса Морзе, заданные как двухэлементный вектор. Первым элементом является параметр симметрии (гамма), который должен быть больше или равен 1. Вторым элементом является параметр временной полосы пропускания, который должен быть строго больше гамма. Отношение параметра временной полосы пропускания к гамма не может превышать 40. Когда гамма равна 3, импульс Морса идеально симметричен в частотной области. Перекос равен 0. Значения гамма более 3 приводят к положительному перекосу, а значения гамма менее 3 - к отрицательному перекосу. WaveletParameters действителен только в том случае, если Wavelet свойство имеет значение 'Morse'.

WaveletParameters и TimeBandwidth не могут быть указаны оба свойства.

Пример: 'WaveletParameters',[4,20]

Число голосов на октаву для использования при определении необходимого разделения между минимальной и максимальной шкалами, указанное как чётное целое число между 4 и 48. Минимальная и максимальная шкалы эквивалентны минимальной и максимальной частотам или максимальным и минимальным периодам соответственно.

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

Типы данных: double

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

свернуть все

Минимальная частота полосы пропускания вейвлетов, возвращаемая в виде скаляра. minfreq находится в циклах/выборке, если SamplingFrequency не указан. В противном случае minfreq в герцах.

Максимальная частота полосы пропускания вейвлетов, возвращаемая как скаляр. maxfreq находится в циклах/выборке, если SamplingFrequency не указан. В противном случае maxfreq в герцах.

Максимальный период полосы пропускания вейвлетов, возвращаемый как скалярная длительность с тем же форматом, что и Ts.

Если количество стандартных отклонений установлено таким образом, что log2(maxperiod/minperiod) < 1/NV где NV - количество голосов на октаву, maxperiod настраивается на minperiod*2^(1/NV).

Минимальный период полосы пропускания вейвлета, возвращаемый как скалярная длительность с тем же форматом, что и Ts.

Если количество стандартных отклонений установлено таким образом, что log2(maxperiod/minperiod) < 1/NV где NV - количество голосов на октаву, maxperiod настраивается на minperiod*2^(1/NV)

Расширенные возможности

.

См. также

|

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

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

Поддержка