exponenta event banner

небольшие волны

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

Описание

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

пример

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

Примеры

свернуть все

Создание семиуровневого набора фильтров DWT с длиной сигнала 1000 выборок с помощью Daubechies db2 вейвлет и частота дискретизации 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-на-N, где L - набор фильтров Level и N представляет собой SignalLength. Вейвлеты упорядочены в psi от наилучшего разрешения шкалы до самого крупного разрешения шкалы.

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

Представлен в R2018a