Получение MODWT сигнала электрокардиограммы (ЭКГ) с использованием значения по умолчанию sym4 вейвлет вниз до максимального уровня. Данные взяты из Percival & Walden (2000), стр. 125 (данные первоначально предоставлены Уильямом Константином и Пером Рейнхаллом, Вашингтонский университет).

load wecg;
wtecg = modwt(wecg);
whos wtecg

  Name        Size               Bytes  Class     Attributes

  wtecg      12x2048            196608  double              

Первые одиннадцать рядов wtecg - вейвлет-коэффициенты для шкал $${\mathrm{}}^{21}$$- $${\mathrm{}}^{\mathrm{}211}$$. Конечная строка содержит масштабные коэффициенты в масштабе $${\mathrm{}}^{\mathrm{}211}$$. Постройте график коэффициентов детализации (вейвлет) для шкалы $${\mathrm{}}^{23}$$.

plot(wtecg(3,:))
title('Level 3 Wavelet Coefficients')

Получение данных MODWT индекса южных колебаний с помощью db2 вейвлет вниз до максимального уровня.

load soi;
wsoi = modwt(soi,'db2');

Получение MODWT данных обменного курса Deutsche Mark - доллар США с помощью масштабирования длины 8 и вейвлет-фильтров Фейера-Коровкина.

load DM_USD;
[Lo,Hi] = wfilters('fk8');
wdm = modwt(DM_USD,Lo,Hi);

Получить MODWT сигнала ЭКГ до шкалы $${\mathrm{}}^{24}$$, соответствующей четвертому уровню. Использовать значение по умолчанию sym4 вейвлет. Данные взяты из Percival & Walden (2000), стр. 125 (данные первоначально предоставлены Уильямом Константином и Пером Рейнхаллом, Вашингтонский университет).

load wecg;
wtecg = modwt(wecg,4);
whos wecg wtecg

  Name          Size              Bytes  Class     Attributes

  wecg       2048x1               16384  double              
  wtecg         5x2048            81920  double              

Размер строки wtecg является L + 1, где в этом случае уровень (L) равен 4. Размер столбца соответствует количеству входных выборок.

Получение MODWT сигнала ЭКГ с помощью обработки границ отражения. Использовать значение по умолчанию sym4 вейвлет и получить преобразование до уровня 4. Данные взяты из Percival & Walden (2000), стр. 125 (данные первоначально предоставлены Уильямом Константином и Пером Рейнхаллом, Вашингтонский университет).

load wecg;
wtecg = modwt(wecg,4,'reflection');
whos wecg wtecg

  Name          Size               Bytes  Class     Attributes

  wecg       2048x1                16384  double              
  wtecg         5x4096            163840  double              

wtecg имеет 4096 столбцов, что в два раза больше длины входного сигнала, wecg.

Загрузка 23-канальных данных ЭЭГ Espiga3 [3]. Каналы расположены столбчато. Данные дискретизируются на частоте 200 Гц.

load Espiga3

Вычислите максимальное перекрывающееся дискретное вейвлет-преобразование до максимального уровня.

wt = modwt(Espiga3);

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

sigN2 = vecnorm(Espiga3).^2;
wtN2 = sum(squeeze(vecnorm(wt,2,2).^2));
bar(1:23,log(sigN2))
hold on
scatter(1:23,log(wtN2),'filled','SizeData',100)
alpha(0.75)
legend('Signal Energy','Energy in Wavelet Coefficients', ...
        'Location','NorthWest')
xlabel('Channel')
ylabel('ln(squared energy)')

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

Выберите sym6 вейвлет. Загрузить и построить график сигнала электрокардиограммы (ЭКГ). Частота дискретизации сигнала ЭКГ составляет 180 герц. Данные взяты из Percival and Walden (2000), стр. 125 (данные первоначально предоставлены Уильямом Константином и Пером Рейнхаллом, Вашингтонский университет).

load wecg
t = (0:numel(wecg)-1)/180;
wv = 'sym6';
plot(t,wecg)
grid on
title(['Signal Length = ',num2str(numel(wecg))])
xlabel('Time (s)')
ylabel('Amplitude')

Возьмите MODWT сигнала.

wtecg = modwt(wecg,wv);

Входные данные являются выборками функции $$f(x$$), оцениваемой в $$$$N-кратных временных точках. Функция может быть выражена в виде линейной комбинации масштабной функции ($$x)$$и $$\mathrm{вейвлет\; (}\mathrm{x)}$$при различных масштабах и $$трансляциях:f\underset{(x\mathrm{)}}{\overset{\mathrm{}}{}}{}_{}{\mathrm{}}^{{}_{\mathrm{}}\mathrm{}}{\mathrm{}}^{{}_{\mathrm{}}}\mathrm{=\sum k=0N-1ck2-J0/2\varphi }(\underset{\mathrm{}}{\overset{{\mathrm{2-J0x-k}}_{\mathrm{)}}}{}}{}_{}$$+∑j=1J0fj ($$x_{\mathrm{),}}\underset{\mathrm{где}\mathrm{}}{\overset{\mathrm{fj}\mathrm{(}}{}}{x}_{)}{\mathrm{}}^{\mathrm{}}{\mathrm{}}^{}$$=∑k=0N-1dj,k2-j/2ψ ($${}_{\mathrm{}}$$2-jx-k) и J0 - количество уровней вейвлет-разложения. Первая сумма является грубой масштабной аппроксимацией $${}_{\mathrm{сигнала}},$$а fj (x) являются деталями на последовательных шкалах. MODWT $$$$возвращает $${}_{}$$коэффициенты числа N $${\{}_{\mathrm{}}ck\}$$и коэффициенты $${}_{\mathrm{детализации\; (J0\; \times \; N)}}$${dj, k} расширения. Каждая строка вwtecg содержит коэффициенты с другой шкалой.

