ИЭ / 3 семестр / графики для курсача в матлабе / stbIk_i_coordinata

clc,close all;

tau=4*10^(-4);   %введи своё тау

Xphi=0:0.01:90;
Xpsi=90:0.01:120;
X1l=90:0.01:120;
X20=120:0.01:150;

i1phi=300*(1-Xphi+Xphi);

%exp(параша) - эту часть нигде не меняйте, потому что она у всех одинаковая

i1psi=-300*(1-exp(-1/tau * (10^(-4)-(120-Xpsi)/300000)));   %введи свою функцию тока i1psi 
i1l=300*exp(-1/tau * (10^(-4)-(120-X1l)/300000));   %введи свою функцию тока i1l
i20=150*exp(-1/tau * (10^(-4)-(X20-120)/300000));   %введи свою функцию тока i20

figure;
hold on;
grid on;

plot(Xphi,i1phi,'-','LineWidth',1.3,'color',[255 255 0]/255);   %строим i1phi
plot([90 120],[300 300],'--','LineWidth',1.3,'color',[255 255 0]/255);
plot(Xpsi,i1psi,'--','LineWidth',1.3,'color',[0 0 255]/255);   %строим i1psi
plot(X1l,i1l,'-','LineWidth',1.3,'color',[255 0 0]/255);   %строим i1l
plot(X20,i20,'-','LineWidth',1.3,'color',[0 255 0]/255);   %строим i20

xlim([0 170]);
ylim([min([min(i1psi) min(i1phi) min(i1l) min(i20)])-75 max([max(i1psi) max(i1phi) max(i1l) max(i20)])+100]);    %если не работает из-за этой строчки, то ставьте сами числа на рандоме и смотрите


ylabel('i(x), А');   %название оси Y

%Строим свою ось Х
x = annotation('arrow');  
x.Parent = gca;           
x.X = [0 165];          
x.Y = [0 0]; 
x.LineWidth  = 1;         
x.HeadWidth  = 8;
x.HeadLength = 8;

%засечки на оси Х
for i=0:3:15
    plot([10*i 10*i], [-3 3],'-','LineWidth',1.3,'color',[0 0 0]);
end
text(28.65, -7,'30', 'FontSize',12);
text(58.65, -7,'60', 'FontSize',12);
text(88.65, -7,'90', 'FontSize',12);
text(118.65, -7,'120', 'FontSize',12);
text(148.65, -7,'150', 'FontSize',12);

%название оси
text(160, 12,'х,км', 'FontSize',12);


%подибрайте значения, чтобы на оси Y были нужные вам засечки
set(gca,'XTick',0:30:160,'YTick',10*floor((min([min(i1psi) min(i1phi) min(i1l) min(i20)]))/10):20:max([max(i1psi) max(i1phi) max(i1l) max(i20)])+100); 

%строим серую штриховку
for i=1:8
    plot([10*i 10*i], [0 i1phi(1)],'-','LineWidth',1.3,'color',[0.7 0.7 0.7]);
end

for i=9:11
    plot([10*i 10*i], [0 i1l(1000*(i-9)+1)],'-','LineWidth',1.3,'color',[0.7 0.7 0.7]);
end
plot([90 90], [0 i1l(1)],'-','LineWidth',1.3,'color',[0.7 0.7 0.7]);

for i=13:15
    plot([10*i 10*i], [0 i20(1000*(i-12)+1)],'-','LineWidth',1.3,'color',[0.7 0.7 0.7]);
end

%строим вертикальную прямую по середине. Меняйте значение высоты, если надо
plot([120 120], [min([min(i1psi) min(i1phi) min(i1l) min(i20)])-75 -13],'-','LineWidth',1.3,'color',[0 0 0]);   
plot([120 120], [0 max([max(i1psi) max(i1phi) max(i1l) max(i20)])+100],'-','LineWidth',1.3,'color',[0 0 0]);
%строим вектора
x1shtrih = annotation('arrow');  
x1shtrih.Parent = gca;           
x1shtrih.X = [120 80];          
x1shtrih.Y = [min([min(i1psi) min(i1phi) min(i1l) min(i20)])-30 min([min(i1psi) min(i1phi) min(i1l) min(i20)])-30]; %меняйте координаты высоты стрелки, если надо  
x1shtrih.LineWidth  = 1;         
x1shtrih.HeadWidth  = 8;
x1shtrih.HeadLength = 8;

