basicGenericPolarisation

 a few test scripts for CBass South

0001 % a few test scripts for CBass South
0002 
0003 %d = read_arcSouth('19-Mar-2012:13:40:20','19-Mar-2012:13:53:00');
0004 %d = read_arcSouth('19-Mar-2012:14:20:20','19-Mar-2012:14:18:00');
0005 %d = read_arcSouth('20-Mar-2012:17:27:00','20-Mar-2012:17:37:30'); %%lid off power supply
0006 d = read_arcSouth('20-Mar-2012:17:58:00','20-Mar-2012:18:01:00'); %%lid on power supply
0007 
0008 time=10; %time in minutes
0009 LLAveraged=mean(d.antenna0.roach1.LL(:,5:62),2);
0010 RRAveraged=mean(d.antenna0.roach1.RR(:,5:62),2);
0011 QAveraged=(mean(d.antenna0.roach1.Q(:,5:62),2));
0012 UAveraged=(mean(d.antenna0.roach1.U(:,5:62),2));
0013 LoadRAveraged=(mean(d.antenna0.roach1.load1(:,5:62),2));
0014 LoadLAveraged=(mean(d.antenna0.roach1.load2(:,5:62),2));
0015 
0016 LLdB=10*log10(LLAveraged);
0017 RRdB=10*log10(RRAveraged);
0018 
0019 LoadRdB=10*log10(LoadRAveraged);
0020 LoadLdB=10*log10(LoadLAveraged);
0021 close all
0022 noiseDiode=d.antenna0.roach1.switchstatus;
0023 
0024 % subplot(5,1,1)
0025 % plot(1440.*(d.antenna0.roach1.utc-d.antenna0.roach1.utc(1)),LLAveraged,'.')
0026 % hold
0027 % plot(1440.*(d.antenna0.roach1.utc-d.antenna0.roach1.utc(1)),RRAveraged,'.r')
0028 % plot(1440.*(d.antenna0.roach1.utc-d.antenna0.roach1.utc(1)),LoadRAveraged,'.c')
0029 % plot(1440.*(d.antenna0.roach1.utc-d.antenna0.roach1.utc(1)),LoadLAveraged,'.k')
0030 %
0031 % legend('LL','RR','LoadR','LoadL');
0032 % xlim([0 time])
0033 % %ylim([93 95])
0034 % grid on
0035 % title('Sky (300K) and Load')
0036 % ylabel('ADC')
0037 % xlabel('Time [m]')
0038 
0039 subplot(5,1,1)
0040 plot(1440./60.*(d.antenna0.roach1.utc-d.antenna0.roach1.utc(1)),10*log10(LLAveraged),'.')
0041 hold
0042 plot(1440./60.*(d.antenna0.roach1.utc-d.antenna0.roach1.utc(1)),10*log10(RRAveraged),'.r')
0043 legend('LL','RR');
0044 xlim([0 time])
0045 %ylim([93 95])
0046 grid on
0047 title('Sky (300K)')
0048 ylabel('Power [dB]')
0049 xlabel('Time [m]')
0050 %ylim([230 290])
0051 
0052 subplot(5,1,2)
0053 plot(1440.*(d.antenna0.roach1.utc-d.antenna0.roach1.utc(1)),10*log10(LoadRAveraged),'.c')
0054 hold
0055 plot(1440.*(d.antenna0.roach1.utc-d.antenna0.roach1.utc(1)),10*log10(LoadLAveraged),'.k')
0056 
0057 legend('LoadR','LoadL');
0058 xlim([0 time])
0059 %ylim([40 50])
0060 %ylim([93 95])
0061 grid on
0062 title('Load Power')
0063 ylabel('Power [dB]')
0064 xlabel('Time [m]')
0065 
0066 subplot(5,1,3)
0067 plot(1440*(d.antenna0.thermal.utc-d.antenna0.roach1.utc(1)),d.antenna0.thermal.ccTemperatureLoad,'r.')
0068 hold
0069 plot(1440*(d.array.frame.utc-d.antenna0.roach1.utc(1)),d.antenna0.thermal.lsTemperatureSensors(:,4),'b.')
0070 xlim([0 time])
0071 ylim([0 20])
0072 title('Temperatures')
0073 ylabel('Kelvin')
0074 xlabel('Time [m]')
0075 legend('CryoCon Bobbin','LakeShore SecondStage')
0076 grid on
0077 
0078 subplot(5,1,4)
0079 plot(1440.*(d.antenna0.roach1.utc-d.antenna0.roach1.utc(1)),noiseDiode,'.k')
0080 title('NoiseDiode')
0081 ylabel('Value')
0082 xlabel('Time [m]')
0083 xlim([0 time])
0084 ylim([-1 2])
0085 
0086 subplot(5,1,5)
0087 Q=d.antenna0.roach1.Qtime;
0088 U=d.antenna0.roach1.Utime;
0089 plot(1440.*(d.array.frame.utc-d.antenna0.roach1.utc(1)),Q(:,31),'.k')
0090 hold
0091 plot(1440.*(d.array.frame.utc-d.antenna0.roach1.utc(1)),U(:,31),'.r')
0092 legend('Q','U')
0093 
0094 title('Polarisation')
0095 ylabel('Value')
0096 xlabel('Time [m]')
0097 xlim([0 time])
0098 
0099 figure
0100 
0101 Hs=spectrum.periodogram;  
0102 Fs=100;
0103 psd(Hs,LLAveraged,'Fs',Fs)
0104 
0105 
0106 
0107 %legend('CryoCon Bobbin','LakeShore SecondStage')