
close all;

% vykresleni
figure      % zalozi novy obrazek
plot(t_v,v,'r',t_v,a_t,'b')   % vykresleni prUbehU
axis([0 T_sum 5/4*min(a_t) 5/4*max(v)])
title('Prubeh rychlosti a zrychleni stroje v case')
xlabel('t [s]'), ylabel('v [km/h], a [m/s^2]')
set(gca,'FontSize',20)
grid on

% vykresleni
figure      % zalozi novy obrazek
plot(t_v,v,'r',t_v,M_pr*0.01,'b')   % vykresleni prUbehU
axis([0 T_sum 5/4*min(M_pr)*0.01 5/4*max(v)])
title('Prubeh rychlosti stroje a krouticiho momentu na vystupu prevodovky v case')
xlabel('t [s]'), ylabel('v [km/h], M_pr*0.02 [Nm]')
set(gca,'FontSize',20)
grid on

P_plus=[];
P_minus=[];
for i=1:length(t_v)
    if P_pr(i)<0
        P_minusi=P_pr(i);
        P_plusi=0;
    else
        P_plusi=P_pr(i);
        P_minusi=0;
    end
    P_plus=[P_plus;P_plusi];
    P_minus=[P_minus;P_minusi];
end
P_2pr=[P_plus P_minus];
figure
plot(t_v,P_2pr)
area(t_v,P_2pr)
title('Prubeh vykonu na vystupu prevodovky v case')
xlabel('t [s]'), ylabel('P_p [kW]')
set(gca,'FontSize',20)
grid on