How to calculate the energy of details and approxmiations coefficients of discret wavelet transform ?

2 visualizaciones (últimos 30 días)
Hi everyone, i calculate the energy of details and approxmiations coefficients of discret wavelet transform .i used DWT(Daub4) Here is the code
for n=1:23
deriv=(val(n,:))/2.559375;
t=(0:length(deriv)-1)/Fs;
result{n} = filtfilt(d1,deriv);
[c,l]=wavedec(result{n},4,'db4');
[cd1,cd2,cd3,cd4]=detcoef(c,l,[1,2,3,4]);
cA4=appcoef(c,l,'db4',4);
[Ea,Ed]=wenergy(c,l)
energie{n} = {Ea, Ed};
end
And when i was looking for the formula of the energy i found this :
and as i understand the energy is the percentage of energy corresponding to the approximation and the details
So when i was trying to write a C code using this formula and the values of cd1,cd2,cd3,cd4 i don't get the same values that i have with matlab So i want to know where i wrong ?
i mean to calcuate the energy i should use [cd1,cd2,cd3,cd4]??
Please i need a clarification,Thanks
  1 comentario
Ewunate Assaye Kassaw
Ewunate Assaye Kassaw el 18 de Jun. de 2020
hello sir
it gives me a hint thanks for that
i have questions
  1. n stands for what?
  2. 2.559375 where this number comes from?
  3. deriv=(val(n,:))/2.559375; i dont undrstand 'val' it gives error when i run it and says fix val with eval and again gives me error
can you help me please?

Iniciar sesión para comentar.

Respuestas (2)

Jiaxin Tu
Jiaxin Tu el 28 de Feb. de 2019
That's because wenergy gives the percentage of power out of 100.
If you divide your calculations with the formula by the sum of energy *100 you will get the same answer

Ewunate Assaye Kassaw
Ewunate Assaye Kassaw el 18 de Jun. de 2020
hello sir
it gives me a hint thanks for that
i have questions
  1. n stands for what?
  2. 2.559375 where this number comes from?
  3. deriv=(val(n,:))/2.559375; i dont undrstand 'val' it gives error when i run it and says fix val with eval and again gives me error
can you help me please?

Categorías

Más información sobre Discrete Multiresolution Analysis en Help Center y File Exchange.

Etiquetas

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by