开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > 波变换 > 查看源码
wtbarrow.m
资源名称：wavelet_toolbox.rar　[点击查看]
上传用户：haiyisale
上传日期：2013-01-09
资源大小：3246k
文件大小：9k
源码类别：

波变换

开发平台：

Matlab

wtbarrow.m：源码内容
  1. function HDL = wtbarrow(option,varargin)
  2. %WTBARROW Create Arrows.
  3. %
  4. %  example:
  5. %       hdl = wtbarrow('create');
  6. %       set(gca,'Xlim',[0,1],'Ylim',[0,1]);
  7. %       pause
  8. %       hdl = wtbarrow('create','Rotation',pi/2,'Scale',[0.5 0.5],'Color','b');
  9. %       set(gca,'Xlim',[-2,2],'Ylim',[-2,2]);
  10. %       pause
  11. %       hdl = wtbarrow('create','Translate',[1,0],'Rotation',pi/4,'Scale',[0.5 0.5],'Color','y');
  13. % Copyright 2003 The MathWorks, Inc.
  15. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 10-Feb-2003.
  16. %   Last Revision: 06-Mar-2003.
  18. if nargin<1 , option = 'create'; end
  20. % Defaults.
  21. %------------
  22. sca = ones(1,2);        
  23. rot = 0;
  24. tra = zeros(1,2);  
  26. nbarg  = length(varargin);
  27. switch option
  28.     case {'create','createBIS','createTER'}
  29.         % ------------------------------------------
  30.         % Arrow construction.
  31.         % ------------------------------------------
  32.         %          x             H = HArrow
  33.         %     <------------> <------->
  34.         %                    =             ^
  35.         %                    ===           |  W = WArrow
  36.         %                    ====          v
  37.         %     ====================      ^
  38.         %     ======================    |
  39.         %     ========================  |  2*E
  40.         %     ======================    |
  41.         %     ====================      v
  42.         %                    ====
  43.         %                    ===
  44.         %                    =
  45.         % ALFA is the angle of the arrow head
  46.         % or W = WArrow.
  47.         % Z = E + W;
  48.         % ------------------------------------------
  49.         
  50.         % defaults.
  51.         %----------
  52.         axe = gca;
  53.         vis = 'On';
  54.         col = [1 0 0]; 
  55.        %-------------
  56.         alf = pi/4;
  57.         W   = NaN;
  58.         xP1 = 3/4;
  59.         H   = 1/4;
  60.         E   = 1/8;
  61.         %-----------
  62.         xP2 = 2*xP1;
  63.         dY1 = 1/2;
  64.         %------------          
  65.         for k = 1:2:nbarg
  66.             argType = varargin{k};
  67.             argVal  = varargin{k+1};
  68.             argType = lower(argType(1:3));
  69.             switch argType
  70.                 case 'axe' , axe = argVal;
  71.                 case 'vis' , vis = argVal;
  72.                 case 'col' , col = argVal; 
  73.                 %------------------------
  74.                 case 'alf' , alf = argVal;
  75.                 case 'xp1' , xP1 = argVal;
  76.                 case 'xp2' , xP2 = argVal;    
  77.                 case 'har' , H   = argVal;
  78.                 case 'war' , W   = argVal;                    
  79.                 case 'wid' , E   = argVal;
  80.                 case 'hy1'  ,HY1 = argVal;
  81.                 %------------------------       
  82.                 case 'sca' , sca = argVal;        
  83.                 case 'rot' , rot = argVal;
  84.                 case 'tra' , tra = argVal;
  85.             end
  86.         end
  87.         if isnan(W)
  88.             Z =  tan(alf)*H;
  89.             W = max(Z-E,0);
  90.         end
  91.         Z  = E + W;
  92.         switch option
  93.             case 'create'
  94.                 ArrowData = [...
  95.                   0  xP1  xP1 xP1+H xP1 xP1 0  0; ...
  96.                  -E  -E   -Z     0   Z   E  E -E ...
  97.                     ];
  98.             case 'createBIS'
  99.                 V  = HY1-H;
  100.                 y1 = -(E+HY1);
  101.                 y2 = y1+V;
  102.                 Xd = [0 2*E  2*E  xP1-E  xP1-E  xP1+E  xP1+E  xP2-E  xP2-E  xP2-Z  xP2   xP2+E+W  xP2+E  xP2+E  0   0];
  103.                 Yd = [y1 y1  -E    -E     y1      y1     -E    -E     y2      y2    y1      y2      y1+V   E    E  y1];
  104.                 ArrowData = [Xd ; Yd];
  105.                 
  106.             case 'createTER'
  107.                 V  = HY1-H;
  108.                 y1 = -(E+HY1);
  109.                 y2 = y1+V;
  110.                 Xd = [0 2*E  2*E  xP1-E  xP1-E  xP1+E  xP1+E  xP2-E  xP2-E  xP2-Z  xP2   xP2+E+W  xP2+E  xP2+E  0   0];
  111.                 Yd = [y1 y1  -E    -E     y1      y1     -E    -E     y2      y2    y1      y2      y1+V   E    E  y1];
  112.                 N = 10;
  113.                 xx = linspace(2*E,xP1-E,N);
  114.                 t  = linspace(0,1,10);
  115.                 yy = y1 +  HY1*4*t.*(1-t);
  116.                 Xd = [Xd(1:2) xx Xd(5:end)];
  117.                 Yd = [Yd(1:2) yy Yd(5:end)];
  118.                 ArrowData = [Xd ; Yd];
  119.               
  120.         end
  121.         %-------------------------------------------
  122.         if ~(isequal(sca,ones(2,1)) | isequal(sca,ones(1,2)))
  123.             MSca = [sca(1) 0 ; 0 sca(2)];
  124.             ArrowData = MSca*ArrowData;
  125.         end
  126.         if ~isequal(rot,0)
  127.             MRot = [cos(rot) sin(rot) ; -sin(rot) cos(rot)];
  128.             ArrowData = MRot*ArrowData;
  129.         end
  130.         if ~(isequal(tra,zeros(2,1)) | isequal(tra,zeros(1,2)))
  131.             VTra = [tra(1) ; tra(2)]; 
  132.             ArrowData = ArrowData + VTra(:,ones(1,size(ArrowData,2)));
  133.         end
  134.         HDL = patch( ...
  135.             'Xdata',ArrowData(1,:),'Ydata',ArrowData(2,:),...
  136.             'FaceColor',col,'Visible',vis);
  138.     case {'special_1'}       
  139.         % defaults.
  140.         %----------
  141.         axe = gca;
  142.         vis = 'On';
  143.         col = [1 0 0]; 
  144.         %-------------
  145.         alf = pi/4;
  146.         W   = NaN;
  147.         H   = 1/4;
  148.         E   = 1/8;
  149.         %-----------
  150.         Xd = [1 2 3];
  151.         Yd = [1 2 3];
  152.         %------------          
  153.         for k = 1:2:nbarg
  154.             argType = varargin{k};
  155.             argVal  = varargin{k+1};
  156.             Larg = min([3,length(argType)]);
  157.             argType = lower(argType(1:Larg));
  158.             switch argType
  159.                 case 'axe' , axe = argVal;
  160.                 case 'vis' , vis = argVal;
  161.                 case 'col' , col = argVal; 
  162.                 %------------------------
  163.                 case 'xva' , Xd  = argVal;
  164.                 case 'yva' , Yd  = argVal;    
  165.                 case 'alf' , alf = argVal;
  166.                 case 'har' , H   = argVal;
  167.                 case 'war' , W   = argVal;                    
  168.                 case 'wid' , E   = argVal;
  169.                 %------------------------       
  170.                 case 'sca' , sca = argVal;        
  171.                 case 'rot' , rot = argVal;
  172.                 case 'tra' , tra = argVal;
  173.             end
  174.         end
  175.         if isnan(W)
  176.             Z =  tan(alf)*H;
  177.             W = max(Z-E,0);
  178.         end
  179.         Ex = E;
  180.         Ey = E;
  181.         x1 = Xd(1); x2 = Xd(2); x3 = Xd(3);
  182.         y1 = Yd(1); y2 = Yd(2); y3 = Yd(3);
  183.         XY = zeros(2,20);
  184.         XY(:,1)  = [x1-Ex ; y1];  
  185.         XY(:,2)  = [x1+Ex ; y1]; 
  186.         XY(:,3)  = [x1+Ex ; y2];  
  187.         XY(:,4)  = [x2-Ex ; y2]; 
  188.         XY(:,5)  = [x2-Ex ; y1];  
  189.         XY(:,6)  = [x2+Ex ; y1]; 
  190.         XY(:,7)  = [x2+Ex ; y2+2*Ey];
  191.         XY(:,8)  = [(x1+x2)/2 + Ex ; y2+2*Ey]; 
  192.         XY(:,9)  = [(x1+x2)/2 + Ex ; y3];
  193.         XY(:,10) = [x3-Ex ; y3]; 
  194.         XY(:,11) = [x3-Ex ; y1+H];
  195.         XY(:,12) = [x3-Ex-W ; y1+H];
  196.         XY(:,13) = [x3 ; y1];
  197.         XY(:,14) = [x3+Ex+W ; y1+H];
  198.         XY(:,15) = [x3+Ex ; y1+H];
  199.         XY(:,16) = [x3+Ex ; y3+2*Ey];
  200.         XY(:,17) = [(x1+x2)/2 - Ex ; y3+2*Ey];
  201.         XY(:,18) = [(x1+x2)/2 - Ex ; y2+2*Ey];
  202.         XY(:,19) = [x1-Ex ; y2+2*Ey];
  203.         XY(:,20) = XY(:,1);
  204.         %-------------------------------------------
  205.         if ~(isequal(sca,ones(2,1)) | isequal(sca,ones(1,2)))
  206.             MSca = [sca(1) 0 ; 0 sca(2)];
  207.             XY = MSca*XY;
  208.         end
  209.         if ~isequal(rot,0)
  210.             MRot = [cos(rot) sin(rot) ; -sin(rot) cos(rot)];
  211.             XY = MRot*XY;
  212.         end
  213.         if ~(isequal(tra,zeros(2,1)) | isequal(tra,zeros(1,2)))
  214.             VTra = [tra(1) ; tra(2)]; 
  215.             XY = XY + VTra(:,ones(1,size(XY,2)));
  216.         end
  217.         HDL = patch( ...
  218.             'Xdata',XY(1,:),'Ydata',XY(2,:),...
  219.             'FaceColor',col,'Visible',vis);
  220.         
  221.         
  222.     case 'set'
  223.         HDL = varargin{1};
  224.         ArrowData = [get(HDL,'Xdata')' ; get(HDL,'Ydata')'];
  225.         vis = get(HDL,'Visible');
  226.         col = get(HDL,'FaceColor');
  227.         for k = 2:2:nbarg
  228.             argType = varargin{k};
  229.             argVal  = varargin{k+1};
  230.             argType = lower(argType(1:3));
  231.             switch argType
  232.                 case 'vis' , vis = argVal;
  233.                 case 'col' , col = argVal; 
  234.                 case 'sca' , sca = argVal;        
  235.                 case 'rot' , rot = argVal;
  236.                 case 'tra' , tra = argVal;
  237.             end
  238.         end
  239.         if ~(isequal(sca,ones(2,1)) | isequal(sca,ones(1,2)))
  240.             MSca = [sca(1) 0 ; 0 sca(2)];
  241.             ArrowData = MSca*ArrowData;
  242.         end
  243.         if ~isequal(rot,0)
  244.             MRot = [cos(rot) sin(rot) ; -sin(rot) cos(rot)];
  245.             ArrowData = MRot*ArrowData;
  246.         end
  247.         if ~(isequal(tra,zeros(2,1)) | isequal(tra,zeros(1,2)))
  248.             VTra = [tra(1) ; tra(2)]; 
  249.             ArrowData = ArrowData + VTra(:,ones(1,size(ArrowData,2)));
  250.         end
  251.         set(HDL, ...
  252.             'Xdata',ArrowData(1,:),'Ydata',ArrowData(2,:), ...
  253.             'FaceColor',col,'Visible',vis);
  254. end