0001 fid = fopen('Feb2012Op.txt'); 0002 data = fscanf(fid, '%f %f %f %f %f %f %f', [7 inf]); 0003 fclose(fid); 0004 0005 data = data'; 0006 0007 SizeD = size(data); 0008 0009 utc = data(:,1); 0010 temp = data(:,2); 0011 relH = data(:,3); 0012 meanRelH = mean(relH); 0013 pipeOp = data(:,4); 0014 meanPipeOp = mean(pipeOp); 0015 myOp = data(:,5); 0016 meanMyOp = mean(myOp); 0017 skydips = data(:,6); 0018 meanSkydips = mean(skydips); 0019 noiseD = data(:,7); 0020 meanNoiseD = mean(noiseD); 0021 0022 time = (utc - utc(1)); 0023 0024 varRelH = relH - meanRelH; 0025 varPipeOp = pipeOp - meanPipeOp; 0026 varMyOp = myOp - meanMyOp; 0027 varSkydips = skydips - meanSkydips; 0028 varNoiseD = noiseD - meanNoiseD; 0029 0030 figure(1) 0031 [A, H1] = plotyy(time, myOp, time, temp); 0032 set(H1,'LineStyle', 'None','Marker', 'o') 0033 hold on 0034 plot(time, pipeOp./3, 'r o'); 0035 hold on 0036 plot(time, skydips, 'k o') 0037 hold off 0038 legend('MyTheoOp', 'PipelineOp/3', 'SkyDips'); 0039 xlabel('Time (days)'); 0040 set(get(A(1), 'YLabel'), 'String', 'Opacity') 0041 set(get(A(2), 'YLabel'), 'String', 'Atmospheric Temperature (*C)'); 0042 title('Opacities Against Temp: 01/02/12') 0043 0044 figure(2) 0045 [A, H1] = plotyy(time, myOp, time, relH); 0046 set(H1, 'LineStyle', 'None','Marker', 'o') 0047 hold on 0048 plot(time, pipeOp./3, 'r o'); 0049 hold on 0050 plot(time, skydips, 'k o') 0051 hold off 0052 legend('MyTheoOp', 'PipelineOp/3', 'SkyDips'); 0053 xlabel('Time (days)'); 0054 set(get(A(1), 'YLabel'), 'String', 'Opacity') 0055 set(get(A(2), 'YLabel'), 'String', 'Relative Humidity (%)'); 0056 title('Opacities Against Humidity: 01/02/12') 0057 0058 figure(3) 0059 plot(time, myOp, 'b o', time, pipeOp, 'r o', time, skydips, 'g o', time, noiseD, 'k o'); 0060 legend('MyTheoOp', 'PipelineOp', 'SkyDips', 'NoiseD'); 0061 xlabel('Time (days)'); 0062 ylabel('Opacity') 0063 title('Opacities: 01/02/12') 0064 0065 figure(4) 0066 plot(time, varPipeOp + 0.001, 'r o') 0067 hold on 0068 plot(time, varMyOp, 'b o') 0069 hold on 0070 plot(time, varSkydips -0.001, 'g o') 0071 hold on 0072 plot(time, varNoiseD -0.002, 'k o') 0073 hold off 0074 legend('PipeTheoOp + 0.001', 'MyTheoOp', 'SkyDips -0.001', 'NosieD -0.002') 0075 title('Mean Variations in Opacities: 01/02/12') 0076 xlabel('Time (days)'); 0077 ylabel('Offset Opacity Mean Variations') 0078 0079