Пропускная способность банка фильтров 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'])
