Cuestion 4 (Sistemas Discretos, regulador)

Comprobaciones y calculos con el Scilab

z=%z;
gz=(z-0.1)/((z+0.1)*(z-0.5));
gzs=syslin('d',gz);
z1=-0.3+%i*0.6
//Angulo en z1 de G
aux=horner(gzs,z1);
aux1=atan(imag(aux),real(aux));
aux2=360*aux1/(2*%pi)
//El angulo del regulador
aux3=-180-aux2
aux4=-2*%pi*aux3 /360
//Vamos a situar el polo del regulador en -0.3 b=0.3;
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)
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])