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

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

波变换

开发平台：

Matlab

wvdtool.m：源码内容
  1. function varargout = wvdtool(option,varargin)
  2. %WVDTOOL Wavelet display tool.
  3. %   VARARGOUT = WVDTOOL(OPTION,VARARGIN)
  5. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 12-Mar-96.
  6. %   Last Revision: 17-Jan-2001.
  7. %   Copyright 1995-2002 The MathWorks, Inc.
  8. %   $Revision: 1.22 $ $Date: 2002/06/17 12:18:39 $
  10. % Test inputs.
  11. %-------------
  12. if nargin==0 , option = 'create'; end
  13. [option,winAttrb] = utguidiv('ini',option,varargin{:});
  15. % MemBloc1 of stored values.
  16. %---------------------------
  17. n_miscella     = 'WvDisp_Misc';
  18. ind_graph_area = 1;
  19. ind_wave_fam   = 2;
  20. ind_wave_nam   = 3;
  21. ind_refinment  = 4;
  22. nb1_stored     = 4;
  24. % Tag property of objects.
  25. %-------------------------
  26. tag_prec_val  = 'Prec_Val';
  27. tag_cmd_frame = 'Cmd_Frame';
  28. tag_display   = 'Display';
  29. tag_pus_inf1  = 'Pus_Inf1';
  30. tag_pus_inf2  = 'Pus_Inf2';
  32. switch option
  33.     case 'create'
  34.         % Get Globals.
  35.         %-------------
  36.         [Def_Txt_Height,Def_Btn_Height,Def_Btn_Width,  ...
  37.          X_Spacing,Y_Spacing,Def_FraBkColor] = ...
  38.             mextglob('get',...
  39.                 'Def_Txt_Height','Def_Btn_Height','Def_Btn_Width', ...
  40.                 'X_Spacing','Y_Spacing','Def_FraBkColor' ...
  41.                 );
  43.         % Window initialization.
  44.         %----------------------
  45.         win_title = 'Wavelet Display';
  46.         [win_loctool,pos_win,win_units,str_numwin,...
  47.                 frame0,pos_frame0,Pos_Graphic_Area,pus_close] = ...
  48.                     wfigmngr('create',win_title,winAttrb,'ExtFig_WDisp',mfilename);
  49.         varargout{1} = win_loctool;
  51. % Add Help for Tool.
  52. %------------------
  53. wfighelp('addHelpTool',win_loctool,'Wavelet Dis&play','WVDI_GUI');
  55. % Add Help Item.
  56. %----------------
  57. wfighelp('addHelpItem',win_loctool,'Available Wavelets','WVDI_FAMILIES');
  58. wfighelp('addHelpItem',win_loctool,'Wavelet Families','WVDI_ADVANCED');
  59. wfighelp('addHelpItem',win_loctool,'Properties of Wavelet Families','UT_WAVELET');
  60. wfighelp('addHelpItem',win_loctool,'Adding Your Own Wavelets','WVDI_NEWWAVE');
  62.         % Begin waiting.
  63.         %---------------
  64.         set(wfindobj('figure'),'Pointer','watch');
  66.         % General parameters initialization.
  67.         %-----------------------------------
  68.         dx = X_Spacing;
  69.         dy = Y_Spacing; dy2 = 2*dy;
  70.         d_txt = (Def_Btn_Height-Def_Txt_Height);
  72.         % Position property of objects.
  73.         %------------------------------
  74.         xlocINI      = pos_frame0([1 3]);
  75.         ybottomINI   = pos_win(4)-1.5*Def_Btn_Height-dy2;
  76.         y_low        = ybottomINI-2*Def_Btn_Height;
  77.         x_left       = pos_frame0(1)+2*dx;
  78.         pos_prec_txt = [x_left, y_low+d_txt/2, Def_Btn_Width, Def_Txt_Height];
  79.         xleftB       = x_left+Def_Btn_Width+dx;
  80.         pos_prec_val = [xleftB , y_low , Def_Btn_Width , Def_Btn_Height];
  81.         wpop         = 3*(pos_frame0(3)-4*dx)/4;
  82. xborder      = (pos_frame0(3)-wpop)/2;
  83.         xpop         = pos_frame0(1)+xborder;
  84.         y_low        = pos_prec_txt(2)-5*Def_Btn_Height;
  85.         pos_display  = [xpop, y_low, wpop, 2*Def_Btn_Height];
  86.         y_low        = pos_display(2)-3*Def_Btn_Height;
  87.         pos_inf_txt  = [x_left, y_low, 2.5*Def_Btn_Width, Def_Btn_Height];
  88.         pos_pus_inf1    = pos_display;
  89. pos_pus_inf1(1) = pos_pus_inf1(1)-xborder/2;
  90. pos_pus_inf1(3) = pos_pus_inf1(3)+xborder;
  91.         pos_pus_inf1(2) = pos_inf_txt(2)-2*Def_Btn_Height-dy;
  92.         pos_pus_inf2    = pos_pus_inf1;
  93.         pos_pus_inf2(2) = pos_pus_inf1(2)-3*Def_Btn_Height;
  95.         % String property of objects.
  96.         %----------------------------
  97.         str_display  = 'Display';
  98.         str_inf_txt  = 'Information on:';
  99.         str_inf1     = 'Haar Family (HAAR)';
  100.         str_inf2     = 'All Wavelet Families';
  101.         str_prec_txt = 'Refinement';
  102.         str_prec_val = ['5 ' ; '6 ' ; '7 ' ; '8 ' ; '9 ' ; '10' ; '11' ; '12' ];
  104.         % Command part construction of the window.
  105.         %-----------------------------------------
  106.         utanapar('create',win_loctool, ...
  107.                  'xloc',xlocINI,'bottom',ybottomINI,...
  108.                  'datflag',0,'levflag',0,...
  109.                  'wtype','all' ...
  110.                  );
  112.         comFigProp = {'Parent',win_loctool,'Unit',win_units};
  113.         comPopProp = {comFigProp{:},'Style','Popup'};
  114.         comPusProp = {comFigProp{:},'Style','Pushbutton'};
  115.         comTxtProp = {comFigProp{:},'Style','Text', ...
  116.            'HorizontalAlignment','left','Backgroundcolor',Def_FraBkColor};
  118.         Tooltip      = 'The wavelet is computed on a 2^Refinement grid';
  119.         txt_prec_txt = uicontrol(comTxtProp{:}, ...
  120.                                  'Position',pos_prec_txt,...
  121.                                  'String',str_prec_txt,...
  122.                                  'Tooltip',Tooltip ...
  123.                                   );
  124.         pop_prec_val = uicontrol(comPopProp{:},...
  125.                                  'Position',pos_prec_val,...
  126.                                  'String',str_prec_val,...
  127.                                  'Value',4,...
  128.                                  'Tooltip',Tooltip,...
  129.                                  'Tag',tag_prec_val...
  130.                                  );
  131.         pus_display  = uicontrol(comPusProp{:},...
  132.                                  'Position',pos_display,...
  133.                                  'String',xlate(str_display),...
  134.                                  'Tag',tag_display...
  135.                                  );
  136.         txt_inf_txt  = uicontrol(comTxtProp{:}, ...
  137.                                  'Position',pos_inf_txt,...
  138.                                  'String',str_inf_txt...
  139.                                  );
  140.         pus_inf1     = uicontrol(comPusProp{:},...
  141.                                  'Position',pos_pus_inf1,...
  142.                                  'String',xlate(str_inf1),...
  143.  'Max',2,...
  144.                                  'Tag',tag_pus_inf1...
  145.                                  );
  146.         pus_inf2     = uicontrol(comPusProp{:},...
  147.                                  'Position',pos_pus_inf2,...
  148.                                  'String',xlate(str_inf2),...
  149.                                  'Tag',tag_pus_inf2...
  150.                                  );
  152.         % Callbacks update.
  153.         %------------------
  154.         utanapar('set_cba_num',win_loctool,[pus_display]);
  155.         [pop_fam,pop_num] = utanapar('handles',win_loctool,'fam','num');
  156.         cb_fam = get(pop_fam,'Callback');
  157.         cb_num = get(pop_num,'Callback');
  158.         cba_upd_fam = [mfilename '(''upd_fam'',' str_numwin ');'];
  159.         cba_update  = [mfilename '(''new'',' str_numwin ');'];
  160.         cba_draw_1d = [mfilename '(''draw1d'',' str_numwin ');'];
  161.         cba_inf1    = [mfilename '(''inf1'',' str_numwin ');'];
  162.         cba_inf2    = [mfilename '(''inf2'',' str_numwin ');'];
  163.         set(pop_fam,'Callback',[cb_fam , cba_upd_fam]);
  164.         set(pop_num,'Callback',[cb_num , cba_update]);
  165.         set(pop_prec_val,'Callback',cba_update);
  166.         set(pus_display,'Callback',cba_draw_1d);
  167.         set(pus_inf1,'Callback',cba_inf1);
  168.         set(pus_inf2,'Callback',cba_inf2);
  170.         % Setting units to normalized.
  171.         %-----------------------------
  172.         Pos_Graphic_Area = Pos_Graphic_Area./[pos_win(3:4),pos_win(3:4)];
  173.         set(findobj(win_loctool,'Units','pixels'),'Units','normalized');
  175. % Add Context Sensitive Help (CSHelp).
  176. %-------------------------------------
  177. hdl_WVDI_GUI = [txt_prec_txt,pop_prec_val];
  178. wfighelp('add_ContextMenu',win_loctool,hdl_WVDI_GUI,'WVDI_GUI');
  179. %-------------------------------------
  181.         % Memory for stored values.
  182.         %--------------------------
  183.         wmemtool('ini',win_loctool,n_miscella,nb1_stored);
  184.         wmemtool('wmb',win_loctool,n_miscella,...
  185.                        ind_graph_area,Pos_Graphic_Area,...
  186.                        ind_wave_fam,'haar', ...
  187.                        ind_wave_nam,'haar', ...
  188.                        ind_refinment,0      ...
  189.                        );
  191.         % End waiting.
  192.         %-------------
  193.         set(wfindobj('figure'),'Pointer','arrow');
  195.     case 'upd_fam'
  196.         %**********************************************************%
  197.         %** OPTION = 'upd_fam' -  UPDATE OF THE WAVELET FAMILY   **%
  198.         %**********************************************************%
  199.         win_loctool = varargin{1};
  200.         new = wvdtool('new',win_loctool);
  201.         if new==0 , return; end
  203.         % Handles of tagged objects.
  204.         %---------------------------
  205.         pus_handles = findobj(win_loctool,'Style','pushbutton');
  206.         pus_inf1    = findobj(pus_handles,'Tag',tag_pus_inf1);
  208.         % Test family and Set visible on or off the wavelet number if exists.
  209.         %--------------------------------------------------------------------
  210.         wav_nam = cbanapar('get',win_loctool,'wav');
  211.         [wav_fn,wav_fsn,tabNums] = wavemngr('fields',wav_nam,'fn','fsn','tabNums');
  212.         if size(tabNums,1)<2 ,add = ' wavelet'; else , add = ' wavelets'; end
  213. strPush = [wav_fn ' Family (' upper(wav_fsn) ')' ];
  214.         set(pus_inf1,'String',strPush);
  216.     case 'inf1'
  217.         %*****************************************%
  218.         %** OPTION = 'inf1' - LOCAL INFORMATION **%
  219.         %*****************************************%
  220.         win_loctool = varargin{1};
  222.         % Getting wavelet.
  223.         %-----------------
  224.         wav_nam = cbanapar('get',win_loctool,'wav');
  225.         wav_fam = wavemngr('fam_num',wav_nam);
  226.         infotxt = [deblankl(wav_fam) 'info.m'];
  227.         [old_info,fig] = whelpfun('getflag');
  228.         if ~isempty(old_info) & strcmp(infotxt,old_info)
  229.             figure(fig); return;
  230.         end
  232.         % Waiting message.
  233.         %-----------------
  234.         wwaiting('msg',win_loctool,'Wait ... loading');
  236.         [str_inf,fid] = wreadinf(infotxt,'noerror');
  237.         if fid==-1
  238.             msg = sprintf('File %s not found !', infotxt);
  239.             errargt(mfilename,msg,'msg');
  240.             wwaiting('off',win_loctool);
  241.             return
  242.         else
  243.             dim     = size(str_inf);
  244.             rowfam  = str_inf(1,:);
  245.             str_inf = str_inf(2:dim(1),:);
  246.             col = 1;
  247.             while all(str_inf(:,col)==' ') , col = col+1; end
  248.             blk  = ' ' ;
  249.             str_inf = [rowfam ; ...
  250.                        str_inf(:,col:dim(2)) blk*ones(dim(1)-1,col-1) ];
  251.         end
  252.         ftnsize = wmachdep('fontsize','winfo');
  253.         whelpfun('create',str_inf,infotxt,ftnsize);
  255.         % End waiting.
  256.         %-------------
  257.         wwaiting('off',win_loctool);
  259.     case 'inf2'
  260.         %*****************************************%
  261.         %** OPTION = 'inf2' - LOCAL INFORMATION **%
  262.         %*****************************************%
  263.         win_loctool = varargin{1};
  265.         % Handles of tagged objects.
  266.         %---------------------------
  267.         infotxt        = 'infowave.m';
  268.         [old_info,fig] = whelpfun('getflag');
  269.         if ~isempty(old_info) & strcmp(infotxt,old_info)
  270.             figure(fig); return;
  271.         end
  273.         % Waiting message.
  274.         %-----------------
  275.         wwaiting('msg',win_loctool,'Wait ... loading');
  277.         [str_inf,fid]   = wreadinf(infotxt,'noerror');
  278.         if fid==-1
  279.             msg = sprintf('File %s not found !', infotxt);
  280.             errargt(mfilename,msg,'msg');
  281.             wwaiting('off',win_loctool);
  282.             return
  283.         end
  284.         ftnsize = wmachdep('fontsize','winfo');
  285.         whelpfun('create',str_inf,infotxt,ftnsize);
  287.         % End waiting.
  288.         %-------------
  289.         wwaiting('off',win_loctool);
  291.     case 'draw1d'
  292.         %************************************************%
  293.         %** OPTION = 'draw1d' - DRAW AXES IN 1D        **%
  294.         %************************************************%
  295.         win_loctool = varargin{1};
  296.         [new,Wave_Fam,Wave_Nam,prec_val] = wvdtool('new',win_loctool);
  297.         if new==0 , return; end
  299.         % Waiting message.
  300.         %-----------------
  301.         wwaiting('msg',win_loctool,'Wait ... computing');
  303.         % Handles of tagged objects.
  304.         %---------------------------
  305.         fra_handles = findobj(win_loctool,'Style','frame');
  306.         hdl_frame0  = findobj(fra_handles,'Tag',tag_cmd_frame);
  308.         % Update parameters selection before drawing.
  309.         %-------------------------------------------
  310.         wmemtool('wmb',win_loctool,n_miscella, ...
  311.                        ind_wave_fam,Wave_Fam,  ...
  312.                        ind_wave_nam,Wave_Nam,  ...
  313.                        ind_refinment,prec_val  ...
  314.                        );
  316.         % General graphical parameters initialization.
  317.         %--------------------------------------------
  318.         pos_g_area     = wmemtool('rmb',win_loctool,n_miscella,ind_graph_area) ;
  319.         pos_hdl_frame0 = get(hdl_frame0,'Position');
  320.         win_units      = 'normalized';
  321.         pos_win        = get(win_loctool,'Position');
  322.         bdx            = 0.08*pos_win(3);
  323.         bdy            = 0.09*pos_win(4);
  324.         bdy_d          = bdy;
  325.         bdy_u          = bdy;
  326.         bdyl           = pos_g_area(2)+bdy_d;
  327.         w_graph        = pos_hdl_frame0(1);
  328.         h_graph        = pos_g_area(4);
  330.         % Computing and displaying wavelets and filters.
  331.         %-----------------------------------------------
  332.         wtype = wavemngr('type',Wave_Nam);
  333.         commonAxesProp = {...
  334.            'Parent',win_loctool, ...
  335.            'Units',win_units,    ...
  336.            'DrawMode','Fast',    ...
  337.            'Box','On'            ...
  338.            };
  339.         stemCOL = wtbutils('colors','wfilters');
  340.         switch wtype
  341.             case 1
  342.                 str_wintitle = [Wave_Fam ' Wavelet --> ' Wave_Nam];
  343.                 [phi,psi,xVal] = wavefun(Wave_Nam,prec_val);
  344.                 [Lo_D,Hi_D,Lo_R,Hi_R] = wfilters(Wave_Nam);
  345.                 xaxis   = [min(xVal)  max(xVal)];
  346.                 xaxis_f = [0          length(Lo_D)-1];
  347.                 xVal_f  = xaxis_f(1):xaxis_f(2);
  348.                 nb      = 1;
  349.                 len     = length(xVal_f);
  350.                 while len>10
  351.                     nb = 2*nb;
  352.                     len = len/nb;
  353.                 end
  354.                 tics    = 0:nb:length(xVal_f);
  355.                 xlabs   = int2str(tics');
  356.                 yaxis   = [min(min(phi),min(psi))     max(max(phi),max(psi))];
  357.                 ecart   = abs(yaxis(2)-yaxis(1))/20;
  358.                 yaxis(1) = yaxis(1)-ecart;
  359.                 yaxis(2) = yaxis(2)+ecart;
  360.                 yaxis_r = [min(min(Lo_R),min(Hi_R))   max(max(Lo_R),max(Hi_R))];
  361.                 yaxis_d = [min(min(Lo_D),min(Hi_D))   max(max(Lo_D),max(Hi_D))];
  362.                 h_axe_wave = (h_graph-2*bdy-bdy_d-bdy_u)/2;
  363.                 h_axe_filt = h_axe_wave/2;
  364.                 w_axe      = (w_graph-3*bdx)/2;
  365.                 pos_phi  = [ bdx,...
  366.                              bdyl+2*bdy+2*h_axe_filt,...
  367.                              w_axe,h_axe_wave ];
  368.                 pos_psi  = [ 2*bdx+w_axe,...
  369.                              bdyl+2*bdy+2*h_axe_filt,...
  370.                              w_axe,h_axe_wave ];
  371.                 pos_Lo_D = [ bdx,...
  372.                              bdyl+bdy+h_axe_filt,...
  373.                              w_axe,h_axe_filt ];
  374.                 pos_Lo_R = [ bdx,...
  375.                               bdyl,w_axe,h_axe_filt ];
  376.                 pos_Hi_D = [ 2*bdx+w_axe,...
  377.                              bdyl+bdy+h_axe_filt,...
  378.                              w_axe,h_axe_filt ];
  379.                 pos_Hi_R = [ 2*bdx+w_axe,bdyl,w_axe,h_axe_filt ];
  381.                 axeProp  = {commonAxesProp{:},'Position',pos_phi};
  382.                 axeTitle = 'Scaling function phi';
  383.                 axe_phi  = plotYval(xVal,phi,axeProp,'r',axeTitle);
  385.                 axeProp  = {commonAxesProp{:},'Position',pos_psi};
  386.                 axeTitle = 'Wavelet function psi';
  387.                 axe_psi  = plotYval(xVal,psi,axeProp,'g',axeTitle);
  389.                 axe_Lo_D = axes(commonAxesProp{:},'Position',pos_Lo_D);
  390.                 wdstem(axe_Lo_D,xVal_f,Lo_D,stemCOL,1);
  391.                 set(axe_Lo_D,...
  392.                     'Ylim',yaxis_d,            ...
  393.                     'Xlim',xaxis_f,            ...
  394.                     'Box','On',                ...
  395.                     'XTicklabelMode','manual', ...
  396.                     'XTick',tics,              ...
  397.                     'XTickLabel',xlabs         ...
  398.                     );
  399.                 wtitle('Decomposition low-pass filter','Parent',axe_Lo_D);
  401.                 axe_Lo_R = axes(commonAxesProp{:},'Position',pos_Lo_R);
  402.                 wdstem(axe_Lo_R,xVal_f,Lo_R,stemCOL,1);
  403.                 set(axe_Lo_R,...
  404.                     'Ylim',yaxis_r,            ...
  405.                     'Xlim',xaxis_f,            ...
  406.                     'Box','On',                ...
  407.                     'XTicklabelMode','manual', ...
  408.                     'XTick',tics,              ...
  409.                     'XTickLabel',xlabs         ...
  410.                     );
  411.                 wtitle('Reconstruction low-pass filter','Parent',axe_Lo_R);
  413.                 axe_Hi_D = axes(commonAxesProp{:},'Position',pos_Hi_D);
  414.                 wdstem(axe_Hi_D,xVal_f,Hi_D,stemCOL,1);
  415.                 set(axe_Hi_D,...
  416.                     'Ylim',yaxis_d,            ...
  417.                     'Xlim',xaxis_f,            ...
  418.                     'Box','On',                ...
  419.                     'XTicklabelMode','manual', ...
  420.                     'XTick',tics,              ...
  421.                     'XTickLabel',xlabs         ...
  422.                     );
  423.                 wtitle('Decomposition high-pass filter','Parent',axe_Hi_D);
  425.                 axe_Hi_R = axes(commonAxesProp{:},'Position',pos_Hi_R);
  426.                 wdstem(axe_Hi_R,xVal_f,Hi_R,stemCOL,1);
  427.                 set(axe_Hi_R,...
  428.                     'Ylim',yaxis_r,            ...
  429.                     'Xlim',xaxis_f,            ...
  430.                     'Box','On',                ...
  431.                     'XTicklabelMode','manual', ...
  432.                     'XTick',tics,              ...
  433.                     'XTickLabel',xlabs         ...
  434.                     );
  435.                 wtitle('Reconstruction high-pass filter','Parent',axe_Hi_R);
  436.                 axe_cmd = [axe_phi axe_psi];
  437.                 axe_act = [axe_Lo_D axe_Hi_D axe_Lo_R axe_Hi_R];
  439.             case 2
  440.                 [phi1,psi1,phi2,psi2,xVal] = wavefun(Wave_Nam,prec_val);
  441.                 [Lo_D,Hi_D,Lo_R,Hi_R]  = wfilters(Wave_Nam);
  442.                 xaxis   = [min(xVal)    max(xVal)];
  443.                 xaxis_f = [0            length(Lo_D)-1];
  444.                 xVal_f  = 0:xaxis_f(2);
  445.                 nb      = 1;
  446.                 len     = length(xVal_f);
  447.                 while len>10
  448.                     nb = 2*nb;
  449.                     len = len/nb;
  450.                 end
  451.                 tics    = 0:nb:length(xVal_f);
  452.                 xlabs   = int2str(tics');
  453.                 yaxis1  = [min(min(phi1),min(psi1))   max(max(phi1),max(psi1))];
  454.                 ecart   = abs(yaxis1(2)-yaxis1(1))/20;
  455.                 yaxis1(1) = yaxis1(1)-ecart;
  456.                 yaxis1(2) = yaxis1(2)+ecart;
  457.                 yaxis2  = [min(min(phi2),min(psi2))   max(max(phi2),max(psi2))];
  458.                 ecart   = abs(yaxis2(2)-yaxis2(1))/20;
  459.                 yaxis2(1) = yaxis2(1)-ecart;
  460.                 yaxis2(2) = yaxis2(2)+ecart;
  461.                 yaxis_r = [min(min(Lo_R),min(Hi_R)) max(max(Lo_R),max(Hi_R))];
  462.                 yaxis_d = [min(min(Lo_D),min(Hi_D)) max(max(Lo_D),max(Hi_D))];
  463.                 str_wintitle = ['Biorthogonal Wavelet --> ' Wave_Nam];
  464.                 h_axe_wave = 2*(h_graph-3*bdy-bdy_d-bdy_u)/6;
  465.                 h_axe_filt = (h_graph-3*bdy-bdy_d-bdy_u)/6;
  466.                 w_axe      = (w_graph-3*bdx)/2;
  468.                 pos_phi1 = [ bdx,...
  469.                              bdyl+3*bdy+2*h_axe_filt+h_axe_wave,...
  470.                              w_axe,h_axe_wave];
  471.                 pos_phi2 = [ bdx,...
  472.                              bdyl+bdy+h_axe_filt,...
  473.                              w_axe,h_axe_wave];
  474.                 pos_psi1 = [ 2*bdx+w_axe,...
  475.                              bdyl+3*bdy+2*h_axe_filt+h_axe_wave,...
  476.                              w_axe,h_axe_wave];
  477.                 pos_psi2 = [ 2*bdx+w_axe,...
  478.                              bdyl+bdy+h_axe_filt,...
  479.                              w_axe,h_axe_wave];
  480.                 pos_Lo_D = [ bdx,...
  481.                              bdyl+2*bdy+h_axe_filt+h_axe_wave,...
  482.                              w_axe,h_axe_filt];
  483.                 pos_Lo_R = [ bdx,bdyl,w_axe,h_axe_filt];
  484.                 pos_Hi_D = [ 2*bdx+w_axe,...
  485.                               bdyl+2*bdy+h_axe_filt+h_axe_wave,...
  486.                               w_axe,h_axe_filt];
  487.                 pos_Hi_R = [ 2*bdx+w_axe,bdyl,w_axe,h_axe_filt];
  489.                 axeProp  = {commonAxesProp{:},'Position',pos_phi1};
  490.                 axeTitle = 'Decomposition scaling function phi';
  491.                 axe_phi1 = plotYval(xVal,phi1,axeProp,'r',axeTitle);
  493.                 axeProp  = {commonAxesProp{:},'Position',pos_phi2};
  494.                 axeTitle = 'Reconstruction scaling function phi';
  495.                 axe_phi2 = plotYval(xVal,phi2,axeProp,'r',axeTitle);
  497.                 axeProp  = {commonAxesProp{:},'Position',pos_psi1};
  498.                 axeTitle = 'Decomposition wavelet function psi';
  499.                 axe_psi1 = plotYval(xVal,psi1,axeProp,'g',axeTitle);
  501.                 axeProp  = {commonAxesProp{:},'Position',pos_psi2};
  502.                 axeTitle = 'Reconstruction wavelet function psi';
  503.                 axe_psi2 = plotYval(xVal,psi2,axeProp,'g',axeTitle);
  505.                 axe_Lo_D = axes(commonAxesProp{:},'Position',pos_Lo_D);
  506.                 wdstem(axe_Lo_D,xVal_f,Lo_D,stemCOL,1);
  507.                 set(axe_Lo_D,...
  508.                     'Ylim',yaxis_d,...
  509.                     'Xlim',xaxis_f,...
  510.                     'Box','On',...
  511.                     'XTicklabelMode','manual',...
  512.                     'XTick',tics,...
  513.                     'XTickLabel',xlabs...
  514.                     );
  515.                 wtitle('Decomposition low-pass filter','Parent',axe_Lo_D);
  517.                 axe_Lo_R = axes(commonAxesProp{:},'Position',pos_Lo_R);
  518.                 wdstem(axe_Lo_R,xVal_f,Lo_R,stemCOL,1);
  519.                 set(axe_Lo_R,...
  520.                     'Ylim',yaxis_r,...
  521.                     'Xlim',xaxis_f,...
  522.                     'Box','On',...
  523.                     'XTicklabelMode','manual',...
  524.                     'XTick',tics,...
  525.                     'XTickLabel',xlabs...
  526.                     );
  527.                 wtitle('Reconstruction low-pass filter','Parent',axe_Lo_R);
  529.                 axe_Hi_D = axes(commonAxesProp{:},'Position',pos_Hi_D);
  530.                 wdstem(axe_Hi_D,xVal_f,Hi_D,stemCOL,1);
  531.                 set(axe_Hi_D,...
  532.                     'Ylim',yaxis_d,...
  533.                     'Xlim',xaxis_f,...
  534.                     'Box','On',...
  535.                     'XTicklabelMode','manual',...
  536.                     'XTick',tics,...
  537.                     'XTickLabel',xlabs...
  538.                     );
  539.                 wtitle('Decomposition high-pass filter','Parent',axe_Hi_D);
  541.                 axe_Hi_R = axes(commonAxesProp{:},'Position',pos_Hi_R);
  542.                 wdstem(axe_Hi_R,xVal_f,Hi_R,stemCOL,1);
  543.                 set(axe_Hi_R,...
  544.                     'Ylim',yaxis_r,...
  545.                     'Xlim',xaxis_f,...
  546.                     'Box','On',...
  547.                     'XTicklabelMode','manual',...
  548.                     'XTick',tics,...
  549.                     'XTickLabel',xlabs...
  550.                     );
  551.                 wtitle('Reconstruction high-pass filter','Parent',axe_Hi_R);
  553.                 axe_cmd = [axe_phi1 axe_phi2 axe_psi1 axe_psi2];
  554.                 axe_act = [axe_Lo_D axe_Hi_D axe_Lo_R axe_Hi_R];
  556.             case 3
  557.                 str_wintitle   = [Wave_Fam '  Wavelet --> ' Wave_Nam];
  558.                 [phi,psi,xVal] = wavefun(Wave_Nam,prec_val);
  559.                 h_axe   = h_graph-bdy_d-bdy_u;
  560.                 w_axe   = (w_graph-3*bdx)/2;
  561.                 pos_phi = [bdx,bdyl,w_axe,h_axe];
  562.                 pos_psi = [2*bdx+w_axe,bdyl,w_axe,h_axe];
  564.                 axeProp  = {commonAxesProp{:},'Position',pos_phi};
  565.                 axeTitle = 'Scaling function phi';
  566.                 axe_phi  = plotYval(xVal,phi,axeProp,'r',axeTitle);
  568.                 axeProp  = {commonAxesProp{:},'Position',pos_psi};
  569.                 axeTitle = 'Wavelet function psi';
  570.                 axe_psi  = plotYval(xVal,psi,axeProp,'g',axeTitle);
  572.                 axe_cmd = [axe_phi axe_psi];
  573.                 axe_act = [];
  575.             case 4
  576.                 str_wintitle = [Wave_Fam '  Wavelet --> ' Wave_Nam];
  577.                 [psi,xVal]   = wavefun(Wave_Nam,prec_val);
  578.                 h_axe   = h_graph-bdy_d-bdy_u;
  579.                 w_axe   = w_graph-2*bdx;
  580.                 pos_psi = [bdx,bdyl,w_axe,h_axe];
  582.                 axeProp  = {commonAxesProp{:},'Position',pos_psi};
  583.                 axeTitle = 'Wavelet function psi';
  584.                 axe_psi  = plotYval(xVal,psi,axeProp,'g',axeTitle);
  586.                 axe_cmd = axe_psi;
  587.                 axe_act = [];
  589.             case 5
  590.                 str_wintitle = [Wave_Fam '  Wavelet --> ' Wave_Nam];
  591.                 [psi,xVal] = wavefun(Wave_Nam,prec_val);
  592.                 psiAbs = abs(psi);
  593.                 psiAng = angle(psi);
  594.                 h_axe  = (h_graph-bdy_d-bdy_u-bdy)/2;
  595.                 w_axe  = (w_graph-3*bdx)/2;
  597.                 pos_abs = [bdx,bdyl,w_axe,h_axe];
  598.                 pos_ang = [2*bdx+w_axe,bdyl,w_axe,h_axe];
  599.                 y_axe   = bdyl+h_axe+bdy;
  600.                 pos_rea = [bdx,y_axe,w_axe,h_axe];
  601.                 pos_ima = [2*bdx+w_axe,y_axe,w_axe,h_axe];
  603.                 yVal     = real(psi);
  604.                 axeProp  = {commonAxesProp{:},'Position',pos_rea};
  605.                 axeTitle = 'Real part of function psi';
  606.                 axe_rea  = plotYval(xVal,yVal,axeProp,'r',axeTitle);
  608.                 yVal     = imag(psi);
  609.                 axeProp  = {commonAxesProp{:},'Position',pos_ima};
  610.                 axeTitle = 'Imaginary part of function psi';
  611.                 axe_ima  = plotYval(xVal,yVal,axeProp,'g',axeTitle);
  613.                 yVal     = abs(psi);
  614.                 axeProp  = {commonAxesProp{:},'Position',pos_abs};
  615.                 axeTitle = 'Modulus of function psi';
  616.                 axe_abs  = plotYval(xVal,yVal,axeProp,'r',axeTitle);
  618.                 yVal     = angle(psi);
  619.                 axeProp  = {commonAxesProp{:},'Position',pos_ang};
  620.                 axeTitle = 'Angle of function psi';
  621.                 axe_ang  = plotYval(xVal,yVal,axeProp,'g',axeTitle);
  622.  
  623.                 axe_cmd = [axe_rea axe_ima axe_abs axe_ang];
  624.                 axe_act = [];
  625.         end
  627.         % Display status line.
  628.         %---------------------
  629.         wfigtitl('string',win_loctool,str_wintitle,'on');
  631.         % Axes attachment.
  632.         %-----------------
  633.         dynvtool('init',win_loctool,[],axe_cmd,axe_act,[0 0]);
  635.         % Setting units to normalized.
  636.         %-----------------------------
  637.         set(findobj(win_loctool,'Units','pixels'),'Units','normalized');
  639.         % End waiting.
  640.         %-------------
  641.         wwaiting('off',win_loctool);
  643.     case 'new'
  644.         %*************************************************%
  645.         %** OPTION = 'new' -  test drawing parameters   **%
  646.         %*************************************************%
  647.         win_loctool = varargin{1};
  649.         % Handles of tagged objects.
  650.         %---------------------------
  651.         pop_handles  = findobj(win_loctool,'Style','popupmenu');
  652.         pop_prec_val = findobj(pop_handles,'Tag',tag_prec_val);
  654.         % Test Main parameters selection before drawing.
  655.         %-----------------------------------------------
  656.         Wave_Nam = cbanapar('get',win_loctool,'wav');
  657.         Wave_Fam = wavemngr('fam_num',Wave_Nam);
  658.         prec_val = get(pop_prec_val,'Value')+4;
  659.         [wfam,wnam,raf] = wmemtool('rmb',win_loctool,n_miscella,...
  660.                                 ind_wave_fam,ind_wave_nam,ind_refinment);
  661.         if raf~=prec_val | ~strcmp(wnam,Wave_Nam) | ~strcmp(wfam,Wave_Fam)
  662.             varargout = {1,Wave_Fam,Wave_Nam,prec_val};
  663.         else
  664.             varargout = {0,Wave_Fam,Wave_Nam,prec_val};
  665.             return
  666.         end
  668.         % Setting refinment to 0 (as flag).
  669.         %----------------------------------
  670.         wmemtool('wmb',win_loctool,n_miscella,ind_refinment,0);
  672.         % Cleaning the graphical part.
  673.         %-----------------------------
  674.         dynvtool('stop',win_loctool);
  675.         axe_handles = findobj(get(win_loctool,'Children'),'flat', ...
  676.                               'type','axes','visible','on');
  677.         delete(axe_handles);
  678.         wfigtitl('vis',win_loctool,'off');
  680.     case 'demo'
  681.         %*******************************************%
  682.         %** OPTION = 'demo' -  for DEMOS or TESTS **%
  683.         %*******************************************%
  684.         win_loctool = varargin{1};
  685.         Wave_Nam    = varargin{2};
  687.         % Handles of tagged objects.
  688.         %---------------------------
  689.         children     = get(win_loctool,'Children');
  690.         uic_handles  = findobj(children,'flat','type','uicontrol');
  691.         pop_handles  = findobj(uic_handles,'Style','popupmenu');
  692.         pop_prec_val = findobj(pop_handles,'Tag',tag_prec_val);
  693.         pus_handles  = findobj(uic_handles,'Style','pushbutton');
  694.         pus_inf1     = findobj(pus_handles,'Tag',tag_pus_inf1);
  696.         cbanapar('set',win_loctool,'wav',Wave_Nam);
  697.         [Wave_Fam,tabNums] = wavemngr('fields',Wave_Nam,'fsn','tabNums');
  699.         if nargin==4
  700.             set(pop_prec_val,'Value',varargin{3});
  701.         end
  702.         new = wvdtool('new',win_loctool);
  703.         str_inf1 = [upper(Wave_Fam) ' wavelet'] ;
  704.         if size(tabNums,1)>1
  705.             str_inf1 = [str_inf1 's'] ;
  706.         end
  707.         set(pus_inf1,'String',str_inf1);
  708.         wvdtool('draw1d',win_loctool);
  710.     case 'close'
  712.     otherwise
  713.         errargt(mfilename,'Unknown Option','msg');
  714.         error('*');
  715. end
  717. %=============================================================================%
  718. % INTERNAL FUNCTIONS
  719. %=============================================================================%
  720. %-----------------------------------------------------------------------------%
  721. function axe = plotYval(xVal,yVal,axeProp,color,axeTitle)
  723. xlim  = [min(xVal) , max(xVal)];
  724. mini  = min(yVal);
  725. maxi  = max(yVal);
  726. dyVal = maxi-mini;
  727. if dyVal<sqrt(eps) , dyVal = sqrt(eps);end
  728. ylim = [mini maxi]+0.02*dyVal*[-1 1];
  729. axe  = axes(axeProp{:},'Xlim',xlim,'Ylim',ylim);
  730. line('Xdata',xVal,'Ydata',yVal,'Color',color,'Parent',axe);
  731. wtitle(axeTitle,'Parent',axe);
  732. %-----------------------------------------------------------------------------%
  733. %=============================================================================%