1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85
| clear all;
clc;
close all;
file1='20120731-full-350nm-slit100ShTi0.1Gain5.asc';
S1=load(file1);
file2='20120731-full-450nm-slit100ShTi0.1Gain5.asc';
S2=load(file2);
file3='20120731-full-550nm-slit100ShTi0.1Gain5.asc';
S3=load(file3);
file4='20120731-full-650nm-slit100ShTi0.1Gain5.asc';
S4=load(file4);
file5='20120731-full-750nm-slit100ShTi0.1Gain5.asc';
S5=load(file5);
file6='20120731-full-850nm-slit100ShTi0.1Gain5.asc';
S6=load(file6);
l1=size(S1,1);
l2=size(S2,1);
l3=size(S3,1);
l4=size(S4,1);
l5=size(S5,1);
l6=size(S6,1);
c1=size(S1,2);
c2=size(S2,2);
c3=size(S3,2);
c4=size(S4,2);
c5=size(S5,2);
c6=size(S6,2);
Q1= S1(129:l1-133,1);
Q2= S2(129:l2-133,1);
Q3= S3(129:l3-133,1);
Q4= S4(129:l4-133,1);
Q5= S5(129:l5-133,1);
Q6= S6(129:l6-133,1);
R1 = sum(S1(129:l1-133,150:c1-150),2);
R2 = sum(S2(129:l2-133,150:c2-150),2);
R3 = sum(S3(129:l3-133,150:c3-150),2);
R4 = sum(S4(129:l4-133,150:c4-150),2);
R5 = sum(S5(129:l5-133,150:c5-150),2);
R6 = sum(S6(129:l6-133,150:c6-150),2);
% figure(1);
% subplot(2,3,1);plot(Q1,R1);title('350nm');
% subplot(2,3,2);plot(Q2,R2);title('450nm');
% subplot(2,3,3);plot(Q3,R3);title('550nm');
% subplot(2,3,4);plot(Q4,R4);title('650nm');
% subplot(2,3,5);plot(Q5,R5);title('750nm');
% subplot(2,3,6);plot(Q6,R6);title('850nm');
%
% saveName = ('Lines');
% saveas(gcf, saveName, 'jpg');
% saveas(gcf, saveName, 'fig');
%
% figure(2);
%
% plot(Q1,R1,'b');
% hold on
% plot(Q2,R2,'r');
% hold on
% plot(Q3,R3,'y');
% hold on
% plot(Q4,R4,'g');
% hold on
% plot(Q5,R5,'c');
% hold on
% plot(Q6,R6,'k');
% title('Full spectra at 100% - at the top (21,5mm from the bottom)');
% xlabel('Wavelength in nm');
%
% hold off
% hold off
% hold off
% hold off
% hold off
%
% saveName = ('Full spectra at 100% - at the top (21,5mm from the bottom)');
% saveas(gcf, saveName, 'jpg');
% saveas(gcf, saveName, 'fig'); |