Двухуровневые биортогональные фильтры первого уровня
Получите фильтры разложения и реконструкции, сопоставленные с биортогональным вейвлет 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., M. Barlaud, P. Mathieu, and I. 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., and A. Tabatabai. «Субдиапазонное кодирование цифровых изображений с использованием симметричных фильтров короткого ядра и методов арифметического кодирования». В ICASSP-88. Международная конференция по акустике, речи и обработке сигналов, 761-64. Нью-Йорк, Нью-Йорк, США: IEEE, 1988. https://doi.org/10.1109/ICASSP.1988.196696.