demomin

 minutes

0001 
0002 % minutes
0003   start=[1990 1 1 5 0 0];    %Gregorian start [yyyy mm dd hh mi sc]
0004   stop=[1990 1 1 5 20 0];
0005   jd=[julian(start):1/5000:julian(stop)];
0006 
0007 % synthesize velocity data
0008   u=sin(.1*jd(:)*500).^2-.5;
0009   v=cos(.1*jd(:)*500);
0010 
0011 % convention: store velocity time series as complex vector:
0012   w=u+i*v;
0013 
0014   h=timeplt(jd,[u v abs(w) w],[1 1 2 3]);
0015   title('Demo of TIMEPLT and STACKLBL:  Minutes')
0016  
0017 % use STACKLBL to label each stack plot panel with title and units:
0018 
0019   stacklbl(h(1),'East + North velocity','m/s');
0020   stacklbl(h(2),'Speed','m/s');
0021   stacklbl(h(3),'Velocity Sticks','m/s');
0022