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

波变换

开发平台:

Matlab

  1. function x = recons(t,node,x,sizes,edges)
  2. %RECONS Reconstruct wavelet coefficients.
  3. %   Y = RECONS(T,N,X,S,E) reconstructs the wavelet
  4. %   packet coefficients X associated with the node N
  5. %   of the wavelet tree T, using sizes S and the edges values E.
  6. %   S contains the size of data associated with
  7. %   each ascendant of N.
  8. %   The children of a node F are numbered from left 
  9. %   to right: [0, ... , ORDER-1].
  10. %   The edge value between F and a child C is the child number.
  11. %
  12. %   This method overloads the DTREE method.
  13. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 09-Sep-1999.
  14. %   Last Revision: 11-May-2003.
  15. %   Copyright 1995-2004 The MathWorks, Inc.
  16. %   $Revision: 1.4.4.2 $  $Date: 2004/03/15 22:36:04 $ 
  17. dwtEXTM = t.dwtMode;
  18. Lo_R  = t.waveInfo.Lo_R;
  19. Hi_R  = t.waveInfo.Hi_R;
  20. mode  = t.dwtMode;
  21. nb_up = length(edges);
  22. f     = zeros(nb_up,length(Lo_R));
  23. shift = zeros(nb_up,1);
  24. K = find(edges==0 | edges==3);
  25. if ~isempty(K)
  26.     f(K,:) = Lo_R(ones(size(K)),:);
  27. end
  28. K = find(edges==1 | edges==2);
  29. if ~isempty(K)
  30.     f(K,:) = Hi_R(ones(size(K)),:);
  31. end
  32. K = find(edges==2 | edges==3);
  33. if ~isempty(K) , shift(K) = 1; end
  34. for k=1:nb_up
  35.     s = max(sizes(k,:));
  36.     x = upsconv1(x,f(k,:),s,dwtEXTM,shift(k));
  37. end