Problema A6.15 OGATA 4edicion pag407

Solapas principales

Vamos a dibujar el lugar de las raices mediante Scilab de un sistema en lazo cerrado con realimentación unitaria cuyo funcion de transferencia en lazo abierto es:

 

Funcion de transferencia, Transformada de Laplace


 

Programa en Scilab

s=%s;

num=(s^2+25)*s;

den=(s^4+404*s^2+1600);

g=syslin('c',num/den);

evans(g)

v=[-20 20 -21 21];

mtlb_axis(v)
Lugar de las raices del sistema con Scilab

El resultado que obtendremos no nos valdra la grafica esta deformada tendremos que programarlo de la siguiente manera

 
Programa en Scilab

s=%s;

num=(s^2+25)*s;

den=(s^4+404*s^2+1600);

g=syslin('c',num/den);

for k1=0.2:0.2:20,

  gl=1+k1*g;

  numl=numer(gl);

  r=roots(numl);

  x=real(r);

  y=imag(r);

  plot(x,y,'o');

end;

for k2=20.02:0.2:30,

  g2=1+k2*g;

  num2=numer(g2);

  r2=roots(num2);

  x2=real(r2);

  y2=imag(r2);

  plot(x2,y2,'o');

end;

for k3=35:5:1000,

  g3=1+k3*g;

  num3=numer(g3);

  r3=roots(num3);

  x3=real(r3);

  y3=imag(r3);

  plot(x3,y3,'o');

end;

v=[-20 20 -21 21];

mtlb_axis(v)

xtitle('lugar de las raices','Eje Real','Eje imaginario')

xgrid;

 

Lugar de las raices con puntos del sistema con Scilab

 

Español

Añadir nuevo comentario

Plain text

  • No se permiten etiquetas HTML.
  • Las direcciones de las páginas web y las de correo se convierten en enlaces automáticamente.
  • Saltos automáticos de líneas y de párrafos.
Pin It