Mensaje de error

Deprecated function: implode(): Passing glue string after array is deprecated. Swap the parameters en drupal_get_feeds() (línea 394 de /home1/montes/public_html/drupal/includes/common.inc).

numer

Español

Ejercicio B.6.15

[adsense:336x280:9156825571]
Vamos a dibujar el lugar de las raíces del sistema en lazo abierto y calcular las raíces en lazo cerrado cuando $K=2$

 

Funcion de transferencia en lazo abierto, Tranformada de Laplace

 

Programa en Scilab:
//Borramos la ventana grafica

clf;
//Declaramos la variable s y la funciones 'g','h' y 'g*h'

s=%s;

g=(s+1)/(s*(s^2+2*s+6));

h=1/(s+1);

gh=g*h;

//Declaramos la funcion g*h lineal y continua

ghs=syslin('c',gh)

//Declaramos ghs2 como la funcion anterior
 pero con K=2

ghs2=2*ghs;
[adsense:336x280:9156825571]
//Calculamos las raices de (1+ghs2)*(s+1), este ultimo es el polo que se anulo
 con el cero

glc=(1+ghs2)*(s+1);

r=roots(numer(glc))

//Dibujamos el lugar de las raices de la funcion ghs

evans(ghs);
mtlb_axis([-5 5 -5  5])

//Dibujamos el polo y el cero que se anularon:

plot(-1,0,'o');

plot(-1,0,'x');

//Dibujamos las raices para K=2;

plot(real(r),imag(r),'+');

xgrid;

xtitle('Lugar de las raices, los polos en lazo cerrado -> + para K=2','Eje real'
,'Eje imaginario');


Resultados
-->r=roots(numer(glc))
 r  =
 
  - 0.3706383              
  - 1.                     
  - 0.8146808 + 2.175406i  
  - 0.8146808 - 2.175406i



Lugar de las raices en lazo cerrado con Scilab
Español

Problema A6.15 OGATA 4edicion pag407

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

Programa 6.8 OGATA 4edicion pag370

Vamos a dibujar el lugar de las raices con indicacines de factor de amortiguamiento $\zeta$=0.5 y se pedira marcar un punto del lugar de las raices y se obtendra la ganancia K y los polos en lazo cerrado mediante Scilab de un sistema en lazo cerrado con realimentación unitaria cuya funcion de transferencia en lazo abierto es:

 

Funcion de transferencia, Transformada de Laplace


Programa en Scilab

num=poly([1 0 0 0],'s','coeff');

den=poly([0 5 4 1],'s','coeff');

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

evans(g);

v=[-3 1 -2 2];

mtlb_axis(v)

sgrid([0.5],[0],32);

//marcar un punto con el raton en el lugar de las raices de la grafica

p=locate(1)

 k=-1/real(horner(g,[1,%i]*p));

 gl=1+k*g;

numgl=numer(gl);

roots(numgl)
k

Resultados:

-->p=locate(1)
 p  =
 
  - 0.6502732  
    1.0550725 
    
-->roots(numgl)
 ans  =
 
  - 0.6412384 + 1.0505487i  
  - 0.6412384 - 1.0505487i  
  - 2.7175231               
 
-->k
 k  =
 
    4.1166111
Español

Programa 6.3 OGATA 4edicion pag362 (Lugar de las Raices)

Vamos a calcular 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 tranferencia del sistema. Laplace


Programa en Scilab

num=poly([1 0 0 0 0],'s','coeff');

den=poly([0 5 10.3 1.1 1],'s','coeff');

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 -25 25];

mtlb_axis(v)

xgrid;

Lugar de las raices por puntos de la trasformada de Laplace con Scilab

Español

Páginas

Pin It
Subscribe to RSS - numer