Programa en Scilab
s=%s;
num=s+3;
den=s*(s+1)*(s^2+4*s+16);
g=syslin('c',num/den);
//Lugar de los raices
evans(g);
//punto de corte en el eje imaginario y ganancia
[Ki,s]=kpure(g)
//Los puntos de corte con el eje imaginario dibujados con un asterisco
plot([real(s) real(s)],[imag(s) -imag(s)],'*r')
v=[-6 6 -6 6]; mtlb_axis(v);
Resultados:
s =
3.140933i
Ki =
33.3273
Añadir nuevo comentario