wavelets

DWT-фильтр, вейвлеты временной области банка

Описание

psi = wavelets(fb) возвращает частотные и центрированные вейвлеты, соответствующие вейвлет-полосовым фильтрам в дискретной группе фильтров вейвлет-преобразования (DWT) fb.

пример

[psi,t] = wavelets(fb) возвращает моменты дискретизации t.

Примеры

свернуть все

Создайте семиуровневую группу фильтров DWT с длиной сигнала 1000 выборок, используя db2 Daubechies вейвлет и частота дискретизации 1 кГц.

wv = "db4";
len = 1000;
lev = 7;
Fs = 1e3;
fb = dwtfilterbank('Wavelet',wv,'SignalLength',len,'Level',lev,'SamplingFrequency',Fs);

Постройте график частотных и центрированных вейвлетов, соответствующих вейвлет-полосно-пропускающим фильтрам.

[psi,t] = wavelets(fb);
plot(t,psi')
grid on
title('Time-domain Wavelets')

Figure contains an axes. The axes with title Time-domain Wavelets contains 7 objects of type line.

Постройте график наивысшей шкалы вейвлета временной области и односторонней величины частотной характеристики соответствующего вейвлет-полосно-пропускающего фильтра.

sc = 1;
[psidft,f] = freqz(fb);
subplot(2,1,1)
plot(t,psi(sc,:))
grid on
xlabel('Time (sec)')
ylabel('Magnitude')
title(['Level ',num2str(sc),' Time-Domain Wavelet'])
subplot(2,1,2)
plot(f(len/2:end),abs(psidft(sc,len/2:end)))
grid on
xlabel('Hz')
ylabel('Magnitude')
title('Magnitude Frequency Response')

Figure contains 2 axes. Axes 1 with title Level 1 Time-Domain Wavelet contains an object of type line. Axes 2 with title Magnitude Frequency Response contains an object of type line.

Входные параметры

свернуть все

Набор фильтров дискретного вейвлета преобразования (DWT), заданный как dwtfilterbank объект.

Выходные аргументы

свернуть все

Ориентированные по времени вейвлеты, соответствующие вейвлет-полосовым фильтрам, возвращаются как L -by - N матрица, где L - банк фильтров Level и N является SignalLength. Вейвлеты упорядочены в psi от самого точного разрешения шкалы до самого грубого разрешения шкалы.

Моменты дискретизации, возвращенные как действительный вектор t длины N, где N - банк фильтров SignalLength. Моменты дискретизации лежат в интервале [½NDT,½NDT), где DT - период дискретизации набора фильтров (обратный частоте дискретизации набора фильтров).

Введенный в R2018a