Этот пример показывает, как переопределить 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 (потому что $${\log }_{2}2=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')