Частотные характеристики банков фильтров DWT
[ возвращает комплексные частотные характеристики для вейвлет-фильтров psidft,f] = freqz(fb)psidft и частотный вектор f для банка фильтров дискретного вейвлет-преобразования (DWT) fb. Частоты находятся в циклах/выборке или в Гц, если частота выборки определена в fb. Частотные характеристики центрированы таким образом, что нулевая частота находится посередине.
freqz( строит график односторонних частотных откликов амплитуды для набора вейвлет-фильтров, fb)fb. Частотные характеристики амплитуды строятся для всех полосовых вейвлет-фильтров и фильтра масштабирования с наибольшим разрешением. Легенда интерактивна. Чтобы переключить видимость отклика величины фильтра, щелкните соответствующую линию в легенде.
Создайте набор фильтров DWT для сигнала длины 4096 и Fejér-Korovkin fk22 вейвлет. Постройте график частотно-амплитудных характеристик вейвлет-фильтров и масштабного фильтра с окончательным разрешением.
len = 4096; fb = dwtfilterbank('Wavelet','fk22','SignalLength',len); freqz(fb)
Получение частотных откликов для вейвлет-фильтров и фильтров масштабирования. Постройте график частотных характеристик масштабных фильтров на всех уровнях разложения.
[psidft,f,phidft] = freqz(fb); plot(f,abs(phidft)') grid on xlabel('Normalized Frequency (cycles/sample)') ylabel('Magnitude') legend('A1','A2','A3','A4','A5','A6','A7')
Постройте график односторонних частотных характеристик амплитуды вейвлет-фильтра и масштабирования на первых двух уровнях разложения. Обратите внимание, как частотные характеристики второго уровня перекрывают амплитудную характеристику фильтра масштабирования первого уровня.
plot(f(len/2:end),abs(psidft(1,len/2:end))') hold on plot(f(len/2:end),abs(phidft(1,len/2:end))') plot(f(len/2:end),abs(psidft(2,len/2:end))') plot(f(len/2:end),abs(phidft(2,len/2:end))') grid on xlabel('Normalized Frequency (cycles/sample)') ylabel('Magnitude') legend('Level 1 Wavelet','Level 1 Scaling','Level 2 Wavelet','Level 2 Scaling')
