ИЭ / 6 семестр (англ) / Курсач / Графики в матлабе+эксель / Hx3_Hx1_cylinder

I=100;
m=600;
Ib=I*(m-1)/(m+1);
I0=2*I/(m+1);
R=0.5;
a=1.25;
b=0.35;
D=0.016;
XA1=a+b;
XA=a;
XB=(R^2)/a;
XB1=(R^2)/(a+b);
x=-1:0.001:3;
h1=x-x;
h2=x-x;
h3=x-x;
h4=x-x;

for i=1:length(x)
if ((x(i)>=(XA1-0.5*D)) && (x(i)<=(XA1+0.5*D)))
h1(i)=I*(XA1-x(i))/(2*pi*(0.25*D^2));
else
if ((x(i)>=-R) && (x(i)<=R))
h1(i)=I0/(2*pi*(XA1-x(i)));
else
h1(i)=I/(2*pi*(XA1-x(i)));
end
end
end

for i=1:length(x)
if ((x(i)>=(XA-0.5*D)) && (x(i)<=(XA+0.5*D)))
h2(i)=-I*(XA-x(i))/(2*pi*(0.25*D^2));
else
if ((x(i)>=-R) && (x(i)<=R))
h2(i)=-I0/(2*pi*(XA-x(i)));
else
h2(i)=-I/(2*pi*(XA-x(i)));
end
end
end

for i=1:length(x)
if ((x(i)>=-R) && (x(i)<=R))
h3(i)=0;
else
h3(i)=-Ib/(2*pi*(XB-x(i)));
end
end

for i=1:length(x)
if ((x(i)>=-R) && (x(i)<=R))
h4(i)=0;
else
h4(i)=Ib/(2*pi*(XB1-x(i)));
end
end

H=h1+h2+h3+h4;

figure
hold on
grid on
set(gca, 'box', 'on');
set(gca,'FontSize',14);

plot(x,log10(abs(H)),'-k','LineWidth',1.3);
xlabel('x, m', 'FontSize', 22);
ylabel('lg(|H_{y}(x)|)', 'Rotation', 0, 'FontSize', 22, 'Position', [-1.12 0 0]);
set(gca,'XTick',(-1:0.1:1));
clear xlim;
xlim([-1 1]);


p=x3(1)/100:0.01:(x3(end)+0.1)/100;
f1=interp1(x3/100, lg3, p, 'pchip');
plot(p,f1,'-k','LineWidth',1.3, 'color',[0.75 0.75 0.75]);
xlabel('x, m', 'FontSize', 22);
ylabel('lg(|H_{y}(x)|)', 'Rotation', 0, 'FontSize', 22, 'Position', [-1.12 0 0]);
set(gca,'XTick',(-1:0.1:1));

LH(1) = plot(nan, nan, '-k','LineWidth',1.3,'color',[0 0 0]);
LH(2) = plot(nan, nan, '-','LineWidth',1.3,'color',[0.75 0.75 0.75]);
L{1} = 'method of images / method of the variable separation';
L{2} = 'Q-FIELD';
legend(LH,L);