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

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

波变换

开发平台：

Matlab

wavemngr.m：源码内容
  1. function varargout = wavemngr(option,varargin)
  2. %WAVEMNGR Wavelet manager.
  3. %   WAVEMNGR is a wavelet manager used to add, delete, restore
  4. %   or read wavelets.
  5. %
  6. %   WAVEMNGR('add',FN,FSN,WT,NUMS,FILE) or
  7. %   WAVEMNGR('add',FN,FSN,WT,NUMS,FILE,B) or 
  8. %   WAVEMNGR('add',FN,FSN,WT,{NUMS,TYPNUMS},FILE) or
  9. %   WAVEMNGR('add',FN,FSN,WT,{NUMS,TYPNUMS},FILE,B)
  10. %   adds a new wavelet family.
  11. %     FN  = family name (string).
  12. %     FSN = family short name (string).
  13. %
  14. %     WT defines the wavelet type: 
  15. %     WT = 1 for orthogonal wavelets.
  16. %     WT = 2 for biorthogonal wavelets.
  17. %     WT = 3 for wavelet with scale function.
  18. %     WT = 4 for wavelet without scale function.
  19. %     WT = 5 for complex wavelet without scale function.
  20. %
  21. %     If the wavelet is a single one, NUMS = ''.
  22. %       examples: mexh, morl.
  23. %     If the wavelet is part of a finite family of wavelets, NUMS
  24. %       is a string containing a blank separated list of items
  25. %       representing wavelet parameters.
  26. %       example: bior, NUMS = '1.1 1.3 ... 4.4 5.5 6.8'.
  27. %     If the wavelet is part of an infinite family of wavelets, 
  28. %       NUMS is a string containing a blank separated list of 
  29. %       items representing wavelet parameters, terminated by the 
  30. %       special sequence **.
  31. %       examples: 
  32. %         db,    NUMS = '1 2 3 4 5 6 7 8 9 10 **'.
  33. %         shan,  NUMS = '1-1.5 1-1 1-0.5 1-0.1 2-3 **'
  34. %     In these last two cases, TYPNUMS specifies the wavelet parameter 
  35. %       input format: 'integer' or 'real' or 'string'; the default 
  36. %       value is 'integer'.
  37. %       examples: db,   TYPNUMS = 'integer'
  38. %                 bior, TYPNUMS = 'real'
  39. %                 shan, TYPNUMS = 'string'
  40. %
  41. %     FILE  = MAT-file or M-file name (string).
  42. %
  43. %     B = [lb ub] specifies lower and upper bounds of
  44. %     effective support for wavelets of type = 3, 4 or 5.
  45. %
  46. %   WAVEMNGR('del',N), deletes a wavelet or a wavelet family where
  47. %     N is the wavelet name or the family short name.
  48. %
  49. %   WAVEMNGR('restore') restores the previous 
  50. %     wavelets.asc ASCII-file.
  51. %   WAVEMNGR('restore',IN2) restores the initial 
  52. %     wavelets.asc ASCII-file.
  53. %
  54. %   OUT1 = WAVEMNGR('read') returns all wavelets family names.
  55. % 
  56. %   OUT1 = WAVEMNGR('read',IN2) returns all wavelet names.
  57. %
  58. %   OUT1 = WAVEMNGR('read_asc') returns all wavelets information
  59. %   retrieved from wavelets.asc ASCII-file.
  61. %----------------------
  62. %   INTERNAL OPTIONS.
  63. %----------------------
  64. %   WAVEMNGR('create')
  65. %   creates wavelets.inf MAT-file using wavelets.asc ASCII-file.
  66. %
  67. %   OUT1 = WAVEMNGR('load') or WAVEMNGR('load')
  68. %   loads Wavelets_Info from wavelets.inf matfile,
  69. %   and puts it in the global variable: Wavelets_Info.
  70. %
  71. %   WAVEMNGR('clear') clear the global Wavelets_Info.
  72. %
  73. %   [OUT1,OUT2,OUT3,OUT4,OUT5,OUT6] = WAVEMNGR('indw',W)
  74. %   returns:
  75. %   family indice, number indice,
  76. %   family, number for the wavelet W.
  77. %   OUT5 is the family table of number.
  78. %   OUT6 is a flag for GUI.
  79. %
  80. %   OUT1 = WAVEMNGR('indf',F) returns indice
  81. %   for wavelet family F (short name).
  82. %
  83. %   varargout = wavemngr('fields',varargin) see below.
  84. %
  85. %   OUT1 = WAVEMNGR('read_struct')
  86. %   OUT1 gives all wavelets structures.
  87. % 
  88. %   OUT1 = WAVEMNGR('read_struct',NUMS)
  89. %   OUT1 gives the wavelets structures.
  90. %   specified by NUMS.
  91. %
  92. %   OUT1 = WAVEMNGR('tfsn') or WAVEMNGR('tfsn',T)
  93. %   returns shortname-table (of type = T : 'dwt' or 'cwt' or 'owt' or 'ccwt').
  94. %
  95. %   OUT1 = WAVEMNGR('tfn') or out1 = WAVEMNGR('tfn',T)
  96. %   returns name-table (of type = T : 'dwt' or 'cwt' or 'owt' or 'ccwt').
  97. %
  98. %   OUT1 = WAVEMNGR('isbior',W)
  99. %   returns 1 for biorthogonal wavelets.
  100. %
  101. %   [FAM,NUM] = WAVEMNGR('fam_num',W) returns
  102. %   family and number for the wavelet W.
  103. %
  104. %   LEN = WAVEMNGR('length',W) returns the
  105. %   length of the support of the wavelet W.
  107. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 12-Mar-96.
  108. %   Last Revision: 26-May-2003.
  109. %   Copyright 1995-2004 The MathWorks, Inc.
  110. % $Revision: 1.17.4.2 $
  112. global Wavelets_Info
  114. %----------------------%
  115. % Wavelets Structures. %
  116. %--------------------------------------------------------%
  117. % Wavelets_Info is a  structure array with size [nb_fam 1]
  118. %
  119. % Wavelet_Struct =
  120. %   struct(...
  121. %     'index'           integer
  122. %     'familyName'      string
  123. %     'familyShortName' string
  124. %     'type'            integer
  125. %     'tabNums'         matrix of string
  126. %     'typNums'         string
  127. %     'file'            string
  128. %     'bounds'          string
  129. %     );
  130. %---------------------------------------------------------%
  132. % Files Names.
  133. %--------------
  134. bin_ini_file = 'wavelets.bin';
  135. asc_ini_file = 'wavelets.ini';
  136. bin_file = 'wavelets.inf';
  137. asc_file = 'wavelets.asc';
  138. sav_file = 'wavelets.prv';
  140. if nargin==0 , option = 'create' ; end
  142. % Miscellaneous Values.
  143. %----------------------
  144. nb_InfoByWave =  7;
  145. NB_FamInWTB   = 11;
  146. WaveTypePOS   = [1:5];
  148. switch option
  149.     case 'load'
  150.       if isempty(Wavelets_Info)
  151.           try
  152.               load(bin_file,'-mat')
  153.           catch
  154.               try
  155.                   load(bin_ini_file,'-mat')
  156.               catch
  157.                   clc
  158.                   disp(' ');
  159.                   disp('---------------------------------------------');
  160.                   sprintf('*** File : %s not found ! ***', bin_ini_file);
  161.                   disp(['***     Using rescue resource ... !     ***']);
  162.                   disp('---------------------------------------------');
  163.                   disp(' ');
  164.                   Wavelets_Info  = wavemngr('rescue');
  165.               end
  166.           end
  167.       end
  168.       if nargout>0 , varargout{1} = Wavelets_Info; end
  170.     case 'clear' , clear global Wavelets_Info 
  172.     case 'indw'
  173.       % in2 : wavelet name
  174.       %-------------------
  175.       % out1 = i_fam
  176.       % out2 = i_num
  177.       % out3 = fam
  178.       % out4 = num_str
  179.       % out5 = tabNums
  180.       % out6 = flag '**' (for GUI)
  181.       %---------------------------
  182.       if isempty(Wavelets_Info) , wavemngr('load'); end
  183.       nb_fam = size(Wavelets_Info,1);
  184.       wname  = deblankl(varargin{1});
  185.       lwna   = length(wname);
  186.       for i_fam=1:nb_fam
  187.           fam = Wavelets_Info(i_fam).familyShortName;
  188.           len = length(fam);
  189.           ok_wave = 0;
  190.           if lwna>=len
  191.               if fam==wname(1:len)
  192.                   tabNums = Wavelets_Info(i_fam).tabNums;
  193.                   for i_num = 1:size(tabNums,1)
  194.                       num_str = noblank(tabNums(i_num,:));
  195.                       if strcmp(num_str,'no') , num_str = '' ; end
  196.                       if strcmp([fam num_str],wname)
  197.                           ok_wave = 1; add_num = 0; break;
  198.                       end
  199.                   end
  201.                   % test for ** number
  202.                   %------------------
  203.                   if ok_wave==0 && strcmp(num_str,'**') && (lwna>len)
  204.                       typNums = Wavelets_Info(i_fam).typNums;
  205.                       num_str = wname(len+1:lwna);
  206.                       switch typNums
  207.                         case 'integer'
  208.                           num = wstr2num(num_str);
  209.                           if ~isempty(num) && (num==fix(num)) && (0<num)
  210.                               ok_wave = 1; add_num = 1;
  211.                           end
  213.                         case 'real'
  214.                           num = wstr2num(num_str);
  215.                           if ~isempty(num)
  216.                               ok_wave = 1; add_num = 1;
  217.                           end
  219.                         case 'string'
  220.                           ok_wave = 1; add_num = 1;
  221.                       end                        
  222.                   end
  223.  
  224.               end
  225.           end
  226.           if ok_wave , break; end
  227.       end
  228.       if ok_wave
  229.           varargout = {i_fam,i_num};
  230.           if nargout<3 , return; end
  231.           varargout = {varargout{:},fam,num_str,tabNums,add_num};
  232.       else
  233.           msg = sprintf('Invalid wavelet name : %s', wname);
  234.           errargt('Wavelet test',msg,'msg');
  235.           error(['***  ' msg '  ***']);
  236.       end
  238.     case 'indf'
  239.       if isempty(Wavelets_Info) , wavemngr('load'); end
  240.       fsnCell  = {Wavelets_Info(:).familyShortName}';
  241.       fsnInput = deblankl(varargin{1});
  242.       varargout{1} = find(strcmp(fsnInput,fsnCell)==1);
  244.     case 'tfsn'
  245.       if isempty(Wavelets_Info) , wavemngr('load'); end
  246.       fsnCell = {Wavelets_Info(:).familyShortName}';
  247.       varargout{1} = strvcat(fsnCell{:});
  248.       if nargin==2
  249.           wtype = lower(varargin{1});
  250.           if isequal(wtype,'all') , return; end
  251.           tab_type = {Wavelets_Info(:).type}';
  252.           tab_type = cat(1,tab_type{:});
  253.           switch wtype
  254.             case {'dwt'}  , ind = find(tab_type==1 | tab_type==2);
  255.             case {'cwt'}  , ind = find(tab_type>0 & tab_type<5);
  256.             case {'owt'}  , ind = find(tab_type==1);
  257.             case {'ccwt'} , ind = find(tab_type==5);
  258.             case {'666'}
  259.                ind = [find(tab_type==5) ; find(tab_type>0 & tab_type<5)];
  260.           end
  261.           varargout{1} = varargout{1}(ind,:);
  262.        end
  264.     case 'tfn'
  265.       if isempty(Wavelets_Info) , wavemngr('load'); end
  266.       fnCell = {Wavelets_Info(:).familyName}';
  267.       varargout{1} = strvcat(fnCell{:});
  268.       if nargin==2
  269.           wtype = lower(varargin{1});
  270.           if isequal(wtype,'all') , return; end
  271.           tab_type = {Wavelets_Info(:).type}';
  272.           tab_type = cat(1,tab_type{:});
  273.           switch wtype
  274.             case {'dwt'}  , ind = find(tab_type==1 | tab_type==2);
  275.             case {'cwt'}  , ind = find(tab_type>0 & tab_type<5);
  276.             case {'owt'}  , ind = find(tab_type==1);
  277.             case {'ccwt'} , ind = find(tab_type==5);
  278.           end
  279.           varargout{1} = varargout{1}(ind,:);
  280.        end
  282.     case 'fields'
  283.       % in2 = {'ind', wavelet index} or
  284.       % in2 = {'fsn', family shortname} or
  285.       % in2 = {'wn' , wavelet name}
  286.       % or
  287.       % in2 = wavelet name
  288.       %----------------------------
  289.       % in3  ... = field(s) name(s)
  290.       %----------------------------
  291.       % out1 ... = field(s) value(s)
  292.       %-----------------------------
  293.       if iscell(varargin{1})
  294.           in_type = varargin{1}{1};
  295.           arg = varargin{1}{2};
  296.       else
  297.           in_type = 'wn';
  298.           arg = varargin{1};
  299.       end
  300.       switch in_type
  301.           case 'ind'
  302.               if isempty(Wavelets_Info) , wavemngr('load'); end
  303.               i_fam = arg;
  304.  
  305.           case 'fsn' , i_fam = wavemngr('indf',arg);
  306.           case 'wn'  , i_fam = wavemngr('indw',arg);
  307.       end
  308.       nb = nargin-2;
  309.       if nb==0
  310.           varargout{1} = Wavelets_Info(i_fam);
  311.           return;
  312.       end
  313.       for k=1:nb
  314.           switch varargin{k+1}
  315.               case 'ind'     , field = 'index';
  316.               case 'fn'      , field = 'familyName';
  317.               case 'fsn'     , field = 'familyShortName';
  318.               case 'type'    , field = 'type';
  319.               case 'tabNums' , field = 'tabNums';
  320.               case 'typNums' , field = 'typNums';
  321.               case 'file'    , field = 'file';
  322.               case 'bounds'  , field = 'bounds';
  323.               otherwise      , field = '';
  324.           end
  325.           if ~isempty(field)
  326.               varargout{k} = Wavelets_Info(i_fam).(field);
  327.           else
  328.               varargout{k} = Wavelets_Info(i_fam);
  329.           end
  330.       end
  332.     case 'type'
  333.         i_fam = wavemngr('indw',varargin{1});
  334.         varargout{1} = Wavelets_Info(i_fam).('type');
  336.     case 'file'
  337.         i_fam = wavemngr('indw',varargin{1});
  338.         varargout{1} = Wavelets_Info(i_fam).('file');
  340.     case 'fn'
  341.         i_fam = wavemngr('indw',varargin{1});
  342.         varargout{1} = Wavelets_Info(i_fam).('familyName');
  344.     case 'fsn'
  345.         i_fam = wavemngr('indw',varargin{1});
  346.         varargout{1} = Wavelets_Info(i_fam).('familyShortName');
  348.     case {'tabnums','tabNums'}
  349.         i_fam = wavemngr('indf',varargin{1});
  350.         varargout{1} = Wavelets_Info(i_fam).('tabNums');
  352.     case {'typNums'}
  353.         i_fam = wavemngr('indf',varargin{1});
  354.         varargout{1} = Wavelets_Info(i_fam).('typNums');
  356.     case 'bounds'
  357.         i_fam = wavemngr('indw',varargin{1});
  358.         varargout{1} = Wavelets_Info(i_fam).('bounds');
  360.     case 'isbior'
  361.         wname = varargin{1};
  362.         if length(wname)>3,
  363.             wname = wname(1:4);
  364.             varargout{1} = isequal(wname,'bior') | isequal(wname,'rbio');
  365.         else
  366.             varargout{1} = 0;
  367.         end
  369.     case 'fam_num'
  370.         [nul,nul,varargout{1},varargout{2}] = wavemngr('indw',varargin{1});
  372.     case 'length'
  373.         wname = varargin{1};
  374.         [wtype,bounds] = wavemngr('fields',{'wn',wname},'type','bounds');
  375.         switch wtype
  376.             case {1,2}
  377.                 Lo_D = wfilters(wname);
  378.                 varargout{1} = length(Lo_D);
  379.             case {3,4}
  380.                 varargout{1} = bounds(2)-bounds(1)+1;
  381.             otherwise
  382.                 errargt(mfilename,'invalid argument','msg'); error('*');
  383.         end
  384.                                        
  385.     case 'create'
  386.         wavemngr('clear')
  387.         beg_fam = 'Family Name : ';
  388.         fid = fopen(asc_file);
  389.         if fid==-1
  390.             fid   = fopen(asc_ini_file); 
  391.             winfo = fread(fid);
  392.             fclose(fid);
  393.             fid   = fopen(asc_file,'w');
  394.             fwrite(fid,winfo); 
  395.             fclose(fid);
  396.         else
  397.             winfo = fread(fid);
  398.             fclose(fid);
  399.         end
  400.         winfo   = setstr(winfo');
  401.         str_NL  = findNL('string',winfo);
  402.         lennewl = length(str_NL);
  403.         ind_NL  = findstr(str_NL,winfo);
  404.         first   = findstr(beg_fam,winfo)+length(beg_fam);
  405.         nb_fam  = length(first);
  407.         %------------------------------%
  408.         % Building Wavelets Structure. %
  409.         %------------------------------%
  410.         nul = cell(nb_fam,1);
  411.         Wavelets_Info = struct(...
  412.                     'index',nul,           ...
  413.                     'familyName',nul,      ...
  414.                     'familyShortName',nul, ...
  415.                     'type',nul,            ...
  416.                     'tabNums',nul,         ...
  417.                     'typNums',nul,         ...
  418.                     'file',nul,            ...
  419.                     'bounds',nul           ...
  420.                     );
  422.         for j = 1:nb_fam
  423.             i_fam   = first(j);
  424.             indexs  = find(ind_NL>i_fam);
  425.             indexs  = ind_NL(indexs(1:nb_InfoByWave));
  426.             fam     = winfo(i_fam:indexs(1)-1);
  427.             sname   = winfo(indexs(1)+lennewl:indexs(2)-1);
  428.             wtype   = winfo(indexs(2)+lennewl:indexs(3)-1);
  429.             nums    = winfo(indexs(3)+lennewl:indexs(4)-1);
  430.             typNums = winfo(indexs(4)+lennewl:indexs(5)-1);                  
  431.             fname   = winfo(indexs(5)+lennewl:indexs(6)-1);
  432.             bounds  = winfo(indexs(6)+lennewl:indexs(7)-1);
  434.             notspace = ~isspace(nums);
  435.             lnot     = length(notspace);
  436.             index1   = find(notspace==1);
  437.             k0       = index1(1);
  438.             k1       = index1(length(index1));
  439.             indnum   = diff(notspace);
  440.             fnum     = find(indnum==1)+1;
  441.             lnum     = find(indnum==-1);
  442.             if k0==1    , fnum = [1 fnum];  end
  443.             if k1==lnot , lnum = [lnum k1]; end
  444.             nb_num   = length(fnum);
  445.             tabNums  = '';
  446.             for p = 1:nb_num
  447.                 tabNums = strvcat(tabNums,nums(fnum(p):lnum(p)));
  448.             end
  450.             Wavelets_Info(j).index = j;
  451.             Wavelets_Info(j).familyName = fam;
  452.             Wavelets_Info(j).familyShortName = sname;
  453.             Wavelets_Info(j).type = wstr2num(wtype);
  454.             Wavelets_Info(j).tabNums = tabNums;
  455.             Wavelets_Info(j).typNums = typNums;
  456.             Wavelets_Info(j).file = fname;
  457.             Wavelets_Info(j).bounds = wstr2num(bounds);
  458.         end
  459.         try
  460.             save(bin_file,'Wavelets_Info')
  461.         catch
  462.             errargt(mfilename,'Changing Wavelets : Save FAILED !','msg');
  463.         end
  465.     case 'read'
  466.         if isempty(Wavelets_Info) , wavemngr('load'); end
  467.         nb_fam  = size(Wavelets_Info,1);
  468.         sep_fam = '=';
  469.         sep_fam = sep_fam(:,ones(1,35));
  470.         sep_num = '-';
  471.         sep_num = sep_num(:,ones(1,30));
  472.         varargout{1} = sep_fam; 
  473.         tab     = setstr(9);
  475.         families = strvcat(Wavelets_Info(:).familyName);
  476.         tabtab   = tab(ones(1,nb_fam),ones(1,2));
  478.         if  nargin==1
  479.             famSName    = strvcat(Wavelets_Info(:).familyShortName);
  480.             famTAB      = [families tabtab famSName];
  481.             varargout{1}= strvcat(varargout{1},famTAB,sep_fam);
  482.         else
  483.             famTAB = [families tabtab];
  484.             for k =1:nb_fam
  485.                  sfname = Wavelets_Info(k).familyShortName;
  486.                  varargout{1} = strvcat(varargout{1},[famTAB(k,:) sfname]);
  487.                  nb     = 0;
  488.                  wnames = [];
  489.                  tabNums = Wavelets_Info(k).tabNums;
  490.                  if size(tabNums,1)>1
  491.                      varargout{1} = strvcat(varargout{1},sep_num);
  492.                  end
  493.                  for j = 1:size(tabNums,1)
  494.                      num_str = noblank(tabNums(j,:));
  495.                      if ~strcmp(num_str,'no')
  496.                          wnames = [wnames sfname noblank(tabNums(j,:)) tab];
  497.                      end
  498.                      if nb<3
  499.                          nb = nb+1;
  500.                      else
  501.                          if ~isempty(wnames)
  502.                              varargout{1} = strvcat(varargout{1},wnames);
  503.                          end
  504.                          nb  = 0;
  505.                          wnames = [];
  506.                      end
  507.                  end
  508.                  if nb>0 && ~isempty(wnames)
  509.                      varargout{1} = strvcat(varargout{1},wnames);
  510.                  end
  511.                  varargout{1} = strvcat(varargout{1},sep_fam);
  512.             end
  513.         end
  515.     case 'read_asc'
  516.         fid = fopen(asc_file);
  517.         if fid==-1 , fid = fopen(asc_ini_file); end
  518.         winfo = fread(fid);
  519.         fclose(fid);
  520.         varargout{1} = setstr(winfo');
  522.     case 'read_struct'
  523.         if isempty(Wavelets_Info) , wavemngr('load'); end
  524.         if nargin==2
  525.             indfam = varargin{1};
  526.         else
  527.             nb_fam = size(Wavelets_Info,1);
  528.             indfam = [1:nb_fam];
  529.         end
  530.         if nargout==0
  531.             sep_fam = '*';
  532.             sep_fam = sep_fam(:,ones(1,35));
  533.             disp(' '); disp(' '); disp(' ');
  534.             disp(sep_fam);
  535.         end
  536.         for k =1:length(indfam)
  537.             tmp{k} = Wavelets_Info(indfam(k));
  538.             if nargout==0
  539.                 disp(tmp{k}); disp(sep_fam);
  540.             end
  541.         end
  542.         if nargout>0 , varargout = tmp; end
  544.     case 'add'
  545.         wavemngr('clear')
  546.         Wavelets_Info = wavemngr('load');
  547.         nb_fam = size(Wavelets_Info,1);
  548.         err = 0;
  549.         family_Name = varargin{1};        
  550.         if isempty(family_Name)
  551.             err = 1; 
  552.             msg = 'Wavelet Family Name is empty !';
  553.         else
  554.             fnIn  = noblank(family_Name);
  555.             tmpCell = {Wavelets_Info(:).familyName}';
  556.             ind   = find(strcmp(fnIn,tmpCell)==1);
  557.             err   = ~isempty(ind);
  558.             if err , msg = 'The Wavelet Family Name is already used !'; end
  559.         end
  561.         if err==0
  562.             family_Short_Name = varargin{2};
  563.             if isempty(family_Short_Name)
  564.                 err = 1; 
  565.                 msg = 'Wavelet Family Short Name is empty !';
  566.             else
  567.                 fsnIn   = deblankl(family_Short_Name);
  568.                 tmpCell = {Wavelets_Info(:).familyShortName}';
  569.                 ind     = find(strcmp(fsnIn,tmpCell)==1);
  570.                 err     = ~isempty(ind);
  571.                 if err
  572.                     msg = 'The Wavelet Family Short Name is already used !';
  573.                 end
  574.             end
  575.         end
  577.         if err==0
  578.             wavelet_Type = varargin{3};
  579.             if isempty(find(wavelet_Type==WaveTypePOS))
  580.                 err = 1;
  581.                 msg = 'Invalid Wavelet Type !';
  582.             end
  583.         end
  585.         if err==0            
  586.             if isempty(varargin{4})
  587.                 wavelet_tabNums = 'no';
  588.                 wavelet_typNums = 'no';
  589.             elseif ischar(varargin{4})
  590.                 wavelet_tabNums = deblank(varargin{4});
  591.                 wavelet_typNums = 'integer';
  592.             elseif iscell(varargin{4})
  593.                 wavelet_tabNums = deblank(varargin{4}{1});
  594.                 wavelet_typNums = deblank(varargin{4}{2});
  595.                 if ischar(wavelet_typNums)
  596.                    switch wavelet_typNums
  597.                      case {'integer','real','string'}
  598.                      otherwise
  599.                        err = 1;
  600.                        msg = 'Invalid Wavelet type of numbers !';
  601.                    end                 
  602.                 else
  603.                     err = 1;
  604.                     msg = 'Invalid Wavelet type of numbers !';
  605.                 end
  606.                 
  607.             else
  608.                 err = 1;
  609.                 msg = 'Invalid Wavelet numbers !';
  610.             end
  611.         end
  613.         if err==0
  614.             wavelet_File = varargin{5};
  615.             if isempty(wavelet_File)
  616.                 err = 1;
  617.             elseif findstr('.mat',wavelet_File)
  619.             else 
  620.                 wavelet_File = deblankl(wavelet_File);
  621.                 ind = findstr('.m',wavelet_File);
  622.                 if ind>0 , wavelet_File = wavelet_File(1:ind-1); end
  623.                 if isempty(wavelet_File) , err = 1; end
  624.             end
  625.             if err==1
  626.                 msg = 'Invalid Wavelet File Name !';
  627.             end
  628.         end
  629.         if err==0
  630.             nbArgIN = length(varargin);
  631.             switch wavelet_Type
  632.               case {1,2}
  633.                 if nbArgIN<6 , wavelet_Bounds = ''; end
  634.                 
  635.               otherwise  
  636.                 if nbArgIN<6
  637.                     err = 1;
  638.                     msg = 'Invalid number of arguments !';
  639.                 else
  640.                     wavelet_Bounds = varargin{6};
  641.                     if length(wavelet_Bounds)~=2,               err = 1;
  642.                     elseif wavelet_Bounds(1)>wavelet_Bounds(2), err = 1;
  643.                     end
  644.                     if err==1
  645.                         msg = 'Invalid value for wavelet bounds !';
  646.                     end
  647.                 end
  648.             end
  649.         end
  650.         if err
  651.             msg = strvcat('Add New Wavelet FAILED !!',msg);
  652.             errargt(mfilename,msg,'msg');
  653.             return
  654.         end
  656.         fid = fopen(asc_file);
  657.         if fid==-1 , fid = fopen(asc_ini_file); end
  658.         winfo = fread(fid);
  659.         fclose(fid);
  660.         fid = fopen(sav_file,'w');
  661.         fwrite(fid,winfo);
  662.         fclose(fid);
  664.         Chrline = findNL('char',winfo);
  665.         beg_fam = 'Family Name : ';
  666.         sep_fam = '------------------------';
  668.         wavelet_Type = sprintf('%.0f',wavelet_Type);
  669.         if ~isempty(wavelet_Bounds)
  670.             wavelet_Bounds = [num2str(wavelet_Bounds(1)) ' ' ...
  671.                               num2str(wavelet_Bounds(2))];
  672.         end
  674.         winfo = [winfo(1:end-1);           Chrline;  ...
  675.                 abs(beg_fam'); ...
  676.                 abs(family_Name(:));       Chrline;  ...
  677.                 abs(family_Short_Name(:)); Chrline;  ...
  678.                 abs(wavelet_Type(:));      Chrline;  ...
  679.                 abs(wavelet_tabNums(:));   Chrline;  ...
  680.                 abs(wavelet_typNums(:));   Chrline;  ...
  681.                 abs(wavelet_File(:));      Chrline;  ...
  682.                 abs(wavelet_Bounds(:));    Chrline;  ...
  683.                 abs(sep_fam');             Chrline   ...
  684.                 ];
  686.         fid = fopen(asc_file,'w');
  687.         fwrite(fid,winfo);
  688.         fclose(fid);
  689.         wavemngr('create');
  691.     case 'del'
  692.         wavemngr('clear')
  693.         Wavelets_Info = wavemngr('load');
  694.         err   = 0;
  695.         i_fam = [];
  696.         if isempty(varargin{1})
  697.             err = 1; 
  698.             msg = 'Wavelet Family (Short) Name is empty !';       
  699.         else
  700.             name = noblank(varargin{1});
  701.             tmpCell = {Wavelets_Info(:).familyName}';
  702.             i_fam = find(strcmp(name,tmpCell)==1);
  703.             if isempty(i_fam)
  704.                 tmpCell = {Wavelets_Info(:).familyShortName}';
  705.                 i_fam = find(strcmp(name,tmpCell)==1);
  706.             end
  707.         end
  709.         if err==0
  710.             if isempty(i_fam)
  711.                 err = 1;
  712.                 msg = 'Invalid Wavelet Family (Short) Name !';
  713.             elseif i_fam<=NB_FamInWTB
  714.                 err = 1;
  715.                 fn  = Wavelets_Info(i_fam).familyName;
  716.                 msg = sprintf('You can''t delete %s Wavelet Family !',fn);
  717.             end
  718.         end
  719.         if err
  720.             errargt(mfilename,msg,'msg');
  721.             return 
  722.         end
  724.         fid = fopen(asc_file);
  725.         if fid==-1 , fid = fopen(asc_ini_file); end
  726.         winfo = fread(fid);
  727.         fclose(fid);
  728.         fid = fopen(sav_file,'w');
  729.         fwrite(fid,winfo);
  730.         fclose(fid);
  732.         str_winfo = setstr(winfo');
  733.         str_NL    = findNL('string',winfo);
  734.         beg_fam   = 'Family Name : ';
  735.         first     = findstr(beg_fam,str_winfo);
  736.         first     = first(i_fam);
  737.         ind_NL    = findstr(str_NL',str_winfo);
  738.         indexs    = find(ind_NL>first);
  739.         indexs    = ind_NL(indexs(1:nb_InfoByWave+1));
  740.         last      = indexs(nb_InfoByWave+1)+length(str_NL)-1;
  742.         winfo(first:last) = [];
  743.         fid = fopen(asc_file,'w');
  744.         fwrite(fid,winfo);
  745.         fclose(fid);
  746.         wavemngr('create');
  748.     case 'restore'
  749.         wavemngr('clear')
  750.         if nargin==1
  751.             fid = fopen(sav_file);
  752.             if fid==-1 , fid = fopen(asc_ini_file); end
  753.         else
  754.             fid = fopen(asc_ini_file);
  755.         end
  756.         winfo = fread(fid);
  757.         fclose(fid);
  758.         fid = fopen(asc_file,'w');
  759.         fwrite(fid,winfo);
  760.         fclose(fid);
  761.         wavemngr('create');
  763.     case 'rescue'
  764.         famtype = ...
  765.            {...
  766.             1, 'Haar',         'haar',  1;
  767.             2, 'Daubechies',   'db',    1;
  768.             3, 'Symlets',      'sym',   1;
  769.             4, 'Coiflets',     'coif',  1;
  770.             5, 'BiorSplines',  'bior',  2;
  771.             6, 'ReverseBior',  'rbio',  2;
  772.             7, 'Meyer',        'meyr',  3;
  773.             8, 'DMeyer',       'dmey',  1;
  774.             9, 'Gaussian',     'gaus',  4;
  775.            10, 'Mexican_hat',  'mexh',  4;
  776.            11, 'Morlet',       'morl',  4;
  777.            10, 'Complex Gaussian',   'cgau',  5;
  778.            11, 'Shannon',            'shan',  5;
  779.            10, 'Frequency B-Spline', 'fbsp',  5;
  780.            11, 'Complex Morlet',     'cmor',  5
  782.             };
  784.         nums = ...
  785.            {...
  786.             '';
  787.             '1 2 3 4 5 6 7 8 9 10 **';
  788.             '2 3 4 5 6 7 8 **';
  789.             '1 2 3 4 5';
  790.             '1.1 1.3 1.5 2.2 2.4 2.6 2.8 3.1 3.3 3.5 3.7 3.9 4.4 5.5 6.8';
  791.             '1.1 1.3 1.5 2.2 2.4 2.6 2.8 3.1 3.3 3.5 3.7 3.9 4.4 5.5 6.8';
  792.             '';
  793.             '';
  794.             '1 2 3 4 5 6 7 8 **';
  795.             '';
  796.             '';
  797.             '1 2 3 4 5 **';
  798.             '1-1.5 1-1 1-0.5 1-0.1 2-3 **';
  799.             '1-1-1.5 1-1-1 1-1-0.5 2-1-1 2-1-0.5 2-1-0.1 **';
  800.             '1-1.5 1-1 1-0.5 1-1 1-0.5 1-0.1 **'
  801.             };
  803.         typNums = ...
  804.            {...
  805.             'no';
  806.             'integer';
  807.             'integer';
  808.             'integer';
  809.             'real';
  810.             'real';
  811.             'no';
  812.             'no';
  813.             'integer';
  814.             'no';
  815.             'no';
  816.             'integer';
  817.             'string';
  818.             'string';
  819.             'string'
  820.             };
  822.         files = ...
  823.            {...
  824.             'dbwavf';
  825.             'dbwavf';
  826.             'symwavf';
  827.             'coifwavf';
  828.             'biorwavf';
  829.             'rbiowavf';
  830.             'meyer';
  831.             'dmey.mat';
  832.             'gauswavf';
  833.             'mexihat';
  834.             'morlet'
  835.             'cgauwavf';
  836.             'shanwavf';
  837.             'fbspwavf';
  838.             'cmorwavf'
  839.             };
  841.         bounds = ...
  842.            {...
  843.             [];
  844.             [];
  845.             [];
  846.             [];
  847.             [];
  848.             [];
  849.             [-8 8];
  850.             [];
  851.             [-5 5];
  852.             [-8 8];
  853.             [-8 8];
  854.             [-5 5];
  855.             [-20 20];
  856.             [-20 20];
  857.             [-8 8]
  858.             };
  860.         nbfam   = size(famtype,1); 
  861.         tabNums = cell(nbfam,1);
  862.         for k = 1:nbfam
  863.             tabNums{k} = '';
  864.             s = deblank(nums{k});
  865.             if ~isempty(s)
  866.                 I = find(isspace(s));
  867.                 i_beg = 1;
  868.                 for j = 1:length(I)
  869.                     i_end = I(j)-1;
  870.                     ss = deblank(s(i_beg:i_end));
  871.                     if ~isempty(ss)
  872.                         tabNums{k} = strvcat(tabNums{k},ss);
  873.                     end
  874.                     i_beg = i_end+2;
  875.                 end
  876.                 ss = deblank(s(i_beg:end));
  877.                 tabNums{k} = strvcat(tabNums{k},ss);
  878.             end
  879.             if isempty(tabNums{k}) , tabNums{k} = 'no'; end
  880.         end
  882.         varargout{1} = ...
  883.            struct(...
  884.                   'index',           famtype(:,1), ...
  885.                   'familyName',      famtype(:,2), ...
  886.                   'familyShortName', famtype(:,3), ...
  887.                   'type',            famtype(:,4), ...
  888.                   'tabNums',         tabNums,      ...
  889.                   'typNums',         typNums,      ...
  890.                   'file',            files,        ...
  891.                   'bounds',          bounds        ...
  892.                   );
  894.     otherwise
  895.         errargt(mfilename,'Unknown Option','msg');
  896.         error('*');
  897. end
  899. %----------------------------------%
  900. % Internal Function(s)             %
  901. %----------------------------------%
  902. function nl = findNL(opt,txt)
  903. %FINDNL Find New Line char or string.
  905. ind10 = find(txt==10);
  906. ind13 = find(txt==13);
  907. if isempty(ind13) ,
  908.     nl = 10;
  909. elseif isempty(ind10)
  910.     nl = 13;
  911. else
  912.     nl = [13;10];
  913. end
  914. if isequal(opt,'string')
  915.     nl = setstr(nl);
  916. end
  917. if size(txt,2)>1 , nl = nl'; end
  918. %----------------------------------%
  919. function s = noblank(x)
  920. %NOBLANK Removes blanks in a string.
  922. if ~isempty(x)
  923.     s = x(find(x~=' ' & x~=0));
  924. else
  925.     s = '';
  926. end
  927. %----------------------------------%