Этот пример показывает, как определить фильтр масштабирования экстремальной фазы Daubechies с заданной суммой. Двумя наиболее распространенными значениями для суммы являются $$\sqrt{2}$$ и 1.

Создайте две версии db4 масштабирующий фильтр. Один масштабирующий фильтр суммируется с $$\sqrt{2}$$ а другая версия равна 1.

NumVanishingMoments = 4;
h = dbaux(NumVanishingMoments,sqrt(2));
m0 = dbaux(NumVanishingMoments,1);

Фильтр с суммой, равной $$\sqrt{2}$$ - фильтр синтеза (реконструкции), возвращаемый wfilters и используется в дискретном вейвлет.

[LoD,HiD,LoR,HiR] = wfilters('db4');
max(abs(LoR-h))

ans = 4.2590e-13

Для ортогональных вейвлеты фильтр анализа (разложения) является обратным по времени фильтром синтеза.

max(abs(LoD-fliplr(h)))

ans = 4.2590e-13

Этот пример показывает, что симлетные и Daubechies масштабирующие фильтры одного и того же порядка являются решениями одного и того же полиномиального уравнения.

Сгенерируйте порядок 4 масштабирующий фильтр Daubechies и постройте его график.

wdb4 = dbaux(4)

wdb4 = 1×8

    0.1629    0.5055    0.4461   -0.0198   -0.1323    0.0218    0.0233   -0.0075

stem(wdb4)
title('Order 4 Daubechies Scaling Filter')

wdb4 является решением уравнения: P = conv (wrev (w), w) * 2, где P является фильтром «Lagrange trous» для N = 4. Оцените P и постройте график. P является симметричным фильтром, и wdb4 является минимальным фазовым решением предыдущего уравнения на основе корней P.

P = conv(wrev(wdb4),wdb4)*2;
stem(P)
title('''Lagrange trous'' filter')

Сгенерируйте wsym4, фильтр симлетного масштабирования порядка 4 и постройте его график. Симлеты являются «наименее асимметричными» вейвлетами Доубехи, полученными из другого выбора между корнями P.

wsym4 = symaux(4)

wsym4 = 1×8

    0.0228   -0.0089   -0.0702    0.2106    0.5683    0.3519   -0.0210   -0.0536

stem(wsym4)
title('Order 4 Symlet Scaling Filter')

Вычислите conv (wrev (wsym4), wsym4) * 2 и подтвердите, что wsym4 является другим решением уравнения P = conv (wrev (w), w) * 2.

P_sym = conv(wrev(wsym4),wsym4)*2;
err = norm(P_sym-P)

err = 1.8677e-15

Этот пример демонстрирует, что для заданной поддержки совокупная сумма квадратов коэффициентов масштабирующего фильтра увеличивается более быстро для экстремального вейвлета фазы, чем другие вейвлеты.

Сгенерируйте коэффициенты масштабирующего фильтра для db15 и sym15 вейвлеты. Оба вейвлета имеют поддержку ширины $2\times 15-1=29$.

[~,~,LoR_db,~] = wfilters('db15');
[~,~,LoR_sym,~] = wfilters('sym15');

Затем сгенерируйте коэффициенты масштабирующего фильтра для coif5 вейвлет. Этот вейвлет также имеет поддержку ширины $6\times 5-1=29$.

[~,~,LoR_coif,~] = wfilters('coif5');

Подтвердите сумму коэффициентов для всех трех вейвлет равной $\sqrt{2}$.

sqrt(2)-sum(LoR_db)

ans = 2.2204e-16

sqrt(2)-sum(LoR_sym)

ans = 0

sqrt(2)-sum(LoR_coif)

ans = 2.2204e-16

Постройте график совокупных сумм квадратов коэффициентов. Обратите внимание, как быстро сумма Daubechies увеличивается. Это связано с тем, что его энергия сосредоточена в небольших абсциссах. Поскольку вейвлет Daubechies имеет экстремальную фазу, совокупная сумма его квадратов коэффициентов увеличивается быстрее, чем два других вейвлета.

plot(cumsum(LoR_db.^2),'rx-')
hold on
plot(cumsum(LoR_sym.^2),'mo-')
plot(cumsum(LoR_coif.^2),'b*-')
legend('Daubechies','Symlet','Coiflet')
title('Cumulative Sum')