% Figure 2. 50 A Type B and 32 A type C MCB trip-time versus time 
% characteristics and associated PEC trip characteristic (yellow curve)

clear
FontSize = 28;
figure(2);
clf

% MCBs Time-Current Characteristics
I_b  = [1.45    1.45    1.5     1.75    2   2.5 3   3.5 4   4.5 5   5]*50;
Time = [10000   4000    2000    400     175 70  40  27  20  15  12  0.1];
semilogy(I_b,Time,'LineWidth',2)
hold on
I_c  = [1.45    1.45    1.5     1.75    2   2.5 3   3.5 4   4.5 5   6   7       8   10  10] * 32;
Time = [10000   4000    2000    400     175 70  40  27  20  15  12  9   7.25   6.25 5   0.1];
semilogy(I_c,Time,'LineWidth',2)

% PEC Coordination
I_PEC = [[1.45    1.45    1.5     1.75    2   2.5 3   3.5 4   4.5 5   5 5 5 0 0 0]*50 + [0 0 0 0 0 0 0 0 0 0 0   0   0 0       8   10  10] * 32] * 1.1 
Time = [10000   4000    2000    400     175 70  40  27  20  15  12  9   7.25   7 6.75 5.5   0.1];
semilogy(I_PEC,Time,'LineWidth',4)

hold off

legend({'50A Type B MCB','32A Type C MCB','PEC', },'FontSize',FontSize,'Location','NorthEast')
xlabel('AC RMS Current [A]','FontSize', FontSize)
ylabel('Time [s]','FontSize', FontSize)

FH = gca;
FH.FontSize = FontSize - 2;
FH.FontName = 'Times New Roman';
grid