%zatěžovací charakteristiky

%% 01 výtahová
figure(1)                       %zolož okno pro graf
om = -1.8:0.1:1.8;              %definice řady úhlových rychlostí (rad/s)
K = 3;                          %konstanta určující zatížení
M = ones(size(om)).*K;          %výtahová = konstatní moment M = k, vynásobíme vektorem jedniček - ones o velikosti vektoru om, aby byly oba vektory stejného rozměru
P = M.*om;                      %výkon na hřídeli je P = M.om
plot(om, M, 'r', om, P, 'b','LineWidth',2.3)    %vykresli graf M(om) červeně, P(om) modře a zvětši tloušťku čar
xlabel('\omega')                %popisky os
ylabel('M, P')
title('Výtohová')               %nadpis

SetAxisParam([1.04 1],'{\itM = K}_1','{\itP = K}_1\omega') %funkce pro změnu osového kříže a jeho popisek  -- je nadefinována dole
% set(gcf,'Position',[1300 750 350 350])

%% hoblovková
figure(2)
om = -1.8:0.01:1.8;
K = 3;
M = K.*sign(om);          %výtahová = konstatní moment
P = M.*om;
plot(om, M, 'r', om, P, 'b','LineWidth',2.3)
xlabel('\omega')
ylabel('M, P')
title('Hoblovková')

SetAxisParam([1.2 1],'{\itM = K}_2{sign}\omega','{\itP = K}_2|\omega|')

% set(gcf,'Position',[1300 750 350 350])

%% kalandrová
figure(3)
om = -1.8:0.01:1.8;
K = 3;
M = K.*om;          %výtahová = konstatní moment
P = M.*om;
plot(om, M, 'r', om, P, 'b','LineWidth',2.3)
xlabel('\omega')
ylabel('M, P')
title('Kalandrová')


SetAxisParam([1.04 1],'{\itM = K}_3\omega','{\itP = K}_3\omega^2')
% set(gcf,'Position',[1300 750 350 350])

%% ventilátorová
figure(4)
om = -1.8:0.01:1.8;
K = 3;
M = K.*om.^2 .* sign(om);          %výtahová = konstatní moment
P = M.*om;
plot(om, M, 'r', om, P, 'b','LineWidth',2.3)
xlabel('\omega')
ylabel('M, P')
title('Ventilátorová')

% set(gcf,'Position',[1300 750 450 335])
SetAxisParam([1.35 0.5],['{\itM = K}_4\omega^2  ... \omega > 0' newline...
    '{\itM = -K}_4\omega^2 ... \omega < 0'],'{\itP = K}_4|\omega^3|')

%% navíječková
figure(5)
om = -1.8:0.1:1.8;
K = 30;
M = K./om;          %výtahová = konstatní moment
P = M.*om;
plot(om, M, 'r', om, P, 'b','LineWidth',2.3)
xlabel('\omega')
ylabel('M, P')
title('Navíječková')


% set(gcf,'Position',[1300 750 350 350])
SetAxisParam([1.04 1],'{\itM = K}_5/\omega','{\itP = K}_5')


%%
function [] = SetAxisParam(legendPosition,legendTxt1,legendTxt2)
axis square
LineWidth = 1.2;
set(gca,'XAxisLocation','origin','YAxisLocation','origin',...   %osový kříž uprostřed
    'XTickLabel',[],'YTickLabel',[],...                         %bez popisek
    'Box','off','FontSize',12,'LineWidth',LineWidth)            %bez rámečku, větší font, a čára

XLim = get(gca,'XLim'); %šipky os
YLim = get(gca,'YLim');
arrowOpen = 0.03;       
line([0.82 1].*XLim(2),[arrowOpen 0].*YLim(2),'LineWidth',LineWidth,'Color','k')
line([0.82 1].*XLim(2),[-arrowOpen 0].*YLim(2),'LineWidth',LineWidth,'Color','k')
line([arrowOpen 0].*XLim(2),[0.82 1].*YLim(2),'LineWidth',LineWidth,'Color','k')
line([-arrowOpen 0].*XLim(2),[0.82 1].*YLim(2),'LineWidth',LineWidth,'Color','k')

h = legend(legendTxt1,legendTxt2,'Location','southeast');
pos = get(h,'Position');
set(h,'FontSize',12,'Position',[pos(1)*legendPosition(1) pos(2)*legendPosition(2) pos(3:4)])

end