Визуализируя вейвлеты, пакеты вейвлета и фильтры вейвлета

Этот пример показывает, как использовать wfilters, wavefun и wpfun, чтобы получить фильтры, вейвлет или пакеты вейвлета, соответствующие конкретному семейству вейвлетов. Можно визуализировать 2D отделимые вейвлеты с 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');

Визуализируйте вейвлет Morlet с действительным знаком. Существует не сопоставленная масштабирующаяся функция.

figure
[psi,xval] = wavefun('morl');
plot(xval,psi,'linewidth',2)
title('$\psi(x) = e^{-x^2/2} \cos{(5x)}$','Interpreter','latex',...
     'fontsize',14);

Получите первые 4 пакета вейвлета для Добечи меньше всего - асимметричный вейвлет с 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

Смотрите также

|