DWT фильтрация bank частотных характеристик
[ возвращает комплексные частотные характеристики для вейвлет psidft,f] = freqz(fb)psidft и вектор частоты f для дискретного вейвлета фильтров вейвлет-преобразования (DWT) fb. Частоты находятся в циклах/выборке или в Гц, если частота дискретизации задана как fb. Частотные характеристики центрированы так, чтобы нулевая частота была посередине.
freqz( строит графики односторонних величин частотных характеристик для вейвлета банка фильтров, fb)fb. Частотные характеристики величины строятся для всех вейвлет фильтров и фильтра масштабирования с наибольшим разрешением. Легенда интерактивна. Чтобы переключить видимость характеристики величины фильтра, щелкните соответствующую линию легенды.
Создайте банк фильтров DWT для сигнала длины 4096 и 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')
