Визуализация вейвлетов, вейвлет-пакетов и вейвлет-фильтров
В этом примере показано, как использовать wfilters, wavefun, и wpfun для получения фильтров, вейвлет или вейвлет пакетов, соответствующих конкретному семейству вейвлет. Вы можете визуализировать 2-D разделяемые вейвлеты с wavefun2.
Получите фильтры разложения (анализа) и реконструкции (синтеза) для биортогональных сплайн вейвлета фильтров с 3 моментами исчезновения в фильтре реконструкции и 5 моментами исчезновения в фильтре разложения.
[LoD,HiD,LoR,HiR] = wfilters('bior3.5'); subplot(2,2,1) stem(LoD,'markerfacecolor',[0 0 1]); title('Lowpass Decomposition Filter'); subplot(2,2,2) stem(LoR,'markerfacecolor',[0 0 1]); title('Lowpass Reconstruction Filter'); subplot(2,2,3) stem(HiD,'markerfacecolor',[0 0 1]); title('Highpass Decomposition Filter'); subplot(2,2,4) stem(HiR,'markerfacecolor',[0 0 1]); title('Highpass Reconstruction Filter');
Визуализируйте реальный вейвлет Морле. Нет связанной функции масштабирования.
figure [psi,xval] = wavefun('morl'); plot(xval,psi,'linewidth',2) title('$\psi(x) = e^{-x^2/2} \cos{(5x)}$','Interpreter','latex',... 'fontsize',14);
Получите первые 4 вейвлет-пакета для наименее асимметричного вейвлета Daubechies с 4 моментами исчезновения sym4.
[wpws,x] = wpfun('sym4',4,10); for nn = 1:size(wpws,1) subplot(3,2,nn) plot(x,wpws(nn,:)) axis tight title(['W',num2str(nn-1)]); end
