filename1 = 'fin125_period250_250_a.xlsx';
filename2 = 'fin125_period250_250_b.xlsx';
filename3 = 'fin125_period250_250_c.xlsx';
filename4 = 'fin125_period250_500_a.xlsx';
filename5 = 'fin125_period250_500_b.xlsx';
filename6 = 'fin125_period250_500_c.xlsx';
filename7 = 'fin125_period250_1000_a.xlsx';
filename8 = 'fin125_period250_1000_b.xlsx';
filename9 = 'fin125_period250_1000_c.xlsx';
filename10 = 'fin125_period250_2000_a.xlsx';
filename11 = 'fin125_period250_2000_b.xlsx';
filename12 = 'fin125_period250_2000_c.xlsx';
filename13 = 'taper_only_a.xlsx';
filename14 = 'taper_only_b.xlsx';
filename15 = 'taper_only_c.xlsx';
sheet1 = 1;
sheet2 = 2;
xlRange = 'A2:B10003';

M1 = min(xlsread(filename1, sheet2, xlRange));
M2 = min(xlsread(filename2, sheet2, xlRange));
M3 = min(xlsread(filename3, sheet2, xlRange));
M4 = min(xlsread(filename4, sheet2, xlRange));
M5 = min(xlsread(filename5, sheet2, xlRange));
M6 = min(xlsread(filename6, sheet2, xlRange));
M7 = min(xlsread(filename7, sheet2, xlRange));
M8 = min(xlsread(filename8, sheet2, xlRange));
M9 = min(xlsread(filename9, sheet2, xlRange));
M10 = min(xlsread(filename10, sheet2, xlRange));
M11 = min(xlsread(filename11, sheet2, xlRange));
M12 = min(xlsread(filename12, sheet2, xlRange));
M13 = min(xlsread(filename13, sheet2, xlRange));
M14 = min(xlsread(filename14, sheet2, xlRange));
M15 = min(xlsread(filename15, sheet2, xlRange));

M13 = cat(1, M1(1,2), M2(1,2), M3(1,2), M4(1,2), M5(1,2), M6(1,2), M7(1,2), M8(1,2), M9(1,2), M10(1,2), M11(1,2), M12(1,2), M13(1,2), M14(1,2), M15(1,2));

%{
plot(X1,Y1,'g');
hold on

plot(X2,Y2,'r');
hold on

plot(X3,Y3,'k');
hold on

plot(X4,Y4,'m');
ylabel ('Total loss (dB)')
hold on
    
yyaxis right
    %ylim ([5 15])
    plot(c,b, 'b');
    ylabel ('dB/mm');
    hold on
    
title('SiWG023-C3: Horizontal-slot, strip-waveguide loss , GC design: 1, devices: c');
    xlabel('Wavelength (nm)');
    ylabel('Loss (dB/mm)');
    legend('WG length = 250um','WG length = 500um', 'WG length = 1000um', 'WG length = 2000um', 'Loss per mm (dB/mm)');
    hold on

    
disp(b);
disp((X1(5)));
disp(X2(5));
disp(X3(5));
%}