ИЭ / 4 семестр / Курсач от Кочетковой / Графики цилиндр в матлабе / UE24

clc, close all
R1=0.01;
R12=0.025;
R2=0.05;
r1=R1:0.0001:R12;
r2=R12:0.0001:R2;

%е1=2е0
U11=35310734.46*r1.^2-30282.04*log(r1)-142985.02; %введи своё U1
U21=27712.75*log(r2)+93019.99; %введи своё U2
E11=-70621468.93*r1 +30282.04./r1; %введи своё E1
E21=-27712.75./r2; %введи своё E2

%е1=е0
U12=70621468.93*r1.^2-54842.2*log(r1)-259619.81; %введи своё U1
U22=33434.64*log(r2)+110161.22; %введи своё U2
E12=-141242937.9*r1 +54842.2./r1; %введи своё E1
E22=-33434.64./r2; %введи своё E2

%е1=6е0
U13=11770244.82*r1.^2-11303.7*log(r1)-53232.49; %введи своё U1
U23=20454.64*log(r2)+71276.61; %введи своё U2
E13=-23540489.64*r1+11303.7./r1; %введи своё E1
E23=-20454.64./r2; %введи своё E2


figure;
subplot(2,1,1);
hold on;
set(gca,'box','on');
grid on;
plot(r1,U11,'-','color',[0 0 0]/255,'LineWidth',1.3);
plot(r2,U21,'-','color',[0 0 0]/255,'LineWidth',1.3);
plot([R12, R12], [-30000, U11(end)],'--','color',[0.5 0.5 0.5],'LineWidth',1.3);
text(R12+0.001, -15000, 'R_{12}=0.025 м', 'FontSize',14);
text(0.015, -6000, ' $\varepsilon_{1}=2\varepsilon_{0}$ ','interpreter','latex', 'FontSize',14);

plot(r1,U12,'-','color',[100 100 100]/255,'LineWidth',1.3);
plot(r2,U22,'-','color',[100 100 100]/255,'LineWidth',1.3);
plot([R12, R12], [-30000, U12(end)],'--','color',[0.5 0.5 0.5],'LineWidth',1.3);
text(0.015, -12000, ' $\varepsilon_{1}=\varepsilon_{0}$ ','interpreter','latex', 'FontSize',14);

plot(r1,U13,'-','color',[185 185 185]/255,'LineWidth',1.3);
plot(r2,U23,'-','color',[185 185 185]/255,'LineWidth',1.3);
plot([R12, R12], [-30000, U13(end)],'--','color',[0.5 0.5 0.5],'LineWidth',1.3);
text(0.015, 0, ' $\varepsilon_{1}=6\varepsilon_{0}$ ','interpreter','latex', 'FontSize',14);

xlim([0.01 0.05])
ylim([-20000 10000]);
xlabel('r, м','fontsize',14);
ylabel('U, В','fontsize',14,'Rotation',0,'Position', [0.006 -7500 0]);

subplot(2,1,2);
hold on;
set(gca,'box','on');
grid on;
plot(r1,E11,'-','color',[0 0 0]/255,'LineWidth',1.3);
plot(r2,E21,'-','color',[0 0 0]/255,'LineWidth',1.3);
plot([R12, R12], [-4000000, E11(end)],'--','color',[0.5 0.5 0.5],'LineWidth',1.3);
text(R12-0.01, -10^6, 'R_{12}=0.025 м', 'FontSize',14);
text(0.0115, 2*10^6, ' $\varepsilon_{1}=2\varepsilon_{0}$ ','interpreter','latex', 'FontSize',14);

plot(r1,E12,'-','color',[100 100 100]/255,'LineWidth',1.3);
plot(r2,E22,'-','color',[100 100 100]/255,'LineWidth',1.3);
plot([R12, R12], [-4000000, E12(end)],'--','color',[0.5 0.5 0.5],'LineWidth',1.3);
text(0.0115, 3.8*10^6, ' $\varepsilon_{1}=\varepsilon_{0}$ ','interpreter','latex', 'FontSize',14);

plot(r1,E13,'-','color',[185 185 185]/255,'LineWidth',1.3);
plot(r2,E23,'-','color',[185 185 185]/255,'LineWidth',1.3);
plot([R12, R12], [-4000000, E13(end)],'--','color',[0.5 0.5 0.5],'LineWidth',1.3);
text(0.0115, 0*10^6, ' $\varepsilon_{1}=6\varepsilon_{0}$ ','interpreter','latex', 'FontSize',14);

xlim([0.01 0.05])
ylim([-1.5*10^6 4.5*10^6])
xlabel('r, м','fontsize',14);
ylabel('E, В/м','fontsize',14,'Rotation',0,'Position', [0.006 0 0]);