При взятии MODWT сигнала длиной $$N$$существуют уровни $${\mathrm{}}_{\mathrm{}}$$разложения floor (log2 (N)) -many (по умолчанию). Коэффициенты детализации формируются на каждом уровне. Коэффициенты масштабирования возвращаются только для конечного уровня. В этом примере, $$поскольку\mathrm{}\mathrm{N}\mathrm{}\mathrm{=}$$$${}_{\mathrm{}}2048\text{, J0}=\mathrm{}\mathrm{}\mathrm{floor}\mathrm{}\mathrm{(}\mathrm{}log2\mathrm{}\mathrm{(}$$2048)) = 11 и количество строк вwtecg является $${}_{\mathrm{J0}}\mathrm{+}1\mathrm{}=\mathrm{}11\mathrm{}$$+ 1 = 12.

MODWT разделяет энергию по различным масштабам и коэффициентам масштабирования: $${}^{\mathrm{}}\underset{\mathrm{}}{\overset{{}_{\mathrm{}}}{}}{{}_{}}^{\mathrm{}}{{}_{{}_{\mathrm{}}}}^{\mathrm{||X||2=\sum j=1J0||Wj||2+||VJ0||2}}$$, где $$X$$ - входные данные, $${}_{\mathrm{Wj}}$$ - коэффициенты детализации в масштабе $$j$$, и $${}_{{}_{\mathrm{VJ0}}}$$ - коэффициенты масштабирования конечного уровня.

Вычислите энергию в каждом масштабе и вычислите их сумму.

energy_by_scales = sum(wtecg.^2,2);
Levels = {'D1';'D2';'D3';'D4';'D5';'D6';'D7';'D8';'D9';'D10';'D11';'A11'};
energy_table = table(Levels,energy_by_scales);
disp(energy_table)

    Levels     energy_by_scales
    _______    ________________

    {'D1' }         14.063     
    {'D2' }         20.612     
    {'D3' }         37.716     
    {'D4' }         25.123     
    {'D5' }         17.437     
    {'D6' }         8.9852     
    {'D7' }         1.2906     
    {'D8' }         4.7278     
    {'D9' }         12.205     
    {'D10'}         76.428     
    {'D11'}         76.268     
    {'A11'}         3.4192     

energy_total = varfun(@sum,energy_table(:,2))

energy_total=table
    sum_energy_by_scales
    ____________________

           298.28       

Подтвердите, что MODWT сохраняет энергию, вычисляя энергию сигнала и сравнивая ее с суммой энергий по всем масштабам.

energy_ecg = sum(wecg.^2);
max(abs(energy_total.sum_energy_by_scales-energy_ecg))

ans = 7.4402e-10

Возьмите MODWTMRA сигнала.

mraecg = modwtmra(wtecg,wv);

MODWTMRA возвращает проекции функции $$f(x$$) на различные вейвлет-подпространства и конечное пространство масштабирования. То есть MODWTMRA возвращает $$\underset{\mathrm{}}{\overset{\mathrm{}}{}}{}_{}{\mathrm{}}^{{}_{\mathrm{}}\mathrm{}}\sum k=0N-1ck2-J0/2\varphi {\mathrm{}}^{({}_{\mathrm{}}}$$2-J0x-k) и $${}_{\mathrm{}}$$$${\mathrm{J0-many}}_{}\{fj$$(x)}, $$\mathrm{оцененные}$$в N-кратных временных точках. Каждая строка вmraecg - проекция $$f(x$$) на другое подпространство. Это означает, что исходный сигнал может быть восстановлен путем сложения всех проекций. Это неверно в случае MODWT. Добавление коэффициентов вwtecg не восстановит исходный сигнал.

Выберите момент времени, добавьте проекции $$f(x$$), вычисленные в этот момент времени, и сравните с исходным сигналом.

time_point = 1000;
abs(sum(mraecg(:,time_point))-wecg(time_point))

ans = 3.0846e-13

Подтвердите, что, в отличие от MODWT, MODWTMRA не является энергосберегающим преобразованием.

energy_ecg = sum(wecg.^2);
energy_mra_scales = sum(mraecg.^2,2);
energy_mra = sum(energy_mra_scales);
max(abs(energy_mra-energy_ecg))

ans = 115.7053

MODWTMRA - это нулевая фазовая фильтрация сигнала. Элементы будут выровнены по времени. Продемонстрируйте это, построив график исходного сигнала и одного из его проекций. Чтобы лучше проиллюстрировать выравнивание, увеличьте изображение.

plot(t,wecg,'b')
hold on
plot(t,mraecg(4,:),'-')
hold off
grid on
xlim([4 8])
legend('Signal','Projection','Location','northwest')
xlabel('Time (s)')
ylabel('Amplitude')

Создайте аналогичный график с использованием коэффициентов MODWT в том же масштабе. Обратите внимание, что элементы не будут выровнены по времени. MODWT не является нулевой фазовой фильтрацией входного сигнала.

plot(t,wecg,'b')
hold on
plot(t,wtecg(4,:),'-')
hold off
grid on
xlim([4 8])
legend('Signal','Coefficients','Location','northwest')
xlabel('Time (s)')
ylabel('Amplitude')