det_intervals

time_array = d.antenna0.servo.utc;
time_array = d.antenna0.servo.utcFast;

0001 function d= det_intervals(d)
0002 
0003 %time_array = d.antenna0.servo.utc;
0004 %time_array = d.antenna0.servo.utcFast;
0005 
0006 time_array1 = (d.antenna0.roach1.utc-d.antenna0.roach1.utc(1))*24*3600;
0007 
0008 for i=1:(size(time_array1)/2-1)
0009       inter1(i) = time_array1(i+1)-time_array1(i);
0010 end;
0011 
0012 time_array2 = (d.antenna0.roach2.utc-d.antenna0.roach2.utc(1))*24*3600;
0013 
0014 for i=1:(size(time_array2)/2-1)
0015       inter2(i) = time_array2(i+1)-time_array2(i);
0016 end;
0017 
0018 figure;
0019 subplot(2,1,1)
0020 plot(inter1-mean(inter1));
0021 title ('Roach 1 timing interval offset (t(i+1)-t(i)) - mean')
0022 ylabel('offset from mean interval (s)');
0023 xlabel('sample number');
0024 
0025 subplot(2,1,2)
0026 
0027 plot(inter2-mean(inter2));
0028 title ('Roach 2 timing interval offset (t(i+1)-t(i)) - mean')
0029 ylabel('offset from mean interval (s)');
0030 xlabel('sample number');
0031 %ylim([-0.000001 0.000001]);
0032 
0033 figure
0034 %plot(time_array2(167000:169000)-time_array1(167000:169000));
0035 plot(time_array2-time_array1);
0036 %plot(time_array2(182000:184000)-time_array1(182000:184000));
0037 grid on
0038 title ('Difference between d.antenna0.roach1.utc and d.antenna0.roach2.utc')
0039 ylabel('Difference (s)');
0040 xlabel('sample number');
0041 
0042 figure
0043 subplot(3,1,1)
0044 plot(time_array2-time_array1);
0045 %plot(time_array2(167000:170000)-time_array1(167000:170000));
0046 %plot(time_array2(182000:184000)-time_array1(182000:184000));
0047 title ('Difference between d.antenna0.roach1.utc and d.antenna0.roach2.utc')
0048 ylabel('Difference (s)');
0049 xlabel('sample number');
0050 grid on;
0051 subplot(3,1,2)
0052 plot(mean(d.antenna0.roach1.RR(:,:),2))
0053 %plot(mean(d.antenna0.roach1.RR(167000:170000,:),2))
0054 %plot(mean(d.antenna0.roach1.RR(182000:184000,:),2))
0055 title ('Roach 1 RR')
0056 ylabel('Roach1 (RR)');
0057 xlabel('sample number');
0058 ylim([4e4 4.4e4])
0059 subplot(3,1,3)
0060 plot(mean(d.antenna0.roach2.RR(:,:),2))
0061 %plot(mean(d.antenna0.roach2.RR(167000:170000,:),2))
0062 %plot(mean(d.antenna0.roach2.RR(182000:184000,:),2))
0063 title ('Roach 2 RR')
0064 ylabel('Roach2  (RR)');
0065 xlabel('sample number');
0066 ylim([3.4e4 3.8e4])
0067 grid on
0068 
0069