Этот пример показывает, как получить матрицу функций рассеяния для сети вейвлет рассеяния и как сравнить матрицу с коэффициентами рассеяния.

Загрузите сигнал ЭКГ, дискретизированный со скоростью 180 Гц. Создайте сеть вейвлет времени, которая может использоваться с сигналом.

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

Вычислите матрицу функции рассеяния с помощью преобразования журнала. Отобразите размерности матрицы.

smat = featureMatrix(sf,wecg,'Transform','Log');
size(smat)

ans = 1×2

   147     8

Теперь вычислите преобразование рассеяния сигнала. Получите коэффициенты рассеяния. Выходы - массив ячеек с тремя элементами. Каждый элемент является таблицей. Подтвердите, что общее количество строк в таблицах равно количеству строк в матрице.

S = scatteringTransform(sf,wecg);
t1rows = size(S{1},1);
t2rows = size(S{2},1);
t3rows = size(S{3},1);
disp(['Total Number of Rows: ',num2str(t1rows+t2rows+t3rows)])

Total Number of Rows: 147

Отобразите основу -2 журналов разрешение коэффициентов рассеяния нулевого порядка.

disp(['Resolution: ',num2str(S{1}.resolution(1))])

Resolution: -8

Получите естественный логарифм коэффициентов рассеяния нулевого порядка. Сравните коэффициенты рассеяния с первой строкой в матрице функций. Количество коэффициентов в каждом равняется абсолютному значению разрешения основы -2 журналов.

logS = log(sf,S);
logScat = logS{1}.signals{1};
[smat(1,:)' logScat]

ans = 8×2

   -1.2914   -1.2914
   -2.4682   -2.4682
   -1.6368   -1.6368
   -1.2716   -1.2716
   -1.6818   -1.6818
   -4.3701   -4.3701
   -1.3199   -1.3199
   -1.0542   -1.0542