В этом примере показано, как выполнить избыточную выборку 1-D преобразования вейвлет-рассеяния.

Загрузите сигнал ЭКГ, дискретизированный на частоте 180 Гц, и создайте сеть вейвлет-временного рассеяния для обработки сигнала. Чтобы выполнить критически пониженное преобразование вейвлет-рассеяния, не изменяйте значение OversamplingFactor свойство в sf. Возвращает количество коэффициентов рассеяния для сети рассеяния.

load wecg
Fs = 180;
sf = waveletScattering('SignalLength',numel(wecg),'SamplingFrequency',Fs);
ncf = numCoefficients(sf)

ncf = 8

Возврат 1-D преобразования вейвлет-рассеяния wecgи постройте график коэффициентов рассеяния нулевого порядка. Подтвердить, что число коэффициентов рассеяния нулевого порядка равно ncf.

s = scatteringTransform(sf,wecg);
display(['Number of zeroth-order scattering coefficients: ',...
  num2str(numel(s{1}.signals{1}))])

Number of zeroth-order scattering coefficients: 8

plot(s{1}.signals{1},'x-')
grid on
axis tight
title('Zeroth-Order Scattering Coefficients')

Чтобы переопределить коэффициенты рассеяния на коэффициент 2, установите значение OversamplingFactor имущество sf равно 1 (потому что $${\mathrm{}}_{\mathrm{}}\mathrm{log22}\mathrm{=}$$1). Возвращает количество коэффициентов рассеяния для отредактированной сети. Подтвердите, что число коэффициентов рассеяния удвоилось.

sf.OversamplingFactor = 1;
ncf = numCoefficients(sf)

ncf = 16

Возврат преобразования вейвлет-рассеяния wecg используя отредактированную сеть, постройте график коэффициентов рассеяния нулевого порядка. Поскольку число коэффициентов в критически дискретизированном преобразовании равно 8, подтвердите, что число коэффициентов нулевого порядка в избыточно дискретизированном преобразовании равно 16.

s = scatteringTransform(sf,wecg);
figure
plot(s{1}.signals{1},'x-')
grid on
axis tight
title('Zeroth-Order Scattering Coefficients')