wenergy.m
上传用户:haiyisale
上传日期:2013-01-09
资源大小:3246k
文件大小:1k
源码类别:

波变换

开发平台:

Matlab

  1. function E = wenergy(t)
  2. %WENERGY Energy for a wavelet packet decomposition.
  3. %   For a wavelet packet tree T, (see WPTREE, WPDEC, WPDEC2) 
  4. %   E = WENERGY(T) returns a vector E, which contains the
  5. %   percentages of energy corresponding to the terminal nodes
  6. %   of the tree T. 
  7. %
  8. %   Examples:
  9. %     load noisbump
  10. %     T = wpdec(noisbump,3,'sym4');
  11. %     E = wenergy(T)
  12. %     ------------------------------
  13. %     load detail
  14. %     T = wpdec2(X,2,'sym4');
  15. %     E = wenergy(T)
  16. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 12-Mar-96.
  17. %   Last Revision: 14-May-2003.
  18. %   Copyright 1995-2004 The MathWorks, Inc.
  19. %   $Revision: 1.2.4.2 $
  20. C = read(t,'allcfs');
  21. Et = sum(C(:).^2);
  22. tn = leaves(t,'s');
  23. for k=1:length(tn)
  24.     C = read(t,'data',tn(k));
  25.     E(k) = sum(C(:).^2);
  26. end
  27. E = 100*E/Et;