thermalNoisePlots

sensititivity

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