order patterns recurrence plot added

crqa.m

@@ -56,7 +56,8 @@ function xout=crqa(varargin)
 %                    (all vectors have the length one).
 %      fan         - Fixed amount of nearest neighbours.
 %      inter       - Interdependent neighbours.
-%      distance    - Distance coded matrix (global CRP, Euclidean norm).
+%      omatrix     - Order matrix.
+%      opattern    - Order patterns recurrence plot.
 %
 %    Normalization of the data series.
 %      normalize   - Normalization of the data.
@@ -134,6 +135,9 @@ function xout=crqa(varargin)
 % $Revision$
 %
 % $Log$
+% Revision 5.5  2004/11/10 07:04:50  marwan
+% initial import
+%
 %
 % This program is part of the new generation XXII series.
 %
@@ -178,7 +182,7 @@ errcode=1;
    w=[];wstep=0; method='max'; method_n=1;
    embed_flag=1;
    nogui=0;
-      check_meth={'ma','eu','mi','nr','fa','in','di'}; 	% maxnorm, euclidean, nrmnorm,  fan, distance
+      check_meth={'ma','eu','mi','nr','fa','in','om','op','di'}; 	% maxnorm, euclidean, nrmnorm,  fan, distance
       check_gui={'gui','nog','sil'};		% gui, nogui, silent
       check_norm={'non','nor'};			% nonormalize, normalize
    
@@ -195,7 +199,7 @@ errcode=1;
      h=findobj('Tag','crqa_eps');
      e=str2num(get(h(1),'String'));
      h=findobj('Tag','crqa_method');
-     method={'Maximum Norm','Euclidean Norm','Minimum Norm','Normalized Norm','FAN','Interdependent'};
+     method={'Maximum Norm','Euclidean Norm','Minimum Norm','Normalized Norm','FAN','Interdependent','Order Matrix','Order Pattern','Distance Plot'};
      method=method{get(h(1),'Value')};
      nonorm=get(h(1),'UserData');
      h=findobj('Tag','crqa_lmin');
@@ -232,7 +236,7 @@ errcode=1;
          if isempty(nogui), nogui=0; end
          if isempty(method_n), method_n=1; end
          if nogui>2, nogui=1; end
-         if method_n>7, method0=7; end
+         if method_n>length(check_meth), method0=length(check_meth); end
          method=check_meth{method_n};
       else
          nogui=0;
@@ -570,7 +574,7 @@ case 'init'
             'Tag','crqa_method',...
 	    'UserData',nonorm,...
 	    'Value',method_n,...
-	    'String','Maximum Norm|Euclidean Norm|Minimum Norm|Normalized Norm|Fixed Amount|Interdependent',...
+	    'String','Maximum Norm|Euclidean Norm|Minimum Norm|Normalized Norm|Fixed Amount|Interdependent|Order Matrix|Order Patterns|Distance Plot',...
             'Position',[89+30 24.6+.2+2.5 22  1.7]);    
 
   h=uicontrol(props.text,...
@@ -777,6 +781,12 @@ case 'store'
 case 'compute'
 
   errcode=11;
+  
+  if length(method)>1 & strcmpi(method(1:3),'dis')
+      disp('Warning: RQA from distance plot not possible!')
+      return
+  end
+  
   if ~nogui
     h_fig=findobj('tag','crqa_Fig');
     setptr(gcf,'watch'), 

crqad.m

@@ -34,7 +34,8 @@ function out=crqad(varargin)
 %                    (all vectors have the length one).
 %      fan         - Fixed amount of nearest neighbours.
 %      inter       - Interdependent neighbours.
-%      distance    - Distance coded matrix (global CRP, Euclidean norm).
+%      omatrix     - Order matrix.
+%      opattern    - Order patterns recurrence plot.
 %
 %    Normalization of the data series.
 %      normalize   - Normalization of the data.
@@ -69,6 +70,9 @@ function out=crqad(varargin)
 % $Revision$
 %
 % $Log$
+% Revision 2.2  2004/11/10 07:05:55  marwan
+% initial import
+%
 
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% programme properties
@@ -100,7 +104,7 @@ splash_gpl('crp');
    if nargin & isnumeric(varargin{1})
 
      % check the text input parameters for method, gui 
-      check_meth={'ma','eu','nr','fa','in','di'}; 	% maxnorm, euclidean, nrmnorm,  fan, distance
+      check_meth={'ma','eu','mi','nr','fa','in','om','op','di'}; 	% maxnorm, euclidean, nrmnorm,  fan, distance
       check_gui={'gui','nog','sil'};		% gui, nogui, silent
       temp_meth=0;
       temp_gui=0;
@@ -117,7 +121,7 @@ splash_gpl('crp');
          if isempty(nogui), nogui=0; end
          if isempty(method_n), method_n=1; end
          if nogui>2, nogui=1; end
-         if method_n>6, method0=6; end
+         if method_n>length(check_meth), method0=length(check_meth); end
          method=check_meth{method_n};
       else
          nogui=0;
@@ -218,6 +222,10 @@ splash_gpl('crp');
 flag=1;
 
 x1=x; x2=y;
+  if length(method)>1 & strcmpi(method(1:2),'di')
+      disp('Warning: RQA from distance plot not possible!')
+      return
+  end
 
 warning off
 if size(x1,1)<size(x1,2), x1=x1'; end