Энергетическое распределение разложения мультисигнала 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 не 'ни один', 'вид', 'возрастание', 'убывание'.
% Load original 1D-multisignal. load thinker % Perform a decomposition at level 2 using wavelet db2. dec = mdwtdec('r',X,2,'db2'); % Compute the energy distribution. [E,PEC,PECFS] = wdecenergy(dec); % Display the total energy and the distribution of energy % for each wavelet component (A2, D2, D1). E31 = E(31) perA2D2D1 = PEC(31,:) % Compare the coefficient energy distribution % for signal 1 and signal 31. PECFS_1 = PECFS(1,:); PECFS_31 = PECFS(31,:); figure; plot(PECFS_1,'--r','linewidth',2); hold on plot(PECFS_31,'b','linewidth',2); grid; set(gca,'Xlim',[1,size(PECFS,2)]) title('PECFS1 dashed line and PECFS31 solid line')
