normalisation issue solved

phasesynchro.m

@@ -64,6 +64,9 @@ function CPRout = phasesynchro(varargin);
 % $Revision$
 %
 % $Log$
+% Revision 5.4  2009/03/24 08:34:24  marwan
+% copyright address changed
+%
 % Revision 5.3  2008/07/01 11:36:05  marwan
 % bug in embedding dimension fixed
 %
@@ -277,6 +280,9 @@ RR(find(isnan(RR)))=0; RR_(find(isnan(RR)))=0;
 rr1 = normalize(RR(w+1:end));
 rr2 = normalize(RR_(w+1:end));
 
+rr1a = (RR(w+1:end));
+rr2a = (RR_(w+1:end));
+
 
 
 % ACF of RR-tau
@@ -293,8 +299,8 @@ CPR = C(1,2);
 
 % plot
 if nogui ~= 2
-    plot(0:i3,rr1(1:i3+1),':',0:i3,rr2(1:i3+1),':'), hold on
-    plot(i3:w,rr1(i3+1:end),i3:w,rr2(i3+1:end)), hold off
+    plot(0:i3,rr1a(1:i3+1),':',0:i3,rr2a(1:i3+1),':'), hold on
+    plot(i3:w,rr1a(i3+1:end),i3:w,rr2a(i3+1:end)), hold off
     title(['CPR: ',sprintf('%2.2f',CPR)], 'fontw','bold')
     xlabel('Lag'), ylabel('RR_{\tau}')
 end