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

波变换

开发平台:

Matlab

  1. function varargout = plot(t,varargin);
  2. %PLOT Plot RWVTREE object.
  3. %   PLOT(T) plots the RWVTREE object T.
  4. %   FIG = PLOT(T) returns the handle of the figure, which
  5. %   contains the tree T.
  6. %   PLOT(T,FIG) plots the tree T in the figure FIG, which
  7. %   already contains a tree.
  8. %
  9. %   PLOT is a graphical tree-management utility. The figure
  10. %   that contains the tree is a GUI tool. It lets you change
  11. %   the Node Label to Depth_Position or Index, and Node Action
  12. %   to Split-Merge or Visualize.
  13. %   The default values are Depth_Position and Visualize.
  14. %
  15. %   You can click the nodes to execute the current Node Action.
  16. %
  17. %   After some split or merge actions you can get the new tree
  18. %   using the handle of the figure, which contains it.
  19. %   You must use the following special syntax:
  20. %       NEWT = PLOT(T,'read',FIG).
  21. %   In fact, the first argument is dummy. Then the most general
  22. %   syntax for this purpose is:
  23. %       NEWT = PLOT(DUMMY,'READ',FIG);
  24. %   where DUMMY is any object parented by an NTREE object.
  25. %
  26. %   DUMMY can be any object constructor name, which returns
  27. %   an object parented by an NTREE object. For example:
  28. %      NEWT = PLOT(ntree,'read',FIG);
  29. %      NEWT = PLOT(dtree,'read',FIG);
  30. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 09-Sep-1999.
  31. %   Last Revision: 21-May-2003.
  32. %   Copyright 1995-2004 The MathWorks, Inc.
  33. %   $Revision: 1.5.4.2 $  $Date: 2004/03/15 22:36:07 $ 
  34. nbin = length(varargin);
  35. fig_tree = NaN;
  36. switch nbin
  37.   case 0    , option = 'create';
  38.   otherwise
  39.     option = varargin{1};
  40.     if isnumeric(option)
  41.         fig_tree = option;
  42.         option = 'create';
  43.     end
  44. end
  45. switch option
  46.   case {'create'}
  47.   case {'length','type'} , fig_tree = varargin{2};  % Added Node Labels
  48. end
  49. switch option
  50.   case 'create'
  51. order = treeord(t);
  52. fig_tree = plot(getwtbo(t,'wtree'),fig_tree);
  53.     if nargout>0 , varargout{1} = fig_tree; end
  54.     set(fig_tree,'NumberTitle','Off','Name',...
  55.         ['Fig ' int2str(fig_tree) ...
  56.          ' - Right Wavelet Tree Object of Order ' int2str(order)]);
  57.     % Store the RWVTREE.
  58.     %-------------------
  59.     plot(dtree,'write',fig_tree,t);
  60.     % Add Node Label menus.
  61.     %----------------------
  62.     plot(ntree,'addNodeLabel',fig_tree,'length');
  63.     plot(ntree,'addNodeLabel',fig_tree,'type');
  64.     % Set default Node Label to 'Index'.
  65.     %-----------------------------------
  66.     plot(ntree,'setNodeLabel',fig_tree,'Index');
  67.   case {'length','type'}
  68.     t = plot(ntree,'read',fig_tree);
  69.     if nbin<3 , nodes = allnodes(t); else , nodes = varargin{3}; end
  70.     n = length(nodes);
  71.     switch option
  72.       case 'length' , labtype = 's';
  73.       case 'type'   , labtype = 't';
  74.     end
  75.     labels = tlabels(t,labtype,nodes);
  76.     err    = ~isequal(n,size(labels,1));
  77.     varargout = {labels,err};
  78.   otherwise
  79.     try
  80.       nbout = nargout;
  81.       varargout{1:nbout} = plot(dtree,varargin{:});       
  82.     end
  83.     
  84. end
  85. %-------------------------------------------------------------------------%
  86. % Internal Functions                                                      %
  87. %-------------------------------------------------------------------------%
  88. function labels = tlabels(t,varargin)
  89. labtype = varargin{1};
  90. if length(varargin)<2
  91.     nodes = allnodes(t);
  92. else
  93.     nodes = varargin{2};
  94. end
  95. nbnodes = length(nodes);
  96. labels  = [];
  97. order = treeord(t);
  98. switch labtype
  99.    case 's'
  100.      sizes = read(t,'sizes',nodes);
  101.      switch order
  102.        case 2         
  103.          for k=1:nbnodes
  104.            labels = strvcat(labels,sprintf('%0.f',max(sizes(k,:))));
  105.          end
  106.        case 4
  107.          for k=1:nbnodes
  108.            labels = strvcat(labels,sprintf('(%0.f,%0.f)',sizes(k,:)));
  109.          end
  110.      end
  111.    case 't'
  112.      [d,p] = ind2depo(order,nodes);
  113.      p = rem(p,order);
  114.      pstr = repLine('a',nbnodes);
  115.      if order==2
  116.          I = find(p==1); pd = repLine('d',length(I)); pstr(I,:) = pd;
  117.      else
  118.          I = find(p==1); pd = repLine('h',length(I)); pstr(I,:) = pd;
  119.          I = find(p==2); pd = repLine('v',length(I)); pstr(I,:) = pd;
  120.          I = find(p==3); pd = repLine('d',length(I)); pstr(I,:) = pd;
  121.      end
  122.      lp = repLine('(',nbnodes);
  123.      rp = repLine(')',nbnodes);
  124.      labels = [lp pstr rp];
  125.      labels(1:end,:) = labels([2:end,1],:);
  126. end
  127. %-------------------------------------------------------------------------%
  128. function m = repLine(c,n)
  129. %REPLINE Replicate Lines.
  130. m = c(ones(n,1),:);
  131. %-------------------------------------------------------------------------%