function [eje_x,eje_menos_y,eje_mas_y]=plot19d(L,zl1u,pzona)
% CIRCULO MHO (datos: L,zl1u,pzona)
modzl1u=abs(zl1u);
Lmodzl1u=L*modzl1u;
%kx=input('Cte.para mas x y menos x(positiva): kx= ');
%kmenosy=input('Cte.para menos y(negativa): menos ky= ');
%kmasy=input('Cte.para mas y: ky= ');
%**************************************
% VALORES TIPICOS DE kx, kmenosy, kmasy
%**************************************
%kx=1;
%kmenosy=-0.5;
%kmasy=1.5;
%**************************************
%**************************************
% OTROS VALORES TIPICOS DE kx, kmenosy, kmasy
%**************************************
%kx=3;
%kmenosy=-3;
%kmasy=3;
%**************************************
kx=1.5;
kmenosy=-1;
kmasy=1.5;
eje_x=kx*Lmodzl1u;
eje_menos_y=kmenosy*Lmodzl1u;
eje_mas_y=kmasy*Lmodzl1u;
zl1=zl1u*L;
zl1pzona=zl1*pzona;
centro=zl1pzona/2;
recentro=real(centro);
imcentro=imag(centro);
radio=sqrt(recentro^2+imcentro^2);
extizq=recentro-radio;
extder=recentro+radio;
intervaloX=extder-extizq;
incrementoX=intervaloX/100;
X=zeros(1,101);
Y1=zeros(1,101);
Y2=zeros(1,101);
X(1)=extizq;
Y1(1)=imcentro;
Y2(1)=imcentro;
for n=2:100
X(n)=X(n-1)+incrementoX;
Y1(n)=imcentro+sqrt(radio^2-(X(n)-recentro)^2);
Y2(n)=imcentro-sqrt(radio^2-(X(n)-recentro)^2);
end
X(101)=extder;
Y1(101)=imcentro;
Y2(101)=imcentro;
Xreves=X(:,101:-1:1);
XX=zeros(1,202);
XX=[X Xreves];
Y2reves=Y2(:,101:-1:1);
Y1Y2=zeros(1,202);
Y1Y2=[Y1 Y2reves];
plot(XX,Y1Y2,'g')
axis('equal')
axis([-kx*Lmodzl1u kx*Lmodzl1u kmenosy*Lmodzl1u kmasy*Lmodzl1u]);
grid
xlabel('R(Ohm)')
ylabel('X(Ohm)')
% LINEA (datos: L,zl1u)
rezl1=real(zl1);
imzl1=imag(zl1);
factor=imzl1/rezl1;
hold on
X3=zeros(1,2);
Y3=zeros(1,2);
X3(1)=0;
Y3(1)=0;
X3(2)=rezl1;
Y3(2)=imzl1;
plot(X3,Y3,'b')
hold off
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
