plotSavedStats

% Load Saved File

0001 %% Load Saved File
0002 date_start = '01-Oct-2012';
0003 date_end = '19-Jan-2014';
0004 load( ['skyStareNoiseStats_' date_start '_' date_end '.mat'] );
0005 
0006 % alpha = Ivals(0)
0007 % f_knee = Ivals(1)
0008 % sigma_w = Ivals(2)
0009 %
0010 % 95% statistical confidence levels in lower and upper
0011 % RMS of fits in rms
0012 % Start MJD of sky stare in MJD
0013 %
0014 
0015 %% Convert Date
0016 [year month day] = mjd2date(MJD);
0017 dateno = datenum( year, month, day );
0018 
0019 %% Side By Side Plotting
0020 % KNEE FREQ
0021 figure
0022 subplot(131)
0023 semilogy( dateno(Ivals(:,2)<10.), Ivals(Ivals(:,2)<10.,2), '.k', 'markerSize',10 )
0024 datetick('x','mmm-YY','keepticks')%,'keeplimits')
0025 xlabel( 'Date' )
0026 ylabel( 'f_{knee}' )
0027 title( 'I' )
0028 grid()
0029 
0030 subplot(132)
0031 semilogy( dateno(Qvals(:,2)<10.), Qvals(Qvals(:,2)<10.,2), '.k', 'markerSize',10 )
0032 datetick('x','mmm-YY','keepticks')%,'keeplimits')
0033 xlabel( 'Date' )
0034 title( 'Q' )
0035 grid()
0036 
0037 subplot(133)
0038 semilogy( dateno(Uvals(:,2)<10.), Uvals(Uvals(:,2)<10.,2), '.k', 'markerSize',10 )
0039 datetick('x','mmm-YY','keepticks')%,'keeplimits')
0040 xlabel( 'Date' )
0041 title( 'U' )
0042 grid()
0043 
0044 % ALPHA
0045 figure
0046 subplot(131)
0047 plot( dateno(Ivals(:,1)<10.), Ivals(Ivals(:,1)<10.,1), '.k', 'markerSize',10 )
0048 datetick('x','mmm-YY','keepticks')%,'keeplimits')
0049 xlabel( 'Date' )
0050 ylabel( '\alpha' )
0051 title( 'I' )
0052 grid()
0053 
0054 subplot(132)
0055 plot( dateno(Qvals(:,1)<10.), Qvals(Qvals(:,1)<10.,1), '.k', 'markerSize',10 )
0056 datetick('x','mmm-YY','keepticks')%,'keeplimits')
0057 xlabel( 'Date' )
0058 title( 'Q' )
0059 grid()
0060 
0061 subplot(133)
0062 plot( dateno(Uvals(:,1)<10.), Uvals(Uvals(:,1)<10.,1), '.k', 'markerSize',10 )
0063 datetick('x','mmm-YY','keepticks')%,'keeplimits')
0064 xlabel( 'Date' )
0065 title( 'U' )
0066 grid()
0067 
0068 % SIGMA_W
0069 figure
0070 subplot(131)
0071 semilogy( dateno(Ivals(:,3)<10.), Ivals(Ivals(:,3)<10.,3), '.k', 'markerSize',10 )
0072 datetick('x','mmm-YY','keepticks')%,'keeplimits')
0073 xlabel( 'Date' )
0074 ylabel( '\sigma_w' )
0075 title( 'I' )
0076 grid()
0077 
0078 subplot(132)
0079 semilogy( dateno(Qvals(:,3)<10.), Qvals(Qvals(:,3)<10.,3), '.k', 'markerSize',10 )
0080 datetick('x','mmm-YY','keepticks')%,'keeplimits')
0081 xlabel( 'Date' )
0082 title( 'Q' )
0083 grid()
0084 
0085 subplot(133)
0086 semilogy( dateno(Uvals(:,3)<10.), Uvals(Uvals(:,3)<10.,3), '.k', 'markerSize',10 )
0087 datetick('x','mmm-YY','keepticks')%,'keeplimits')
0088 xlabel( 'Date' )
0089 title( 'U' )
0090 grid()
0091 
0092 
0093 %% On Top Plotting
0094 figure
0095 semilogy( dateno(Ivals(:,2)<10.), Ivals(Ivals(:,2)<10.,2), '.r', 'markerSize',10 )
0096 hold
0097 semilogy( dateno(Qvals(:,2)<10.)+1, Qvals(Qvals(:,2)<10.,2), '.g', 'markerSize',10 )
0098 semilogy( dateno(Uvals(:,2)<10.)+2, Uvals(Uvals(:,2)<10.,2), '.b', 'markerSize',10 )
0099 
0100 datetick('x','mmm-YY','keepticks')%,'keeplimits')
0101 xlabel( 'Date' )
0102 ylabel( 'f_{knee}' )
0103 title( 'I-red, Q-green, U-blue' )
0104 grid()
0105 
0106 % ALPHA
0107 figure
0108 plot( dateno(Ivals(:,1)<10.), Ivals(Ivals(:,1)<10.,1), '.r', 'markerSize',10 )
0109 hold
0110 plot( dateno(Qvals(:,1)<10.)+1, Qvals(Qvals(:,1)<10.,1), '.g', 'markerSize',10 )
0111 plot( dateno(Uvals(:,1)<10.)+2, Uvals(Uvals(:,1)<10.,1), '.b', 'markerSize',10 )
0112 
0113 datetick('x','mmm-YY','keepticks')%,'keeplimits')
0114 xlabel( 'Date' )
0115 ylabel( '\alpha' )
0116 title( 'I-red, Q-green, U-blue' )
0117 grid()
0118 
0119 % SIGMA_W
0120 figure
0121 semilogy( dateno(Ivals(:,3)<10.), Ivals(Ivals(:,3)<10.,3), '.r', 'markerSize',10 )
0122 hold
0123 semilogy( dateno(Qvals(:,3)<10.)+1, Qvals(Qvals(:,3)<10.,3), '.g', 'markerSize',10 )
0124 semilogy( dateno(Uvals(:,3)<10.)+2, Uvals(Uvals(:,3)<10.,3), '.b', 'markerSize',10 )
0125 
0126 datetick('x','mmm-YY','keepticks')%,'keeplimits')
0127 xlabel( 'Date' )
0128 ylabel( '\sigma_w' )
0129 title( 'I-red, Q-green, U-blue' )
0130 grid()