wdecenergy

Энергетическое распределение разложения мультисигнала 1-D

Синтаксис

[E,PEC,PECFS] = wdecenergy(DEC) [E,PEC,PECFS,IDXSORT,LONGS] = wdecenergy(DEC,'sort') [E,PEC,PECFS] = wdecenergy(DEC,OPTSORT,IDXSIG) [E,PEC,PECFS,IDXSORT,LONGS] = wdecenergy(DEC,OPTSORT,IDXSIG)

Описание

[E,PEC,PECFS] = wdecenergy(DEC) вычисляет вектор E, который содержит энергию (L2-норма) каждого анализируемого сигнала, матричный PEC, который содержит процент энергии для каждого компонента вейвлета (приближение и детали) каждого сигнала и матричного PECFS, который содержит процент энергии для каждого коэффициента.

E (i) является энергией (L2-норма) сигнала ith.
PEC (i, 1) является процентом энергии для приближения уровня МАКСЛЕВ = DEC.level сигнала ith.
PEC (i, j), j=2..., MAXLEV+1 является процентом энергии для детали уровня (MAXLEV+1-j) сигнала ith.
PECFS (i, j), процент энергии для jth коэффициентов сигнала ith.

[E,PEC,PECFS,IDXSORT,LONGS] = wdecenergy(DEC,'sort') возвращает отсортированный PECFS (строкой) в порядке возрастания и вектор индекса IDXSORT.

Заменение 'вида' 'возрастанием' возвращает тот же результат.
Заменение 'вида' 'убыванием' возвращает PECFS, отсортированный в порядке убывания.

LONGS является вектором, содержащим длины каждого семейства коэффициентов.

[E,PEC,PECFS] = wdecenergy(DEC,OPTSORT,IDXSIG) возвращает значения для сигналов, индексы которых даны IDXSIG вектор.

[E,PEC,PECFS,IDXSORT,LONGS] = wdecenergy(DEC,OPTSORT,IDXSIG) возвращает значения для сигналов, индексы которых даны вектором IDXSIG, вектор индекса IDXSORT и LONGS, который является вектором, содержащим длины каждого семейства коэффициентов. Допустимые значения для OPTSORT не 'ни один', 'вид', 'возрастание', 'убывание'.

Примеры

свернуть все

Мультисигнал 1-D Энергетическое распределение разложения

Скрипт Open Live Script

Загрузите 23 канала данные EEG Espiga3 [1]. Каналы располагаются по столбцам. Данные производятся на уровне 200 Гц.

load Espiga3

Выполните разложение на уровне 2 с помощью db2 вейвлет.

dec = mdwtdec('c',Espiga3,2,'db2')

dec = struct with fields:
        dirDec: 'c'
         level: 2
         wname: 'db2'
    dwtFilters: [1x1 struct]
       dwtEXTM: 'per'
      dwtShift: 0
      dataSize: [995 23]
            ca: [249x23 double]
            cd: {[498x23 double]  [249x23 double]}

Вычислите энергетическое распределение.

[e,pec,pecfs] = wdecenergy(dec);

Отобразите полную энергию и распределение энергии для каждого компонента вейвлета (A2, D2, D1) во втором канале.

idx = 2;
e(idx)

ans = 8.0445e+05

perA2D2D1 = pec(idx,:)

perA2D2D1 = 1×3

   99.0126    0.8986    0.0887

Сравните содействующее энергетическое распределение для сигнала 1 и сигнала 10. Поскольку большая часть энергии находится в коэффициентах приближения, увеличение оси X количеством коэффициентов приближения.

sigA = 1;
sigB = 10;
pecfsA = pecfs(sigA,:);
pecfsB = pecfs(sigB,:);
plot(pecfsA,'r--')
hold on
plot(pecfsB,'b')
grid on
legend('pecfsA','pecfsB')
xlim([0 size(dec.ca,1)])

Ссылки

[1] Мезаструктура, Гектор. “Адаптированные Вейвлеты для Обнаружения Шаблона”. Происходящий в Распознавании образов, Анализе изображения и Приложениях, отредактированных Альберто Санфелиу и Мануелем Лазо Кортесом, 3773:933–44. Берлин, Гейдельберг: Спрингер Берлин Гейдельберг, 2005. https://doi.org/10.1007/11578079_96.

Смотрите также

mdwtdec | mdwtrec

Документация