SISTEMAS DE CONTROL AVANZADO
LUIS EDO GARCÍA JAIMES
POLITÉCNICO COLOMBIANO J.I.C
PRIMERA PARTE
Luis Edo García Jaimes
ANÁLISIS DE SISTEMAS DE CONTROL EN EL
ESPACIO DE ESTADO
Este método tiene como objetivo la descripción de
un sistema en función de 𝑛 ecuaciones en
diferencias o diferenciales de primer orden, las
cuales pueden combinarse para formar una
ecuación matricial en diferencias o una diferencial
de primer orden.
Luis Edo García Jaimes
CONCEPTOS BÁSICOS
• Estado: el estado de un sistema dinámico es el conjuntomás pequeño de variables, tales que el conocimiento dedichas variables en 𝑡 = 𝑡𝑜 junto con el conocimiento de laentrada para 𝑡 ≥ 𝑡𝑜 , determinan completamente elcomportamiento dinámico del sistema para 𝑡 ≥ 𝑡𝑜.
• Variables de Estado: son las variables que conforman elconjunto más pequeño de variables que determinan elestado de un sistema dinámico: 𝑥1, 𝑥2 . . . 𝑥𝑛.
• Vector de Estado: es el vector formado por el conjunto de
las 𝑛 variables de estado 𝑥1, 𝑥2 . . . 𝑥𝑛 que sean necesarias
para determinar completamente el comportamiento del
sistema.
Luis Edo García Jaimes
REPRESENTACIÓN EN EL ESPACIO DE ESTADO DE UN SISTEMA CONTINUO
Las variables 𝑢𝑖(𝑡), 𝑖 = 1. . . 𝑟, representan las entradas
Las variables 𝑦𝑖(𝑡), 𝑖 = 1. . . 𝑚, representan las salidas
Las variables 𝑥𝑖(𝑡), 𝑖 = 1. . . 𝑛, representan las variables de estado
𝑢 𝑡 =
)𝑢1(𝑡
𝑢2 𝑡∙∙
𝑢𝑟 𝑡
y 𝑡 =
)𝑦1(𝑡
)𝑦2(𝑡∙∙
)𝑦𝑚(𝑡
𝑥 𝑘𝑡 =
)𝑥1(𝑡
)𝑥2(𝑡∙∙
)𝑥n(𝑡
Luis Edo García Jaimes
y2(t)
y3(t)
ym(t)
Variables de
estadou1(t)
u2(t)
u3(t)
ur(t)
y1(t)
x1(t)
x2(t)
x3(t)
xn(t)
u(t) y(t)Vector de
Estado
a)
ECUACIÓN DE ESTADO
En general, la ecuación que describe el estado de un sistema de tiempo
continuo, se puede escribir en la forma:
𝑥 𝑡 = 𝑓 𝑥 𝑡 , 𝑢(𝑡)
La salida del sistema se puede dar como:
𝑦 𝑡 = 𝑔 𝑥 𝑡 , 𝑢(𝑡)
Para los sistemas lineales, de tiempo continuo e invariantes en el tiempo, la
ecuación de entrada y la ecuación de salida se pueden escribir así
𝑥 𝑡 = 𝐴𝑥 𝑡 + 𝐵𝑢 𝑡
𝑦 𝑡 = 𝐶𝑥 𝑡 + 𝐷𝑢 𝑡
Luis Edo García Jaimes
DEFINICIÓN DE TÉRMINOS
En donde:
𝑥(𝑡) = Vector de estado (vector 𝑛)
𝑦(𝑡) =Vector de salida (vector 𝑚)
𝑢(𝑡) = Vector de entrada (vector 𝑟)
𝐴 =Matriz de estado (matriz 𝑛 × 𝑛)
𝐵 = Matriz de entrada (matriz 𝑛 × 𝑟)
𝐶 = Matriz de salida (matriz 𝑚 × 𝑛)
𝐷 = Matriz de transmisión directa (matriz 𝑚 × 𝑟)
Luis Edo García Jaimes
DIAGRAMA DE BLOQUES SISTEMA CONTINUO
𝑥 𝑡 = 𝐴𝑥 𝑡 + 𝐵𝑢 𝑡
𝑦 𝑡 = 𝐶𝑥 𝑡 + 𝐷𝑢 𝑡
Luis Edo García Jaimes
++
++
B
A
D
Cu(t) y(t)
x(t)x(t)
REPRESENTACIÓN EN EL ESPACIO DE ESTADO DE UN SISTEMA DISCRETO
Las variables 𝑢𝑖(𝑘), 𝑖 = 1. . . 𝑟, representan las entradas
Las variables 𝑦𝑖(𝑘), 𝑖 = 1. . . 𝑚, representan las salidas
Las variables 𝑥𝑖(𝑘), 𝑖 = 1. . . 𝑛, representan las variables de estado
𝑢 𝑘 =
)𝑢1(𝑘
𝑢2 𝑘∙∙
𝑢𝑟 𝑘
𝑦 𝑘 =
)𝑦1(𝑘
)𝑦2(𝑘∙∙
)𝑦𝑚(𝑘
𝑥 𝑘 =
)𝑥1(𝑘
)𝑥2(𝑘∙∙
)𝑥n(𝑘
Luis Edo García Jaimes
y2(k)
y3(k)
ym(k)
Variables de
estadou1(k)
u2(k)
u3(k)
ur(k)
y1(k)
x1(k)
x2(k)
x3(k)
xn(k)
u(k) y(k)Vector de
Estado
b)
ECUACIÓN DE ESTADO SISTEMA DISCRETO
En general, la ecuación que describe el estado de un sistema de tiempo discreto,
en el instante 𝑘 + 1 se puede escribir en la forma:
𝑥 𝑘 + 1 = 𝑓 𝑥 𝑘 , 𝑢(𝑘)
Así mismo, la salida del sistema se puede dar como:
𝑦 𝑘 = 𝑔 𝑥 𝑘 , 𝑢(𝑘)
Para los sistemas lineales de tiempo discreto invariantes en el tiempo, la ecuación
de entrada y la ecuación de salida se pueden escribir así:
𝑥 𝑘 + 1 = 𝐴𝑥 𝑘 + 𝐵𝑢 𝑘
𝑦 𝑘 = 𝐶𝑥 𝑘 + 𝐷𝑢 𝑘
Luis Edo García Jaimes
DEFINICIÓN DE TÉRMINOS
En donde:
𝑥(𝑘) = Vector de estado (vector 𝑛)
𝑦(𝑘) =Vector de salida (vector 𝑚)
𝑢(𝑘) = Vector de entrada (vector 𝑟)
𝐴 =Matriz de estado (matriz 𝑛 × 𝑛)
𝐵 = Matriz de entrada (matriz 𝑛 × 𝑟)
𝐶 = Matriz de salida (matriz 𝑚 × 𝑛)
𝐷 = Matriz de transmisión directa (matriz 𝑚 × 𝑟)
Luis Edo García Jaimes
DIAGRAMA DE BLOQUES SISTEMA DISCRETO
𝑥 𝑘 + 1 = 𝐴𝑥 𝑘 + 𝐵𝑢 𝑘
𝑦 𝑘 = 𝐶𝑥 𝑘 + 𝐷𝑢 𝑘
Luis Edo García Jaimes
++
++
B
A
D
Cz-1
u(k)
x(k+1) x(k)
y(k)
FORMAS CANÓNICAS PARA ECUACIONES EN EL ESPACIO DE ESTADO EN TIEMPO DISCRETO
Sea el sistema en tiempo discreto definido por la ecuación de diferencias:
𝑦 𝑘 + 𝑎1𝑦 𝑘 − 1 + 𝑎2𝑦 𝑘 − 2 + ⋯𝑎𝑛𝑦 𝑘 − 𝑛 = 𝑏𝑜𝑢 𝑘 + 𝑏1𝑢 𝑘 − 1 ⋯+ 𝑏𝑛𝑢 𝑘 − 𝑛
En donde 𝑢(𝑘) es la entrada y 𝑦(𝑘) es la salida del sistema.
La función de transferencia de pulso correspondiente a la ecuación anterior está
dada por:
𝐺 𝑧 =𝑌(𝑧)
𝑈(𝑧)=
𝑏𝑜 + 𝑏1𝑧−1 + 𝑏2𝑧
−2 ⋯+ 𝑏𝑛𝑧−𝑛
1 + 𝑎1𝑧−1 + 𝑎2𝑧
−2 ⋯+ 𝑎𝑛𝑧−𝑛
𝐺 𝑧 =𝑌(𝑧)
𝑈(𝑧)=
𝑏𝑜𝑧𝑛 + 𝑏1𝑧𝑛−1 + 𝑏2𝑧
𝑛−2 ⋯+ 𝑏𝑛
𝑧𝑛 + 𝑎1𝑧𝑛−1 + 𝑎2𝑧
𝑛−2 ⋯+ 𝑎𝑛
Luis Edo García Jaimes
TIPOS DE FORMAS CANÓNICAS
Existen diferentes formas para representar el sistema discreto
definido por las ecuaciones dadas.
Las más utilizadas son las llamadas formas canónicas a saber:
• Forma canónica controlable (FCC).
• Forma canónica observable (FCO).
• Forma canónica diagonal (FCD).
• Forma canónica de Jordan (FCJ).
Luis Edo García Jaimes
FORMA CANÓNICA CONTROLABLE
La representación en el espacio de estado del sistema en tiempo discreto
definido por las funciones de transferencia dadas se puede expresar en
la forma:
𝑥1(𝑘 + 1)𝑥2(𝑘 + 1)𝑥3(𝑘 + 1)
∙∙
𝑥𝑛(𝑘 + 1)
=
−𝑎1 −𝑎2 ∙ −𝑎𝑛−1 −𝑎𝑛
1 0 ∙ 0 00 1 ∙ 0 0∙ ∙ ∙ ∙ ∙∙ ∙ ∙ ∙ ∙0 0 ∙ 1 0
𝑥1(𝑘)𝑥2(𝑘)𝑥3(𝑘)
∙∙
𝑥𝑛(𝑘)
+
10∙∙00
𝑢 𝑘
𝑦 𝑘 = 𝑏1 − 𝑎1𝑏𝑜 𝑏2 − 𝑎2𝑏𝑜 ∙ 𝑏𝑛 − 𝑎𝑛𝑏𝑜
𝑥1(𝑘)𝑥2(𝑘)
∙∙
𝑥𝑛−1(𝑘)𝑥𝑛(𝑘)
+ 𝑏𝑜𝑢(𝑘)
Luis Edo García Jaimes
FORMA CANÓNICA OBSERVABLE
La representación en el espacio de estado del sistema en tiempo discreto,
definido por las funciones de transferencia dadas, se puede expresar en la
forma:
𝑥1(𝑘 + 1)𝑥2(𝑘 + 1)
∙∙
𝑥𝑛−1(𝑘 + 1)𝑥𝑛(𝑘 + 1)
=
−𝑎1 1 0 ∙ 0−𝑎2 0 1 ∙ 0
∙ ∙ ∙ ∙ ∙∙ ∙ ∙ ∙ ∙
−𝑎𝑛−1 0 0 ∙ 1−𝑎𝑛 0 0 ∙ 0
𝑥1(𝑘)𝑥2(𝑘)
∙∙
𝑥𝑛−1(𝑘)𝑥𝑛(𝑘)
+
𝑏1 − 𝑎1𝑏𝑜
𝑏2 − 𝑎2𝑏𝑜
∙∙
𝑏𝑛−1 − 𝑎𝑛−1𝑏𝑜
𝑏𝑛 − 𝑎𝑛𝑏𝑜
𝑢 𝑘
𝑦 𝑘 = 1 0 0 ∙ 0
𝑥1(𝑘)𝑥2(𝑘)
∙∙
𝑥𝑛−1(𝑘)𝑥𝑛(𝑘)
+ 𝑏𝑜𝑢(𝑘
Luis Edo García Jaimes
FORMA CANÓNICA DIAGONAL
Si los polos de la función de transferencia de pulso dada son todos distintos, es decir, si ella
se puede expandir en fracciones parciales en la forma:
𝐺 𝑧 =𝐶1
𝑧 − 𝑝1+
𝐶2
𝑧 − 𝑝2+ ⋯
𝐶𝑛
𝑧 − 𝑝𝑛+ 𝑏𝑜
La representación en el espacio de estado definido se puede expresar en la forma:
𝑥1(𝑘 + 1)𝑥2(𝑘 + 1)
∙∙∙
𝑥𝑛(𝑘 + 1)
=
𝑝1 0 ∙ 0 00 𝑝2 ∙ 0 00 0 ∙ 0 0∙ ∙ ∙ ∙ ∙∙ ∙ ∙ ∙ ∙0 0 ∙ 0 𝑝𝑛
𝑥1(𝑘)𝑥2(𝑘)
∙∙∙
𝑥𝑛(𝑘)
+
11∙∙∙1
𝑢 𝑘
𝑦 𝑘 = 𝐶1 𝐶2 ∙ ∙ 𝐶𝑛
𝑥1 𝑘
𝑥2 𝑘∙∙∙
𝑥𝑛 𝑘
+ 𝑏𝑜𝑢 𝑘
Luis Edo García Jaimes
FORMA CANONICA DE JORDAN
Si al descomponer en fracciones parciales la función de transferencia dada se obtiene un polo
múltiple de orden 𝑟 en 𝑧 = 𝑝1 y todos los demás polos son distintos, es decir:
𝐺 𝑧 =𝐶1
(𝑧 − 𝑝1)𝑟+
𝐶2
(𝑧 − 𝑝1)𝑟−1
+ ⋯𝐶𝑟
𝑧 − 𝑝1+
𝐶𝑟+1
𝑧 − 𝑝𝑟+1+
𝐶𝑟+2
𝑧 − 𝑝𝑟+2+ ⋯
𝐶𝑛
𝑧 − 𝑝𝑛+ 𝑏0
La representación en el espacio de estado se puede expresar en la forma
𝑥1(𝑘 + 1)𝑥2(𝑘 + 1)
∙𝑥𝑟(𝑘 + 1)
𝑥𝑟+1(𝑘 + 1)∙∙
𝑥𝑛(𝑘 + 1)
=
𝑝1 1 0 00 𝑝1 1 0∙ ∙ ∙ ∙0 0 0 𝑝1
0 0 ∙ 0∙ ∙ ∙ ∙∙ ∙ ∙ ∙0 0 ∙ 0
0 ∙ 00 ∙ 0∙ ∙ 00 ∙ 0
𝑝𝑟+1 ∙ 00 ∙ 0∙ ∙ ∙0 ∙ 𝑝𝑛
𝑥1(𝑘)𝑥2(𝑘)
∙𝑥𝑟(𝑘)
𝑥𝑟+1(𝑘)∙∙
𝑥𝑛(𝑘)
+
00∙11∙∙1
𝑢
𝑦 𝑘 = 𝐶1 𝐶2 ∙ ∙ 𝐶𝑛
𝑥1 𝑘
𝑥2 𝑘∙∙∙
𝑥𝑛 𝑘
+ 𝑏𝑜𝑢 𝑘
Luis Edo García Jaimes
EJEMPLO
Obtenga tres representaciones en el espacio de estado para el sistema descrito
mediante la ecuación en diferencias dada. Considere condiciones iniciales
iguales a cero.
𝑦 𝑘 − 𝑦 𝑘 − 1 + 0.24𝑦 𝑘 − 2 = 0.4𝑢 𝑘 − 1 + 0.1𝑢(𝑘 − 2)
Tomando la transformada z:
1 − 𝑧−1 + 0.24𝑧−2 𝑌 𝑧 = 0.4𝑧−1 + 0.1𝑧−2 𝑈(𝑧)
La función de transferencia es:
𝐺 𝑧 =𝑌(𝑧)
𝑈(𝑧)=
0.4𝑧−1 + 0.1𝑧−2
1 − 𝑧−1 + 0.24𝑧−2 ∴𝑏0 = 0 𝑏1 = 0.4 𝑏2 = 0.1
𝑎1= −1 𝑎2 = 0.24
Luis Edo García Jaimes
SOLUCIÓN
a) FORMA CANÓNICA CONTROLABLE (FCC)
𝑥1(𝑘 + 1)𝑥2(𝑘 + 1)
=1 −0.241 0
𝑥1(𝑘)𝑥2(𝑘)
+10
𝑢 𝑘 𝑦 𝑘 = 0.4 0.1𝑥1(𝑘)𝑥2(𝑘)
b) FORMA CANÓNICA OBSERVABLE (FCO)
𝑥1(𝑘 + 1)𝑥2(𝑘 + 1)
=1 1
−0.24 0
𝑥1(𝑘)𝑥2(𝑘)
+0.40.1
𝑢 𝑘 𝑦 𝑘 = 1 0𝑥1(𝑘)𝑥2(𝑘)
c) FORMA CANÓNICA DIAGONAL (FCD)
𝐺 𝑧 =𝑌(𝑧)
𝑈(𝑧)=
0.4𝑧−1 + 0.1𝑧−2
1 − 𝑧−1 + 0.24𝑧−2=
1.7
𝑧 − 0.6−
1.3
𝑧 − 0.4
𝑥1(𝑘 + 1)𝑥2(𝑘 + 1)
=0.6 00 0.4
𝑥1(𝑘)𝑥2(𝑘)
+11
𝑢 𝑘 𝑦 𝑘 = 1.7 −1.3𝑥1(𝑘)𝑥2(𝑘)
Luis Edo García Jaimes
FUNCIÓN DE TRANSFERENCIA Y REPRESENTACIÓN EN EL ESPACIO DE ESTADO PARA UN SISTEMA CONTINUO
La representación en el espacio de estado de un sistema continuo es:
𝑥 𝑡 = 𝐴𝑥 𝑡 + 𝐵𝑢 𝑡𝑦 𝑡 = 𝐶𝑥 𝑡 + 𝐷𝑢 𝑡
Tomando la transformada de Laplace con CI=0 se obtiene𝑆𝑋 𝑆 = 𝐴𝑋 𝑆 + 𝐵𝑈(𝑆)
𝑆𝐼 − 𝐴 𝑋 𝑆 = 𝐵𝑈(𝑆)
Despejando 𝑋(𝑆) :𝑋 𝑆 = 𝑆𝐼 − 𝐴 −1𝐵𝑈(𝑆)
Es decir:𝑌 𝑆 = 𝐶 𝑆𝐼 − 𝐴 −1𝐵𝑈 𝑆 + 𝐷𝑈 𝑆
Como 𝑢(𝑡) es la entrada al sistema e 𝑦(𝑡) su salida, la función de transferencia es:
𝐺 𝑆 =𝑌(𝑆)
𝑈(𝑆)= 𝐶 𝑆𝐼 − 𝐴 −1𝐵 + 𝐷
Luis Edo García Jaimes
ECUACIÓN CARACTERÍSTICA SISTEMA CONTINUO
Por definición:
𝐴−1 =𝐴𝑑𝑗(𝐴)
𝑑𝑒𝑡(𝐴)=
𝐴𝑑𝑗(𝐴)
𝐴Para la expresión:
𝐺 𝑆 =𝑌(𝑆)
𝑈(𝑆)= 𝐶 𝑆𝐼 − 𝐴 −1𝐵
Se obtiene:
𝐺 𝑆 =𝑌(𝑆)
𝑈(𝑆)= 𝐶 𝑆𝐼 − 𝐴 −1𝐵 =
𝐶 ∗ 𝐴𝑑𝑗 𝑆𝐼 − 𝐴 ∗ 𝐵
𝑆𝐼 − 𝐴
La ecuación característica del sistema es:
𝑆𝐼 − 𝐴 = 𝑆𝑛 + 𝛼1𝑆𝑛−1 + 𝛼2𝑆
𝑛−2 + ⋯𝛼𝑛 = 0
Luis Edo García Jaimes
FUNCIÓN DE TRANSFERENCIA DE PULSO Y REPRESENTACIÓN EN EL ESPACIO DE ESTADO DISCRETO
La representación en el espacio de estado de un sistema discreto es:
𝑥 𝑘 + 1 = 𝐴𝑥 𝑘 + 𝐵𝑢 𝑘𝑦 𝑘 = 𝐶𝑥 𝑘 + 𝐷𝑢 𝑘
Tomando la transformada z con CI=0 se obtiene:
𝑧𝑋 𝑧 = 𝐴𝑋 𝑧 + 𝐵𝑈(𝑧)𝑧𝐼 − 𝐴 𝑋 𝑧 = 𝐵𝑈(𝑧)
Despejando 𝑋 𝑧 :
𝑋 𝑧 = 𝑧𝐼 − 𝐴 −1𝐵𝑈(𝑧)Es decir:
𝑌 𝑧 = 𝐶 𝑧𝐼 − 𝐴 −1𝐵𝑈 𝑧 + 𝐷𝑈 𝑧Como 𝑢(𝑘) es la entrada al sistema e 𝑦(𝑘) su salida, la función de
transferencia es:
𝐺 𝑧 =𝑌(𝑧)
𝑈(𝑧)= 𝐶 𝑧𝐼 − 𝐴 −1𝐵 + 𝐷
Luis Edo García Jaimes
ECUACIÓN CARACTERÍSTICA SISTEMA DISCRETO
Por definición:
𝐴−1 =𝐴𝑑𝑗(𝐴)
𝐷𝑒𝑡(𝐴)=
𝐴𝑑𝑗(𝐴)
𝐴Para la expresión:
𝐺 𝑧 =𝑌(𝑧)
𝑈(𝑧)= 𝐶 𝑧𝐼 − 𝐴 −1𝐵
Se obtiene:
𝐺 𝑧 =𝑌(𝑧)
𝑈(𝑧)= 𝐶 𝑧𝐼 − 𝐴 −1𝐵 =
𝐶 ∗ 𝐴𝑑𝑗 𝑧𝐼 − 𝐴 ∗ 𝐵
𝑧𝐼 − 𝐴
La ecuación característica del sistema es:
𝑧𝐼 − 𝐴 = 𝑧𝑛 + 𝛼1𝑧𝑛−1 + 𝛼2𝑧
𝑛−2 + ⋯𝛼𝑛 = 0
Luis Edo García Jaimes
EJEMPLOHallar la función de transferencia de pulso del sistema cuyo comportamiento
dinámico está descrito mediante la ecuación:
𝑥 𝑘 + 1 =0.6 0.20.5 0.3
𝑥 𝑘 +0.40.7
𝑢(𝑘) 𝑦 𝑘 = 0.5 −0.4 𝑥(𝑘)
La función de transferencia del sistema es:
𝐺 𝑧 =𝑌(𝑧)
𝑈(𝑧)= 𝐶 𝑧𝐼 − 𝐴 −1𝐵
𝑧𝐼 − 𝐴 =𝑧 00 𝑧
−0.6 0.20.5 0.3
=𝑧 − 0.6 −0.2−0.5 𝑧 − 0.3
𝑧𝐼 − 𝐴 −1 =𝑎𝑑𝑗(𝐴)
𝑧𝐼 − 𝐴=
𝑧 − 0.3 0.20.5 𝑧 − 0.6
𝑧 − 0.6 𝑧 − 0.3 − 0.1=
𝑧 − 0.3 0.20.5 𝑧 − 0.6
𝑧2 − 0.9𝑧 + 0.08
𝐺 𝑧 = 𝐶(𝑧𝐼 − 𝐴)−1𝐵 =0.5 −0.4
𝑧 − 0.3 0.20.5 𝑧 − 0.6
0.40.7
𝑧2 − 0.9𝑧 + 0.08
𝐺 𝑧 =−0.08𝑧 + 0.098
𝑧2 − 0.9𝑧 + 0.08
Luis Edo García Jaimes
DISEÑO DE SISTEMAS DE CONTROL EN EL ESPACIO DE ESTADO
• El diseño de un sistema de control digital consiste en determinar un
algoritmo que permita generar una secuencia de valores de las
variables de control de la planta 𝑢(𝑘), de manera que las salidas
𝑦(𝑘) cumplan con las especificaciones de funcionamiento
establecidas.
• En esta sección se presenta el diseño de controladores en el
espacio de estado, utilizando el método de asignación de polos.
Para su aplicación se requiere que el sistema sea completamente
controlable y completamente observable.
Luis Edo García Jaimes
CONTROLABILIDAD
Sea el sistema discreto:
𝑥 𝑘 + 1 = 𝐴𝑥 𝑘 + 𝐵𝑢 𝑘
𝑦 𝑘 = 𝐶𝑥 𝑘 + 𝐷𝑢 𝑘
Se dice que dicho sistema es de “estado completamente
controlable”, si es posible transferir el sistema desde un estado
inicial arbitrario a cualquier otro estado deseado en un intervalo
de tiempo finito mediante la aplicación de una señal de control,
no restringida, 𝑢(𝑘𝑇).
Luis Edo García Jaimes
CONDICIÓN DE CONTROLABILIDAD
El sistema descrito por la ecuación:
𝑥 𝑘 + 1 = 𝐴𝑥 𝑘 + 𝐵𝑢 𝑘
𝑦 𝑘 = 𝐶𝑥 𝑘 + 𝐷𝑢 𝑘
Es controlable si:
𝑅𝑎𝑛𝑔𝑜 𝐵 𝐴𝐵 𝐴2𝐵 ⋯ 𝐴𝑛−1𝐵 = 𝑛
Siendo 𝑛 × 𝑛 el orden de la matriz A.
Una condición suficiente y necesaria para la controlabilidad completa del
estado, es que no se presente cancelación de ceros y polos en la función
de transferencia del sistema.
Luis Edo García Jaimes
OBSERVABILIDAD
Sea el sistema discreto definido por:
𝑥 𝑘 + 1 = 𝐴𝑥 𝑘 + 𝐵𝑢 𝑘
𝑦 𝑘 = 𝐶𝑥 𝑘 + 𝐷𝑢 𝑘
Se dice que dicho sistema es de “estado complemente
observable” si cualquier estado inicial 𝑥(0) puede determinarse
a partir de la observación de 𝑦(𝑘) en 𝑛 períodos de muestreo
como máximo.
Luis Edo García Jaimes
CONDICIÓN DE OBSERVABILIDAD
El sistema discreto definido por la ecuación:𝑥 𝑘 + 1 = 𝐴𝑥 𝑘 + 𝐵𝑢 𝑘
𝑦 𝑘 = 𝐶𝑥 𝑘 + 𝐷𝑢 𝑘Es de estado completamente observable sí:
𝑅𝑎𝑛𝑔𝑜
𝐶𝐶𝐴𝐶𝐴2
⋮𝐶𝐴𝑛−1
= 𝑛
Una condición suficiente y necesaria para la observabilidad completa
del estado es que no se presente cancelación de ceros y polos en la
función de transferencia de pulso.
Luis Edo García Jaimes
EJEMPLO DE CONTROLABILIDAD Y OBSERVABILIDAD
Dado el sistema en tiempo discreto definido por:
𝑥 𝑘 + 1 =0 1 00 0.2 1
−0.6 −0.5 −0.3𝑥 𝑘 +
110
𝑢 𝑘 𝑦 𝑘 = 1 0 1 𝑥 𝑘
a) Es el sistema completamente controlable?
b) Es el sistema completamente observable?
Solución:
a) La matriz de controlabilidad es: 𝐶𝑜 = 𝐵 𝐴. 𝐵 𝐴2. 𝐵
𝐶𝑜 =1 1 0.21 0.2 −1.060 −1.1 −0.37
det 𝐶𝑜 = −1.09 𝑅𝑎𝑛𝑔𝑜 𝐶𝑜 = 3
𝐸𝑙 𝑠𝑖𝑠𝑡𝑒𝑚𝑎 𝑒𝑠 𝐶𝑜𝑛𝑡𝑟𝑜𝑙𝑎𝑏𝑙𝑒
Luis Edo García Jaimes
MATRIZ DE OBSERVBILIDAD
b) La matriz de observabilidad para el sistema dado es:
𝑂𝑏=𝐶𝐶𝐴𝐶𝐴2
𝑂𝑏 =1 0 1
−0.6 0.5 −0.30.18 −0.35 0.59
det 𝑂𝑏 = 0.31 𝑅𝑎𝑛𝑔𝑜 𝑂𝑏 = 3
𝐸𝑙 𝑠𝑖𝑠𝑡𝑒𝑚𝑎 𝑒𝑠 𝑜𝑏𝑠𝑒𝑟𝑣𝑎𝑏𝑙𝑒.
NOTA: Tener en cuenta que las matrices son de orden 3.
Luis Edo García Jaimes
CONTROL POR REALIMENTACIÓN DEL ESTADO Y ASIGNACIÓN DE POLOS
• El método de asignación de polos, comienza con ladeterminación de los polos de lazo cerrado deseados,utilizando especificaciones basadas en la respuestatransitoria y/o en los requerimientos de respuesta enfrecuencia.
• Si se desea ubicar los polos de lazo cerrado en 𝑧 = 𝑧1, 𝑧 =𝑧2, … 𝑧 = 𝑧𝑛 es posible elegir una matriz de ganancia derealimentación K adecuada, que force al sistema a tener lospolos de lazo cerrado en el lugar deseado siempre y cuandoel sistema sea de estado completamente controlable ycompletamente observable.
Luis Edo García Jaimes
CONTROL POR REALIMENTACIÓN DEL ESTADO
Sea el sistema de control en lazo abierto de la fig a, definido por laecuación de estado:
𝑥 𝑘 + 1 = 𝐴𝑥 𝑘 + 𝐵𝑢 𝑘
Si se elige como ley de control:
𝑢 𝑘 = −𝐾𝑥 𝑘
Se obtiene el sistema de control realimentado de la fig b. A esteesquema se le denomina “sistema con realimentación de estado”.
Sistema en lazo abierto Sistema en lazo cerrado
u(k)B
+
+
x(k+1) x(k)z-1I
A
-K
b.
Luis Edo García Jaimes
MATRIZ DE REALIMENTACIÓN KLa matriz 𝐾 = [𝑘1 𝑘2 ⋯
𝑘𝑛] se llama “matriz de ganancia de
realimentación” y convierte al sistema en un sistema de control en lazo
cerrado, con sus polos ubicados en el lugar deseado.Reemplazando la ley de control: 𝑢 𝑘 = −𝐾𝑥(𝑘):
𝑥 𝑘 + 1 = 𝐴𝑥 𝑘 − 𝐵𝐾𝑥(𝑘)
𝑥 𝑘 + 1 = 𝐴 − 𝐵𝐾 𝑥 𝑘
Tomando la transformada z a la ecuación se obtiene:
𝑧𝑋 𝑧 − 𝑧𝑥 0 = 𝐴 − 𝐵𝐾 𝑋(𝑧)
𝑧𝐼 − 𝐴 + 𝐵𝐾 𝑋 𝑧 = 𝑧𝑥 0
𝑋 𝑧 =𝑧. 𝑎𝑑𝑗 𝑧𝐼 – 𝐴 + 𝐵𝐾 . 𝑥(0)
𝑧𝐼 – 𝐴 + 𝐵𝐾Así, la ecuación característica del sistema en lazo cerrado es:
𝑧𝐼 – 𝐴 + 𝐵𝐾 = 𝑧𝑛 + 𝛼1𝑧𝑛−1 + 𝛼2𝑧
𝑛−2 ⋯+ 𝛼𝑛−1𝑧 + 𝛼𝑛 = 0
En donde 𝛼1, 𝛼2, ⋯𝛼𝑛 son los coeficientes de la ecuación característica
deseada.Luis Edo García Jaimes
CÁLCULO DE LA MATRIZ DE GANACIA DE REALIMENTACIÓN KLa matriz de ganancia de realimentación 𝐾 se puede obtener por
diferentes métodos. Uno de los más utilizados es el de la Formula de
Ackerman el cual permite calcular directamente la matriz de ganancia
de realimentación, a partir de la ecuación:
𝐾 = 0 0 ⋯ 1 𝐵 𝐴𝐵 𝐴2𝐵 ⋯ 𝐴𝑛−1𝐵 −1𝜙 𝐴
En donde:
𝜙 𝐴 = 𝐴𝑛+𝛼1𝐴𝑛−1 + 𝛼2𝐴
𝑛−2 ⋯+ 𝛼𝑛−1𝐴 + 𝛼𝑛𝐼
Siendo 𝛼1, 𝛼2 ⋯𝛼𝑛 los coeficientes de la ecuación característica
deseada:
𝑧 − 𝑧1 𝑧 − 𝑧2 . . 𝑧 − 𝑧𝑛 = 𝑧𝑛 + 𝛼1𝑧𝑛−1 + 𝛼2𝑧
𝑛−2 ⋯+ 𝛼𝑛−1𝑧 + 𝛼𝑛 = 0
Luis Edo García Jaimes
SISTEMA CON REALIMENTACIÓN DEL ESTADO Y ENTRADA DE REFERENCIA
El sistema de control anterior no tiene entrada de referencia. Este tipo de
control se denomina “sistema de control tipo regulador”. En la mayoría
de los casos, es necesario que la salida 𝑦(𝑘) siga a una entrada de
referencia 𝑟(𝑘) , este sistema se denomina “sistema de control tipo
Servo” y su configuración básica se muestra en la figura
u(k)B
++
x(k+1) x(k)z-1I
A
-K
++
r(k)Ko
v(k) y(k)C
Luis Edo García Jaimes
CÁLCULO DEL FACTOR DE CORRECCIÓN DE ERROR Ko
Luis Edo García Jaimes
Considerando el sistema de la figura anterior, se obtiene:
𝑥(𝑘 + 1) = 𝐴𝑥(𝑘) + 𝐵𝑢(𝑘)
𝑦(𝑘) = 𝐶𝑥(𝑘)
La señal de control 𝑢(𝑘) está dada por:
𝑢(𝑘) = 𝐾𝑜𝑟(𝑘) − 𝐾𝑥(𝑘)
La función de transferencia del sistema en lazo cerrado es:
𝐺𝑤(𝑧) = 𝑌(𝑧)
𝑅(𝑧)= 𝐶[𝑧𝐼 − 𝐴 + 𝐵𝐾]−1𝐵𝐾𝑜 = 𝐶
𝑎𝑑𝑗[𝑧𝐼 − 𝐴 + 𝐵𝐾]
𝑧𝐼 − 𝐴 + 𝐵𝐾 𝐵𝐾𝑜
Para tener error de estado estable igual a cero se debe cumplir que: 𝑦𝑠𝑠 = 𝑟
Aplicando el teorema del valor final resulta:
𝐾𝑜 lim𝑧→1
𝐶[𝑧𝐼 − 𝐴 + 𝐵𝐾]−1𝐵 = 1
EJEMPLO
Luis Edo García Jaimes
Obtener para el proceso a) La matriz de ganancia de realimentación 𝐾 de modo que
el sistema en lazo cerrado tenga constante de tiempo igual a 180 s y coeficiente de
amortiguamiento igual a 0.8 b) El coeficiente de corrección de error 𝐾0 para que el
error de estado estable sea igual a cero
La figura muestra el sistema de control de temperatura de una columna de
destilación. La dinámica del sistema se aproximó a un sistema de primer orden y
resultó ser:
𝐺𝑝(𝑆) =𝑇(𝑆)
𝐹(𝑆)=
0.5𝑒−60𝑆
250𝑆 + 1 𝑇 = 50 𝑠.
SOLUCIÓN
Luis Edo García Jaimes
La función de transferencia de pulso para el sistema está dada por:
𝐻𝐺(𝑧) = (1 − 𝑧−1)𝑧−𝑁ℑ𝑚 𝐺(𝑆)
𝑆 𝐺(𝑆) =
0.5𝑒−60𝑆
250𝑆 + 1
𝐻𝐺(𝑧) =0.0739𝑧 + 0.0167
𝑧3 − 0.8187𝑧2
La representación del sistema en FCC es:
𝑥(𝑘 + 1) = 0.8187 0 0
1 0 00 1 0
𝑥(𝑘) + 100 𝑢(𝑘)
𝑦(𝑘) = [0 0.0739 0.0167]
La ubicación de los polos de lazo cerrado se obtiene a partir de las especificaciones
constante de tiempo y del coeficiente de amortiguamiento requeridos.
𝜏 =1
𝜉. 𝜔𝑛 𝜔𝑛 =
1
𝜏. 𝜉=
1
180 ∗ 0.8 𝜔𝑛 = 0.00694 𝑟𝑎𝑑/𝑠
SOLUCIÓN
Luis Edo García Jaimes
La ubicación de los polos se obtiene a partir de las ecuaciones:
𝑧 = 𝑒−𝜉𝜔𝑛𝑇 = 𝑒−0.8∗0.00694∗50 = 0.7575 Para diseño 𝑧 < 1
𝜃 = 57.3𝜔𝑛𝑇 1 − 𝜉2 = 57.3 ∗ 0.00694 ∗ 50 ∗ 1 − 0.82 = 11.92° Para diseño 0° ≤ 𝜃 ≤ 80°
𝑧 = 𝑧 [𝑐𝑜𝑠𝜃 ± 𝑗𝑠𝑒𝑛𝜃]
𝑧 = 0.7575(𝑐𝑜𝑠11.92 ± 𝑗𝑠𝑒𝑛11.92) = 0.741 ± 𝑗0.156
Los polos de lazo cerrado quedan ubicados en:
𝑧 = 0.741 + 𝑗0.156 y 𝑧 = 0.741 + 𝑗0.156
El tercer polo se asigna en 𝑧 = 0.02 de modo que no sea polo dominante
La ecuación característica deseada es:
(𝑧 − 0.741 − 𝑗0.156)(𝑧 − 0.741 + 𝑗0.156)(𝑧 − 0.02) = 0
𝑧3 − 1.502𝑧2 + 0.603𝑧 − 0.01146 = 0
CÁLCULO DE LA MATRIZ K
Luis Edo García Jaimes
Utilizando la fórmula de Ackerman:
𝐾 = [0 0 0][𝐵 𝐴𝐵 𝐴2𝐵]−1𝜙(𝐴)
La ecuación característica deseada dio: 𝑧3 − 1.502𝑧2 + 0.603𝑧 − 0.01146 = 0
Por lo tanto:
𝜙(𝐴) = 𝐴3 − 1.502𝑧2 + 0.603𝐴 − 0.01146𝐼 = 0.02422 0 00.04358 −0.01146 0−0.6833 0.603 −0.01146
𝐶𝑂 = [𝐵 𝐴𝐵 𝐴2𝐵] = 1 0.8187 0.67020 1 0.81870 0 1
𝐾 = [0 0 1] 1 0.8187 0.67020 1 0.81870 0 1
−1
0.02422 0 00.04358 −0.01146 0−0.6833 0.603 −0.01146
𝐾 = [−0.6833 0.603 −0.01146]
CÁLCULO DEL FACTOR DE CORRECCIÓN DE ERROR Ko
Luis Edo García Jaimes
𝐾𝑜 lim𝑧→1
𝐶[𝑧𝐼 − 𝐴 + 𝐵𝐾]−1𝐵 = 1
𝐺𝑤(𝑧) = 𝐶[𝑧𝐼 − 𝐴 + 𝐵𝐾]−1𝐵 =0.0739𝑧 + 0.0167
𝑧3 − 1.502𝑧2 + 0.603𝑧 − 0.01146
𝐾𝑜 lim𝑧→1
0.0739𝑧 + 0.0167
𝑧3 − 1.502𝑧2 + 0.603𝑧 − 0.01146 = 1 𝐾𝑜 = 0.988
Sin factor de corrección Con factor de corrección Ko=0.988
0BSERVADORES O ESTIMADORES DE ESTADO
En la práctica, no todas las variables de estado de un sistema se pueden medir en
forma directa. Este hecho hace necesario estimar el valor de aquellas variables de
estado cuya medición directa no es posible.
El sistema que posibilita la estimación se denomina “Observador o estimador de
estado”.
El observador de un sistema dinámico lineal en tiempo discreto es otro sistema
dinámico lineal en tiempo discreto que tiene como entradas la entrada y la salida
del sistema discreto y como salida, los valores de las variables de estado.
Luis Edo García Jaimes
TIPOS DE OBSERVADORESPara resolver el problema de la observación son posibles dos soluciones:
a. Utilizar un observador tipo predictor que permite obtener el estado del
sistema en el instante (𝑘 + 1), estimando 𝑥(𝑘 + 1) a partir de la entrada 𝑢(𝑘) y
de la salida 𝑦(𝑘).
b. Utilizar un observador corriente que permite obtener el estado del sistema
en el instante (𝑘 + 1) estimando 𝑥(𝑘 + 1) a partir de la entrada 𝑢(𝑘) y de la
salida 𝑦(𝑘 + 1)
Las figuras representan, los dos tipos de observadores.
Observadores de estado a) Tipo Predictor b) Tipo Corriente Luis Edo García Jaimes
OBSERVADOR DE ESTADO TIPO PREDICTOR
Para obtener las ecuaciones que describen a este observador, se supone que el
estado real del sistema 𝑥(𝑘) no puede medirse directamente.
Si el estado 𝑥(𝑘) debe estimarse, es necesario que el estado estimado 𝑞(𝑘) y el
estado real 𝑥(𝑘) sean iguales.
La figura ilustra cómo se realiza la estimación de los estados.
Observador de estado tipo predictor
Luis Edo García Jaimes
ECUACIÓN DEL OBSERVADOR TIPO PREDICTOR
La planta está descrita mediante la ecuación:
𝑥(𝑘 + 1) = 𝐴𝑥(𝑘) + 𝐵𝑢(𝑘)
𝑦(𝑘) = 𝐶𝑥(𝑘)
De la figura del observador se deduce que el sistema tiene dos entradas 𝑢(𝑘) e
𝑦(𝑘), entonces, su ecuación se puede escribir en la forma:
𝑞(𝑘 + 1) = 𝐹𝑞(𝑘) + 𝐿𝑦(𝑘) + 𝐻𝑢(𝑘)
En donde 𝐹, 𝐿 y 𝐻 son matrices desconocidas.
Para que 𝑞(𝑘) = 𝑥(𝑘). Las matrices 𝐹, 𝐿 y 𝐻 deben cumplir que:
𝐻 = 𝐵 y 𝐴 = 𝐹 + 𝐿𝐶.
Entonces, la ecuación del observador predictor, se puede escribir en la forma:
𝑞(𝑘 + 1) = (𝐴 − 𝐿𝐶)𝑞(𝑘) + 𝐿𝑦(𝑘) + 𝐵𝑢(𝑘)
La matriz 𝐿 se denomina Matriz de ganancia de realimentación del observador.
La ecuación característica del observador es:
𝑧𝐼 − 𝐴 + 𝐿. 𝐶 = 0
Luis Edo García Jaimes
DIAGRAMA DE BLOQUES DEL SISTEMA CON EL OBSERVADOR TIPO PREDICTOR
La figura representa el sistema de control con la matriz de ganancia de
realimentación 𝐾 y el observador de estado incluidos.
Bx(k+1)
z-1x(k)
Cy(k)
++
A
-K
Cz-1
A
L
B+
++
+
+
-
u(k)
q(k+1) q(k) y(k)^
La ley de control es 𝑢(𝑘) = − 𝑘𝑞(𝑘) así la ecuación del observador tipo predictor
de orden completo se puede escribir en la forma:
𝑞(𝑘 + 1) = (𝐴 − 𝐿𝐶 − 𝐵𝐾)𝑞(𝑘) + 𝐿𝑦(𝑘) Luis Edo García Jaimes
FUNCIÓN DE TRANSFERENCIA DE PULSO DEL CONTROLADORUna vez obtenida la matriz de ganancia de realimentación 𝐾 y la matriz de ganancia
del observador 𝐿, es posible obtener la función de transferencia de pulso del
controlador. Para este controlador, la entrada es −𝑌(𝑧) y la salida 𝑈(𝑧).
La ecuación del observador de estado es:
𝑞(𝑘 + 1) = (𝐴 − 𝐿𝐶 − 𝐵𝐾)𝑞(𝑘) + 𝐿𝑦(𝑘)
Tomando la transformada z a esta ecuación del observador con CI=0
𝑧𝑄(𝑧) = [𝐴 − 𝐿𝐶 − 𝐵𝐾]𝑄(𝑧) + 𝐿𝑌(𝑧)
[𝑧𝐼 − 𝐴 + 𝐵𝐾 + 𝐿𝐶]𝑄(𝑧) = 𝐿𝑌(𝑧)
𝑄(𝑧) = [𝑧𝐼 − 𝐴 + 𝐵𝐾 + 𝐿𝐶]−1𝐿𝑌(𝑧)
La ley de control es: 𝑢(𝑘) = −𝐾𝑞(𝑘) 𝑈(𝑧) = −𝐾𝑄(𝑧)
Entonces: 𝑈(𝑧) = −𝐾𝑄(𝑧) = −𝐾[𝑧𝐼 − 𝐴 + 𝐵𝐾 + 𝐿𝐶]−1𝐿𝑌(𝑧)
𝐷(𝑧) = −𝑈(𝑧)
𝑌(𝑧)= 𝐾[𝑧𝐼 − 𝐴 + 𝐵𝐾 + 𝐿𝐶]−1𝐿
Esta ecuación permite estimar la función de transferencia de pulso del controlador. Luis Edo García Jaimes
FACTOR DE CORRECCIÓN DE ERROR
Si se desea tener error de estado estable igual a cero, ante una entrada en escalón,
es necesario adicionar un factor de corrección de error 𝐾𝑜 en el circuito del set-point
+-
R(z)Ko HG(z)
D(z)
Y(z)
Sistema de control por realimentación de estados con factor de corrección
de error en el circuito del set-point
De la figura se obtiene
𝐺𝑤(𝑧) =𝑌(𝑧)
𝑅(𝑧)=
𝐾0𝐻𝐺(𝑧)
1 + 𝐷(𝑧)𝐻𝐺(𝑧)
Para que el error sea cero debe cumplirse que 𝑦(∞) = 𝑟(𝑡), por lo tanto:
𝐾0 ∗ lim𝑧→1
𝐻𝐺(𝑧)
1 + 𝐷(𝑧)𝐻𝐺(𝑧)= 1
Luis Edo García Jaimes
EJEMPLO
Luis Edo García Jaimes
La figura representa el sistema de control de temperatura de una marmita. Se trata
de controlar la temperatura del producto regulando la entrada de vapor a la camisa
de la marmita. Una vez linealizado el sistema, su dinámica se puede aproximar a la
de un modelo de segundo orden así:
𝐺𝑝(𝑆) =𝑇𝑝(𝑠)
𝐹𝑣(𝑠)=
1.75
(100𝑆 + 1)(20𝑆 + 1) 𝑇 = 10𝑠
b) Diseñar un observador adecuado para el sistema. c) Obtener la ecuación del
controlador y el factor de corrección de error 𝐾𝑜
a) Determinar la matriz de ganancia de
realimentación K de modo que el sistema en lazo
cerrado tenga máximo sobreimpulso del 10% y tiempo
de establecimiento de 300 s.
SOLUCIÓNPrimero que todo se discretiza el sistema con T=10 s.
𝐻𝐺(𝑧) = (1 − 𝑧−1)ℑ 𝐺(𝑆)
𝑆 𝐺(𝑆) =
1.75
(100𝑆 + 1)(20𝑆 + 1)
𝐻𝐺(𝑧) =0.03602𝑧 + 0.0295
𝑧2 − 1.5114𝑧 + 0.5488
La representación del sistema en su forma canónica controlable es:
𝑥(𝑘 + 1) = 1.5114 −0.5488
1 0 𝑥(𝑘) +
10 𝑢(𝑘) 𝑦(𝑘) = [0.03602 0.0295]𝑥(𝑘)
a) Ubicación de los polos de lazo cerrado para calcular la matriz de ganancia de
realimentación 𝐾
𝑀𝑝 = 𝑒−
𝜋𝜉
1−𝜉2 𝜉 =
ln(𝑀𝑝)
𝜋2 + (ln(𝑀𝑝))2 𝜉 = 0.591
𝑡𝑠 =4
𝜉. 𝜔𝑛 𝜔𝑛 =
4
𝜉. 𝑡𝑠 𝜔𝑛 = 0.0225 𝑟𝑎𝑑/𝑠
Luis Edo García Jaimes
SOLUCIÓN: UBICACIÓN DE POLOS
Luis Edo García Jaimes
Para ubicar los polos deseados se utilizan las siguientes ecuaciones:
𝑧 = 𝑒−𝜉𝜔𝑛𝑇 = 𝑒−0.591∗0.0225∗10 = 0.875
𝜃 = 57.3𝜔𝑛𝑇 1 − 𝜉2 = 57.3 ∗ 0.0225 ∗ 10 ∗ 1 − 0.5912 = 10.4°
𝑧 = 𝑧 [𝑐𝑜𝑠𝜃 ± 𝑗𝑠𝑒𝑛𝜃] 𝑧 = 0.86 ± 𝑗0.158
La ecuación característica deseada debe ser de orden dos.
(𝑧 − 0.86 − 𝑗0.158)(𝑧 − 0.86 + 𝑗0.158) = 0 𝑧2 − 1.72𝑧 + 0.7645 = 0
La matriz 𝐾 se calcula con la fórmula de Ackerman:
𝐾 = [0 1][𝐵 𝐴𝐵]−1𝜙(𝐴)
𝜙(𝐴) = 𝐴2 − 1.72𝐴 + 0.7645 = −0.0995 0.1144−0.2086 0.2157
𝐶𝑜 = [𝐵 𝐴𝐵] = 1 1.51140 1
CÁLCULO DE LA MATRIZ K Y POLOS DEL OBSERVADOR
Luis Edo García Jaimes
𝐾 = [0 1] 1 1.51140 1
−1
−0.0995 0.1144−0.2086 0.2157
𝐾 = [−0.2086 0.2157]
Para el diseño del observador se recomienda que su coeficiente de
amortiguamiento sea igual al seleccionado para el cálculo de la matriz 𝐾 y que su
velocidad de respuesta sea mayor que la estimada para el sistema.
Se asume 𝜉 = 0.591 y 𝜔𝑛 = 0.04 𝑟𝑎𝑑/𝑠. (Mayor que la seleccionada para 𝐾)
𝑧 = 𝑒−𝜉𝜔𝑛𝑇 = 𝑒−0.591∗0.04∗10 = 0.789
𝜃 = 57.3𝜔𝑛𝑇 1 − 𝜉2 = 57.3 ∗ 0.04 ∗ 10 ∗ 1 − 0.5912 = 14.48°
Polos deseados: 𝑧 = 𝑧 [𝑐𝑜𝑠𝜃 ± 𝑗𝑠𝑒𝑛𝜃] 𝑧 = 0.763 ± 𝑗0.197
Ecuación característica para el observador: 𝑧2 − 1.526𝑧 + 0.6209 = 0
MATRIZ L Y ECUACIÓN DEL OBSERVADOR
Luis Edo García Jaimes
La matriz L del observador se calcula utilizando la fórmula de Ackerman:
𝐿 = 𝜙(𝐴) 𝐶𝐶𝐴
−1
01
𝜙(𝐴) = 𝐴2 − 1.526𝐴 + 0.6209𝐼 = 0.05 0.008
−0.0146 0.0721
𝑂𝑏 = 𝐶𝐶𝐴
= 0.036 0.0295
0.0839 −0.01976
𝐿 = 0.05 0.008
−0.0146 0.0721
0.036 0.02950.0839 −0.01976
−1
01 𝐿 =
0.3724−0.9496
Ecuación del observador:
𝑞(𝑘 + 1) = [𝐴 − 𝐿𝐶 − 𝐵𝐾]𝑞(𝑘) + 𝐿𝑦(𝑘)
𝑞(𝑘 + 1) = 1.7065 −0.77541.0342 0.0280
𝑞(𝑘) + 0.3724
−0.9496 𝑦(𝑘)
CÁLCULO DE LA FTP DEL CONTROLADOR
Luis Edo García Jaimes
La ecuación del controlador con el observador es:
𝐷(𝑧) = −𝑈(𝑧)
𝑌(𝑧)= 𝐾[𝑧𝐼 − 𝐴 + 𝐵𝐾 + 𝐿𝐶]−1𝐿
[𝑧𝐼 − 𝐴 + 𝐵𝐾 + 𝐿𝐶] = 𝑧 − 1.7065 0.7754−1.0342 𝑧 − 0.028
[𝑧𝐼 − 𝐴 + 𝐵𝐾 + 𝐿𝐶]−1 = 𝑧 − 0.028 −0.7754
1.0342 𝑧 − 1.7065
𝑧2 − 1.734𝑧 + 0.8497
𝐷(𝑧) = −𝑈(𝑧)
𝑌(𝑧)=
[−0.2086 0.2157] ∗ 𝑧 − 0.028 −0.7754
1.0342 𝑧 − 1.7065 ∗
0.3724−0.9496
𝑧2 − 1.734𝑧 + 0.8497
𝐷(𝑧) = −𝑈(𝑧)
𝑌(𝑧)=
−0.28279𝑧 + 0.28147
𝑧2 − 1.734𝑧 + 0.8497
CÁLCULO DEL FACTOR DE CORRECCIÓN DE ERROR Ko
Luis Edo García Jaimes
Para que el error de estado estable sea cero debe cumplirse que:
𝐾𝑜 lim𝑧→1
𝐺𝑤 (𝑧) = 1 𝐾0 lim𝑧→1
𝐻𝐺(𝑧)
1 + 𝐷(𝑧). 𝐻𝐺(𝑧)= 1
Con
𝐻𝐺(𝑧) =0.03602𝑧 + 0.0295
𝑧2 − 1.5114𝑧 + 0.5488 𝑦 𝐷(𝑧) = −
𝑈(𝑧)
𝑌(𝑧)=
−0.28279𝑧 + 0.28147
𝑧2 − 1.734𝑧 + 0.8497
Resulta:
𝐺𝑤(𝑧) =0.03602𝑧3 − 0.03299𝑧2 − 0.02055𝑧 + 0.02507
𝑧4 − 3.246𝑧3 + 4.0103𝑧2 − 2.2348𝑧 + 0.4747
Por lo tanto:
𝐾𝑜 lim𝑧→1
0.03602𝑧3 − 0.03299𝑧2 − 0.02055𝑧 + 0.02507
𝑧4 − 3.246𝑧3 + 4.0103𝑧2 − 2.2348𝑧 + 0.4747 = 1
1.7857𝐾𝑜 𝐾𝑜 = 0.56
DIAGRAMA DE BLOQUES PARA SIMULAR EL SISTEMA
Luis Edo García Jaimes
RESPUESTAS DEL SISTEMA CONTROLADO CON Y SIN FACTOR DE CORRECCIÓN DE ERROR
Respuesta sin factor de corrección Respuesta con factor de corrección
Luis Edo García Jaimes
SISTEMA TIPO SERVO CON INTEGRADORLa figura muestra un sistema de control por realimentación del estado en el cual se
utiliza un integrador adicional para estabilizar adecuadamente el sistema y mejorar
su exactitud
Bx(k+1)
z-1x(k)
Cy(k)
++
A
K1
Cz-1
A
L
B+
++
+
+-
u(k)
q(k+1) q(k) y(k)^
K1q(k)
+-
K i
z-1
++
+-
r(k) v(k)
Sistema tipo servo con integrador Luis Edo García Jaimes
DISEÑO DEL SISTEMA TIPO SERVO CON INTEGRADOR
La ecuación de estado de la planta y su correspondiente ecuación de salida son,
respectivamente:
𝑥(𝑘 + 1) = 𝐴𝑥(𝑘) + 𝐵𝑢(𝑘)
𝑦(𝑘) = 𝐶𝑥(𝑘)
La ley de control para el sistema es:
𝑢(𝑘) = −𝐾1𝑥(𝑘) + 𝐾𝑖𝑣(𝑘)
𝑣(𝑘) = 𝑟(𝑘) − 𝑦(𝑘) + 𝑣(𝑘 − 1)
Para realizar el diseño, utilizando la técnica de asignación de polos, se debe estimar
la matriz 𝐾𝑖 correspondiente al integrador y la matriz 𝐾1 correspondiente a la matriz
de ganancia de realimentación.
Luis Edo García Jaimes
ECUACIONES DE DISEÑO PARA EL SISTEMA TIPO SERVO CON INTEGRADOR
Para el cálculo de las matrices 𝐾1 y 𝐾𝑖 se utiliza la ecuación:
[𝐾1 ⋮ 𝐾𝑖] = 𝐾 + [0 ⋮ 𝐼𝑚 ] 𝐴 − 𝐼𝑛
𝐶𝐴
𝐵
𝐶𝐵 −1
𝐾 = [0 0 ⋯ 1] 𝐵 𝐴 𝐵 𝐴 2𝐵 ⋯ 𝐴 𝑛−1𝐵 −1
𝜙 𝐴
∅ 𝐴 = 𝐴 𝑛 + 𝛼1𝐴 𝑛−1 + 𝛼2𝐴
𝑛−2 + ⋯ 𝛼𝑛𝐼
𝛼1 , 𝛼2 …𝛼𝑛 : son los coeficientes de la ecuación característica deseada.
𝐴 = 𝐴
0 𝐵
0 (𝑛+𝑚)×(𝑛+𝑚)
𝐵 = 0
𝐼𝑚
(𝑛+𝑚)×𝑚
Conocidas las matrices 𝐾1 y 𝐾𝑖 , la ley de control para el sistema está dada por:
𝑈(𝑧) =[1 + 𝐾1[𝑧𝐼 − 𝐴 + 𝐿𝐶]−1𝐵]−1[𝐾𝑖𝑧[𝑅(𝑧) − 𝑌(𝑧)] − 𝐾1(𝑧 − 1)[𝑧𝐼 − 𝐴 + 𝐿𝐶]−1𝐿𝑌(𝑧)]
𝑧 − 1
La matriz 𝐿, correspondiente a la matriz de ganancia del observador, se calcula en
la misma como la del observador de orden completo con la fórmula de Ackerman. Luis Edo García Jaimes
EJEMPLO DISEÑO SISTEMA TIPO SERVO CON INTEGRADORSea el tanque con agitador representado en la figura.
El objetivo es controlar la temperatura 𝑇𝑜 del fluido de salida 𝑓𝑜 , manipulando el
caudal de vapor 𝑞𝑖 que pasa a través del serpentín. Mediante la aplicación de una
entrada en forma de escalón, se obtuvo la función de transferencia:
𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑒𝑙 𝑡𝑎𝑛𝑞𝑢𝑒 (%)
Flujo de vapor (%)=
𝑇0(𝑆)
𝑄𝑖(𝑆)=
2.5𝑒−20.3𝑆
75.4𝑆 + 1
Diseñe para el sistema un controlador tipo servo con integrador para regular la
temperatura del tanque. (Tiempo en s) Luis Edo García Jaimes
SOLUCIÓN: DISCRETIZACIÓN DEL SISTEMA
Selección del periodo de muestreo:
0.2(𝜏𝑒𝑞 + 𝜃′) ≤ 𝑇 ≤ 0.6(𝜏𝑒𝑞 + 𝜃′ )
𝜏𝑒𝑞 = Constante de tiempo equivalente del sistema en lazo cerrado sin el retardo.
𝜏𝑒𝑞 = 21.54𝑠. 8.36 ≤ 𝑇 ≤ 25.1 Se asume 𝑇 = 21 𝑠
𝐺𝑝(𝑧) = (1 − 𝑧−1)𝑧−𝑁ℑ𝑚 𝐺𝑝(𝑆)
𝑆 𝐺𝑝(𝑧) = (1 − 𝑧−1)ℑ𝑚
2.5
𝑆(75.4𝑆 + 1)
𝐺𝑝(𝑧) =0.0329𝑧 + 0.5747
𝑧(𝑧 − 0.7569)=
0.0329𝑧 + 0.5747
𝑧2 − 0.7569𝑧
Se representa el sistema en FCO:
𝑥(𝑘 + 1) = 0.7569 1
0 0 𝑥(𝑘) +
0.03290.5747
𝑢(𝑘) 𝑦(𝑘) = [1 0] 𝑥(𝑘)
Luis Edo García Jaimes
CÁLCULO DE LA MATRIZ
La matriz 𝐾 , que permite calcular a la matriz de realimentación 𝐾1 y a la matriz del
integrador 𝐾𝑖 se calcula a partir de la ecuación:
𝐾 = [0 0 ⋯ 1][𝐵 𝐴 𝐵 𝐴 2𝐵 ⋯ 𝐴𝑛−1𝐵]−1𝜙(𝐴 )
𝐴 = 𝐴
0 𝐵
0 =
0.7569 10 0
0 0
0.03290.5747
0 𝐵 =
0𝐼𝑚
= 001
[𝐵 𝐴 𝐵 𝐴 2𝐵 ] = 0 0.0329 0.59960 0.5147 01 0 0
𝜙 𝐴 = 𝐴 3 + 𝛼1𝐴 2 + 𝛼2𝐴 + 𝛼3𝐼
La ubicación de los polos para la matriz 𝐾 , se obtiene con las ecuaciones:
𝑧 = 𝑒−𝜉𝑤𝑛𝑇 𝜃 = 57.3𝑤𝑛𝑇 1 − 𝜉2 𝑧 = 𝑧 [𝑐𝑜𝑠𝜃 ± 𝑗𝑠𝑒𝑛𝜃]
Estos polos permiten el obtener la ecuación característica y por lo tanto el de los
coeficientes 𝛼1, 𝛼2, 𝛼3 Luis Edo García Jaimes
UBICACIÓN DE LOS POLOS PARA LA MATRIZ K
El tiempo de establecimiento del sistema en lazo abierto es: 𝑡𝑠𝐿𝐴 = 4𝜏 = 301.6 𝑠
Se selecciona el tiempo de establecimiento: 0.3𝑡𝑠𝐿𝐴 ≤ 𝑡𝑠 ≤ 0.8𝑡𝑠𝐿𝐴 𝑡𝑠 = 225 𝑠
Se selecciona el coeficiente de amortiguamiento: 0.6 ≤ 𝜉 ≤ 0.9 𝜉 = 0.8
La ubicación de los polos es:
𝑧 = 𝑒−𝜉𝑤𝑛𝑇 𝜃 = 57.3𝑤𝑛𝑇 1 − 𝜉2 𝑧 = 𝑧 [𝑐𝑜𝑠𝜃 ± 𝑗𝑠𝑒𝑛𝜃]
𝑡𝑠 =4
𝜉𝑤𝑛 𝑤𝑛 =
4
𝜉𝑡𝑠=
4
0.8 ∗ 225 𝑤𝑛 = 0.0222 𝑟𝑎𝑑/𝑠
𝑧 = 𝑒−0.8∗0.0222∗21 = 0.688 𝜃 = 57.3 ∗ 0.0222 ∗ 21 ∗ 1 − 0.82 = 16.02°
Polos están ubicados en 𝑧 = 0.661 ± 𝑗0.1898. Se adiciona un tercer polo en 𝑧 = 0
La ecuación característica es:
(𝑧 − 0.661 − 𝑗0.1898)(𝑧 − 0.661 + 𝑗0.1898)𝑧 = 𝑧3 − 1.322𝑧2 + 0.4729𝑧 = 0
Luis Edo García Jaimes
CÁLCULO DE LAS MATRICES K1 Y Ki
Con la ecuación característica obtenida resulta:
𝜙 𝐴 = 𝐴 3 − 1.322𝐴 2 + 0.4729𝐴 = 0.0342 0.0451 −0.3232
0 0 0.27370 0 0
𝐾 = [0 0 1] 0 0.0329 0.59960 0.5147 01 0 0
−1
0.0342 0.0451 −0.3232
0 0 0.27370 0 0
𝐾 = [0.05702 0.0753 −0.5651]
[𝐾1 𝐾𝑖] = 𝐾 + [0 ⋮ 𝐼𝑚 ] 𝐴 − 𝐼𝑛 𝐵𝐶. 𝐴 𝐶. 𝐵
−1
[𝐾1 𝐾𝑖] = [0.05702 0.0753 −0.5651] + [0 0 1] −0.2431 1 0.0329
0 −1 0.57470.7569 1 0.0329
−1
[𝐾1 𝐾𝑖] = [0.5386 0.7117 0.2483]
𝐾1 = [0.5386 0.7117] 𝐾𝑖 = [0.2483]
Luis Edo García Jaimes
CÁLCULO DEL OBSERVADOREl diseño del observador se realiza utilizando la fórmula de Ackerman:
𝐿 = 𝜙(𝐴) 𝐶𝐶𝐴
−1
01
Ubicación de los polos para el observador: se toma 𝑡𝑠𝑜 < 𝑡𝑠𝑘 . En este caso se asume
𝑡𝑠𝑜 = 200 𝑠 𝑦 𝜉 = 0.8
𝑤𝑛 =4
𝜉. 𝑡𝑠=
4
0.8 ∗ 200= 0.025 𝑟𝑎𝑑/𝑠
𝑧 = 𝑒−0.8∗0.025∗21 = 0.657 𝜃 = 57.3 ∗ 0.025 ∗ 21 ∗ 1 − 0.82 = 18.04°
𝑧 = 𝑧 [𝑐𝑜𝑠𝜃 ± 𝑗𝑠𝑒𝑛𝜃] = 0.624 ± 𝑗0.203
Ecuación característica: 𝑧2 − 1.248𝑧 + 0.4305 = 0
𝜙(𝐴) = 𝐴2 − 1.248𝐴 + 0.4305𝐼 = 0.0587 −0.4917
0 0.4305
𝐶𝐶𝐴
= 0.0587 −0.4917
0 0.4305
𝐿 = 0.0587 −0.4917
0 0.4305
0.0587 −0.49170 0.4305
−1
01 𝐿 =
−0.49110.4305
Luis Edo García Jaimes
CÁLCULO DE LA LEY DE CONTROL
La ley de control es:
𝑈(𝑧) =[1 + 𝐾1[𝑧𝐼 − 𝐴 + 𝐿𝐶]−1𝐵]−1[𝐾𝑖𝑧[𝑅(𝑧) − 𝑌(𝑧)] − 𝐾1(𝑧 − 1)[𝑧𝐼 − 𝐴 + 𝐿𝐶]−1𝐿𝑌(𝑧)]
𝑧 − 1
𝐾1[𝑧𝐼 − 𝐴 + 𝐿𝐶]−1𝐵 =0.4254𝑧 − 0.2103
𝑧2 − 1.248𝑧 + 0.4305
𝐾1 [𝑧𝐼 − 𝐴 + 𝐿𝐶]−1𝐿 =0.04171𝑧
𝑧2 − 1.248𝑧 + 0.4305
𝑈(𝑧) =[0.248𝑧3 − 0.3098𝑧2 + 0.10689𝑧]𝑅(𝑧) − [0.248𝑧3 − 0.2681𝑧2 + 0.0652𝑧]𝑌(𝑧)
𝑧3 − 1.8226𝑧2 + 1.0428𝑧 − 0.2202
Ecuación en diferencias correspondiente a la ley de control:
𝑢(𝑘) = 0.248𝑟(𝑘) − 0.3098𝑟(𝑘 − 1) + 0.10689𝑟(𝑘 − 2) − 0.248𝑦(𝑘) + 0.2681𝑦(𝑘 − 1)
− 0.0652𝑦(𝑘 − 2) + 1.8226𝑢(𝑘 − 1) − 1.0428𝑢(𝑘 − 2) + 0.2202𝑢(𝑘 − 3) Luis Edo García Jaimes
RESPUESTA DEL SISTEMA TIPO SERVO CON INTEGRADOR
Luis Edo García Jaimes
SISTEMAS NO LINEALES
Los sistemas no lineales representan sistemas cuyo
comportamiento no se puede expresar como la suma de los
comportamientos de sus descriptores, es decir son sistemas que
no cumplen el principio de superposición.
En los sistemas no lineales, las ecuaciones de movimiento,
evolución o comportamiento que regulan su comportamiento
dinámico son no lineales.
Luis Edo García Jaimes
EJEMPLO DE SISTEMAS NO LINEALES (2)
Luis Edo García Jaimes
REPRESENTACIÓN DE UN SISTEMA NO LINEAL
Un sistema no lineal se puede representar mediante ecuaciones de estado en la
siguiente forma:
𝑥 1 = 𝑓1(𝑥, 𝑢, 𝑡)
𝑥 2 = 𝑓2(𝑥, 𝑢, 𝑡)
⋯ ⋯
𝑥 𝑛 = 𝑓𝑛(𝑥, 𝑢, 𝑡)
𝑦 = ℎ(𝑥, 𝑢, 𝑡)
Estas ecuaciones se pueden escribir en la forma matricial así:
𝒙 (𝑡) = 𝒇[𝑥(𝑡), 𝑢(𝑡)]
𝑦(𝑡) = 𝒉[𝒙(𝑡), 𝒖(𝑡)]
En donde 𝒙(𝑡) es el vector de estado (𝑛 × 1), 𝒖(𝑡) es el vector de entradas (𝑟 × 1)
y 𝒇[𝑥(𝑡), 𝑢(𝑡)] es un vector que es función del vector de estado y del vector de
entrada. Luis Edo García Jaimes
LINEALIZACIÓN DE UN SISTEMA NO LINEALPara linealizar un sistema no lineal existen diferentes métodos: uno de ellos
consiste en la expansión de las ecuaciones de estado no lineales en series de Taylor
alrededor de un punto o de una trayectoria de operación nominal del sistema,
despreciando los términos de orden superior al primero, con lo cual resulta una
aproximación lineal de las ecuaciones de estado en un punto determinado.
𝑥 (𝑡) = 𝐴 𝑥(𝑡) + 𝐵 𝑢(𝑡)
En donde:
𝐴 =
𝜕𝑓1
𝜕𝑥1
𝜕𝑓1
𝜕𝑥2⋮
𝜕𝑓1
𝜕𝑥𝑛
𝜕𝑓2
𝜕𝑥1
𝜕𝑓2
𝜕𝑥2⋮
𝜕𝑓2
𝜕𝑥𝑛
⋯ ⋯ ⋮ ⋯𝜕𝑓𝑛𝜕𝑥1
𝜕𝑓𝑛𝜕𝑥2
⋮𝜕𝑓𝑛𝜕𝑥𝑛
𝑃𝑜
𝐵 =
𝜕𝑓1
𝜕𝑢1
𝜕𝑓1
𝜕𝑢2⋮
𝜕𝑓1
𝜕𝑢𝑛
𝜕𝑓2
𝜕𝑢1
𝜕𝑓2
𝜕𝑢2⋮
𝜕𝑓2
𝜕𝑢𝑛
⋯ ⋯ ⋮ ⋯𝜕𝑓𝑛𝜕𝑢1
𝜕𝑓𝑛𝜕𝑢2
⋮𝜕𝑓𝑛𝜕𝑢𝑛
𝑃𝑜
𝑃𝑜: corresponde al punto de equilibrio alrededor del cual se va a linealizar el
sistema. Los valores de 𝑥(𝑡) y de 𝑢(𝑡) deben mantenerse siempre lo más cerca
posible a los valores de referencia 𝑥𝑜 y 𝑢𝑜 respectivamente Luis Edo García Jaimes
PROPIEDADES DEL SISTEMA LINEALIZADO
Luis Edo García Jaimes
La línea continua representa la función no lineal 𝑓 (𝑥, 𝑢). El cırculo negro es el punto
[𝑥𝑜 , 𝑢𝑜] alrededor del cual se realiza la linealización de la función no lineal. Como
se observa, linealización solo es válida en el interior de una región, demarcada por
el círculo externo. En términos generales no se puede asegurar decir de qué
tamaño es la región donde es válida la linealización; sólo se puede decir es que es
pequeña. La línea discontinua representa la función linealizada.
EJEMPLO DE LINEALIZACIÓN
Luis Edo García Jaimes
𝐹 [𝑙 𝑚𝑖𝑛 ] 𝑉 [𝑙] 𝐶𝐴𝑂 [𝑚𝑜𝑙 𝑙 ] 𝐶𝐵𝑂 [𝑚𝑜𝑙 𝑙 ] 𝐾1[𝑙 𝑚𝑜𝑙. 𝑚𝑖𝑛 ] 𝐾2 [𝑙 𝑚𝑖𝑛 ]
100 1000 200 0 0.01 0.5
Reactivos
FIC
TIC
LIC
Producto
T
Tj
CWS
CWR
F
CAo
CA
FJ
V
To
Tjo
CBO
CB
C
Se tiene un reactor continuamente agitado en el cual ocurre el siguiente sistema de
reacciones en serie A → B → C. El modelo matemático del reactor, que describe la
variación de las concentraciones de los componentes A y B, está dado por:
𝑑𝐶𝐴
𝑑𝑡=
𝐹
𝑉(𝐶𝐴𝑂 − 𝐶𝐴) − 𝐾1𝐶𝐴
2 1.
𝑑𝐶𝐵
𝑑𝑡=
𝐹
𝑉(𝐶𝐵𝑂 − 𝐶𝐵) + 𝐾1𝐶𝐴
2 − 𝐾2𝐶𝐵 2.
Linealice el sistema alrededor de la
zona de trabajo dada y diseñe un
controlador tipo servo con integrador de
modo que el sistema tenga polos en
𝑧 = 0.8, 𝑧 = 0.6, 𝑧 = 0 y los polos del
observador en 𝑧 = 0.6, 𝑧 = 0.5
La zona de trabajo del reactor está definida por los siguientes parámetros:
LINEALIZACIÓN DEL MODELOLa linealización se inicia con el cálculo de los puntos de equilibrio del sistema en la
zona de trabajo especificada 𝐹 = 100 𝑙 𝑚𝑖𝑛, 𝐶𝐴𝑂 = 200 𝑚𝑜𝑙 𝑙, 𝐶𝐵𝑂 = 0
Los puntos de equilibrio se calculan haciendo cero las derivadas y resolviendo las
ecuaciones resultantes con respecto a las variables de entrada y de salida.
Se consideran como salidas del sistema las concentraciones 𝐶𝐴 𝑦 𝐶𝐵 y como
entrada el flujo 𝐹, por lo tanto:
20 − 0.1𝐶𝐴 − 0.01𝐶𝐴2 = 0 𝐶𝐴 = 40 𝑚𝑜𝑙 𝑙
−0.1𝐶𝐵 + 0.01𝐶𝐴2 − 0.5𝐶𝐵 = 0 𝐶𝐵 = 26.66 𝑚𝑜𝑙 𝑙
Punto de equilibrio: 𝑃0 = [𝐶𝐴 𝐶𝐵 𝐹] = [40 26.66 100]
El sistema linealizado se puede escribir en la forma:
𝑥 = 𝐴𝑥 + 𝐵𝑢 𝑦 = 𝐶𝑥 𝐴 =
𝜕𝑓1
𝜕𝐶𝐴
𝜕𝑓1
𝜕𝐶𝐵
𝜕𝑓2
𝜕𝐶𝐴
𝜕𝑓2
𝜕𝐶𝐵
𝑃0
𝐵 =
𝜕𝑓1
𝜕𝐹𝜕𝑓2
𝜕𝐹
𝑃0
Luis Edo García Jaimes
MODELO LINEALIZADO
Luis Edo García Jaimes
El cálculo de las derivadas parciales con respecto a las variables de entrada y de
salida evaluado en el punto de equilibrio establecido da:
𝜕𝑓1
𝜕𝐶𝐴= −0.001𝐹 − 2𝐾1𝐶𝐴 = −0.9
𝜕𝑓1
𝜕𝐶𝐵= 0
𝜕𝑓1
𝜕𝐹= 0.001𝐶𝐴0 − 0.001𝐶𝐴 = 0.16
𝜕𝑓2
𝜕𝐶𝐴= 2𝐾1𝐶𝐴 = 0.8
𝜕𝑓2
𝜕𝐶𝐵= −0.001𝐹 − 𝐾2 = −0.6
𝜕𝑓2
𝜕𝐹= 0.001𝐶𝐵0 − 0.001𝐶𝐵 = −0.0266
Por lo tanto, la representación del sistema linealizado en el espacio de estado es:
𝐶
𝐴
𝐶 𝐵
= −0.9 00.8 −0.6
𝐶𝐴
𝐶𝐵 +
0.16−0.0266
𝐹 𝐶𝐴 = [1 0] 𝐶𝐴
𝐶𝐵
La función de transferencia del sistema continuo es:
𝐺𝑃(𝑆) = 𝐶[𝑆𝐼 − 𝐴]−1𝐵
𝐺𝑃(𝑆) = [1 0] 𝑆 + 0.9 0−0.8 𝑆 + 0.6
−1
0.16
−0.0266 𝐺𝑃(𝑆) =
𝐶𝐴(𝑆)
𝐹(𝑆)=
0.16𝑆 − 0.0477
(𝑆 + 0.9)(𝑆 + 0.6)
La constante de tiempo del sistema es: 𝜏 = 2.77 𝑚𝑖𝑛. Para discretizar el proceso se toma
como periodo de muestreo 𝑇 = 0.5 𝑚𝑖𝑛
ECUACIÓN DE ESTADO DISCRETA Y CÁLCULO DE LA MATRIZ 𝐾
Luis Edo García Jaimes
𝐻𝐺(𝑧) = (1 − 𝑧−1)ℑ 𝐺(𝑆)
𝑆 𝐻𝐺(𝑧) = (1 − 𝑧−1)ℑ
0.16𝑆 + 0.096
𝑆(𝑆 + 0.9)(𝑆 + 0.6)
𝐻𝐺(𝑧) =𝐶𝐴(𝑧)
𝐹(𝑧)=
0.0644𝑧 − 0.0477
𝑧2 − 1.3784𝑧 + 0.4723
La ecuación de estado en FCO es:
𝐶𝐴(𝑘 + 1)
𝐶𝐵(𝑘 + 1) =
1.3784 1−0.4723 0
𝐶𝐴(𝑘)
𝐶𝐵(𝑘) +
0.0644−0.0477
𝐹(𝑘) 𝑦(𝑘) = [1 0] 𝐶𝐴(𝑘)
𝐶𝐵(𝑘)
La ecuación característica deseada para el sistema en lazo cerrado es:
(𝑧 − 0.8)(𝑧 − 0.6)𝑧 = 0 𝑧3 − 1.4𝑧2 + 0.48𝑧 = 0
𝐾 = [0 0 1][𝐵 𝐴 𝐵 𝐴 2𝐵 ]−1𝜙(𝐴 )
𝜙 𝐴 = 𝐴 3 − 1.4𝐴 2 + 0.48𝐴 = −0.0202 −0.02207 −0.000390.01042 0.01020 0.00028
0 0 0
𝐶𝑂 = [𝐵 𝐴 𝐵 𝐴 2𝐵 ] = 0 0.0644 0.041060 −0.0477 −0.030411 0 0
CÁLCULO DE LA MATRIZ [K1 Ki]
𝐾 = [0 0 1] 0 0.0644 0.041060 −0.0477 −0.030411 0 0
−1
−0.0202 −0.02207 −0.000390.01042 0.01020 0.00028
0 0 0
𝐾 = [−1533.55 −2069.99 −0.0000216]
[𝐾1 𝐾𝑖] = 𝐾 + [0 ∶ 𝐼𝑚 ] 𝐴 − 𝐼𝑛 𝐵
𝐶𝐴 𝐶𝐵 −1
[𝐾1 𝐾𝑖] = [−1533.55 −2069.99 −0.0000216] + [0 0 1] 0.3784 1 0.00644
−0.4723 −1 −0.04771.3784 1 0.0644
−1
[𝐾1 𝐾𝑖] = [5966.28 8041.06 4.79] [𝐾1] = [5966.28 8041.06] [𝐾𝑖] = [4.79]
Cálculo de la matriz 𝐿 del observador:
𝐿 = 𝜙(𝐴) 𝐶𝐶𝐴
−1
01
Polos deseados para el observador: 𝑧 = 0.6 , 𝑧 = 0.5
Ecuación característica para el observador: 𝑧2 − 1.1𝑧 + 0.3 = 0 Luis Edo García Jaimes
CALCULO DEL OBSERVADOR Y DE LA LEY DE CONTROL
𝜙(𝐴) = 𝐴2 − 1.1𝐴 + 0.3𝐼 = 0.2114 0.2784
−0.1314 −0.1723
𝐶𝐶𝐴
= 1 0
1.3784 1
𝐿 = 0.2114 0.2784
−0.1314 −0.1723
1 01.3784 1
−1
01 𝐿 =
0.2787−0.1723
La ecuación del observador es:
𝑞(𝑘 + 1) = [𝐴 − 𝐿𝐶]𝑞(𝑘) + 𝐿𝑦(𝑘) + 𝐵𝑢(𝑘)
𝑞(𝑘 + 1) = 1.1 1
−0.3 0 𝑞(𝑘) +
0.2787−0.1723
𝑦(𝑘) + 0.0644
−0.0477 𝑢(𝑘)
La ley de control está dada por:
𝑈(𝑧) =[1 + 𝐾1[𝑧𝐼 − 𝐴 + 𝐿𝐶]−1𝐵]−1[𝐾𝑖𝑧[𝑅(𝑧) − 𝑌(𝑧)] − (𝑧 − 1)𝐾1[𝑧𝐼 − 𝐴 + 𝐿𝐶]−1𝐿𝑌(𝑧)]
𝑧 − 1
𝑈(𝑧)
=[4.7904𝑧3 − 5.2694𝑧2 + 1.4371𝑧]𝑅(𝑧) − [4.7904𝑧3 + 270.268𝑧2 − 449.658𝑧 + 175.557𝑧]𝑌(𝑧)
𝑧3 − 1.4301𝑧2 − 17.3𝑧 + 17.7304
Luis Edo García Jaimes
DIAGRAMA DE BLOQUES PARA EL CONTROL DEL REACTOR
Luis Edo García Jaimes
RESPUESTA DEL SISTEMA CON EL CONTROLADOR DISEÑADO
Luis Edo García Jaimes
EJEMPLO DE CONTROL PI PARA EL REACTOR
Luis Edo García Jaimes