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

Загрузите сигнал ECG, произведенный на уровне 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