Cuestion 4 (Sistema Discreto, regulador, polo dominante)

z=%z;
gz=2*(z-0.3)/((z+0.3)*(z-0.5));
gzs=syslin('d',gz);
z1=-0.6+%i*0.5
//Angulo en z1 de G
aux=horner(gzs,z1);
aux1=atan(imag(aux),real(aux));
aux2=360*aux1/(2*%pi)
//El angulo del regulador de atras
aux3=-180-aux2
aux4=-2*%pi*aux3 /360;
//Vamos a situar el cero del regulador en -0.6
b=0.6;
aux5=imag(z1)*tan(aux4)
a=-(real(z1)+aux5)
//Parte del regulador nos quedara
gc=(z+b)/(z+a)
gcs=syslin('d',gc);
aux6=horner(gcs,z1);
aux7=atan(imag(aux6),real(aux6));
aux8=360*aux7/(2*%pi)
//Vamos a calcular la kc del regulador
aux9=horner(gz*gc,z1);
aux10=abs(aux9)
kc=1/aux10
//Vamos a comprobar los resultados
gt=kc*gc*gz;
gts=syslin('d',gt);
aux11=horner(gts,z1);
aux12=atan(imag(aux11),real(aux11));
aux13=360*aux12/(2*%pi)
aux14=abs(aux11)
//Representacion del lugar de las raices
clf;
subplot(2,1,1)
evans(gzs)
plot(real(z1),imag(z1),'*');
plot(real(z1),-imag(z1),'*');
mtlb_axis([-3 1.5 -2 2])
subplot(2,1,2)
evans(gts)
plot(real(z1),imag(z1),'*');
plot(real(z1),-imag(z1),'*');
mtlb_axis([-3 1.5 -2 2])