function xout=jrp(varargin)
%JRP Creates a joint recurrence plot.
% JRP(X [,Y] [,param1,param2,...]) creates a simple
% joint recurrence plot/ recurrence plot of maximal
% two data series. Embedding parameters will be the
% same for both systems. Results can be stored into
% the workspace.
%
% R=JRP(X,M,T,E) uses the dimension M, delay T
% and the size of neighbourhood E and creates a recurrence
% plot of X.
%
% R=JRP(X,Y,'order') creates an order matrix plot
% with normalization of the data.
%
% R=JRP(X,Y,'distance','nonormalize') creates a
% distance coded matrix plot without normalization
% of the data.
%
% Allows to change the parameters interactively by
% using a GUI.
%
% The source-data X and test-data Y can be one- or
% a two-coloumn vectors (then, in the first column
% have to be the time); if the test-data Y is not
% specified, a simple recurrence plot is created.
%
% Parameters: dimension M, delay T and the size of
% neighbourhood E are the first three numbers after
% the data series; further parameters can be used
% to switch between various methods of finding the
% neighbours of the phasespace trajectory, to suppress
% the normalization of the data and to suppress the
% GUI (useful in order to use this programme by other
% programmes).
%
% Methods of finding the neighbours/ of plot.
% maxnorm - Maximum norm.
% euclidean - Euclidean norm.
% minnorm - Minimum norm.
% nrmnorm - Euclidean norm between normalized vectors
% (all vectors have the length one).
% fan - Fixed amount of nearest neighbours.
% inter - Interdependent neighbours.
% omatrix - Order matrix.
% opattern - Order patterns recurrence plot.
% distance - Distance coded matrix (global JRP, Euclidean norm).
%
% Normalization of the data series.
% normalize - Normalization of the data.
% nonormalize - No normalization of the data.
%
% Suppressing the GUI.
% gui - Creates the GUI and the output plot.
% nogui - Suppresses the GUI and the output plot.
% silent - Suppresses all output.
%
% Parameters not needed to be specified.
%
% Current limitation: for higher speed in
% output the whole matrix of the recurrence
% plot is in the work space - this limits
% the application of long data series.
%
% Examples: a = sin((1:1000) * 2 * pi/67);
% b = sin(.01 * ([1:1000] * 2 * pi/67) .^ 2);
% jrp(a,b,3,12,'fan')
%
% See also CRP, JRQA.
%
% References:
% Romano, M., Thiel, M., Kurths, J., von Bloh, W.:
% Multivariate Recurrence Plots, Phys. Lett. A , 330, 2004.
% Copyright (c) 2004-2005 by AMRON
% Norbert Marwan, Potsdam University, Germany
% http://www.agnld.uni-potsdam.de
%
% $Date$
% $Revision$
%
% $Log$
% Revision 2.2 2005/04/15 09:02:32 marwan
% minor bugfix in plugin section
%
% Revision 2.1 2005/04/08 09:54:19 marwan
% plugin added
%
% Revision 1.3 2005/04/06 13:00:56 marwan
% *** empty log message ***
%
% Revision 1.2 2005/04/04 09:52:53 marwan
% bug in distance plot and colormap selection fixed
%
% Revision 1.1 2005/03/16 16:22:29 marwan
% support for joint recurrence plots added
%
%
%
% This program is part of the new generation XXII series.
%
% This program is free software; you can redistribute it and/or
% modify it under the terms of the GNU General Public License
% as published by the Free Software Foundation; either version 2
% of the License, or any later version.
warning off
global errcode props
errcode=0;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% programme properties
init_properties
hCRP=[];hCtrl=[];nogui=[];obj=[];mflag=[];
set(0,'ShowHidden','On')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% check and read the input
error(nargchk(1,8,nargin));
if nargout>1, error('Too many output arguments'), end
check_meth={'ma','eu','mi','nr','fa','in','om','op','di'}; % maxnorm, euclidean, nrmnorm, fan, distance
check_norm={'non','nor'}; % nonormalize, normalize
check_gui={'gui','nog','sil'}; % gui, nogui, silent
if isnumeric(varargin{1})==1 % read commandline input
varargin{9}=[];
i_double=find(cellfun('isclass',varargin,'double'));
i_char=find(cellfun('isclass',varargin,'char'));
% check the text input parameters for method, gui and normalization
temp_meth=0;
temp_norm=0;
temp_gui=0;
if ~isempty(i_char)
for i=1:length(i_char),
varargin{i_char(i)}(4)='0';
temp_meth=temp_meth+strcmpi(varargin{i_char(i)}(1:2),check_meth');
temp_norm=temp_norm+strcmpi(varargin{i_char(i)}(1:3),check_norm');
temp_gui=temp_gui+strcmpi(varargin{i_char(i)}(1:3),check_gui');
end
method=min(find(temp_meth));
nonorm=min(find(temp_norm))-1;
nogui=min(find(temp_gui))-1;
if isempty(method), method=1; end
if isempty(nonorm), nonorm=1; end
if isempty(nogui), nogui=0; end
if method>length(check_meth), method=length(check_meth); end
if nonorm>1, nonorm=1; end
if nogui>2, nogui=2; end
else
method=1; nonorm=1; nogui=0;
end
if nogui==0 & nargout>0, nogui=1; end
% get the parameters for creating RP
if max(size(varargin{1}))<=3
error('To less values in data X.')
end
x=double(varargin{1});
% if isempty(varargin{2}) | ~isnumeric(varargin{2}), y=x; end
% if ~isempty(varargin{2}) & isnumeric(varargin{2}), y=double(varargin{2}); end
if isempty(varargin{2}) | ~isnumeric(varargin{2}), y=x;
else, y=double(varargin{2}); end
if (isnumeric(varargin{2}) & max(size(varargin{2}))==1) | ~isnumeric(varargin{2})
y=x;
if ~isempty(varargin{i_double(2)}), m=varargin{i_double(2)}(1); else m=1; end
if ~isempty(varargin{i_double(3)}), t=varargin{i_double(3)}(1); else t=1; end
if ~isempty(varargin{i_double(4)}), e=varargin{i_double(4)}(1); else e=.1; end
else
if ~isempty(varargin{i_double(3)}), m=varargin{i_double(3)}(1); else m=1; end
if ~isempty(varargin{i_double(4)}), t=varargin{i_double(4)}(1); else t=1; end
if ~isempty(varargin{i_double(5)}), e=varargin{i_double(5)}(1); else e=.1; end
end
if e<0, e=1; disp('Warning: The threshold size E cannot be negative and is now set to 1.'), end
if t<1, t=1; disp('Warning: The delay T cannot be smaller than one and is now set to 1.'), end
t=round(t); m=round(m); mflag=method;
action='init';
if max(size(x))~=max(size(y)),
if ~nogui, errordlg('Data must have the same length.','Check Data'), waitforbuttonpress, return, else error('Data must have the same length.'), end
end
Nx=length(x); Ny=length(y);
NX=Nx-t*(m-1);NY=Ny-t*(m-1);
if (NX<1 | NY<1) & nogui==0;
errordlg('The embedding vectors cannot be created. Dimension M and/ or delay T are to big. Please use smaller values.','Dimension/ delay to big')
waitforbuttonpress
end
if size(x,1)<size(x,2), x=x'; end
if size(y,1)<size(y,2), y=y'; end
if ~isempty(find(isnan(x)))
disp('Warning: NaN detected (in first variable) - will be cleared.')
for k=1:size(x,2), x(find(isnan(x(:,k))),:)=[]; end
end
if ~isempty(find(isnan(y)))
disp('Warning: NaN detected (in second variable) - will be cleared.')
for k=1:size(y,2), y(find(isnan(y(:,k))),:)=[]; end
end
if size(x,2)>=2
xscale=x(:,1);
if ~isempty(find(diff(xscale)<0)), error('First column of the first vector must be monotonically non-decreasing.'),end
if nonorm==1, x=(x(:,2)-mean(x(:,2)))/std(x(:,2)); else x=x(:,2); end
else
if nonorm==1, x=(x-mean(x))/std(x); end
xscale=(1:length(x))';
end
if size(y,2)>=2
yscale=y(:,1);
if ~isempty(find(diff(yscale)<0)), error('First column of the second vector must be monotonically non-decreasing.'),end
if nonorm==1, y=(y(:,2)-mean(y(:,2)))/std(y(:,2)); else y=y(:,2); end
else
if nonorm==1, y=(y-mean(y))/std(y); end
yscale=(1:length(y))';
end
if method==7 & m > 1,
m=1;
disp('Warning: For order matrix a dimension of one is used.')
end
ds=eye(m);
else % read input from the GUI
action=varargin{1};
nogui=0;
h=get(gcf,'Name');h=h(findstr(h,'(')+1:findstr(h,')')-1);
hCRP=findobj('Name',['Joint Recurrence Plot (' h ')']);
hCtrl=findobj('Name',['Control (' h ')']);
h=str2num(h);
xshuttle=get(findobj('Parent',hCRP,'Tag','DataPlot2'),'UserData');
if ~isempty(xshuttle)
xscale=xshuttle(:,1);
x=xshuttle(:,2);
yshuttle=get(findobj('Parent',hCRP,'Tag','DataPlot1'),'UserData');
yscale=yshuttle(:,1);
y=yshuttle(:,2);
if ~isempty(hCtrl)
if get(findobj('Tag','Unthresh','Parent',hCtrl),'Value')
mflag=length(check_meth);
else
mflag=get(findobj('Tag','Method','Parent',hCtrl),'Value');
end
m=get(findobj('Tag','Dim','Parent',hCtrl),'Value');
t=str2num(get(findobj('Tag','Delay','Parent',hCtrl),'String'));
e=str2num(get(findobj('Tag','Size','Parent',hCtrl),'String'));
if e<0, e=1;
errordlg('The threshold size E cannot be negative.','Threshold size to small')
waitforbuttonpress
set(findobj('Tag','Size','Parent',hCtrl),'String','1')
action='';
end
if t<1, t=1;
errordlg('The delay T cannot be smaller than one.','Delay to small')
waitforbuttonpress
set(findobj('Tag','Delay','Parent',hCtrl),'String','1')
action='';
end
ds = get(findobj('Tag','Dim','Parent',hCtrl),'UserData');
if mflag==7 | mflag==8 | mflag==length(check_meth)
set(findobj('Tag','Size','Parent',hCtrl),'enable','off');
set(findobj('Tag','Sizetext','Parent',hCtrl),'enable','off');
else
set(findobj('Tag','Size','Parent',hCtrl),'enable','on');
set(findobj('Tag','Sizetext','Parent',hCtrl),'enable','on');
end
if mflag==7 & m > 1,
m=1; ds = 1;
disp('Warning: For order matrix a dimension of one is used.')
end
Nx=length(x); Ny=length(y);
NX=Nx-t*(m-1);NY=Ny-t*(m-1);
if (NX<1 | NY<1) & strcmpi(action,'apply');
errordlg('The embedding vectors cannot be created. Dimension M and/ or delay T are to big. Please use smaller values.','Dimension/ delay to big')
waitforbuttonpress
action='';
end
end
end
clear xshuttle yshuttle
cm_old=get(hCRP,'Colormap');
cm=[{hsv(256)}; {hot(256)}; {gray(256)};...
{french(256)}; {bone(256)}; {copper(256)}; {pink(256)};...
{flag(256)}; {lines(256)}; {colorcube(256)};...
{jet(256)}; {prism(256)}; {cool(256)};...
{autumn(256)}; {spring(256)}; {winter(256)};...
{summer(256)}; {flipud(gray(2))}; {flipud(cm_old)}];
if isempty(findobj('Tag','CRPFig')) & nogui==0
action='init';
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% splash the GPL
splash_gpl('jrp')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% nogui
if nogui>0
hCRP=9999;
if nogui~=2
tx(1)={'maximum norm'};
tx(2)={'Euclidean norm'};
tx(3)={'minimum norm'};
tx(4)={'Euclidean norm of normalized distance'};
tx(5)={'fixed amount of nearest neighbours'};
tx(6)={'interdependent neighbours'};
tx(7)={'order matrix'};
tx(8)={'order pattern'};
tx(9)={'distance plot'};
disp(['use method: ', char(tx(method))]);
if nonorm==1, disp('normalize data'); else disp('do not normalize data'); end
end
action='compute';
if (NX<1 | NY<1)
disp('Warning: The embedding vectors cannot be created.')
disp('Dimension M and/ or delay T are to big. Please use smaller values.')
action='';
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% switch routines
try
switch(action)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% initialization
case 'init'
errcode=1;
xshuttle(:,1)=xscale;
yshuttle(:,1)=yscale;
xshuttle(:,2)=x;
yshuttle(:,2)=y;
ds=eye(m);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% create GUI
errcode=2;
root_ud=get(0,'UserData');
if isstruct(root_ud)
if isfield(root_ud,'crp')
if ~isempty(root_ud.crp)
root_ud.crp=[root_ud.crp max(root_ud.crp)+1];
else
root_ud.crp=1;
end
h=num2str(root_ud.crp(end));
else
root_ud.crp=1;
h=num2str(1);
end
else
root_ud.old=root_ud;
root_ud.crp=1;
h=num2str(1);
end
set(0,'UserData',root_ud)
%%%%%%%%%%%%%%%%% JRP Figure
h_axes=create_JRPfig(h,xshuttle,yshuttle);
if nargout>0, xout=haxes; end
%%%%%%%%%%%%%%%%% Control Figure
errcode=3;
create_Ctrlfig('jrp',h,ds,m,t,e,method)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% vectorswitch
case 'vectorswitch'
errcode=4;
vectorswitch
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% fit dimension display
case 'fitdim'
errcode=5;
ds2=eye(m);
if m>length(ds)
ds2(1:length(ds),1:length(ds))=ds;
else
ds2=ds(1:m,1:m);
end
ds=ds2; clear ds2
set(findobj('Tag','Dim','Parent',gcf),'UserData', ds);
for i=1:20
if i>m, set(findobj('Tag',['DimShift' num2str(i)],'Parent',gcf), 'Enable', 'off');
else, set(findobj('Tag',['DimShift' num2str(i)],'Parent',gcf), 'Enable', 'on'); end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% unthresh
case 'unthresh'
errcode=6;
switch_unthresholded(hCRP)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% stretch
case 'stretch'
errcode=7;
stretch(hCRP,xscale,yscale,Nx,Ny)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% change colormap
case 'log'
errcode=82;
change_colormapscale(hCRP,cm)
case 'colormap'
errcode=81;
change_colormap(hCtrl,hCRP,cm)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% store
case 'store'
errcode=9;
X=get(findobj('Tag','CRPData','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'UserData');
if isempty(X)
warndlg('The JRP matrix is still empty. Please start the computation before storing.','Empty JRP')
waitforbuttonpress
else
assignin('base','X', double(X))
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% handlevisON
case 'handlevisON'
set(hCRP, 'HandleVis','on')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% close
case 'close'
errcode=101;
close_all
case 'smartclose'
errcode=102;
if ishandle(hCRP) & ishandle(hCtrl)
smart_close(hCRP,hCtrl)
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% computation
case 'compute'
errcode=11;
txt_cross='Joint ';
if size(x)==size(y)
if x==y, txt_cross=''; end
end
if nogui==0
setptr([hCRP,hCtrl],'watch')
h=findobj('Tag','Apply','Parent',hCtrl);
obj.children=get(hCtrl,'Children');
obj.enable=get(obj.children,'Enable');
set(obj.children,'Enable','off')
set(h(1),'String','Stop',...
'ToolTip','Stops the computation.',...
'Enable','on',...
'Callback','set(gcbo,''String'',''Stopped'')')
end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'Visible','on')
set(findobj('Tag','CRPData','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'Visible','off')
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Building Embedding Vectors'),drawnow
X=uint8(zeros(NY,NX));
if check_stop(hCRP,hCtrl,nogui,obj), return, end
% embedding vectors
i=(1:NX)';j=0:t:0+(m-1)*t;
i=reshape(i(:,ones(1,m))+j(ones(NX,1),:),m*NX,1);
x1=x(i);
x2=reshape(x1,NX,m);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
i=(1:NY)';j=0:t:0+(m-1)*t;
i=reshape(i(:,ones(1,m))+j(ones(NY,1),:),m*NY,1);
y1=y(i);
y2=reshape(y1,NY,m);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
y2=y2*ds; % switch vectors
% check if plugin exist and is executable
[plugin_exist, plugin_name, plugin_path] = is_crp_plugin;
if plugin_exist & ( mflag < 4 | mflag == 9 ) & length(x) == length(y) % if plugin exist and method is MAX, MIN, EUC ord DIS
if nogui == 1, disp('(plugin used)'), end
[X1 matext] = crp_plugin(x2, x2, 1, 1, e, mflag, hCRP, plugin_path, 0);
[X2 matext] = crp_plugin(y2, y2, 1, 1, e, mflag, hCRP, plugin_path, 0);
X = uint8(X1' .* X2)'; clear X1 X2
else
% use builtin implementation
[NX, mx] = size(x2);
[NY, my] = size(y2);
clear jx jy x1 y1
switch(mflag)
%%%%%%%%%%%%%%%%% local JRP, fixed distance
case {1,2,3}
errcode=111;
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Reshape Embedding Vectors (1)'),drawnow
px = permute(x2, [ 1 3 2 ]);
py = permute(x2, [ 3 1 2 ]);
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Compute Distance Matrix (1)'),drawnow
if check_stop(hCRP,hCtrl,nogui,obj), return, end
s1 = px(:,ones(1,NX),:) - py(ones(1,NX),:,:);
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Reshape Embedding Vectors (2)'),drawnow
px = permute(y2, [ 1 3 2 ]);
py = permute(y2, [ 3 1 2 ]);
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Compute Distance Matrix (2)'),drawnow
if check_stop(hCRP,hCtrl,nogui,obj), return, end
s2 = px(:,ones(1,NY),:) - py(ones(1,NY),:,:);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
switch mflag
case 1
%%%%%%%%%%%%%%%%% maximum norm
sx = max(abs(s1),[],3);
sy = max(abs(s2),[],3);
matext=[num2str(round(100*e)/100) '\sigma (fixed distance maximum norm)'];
case 2
%%%%%%%%%%%%%%%%% euclidean norm
errcode=112;
sx = sqrt(sum(s1.^2, 3));
sy = sqrt(sum(s2.^2, 3));
matext=[num2str(round(100*e)/100) '\sigma (fixed distance euclidean norm)'];
case 3
%%%%%%%%%%%%%%%%% minimum norm
errcode=113;
sx = sum(abs(s1), 3);
sy = sum(abs(s2), 3);
matext=[num2str(round(100*e)/100) '\sigma (fixed distance minimum norm)'];
end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')'),'String','Building JRP Matrix'),drawnow
if check_stop(hCRP,hCtrl,nogui,obj), return, end
X1=sx<e;
X2=sy<e;
if check_stop(hCRP,hCtrl,nogui,obj), return, end
X = uint8(X1' .* X2)'; clear s s1 s2 sx sy x1 y1 px py X1 X2
%%%%%%%%%%%%%%%%% local JRP, normalized distance euclidean norm
case 4
errcode=114;
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Normalize Embedding Vectors'),drawnow
Dx=sqrt(sum(((x2(:,:)).^2)'))';
if check_stop(hCRP,hCtrl,nogui,obj), return, end
Dy=sqrt(sum(((y2(:,:)).^2)'))';
if check_stop(hCRP,hCtrl,nogui,obj), return, end
x1=x2./repmat(Dx,1,m);x2=x1;
if check_stop(hCRP,hCtrl,nogui,obj), return, end
y1=y2./repmat(Dy,1,m);y2=y1; clear Dx Dy y1 x1
if check_stop(hCRP,hCtrl,nogui,obj), return, end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Reshape Embedding Vectors (1)'),drawnow
px = permute(x2, [ 1 3 2 ]);
py = permute(x2, [ 3 1 2 ]);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Compute Distance Matrix (1)'),drawnow
s1 = px(:,ones(1,NX),:) - py(ones(1,NX),:,:);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
sx = sqrt(sum(s1.^2, 3));
if check_stop(hCRP,hCtrl,nogui,obj), return, end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Reshape Embedding Vectors (2)'),drawnow
px = permute(y2, [ 1 3 2 ]);
py = permute(y2, [ 3 1 2 ]);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Compute Distance Matrix (2)'),drawnow
s1 = px(:,ones(1,NY),:) - py(ones(1,NY),:,:);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
sy = sqrt(sum(s1.^2, 3));
if check_stop(hCRP,hCtrl,nogui,obj), return, end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')'),'String','Building JRP Matrix'),drawnow
X=uint8( (255*sx/max(sx(:)))<(255*e/max(sx(:)))' .* (255*sy/max(sy(:)))<(255*e/max(sy(:)))' );
clear sx sy s1 x1 y1 px py
matext=[num2str(round(100*e)/100) '\sigma (normalized distance euclidean norm)'];
%%%%%%%%%%%%%%%%% local JRP, fixed neigbours amount
case 5
errcode=115;
if e>=1
e=round(e)/100;
txt=['The value for fixed neigbours amount has to be smaller '...
'than one. Continue the computation with a value of ' ...
num2str(e)];
if nogui==0
warndlg(txt,'Threshold value mismatch');
drawnow
waitforbuttonpress
set(findobj('Tag','Size','Parent',gcf),'String',num2str(e))
end
end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Reshape Embedding Vectors (1)'),drawnow
px = permute(x2, [ 1 3 2 ]);
py = permute(x2, [ 3 1 2 ]);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Compute Distance Matrix (1)'),drawnow
s1 = px(:,ones(1,NX),:) - py(ones(1,NX),:,:);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
% s = sqrt(sum(s1.^2, 3));
s = (sum(s1.^2, 3));
if check_stop(hCRP,hCtrl,nogui,obj), return, end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Sort Distance Matrix (1)'),drawnow
mine=round(NY*e);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
[SS, JJ]=sort(s');
JJ=JJ';
if check_stop(hCRP,hCtrl,nogui,obj), return, end
X1(NX*NX)=0; X1(JJ(:,1:mine)+repmat([0:NX:NX*NX-1]',1,mine))=1;
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Reshape Embedding Vectors (2)'),drawnow
px = permute(y2, [ 1 3 2 ]);
py = permute(y2, [ 3 1 2 ]);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Compute Distance Matrix (2)'),drawnow
s1 = px(:,ones(1,NY),:) - py(ones(1,NY),:,:);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
% s = sqrt(sum(s1.^2, 3));
s = (sum(s1.^2, 3));
if check_stop(hCRP,hCtrl,nogui,obj), return, end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Sort Distance Matrix (2)'),drawnow
mine=round(NY*e);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
[SS, JJ]=sort(s');
JJ=JJ';
if check_stop(hCRP,hCtrl,nogui,obj), return, end
X2(NY*NY)=0; X2(JJ(:,1:mine)+repmat([0:NY:NY*NY-1]',1,mine))=1;
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Building JRP Matrix'),drawnow
if check_stop(hCRP,hCtrl,nogui,obj), return, end
X=uint8(reshape(X1,NX,NX) .* reshape(X2,NY,NY)); clear X1 X2 SS JJ s px py
matext=[num2str(round(1000*mine/NY)/10) '% (fixed neighbours amount)'];
%%%%%%%%%%%%%%%%% local JRP, interdependent neigbours
case 6
errcode=116;
if e>=1
e=round(e)/100;
txt=['The value for fixed neigbours amount has to be smaller '...
'than one. Continue the computation with a value of ' ...
num2str(e)];
if nogui==0
warndlg(txt,'Threshold value mismatch');
drawnow
waitforbuttonpress
set(findobj('Tag','Size','Parent',gcf),'String',num2str(e))
end
end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Reshape Embedding Vectors (1)'),drawnow
px = permute(x2, [ 1 3 2 ]);
py = permute(x2, [ 1 3 2 ]);
px2 = permute(x2, [ 3 1 2 ]);
py2 = permute(x2, [ 3 1 2 ]);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Compute Distance Matrix (1)'),drawnow
s1 = px(:,ones(1,NX),:) - px2(ones(1,NX),:,:);
sx = sqrt(sum(s1.^2, 3));
if check_stop(hCRP,hCtrl,nogui,obj), return, end
s1 = py(:,ones(1,NX),:) - py2(ones(1,NX),:,:);
sy = sqrt(sum(s1.^2, 3));
if check_stop(hCRP,hCtrl,nogui,obj), return, end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Sort Distance Matrix (1)'),drawnow
mine=round(min(NX,NX)*e);
[SSx, JJx]=sort(sx);%SSx(1,:)=[]; JJx(1,:)=[];
if check_stop(hCRP,hCtrl,nogui,obj), return, end
[SSy, JJy]=sort(sy);%SSy(1,:)=[]; JJy(1,:)=[];
if check_stop(hCRP,hCtrl,nogui,obj), return, end
e=mean(SSy(mine:mine+1,:));
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Building RP Matrix (1)'),drawnow
for i=1:min(NX,NX),
if check_stop(hCRP,hCtrl,nogui,obj), return, end
JJx((JJx(1:mine,i)>min(NX,NX)),i)=1;
X1(i,JJx(1:mine,i))=(sy(i,JJx(1:mine,i))<=e(i))';
end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Reshape Embedding Vectors (2)'),drawnow
px = permute(y2, [ 1 3 2 ]);
py = permute(y2, [ 1 3 2 ]);
px2 = permute(y2, [ 3 1 2 ]);
py2 = permute(y2, [ 3 1 2 ]);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Compute Distance Matrix (2)'),drawnow
s1 = px(:,ones(1,NY),:) - px2(ones(1,NY),:,:);
sx = sqrt(sum(s1.^2, 3));
if check_stop(hCRP,hCtrl,nogui,obj), return, end
s1 = py(:,ones(1,NY),:) - py2(ones(1,NY),:,:);
sy = sqrt(sum(s1.^2, 3));
if check_stop(hCRP,hCtrl,nogui,obj), return, end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Sort Distance Matrix (2)'),drawnow
mine=round(min(NY,NY)*e);
mine(~mine) = 1;
[SSx, JJx]=sort(sx);%SSx(1,:)=[]; JJx(1,:)=[];
if check_stop(hCRP,hCtrl,nogui,obj), return, end
[SSy, JJy]=sort(sy);%SSy(1,:)=[]; JJy(1,:)=[];
if check_stop(hCRP,hCtrl,nogui,obj), return, end
e=mean(SSy(mine:mine+1,:));
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Building RP Matrix (2)'),drawnow
for i=1:min(NY,NY),
if check_stop(hCRP,hCtrl,nogui,obj), return, end
JJx((JJx(1:mine,i)>min(NY,NY)),i)=1;
X2(i,JJx(1:mine,i))=sy(i,JJx(1:mine,i))<=e(i);
end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Building JRP Matrix'),drawnow
X = X1.* X2;
clear X1 X2 SS* JJ* s sx sy s1 px py px2 py2
matext=[num2str(round(1000*mine/NY)/10) '% (interdependent neighbours)'];
%%%%%%%%%%%%%%%%% order matrix
case 7
errcode=117;
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Compute Order Matrix (1)'),drawnow
px = permute(x2, [ 1 3 2 ]);
py = permute(x2, [ 3 1 2 ]);
X1 = px(:,ones(1,NX),:) >= py(ones(1,NX),:,:) - e;
if check_stop(hCRP,hCtrl,nogui,obj), return, end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Compute Order Matrix (2)'),drawnow
px = permute(y2, [ 1 3 2 ]);
py = permute(y2, [ 3 1 2 ]);
X2 = px(:,ones(1,NY),:) >= py(ones(1,NY),:,:) - e;
X = X1' .* X2';
if check_stop(hCRP,hCtrl,nogui,obj), return, end
clear X1 X2 s sx sy s1 px py
matext='';
%%%%%%%%%%%%%%%%% order pattern recurrence plot
case 8
errcode=118;
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Compute Order Patterns'),drawnow
% create a order pattern test
cmdStr = '';
for i=2:m;
cmdStr = [cmdStr, ' permX(:,', num2str(i-1) ,') < permX(:,', num2str(i), ') + eps'];
if i < m, cmdStr = [cmdStr, ' &']; end
end
if m==1
cmdStr = '1';
disp('Warning: No order patterns for dimension one; please use higher dimension!')
end
% order patterns by permutation of the set of values
clear patt*
for i=1:length(x2)
permX=perms(x2(i,:));
orderPattern = find(eval(cmdStr));
if isempty(orderPattern) orderPattern = 0; end
pattX(i) = orderPattern(1);
end
for i=1:length(y2)
permX=perms(y2(i,:));
orderPattern = find(eval(cmdStr));
if isempty(orderPattern) orderPattern = 0; end
pattY(i) = orderPattern(1);
end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Create Order Patterns Matrix (1)'),drawnow
px = permute(pattX', [ 1 3 2 ]);
py = permute(pattX', [ 3 1 2 ]);
X1 = px(:,ones(1,length(x2)),:) == py(ones(1,length(x2)),:,:);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Create Order Patterns Matrix (2)'),drawnow
px = permute(pattY', [ 1 3 2 ]);
py = permute(pattY', [ 3 1 2 ]);
X2 = px(:,ones(1,length(y2)),:) == py(ones(1,length(y2)),:,:);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Building JRP Matrix'),drawnow
X = X1' .* X2';
if check_stop(hCRP,hCtrl,nogui,obj), return, end
clear X1 X2 sx sy s1 px py pattY pattX orderPattern permX permY
matext='';
%%%%%%%%%%%%%%%%% global JRP
case length(check_meth)
errcode=110+length(check_meth);
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Compute Distance Matrix (1)'),drawnow
px = permute(x2, [ 1 3 2 ]);
py = permute(x2, [ 3 1 2 ]);
s1 = px(:,ones(1,NX),:) - py(ones(1,NX),:,:);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
sx = sqrt(sum(s1.^2, 3))';
if check_stop(hCRP,hCtrl,nogui,obj), return, end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Compute Distance Matrix (2)'),drawnow
px = permute(y2, [ 1 3 2 ]);
py = permute(y2, [ 3 1 2 ]);
s1 = px(:,ones(1,NY),:) - py(ones(1,NY),:,:);
if check_stop(hCRP,hCtrl,nogui,obj), return, end
sy = sqrt(sum(s1.^2, 3))'
if check_stop(hCRP,hCtrl,nogui,obj), return, end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')),'String','Building JRP Matrix'),drawnow
% s=exp(-(s^2)/(4*.1^2));
X = double(sqrt(sx.^2 + sy.^2));
size(X)
matext='';
clear sx sy s1 px py
end
end % end plugin
if nogui==0
for i=1:length(obj.enable), set(obj.children(i),'Enable',obj.enable{i}); end
set(h(1),'String','Apply',...
'ToolTip','Starts the computation - be patient.',...
'Callback','jrp compute')
setptr([hCRP,hCtrl],'arrow')
end
%%%%%%%%%%%%%%%%% show JRP
if nogui==0
Shuttle.hCRP=hCRP;
Shuttle.hCtrl=hCtrl;
Shuttle.matext=matext;
Shuttle.xscale=xscale;
Shuttle.yscale=yscale;
Shuttle.mflag=mflag;
Shuttle.m=m;
Shuttle.t=t;
Shuttle.cm=cm;
Shuttle.txt_cross=txt_cross;
show_crp(X,Shuttle)
else
if nargout==1, xout=X;end
end
end
warning on
try, set(0,props.root), end
set(0,'ShowHidden','Off')
%%%%%%% error handling
%if 0
catch
try, if nogui==0
for i=1:length(obj.enable), set(obj.children(i),'Enable',obj.enable{i}); end
set(h(1),'String','Apply',...
'ToolTip','Starts the computation - be patient.',...
'Callback','jrp compute')
setptr([hCRP,hCtrl],'arrow')
end, end
z=whos;x=lasterr;y=lastwarn;in=varargin{1};
print_error('jrp',z,x,y,in,mflag,action)
try, set(0,props.root), end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%