Sistemas de Control Digital: University of Illinois at Champaign
Sistemas de Control Digital: University of Illinois at Champaign
BENJAMIN C. KUO
University of Illinois at Champaign
PRIMERA EDICIÓN
MÉXICO, 1997
Traducción:
Gerardo Urbina Medal
Departamento de Ingeniería Eléctrica
UAM, Iztapalapa
Revisión técnica:
Gerardo René Espinoza Pérez
Investigador asociado,
Instituto de Ingeniería, UNAM
Diseño de portada:
Estrategia Visual
ISBN 968-26-1292-6
Impreso en México
Printed in Mexico
Este libro trata la teoría y aplicaciones de los sistemas de control digital y constituye
una extensa revisión de la primera edición, publicada en 1980. Desde 1963 el autor
ha escrito varios libros sobre sistemas de control digital y de datos muestreados.
El primero de ellos, Analysis and Synthesis of Sample-Data Control Systems, fue
publicado en 1963 por Prentice-Hall. El segundo, Discrete-Data Control Systems,
apareció en 1970 publicado por la misma casa editorial y presentó un enfoque
totalmente nuevo con respecto al ofrecido en la edición de 1963, ya que en él se hizo
hincapié en variables de estado y teoría de control moderna. La primera edición de
Digital Control Systems fue publicada por Holt, Reinhart and Winston en 1980, y es
una revisión completa del libro editado en 1970.
CONTENIDO SINTÉTICO
El capítulo 1 proporciona una introducción general ala materia de los sistemas de control
digital, al presentar varios ejemplos comunes de sistemas, de control de este tipo. El
capítulo 2 abarca la conversión de señales, el procesamiento de las mismas y la operación
de muestreo. El contenido de este capítulo sirve como motivación del tratamiento
matemático y analítico de datos y señales digitales en sistemas de control. El capítulo
3 presenta las transformadas z y z modificada, incluyendo los sistemas con muestreo no
uniforme y aquellos donde las señales se muestrean con frecuencias distintas. El capítulo
PREFACIO VII
Cada capítulo comienza con una lista de palabras clave y temas que dan un panorama
de los plintos más importantes estudiados en él. Hasta donde es posible, cada uno
comienza con un introducción que sirve para motivar los temas tratados en el capítulo,
además de indicar los conocimientos previos necesarios para entender el material. En
todo el libro aparecen muchos ejemplos ilustrativos, la mayor parte tomados de
sistemas prácticos. Al término de cada capítulo, comenzando con el 2, se encuentran
muchos ejercicios y problemas.
VIII SISTEMAS DE CONTROL DIGITAL
RECONOCIMIENTOS
B.C. Kuo
Champaign, Illinois
Contenido
Capítulo 1 Introducción 1
1-1 Introducción 1
1-1-1 Elementos básicos de un sistema de control de datos discretos 2
1-1-2 Ventajas de los sistemas de control de datos discretos 3
1-2 Ejemplos de sistemas de control digital y discretos 6
1-2-1 Sistema de control simplificado de un solo eje para un piloto
automático 6
1-2-2 Sistema de regulación de laminado controlado por computadora
digital 8
1-2-3 Controlador digital para una turbina y generador 8
1-2-4 Sistema de control de un motor de pasos 9
1-2-5 Sistemas controlados por microprocesador 10
1-2-6 Modelo discreto de un problema de pago de intereses 11
B ibl iografia 12
lx
Capítulo 3 Transformada z 75
CONTENIDO XI
CONTENIDO XIII
CONTENIDO XV
Apéndice A Números de punto fijo y de punto flotante 725
A-1 Representación de números en punto fijo 725
A-2 Representación de números en punto flotante 728
Apéndice Modelado matemático del muestreo con la Integral
de convolución 731
B-1 Expresiones para la salida del muestreador de ancho de pulso finito
Fp* (s) 731
B-2 Expresiones para la salida del muestreador ideal Fe (s) 736
B-3 Expresiones alternas para la transformada z modificada 737
Sistemas de datos discretos • sistemas de control digital • datos muestreados • ventajas del control
digital y discreto • ejemplos de sistemas de control de datos discretos y digitales
1-1 INTRODUCCIÓN
1
2 SISTEMAS DE CONTROL DIGITAL
discretos en un sentido amplio para describir todos los sistemas donde aparecen
señales digitales o muestreadas. Por otra parte, los métodos de análisis y diseño son
esencialmente los mismos, sin importar que el sistema contenga datos muestreados o
codificados de manera digital.
La figura 1-1 presenta los elementos básicos de un sistema común de control de lazo
cerrado con datos muestreados. El muestreador simplemente representa un dispositivo
u operación cuya salida es un tren de pulsos. Ninguna clase de información se
transmite entre dos pulsos consecutivos. La figura 1-2 ilustra la entrada y salida
comunes de un muestreador. La señal de entrada e(t) continua queda muestreada por
el dispositivo, y la salida es una secuencia de pulsos. En el caso ilustrado, se supone
que la frecuencia de muestreo es uniforme. Las magnitudes que tienen los pulsos en
los instantes de muestreo representan los valores de la señal de entrada e(t) justo en
esos momentos. En general, el esquema de muestreo tiene muchas variaciones:
algunos son periódicos, de frecuencia cíclica, con varias frecuencias, aleatorios y
modulados por el ancho de pulso. El tipo más común de muestreo encontrado en
sistemas prácticos es el muestreo uniforme con una sola frecuencia, tal como se
muestra en la figura 1-2.
El filtro ubicado entre el muestreador y el proceso controlado se utiliza para
suavizar la salida del primero debido a que la mayor parte de los procesos controlados,
como aquellos donde intervienen motores convencionales de ca o cd, están diseñados
y construidos para recibir señales analógicas.
La figura 1-3 presentí el diagrama de bloques de un sistema de control digital
representativo. La existencia de señales codificadas de manera digital, como las
señales codificadas en binario, en ciertas partes del sistema requiere convertidores
digital-analógicos (D/A) y analógico-digitales (A/D). El bloque de la figura 1-3 que
e(t)
(a) t (b)
Figura 1-2. a) Entrada continua al muestreador. b) Salida discreta del muestreador.
INTRODUCCIÓN 3
Código
digital de Computadora Convertidor, Proceso Salida
entrada digital digital-analógico controlado
Convertidor
analógico-digital
Barra relevadora
Galvanómetro
Leva accionada
Hacia por un motor
la carga
Polea
accionada por Contactos del relevador
un motor (normalmente cerrados)
Varilla de contacto
Termómetro de mercurio
IZZ
Bobina de Homo
calentamiento
Corriente de
Temperatura del horno h(t) calentamiento El ancho del pulso es proporcional a la amplitud
del horno de h(t) en 1 = nT
4(0
t
kT
fi 11H
(k+1)T (k4-2)T (k+3)T (k+4)T
►t
Proceso Salida
Controlador controlado
que plantear la siguiente pregunta: ¿por qué los datos muestreados y el control digital?
En otras palabras, ¿cuáles son las ventajas y características de los datos muestreados
y del control digital? Primero es necesario reconocer que el muestreo es una caracte-
rística inherente de muchos sistemas físicos, y que el comportamiento de éstos puede
describirse mediante datos discretos o modelos digitales; por ejemplo, en un sistema
de radar, las señales transmitidas y recibidas por éste tienen la forma de trenes de
pulsos. La operación de rastreo del radar realiza la función de muestreador, que
convierte en datos muestreados la información relacionada con el azimut y la
elevación. En el mundo existen muchos fenómenos y sistemas sociales, económicos
o biológicos cuya dinámica puede representarse con modelos de datos discretos.
Para abordar las ventajas y características de los sistemas de datos discretos,
primero es necesario examinar las de los sistemas de control de datos continuos, los
cuales tienen muchos años de existencia. La figura 1-7 muestra un diagrama de blo-
ques que contiene los elementos básicos de un sistema de control de datos continuos
formado por un solo lazo. En general, el controlador es un circuito electrónico que
trabaja con una señal analógica y cuya salida es otra señal del mismo tipo. La ventaja
del controlador analógico es que el sistema trabaja en tiempo real y puede tener un
ancho de banda muy grande. Esto es equivalente a tener una frecuencia de muestreo
infinita, de modo que el efecto del controlador siempre está presente. Las desventajas
del controlador analógico es que sus elementos se encuentran montados de una manera
rígida, de modo que las características del controlador están fijas de antemano, lo que
dificulta modificar el diseño; por otra parte, el envejecimiento de los componentes y
la sensibilidad de éstos a cambios en el medio pueden llegar a ser muy importantes.
Asimismo, los componentes analógicos exhiben una susceptibilidad mayor a proble-
mas de ruido. Algunas de las ventajas del control digital son las siguientes:
Los componentes digitales tienen menor susceptibilidad al envejecimiento y a
las variaciones de las condiciones ambientales.
Los componentes digitales son menos sensibles al ruido y las perturbaciones.
Los procesadores digitales tienen un tamaño y peso menor. Pueden fabricarse MP
y PDS de un solo circuito integrado muy versátiles y de gran capacidad para
aplicaciones de control.
Los costos de los MP y PDS son cada vez menores.
Los MP y PDS permiten una flexibilidad mayor en la programación; un cambio
en el diseño no requiere cambios en el hardware.
Son más confiables.
7. Proporcionan una mejor sensibilidad a variaciones en los parámetros.
El galvanómetro de barra relevadora de la figura 1-4 es un ejemplo de la manera
en que el muestreo puede mejorar la sensibilidad. En este caso se tiene la amplifica-
ción, para fines de control, de una señal pequeña mediante la operación de muestreo.
6 SISTEMAS DE CONTROL DIGITAL
INTRODUCCIÓN 7
Comando Altitud
de altitud Controlador Dinámica e e
Prefiltro analógico del fuselaje
Giróscopo
Transductor
de posición
(a)
Comando
de altitud Controla - Dinámica
Prefiltro D/A —..«del fuselaje -
dor digital
Giroscopio
lo
Transductor
de posición
(b)
Figura 1-8. a) Sistema de control simplificado de un solo eje para un piloto automático con
datos analógicos. b) Sistema de control simplificado de un solo eje para un piloto automático
con datos digitales.
pulso hasta que llegue otra muestra. La figura 1-9 ilustra la situación en que se tienen
dos muestreadores con diferentes periodos de muestreo Ti y T2. En general, si la
rapidez de variación de la señal en un lazo es mucho menor que la de otro lazo,
entonces puede aumentarse la duración del periodo de muestreo del lazo más lento.
El sistema de la figura 1-9 tiene muestreadores que trabajan con frecuencias de
muestreo diferentes. Este tipo de sistemas se conoce como sistemas de datos con
muestreo múltiple.
Una de las ventajas del uso del muestreo y del muestreo múltiple es que pueden
emplearse ciertos componentes de alto costo sobre una base de tiempo compartido.
Comando
de altitud
Prefiltro _91 dorControla-
digital D/A
Dinámica
lel fuselaje
•••••••1111.
-•
Retenedor Giroscopio
Tt
Muestreador
Transductor
Retenedor 1--
T2 de posición
Figura 1-9. Sistema de control digital simplificado de un solo eje para un piloto automático con
muestreo múltiple.
8 SISTEMAS DE CONTROL DIGITAL
Rodillos
Hoja
o
Control de • ransclu Sensor de
Rodillo de espesor Control Sensor espesor Rodillo de
rebobinado de tensión de tensión desenrollado
Velocidad 111.•n•
Comandos Tensión
de entrada Computadora digital TCR de visualización
Espesor
Sensor de
ND espesor
Del
calentador
(Control de velocidad) 4
Transductor
Comando H D/A Válvula -a. Turbina ,-410 Variables
-41.electrohidráulico
codificado de salida
Mini - Generador
computadora
--41 D/A Excitador
(Control de voltaje)
Adquisición
de datos
Variables
de salida Transducto- Multiplexor Salida digital hacia
res digitales digital la computadora
En ocasiones se tienen sistemas en que todos los elementos son digitales, de modo
que no es necesario usar convertidores A/D y D/A a fin de acoplar las señales. La
figura 1-15 muestra uno de estos sistemas, el cual se utiliza para el control de la cabeza
de lectura/escritura de una unidad de disco de computadora. El elemento principal
empleado por este sistema es un motor de pasos cuya respuesta está dictada por una
serie de pulsos. Por cada pulso de entrada, el motor responde con un desplazamiento
fijo. Dadas estas características, el sistema puede considerarse como un sistema de
datos completamente digital.
Variables Salida
de salida Transductor Multiplexor Muestreador digital hacia
analógico y retenedor A/D
analógico la computadora
Codificador
Computadora Circuito Motor digital de posición
digital ••11> Secuencia de mando de pasos
Cabeza
de lectura
escritura
Figura 1-15. Sistema de control del motor de pasos para el posicionamiento de la cabeza de
lectura-escritura de una unidad de disco.
GALIL
Circuito integrado
GL-1200
Dirección Comando Posición
de posición
Posición Codificador
del motor incremental
Figura 1-16. Sistema de control de posición con el controlador de un solo circuito integrado
GL 1200.
INTRODUCCIÓN 11
Procesador Huésped
huésped E/S
Figura 1-17. Sistema de control de posición con el microprocesador de un solo circuito integrado
LM628.
mejora mayor en el desempeño del sistema. La razón de esto es que los PDS tienen
la capacidad de realizar cálculos con una velocidad mayor que la de los microproce-
sadores. Como consecuencia de lo anterior, los PDS pueden programarse con la
finalidad de dotar de características de filtrado con pendientes de corte muy grandes,
imposibles de obtener con los microprocesadores convencionales.
Uk =U O p(k + 1)
• Retraso p k) •
1+R •
Figura 1-18. Modelo de datos discretos para el problema de pago de intereses.
12 SISTEMAS DE CONTROL DIGITAL
donde p(k + 1) denota la cantidad adeudada al término del (k + 1)-ésimo periodo. Las
condiciones de frontera son p(0), la cual está dada, y p(N) = O. La ecuación de
diferencias 1-1 puede resolverse para U mediante el método recursivo o la transfor-
mada z (Cap. 3). Esta ecuación también se conoce como ecuación de diferencias de
primer orden. Por otra parte, el modelo del sistema de datos discretos de la figura 1-18
está dado en forma de diagrama de bloques.
Bibliografía
Conversión y
procesamiento de
señales
2-1 INTRODUCCIÓN
13
14 SISTEMAS DE CONTROL DIGITAL
Convertidor D/A
Un convertidor digital-analógico decodifica una señal digital de entrada. La salida
del convertidor D/A es una señal analógica, usualmente en forma de un voltaje o
corriente. n
Convertidor A/D
Un convertidor analógico-digital es un dispositivo que convierte una señal analógica
en una señal codificada digitalmente. n
CONVERSIÓN Y PROCESAMIENTO DE SEÑALES 15
MultIpl«or
El multiplexor se emplea para acoplar señales que provienen de varias fuentes, de
modo que las pueda procesar el mismo procesador o canal de comunicación. La
finalidad de este dispositivo es compartir el tiempo del procesador entre todas las
señales que llegan a él. n
7/8
3/4
5/8
1/2
"
1/4
LSB
1/8
I 1 1 1 I 1 a
0.000 0.010 0.100 0.110
0.001 0.011 0.101 0.111
Palabra digital (fracción)
- Relación entre el código binario fraccionario y números decimales.
Rozara 2 1.
En los procesos de conversión D/A y A/D, los bits más significativo (MSB, por sus
siglas en inglés) y menos significativo (LSB, por sus siglas en inglés), así como el
peso de cada uno en la palabra, tienen una gran importancia en la comprensión del
proceso de conversión. Los convertidores D/A y A/D basados en el código binario
natural usan un código fracciona'. Según se muestra en la tabla A-2 del apéndice A,
para un código fracciona' de 3 bits, el MSB tiene un peso igual a la mitad de la escala
completa (1/2 FS, por sus siglas en inglés); el segundo bit, un peso de 1/4 FS, mientras
que el último, el LSB, un peso de 1/8 FS. En general, para un código fracciona' de n
bits, el MSB sigue teniendo un peso igual a 1/2 FS y el LSB, igual a FS. Ahora bien,
sin importar cuál sea el entero o código fracciona' utilizado, una palabra de n bits
siempre define 2" estados distintos, con lo cual proporciona una resolución de una
parte en 2"; por ejemplo, en la figura 2-1 la resolución es 2-3 = 1/8 para un código
fraccional de 3 bits. En este caso, el LSB es 1/8 y FS = 1. La resolución puede
mejorarse con más bits. Si se utiliza el código fracciona' con una palabra de cuatro
bits y se mantiene la misma escala completa, entonces el valor del LSB es 2 -4, o 1/16,
y la resolución mejora a una parte en 16. Con este cambio, el diagrama de la figura
2-1 tendría 16 niveles distintos. Nótese que, en general, el código digital utilizado no
tiene correspondencia con la señal analógica y viceversa; por consiguiente, los
números analógicos representados en la escala vertical de la figura 2-1 deben consi-
derarse como una parte de la escala completa. También se aplica el mismo principio
al código entero. Siempre debe tenerse en mente que para mejorar la resolución
mediante el incremento del número de bits, es necesario conservar la misma escala
completa para la señal analógica o digital.
La operación de conversión de una señal analógica en un número digital se conoce
como conversión A/D. Si el número de bits en la palabra digital es finito, como es el
caso de los sistemas reales, entonces el proceso de conversión A/D tiene una resolu-
ción finita; dado que la salida digital sólo cuenta con un número finito de niveles, el
número analógico queda redondeado o truncado en un número digital. Este procedi-
CONVERSIÓN Y PROCESAMIENTO DE SEÑALES 17
Señal codificada
digitalmente D
7q— FS-LSB = g FS
6q— FS
5q— TiFS
447— IPS
39— in
-§-FS 2q :$-FS
2 2
A2 di2 A2 2
q
q 3:7 4:7 54q i
1 1 1 1 I 1 nab
tlQ 13q
_'2 31 :71
Señal analógica A
q = nivel de cuantización
-4q
-5q
--6q
--7q
(a)
Figura 2-2. (a) Relación entrada-salida de un A/D de tres bits con cuantización por redondeo.
Señal codificada
digitalmente D
7 q — FS-LSB = FS
6 9— 71 FS
5q— FS
1
49— Fs
3
39 — "g FS
8 Fs 24- 1 FS
(/
-7q -6q -5q -4q -3q -2q -q /
1 . 1 if eitftie e/ I I 1 1
1 e 1
/
-q
q 2 9 39 4q 5 9 69 7q
/ Señal analógica
'29
-39 q nivel de cuantización
--4q
— -5q
--6q
— -7 q
AIA/A44/
PrIlr/P9rI
Error de cuantización
(b)
Figura 2-2. (b) Relación entrada-salida de un A/D de tres bits con cuantización por truncamiento.
Tabla 2-1 Relación entre números analógicos y digitales para una conversión
analógico-digital de tres bits (por redondeo)
Tabla 2-2 Relación entre números analógicos y digitales para una conversión
analógico-digital (por truncación)
MSB LSB
(x 4q) (x 2q) (x q)
jA(S q O O O O
q SIA <2 0 0 1 1/8 FS = q =LSB
2q51A <3q 0 1 0 1/4 FS - 2q
3q SIA <4q 0 1 1 3/8 FS - 3q
4qSIA <5q 1 0 0 1/2 FS ,- 4q
5qSIA <6q 1 0 1 5/8 FS = 5q
6qS1A <7q 1 1 0 3/4 FS = 6q
7q1lill< 00 1 1 1 7/8 FS - FS -LSB - 7q
Nótese que las cuantizaciones de las figuras 2-2(a) y 2-2(b) no son relaciones uno
a uno. El parámetro q, que es igual al valor del bit menos significativo, se conoce como
nivel de cuantización. Para el caso de la palabra de tres bits, el LSB es igual a 1/8 FS,
como se indica en las tablas 2-1 y 2-2. En las figuras 2-2(a) y 2-2(b), las líneas punteadas
representarían la salida ideal si no existiera el fenómeno de la cuantización. La diferencia
entre la línea recta y la característica de cuantización se conoce como error de cuanti-
por tanto, el valor máximo de la señal analógica que el cuantizador puede convertir
sin exceder el error de cuantización de ±0.625 es
FS + = 9.375 V (2-2)
8
Reloj
o+
Interruptor
es(t) eo (t)
Figura 2 -3. Circuito sencillo que ilustra el principio de operación del muestreador-retenedor.
carga con una constante de tiempo R.,C hasta alcanzar el valor de la señal muestreada.
Por otra parte, el funcionamiento del muestreador no es instantáneo, ya que le lleva cierto
tiempo responder a los comandos de muestreo y retención.
En la práctica el funcionamiento del S/H puede tener muchos errores e imperfec-
ciones, y la salida del dispositivo puede desviarse de manera considerable de la forma
de onda ideal de la figura 2-4.
La figura 2-5 ilustra una señal de entrada es(t) y la correspondiente señal de salida
del S/H. En general, la salida de un S/H se caracteriza por varias fuentes de retrasos
y de retención imperfecta durante el modo de retención. Estas características se
ilustran en la figura 2-5 y a continuación se proporcionan sus defmiciones.
es(t)
Seguimiento
e0(1)
e o(0
p p=0
t
Pulso de reloj
Muestreo
Retención
Figura 2-4. Señales simplificadas en el proceso de muestreo y retención; R, 0.
22 SISTEMAS DE CONTROL DIGITAL
r Atenuación
Señal de
entrada es(t)
Salida
ideal eo(t)
—*1 T
Muestreo Retención
Interruptor
t 'abierto
Comando Interruptor Comando
de muestreo cerrado de retención
Figura 2-5. Señales de entrada y salida del S/H con retrasos finitos.
L 0+
Entrada V Muestreador Salida
analógica e3(t) Amplificador Amplificador analógica eo(t)
de aislamiento de aislamiento
de la entrada de la salida _2_
Comando de muestreo
y retención
o T 2T 3T
Duración Periodo
del muestreo de muestreo
(seguimiento) F.o— Periodo Retención —.I
p de muestreo T
Figura 2-7. Señales de entrada y salida del dispositivo S/11 con frecuencia de muestreo uniforme.
24 SISTEMAS DE CONTROL DIGITAL
e d(t) e-Tds
7S Retenedor -4» Filtro Salida éo(t)
Entrada T (p)
Muestreador Retraso en el tiempo
Figura 2-8. Diagrama de bloques aproximado del comportamiento del dispositivo de muestreo
y retención.
ideal del S/H para una entrada analógica dada. El tiempo transcurrido entre las órdenes
de muestreo se conoce como periodo de muestreo T.
ideal
e5(t)
t
T 2T 3T 4T 5T
Comando de muestreo
La figura 2.10 presenta las formas de onda características para la entrada y salida de un
dispositivo S/H ideal. Se mostrará más adelante que el modelo del muestreador ideal
establece un modelo matemático conveniente para la operación S/H.
—o Corriente
—NsAr--
R1
Entradas —II'. Circuitería • Interruptores Red
digitales —11' lógica 11"
• electrónicos
1n•nnn•n•
resistiva
E09o deVoltaje
salida
t _L.
Referencia
de precisión
Voltaje
de referencia
Figura 2-12. Convertidor D/A de tres bits con red resistiva ponderada.
de cada bit de entrada. La figura 2-12 presenta un convertidor D/A sencillo de tres
bits. Los valores de los resistores que forman la red de adición están ponderados de
manera binaria. Cada uno está conectado mediante un interruptor electrónico al
voltaje de referencia o a tierra. Cuando aparece un uno binario en la circuitería lógica
de control del interruptor, éste se cierra y conecta el resistor al voltaje de referencia.
Por otra parte, un valor binario cero lo conecta a tierra. Para un amplificador
operacional de alta ganancia, la impedancia de entrada es tan grande que el voltaje en
el punto de suma es prácticamente cero; por tanto, si el interruptor que corresponde a
la rama del MSB está conectado al voltaje de referencia —Er, y los otros dos
interruptores lo están a tierra, lo que corresponde a la palabra digital 100, el voltaje
de salida E0 es
E, R f l f (2-4)
Dado que 1f = 1i y que Ji = &IR, entonces
R E,
= f (2-5)
El valor del LSB que corresponde a la palabra digital 001 es RfEr /4R. Es así como
la salida a escala completa es 8RfEr /4R = 2RfEr /R. La tabla 2-3 contiene los voltajes
de salida del proceso de conversión D/A correspondientes a todas las palabras de tres
bits para una FS de 10 V.
CONVERSIÓN Y PROCESAMIENTO DE SEÑALES 27
Tabla 2-3 Voltajes de salida del convertidor D/A de la figura 2-12 correspondientes a todas
las entradas digitales posibles
RfEr FS 2.50
010
2R 4
3RJEr 3
011 — FS 3.75
4R 8
RfEr 1
100 — FS 5.00
R 2
5RfEr 5
101 — FS 6.25
4R 8
6RJEr 3
110 — FS 7.50
4R 4
r 7
111 FS — LSB FS 8.75
7 4R 8
donde ao, al, ..., an..1 son uno o cero, lo que depende de la palabra digital de entrada
al convertidor.
Dado que un dispositivo D/A convierte una señal digital en una analógica con
una magnitud equivalente al valor de la palabra digital, entonces desde un punto de
vista funcional, el dispositivo puede considerarse como formado por un decodificador
Entrada Salida
digital analógica
D/A •
Entrada Salida
digital analógica
Decodificador
Entrada Salida
analógica digital
Entrada Salida
analógica digital
S/H Q n•n•••n• Decodificador
Cuantizador
Figura 2-14. Representación en diagrama de bloques de un convertidor A/D.
CONVERSIÓN Y PROCESAMIENTO DE SEÑALES 29
Cada uno de estos tipos tiene sus propias ventajas y limitaciones. Cada uno es útil
para una clase específica de aplicaciones de acuerdo con la velocidad de conver-
sión, costo, exactitud y tamaño. En el mercado pueden hallarse unidades A/D des-
de seis hasta 16 bits. La exactitud de conversión más común es de 0.01% de la FS
±1/2 LSB.
Con la finalidad de ilustrar el proceso básico de conversión A/D, a continuación
se describe de manera breve el convertidor A/D por aproximaciones sucesivas. La
Salida de reloj
Salida en
t Estado formato serie
Reloj
•••
Referencia
analógica EVA
fr
Entrada
analógica ° Comparador
Figuro 2-15. Diagrama de bloques simplificado de un convertidor A/D por aproximaciones
sucesivas.
30 SISTEMAS DE CONTROL DIGITAL
Inicio de
la conversión
o
Fin de
la conversión II Lp
Salida
del reloj
o
Bit n
LSB o
Figura 2-16. Diagrama de temporización de un convertidor A/D por aproximaciones sucesivas.
CONVERSIÓN Y PROCESAMIENTO DE SEÑALES 31
es (r)
vs A
— V
esta señal primero pasa por el dispositivo S/H, el cual muestrea y mantiene la entrada
hasta que la conversión ha terminado.
Para ilustrar las implicaciones del tiempo de conversión, considérese la señal
analógica de la figura 2-17. El objetivo del uso de un convertidor A/D es convertir el
valor Vs en un número digital. Si el tiempo necesario para hacer la medición o
conversión es n, entonces existe una incertidumbre en la amplitud medida si la señal
cambia en el intervalo definido por n. Sea VFS el valor de escala completa de la señal
de la figura 2-17, y supóngase que el dispositivo A/D es un convertidor binario de n
bits. Lo anterior significa que la señal será digitalizada con una resolución de n bits;
esto es, de una parte en 2". El cambio de amplitud en la señal durante el tiempo de
conversión n puede aproximarse con la siguiente expresión:
des(t)
A V x (2-9)
dt t=t<
Para una resolución de 2-1%
AV
< 2 -~" (2-10)
VFS
Tabla 2-4 Limitac ones en frecuencia del A/D debidas al tiempo de conversión Tc
Longitud de
palabra n
bits
Tiempo de 4 6 8 10 12 14 16
conversión T (s)
Frecuencia máxima (rad/s)
10-3 MS) 125 31.25 7.8125 1.953 0.488 0.122 0.0305
10 -4 (100 ps) 1250 312.5 78.125 19.53 4.88 1.22 0.305
10-5 (10 AIS) 12 500 3125 781.25 195.3 48.8 12.2 3.05
10-6 (1 ps) 125 000 31 250 7812.5 1953 488 122 30.5
10-7 (100 ns) 1 250 000 312 500 78 125 19 530 4880 1220 305
CONVERSIÓN Y PROCESAMIENTO DE SEÑALES 33
ejemplo del convertidor A/D y el S/H de 10 bits dado anteriormente, el retraso total
es 1 ps. Si se emplea una frecuencia de muestreo de 1 MHz (que en teoría es la
frecuencia de muestreo máxima desde el punto de vista del S/H y el A/D), la señal de
entrada no puede tener componentes en frecuencia mayores de 0.5 MHz.
En los capítulos finales del libro se demostrará que la estabilidad de lazo cerrado
de un sistema de control digital se encuentra muy relacionada con el periodo de
muestreo. En muchos casos, una frecuencia de muestreo baja tiene efectos degrada-
torios sobre la estabilidad de los sistemas de lazo cerrado; por consiguiente, la
selección de una frecuencia de muestreo adecuada para un sistema de control digital
ha de basarse en consideraciones de estabilidad y de procesamiento de la señal.
Es interesante comparar los diagramas de bloque de los convertidores D/A (Fig. 2-13)
y A/D (Fig. 2-14). Si la resolución del convertidor A/D es muy grande, puede
despreciarse el efecto no lineal de la cuantización, y ya que las relaciones de
transferencia del codificador y del decodificador pueden representarse con ganancias
constantes, los dos diagramas se reducen esencialmente a una operación S/H. Por otra
parte, éste es el caso en los estudios analíticos de los sistemas de control digital; por
ejemplo, para fines de análisis, el sistema digital del piloto automático de la figura
1-8 puede representarse con el diagrama de bloques de la figura 2-18, donde los
convertidores A/D y D/A están simbolizados con dispositivos S/H. En este caso, el
controlador digital está representado por la función de transferencia D(z), mientras
que los elementos analógicos están modelados por sus correspondientes funciones de
transferencia analógicas. Quizá otra buena razón para despreciar los efectos de la
cuantización es que el cuantizador es un elemento no lineal. Si bien es posible tratar
de manera matemática la no linealidad de la cuantización y estimar el error producido
por ésta, en realidad no hay un método general que permita el diseño de una clase de
sistemas de control digital con cuantizadores.
Controlador Dinámica
digital del fuselaje
e e
Gf(s) S/I1 D(z) S/H GA (s) 1/s
Prefiltro
G. (s)
Giroscop'o
G (s)
Transductor
de posición
Figura 2-18. Diagrama de bloques analítico del piloto automático digital de la figura 1-18.
CONVERSIÓN Y PROCESAMIENTO DE SEÑALES 35
f(t) ft, (0
T (p)
Figura 2-19. Muestreador de frecuencia de muestreo uniforme con duración de muestreo finita.
36 SISTEMAS DE CONTROL DIGITAL
modulador de amplitud de pulso. Para ello, se considera que la entrada f(t) está
multiplicada por una señal portadora p(t) formada por un tren de pulsos periódicos,
cada uno con amplitud unitaria. La figura 2-21 ilustra formas de onda representativas
para la señal de entrada f(t), la portadora p(t) y la salidaff,i(t).
Algunas operaciones de muestreo pueden describirse con un muestreador con
modulación de ancho de pulso. Es característico que la salida de este tipo de
muestreador esté formada por pulsos cuya duración cambia en función de la amplitud
de la entrada en los instantes de muestreo, en tanto que la amplitud de los pulsos es
constante. La figura 2-22 ilustra las señales de entrada y salida de un modulador de
ancho de pulso. También existe un esquema aún más elaborado en que la amplitud y
t
p(r)11 I
El -11-P EI El 11
hs-T-411
filr-lr111(.fir
HP,-T-01
Figura 2-21. Formas de onda características para la entrada y la salida del muestreador uniforme.
f(0
fw(t)
11 II • t
he-T-01
Figura 2-22. Señales de entrada y salida características del modulador de ancho de pulso.
CONVERSIÓN Y PROCESAMIENTO DE SEÑALES 37
P(t)
1 1-1
•• ••
en frecuencia mucho mayores que f(t); así pues, el muestreador puede considerarse
como un generador de armónicos.
Dado que el tren de pulsos unitarios p(t) es una función periódica con periodo T,
puede representarse con la serie de Fourier,
00
Dado que n se extiende desde —°° hasta co, la última ecuación también puede
expresarse de la siguiente manera
CO
P;(jco) = C„F(jco + jnco.,) (2-29)
n = — co
Co = lím C„ (2-30)
n-'0O
4w 2x - r 2n 4x
P
(.0
-6w, -4w, -2w, -w, o w, 2w, 4w, 6w,
IF(jw)I
• In (jw)I
•••
---
4x 2w 27C 47C
.. Tv... ./T\ P
-6w,
7N-..-.-`191(
-4w, o
I
w,
ANtf>13-7;r
i
2w, 4w,
-111 6w,
- 0)s - (0c W c 10i
2 2
I In (i(0)1
la cual puede emplearse como medida cualitativa para ilustrar el espectro de amplitud
de FINO. Si el espectro de amplitud de la entrada de tiempo continuo f(t) está
limitado en banda, como se muestra en la figura 2-24(b), el espectro de amplitud de
Fp*(10) tendrá la forma que aparece en la figura 2-24(c). Nótese que la gráfica
de la figura 2-24(c) no es exactamente la de IF;(4))1, pero indica la forma general del
espectro.
El requisito de que ws, sea por lo menos dos veces mayor que el componente de
mayor frecuencia contenido en la señal f(t) (frecuencia de Nyquist) se conoce como
teorema del muestreo, y se plantea de manera formal en la sección 2-8.
Claro está que en la realidad muchas de las señales que aparecen en los sistemas
no exhiben una frecuencia de corte abrupta en el espectro de amplitud como la que se
muestra en la figura 2-24(b), de modo que siempre existe un efecto de traslapamiento
a pesar incluso de que la frecuencia de muestreo sea mayor que esta última. Adelante
se mencionan más aspectos e implicaciones físicas de la frecuencia de traslapamiento
sobre el desempeño de los sistemas de control digital.
,r
Lo anterior significa que cuando se conocen los valores de la primera derivada def(t)
en t = kTpara k = 0, 1, 2, ..., además de los valores de f(k7), la frecuencia de muestreo
mínima es coc, que es la mitad de la requierida cuando sólo se mide f(kT). Por otra
parte, la adición de cada derivada subsecuente de mayor orden permite alargar el
tiempo entre las muestras de acuerdo con (1/2)(n + 1)(2x/co,), donde n es el orden de
la derivada de mayor orden cuando se observan las demás derivadas de menor orden
en todas las muestras.
dD(Z)
IY = (2-40)
4= 4,,
donde 1„ es el n-ésimo polo simple de F(1), con n = 1, 2, ..., k.
CONVERSIÓN Y PROCESAMIENTO DE SEÑALES 45
1 al - 1
. Snr"F(S) (2-42)
K"` (i — 1)! Os' s=
F7,(s) = E f (kT)
°
G [1 — e "1 - kTs
e (2-45)
k=0
Dado que la hipótesis es que la duración del muestreo p es muy pequeña, el
término e-P's puede aproximarse tomando sólo los primeros dos términos del corres-
pondiente desarrollo en serie de potencias. En consecuencia,
(ps)2
1 — e -Ps = 1 — -[1 ps + ps (2-46)
2!
Por consiguiente, la ecuación 2-45 puede escribirse como
CO
Figura 2 - 26. (a) Muestreador con ancho de pulso finito. (b) Muestreador ideal conectado en
serie con un atenuador p. (a) y (b) equivalen si p es mucho más pequeño que T y que la menor
constante de tiempo de f(t).
Por tanto, una aproximación del muestreador de ancho de pulso finito es el muestrea-
dor ideal conectado en serie con un atenuador con atenuación p. Las figuras 2-26(a)
y (b) ilustran esta situación.
El muestreador ideal está definido como un muestreador que abre y cierra de
manera instantánea, en tiempo cero, cada T segundos.
La salida del muestreador ideal se expresa como
fió
8r (t)
t •••
t t t t t t •t
o T 2T 3T 4T 5T 6T 7T
f(27)
«AM it
f(37) f(7T)
► t
o ll +
f(67) 4
f(4T) f(57)
(c)
Figura 2-27. (a) Señal de entrada al muestreador ideal. (b) Tren de impulsos unitarios. (c) Salida
del muestreador ideal. f* (t) f(1)67(1).
muestreo correspondientes. Dado que la función impulso está definida como un pulso
de duración cero y amplitud infinita, en la figura 2-27 los impulsos aparecen indicados
con flechas cuyas amplitudes representan la intensidad de los impulsos.
Es posible derivar otras expresiones para F*(s) comparables con las de Fp*(s) en las
ecuaciones 2-37 a 2-42. A continuación se resumen estas expresiones para el proceso
de muestreo ideal. Los detalles de las deducciones se encuentran en el apéndice B.
f *(t) = i
cc f (kT)8(t — kT) (2-52)
kO
Dominio de la frecuencia
F*(s) =
i(0 +) + — 1 F(s + jnw) (2-53)
2 T„=- ao
48 SISTEMAS DE CONTROL DIGITAL
donde
=N DW (2-56)
dD(1)
= (2-57)
4 z„
y en es el n-ésimo polo de F (e), n = 1, 2, ..., k.
11
donde el residuo en el polo e = e., el cual tiene una multiplicidad mn (>1) es
1)! 1m" -
[
- 1.)-"Fw 1 — 014=1.
(2-60)
Ejemplo 2-1
F*(s) — (2-68)
1 — e - Ts s=s+a 1 — e—(s+a""
Ejemplo 2-2
Como ejemplo de la aplicación de las ecuaciones 2-58 y 2-60 a una F(s) con polos
de orden múltiple, considere el proceso de muestreo de una función rampa unitaria
con un muestreador ideal. Sea
f (t) = tu s(t) (2-69)
a
K12 = s2 F( s) =0 (2-72)
as s=0
p2 1 1 1 Te'
F*(s) = z o = (1 e- Ts)z (2-74)
V 1— e - T(s —
En principio, puede emplearse la ecuación 2-53 para obtener los mismos resul-
tados en los dos últimos ejemplos, pero los cálculos son más complicados. Por esto
la ecuación 2-53 se utiliza muy raras veces para determinar la solución en forma
cerrada de F * (s). Si f(t) está dada y la señal es muestreada con un muestreador ideal,
IF*(jco)1
••• •••
• 0)
ms ws
—w s — 2°2-1 —wc wc
2
2w,
(b)
Figura 2-28. Espectro de amplitud de las señales de entrada y salida de un muestreador ideal.
(a) Espectro de amplitud de la señal continua de entrada f(t). (b) Espectro de amplitud de la
salida del muestreador ((os > 2ffic).
CONVERSIÓN Y PROCESAMIENTO DE SEÑALES 51
entonces puede usarse la ecuación 2-54 a fin de hallar F*(s), aunque la desventaja es
que la respuesta se obtiene como una serie infinita, a partir de la cual debe encontrarse
la solución en forma cerrada. Como alternativa, cabe tomar la transformada de
Laplace de f(t) y utilizar las ecuaciones 2-55 o 2-58, según el tipo de polos de F(s).
Si F(s) es dada, se pueden utilizar directamente las ecuaciones 2-55 o 2-58; no hay
necesidad de optar por la transformada inversa de Laplace de F(s) para hallar primero
f(t) y después utilizar la ecuación 2-54.
La ecuación 2-53 es útil al interpretar la operación del muestreador ideal como un
generador de armónicos. Según indica esta ecuación, el muestreador ideal reproduce el
espectro de la entrada f(t) en la salida continua de éste, así como los componentes
armónicos que son múltiplos enteros de la frecuencia de muestreo ah y todos los que
tienen la amplitud 1/T. Ahora suponga que el espectro de amplitud de f(t) está limitado
en banda, como se ilustra en la figura 2-28(a), donde oh es la frecuencia de muestreo en
rad/s y wc es la frecuencia del componente de mayor frecuencia contenido en f(t). De
nuevo, si la frecuencia de muestreo es menor que 2coc, el espectro de la salida estará
distorsionado debido al traslapamiento de las bandas laterales en iF*(j col También es
interesante notar que, para el muestreo con ancho de pulso finito, la amplitud del espectro
de ¡Fp* (jco)i disminuye a medida que co aumenta, en tanto que en el caso del muestreo
ideal, la amplitud es constante para toda cf.).
Bandas
complementarias
Banda
primaria
Bandas
complementarias
P0)5/2
j5o)312
ficosa
jco,12
-.10s12
-fico,12
-Mos12
-1%s/2
Musa
X X
j5cos/2
X X
j3ws12
- --#
a
o
-"12
x- - -- -XE *
-
Polos de F(s) X X
-j5cosi2
X X Polos de F*(s)
ws12
(a) (b)
Figura 2-31. Polos de F(s) y F*(s) que ilustran los efectos del traslapamiento.
54 SISTEMAS DE CONTROL DIGITAL
d'Y(t)
f(")(kT) = (2-79)
dt" t=kT
para n = 1, 2,...
A fin de evaluar los coeficientes de la serie de potencias que aparecen en la
ecuación 2-77, es necesario obtener el valor dé las derivadas de f(t) en los instantes
de muestreo. Dado que sólo se tiene información sobre f(t) y su magnitud en los
instantes de muestreo, habrá que estimar las derivadas de f(t) a partir de los valores
CONVERSIÓN Y PROCESAMIENTO DE SEÑALES 55
def(kT). Una expresión sencilla para hacerlo donde sólo aparecen dos pulsos de datos
y que da una estimación de la primera derivada de f(t) en t = kT es
f (1) (kT) = 7. [f ( kT) — f [(k — 1)7]] (2-80)
donde
f ( 1) [(k — 1 ) T] = — [f[( k — 1)7] — f [(k — 2)7]] (2-82)
Por las razones mencionadas, a menudo sólo se emplea el primer término de la ecua-
ción 2-77 a fin de aproximar f(t) en el intervalo comprendido entre kT S t < (k + 1)T.
Esto es,
fk(t) = f (kT) (2-84)
1 Impulso unitario 1
o T
(a) (b)
Figura 2-32. (a) Entrada impulso unitario al retenedor de orden cero. (b) Respuesta al impulso
del retenedor de orden cero.
Este dispositivo se conoce más como retenedor de orden cero (ROC), pues mantiene
el valor de f(kT) en el intervalo kT S t < (k+ 1)T hasta que llegue la siguiente muestra,
jr(k + 1)7]. Este retenedor tiene las mismas características que el dispositivo S/H
estudiado en la sección 2-4.
Puede demostrarse con facilidad que el ROC es un dispositivo lineal, ya que
satisface el principio de superposición. Para obtener la función de transferencia del
ROC, considérese la respuesta al impulso de la figura 2-32(b), la cual es resultado
evidente de la ecuación 2-84. La respuesta al impulso puede expresarse como
1 — e - Ts
G h o(S) = (2-86)
s
La figura 2-33(b) muestra la salida h(t) del ROC cuando la entrada a éste es una
secuencia de pulsos f(k7), k = 0, 1, 2,..., o una señal muestreadaf*(0. Las señales de
la figura 2-33 indican con claridad que la exactitud del retenedor de orden cero como
dispositivo de extrapolación &pende, en gran medida, de la magnitud de la frecuencia
de muestreo ois. Conforme esta frecuencia aumenta hasta el infinito, T tiende a cero
y la salida h(t) del ROC tiende hacia la señal continua en el tiempof(t). En el siguiente
capítulo se verifica este hecho de manera analítica.
"Oi
n
\ 7T 8T
•
o T 2T 3T 4T 5T 6? "- ,L V 9T 10T
h(0 4
7T 8T t. 5
O T 2T 3T 4T 5T 6T n 10T
9T
Figura 2 - 33. Funcionamiento del retenedor de orden cero en el dominio del tiempo. (a) Señal de
entradaf(t) y señal muestreadap(t). (b) Forma de onda de salida del retenedor de orden cero.
o
sen (coT /2) - i.772
Gh000)) = T e (2-89)
(2)772
•
Dado que T es el periodo de muestreo en segundos y T = 2n/cos, donde 0, es la
frecuencia de muestreo en rad/s, la ecuación 2-89 se convierte en
(nicowslcos) _
0 sennto
27
Gho( jw) — e (2-90)
La magnitud de Gho(jw) es
2n sen (nco/tos)
G( ico s)i = (2-91)
cos itcplcos
El signo de sen(nTo/cos) cambia para cada múltiplo entero de Ira 3axf. Este cambio de
signo de + a — puede considerarse como un cambio de fase de —180'.
58 SISTEMAS DE CONTROL DIGITAL
10.00
0.00
20.00
40.00
-60.00
-80.00
100.00
120.00
0.10E+00 0.10E+01 0.10E+02
0.00
360.00
720.00
-1080.00
o
-1440.00
-1800.00
3
ZD -2160.00
-2520.00
-2880.00
-3240.00
-3600.00
0.10E+00 0.10E+01 0.10E+02
cok%
Figura 2-34. Gráficas de Bode de las características de ganancia y fase del retenedor de
orden cero.
z.uu
1 .80 ---- -- --
--- -------- - - - . - -- --- ---------
1 .40 ammiarniumiamiaman
1.0:8000 1.1.1111111111111111111111111111111111111111
10111111.1111101111111111111111111111111111111. I.
0 11111111.1111111111111111111111111111111111111
0:6° 111111111111111111111111111111111111111 .
° 4 1.1011111111111111111111111111111111111
w er,wilig_ItIllffl_11111111111111111111 -----
0110 1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 9.00 10.
*A%
0.00
-360.00
-720.00
-1080.00
-1440.00
g -1800.00
j
,1
- -2160.00
-2520.00
-2880.00
-3240.00
-3600.00
0.00
(t — 2T)
+ us(t — 2T)us( — 2T) + us(t — 2T) (2-100)
G h 1 (S) =
1 + T s [1 -- e-12 (2-101)
T s
o
1+Ts
G hl( s) — T [ Gho(s)] 2 (2-102)
CONVERSIÓN Y PROCESAMIENTO DE SEÑALES 61
Impulso unitario
(a)
Figura 2-36. (a) Impulso unitario. (b) Respuesta al impulso del retenedor de primer orden.
111n41::lmit nuaf(1)
2T 3T 4T 5T 6T
111115
Figura 2-37. Reconstrucción de una señal continua con un retenedor de primer orden.
62 SISTEMAS DE CONTROL DIGITAL
0.00
—360.00
—720.00
75, —1080.00
—1440.00
—1800.00
.1-2160.00
Z —2520.00
—2880.00
—3240.00
—3600.00
0.10E+00 0.10E+01 0.10E+02
(0/(os
Figura 2-38. Gráficas de Bode de las características de ganancia y fase del retenedor de
primer orden.
•
2.00
1.80
1.60 xmAlu
WAVAVA .W.V.W.TYMWV.W.W.W/W*
1.40 Lw:*
1.20 IrtA
1.00
0.80
0.60 ano
111111111
0.40
0.20
0.00
0.00 1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 9.00. 10.00
oilw
0.00
-360.00
-720.00
-1080.00
1 -1440.00
-1800.00
-2160.00
4.5
-2520.00
-2880.00
3240.00
3600.00
0.00 1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 9.00 10.00
0)/(os
caso del ROC, pero con una pendiente de descenso más pronunciada. Sin embargo,
los valores pico del RPO en cada intervalo definido por micos poseen Una magnitud
mayor que los del ROC, lo cual explica la presencia de componentes tipo rampa en
la salida del RPO. La característica en fase del RPO, tal como está dada por la.ecuación
2-105, es una curva suave para todos los valores de u).
La figura 2-40 presenta una comparación entre las respuestas del ROC y el RPO
en el dominio de la frecuencia, desde co/cos = O hasta Oros = 5. Nótese que para
frecuencias pequeñas sobre cada intervalo co/cos, el retraso en fase prOducido por el
ROC es mayor que el del RPO, pero para frecuenciat grandes sucede lo contrario. En
los sistemas de control retroalimentados, las características de la fimción de transfe-
rencia de lazo abierto a frecuencias grandes son las que, en general, determinan la
estabilidad del sistema de lazo cerrado; por tanto, en general, la aplicación de un RPO
crea más problemas de estabilidad que el ROC.
64 SISTEMAS DE CONTROL DIGITAL
2.00
1.80
1.60
Wilil 11111.11111111111111111111111111111111111111
1.40
1.20 1111111111111111111111111
1111
1.00
0.80
ROC
111111111111111111111
0.60
ummumemema
MIMO
0.40
0.20 11111&111~11111.1111111
wwwww".8115~.~.----
0.00
0.00 0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.00
ca3ws
0.00
-180.00
-360.00
-540.00
-720.00
-900.00
-1080.00
-1260.00
-1440.00
-1620.00
-1800.00
0.00 0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.00
tolw,
Figura 2-40. Comparación de las características de los retenedores ROC, RPO y con retraso.
AL f(r)
h t)
3T
3 1( 4T 5T
T 2T
Nuw." 6T 7T ST 9T
f(t)
T 2T
3T 4T 5T../.411r,
6T 7T 8T 9T
Figura 2-42. Interpretación de la salida del retenedor poligonal como la suma de la salida de
dos retenedores triangulares h i(t) y h 2(0 .
SISTEMAS DE CONTROL DIGITAL
Te — 2 + e—Ts
G h p(S) = Tse (2-107)
(1 — e - 'V Ts
Ghp(s) = e (2-108)
Tse
eTs
[Gho(s)]2
Ghp(s) (2-109)
f(t) h(t)
4T 5T 6T
2T „ 7T 8T 9T 10T
Las características de amplitud del RP indican que es un filtro pasa bajas mejor que
el ROC y el RPO. La amplitud de Ghp(jco) disminuye con mayor rapidez después de
= 1.
2nco
L Ghs(jo) - rad (2-113)
Co s
Por tanto, la característica de amplitud del RR es igual a la del RP, pero la fase es una
línea recta, tal como se indica en la figura 2-40.
PROBLEMAS
2-1 Suponga que adquiere un automóvil nuevo con un préstamo p(0) a pagar en un
periodo de Nmeses, con una tasa de interés mensual igual a R por ciento. Los intereses
y el principal deben pagarse en N pagos iguales. Determine el monto U de cada pago
mensual. La ecuación de diferencias para este problema está dada por la expresión
68 SISTEMAS DE CONTROL DIGITAL
1-1. Resuelva el problema cuando p(0) = US$15 000, N = 36 meses y R = 1%; esto
es,R = 0.01.
2-2 Considere la compra de un auto nuevo con un préstamo a plazos igual ap(0), el
cual debe pagarse en N meses. La tasa de interés es R% mensual. El principal y los
intereses han de pagarse en Npagos mensuales del mismo monto U. Suponga que el
comprador sólo puede efectuar pagos mensuales de US$300 durante 48 meses.
Determine el monto máximo del préstamo que puede otorgársele si R = 1.2%; esto
es, R = 0.012.
2-3 Considere un préstamo de montop(0) a pagar en Nmeses, con una tasa de interés
mensual igtial a R% sobre el saldo insoluto. El principal y los intereses se pagan al
final de cada periodo, donde cada pago tiene un monto dos veces mayor que el del
periodo anterior. Escriba una ecuación de diferencias que exprese el monto de la deuda
al término del (k + 1)-ésimo periodo en función del lapso anterior y del pago inicial
u(0). Exprese la solución p(N) en términos de p(0) y u(0). Si N = 5, R = 0.01, p(0)
= US$50 000, encuentre u(0),"p(k), k = 1, 2, 3, 4, y el monto total (principal más
intereses) pagado en todo el periodo.
2-4 Se quiere obtener una expresión recursiva (ecuación de diferencias) para hallar
la r-ésima raíz de un número real N. Para ello, primero se escribe el desarrollo en serie
de Taylor de la función f(x) alrededor del punto xn:
(x — x„)2
f(x) = f (x„) + (x — x„)f(1)(x„) + f (2)(x„) + • • • (P2-4-1)
2!
donde
(x)
fl)(xn) =
dtk
Denote con x el valor de la siguiente iteración, x,,+1. Por otra parte, x también representa
la solución de la ecuaciónf(x) = 0. En consecuencia, la ecuación P2-4-2 se convierte en
O = f(xn) + (xn+ 1— xn)f(1)(x.) (P2-4-3)
o
f(x0)
Xn+ Xn fkxn) (P2-4-4)
2-7
o. Halle el tiempo de conversión máximo necesario para digitalizar una señal
senoidal de 1 kHz, v(t) = V sen no, con una resolución de 10 bits.
b. Demuestre que, para una resolución fija, la relacióri entre el tiempo de conversión
y la frecuencia de la señal senoidal de la parte a es, en coordenadas logarítmicas,
una línea recta. Ilustre los casos correspondientes a resoluciones de 10 y 16 bits
para una escala de tiempo de conversión de 10-9 a 10-3 sobre el eje vertical y
una frecuencia de 1 Hz a 100 kHz en el eje horizontal. Utilice coordenadas
logarítmicas.
2-8 Encuentre el voltaje que aparece en la salida de un convertidor D/A de ocho bits
con una escala completa de 10 V si la entrada es la lista de palabras digitales de la
siguiente tabla. Proporcione también el factor de ponderación de cada bit en términos
de una fracción de la escala completa.
00000001
00000010
00000100
00001000
00100000
10000000
11111111
Figura P2-9
2-13 Se considera que un muestreador con "retraso" está cerrado durante un lapso
muy breve p en los instantes t A, A -F T, A + 2T,..., A + nT, donde T> A >0yT es
el periodo de muestreo en segundos. La entrada al muestreador es una señal f(t)
continua. Se supone que la salida del muestreador, fi* (O& es un tren de pulsos
máximamente plano, dado que A 4 7'. Deduzca la transformada maestreada de
f/,` (t)e, Fp* (s, A). [Note que también puede hacerse referencia a esta transformada.como
"transformada muestreada con retraso" de "0.] Para p pequeño, suponga que 1 —
ps. Determine la tranformada muestreada con retraso de f(t) = e-2t.
2-14 La transformada muestreada con retraso del problema 2-13 puede aplicarse a
la consideración de ancho de pulso finito en sistemas de datos tnuestreados. Según se
indica en la figura P2-14(a)", un pulso de duración p puede considerarse como el
resultado de Npulsos elementales de duración A, tal que NA = p con A muy pequeño;
por tanto, un muestreador real puede representarse con N muestreadores So, SA, S,.,
SW-1», conectados en paralelo, como se aprecia en la figura P2-14(b). Lo anterior
implica que el muestreador So toma muestras en los instantes 0, T, 2T,.., nT,...; el
muestreador Se hace lo propio en los instantes A, A + T, A+ 2T,..., A+ nT, etc. Obtenga
la transformada muestreada de fp* (t).
CONVERSIÓN Y PROCESAMIENTO DE SEÑALES 71
o
p= T+p
fo(t)
_25 fa(t)
Se
f -1 (1)
S (N_ i» Figura P2-14 a) y b)
f(r)
-T Figura P2-15
f(t) Al)
T 2T 3T 4T ST 6T t 0 T 2T 3T 4T 5T 6T
(a) (b) Figura P2-16
Posición
Determine cuántos cruces por cero en las salidas del codificador corresponden a
un carácter del disco de impresión.
Suponga que se emplea un reloj de 500 kHz para medir la velocidad del disco de
impresión, en conjunción con las señales del codificador. La velocidad se mide
contando el número de pulsos entre dos cruces por cero consecutivos de la salida
del codificador. ¿Cuál es la velocidad mínima, en revoluciones por minuto (RPM),
que puede medirse con esta configuración sin sobrepasar la capacidad de un
contador de 12 bits?
c. El ángulo de giro máximo del disco de impresión es 180 grados en cualquier
dirección, o 48 caracteres. Este giro debe hacerse en un tiempo cuando mucho
igual a 30 ms. Si se supone que el disco representa una carga inercial pura y que
el movimiento de 48 caracteres se logra con un perfil de aceleración constante
durante el recorrido de la primera mitad de la distancia, con una desaceleración
constante en lo que resta de ésta, encuentre la velocidad máxima en RPM que el
codificador debe medir, en número de cruces por cero por segundo.
CONVERSIÓN Y PROCESAMIENTO DE SEÑALES 73
w(t) rad/s Gs
Pulsos por periodo
de muestreo T Figura P2 - 19
Bibllografia
Procesamiento digital dé señales
Rabiner, L.R. y C.M. Rader, Editores, Digital Signal Processing, IEEE Press,
Nueva York, 1972.
Helms, H.D. y L.R. Rabiner, Editores, Literature in Digital Processing.. Termi-
nology and Pennuted Iltle Index, IEEE Press, Nueva'York, 1972.
Openheim, A.V. y R.W. Schafer, Digital Signal Processing, Prentice-Hall, En-
glewood Cliffs, N.J., 1975.
75
76 SISTEMAS DE CONTROL DIGITAL
Y
Im z = re sen coT (3-5)
con
S = o ± jco (3-6)
La relación entre s y z, tal y como está dada por la ecuación 3-2, puede definirse
como la transformación z. Al sustituir la ecuación 3-2 en la 3-1, se tiene que
1
F*[s = — ln zj= F(z) = f (kT)z' (3-7)
k= o
F(z) = 21f(t)] = f(kT)z-k (3-14)
k=0
respectivamente. Dado que F(z) se obtiene al muestrear primero f(t) con un mues-
treador ideal con un periodo de muestreo T, existe la idea equivocada de que cuando
78 SISTEMAS DE CONTROL DIGITAL
Como f*(t) representa un tren de impulsos espaciados T segundos entre sí, a medida
que T se reduce, el tren de impulsos se convierte en un "racimo" de impulsos en t =
O, lo que de ninguna manera se parece a fit). El equivalente en el dominio del tiempo
de la ecuación 3-15 es
lím f *(t) f (t) (3-16)
T-+ 0
Entonces, ¿cómo explicar que, en la realidad, la serial muestreada tienda hacia la señal
continua cuando el periodo de muestreo tiende a cero? Este problema se trata en el
capítulo 4.
F(Z) = (3-18)
D(Z)
tiene k polos distintos en e_ fin, n = 1, 2, ..., k. D'(e,,) está defmida como
dD()
1) 11„) = (3-19)
clZ 474,,
F(z)
-1 r" - 'Km [a m. -- 1
(3-20)
„= 1 i 1 (Mn
as,„„-i _ e — Ts s=s_s
9. z7
TRANSFORMADA Z 79
donde
[Di-
K ,,. snr'nF(s)] (3-21)
(i — 1)! Lasi — i (s 5=5„
Ejemplo 3-1
zis(t)
1 I I I I 1 I t
2T 3T 4T sr 6T 7T 8T
Figura 3-1. (a) Función escalón unitario. (b) Función escalón unitario muestreada.
El desarrollo anterior se detalló con la finalidad de ilustrar los tres pasos para
evaluar F (z) a partir de f(t). En la práctica se sustituye f(kT) en la ecuación 3-11, sin
tratar con expresiones donde intervienen transformadas de Laplace.
TRANSFORMADA Z 81
Ejemplo 3-2
Esta serie infinita converge para todos los valores de z que satisfacen arz-l i < ic
1
= 1 e-aTz- 1 z e-aT
(3-33)
82 SISTEMAS DE CONTROL DIGITAL
Ejemplo 3-3
En este caso resulta más conveniente expresar sen con' en forma exponencial;
entonces, la ecuación 3-34 toma la forma
00 ejeukT e - jcokT
F(z) --= E- (3-35)
k = 0 2j
Esta serie infinita es convergente en lz-1 1< 1, y puede escribirse en forma cerrada
utilizando el mismo método descrito en los ejemplos 3-1 y 3-2. El resultado es
1 1
F(z) = 1 — -jcok Tz — 1.1 (3-36)
2j[l — eiwkr 1
N(1 ) = N( 2)_
ya que N( 4 ) no es función de
D'(1) = 2jw
Dn 2) = —2jw
TRANSFORMADA Z 83
2j (3-39)
j 1 — eiwT Z - 1 1 — e - iw 1. z - I
que aparentemente conduce al mismo resultado que el dado por la ecuación 3-37.
Ejemplo 3-4
la cual tiene un polo repetido dos veces en s = 0. Para obtener F(z) en forma cerrada
directamente a partir de F(s) es necesario emplear la ecuación 3-20 ya que, en este
caso, F(s) tiene un polo de segundo orden. En consecuencia, de las ecuaciones 3-20
y 3-21, resulta k = 1, s1 = O y mi = 2. La ecuación 3-20 da
1 rai- 1 2 1]
(i — 1) !
S — (3-46)
s2 s= o
SISTEMAS DE CONTROL DIGITAL
[ 1 Tz
F(z) = s2 —, eTsZ — 1 = 2
(3-48)
as s1- S 1)
Ejemplo 3-5
1 1 1
= 1 ± - Z -1 + - Z -2 ± - Z -3 + • • (3-50)
2 4 8
Note que, en este caso, no aparece el periodo de muestreo T, por lo que éste puede
considerarse igual a uno.
Con objeto de expresar F(z) en forma cerrada, se multiplican los dos miembros
de la ecuación 3-50 por (1/2)f' y el resultado se resta de la misma ecuación.
Después de resolver para F(z), se llega a la siguiente expresión
1 z
F(z) = (3-51)
( 1— z (z — 1)
iw
Plano s
cok
2
Banda primaria
o
•
2
(a)
Re z
(b)
Figura 3-2. Mapeo de la banda primaria del semiplano izquierdo del plano s en el plano z
realizado por la transformada z.
86 SISTEMAS DE CONTROL DIGITAL
círculo unitario (círculo con radio uno), con centro en el origen del plano z, mediante
la transformación z = en, según se ilustra en la figura 3-2(b). Dado que
e(s+inuo.,)T e TsejUme Ts = z (3-52)
con m entero, todas las bandas complementarias del semiplano izquierdo del plano
serán mapeadas también en el mismo circulo unitario del plano z; por tanto, a partir
del mapeo realizado por la transformada z entre los planos s y z cabe concluir.
Todos los puntos que están en el semiplano izquierdo del planos corresponden
a puntos dentro del círculo unitario del plano z.
Todos los puntos en el semiplano derecho del plano s corresponden a puntos fuera
del círculo unitario del plano z.
Todos los puntos sobre el eje j co del plano s, corresponden a puntos sobre el
círculo unitario izi = 1 del plano z.
A continuación se investiga el mapeo del plano s al plano z de los lugares de
amortiguamiento constante, defrecuencia constante y de relación de amortiguamien-
to constante.
•/`o
Pianos
Lugares de
amortiguamiento
constante
(a) (b)
Figura 3-3. Lugares de amortiguamiento constante en los planos s y z.
la cual representa una línea recta que parte del origen con un ángulo igual a O = anT
rad, medido a partir del eje positivo real, como se indica en la figura 3-4. En la misma
figura aparecen otros lugares de frecuencia constante en los planos s y z; por tanto,
los polos ubicados en los cuadrantes primero y cuarto del plano z tienen frecuencias
menores que las de los que están en los cuadrantes segundo y tercero.
fto
Plano s
ws
./(02
a Re z
(a) (b)
donde
fi = sen -1 C (3-56)
con fi medida en grados a partir del eje jco, como se ve en la figura 3-5(a). El mapeo
en el plano z de esta relación es
z e Ts e T(-w tan fi+ jw) e - 24m tan filo), brco/cos (3-57)
(a) (b)
Figura 3-5. Lugares de la relación de amortiguamiento constante para 13 = 30* (1= 50%) en los
planos s y z.
TRANSFORMADA Z 89
(a) (b)
Figura 3-6. Lugares de la relación de amortiguamiento constante acotados por la banda primaria
en el plano s y lugares correspondientes en el plano z.
Dado que muchos sistemas de control tienen características pasa bajas, lo usual
es que sólo tengan importancia los polos próximos a la banda primaria. Lo anterior
significa que, para fines prácticos, nada más tiene interés la parte del lugar de relación
de amortiguamiento constante que corresponde a frecuencias que van desde co = O
hasta co = cos/2; por ejemplo, si se desea una relación de amortiguamiento de 70.7%
o mejor (fi 45'), basta considerar las regiones que abarca la banda primaria en los
planos s y z (Fig. 3-6).
O T 2T 3T 4T 5T 6T 7T 8T 9T
muestreada por un rnuestreador ideal, se pierde por completo la información entre los
instantes de muestreo; por esta razón no es posible esperar la recuperación de f(t) a
partir de F(z). La figura 3-7 ilustra el caso sencillo en que se toma la transformada z
de la función escalón unitario; en este caso, la operación de muestreo genera un tren de
impulsos unitarios con un periodo de muestreo igual a T s. La transformada z inversa
de F(z) no es única, ya que f(t) puede ser cualquier función de t con un valor unitario
en los instantes de muestreo.
El problema de la no unicidad de la transformada z inversa es una de las
limitaciones del método de dicha transformada. Sin embargo, su simplicidad y
ventajas superan con mucho tal desventaja. Aún más, si el periodo de muestreo se
elige con cuidado, tomando en cuenta la variación en el tiempo de la función f(t), la
secuencia f(k7), e incluso la transformada z F(z), constituyen una representación
muy exacta de la función continua f(t).
f (kT) = z )]
[F(z)] (3-58)
Tabla de transformadas
Para funciones sencillas, la transformada z inversa puede buscarse en una tabla de
transformadas z, como la que aparece en el apéndice C.
Función escalón
z
F(z) = (3-61)
z—1
Función rampa
Tz2
F(z) — (z _ 1)2 (3-62)
Función parábola
2z(z + 1)
F(z) — T (3-63)
(z — 1)3
Un solución muy sencilla para esta situación es "guardar" la z del numerador de
F(z)(si es que la hay) y desarrollar en fracciones parciales F(z)/z. Después de efectuar
el desarrollo, la z en el denominador de F(z) es multiplicada del otro lado de la
expresión por cada una de las fracciones parciales. Si F(z) no tiene un cero en z = O,
como sucede cuando el sistema presenta un retraso puro que es múltiplo entero del
periodo de muestreo T, el/desarrollo en fracciones parciales de F(z) se realiza como
en la ecuación 3-59; por ejemplo,
2,- ir 1
L z - (kT) (344)
donde k= 1, 2, 3, ... Otro enfoque es ver que la ecuación 3-64 denota la secuencia 1,
1, 1, ..., la cual comienza en k = 1.
Ejemplo 3-6
Considere la transformada z
(1 — ar)z
F(z) = (3-65)
(z - 1)(z - ar)
Un programa para computadora digital como PFEZ del paquete DCSP puede
tener mucha utilidad para el desarrollo de transformadas z en fracciones parciales.
TRANSFORMADA Z 93
Ejemplo 3-7
Para la función F(z) dada por la ecuación 3-65, al dividir los polinomios del numerador
y denominador se tiene
F(z)= ( 1 _ e -ar)z - + (1 _ e -2a1)z -2 (1 — e 3° T)z7 3 + • •• (3-71)
donde, en este caso, el valor inicial de k puede ser cero o uno y el resultado de la
ecuación 3-72 es el mismo. Es evidente que la ecuación 3-71 indica quef(0) = O.
.ft0
C tic»
c —ioo
(a) (b)
Se observa que la trayectoria de integración pasa verticalmente por todas las bandas
periódicas del plano s; por tanto, la integral de la ecuación 3-75 puede separarse en
una suma de integrales, cada una evaluada dentro de una banda periódica. En
consecuencia, la ecuación 3-85 puede escribirse como
1 a) c + joh( + 1/2)
f(kT) = — E F(s)ekTs ds (3-76)
bri i = — oo c + jeo,(i — 1 /2)
donde (os = 2ir/T. Al reemplazar s por s +jicos, donde i es un entero, la ecuación 3-76
se convierte en
1 co c + jahl 2
f(kT)= — F(s + jicos)ekT(s+i ") d(s + jicos) (3- 77)
27rii = — co c — jw,I2
Como exp[kT(s +fico.,)]= exp (kTs) para valores enteros de j e i y d(s + ficos) = ds, la
ecuación 3-77 puede escribirse como
1 co c + jw,I2
f(kT) = — E F(s + jiws)ekTs ds (3- 78)
bri i = — oo c — jw,I2
=
1
T
ln zj= T z dz 3-84)
Ejemplo 3-8
1 z(1 - e - aT)
f(kT) = 2irj dz (3-86)
(z - 1Xz - e- aT) Zk— 1
donde 1' es un círculo con centro en z = O y radio suficiente para abarcar los polos de
F(z), los cuales están en z = 1 y z = cal.. La integral de contorno de la ecuación 3-86
suele evaluarse con el teorema de los residuos de la teoría de variable compleja. Si el
integrando de una integral contiene una función racional —esto es, una función con
un número finito de polos y ceros— el método de los residuos es igual, en esencia, al
para k = 0, 1, 2, Este resultado concuerda con los que aparecen en las ecuaciones
3-68 y 3-72, los cuales se obtuvieron empleando los dos métodos descritos.
. o 1. Adición y sustracción
Si las transformadas z deji(0 yf2(t) son Fi(z) y F2(z) respectivamente, y tenemos
= E fi(k nz-k ± E Mknz-k (3-92)
k = 4 k=O
F1(z) F2(z)
TRANSFORMADA / 97
f[f(t — nT)us(t — nT)] = f(kT — nT)z - (3-96)
k=0
Es necesario señalar que el teorema anterior sólo puede aplicarse a una función
real con un retraso que sea múltiplo entero del periodo de muestreo. En la práctica,
es frecuente encontrar situaciones donde el retraso no es un múltiplo entero del periodo
de muestreo y que, por tanto, deben tratarse de manera diferente.
Ejemplo 3-9
Ejemplo 3-10'
z sencoT
F(z) = f[sen cot] = (3-109)
— 2z cos coT + 1
100 SISTEMAS DE CONTROL DIGITAL
si el límite existe.
El teorema simplemente afirma que el valor de f(0) es igual al valor de F(z) en
z = co. n
e
Teorema del valor final
Si la transformada z de la funciónf(t) es F(z), y si la función
(1 — z-1)F(z) no tiene polos sobre el círculo unitario IzI = 1 en
el plano z o fuera de éste, entonces
lím f(kT) = lím (1 — z')F(z) (3-114)
k—.co
E f[(k — 1)T]z -k = f(0)z' + f(T)z - 2 + • • • + f[(n — 1)1]z - " (3-116)
k=0
TRANSFORMADA 2' 101
donde n es un entero positivo. Como f(t) O para toda t < O, el término f(-7) que
debería aparecer en la ecuación 3-116 es cero. Ahora, al comparar las ecuaciones 3-115
y 3-116, se advierte que ésta puede escribirse como
n-1
f[(k 1)T]z -k = z - 1 E f(kT)z -k (3-117)
k=0 k=0
Al efectuar la diferencia entre las ecuaciones 3-115 y 3-117, y dejar que z tienda a
uno, se tiene
n-1 n—1
lím [E f(kT)z -k — z - E f(kT)z E f(kT) E f(kT) = f(nT)
z •-• 1 k=0 k = 0 k=0
(3-118)
. (3-119)
n4. 00 n —n 1
z—P =O k=0
Al intercambiar los límites del miembro derecho de la ecuación 3-119, y puesto que
n—1
lím E f (kT)z k = him E f (k'T)z -k = F(z) (3-120)
n—, co k=0 n—.co k=0
Ejemplo 3-11
el valor final def(k7) se encuentra a partir de F(z) mediante el teorema del valor final,
ecuación 3-114. Como
0.792z
(1 — z')F(z) = (3-123)
z2 — 0.416z + 0.208
no tiene polos o ceros sobre el círculo unitario JzJ =1 o fuera de éste, cabe aplicar el
teorema del valor final; por tanto, de la ecuación 3-114,
0.792z
lím f(kT) = lím — 1 (3-124)
k—>co
y__,1 Z 2 — 0.416z + 0.208
102 SISTEMAS DE CONTROL DIGITAL
Este resultado puede comprobarse con rapidez si se desarrolla F(z) en una serie de
potencias en z-1.
F(z) = 0.792z -1 + 1.12z -2 + 1.091z -3 + 1.010z -4 + 0.983z-5
+ 0.989z -6 + 0.990z + • • (3-125)
Ejemplo 3-12
a
.2[te -at] = e-ati (3-128)
da
La transformada z de ca' es z/(z —e—a7); por consiguiente, a partir de la ecuación
3-126 la transformada z de te' es
= a 1
797. _
O z Tze'T (3-129)
a a Lz — aT (z e — aT)2
TRANSFORMADA Z 103
donde
Demostración: por definición, la transformada z del producto defi (t) por f2(t) es
Puesto que
-k
F2(Zr 1 ) -= f2(kTXZC (3-144)
k=O
TRANSFORMADA Z 105
expresión que converge de manera absoluta para lzr i l> Q2, o l l < 1z1/ a2; con esto, la
ecuación 3-141 se convierte en
1 Fig)F2(zr
[A(t)f2(t)] ="j7- 1) (3-145)
ri C
con
Ejemplo 3-13
donde I" es una trayectoria circular que se encuentra en el interior del anillo descrito
por
Y IzI > 1.
, Por consiguiente, la trayectoria de integración de la ecuación 3-153 contiene sólo
los polos del integrando en Ç = 1. La forma usual de evaluar la integral de contorno
de la ecuación 3-153 es aplicar el teorema de los residuos de variables complejas. En
este caso, dado que el integrando de la integral de la ecuación 3-153 es una función
racional, puede emplearse el método de desarrollo en fracciones parciales. El coefi-
ciente del desarrollo que corresponde al término (1-- 1)2 es
a[ nc-i l Tze-°T
(3-155)
1f1(0f2(0] = (zu 1 e— ar
) 1 — e' T)2
que concuerda con el resultado obtenido en la ecuación 3-129 del ejemplo 3-12.
donde
p(k) = monto adeudado al final del k-ésimo periodo
R = tasa de interés mensual (por ciento)
U = pago mensual (constante)
La cantidad abonada al inicio esp(0). Por otra parte, el plazo de pago es Nmeses.
La ecuación 3-158 es una ecuación de diferencias lineal de primer orden con
coeficientes constantes, la cual puede resolverse con el método de la transformada z.
Sin embargo, antes de usar este método, es conveniente examinar la sólución en el
dominio real utilizando para ello el método recursivo.
Al sustituir k= O en la ecuación 3-158 se tiene
p(1) = (1 + R)p(0) — U (3-159)
108 SISTEMAS DE CONTROL DIGITAL
donde
X = (1 + R)N-1 U — (1 + R)N - 2 U — • ••— (1 + R)U — U (3-162)
La suma de la última serie puede obtenerse multiplicando ambos miembros de la
ecuación por (1 + R), para después restar la ecuación 3-162 de la expresión resultante,
con lo que se obtiene
RX = (1 + R) N U — U (3-163)
Al resolver para X la última ecuación, se tiene que
(1 + R)N
1
X= (3-164)
R R
La sustitución de X en la ecuación 3-161 permite obtener la solución para p(N).
(3-183)
Por tanto,
y(k) = —0.909(0.5) k + 2.754(0.316) k cos (2.255k — 1.234) (k 0) (3-184)
Muestreador ficticio
--- 75 - — — c•(t)
S1 T
r(t) r*(t)
• nt G(s) c(t)
Si T
I Muestreador ficticio
Muestreador
I — — — — - — — ► c*(t)N
SN TIN
Figura 3-11. Respuesta de la salida del muestreador múltiple con un periodo de muestreo de T/4.
Si con un periodo de muestreo T, de modo que c(t) esté caracterizada por la señal
muestreada c*(t). Sin embargo, si c(t) contiene componentes de alta frecuencia
demasiado rápidas para el periodo de muestreo T, se puede introducir otro muestreador
ficticio SN en la salida con un-Periodo de muestreo T/N, donde Nes un entero positivo
mayor que uno. La salida de este muestrador ficticio SN puede escribirse como
c*(t)N — c(kT N)5(t — kT IN) (3-185)
k=O
•
la cual contiene N — 1 pulsos adicionales entre los instantes de muestreo kT y (k+ 1)T
del muestreador Si. La figura 3-11 ilustra una serial muestreada con este muestreador,
con un periodo de muestreo igual a T/4.
Al tomar la transformada z de los dos miembros de la ecuación 3-185 se tiene
c*(t) r c*(t—áT)
Otro modo de representar los detalles de una señal entre los instantes de muestreo, si
el periodo de muestreo es T, es retrasar la señal analógica un tiempo AT, donde
O < A < 1; hecho esto, la señal retrasada se muestrea con un muestreador convencional
en t = kT, k= 0, 1, 2, ... Al modificar la magnitud del retraso AT, es posible recuperar
toda la información que se quiera sobre la señal entre los instantes de muestreo. La
figura 3-12 presenta un sistema con un retraso ficticio AT en la salida, donde el
muestreo de la señal retrasada se efectúa con un muestreador con periodo T.
La salida retrasada del sistema de la figura 3-12 es c(t — A7)14(t — Al). La salida
muestreada del retraso ficticio puede expresarse como
Figura 3 - 13. Formas de onda que ilustran las operaciones que hay detrás de la transformada z
con retraso.
1 14 SISTEMAS DE CONTROL DIGITAL
Ejemplo 3-14
La transformada z con retraso def(t) = e'us(t) se obtiene con la ecuación 3-189 y resulta
00
F(z, A) =
z - e-a
= eaAT e -a(T Anz - 1 + e - a(2T-áT)z-2
(3-191)
Para superar esta dificultad con la transformada z, se introduce un factor m tal que
m = 1 - A (3-192)
Dado que A- tiene un valor entre cero y uno, m también se encuentra en el mismo
intervalo de valores.
Al sustituir la ecuación 3-192 en la 3-189, se tiene que
.00
Z e - amT
= Z — 1 —amT
(3-199)
z — e - aT= Z — °T
(a)
t
—mT
o".
, • c*(t + mT)
n ir.,
r
i
..... -...-,,,
..-
e .4.
e .111 Mi; 1.1
e
c*(t + —7)
1
2r r 4r sr ir
4
7T
+
8T
f
9T
•t
(d)
Figura 3-14. Ilustración de los pasos para obtener la transformada z modificada. (a) Función
del tiempo c(t). c(t) se desplaza a la izquierda m 7' unidades, O < m < 1. (c) El muestreo de
c(t + m7) comienza en t O. (d) La secuencia muestreada se desplaza 7" unidades a la derecha.
En este caso, puesto que la ecuación 3-201 implica que c(0) 0, az,m) = C(z) cuando
m 1.
Las ecuaciones 3 - 194, 3 -200 y 3 -201 presentan tres expresiones alternativas para
la transformada k modificada. Estas ecuaciones están sujetas a diferentes condiciones
de validez y son Útiles para varios fines. La ecuación 3-194 es la más general, ya que
es válida para cualquier función c(t). La ecuación 3-200 lo es para cualquier función
TRANSFORMADA Z 117
c(t) con una transformada de Laplace C(s). La ecuación 3-201 sólo es válida cuando
c(0) = 0, razón por la que no puede emplearse con funciones que tienen un disconti-
nuidad en t = O.
En el siguiente ejemplo se ilustra el empleo de la ecuación 3-200 para evaluar la
transformada z modificada. En d'apéndice C se encuentra una tabla de transformadas
z modificadas, con transformadas z regulares.
Ejemplo 3-15
Considere la funciónf(t) = e-alus(t) del ejemplo 3-14. Al usar la ecuación 3-200 se tiene
1 emT4 j
F(z, m) = z - 1 1[ residuo de at Z = —a
Z + a 1 — e T4z - 1
[ emr4 1 e -otra
= z- 1— (3-203)
1 — e rlz -1 1, z — e-ar
(3-204)
Ejemplo 3-16
El coeficiente del término z-k en la serie infinita representa el valor Uf(t) entre los
instantes de muestreo t (k — 1)T y t = k7; donde k = 1,2, ..., y 0 < m < 1. Durante
118 SISTEMAS DE CONTROL DIGITAL
el, primer periodo de muestreo, la funciónf(t) está descrita por el coeficiente e-amT de
z-1 en la ecuación 3-205. Cuando m = 0, se obtiene el valor def(t) en t 0+; cuando
m= 1, e'T proporciona el valor de c(t) en t = T. De manera similar, para el k-ésimo
periodo de muestreo, k = 1, 2, ...,
m= O
m= 1
Fórmula de inversión
Al igual que con el método de la transformada z, la transformada z modificada inversa
puede obtenerse con la integral de inversión
TRANSFORMADA Z 119
Ejemplo 3-17
Considere la función del ejemplo 3-15, la cual se evaluará con la fórmula de inversión.
Después de sustituir la F(z,m) de la ecuación 3-203 en la 3-206, se tiene que
1 e —mT zk -1 dz
f(kT, m) = — O. (3-207)
2nj r z — e'
donde I" puede ser un círculo en el plano z con centro en el origen y con un radio
suficiente para abarcar el polo del integrando de la ecuación 3-207, que es z =
La evaluación de la integral de inversión (Ec. 3-207) consiste, básicamente, en
encontrar el residuo del integrando en el polo. En consecuencia,
e - mra k— 1
f(kT, m) = E residuo de Z at z = e- aT
z - aT — e
= e —mTa z k — 1 i z=e_.r = e - (k — 1+ oti)aT
(3-208)
la cual puede demostrarse que es la solución en forma cerrada del problema 3-16.
PROBLEMAS
3-1 Determine las transformadas z de las siguientes funciones. Puede utilizar una
tabla de transformadas z o el método de desarrollo en fracciones parciales.
5
F(s)
s(s2 + 4)
4
F(s) —
s2(s + 2)
2
F(s) =
s2 + s + 2
2(s + 1)
F(s) =
s(s + 5)
10
e. F(s) =
s(s2 + s + 2)
3-2 La transformada z de g(t) para T = 1 s es
z(z — 0.2)
G(z) =
4(z — 0.8)(z — 1)
encuentre la secuencia g(kT) para k = 0, 1, 2, ..., 40. Utilice un programa de compu-
tadora como INVZ del paquete DCSP, si es que lo tiene a la mano. ¿Cuál es el valor
final de g(kT) cuando k --> co?
3-3 Dadas las funciones fi(t) = t2us(t) y f2(t) = C2'14(0, utilice el teorema de la
derivada parcial de la transformada z para encontrar la transformada z de ji(t)fi(t).
3-4 Dadas las funciones fi (t) = t2 us(t) y f2(t) = e-2'us(t), aplique el teorema de la
convolución compleja de la transformada z a fin de hallar la transformada z de
fi(t) f2(t).
3-5 Detetrnine la transformada z de las siguientes secuencias. Exprese los resultados
en forma cerrada.
2z + 1
F(z) —
(Z — 0.1)2
3-8 Halle la transformada z inversa f(k) de la siguiente función.
2z
F(z) = z2
1.2z + 0.5
3-9 Dada la función
f (k) = (0.1)kus(k) + 0.5k(0.1)kus(k)
halle la transformada z de f(k), F(z).
3-10 Dada la función
f(k) = (0.1)k u.,(k) + 0.5k(0.1)k- Ius(k — 1)
encuentre la transformada z de f(k), F(z).
3-11 Encuentre la transformada z inversa f(kT) de las siguientes funciones.
z
F(z) =
z2 + 1
10z
F(z) z2 1-
1
c. F(z) =
(z — 0.2)
TRANSFORMADA Z 121
F1(z) = .2[F1(s)]
Si Fi(s) y F2(s) están limitadas en banda con IFi(ico)i = O para 1w1 >_ coi y1 F2(1 (o)! = O
para! a, I co2 y la frecuencia de muestreo es cos >_ 2 máx(col , W2), demuestre que
Y[Fl (s)F2(s)] = TFi(z)F2(z)
3-17 Considere el sistema de datos muestreados de la figura P3-17, el cúal tiene un
retraso que no es un múltiplo entero de la frecuencia de muestreo. Pruebe que la
transformada z de la salida es
00
c[t-(1-m)T] cnt-(1-m)T]
C(I). e—(1—m)Ts
0<m51 T Cm(z) Figura P3-17.
Bibliografía
Muestreo y transformadas z
"Additions to Sampled-Data Theory " de G. V. Lago en Proc. National Electro-
nics Conf., vol. 10, pp. 758-766, 1954.
"The z-transformation " de H. A. Helm en Bell System Technical Journal, vol.
38, pp. 177-196, enero 1959.
"The Analysis of Sampled-Data Systems " de J. R. Ragazzini y L. H. Zadeh en
Trans. ALEE, vol. 71, parte 2, pp. 225-234, 1952.
Kuo, B.C., Analysis and Synthesis of Sampled-Data Control Systems, Prentice-
Hall, Englewood Cliffs, N.J., 1963.
5. "A general z-Transform Formula for Sampled-Data Systems " de E. I. Jury en
IEEE Trans. Automatic Control, vol. AC-12, pp. 606-608, octubre 1967.
TRANSFORMADA Z 123
Jury, E.I., Theory and Applications of the z-Transform Method, John Wiley &
Sons, Nueva York, 1964.
"The z-Transform of a Product of Two Functions" de I. Kliger y W. C. Lipinski
en IEEE Trans. Automatic Control, vol. AC-9, pp. 582-583, octubre 1964,
"Variable Frequency Sampling" de P. M. Will en IRE Trans. Automatic Control,
vol. AC-7, p. 126, octubre 1962.
"Adaptive Sampling Frecuency for Sampled-Data Systems" de R. C. Dorf,
M. C. Farren y C. A. Phillips en IRE Trans. Automatic Control, vol. AC-7, pp.
38-47, enero 1962.
"The z-transform of an impulse function" de R. K. III, Cavin, D. L. Chenoweth
y C. L. Phillips en IEEE Trans. Automatic Control, vol. AC-12, p. 113, febrero
1967.
"A note on Sampled-Data Control Systems" de C. L. Phillips en IEEE Trans.
Automatic Control, vol. AC-10, pp. 489-490, octubre 1965.
"Additions to the z-Transformation Theory for Sarbpled-Data Systems" de G. V.
Lago en Trans. ALEE, vol. 74, parte 2, pp. 403-407, enero 1955.
"Analysis of Cyclic-Rate Sampled-Data Feedback-Control Systems" de R. E.
Hufnagel en Trans. ALEE (Applications and Industry), vol. 77, pp. 421.425,
noviembre 1958.
"The Analysis of Sampled-Data control Systems with a Periodically Time-Var-
ying Sampling Rate" de E. I. Jury y F. J. Mullin en IRE Trans. Automatic Control,
vol. AC-4, pp. 15-21, mayo 1959.
Gabel, R.A. y R.A. Roberts, Signals and Linear Systems, 3a. ed., John Wiley &
Sons, Nueva York, 1987.
Transformada z modificada
"Additions to the Modified z-Transforni Method" de E..I. Jury en IRE WESCON
Convention Record, parte 4, pp. 136-156, 1957.
"Tables of z-Transforms and Módified z-Transforms of Various Sampled-Data
Systems Configurations" de E. I. .11ry y Farmanfarma en Univ. of California,
Berkeley, Electronics Research Lab., Report 136A, Ser. 60, 1955.
18. "A Theorem on the Modified z-Transform" de W. Mesa y C. L. Phillips en IEEE
Trans. Automatic Control; vol. AC-10, p. 489, octubre 1965.
1,5:5 .1.510 1111 I.:11:y
Funciones de
transferencia,
diagramas de bloques
y gráficas de flujo de
señal
TÉRMINOS Y TEMAS IMPORTANTES
Función de transferencia muestreada • función de transferencia en z • ecuación característica •
diagramas de bloqué • gráficas de flujo de señal • sistemas con muestreo múltiple.
4-1 INTRODUCCIÓN
Hasta el momento el estudio de los sistemas de datos discretos y de control digital se
ha concentrado en el tratamiento matemático de las señales maestreadas y de los datos
digitales. Ahora es momento de considerar el modelado y análisis de sistemas sujetos
a entradas discretas y datos digitales.
Todos los conceptos relacionados con el análisis y diseño de sistemas de control
analógicos lineales, como funciones de transferencia, diagramas de bloques y gráficas
de flujo de. señal, pueden aplicarse al estudio de los sistemas lineales de control
digitales y discretos.
124
FUNCIONES DE TRANSFERENCIA, DIAGRAMAS DE BLOQUES Y GRÁFICAS DE FLUJO DE SEÑAL 125
Muestreador ficticio
c* (ti
s2 (s)
Figura 4-1. (a) Sistema lineal con entrada continua. (b) Sistema lineal con entrada muestreada.
que las condiciones iniciales del sistema son todas cero. Si la entrada al sistema es
una señal muestreada o digital r*(t), donde esta última se encuentra simulada por el
muestreador ideal Si de la figura 4-1(b), el objetivo ahora es encontrar las relaciones
entrada- salida del sistema de datos discretos y, si es posible, una función de transfe-
rencia entrada-salida en los dominios s y z. La transfórmada de Laplace de la salida
del sistema de la figura 4-1(b) es
C(s) = R*(s)G(s) • (4-2)
donde R*(s) denota la transformada de Laplace de r*(t), y G(s) está definida por la
ecuación 4-1. Aunque esta función de transferencia es fácil de obtener, contiene
una mezcla de señales discretas y continuas con la que es dificil de trabajar analítica-
mente; por tanto, es deseable expresar las características del sistema mediante una
expresión de transferencia que relacione r*(t) con c*(t), que es la forma muestreada
de la salida c(t). Para este fin se introduce un muestreador ficticio 52 en la salida del
sistema, como se indica en la figura 4-1(b), sincronizado con Si; esto es, los dos
muestreadores tienen el mismo periodo de muestreo T, y realizan el muestreo en los
mismos instantes. Con esto se tiene a disposición tanto a c(t) como a c*(t). De las
ecuaciones,2-53 y 2-54 pueden obtenerse dos expresiones útiles para C*(s) que son
las siguientes;
1
C*(s) --= — E C(s + jnos) c(0) = 0 (4-3)
CO
1 a° 1 °°
C*(s) = — C(s + jntos) = E R*(s + jncos)G(s + jncos) (4-5)
de la figura 4-1(b). La salida del sistema es la respuesta al impulso g(t). La salida del
muestreador ficticio S2 es
00
Muestreador ficticio
— — c*(t)
—
s2 C*(s)
en cascada, G1 y G2. Los dos elementos están separados por un segundo muestreador
52, que es idéntico al de la entrada y que funciona sincronizado con éste. Para facilitar
la aplicación del método de la transformada z, se introduce un muestreador ficticio 52
en la salida del sistema, como se indica en la figura 4-2. Una vez hecho esto, la
transformada z de todo el sistema se obtiene de la siguiente manera.
Las relaciones de transferencia de los sistemas Gi y G2 son
D(z) G 1 (z)R(z) (4-18)
(4-19)
C(z) G2(Z)D(Z)
donde D(z) es la transformada z de d(t), que es la salida de G i . Al sustituir la ecuación
4-18 en la 4-19 se tiene que
C(z) G i(z)G2(z)R(z) (4-20)
Por consiguiente, puede afirmarse que la función de transferencia en z de dos sistemas
lineales separados por un muestreador es igual al producto de las transformadas z
de los dos sistemas.
I.
d(r)
G 1 (s) G2(s)
NY)
Dado que el retenedor de orden cero (ROC) es un elemento importante en los sistemas
de control de datos rnuestreados, en esta sección se estudia con interés especial la
función de transferencia de éste.
La función de transferencia del ROC, dada por la ecuación 2-86, es
1—e '
Gho(s) — (4-25)
(1 — 1) z = 1 (4-27)
z—1
donde
—
.11G1,0(s)G(s)] = e G(s)]
-
= (1 — z 1)11G‘(s)-1 (4-29)
sj
entonces
así como
En consecuencia, los resultados que aparecen en las dos últimas ecuaciones significan
que si se envía una serial continua a un dispositivo S/H, con periodo de muestreo T,
la salida de este último se convierte en r(t) si el periodo de muestreo se hace cero.
Ejemplo 4-1
Como ilustración del proceso de límite en la unidad S/H, suponga que la entrada r(t)
del sistema de la figura 4-4 es
r(t) = e - atus(t) (4-34)
La serie infinita del miembro derecho de la ecuación 4-45 se convierte en G(s) cuando
T tiende a cero ya que, como se indica en la figura 2-34, todos los componentes
armónicos se desplazan hacia el infinito a medida que oh tiene a infinito; por tanto,
la ecuación 4 -45 se convierte en
lím [Gho(s)G(s)] * = G(s) (4-46)
T.-+0
Con esto se demuestra la identidad que aparece en la ecuación 4-40, ya que [Gho(s)G(s)]*
es, en esencia, igual que .1°[Gho(s)G(s)] con excepción de un cambio de variable, z = e7.
Ejemplo 4-2
G(s) (4-47)
s(s + a)
donde K y a son constantes. La transformada z de la combinación retenedor de orden
cero G(s) es
K
2IGho(s)G(s)] = ( 1 — Z—
le[ .5 2 (s + a)]
KT K(1 — e-a')
(4-48)
a(z — 1) a2 (z — °T)
FUNCIONES DE TRANSFERENCIA, DIAGRAMAS DE BLOQUES Y GRÁFICAS DE FLUJO DE SEÑAL 133
E(s) ?S E*(s) el
G(s)
C(s) e».
r(t) 47%.......) e(t) T e*(t) c(t)
-
H(s)
Figura 4-5. Sistema de control de datos discretos con un inuestreador en la trayectoria directa.
134 SISTEMAS DE CONTROL DIGITAL
Nótese que en el miembro izquierdo de las ecuaciones anteriores aparecen sólo las
variables de salida, mientras que las de entrada nada más se encuentran en los
miembros derechos de las expresiones. Ahora, la transformada muestreada de la
ecuación 4-51 es
E (s) = R*(s) — GH*(s)E*(s) (4-53)
donde
GH*(s) [G(s)H(s)]*
-
_ E G(s + jncos)H(s + jncos) (4-54)
T
Al resolver para E*(s) la ecuación 4-53, se tiene que
H(s)
c*(t) T
El lector debe examinar las expresiones anteriores con . cuidado y notar que, además
de las funciones de transferencia del sistema, los miembros derechos de ambas
ecuaciones contienen sólo la señal de entrada, mientras que las variables de salida
aparecen únicamente en los miembros izquierdos. La transformada muestreada de las
ecuaciones 4-61 y 4-62 son
E*(s) = R*(s) — H*(s)C*(s) - (4-63)
C*(s) = GR*(s) — GH*(s)C*(s) (4-64)
donde
GR*(s) = [G(s)R(s)]* (4-65)
136 SISTEMAS DE CONTROL DIGITAL
Y
GH*(s) = [G(s)11(s)]* (4-66)
De la ecuación 4-64, se tiene
C*(s) = (4-67)
1 + GH*(s) GR*(s)
La transformada z de la salida c*(t) se obtiene de manera directa de la ecuación
anterior:
Ejemplo 4-3
Suponga que la función de transferencia directa G(s) del sistema de la figura 4-5 es
10
G(s) — (4-73)
s(s + 5)
10 2(1 — e - 5 T)z
G(z) = GH(z) = 21G(s)H(s)] = .1°(
s.(s + 5)) (z — 1Xz — e_5T)
0.7869z
= (4-74)
(z — 1Xz — 0.'6065)
Al sustituir la ecuación anterior en la 4-60, la función de transferencia de lazo
cerrado es
C(z) 0.7869z
= (4-75)
R(z) z2 — 0.8196z + 0.6065
Por tanto, la ecuación característica del sistema es
z2 — 0.8196z + 0.6065 = O (4-76)
H(s)
(a)
-H(s)
(b)
Figura 4-7. (a) Diagrama de bloques de un sistema de control digital de lazo cerrado. (b) Gráfica
de flujo de señal equivalente del sistema.
—GH*(s)
G*(s)
R*(s) o •1 E*(s)
• 0 C*(s)
Figura 4-8. Gráfica de flujo de señal muestreada del sistema de la figura 4-7.
140 SISTEMAS DE CONTROL DIGITAL
G(s)
C(s) = R (s) (4-80)
1 + GH*(s)
Las relaciones para las salidas discretas dadas por las ecuaciones 4-78 y 4-79 también
pueden obtenerse de la gráfica de flujo de señal compuesta de la figura 4-9 mediante
la fórmula de la ganancia.
El método de la gráfica de flujo de señal muestreada también es aplicable en
sistemas lineales con muestreo múltiple y varios lazos, siempre y cuando los mues-
treadores estén sincronizados y tengan la misma frecuencia de muestreo. El método
también puede emplearse a fin de obtener transformadas z modificadas, de modo que
sea posible recuperar las respuestas entre los instantes de muestreo.
—GH*(s)
G*(s)
R*(s) C*(s)
h"*(s)
G(5)
R(s) C(s)
E(s)
—H(s)
Figura 4-9. Gráfica compuesta de flujo de señal del sistema digital de la figura 4-7.
FUNCIONES DE TRANSFERENCIA, DIAGRAMAS DE BLOQUES Y GRÁFICAS DE FLUJO DE SEÑAL 141
Ejemplo 4-4
H(s)
Gi(s) G2(s)
R(s) C(s)
Figuro 4-10. (a) Diagrama de bloques de un sistema de control de datos discretos con varios
lazos. (b) Gráfica de flujo de señal del sistema.
142 SISTEMAS DE CONTROL DIGITAL
R*(s) C*(s)
R(s) C(s)
•
Figura 4-11. Gráficá compuesta de flujo de señal del sistema de la figura 4-10.
1 + G2H*(s) R*(s)
E*(s) = (4-84)
1 + G i Gt(s) + G2H*(s)
GiGi(s)
C*(s) = R*(s) (4-85)
1 + G i Gt(s) G 2 H*(s)
G i (s)G2 (s)Ci + G2 H*(s)] — G2 (s)H(s)G iGt(s)
C(s) = R*(s) (a-86)
1 + G i Gt(s) + G2H*(s)
1+ G2 H(z)
E(z) = R(z) (4-87)
1 + G 1 G2(z) + G2H(z)
G i G 2(z)
C(z) = R(z) (4-88)
1 + G i G 2(z) + G2 H(z)
Ejemplo 4-5
En este ejemplo se ilustra un sistema de datos discretos que contiene un lazo que no
está separado por muestreadores. La figura 4-12(a) presenta el diagrama de bloques
de este sistema. Los pasos son los siguientes.
Se dibuja la gráfica de flujo de señal del sistema original (Fig. 12b). Las variables
de entrada y salida son
variables de entrada: r(t), e*(t)
variables de salida: e(t), c(t).
Note que la gráfica de flujo de señal tiene un lazo analógico con una ganancia de
lazo igual a -G2(s). El bloque D(s), el cual cuenta con muestreadores en la
entrada y la salida, es la representación en diagrama de bloque de un controlador
digital con entrada y salida discreta. La función do transferencia D(s) junto con
el muestreador de la salida, pueden reemplazarse con una función de transferencia
muestreada D*(s), según se indica en la figura 4-12.
Las siguientes ecuaciones se obtienen a partir del diagrama de la figura 4-12(a)
con la fórmula de Mason para la ganancia, donde r(t) y e*(t) son las entradas y
e(t) y c(t)las salidas, tomando al mismo tiempo la transformada muestreada:
E*(s) ( 1 +R(Gs) (s) )* D*(s)( G i(s)G 2(s))* E*(s)
(4-89)
1 + G 2(s)
c*(s) GR(±s)GG2((ss)))* D*(s)( G i (s)G 2(s)\*
(4-90)
1 + G2(s)j E.
Note que en este ejemplo, que es similar al del sistema de la figura 4-6, la
señal de entrada no aparece como una variable independiente, ya que ésta no se
muestrea de manera directa.
C(s)
G2(s)
(b)
Figura 4-12. (a) Diagrama de bloques de un sistema de control de datos discretos. (b) Gráfica
de flujo de señal.
144 SISTEMAS DE CONTROL DIGITAL
-1
Figura 4-13. Gráfica compuesta de flujo de señal del sistema digital de la figura 4-12.
(G i (s)G 2(s)*
D*(s)
(
1(+s)GG2(s))*
(R 1 + G 2(s) ) R(s) \*
C* (s) — (4-91)
(G i(s)G 2(s)* + G 2(s))
)1 + 1 + D* (s)
1 + G2(s))
o, en el dominio de la transformada z,
1 1(-1-s):22(:))(z) G 2(s)
donde
R(s)G 2(s)) (R(s)G 2(s))
, (z) (4-93)
( 1 + G 2(s) 1 + G2(s))
(Gi(s)Gz(s))D*(s)( R(s) )*
2( s)
R(s)G 2(s) 1 + G 2 (s) ) - 1 G2(s)
C(s) = + (4-94)
1 + G 2(s) G i(S)G2(S))*
1 + D*(S)(
1 + G 2(s)
--
-x c*(1)
- —C(z)
T --0k
E(s) E*(s) c(t)
G(s) •
e(t) T e*(t) C(z)
—
"X c*(t)
T
- —C(z)
— —► T
c*(t—AT)
- -7;
1-(z/n)
E(s) E*(s)♦ e—áTs
G(s)
e(t) T e*(1) C(s) c(t—M)
c(1) Retraso
ficticio
Figura 4-14. (a) Sistema de datos discretos. (b) Sistema de datos discretos con retraso ficticio.
146 SISTEMAS DE CONTROL DIGITAL
C(z, 1.= en
T, ,
= E(z)G(z, m) (4-97)
Por otra parte, la transformada z modificada de una señal muestreada e*(t) es igual
a la transformada z de la señal, o
„,[E*(s)] = E(z) (4-100)
Ejemplo 4-6
1
G(s) = (4-101)
s+a
donde a es constante. La entrada al sistema es una función escalón unitario, e(t) =
u.s(t). El problema consiste en evaluar la salida del sistema con el método de la
transformada z modificada.
E(z) = (4-103)
z—1
mientras que la transformada z modificada de G(s) está dada por la ecuación 3-203,
e'
G(z, m) (4.104)
z — e - aT
En consecuencia,
e - maT
C(z, m) — (4-105)
z — e —ar Z — 1
m=1 m= i
m= 1
T 2T 3T 4T 5T 6T 7T t
Figura 4-15. Respuesta del sistema de datos discretos del ejemplo 4-5.
c*(t) c*(t--.117)
r -?5 - "-C-7(Znil
T C(z)
R(s) E(s) E*(s) C(s) i
s) e-ATs
+ T c(t) o-67)
Retraso
ficticio
H(s)
Figura 4-16. Sistema de datos discretos de lazo cernido con retraso ficticio.
FUNCIONES DE TRANSFERENCIA, DIAGRAMAS DE BLOQUES Y GRÁFICAS DE FLUJO DE SEÑAL 149
Los sistemas de datos discretos considerados hasta el momento tienen una frecuencia
de muestreo uniforme. Cuando los datos son digitales, la rapidez con que fluyen en
las entradas y salidas de los procesadores y controladores digitales de todo el sistema
es idéntica y uniforme. Para sistemas de datos muestreados, los muestreadores que
hay en todo el sistema tienen la misma frecuencia de muestreo.
En la práctica, el ancho de banda de cada una de las partes que forma el sistema
puede ser muy diferente. En consecuencia, resulta más apropiado muestrear una señal
lenta con una frecuencia de muestreo baja, mientras que las señales que tienen
componentes en frecuencia altos, deben muestrearse con frecuencias de muestreo
altas. Cuando un sistema de datos discretos tiene frecuencias de muestreo diferentes
en distintos puntos del mismo, el sistema recibe el nombre de sistema con muestreo
múltiple. La figura 4-17 presenta el diagrama de bloques de un sistema de control
automático para proyectiles. Como consecuencia de la diferencia entre las dinámicas
de los lazos de posición y retroalimentación, los muestreadores del sistema tienen dos
periodos de muestreo, Ti y T2. Nótese que los dos retenedores de datos que hay en el
sistema ofrecen características distintas, ya que ambos retienen datos muestreados
lapsos Ti y T2, respectivamente.
Dada la gran variedad de esquemas de muestreo, el análisis y diseño de sistemas
con muestreo múltiple es, en general, muy complicado. Para comenzar, primero se
establecerán los esquemas fundamentales de muestreo múltiple. La figura 4-18
presenta tres tipos básicos de sistemas discretos con muestreo múltiple. El sistema de
la figura 4-18(a) tiene un muestreador lento-rápido en la entrada y otro rápido-lento
en la salida. Este tipo de sistema se conoce como sistema con muestreo múltiple
lento-rápido. La figura 4-18(b) presenta un sistema donde el muestreador de la entrada
cuenta con una frecuencia de muestreo mayor que la del muestreador de la salida; este
sistema recibe el nombre de sistema con muestreo múltiple rápido-lento. Un sistema
de control digital que tiene controladores con diferentes frecuencias de muestreo en
la entrada y la salida, puede modelarse con el diagrama de bloques que aparece en la
figura 4-18(c). En general, resulta imposible relacionar los periodos de muestreo Ti
y T2 mediante un factor multiplicativo entero. Para fines de análisis, las matemáticas
se simplifican mucho si se supone que Ti y T2 son múltiplos enteros; esto es, con
respecto a la figura 4-18(a), Tl = N7'2, donde N es un entero positivo mayor que uno.
Dinámica
HRetenedor Retenedor Actuador 1--:2 del
T2 proyectil
Controlador Sensor
digital T2
Controlador
digital Sensor
Figura 4-17. Sistema de control digital del.piloto automático de un proyectil con muestreadores
que trabajan con frecuencias distintas.
150 SISTEMAS DE CONTROL DIGITAL
Entrada ?S Salida
G(s)
(a)
Entrada
G(s)
s
? Salida
T2
(b)
Entrada ?S Salida
G(z)
T2
(c)
Figura 4-18. Sistemas de datos discretos con muestreo múltiple y de lazo abierto. (a) Sistema
con muestreo lento-rápido (Ti > T2). (b) .Sistema con muestreo rápido-lento (Ti > T2). (c)
Sistema completamente digital.
G(z)N
k=0
E0. g(kT/ N)z (4-114)
R(.․)?5 s G(s)
c(t) c*(ON
T R(z) T/N I TIN C(z)N
Muestreador
ficticio
Figura 4-19. Sistema con muestreo múltiple lento-rápido con muestreador ficticio rápido.
FUNCIONES DE TRANSFERENCIA, DIAGRAMAS DE BLOQUES Y GRÁFICAS DE FLUJO DE SEÑAL 151
Ti
ersfiV
(t)
G(s)
T R(z) e kTsIN
(N-1)7'0N
Figura 4-20. Sistema con muestreo múltiple lento-rápido donde el muestreador rápido se
descompone en N muestreadores lentos con unidades de retraso y adelanto.
y C(z)N se define del mismo modo. Al igual que en la ecuación 3-187, G(z)N puede
obtenerse a partir de G(z) con la siguiente transformación
G(z)N = G(z)I
•
(4-115)
= T IN
2TekTsIN G(s)] =
n = O
gRn + —Nk )Tiz (4-119)
Dado que en esta aplicación kIN siempre es menor que uno, al comparar el
miembro derecho de la última expresión con la ecuación que define a la transformada
z modificada, ecuación 3-194, se llega a la siguiente relación:
152 SISTEMAS DE CONTROL DIGITAL
X[ekTsING(s)] =-- [zG(z, n)]im=kIN (4-120)
Ejemplo 4-7
Suponga que la función de transferencia del sistema con muestreo múltiple lento-rá-
pido de 1a figura 4-18(a) es
1
G(s) = (4-124)
s(s + 1)
La entrada del sistema es una función escalón unitario 140. Los periodos de muestreo
son T1 = 1 s y T2 = 1/2 s.
La transformada en z modificada de G(s) para un periodo de muestreo de 1 s es
1 e -m
G(z, m) = (4-125)
z- 1 z - e-1
Tabla 4-1 c *(t)2 y c*(t) del sistema con muestreo múltiple del ejemplo 4-7
c*(02 c*(t)
Periodos de muestreo Valores Periodos de muestreo Valores
O 0.000000 0 0.000000
1 0.393000 1 0.632000
2 0.631551 2 1.496576
3 1.169351 3 2.446740
4 1.495796 4 3.428401
5 2.086948 5 4.421651
6 2.445778 6 5.419168
7 3.056587 7 6.418253
8 3.427348 8 7.417916
9 4.045400 9 8.417791
10 4.420558 10 9.417746
11 5.041279 11 10.417728
12 5.418056 12 11.411722
13 6.039761 13 12.417719
14 6.417135 14 13.417717
15 7.039202 15 14.417716
16 7.416796 16 15.417716
17 8.038996 17 15.417715
18 8.416671 18 17.417713
19 9.038919 19 18.417711
20 9.416624
21 10.038890
22 10.416605
23 11.038877
24 11.416596
25 12.038872
26 12.416593
27 13.038869
28 13.416591
29 U.038868
30 14.416591
31 15.038868
32 15.416590
33 16.038866
34 16.416588
35 17.038864
36 17.416586
37 18.038862
38 18.416584
39 19.038860
Después de cancelar los factores comunes en la última ecuación, C(z)2 queda como
,2
C(2)2 = 0.393z2
L(z 2 - ixz, - 0.607)12j - 1
(4-129)
La expresión entre paréntesis es G(z)2, misma que puede obtenerse a partir de G(z)
con la ecuación 4-115.
El desarrollo en serie de potencias en zz l mediante división larga, o con ayuda de
un paquete de software, permite obtener c(t) en los instantes de muestreo kT2. La tabla
4-1 contiene la respuesta para k = O hasta 39. El lector puede verificar que si la
frecuencia de muestreo del muestreador de salida de la figura 4-18(a) tiene un periodo
= 1 s, la transformada z de la salida del sistema es
1 z
C(z) =
0.632z
1.368z + 0.3681 — 1
L2 (4-130)
Por tanto, C*(s)N puede expresarse del mismo modo. Las relaciones anteriores son
útiles para el análisis en el dominio de la frecuencia. En el siguiente desarrollo se
presenta otra aplicación de la ecuación 4-132.
donde G(z)N está definida en las ecuaciones 4-114 y 4-115, mientras que R(z)N está
definida de manera similar. Una vez que se determina C(z)N, que es la transformada.
Z de la salida con muestreo rápido c*(t)N, la salida del muestreador ficticio lento c*(t)
se obtiene al "extraer" muestras de c*(ON en t = kT, k= 0, 1, 2,..., o al obtener c(k7)
a partir del conjunto c(kT/N). Luego entonces, C(z) es simplemente la transformada z
de c*(t).
A continuación se presentan tres métodos distintos para obtener C(z) a partir de
C(z)N.
(a) Método del dominio de la frecuencia
Primero se escribe
*(S) =
T
C[s j 21 1
Tj
(4-134)
Sea n = pN+ k, donde —00 .p 00, y k = 0, 1, 2, ..., N— 1, entonces la última ecuación
se convierte
1 2irpN 2nk
C*(s) = C[s +j +j (4-135)
I T Ti
la cual puede escribirse como
1 1 T N j27rP N +. 2;k
1
C*(s) (4-136)
En vista de la relación para G*(s)N dada por la ecuación 4-132, la última expresión
queda como
1 N - 1 [ 27rk
C*(s) = C* [s + -- (4-137)
N k = 6 T
156 SISTEMAS DE CONTROL DIGITAL
Como k es entero, R*(s + firkIT)= R*(s); por tanto, la ecuación 4-138 se convierte
en
C*(s) = G*(s)R*(s) (4-139)
C(z) G(z)R(z) (4-140)
donde
1 03
dzN
= C(z N) E zÑ z " (4-146)
r
n=o ZN
Ejemplo 4-8
Suponga que la función de transferencia G(s) del sistema de la figura 4-21 es la dada
por la ecuación 4-124, y que T = 1 s, N = 2. La entrada r(t) es una función escalón
unitario. La transfonnada z de la salida del muestreador rápido es C(z)N = G(z)NR(z)N;
por tanto,
0.393z2 z2
C(z) 2 — (4-149)
[(z2 — 1)(z2 — 0.607) z2 — 1
En este ejemplo se utilizan los tres métodos delineados con anterioridad para encon-
trar la transformada z de la salida, C(z). La transformada z de la salida del muestreador
lento c*(t) se obtiene al sustituir la última ecuación en la 4-141. Con ello se tiene
1 1 c( z2 ej2nk/2)2
C(z) = -
2 kír0
1 0.393z5 0.393z1
(4-150)
2 (z 2 — 1) 2 (z 2 — 0.607) (—z 2 — 1)2( — z 2 — 0.607)
158 SISTEMAS DE CONTROL DIGITAL
Por tanto,
1.5445z 1.5445z 2 - 3.5445z (1.0246z + 0.2386)z
C(z) = (4-159)
z - 0.368 (z - 1)2 (z - 1)2(z - 0.368)
que es el resultado esperado.
Las respuestas c*(t)2 y c*(t) se obtienen mediante el desarrollo de C(z) 2 y C(z)
en series de potencias en z1-1 y z-1 , respectivamente. Los resultados aparecen en la
tabla 4-2.
FUNCIONES DE TRANSFERENCIA, DIAGRAMAS DE BLOQUES Y GRÁFICAS DE FLUJO DE SEÑAL 159
Tabla 4-2 c*(t)2 y c*(t) del sistema con muestreo múltiple del ejemplo 4-8
c*(t)2 c*(t)
Periodos de muestreo Valores Periodos de muestreo Valores
0 0.000000 0 0.000000
1 0.393000 1 1.024600
2 1.Ó24551 2 2.664853
3 1.800902 3 4.531666
4 2.665147 4 6.481853
5 3.582742 5 8.462722
6 4.532722 6 10.454881
7 5.502357 7 12.451195
8 6.483922 8 14.449038
9 7.472728 9 16.447441
10 8.465928 10 18.446051
11 9.461794 11 20.444735
12 10.459277 12 22.443445
13 11.457741 13 24.442163
14 12.456799 14 26.440884
15 13.456215 15 28.439604
16 14.455849 16 30.438324
17 15.455612 17 32.437046
18 16.455454 18 34.435768
19 17.455341 19 36.434490
20 18.455256
21 19.455183
22 20.455116
23 21.455050
24 22.454979
25 23.454903
26 24.454821
27 25.454733
28 26.454639
29 27.454538
30 28.454430
31 29.454315
32 30.454193
33 31.454063
34 32.453926
35 33.453777
36 34.453617
37 35.453445
38_ 36.453262
39 37.453068
160 SISTEMAS DE CONTROL DIGITAL
Ti
enIN
r(t)
R(s) ekTs/N
(N-1)Ts1N
Figura 4-22. Sistema con muestreo múltiple rápido-lento donde el muestreador rápido se
descompone en N muestreadores lentos con unidades de retraso y adelanto.
Por tanto,
[G(s)e kTsi n = G(z, m)I„, = 1 - /o/ (4-163)
donde kIN < 1. Esto siginifica que el segundo término de la ecuación 4-160 puede
obtenerse al hacer m = 1— k/N en la transformada z modificada de G(s).
La ecuación 4-160 queda ahora como
- k
C(z) = R(z)G(z) + NE zR(z, —)G(z, 1 — (4-164)
k = O N
Ejemplo 4-9
Suponga que la función de transferencia G(s) del sistema de la figura 4-21 es igual a
la dada por la ecuación 4-124, y que T = 1 s, N = 2. La entrada r(t) es un escalón
FUNCIONES DE TRANSFERENCIA, DIAGRAMAS DE BLOQUES Y GRÁFICAS DE FLUJO DE SEÑAL 161
1 e - mT
e G(z, m) = (4-167)
z — 1 z — e -T
La sustitución de R(z), G(z) y de las dos últimas expresiones en la ecuación 4-165,
proporciona el siguiente resultado
e - T12
1 C(z) = z--T z1 z—
z— 1 z— 1 z—e 1 z—e-Ti
(1.0246z + 0.2386)z
(4-168)
(z — 1)2(z — 0.368)
que concuerda con el obtenido en la ecuación 4-159.
4-7-3 Sistemas con muestreo múltiple donde todos los elementos son
digitales
Es común tener sistemas digitales con muestreo múltiple en que varios componentes son
completamente digitales, donde las entradas y salidas de éstos se encuentran codificadas
de manera digital, razón por la que no hay muestreo fisico en ellas. A continuación se
demuestra que las relaciones matemáticas obtenidas en las secciones anteriores para
sistemas con muestreo múltiple también pueden aplicarse en este caso.
La figura 4-23 representa el diagranía de bloques de un sistema con muestreo
múltiple completamente digital con un muestreo de entrada lento y. otro de salida
rápido. En este caso, la entrada al sistema está formada por datos digitales, repre-
sentados por trenes de pulsos con periodo T s. De manera similar, los datos dé salida
también están representados por un tren de pulsos con periodo T/N. La función de
transferencia en z entrada-salida es
C(z)N = G(z)N R(z) (4-169)
donde C(z)N y G(z)N están definidos sólo de la siguiente manera
R(z)_01
G(z)N
C(z)N
►
Figura 4-23. Sistema completamente digital con muestreo múltiple lento-rápido.
162 SISTEMAS DE CONTROL DIGITAL
R(z)N C(z)
G(z)N
donde C(z) también puede obtenerse a partir de C(z)N con la ecuación 4-141.
La figura 4-24 contiene la representación en diagrama de bloques de un sistema
con muestreo múltiple completamente digital del tipo rápido-lento. En este caso,
con lo que la salida, que tiene una frecuencia baja, está dada por
1 N-1
C(Z) = — E G(zN ei2nmIN)NR(zNei2xmlN)1N (4-173)
N „,=.0
En esta sección se aplican algunos de los métodos de análisis de sistemas con muestreo
múltiple introducidos en las secciones anteriores a los sistemas de control de lazo
cerrado.
Ejemplo 4-10
Considere el sistema de control digital con muestreo múltiple de lazo cerrado que
aparece en la figura 4-25(a). El muestreador rápido se descompone en N muestteado-
res lentos con unidades de retraso y adelanto, como se indica en la figura 4-25(b).
Note que el retenedor de orden cero (ROC) del muestreador rápido tiene una función
de transferencia diferente a la del muestreador lento. Sin importar cuál sea el caso, el
periodo de retención del ROC debe ser igual al periodo de muestreo.
La función de transferencia del ROC del muestreador con un periodo de mues-
treo T es
1 — e - Ts
GhO(S) = (4-174)
s
La función de transferencia del ROC asociado con el muestreador que tiene un periodo
de muestreo T/2 es
1 — e - Ts12
Gh02(S)1"-- (4-175)
FUNCIONES DE TRANSFERENCIA, DIAGRAMAS DE BLOQUES Y GRÁFICAS DE FLUJO DE. SEÑAL 163
s+1
zoh
T!2
C(z)
R(z)
Y[Gho(s)G(s)]
1 + [ G ho 2 (s)11(s)].21G ho(s)G(s)] + Ile rs12 Gh o(s)GWine Ts12Gh02(s)11(S)]
(4-176)
donde
1
[G ho(s)G( s)] = ( 1 —1)511=
s
(4-177)1
—
—TsI2 1
[G hoz(S)HCSU = 111
s s+ 1
_ 1 re— s12
— (4-178)
s(s + 1)] s(s + 1)]
El último término de la expresión anterior se evalúa aplicando la fórmula 4-161; en
consecuencia,
[ e - TsI2 1 1 1 1 e - T/2
2' — «1`,„ (4-179)
s(s + 1)] [s(s + 1) j„,= 1/2 z — 1 z — e - T
Después de sustituir la ecuación 4-179 en la 4-178 y efectuar las simplificaciones
respectivas, se tiene
i oTs/2
►
Y[eTsI2G o(s)G(s)] eTsI2 1 — e—rs
— (1 z 1)2 1I--] (4-180)
S2
164 SISTEMAS DE CONTROL DIGITAL
e Ts12 1 zT/2 Tz
(4-181)
s 2 ]. = 1
"[ 2]
s m=1/2 Z — 1 ( z — 1)2
1 — Ts12 1
Ts12 G hi (s)H(s)] = e- TsI2
s s+ lj
e - Ts12 [ e - Ts 1
=
Ls(s + 1) .s(s + 1)]
1— e - 12
(4-182)
Z—e
—T
C(z) T(z — e - T)
R(z) 2 T T e - T12 z e - T T T
z —1+ Te - T e - 172 —
2 2 2 2
(4-183)
Ejemplo 4-11
1 — e-Ts/3
G h0 3( S) = (4-185)
s
Las salidas de los muestreadores se denotan como X( *) , Yik e 1T, como
se indica en la figura 4-26(c). Es importante considerar a todas como "entradas," ya
que así pueden descartarse los muestreadores. La salida c* (02 junto con la transfor-
mada z de ésta, C(z)2, describen a c(t) en t = O, T/2, 3T/2,...
FUNCIONES DE TRANSFERENCIA, DIAGRAMAS DE BLOQUES Y GRÁFICAS DE FLUJO DE SEÑAL 165
1- 25 c*(t)Ni
1 Si
c(t)
C(s)
3T 2T ST 3T
t
2 2
R(s) • E(s)
(c)
Figura 4-26. (a) Sistema de control con muestreo múltiple. (b) Esquemas de muestreo de los
muestreadores. (c) Sistema con descomposición de muestreadores equivalente. .
donde
Ft(s) = [e- TsI2 Gh02 (s)G(s)]* [Gh o 3 H(s)]*
+ [e - Ts/6 Gh02(s) G(s)]* [ e Ts13 Gh03(S)H(S)J*
eTsI6 2TsI3G (4-195)
L Gho2( s) G(s)] * Ce- h03(s) 11 (s)] *
La C(z) de la ecuación 4-197 proporciona información sólo sobre c(t) en los instantes
de muestreo t kT, k = 0, 1, 2, 3,... Para describir c(t) en los instantes de muestreo
correspondientes a Si, en la ecuación 4-187 se sustituyen las expresiones correspondien-
tes aXd ydri*, dadas por las ecuaciones 4-193 y 4-194, y luego se obtiene la transformada
z; con esto se llega a
1/2 2(eTs12,-,h02 G)] [1 + F
(G ho2 G) + z- 1(z)] — Ts12 G h0 2 G)
+z 1 1(e
rs12 G hoz G)V 2(z)
C(z)2 = R(z)
1 + F1 (z) + F2(z)
10(eTs/2Gh02 G [1 F2(z)]
[2 ( e- Tsi 2 G hoz G) + z )] +
1/2.2.( _TsI2
C2' (G1,0 2 G) + z- Gh02 G)Vi(z) "eTsI2R)
1 + Fi (z) + F2(z)
(4-209)
donde
— e - Ts/IVI
G o(s) = G(s) (4-213)
s
e- Ts( 1 — 1 iN„)
G„(s). e - Tsfiv„ 1 — G(s) (4-215)
s
Con la C(s) dada por la ecuación 4-210, se obtiene C(z) al tomar la transformada
z de ambos miembros de la ecuación. Sin embargo, C(z) brinda información sobre c(t)
sólo en t = kT, k = 0, 1, 2,... Para obtener información adicional sobre c(t) entre los
instantes básicos de muestreo, puede colocarse un muestreador ficticio S, en la salida,
o aplicar a C(s) el método de la transformada z modificada para obtener C(z, m).
FUNCIONES DE TRANSFERENCIA, DIAGRAMAS DE BLOQUES Y GRÁFICAS DE FLUJO DE SEÑAL 169
r
T c*(t)
R(s) E(s) E*(s)
zoh n—► D(s) zoh --O> G(s)
+ Sc c(t)
C(s)
(a)
• •
o t
•
T T T T T T
O Ni _
N2 N3 T
T+— T+—
NI N2
T+
Ñ3 Zr
(b)
ersiNt
ers/N2
eTs/N3
(c)
(d)
Figura 4-27. (a) Sistema de control con muestreo cíclico. (b) Esquemas de muestreo de los
muestreadores. (c) Descomposición de los muestreadores con retrasos y adelantos. (d) Sistema
equivalente con la descomposición de muestreadores.
170 SISTEMAS DE CONTROL DIGITAL
Ejemplo 4-12
c(t)
ROC G(s)
C(s)
G(s)= s
T=0.1s
1 1 t
—
T+N T 3T
2T 2T + —
Ñ. N
eTs/N
T T.
R(s E(s)
G(s) c(t)
C(s)
'CID( XT
.► Ts/N y — e—TsIN
T=0.18
G(s) = 77,1
(c)
Figura 4-28. (a) Sistema de control con muestreo cíclico. (b) Esquema de muestreo. (c) Sistema
equivalente con la descomposición de muestreadores.
FUNCIONES DE TRANSFERENCIA, DIAGRAMAS DE BLOQUES Y GRÁFICAS DE FLUJO DE SEÑAL 171
[e - T'IN G i(s)]
2[ersIN R(s)] (4-222)
A(z)
G 0 (z)[1 + G 1(z)] + [eTsIN Go(s)]{z - 'IN — G i(s)] R(z)
C(z)2 A(z)
t[e- Ts /NG1(s)](1 + z ' IN {1 + G 0(z) + I[eT" Go(s)]}) 2,[er R(s)]
A(z)
(4-223)
G4,5
(e-TsIN)*(eTsING0)s
Figura 4-29. Gráfica de flujo de señal muestreada del sistema con muestreo cíclico de la
figura 4-28.
172 SISTEMAS DE CONTROL DIGITAL
donde
L(z) = 1 + Go(z) + G 1 (z) + G 0(z)G i (z) — [eT G 0(s)P[e - sl N G i(s)]
(4-224)
Las funciones C(z) y C(z)2 se obtienen con las ecuaciones 4-161 y 4-163,
_ e - TsIN e- T(1 — 1IN) e - T 0.0307
G°(z ) = (4-225)
s(s + 1)) z — e - T z — 0.9048
Las transformadas de las salidas se obtienen al sustituir las últimas cinco expre-
siones en las ecuaciones 4-222 y 4-223 y efectuar las simplificaciones respectivas,
0.0931z2 — 0.0842z
C(z) = z3 (4-230)
2.7186z 2 + 2.455z — 0.7364
Tabla 4-3 c*(02 del sistema con muestreo múltiple del ejemplo 4-12
PROBLEMAS
4-4 Encuentre la salida c(k7) en los instantes de muestreo del sistema de la figu-
ra P4-4.
r(t) ?s r*(t) h(t) c(t)
Retenedor de
primer orden s+1 Figura P44.
(a)
c(kT)
C(z)
(b)
Figura P4-5 (a), (b).
FUNCIONES DE TRANSFERENCIA, DIAGRAMAS DE BLOQUES Y GRÁFICAS DE FLUJO DE SEÑAL 175
K,
Retroalimen Retroalimentación
tación de velocidad de corriente
K,
Figura P4-6.
T= 0.01 s.
T = 0.1 s.
Los parámetros del sistema son:
Resistencia del motor R = 5 -2
Inductancia del motor L OH
Inercia del motor y la carga J= 0.001 oz-in-s2
Relación de engranes N = 100
Constante de par del motor Ka = 10 oz-in/A
Constante de fuerza contraelectromotriz del motor Kb = 0.07 V/rad/s
Ganancia del tacómetro en retroalimentación Kt = 0.0005 V/rad/s
Ganancia en corriente en retroalimentación = 1 V/A
Ganancia del amplificador K= 20
4-7 Repita el problema 4-6 con L = 0.05_11. Los demás parámetros del sistema son
los mismos del problema 4-6.
4-8 En la figura P4-8 aparece una variación del sistema de control de vuelo descrito
en el problema 4-6. Los parámetros del sistema son los del problema 4-6. Obtenga la
FCEM
del motor
Kb
Amplificador
T,,, eo
e 0* 00-n
-441 K Ls+R a Ja
Ki
Retroalimen- Retroalimentación
tación de velocidad de corriente
K
Figura P4-8.
176 SISTEMAS DE CONTROL DIGITAL
4-9 Repita el problemá 4-8 con L = 0.05 H. Los parámetros del sistema aparecen en
el problema 4-6.
Retroalimentación de velocidad
Figura P4-10.
Td
r(s) x e',0(t)
—O» ROC
s(1 +0.023)
cc
n•► e-0.1s
Te • 0.4
1 + 0.1s
_
Motor-actuador Retraso Máquina
Figura P4-11.
FUNCIONES DE TRANSFERENCIA, DIAGRAMAS DE BLOQUES Y GRÁFICAS DE FLUJO DE SEÑAL 177
Figura P4-12.
c(t)
r(t) -X t r(t) Dispositivo rd(t) G(s)
T retenedor
L Figura P4-13.
178 SISTEMAS DE CONTROL DIGITAL
4-14 Para los sistemas de control digital de la figura P4-14, determine C(s) y C(z)
con el método de la gráfica de flujo de señal muestreada. No intente resolver los
problemas de manera algebraica.
R(s) C(s)!
G(s)
H1(s) D*(s)
H2(s) .4
(b)
C(s)
R(s)
G 1 (s) G2 (s)
H(s)
(d) Figura P4-14 (a), (b), (c), (d), (e).
FUNCIONES DE TRANSFERENCIA, DIAGRAMAS DE BLOQUES Y GRÁFICAS DE FLUJO DE SEÑAL 179
H(s)
e*
«I> ROC D(s) —?5 — ROC G(s) •••n•n•n ••
TIN
Figuro P4-15.
s + 3 5
G2 (s) —
Gi(s) s + 5 s(s + 2)
Figura P4-16.
Figura P4-17.
180 SISTEMAS DE CONTROL DIGITAL
zN = Z1IN =
WN +1
WN — 1
2C(w N) 1
C(w) = E residuos de
WN + 1 NW-14—
(wN - 1 w+1
en los polos de C(wN)/(wN + 1). Dada C(zN) = ZAT/(ZN — e-9, halle C(w) con el método
de transformación descrito líneas arriba.
c
ROC G(s)
Muestreador
no uniforme
t tt ft tt tt••• •t
IN / 2T T 3T \47.
+ 1 IN)T + 1 IN)T (3 + 1IN)T Figura P4-20.
FUNCIONES DE TRANSFERENCIA, DIAGRAMAS DE BLOQUES Y GRÁFICAS DE FLUJO DE SEÑAL 181
T, (1 + 1/N)T, 2T, (2 + 11IV)T, 3T, ..., donde N es un entero positivo mayor que uno.
Demuestre que la transformada z de la salida puede escribirse como
1
C(z) = {[G 0(z)[1 + G i(z)] --- [ersING0(s)] lle-niNG i(s)nR(z)
donde
e e - Ts(1 - IN))
G 1 (z) = G(s).
Encuentre C(z) para G(s) = K/(s + 1), r(t) = 74(0 y T=1 s. No se especifican los valores
de K y N.
4-21 Existe un tipo especial de esquema con muestreo múltiple donde se supone que
la entrada y la salida del sistema son trenes de pulsos, donde los pulsos aparecen en
t= kT/N, k = 0, 1, 2, ..., con excepción de que el valor de los pulsos de entrada cambia
sólo en t = kT, k = 0, 1, 2,... La figura P4-21 presenta el diagrama de bloques del
sistema, donde Si representa el muestreador especial de la entrada. El muestreador de
la salida es un muestreador regular con periodo de muestreo TIN. Demuestre que el
muestreador Si puede descomponerse en una combinación de muestreador con
periodo T y N retrasos conectados en p-aralelo. Dibuje el diagrama de bloques
equivalente del sistema. Obtenga la transformada z de lá salida c*(t)N, C(zN) en
términos de la función de transferencia del sistema y de la transformada de la
entrada R(z).
r(t)
G(s) TIN--o' c*(01,1 Figura P4-21.
Sl
r(t)
G(zN) C(zN)
Si TIN Figura P4-22.
182 SISTEMAS DE CONTROL DIGITAL
Figuro P4-23.
r(t) c(s)
-I> ROC D(s) ROC G(s)
Figura P4-24.
Bibliografía
5-1 INTRODUCCIÓN
En el capítulo anterior se presentaron métodos convencionales para el análisis de
sistemas de control de datos discretos mediante transformadas z, funciones de trans-
ferencia, diagramas de bloque o gráficas de flujo de señal. Como alternativa también
puede utilizarse el enfoque calificado como moderno, el cual emplea la formulación
con variables de estado. En la práctica, muchos ingenieros se sienten más cómodos
con los métodos convencionales, ya que son más fáciles de comprender y tienen el
apoyo de una buena cantidad de técnicas gráficas. Por otra parte, las técnicas con
variables de estado poseen las siguientes ventajas sobre los métodos convencionales,
al menos en lo que respecta al estudio de sistemas de control de datos discretos.
184
TÉCNICAS CON VARIABLES DE ESTADO 185
Ecuación de salida
c(t) = g[x(t), u(t), t] (5-4)
donde x(t) es una matriz columna de dimensión n x 1 que se conoce como vector de
estados; esto es,
x 1(0
x(t) =[x2(t) (5-5)
x,.,(t)
u p(t)
ct)
que es una matriz columna o vector de dimensión q x 1.
Si el sistema es lineal e invariable en el tiempo, las ecuaciones dinámicas pueden
escribirse como
dx(t)
— Ax(t) + Bu(t) (5-8)
dt
y
c(t) = Dx(t) + Eu(t) (5-9)
Sea x(0) el estado inicial en t = 0. Entonces, 4(t) también está definida por la siguiente
ecuación matricial
x(t) = ol)(07[(0) (5-12)
Otro método para resolver la ecuación de estados homogénea es suponer una solución,
igual que con el método clásico de solución de ecuaciones diferenciales ordinarias
lineales, y sustituirla en la ecuación de estados a fin de comprobar si es correcta.
Supóngase que la solución de la ecuación 5-10 es
x(t) = eAtx(0) (5-17)
para t O, donde eM representa un desarrollo en serie de potencias de la matriz At,
e At = 1 1
At + — A2t 2 + — A 3t3 + ••• (5-18)
2! 3!
La derivada de la última ecuación con respecto a t es
d(eAr) A A,
, (5-19)
dt
que tiene la misma forma que la ecuación 5-11.
Además de la ecuación 5-16, se tiene otra expresión útil para la matriz de
transición de estados,
«o = eAt = I ± At + 1
ít A 2t2 + 1
—
31 A3t3 + • • • (5-20)
para kT t S (k + 1)T. La ultima ecuación describe el vector de estados x(t) para todo
tiempo t entre los instantes de muestreo kt y (k + 1) T, con k = 0, 1, 2, ...
l Se
Las ecuaciones de estado discretas dadas por la expresión 5-39, junto con la ecuación
de salida de la 5-43, se definen como las ecuaciones dinámicas del sistema de la
figura 5-2.
En ciertas situaciones la operación de muestreo no se efectúa con respecto al
tiempo y éste no es la variable independiente. En estas circunstancias las ecuaciones
dinámicas discretas describen "eventos" discretos y se expresan como
TÉCNICAS CON VARIABLES DE ESTADO 191
x(k + 1) = 4(1)x(k) + 0(1)u(k) (5-44)
c(k) = Dx(k) + Eu(k) (5-45)
la cual se obtiene de las ecuaciones 5-39 y 5-43 si se fija T = 1 en ellas. Otra manera
de considerar esta notación es que el periodo de muestreo ha sido normalizado con
un valor igual a uno en las expresiones anteriores.
=kT [x[(k + 1)T — x(kT)]] (5-50)
T
[[x(k + 1)T — x(kT)]] = Ax(kT) + Bu(kT) (5-51)
192 SISTEMAS DE CONTROL DIGITAL
Al reacomodar la ecuación 5-51 de modo que tenga la forma de una ecuación de datos
discretas, se tiene
x[(k + 1)T] = (1 + T A)x(kT) + TBu(kT) (5-53)
Las soluciones de las ecuaciones de estado discretas 5-39, 5-44 o 5-46, son las ecua-
ciones de transición de estados del sistema de datos discretos. Dado que todas estas
ecuaciones tienen la misma forma, las soluciones también son muy parecidas, empero
las ecuaciones 5-39 y 5-44 contienen la matriz de transición de estados 4(7) de la
matriz A de la ecuación 5-8, la cual es no singular; por tanto, la transición de estados
para los sistemas descritos por las ecuaciones 5-39 y 5-44 puede llevarse a cabo hacia
adelante o hacia atrás. En general, para ecuaciones de estado discretas dadas por
la expresión 5-46, no hay restricción alguna sobre la matriz A, no obstante, si la
transición de estados se da hacia atrás, entonces A debe ser no singular.
A continuación se derivan las ecuaciones de transición de estados de las ecuacio-
nes de estado 5-39, resultado que después se extiende a las ecuaciones 5-44 y 5-46.
El método más directo para resolver la ecuación 5-39 es la recursión. Las ecuacio-
nes siguientes se obtienen al sustituir de manera progresiva k = 0, 1, 2, ... en la
ecuación 5-39
k=O
x(T) = (f)( T)x(0) + 0(T)u(0) (5-54)
k=1
x(2T) = 4('T)x(T) + 0(T)u(T)
= (1)(T)[4)(T)x(0) + 0(T)u(0)] + 0(T)u(T)
= 4(7')1( T)x(0) + 4(T)0(T)u(0) + O( T)u( T) (5-55)
Ecuación de estados
x(k + 1) = 4)(1)x(k) + 0(1)u(k) (5-44)
194 SISTEMAS DE CONTROL DIGITAL
N-1
x(N) = (f►(N));(0) + E4)(N
1 — i — 1)0(1)u(i) (5-65)
¿ =o
Ecuación de estados
x(k + 1) = Ax(k) + Bu(k) (5-66)
donde
N -->
que tiene la misma forma que la ecuación 5-59. La solución de las ecuaciones de
estado 5-44 y 5-46 se obtiene de manera similar con el método dé la transformada z.
Nótese que tomar la transformada z de AA implica que
Y(Ak) = E
c° AkZ—k (5-79)
k=O
Ejemplo 5-1
Aquí se ilustra el análisis de un sistema de datos discretos de lazo abierto con el método
de variables de estado. La figura 5-3 presenta el diagrama de bloques del sistema en
consideración. Las ecuaciones dinámicas que describen el proceso lineal son
-
d x 1(0
dt 311[xx:Ittl [01]u(0
(5-80)
dx 2(t)
dt
c(t) = x1(t) (5-81)
196 SISTEMAS DE CONTROL DIGITAL
T
Figura 5-3. Sistema digital de lazo abierto.
donde xi(t) y x2(t) son las variables de estado, c(t) la salida escalar y u(t) la entrada
escalar. También como u(t) es la salida de un retenedor de orden cero,
u(t) = u(kT) = r(kT) (5-82)
(5-89)
donde N es cualquier entero positivo.
donde
U1
U(z) = U i
u 2(z) (5-91)
U p(z)
Cq(z)
es la transformada del vector de salida de dimensión q x 1;
G i i (Z) G 12(z) - - - Gip(z)
[ G2 (z) G22(z) • ..z G2P(z)
CY(z ) = ! (5-93)
G. i (z) G. 2(z) • - - Gqp(z)
donde 1 i. qy 1 j p.
Ejemplo 5-2
C(z) — e'
G(z) = Gi(s))
R(z) s
1
= (1 — (5-112)
s(s + 1Xs + 2))
orden está dada por la ecuación matricial de estados de la expresión 5-39. La ecuación
característica del sistema, a menudo también conocida como ecuación característica
de «II está definida como el determinante de zI — (1)(t) igualado con cero; esto es,
zI — 4)(T)1 = 0 (5-113)
Ejemplo 5-3
Algunas de las propiedades más importantes de los valores propios son [la
ecuación de estado de la expresión 5-46 se emplea en el entendido de que A se
reemplaza con (1)(7) si la ecuación de estado tiene la forma de la expresión 5-39]:
Ejemplo 5-4
Considérese el sistema de datos discretos del ejemplo 5-2. Los vectores propios
asociados con los valores propios zi = e-2T y z2 = e- -T son p i y p2, respectivamente.
Por definición, los vectores propios satisfacen la siguiente ecuación:
Ez i I — 4)(T)]p i = 0 (5-121)
Es fácil demostrar que ambas son linealmente dependientes,