Двойные древовидные биоортогональные фильтры первого уровня
Получите фильтры разложения и реконструкции, сопоставленные с биоортогональным вейвлетом nearsym5_7.
fname = 'nearsym5_7';
[LoD,HiD,LoR,HiR] = qbiorthfilt(fname);Используйте dwtfilterbank функция, чтобы создать 7-уровневый дискретный банк фильтра преобразований вейвлета с биоортогональными фильтрами. Задайте тип фильтра вейвлета как analysis. Поскольку фильтры не имеют даже длин, расширяют фильтры соответственно, чтобы совпадать со степенями их z-преобразований.
scal(:,1) = [0 0 LoD' 0]; scal(:,2) = [0 LoR']; wavf(:,1) = [0 HiD']; wavf(:,2) = [0 0 HiR' 0]; fb = dwtfilterbank('Wavelet','Custom',... 'CustomScalingFilter',scal,... 'CustomWaveletFilter',wavf,... 'Level',7,... 'FilterType','analysis');
Получите вейвлеты временного интервала, соответствующие фильтрам полосы пропускания вейвлета. Постройте вейвлет самой грубой шкалы.
[psi,t] = wavelets(fb); plot(t,psi(end,:)) grid on xlabel('Time') ylabel('Amplitude')
[1] Antonini, M., М. Барло, П. Мэтью и я. Daubechies. “Отобразите Кодирование Используя Преобразование Вейвлета”. Транзакции IEEE на Обработке изображений 1, № 2 (апрель 1992): 205–20. https://doi.org/10.1109/83.136597.
[2] Кингсбери, Ник. “Комплексные Вейвлеты для Анализа Инварианта Сдвига и Фильтрации Сигналов”. Примененный и Вычислительный Гармонический Анализ 10, № 3 (май 2001): 234–53. https://doi.org/10.1006/acha.2000.0343.
[3] Le Gall, D. и А. Тэбэйтабай. “Кодирование поддиапазона Цифровых изображений Используя Симметричные Короткие Фильтры Ядра и Методы Кодирования Арифметики”. В ICASSP-88., Международная конференция по вопросам Акустики, Речи и Обработки сигналов, 761–64. Нью-Йорк, Нью-Йорк, США: IEEE, 1988. https://doi.org/10.1109/ICASSP.1988.196696.