Teoria Control
Teoria Control
Introducción al
Control de Procesos
para Ingenieros
Jean-François DULHOSTE
Escuela de Ingeniería Mecánica - ULA
2 Introducción al Control de Procesos para Ingenieros
Jean-François DULHOSTE
Tabla de Contenido
Tala de Contenido..................................................................................................................3
Introducción...........................................................................................................9
Tema 1. Introducción a los Sistemas de Control..........................................................................11
Sistema de Control..........................................................................................................11
Definición...................................................................................................................... 11
Ejemplo de un Sistema de Control.....................................................................................11
Ventajas de un control automático......................................................................................11
Términos Básicos............................................................................................................11
Elementos de un Sistema de Control....................................................................................11
Clases de sistemas de control...................................................................................12
Sistema de control Pasivo, sistema de control de Lazo Abierto y sistema de control Retroalimentado.
12 Sistema de Control Continuo vs. Sistema de Control Discontinuo...........................................14
Sistemas de Control Analógicos vs. Sistemas de Control Digitales.............................................17
Sistemas de Control Reguladores vs. Sistemas de Control Seguidores........................................17
Sistemas de control invariantes en el tiempo vs. Sistemas de control variable en el tiempo.............18
Sistemas de control con parámetros concentrados vs. Sistemas de control con parámetros distribuidos. .19
Sistemas de control lineales vs. Sistemas de control no lineales..............................................19
Sistemas de control una entrada y una salida (SISO) vs. Sistemas de control de múltiples entradas
y salidas (MIMO).................................................................................................19
Sistemas de control determinísticos vs. Sistemas de control estocásticos............................20
Características de los sistemas de control....................................................................20
Estabilidad...................................................................................................20
Exactitud......................................................................................................21
Velocidad de Respuesta.................................................................................................21
Lazos de control comúnmente utilizados en procesos industriales................................................21
Control en cascada.....................................................................................................22
Control de relación.........................................................................................................22
Control de rango partido...............................................................................................23
Ejercicios..........................................................................................................23
Tema 2. Modelado Matemático.............................................................................................25
Introducción.......................................................................................................25
Sistemas Mecánicos..............................................................................................25
Sistemas Mecánicos Rotativos..................................................................................26
Sistemas Eléctricos...............................................................................................27
Analogía Electromecánica................................................................................................29
Diagrama de impedancias.............................................................................................29
Sistemas Térmicos................................................................................................30
Sistemas Hidráulicos.............................................................................................32
Sistemas Neumáticos.............................................................................................34
Ejercicios..........................................................................................................36
Tema 3. Formas de representación del modelo matemático.........................................................38
Introducción.......................................................................................................38
Representación de un modelo matemático con la Función de Transferencia...................................39
La Matriz de Transferencia............................................................................................40
Polos y ceros de un sistema LTI, SISO................................................................................41
Representación de un modelo matemático en Espacio de Estado.................................................43
Definiciones..............................................................................................................44
Representación de sistemas dinámicos en el espacio de estado...............................................45
Relación entre función de transferencia y espacio de estado...................................................50
No unicidad del conjunto de variables de estado.................................................................52
Tema 4. Diagramas de Bloque..............................................................................................59
Introducción.......................................................................................................59
Elementos de un diagrama de bloques.................................................................................59
Ejemplos de diagrama de bloques...................................................................................60
Forma canónica de un sistema de control retroalimentado.....................................................60
Álgebra de bloques......................................................................................................61
Simplificación de diagramas de bloque.................................................................................63
Diagramas de bloque de varias entradas...........................................................................64
Entradas y salidas múltiples..................................................................................65
Construcción de un diagrama de bloques a partir de las funciones de transferencia de los elementos del
sistema.......................................................................................................68
Tema 5. Respuesta de Sistemas...................................................................................74
Introducción.......................................................................................................74
Tipos de excitación............................................................................................................74
Entrada en escalón......................................................................................................74
Entrada en rampa.......................................................................................................75
Entrada sinusoidal..........................................................................................75
Entrada en pulso............................................................................................75
Entrada en impulso.........................................................................................76
Entrada en onda cuadrada................................................................................76
Respuestas de sistemas de primer y segundo orden.........................................................77
Respuesta de sistemas de Primer Orden................................................................77
Respuesta de sistemas de Segundo Orden..............................................................83
Solución de las ecuaciones en representación de estado.................................................91
Tema 6. Estabilidad de Sistemas..............................................................................98
Introducción............................................................................................98
Definiciones de estabilidad BIBO.........................................................................98
Estabilidad limitada........................................................................................99
Conclusión...................................................................................................99
Criterio de Estabilidad Routh..............................................................................100
Observaciones al criterio de Routh.....................................................................101
Criterio de estabilidad de Hurwitz........................................................................103
Análisis de estabilidad de sistemas en Espacio de Estado............................................105
Determinación de los valores propios de una matriz....................................................105
Ejercicios...................................................................................................106
Tema 7. Análisis de la respuesta en frecuencia...........................................................109
Introducción..................................................................................................109
Calculo de la respuesta en frecuencia...................................................................109
Procedimiento para la obtención de la respuesta en frecuencia..............................................110
Diagramas de BODE..........................................................................................110
Diagramas de Bode de Funciones comunes................................................................111
Procedimiento general para trazar diagramas de Bode....................................................114
Análisis de estabilidad utilizando los diagramas de Bode...............................................114
Diagrama de Nyquist........................................................................................122
Diagramas de Nyquist de Funciones comunes........................................................123
Formas generales de los diagramas polares...........................................................125
Ejemplos....................................................................................................126
Análisis de estabilidad utilizando diagramas de Nyquist................................................126
Ejemplos....................................................................................................128
Tema 8. Controlabilidad y Observabilidad.................................................................133
Controlabilidad...............................................................................................133
Observabilidad...............................................................................................136
Dualidad.......................................................................................................137
Teoría de la realización........................................................................................137
Formas canónicas de representación de estado........................................................138
Ejercicios.......................................................................................................141
Tema 9. Acciones de Control.................................................................................143
Introducción..................................................................................................143
Tipos de acción de control.................................................................................143
1............................................Acción de control discontinua o de dos posiciones (ON-OFF) 143
2. Acción de control proporcional (P)..................................................................144
3. Acción de control derivativa (D).....................................................................145
4. Acción de control integral (I).........................................................................146
5. Control proporcional más derivativa (PD)..........................................................147
6. Control proporcional más integral (PI)..............................................................147
7.Control proporcional más integral más derivativo (PID)..............................................148
Resumen....................................................................................................148
Ejemplo de estudio del efecto de una acción de control en un sistema........................................149
Si el controlador es proporcional.......................................................................149
Si el controlador es proporcional mas derivativo.........................................................151
Si el controlador es proporcional más integral.......................................................153
Si el controlador es proporcional más integral más derivativo.........................................155
Tema 10. Ajuste de Controladores PID.....................................................................159
Introducción..................................................................................................159
Métodos de Ziegler - Nichols...............................................................................159
Método de Ziegler - Nichols basado en la respuesta al escalón........................................159
Método de Ziegler - Nichols basado en la respuesta frecuencial.......................................160
Métodos por asignación de Polos Dominantes...........................................................161
Polos Dominantes..........................................................................................161
Método aproximado para designar polos dominantes....................................................164
Ejercicio....................................................................................................166
Tema 11. Control por Retorno de Estado...................................................................169
Retorno de estado lineal...................................................................................169
El control modal...........................................................................................170
Introducción
El control de procesos, que en sus inicios estaba restringido a máquinas sofisticadas y procesos muy
complejos y costosos, está hoy en día en prácticamente todas las actividades humanas. Entre estas
actividades las actividades de la Ingeniería Mecánica tienen una importancia primordial, pues estas van
desde las maquinarias industriales, pasando por los automóviles hasta llegar a los equipos simples de
los hogares como hornos, neveras, calentadores, etc.
A la vez de si difusión en las actividades humanas las técnicas del control de procesos han evolucionado, y
se ha vuelto una ciencia que para ser manejado en su globalidad requiere de estudios especiales de
este dominio, sin embargo casi todos los ingenieros mecánicos se topan en su carrera con
sistemas automatizados a los cuales deberán operar, mantener o incluso modificar. En vista de esto la
carrera de Ingeniero mecánico incluye en su pensum dos asignaturas referidas a esta ciencia, la
primera Instrumentación cuyo objetivo es el aprendizaje de la ciencia de las mediciones de variables y la
segunda Teoría de Control cuyo objetivo es el aprendizaje de los fundamentos primordiales del control de
procesos.
Este documento se presenta como una ayuda a la enseñanza de la asignatura Teoría de Control para la
carrera de Ingeniería Mecánica, para lo cual se ha hecho un intento de presentar de una forma simple
los conocimientos primordiales que debería poseer un Ingeniero Mecánico para trabajar con sistemas de
control de procesos en todos sus ámbitos de trabajo, y que deben ser adquiridos en el transcurso de un
semestre académico. No pretende por lo tanto ser un documento extenso sobre el tema, como lo son
libros fundamentales sobre el tema.
Con este objetivo el documento se desarrolla en tres unidades que engloban estos conocimientos
primordiales del control de procesos.
Una primera unidad referida a la Representación Matemática de Sistemas Lineales, la cual presenta las
diversas formas de representación de sistemas lineales desde el punto de vista matemático y gráfico, útiles
para las diversas técnicas existentes en el control de procesos, y además comienza con una somera
introducción al control de procesos, esta se divide en cuatro temas: Tema 1 Introducción a los sistemas de
control, Tema 2 Modelado Matemático, Tema 3 Formas de representación de modelos matemáticos y Tema
4 Diagramas de Bloque.
Una segunda unidad referida al Análisis de la Dinámica de los Sistemas, en la cual se presentan los
principales criterios de estudio en referencia a los sistemas desde el punto de vista de sus diversas
representaciones. Esta incluye los siguientes temas: Tema 5 Respuesta de sistemas, Tema 6 Estabilidad de
sistemas, Tema 7 Análisis de la respuesta en frecuencia y Tema 8 Observabilidad y Controlabilidad .
Una tercera unidad referida a las Técnicas de Control, en la cual se presentan las principales técnicas de
control desde el punto de vista del control clásico o del control moderno, y culmina con una introducción a
los sistemas digitales. Esta unidad está organizada en los siguientes temas: Tema 9 Acciones de
control, Tema 10 Ajuste de Controladores PID, Tema 11 Control por retorno de estado y Tema 12
Introducción a los Sistemas Digitales.
Definición
Una definición de un sistema de control puede ser:
“Es un arreglo de componentes físicos conectados de tal manera, que el arreglo pueda comandar, dirigir o
regular a sí mismo o a otro sistema”.
Pivote
Flotador
Válvula
Entra agua
Agua
Sale agua
Términos Básicos
Planta. Es un equipo o conjunto de equipos que permiten realizar una operación determinada. Cuando se
tiene un conjunto de equipos interactuando para generar un producto se tiene una planta industrial.
Proceso. Está constituido por una serie de operaciones coordinadas sistemáticamente para producir
un resultado final que puede ser un producto.
Sistema. Es una combinación de componentes físicos que actúan conjuntamente para cumplir un
determinado objetivo.
Esto se puede representar en forma de un diagrama de bloques que muestra la interacción lógica de
los elementos que conforman un sistema de control retroalimentado.
Perturbación
Variable manipulada Variable
Señal de Elemento final de Proceso a
referencia Error Controlador Controlada
Elemento de medición
Señal de retroalimentación o
medición
Jean-François DULHOSTE
Sistema de control continuo
Es aquel en que la modificación de la
variable manipulada se efectúa Transmisor
Controlador
continuamente. de
Por ejemplo un controlador de temperatura
con vapor.
Este instrumento calienta el agua con la
ayuda de un serpentín por el cual pasa Sale agua
un flujo de vapor. Este flujo de vapor caliente
es aumentado o disminuido por una válvula
de control, la cual se abre o cierra
suavemente en una cierta proporción
Vapor
según el mensaje enviado por el
controlador, el cual trabaja con la Intercambiad
diferencia entre una señal de Entra or de calor
referencia y el valor de la temperatura del agua fría
agua de salida. La temperatura del agua es Drenaje
medida por un transmisor de temperatura.
Como la apertura de la válvula puede
tomar
una infinidad de valores entre totalmente Esquema de un calentador de agua a vapor
cerrada y totalmente abierta, entonces el
control es continuo.
Máximo
Mínimo
t t
Sistema continuo Sistema discontinuo
Sistemas de Control Analógicos vs. Sistemas de Control Digitales
Un sistema de control Analógico es aquel en que todas las variables son continuas y en función de un tiempo
continuo. En cambio, en un sistema de control digital las variables son discretas, es decir solo pueden tomar
valores predeterminados en función de la precisión del sistema, y estas variables se conocen solo en algunos
instantes de tiempo. Los sistemas de control que usan una computadora son en esencia sistemas digitales.
Variable
Continua Discreta
x(t) x(t)
Continuo
Tiempo
t t
Señal Analógica Señal discreta tiempo continuo
x(t) x(t)
Discreto
t t
Señal Continua tiempo discreto Señal Digital
Sistema Regulador
Es aquel cuya función es mantener la variable controlada en un valor constante
Un ejemplo de este tipo de control es el regulador de Watt, que permite controlar la velocidad de motor,
manteniéndola en un valor constante cualquiera sea la carga.
En este sistema el motor está conectado mecánicamente a un sistema de centrífugas. Si la velocidad
del motor aumenta entonces la fuerza centrifuga ejerce una fuerza que hace cerrar la válvula de paso
de combustible, esto hasta que la velocidad se estabilice nuevamente en el valor deseado. Lo contrario
ocurre si
Lector de Controlador
posición
Centrifuga
Motor
Abre
Cierra
Combustible C
Válvula de
control
Esquema de un regulador de velocidad de Watt Control de dirección de un barco
la velocidad del motor disminuye.
Sistema seguidor
Es cuando la señal de referencia varía constantemente y el sistema de control trata de mantener la
igualdad entre la variable controlada y la señal de referencia. El control trata de seguir el punto de ajuste.
Por ejemplo el control de dirección de un barco.
En un barco grande es imposible llevar el timón en forma manual. Luego el movimiento se realiza mediante un
sistema hidráulico. La rueda del timón indica la posición deseada (señal de referencia que en este
caso cambia continuamente) y es el motor hidráulico el que se encarga de mover la pala.
Se dice que una ecuación es lineal si a esta se le puede aplicar el principio de superposición.
El principio de superposición establece que la respuesta producida por dos funciones excitadoras distintas
es la suma de las respuestas individuales.
Si s = f(e) entonces e = a + b → s = f(a) + f(b)
Es de notar que los sistemas reales son todos no lineales, sin embargo en muchos casos la extensión de
la variación de las variables del sistema no es amplia y se puede linealizar el sistema dentro de un
rango relativamente estrecho de valores de variables.
Sistemas de control una entrada y una salida (SISO) vs. Sistemas de control de múltiples
entradas y salidas (MIMO)
Un sistema puede tener una entrada y una salida, por ejemplo los controles de temperatura vistos
anteriormente. Pero existen también sistemas con múltiples entradas y salidas, por ejemplo el control de
una caldera en donde las variables controladas (salidas) son la temperatura y la presión, las cuales son
controladas a través de un flujo y del calentamiento.
U Y
A
C
U Y V Z
G D
SISO MIMO
Estabilidad
Se dice que un sistema
de control estable es
aquel que responde en
t t
forma limitada a
Estable Inestable
cambios limitados en
la variable controlada.
Es decir si ocurre un cambio en la señal de referencia o se produce una perturbación el sistema al principio
se desviara de su valor y volverá luego a alcanzar el valor correcto.
Un sistema inestable en cambio producirá oscilaciones persistentes o de gran amplitud de la variable
controlada.
Exactitud
Un sistema exacto es aquel VC VC
capaz de mantener el error
en un valor mínimo, o en
todo caso aceptable.
En la realidad no existen
sistemas absolutamente
exactos debido a las
pequeñas imperfecciones
de sus componentes, pero
se consideran sistemas
exactos aquellos que
t t
satisfacen los
Exacto Inexacto
requerimientos del sistema.
Por lo general el costo de un sistema tiende a aumentar con la exactitud.
Velocidad de Respuesta
Es la rapidez con que la variable controlada se aproxima a la señal de referencia.
Un sistema debe
responder a cualquier
entrada en un tiempo VC VC
aceptable, ya que
aunque un sistema sea
estable y tenga la
exactitud requerida, si
este es demasiado
lento no tiene ningún
valor.
Por lo general la
estabilidad y la t t
velocidad de respuesta Alta velocidad de Baja velocidad de
son características que
se contraponen, es decir mientras más rápido sea un sistema mayor será la tendencia a la
inestabilidad y viceversa.
Señal de Controlador
referencia maestro
TC TT
REACTOR
Controlador LC
esclavo TC TT
Agua Condensad
enfriamiento
Control de relación
Es aquel en el que se controla una variable en relación a otra. El objetivo es mantener la relación entre dos
variables en un valor específico. El control se efectúa manipulando una válvula que afecta a una de las
variables mientras que la otra
permanece constante. Producto
Por ejemplo en un control de una Señal de FT
mezcla de componentes se controla referencia
la proporción de uno (o varios) de RS Mezcla
ellos en función del componente
principal para así mantener
constante la proporción. FC
FT
Producto Flujo
En este caso una de las variables es manipulado
la controlada mientras que la otra
se utiliza para generar el valor Control de una mezcla de componentes
deseado. La no controlada es
multiplicada por un coeficiente ajustable en un instrumento llamado estación de relación (RS). La señal de
salida de la estación de relación es la señal de referencia del controlador de flujo.
Control de rango partido
En este caso el controlador tiene una entrada y dos salidas. Una de las salidas actuará desde un valor
A hasta uno B y el otro desde uno B
hasta uno C.
TC TT
Por ejemplo en un control de rango Señal de
partido de temperatura. referencia
Ejercicios
2) Las neveras para mantener el producto a una temperatura dada poseen un sistema de control.
a) Explique como es o como cree que debe que debe ser el sistema de control (Haga un esquema).
b) Haga el diagrama de bloque típico de un sistema de control e identifique los componentes en el
sistema de la nevera.
c) Clasifique el sistema de control.
3) Los motores d ecombustión interna requieren, para funcionar correctamente, de un sistema de
enfriamiento, que mantenga la temperatura del motor en un valor adecuado. Ni muy alto ni muy bajo.
Existen varios sistemas para controlar esa temperatura, los dos mas conocidos son:
El clásico que usa un termostato que regula el flujo de refrigerante (agua), y donde el ventilador
esta acoplado mecánicamente al motor.
En los vehiculos más recientes se usa, además del termostato, un controlador de temperatura que
prende y apaga un electroventilador.
Seleccione uno de los dos sistemas antes mensionados y:
a) Explique su funcionamiento, utilizando figuras.
b) Haga el diagrama de bloque típico de un sistema de control retroalimentado identificando
los componentes del sistema.
c) Clasifique el sistema de control.
Tema 2. Modelado Matemático 25
Sistemas Mecánicos
Un sistema mecánico está conformado por los elementos siguientes:
Elementos Representación gráfica Ecuación fundamental
Resorte F Kx
dx
Amortiguador F CV C
dt
dx
Fricción F BV B
dt
Masa d2x
F Ma M
M dt 2
Donde:
F : Fuerza x : Desplazamiento V : Velocidad
a : Aceleración K : Constante del resorte C : Constante del amortiguador
B : Coeficiente de fricción M : Masa
El modelo matemático se obtiene haciendo un diagrama de cuerpo libre sobre cada masa del sistema.
dx
M F Kx C M d2x
dt dt 2
Fr Fa O escrito ene. Orden común de una ecuación diferencial
ordinaria:
d2
x C Kx F
M dx
dt 2 dt
Esta ecuación es una relación del desplazamiento de la masa (salida) en función de la fuerza aplicada
(entrada).
Para simplificar la escritura de la ecuación diferencial se puede utilizar el operador matemático de derivada:
d
D
dt
d2
Que para una derivada de segundo orden es: D2 2
dt
Con esta representación la ecuación de nuestro sistema mecánico se escribe: MD2 x CDx Kx F
Ejes T G
d
Cojinete T C C
dt
d 2
Masa o Volante de T I I
I
inercia dt 2
1 1 N2
relación de velocidad
Tren de engranes 2 2 N1
T11 T22 relación de trabajo
Donde:
T : Torque o momento : Desplazamiento angular o deformación angular
: Velocidad angular : Aceleración angular
G:
Coeficiente
de
deformación
de ejes C :
Coeficiente
de fricción
viscosa
I : Momento de inercia de masas N : Numero de dientes de engrane
El modelo matemático se obtiene haciendo un diagrama de cuerpo libre sobre cada volante de inercia del
sistema.
Ejemplo 2:
Se hace el diagrama de cuerpo libre sobre el volante de inercia:
G1 1 d2 d 2 2
G22 C I
T dt dt 2
G2
C
G1 2
I Y se escribe adicionalmente la ecuación que relaciona el
momento aplicado con al extremo del eje con el que momento que
recibe el volante de inercia:
T G1 1 2
2 1 Con estas dos ecuaciones se puede hallar una expresión entre el
momento aplicado al sistema (entrada) y el movimiento angular
del momento de inercia:
d 2 2 d 2
I C G22 T
dt 2 dt
O escrito utilizando el operador matemático: ID 22 CD2 G22 T
También se puede definir una salida diferente por ejemplo el desplazamiento angular en el extremo del eje. En
este caso se combinan las dos ecuaciones de una forma diferente:
Sistemas Eléctricos
Un sistema eléctrico está conformado por los elementos siguientes:
Elementos Representación gráfica Ecuación fundamental
Resistencia V RI , ZR R
1 1
Capacitor V t Idt; Z
C
C 0 CD
dI
Bobina V L , Z LD
L
dt
Elemento cualquiera V ZI
Z
En un nodo
I 0
En una malla
V 0
Elementos en serie ZT Z i
Z Z
Elementos 1
Z Z
T
1Zi
en paralelo
Z
Donde:
V : Voltaje o diferencia de potencial
I : Intensidad
Z : Impedancia
R : Resistencia
C : Capacitancia
L : Inductancia
Ejemplo 3: Hallar V f I y V f V
3
R2
R1
Primera parte V f I
V1 C2
Sabemos inicialmente que:
V2
V ZT I
Donde
ZT Z1 Z 2 Z3 Z 4 V C R
V3 L 3
Z R; 1
Z ; V4
1 1 2 1 C 2D L4
R2
1
Z3 ; Z4 L4 D
1 C 3D 1
R3 L3 D
Luego
V R 1 1 L D I
1 4
1 C 2D 1 C 3D 1
R2 R3 L3 D
Segunda parte V f V3
Hallamos primero I f V3
V3 Z3 V3
I
I Z
3
Luego
V ZT V
Z3 3
Analogía Electromecánica
Este es un método que permite resolver en forma relativamente más sencilla problemas mecánicos, como i
se tratase de sistemas eléctricos. En este caso hacemos:
V análogo a F, e I análogo a x
Resorte ZR K
Amortiguador Z A CD
Fricción ZF BD
Z M MD 2
Masa
M
Elemento cualquiera F Zx
Z
Elementos en serie ZT Z i
Z1 Z
Elementos 1
Z
Z1 T
1Zi
en paralelo
Z
El método sirve para sistemas con una sola fuerza y se resuelven los problemas haciendo primero el
diagrama de impedancias.
Diagrama de impedancias
Para realizar el diagrama de impedancias:
Se coloca en la parte superior una línea horizontal que representa la coordenada donde
está aplicada la fuerza.
Se coloca en la parte inferior una línea que representa la tierra, o referencia.
Se colocan entre las dos anteriores líneas que representen las otras coordenadas existentes.
Se colocan las impedancias correspondientes a cada elemento y se hace la conexión de este a las
coordenadas correspondientes. Nótese que cada elemento estará conectando siempre dos
coordenadas. En el caso de las masas estas siempre irán conectando la tierra y la coordenada
donde se encuentran, mientras que los otros elementos pueden conectar dos coordenadas
diferentes a la tierra.
Ejemplo 4: F F
M
La ecuación del sistema obtenida por x
éste método será:
F Z R Z A Z M x ; ZM ZA ZR
F K CD MD 2 x
F Kx CDx MD2 x
T1 T2
Si T1 T2 : Q T1 T2
Pared delgada Q
(no absorbe calor) Rt
Rt
dT
Pared gruesa
T1 Tp T2 Q C T
dt
(con Q1
almacenamiento de Q2 T1 Tp Tp
CT Q T;2 Q
calor)
Rt1 Rt2 1
Rt1
2
Rt 2
Q1 Q2 CT DTp
Donde:
Q : Flujo de calor T : Temperatura
RT : Resistencia térmica CT : Capacitancia térmica (masa por calor
específico)
Ejemplo 5: Termómetro de mercurio con pozo térmico de cobre.
El termómetro está formado de tres paredes que absorben calor, más un
elemento receptor que también absorbe calor:
Hg: Mercurio THg, CHg
V: Vidrio TV, CV
C: Cobre TC, CC
T a D3T a2 D a a T
T DT
E 1 Hg 2 Hg 3 4 Hg
Hg
Con:
a1 R1R3CHg CV
CC R3CHg R3CHg CC
CV
R3CHg CV
a2 R1 CHg CC CV CC
R1 R2 R2
R3CHg C R3CHg CHg
Hg CC CC
CV CV CC
a3 R1 R R R R R2 R R
R3CHg
1 2 1 1 2 2 2 2
R3 R3
1 1
a4 R1 R R 2 1
1 2 R2
Sistemas Hidráulicos
Un sistema hidráulico está conformado por los elementos siguientes:
Elementos Representación gráfica Ecuación fundamental
Escuela de Ingeniería Mecánica - ULA
Tema 2. Modelado Matemático 33
dP
Tanques
Q C ; Q Q C DP
Qe Ch Qs
h
dt
e s h
h
P
P h
Q
Pe Ps
Ductos P P Ps
Q ;Q e
Rh Rh Rh
Donde:
Q : Flujo o caudal
P : Presión
h : Nivel
Rh : Resistencia hidráulica (perdidas que se producen en tuberías y accesorios)
Ch : Capacitancia hidráulica (volumen que es capaz de absorber)
Ejemplo 6: Hallar h2
f Qe
Qe
Ch1
h1 Ch2
R1 h2 R2 Qs
P1 P2
Q
El sistema hidráulico está conformado por dos tanques conectados entre sí. Estos tienen una entrada
de agua por el primer tanque y una salida por el segundo. Las ecuaciones fundamentales del sistema,
considerando presiones manométricas ( Patm 0 ) serán:
(1) P h
1 1
(2) P2 h2
(3) Qe Q
(4) C DP Q Q
h1 1 s
(5) Ch 2 DP2
P1 P2
Q
R1
(6) Qs P2
R2
1 1
1 1
Con 6 en 8 obtenemos P1 P2 : (9) P R C D P2
1 1 h2
f R R
1 2
1 1 1 P2
Con 9 en 7 obtenemos Qe f P2 : (10) Qe R C D R C D P2
h1 1 h2
R R
1 1 2 1
2 R1Ch1 1
De R1Ch1Ch 2
Q P2 Ch1 Ch 2 DP2 P2
R2 R2
2 R1Ch1
Con 2 en 10 obtenemos Qe f h2 : Qe R1Ch1Ch 2D h2 Ch1 Ch 2 Dh2 h2
R
2 2
Sistemas Neumáticos
Un sistema neumático está conformado por los elementos siguientes:
Elementos Representación gráfica Ecuación fundamental
dP
Tanques m e Cn P m s m C ; m s Cn DP
m
n e
dt
m
Pe Ps P P
Ductos m ; m e
Ps
Rn Rn
R
n
Donde:
m : Flujo másico
P : Presión
Rn : Resistencia neumática (perdidas que se producen en tuberías y accesorios)
Cn : Capacitancia neumática (V RT )
Jean-François DULHOSTE
Tema 2. Modelado Matemático 35
Nota: m 2 debe ser conocido (entrada), o en su defecto la presión de entrada debe ser conocida.
R
C1C2 1 1 3 1
2
C1 2 1
C1
m 2 D P1 C2 DP1 R P1
2
R R R
1 1 1 3 1 3
1. Sistema mecánico
Hallar:
K1 x2 f F y
K2 C2 x1 f F
F M1
Por los dos métodos
C1 x1 x2
M2
2. Sistema termo-neumático
R2 R3 Hallar:
R1 P2 Patm
P3 m 3 m 3 f TE , P4
C2 R4
m1 m 2
C3
P4 = Cte Nota: la ecuación de
m 4 relación entre los dos
P1 sistemas:
Gas Pv mRT
P RT
Donde se supone:
Vidrio
R Constante
Cobre TE
3. Sistema termo-eléctrico
R3
QR VI I 2 R V 2 R
4. Sistema hidráulico
QE1
Ch3
Ch1
h1 Ch2
R1 R2 h2 R3
P1 P2
h3
QS a
b R4
QE2
F
P3 P2
Pistón
de área m1
C3 C2
m
R1 Patm
P1 m
R3
2
R2
Hallar m f F
Nota la ecuación que relaciona el sistema neumático
con el pistón es: P F A QE
n1
X s a0 s an1s an
a1s
Utilizando este concepto de función de transferencia se puede representar la dinámica de un sistema por
ecuaciones algebraicas en s. Si la potencia más alta de s en el denominador es n se dice que el sistema es
de orden n.
Sea un sistema LTI, SISO sometido a una entrada u(t) y representado U(s) Y(s)
G(s)
por su función de transferencia G(s).
Definición de la respuesta al impulso: Un sistema que tiene como función de transferencia G(s), tiene como
respuesta al impulso la función:
g(t) = f–1[G(s)]
La respuesta de este sistema a una entrada cualquiera u(t) se puede calcular utilizando el teorema de
convolución: La respuesta de un sistema cual función de transferencia G(s) esta dado por la siguiente
integral de convolución:
t t
y(t) = ƒ g(r)u(t — r)dr = ƒ g(t — r)u(r)dr
0 0
El producto de convolución se expresa en general como y(t) = g(t) ∗ u(t) = u(t) ∗ g(t)
La Matriz de Transferencia
Uj(s) Yi(s)
El concepto de Matriz de Transferencia es una extensión a sistemas MIMO G(s)
de la función de transferencia.
Definición: la matriz G(s) ∈ Œr×m se denomina Matriz de Transferencia, y relaciona la entrada U(s) con la
salida Y(s).
Y(s) = G(s)U(s)
Con:
m : Número de entradas; r : número de salidas
Que también puede expresarse en notación matricial explícita por elemento:
Yi(s) = Gij(s)Uj(s)
Ejemplo 2: Se tiene el sistema mecánico MIMO con dos entradas (F1 et F2) F1
y dos salidas (y1 et y2): M1
Los polos y los ceros permiten la caracterización dinámica de un sistema. Estos se pueden definir a partir de
funciones o matrices de transferencia (mas fácil para los sistemas SISO) o a partir de modelos de
estado (mas practico en modelos MIMO).
Ejemplo 3:
Para la función de transferencia del sistema mecánico
N(s) del ejemplo
1 1:
G(s) = =
D(s) MS2 + Cs + K
G(s) = [s+1
1
s+2
1
]
s+1 s+1
Los polos del sistema serán las raíces de esta ecuación característica:
s1 = —1; s2 = —1; s3 = —2
Los ceros
En el caso de los sistemas denominados cuadrados, en donde el número de entradas es igual al número de
salidas, los ceros se pueden determinar mediante la matriz o función de transferencia.
Se define el polinomio o ecuación de los ceros N(s) como el más grande común divisor de los numeradores
de los menores de orden máximo de G(s) normalizado, para tener la ecuación característica D(s) como
denominador. Este polinomio se obtiene con:
N(s)
= |G(s)|
D(s)
Les ceros (z0) son las raíces de este polinomio de orden m, y se obtienen para: N(z0) = 0
2
3
Ejemplo 5: para la misma matriz de transferencia del ejemplo 4: G(s) = [ s+1
1 s+2
]
s+1 1
s+1
Existe un menor máximo de orden 2 que
2 es: 1 3
|G(s)| = — 1 —s + 1
=
(s + 1) (s + (s + 2) (s + (s + 1)2(s + 2)
El polinomio de los ceros es: 1) 1)
N(s) = —s + 1
Existe un solo cero: z0 = 1
NOTA:
Todas las definiciones son aplicables al caso más simple de un sistema SISO para las cuales G(s) es una
fracción racional donde:
El numerador es πz(s) y sus raíces son los ceros
El denominador es πs(s) y sus raíces son los polos.
Definiciones:
La diferencia de grados entre D(s) y N(s) (n y m) se denomina el grado relativo.
Si n — m > 0 el modelo es estrictamente propio (grado relativo positivo)
Si n — m = 0 el modelo es bipropio (grado relativo cero).
Si n — m ≥ 0 el modelo es propio.
Si n — m < 0 el modelo es impropio (grado relativo negativo).
Los sistemas reales son estrictamente propio.
Los controladores pueden ser propios o impropios. Los impropios se modifican para
poder construirlos.
La representación en espacio de estado, también conocida como representación interna, fue utilizada en
otras disciplinas como la mecánica o termodinámica desde hace largo tiempo. Por ejemplo, para el
comportamiento macroscópico de un gas puede describirse y predecirse con un número finito de variables
físicas: el volumen V de ese gas, su presión P y su temperatura T. El conjunto P, V, T representa el estado
termodinámico del gas. Su evolución en el tiempo dependerá del entorno exterior (aporte de calor por
ejemplo) pudiéndose caracterizar su comportamiento dinámico con el conocimiento de ese entorno, que en
control denominamos entrada del sistema.
En conclusión el estado dinámico de un sistema puede ser representado por un conjunto de variables
denominadas variables de estado. Este conjunto de variables caracteriza completamente la configuración
dinámica actual del sistema. Para esto se requiere de un número mínimo de variables de estado necesarias
y suficientes que permiten la descripción dinámica del sistema.
Los sistemas automáticos modernos, a partir de los cuales se desarrollo la representación de estado
para el control de procesos, aparecen en los años 60 para permitir el control de sistemas complejos tales
como las aplicaciones espaciales Apolo y Polaris, las cuales tienen múltiples entradas y salidas (MIMO), y
criterios de funcionamiento cada vez más severos. El uso del espacio de estado para representación de
sistemas de control proviene de la capacidad que tiene esta representación de representar sistemas
multivariables complejos. Su desarrollo y aplicación crece luego con el uso de los computadores.
El conjunto de variables de estado no es único, pero debe estar conformado para cada sistema por un
número idéntico de variables de estado independientes. Esto significa que la selección de estas variables,
así como de sus condiciones iniciales, constituye un conjunto que se puede fijar de forma arbitraria.
El estado inicial del sistema constituye su memoria: dado un estado inicial a un instante dado el
conocimiento del pasado no permite el conocimiento del futuro del sistema, se requiero por lo tanto de
unas funciones (ecuaciones de estado) que permiten la predicción del futuro, las funciones comúnmente
utilizadas son las resultantes de una integración.
Para comprender correctamente el funcionamiento de esta representación se estudiaran las definiciones
básicas de estado, variable de estado, vector de estado y espacio de estado. Luego se presentará la forma
de las ecuaciones en espacio de estado, su relación con las funciones de transferencia y la forma de
representar sistemas lineales en espacio de estado.
Definiciones
Estado
El estado de un sistema dinámico es el conjunto más pequeño de variables x(t) (denominadas variables de
estado) tales que el conocimiento de esas variables en t t0 , conjuntamente con el conocimiento de la
entrada u(t) para todo tiempo t t , y las ecuaciones que describen la dinámica f(x, t, u) y g(x, t, u),
0
determinan completamente el comportamiento futuro de los estados x(t) y salidas y(t) del sistema
para cualquier tiempo t t0 .
Variables de estado
Las variables de estado de un sistema dinámico x(t) son las variables que constituyen el conjunto más
pequeño de variables que determinan el estado de un sistema dinámico.
Nótese que las variables de estado no deben ser necesariamente cantidades físicas mensurables u
observables. Sin embargo es conveniente escoger como variables de estado de un sistema magnitudes.
Vector de estado
Si se requieren n variables de estado para describir completamente el comportamiento de un sistema dado,
se puede considerar a esas n variables como los n componentes de un vector x. Vector que recibe el nombre
de vector de estado.
y a1 y an1 y an y
n1
u
Suponiendo que las condiciones iniciales
y 0 , y 0 , y y la entrada ut para un tiempo t 0 son
0
conocidas, entonces las variables de estado deben ser tales que definan completamente el comportamiento
futuro del sistema. Bajo esta premisa se puede entonces escoger como variables de estado:
n1
x1 y; x2 y ; xn y
Entonces la ecuación diferencial se puede escribir como:
x1
x2 x2
x3
⁝
xn1 xn
xn an x1 a1 xn
u y x1
O en forma matricial:
x Ax Bu
y Cx
Donde:
x 0 1 0 0 0
1
x 0 0 1 0 0
2
x ; A ⁝ ⁝ ⁝ ⁝ ; B ⁝ ; C 0 0
⁝
1
0 0 0 1 0
x
n an an2 a1 1
an 1
Esta forma de representación se denomina comúnmente forma canónica controlador.
Mx2
Cx2 Kx1 u x u Cx2 Kx1
2 M
Se obtiene entonces el sistema de ecuaciones de estado:
x1 x2
K
x C 1
x x u
2 1 2
M M
y x1
M
El cual puede expresarse matricialmente como:
x Ax Bu
y Cx
Donde:
0 1 0
A=[ K C ] ; B = [ 1 ] ; C = [ 1 0]
— —
M M M
Caso de ecuaciones diferenciales de orden n en donde la función excitadora incluye términos derivativos
Sea el siguiente sistema de orden n:
n n1 n n1
Suponiendo que las condiciones iniciales y 0 , y 0 , y 0 y la entrada ut para un tiempo t 0 son
conocidas, entonces las variables de estado deben ser tales que definan completamente el comportamiento
futuro del sistema. En este caso en particular las variables de estado deberán además ser tales que
eliminen las derivadas de u en la ecuación de estado. Bajo esta premisa se pueden escoger como variables
de estado:
x1 y 0u
x2 y 0u 1u x1 1u
x3 y 0u 1u 2u x2 2u
⁝
n1
n n2
1
xn y 0 u 1 u n2u n1u xn1 n1u
Donde 0 , 1, n son coeficientes que se determinan como:
0 b0
1 b1 a10
2 b2 a11 a2 0
3 b3 a12 a2 1 a30
⁝
n bn a1n1 an11 an0
Con esta escogencia de variables de estado se obtiene el sistema de ecuaciones de estado:
x1 x2
1u x2 x3
2u
⁝
xn1 xn n1u
xn an x1 an1 x2 a1 xn
nu y x1 0u
O en forma matricial:
x Ax Bu
y Cx
Donde:
x 0 1 0 0 1
2
0
1
x 0
1 0
x ; A ⁝ ⁝ ; B ⁝ ; C 1 0 0 ; D 0 b
2
⁝ ⁝
⁝ 0
0 1
n1
x 0 0
n an an2 a1
n
an 1
x2 1 x1 1 x2 1 x2 1 x3
R2 CV R2 CV R3 CV R3 CV
1 1
x x x
3 2
R3CHg R3CHg 3
y x3
El cual puede escribirse en forma matricial como:
x˙ = Ax + Bu
y = Cx
Donde:
1 1 1
RC RC R 0 1
C
1 C 2
1 2 C 1 1
C
1
A ; B 0 ; C 0 1
0
R RC R RC 0
C C
2 V 3
2 V 3 V
V 1 1
0
R3CHg R3CHg
Ejemplo 9: Consideremos el sistema mecánico MIMO con dos entradas (F1 y F2) y dos salidas (y1 y y2),
mostrado en la figura:
Las ecuaciones del sistema son:
F1 — K1(y1 — y2) — C1D(y1 — y2) = M1D2y1
F1
F2 + K1(y1 — y2) + C1D(y1 — y2) — K2y2 — C2Dy2 = M2D2y2
M1
En este caso tenemos:
u1 y1
F1 y1
u = [u ] = [ ] ; y = [y ] K1 C1
2 F2 2 F2
Por lo
tanto:
Gs C sI A B D
1
Ejemplo 10: para el sistema mecánico teníamos que el modelo matemático expresado en espacio de estado
es:
x Ax Bu
y Cx
Donde:
0 1
0
A K
B
C; 1 ; C 1 0
M M M
Si queremos obtener la función de transferencia a partir de esta expresión del modelo debemos entonces
usar la expresión:
1
0 0 1 1 0
G s C sI A B D 1 0 s C 1 0
1
K
0 1
M M
Al resolver esta ecuación obtenemos:
s 1 1 0 M
Gs 1 0 K C 1
s M M
M
Recordatorio : Donde:
A–1 = 1 adj(A)T dij = (—1i+j )|A˜ij | elementos de la matriz adj(A)
|Æ|
A˜ij menores principales
C 1 0
s 0
Gs Gs
1 C
1 0
1 1 s 11
M
2 C K K s 2 C K
s M
s M s s M
M M M
M ;
Gs
1 1 Gs M
1
C K M M
s
2
s ;
M Ms2 Cs K
Que es exactamente la función de transferencia encontrada a partir de la ecuación diferencial.
Ejemplo 11: A partir del modelo matemático en representación de estado obtenido en el ejemplo 8, queremos
obtener la función de transferencia de éste sistema se puede entonces obtener con la expresión:
G s C sI A B D
1
1 1 1
s 0
1
RC RC RC
Gs 1
0 1 C 2 1 2 C 1
0 1 C s 1 0
1
R RC R
C C RC
2
2
V
V 3 3 V
0
V 1 1
0 s
Con: R3CHg
R3CHg
s 1
1
1 s
RC
R2C 3 V R3 CHg
1 s
1
1
1
V
RC R C
1 1 2 V
2 V
R3 R3
Hg
C
Hg
C R R C C
3 V 3 Hg
1 1 1 1 1 1 1
1
s s
s s
R C R C
R C
3 Hg R1CC R C R3 CHg
R C
1 C
R C
2 C 3 Hg
2 C 2 C
1
1 1 1
C s R RC
R C C s R
1
1 s 1 1 1 1 C 2 C 2 V 3 V
R C
R2CC 3 V R2CC 3 V 1 1
RC
R1CC
R2CC R2CV
sI A1 1
1 1
s 1 1
s R C C R C s R C
1 C 2 C 2 V 3 V 3 Hg
1 1 1 1 1 1 1
R C R C s C RC C s RC RC
2 C 2 V 3 Hg 3 V 3 Hg 1 C 2 C
En este caso se observa que el cálculo algebraico se vuelve relativamente largo, por lo cual es más
fácil obtener la función de transferencia directamente de la ecuación diferencial del sistema:
T 2 a D3T a a aT
DT DT
E 1 Hg 2 Hg 3 4 Hg
Hg
Con:
a1 R1 R3CHg CV
CC R3CHg R3CHgCC
CV
R3CHgCV
a2 R1 CHg CC CV CC
R1 R2 R2
R3CHg CHg R3CHg CHg CC CC
CV CV CC
R R R2 R R
a3 R1 R R R3CHg
1 2 1 1 2 2 2 2
R3 R3
1 1
a4 R1 2 1
R
R1R 2 2
La no unicidad del conjunto de variables de estado significa que para un sistema cualquiera existen diversas
representaciones de estado posibles. De forma recíproca matrices A diferentes pueden representar un
mismo sistema y por ende una misma ecuación característica.
Para probar que esto es posible utilicemos el ejemplo siguiente:
Ejemplo 12: Supongamos que inicialmente se tiene el sistema en forma de ecuación diferencial:
y 6 y 11y 6 y 6u
Para obtener una representación en forma de espacio de estado se pueden tomar los siguientes estados:
x1 y; x2 x3 y
y;
Con estos estados se obtiene la representación de estado siguiente:
x1
x2 x2
x3
x3 6x1 11x2 6x3 6u
Que puede expresarse en forma matricial como:
x Ax Bu
y Cx
Donde:
0 1 0 0
A=[ 0 0 1 ]; B = [0]; C = [1 0 0]
—6 —11 —6 6
1 0 0 1
F(s + 1) I 3
s+1 0
0 –1 I
3 I 1 I
([(s) = [1 1 1]
G
0 s+2 0 ]) [—6] = [1 1 1] 0 0 [—6]
0 0 s+3 3 I (s + 2) I
I 3
1
1 1 I
G(s) = [ 0 0 I
L (s + 3)ے
1 3 =
] [—6] 3 6
— 3 +
(s + 1) (s + 2)
(s + 3) (s + 1) (s + 2) (s + 3)
3
6
G(s) =
(s + 1)(s + 2)(s + 3)
Se observa que esta segunda representación de estado corresponde exactamente al mismo sistema
pues posee la misma función de transferencia.
Forma canónica de Jordan
La forme canónica de Jordan o forma Modal, es la correspondiente a la segunda representación del ejemplo
anterior, en la cual la matriz A solo posee elementos en la diagonal, es una matriz llamada diagonal. Estos
elementos de la diagonal corresponden directamente a los valores propios de la matriz, los cuales son los
polos del sistema o raíces de su ecuación característica.
1 0 0
AJ 0 2 0
0 0 3
La ventaja de la representación canónica de Jordan es que muestra directamente la estabilidad del sistema
y además las operaciones matemáticas con las matrices diagonales son más sencillas.
La representación canónica de Jordan es solo otra representación en espacio de estado posible para
un sistema lineal.
Diagonalización de matrices
De hecho la representación de Jordan es un caso particular de cambio de variable, donde la matriz P que
permite la diagonalización de la matriz A. La forma de la matriz P necesaria para la obtención de la
representación canónica de Jordan se puede generalizar par una matriz A cuya forma es de tipo canónica
controlador:
0 1 0 0
0 1 0
⁝ ⁝
0
A ⁝ ⁝ 0 1
0
0
an an2 a1
1
Como: an
1 1 1 1
2
1
3 n
P2 2 2 2
1 2 3 n
⁝ ⁝ ⁝ ⋱ ⁝
1 n1 2 n1 3 n1 n
n1
Para nuestro ejemplo la matriz P que permite la obtención de la forma canónica de Jordan es:
1 1 1
P 1 2 3
1 4 9
Cuya inversa
es:
3 2.5 0.5
1
P 3 4 1
1 1.5 0.5
De hecho:
3
2.5 0.5 0 1 0 1 1 1 0 0
1
P1 AP 3 4 1 0 0 1 1 2 3 0 2 0 A
J
Nota: Para que una representación de estado sea equivalente a otra obtenida por cambio de variable
es necesario la transformación de todas las matrices que representan el sistema.
Ejercicios
1. Para los ejercicios del tema 2 obtener una representación en espacio de estado de las ecuaciones
del modelo, utilizando la ecuación diferencial obtenida y otra utilizando el sistema de ecuaciones que
representa el sistema.
2. Obtener la función o matriz de transferencia de los sistemas utilizados para el ejercicio 1.
3. Para los siguientes sistemas en representación de estado obtener la función de transferencia, la
ecuación diferencial del sistema y la forma canónica de Jordan.
0 1 0 0
A 0 0 1 B 0
1 8 1 C 1 0 0
2
4
, ,
0 1 0 0
A
0 0 1 B 0
2 1
2 2 C 1 0 0
, ,
0 2 0 0
A 0
0 1 B 0
23 C 1 0 0
7 7
,
3
,
Tema 4. Diagramas de Bloque 59
1. B
loques
Entrada Salida U(s) Y(s)
Sistema G(s)
Es el elemento principal de un diagrama de bloques ya que representa los componentes del sistema. Es el
elemento que contiene la relación entre la entrada y la salida del componente del sistema, ya sea con un
nombre o directamente la función de transferencia del elemento.
2.Líneas
Son las que representan la trayectoria de las señales que transportan la información y energía. Dicho
de otra manera indican la conexión de una variable entre diversos bloques. Se usan por lo general flechas
para indicar el sentido de la señal.
3.Punto de suma
U Y U Y
V
V
Este se representa por un círculo e indica que la salida es igual a la suma de las dos señales que entran, o la
resta en caso de que aparezca el signo negativo:
Z X Y
Existen dos nomenclaturas posibles, ya sea con los signos al exterior, en cuyo caso la ausencia de signo
indica una suma, o separando el círculo en segmentos mediante una X y asignando el signo dentro del
segmento correspondiente.
X X X
X X
X
X
Es el lugar donde la señal se separa para seguir dos o más trayectorias. En otras palabras indican que
una variable se usará en varios bloques.
Elemento de medición
Señal de retroalimentación o
medición
sint
U 1 Y
Ki KV
p
K K s sd Ms2 Cs K
KM
X H
G Y
U
1
GH
Álgebra de bloques
Los diagramas de bloques de sistemas de control complicados se pueden simplificar usando una serie
de teoremas de transformación, las cuales se obtienen fácilmente por deducción del significado de los
elementos.
Teoremas de transformación
Diagrama de Bloque Diagrama de bloque equivalente Ecuación
1 Combinación de bloques en cascada
U
U V Y GG Y V G1U ; Y G2V
G1 G2
Y G1G2 U
1 2
U V Y V G1U ; W G2U
G1 U
G ±G Y Y V W
W 1 2
G2 Y G1 G2 U
3 Retroalimentación negativa
U X Y Z G2Y ; Y G1X
G1
U G1 Y X UZ
Z 1 G1G2 G1
G2 Y U
1 G1G2
4 Retroalimentación positiva
U X Y Z G2Y ; Y G1X
G1
U G1 Y X UZ
G1
Z G2 1 G1G2 Y U
1 G1G2
5 Movimiento de un punto de suma después de un bloque
U X Y U X Y
G1 G1 X G 1U ; Z G 1V
V V Z Y G1(U V )
G1
Y U X Y
U G1 X G1 X G1U
V 1
V Y G 1U V
G1
U Y
U Y G1 Y G1U
G1
1
U 1 U U Y
G1 G1
Y U Y
U G1
G1
Y Y Y G1U
G1
U X Y U Z Y X U V; Z U
WY X W ZV
V W W V
Y U V W U W V
10 Redistribución de puntos de ramificación
Y Y
H1 G2
H2
Para reducir el diagrama se pueden seguir los siguientes pasos:
1. Numerar todos los puntos de suma y ramificación:
U 1 2 56 Y
G1 G43 4 G3
H1 G2
H2
2. Reducir desde lo más interno, por ejemplo entre 2 y 3, y entre 4 y 5:
U 1 2 3 4 5 6 Y
G1 G4 G3 + G2
H1
H2
3. Llevar el diagrama a la forma canónica de un sistema de control retroalimentado:
U 1 G1G4 G2 G3 6 Y
1 G1G4 H1
H2
M
Igualamos primero W = 0 y reducimos el diagrama resultante:
U Y1
C E Q
U CEQ Y1
M M
U CEQ Y1
1
MCEQ
Igualamos primero U = 0 y reducimos el diagrama resultante:
W
Y2
C E Q
W Y2
Q
M CEM
W Q Y2
1
MCEQ
Finalmente el diagrama resultante será la suma de los dos diagramas obtenidos:
CEQ Y1
U
1 MCEQ
Y
Q Y2
W
1 MCEQ
U CE
Q Y
1 MCEQ
W
C
V Z
D
Primero determinamos a Y ignorando a Z:
U Y
A
C
V
D
Igualamos V a cero y determinamos la función de transferencia entre U y Y1 :
U Y1
A
BC
D Y1
U
A
BCD
U Y1 A
1
ABCD
Igualamos U a cero y determinamos la función de transferencia entre V y Y2 :
Y2
A
B
V Y2
C —ACD
V
D B
V Y2 - ACD
1- ABCD
V - ACD 1- ABCD Y2
Segundo determinamos a Z ignorando a Y:
U A
BC
V Z
D
U —ABD Z1
Z1
C
U Z1 - ABD
1 ABCD
A
BC
—ABC
V Z2 V Z2
D D
V Z2 D
1
ABCD
Para la salida Z tenemos entonces:
U - ABD 1 ABCD Z1
Z
V D Z2
1 ABCD
Ejemplo 4. Representación en
diagramas de bloque a partir del
Válvula de
esquema del sistema de control
control
retroalimentado.
H 2O KT
Elemento de medición
Donde:
Qe : caudal de entrada al
tanque Qs : caudal de salida del
tanque h: nivel del tanque
Pm : presión medida que corresponde directamente al valor del nivel
Pr : presión de referencia, la deseada en el tanque, que tiene una correspondencia directa con el nivel
deseado
xv : posición del vástago de la válvula que regula el paso del flujo de entrada.
Observamos aquí que necesitamos conocer la función de transferencia del proceso para completar el
diagrama de bloque.
Esta función de transferencia corresponde al de un sistema hidráulico, donde la entrada corresponde a la
diferencia entre el caudal de entrada y el de salida (Qe — Qs), y la salida del sistema será el nivel (h).
Se tienen las siguientes ecuaciones:
Qe Qs (1) PH (2)
CDP 2Oh
Con estas dos ecuaciones podemos obtener la ecuación de la dinámica del sistema que es:
Qe Qs H 2OCDh
A partir de esta ecuación podemos obtener la función de transferencia:
hs 1
Qe Qs s H 2OCs
El diagrama de bloques completo del sistema será entonces:
Qs
P 1 h
K P KV
r
H 2OCs
H 2O KT
Para obtener la función de transferencia de lazo cerrado debemos reducir el diagrama de bloques del
sistema, que es un diagrama de dos entradas (Pr, Qs) y una salida (h).
Igualamos primero Qs a cero y obtenemos:
Pr K P KV h’
KP KV
H 2OCs Pr H 2OCs h’
KP KV
H 2O KT 1 K
H 2O T
H 2O Cs
Preferencia KP KV h’
H 2O Cs KT KP KV
Igualamos luego Pr a cero y obtenemos:
Qs 1 h’’
H 2OCs 1
Qs H 2OCs h’’
1 H 2O KT K P KV
H 2O KT K P KV H 2OCs
Qs 1 h’’
H 2O Cs KT KP KV
El diagrama de bloques simplificado será:
Pr K P KV
1 h
H 2O Cs KT KP KV
Qs
Este diagrama de bloque se puede entonces expresar como una función de transferencia o ecuación en
transformada de Laplace que relaciona la entrada con la salida:
h
KP
P
1
Q
KV s
KT
Cs K
KVP r KT
Cs KP KV
H 2O H 2O
Ejemplo 5: Representación del diagrama de bloques a partir de las ecuaciones de los componentes
Las siguientes ecuaciones representan un sistema de control:
C1 C6
xv Kz C4 y
1D q C8w
v e; z C2u; q y
e AD x
; ; C5v;
2 KC C
7 3
Donde:
e, q, v, w, x, z: son señales
y: es la variable controlada
u: es la señal de referencia
D: es el operador diferencial (derivada respecto de t)
A, C1, C2, C3, C4, C5, C6, C7, C8, K: son constantes (ganancias)
Determinar la función de transferencia de: Identificar:
Proceso Variable manipulada
Controlador Perturbación
Elemento final de control
Elemento de medición
u C2 z K K C7C3
w C6 y C4
C8 1 D K C7C3
q
Luego se vuelve a dibujar el diagrama de bloques pero en una forma más ordenada que permita identificar
más fácilmente los componentes, con u como entrada y y como salida:
w
C8
K q
u z x 1 e Cv C6
C2 C5 y
K C7C3 2
1
1 D
AD
v
C4
K C7C3
Comparando con el diagrama de bloque típico de un sistema de control retroalimentado se pueden
fácilmente identificar los componentes, se hacen determinan entonces las funciones de transferencia para
cada uno de ellos:
w
C8
u C2 K x C1 C6
v Cq y
K C7C3 2 AD C1 5 1 D
C4
K C7C3
w
u C2 K x C1
v q 1 C6
K C7C3 2 AD C1 C5 C8 C8
y
1 D
K C7C3 C2 K C4
K C7C3
w
u C2 KC1 v C5 q C8C6 y
1 D
K C7C3 2AD C1 C8
Controlador Elemento final de Proceso
control
C4
C2 K
Elemento de medición
De esta última forma del diagrama observamos entonces que la variable manipulada es q y que la
perturbación es w.
Tema 5. Respuesta de Sistemas
Introducción
El modelo matemático dinámico de un sistema físico está representado por una ecuación diferencial.
La ecuación diferencial suele ser en general de cierta complejidad, pero en muchos casos prácticos el
modelo puede simplificarse y transformarse en una ecuación diferencial ordinaria lineal de la forma:
n n1
Tipos de excitación
La función excitación de un sistema puede tomar cualquier forma en la realidad, y por lo general no puede
conocerse con anticipación ya que suele ser de carácter aleatorio. Sin embargo en muchos casos prácticos
el estudio de algunas funciones de excitación específicas suele ser suficiente, ya que puede aportar
conocimientos importantes sobre el sistema.
Las funciones de excitación más representativas son:
Entrada en escalón
En este caso la entrada sufre un cambio instantáneo y u
finito.
Esta se puede expresar de la siguiente forma:
0 para t 0
ut
para t 0
H
O expresada en transformada de Laplace:
H
U s H
s
Donde:
H: es una constante
Cuando H es igual a 1 se dice que la entrada es un
0 t
escalón unitario.
Nótese que la función escalón queda indefinida para t = 0.
La función escalón es la más simple de las entradas que se le puede aplicar a un sistema, y por esto
suele ser la más utilizada porque permite hacer análisis simples de la respuesta del sistema. Esta
permite
Tema 5. Respuesta de Sistemas 75
principalmente definir la respuesta de sistemas en donde los cambios son instantáneos y que luego se
mantienen en el tiempo, por ejemplo la respuesta de un termómetro clínico al ser introducido en un baño a
una temperatura dada.
Entrada en rampa
En este caso la entrada varía linealmente con el u
tiempo.
Esta se puede expresar de la siguiente forma:
0 para t 0
ut
H para t a
0
O expresada en transformada de Laplace: b
H
U s
s2
En este caso H define la pendiente de la rampa:
a
H
b 0 t
La función rampa suele ser también muy utilizada ya
que permite el estudio de la respuesta de un sistema sometido a una entrada que varía
constantemente. Por ejemplo la respuesta de un motor al aumentar progresivamente la inyección de
combustible.
Entrada sinusoidal
En este caso la entrada cambia su magnitud de u
acuerdo a una función sinusoidal de amplitud
constante 2H y frecuencia ω.
Esta se puede expresar de la siguiente forma: H
ut 0 para t 0
H sint para t 0
O expresada en transformada de Laplace: H
H
U s
s2
2
Entrada en pulso
u
En este caso la entrada sufre un cambio instantáneo y
finito, pero de duración finita en el tiempo. Esta se
puede considerar como una entrada en escalón de
altura H/t0 pero de duración finita en el tiempo,
comenzando en el tiempo t = 0 y terminando en el
tiempo t = t0 .
Esta se puede expresar de la siguiente forma: H/t0
76 Introducción al Control de Procesos para Ingenieros
0
para t 0, t t0
ut H
t para 0 t t0
0
O expresada en transformada de Laplace:
U s
H
1 e st0
t s 0
Entrada en impulso
Cuando el tiempo t0 tiende a cero entonces se trata de un caso especial de la función pulso
denominada función impulso:
Esta se puede expresar de la siguiente forma:
0 para t 0, t
ut
t0 H
lim para 0 t t0
t00
t0
O expresada en transformada de Laplace:
U s H
En este caso como la altura de la función impulso es igual a H/t0 y la duración es t0 entonces el área bajo la
curva será igual a H. Cuando la el área de la función impulso es igual a la unidad entonces esta función
recibe el nombre de impulso unitario o delta de Dirac. Esta se designa comúnmente como:
t t0
Y cumple las siguientes condiciones:
t t0 0 para t t0
t t0 para t t0
t t
0 dt 1
Se debe notar que un impulso de magnitud infinita y duración cero es un término netamente matemático
que no existe en la realidad. Sin embargo, cuando la magnitud de un pulso es muy grande y su
duración muy pequeña, entonces el pulso de entrada se puede aproximar a una función impulso que es más
sencilla de tratar matemáticamente.
La función impulso unitario es muy utilizada para el estudio de respuesta de sistemas ya que permite entre
otras cosas estudiar la respuesta de sistemas en presencia de una entrada discontinua, como el estudio de
impactos por ejemplo.
0 para t 0
H
para 0 t t0
t
ut 0
H para t0 t 2t0
t0
⁝
Esta función suele ser muy utilizada para definir la respuesta de sistemas en donde la entrada varía
muy rápidamente y oscilar entre valores máximos y mínimos. Por ejemplo la respuesta de sistemas
electrónicos y de audio. El objetivo de esta función de entrada es muy similar al de la función
sinusoidal pero permite observar más fácilmente los resultados obtenidos siendo su cómputo más
complejo por tratarse de una función discontinua.
Respuesta completa
La respuesta completa del sistema será entonces la suma de la respuesta transitoria más la respuesta
en estado estable:
y yT yE
y Cet KH
Falta entonces obtener el valor de C con las condiciones iniciales. Para un sistema de primer orden se
requiere de una sola condición inicial, la cual depende del estado en que se encuentra el sistema en el
momento inicial en que se aplica la entrada (t = 0). En este caso tomaremos por simplicidad y sin pérdida de
generalidad las condiciones iniciales siguientes:
t0 y0 0
Para determinar el valor de la constante C se debe sustituir el valor de las condiciones iniciales y del tiempo
en el cual estas se aplican en la ecuación de la respuesta del sistema:
0 Ce0 KH
C KH
Finalmente la respuesta del sistema será:
y KHet KH
Con esta ecuación que define la respuesta del sistema se puede determinar el valor de la respuesta y para
cualquier instante de tiempo, es así como por ejemplo, y en el caso que la ganancia K = 1:
Para t
y H 1 e 0.6321H
Para t
4 y H 1 e4 0.9817H
Esta respuesta se puede representar gráficamente como sigue:
u
98.17%
63.21% y
H
0 r 4r t
A partir de estos resultados se define la constante de tiempo (r) como el tiempo requerido para que la
respuesta de un sistema de de primer orden alcance el 63.21% del valor del escalón al cual fue sometido.
Otra definición de la constante de tiempo es el tiempo necesario para que un sistema de primer orden
sometido a una entrada en escalón alcance su valor en estado estable si mantuviera la velocidad de
respuesta inicial.
Esta constante de tiempo es entonces un parámetro que define la velocidad de respuesta de los
sistemas de primer orden, mientras más grande es la constante de tiempo más lenta será la respuesta del
sistema. Se considera que un sistema de primer orden se encuentra en estado estable cuando el tiempo
transcurrido a partir de la aplicación de la entrada es superior a 4r.
Respuesta completa
La respuesta completa del sistema será entonces la suma de la respuesta transitoria más la respuesta
en estado estable:
y yT yE
y Cet KH KHt
Falta entonces obtener el valor de C con las condiciones iniciales. En este caso tomaremos por simplicidad y
sin pérdida de generalidad las condiciones iniciales siguientes:
t0 y0 0
Para determinar el valor de la constante C se debe sustituir el valor de las condiciones iniciales y del tiempo
en el cual estas se aplican en la ecuación de la respuesta del sistema:
0 Ce0 KH KH0
C KH
Finalmente la respuesta del sistema será:
y KHet KH KHt
Esta respuesta se puede representar gráficamente, en el caso de K = 1, como sigue:
u = Ht
a
y
b
0 t
—rH
En este caso la respuesta en estado estable tiende a una recta paralela a la entrada pero desfasada en un
valor —rH.
yE H H
sin t cos t
1 1
2 2
2 2
Respuesta completa
La respuesta completa del sistema será entonces la suma de la respuesta transitoria más la respuesta
en estado estable:
y yT yE
H H
y Cet sin t cos t
1 1
2 2
2 2
Falta entonces obtener el valor de C con las condiciones iniciales. En este caso tomaremos por simplicidad y
sin perdida de generalidad las condiciones iniciales siguientes:
t 0 y0 0
Para determinar el valor de la constante C se debe sustituir el valor de las condiciones iniciales y del tiempo
en el cual estas se aplican en la ecuación de la respuesta del sistema:
0 Ce 0 H H
sin 0 cos 0
1 1
2 2
2 2
H
C 1 2 2
Finalmente la respuesta del sistema será:
H t H H
y e sin t cos t
1 1 1
2 2 2 2 2
2
Para tener una forma de respuesta más fácil de representar gráficamente utilizaremos la relación
trigonométrica:
M sin t N cos t M 2 N 2 1 2 sin t M
tan 1 N
Obtenemos entonces que la solución en estado estable se puede expresar como:
H H
y et sint
1
2 2
1 22
Donde:
2H2
1
tan 1
tan1
H2 2
1
H
H y 1 2 2
0 t
1/ω
La respuesta en estado estable tiende a una función sinusoidal de menor amplitud y desfasada un ángulo .
La respuesta para un K G 1 puede obtenerse con la mismo procedimiento expresado para los sistemas
de primer orden.
Respuesta de un sistema de segundo orden a una entrada en escalón
En este caso la entrada para t 0 será: u H , luego la ecuación diferencial del sistema se puede escribir
como:
D 2 y 2nDy n2 y 2 H
La respuesta de este sistema corresponde entonces a la solución de la ecuación diferencial, la cual se
obtiene en dos partes:
2 4 2 n 4 n
2 2
D1,2
n
2
D1,2 2 1
n n
Luego la solución homogénea depende de la forma de las raíces de la ecuación característica, en este caso
la forma de las raíces depende del radio de amortiguamiento (£), la solución homogénea puede entonces
tomar varias formas como se presenta a continuación:
Si 1 aparecen dos raíces reales positivas.
La solución transitoria será entonces de la forma
n 2 1 n 2 1 t
t
yT C1e
n
n C2 e
Si 1 aparecen dos raíces reales iguales.
La solución transitoria será entonces de la forma
yT C 1e t C 2te t
n n
Donde C1 y C2 son constantes que dependen de las condiciones iniciales del sistema.
Respuesta completa
La respuesta completa del sistema será entonces la suma de la respuesta transitoria más la respuesta
en estado estable:
y yT yE
Donde yT es cualquiera de las cuatro opciones determinadas anteriormente.
Si 1 :
2 1 n 2 1 t
n
t
y Ce
n
H
Si 1 n C2 e
T 1
:
t t
y C e n C te n H
Si 0 1
T 1 2
:
y e n C sin 1 2 t
T 1
t
n
cos 1 2 t H
Si 0 : C2 n
Falta entonces obtener el valor de las constantes C1 y C2 con las condiciones iniciales. Para un sistema
de segundo orden se requiere de dos condiciones iniciales, la cuales depende del estado en que se
encuentra el sistema en el momento inicial en que se aplica la entrada (t = 0). En este caso tomaremos por
simplicidad y sin perdida de generalidad las condiciones iniciales siguientes:
t 0 y0 0; Dy0 0
Para determinar el valor de las constantes C1 y C2 se debe sustituir el valor de las condiciones iniciales y del
tiempo en la ecuación de la respuesta del sistema y en su derivada:
Si 1 :
2 1 n 2 1 0
n
0 C1e
n
0
n H
C2e
1
n 2 1 0 n n 2 1 0
n 2
2 1
0 C1 n n e
C2 n n e
C1 C2 H
C1 n n 1 C
2
2 n
2
n 1 0
C H H 21
n
n
1
2 2 1
C H
n n n 1
2
2
2 2 1
n
Si 1 :
0
0 C1 e n C 2
0e 0 H n
C 0e n 0
n 1 2 n 2
0 C1 H
C1 H
0 C C C H
n 1 2 2 n
Si 0 1 :
0 e n C sin 1
0
2
0 C 1
2
0 H
1 2
cos
n
0 n
e n 0
sin 1
2
0 2
cos 1 2
0
C 1
n C n
e n C1 n 1 2 cos 1 2 0
0
1 sin 1 0
2 2
n C2 n n
0 C2 H
C2 H H
0 n C 2 C1
n
12 C1
1
2
n
n
Si 0
:
0 C1 sin n 0 C2 cos n 0 H
0 C1n cos n 0 C2n sin n 0
0 C2 H C2 H
0C C 0
1 n 1
Si 1 : n 1
yT He n n Hte n H
t t
Si 0 1 :
H
e n H
n 1 2 t cos n 1 2 t
nt
yT
sin n 1 H
2
Si 0
:
yT H cos nt H
Con esta ecuación que define la respuesta del sistema se puede determinar el valor de la respuesta y para
cualquier instante de tiempo. Esta respuesta se puede representar gráficamente como sigue:
y0 1
y 0
H y 1
y 1
0 t
En función de estos resultados las respuestas de los sistemas de segundo orden se clasifican según el valor
del coeficiente de amortiguamiento ( ):
Si 1: Respuesta Sobre-Amortiguada. El sistema responde a una entrada en escalón lentamente.
Si 1: Respuesta Críticamente Amortiguada. El sistema responde a una entrada en escalón con
máxima velocidad de respuesta sin oscilar.
Si 0 1 : Respuesta Sub-Amortiguada. El sistema responde rápidamente pero comienza a oscilar
tratando de alcanzar el valor del escalón.
Si 0 : Respuesta sin amortiguamiento. El sistema oscila constantemente entre 0 y 2H.
2B 2 BH
Hn n 2H
2 B 2 A A
n
0 n
Respuesta completa
La respuesta completa del sistema será entonces la suma de la respuesta transitoria más la respuesta
en estado estable:
y yT yE
yE yT 2H Ht
n
Donde ( yT ) es cualquiera de las cuatro posibilidades vistas en la entrada en escalón. Falta
entonces obtener el valor de C1 y C2 con las condiciones iniciales, para obtener la solución completa.
y 1
0 t
y0 1
2H
n
En este caso la respuesta en estado estable tiende a una recta paralela a la entrada pero desfasada en un
2H
valor .
n
Se sustituye la solución supuesta y su derivada en la ecuación original, para obtener el valor final de la
solución particular:
A 2
sin t B 2 cost 2
n A cost B sin t n
2
Asin t B cost 2 H sin t
Se obtienen dos ecuaciones, una con los términos con sinωt y otra con los términos con cosωt:
A 2 2n B n 2 A 2 H
n
B 2 2 A n 2 B
n
0 2
2
2
2 2
2 A n
H
2nn
2 2
2 2
n A 2nB H n
n n
2 A n 2 2 B 0 2 3H
B n
22 2
2
n 2
n
2 2
n
2H
n 2
2 sin t
cos t
2
n
2n 2
2 2
yEas
2
n
n
Respuesta completa
La respuesta completa del sistema será entonces la suma de la respuesta transitoria más la respuesta
en estado estable:
y yT yE
Donde ( yT ) es cualquiera de las cuatro posibilidades vistas en la entrada en escalón. Falta
entonces obtener el valor de C1 y C2 con las condiciones iniciales, para obtener la solución completa.
y 1
u
H y <H
t
y0 1
0 1/ω
La respuesta en estado estable tiende a una función sinusoidal de menor amplitud y desfasada un ángulo .
Solución de las ecuaciones en representación de estado.
En el caso de los sistemas representados en espacio de estado para obtener la respuesta del sistema
es necesario primero resolver la ecuación de estado (ecuación dinámica):
x˙(t) = A(t)x(t) + B(t)u(t)
Una vez que se tenga esta solución la respuesta del sistema corresponde a una operación algebraica
establecida en la ecuación de salida:
y = Cx
Como se trata simplemente de una representación diferente de un mismo sistema tendremos también una
respuesta transitoria (también llamada respuesta libre) y una respuesta en estado estable (también
llamada respuesta forzada).
Respuesta transitoria
En el caso general en donde el modelo es variable en el tiempo, se requiere definir una matriz de transición de
estado Φ que interviene en la solución de la ecuación de estado.
Donde A(t) es continua respecto de t, está dada por la respuesta libre del sistema, que es la respuesta
del sistema producto de sus condiciones iniciales:
x(t) = Φ(t, t0)x0 ∀t ≥ t0
En el caso general, de un sistema variable en el tiempo, es raramente posible calcular la solución exacta de
la matriz de transición, por lo cual su solución suele calcularse de forma numérica. Mientras que para
sistemas LTI, la matriz de transición si se puede calcular de forma explícita, solución que será
mostrada más adelante.
Respuesta completa
La respuesta completa del sistema representado por la solución de la ecuación dinámica de estado:
x˙(t) = A(t)x(t) + B(t)u(t) x(t0) = x0
Se puede obtener entonces, aplicando el principio de superposición, a partir de matriz de transición, este
será la suma de la respuesta transitoria más la respuesta en estado estable, de la cual resulta la expresión:
t
x(t) = Φ(t, t0)x0 + ƒ Φ(t, r)B(r)u(r)dr
t0
La salida del sistema se puede entonces calcular con la expresión:
t
y(t) = C(t)Φ(t, t0)x0 + C(t) ƒ Φ(t, r )B(r )u(r )dr + D(t)u(t)
t0
Caso particular de los sistemas Lineales Invariantes en el Tiempo (LTI). Solución en el dominio del
tiempo.
Consideremos el modelo LTI:
x˙(t) = Ax(t) + Bu(t)
y(t) = Cx(t) + Du(t)
En el caso de un modelo estrictamente propio, para el cual la matriz de transición directa D = 0 , se tiene
el modelo LTI homogéneo (sin entrada) y con las condiciones iniciales x(t0 ) = x0 ,
x˙(t) = Ax(t)
y(t) = Cx(t)
Nota: Se debe tomar en cuenta que a diferencia de los números escalares (a. b = b. a) las operaciones con
matrices no son conmutativas (A. B G B. A), por lo tanto:
e(a+b)t G eat+ebt
Teorema de Cayley-Hamilton.
Este teorema permite simplificar el cálculo de la aproximación del exponencial por la serie de Taylor,
limitándolo a un número finito de términos.
Toda matriz cuadrada A es solución de su ecuación característica. Por lo tanto, si QA(s) es el
polinomio característico de A ∶ QA(s) = det(sI — A)., entonces para una matriz cuadrada A de
dimensión n, QA(s) es un polinomio de grado n.
QA(s) = det(sI — A) = sn + an–1sn–1 + . . . + a1s + a0
Es decir, An se puede expresar como una combinación lineal de las potencia inferiores de A :
An = —an–1An–1 — . . . — a1A – a0I
Caso particular de los sistemas Lineales Invariantes en el Tiempo (LTI). Solución en el dominio de
Laplace.
También se puede determinar la matriz de transición utilizando la transformada de Laplace del modelo del
sistema, para ello partimos del modelo LTI:
x˙(t) = Ax(t) + Bu(t)
y(t) = Cx(t) + Du(t)
Considerando las condiciones iniciales iguales a cero, condición que se aplica a todas las transformadas de
Laplace, le solución de la ecuación de estado es de la forma:
t
x(t) = Φ(t)x(0) + ƒ Φ(t — r)Bu(r)dr
0
Y se puede demostrar que:
Φ(s) = f{Ф(t)} = [sI — A]–1
Ejemplo: se quiere determinar la matriz de transición del sistema LTI representado por la matriz:
0 —2
A=[ ]
1 —3
La matriz de transición en transformada de Laplace se determina con la expresión:
1 0 0 —2 –1 s 2 –1
Φ(s) = [sI — A]–1 = [s ]— ]] = [ ]
[ [ 1 —3 —1 s + 3
1 s + 3 —2 0 1 1
Φ(s) = [ ]= s + 3 —2
[
Δs 1 s (s2 + 3s + 2) 1 s ]
(s + 3) 2
F(s + 1)(s + 2) —(s + 1)(s + 2)1
Φ(s) = I I
I 1 s I
L(s + 1)(s + 2) (s + 1)(s + 2) ے
Si consideramos unas condiciones iniciales diferentes de conocidas, por ejemplo x1(0) = x2(0) = 1, se
pueden determinar entonces los valores de los estados para cualquier tiempo t:
x1(t) 1 e
[ ] = Φ(t) [ ] = [ –2t]
x2(t) 1 e–2t
Ejercicios F
1. Cuál es el valor o el rango de valores de C para que el sistema de
M
la figura responda en forma sobreamortiguada, críticamente
amortiguada, subamortiguada y sin amortiguamiento. Hallar las
expresiones para ωn y £. Para el mismo sistema trate de hacer y
una gráfica aproximada de la respuesta (y) para cada uno de los
cuatro casos. Determine adicionalmente el valor de la frecuencia
k C
natural amortiguada d n 1 2 .
5. Se tiene una resistencia eléctrica introducida en un baño, el cual se puede considerar aislado
térmicamente del exterior, y que se encuentra inicialmente a temperatura ambiente. La resistencia
se prende, y a los 5 minutos se introduce un termómetro de primer orden. ¿Cuál es la lectura del
termómetro a los 2 minutos de haberse introducido? Si: Ta = 20 ºC, I = 2A, R = 2000 fi,
masa = 10Kg, Cp = 8 KJ/KgºC (KW. s/KgºC), constante de tiempo del termómetro de r =
30 s. Haga una gráfica aproximada de la respuesta del proceso.
7. Se tiene un envase con mezcla de agua y hielo sobre un mechero para calentarla, y fuera del envase
se tiene un termómetro de primer orden con una constante de tiempo de 0.5 minutos. Se prende
el mechero, con lo cual el agua comienza a calentarse a razón de 1 ºC por minuto. Treinta
segundos después se introduce en el envase el termómetro. Un minuto después de introducido el
termómetro se reduce la llama del mechero, con lo cual la temperatura del agua se mantiene
constante. Cuál es el valor indicado por el termómetro al cabo de 2 minutos de haberse prendido el
mechero.
10. Se tiene un voltímetro de corriente continua, el cual debido a su estructura interna se comporta
como un instrumento de segundo orden con relación de amortiguamiento igual a £ = 0.9 y
frecuencia natural igual a ωn = 1 rad/s. El instrumento se encuentra inicialmente desconectado y
se conecta a una batería de V = 12 v. Al cabo de un minuto se le conecta a la batería un cargador
por lo cual su voltaje comienza a aumentar a razón de 0.1 voltios por segundo hasta llegar a V = 14
v, a partir de este momento el voltaje de la batería permanece constante. ¿Cuál será el voltaje
indicado por el instrumento al cabo de 2 minutos de haberse conectado?
11. Se tiene un medidor de flujo para gasolina, con el cual se quiere medir la cantidad de combustible que
gasta un vehículo en un momento determinado. Inicialmente el vehículo se encuentra apagado por lo
que el flujo de combustible es nulo, y se enciende el motor dejándolo en velocidad de marcha mínima,
con lo cual empieza a consumir un caudal de combustible de 5 litros por hora. Al cabo de 5 minutos
se comienza a acelerar el motor progresivamente por lo que el consumo de combustible aumenta a
razón de 1 litro por hora cada minuto hasta alcanzar la velocidad máxima del motor, en donde éste
consume 20 litros por hora. Tres minutos después de haber alcanzado la velocidad máxima se apaga
el motor. Si el instrumento es de segundo orden con radio de amortiguamiento de 1.2 y frecuencia
natural de 60 radianes por segundo ¿Cuál será el valor de flujo indicado por el instrumento 24
minutos después de de haberse encendido el motor?
12. En la figura se observa un tanque de agua en donde se quiere medir el nivel mediante un manómetro
conectado en el fondo a través de la válvula V1, dicho tanque tiene un área A = 10 m2 posee una
entrada de agua Qe a través de la válvula V3 y una salida Qs a travez de la válvula V2. El
manómetro
1
n= rad
funciona como un instrumento de segundo orden con £ = 0.8 y 60
s
, e inicialmente la válvula V1
se encuentra cerrada, indicando el manómetro presión atmosférica; el nivel de agua en el tanque se
encuentra en 2 metros y los flujos de entrada y salida son iguales (Qe = Qs = 5m3/min).
Repentinamente se abre la válvula V1 y por tanto el
manómetro empieza a medir la presión en el fondo del V3
tanque que es proporcional al nivel, al cabo de 2 Qe
minutos se abre un poco más la válvula V2 por lo que 2m
el flujo de salida aumenta repentinamente al doble del
flujo de entrada (Qs = 2Qe), dos minutos después se h Ch =
vuelve a cerrar un poco la válvula V2 hasta que los A Qs
flujos de entrada y salida vuelven a ser iguales (Qs =
Qe), esto se mantiene durante 1 minuto al cabo del V1 V2
cual se abre un poco la válvula V3 por lo que el flujo de
entrada pasa a ser el doble del flujo de salida (Qe = 2Qs). a) Haga la gráfica aproximada del
proceso, entrada (nivel del tanque “h”vs tiempo), salida (nivel indicado por el manómetro
“hm” vs Tiempo). b) Cual será el valor del nivel en el cual se encuentra el tanque (h) al cabo de
7 minutos de comenzar el proceso. c) Cual es el valor de presión (Pm), en psi., indicado por el
manómetro el cabo de 7 minutos de comenzar el proceso. d) Cual es el error que se comete en
la lectura del nivel a través del manómetro.
Datos adicionales: agua = 9.800 N/m3 . 1 Pa (N/m2 ) = 1.45 x10 — 4 psi (lbf/plg2 ). P =
h
Tema 6. Estabilidad de Sistemas
Introducción
La noción de estabilidad es fundamental en el desarrollo de sistemas de control y en particular para
los sistemas retroalimentados. La ausencia de esta propiedad vuelve inútil en la práctica a cualquier
sistema.
Existen diversas formas de definir la estabilidad. Por ejemplo se puede hablar de la noción de estabilidad de
un sistema autónomo que no es idéntica a la utilizada en sistemas sometidos a entradas y salidas (en
donde la energía puede tener ciertos límites). También podemos referir que la estabilidad entre la entrada
y la salida no necesariamente implica una estabilidad interna a los sistemas. Se puede hablar entonces
de estabilidad local, global o semiglobal según el sistema no lineal en consideración.
Aquí nos vamos a interesar solamente a un tipo de estabilidad, utilizada para una clase de sistemas
específicos, entre los que se encuentran los sistemas LTI, la cual es la descripción de estabilidad
entrada salida, la cual nos lleva a la definición de estabilidad de entradas y salidas limitadas,
comúnmente denominada estabilidad BIBO (Bounded Input Bounded Output).
VC VC VC
t t t
Estable Estabilidad Inestable
Tal consideración sugiere que los coeficientes de t en los términos exponenciales de la respuesta transitoria
sean números reales negativos o números complejos con parte real negativa.
Esto implica que para que un sistema sea estable las raíces de la ecuación característica deben ser
negativas o con parte real negativa. Esto es ya que la ecuación característica representa la parte
transitoria (homogénea) de la ecuación que rige el sistema.
De lo anterior podemos observar que la estabilidad no depende de la entrada sino que es una característica
propia del sistema.
yT Cet
Este será estable siempre que la constante de tiempo sea positiva.
y e n C sin 1 2 t C cos 1 2 t
T 1
t
n 2
n
yT en C1 sin n 1 2 t C
t
2 n
12 t
cos
Solo en este último caso el sistema será inestable, ya que el coeficiente de t en el
exponencial es positivo, luego los transitorios aumentan en el tiempo.
Estabilidad limitada
Es el caso que sirve de frontera entre l estabilidad absoluta y la inestabilidad, y se presenta cuando
las raíces de la ecuación característica tienen partes reales iguales a cero.
La respuesta resulta ser una oscilación permanente cuya amplitud ni crece ni decae en el
tiempo. Esto es el caso de 0 para un sistema de segundo orden por ejemplo.
Conclusión
De todo lo anterior podemos obtener las siguientes conclusiones:
1. Un sistema es estable si todas la raíces de la ecuación característica son negativas o con parte
real negativa.
2. Un sistema es inestable, si tiene en su ecuación característica alguna raíz positiva o con parte real
positiva.
3. Un sistema tiene estabilidad limitada si alguna de sus raíces son pares de imaginarios puros. Una
raíz cero no influye sobre la estabilidad porque la respuesta no es oscilatoria.
4. Si el polinomio de la ecuación característica tiene algún coeficiente ausente o negativo entonces el
sistema es inestable.
Jean-François DULHOSTE
4. s
4 s2 16 0
2
Estabilidad limitada ya que tiene dos pares de raíces imaginarias
puras
5. s2 16 s 1 0 Estabilidad limitada. Ya que el exponencial domina sobre la otra
función y no hay ninguna raíz positiva, ni con parte real positiva. Los exponenciales siempre dominan
sobre los cosenos y senos.
6. s
2s 8 s 2 s2 1 0 Que es lo mismo que: s 4s 2s 2 s2 1 0 Es inestable
2
Para hallar la estabilidad de cualquier sistema, se pueden utilizar métodos para facilitar la operación, Ya
que en algunos casos si el polinomio es de orden elevado será poco práctico encontrar las raíces. Entre
los métodos más sencillos para determinar la estabilidad están:
Criterio de Routh
Criterio de Hurwitz
Criterio de la fracción continuada
Estos métodos sirven para hallar la estabilidad de ecuaciones características en forma de polinomios.
Existirán además métodos más complejos para determinar la estabilidad a otros sistemas y obtener
mayor información.
Donde:
an , an1,, a0 : son los coeficientes de l ecuación característica.
a a anan3 a a anan5 a a anan7
b1 n1 n2 ; 2 n1 n4 ; 1 n1 n6
b b
an1 an1 an1
ba a b ba a b
c1 1 n3 n1 2 ; c 1 1 n5 n1 3 ; etc.
b1 b1
La tabla se continúa horizontal y verticalmente hasta que solo se obtengan ceros.
Tema 6. Estabilidad de Sistemas 101
Ejemplo 3.
Para el sistema: La tabla de Routh será:
s 6s 12s 8 0
3 2
1 12 0
s3
6 8 0
s2 32 3 0 0
s
8 0 0
s0
Como en la primera columna no hay cambios de signo entonces el sistema tiene 3 raíces negativas o
con parte real negativa y por lo tanto es estable.
Ejemplo 4.
Si se tiene la siguiente primera columna de Routh
determinar todas las conclusiones posibles:
s5 2
El sistema tiene dos raíces positivas o con parte
s4 3
real positiva porque hay dos cambios de signo.
s3 4 El sistema tiene una raíz cero por el cero terminal.
s2 5 Hay dos raíces negativas o con parte real negativa.
El sistema es de quinto orden.
s 2
Por lo tanto es inestable.
s0 0
Ejemplo 5.
Determinar todas las conclusiones posibles de la
siguiente primera columna de Routh:
El sistema es de orden 4, posee 4 raíces.
No tiene cambios de signo, luego no tiene raíces
positivas.
Escuela de Ingeniería Mecánica - ULA
102 Introducción al Control de Procesos para Ingenieros
Posee u par de raíces imaginarias puras.
Jean-François DULHOSTE
s4 3 Posee un par de raíces negativas o con parte real
s3 2 negativa.
El sistema tiene estabilidad limitada por el cero no
s2 5 terminal.
s0
s0 2
Ejemplo 6.
Haga un estudio de la estabilidad del siguiente El primer paso es entonces el construir la tabla
sistema usando el criterio de Routh de Routh:
D5x 3D4x 7D3x 20D2x 6Dx 15x u
s5 1 7 6 0
3 20 15 0
s4 1 3 1 0
s3 11 15 0
6 11 0
s2
15 0
s
s0
La primera columna de la tabla de Routh no posee cambios de signo luego todas sus raíces son negativas o
con parte real negativa, por lo tanto el sistema es estable.
Ejemplo 7.
Determinar todas las conclusiones posibles de la
siguiente primera columna de Routh:
s10 3
s9 2 El sistema es de orden 10
s8 5 Tiene una raíz cero
Tiene dos raíces imaginarias puras
s7 4
Tiene un par de raíces positivas o con parte real
s9 3 positiva
s5 2 Es inestable por poseer un cambio De signo
s4 0
s3 1
s2 0
s 2
s0 0
Ejemplo 8. 0
Determinar, con el criterio de Routh, la estabilidad
del sistema cuya ecuación característica es:
s4 s3 s 1
La tabla de Routh será:
Tema 6. Estabilidad de Sistemas 103
s4 1 0 1
s3 1 1 0
s2
1 1 0
s 0
s0 1 0
d lim 1
1
0
Por lo tanto el sistema es inestable.
Sin embargo el cálculo de la tabla no era necesario ya que la ecuación característica posee dos coeficientes
negativos y hay uno ausente, se podría concluir directamente que el sistema es inestable.
Ejemplo 9.
Determinar la estabilidad del siguiente sistema: Le tabla de Routh es:
s5 + 2s4 + 3s3 + 6s2 + 5s + 3 = 0 s5 1 3 5
s4 2 6 3
s3 0 → ε 7/2 0 Le sistema es inestable.
—∞ 3 0
s2 | 7/2 0
s
s 0 3 0
Jean-François DULHOSTE
an1 an3 an5 an7
an1 an3 an1 an3 an5
an an2 an4 an6
1 an1 ; 2 ; 3 an an2 an4 ; 4
an an2 0 a a an1 an3 an5
n1 n3
0 an2 an4
an
Y así sucesivamente hasta llegar a Δn donde n es el orden de la ecuación característica del sistema.
El criterio de Hurwitz dice que dice que todas las raíces de la ecuación característica serán negativas o con
parte real negativa si los determinantes son positivos ( i 0 para i 1,2,n ).
Por supuesto al igual que para el criterio de Routh si existe algún coeficiente ausente o negativo en la
ecuación característica se puede concluir que el sistema es inestable y no es necesario aplicar el método.
Ejemplo 10. Determinar la estabilidad del sistema siguiente usando el criterio de Hurwitz:
a s 4 a s3 a s 2 a s a 0
4 3 2 1 0
Obtenemos entonces los determinantes de Hurwitz, que en este caso serán 4 por ser el sistema de
4to orden:
1 a3
a3 a1
a a a a
2 3 2 4 1
a 4 a2
a3 a1 0
3 a4 a2 a0 a a a a a 2
a a 2
3 2 1 0 3 4 1
0 a3 a1
a3 a1 0 0
a4 a2 a0 0
4
a3 a1 0
a a a a a a a
2
a4 a1
2
0 a3 a0 a4 a2 0 0 3 2 1 0 3
0 0 a3 a1
0 a0
a1
a4 a2
Por lo tanto el sistema será estable si:
a3 0
a3 a2 a4 a1 0
a3 a2 a1 a 0a 32 a 4a12 0
a
0
a a a a a
3 2 1
2
0 3
a a2 0
4 1
Ejemplo 11. Determinar si la siguiente ecuación característica representa un sistema estable o inestable:
s3 8s2 14s 24 0
Obtenemos entonces los determinantes de Hurwitz, que en este caso serán 3 por ser el sistema de
3er orden:
1 8
8 24
88
2
1 14
8 24 0
3 1 14 0 2112
0 8 24
Como todos los determinantes son positivos entonces el sistema es estable.
8 6 5 0 5
4 0 5 3
535
3 0
0 8 6 5
El sistema es inestable por tener determinantes de Hurwitz negativos.
Por ejemplo:
Si se tiene un sistema lineal cuya matriz A es:
0 1 0
A 0 0 1
6 11 6
La ecuación característica del sistema será:
1 0
sI A 0 1 3 62 11 6
6 11 6
1 2 3 0
Por lo tanto los valores propios de la matriz A serán:
1 1
2 2
3 3
Lo cual indica que este sistema es estable.
Ejercicios
1. Determinar si los siguientes sistemas son estables (utilice los dos criterios):
a. s3 4s2 8s 12 0
b. s3 2s2 s 2 0
c. s3 7s2 7s 46 0
d. s 4 3s3 6s2 9s 12 0
e. s4 8s3 24s2 32s 16 0
f. D6 x 2D5x 4D4x 6D3x 8D2x 4Dx 5x u
2. Citar razonablemente toda la información suministrada por las siguientes primeras columnas de
Routh:
2 2 1 1 2
3 3 2 2 3
4 4 3 3 4
e.
5 5 4 d. 4 5
a.
0 b. 0 c. 0 5 6
6 6 5 6 7
7 7 0 7
0 0 6
0 8
3. Hallar el valor o el rango de valores de K para que los siguientes sistemas sean estables (utilice los
dos criterios):
a. 3s3 3s2 s K 2 0
b. s 4 6s3 11s2 6s K 0
c. s 3 4 K s 2 6s 16 8K 0
d. 3s 4 2s3 Ks2 s 2 s 2 K 0
e. s5 s 4 s3 s2 Ks3 s K 0
4. Determine la estabilidad de los sistemas representados por las siguientes matrices de estado:
0 1 0 0
a. A 0 0 1, B 0 , C 1 0 0
1 8 1
2 4 0
0 1 0 B 0 , C 1 0 0
b. A 0 0 1 ,
2 1 2 2
0 2 0 0
c. A 0 0 1, B 0 , C 1 0 0
23 7 3
7
5. Obtenga la representación canónica de Jordan para los sistemas anteriores.
6. Determine el valor de K para que los siguientes sistemas sean estables
K2
0 0 0
d. A 0 K 1 0 , B 0 , C 1 0 0
0 K 1
0
0
0 1 0 , B 0 C 1 0 0
e. A 0 0 1 ,
1 K2
1 6
3
3
0 1 0 1
A 0 1 B 1
0
4 K 8 2 K 1
, C 1 0 0
6
,
110 Introducción al Control de Procesos para Ingenieros
Asint B sint
Sistema
La respuesta en estado estable a una entrada de este tipo será de la forma: Bsint
La representación gráfica de la entrada y respuesta en estado estable está dada en la siguiente figura.
Asint
B sint
2π
π
ωt
En la práctica raras veces los sistemas de control están sometidos a señales sinusoidales, pero la
información que se obtiene por el análisis sinusoidal se puede usar para establecer la naturaleza de la
respuesta a una gran variedad de señales. Además el análisis es conveniente para el manejo analítico
y experimental.
Gs 1
s 1
La respuesta en estado estable para una entrada sinusoidal de la forma Hsin(ωt) se puede hallar
sustituyendo a s por jω como sigue:
Gs 1
j 1
Si se multiplica y divide por el conjugado para no tener términos imaginarios en el denominador:
1 j 1 1 j
Gs
j 1 j 1 2 2 1
En este caso la relación de amplitud y el ángulo de fase serán:
2 2
2 1 2
1 1
2 2 2 2
2 2 12 2 2 12
M G j 1 1 1 2 2
Im 2 2 1
1 tan1
G j tan
1
tan
Re 2 12
1
La respuesta en estado estable (respuesta en frecuencia) será:
y t EE H
sint tan1t
1 2 2
Diagramas de BODE
Los diagramas de Bode son una forma de representar la respuesta en frecuencia de un sistema de
control en coordenadas rectangulares, este consiste de dos gráficos:
1. Diagrama de amplitud o atenuación: en donde se grafica la relación de amplitud M de G(jω) en
decibelios (db) contra el logaritmo de la frecuencia.
El valor M expresado en decibelios se obtiene como: M db 20 log10 M
2. Diagrama de fase: en donde se grafica el ángulo de fase de G(jω) contra el logaritmo de la
frecuencia.
Los diagramas de Bode tienen las siguientes características:
Como se grafica el logaritmo de ω, los factores producto y cociente de G(jω) se convierten en
sumas de rectas.
Los ángulos de fase también se suman y restan en forma natural.
Tema 7. Análisis de la respuesta en frecuencia 111
El diagrama de Bode en la mayoría de los casos se puede aproximar mediante segmentos de recta lo
cual simplifica la construcción.
Ganancia ( Gs K )
Amplitud:
M db 20 log K ctte
0
Fase: tan1 0º si K 0
K 180º si K 0
M db
20log K
Fase:
tan1j P 90º P
La amplitud en este caso es una recta con pendiente ±20P y pasa por Mdb = 0 para ω = 1. Es así como
la pendiente será de 20 db por década de frecuencia. Es decir cambia de 20 db cada vez que el valor
de la frecuencia es multiplicado por 10. El ángulo de fase es una recta horizontal que pasa por 90P.
0
0.1 1 10 100 1000 log
-20 P 1
P 2
-40
P 1
90
0º
0.1 1 10 100 1000
P 1 log
-90º
P 2
-180º
20 log1 1
Amplitud: M 20log
db j 1 j
Fase: tan1
Se obtienen en este caso curvas idénticas pero con signo invertido
Si se requiere el valor exacto del diagrama de Bode basta con corregir el valor de la amplitud y esto se hace
calculando el error en algunos puntos.
El error máximo se produce a la frecuencia de cruce y es aproximadamente igual a 3db ya que:
Emax 20 11 20 log 10 log 2 3.03db
log
M db
1s
20
Asíntotas
0 10100
0.1 1 1000
log
-20 1
Curva exacta
1 s
-40
1s
45º
0º
0.1 1 10 100 1000
log
-45º
1
- 1
s
2
1
Factores de cuadráticos ( Gs s 2 s ) n
2 2
n n
j 2n j n j n 2j n 1 n n
Se estudiará aquí el caso en que £ ≤ 1, ya que en los otros casos G(s) puede descomponerse en el producto
de polos simples.
Amplitud:
M 20 log
1 n
2 2
4 2 n
2
db
2 n
1 n 2
el cual depende del valor de £ y tiende a crecer cuando £ decrece.
Fase:
tan1
0 0º
n 90º
180º
La curva del ángulo de fase es antisimétrica respecto al punto de inflexión y también depende del
valor de £.
M db 0.05
4 0.1
0
0.1
1
1 10 100 1000
log n
-40
-80
0 0.1 1 10 100 1000
0.05 log n
0.1
-90º 1
-180º
s2 2 s 2
n
El diagrama de Bode para G s
n se puede obtener invirtiendo las curvas anteriores.
2
n
Estabilidad relativa
Los conceptos de margen de ganancia y margen permiten determinar la estabilidad de un sistema según el
siguiente criterio:
Un sistema es estable si MG 0 y m 0
Nótese aquí que estos conceptos no solo indican la estabilidad en términos absolutos sino que permiten dar
un margen de que tan lejos está un sistema de la estabilidad o inestabilidad, esto ya que entre más
pequeños se hagan los valores del margen de ganancia y de fase más tendera el sistema hacia la
inestabilidad y viceversa.
1. Se determina la frecuencia crítica (ωC) trazando para ello una horizontal que pase por -180º hasta
cortar el diagrama de fase luego se traza una vertical hasta el eje de las frecuencias.
2. Se traza una vertical que pase por ωC y que corte el diagrama de amplitud. La diferencia entre
el punto de corte con 0db da el margen de ganancia (MG).
3. Se determina la frecuencia de ganancia crítica (ωg ), trazando para ello una horizontal que pase por
0db hasta cortar el diagrama de amplitud luego se baja verticalmente hasta el eje de las
frecuencias.
4. Se traza una vertical que pasa por ωg y corte el diagrama de fase. La diferencia entre -180º y
el punto de corte da el margen de fase.
5. Si MG > 0 y mф > 0 se dice que el sistema es estable.
M db
4
g
0
-40
MG 0 log
-80
0º
log
C
m 0
-180º
6.
Gs 10 j 3
j j 2 j j 2
2
Por comodidad para el trazado de las graficas se rescribe la ecuación en la forma normalizada con los
términos independientes iguales a 1:
Gs
103 j 3
1
7.5 j 3 1
j 2 j 2 12
2
2 j 2 1
j j 2 1 j 2 j 2 1
2
Esta función se descompone en factores conocidos, que son en este caso:
1 1 1
7⏟.5
1 j
⏟
––
1 j 3
3
1 j 2
––
1 2 j 2
––––––
2 4 j 2 5
Gráfica asintótica
Seguidamente se hace la gráfica con las asíntotas identificando las frecuencias de cruce para los factores
3, 4 y 5:
1. Amplitud: M db 20 log 7.5 17.5 : recta horizontal que pasa por 17.5
Fase: 0
tan 1 0º : recta horizontal que pasa por 0º
7.5
2. Amplitud: M 20log j 20log : con pendiente -20P y pasa por M 0 para 1.
1
db
db
Fase:
tan1 j 1 90º : recta horizontal que pasa por 90.
3. Frecuencia de cruce 1 como 1 3
3 1
Amplitud: M db 20 log1 3 20 log
j 1 j 3
Fase:
tan1 3
0 0º
1 45º
90º
4. Frecuencia de cruce 2
1
Amplitud: M 20log 20log1 2 1 22
db
1 j
n
2 2
5. Frecuencia de cruce como 2 ,y 1 2
n 2
41 2 2 2
1 2 2 2 2 2
Amplitud:
M 20log
db
Si 2
se puede aproximar a M db 20 log1 línea recta horizontal
0
4
2 2
Si 2 se puede aproximar a M 20 40 log línea recta con
log db
pendiente -40db/década
2 2
Fase: tan1
1 2 2
0 0º
n 90º
180º
Grafica exacta
Para la gráfica exacta se debe calcular el valor exacto de la curva de magnitud y fase en algunos
puntos adicionales, los puntos mínimos necesarios son en general los correspondientes a las frecuencias de
cruce.
Magnitud
En 2
2 2
2 2 2 2 12 2 2
M
2
2 20log 12 20log 12 20log
20log 7.5 3 2 2 2
db
20log
M db 16.61
En 2
M db 7.06
En 3
3 2 3 2 3 2 2 3 2
32 20 log 12 20 log 12 20 log 1
M db 20 log 7.5 20 3 2 2 2
log
M db 5.75
1
Fase
1
En 2
1 0
2 3
2 2 1
2 2 190º
1
tan tan tan tan
tan 1 2
2 2
7.5 0 1 1
En 2
1 0 2
2 2
1 1 23 1 1 22
tan 122 2 143
tan
7.5 tan
0 tan1 tan1
En 3
0 3
1
1 1 33 1 32 1 32
tan 312 2 258.11
tan
7.5 tan
0 tan1 tan1
Diagrama de Bode
Análisis de estabilidad
MG 18
m 50º
Gs s 2 s 2
s 3 s 4s 1s 2 3s 1
s5 11s 4 44s3 77s 2 55s 12
Determine la estabilidad e las siguientes funciones mediante los diagramas de Bode.
Gs s 2 s2
s s 4s 1s 2 3s s5 11s 4 44s3 77s 2 55s
12
3 1
Diagrama de Nyquist
El diagrama polar comúnmente denominado diagrama de Nyquist de una función de transferencia sinusoidal
G(jω) es un diagrama de la magnitud o módulo de G(jω) en función
del ángulo de fase de G(jω) en coordenadas polares al variar el valor
de ω de menos infinito a infinito. Entonces, el diagrama polar es el
lugar de los vectores |G(jω)|∠G(jω) cuando ω varía de —∞ a ∞.
La figura presenta un ejemplo de este diagrama. Cada punto en el
diagrama polar de G(jω) representa el extremo terminal de un
vector para un valor determinado de ω.
Para construir el diagrama polar se calcula directamente la
magnitud |G(jω)| y el ángulo de fase ∠G(jω) para cada frecuencia
ω.
Los datos también se pueden obtener a partir del diagrama de Bode.
Si se requiere la multiplicación de dos funciones de transferencia
sinusoidales, se puede obtener multiplicando, para cada frecuencia,
las funciones de transferencia sinusoidales individuales por medio de
una multiplicación algebraica compleja. Es decir:
Si
G j G1 jG2 j
Entonces
G j G1 j G2 j
G j G1 j G2 j
En general, si se desea el diagrama polar de G(jω) = G1 (jω) ∙ G2 (jω), es conveniente trazar primero el
diagrama logarítmico de G(jω) y luego convertirlo en un diagrama polar, en lugar de dibujar los diagramas
polares de G1(jω) y G2(jω) y multiplicar ambos en el plano complejo para obtener el diagrama polar de
G(jω).
Una ventaja al utilizar un diagrama polar es que presenta las características de respuesta en frecuencia de
un sistema en todo el rango de frecuencias, en un solo diagrama, la desventaja es que el diagrama no indica
claramente las contribuciones de cada factor individual de la función de transferencia.
G j 1 1
tan1 T
1 jT
1 2T 2
0 : G j0 10º
Para
1T:
1 1 45º
2
G j
T
Si ω tiende a infinito, la
magnitud de G(jω) tiende a cero
y el ángulo de fase tiende a —
90º. El diagrama polar de esta
función de transferencia es un
semicírculo, con centro en 0.5
sobre el eje real
y radio igual a 0.5, cuando la frecuencia ω varía de cero a infinito, como se puede ver en la figura.
Para el caso sobreamortiguado, al incrementar £ más allá de la unidad, el lugar de G(jω) tiende a ser
un semicírculo, lo que resulta del hecho de que para un sistema fuertemente amortiguado las raíces
características son reales, y una de ellas es mucho más pequeña que la
otra. Como para un valor de £ suficientemente alto la raíz más grande de
la respuesta se hace muy pequeña, el sistema se comporta como uno de
primer orden.
n n n n
La porción de baja frecuencia de la curva es: limG j 10º
0
La porción de alta frecuencia es: lim G j 180º
Como la parte imaginaria de G(jω) es positiva y crece en forma monótona para ω > 0, y la parte real
de G(jω) decrece en forma monótona a partir de la unidad, la forma general del diagrama polar de
G(jω) es como aparece en la figura. El ángulo de fase está entre 0° y 180°.
Retardo de transporte.
El retardo de transporte, o tiempo muerto G j e j T , se puede escribir
como:
G j 1cosT j sinT
Como la magnitud de G(jω) es siempre la unidad, y el ángulo de fase varía
linealmente con ω el diagrama polar del retardo de transporte es un círculo
unitario, como se puede ver en la figura.
Tema 7. Análisis de la respuesta en frecuencia 125
G j
a 0 j a1 j j 1 jT1 1 j2T
n n1
Donde n > m o el grado del polinomio denominador es mayor que el del numerador, tendrán las siguientes
formas generales:
1. Para λ = 0 o sistemas tipo 0: el punto de inicio del diagrama polar (que corresponde a ω = 0) es
finito y está sobre el eje real positivo. La tangente al diagrama polar en ω = 0 es perpendicular al
eje real. El punto terminal, que corresponde a ω = λ, está en el origen y la curva es tangente a uno
de los ejes.
2. Para λ = 1 o sistemas de tipo 1: el término jω en el denominador contribuye con - 90° al ángulo de
fase total de G(jω) para 0 ≤ ω ≤ ∞. Para ω = 0, la magnitud de G(jω) es infinita, y el ángulo de
fase es igual a menos 90°. A frecuencias bajas, el diagrama polar es asintótico a una línea
paralela al eje imaginario negativo. En ω = λ, la magnitud es cero, la curva converge hacia el
origen y es tangente a uno de los ejes.
3. Para λ = 2 o sistemas de tipo 2: el término (jω)2 en el denominador contribuye con -180° al ángulo
de fase total de G(jω) para 0 ≤ ω ≤ ∞. En ω = 0 la magnitud de G(jω) es infinita, y el ángulo de
fase es igual a menos 180°. A frecuencias bajas, el diagrama polar es asintótico a una línea paralela
al eje real negativo. En ω = λ, la magnitud se hace cero, y la curva es tangente a uno de los ejes.
En la figura se pueden ver las formas generales de las porciones de baja frecuencia de los diagramas polares
de los sistemas de tipo 0,
tipo 1 y tipo 2.
e
G j
jL 1
e
jL
1 j T 1 jT
La magnitud y el ángulo de fase son, respectivamente,
1 1
G j jL
e 1 jT 1 2T 2
y
G j e
jL 1
L tan 1 T
1 jT
Como la magnitud decrece en forma monótona a partir de la unidad y el ángulo de fase
también decrece monótona e indefinidamente, el diagrama polar de la función de transferencia
dada es una espiral, como aparece en la figura.
Nótese que en los ejemplos mostrados anteriormente se desarrollo el diagrama polar para una frecuencia ω
entre 0 e infinito. El diagrama polar para una frecuencia entre menos infinito y cero será la imagen
espejo respecto al eje real del diagrama para frecuencias de cero a infinito. Y el diagrama entre
infinito y menos infinito será siempre el origen de coordenadas.
Al examinar la estabilidad de los sistemas de control lineales utilizando el criterio de estabilidad de Nyquist,
se pueden presentar tres posibilidades:
1. No hay rodeo del punto — 1 + j0. Esto implica que el sistema es estable si no hay polos de
G(s)H(s) en el semiplano derecho del plano s; en caso contrario, el sistema es inestable.
2. Hay un rodeo en sentido anti horario o rodeos del punto —1 + j0. En este caso el sistema es
estable si la cantidad de rodeos antihorarios es la misma que la cantidad de polos de G(s)H(s)
en el semiplano derecho del plano s; en caso contrario, el sistema es inestable.
3. Hay un rodeo o rodeos del punto —1 + j0 en sentido horario. En este caso, el sistema es inestable.
Ejemplos
En los ejemplos siguientes, se supone que los valores de la ganancia K y de las constantes de tiempo (como
T, T1, y T2) son todos positivos.
Para valores pequeños de K, no hay rodeo al punto —1 + j0. Por lo tanto, el sistema es estable.
Para valores grandes de K, el lugar de G(s)H(s) rodea al punto —1 + j0 dos veces en
sentido horario, lo que indica dos polos en el semiplano derecho del plano s y el sistema es
inestable.
Para lograr buena exactitud, K debe ser grande. Sin embargo, desde el punto de vista de la
estabilidad un valor elevado de K brinda estabilidad pobre, incluso inestabilidad. Se debe hallar una
solución de compromiso entre exactitud y estabilidad.
T<T
1 2
T=T T1 > T 2
1 2
Estable Inestable
Esto significa que el sistema de lazo cerrado tiene dos polos de lazo cerrado en el semiplano
derecho del plano s y es inestable.
5. Determine la estabilidad de un sistema de lazo cerrado que tiene la siguiente función de transfe-
rencia de lazo abierto:
K s 3
GsH s
s s 1
Se pide determinar el rango de ganancia K que hace al sistema estable utilizando el criterio de
estabilidad de Nyquist. (La ganancia K es positiva).
Para examinar la estabilidad del sistema de control, hay que trazar el lugar de Nyquist de G(s),
donde
Gs G1sG2s
Sin embargo, en este punto no se conocen los polos de G(s). Por lo tanto, hay que examinar el lazo
menor por posible presencia de polos en el semiplano positivo del plano s. Esto es fácil de realizar
utilizando el criterio de estabilidad de Routh. Como:
G2 s 1
s s2 1
3
s3 s12 00
s s01 10
1 1 00
00
Como G(s) tiene dos polos en el semiplano derecho del plano s, se tiene que P1 = 2. Como Z1 = N1 +
P1, para la estabilidad se requiere que Z1 = 0 o N1 = —2.
Es decir, el diagrama de Nyquist de G(jω)/K debe rodear al punto — 1 + j0 dos veces en sentido
anti horario. De la figura se ve que, si el punto crítico queda entre 0 y — 0.5, entonces el diagrama
de G(jω)/K rodea al punto crítico dos veces en sentido antihorario. Por lo tanto, se requiere
que
—0.5K < —1
Entonces el rango de ganancia de K por condición de estabilidad es 2 < K.
Jean-François DULHOSTE
Tema 8. Controlabilidad y Observabilidad
Los métodos para la determinación de estas características de los sistemas varían según el tipo de
sistema. En este se estudia solo lo correspondiente a los Sistemas Lineales Invariantes en el tiempo (LTI).
Controlabilidad
La controlabilidad es una característica de un sistema, generalmente representado por un modelo en
espacio de estado, que nos indica si la evolución de una o varias de sus dinámicas (estados) pueden ser
modificadas por las entradas del sistema (control).
Definición
Un estado xi es controlable en t0 siempre que se pueda determinar una entrada u(t) que
conduzca todo estado inicial xi(t0) hacia 0 en un tiempo t0 ≤ t1 ≤ tf .
Si esta propiedad se cumple ∀ t y ∀ i = 1,· · · n entonces se dice que el sistema es completamente
controlable.
Notas:
Si un sistema no es completamente controlable entonces para algunas condiciones iniciales no
existe ninguna entrada capaz de llevar el sistema al origen.
La controlabilidad es una noción de suma importancia puesto que nos permite determinar si es
posible controlar un sistema para modificar su comportamiento (estabilización de un sistema
inestable, modificación de las dinámicas propias del sistema). Y es por ello que es fundamental en
la teoría de la síntesis de controladores en espacio de estado.
La matriz de controlabilidad es
entonces: 0 0 1 0 2 0
C = [1 0 1 0 1 0]
1 0 1 0 1 0
La matriz de controlabilidad tiene rango 2, puesto que las dos últimas líneas son idénticas, y por tanto solo
se pueden determinar menores principales diferentes de cero de dimensión 2. Por lo tanto esta
representación de estado No es Controlable.
Definición
Un estado xi es observable en t0 siempre que sea posible determinar xi(t0) conociendo a
y(t).
Nota: La noción de observabilidad es crucial para los sistemas donde es imposible medir la totalidad
del vector de estado, y en consecuencia se requiere de la estimación de este a partir de las variables de
salida.
Dualidad
Se dice que las propiedades de Observabilidad y Controlabilidad son dos nociones duales, puesto que
son propiedades intercambiables para sistemas duales adjuntos. Para explicar esto consideremos dos
sistemas S y S∗ definidos de la siguiente manera:
x˙ = Ax(t) + Bu(t) x˙∗ = ATx∗(t) + CTu∗(t)
y = Cx(t) y∗ = BTx∗(t)
S∗ se denomina sistema dual adjunto de S, puesto que está hecho con las mismas matrices pero
transpuestas e intercambiando la posición e C y B.
Teoría de la realización
Se denomina realización de una matriz de transferencia G(s) ∈ Œr×m , a toda representación de estado
(A, B, C, D) obtenida con la expresión:
G(s) = C(s1 — A)–1B + D
Cuando la realización (A, B, C, D) tiene un orden mínimo n (x ∈ Rn) esta se denomina realización mínima o
irreductible.
El ejemplo siguiente muestra como la obtención de una realización influye notoriamente en las propiedades
de observabilidad y controlabilidad del modelo en espacio de estado obtenido.
- No observable
—1y no0 controlable:0 - Observable y0 no—2
controlable: 1
A1 = [ ] B1 = [ ] C1 = [0 1] A2 = [ ] B2 = [ ] C2 = [0 1]
0 —2 1 1 —3 1
- Controlable y0 no observable
1 : 0 - Observable y controlable:
A3 = [ ] B3 = [ ] C3 = [1 1] A4 = [—2] B4 = [1] C4 = [1]
—2 —3 1
Se observa que la cuarta realización es mínima mientras que las otras tres no lo son.
El ejemplo muestra que por un lado las representaciones de estado no son equivalentes entre sí, y por otro
lado que la simplificación de un polo con un cero está estrechamente a las propiedades de observabilidad y
controlabilidad de la realización obtenida.
De esto se deduce que: Une realización de estado (A, B, C, D) de G(s) es mínima si y solo si esta es
observable y controlable.
Y(s) n
Y(s) c
i
G(s) = = = c0 + Σ
U(s) ∏n (s — (s — λi)
λi)i=1 i=1
Y en este caso se puede obtener una realización, donde la matriz A tenga una forma diagonal:
λ1 0 0 ⋯ 0 0 λ2 0
F I 0
1 I
A=I0 0 ⋱
⋮I 1
F11
I I
B = I⋮ I
I⋮ ⋱ ⋮I I⋮
L0 0 ⋯ ⋯ λn ے I
L1 ے
C = [c1 c2 ⋯ ⋯ c n] D = [c0]
En el caso particular de polos reales
conjugados, se puede realizar la α + iβ 0 1/2 – i/2 α β
[ =[ ][ ] [1 1 ]
siguiente transformación para obtener 0 α+ ] 1/2 i/2 —β α i —i
una matriz A de términos reales: iβ
nF n
ci c
Y(s) N(s)
G(s) = U(s) = (s — λ1)nF ∏n F (s — λi) = Σ (s — λ1)i + Σ (s —i λi)
i=n +1 i=1 i=n´+1
Se puede por lo tanto obtener una realización con la matriz A de la forma denominada cuasi-diagonal
de Jordan:
λ1 1 0 0 ⋯ 0 0
F0 ⋱ 0⋮ 1
II 0 0 λ1 0
0 F⋮1
1
II
I I
1
A= B=I I
I λnF+1 ⋮I I 1I
I⋮ ⋱ 0I I⋮
L0 0 ⋯ ⋯ λn ے I
L1 ے
C = [c1 c2 ⋯ ⋯ cn]
Para obtener esta representación multiplicamos y dividimos la función de transferencia por una variable
intermediaria Z(s) :
Y(s) (b1sn–1 + ⋯ + bn)Z(s)
G(s) = = n
U(s) (s + a1sn–1 + ⋯ + an)Z(s)
Estas operaciones no modifican las relaciones en el sistema, podemos escribir entonces las dos relaciones:
Y(s) = (b1sn–1 + ⋯ + bn)Z(s)
U(s) = (sn + a1sn–1 + ⋯ + an)Z(s)
La transformada inversa de Laplace de las dos relaciones es:
y = b1Dn–1z + ⋯ + bnz
u = Dnz + a1Dn–1z + ⋯ + anz
Para obtener esta realización se considera la ecuación diferencial original del sistema:
Dny + a1Dn–1y + ⋯ + an–1Dy + a n y = b1Dn–1u + ⋯ + bn–1Du + bnu
También se puede obtener una forma canónica a partir de una realización cualquiera realizando un cambio de
variable lineal.
Ejemplo 7: Obtener las realizaciones en las tres formas canónicas para el sistema representado por la
ecuación diferencial siguiente:
D3y + 6D2y + 11Dy + 6y = 6u
Ejercicios
Para las siguientes funciones de transferencia obtener las realizaciones, determinar la observabilidad y
controlabilidad de las mismas y presentar las realizaciones en las tres formas canónicas estudiadas.
–s–2
1. G(s) = s
2
s 3+10s 2+29s+20
s3+10s2+29s+20
2. G(s) = s4+8s 3+23s 2+28s+12
3. G(s) = s+1
s4–5s 2+4
4. G(s) = s+1
s 3+11s 2+36s+36
–s–2
5. G(s) = s
2
s4–5s 2+4
6. G(s) = s+1
s 4+8s 3+23s 2+28s+12
Tema 9. Acciones de Control 143
KM
et
Jean-François DULHOSTE
En la práctica un controlador discontinuo debe tener una zona muerta o histéresis (brecha diferencial). Esta
zona muerta se debe minimizar para mantener el error pequeño pero debe estar para evitar que se
produzcan oscilaciones.
Solenoide 115 V
Flotador
Válvula
Entra
agua
h Sale agua
La presencia de la brecha diferencial produce un error entre el valor deseado (nivel deseado en el ejemplo) y el
valor real de la variable (nivel real en el ejemplo). Pero esta es necesaria para evitar conexiones y
desconexiones muy cercanas en tiempo.
h
Brecha diferencial
t
Es por esto que a brecha diferencial se debe ajustar dependiendo de la exactitud deseada, de la frecuencia
de conexión y desconexión del elemento final de control (válvula solenoide) y de los valores a obtener.
Pivote
Flotador
Válvula
Entra
agua
h
Sale agua
Es
D
dt D
En otras palabras la salida del controlador varía en proporción a la velocidad de cambio del error, si el error
es constante no se produce ninguna acción de control.
Pivote
Flotador
Válvula
Entra
agua
h
Sale agua
m t K M s KI
t
e t ó F s
E s s
dt
I 0
En otras palabras la velocidad de la salida del controlador es proporcional al error o la rata de cambio de la
salida del controlador es proporcional al error. Por lo cual la tendencia es a minimizar el error.
Pivote
Aceite a
presión
Flotador
Válvula
Entra
agua h
Sale agua
Esta se define como el cambio lineal en el error, cuando la respuesta proporcional iguale a la derivativa.
Se acostumbra expresar la acción derivativa en minutos de adelanto, que representa el tiempo en minutos
con que la acción derivativa se anticipa al efecto de acción proporcional.
La principal ventaja de este control es que produce señales de adelanto que actúan rápidamente cuando la
variable controlada cambia bruscamente.
m t K M s
Pet
t
et ó F s K KI
K dt E s s
I 0
P
Este es un esfuerzo para combinar las ventajas del controlador proporcional más el controlador integral, es
decir tiene una buena respuesta transitoria por el efecto proporcional y corrige el error en estado establo
por el efecto integral.
I
KI
Esta se define como el tiempo necesario para que la respuesta integral iguale a la proporcional después de
un cambio en escalón del error.
m t K M s
e t
t
et dt KD det ó F s K KI K s
s
E s
P I D
K
0
dt P
Este controlador ofrece rápida respuesta proporcional al error, mientras que tiene un reajuste automático
desde la parte integral que elimina el error en estado estable. La acción derivativa permite que el controlador
responda rápidamente a cambios en el error.
Generalmente en los controladores PID se puede ajustar:
KP en % de la banda proporcional
TI en minutos de acción integral
TD en minutos de acción derivativa
Resumen
Control Función de Velocidad de Error en estado Uso solo Costo
transferencia respuesta estable
Proporcional KP Media Existe Si Bajo
(P)
Derivativo (D) KDs Alta Existe No Medio
Integral (I) KI Baja No hay No Medio
s
PI KI Media No hay Si Alto
K
P
s
PD KP KDs Alta Existe Poco Alto
PID KI Alta No hay Si Alto
K K s
P D
s
KM
Si se supone la referencia R(s) es constante, el elemento que puede afectar al sistema es la perturbación.
Por tanto por comodidad podemos tomar R(s) = 0, esto implica que la respuesta del sistema C(s) debe
tender a cero para cualquier perturbación P(s) en el transcurso del tiempo.
En estas condiciones puedo transformar el diagrama de bloque del sistema completo en:
Ps 1 Y s
s 1
KM KV F s
Ps 1 Y s
s 1 KM KV F s
Si el controlador es proporcional
F s KP
La función de transferencia del sistema completo será:
Y s 1
Ps s 1 KM KV KP
Donde para simplificar se puede escribir Luego la función de transferencia será:
KM KV K P 1 K
Y s 1
Ps s K
H H yt
K Error
t
0
et ut K
Myt
mt K P et
Si suponemos que la perturbación cambia bruscamente y luego se mantiene constante en el tiempo, lo cual
está representado por un escalón:
pt H Ps
H
s
El valor en estado estable será:
YEE sH H
lim
s0 K s K s K
Según lo supuesto al principio el valor deseado para la salida en estado estable debería ser cero, puesto que
nuestra señal de referencia es cero, luego este valor en estado estable indica que con este tipo de
controlador se produce un error en estado estable que es igual a H/K.
p t
H H yt
K Error
t
e t u t K
Myt
0
La diferencia entre una acción de control proporcional y una proporcional derivativa, se puede apreciar más
claramente si suponemos que el error varía en forma de rampa (esta suposición es meramente ilustrativa ya
que en un caso real es muy difícil obtener un error variando de esta forma, en todo caso no corresponde a
una perturbación variando en forma de rampa). Podemos observar si graficamos el error y la salida del
controlador lo siguiente:
mt K P et K D
Det e t u t K M y t
0
t
F s K P K I
s
La función de transferencia del sistema completo será:
Y s 1 1
Ps Ks 1 K KI Ks 1 KM KV KI
K K K
M V P
s M V P
s
Y s s
2
Ps s 1 K K K s K K K
M V P M V I
En este caso el sistema se vuelve de segundo orden, la función de transferencia será de la forma:
Y s s
Ps s2 2 s
2
n n
Donde:
2n
1 KM KV K P ; 2 K M KV K I
n
Donde: D2 y 2 n Dy n 2 y 2
p'
KM KV KI
p' 4 Dp
n
Por lo tanto el sistema es un sistema de segundo orden y responde como tal. Suponiendo que en las
condiciones iniciales la salida toma el valor de la señal de referencia:
t 0 Y0 0
La respuesta del sistema podrá tomar cualquiera de las cuatro formas posibles en sistemas de segundo
orden
La gráfica de la respuesta de este sistema será, para el caso de una perturbación (p) en escalón:
1
p't
Si
pt
H yt
0
p't
Si
1
p t
y t
H
0 t
e t u t K M y t
m t K etdt
et K
I t
P
0
F s K P K I
s KDs
La función de transferencia del sistema completo será:
Y s 1 K K1M KV KI
Ps Ks 1 K K I K s Ks 1 K K K K s
K
M V P D
s M V P
s
M V D
Y s s
Ps KM KV KD s 2 1 K K K s K K K
M V P M V I
En este caso el sistema se vuelve de segundo orden, la función de transferencia será de la forma:
Y s KM KV KD s
Ps s2 2 s 2
n n
Donde:
2
1 KM KV KP ; 2 KM KV KI
KM
n n
K K
V D
KM KV KD
Estudio de estabilidad del sistema
Como el sistema es de segundo orden su ecuación característica es de la forma:
s2 2 ns 2 0
Este sistema tiene dos raíces, y la estabilidad y forma de la respuesta dependen del valor de £.
Se observa entonces aquí que la principal diferencia que aparece con respecto a la acción de control
proporcional más integral es que el radio de amortiguamiento y la frecuencia natural son diferentes, por lo
cual la parte transitoria será probablemente más corta y la respuesta estará más cerca del valor deseado.
Sin embargo al afectar el valor del radio de amortiguamiento esta acción podría hacer tender al
sistema hacia la inestabilidad si este valor se acerca mucho a cero, por lo cual se debe tener
cuidado con la calibración del tiempo de acción integral.
La gráfica de la respuesta de este sistema será, para el caso de una perturbación (p) en escalón:
1
p't
Si
p t
H yt
t
0
et ut K
Myt
0
0
p't
Si
1
pt
yt
H
0 t
e t u t K M y t
m t K etdt K D Det
e t
I t
P
0
K
Tema 10. Ajuste de Controladores PID 159
Introducción
En este tema se describen algunos métodos de determinación los valores de los parámetros de un
controlador PID que son la ganancia proporcional (Kp ), el tiempo de acción derivativo (Td ) o el tiempo de
acción integral (Ti). A este proceso de selección de los parámetros del controlador, para que cumpla
las especificaciones de operación, se le conoce como ajuste o entonación del controlador. Este ajuste se
hace con el objetivo que el sistema controlado cumpla con las especificaciones de respuesta
transitoria y de estado estacionario que se requieran.
Existen diversos métodos para ajuste de controladores que difieren en función al conocimiento que se
requieren sobre la dinámica del proceso a controlar.
Cuando se conoce un modelo matemático del proceso, entonces se pueden aplicar métodos analíticos para
determinar los parámetros del controlador.
Pero si el proceso es muy complejo, suele ser difícil la obtención de un modelo matemático, en cuyo caso
no es posible la utilización de un método analítico de cálculo de los parámetros del controlador PID y se
debe recurrir a procedimientos experimentales para el ajuste.
Los procedimientos convencionales analíticos, por lo general están restringidos a proceso que se comportan
como sistemas lineales, de una entrada y una salida, e invariantes en el tiempo.
Estos valores para los parámetros del controlador permiten obtener un control que regula en forma estable
y con una buena velocidad de respuesta al proceso a controlar. Sin embargo como estos valores se obtienen
partiendo de la aproximación antes mencionada y el proceso real no corresponde exactamente a ese modelo
entonces generalmente es conveniente realizar luego un ajuste fino. Con esto se busca mejorar
la estabilidad relativa que se traduce en el aumento del margen de fase y ganancia.
El método diseñado por Ziegler y Nichols conduce a fórmulas simples para obtener los parámetros del
controlador, en términos de la ganancia critica y el periodo crítico, como se puede ver en la tabla siguiente.
La tabla también muestra una estimación del periodo de la dinámica dominante (Tp) del sistema en
lazo cerrado.
Parámetros de controladores PID según método de la respuesta frecuencial de Ziegler - Nichols
Controlador Kp Ti Td Tp
P 0.5 Kc - - tc
Pl 0.4 Kc 0.8 tc - 1.4 tc
PID 0.5 Kc 0.5 tc 0.12 tc 0.85 tc
Polos Dominantes
U(s) H(s) G(s) Y(s)
Considere el sistema mostrado en la figura, donde
G(s) es la función de transferencia del proceso y
H(s) es la función de transferencia del controlador.
La función de transferencia del sistema en lazo
cerrado (Gc(s)) está dada por:
GsH
G s s
c
1 GsH s
Muchas de las propiedades de los sistemas en lazo cerrado pueden ser deducidas a partir de los polos y
ceros de Gc (s), estos últimos son también los ceros de G(s)H(s). Los polos en lazo cerrado se pueden
obtener a partir de la raíz de la siguiente ecuación:
1 GsH s 0 Im
La ubicación en el plano imaginario de los polos y ceros de un P1
sistema en lazo cerrado puede variar considerablemente con
respecto al lazo abierto. Muchos lazos de retroalimentación
simple, tendrán una configuración del tipo mostrado en la figura P3 Z1
siguiente, donde las principales características de la respuesta Re
están dadas por un par de polos complejos, P1 y P2 ,
denominados polos dominantes. La respuesta es también
influenciada por los otros polos y los ceros P3 y Z1 , P2
respectivamente. También pueden existir otros polos y ceros
alejados del origen. Los polos y ceros cuyas partes reales son mucho
más pequeñas que la parte real de los polos dominantes, tienen
pequeña influencia en la respuesta transitoria, por lo cual este efecto puede despreciarse.
Control PI
El método de diseño por asignación de polos dominantes será aplicado primero a un control PI. Dos polos en
lazo cerrado se pueden especificar en este caso, ya que un controlador PI tiene dos parámetros ajustables.
Con el controlador PI dado por la ecuación:
H s K p Ki
s
Donde Kp es la ganancia proporcional y Ki es la ganancia integral. Los parámetros Kp y Ki serán
determinados de tal forma que el sistema en lazo cerrado tenga los polos dominantes en s = P1 y s = P2 ,
tendremos entonces:
P1 0 1 2 i
i0
1 2 i
P i
1 0 0
Esto implica que:
Ki
1 K p G P1 0
P1
K
1 K i G P 0
p 2
P2
La ecuación anterior es lineal en Kp y Ki . Esta tiene una solución si |G(P1 )| G 0 y la solución en este
caso es:
1 2 A B
KP0 0 0
1 2 A 2 2
B
K 0 0
0B0 0
Donde:
2 B 2
i 0 1
2
A
0
A0 ReGP1
B0 ImGP1
El parámetro ω0 puede ser visto como un parámetro que permite determinar la velocidad de respuesta del
sistema de lazo cerrado obtenido. Para valores pequeños de ω0 se obtiene un sistema lento, y para valores
altos se obtiene un sistema rápido.
Si la dinámica de un proceso es de primer orden, el sistema en lazo cerrado tiene solamente
dos polos y por lo tanto el parámetro ω0 puede ser elegido arbitrariamente.
Para sistemas de orden elevado, en lazo cerrado se tendrá un número elevado de polos, pero se
escogerán de tal forma que dos sean dominantes.
Para obtener un sistema estable con polos no dominantes sobre el eje real, estos polos deberán ser
de parte real más pequeños que – £ω0 . La condición para obtener dos polos dominantes es que ω0
debe ser escogido suficientemente pequeño. El límite superior de ω0 puede ser determinado usando
la condición de que el polo más alejado sobre el eje real se encuentre en el plano imaginario en s =
—αω0 .
Para un proceso estable, la función A(ω0) es positiva y B(ω0) es pequeño para un ω0 pequeño.
Entonces resulta que la ganancia proporcional Kp(ω0) es negativa para un ω0 pequeño. Y como
es normalmente deseable obtener las ganancias del controlador positivas, un límite inferior para el
parámetro de diseño se puede obtener de la condición KP(ωoI) = 0. El valor de ω0I corresponde
en este caso a un control integral puro. Una alternativa para elegir a ω0 , basado en los
polos dominantes, es seleccionar un ω0 que genere el valor máximo de la ganancia integral. Los
valores obtenidos en este caso serán muy cercanos a los obtenidos para la condición de polos
dominantes.
Control PD
El método de diseño por asignación de polos dominantes puede ser aplicado al control PD. La función
de transferencia de este controlador es:
H s K p Kd s
Y de igual forma que para el caso del control PI se requiere que el sistema en lazo cerrado tenga polos en P1
y P2 . Con cálculos análogos al caso del controlador PI los parámetros del controlador se pueden obtener
mediante las expresiones:
1 2 A B
0 0
KP 0
1 2
A
2
2
B
K 0 0
B 0
2
d 0
0 1 2 A 2
B 0 0
El valor del parámetro de diseño ω0 debe ser aquí aún mayor que para el controlador PI. Pudiéndose estimar
su valor inferior mediante la expresión KP(ωOD) = 0. El valor de ω0D corresponde en este caso a un control
derivativo puro. Un valor razonable para este parámetro es escoger el que proporciona el mayor valor para
la ganancia proporcional.
Control PID
Con un control PID, es posible posicionar tres polos en lazo cerrado. Dada la función de transferencia
del controlador PID como:
K'
H s K i K ' d s
'p s
Donde K' es la ganancia proporcional., K' es la ganancia para el término integral, y K' es la ganancia para el
p i d
término derivativo. Dos polos en lazo cerrado serán posicionados de acuerdo a la ecuación obtenida para el
control PI. Y se asume que con esto el problema de diseño del controlador PI está resuelto, es decir, que se
conocen los valores de KP (ω0 ) y Ki (ω0 ). El valor de la función de transferencia del controlador H’(s) en
P1 = —σ + iω es:
Suponiendo que la función de transferencia tiene el mismo valor que la función de transferencia para el
controlador PI, se tiene:
d
K ' p K 'i
2
K
0
Kp
'
0
Ki 2
K 'i Ki
2 K 'd
0 0 2
Por lo tanto se tiene que:
K 'p 0 K p 0 20K
'd K 'i 0 Ki 0 0 K
2
'd
Así, se obtiene una familia de ganancias, con dos parámetros (ω0 , K 'd), para un controlador PID, lo cual
genera un sistema de lazo cerrado con polos en s = P1 y s = P2 . El parámetro Kd' será ahora determinado
de manera que el sistema de lazo cerrado también tenga un polo en s = —
ω0 . Con esto se obtiene:
K 'i
1 K ' p K 'd 0 G 0 0
0
Introduciendo las expresiones de KP(ω0) y Ki(ω0) se obtiene:
Ki
1 K p 2K 'd 20K 'd G 0 0
0
De aquí, si G(—ω0) G 0, se obtiene:
K
i 0
1 K p 0 G 0
K 'd 0
0
20 1 G 0
Con esta ecuación queda definida una familia de controladores PID con un parámetro de diseño (ω0), con los
cuales se obtiene un sistema de lazo cerrado con tres polos dominantes ubicados en – £ω0 ± iω0√1 — £2 y
– ω0 . Valores pequeños de ω0 producen un sistema con velocidad de respuesta lenta y valores grandes de
ω0 producen un sistema con velocidad de respuesta rápida.
Ejercicio
Se tiene un proceso con la siguiente función de transferencia:
s 1
G
1 1 1
1 1.26s 2
3
4
0.8
0.6
0.4
0.2
0 0.5 1 1.5 2
GsH s Kp
1 s1 0.2s1 0.05s1 0.01s
La figura siguiente muestra los diagramas de Nyquist para el sistema con una ganancia proporcional igual a
10 y a 25.
Las figuras muestran que:
La ganancia crítica esta alrededor de Kp = 25
6 15
4 10
2 5
Imaginary Axis
Imaginary Axis
To: Y(1)
To: Y(1)
0 0
-2 -5
-4 -10
-6 -15
-8 -20
-2 0 2 4 6 8 10 -5 0 5 10 15 20 25
Para Kp = 10 Para Kp = 25
Con estos valores se busca en la tabla los valores para las ganancias del controlador obteniéndose los
siguientes:
Para un controlador PI:
Kp = 10 y Ti = 0.5
Para un controlador PID
Kp = 15, Ti = 0.31 y Td = 0.08
Obsérvese que los valores no son idénticos en los dos métodos, esto es debido a que son fruto de
aproximaciones distintas, y se trata en todo caso de una aproximación en donde se puede realizar
posteriormente un ajuste más fino. Sin embargo la magnitud de los valores es relativamente próxima en los
dos casos.
170 Introducción al Control de Procesos para Ingenieros
El principio básico del control en lazo cerrado clásico se muestra en la siguiente figura:
Perturbación
Acciones de Control
Sistema de Control Error
Referenci Controlador
Sistema
E.F.C.
dinámico a
controlar
Elementos
Medidas de Salidas
medición
En este las acciones de control son calculadas por la medición de las salidas (variables controladas),
este tipo de control suele denominarse control por retorno de salida.
Mientras que por su lado el principio básico del control por retorno de estado se basa en realizar las
acciones de control en función de los valores medidos de los estados. Este se puede apreciar en la siguiente
figura:
Jean-François DULHOSTE
Como ya estudiamos anteriormente, en un sistema expresado en espacio de estado, los polos de la función
de transferencia son los valores propios de la matriz de estado A, por lo tanto una retroalimentación sobre
los estados con una matriz K, es equivalente a realizar una modificación de la matriz de estado del
sistema, tal como se explica a continuación.
En consecuencia la matriz de estado del sistema en lazo cerrado es: ALC = A — BK.
La dinámica del sistema en lazo cerrado estará determinada entonces por de los valores propios de la
matriz (ALC)
Estos valores propios como bien sabemos son las raíces de la ecuación característica λÆ–BK(s) , que se
obtienen con:
|λI – ALC| = 0
El control modal
También denominado control por ubicación de polos. Este se puede realizar en espacio de estado, o en
forma algebraica con las funciones de transferencia.
Se denomina control modal al control que consiste en determinar una matriz de retroalimentación de
estado K tal que los valores propios de la matriz ALC se ubiquen en posiciones prefijadas
arbitrariamente (λ0, λ 1 ,⋯, λn–1) (valores reales o complejos conjugados). La existencia de una solución
depende de la controlabilidad del sistema y de la posibilidad de acceder a los estados del mismo.
Donde (λ0, λ 1 ,⋯, λn–1)) son los polos que se quieren imponer al sistema en lazo cerrado.
Otra forma de calcular la matriz K consiste en efectuar un cambio de variable lineal al sistema para obtener
una forma canónica de controlabilidad.
Notas:
Para que el control se pueda realizar físicamente se deben escoger valores propios reales o
complejos conjugados por pares, lo cual garantiza una función de transferencia de
coeficientes reales.
Como la estabilidad es la primera característica que se debe asegurar estos coeficientes deben
tener su parte real estrictamente negativa.
Ejemplo 1:
Consideremos el sistema LTI definido por las matrices:
1 4 —1 2
A = [6 —1 3 ] ; B = [ 3 ] ; C = [1 0 0]
2 2 —5 —1
Queremos realizar un control por retorno de estado para obtener un sistema en lazo cerrado cuyos polos se
ubiquen en λ1 = —1, λ2 = —1 — 2i, λ3 = —1 + 2i.
s3 + (2k1 + 3k2 — k3 + 5)s2 + (25k1 + 21k2 + 10k3 — 29)s + (41k1 + 72k2 + 71k3 — 129)
= s3 + 3s2 + 7s + 5
k1 2 3 —1 –1 3 5
[k2] = [25 21 10 ] ([7]— [ —29 ])
k3 41 72 5 —129
71
Cuando el sistema se encuentra en la forma canónica de controlabilidad, las matrices A y B tienen formas
particulares que permiten facilitar el cálculo de la matriz K:
0 1 0 ⋯ 0 0
F
I 0 0 1 01 IF 1I
0 0 ⋱ ⋱ ⋮ I B = 0⋮
A= I I II
I ⋮ ⋱ 1I I0I
L—an —an–1 ⋯ ⋯ a1 ے L1ے
C = [bn bn–1 ⋯ ⋯ b1]
Como lo que se quiere determinar es una matriz K = [k1 k2 ⋯ kn–1], tal que:
0 1 0 ⋯ 0 0 F 0 1 0 ⋯ 0 1
F 0 1I FI01I 0 0 1 0
I 0 0 1 I I
ALC = 0 0 ⋱ ⋱ ⋮ — ⋮ [k1 k2 ⋯ kn–1] = 0 0 ⋱ ⋱ ⋮
I I I I I I
I ⋮ ⋱ 1 I I0I I ⋮ ⋱ 1 I
L—an —an–1 ⋯ ⋯ —a1 ےL1ے L—an — k1 —an–1 — k2 ⋯ ⋯ —a1 — kn–1 ے
El control modal impone que el denominador de la función de transferencia del sistema en lazo cerrado sea:
den(HBF (s)) = (s– λ0)(s — λ1).. . (s — λn–1) = sn + a' sn–1 + . . . +a' s + a′n
1 n–1
La ubicación de los polos en este caso no modifica la forma de la representación (sigue siendo en lazo
cerrado de forma canónica de controlabilidad). En consecuencia, obtenemos dos escrituras distintas para la
matriz de estado del sistema en lazo cerrado (A — BK):
0 1 0 ⋯ 0
0 1 0 ⋯ 0
IF 0 0 1 0 I1 IF 0 0 1 0 1I
ALC = I 0 0 ⋱ ⋱ ⋮ =
I I 0 0 ⋱ ⋱ ⋮ I
I ⋮ ⋱ 1 I I ⋮ ⋱ 1 I
L—an — k1 —an–1 — k2 ⋯ ⋯ —a1 — kn–1 ے L—a′n —a′n–1 ⋯ ⋯ —a′1 ے
sencilla. Ejemplo 2:
Para el sistema del ejemplo 1 la función de transferencia es:
2s2 + 25s + 41
G(s) = 3
s + 5 s2 — 29 s — 129
Se quiere un control por retorno de estado para que el sistema en lazo cerrado tenga los polos en λ1 =
—1, λ2 = —1 — 2i, λ3 = —1 + 2i ; lo cual es equivalente a tener un sistema cuya ecuación característica
sea:
Ec = (s + 1)(s + 1 + 2i)(s + 1 — 2i) = (s + 1)(s2 + 2s + 5) = s3 + 3s2 + 7s + 5
Por lo tanto la matriz K se obtiene con:
0 1 0 0 1 0
ALC = [ 0 0 1 ]=[0 0 1]
129 — k1 29 — k2 —5 — k3 —5 —7 —3
K = [134 36 —2]
Ejemplo 3:
Para el sistema de segundo orden siguiente:
Y(s) 1
G(s) = =
U(s) s2
Se quiere obtener mediante un control por ubicación de polos, un sistema en lazo cerrado cuyos polos se
encuentren en:
λ = —1 ± i , es decir que su ecuación característica sea:
q(λ) = λ2 + 2λ + 2
Se puede obtener una representación de estado del sistema escogiendo como estados a: x1 = y; x2 = y˙,
con los cuales obtenemos para la representación
0 1 de estado
0 las matrices siguientes:
A=[ ] ;B = [ ] ;C = [1 0]
0 0 1
La matriz de controlabilidad es:
0 1
C = [B AB] = [ ]
1 0
Y su inversa:
1
C–1 = 0 —1 0 1
[ ]=[ ]
—1 —1 0 1 0
Por otro lado:
0 1 0 1 1 0 2 2
q(A) = [ ]2 + 2 [ ]+2[ ]=[ ]
0 0 0 0 0 1 0 2
Inserción de un pre-compensador
Cuando el objetivo de control es obtener: lim y(t) = ycdonde y ces el valor de la referencia (consigna o
t→∞
valor deseado), que supondremos constante en el tiempo (Escalón). Se requiere calcular el valor de la
entrada necesaria v para obtener como salida ese valor deseado yc .
Este resultado nos muestra que solo requerimos de agregar un bloque con una matriz que multiplique a la
entrada real del sistema, es decir simplemente una corrección estática como la mostrada en la siguiente
figura:
yc v u
N x˙ x C y
_
w
v u x˙ x C y
—w —yc
K1
K2 z z˙
ƒ
En este caso el objetivo es obtener: lim y(t) = y cdonde y ces el valor de referencia considerado constante
t→∞
en este caso, y v se comporta como una perturbación.
x
Lo que se busca es un control u = —[K1 K2] [ z] que anule el error estático de la salida del sistema.
El sistema en lazo cerrado Ben régimen estacionario
xs será:
A 0 xs [K K ] B 0 v A — BK1 —BK2 xs B 0 v
0= [ ] [z ] — [ ] 1 2 [
z ] + [ ] [y ] = [ ] [ ] + ] [y ]
C 0 s 0 [ z
xs s 0 —1 c C 0 s 0 —1 c
{y s = [C 0] [ ]
zs
Para que este sistema aumentado tienda a un estado de equilibrio es suficiente y necesario que la matriz
A — BK1 —BK2
del sistema en lazo cerrado ALC = [ ] sea estable. Razonamiento que es válido para toda
C
entrada v.
Exemple 4:
Sea el sistema LTI definido por las matrices:
1 4 —1 2
A = [6 —1 3 ] ; B = [ 3 ] ; C = [1 0 0]
2 2 —5 —1
Se quiere realizar un control por retorno de estado para obtener un sistema en lazo cerrado cuyos polos se
ubiquen en: λ1 = —1, λ2 = —1 — 2i, λ3 = —1 + 2i.
La matriz del sistema en lazo cerrado es:
1 4 —1 2 2
A — BK1 —BK2 [6 —1 3 ] — [ 3 ] [k11 k12 k13] —[ 3 ] k2
ALC = [ C 0 ]=[ 2 2 —5 —1 1 ]
[1 0 0] 0
1 — 2k11 4 — 2k12 —1 — 2k13 —2k2
6 — 3k —1 — 3k12 3 — 3k13 —3k2
ALC = [ 2 + k1111 2 + k12 —5 + k13 k2 ]
1 0 0 0
Para que esta matriz sea estable es necesario y suficiente con escoger los Kij tales que los valores propios
de la matriz sean negativos o con parte real negativa, tres de estos valores propios están definidos por
el enunciado del problema:
|sI — ALC| = (s + 1)(s + 1 + 2i)(s + 1 — 2i)(s — λ4) = (s3 + 3s2 + 7s + 5)(s — λ4)
|sI — ALC| = s4 + (3 — λ4)s3 + (7 — 3λ4)s2 + (5 — 7λ4)s + 5λ4
Si por ejemplo escogemos como cuarto valor propio (el correspondiente a la parte integral) λ4 = —4:
|sI — ALC| = s4 + 7s3 + 19s2 + 33s + 20
Control Parcial
Consideremos ahora un sistema que es solo parcialmente controlable, es decir rango(A) = rÆ < n. En
este caso rÆ designa el grado de controlabilidad del sistema, es decir el número de estados controlables que
este posee. Para el sistema:
x˙ = Ax + Bu
Seleccionamos como vector de estado x un vector propio del sistema. Y separamos el vector de estado en
dos partes: la parte controlable con rÆ variables de estado: x1 y la parte no controlable con (n – rÆ)
variables de estado: x2 .
Podemos entonces escribir la ecuación dinámica del sistema como:
x˙1 A11 A12 x1 B1
[ ]= [ ] [x ] + [ ] u
x˙2 0 A22 2 0
Los valores propios A (modos del sistema) son los valores propios de A11 junto a los valores propios de A22 .
|sI — A| = |sI — A11| × |sI — A22|
Notas:
Obsérvese que el control no tiene ningún efecto sobre la evolución de la parte no controlable del
sistema pues:
x˙2 = A22x2
Para un sistema con partes no controlables, solo los modos de la parte controlable (A11) se pueden
modificar mediante un control por retorno de estado.
La parte K2 de la matriz K se determina de una forma distinta. Se pueden considerar para ello
criterios adicionales como por ejemplo el hecho de hacer que la salida del sistema sea independiente
de la parte no controlable x2 .
Ejercicios:
Para los siguientes sistemas determine el control por retorno de estado requerido para estabilizarlos,
un control directo para anular el error estático y un control integral.
0 1 0
A=[ ] B=[ ] 0 1 0 0
A = [0 0 1] B = [0]
1. 2 —3 1; ;
2 —7 4 1
C = [ 1 2]
C = [0 1 2]
2.
Observadores
En algunos casos no es posible acceder a la totalidad de los estados, ya sea porque es imposible
realizar una medición física de los mismos o porque colocar medidores para cada estado haría muy
costoso o complejo el sistema. En este caso la implementación del control u = —Kx no se puede
realizar directamente. Y como la matriz C no es invertible en la mayoría de los casos, la ecuación de salida
y = Cx no permite determinar el valor de x.
Se requiere en este caso la reconstrucción de los estados x a partir de la información disponibles para
la entrada u y la salida y para ello se utiliza un observador, el cual está constituido por un sistema dinámico
que permite obtener una estimación aproximada x^ de x. Este observador suele también
denominarse reconstructor, estimador o filtro.
— Si z y x tienen la misma dimensión, entonces se dice que el observador es completo (se estima todo el
vector de estado). En ese caso T = I ; por tanto z = x y z^ = x^.
—Si dim(z) < dim(x), entonces se dice que el observador es de orden reducido.
Observador Identidad
Un observador identidad es un observador completo sin desvío ( z^ (t) → z si t → ∞) donde z = x. Con
este observador obtenemos las ecuaciones de la pareja sistema-observador siguientes:
x^˙ = Fx^ + Ly + Ju
{x˙ = Ax + Bu
y = Cx
Esta relación que caracteriza al observador puede representarse por el siguiente diagrama de boques:
El cálculo de un observador para un sistema LTI consiste entonces en calcular la matriz L, para obtener unos
valores propios del sistema retroalimentado (observador) en una posición seleccionada arbitrariamente. Se
dice entonces que se trata de un problema de control modal del sistema dual.
Ejemplo 5:
Consideremos el sistema lineal x˙ = Ax + Bu; y = Cx definido por las matrices:
2 3 0
A=[ ] B=[ ]
—1 4 1
[
C= 1 0 ]
Si suponemos que solo la salida y = x1 puede medirse, requerimos de un observador para determinar el
estado no medible x2 . En este caso aunque solo requerimos determinar uno de los dos estados, vamos a
construir un observador que permite determinar todos los estados, es decir un observador completo. Sería
posible sin embargo la construcción de un observador reducido, sin embargo con frecuencia es preferible
construir un observador completo y utilizar los estados redundantes para corregir eventuales defectos de
estimación causados por ruido u otras perturbaciones. Un ejemplo de esto sería el filtro de Kalman
(observador optimo de tiempo variable) el cual permite la resolución de problemas de observación con ruido
en la medición.
por lo tanto el sistema es observable, y podemos entonces construir un observador para determinar los
estados del sistema.
Para la construcción del observador vamos a seleccionar una dinámica tal que el comportamiento de éste
sea equivalente a un sistema de segundo orden cuya ecuación característica sea:
p(λ) = λ2 + 2ξωnλ + ω2n
Donde: ξ = 0.8 y ωn = 10.
Estos valores nos permiten obtener un tiempo de estabilización del observador de Ts = 0.5 s
Repaso: El tiempo de estabilización es el tiempo necesario para que la salida de un sistema sea igual al 98%
del valor de la respuesta en estado estable (2% de error) después de un cambio en escalón de la
entrada, este se puede calcular para un sistema de primer orden con: Ts = 4τ y para un sistema de
segundo orden
4
con: Ts = .
(ωn
La ecuación característica del observador será: 1 0 2 3 L1
p(λ) = |λI — (A — LC)| = |λ [ ] — ([ ] — [ ] [1 0])|
0 1 —1 4 L2
λ — 2 + L1 3 | = (λ — 2 + L1)(λ — 4) — (—3)(1 + L2)
p(λ) = |
Por lo tanto: 1 + L2 λ—4
Control por retorno de estado. Se calcula la ganancia del controlador suponiendo que todo el
vector de estado x puede medirse. El control lineal es de la forma u = —Kx, donde K se determina
imponiendo los polos que se quieren obtener en lazo cerrado.
Reconstrucción de los estados. Si solo es posible medir y, se necesita determinar un observador, es
decir se necesita calcular la ganancia L que asegure la estabilidad del observador y una
observación sin desvío.
Control por retorno de salida. El control del sistema se realiza finalmente a partir del
estado estimado u = —Kx^.
A pesar de que la ganancia K se calcula para garantizar la estabilidad y rendimiento del sistema, y para ello
se escogieron de forma que las raíces de la ecuación característica |λI – (A — BK)| = 0 sean negativas o
con parte real negativa. Y por lo tanto x(t) → 0 cuando t → ∞. Es necesario verificar que ésta condición se
garantiza aún con los estados estimados x^.
Para mostrar la estabilidad del sistema vamos a estudiar las dinámicas del error del sistema que se
introduce con la sustitución de x por x^ en el cálculo del control.
Si reemplazamos el valor de la entrada u en la ecuación del observador obtenemos un sistema con entrada y
y con salida u :
˙ = (A — LC — BK)x^ +
{x^
Ly u = —Kx^
Lo valores propios de la nueva matriz dinámica del sistema son por lo tanto las raíces de:
|λI — (A — BK)||λI — (A — LC)| = 0
Por lo tanto:
si las raíces de |λI — (A — BK)| = 0 son negativas o con parte real negativa, lo cual está
garantizado pues se calculó el control para garantizar esa condición,
y las raíces de |λI — (A — LC)| = 0 son negativas o con parte real negativa, lo cual está
garantizado pues el observador fue calculado para garantizar esa condición,
entonces los valores propios de la matriz dinámica serán negativos o con parte reales negativas.
En conclusión los valores propios del sistema en lazo cerrado son los valores propios de (A — BK), es decir
las relativas al control del sistema, mas los valores propios de (A — LC), es decir las del observador, y en
consecuencia se observa que la sustitución de x por x^ no modifica los valores propios obtenidos en el
cálculo del control. Solamente se superponen valores propios adicionales que provienen del
observador. La estabilidad del sistema en lazo cerrado no es modificada por la presencia del observador,
siempre que este sea estable y no tenga desvío. Esto es conocido bajo el nombre de principio de
separación el cual se satisface para todos los sistemas lineales.
Para que el comportamiento del sistema en lazo cerrado no sea modificado de manera notable por la
presencia del observador, es necesario que la reconstrucción del estado sea rápida respecto de la dinámica
del sistema de lazo cerrado (polos de (A — LC) de gran modulo respecto de (A — BK)).
Inserción de un pre-compensador
De igual forma que para el control por retorno de estado, el control por retorno de salida no garantiza que
el vector de la señal de referencia v sea nulo. El estado del sistema x va a converger hacia 0 con una
dinámica determinada por los
polos ya ubicados. x˙ = Ax + Bu y = Cx y
yc N v u
Cuando v ya no sea nulo,
el estado converge hacia
un valor que no es -
necesariamente nulo. Un
w x^˙ = (A — LC)x^ + Ly + Bu
pre-compensador es en
este caso una matriz
cuadrada N, que se K x^
ubica justo después del Regulador
vector de la referencia,
tal como se observa en
la figura.
Este pre-compensador no cambia los polos del sistema en lazo cerrado. Simplemente permite hacer
corresponder ciertos componentes de la referencia con ciertas variables de estado previamente
seleccionadas.
En este caso el control es:
u = v — Kx^
Por otro lado si vamos al sistema que queremos controlar x˙ = Ax + Bu; y = Cx y se sustituye la ley
de control en la ecuación dinámica obtenemos:
x˙ = Ax + Bv — BKx^
Como lo que queremos es estudiar la dinámica del error sustituimos x^ = x — e :
x˙ = Ax — BKx + BKe + Bv = (A — BK)x + BKe + Bv
Como |A — LC| G 0 pues la matriz L fue seleccionada para que el observador sea estable, es decir con polos
estrictamente inferiores a cero, entonces obligatoriamente e = 0 (error de observación nulo en régimen
estacionario).
La cual es idéntica a la ecuación de control por retorno de estado directo. Por lo tanto la salida y
tiende hacia la referencia yc si calculamos la matriz N con la expresión:
N = [C(—(A — BK)–1B)]–1
Ejemplo 7:
Para el sistema presentado en el ejemplo 6, representado por las matrices:
2 3 0
A=[ ] B=[ ]
—1 4 1
C = [1 0]
Se quiere realizar un control por retorno de salida, ya que el estado x2 no se puede medir físicamente, con
el objetivo de obtener un sistema en lazo cerrado con polos en λ1 = —1 + 2i, λ2 = —1 — 2i.
Para este sistema en el ejercicio 6 fue calculada la ganancia de un observador que recordamos aquí:
L1 22
L=[ ]=[ ]
L2 59
Para el control por retorno de salida se requiere calcular la ganancia del controlador, igual como si fuese un
retorno de estado simple: s—2 —3
1 0 2 3 0
|sI — A + BK| = |s [ ]—[ ] + [ ] [k1 k2]| = |[ ]|
0 1 —1 4 1 1 + k1 s — 4 + k2
s—2 3
|[ 1 + k1 s—4+
k2 ]| = (s — 2)(s — 4 + k2) — (—3)(1 + k1) = (s + 1 + 2i)(s + 1 — 2i)
Paras realizar el control por retorno de salida es suficiente con realizar un retorno de estado:
˙ = (A — LC)x^ + Bu +
{x^
Ly u = v — Kx^
–1
2 3 –1 0 —0.8 —0.6 0 –1
N = [[1 0] (— ([ ] ) [ ])] = [[1 0] (— [ ][ ])]
—13/3 —4 1 0.8667 0.4 1
N = —1,6667
Ejercicios.
Para los siguientes sistemas calcular un observador para la estimación de los estados que se suponen
no medidos, y determinar un control por retorno de salida para estabilizar los sistemas, agregando
adicionalmente un pre-compensador para anular el error.
0 1 0
A=[ ] B=[ ]
1. 2 —3 1;
C = [1 2]
0 1 0 0
A = [0 0 1] B = [0]
2. 2 —7 4 ;
1
C = [0 1 2]
Tema 12. Introducción a los Sistemas Digitales 187
El control mediante un computador digital requiere de elementos adicionales en el sistema de control que
modifican parcialmente el lazo de control y requieren de conocimientos adicionales tradicionalmente
ligados al tratamiento de señales.
Elementos
Medidas de Salidas
medición
Control Retroalimentado Continuo
Controles Perturbación
Controles digitales
Sistema de Control analógicos
Sistema dinámico a
Referenc Error Computador
CDA E.F.C.
digital
Sensores / transmisores
Medidas CAD
Medidas Salidas
digitales analógicas
Sistema de control
u(t)
r(kT) e(kT) v(kT) v(t) y(t)
Kv G(s)
C(z) CDA
CAD
m(kT) m(t) H(s)
Convertidor Análogo-Digital
Un convertidor análogo digital realiza principalmente dos operaciones:
El muestreo de la señal, es decir su discretización temporal, con un período fijo que se denota
generalmente como Te .
Una señal discreta es una señal que posee un valor solo a unos instantes precisos. Estos
instantes de tiempo están definidos por el período de de muestreo Te , que es un valor de tiempo
constante que separa dos de estos instantes.
Un discretizador ideal se puede ver como un interruptor que se conecta cada Te segundos, durante
un instante muy corto de tiempo y esta desconectado el resto del tiempo. Si consideramos
como variable de entrada del discretizador a f(t) y como salida a
fe(t) tendremos como función que los relaciona a fe(t) =
f(t)ðTe o ðTe = ∞ fe(t)
k=–∞ ð(t — kTe) que no es más que la f(t)
función impulso ∑ o delta de Dirac. A cada paso de
discretización fe (t) = f(kTe )ð(t — kTe ). Continuo Discreto
Variable
Continua Discreta
x(t) x(t)
Continuo
Tiempo
t t
Señal Analógico Señal discreta tiempo continuo
x(t) x(t)
Discreto
t t
Señal Continua tiempo discreto Señal Digital
En la práctica un convertidor análogo digital no puede adquirir la señal de forma instantánea, por lo cual se
produce una diferencia entre la discretización ideal y la real,
puesto que el segundo se produce generalmente haciendo un
promedio del valor de la señal en un tiempo dado Tr . Esta
diferencia produce una distorsión en la señal discretizada,
la cual se demuestra a continuación. Consideremos el
Valordela
siguiente ejemplo:
Valor de la f(t) = A sin(ωt)
La señal discretizada ideal correspondiente es:
fe(kTe) = A sin(ωkTe)
Convertidor Digital-Análogo
Un convertidor digital-análogo por su parte tiene como objetivo obtener de nuevo la señal analógica
previamente discretizada equivalente de forma exacta para todo tiempo, es decir fr(t) = f(t) ∀t.
La respuesta al impulso de este tipo de sistema es una función rectángulo unitario de ancho Te .
Se habla de “bloqueador de orden cero” (BOZ), pues de hecho realiza una interpolación de orden cero de
los valores de la señal discreta (el único orden posible que sea causal, en consecuencia el único posible de
realizar en tiempo real).
La frecuencia de muestreo, la cual tiene una influencia sobre la precisión de la medida y el retraso
que agrega el bloqueador de orden cero.
La precisión del sistema numérico, el cual puede adoptar solo valores numéricos precisos, por
tanto tiene una influencia sobre el error de cuantificación (Amplitude Quantization Error).
Ecuaciones recurrentes
Un sistema numérico lineal invariante en el tiempo se define mediante una relación de la siguiente forma:
n m
Σ aiy(k — = Σ bj u(k — j)
i) j=0
i=0
Ejemplo 1
Consideremos las siguientes ecuaciones recurrentes de orden 2 (en ellas se consideran dos pasos de
muestreo previos):
2y(k) + y(k — 1) — y(k — 2) = u(k) — u(k — 1)
Para que se puedan resolver esta ecuación se requiere conocer dos condiciones iniciales, y la entrada para
todo paso de muestreo.
Si tenemos como condiciones iniciales: y(—1) = 1 e y(—2) = 1, y se tiene como entrada del sistema un
escalón unitario, es decir:
0 si k < 0
u(k) = { 1 si k ≥ 0
Esto es equivalente a la
expresión: y(k — 1 k+1
1)
y(k) = = ( )
2 2
Es de notar que la solución de esta ecuación recurrente depende de las condiciones iniciales.
La transformada en Z
Los sistemas analógicos lineales se pueden representar mediante la función de transferencia, que no es más
sino la transformada de Laplace de su respuesta al impulso.
0 0 k=–∞
Mediante algunas consideraciones y operaciones matemáticas, así como propiedades de la función impulso
obtenemos la siguiente expresión:
∞
Propiedades de la transformada en Z
Consideremos dos señales discretas f(k) y g(k). Notamos F(z) y G(z) las transformadas en Z
respectivas, α y β números reales cualesquiera.
Linealidad:
S(αf(k) + βg(k)) = αF(z) + βG(z)
Retraso:
S(f(k — n)) = z –n F(z) ∀n ∈ Գ
Avance:
194 Introducción al Control de Procesos para Ingenieros
n–1
Convolución:
S(f(k) ∗ g(k)) = F(z) ∗ G(z)
La transformada es:
∞ ∞ 1 z
F (z) = Σ f(k) z–k = Σ z–k = 1 + z–1 + z–2 + ⋯ = =
e
k=0 k=0 1 — z–1 z — 1
Jean-François DULHOSTE
Tema 12. Introducción a los Sistemas Digitales 195
1 t = kTe
ð(t — kTe) = { –kTes
e ð(k — kTe) z–k
0 t G kTe
Ejemplo 2:
Queremos obtener la representación discreta del
u(t) ue(t) y(t)
sistema BOZ H(s)
1
H(s) =
s(s + 1)
Tomando en cuenta para ello que en la discretización se
debe agregar un bloqueador de orden cero, con el cual U(z) Y(z)
G(z)
el sistema se convierte en:
( ) ( )
1 — e–sTe 1 1 — e–sTe
G0 s H s = =
s s(s + 1) s2(s + 1) 1 1 1
G(s) = G (s)H(s) = (1 — e–sTe) ( )
0 — +
s2 s s+1
Utilizando las propiedades de la transformada en Z y la tabla–sTde transformadas
1 1 tenemos:
1
G(z) = S(G (s))S(H(s)) = S(1 — e e) ∗ S ( )
0 — +
s2 s s+1
zTe z z
G(z) = (1 — z–1) ( — + )
2 –Te
z—1 (z —
Te 1) z —1 1 z — 1e
G(z) = ( )z( — + )
z (z — 1)2 z — 1 z — e–Te
Como la función de transferencia se define para condiciones iniciales iguales a cero, entonces con el
teorema del retraso obtenemos:
n m
Σ ai y(k — i) = Σ bj u(k — j)
También se puede obtener una relación equivalente con una función de transferencia discreta de la forma:
G(z) = Y(z) ∑m bj z j
= j=0
U(z) ∑ni=0 aizi
Ejemplo 4. Se tiene:
Y(z) b1z + b0
G(z) = U(z) = z2 + a1z + a0
Transformada inversa de Z
Existen 4 métodos que permiten determinar la inversa de una transformada en Z.
Donde Γ es un contorno cerrado del plano complejo que contiene todas las singularidades de F(z). Esta
técnica de inversión se utiliza raramente en casos prácticos, pero permite de poner en evidencia que existe
una relación única entre la señal discreta y su transformada en Z.
Se debe recordar sin embargo que la transformada en Z inversa solo proporciona información de la
señal temporal en los instantes de muestreo, por lo cual se requiere que las condiciones de Shannon hayan
sido respetadas en la discretización para poder obtener la señal temporal.
Tablas de transformadas.
Estas tablas son un repertorio de transformadas en Z y transformadas de Laplace de diversas
funciones, donde la transformada en Z se coloca en general bajo la forma de una función racional:
Num(z)
F(z) = Den(z)
Basta entonces con utilizar las tablas para obtener la inversa de la transformada de cada uno de los
términos, obteniéndose entonces la señal temporal completa.
División polinomial.
Esta técnica se basa en la definición de la transformada en Z de una señal, que es un polinomio de grado
infinito en z –1 . La idea es obtener el polinomio término a término a partir de la fracción racional F(z).
Se comienza por escribir a F(z) bajo la forma de una fracción racional en z–1 :
Se calcula a continuación el resultado de la división polinomial de N(z–1) entre D(z–1) siguiendo los
términos crecientes de z –1 . Obtenemos de esta manera a F(z) bajo la forma de un polinómio de orden
infinito:
Este método no permite determinar la expresión general ni obtener el enésimo término de la secuencia
temporal. Pero es útil para una programación utilizando una computadora o para la verificación de los
primeros términos obtenidos mediante algún otro método de cálculo.
Ejemplo 5.
Cálculo de la transformada inversa de:
e–1z + 1 — 2e–1 0.3679z + 0.2642
( )
G z = 2 =
z — (1 + e–1)z + e–1 1z2 — 1.3679z + 0.3679
Utilizando la descomposición en fracciones parciales, descomponemos:
G(z) e–1 + (1 — 2e–1)z–1 c1 c2
= = +
z (z — 1)(z — e–1) (z — 1) (z — e–1)
(z–1)(e—1+(1–2e—1)z—1)
(z–e—1)(e—1+(1–2e—1)z—1) e—1–1
c1 = (z–1)(z–e—1) = 1 ; c2 = | =
z=1 (z–1)(z–e—1) z=e—1 e—1
0 0.9145z–3 —0.2824z–4
Nota: la factorización en función del orden más elevado no es indispensable, puesto que la división con la
función original produce el mismo resultado.
Donde u(k) es conocido, y denotamos Y(z) = S(y(k)) y U(z) = S(u(k)), y además según el teorema del
retraso para condiciones iniciales nulas: S(y(k — i)) = z–iY(z) + ∑i–1 –l
l=0 z y(l — i)
Para determinar la secuencia de la solución y(k), basta con calcular la transformada en Z inversa de Y (z).
Obteniéndose una ecuación polinomial que puede resolverse manualmente.
Sistemas Digitales en Lazo Cerrado
La figura muestra la relación entre el modelo continuo y el discreto para el mismo sistema
0 10 20 30 40
60
1.5
0.5
Continuo
Discreto
0
Nót ema digital (40% en lugar de 16%).
ese
aqu un tiempo de estabilización dos
í veces mayor para el sistema digital.
que
:
Estos valores por supuesto
dependen del período de muestreo utilizado.
u
La
n
s
o
b
r
e
i
m
p
u
l
s
o
m
á
s
i
m
p
o
r
t
a
n
t
e
p
a
r
a
e
l
s
i
s
t
Tema 12. Introducción a los Sistemas Digitales 201
diferencia será más pequeña entre más pequeño sea este tiempo.
Donde Ad, Bd, Cd et Dd son las matrices discretas del sistema lineal, y Te es el período de muestreo. De igual
forma que para los sistemas continuos el modelo se representa mediante una serie de ecuaciones de primer
orden, en este caso en adelanto.
En el semiplano izquierdo del plano complejo σ < 0 (parte real del polo), por lo tanto la magnitud de z < 1,
y no existe influencia del ángulo de fase.
En consecuencia: Un sistema
digital será estable si los I I
polos de G(z) están dentro
del círculo unidad del plano r =1
complejo z, es decir si su R
magnitud es estrictamente
R
inferior a 1.
Ejemplo 7.
Para el sistema cuya función de transferencia continua es:
K
G(s) =
s(s + 1)
K(0.3679z + 0.2642)
T(z) = z2 + (0.3679K — 1.3679)z + 0.2642K + 0.3679
La estabilidad del sistema en lazo cerrado depende de las raíces de los polinomios:
s2 + s + K = 0
12
2 2
El sistema es estable si se cumple la condición: 0 < K < 2.39 .
A partir de la relación entre les polos para un sistema digital y continuo observamos que:
los polos reales estables en s (< 0) se convierten en polos reales positivos inferiores a 1.
los polos reales inestables se convierten en polos reales positivos superiores a 1.
los polos complejos conjugados estables se convierten en polos complejos con norma inferior a 1.
los polos complejos conjugados
inestables se convierten en
polos complejos conjugados de
norma superior a 1.
Jean-François DULHOSTE
los polos en el eje imaginario se convierten en polos complejos en el círculo unidad.
Los polos simples reales negativos en Z no tienen equivalente en s.
z
Polo simple real: Gi =
z–pi
Si pi > 1 el modo es divergente.
Si pi < —1 el modo es oscilatorio divergente.
Si —1 < pi < 0 el modo es oscilatorio amortiguado.
Si 0 < pi < 1 el modo es amortiguado.
Si pi = 1 el modo es permanente.
Si pi = —1 el modo es oscilatorio permanente.
z
Polo real múltiple: Gi = i
(z–p )di
Si pi > 1 el modo es divergente.
Si pi < —1 el modo es oscilatorio divergente.
Si —1 < pi < 0 el modo es oscilatorio amortiguado.
Si 0 < pi < 1 el modo es amortiguado.
Si pi = 1 el modo es divergente.
Si pi = —1 el modo es oscilatorio divergente.
Se nota que un polo doble en 1 conduce a un modo divergente mientras que un polo simple en 1 es
convergente.
Además, se constata que la convergencia de los modos hacia cero es más rápida mientras el módulo del polo
pi sea más cercano a cero. Se deduce entonces que los polos más rápidos del sistema digital son los polos
más cercanos a punto 0.
— Sistemas de orden 1
K K′z
G(s) = → G(z) =
s– p z – eTep
2
Con: p = —£ωn + iωn√1 — £2 y z1,2 = e–Te£ωne±jTeωn√1–£
Nota: Un sistema continuo sin ceros puede producir un sistema discreto con ceros. No existe una relación
directa entre los ceros del modelo continuo y los ceros del modelo digital.
La siguiente figura muestra las curvas de isoamortiguamiento (azul) y de isopulsación (rojo) para sistemas
numéricos de orden 2.
Jean-François DULHOSTE
Criterio algebraico de Jury
Este es un criterio que permite determinar si los polos del sistema en Z están o no en interior del
círculo unitario, sin necesidad de calcular el valor de estos polos. Es el equivalente a los criterios de Routh
o Hurwitz en sistemas continuos.
Para una función de transferencia en Z de la forma:
N(z)
G(z) =
D(z)
Con la ecuación característica: D(z) = a0 + a1z + a2z2 + ⋯ + anzn
Se construye la tabla del criterio de Jury siguiente:
a1
El criterio de Jury establece que la ecuación D(z) tiene sus polos en el interior del círculo unitario si y solo
si se cumplen todas las desigualdades siguientes:
1. |a0| — an < 0
2. D(1) > 0
3. (—1)nD(—1) > 0
4. |bj,0 | — |bj,n–j | > 0, ∀j ∈ {1, n — 2}
Nota:
La condición 1 implica que an > 0
La condición 4, que debe verificarse para sistemas de orden 3 o superior, se cumple si para cada
línea impar de la tabla, a excepción de la primera, el coeficiente de la primera columna tiene un valor
absoluto mayor que el de la última columna no nula.
|bj,0|— |bj,n–j| > 0
Para un sistema de orden 3, la cuarta condición puede simplificarse como sigue:
|a2 – a 2 |— |a0a2 – a1a3| > 0
0 3
Ejemplo 8.
Determinar la estabilidad del sistema representado por su ecuación característica discreta, en función de la
ganancia K > 0 del controlador, utilizando el criterio de Jury:
D(z) = z3 — 0.5z2 — 0.01z + 0.005 + K
Solución:
1. |a0| — an = |0.005 + K| — 1 < 0 ⇒ |0.005 + K| < 1 ⇒
—1.005 < K < 0.995
Método Directo:
Se determina inicialmente el modelo digital del sistema discretizando el modelo continuo y utilizando
la transformada en Z.
Se determina el controlador digital adaptado utilizando para ello técnicas de análisis y
concepción digitales.
Se realiza un análisis del comportamiento del sistema en lazo cerrado con las ecuaciones discretas.
Métodode transposición:
Se determina el modelo continuo del sistema.
Se determina un controlador en tiempo continuo para el sistema.
Se realiza una conversión del controlador en tiempo continuo en un controlador en tiempo
discreto por transposición (esto se puede realizar con la función c2d de Matlab).
Se analiza el comportamiento del sistema en lazo cerrado.
El método directo se utiliza menos en la práctica, puesto que requiere de herramientas digitales con
frecuencia poco dominadas por los ingenieros, pero los resultados obtenidos son más precisos que los que se
obtienen por transposición, debido a que no se realizan aproximaciones del comportamiento continuo de los
controladores digitales.
El método de Transposición
Una de las dificultades que presenta la síntesis de controladores por transposición es la selección del
método de transposición adecuado en cada caso. El objetivo de la transposición es el de determinar
un controlador digital para el cual el comportamiento del control digital se aproxime lo más posible al
comportamiento del control analógico.
Es imposible obtener un comportamiento idéntico debido al efecto del muestreo y de la cuantificación.
El control digital obtenido será en el mejor de los casos equivalente al control analógico pero en ningún
caso mejor.
Existen varias técnicas de transposición bien conocidas. Cada una de ellas presenta sus ventajas e
inconvenientes y por tanto se adaptan más o menos a cada problema en particular, ninguna de ellas puede
considerarse mejor que las otras de forma general.
Podemos citar entre los métodos existentes a los siguientes:
Transposición por muestreo — bloqueo de orden cero.
Transposición por aproximación bilineal.
Transposición por aproximación de Euler.
Transposición por conservación de polos y ceros.
Transposición por muestreo — bloqueo de orden uno (argumento ’foh’ de la función c2d de Matlab).
Con este método de transposición la salida del controlador digital es igual a la salida del controlador
analógico en los instantes de muestreo, y este valor se mantiene durante todo el período de muestreo. La
salida digital presentará siempre un retraso respecto a la señal analógica.
La siguiente figura muestra el efecto del período de muestreo para este método
En la práctica con esta técnica de transposición se obtienen buenos resultados solamente si el período de
muestreo es muy pequeño en relación al tiempo de respuesta de los sistemas controlados.
Si no se puede realizar un muestreo muy rápido, se recomienda el uso de alguna otra técnica de
transposición menos sensible al período de muestreo.
Sabiendo que para los sistemas continuos la integración es equivalente a una multiplicación por 1/s en el
dominio de la frecuencia, se quiere hacer una aproximación de 1/s con la función de transferencia digital:
Te (1 + z—1)
. Esto es equivalente a hacer el cambio de variable:
2 (1 – z—1)
2 (z — 1
s →) Te (z + 1
) 2+Tes
En la función de transferencia a transponer o lo que es equivalente imponer la igualdad z = .
2–Tes
Efecto en dominio de la frecuencia. A diferencia de la transposición por muestreo — bloqueo de orden cero, la
transformación bilineal no introduce retraso. El comportamiento de la fase será mejor en términos globales.
En cuanto a la amplitud, se observa que el filtrado introduce una distorsión muy importante para
frecuencias cercanas a la frecuencia de Nyquist.
Efecto sobre un control integral o derivativo. La amplificación de las frecuencias cercanas a la frecuencia de
Nyquist conduce a la amplificación de los ruidos a alta frecuencia, lo que no es para nada deseable, y menos
aún para una acción de control derivativa. Por el contrario el efecto sobre el integrador es beneficioso, debido
a que las frecuencias altas se atenúan.
Efecto del período de muestreo. Para el mismo período de muestreo la aproximación bilineal es mejor que
el bloqueo de orden cero.
Corrección de la distorsión. Para corregir la distorsión que introduce la aproximación bilineal, se puede añadir
un término de deformación previa “prewarping”, el cual permite corregir la distorsión para una frecuencia
específica ω0 . Para ello modificamos el cambio de variable como sigue:
ω0Te 2 z — 1
s→ 2
ω T Te z +
tan 0 e
1 2
El comportamiento del controlador digital se vuelve en este caso casi idéntico al del controlador analógico,
pero solo a la frecuencia específica ω0 .
Conclusión
La aproximación bilineal permite la transposición de forma eficaz de los controladores de tipo paso bajo.
Asociados a un prewarping, se vuelve también eficaz para filtros selectivos. Su buen comportamiento en
ángulo de fase la convierte en la primera opción para controladores de tipo adelanto o atraso.
Por el contrario no está adaptada para filtros con un banda pasante ancha debido a la amplificación
muy importante que produce a altas frecuencias. Se deberá entonces tener mucho cuidado en la
transposición de controladores de tipo PD o PID cuando los polos del controlador trabajan a latas
frecuencias.
A pesar de esto, como la mayoría de los controladores suelen ser de tipo paso bajo o pasa banda, con
frecuencia es la opción por defecto.
Período de muestreo
En la práctica el período de muestreo (Te) depende del tipo de proceso (químico, térmico, mecánico, etc.) y
debe seleccionarse pequeño en relación con el tiempo de respuesta (Tr) del sistema en lazo cerrado. Un
valor práctico es que se debe realizar al menos 4 a 10 muestras durante el tiempo de respuesta Nech =
Tr/Te = 4~10.
Retraso
El retraso en los sistemas de control digital debe siempre tomarse en cuenta, pues agrega una dificultad al
sistema. Este retraso puede tener varios orígenes a saber:
Muestreo.
Filtros y retrasos de comunicación.
Tiempo de ejecución de los códigos.
Interrupciones.
Dependencia de datos.
Preeminencia de procesos de orden superior.
El retraso introduce un desfasaje en las señales, el cual puede producir una inestabilidad en el sistema.
Según Bode el mayor retraso admisible puede calcularse como la siguiente expresión:
Da Mфπ [s]
= 180ω0
d
Donde:
Mф : Margen de fase
ω0 : Frecuencia de cruce con |G(jω0)| = 1
Bibliografía: 213
Bibliografía:
Katsuhito OGATA. Ingeniería de control moderna. Prentice Hall.
Benjamín Kuo. Sistemas de Control Automático. Pearson Education.
Smith y Corripio. Control Automático de Procesos. Limusa.
Distéfano, Stubberud y Williams. Retroalimentación y sistemas de Control. Serie Shaum. Mc Graw
Hill.
Richard C. Dorf & Robert H. Bishop. Modern Control Systems. Pearson, 11ava Edición, 2008.
Philippe de Laminart. Automatique, Commande de systèmes linéaires. Hermes 1993.
Florent Nageotte. Analyse et Commande des Systèmes Numériques. https://lsiit-
cnrs.unistra.fr/avr-fr/index.php/Page_personnelle_de_Florent_Nageotte.
Yves Brière. Commande des systèmes par calculateur. http://personnel.isae.fr/yves-briere/?lang=fr