Вейвлеты фильтров рассеяния вейвлет-времени
filters = filterbank(sf)filters - массив ячеек из массивов структур с norder элементами, где norder - количество порядков рассеяния. Первый элемент filters содержит фильтр масштабирования, phift, используемый в расчете коэффициентов рассеяния 0-го порядка. Последующие элементы filters содержат вейвлет, psift, и масштабирующий фильтр, phift, для соответствующих блоков фильтров рассеивающего разложения.
Точность phift и psift зависит от точности, заданной в сети рассеяния sf.
[ возвращает параметры фильтра для каждого элемента filters,f,filterparams] = filterbank(sf)filters. filterparams - массив ячеек с norder элементами. Каждый элемент filterparams является MATLAB® таблица.
[___] = filterbank( возвращает банки фильтров, используемые для вычисления заданной sf,order)order коэффициенты рассеяния. order - целое число от 0 до nfilters включительно, где nfilters - количество банков фильтров в сети рассеяния. Эти входные параметры могут использоваться с любым из выходных синтаксисов, показанных ранее.
Создайте сеть вейвлет для сигнала, дискретизированного с частотой 25 Гц.
sf = waveletScattering('SamplingFrequency',25)sf =
waveletScattering with properties:
SignalLength: 1024
InvarianceScale: 20.4800
QualityFactors: [8 1]
Boundary: 'periodic'
SamplingFrequency: 25
Precision: 'double'
OversamplingFactor: 0
OptimizePath: 0
Получите блоки фильтров, интервалы частоты ДПФ и параметры банка фильтров.
[filters,f,fparams] = filterbank(sf);
Постройте график вейвлет, используемых при вычислении коэффициентов первого порядка. Постройте график и центральных частот вейвлета.
coefOrder = 1;
wvFilters = filters{coefOrder+1}.psift;
wvcenFrq = fparams{coefOrder+1}.omegapsi;
plot(f,wvFilters)
hold on
cf = plot(wvcenFrq,max(wvFilters),'rx');
grid on
title('Wavelet Filters')
xlabel('Hz')
ylabel('Magnitude')
legend(cf,'Center Frequencies')