Банки фильтров вейвлет-временного рассеяния
filters = filterbank(sf)filters - массив ячеек структурных массивов с элементами нордера, где 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 до nфильтров включительно, где nфильтры - количество банков фильтров в сети рассеяния. Эти входные аргументы могут использоваться с любым из выходных синтаксисов, показанных ранее.
Создайте сеть вейвлет-временного рассеяния для сигнала, дискретизированного на частоте 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
Получите банки фильтров, ячейки частот DFT и параметры банков фильтров.
[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')