Mensaje de error

  • Deprecated function: The each() function is deprecated. This message will be suppressed on further calls en book_prev() (línea 775 de /home1/montes/public_html/drupal/modules/book/book.module).
  • Notice: Trying to access array offset on value of type int en element_children() (línea 6422 de /home1/montes/public_html/drupal/includes/common.inc).
  • Notice: Trying to access array offset on value of type int en element_children() (línea 6422 de /home1/montes/public_html/drupal/includes/common.inc).
  • Notice: Trying to access array offset on value of type int en element_children() (línea 6422 de /home1/montes/public_html/drupal/includes/common.inc).
  • 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).

Cuestion 3 (Sistema Discreto, estabilidad, criterio de Jury)

Solapas principales

Enunciado de la cuestion 3 del examen de la 2 semana de Septiembre de Regulacion Automatica II

Solucion:

Dado el siguiente sistema discreto que se realimenta negativa y unitariamente. Vamos a comprobar su estabilidad.

 

Funcion de transferencia del sistema discreto, Transformada Z


 

Vamos a utilizar el criterio de estabilidad de Jury para verificar la estabilidad. La Ecuacion caracteristica sera:

 

Ecuacion del criterio de estabilidad de Jury


 

 

Ecuacion del criterio de estabilidad de Jury con los coeficientes


 

El criterio tiene que verificar:
  • Condicion del criterio de estabilidad de Jury

  • Condicion del criterio de estabilidad de Jury parte2
  • Condicion del criterio de estabilidad de Jury parte3 si n par y

    Condicion del criterio de estabilidad de Jury parte4 si n impar

  • Condicion del criterio de estabilidad de Jury parte5

  • ....

 

La primera condicion no se verifica $3.5<1$ con lo que no es estable

 

Calculos y comprobacion con scilab:

 

z=%z;
g=2*(z-2)/(z^4+0.5*z^2+0.5)
gz=syslin('d',g)

//Ecuacion caracteristica
ec=numer(gz)+denom(gz)


//Comprobacion mediante transformacion bilineal y Routh
s=%s
z1=(s+1)/(s-1)
g1=2*(z1-2)/(z1^4+0.5*z1^2+0.5)
gs1=syslin('c',g1)
axec2=numer(gs1)+denom(gs1);
ec2=numer(axec2)
routh_t(ec2)

 

 

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