reduce_on_off

0001 clear
0002 figure
0003 d = read_arc('26-jan-2010:02:13:00','26-jan-2010:02:36:00'); %TauA
0004 %d = read_arc('26-jan-2010:02:43:00','26-jan-2010:03:00:00'); %CasA
0005  d_on = framecut(d, d.array.frame.features==1);
0006  d_off = framecut(d, d.array.frame.features==2);
0007  
0008  d_on = reshapeRx(d_on);
0009 d_on = reshapeAz(d_on);
0010 d_on = interpAz(d_on);
0011 
0012 d_off = reshapeRx(d_off);
0013 d_off = reshapeAz(d_off);
0014 d_off = interpAz(d_off);
0015 
0016 plot(d_on.antenna0.receiver.utc, d_on.antenna0.receiver.dataNew(:,1),'r');
0017 hold on;
0018 plot(d_off.antenna0.receiver.utc, d_off.antenna0.receiver.dataNew(:,1),'b');
0019 
0020 %Intensity averages
0021 mean_on=mean(d_on.antenna0.receiver.dataNew(:,1))
0022 mean_off=mean(d_off.antenna0.receiver.dataNew(:,1))
0023 error_on = rms(d_on.antenna0.receiver.dataNew(:,1))/(sqrt(length(d_on.antenna0.receiver.dataNew(:,1))))
0024 error_off = rms(d_off.antenna0.receiver.dataNew(:,1))/(sqrt(length(d_off.antenna0.receiver.dataNew(:,1))))
0025 
0026 mean_on_a=mean(d_on.antenna0.receiver.dataNew(:,6))
0027 mean_off_a=mean(d_off.antenna0.receiver.dataNew(:,6))
0028 error_on_a = rms(d_on.antenna0.receiver.dataNew(:,6))/(sqrt(length(d_on.antenna0.receiver.dataNew(:,6))))
0029 error_off_a = rms(d_off.antenna0.receiver.dataNew(:,6))/(sqrt(length(d_off.antenna0.receiver.dataNew(:,6))))
0030 
0031 
0032 
0033 figure
0034 
0035 plot(d_on.antenna0.receiver.utc, d_on.antenna0.receiver.dataNew(:,2),'r');
0036 hold on;
0037 plot(d_off.antenna0.receiver.utc, d_off.antenna0.receiver.dataNew(:,2),'b');
0038 mean_on_pol=mean(d_on.antenna0.receiver.dataNew(:,2))
0039 mean_off_pol=mean(d_off.antenna0.receiver.dataNew(:,2))
0040 error_on_pol = rms(d_on.antenna0.receiver.dataNew(:,2))/(sqrt(length(d_on.antenna0.receiver.dataNew(:,2))))
0041 error_off_pol = rms(d_off.antenna0.receiver.dataNew(:,2))/(sqrt(length(d_off.antenna0.receiver.dataNew(:,2))))
0042 mean_on_pol-mean_off_pol