Вейвлет, рассеивающий полосу пропускания, сосредотачивает частоты
F = centerFrequencies(sf)F = centerFrequencies(sf,filterbanks)F = centerFrequencies(sf)sf. Вывод F является массивом ячеек с элементами Nfb, где Nfb является количеством рассеивания наборов фильтров. Каждый элемент F является вектором с действительным знаком. Если вы задаете частоту дискретизации в sf, F находится в герц. Если вы не задаете частоту дискретизации, F находится в циклах/выборке.
Если существует только один набор фильтров в рассеивающейся среде, F является вектором с действительным знаком, содержащим частоты центра полосы пропускания вейвлета.
F = centerFrequencies(sf,filterbanks)filterbanks. Аргумент filterbanks является скаляром или вектором со всеми элементами между 1 и Nfb включительно, где Nfb является количеством рассеивания наборов фильтров.
Создайте рассеивающуюся среду разложения с частотой дискретизации 50 Гц.
sf = waveletScattering('SamplingFrequency',50)sf =
waveletScattering with properties:
SignalLength: 1024
InvarianceScale: 10.2400
QualityFactors: [8 1]
Boundary: "periodic"
SamplingFrequency: 50
Precision: "double"
OversamplingFactor: 0
Постройте частоты центра полосы пропускания вейвлета для всех наборов фильтров.
bpcf = centerFrequencies(sf);
plot(bpcf{1},'rx-')
hold on
plot(bpcf{2},'bo-')
grid on
title('Wavelet Bandpass Center Frequencies')
legend('Filter Bank 1','Filter Bank 2')
ylabel('Hz')
Постройте фильтры вейвлетов, используемые в вычислении коэффициентов рассеивания второго порядка.
orderCoef = 2;
[filters,f] = filterbank(sf);
figure
plot(f,filters{orderCoef+1}.psift)
grid on
title('Wavelet Filters with Q = 1')
xlabel('Hz')
ylabel('Magnitude')
hold on
pl = plot(bpcf{orderCoef},max(filters{orderCoef+1}.psift),'bo');
legend(pl,'Center Frequencies')
.