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

clc, close all
R1=0.01;
R121=0.025;
R2=0.05;
r11=R1:0.0001:R121;
r21=R121:0.0001:R2;
U11=35310734.46*r11.^2-30282.04*log(r11)-142985.02; %введи своё U1
U21=27712.75*log(r21)+93019.99; %введи своё U2
E11=-70621468.93*r11 +30282.04./r11; %введи своё E1
E21=-27712.75./r21; %введи своё E2

R122=0.012;
r12=R1:0.0001:R122;
r22=R122:0.0001:R2;
U12=35310734.46*r12.^2-6777.34*log(r12)-34741.88; %введи своё U1
U22=6784.3*log(r22)+30323.96; %введи своё U2
E12=-70621468.93*r12 +6777.34./r12; %введи своё E1
E22=-6784.3./r22; %введи своё E2

R123=0.04;
r13=R1:0.0001:R123;
r23=R123:0.0001:R2;
U13=35310734.46*r13.^2-50963.09*log(r13)-238224.78; %введи своё U1
U23=124062.52*log(r23)+381658.1; %введи своё U2
E13=-70621468.93*r13+50963.09./r13; %введи своё E1
E23=-124062.52./r23; %введи своё E2


figure;
subplot(2,1,1);
hold on;
set(gca,'box','on');
grid on;
plot(r11,U11,'-','color',[0 0 0]/255,'LineWidth',1.3);
plot(r21,U21,'-','color',[0 0 0]/255,'LineWidth',1.3);
plot([R121, R121], [-30000, U11(end)],'--','color',[0.5 0.5 0.5],'LineWidth',1.3);
text(R121+0.001, -11000, 'R_{12}=0.025 м', 'FontSize',14);

plot(r12,U12,'-','color',[100 100 100]/255,'LineWidth',1.3);
plot(r22,U22,'-','color',[100 100 100]/255,'LineWidth',1.3);
plot([R122, R122], [-30000, U12(end)],'--','color',[0.5 0.5 0.5],'LineWidth',1.3);
text(R122+0.0005, 5000, 'R_{12}=0.012 м', 'FontSize',14);

plot(r13,U13,'-','color',[185 185 185]/255,'LineWidth',1.3);
plot(r23,U23,'-','color',[185 185 185]/255,'LineWidth',1.3);
plot([R123, R123], [-30000, U13(end)],'--','color',[0.5 0.5 0.5],'LineWidth',1.3);
text(R123+0.0005, -18500, 'R_{12}=0.04 м', 'FontSize',14);

xlim([0.01 0.05])
ylim([-30000 10000]);
xlabel('r, м','fontsize',14);
ylabel('U, В','fontsize',14,'Rotation',0,'Position', [0.007 -7500 0]);
subplot(2,1,2);
hold on;
set(gca,'box','on');
grid on;
plot(r11,E11,'-','color',[0 0 0]/255,'LineWidth',1.3);
plot(r21,E21,'-','color',[0 0 0]/255,'LineWidth',1.3);
plot([R121, R121], [-4000000, E11(end)],'--','color',[0.5 0.5 0.5],'LineWidth',1.3);
text(R121+0.001, -2000000, 'R_{12}=0.025 м', 'FontSize',14);

plot(r12,E12,'-','color',[100 100 100]/255,'LineWidth',1.3);
plot(r22,E22,'-','color',[100 100 100]/255,'LineWidth',1.3);
plot([R122, R122], [-4000000, E12(end)],'--','color',[0.5 0.5 0.5],'LineWidth',1.3);
text(R122+0.0005, -3000000, 'R_{12}=0.012 м', 'FontSize',14);

plot(r13,E13,'-','color',[185 185 185]/255,'LineWidth',1.3);
plot(r23,E23,'-','color',[185 185 185]/255,'LineWidth',1.3);
plot([R123, R123], [-4000000, E13(end)],'--','color',[0.5 0.5 0.5],'LineWidth',1.3);
text(R123+0.0005, -2000000, 'R_{12}=0.04 м', 'FontSize',14);

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