x2 = annotation('arrow');  
x2.Parent = gca;           
x2.X = [120 160];          
x2.Y = [min([min(i1psi) min(i1phi) min(i1l) min(i20)])-40 min([min(i1psi) min(i1phi) min(i1l) min(i20)])-40];   %меняйте координаты высоты стрелки, если надо   
x2.LineWidth  = 1;         
x2.HeadWidth  = 8;
x2.HeadLength = 8;

%надписи векторов
text(80, min([min(i1psi) min(i1phi) min(i1l) min(i20)])-15,'x_{1}''', 'FontSize',12);   
text(160, min([min(i1psi) min(i1phi) min(i1l) min(i20)])-30,'x_{2}', 'FontSize',12);

%засечки
for i=9:11
    plot([10*i 10*i], [min([min(i1psi) min(i1phi) min(i1l) min(i20)])-33 min([min(i1psi) min(i1phi) min(i1l) min(i20)])-27],'-','LineWidth',1.3,'color',[0 0 0]);
end
text(88.65, min([min(i1psi) min(i1phi) min(i1l) min(i20)])-40,'30', 'FontSize',12);
text(98.65, min([min(i1psi) min(i1phi) min(i1l) min(i20)])-40,'20', 'FontSize',12);
text(108.65, min([min(i1psi) min(i1phi) min(i1l) min(i20)])-40,'10', 'FontSize',12);

for i=13:15
    plot([10*i 10*i], [min([min(i1psi) min(i1phi) min(i1l) min(i20)])-43 min([min(i1psi) min(i1phi) min(i1l) min(i20)])-37],'-','LineWidth',1.3,'color',[0 0 0]);
end
text(128.65, min([min(i1psi) min(i1phi) min(i1l) min(i20)])-50,'10', 'FontSize',12);
text(138.65, min([min(i1psi) min(i1phi) min(i1l) min(i20)])-50,'20', 'FontSize',12);
text(148.65, min([min(i1psi) min(i1phi) min(i1l) min(i20)])-50,'30', 'FontSize',12);

%Стрелки возле графиков и подписи к ним 

%i1psi
ii1psi = annotation('arrow');  
ii1psi.Parent = gca;           
ii1psi.X = [Xpsi(2000) Xpsi(2000)-15];          
ii1psi.Y = [i1psi(2000) i1psi(2000)];   %меняйте координаты высоты стрелки, если надо   
ii1psi.LineWidth  = 0.5;         
ii1psi.HeadWidth  = 4;
ii1psi.HeadLength = 4;
text(Xpsi(2000)-15, i1psi(2000)+12,'i_{\psi1}(x_{1}'')', 'FontSize',12);

%i1l
ii1l = annotation('arrow');  
ii1l.Parent = gca;           
ii1l.X = [X1l(2000) X1l(2000)-15];          
ii1l.Y = [i1l(2000) i1l(2000)];   
ii1l.LineWidth  = 0.5;         
ii1l.HeadWidth  = 4;
ii1l.HeadLength = 4;
text(X1l(2000)-15, i1l(2000)+11,'i_{1l}(x_{1}'')', 'FontSize',12);

%i20
ii20 = annotation('arrow');  
ii20.Parent = gca;           
ii20.X = [X20(1500) X20(1500)+15];          
ii20.Y = [i20(1500) i20(1500)];   
ii20.LineWidth  = 0.5;         
ii20.HeadWidth  = 4;
ii20.HeadLength = 4;
text(X20(1500)+15-2, i20(1500)+30,'i_{20}(x_{2})', 'FontSize',12);

%изменяйте координаты расположения текста как вам надо, чтоб было красиво
text(12, i1phi(1)+15,'i_{\phi1}', 'FontSize',12);
% ax=gca
% ax.XColor=[240 240 240]/255;
% ax.YColor=[240 240 240]/255;