ANÁLISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO INTRODUCCIÓN En los capítulos anteriores se desarrollaron las herramientas necesarias para el análisis y el diseño de un sistema de control en tiempo continuo. En este capítulo se desarrollará el modelo necesario para el análisis de un sistema de control en tiempo discreto, conocido como sistema de datos muestreados, en razón de que las señales utilizadas se obtienen mediante el proceso de muestreo de señales de campo en tiempo continuo. El nuevo modelo será reconocido como la Función de Transferencia de Pulsos (FTP), cuyo nombre se debe al modelo ideal utilizado para representar el proceso de muestreo la señal continua x(t). Una vez definido el concepto de FTP, se aplicará a esquemas prácticos que incluyen varios lazos de control, desarrollando un procedimiento general y sistemático para encontrar la FTP equivalente de lazo cerrado. En este análisis comprobaremos una limitación del modelo FTP en el sentido de que solo permite obtener la respuesta () yt para cada instante de muestreo t kT = . Sin embargo, se demostrará que es posible desarrollar herramientas para evaluar la respuesta entre intervalos de muestreo, utilizando la Transformada de Laplace (TL) y la Transformada Z Modificada (TZM). Esta última es de gran interés para el tratamiento de sistemas de datos muestreados con atrasos. Utilizando la FTP de lazo cerrado se analizará el efecto del período de muestreo T en la respuesta dinámica del sistema de control. El concepto de ecuación característica será fundamental para evaluar una de las componentes de la respuesta dinámica: la repuesta transitoria. La respuesta permanente, será evaluada a través del error estacionario como una medida de la calidad del sistema de control, utilizando un modelo similar al que fue desarrollado en el capítulo 2 para el sistema de control en tiempo continuo. Para facilitar el cálculo de los valores característicos de la respuesta transitoria, será necesario analizar la correlación que existe entre los puntos del plano-s y puntos del plano-z, utilizando la regla de transformación sT z e = que es el fundamento de la transformada Z. Una vez analizada la respuesta dinámica, se abordaran temas relacionados con la estabilidad del sistema de control digital, aplicando herramientas similares a las utilizadas en el capítulo 3 para los sistemas de control en tiempo continuo y desarrollando otras que solo son aplicables los modelos discretos. Finalmente, se tratarán aspectos relacionados con el método del lugar de las raíces y la respuesta de frecuencia, que al igual que en los sistemas continuos o analógicos, son herramientas básicas para el análisis y diseño de sistemas de control digital. 5
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
ANÁLISIS DE SISTEMAS DE CONTROL
EN TIEMPO DISCRETO
INTRODUCCIÓN
En los capítulos anteriores se desarrollaron las herramientas necesarias para el análisis y
el diseño de un sistema de control en tiempo continuo. En este capítulo se desarrollará el modelo
necesario para el análisis de un sistema de control en tiempo discreto, conocido como sistema de
datos muestreados, en razón de que las señales utilizadas se obtienen mediante el proceso de
muestreo de señales de campo en tiempo continuo. El nuevo modelo será reconocido como la
Función de Transferencia de Pulsos (FTP), cuyo nombre se debe al modelo ideal utilizado para
representar el proceso de muestreo la señal continua x(t).
Una vez definido el concepto de FTP, se aplicará a esquemas prácticos que incluyen
varios lazos de control, desarrollando un procedimiento general y sistemático para encontrar la
FTP equivalente de lazo cerrado. En este análisis comprobaremos una limitación del modelo FTP
en el sentido de que solo permite obtener la respuesta ( )y t para cada instante de muestreo
t kT= . Sin embargo, se demostrará que es posible desarrollar herramientas para evaluar la
respuesta entre intervalos de muestreo, utilizando la Transformada de Laplace (TL) y la
Transformada Z Modificada (TZM). Esta última es de gran interés para el tratamiento de
sistemas de datos muestreados con atrasos. Utilizando la FTP de lazo cerrado se analizará el efecto del período de muestreo T en la
respuesta dinámica del sistema de control. El concepto de ecuación característica será
fundamental para evaluar una de las componentes de la respuesta dinámica: la repuesta
transitoria. La respuesta permanente, será evaluada a través del error estacionario como una
medida de la calidad del sistema de control, utilizando un modelo similar al que fue desarrollado
en el capítulo 2 para el sistema de control en tiempo continuo. Para facilitar el cálculo de los
valores característicos de la respuesta transitoria, será necesario analizar la correlación que existe
entre los puntos del plano-s y puntos del plano-z, utilizando la regla de transformación sTz e=
que es el fundamento de la transformada Z. Una vez analizada la respuesta dinámica, se abordaran temas relacionados con la
estabilidad del sistema de control digital, aplicando herramientas similares a las utilizadas en el
capítulo 3 para los sistemas de control en tiempo continuo y desarrollando otras que solo son
aplicables los modelos discretos. Finalmente, se tratarán aspectos relacionados con el método del
lugar de las raíces y la respuesta de frecuencia, que al igual que en los sistemas continuos o
analógicos, son herramientas básicas para el análisis y diseño de sistemas de control digital.
5
5 - 2 Capítulo 5 – ANALISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO
5.1 SISTEMAS DE CONTROL EN TIEMPO DISCRETO
En esta sección se hará una revisión breve de los conceptos fundamentales, necesarios para
el análisis de un sistema de control en tiempo discreto, cuyo modelo básico es la función de
transferencia de discreta (FTD) que se obtiene [ReySoto06] evaluando la ecuación en
diferencias (EED) en el dominio-z y utilizando como herramienta la transformada Z.
Modelos clásicos de sistemas control
En los capítulos 2, 3 y 4 se utilizó el diagrama en forma canónica mostrado en la figura 5.1,
para el análisis y diseño del sistema de control en tiempo continuo, el cual se desarrolló
aplicando el concepto de función de transferencia (FT) para modelar el controlador ( )cG s y
el proceso ( )pG s . De modo similar, utilizando la transformada Z es posible desarrollar la forma canónica
mostrada en la figura 5.2 que será utilizado para el análisis y diseño del sistema de control
en tiempo discreto, donde ( )D z es el modelo del controlador o compensador el cual será
desarrollado usando el concepto de función de transferencia discreta (FTD) [ReySoto06].
Por otro lado, ( )G z es el modelo discreto equivalente del proceso ( )pG s que será obtenido
aplicando el concepto de función de transferencia de pulsos (FTP).
En aplicaciones recientes, la función que realiza el controlador o compensador ( )D z en el
sistema de la figura 5.2, se puede lograr mediante el desarrollo e implantación de un
algoritmo de control en un microcontrolador ( )Cµ , tal como se muestra en la figura 5.3.
Sin embargo, como el Cµ solo reconoce magnitudes digitales (valores binarios) es
Figura 5.1 Diagrama en forma canónica para el análisis y diseño de un sistema de control en tiempo continuo.
Figura 5.2 Diagrama en forma canónica para el análisis y diseño de un sistema de control en tiempo discreto.
−
+ ( )R s ( )M s( )cG s
( )Y s ( )E s( )pG s
( )H s
−
+ ( )R z ( )M z( )D z
( )Y z( )E z( )G z
( )H z
5.1 – SISTEMAS DE CONTROL EN TIEMPO DISCRETO 5 - 3
necesario utilizar un convertidor analógico digital (A/D) y un convertidor digital analógico
(D/A), como dispositivos de interfaz del sistema de control. El convertidor A/D transforma la señal de error continua ( )e t en una señal que pueda ser
procesada digitalmente por el µC. Por otro lado, el convertidor D/A traduce la señal de
control digital en una señal continua ( )m t que sea capaz de modificar la planta o proceso,
para mantener la variable controlada ( )y t en un valor cercano al valor deseado ( )r t .
Algoritmo básico del controlador digital
Como una primera aproximación del problema a considerar, asumiremos que se trata de
modelar un controlador PI cuya ecuación característica en el dominio-s según ecuación
(4.9), viene dada por:
( ) ( )ip
kM s k E s
s
= + ⋅
(5.1)
La expresión equivalente de (5.1) en el dominio del tiempo es:
0
( ) ( ) ( )t
p im t k e t k e t dt= ⋅ + ∫ (5.2)
donde ( )e t es la señal de error y ( )m t la señal de control, mostradas en la figura 5.3. Por
otro lado, pk es la ganancia de la acción proporcional y ik la ganancia de la acción
integral, cuyos valores dependen del diseño específico del controlador y son ajustados
como parte del proceso de entonación del lazo de control. Como el microcontrolador de la figura 5.3 está programado para efectuar operaciones de
suma y multiplicación, debe ser posible simular la expresión anterior. Sin embargo su
evaluación deberá realizarse mediante aproximación numérica de la integral, tal como se
muestra en la figura 4.4, en la cual se asume el método del rectángulo del lado derecho o
de la diferencia posterior [ReySoto06].
−
+ r(t) m(t)µC
Gp(s)
y(t) e(t)
Sensor o Transmisor
A/D
D/A
Controlador/Compensador Planta/Proceso Figura 5.3 Esquema típico de un sistema de control digital directo.
5 - 4 Capítulo 5 – ANALISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO
Si ( )x t representa la integración numérica de ( )e t , haciendo t kT= en (5.2) obtenemos:
( ) ( ) ( )p im kT k kT k x kT= + (5.3)
De la figura 5.4 podemos expresar el valor actual de la integral ( )x kT , como
( ) [( 1) ] ( )x kT x k T e kT T= − + × (5.4)
siendo T el intervalo de integración. Sustituyendo (5.4) en (5.3) obtenemos:
( ) ( ) ( ) [ 1) ]p i im kT k kT e kT k x k T= + + − (5.5)
La expresión anterior constituye el algoritmo básico para desarrollar el controlador PI, que
utilizando notación comprimida, puede formularse como: 1( )k p i k i km k kT e k x −= + + (5.6)
La expresión (5.6) se reconoce como el algoritmo de control y corresponde a una ecuación en diferencias (EED) de primer orden, cuya forma general [ReySoto07] es 1 0 1( ) ( 1) ( ) ( 1)y k a y k b x k b x k+ − = + − (5.7)
la cual establece la relación entrada-salida ( ) ( )x k y k↔ . En la expresión anterior se ha
normalizado el coeficiente de ( )y k para que sea 1; esta EED se refiere como de diferencias
hacia atrás, por cuanto hace referencia a valores anteriores de ( )y k y ( )x k . En un caso
más general, la EED hacia atrás de un sistema discreto LIT de orden - n viene dada por:
1 0
( ) ( ) ( )n m
i j
i j
y k a y k i b x k j= =
+ − = −∑ ∑ (5.8)
En la expresión anterior asumimos que m n≤ para garantizar la causalidad del sistema de
control. Como se demostrará en el capítulo 6 la causalidad garantiza la realización del
algoritmo de control. Una forma alterna de (5.8) se consigue formulando la ecuación en
términos de diferencias hacia delante, que se consigue cambiando k por k n+ en (5.8)
1 0
( ) ( ) ( )n m
i j
i j
y k n a y k n i b x k n j= =
+ + + − = + −∑ ∑ (5.9)
( )e t
t
kT ( 1)k T−
Figura 5.4 Aproximación numérica de la integral en el controlador PI.
5.1 – SISTEMAS DE CONTROL EN TIEMPO DISCRETO 5 - 5
Aunque algebraicamente (5.8) y (5.9) son equivalentes, en la implementación del algoritmo
de control se prefiere usar la representación de diferencias hacia atrás, para garantizar su
causalidad al hacer referencia a valores que han ocurrido anteriormente. La notación de
diferencias hacia delante es útil en el desarrollo del modelo de estado, tal como se verá en
el capítulo 7. Finalmente, para efectos de notación en sistemas de datos muestreados, si el
período de muestreo T es constante, se considera que ( ) ( )y k y kT≡ y ( ) ( )x k x kT≡ .
Estrategias de diseño del controlador digital
El propósito final del diseño del controlador digital es calcular la FT ( )D z mostrada en la
figura 5.2, para satisfacer requerimientos de diseño del sistema de control, asumiendo que
se conoce el modelo equivalente discreto G(z) del proceso. A diferencia del diseño del
controlador ( )cG s en el dominio continuo, para el diseño del controlador digital existen dos
estrategias posibles, conocidas como método directo y método indirecto.
1. Método directo: Esta estrategia utiliza el diagrama de la figura 5.2 para diseñar el
controlador ( )D z a partir del modelo discreto equivalente ( )G z del proceso, usando
métodos clásicos en el plano-z. En una primera fase es necesario determinar ( )G z a
partir del modelo continuo del proceso: ( )pG s , utilizando métodos de transformación
o de discretización [ReySoto06].
2. Método indirecto: En el diagrama de la figura 5.1, utilizando métodos clásicos en el
plano-s, es posible diseñar ( )cG s a partir del modelo continuo del proceso ( )pG s .
Una vez obtenido ( )pG s se desarrolla su equivalente discreto ( )D z usando métodos
de transformación o de discretización [ReySoto06].
El método directo se aplica generalmente cuando el controlador no existe y se debe diseñar
por primera vez. Los métodos de diseño utilizando esta estrategia se presentarán en el
capítulo 6; sin embargo, antes será necesario desarrollar el modelo equivalente discreto del
proceso ( )G z , de lo cual nos ocuparemos con detalle en este capítulo. El método indirecto
tiene aplicación práctica cuando ya existe el controlador analógico y se desea sustituirlo por
uno digital. El diseño del controlador usando esta estrategia se presentará en el capítulo 6,
donde se evaluará el efecto del convertidor D/A en la forma final del controlador ( )D z .
El problema de diseño del controlador digital
La EED presentada en (5.8) puede considerarse como el modelo elemental de un filtro
discreto LIT, que usualmente se reconoce como un filtro digital . El diseño de este filtro
implica algo más que la simple determinación de los coeficientes ia y ib de la EED, tal
como sucede en el diseño del filtro continuo. En nuestro caso, existen un conjunto de
5 - 6 Capítulo 5 – ANALISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO
elementos adicionales que deben ser tomados en cuenta. En este sentido podemos
establecer dos fases a considerar en el diseño del controlador digital: Aspectos básicos del diseño:
1. Seleccionar el orden-n del sistema que corresponde al orden de la EED, el cual depende
del grado de exactitud y la complejidad del modelo a lograr.
2. Determinar los coeficientes ia y ib de la EED, que satisfagan las especificaciones de
diseño del filtro digital.
3. Seleccionar el intervalo de integración T para lograr una adecuada representación del
modelo discreto del sistema. Este valor se reconocerá más adelante como el período de
muestreo del sistema de control digital y su efecto en el comportamiento del sistema
será analizado con detalle en este capítulo.
Aspectos complementarios:
1. Seleccionar el tamaño adecuado de la palabra digital para minimizar los errores de
redondeo y truncamiento, en el proceso de almacenamiento interno del µC.
2. Controlar y eliminar el ruido digital generado durante la fase de implementación del
algoritmo de control (errores de redondeo y cuantización). Un caso especial ocurre en
la implementación de la acción derivativa del controlador PID.
5.2 SISTEMAS DE CONTROL DE DATOS MUESTREADOS
En la sección 5.1 se demostró que la función de transferencia discreta (FTD) es un modelo
adecuado para evaluar la respuesta dinámica de un sistema discreto, utilizando como
herramienta básica la transformada Z. En esta sección se desarrollará el modelo del
sistema de control de la figura 5.3, en el que intervienen señales obtenidas por muestreo de
una señal continua, reconocido como sistema de datos muestreados. Aunque en esencia
este sistema es discreto, la evaluación del efecto del muestreo de la señal de error ( )e t ,
realizado por el convertidor A/D y de la reconstrucción de la señal de control ( )m kT
realizada por el convertidor D/A, permitirá demostrar que no es posible aplicar el modelo
de FTD. Este análisis conducirá a la transformada estrella de Laplace (TEL) como modelo
del sistema de datos muestreados.
Muestreo y reconstrucción de una señal
Las aplicaciones prácticas del control digital utilizan el procesamiento digital de una señal
continua ( )x t para someterla a dos acciones fundamentales: muestreo y reconstrucción. El
muestreo es realizado por el convertidor A/D de la figura 5.3, para transformar la señal
continua ( )e t en una señal muestreada ( )e kT , que luego es sometida a un proceso de
5.3 – SISTEMAS DE CONTROL EN TIEMPO DISCRETO 5 - 7
cuantización para obtener la señal digital ( ) ( )e k e kT≡ . Esta señal digital al ser procesada
por el microcontrolador genera como respuesta otra señal digital ( ) ( )m k m kT≡ . La salida del procesador es la señal de muestreada ( )m kT , que es reconstruida por el
convertidor D/A de la figura 5.3 para transformarla nuevamente en una señal continua
( )m t , de modo que pueda ser aplicada a los dispositivos analógicos finales. En esencia, el
sistema de control de datos muestreados se pueden identificar tres tareas fundamentales:
- muestreo de la señal de error ( )e t para obtener la señal muestreada ( )e kT
- procesamiento de la señal ( )e kT para obtener la señal ( )m kT
- reconstrucción de la señal de control ( )m t a partir de la señal muestrea ( )m kT Para modelar el muestreo y reconstrucción de una señal continua se recurre al dispositivo
M-R (muestreador-retensor) mostrado en la figura 5.7, en la cual se ha omitido por el
momento el microcontrolador de la figura 5.3. Para simular el muestreo de la señal
continua ( )x t realizado por el convertidor A/D, se utiliza un interruptor lógico que
permanece abierto un tiempo T y cerrado un tiempo 0t . El convertidor D/A es simulado por
el retensor que procesa la señal muestreada *( ) ( )x t x kT≡ para obtener la señal
reconstruida ɵ( )x t , donde se espera que ɵ( ) ( )x t x t≈ .
La figura 5.8 presenta la salida probable del muestreador donde la señal *( )x t es ahora un
tren de pulsos de duración 0t y período T. Se logra así una representación aproximada
( ) *( )x kT x t≈ de la señal muestreada, que se puede mejorar en la medida en que 0t T<< .
Figura 5.7 Dispositivo M-R para simular el muestreo y reconstrucción de una señal.
Figura 5.8 Señal muestreada como un tren de pulsos.
*( )x t ɵ( )x t RETENSOR
( )x t
T
Dispositivo M-R
0t
5 - 8 Capítulo 5 – ANALISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO
Una primera aproximación de la característica operacional del retensor se muestra en la
figura 5.9, donde la reconstrucción se consigue manteniendo constante el último valor
muestreado mientras dure abierto el interruptor lógico de la figura 5.8. La señal
reconstruida ɵ( )x t se aproximará más a la señal ( )x t en la medida en que el período de
muestreo T sea pequeño. Sin embargo, esto afecta la aproximación de la señal muestreada
de la figura 5.8, ya que exige un valor aún más pequeño de 0t para que *( ) ( )x t x kT≈ .
Transformada estrella de Laplace
Como la señal reconstruida ɵ( )x t es continua por intervalos, es necesario recurrir al uso de
la transformada de Laplace. Para esto expresamos ɵ( )x t como
ɵ( ) (0)[ ( ) ( )] ( )[ ( ) ( 2 )]x t x u t u t T x T u t T u t T= − − + − − − +⋯ (5.10)
Llevando la ecuación (5.10) al dominio-s, obtenemos
21( ) (0) ( ) (2 )
sTsT sTe
X s x x T e x T es
−− −−
= + + + ⋯
Expresando el término entre corchetes como una sumatoria, la salida del retensor es
0
1( ) ( )
sTksT
k
eX s x kT e
s
− ∞−
=
−= ∑ (5.11)
Como la sumatoria solo depende de los valores de *( ) ( )x t x kT= y al evaluarla resulta un
función en el dominio-s, puede expresarse como
0
*( ) ( ) ksT
k
X s x kT e∞
−
=∑≜ (5.12)
Esta expresión se define como la transformada estrella de Laplace (TEL) de la señal
muestreada *( ) ( )x t x nT= . Sustituyendo (5.12) en (5.11), obtenemos
1( ) *( )
sTeX s X s
s
−−= (5.13)
Figura 5.9 Señal muestreada como un tren de pulsos.
5.3 – SISTEMAS DE CONTROL EN TIEMPO DISCRETO 5 - 9
Asociando (5.13) con la figura 5.8, la función de transferencia del retensor es
0
1( )
sTeG s
s
−−= (5.14)
conocido como retensor de orden-0 (ZOH). Las expresiones anteriores son suficientes para
lograr un primer modelo del dispositivo M-R en el dominio-s, mostrado en la figura 5.10.
Para desarrollar el modelo del muestreador asumimos que la señal ( )x t es sometida a un
proceso de modulación por impulsos, como se muestra en la figura 5.11. Este modelo es
solo una aproximación ideal del modelo original considerado en la figura 5.9, ya que ahora
la señal muestreada se representa por un tren de impulsos, en lugar de un tren de pulsos. Sin
embargo, a continuación se demostrará que este modelo ideal del muestreador es suficiente
para lograr el modelo del dispositivo M-R.
A partir de este modelo ideal del muestreador, podemos expresar *( )x t como
[ ]0
*( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( 2 )T
k
x t x t t x t t kT x t t t T t T∞
== ⋅δ = ⋅ δ − = ⋅ δ + δ − + δ − +∑ ⋯ (5.15)
Aplicando la propiedad de muestreo de la señal impulso [Strum00] o [ReySoto07]
*( ) (0) ( ) ( ) ( ) (2 ) ( 2 )x t x t x T t T x T t T= δ + δ − + δ − +⋯
Llevando la expresión anterior al dominio-s, obtenemos
2
0
*( ) (0) ( ) (2 ) ( )sT sT ksT
k
X s x x T e x T e x kT e∞
− − −
== + + + =∑⋯ (5.16)
*( )X s ( )X s 1 sTe
s
−−
( )X s
T
Figura 5.10 Modelo del dispositivo M-R.
( )x t *( )x t
( )T tδ
X
Figura 5.11 Modelo ideal del muestreador.
5 - 10 Capítulo 5 – ANALISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO
Se logra así la misma expresión (5.12) como señal de salida del muestreador y como
entrada del retensor. Luego, el modelo ideal del muestreador como un tren de impulsos, es
compatible con el modelo aproximado del retensor propuesto en la figura 5.9.
Resumiendo, el modelo de la figura 5.10 representa con propiedad las operaciones de
muestreo y reconstrucción del sistema de datos muestreados, donde *( ) *( )X s x t↔ se
define como la transformada estrella de Laplace. Si *( )x t tiene discontinuidades en
t kT= se deberá tomar el valor *( )x t+ en la evaluación de *( )X s .
Métodos para evaluar la TEL
Aunque es posible aplicar la definición (5.12) para evaluar la transformada estrella de
Laplace (TEL) de una señal muestreada, existen formas más prácticas basadas en
- método de equivalencia con la transformada Z
- método de residuos modificado
El método de equivalencia se basa comparar la definición de la TEL dada en (5.12) con la
definición (C.3) de la transformada Z. En efecto, se observa que existe la relación
*( ) ( ) sTz eX s X z == (5.17)
donde la expresión sTz e= se reconoce como la regla de transformación en el sentido de
que permite establecer la relación *( ) ( )X s X z↔ entre la señal muestreada *( )x t en el
dominio-s y su equivalente en el dominio-z. Como ( )X z se logra por transformación de *( )X s , se reconocerá en adelante como el
modelo equivalente en el dominio-z de *( )x t . Esto permite formular la relación
*( ) ( )x t X z↔ para una señal muestreada y establecer una diferencia conceptual con la
relación ( ) ( )x k X z↔ , donde ( )X z es simplemente la TZ de la señal discreta ( )x k .
Ejemplo 5.1: Evaluar la TEL de señal 2( ) 5 tx t e−= , asumiendo 0.5T s= .
Solución: Muestreando la señal ( )x t , obtenemos 2( ) 5 nTx nT e−= . Aplicando la tabla C.1
2 0.5 1
5 5( )
1 0.3679
zX z
z e z− × −= =− −
Aplicando la transformación sTz e= para 0.5T s= , obtenemos
TEL de una señal continua arbitraria.
5.3 – SISTEMAS DE CONTROL EN TIEMPO DISCRETO 5 - 11
0.5
1*( )
1 0.3679 sX s
e−=−
El resultado del ejemplo 5.2 muestra que la TEL de una señal muestreada es una fracción
no racional, que puede presentar dificultades en su manipulación algebraica. El método de residuos modificado [Kuo92] permite obtener directamente el modelo
equivalente discreto ( )X z de la señal muestreada *( )x t a partir del modelo ( )X s de la
señal continua ( )x t . En este método ( )X z se expresa como la suma de los residuos
evaluados a partir de lo polos de ( )X s , como:
1 ( )
( ) ( )n
sTi polos X s
zX z X s
z e==
−∑Residuos (5.18)
Según la forma de estos polos pueden ocurrir 2 casos: Caso 1: Polos simples, reales o complejos
( ) ( ) ( )i
i i sTs s
zR z s s X s
z e =
= −−
(5.19)
Caso 2: Polos múltiples, reales o complejos
1
1
1( ) ( ) ( )
( 1)!i
mm
i im sT
s s
d zR z s s X s
m ds z e
−
−=
= − − − (5.20)
Ejemplo 5.2: Obtener la TZ equivalente de una señal muestreada *( )y t con 0.1T s= , si
1( )
( 1)( 2)Y s
s s=
+ +
Solución: Aplicando (5.19) para los 2 polos simples, obtenemos
110.1
1( )
( 2) 0.9048sTsT
z zR z
s z e z=−=
= =+ − −
220.2
1( )
( 1) 0.8187sTsT
z zR z
s z e z=−=
−= =+ − −
Sumando estos residuos obtenemos finalmente,
Método de residuos modificado con polos simples.
5 - 12 Capítulo 5 – ANALISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO
0.0861
( )0.9048 0.8187 ( 0.9048)( 0.8197)
z z zY z
z z z z= − =
− − − −
Ejemplo 5.3: Obtener el modelo equivalente discreto de una señal muestreada *( )x t , para
0.1T s= , asumiendo que la TL de ( )x t viene dada por
2
1( )
( 1)X s
s s=
+
Solución: Existen 2 polos, uno simple en 1 1p = − y otro en 2 0p = con multiplicidad
2m = , que generan 2 residuos: 1( )R z y 2( )R z . Aplicando (5.19) para el polo
real simple obtenemos
1 210.1
1( )
0.9048sTsT
z zR z
s z e z=−=
= =− −
Para el polo múltiple aplicamos (5.20)
2 2 200.1
1 0.1 ( 1.1)( )
1 1 ( 1) ( 1)sTsT
d z z z z zR z
ds s z e z z z==
− − = = + = − + − − − −
Sumando los 2 residuos
2
0.0048374 ( 0.9672)( )
( 1) ( 0.9048)
z zX z
z z
+=− −
De este modo se puede decir que existen dos modelos para el análisis de una señal
muestreada *( )x t : la TEL *( )X s en el dominio-s y la TZE ( )X z en el dominio-z. Para
facilitar el cálculo de ( )X z , se desarrolló la función especial residuosm() en MATLAB ®,
cuya descripción se presenta en el apéndice B. Su sintaxis es:
Xz = residuosm(Xs,T)
Introduciendo el nombre de esta función se consigue ayuda y un ejemplo demostrativo.
0.0048374 z (z+0.9672) ---------------------- (z-0.9048) (z-1)^2
Método de residuos modificado con polos múltiples
5.3 – SISTEMAS DE CONTROL EN TIEMPO DISCRETO 5 - 13
Sampling time: 0.1
Aplicabilidad y limitaciones del modelo de la TEL y la TZE
El modelo *( )X s y su equivalente ( )X z son suficientes para modelar la relación entre la
señal continua ( )e t y la señal reconstruida ( )e tɵ , formulada a través del dispositivo M-R.
Sin embargo, presentan las siguientes limitaciones:
1. Los dos componentes del dispositivo M-R mostrados en la figura 5.10 no modelan
ningún elemento físico. Simplemente permiten formular con un cierto grado de
aproximación la relación entre la señal continua ( )X s y la señal reconstruida ( )X s .
2. Aunque la señal *( )X s permite interconectar los modelos del muestreador ideal y del
retensor, tiene el inconveniente de que su expresión en el dominio-s no es una
fracción racional y por lo tanto su manipulación algebraica es compleja.
3. Este problema se supera si en lugar del modelo *( )X s se utiliza su modelo
equivalente ( )X z , que sí es una fracción racional.
4. Sin embargo el modelo equivalente ( )X z solo permite obtener valores de la señal
reconstruida ɵ( )x t en los instantes de muestreo, es decir ɵ( )x kT .
Para obtener valores entre intervalos de muestreo se puede utilizar el método de la
transformada de Laplace que será presentado en la sección 5.4. Existen otros métodos para
evaluar la respuesta entre intervalos de muestreo [ReySoto06] como la transformada Z
modificada y el modelo de estado, que están fuera del alcance de este libro.
Espectro de frecuencia de la señal muestreada
Aplicando la propiedad de convolución de la transformada de Fourier a la señal
*( ) ( ) ( )Tx t x t t= ⋅δ de la figura 5.11, se puede demostrar [ReySoto07] que su espectro es
1
*( ) ( )m
k
X j X j jkT
∞
=−∞ω = ω − ω∑ (5.21)
donde 2 /m Tω = π es la frecuencia de muestreo en [rad/s]. Interpretando (5.21) se
desarrolló la figura 5.12, la cual permite reconocer que el muestreo de una señal continua
( )x t de banda limitada con espectro aperiódico ( )X jω , genera la señal ( )x kT con
espectro periódico con período 2 /m Tω = π , siendo mω la frecuencia de muestreo en rad/s
y T el período de muestreo en segundos.
Figura 5.12 Espectro de frecuencia de una señal muestreada.
ω
( )X ω
B−ωBω
A
mω
m−ω
ω
*( )X ω
B−ω Bω
A/T
½ mω
½ m− ω
T
5 - 14 Capítulo 5 – ANALISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO
De acuerdo con la figura 5.12, la reconstrucción de una señal muestreada en el dominio-ω
se reduce a un proceso de filtrado para eliminar los componentes armónicos por encima de
la frecuencia ½ mω . En efecto utilizando el filtro pasa-bajo ideal mostrado en la figura 5.12,
con amplitud T y ancho de banda ½ mω debería ser posible recuperar el espectro ( )X jω de
la señal continua, a partir del espectro *( )X jω de la señal muestreada. Sin embargo, para lograr la recuperación de la señal muestreada el ancho de banda Bω de
la señal a muestrear debe ser inferior al ancho de banda ½ mω del filtro ideal. Esta
condición se reconoce como el teorema de muestreo, que se formula así: Definición 5.1: Teorema de muestreo Una señal continua ( )x t de banda limitada con ancho de banda Bω , puede
ser reconstruida a partir de su versión muestreada *( ) ( )x t x nT= , si la
frecuencia de muestreo mω es superior al doble de Bω , es decir: 2m Bω > ω (5.22)
Definición 5.2: Frecuencia de Nyquist La frecuencia de Nyquist se define como
½N mω ω≜ (5.23)
y representa el ancho de banda del filtro ideal de la figura 5.12. En relación
con el teorema de muestreo, es la máxima frecuencia que puede estar
contenida en una señal muestreada *( )x t de ancho de banda Bω , para
poder recuperar la señal continua ( )x t . Luego, B Nω < ω .
Si B Nω > ω , se presenta solapamiento entre los espectros de los armónicos, tal como se
muestra en la figura 5.13, siendo imposible recuperar la señal muestreada. Esta situación se
reconoce como el efecto de aliasing de la señal reconstruida en el dominio del tiempo.
De las figuras 5.12 y 5.13 se puede concluir que para reconstruir la señal ɵ( )x t a partir de la
señal muestreada *( ) ( )x t x nT= , el retensor de la figura 5.12 debe ser un filtro pasa-bajo,
con ancho de banda N Bω > ω , siendo Bω es el ancho de banda de la señal continua ( )x t .
ω
*( )X ω
A/T
½ mω½ m− ω
Figura 5.13 Efecto de solapamiento en el espectro de una señal muestreada.
5.3 – SISTEMAS DE CONTROL EN TIEMPO DISCRETO 5 - 15
En aplicaciones prácticas se recomienda la selección de la frecuencia de muestreo como (10 20)m Baω ≈ ×ω (5.24)
La magnitud de la frecuencia mω puede estar limitada por los componentes del sistema de
control, particularmente por el efecto de calentamiento de los convertidores A/D y D/A
mostrados en la figura 5.3. En estos casos se puede recurrir a un filtro antialiasing o filtro
guardián con el objeto de reducir el ancho de banda de la señal continua ( )x t , antes de
someterla al proceso de muestreo.
Respuesta de frecuencia del retensor
Evaluando (5.14) para s j= ω es posible comprobar que el retensor de orden-0 (ZOH) es
un filtro pasa-bajo con amplitud T y ancho de banda ½N mω = ω . Evaluando
/ 20
1 ( / 2)( )
/ 2
j Tj Te sen T
G j T ej T
− ω− ω− ωω = =
ω ω
Sustituyendo / 2 / mTω = πω ω , obtenemos para la respuesta de magnitud
0
( / )| ( )|
/m
m
senG j T
πω ωω =πω ω
(5.25)
y para la respuesta de fase
0
0, ( / ) 0( ) ,
, ( / ) 0m
mm
senH j
sen
πω ω >πω∠ ω = − + ϕ ϕ = π πω ω <ω (5.26)
A partir de (5.25) y (5.26), para T=0.1s, se obtuvo la RDF mostrada en la figura 5.29.
Figura 5.15 Respuesta de frecuencia del retensor ZOH.
5 - 16 Capítulo 5 – ANALISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO
Para efecto de comparación, la figura 5.15 incluye la característica del filtro ideal de la
figura 5.12. A continuación se presenta un análisis de las características más importantes
que se derivan de esta respuesta de frecuencia:
1. La respuesta de magnitud del ZOH corresponde a un filtro pasa-bajo real con ancho
de banda 10Nω = π y amplitud 0.1A = , característica requerida para que el retensor
pueda recuperar la señal muestreada.
2. La fase negativa del ZOH reduce la estabilidad del sistema.
3. Disminuyendo el período de muestreo T se puede reducir el efecto desestabilizador.
La selección de mω usando el criterio formulado en (5.24) favorece esta condición.
4. Desde el punto de vista práctico la mayor parte de los convertidores D/A se
comportan como un ZOH.
Correspondencia entre el plano-s y el plano-z
En los párrafos anteriores se demostró que una señal continua ( )x t con espectro de
frecuencia aperiódico, al ser muestreada se produce una señal ( )x nT cuyo espectro es
periódico, con período 2 /m Tω = π , siendo mω la frecuencia de muestreo en rad/s. En este
párrafo se analizará el efecto de la periodicidad de *( )X s , desde el punto de vista de la
correlación entre los puntos del plano-s y el plano-z. Si partimos de la regla de transformación sTz e= que se aplicó en (5.17) para pasar del
plano-s al plano-z y considerando en principio los puntos sobre el eje imaginario, s j= ω
( )
01Ts j T j T jz e e e eσ+ ω ω Ω
σ== = = = = ∠Ω (5.27)
Luego los puntos sobre el eje-jω se transformarían en puntos sobre el círculo unitario, tal
como se muestra en la figura 5.16.
c
b
a x
j y
π/T
− π/T
3π/T
− 3π/T
σ
jω
Transformación varios-a-uno
Tsz ε=
Figura 5.16 Correspondencia entre el plano-s y el plano-z.
5.3 – SISTEMAS DE CONTROL EN TIEMPO DISCRETO 5 - 17
Sin embrago, evaluando la ecuación (5.27) para TΩ = ω en el intervalo [0, ]π , obtenemos
=0
: 1 0º 0,
: 1 90º /2 / 2 /4
: 1 180º / /2m
m
punto a z
punto b z T
punto c z T
= ∠ Ω = ω= ∠ Ω = π ω = π = ω= ∠ Ω = π ω = π = ω
(5.28)
Por lo tanto, los puntos de *( )X s en el intervalo [0, / ]Tπ se transforman sobre la parte
superior del círculo unitario. Los puntos para / 2m Nω > ω = ω se solapan sobre el círculo
unitario, debido a la periodicidad de je Ω , tal como se muestra en la figura 5.16.
Propiedades de la transformada estrella de Laplace
El resultado anterior permite identificar un conjunto de propiedades de la TEL, asociadas
con la periodicidad del espectro de la señal muestreada. Sustituyendo j sω = en (5.21)
obtenemos una expresión equivalente de *( )X jω en el dominio-s
1
*( ) ( )m
k
X s X s j kT
∞
=−∞= + ω∑ (5.29)
La ecuación (5.29) permite reconocer dos propiedades de *( )X s :
P1. *( )X s es periódica con período mjω , es decir *( ) *( )mX s X s j k= ± ω
Esta propiedad se puede demostrar usando la definición (5.12) de la TEL
( )
0
*( ) ( ) mnT s j k
m
n
X s j k x nT e∞
− ± ω
=± ω =∑ (5.30)
Evaluando la expresión anterior para 2 /m Tω = π
( ) 2 1 2mnT s j k nTs j kn nTs nTse e e e kn e− ± ω − π − −= ⋅ = ⋅ ∠ π =∓ ∓
Sustituyendo esta expresión en (5.30) obtenemos
0
*( ) ( ) *( )nTs
m
n
X s j k e nT e X s∞
−
=± ω = =∑
De acuerdo con este resultado, *( )X s es periódica en el plano-s, con período ωm.
P2. Si ( )X s tiene un polo en oss = , *( )X s tendrá polos en o ms s j k= ± ω , 0, 1, 2,k = …
Esta propiedad puede demostrarse desarrollando (5.29)
1
*( ) [ ( ) ( ) ( ) ( 2 ) ]m m mX s X s j X s X s j X s jT
= + − ω + + + ω + + ω +⋯ ⋯ (5.31)
Si ( )X s tiene un polo en 1s p= , cada término de *( )X s en (5.31) contribuirá con un
polo en 1 ms p j k= ± ω , para 0,1,2,k = ∞… .
5 - 18 Capítulo 5 – ANALISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO
5.3 MODELO DEL SISTEMA DE DATOS MUESTREADOS
En la sección anterior se demostró que utilizando el dispositivo M-R es suficiente para
analizar las características que se derivan del muestreo y reconstrucción de una señal. En esta
sección se utilizará este dispositivo para desarrollar un modelo para la representación del
sistema de datos muestreados en el dominio-z, a través del cual es posible evaluar su
comportamiento dinámico en lazo abierto y en lazo cerrado, ante una entrada arbitraria.
Función de transferencia de pulsos
Consideremos el sistema de la figura 5.17, en la cual se muestra el dispositivo M-R, en
cascada con el modelo continuo del proceso ( )pG s de un sistema de control de datos
muestreados, donde ( )E s es la señal de error.
Considerando el bloque ZOH y ( )pG s como ( )G s , la salida ( )Y s puede expresarse como
( ) ( ) *( )Y s G s E s= ⋅ (5.32)
Esta expresión es especial porque es el producto de una función aperiódica ( )G s y una
función periódica *( )E s . Tomando la transformada estrella de Laplace (TEL) de los dos
miembros y aplicando la integral de convolución (5.29)
* * * *1( ) ( ) ( ) ( ) ( )m m
n
Y s G s E s G s jn E s jnT
∞
=−∞
= ⋅ = + ω ⋅ + ω ∑
Aplicando la propiedad P1 de periodicidad de la TEL, *( ) *( )mE s jn E s+ ω = y como
*( )E s no depende de la sumatoria, obtenemos
1
*( ) *( ) ( ) *( ) *( )m
n
Y s E s G s jn E s G sT
∞
=−∞
= ⋅ + ω = ⋅∑ (5.33)
Comparando (5.32) con (5.33) se observa que si en un producto de dos funciones uno de
los factores está estrellado, es posible llevar toda la expresión al dominio-s*. De este modo
la función *( )G s permite establecer la relación entre la entrada-salida *( ) *( )E s Y s→ del
bloque punteado de la figura 5.17. A partir de (5.33) definimos
*
*
*
( )( )
( )
Y sG s
E s≜ (5.34)
*( )E s ( )E sZOH
( )E s
T ( )pG s
( )Y s
( )G s
Figura 5.17 Sistema de control de datos muestreados.
5.4 – MODELO DE SISTEMAS DE CONTROL DISCRETOS 5 - 19
como la función de transferencia de pulsos (FTP) del sistema mostrado la Figura 10.14.
Definición 5.3: Función de Transferencia de Pulsos La función de transferencia de pulsos (FTP) es la relación entre la TEL de
la salida y la TEL de entrada, de un sistema cuya entrada es muestreada.
Usando la transformación (5.17) puede expresarse en el dominio-z, como
( )
( )( )
Y zG z
E z≜ (5.35)
que se reconoce como el modelo discreto equivalente de *( )G s .
El nombre dado a *( )G s es para hacer referencia a que el modelo la señal de entrada
*( )E s se logra mediante la modulación por pulsos. Usando (5.34) y (5.35) se llega al
modelo equivalente del sistema de datos muestreados mostrado en la figura 5.18.
El muestreador ficticio mostrado en la salida del bloque ( )G s se usa para reconocer que el
modelo FTP permite obtener solo los valores de *( ) ( )y t y nT= . La expresión equivalente
de la salida en el dominio-z, de acuerdo con (5.35) es ( ) ( ) ( )Y z G z E z= , donde
y *( ) ( ) ( ) ( )sTz e
G z G s E z E s=
= =Z (5.36)
Considerando que G(s) en (5.36) incluye la FT del ZOH, usando (5.14) podemos escribir
1( ) ( ) ( )
sT
p
eG z G s G s
s
− −= = ⋅
Z Z
Interpretando el término sTe- como el atraso de una muestra en el dominio-z, obtenemos
1( )
( ) (1 )pG s
G z zs
− = − ⋅
Z (5.37)
Según la figura 5.18, la aplicación de (5.37) establece una condición implícita, en el sentido
de que la señal de entrada al bloque G(s) debe estar muestreada para que existe la FTP. La
expresión Z puede evaluarse utilizando uno de los siguientes métodos:
1. Método de residuos modificado.
2. Función especial residuosm de MATLAB .
( )G z ( )Y z ( )E z ( )Y s
( )G s
*( )E s ( )E s
T *( )Y s
T
Figura 5.18 Modelo discreto equivalente del sistema de datos muestreados.
5 - 20 Capítulo 5 – ANALISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO
Ejemplo 5.4: Obtener modelo discreto equivalente del sistema mostrado en la figura 5.16,
asumiendo 0.1T s= y que el proceso está modelado por
1
1)(
+=
ssG p
Solución: Aplicando (5.37)
)()1()1(
1)1()( 1
11 zGzss
zzG ⋅−=
+⋅−= −−Z
Para evaluar 1( )G z utilizamos (5.19) del método de residuos modificado
1 20 10.1 0.11
1 1( ) ( )
1 1 0.9048sT sTs sT T
z z z zR z R z
s z e z s z e z= =−= =
−= ⋅ = = ⋅ =+ − − − −
Sumando los dos residuos y sustituyendo en la expresión de ( )G z
1 0.095163 0.095163( ) (1 )
( 1)( 0.9048) 0.9048
zG z z
z z z
−= − ⋅ =− − −
Utilizando la función especial residuosm() obtenemos
5.4 – MODELO DE SISTEMAS DE CONTROL DISCRETOS 5 - 25
Sistemas de control de lazo abierto incluyendo el c ontrolador digital
El análisis hecho hasta el momento se ha dedicado a evaluar el efecto del dispositivo
muestreador-retensor (M-R) en la respuesta de un sistema de datos muestreados, sin incluir
el controlador. La figura 5.21 muestra un sistema de lazo abierto donde se incluye el
controlado digital, conocido como sistema de control digital directo.
El diagrama anterior representa la rama directa del diagrama de la figura 5.3, donde el
controlador digital es implementado en el µC mediante un algoritmo de control. Como se
demostró en la ecuación (5.6) el algoritmo de control es básicamente una EED que
establece la relación entrada-salida entre dos señales discretas: *( ) *( )e t m t→ o
( ) ( )e k m k→ . Por lo tanto llevando al dominio-z esta relación se obtiene ( ) ( )E z M z→ , a
partir de la cual es posible calcular la función de transferencia discreta (FTD):
( ) ( ) / ( )D z M z E z= , como modelo equivalente del controlador digital. La salida del sistema de la figura 5.19, considerando que el convertidor A/D es un ZOH es
1( ) ( ) ( ) ( ) *( ) ( ) *( )
sT
p p
eY s G s M s G s M s G s M s
s
−−= ⋅ = ⋅ ⋅ = ⋅ (5.41)
donde *( )M s es la salida del controlador. Asumiendo que ( )D z representa su FTD
)()()( zEzDzM ⋅=
Usando la transformación sTz e= podemos expresar en el dominio-s*, como
*( ) *( ) *( )M s D s E s= ⋅ (5.42)
Sustituyendo (5.42) en (5.41)
( ) ( ) *( ) *( )Y s G s D s E s= ⋅ ⋅
Como existe un solo factor no estrellado, podemos llevar al dominio-s*
*( ) *( ) *( ) *( )Y s G s D s E s= ⋅ ⋅
que finalmente puede representarse en el dominio-z como
( ) ( ) ( ) ( )Y z G z D z E z= ⋅ ⋅ (5.43)
donde ( )G z es el modelo equivalente discreto del proceso, calculado usando (5.37). Con
base en este resultado, la figura 5.22 muestra el modelo equivalente del sistema de control
directo en lazo abierto.
Figura 5.21 Sistema de control digital directo de lazo abierto.
M*(s)
e(t) A/D
y(t) CONTROLADOR
DIGITAL e*(t)
D/A m*(t)
E(s) E*(s) Y(s) )(sM
)(tm GP(s)
5 - 26 Capítulo 5 – ANALISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO
En el desarrollo de este modelo, conviene reconocer que
( )G z : Modelo discreto equivalente del proceso
Se calcula a partir de (5.37) con base en el concepto de FTP.
( )D z : Modelo discreto equivalente del controlador digital
Se calcula como una FTD, a partir del algoritmo de control (EED).
La interpretación de este modelo en el dominio-s se muestra en la figura 5.23, donde ( )G s
representa la combinación del ZOH con la función continua del proceso. Por la naturaleza
misma de las señales mixtas (continuas y discretas) que se manejan en este diagrama, se
refiere generalmente como un modelo seudo-continuo del sistema de control de datos
muestreados en lazo abierto.
La única limitación del modelo discreto de la figura 5.22 está en que no es posible llegar a
una expresión de la señal reconstruida ( )m t que se muestra en la figura 5.23. Se puede
demostrar [ReySoto06] que usando el modelo de estado es posible obtener esta señal ( )m t .
Ejemplo 5.7: Obtener la respuesta escalón del sistema de control digital mostrado en la
figura 5.21, asumiendo
1( ) , 0.5
1pG s ZOH con T s
s= =
+
Asumir que el algoritmo del controlador está dado por la siguiente EED
( ) ( ) 0.5 ( 1)m k e k e k= − −
Solución: Aplicando el concepto de FTP el modelo equivalente discreto de ( )pG s es
E(z) D(z)
M(z) G(z)
Y(z) Figura 5.22 Modelo equivalente del sistema de control digital directo de lazo abierto.
e*(t)
T
e(t) D*(s)
m*(t) GP(s)
s
e sT−−1
m(t)
G(s)
y(t)
y*(t)
Figura 5.23 Modelo seudo-continuo del sistema de control digital directo de lazo abierto.
FTP equivalente en sistemas en lazo abierto.
5.4 – MODELO DE SISTEMAS DE CONTROL DISCRETOS 5 - 27
1
2
1 0.3935( ) (1 )
( 1) 0.6065G z z
s s z
− = − ⋅ = + −
Z
Este resultado puede verificarse usando c2d() del TBC
Función de transferencia de pulsos de lazo cerrado
En este párrafo se demostrará que la ubicación del muestreador es determinante en la
obtención del modelo de un sistema de datos muestreados en lazo cerrado, lo cual
conducirá a la necesidad de desarrollar un método sistemático para el cálculo de la función
de transferencia de pulsos de lazo cerrado (FTPLC) del sistema: ( )T z . Consideremos en principio la forma canónica de control de lazo cerrado mostrado en la
figura 5.24, donde se asume que ( )G s es el modelo continuo de la rama directa, que de
acuerdo a lo señalado en la figura 5.17, es la combinación del ZOH y del proceso ( )pG s .
Por otro lado, ( )H s es el modelo continuo de la rama inversa asociada con el sistema de
medición. Sin pérdida de generalidad en este diagrama se ha omitido nuevamente la
presencia del controlador, que será incluido posteriormente.
En la figura 5.24 la señal de salida del sumador E(s) viene dada por
( ) ( ) ( ) ( )E s R s H s Y s= − ⋅ (5.44)
donde la salida ( )Y s puede expresarse como
( ) ( ) *( )Y s G s E s= ⋅ (5.45)
Sustituyendo (5.45) en (5.44), obtenemos
*( ) ( ) ( ) ( ) ( )E s R s H s G s E s= − ⋅ ⋅ (5.46)
Llevando al dominio estrellado y agrupando el producto *( ) ( )G s E s⋅
*( )E s
T
E(s) G(s)
Y(s)
*( )Y s
H(s)
+ −
R(s)
− Figura 5.24 Forma canónica del sistema de control de datos muestreados de lazo cerrado.
5.4 – MODELO DE SISTEMAS DE CONTROL DISCRETOS 5 - 29
** * *( ) ( ) ( ) ( )E s R s GH s E s= − ⋅ (5.47)
Despejando *( )E s
*
*
*
( )( )
1 ( )
R sE s
GH s=
+ (5.48)
Llevando (5.45) al dominio estrellado y sustituyendo *( )E s por (5.48), obtenemos
*
* *
*
( )( ) ( )
1 ( )
G sY s R s
GH s=
+ (5.49)
A partir de esta expresión es posible identificar a
*
*
*
( )( )
1 ( )
G sT s
GH s=
+ (5.50)
como la función de transferencia de pulsos (FTP) de lazo cerrado del sistema de la figura
5.24. Como todos los términos en (5.50) están estrellados es posible obtener una expresión
equivalente en el dominio-z, como
( )
( )1 ( )
G zT z
GH z=
+ (5.51)
que se reconoce como la FTP equivalente de lazo cerrado del sistema de la figura 5.24,
donde ( ) ( )G z G s=Z y ( ) ( ) ( )GH z G s H s= ⋅Z . Utilizando (5.51) es posible evaluar la
salida del sistema de lazo cerrado como ( ) ( ) ( )Y z T z R z= ⋅ . Sin embargo, esta expresión
restringe la respuesta a valores en cada instante de muestreo, lo cual se identifica con el
muestreador ficticio de la figura 5.24. El orden seguido anteriormente en el desarrollo de las expresiones algebraicas, facilitó la
obtención del modelo de *( )T s . Sin embargo, en diagramas más complejos es posible que
el resultado no sea tan explícito. En general, si la señal de error ( )E s no está estrellada, tal
como se muestra en la figura 5.25, no es posible obtener la función de transferencia de
pulsos de lazo cerrado (FTPLC).
T
E(s) G(s)
Y(s)
H(s)
+ −
R(s)
*( )Y s
Figura 5.25 Sistema de control de datos muestreados donde no es posible evaluar la FTPLC.
5 - 30 Capítulo 5 – ANALISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO
Siguiendo el procedimiento anterior, se puede demostrar [ReySoto06] que
*
*
*
( ) ( )( ) ( )
1 ( )1 ( )
GR s GR zY s Y z
GH zGH s= =
++ (5.52)
A pesar de que esta expresión permite obtener *( )Y s y ( )Y z , no existe forma de calcular la
FTPLC. De modo similar, se existen otras configuraciones típicas [Ogata95] algunas de las
cuales serán consideradas más adelante como parte de los ejemplos a desarrollar.
Procedimiento general para obtener la FTPLC: ( )T z
Los sistemas de datos muestreados en lazo cerrado manejan señales mixtas (continuas y
discretas) y el muestreador puede ser colocado en cualquier lugar del diagrama. Esto puede
resultar en un tratamiento algebraico de cierta complejidad, que generalmente no ocurre en
los sistemas continuos. Lo anterior sugiere una estrategia para garantizar que se llegue a
una expresión de la FTPLC, mediante un procedimiento general. El procedimiento exige la formulación de las relaciones causa-efecto en cada componente
del diagrama de bloques (DB) o del gráfico de flujo de señales (GFS) del sistema de
control, se desarrollen como ecuaciones en forma estándar, dejando del lado izquierdo las
señales de salida y del lado derecho las señales de entrada.
El cálculo de la FTPLC se hará aplicando la fórmula de ganancia de Mason (FGM)
presentada en la sección 1.5. El procedimiento consiste en:
1. Construir el GFS original (GFSO) identificando las señales de entrada y salida. Para esto se considera cada muestreador como un dispositivo abierto, tal como se
muestra en la figura 5.23. El desarrollo del GFSO se logra fácilmente si se consideran
solo las señales de entrada a cada bloque del DB. Como se observa en la figura 5.27,
el GFSO deberá incluir señales mixtas.
Una vez construido el GFSM se deben identificar las señales de entrada y salida.
Como cada muestreador es un dispositivo abierto, la señal muestreada *( )E s se
convierte en una entrada adicional y la señal no muestreada ( )E s es una salida
*( )E s E(s) G(s)
Y(s)
+ −
R(s)
X(s)
−
E(s) *( )E s
X(s)
R(s) Y(s)
−1
Figura 5.27 Representación del muestreador como un dispositivo abierto.
5.4 – MODELO DE SISTEMAS DE CONTROL DISCRETOS 5 - 31
adicional. Por lo tanto, cada muestreador genera en el GFSO un par de nodos
entrada-salida adicionales a los nodos propios del sistema.
2. Formular una ecuación en la forma estándar para cada salida.
Cada ecuación debe expresarse únicamente en función de las entradas del grafo.
3. Obtener la TEL de cada ecuación en la forma estándar. Se aplican las reglas desarrolladas sobre factores estrellados y no estrellados. Al final de esta etapa, todas las ecuaciones deben incluir solo señales estrelladas.
4. Construir el GFS modificado (GFSM) Interpretar cada ecuación estándar estrellada como una relación causa-efecto.
5. Aplicar la FGM para obtener las relaciones E-S del sistema en lazo cerrado.
Ejemplo 5.8: Aplicando el procedimiento anterior, obtener ( )T z para la forma canónica
de la figura 5.21, asumiendo que cada bloque está dado por
1 1( ) , 0.5 ( )
1 2G s ZOH con T s y H s
s s= = =
+ +
Solución: Aplicando el procedimiento propuesto, obtenemos
1. GFSO a partir del DB.
Salidas: ( )Y s , ( )E s
Entradas: ( )R s , *( )E s
2. Ecuación estándar para cada salida.
( ) ( ) *( )
( ) ( ) ( ) ( ) *( )
Y s G s E s
E s R s G s H s E s
= ⋅= − ⋅ ⋅
3. TEL de cada ecuación estándar.
Aplicando la regla sobre factores estrellados, obtenemos:
*( ) *( ) *( )
*( ) *( ) *( ) *( )
Y s G s E s
E s R s GH s E s
= ⋅
= − ⋅
4. GFSM a partir de las ecuaciones estándar estrelladas.
Interpretando gráficamente las dos ecuaciones anteriores
FTP de LC de forma canónica usando procedimiento general.
E(s) *( )E s G(s)
− H(s)
Y(s) R(s)
*( )GH s−
*( )E s *( )G s *( )R s *( )Y s
G(s)
Y(s)
5 - 32 Capítulo 5 – ANALISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO
5. FGM para obtener la relación *( ) *( )R s Y s→ A partir del GFSM, obtenemos
*( )*( ) *( )
1 *( )
G sY s R s
GH s= ⋅
+
Expresando en el dominio-z la FTPLC es
)(1
)()(
zGH
zGzT
+=
Se logra así el mismo resultado obtenido en la ecuación (5.50) por el método
algebraico. A partir de (5.37) evaluamos )(zGH para los valores dados
1 1 1 1( ) (1 )
1 2GH z z
s s s
− = − × × × + + Z
Evaluando el término entre corchetes usando residuos modificado
1 0.0774 ( 0.6065)
( 1)( 2) ( 1)( 0.6065)( 0.3679)
z z
s s s z z z
+= + + − − − Z
Sustituyendo en la expresión de ( )GH z , obtenemos
0.0774( 0.6065)( )
( 0.6065)( 0.3679)
zGH z
z z
+=− −
Este resultado se puede verificar usando c2d() del TBC
En el ejemplo anterior se observa que no es posible aplicar la función feedback() del TBC
para evaluar ( )T z del diagrama de la figura 5.24. Esta función se utilizó en el ejemplo 1.7
para calcular la FT de lazo cerrado ( )T s de la forma canónica del modelo continuo. La
razón se debe a la forma típica del denominador de ( )T z en (5.51).
Se puede demostrar [ReySoto06] que en el caso de la figura 5.28, donde la señal de entrada
al bloque de realimentación está muestreada, la FTPLC viene dada por
( )
( )1 ( ) ( )
G zT z
G z H z=
+ (5.53)
donde sí es posible aplicar la función feedback() para evaluar ( )T z , dado que están
definidas explícitamente ( )G z y ( )H z . La sintaxis es
Tz=feedback(Gz,Hz)
Un caso particular ocurre cuando el sistema tiene realimentación estática o realimentación
unitaria, donde también se puede aplicar esta función para el cálculo de la FTPLC. En el ejemplo 5.12 se demuestra que el procedimiento propuesto es efectivo para obtener la
respuesta *( )Y s y en algunos casos la FTP equivalente ( )T z del sistema de lazo cerrado,
con la limitación de que a partir de *( )Y s se obtendría la respuesta en cada instante de
Figura 5.28 Forma canónica donde es posible aplicar la función feedback() para evaluar ( )T z .
*( )E s
T
E(s) G(s)
Y(s)
H(s)
+ −
R(s)
−
T
5 - 34 Capítulo 5 – ANALISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO
muestreo. Sin embargo, el GFSO incluye el nodo ( )Y s que puede usarse para calcular
( )y t , lo cual permitiría obtener valores entre intervalos de muestreo. Esta señal ( )Y s se se
obtiene transmitiendo la señal *( )E s con una ganancia ( )G s . Como la señal *( )E s también
aparece en el GFSM, bastaría incluir una rama auxiliar con ganancia ( )G s para lograr el
valor de ( )Y s , tal como lo muestra en el ejemplo anterior.
Una solución más general consiste en efectuar la interconexión del GFSO y el GFSM para
obtener el modelo general del sistema [Kuo92], tal como se muestra en la figura 5.29, a
través del cual se pueden formulas las siguientes relaciones
*
*
( ) ( )( ) ( ) ( ) ( )
1 ( ) 1 ( )
G z G sY z R z Y s R s
GH z GH s= ⋅ = ⋅
+ + (5.54)
De este modo es posible obtener la respuesta discreta ( ) ( )Y z y kT↔ y la respuesta
continua ( ) ( )Y s y t↔ , esta última para evaluar la salida entre intervalos de muestreo.
En los siguientes ejemplos se aplicará el procedimiento propuesto, donde se demostrará que
es posible identificar la rama auxiliar en el GFSM para la evaluación de la respuesta ( )Y s .
Ejemplo 5.9: Determinar la respuesta en LC del siguiente sistema de control en cascada
1 2
1 1 5( ) , 0.1 ( ) , 0.1 ( )
( 1) 2 5G s ZOH con T s G s ZOH con T s H s
s s s s= = = = =
+ + +
Respuesta de un sistema de control de datos muestreados en cascada.
Figura 5.29 Modelo general del sistema de la figura 5.24
*( )GH s−
*( )E s
( )E s *( )E s
*( )G s
( )G s
( )H s−
*( )R s *( )Y s
( )Y s ( )R s
2( )E s
+ − −
2 *( )E s
T 2( )G s
( )Y s
( )H s
+
( )R s 1 *( )E s 1( )E s
1( )G s T
5.4 – MODELO DE SISTEMAS DE CONTROL DISCRETOS 5 - 35
Solución: Aplicando el procedimiento propuesto, obtenemos
1. GFSO a partir del DB.
Considerando la señal de entrada a cada bloque
Salidas: ( )Y s , 1( )E s , 2( )E s
Entradas: ( )R s , *1 ( )E s , *
2 ( )E s
2. Ecuaciones estándar.
Para las 3 salidas del sistema y a partir del GFSO, obtenemos
2 2
1 2 2
2 1 1 2 2
( ) *( )
( ) ( ) *( )
( ) *( ) *( )
Y s G E s
E s R s G E s
E s G E s G H E s
= ⋅= − ⋅= ⋅ − ⋅ ⋅
3. TEL de cada ecuación estándar.
Aplicando la regla de factores estrellados y no estrellados
2 2
*2 2 2
2 2 2 2 2
*( ) * *( )
*( ) ( ) * *( )
*( ) * *( ) * *( )
Y s G E s
E s R s G E s
E s G E s G H E s
= ⋅
= − ⋅
= ⋅ − ⋅
4. GFSM a partir de las ecuaciones estándar estrelladas.
La rama auxiliar se obtiene interpretando la primera ecuación no estrellada.
5. FGM para calcular la respuesta *( )Y s del sistema.
1
1
* *2* *
** *2 2
( ) ( )( ) ( )
1 ( ) ( ) ( )
G s G sY s R s
G s G s G H s
⋅=
+ ⋅ +
*2G H−
*2( )E s *
2G *( )R s *( )Y s
2G
( )Y s *2G−
*1G
*1( )E s
H−
1( )E s 2 *( )E s 2G
1−
( )Y s ( )R s 1 *( )E s 1G
2( )E s
5 - 36 Capítulo 5 – ANALISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO
Como todos los términos están estrellados podemos obtener el modelo
equivalente en el dominio-z, como
)()()(1
)()()(
221
21
zHGzGzG
zGzGzT
++=
Con ayuda del módulo scejem5p15.m evaluamos )(2 zHG :
)8187.0)(6065.0(
)7919.0(020.0)(2 −−
+=zz
zzHG
Asimismo, para 1( )G z y 2( )G z , obtenemos:
8187.0
0906.0)(
)9048.0)(1(
)9672.0(0048.0)( 21 −
=−−+=
zzG
zz
zzG
Sustituyendo estos valores en la expresión de ( )T z obtenemos:
)5049.0396.1)(9176.0915.1(
)6065.0)(9672.0(00043844.0)(
22 +−+−−+=
zzzz
zzzT
La figura 5.26 muestra la respuesta escalón del sistema de control en cascada. Utilizando la rama auxiliar del GFSM se puede obtener
** 1 2
** *1 2 2
( ) ( )( ) ( )
1 ( ) ( ) ( )
G s G sY s R s
G s G s G H s
⋅=+ ⋅ +
que puede usarse para calcular ( )y t mediante la transformada inversa de
Laplace. El calculo de ( )y t es complejo y será mostrado en el ejemplo 5.17. ------------------------- scejem5p15.m -------------------------
Tabla 5.2 Valores característicos de la respuesta e scalón
Modelo ptM ssT DCK
Continuo 1.1630 8.1718 1.0000
0.2T s= 1.2061 8.4000 1.0000
1T s= 1.3996 16.0000 1.0000
Los resultados de la tabla 5.2 demuestran el efecto del período de muestreo en el
comportamiento dinámico de un sistema de datos muestreados. Según esto, un valor grande
de T puede afectar considerablemente la estabilidad relativa del sistema (respuesta
transitoria). La ganancia DC (respuesta permanente) no se afecta, por la presencia del
retensor (ZOH). Como se demostrará en las secciones 5.6 y 5.7, en un caso más extremo
podría afectar la estabilidad absoluta del sistema.
Ecuación característica del sistema de datos muestr eados
En la sección 1.5 se introdujo el concepto de ecuación característica (EC) como la
expresión que permite calcular los polos del sistema de control en lazo cerrado, que para el
caso de la forma canónica de lazo cerrado del sistema continuo de la figura 1.25, se
presentó en (1.20) como: 1 ( ) ( ) 0G s H s+ = . Las raíces de esta ecuación determinan la forma
característica de la respuesta transitoria y se reconocen como los polos del sistema en lazo
cerrado. De modo similar, la ecuación característica del sistema de datos muestreados para la forma
canónica de la figura 5.16, tomada como el denominador de la FTLC (5.51), es
1 ( ) 0GH z+ = (5.55)
Para un caso más general de un sistema diferente a la forma canónica, la EC puede ser
obtenida a partir del gráfico de flujo de señales modificado (GFSM), que fue desarrollado
como parte del procedimiento propuesto en la sección 5.4 para evaluar la FT de lazo
cerrado del sistema. En efecto, aplicando la fórmula de ganancia de Mason al GFSM
1
( ) ( )
( )( )
N
k k
k
T z z
T zz
=
⋅ ∆=
∆
∑ (5.56)
5.5 – RESPUESTA DINAMICA DEL SISTEMA DE CONTROL DISCRETO 5 - 47
Según (5.56) el denominador de ( )T z es el determinante del sistema: ( )z∆ . Luego la
ecuación característica equivalente en el dominio-z de un sistema arbitrario de datos
muestreados, puede evaluarse como ( ) 0z∆ = (5.57)
En el capítulo 7 se podrá identificar una tercera forma de obtener la EC de un sistema de
control de datos muestreados, a partir de su modelo de estado. Los polos del sistema serán
asociados con los autovalores de la matriz de estado A del sistema. De este modo, (5.55) y (5.57) son expresiones alternas para evaluar la EC del sistema,
dependiendo del modelo utilizado para su representación. El procedimiento presentado en
la sección 5.4 es una herramienta válida para obtener la EC de un sistema de datos
muestreados, a partir de (5.57). Los siguientes ejemplos muestran su aplicación.
Ejemplo 5.13: Obtener la EC del sistema de datos muestreados del ejemplo 5.15, cuyo
GFSM se presenta a continuación.
Solución: Aplicando la FGM se puede obtener el determinante del sistema como
1 2 2( ) 1 ( ) ( ) ( )z G z G z G H z∆ = + +
Aplicando (5.57) obtenemos la EC como
1 2 2( ) 1 ( ) ( ) ( ) 0Q z G z G z G H z= + + =
que es el denominador de ( )T z que se obtuvo en el ejemplo 5.15. Para los
valores dados de 1( )G s , 2( )G s y ( )H s
0.0048( 0.9672) 0.0906 0.020( 0.7919)( ) 1 0
( 1)( 0.9048) 0.8187 ( 0.6065)( 0.8187)
z zQ z
z z z z z
+ += + × + =− − − − −
Simplificando, obtenemos
4 3 2( ) 3.3102 4.0946 2.2474 0.4633 0Q z z z z z= − + − + =
cuyas raíces son 1,2 0.9579 2.14ºp = ±∡ , 3,4 0.7106 10.87ºp = ±∡ . Como todos
los polos están dentro del círculo unitario, el sistema es absolutamente estable.
Ecuación característica a partir del GFSM de un modelo donde existe T(z).
*2G H−
*2( )E s *
2G *( )R s *( )Y s
*2G−
*1G
*1( )E s
5 - 48 Capítulo 5 – ANALISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO
En el ejemplo 5.22 se utilizó la función damp() del Toolbox de Control de MATLAB ® para
verificar el cálculo de los parámetros ζ y ωn de la FT del modelo continuo ( )T s y del
modelo discreto ( )T z . Su sintaxis general es EDU» [wn,zeta,polos]=damp(sis)
donde sis es la FT del modelo continuo o discreto, creado como un objeto LIT. Los
parámetros ζ y ωn se evalúan para cada polo complejo conjugado, a partir del triángulo
característico. Si el modelo es discreto, primero se calculan los polos equivalentes en el
5.5 – RESPUESTA DINAMICA DEL SISTEMA DE CONTROL DISCRETO 5 - 53
dominio-s a partir de (5.60). Si el polo es real devuelve ζ=1 y un valor de ωn que no tiene
significado. La función se puede usar también para obtener los polos del sistema. Los dos ejemplos anteriores permiten verificar el efecto del período de muestreo en el
desarrollo del modelo discreto equivalente de un sistema de datos muestreados. Utilizando
la función especial valrespt() presentada en la sección 2.2, se desarrolló la tabla 5.3
tomando como base el sistema del ejemplo 5.22. Los resultados mostrados son suficientes
para demostrar el efecto del período de muestreo T en los valores característicos de la
respuesta escalón y en los parámetros de la FT del modelo discreto.
Tabla 5.3 - Efecto del período de muestreo T en sistema de 2º orden
Modelo ptM ssT ζ
nω
Continuo 1.1630 8.1718 0.5000 1.0000
0.1T s= 1.1837 8.3000 0.4746 0.9991
0.5T s= 1.2864 11.0000 0.3695 0.9780
1T s= 1.3996 16.0000 0.2494 0.9197
Resolución del modelo discreto: efecto de T
El resultado de la tabla 5.3, se debe a que la correlación entre los polos del modelo continuo
( )T s y de su equivalente discreto ( )T z está influenciada por el período de muestreo T. Una
forma más de evaluar el efecto del período de muestreo en el modelo discreto equivalente,
consiste en determinar la resolución del modelo discreto, expresada en términos del número
de muestras ( )mN por unidad referencial de tiempo (URT). El valor esperado mN puede
usarse como un criterio práctico para la selección del período de muestreo. Si el modelo continuo ( )T s tiene un polo real en 0 1/s = τ , se puede establecer como
URT la constante de tiempo ( )τ , con lo cual el número esperado de muestras sería de
/mN T= τ . El polo de ( )T z aparece en /0
Tz e− τ= . Despejando la relación T/τ , y
considerando el caso de 0 0z < obtenemos
0
1/
( )mN T
ln z= τ = (5.61)
Si el polo de ( )T s es complejo conjugado 0 ds j= −α ± ω , se puede establecer como URT el
período dT asociado con la frecuencia natural amortiguada 2 /d dTω = π y el número
esperado de muestra sería de /m dN T T= . En este caso el polo continuo de transforma en
( )0
d dj T T j Tz e e e r−α± ω −α ± ω= = = ∠ ± θ (5.62)
5 - 54 Capítulo 5 – ANALISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO
donde Tr e−α= y dTθ = ω . Considerando que 2 /d Tω = π , obtenemos
2d
m
TN
T
π= =θ
(5.63)
De este modo el número de muestras logrado, depende de la posición final del polo 0z en el
plano-z. El criterio práctico se basa en seleccionar el período de muestreo T, de modo que
el número mínimo de muestras logradas sea de (8 10)m dN a muestras por o T= τ (5.64)
Este criterio será utilizado en el capítulo 7 para la selección del período de muestro, como
fase preliminar del diseño de controlador. A partir de (5.62) se pueden desarrollar las
siguientes expresiones
21nT
nr e T−ζω= θ = ω − ζ ⋅ (5.65)
la cuales permiten obtener gráficas características de ζ y ωn constantes [ReySoto06], que
son de utilidad en aplicaciones prácticas de análisis y diseño del sistema de control en
tiempo discreto. El Toolbox de Control de MATLAB ® incorpora las funciones sgrid() y
zgrid() que se pueden utilizar con este propósito y fueron utilizadas para generar la
gráfica de la figura 5.29.
A partir de (3.22) se puede verificar que cada línea guarda relación directa con el número
de muestras (Nm). En efecto, para la primera línea θ=π/6 radianes, que equivale a
dm TmuestrasN /126/
22 =π
π=θπ=
Las siguientes líneas corresponden a , , 6 4mN = … ; este resultado es importante, porque
muestra que para lograr un valor alto de mN que garantice una adecuada resolución del
Figura 5.29 Curvas características de ζ y nω constantes.
5.5 – RESPUESTA DINAMICA DEL SISTEMA DE CONTROL DISCRETO 5 - 55
modelo discreto, el polo complejo conjugado dominante de ( )T z deberá estar ubicado en el
primer cuadrante del plano-z.
Desplazamiento de polos de ( )T z : efecto de T
De acuerdo con la figura 5.29, la posición final del polo del modelo discreto equivalente
( )T z depende del valor seleccionado para el período de muestro T. De acuerdo con esto, si
se cambia el valor de T debe ocurrir un desplazamiento del polo de T(z) en el plano-z,
modificando substancialmente el comportamiento dinámico del sistema discreto. Para evaluar el efecto de desplazamiento que sufre el polo del modelo discreto equivalente ( )T z en el plano-z al cambiar el período de muestreo T, consideramos dos casos:
Caso 1: Polo real de ( )T s en 0 6s = − , del ejemplo 5.21.
Caso 2: Polo complejo conjugado de ( )T s en 0 0.5 0.8660s j= − ± , del ejemplo 5.22.
La tabla 5.4 presenta un resumen de los resultados, donde se observa el efecto del período
de muestro T en la posición final de los polos del modelo discreto equivalente ( )T z .
Tabla 5.4 – Efecto de T en la resolución y en la posición final de los pol os de T(z)
T polo real 0 6s = − polo complejo 0 0.5 0.8660s j= − ±
El modelo del controlador discreto PID será desarrollado en el capítulo 7 y puede diferir de la
forma del ejemplo 5.24, dependiendo del algoritmo utilizado en su desarrollo [Ogata95]. El
siguiente ejemplo muestra el uso de la función especial errores() en el cálculo del error
estacionario de un sistema de control de datos muestreados con realimentación dinámica.
Ejemplo 5.19: Calcular el error estacionario del sistema mostrado a continuación.
Solución: Se observa que este sistema es la versión de datos muestreados del ejemplo
2.12. Como el sistema tiene realimentación dinámica, no se podrían aplicar las
expresiones de la tabla 5.5. Para comprobar lo anterior, obtenemos ( )GH z
1
2
2( 2) 0.0088105( 0.8756)( 0.8187)( ) (1 )
( 1)( 5) ( 1)( 0.6065)( 0.9048)
s z zGH z z
s s s z z z
− + + −= − = + + − − − Z
Como el sistema es tipo 1N = , aplicando la tabla 5.5 obtenemos
Error estacionario en sistema de datos muestreados con realimentación dinámica.
1
( 1)s s +
Y(s)
+
R(s) E*(s) E(s)
−
ZOH
T=0.5
2( 2)
5
s
s
++
5 - 62 Capítulo 5 – ANALISIS DE SISTEMAS DE CONTROL EN TIEMPO DISCRETO
1
0.0088105( 0.8756)( 0.8187)( 1) 0.080
( 1)( 0.6065)( 0.9048)v
z
z zK z
z z z
∗
=
+ −= − =− − −
Luego el error estacionario discreto de velocidad sería 1/0.08 12.5sve∗ = = . Sin
embargo, como la realimentación es unitaria es necesario calcular el error de
velocidad usando el teorema del valor final a *( )e t del diagrama anterior. Sin embargo, como el error estacionario no se ve afectado por el período de
muestreo T, podemos evaluarlo a partir del modelo continuo:
[ ]( ) ( ) ( ) 1 ( )He t E s R s k T s↔ = −
donde, según (2.48), 0
( ) 0.8H sk H s == = . Para ( )T s obtenemos
3 2
( ) 5( )
1 ( ) ( ) 6 7 4
G s sT s
G s H s s s s
+= =+ + + +
Por lo tanto, para 2( ) 1/R s s=
2
2 3 2 3 2
1 5 6 6.2( ) 1 0.8
6 7 4 ( 6 7 4)
s s sE s
s s s s s s s s
+ + + = − × = + + + + + +
Aplicando el teorema del valor final
2
3 20
0
6 6.2( ) 1.55
6 7 4sv s
s
s se sE s
s s s==
+ += = =+ + +
Se demuestra así que en este caso no se pueden aplicar las expresiones de la
tabla 5.5. Otra alternativa es aplicar la función especial errores() que se
fundamenta en el algoritmo especial desarrollado en la sección 2.3 ------------------------ SCEjem5p25.m --------------------------
%modelo discreto de la FTLA: GH(z) Gs=zpk([],[0 -1],1), Hs=zpk(-2,-5,2), T=0.1; GHz=c2 d(Gs*Hs,T)
%error estacionario de velocidad utilizando la tabl a 5.5 Kp=dcgain(GHz), esp=1/(1+Kp) Kv=dcgain(minreal(tf([1 -1],1,T)*GHz)), esv=1/Kv