thermalNoise

0001 subplot(211)
0002 plot(d.antenna0.roach1.LL)
0003 subplot(212)
0004 plot(d.antenna0.roach1.switchstatus)
0005 
0006 
0007 %sensititivity
0008 d=read_arcSouth('15-Jan-2013:12:30:00','15-Jan-2013:13:00:00') % @ 12K both sky and load
0009 d=read_arcSouth('15-Jan-2013:14:20:00','15-Jan-2013:14:50:00') % @ 16K both sky and load
0010 d=read_arcSouth('15-Jan-2013:15:00:00','15-Jan-2013:15:20:00') % @ 18K both sky and load
0011 %d=read_arcSouth('16-Jan-2013:20:00:00','16-Jan-2013:20:10:00') % @ 14K both sky and load
0012 %d=read_arcSouth('17-Jan-2013:01:00:00','17-Jan-2013:01:20:00') % @ 24K both sky and load
0013 d=read_arcSouth('19-Jan-2013:12:00:00','19-Jan-2013:12:20:00') % @ 26K both sky and load
0014 %16-Jan-2013:18:13:16      17-Jan-2013:13:39:08
0015 %19-Jan-2013:04:19:05      19-Jan-2013:17:32:17
0016 BU2Kelvin=2700; %approximate conversion from BU to Kelvin from noise steps and noise diode
0017 %d=read_arcSouth('09-Dec-2012:21:00:00','09-Dec-2012:21:30:00') % @ sky @10.5K and load at 18K
0018 %d=read_arcSouth('09-Dec-2012:20:00:00','09-Dec-2012:20:30:00') % @ sky @10.5K and load at 15K
0019 %BU2Kelvin=1400; %approximate conversion from BU to Kelvin from noise steps and noise diode
0020 
0021 %balanced load
0022 LL1=d.antenna0.roach1.LL-d.antenna0.roach1.load2;
0023 RR1=d.antenna0.roach1.RR-d.antenna0.roach1.load1;
0024 LL2=d.antenna0.roach2.LL-d.antenna0.roach2.load2;
0025 RR2=d.antenna0.roach2.RR-d.antenna0.roach2.load1;
0026 
0027 
0028 close all
0029 
0030 
0031 % to combine both use this
0032 LLcomb = [LL1 LL2];
0033 RRcomb = [RR1 RR2];
0034 
0035 %LLcomb = [LL1];
0036 %RRcomb = [RR1];
0037 
0038 LLmean=mean(LLcomb,2);
0039 RRmean=mean(RRcomb,2);
0040 
0041 indextt=[1:12000]; % index of the 100Hz Samples
0042 indextt2=[1000:1200]; % index of the downsampled 1Hz Samples
0043 
0044 subplot(211)
0045 plot(RRmean(indextt))
0046 title('RR')
0047 grid on
0048 
0049 subplot(212)
0050 plot(LLmean(indextt))
0051 title('LL')
0052 grid on
0053 
0054 
0055 
0056 
0057 
0058 
0059 ndOnindex=logical(bitand(d.antenna0.roach1.switchstatus,1));
0060 ndOffindex=logical(~(bitand(d.antenna0.roach1.switchstatus,1)));
0061 
0062 figure
0063 plot(RRmean(ndOnindex,:))
0064 
0065 LLmeanK=LLmean./BU2Kelvin;
0066 RRmeanK=RRmean./BU2Kelvin;
0067 
0068 binLLmean=bin(LLmeanK,100);
0069 binRRmean=bin(RRmeanK,100);
0070 
0071 subplot(211)
0072 plot(binRRmean(1:1.750e3))
0073 title('RR')
0074 grid on
0075 
0076 subplot(212)
0077 plot(binLLmean(1:1.750e3))
0078 title('LL')
0079 grid on
0080 
0081 noiseOnL=mean(LLmeanK(ndOnindex));
0082 noiseOffL=mean(LLmeanK(ndOffindex));
0083 
0084 noiseOnR=mean(RRmeanK(ndOnindex));
0085 noiseOffR=mean(RRmeanK(ndOffindex));
0086 
0087 figure
0088 subplot(211)
0089 plot(binRRmean(1:1.650e3))
0090 title('RR')
0091 grid on
0092 
0093 subplot(212)
0094 plot(binLLmean(1:1.650e3))
0095 title('LL')
0096 grid on
0097 
0098 
0099 figure
0100 plot(binRRmean(indextt2))
0101 hold all
0102 plot(binLLmean(indextt2))
0103 
0104 
0105 stdR= std(binRRmean(indextt2))
0106 stdL = std(binLLmean(indextt2))
0107 
0108 
0109 freq=generatePowerSpectrav2(utc,LLmean(indextt))
0110 
0111 loglog(freq(:,2),freq(:,1))
0112 
0113 1.8
0114 
0115 
0116 stdRa=std(RRmeanK(indextt))
0117 stdLa=std(LLmeanK(indextt))
0118 
0119 subplot(211)
0120 plot(1000.*(RRmeanK(indextt)-mean(RRmeanK(indextt))))
0121 title(['RR :' texlabel('sigma') '=' num2str(1000.*stdRa) ' mK'])
0122 ylabel('mK')
0123 xlabel('Sample 100Hz')
0124 grid on
0125 
0126 subplot(212)
0127 plot(1000.*(LLmeanK(indextt)-mean(LLmeanK(indextt))))
0128 title(['LL :' texlabel('sigma') '=' num2str(1000.*stdLa) ' mK'])
0129 ylabel('mK')
0130 xlabel('Sample 100Hz')
0131 grid on
0132 
0133 
0134 print('-depsc','ThermalNoised_26K.eps')
0135 
0136 
0137 
0138 stdRa=std(RRmeanK(indextt))
0139 stdLa=std(LLmeanK(indextt))