Полосы пропускания набора фильтров CWT 3 дБ
bw = powerbw(fb)fb. bw Ns-by-4 MATLAB®
table, где Ns является количеством частот полосы пропускания вейвлета (равный количеству шкал). Для каждого фильтра в fb, таблица содержит соответствующую полосовую частоту, полосу пропускания на 3 дБ, и более низкую частоту и верхние пределы частоты полосы пропускания на 3 дБ.
Полоса пропускания на 3 дБ ограничивает метку, где степень фильтра является половиной ее пикового значения. Частотная характеристика величины в пределах равна 1 / √ 2 раза пиковая величина. Начиная с полос пропускания в fb нормированы с пиковыми величинами, приблизительно равняются 2, частотная характеристика величины в каждом пределе приблизительно равна 2 / √ 2. Полоса пропускания на 3 дБ также известна как полосу пропускания на уровне половинной мощности потому что
Создайте набор фильтров CWT.
fb = cwtfilterbank;
Получите полосы пропускания (на уровне половинной мощности) на 3 дБ набора фильтров. Получите частотные характеристики вейвлетов.
bw = powerbw(fb); [psidft,f] = freqz(fb);
Выберите полосовой фильтр вейвлета из набора фильтров. Извлеките из таблицы bw пределы на 3 дБ полосового фильтра.
wv = 5; frq = bw.Frequencies(wv); lfb = bw.LowFrequencyBorder(wv); hfb = bw.HighFrequencyBorder(wv);
Постройте частотную характеристику и пределы на 3 дБ. Поскольку частотная характеристика масштабируется, чтобы иметь максимальное значение, равное 2, смотреть график подтвердить, что более низкие и верхние границы частоты пересекают частотную характеристику в .
plot(f,psidft(wv,:)) grid on hold on plot([lfb lfb],[0 2],'r') plot([hfb hfb],[0 2],'r') xlabel('Normalized Frequency (cycles/sample)') ylabel('Magnitude') title(['Bandpass Frequency: ' num2str(frq) ' cycles/sample'])
