Банк фильтров 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'])
