Sistemas Discretos de Control - Basil M Al
Sistemas Discretos de Control - Basil M Al
Sistemas Discretos de Control - Basil M Al
CONTROL
Un Enfoque Práctico
BASIL M. AL-HADITHI
DOCTOR INGENIERO INDUSTRIAL
SISTEMAS DISCRETOS DE
CONTROL
Un Enfoque Práctico
2007
SISTEMAS DISCRETOS DE CONTROL
Un Enfoque Práctico
ISBN: 978-84-9821-872-5
Depósito legal: M-42663-2007
Impreso en España
Printed by Publidisa
Contenido
Prólogo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . vii
2 Herramientas Matemáticas 5
2.1 La Transformada Z . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.1.1 Secuencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.1.2 Ecuaciones en Diferencias . . . . . . . . . . . . . . . . . . . 6
2.2 Propiedades y Teoremas . . . . . . . . . . . . . . . . . . . . . . . . 7
2.2.1 Adición . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2.2.2 Multiplicación por una Constante . . . . . . . . . . . . . . . 8
2.2.3 Linealidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2.2.4 Multiplicación por ak . . . . . . . . . . . . . . . . . . . . . 8
2.2.5 Teorema del Desplazamiento . . . . . . . . . . . . . . . . . 9
2.2.6 Teorema de la Traslación Compleja . . . . . . . . . . . . . . 9
2.2.7 Teorema de la Diferenciación Compleja . . . . . . . . . . . 9
2.2.8 Convolución de Secuencias . . . . . . . . . . . . . . . . . . . 9
2.2.9 Teorema del Valor Inicial . . . . . . . . . . . . . . . . . . . 10
2.2.10 Teorema del Valor Final . . . . . . . . . . . . . . . . . . . . 10
2.3 Transformada Z de Funciones Básicas . . . . . . . . . . . . . . . . 10
2.4 La Transformada Z Inversa . . . . . . . . . . . . . . . . . . . . . . 11
2.4.1 Método de la División Directa . . . . . . . . . . . . . . . . 11
2.4.2 Método de la Expansión en Fracciones Parciales . . . . . . 11
2.4.3 Método de la Integral de Inversión . . . . . . . . . . . . . . 14
2.5 Resolución de Ecuaciones en Diferencias . . . . . . . . . . . . . . . 14
2.6 Función de Transferencia Pulso . . . . . . . . . . . . . . . . . . . . 16
iv Sistemas Discretos de Control
3 Sistemas Muestreados 41
3.1 Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3.2 Sistemas Muestreados . . . . . . . . . . . . . . . . . . . . . . . . . 41
3.2.1 Muestreo mediante Impulsos . . . . . . . . . . . . . . . . . 42
3.2.2 Teorema del Muestreo . . . . . . . . . . . . . . . . . . . . . 43
3.2.3 Reconstrucción de Señales . . . . . . . . . . . . . . . . . . 44
3.3 Análisis de Sistemas Discretos . . . . . . . . . . . . . . . . . . . . 47
3.3.1 Función de Transferencia Equivalente . . . . . . . . . . . . 49
3.3.2 Sistemas con Elementos en Serie . . . . . . . . . . . . . . . 50
3.3.3 Sistema con Bloqueador de Orden Cero . . . . . . . . . . . 51
3.3.4 Sistema con Bloqueador de Orden Uno . . . . . . . . . . . 52
3.3.5 Otras Configuraciones de Sistemas Realimentados . . . . . 52
3.3.6 Sistema Realimentado con Controlador Discreto . . . . . . 57
3.4 Resolución mediante la Integral de Convolución . . . . . . . . . . 58
Problemas Resueltos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
Ejercicios con MATLAB . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
Prólogo
El presente libro se enmarca dentro del área del control de sistemas discretos.
Respecto a los contenidos, se ha procurado simplificar al máximo la teorı́a, sin
profundizar en los desarrollos matemáticos, proporcionando al lector involucrado
en esta disciplina las ideas y conceptos necesarios para abordar el análisis y diseño
de los sistemas discretos. El libro contiene el material necesario para ampliar
las técnicas de control de sistemas, aplicadas a los sistemas continuos en cursos
anteriores, para enlazar ahora con aquellos sistemas que van a ser tratados con
técnicas computacionales.
La bibiliografia disponible en el campo del control, tanto de sistemas conti-
nuos como discretos, es abundante. Este libro proporciona al alumno/lector una
serie de problemas acompañados de ejercicios resueltos mediante la aplicación
MATLAB1 . De esta manera, el alumno/lector tiene la posibilidad de afianzar
los conceptos teóricos y verificar por si mismo las soluciones de los ejercicios del
libro y de todos aquellos que le vayan surgiendo. El uso de MATLAB puede acer-
car al alumno que cursa alguna asignatura de control al campo del análisis y del
diseño de los sistemas, de una forma más participativa y, a su vez, tremendamente
atractiva.
Se trata de un libro de control de sistemas discretos, estructurado en seis
grandes capı́tulos, cada uno de los cuales comienza con una introducción que
resume los fundamentos teóricos que posteriormente se aplicarán. Puede consi-
derarse un segundo curso de control, donde se asume que el lector ya es conocedor
del funcionamiento básico de teorı́a de sistemas continuos.
En el Capı́tulo 1 se encontrará una introducción a los sistemas discretos de
control. Se trata principalmente de definir los distintos tipos de señales como la
señal muestreada y la señal digital. También se explican los bloques principales de
un sistema de control digital como el muestreador, el bloqueador, el convertidor
1
r es una marca registrada de The MathWorks Inc.
MATLAB
viii Sistemas Discretos de Control
Basil M. Al-Hadithi
Madrid, Septiembre de 2007
1
Capı́tulo 1
1.1 Introducción
A medida que el computador y los microcontroladores se fueron introduciendo
en los esquemas de control, se hizo necesario el estudio de nuevas técnicas de
análisis y diseño para controlar los sistemas. Aparece ası́ una nueva disciplina
bajo el nombre de sistemas discretos. En la teorı́a de control, el concepto de
señal y sistema continuos se asocian al mundo analógico, mientras que la señal y
sistema discretos lo hacen al mundo digital.
Una de las ventajas principales que presentan los nuevos sistemas discretos
de control frente los sistemas continuos es su flexibilidad, ya que permiten la
implementación de distintas estrategias de control a muy bajo coste. Los contro-
ladores o reguladores analógicos están rı́gidamente unidos a su diseño hardware
haciéndolos inviables a cualquier cambio en su algortimo de control sin que haya
que rediseñarlo completamente. Este inconveniente queda subsanado con el uso
de microcontroladores y programas de ordenador para implementar controladores
o reguladores digitales. El notable avance de la tecnologı́a hace que estos nuevos
dispositivos digitales sean cada vez más rápidos y precisos.
La realidad de nuestro entorno es analógica, los sistemas y procesos de con-
trol manipulan señales continuas, esto es, toman distintos valores de forma con-
tinua a lo largo de tiempo. La inclusión del computador hace imposible el uso
directo de señales continuas (analógicas) si previamente no se convierten a mag-
nitudes digitales (bits). Esto conlleva la incorporación de nuevos dispositivos
(bloqueadores, muestreadores, convertidores...etc) a los esquemas de control y
2 Capı́tulo 1: Introducción a los Sistemas Discretos de Control
• Señal Muestreada: es una señal cuya amplitud tiene valor sólo en ins-
tantes discretos de tiempo, tal y como se muestra en la figura 1.1.
• Señal Digital: suele ser la señal muestreada y codificada para ser inter-
pretada por un computador.
Figura 1.1: Tipos de señales: (a) señal analógica, (b) y (c) señal muestreada
• Convertidor A/D: convierte una señal analógica en una señal digital con
datos codificados (normalmente en binario).
Capı́tulo 2
Herramientas Matemáticas
2.1 La Transformada Z
En los sistemas continuos se emplea la transformada de Laplace para el análisis
y diseño de aquellos, facilitando el tratamiento de los modelos que suelen venir
dados a través de ecuaciones diferenciales.
El hecho de trabajar con señales que toman valores sólo en los instantes de
tiempo obtenidos por el proceso de muestreo, obliga a introducir nuevas herra-
mientas, conceptos y teoremas que faciliten el análisis de los sistemas discretos.
Estas señales se obtienen al muestrear un señal continua x(t), sólo en los ins-
tantes de muestreo, para obtener la señal muestreada x(kT ). La dinámica de
los sistemas discretos se describen mediante ecuaciones en diferencias que son el
equivalente a las ecuaciones diferenciales de los sistemas continuos.
2.1.1 Secuencias
Una secuencia es un conjunto numerado de elementos que se denomina x(k)
siendo “k” el ı́ndice de la secuencia con un rango de variación −∞ ≤ k ≤ ∞.
{x(k)} = {0, 1, 2, 3, 4, }
6 Herramientas Matemáticas
Figura 2.1:
Figura 2.2:
2.2.1 Adición
La transformada Z de la suma de dos secuencias es igual a la suma de la trans-
formada Z de cada una de las secuencias, es decir,
2.2.3 Linealidad
Aplicando las dos propiedades anteriores, se puede obtener la siguiente propiedad
de linealidad como sigue:
dX(z)
T = −z −1 Z [tx(t)]
dz
Z [x(k)] = X(z)
Z [y(k)] = Y (z)
∞
∞
1
X(z) = Z [x(t)] = 1.z −k = 1 + z −1 + z −2 + · · · =
k=0 k=0
1 − z −1
- Rampa unitaria
t ∀t ≥ 0
x(t) =
0 t<0
∞
∞
X(z) = Z [x(t)] = kT.z −k = T z −1 + 2T z −2 + 3T z −3 + · · ·
k=0 k=0
T z −1
= (2.5)
(1 − z −1 )2
2.4 La Transformada Z Inversa 11
- Exponencial
x(t) = e−at ∀t ≥ 0
∞
∞
X(z) = Z [x(t)] = e−akT .z −k = 1 + e−aT z −1 + e−2aT z −2 + · · ·
k=0 k=0
1
= (2.6)
1 − e−aT z −1
bo z m + b1 z m−1 + · · · + bm−1 z + bm
X(z) =
z n + a1 z n−1 + · · · + an−1 z + an
Primero hay que factorizar el polinomio del denominador de X(z) de la forma:
12 Herramientas Matemáticas
bo z m + b1 z m−1 + · · · + bm−1 z + bm
X(z) =
(z − p1 )(z − p2 ) · · · (z − pn )
Luego se expande X(z) en fracciones simples de manera que puedan encon-
trarse sus transformadas fácilmente en la siguiente tabla 2.1.
a1 a2 an
X(z) = + + ··· +
z − p1 z − p2 z − pn
siendo cada coeficiente,
a1 = [(z − p1 )X(z)]z=p1
..
.
an = [(z − pn )X(z)]z=pn
2.4 La Transformada Z Inversa 13
• Polos múltiples:
Si X(z) tiene polos múltiples p1 de orden r, junto a otras raı́ces simples:
br br−1 b1 ar+1
X(z) = + + ··· + + ··· +
(z − p1 )r (z − p1 )r−1 z − p1 z − pr+1
ar+2 an
+ + ··· + (2.7)
z − pr+2 z − pn
br = [(z − p1 )r X(z)]z=p1
d
br−1 = [(z − p1 )r X(z)]
dz z=p1
1 dj
br−j = [(z − p1 )r X(z)]
j! dz j z=p1
1 dr−1
b1 = [(z − p1 )r X(z)]
(r − 1)! dz r−1 z=p1
α1 z + α2 a3 an
X(z) = + + ··· +
(z − p1 )(z − p2 ) z − p3 z − pn
1
Z −1 [X(z)] = x(kT ) = x(k) = X(z)z k−1 dz
2πj C
siendo C un cı́rculo en el plano Z, con centro en el origen, dentro del cual
quedan todos los polos de X(z)z k−1 .
La solución de la integral de inversión aplicando el teorema de los residuos es:
m
Z −1 [X(z)] = x(k) = residuo de X(z)z k−1 en el polo z = zi de X(z)z k−1
i=1
= K1 + K2 + · · · + Km (2.8)
U(z) Y(z)
G(z)
Figura 2.3:
Y (z) b0 + b1 z −1 + · · · + bm z −m
=
U (z) 1 + a1 z −1 + a2 z −2 + · · · + an z −n
La función de transferencia G(z) coincide con la respuesta del sistema ante
entrada impulso δk .
u(k) = δk
Z [u(k)] = Z [δk ]
Figura 2.4:
Problemas Resueltos
Problemas Resueltos 19
Problema 2.1
xk = x0 δk + x1 δk−1 + x2 δk−2
luego,
Problema 2.2
∞
y0 = x0−n gn = x0 g0 + x−1 g1 + x−2 g2 = 0 · 1 + 0 · 2 + 0 · 1 = 0
n=0
∞
y1 = x1−n gn = x1 g0 + x0 g1 + x−1 g2 = 3 · 1 + 0 · 2 + 0 · 1 = 3
n=0
∞
y2 = x2−n gn = x2 g0 + x1 g1 + x0 g2 = 4 · 1 + 3 · 2 + 0 · 1 = 10
n=0
∞
y3 = x3−n gn = x3 g0 + x2 g1 + x1 g2 + x0 g3 = 0 · 1 + 4 · 2 + 3 · 1 + 0 · 0 = 11
n=0
∞
y4 = x4−n gn = x4 g0 +x3 g1 +x2 g2 +x1 g3 +x0 g4 = 0·1+0·2+4·1+3·0+0·0 = 4
n=0
20 Problemas Resueltos
∞
y5 = x5−n gn = x5 g0 +x4 g1 +x3 g2 +x2 g3 +x1 g4 = 0·1+0·2+0·1+4·0+3·0 = 0
n=0
G(z) = Z [gk ] = 1z 0 + 2z −1 + 1z −2 = (1 + 2z −1 + z −2 )
Problema 2.3
1. Escalón unitario
0, ∀t < 0
x(t) =
1, ∀t ≥ 0
∞
1 z
X(z) = Z [x(t)] = 1 z −k = 1 + z −1 + z −2 + . . . = =
k=0
1 − z −1 z−1
2. Rampa unitaria
Problemas Resueltos 21
0, ∀t < 0
x(t) =
t, ∀t ≥ 0
∞
∞
−k
X(z) = Z [x(t)] = kT z =T k z −k = T z −1 + 2T z −2 + 3T z −3 + . . . =
k=0 k=0
T z −1 Tz
= −1 2
=
(1 − z ) (z − 1)2
3. Exponencial
x(t) = e−at ∀t ≥ 0
∞
X(z) = Z [x(t)] = e−akT z −k = 1 + e−aT z −1 + e−2aT z −2 + e−3aT z −3 + . . . =
k=0
1 z
= =
1− e−aT z −1 z − e−aT
Problema 2.4
∞
X(z) = Z[x(t)] = ak kT z −k = T (az −1 + 2a2 z −2 + 3a3 z −3 + . . . =
k=0
az −1 az
=T =T
(1 − az −1 )2 (z − a)2
Tz
Z[t] = Z[kT ] =
(z − 1)2
y sustituyendo z por (a−1 z) queda finalmente,
a−1 z az
X(z) = T −1 2
=T
(a z − 1) (z − a)2
dX(z)
−z = Z[kx(k)]
dz
En la expresión a transformar ak t = ak kT , se define
x(k) = ak T
de manera que la transformada a calcular se convierte en,
Z[ak t] = Z[kx(k)]
De la tabla de transformadas se conoce,
z
X(z) = Z[x(k)] = Z[ak T ] = T
z−a
Calculando la derivada,
dX(z) (z − a) − z −a
=T 2
=T
dz (z − a) (z − a)2
Finalmente, por aplicación de la propiedad de la derivada se obtiene,
dX(z) −a az
Z[ak t] = Z[kx(k)] = −z = −zT =T
dz (z − a)2 (z − a)2
Problema 2.5
∞
X(z) = Z[x(t)] = (kT )2 z −k = 0 + T 2 z −1 + (2T )2 z −2 + (3T )2 z −3 + . . . =
k=0
Problemas Resueltos 23
z −1 (1 + z −1 ) z(z + 1)
= T2 −1 3
= T2
(1 − z ) (z − 1)3
dX(z)
Z[tx(t)] = −T z
dz
siendo x(t) = t y su transformada X(z) = Z[x(t)] = T (z−1)
z
2
−T z − T T 2z2 + T 2z T 2 z(1 + z) T 2 z −1 (z −1 + 1)
Z[t2 ] = −T z = = =
(z − 1)3 (z − 1)3 (z − 1)3 (1 − z −1 )3
Problema 2.6
0.5z 2
X1 (z) =
(z + 0.5)(0.5z − 0.5)
10
X2 (z) =
z2 − 1.5z + 0.5
0.5z 2 z2 1
X1 (z) = = =
(z + 0.5)0.5(z − 1) (z + 0.5)(z − 1) (1 + 0.5z −1 )(1 − z −1 )
10 10z −2
X2 (z) = =
(z − 1)(z − 0.5) (1 − z −1 )(1 − 0.5z −1 )
1
x1 (∞) = lim (1 − z −1 )X1 (z) = lim (1 − z −1 ) = 0.6667
z→1 z→1 (1 + 0.5z −1 )(1 − z −1 )
10z −2
x2 (∞) = lim (1 − z −1 )X2 (z) = lim (1 − z −1 ) = 20
z→1 z→1 (1 − z −1 )(1 − 0.5z −1 )
Problema 2.7
z2 A B
= +
(z − 2)(z − 1) z−2 z−1
z2
A = (z − 2) =4
(z − 2)(z − 1) z=2
z2
B = (z − 1) = −1
(z − 2)(z − 1) z=1
Ası́ pues,
4 −1 4z −1 z −1
Y [z] = + = −1
−
z−2 z−1 1 − 2z 1 − z −1
siendo su antitransformada,
{yk } = {1, 3, 7, 15 . . .}
Problemas Resueltos 25
z2 z −2 1
Y (z) = · −2 =
(z − 2)(z − 1) z (1 − 2z −1 )(1 − z −1 )
Descomponiendo en fracciones simples,
1 A B
= +
(1 − 2z −1 )(1 − z −1 ) 1 − 2z −1 1 − z −1
−1 1
A = (1 − 2z ) =2
(1 − 2z )(1 − z −1 )
−1
z=2
1
B = (1 − z −1 ) = −1
(1 − 2z −1 )(1 − z −1 ) z=1
Ası́ pues,
2 −1
Y (z) = −1
+
1 − 2z 1 − z −1
siendo su antitransformada,
yk = 2 · 2k − 1k = 2(k+1) − 1
Por último, dando valores a k se obtiene:
{yk } = {1, 3, 7, 15 . . .}
z2 1 1
Y (z) = = −1 −1
= −1
(z − 2)(z − 1) (1 − 2z )(1 − z ) 1 − 3z + 2z−2
1 1 − 3z −1 + 2z −2
−1
−1 + 3z − 2z −2 1 + 3z −1 + 7z −2 + 15z −3 + . . .
3z −1 − 2z −2
−3z −1 + 9z −2 − 6z −3
7z −2 − 6z −3
−7z −2 + 21z −3 − 14z −4
15z −3 − 14z −4
..
.
{yk } = {1, 3, 7, 15 . . .}
Y (z) 1
=
X(z) 1 − 3z −1 + 2z −2
Y (z)(1 − 3z −1 + 2z −2 ) = X(z)
Calculando la transformada inversa se obtiene,
yk − 3yk−1 + 2yk−2 = xk
Como xk es la secuencia impulso 1 0 0 . . . , se pueden obtener los valores de
yk :
Problema 2.8
0.5 − e−aT
X(z) =
(z − 0.5)(z − e−aT )
0.5 − e−aT A B
−aT
= +
(z − 0.5)(z − e ) z − 0.5 z − e−aT
0.5 − e−aT 0.5 − e−aT
A = (z − 0.5) = =1
(z − 0.5)(z − e−aT ) z=0.5
0.5 − e−aT
−aT 0.5 − e−aT 0.5 − e−aT
B = (z − e ) = = −1
(z − 0.5)(z − e−aT ) z=e−aT
e−aT − 0.5
Ası́ pues,
1 1 z −1 z −1
X(z) = − = −
z − 0.5 z − e−aT 1 − 0.5z −1 1 − z −1 e−aT
siendo su antitransformada,
Problema 2.9
1.5z −2 1.5z −1
F (z) = −1 −2
= z −1 = F1 (z)z −1
1 − 0.2z + 0.1z 1 − 0.2z −1 + 0.1z −2
28 Problemas Resueltos
e−aT sen(wT )z −1
X(z) = → x(t) = e−at sen(wt)
1 − 2e−aT cos(wT )z −1 + e−2aT z −2
1.5 0.3z −1
· · z −1 = F (z)
0.3 1 − 0.2z −1 + 0.1z −2
con antitransformada:
1.5 −1.15(t−1)
f (t) = e sen(1.25(t − 1))
0.3
Problema 2.10
z −1
G(z) =
(1 + 2z −1 )n
cuando n toma los valores de 1, 2, 3.
z −1 1 + 2z −1
−z −1 − 2z −2 z −1 − 2z −2 + 4z −3 − 8z −4 + . . .
−2z −2
2z −2 + 4z −3
..
.
Problemas Resueltos 29
gk = k(−2)k−1
Y (z) z −1
=
X(z) 1 + 4z −1 + 4z −2
Y (z)(1 + 4z −1 + 4z −2 ) = z −1 X(z)
Calculando la transformada inversa se obtiene,
z −1 z3 z2
G(z) = · =
(1 + 2z −1 )3 z 3 (z + 2)3
z2
gk = Z −1 [G(z)] = residuos z k−1 = [k]z=−2
(z + 2)3 triple
1 d2 3 z2 k−1 1 d2 k+1
[k]z=−2 = lim (z + 2) z = lim z =
(3 − 1)! z→−2 dz 2 (z + 2)3 2 z→−2 dz 2
1 k(k + 1)
= lim k(k + 1)z k−1 = (−2)k−1 = gk
2 z→−2 2
Dando valores a k se obtiene la secuencia de ponderación pedida:
Problema 2.11
T z −1
W (z) =
(1 − z −1 )2
La señal que resulta de convolucionar ambas:
T z −1
Y (z) = X(z)W (z) =
(1 − z −1 )3
Realizando la división directa,
Problemas Resueltos 31
T z −1 1 − 3z −1 + 3z −2 − z −3
−T z −1 +3T z −2 −3T z −3 +T z −4 T z −1 +3T z −2 +6T z −3 +10T z −4 . . .
−2 −3
3T z − 3T z + T z −4
−3T z −2 + 9T z −3 − 9T z −4 + 3T z −5
6T z −3 − 8T z −4 + 3T z −5
−6T z −3 + 18T z −4 − 18T z −5 + 6T z −6
..
.
T z −1 T z2
Y (z) = =
(1 − z −1 )3 (z − 1)3
T z2
yk = Z −1 [Y (z)] = residuos z k−1 = [k]z=1
(z − 1)3 triple
1 d2 3 Tz
2
[k]z=1 = lim (z − 1) z k−1 =
(3 − 1)! z→1 dz 2 (z − 1)3
1 d2 k+1 1
= lim Tz = T k(k + 1)
2 z→1 dz 2 2
∞
{yk } = {xk } ∗ {wk } = xn wk−n
n=−∞
y0 = x0 w0 + x1 w−1 = 0
y1 = x0 w1 + x1 w0 + x2 w−1 = T + 0 = T
y2 = x0 w2 + x1 w1 + x2 w0 + x3 w−1 = 2T + T + 0 = T = 3T
y3 = x0 w3 + x1 w2 + x2 w1 + x3 w0 + x4 w−1 = 3T + 2T + T + 0 = T = 6T
..
.
Obteniendo la misma secuencia:
Problema 2.12
Dibujar la secuencia de salida yk+2 del sistema G(z) dado, cuando se le somete
a una entrada escalón unitario.
1 + z −1 + 0.5z −2
G(z) =
1 + 0.5z −1 + z −2 + 2z −3
2.5
1.5
0.5
−0.5
−1
−1.5
0 1 2 3 4 5
−2 −1 0 1 2 3
34 Problemas Resueltos
Amplitud
Indice
6
Amplitud
0
0 2 4 6 8 10
Indice
Secuencia exponencial
1
0.9
0.8
0.7
0.6
Amplitud
0.5
0.4
0.3
0.2
0.1
0
0 2 4 6 8 10
Indice
Convolucion x con g
12
10
8
Amplitud
0
0 2 4 6 8
Indice
Ejercicio 2.5 Obtener la transformada Z del escalón unitario como función simbólica:
z2
Y (z) =
(z − 2)(z − 1)
0.5 − e−aT
X(z) =
(z − 0.5)(z − e−aT )
>> syms a T z;
>> X = sym((0.5-exp(-a*T))/(z-0.5)/(z-exp(-a*T)));
>> x = iztrans(X)
x =
2*(1/2)∧ n-(1/exp(a*T))∧ n*exp(a*T)
>> pretty(x)
n / 1 \ n
2 (1/2) - | ----------- | exp (a T)
\ exp (a T) /
z −1
G(z) =
(1 + 2z −1 )n
cuando n toma los valores de 1, 2, 3.
>> G = tf([1 1 0.5 0],[1 0.5 1 2],-1); % ’G’ como objeto ’tf’
>> step(G); % respuesta al escalón
>> xlabel(’Indice’) % insertar texto en eje X
>> ylabel(’Amplitud’) % insertar texto en eje Y
>> title(’Respuesta al escalón’) % insertar título
Respuesta al escalón
6
2
Amplitud
−2
−4
−6
−8
0 2 4 6 8
Indice
41
Capı́tulo 3
Sistemas Muestreados
3.1 Introducción
El análisis de los sistemas discretos de control implica el uso de la transformada
Z. Como se ha mencionado anteriormente, que gracias a la transformada Z se
podrán aplicar los conocimientos de análisis y diseño de los sistemas continuos a
los sistemas discretos.
Las señales en tiempo discreto se obtienen mediante el proceso de muestreo
(conversión A/D) que después de procesarlas en el computador se convierten a
señales en tiempo continuo mediante el proceso de reconstrucción (conversión
D/A).
El proceso de muestreo es necesario debido a que el computador realiza los
cálculos en determinados instantes de tiempo.
Es importante especificar que el periodo de muestreo T va a ser uniforme y
constante, esto es, cada dato muestreado se obtiene a intervalos iguales de tiempo.
Después del muestreo y la reconstrucción va a haber pérdida de información
debido a que primero, el muestreador obvia el valor de la señal continua entre ins-
tantes consecutivos y segundo, el bloqueador, dependiendo de su tipo, aproxima
con mayor o menor exactitud la señal muestreada.
x(t) x*(t)
d (t)
T
t t
Modulador
∞
x∗ (t) = x(t)δT (t − kT )
k=0
= x(0)δ(t) + x(T )δ(t − T ) + x(2T )δ(t − 2T ) + · · · (3.1)
3.2 Sistemas Muestreados 43
z = eT s
o
1
s = ln z (3.3)
T
Por lo tanto la ecuación (3.2) se convierte en:
∞
∗ ∗
L [x (t)]s= 1 lnz = X (s) |s= 1 ln z = x(kT )z −k = X(z) = Z [x(kT )] = Z [x(t)]
T T
k=0
y también:
fs > 2f0
siendo fs la frecuencia de muestreo y f0 la componente de frecuencia más
elevada de la señal x(t).
44 Capı́tulo 3: Sistemas Muestreados
ωs > 2ω0
o bien,
1 1
>2
Ts T0
T0
Ts <
2
Muestreador
d (t)
T
t t t
Modulador Bloqueador
orden cero
Muestreador
a- Modelo real
b- Modelo matemático
h(t) = x(kT )
1 − e−T s
B0 (s) =
s
• Bloqueador Orden Uno
Aproxima la salida mediante un polinomio de orden 1. La reconstrucción
se realiza considerando el valor en el instante (k-1) y proyectando desde
3.3 Análisis de Sistemas Discretos 47
el instante k hasta (k+1) una pendiente igual a la que une los valores de
los instantes [(k-1), k] como se ve en la figura 3.8. Matemáticamente esta
interpolación del bloqueador se expresa como sigue:
1 − e−T s 2 T s + 1
B1 (s) = ( )
s T
x(t) x*(t)
h(t)
T (t)
t t t
Modulador Bloqueador
orden uno
• Bloqueador Poligonal
Es también de orden 1 pero más preciso y exacto. genera una pendiente
igual a la que une, la muestra actual y la anterior pero proyectándola desde
la muestra actual en el punto de amplitud igual a la de la muestra anterior.
Los bloqueadores de orden superior a 1 van extrapolando con funciones que
incorporan un mayor número de muestras anteriores, de manera que van
consiguiendo una reconstrucción cada vez mejor. El problema es que todo
este proceso introduce retardos en la actuación final, pudiendo inestabilizar
los sistemas.
• Análisis del sistema continuo equivalente que es propio del análisis de los
sistemas de control continuo como se muestra en la figura 3.9.
Muestreador
x(t) x(kT) Sistema y(kT) Bloqueador y(t)
discreto orden cero
T
Y (z)
= G(z) (3.6)
X(z)
x(t) y(t)
G(s)
X(s) Y(s)
Figura 3.12:
Se cumple que
50 Capı́tulo 3: Sistemas Muestreados
Y ∗ (s) = [G(s)X(s)]∗
Y ∗ (s) = [GX(s)]∗ (3.8)
En términos de la transformada Z:
Y (z)
= G2 (z)G1 (z) (3.11)
X(z)
3.3 Análisis de Sistemas Discretos 51
En términos de la transformada Z:
Y (z) = G1 G2 (z)X(z)
Se debe mencionar que G1 G2 (z) = G1 (z)G2 (z)
1 − e−T s
B0 (s)G(s) = G(s)
s
desde la ecuación (3.2) se sabe que z = eT s , entonces e−T s representa un
retardo en el dominio continuo equivalente a un z −1 en dominio Z.
Por lo tanto,
1 − e−T s G(s)
Z [B0 (s)G(s)] = Z G(s) = (1 − z −1 )Z
s s
52 Capı́tulo 3: Sistemas Muestreados
1 − e−T s 2 T s + 1
B1 (s)G(s) = ( ) G(s)
s T
1 − e−T s 2 T s + 1 Ts + 1
Z [B1 (s)G(s)] = Z ( ) G(s) = (1 − z −1 )2 Z G(s)
s T T s2
H(s)
por tanto
Como
G∗ (s)X ∗ (s)
Y ∗ (s) =
1 + GH ∗ (s)
En términos de la notación de la transformada Z, se obtiene:
G(z)X(z)
Y (z) =
1 + GH(z)
La función de transferencia pulso es:
Y (z) G(z)
=
X(z) 1 + GH(z)
Ejemplo 3.6 Calcular la función de transferencia del lazo cerrado del sistema
de la figura 3.18
H(s)
Figura 3.18:
Y ∗ (s)
= X ∗ (s) − H ∗ (s)Y ∗ (s)
G∗ (s)
G∗ (s)X ∗ (s)
Y ∗ (s) =
1 + G∗ (s)H ∗ (s)
En términos de la notación de la transformada Z, se obtiene:
G(z)X(z)
Y (z) =
1 + G(z)H(z)
La función de transferencia pulso es:
Y (z) G(z)
=
X(z) 1 + G(z)H(z)
Ejemplo 3.7 Calcular la función de transferencia del lazo cerrado del sistema
de la figura 3.19
Y*(s)
H(s)
d (t)
T
Figura 3.19:
Y (s)
= X(s) − H(s)Y ∗ (s)
G(s)
GX(z)
Y (z) =
1 + GH(z)
Ejemplo 3.8 Calcular la función de transferencia del lazo cerrado del sistema
de la figura 3.20
M*(s) M(s)
H 2 (s) H1 (s)
d (t)
T
Figura 3.20:
Sustituyendo en (3.19)
Y ∗ (s) X ∗ (s)
=
G∗ (s) 1 + GH1∗ (s).H2∗ (s)
G∗ (s)X ∗ (s)
Y ∗ (s) =
1 + GH1∗ (s).H2∗ (s)
(3.23)
G(z)X(z)
Y (z) =
1 + GH1 (z).H2 (z)
La función de transferencia pulso es:
Y (z) G(z)
=
X(z) 1 + GH1 (z).H2 (z)
H(s)
U ∗ (s) = RD
∗
(s)E ∗ (s) (3.26)
por tanto
U ∗ (s) = RD
∗
(s)E ∗ (s)
entonces
U ∗ (s) = RD∗
(s)X ∗ (s) − RD ∗
(s)BGH ∗ (s)U ∗ (s)
RD∗ (s)X ∗ (s)
U ∗ (s) = ∗ (s)BGH ∗ (s) (3.29)
1 + RD
BG(z)RD (z)X(z)
Y (z) =
1 + RD (z)BGH(z)
La función de transferencia pulso es:
X(s)z
X(z) = residuos de en polos de X(s)
z − eTs
Al igual que se estudió en el apartado (2.4) de la transformada Z inversa, los
residuos se obtienen de la siguiente manera:
Problemas Resueltos
Problemas Resueltos 61
Problema 3.1
(s + 1) A B
= +
s(s + 2) s s+2
(s + 1) 1
A= s =
s(s + 2) s=0 2
(s + 1) 1
B = (s + 2) =
s(s + 2) s=−2 2
Ası́ pues,
0.5 0.5
+ G(s) =
s s+2
El primer sumando representa la transformada de Laplace del escalón unitario;
el segundo es la transformada de la exponencial (e−at ). Consultando directamente
en la tablas la tranformada Z de estos términos se obtiene:
(s + 1) z
G(z) = Z[G(s)] = residuos = [k]s=0 + [k]s=−2
s(s + 2) (z − eT s )
(s + 1) z 1 z 0.5
[k]s=0 = lim s = =
s→0 s(s + 2) (z − eT s ) 2 (z − 1) 1 − z −1
(s + 1) z 1 z 0.5
[k]s=−2 = lim (s + 2) = −2T
=
s→−2 s(s + 2) (z − e )
T s 2 (z − e ) 1 − e−2T z −1
La suma de ambos residuos da la solución final:
62 Problemas Resueltos
Problema 3.2
(s + 1)
G(s) =
s2 (s + 4)
(s + 1) z
G(z) = Z[G(s)] = residuos = [k]s=0 + [k]s=−4
s2 (s + 4) (z − eT s ) doble
1 d (s + 1) z d zs + z
[k]s=0 = lim s2 2 = lim =
(2 − 1)! s→0 ds s (s + 4) (z − e )
T s s→0 ds zs + 4z − seT s − 4eT s
(s + 1) z (s + 1) z 3 z
[k]s=−4 = lim (s + 4) 2 = lim 2
=−
s→−4 s (s + 4) (z − e )
T s s→−4 s (z − e )
T s 16 (z − e−4T )
Problema 3.3
G(s)
x(t) {xk} 1 y(t) {yk}
B0(s)
T s(s+a) T
(1 − e−T s ) 1 1
B0 (s)G(s) = = (1 − e−T s ) 2
s s(s + a) s (s + a)
−1 1
Z[B0 G(s)] = B0 G(z) = (1 − z )Z 2
s (s + a)
Método 1. Por cálculo de los residuos
1 1 z
Z 2 = residuos 2 = [k]s=0 + [k]s=−a
s (s + a) s (s + a) (z − eT s ) doble
1 d 1 z
[k]s=0 = lim s2 =
(2 − 1)! s→0 ds s2 (s + a) (z − eT s )
d z −z(z − 1 − aT )
= lim =
s→0 ds (s + a)(z − eT s ) a2 (z − 1)2
1 z z
[k]s=−a = lim (s + a) 2
= 2
s→−a s (s + a) (z − e )
T s a (z − e−aT )
Sumando ambos residuos queda,
1 z 2 (aT + e−aT − 1) + z(1 − e−aT − aT e−aT )
Z =
s2 (s + a) a2 (z − 1)2 (z − e−aT )
Multiplicando numerador y denominador por z −3 ,
1 z −1 (aT + e−aT − 1) + z −2 (1 − e−aT − aT e−aT )
Z 2 =
s (s + a) a2 (1 − z −1 )2 (1 − z −1 e−aT )
−1 z −1 (aT + e−aT − 1) + z −2 (1 − e−aT − aT e−aT )
B0 G(z) = (1 − z ) =
a2 (1 − z −1 )2 (1 − z −1 e−aT )
64 Problemas Resueltos
1 T z −1 1 1 1 1
B0 G(z) = (1 − z −1 ) −1 2
− 2 −1
+ 2 −aT
=
a (1 − z ) a (1 − z ) a (1 − e z −1 )
Problema 3.4
G(s)
x(t) {xk} (s+1) y(t) {yk}
B0(s) 2
T (s+4) T
(1 − e−T s ) (s + 1) (s + 1)
B0 (s)G(s) = = (1 − e−T s )
s (s + 4)2 s(s + 4)2
(s + 1)
Z[B0 G(s)] = B0 G(z) = (1 − z −1 )Z
s(s + 4)2
Resolviendo a través del cálculo de los residuos,
(s + 1) (s + 1) z
Z 2
= residuos 2
= [k]s=0 + [k]s=−4
s(s + 4) s(s + 4) (z − eT s ) doble
(s + 1) z 1 z 1 1
[k]s=0 = lim s 2
= =
s→0 s(s + 4) (z − e )
T s 16 (z − 1) 16 (1 − z −1 )
1 d (s + 1) z −z 2 + z(12T e−4T + e−4T )
[k]s=−4 = lim (s + 4)2 = =
(2 − 1)! s→−4 ds s(s + 4)2 (z − eT s ) 16(z − e−4T )2
−z 2 + 1.4716z 1.4716z −1 − 1
= =
16(z − 0.3679)2 16(1 − 0.3679z −1 )2
Finalmente, el sistema discreto equivalente buscado es:
−1 −1 1 1 1.4716z −1 − 1
B0 G(z) = (1−z )([k]s=0 +[k]s=−4 ) = (1−z ) + =
16 (1 − z −1 ) 16(1 − 0.3679z −1 )2
−1.3362z −2 + 1.7358z −1
=
16(1 − 0.3679z −1 )2
66 Problemas Resueltos
Problema 3.5
G(s)
Y (z)
1. La función de transferencia M (z) = X(z) .
−1 0.5
B0 G(z) = (1 − z )Z
s(s + 1)
Descomponiendo en fracciones simples,
0.5 A B
= +
s(s + 1) s (s + 1)
0.5
A= s = 0.5
s(s + 1) s=0
0.5
B = (s + 1) = −0.5
s(s + 1) s=−1
Consultando la tabla de transformadas Z,
0.5 0.5 0.5
Z = −
s(s + 1) (1 − z ) (1 − e−T z −1 )
−1
Problemas Resueltos 67
luego,
−1 0.5 0.5 0.25z −1
B0 G(z) = (1 − z ) − =
(1 − z −1 ) (1 − e−T z −1 ) T =0.7 1 − 0.5z −1
0.25xk−1 = yk − 0.25yk−1
yk = 0.25xk−1 + 0.25yk−1
Problema 3.6
G1(s) G2(s)
−1 2z −1 2 2
B0 Gdirecta (z) = (1 − z ) − + =
(1 − z −1 )2 (1 − z −1 ) (1 − 0.37z −1 )
Problema 3.7
1. H(s) = 1
2. H(s) = 1/s
3. H(s) = 2s/(s + 1)
G(s)
H(s)
1. Para H(s) = 1
−1 1 A B
B0 G(z) = (1 − z )Z = (1 − z −1 )Z + =
s(s + 2) s (s + 2)
−1 0.5 0.5 0.43z −1 0.43
= (1 − z ) − = =
(1 − z −1 ) (1 − 0.14z −1 ) (1 − 0.14z −1 ) (z − 0.14)
Por lo tanto,
B0 G(z) 0.43
M1 (z) = =
1 + B0 G(z) z + 0.29
2. Para H(s) = 1/s
−1 1 0.43
B0 G(z) = (1 − z )Z =
s(s + 2) (z − 0.14)
Problemas Resueltos 71
−1 1 A B C
B0 GH(z) = (1 − z )Z 2 = (1 − z −1 )Z 2 + + =
s (s + 2) s s (s + 2)
−1 0.5T z −1 0.25 0.25
= (1 − z ) −1 2
− + =
(1 − z ) (1 − z ) (1 − e−2T z −1 )
−1
T =1
−1 1 2s 2
B0 GH(z) = (1 − z )Z = (1 − z −1 )Z =
s(s + 2) (s + 1) (s + 2)(s + 1)
A B −2 2
= (1 − z −1 )Z + = (1 − z −1 ) + =
(s + 2) (s + 1) (1 − 0.14z ) (1 − 0.37z −1 )
−1
Problema 3.8
G(s)
-
H(s)
10
B0(s)
s
−1 1 0.25 0.25
B0 G(z) = (1 − z )Z = (1 − z −1 )Z − =
s(s + 4) s (s + 4)
−1 0.25 0.25 0.25z −1
= (1 − z ) − =
(1 − z −1 ) (1 − 0.02z −1 ) 1 − 0.02z −1
Problema 3.9
74 Problemas Resueltos
Problemas Resueltos 75
Problema 3.10
- -
1
s
G(s)
4
B0(s)
s T = 0.7
H(s)
Gdirecta (z) 1 − 2z −1 + z −2
M (z) = = =
1 + Gdirecta (z)B0 H(z) 1 + 1.3z −1 − 2.3z −2
z 2 − 2z + 1 z−1
M (z) = 2
=
z + 1.3z − 2.3 z + 2.3
Problema 3.11
G(s)
0.5
H(s)
25
B0 G(z) = (1 − z −1 )Z 2
s (s + 5)
Descomponiendo en fracciones simples,
25 A B C
2
= 2+ +
s (s + 5) s s (s + 5)
25 2
A= s 2 =5
s (s + 5) s=0
1 d 2 25
B= s 2 = −1
(2 − 1)! ds s (s + 5) s=0
25
C = (s + 5) =1
s2 (s + 5) s=−5
Consultando la tabla de transformadas Z,
25 5T z −1 1 1
Z 2 = − +
s (s + 5) (1 − z −1 )2 (1 − z −1 ) (1 − e−5T z −1 ) T =0.4
−1 2z −1 1 1
B0 G(z) = (1 − z ) −1 2
− + =
(1 − z ) (1 − z ) (1 − 0.135z −1 )
−1
Problema 3.12
1
B0(s)
s T=1
H(s)
−1 1 0.5 0.5
B0 G(z) = (1 − z )Z = (1 − z −1 )Z − =
s(s + 2) s (s + 2)
−1 0.5 0.5 0.43z −1
= (1 − z ) −1
− =
1−z 1 − e−2T z −1 T =1 1 − 0.135z −1
Y en el lazo de realimentación,
−1 1 T z −1 z −1
B0 H(z) = (1 − z )Z 2 = (1 − z −1 ) =
s (1 − z −1 )2 1 − z −1
La función de transferencia equivalente será:
G(s)
x(t) {xk} (s+1) y(t) {yk}
B0(s) 2
T (s+4) T
G(s)
-
Ejercicios con MATLAB 79
H(s)
Sampling time: 1
80 Problemas Resueltos
-
H(s)
10
B0(s)
s
Sampling time: 1
81
Capı́tulo 4
4.1 Introducción
Un aspecto fundamental en el análisis de los sistemas continuos y discretos de
control es el estudio de la estabilidad de los mismos. Pero además de exigir a los
sistemas una estabilidad absoluta, se pretende que tengan unas caracterı́sticas
dinámicas y estáticas que se reflejan en el régimen transitorio y permanente de
la respuesta.
En los sistemas continuos de control, la estabilidad se analizaba en el plano
complejo s, de manera que si los polos del sistema en lazo cerrado caı́an dentro del
semiplano negativo, el sistema era estable. Cuando los polos tomaban posiciones
sobre el eje real, el comportamiento del sistema era sobreamortiguado; si los polos
resultaban ser complejos conjugados, el sistema pasaba a ser subamortiguado,
siendo cada vez más oscilatorio cuanto más próximos estaban los polos del eje
imaginario jω. Cuando los polos se situaban sobre este el mismo eje imaginario
se tenı́an sistemas oscilatorios o marginalmente estables.
En el caso de los sistemas discretos, la estabilidad de los mismos se resolverá
también por la localización de los polos en el plano z. Para abordar este estudio
habrá que buscar la correspondencia de puntos del plano s con puntos del plano
z.
82 Capı́tulo 4: Análisis de los Sistemas Discretos de Control
z = eT s (4.1)
Siendo s la variable compleja y que se expresa como:
s = σ + jω (4.2)
Sustituyendo el valor de s de 4.2 en 4.1
z = eT s = eT σ+jω = eT σ · ejT ω
En el semiplano izquierdo del plano s el valor de σ es negativo, haciendo que
el módulo de z sea siempre una magnitud menor que la unidad, de forma que se
cumple |z| = eT σ < 1. El término ejT ω es un valor que decribe una circunferencia
cuando σ varı́a entre −∞ y ∞.
z = eT σ · ejT ω = 1ejT ω
obteniendo una circunferencia de radio unidad tal y como se ve en la figura
4.1
z = e−T σ1 · ejT ω
se obtiene una circunferencia interior a la de radio unidad tal y como se
muestra en la figura 4.2 y 4.3. De esta manera, todos los puntos del semiplano
negativo en s se transformarán en puntos interiores al cı́rculo de radio unidad.
En el caso de σ = σ1 se obtendrá una circunferencia de radio eT σ1 que será
exterior al cı́rculo de radio unidad. Se puede observar que al alejarse del eje
imaginario en el plano s, las cirunferencias en el plano z se aproximan al origen.
4.2 Correspondencia entre el Plano s y el Plano z 83
Figura 4.1:
Figura 4.2:
Im
j
Plano s Plano z
r=1
Re
Figura 4.3:
Im
j
Plano s Plano z
j
T T T
z=e e 1
= e
1
1
Re
Figura 4.4:
Ejemplo 4.1 El área representada en la figura 4.5 está formado por las siguientes
partes:
4.3 Análisis de la Estabilidad 85
−π
I- eje imaginario. σ = 0 y T ≤ω≤ π
T (s = 0 + jω)
I- z = ejωT , vector giratorio desde e−jπ hasta ejπ describiendo una circunfer-
encia en sentido antihorario.
II- z = eσT · ejπ = −eσT , recta a lo largo del eje real desde -1 hasta −e−cT .
IV- z = eσT · e−jπ = −eσT , recta a lo largo del eje real desde −e−cT hasta -1.
jw Im
Plano s Plano z
I
II P
T
III II III
s
I IV Re
-
P
IV T
Figura 4.5:
H(s)
Figura 4.6:
Y (z) B0 G(z)
= (4.3)
R(z) 1 + B0 GH(z)
Se denomina ecuación caracterı́stica a:
Una vez obtenida la tabla y calculados los coeficientes que en ella aparecen, el
criterio de estabilidad de Jury dice que un sistema con la ecuación caracterı́stica
P(z), donde a0 > 0, será estable si se cumplen las condiciones siguientes:
1. |an | < a0
2. p(1) > 0
> 0 si n par
3. p(−1)
< 0 si n impar
b
b
b
b2
0 b4
0 b6
c1 = − b11
0
, c2 = − b11
, c3 = − b11
b1 b3
b1 b5
b1 b7
b
1 b3
1
b1 b5
d1 = − c11
, d2 = − c 1
c1 c2
c1 c3
Y (z) b
=
X(z) 1 − az −1
y(k) Im
Plano z
b
0<a<1 Re
Figura 4.7:
y(k)
Im
Plano z
a>1 Re
b
Figura 4.8:
y(k)
Im
Plano z
k -1<a<0 Re
-b
Figura 4.9:
bz z
y(0) = lim Y (z) = lim . =b
z→∞ z→∞ z − a z − 1
y(k)
Im
Plano z
a<-1 Re
b
Figura 4.10:
y(k) Im
Plano z
b
1-a
0<a<1 Re
Figura 4.11:
bz z b
y(∞) = lim 1 − z −1 Y (z) = lim 1 − z −1 . =
z→1 z→1 z−a z−1 1−a
94 Capı́tulo 4: Análisis de los Sistemas Discretos de Control
Im
Plano z
b
1-a -1<a<0 Re
Figura 4.12:
La respuesta ante escalón en los sistemas de primer orden será más rápida
cuanto más cerca esté el polo del origen del plano z (esto se corresponde en el
plano s con los polos más alejados del eje imaginario).
Y (z) b
=
X(z) 1 + a1 z + a2 z −2
−1
4.4 Análisis de la Respuesta de los Sistemas Discretos 95
Figura 4.13:
ω2 n
G(s) =
s(s + 2ζωn )
Y (s) G(s) ω2 n
= = 2
X(s) 1 + G(s) s + 2ζωn + ω 2 n
96 Capı́tulo 4: Análisis de los Sistemas Discretos de Control
y(k)
Im
Plano z
Re
Figura 4.14:
Figura 4.15:
Las caracterı́sticas dinámicas que se pueden medir sobre la respuesta del sis-
tema se calculan a través de las siguientes ecuaciones:
4.4 Análisis de la Respuesta de los Sistemas Discretos 97
π−β
tr =
ωn 1 − ζ 2
π
tp =
ωn 1 − ζ 2
√−ζπ
µp = e 1−ζ 2
4
ts = (criterio 2%)
ζωn
3
ts = (criterio 5%) (4.7)
ζωn
Figura 4.16:
Figura 4.17:
98 Capı́tulo 4: Análisis de los Sistemas Discretos de Control
y(t)
Mp
0,05
1 o
0,02
0,5
td t
tr
tp
ts
Figura 4.18:
• Intervalo de subida nr
• Intervalo de pico np
• Intervalo de retardo nd
• Pico de sobreoscilación Mp
• Intervalo de establecimiento ns
y(k)
Mp
0,05
1 o
0,02
0,5
nd k
nr
np
ns
Figura 4.19:
tiguado. Este rango puede establecerse del 5% al 95%, o del 0% al 100% para un
sistema subamortiguado.
2. Intervalo de Pico ( np ): es el número de instantes que transcurren hasta
que la salida alcanza su máximo por encima del valor final.
3. Intervalo de Retardo ( nd ): es el número de instantes que transcurren hasta
alcanzar por primera vez la salida la mitad del valor final.
4. Pico de Sobreoscilación Mp : Es la amplitud sobrepasada con respecto
al valor final de la salida. Este parámetro se espresa en % y se puede obtener
mediante:
100 Capı́tulo 4: Análisis de los Sistemas Discretos de Control
y(np ) − y(∞)
Mp =
y(∞)
5. Intervalo de Establecimiento ns : es el número de instantes que transcurren
hasta alcanzar la salida el primer valor que queda acotado dentro de una franja
en torno al valor final. Esta franja se suele establecer en el 2% ó el 5%.
Figura 4.20:
Y (z) c
=
X(z) (z + a + jb)(z + a − jb)
Se supone que la entrada es un escalón unitario.
|P1 | = e−σT
θ = ωd .T (4.8)
4.4 Análisis de la Respuesta de los Sistemas Discretos 101
1. Intervalo de Subida ( nr ):
γ γ
nr = =
θ ωd T
2. Intervalo de Pico ( np ):
π π
np = =
θ ωd T
3. Pico de Sobreoscilación Mp :
πσ
Mp = e− tan θ .100% = e ωd .100% = |P1 |np .100%
π
4. Intervalo de Establecimiento ns :
y(ns ) − y(∞)
= ±0.05
y(∞)
π
ns =
σ.T
5. Ganancia Estática: es el valor de la salida en régimen permanente cuando
la entrada es la secuencia escalón unitario.
Aplicando el teorema del valor final, se obtiene:
1
y∞ = lim (1 − z −1 )G(z)( )=1
z→1 1 − z −1
P (z)
G(z) =
(z − 1)n Q(z)
P (z)
Nota: Q(z) no tiene polos en z=1
Se tiene entonces:
P (z)
Sistema tipo 0 → n = 0 → G(z) = Q(z)
P (z)
Sistema tipo 1 → n = 1 → G(z) = (z−1)Q(z)
P (z)
Sistema tipo 2 → n = 2 → G(z) = (z−1) 2 Q(z)
Figura 4.21:
Y (z) B0 G(z)
=
X(z) 1 + B0 GH(z)
La señal B(z) vale:
x(t) = 1(t)
x(k) = 1(k)
1
X(z) = (4.13)
1 − z −1
−1 1 1
e∞ = lim (1 − z )
z→1 1 + B0 GH(z) 1 − z −1
1
e∞ = lim (4.14)
z→1 1 + B0 GH(z)
Kp = lim B0 GH(z)
z→1
x(t) = t
x(k) = kT
T z −1
X(z) = (4.15)
(1 − z −1 )2
4.5 Error en Régimen Permanente 105
−1 1 T z −1
e∞ = lim (1 − z )
z→1 1 + B0 GH(z) (1 − z −1 )2
T z −1
e∞ = lim
z→1 (1 − z −1 ) + (1 − z −1 )B0 GH(z)
T
e∞ = lim (4.16)
z→1 (1 − z −1 )B0 GH(z)
(1 − z −1 )B0 GH(z)
Kv = lim
z→1 T
pudiéndose expresar el error de la forma:
1
e∞ =
Kv
Observar que el error ante entrada rampa se anulará cuando Kv = ∞, cosa
que ocurrirá cuando B0 GH(z) tenga un polo doble en z=1.
t2
x(t) =
2
(kT )2
x(k) =
2
T 2 (1 + z −1 )z −1
X(z) = (4.17)
2(1 − z −1 )3
1 T 2 (1 + z −1 )z −1
e∞ = lim (1 − z −1 )
z→1 1 + B0 GH(z) 2(1 − z −1 )3
T2
e∞ = lim (4.18)
z→1 (1 − z −1 )2 B0 GH(z)
(1 − z −1 )2 B0 GH(z)
Ka = lim
z→1 T2
106 Capı́tulo 4: Análisis de los Sistemas Discretos de Control
H(s)
Figura 4.22:
X(z)
E(z) =
1 + G(z)H(z)
X(z)
e∞ = lim (1 − z −1 )
z→1 1 + G(z)H(z)
1
e∞ =
1 + Kp
Kp = lim G(z)H(z)
z→1
1
e∞ =
Kv
(1 − z −1 )G(z)H(z)
Kv = lim
z→1 T
1
e∞ =
Ka
(1 − z −1 )2 G(z)H(z)
Ka = lim
z→1 T2
Ejemplo 4.3 Calcular el error en régimen permanente del sistema de la figura
4.23
H(s)
Figura 4.23:
X(z)
E(z) =
1 + G1 (z)G2 H(z)
−1 X(z)
e∞ = lim (1 − z )
z→1 1 + G1 (z)G2 H(z)
1
e∞ =
1 + Kp
1
e∞ =
Kv
4.5 Error en Régimen Permanente 109
(1 − z −1 )2 G1 (z)G2 H(z)
Ka = lim
z→1 T2
Ejemplo 4.4 Calcular el error en régimen permanente del sistema de la figura
4.24
Figura 4.24:
X(z)
E(z) =
1 + G1 (z)G2 (z)H(z)
X(z)
e∞ = lim (1 − z −1 )
z→1 1 + G1 (z)G2 (z)H(z)
1
e∞ =
1 + Kp
1
e∞ =
Kv
(1 − z −1 )2 G1 (z)G2 (z)H(z)
Ka = lim
z→1 T2
X(z) Y(z)
R(z) G(z)
+ +
- -
Figura 4.25:
X(z) Y1(z)
R(z) G(z)
+
-
Figura 4.26:
Y1 (z) G(z).R(z)
=
X(z) 1 + G(z).R(z)
El error debido a la entrada X(z) es:
112 Capı́tulo 4: Análisis de los Sistemas Discretos de Control
N(z) Y (z)
2
G(z)
+
-
R(z)
Figura 4.27:
Y2 (z) G(z)
=
N (z) 1 + G(z).R(z)
El error debido a la perturbación es:
H(s)
Figura 4.28:
B0 G(z)
M (z) =
1 + B0 GH(z)
Se sabe que la resolución de la ecuación caracterı́stica proporciona los polos
del sistema en lazo cerrado, pudiéndose expresar de la siguiente forma:
m
(z − zi )
1 + B0 GH(z) = 1 + K ni=1 (4.25)
j=1 (z − pj )
Resolviendo la ecuación caracterı́stica se obtiene la nueva posición de los polos
del sistema realimentado al variar el parámetro K. El lugar de las raı́ces propor-
ciona un método gráfico para la resolución de dicha ecuación y permite analizar
la evolución de los polos al variar la ganancia K o cualquier otro parámetro del
sistema.
Todo punto del plano z que cumpla en la ecuación caracterı́stica será un polo
del sistema realimentado, además cumplirá los llamados criterio del módulo y del
argumento:
n
j=1 (z − pj )
K = − m
i=1 (z − zi )
n
j=1 |(z − pj )|
|K| = m
i=1 |(z − zi )|
d1 .d2 .d3
K|z=p =
d4
Figura 4.29:
Este punto del plano deberá cumplir también el criterio del argumento:
n
m
(z − pj ) − (z − zi ) = ±180(2k + 1) k = 0, 1, 2, · · ·
j=1 i=1
que representa la diferencia entre los ángulos a los polos y los ceros desde el
punto del plano considerado, tal y como se muestra en la figura 4.29.
4.6 Lugar de las Raı́ces 115
α1 + α2 + α3 − β1 = ±180(2k + 1) k = 0, 1, 2, · · ·
2. El lugar de las raı́ces tendrá tantas ramas como polos tenga la función de
transferencia en lazo abierto.
4. Un punto del eje real pertence al lugar de las raı́ces si deja a la derecha un
número impar de polos y ceros
Teniendo que derivar esta expresión con respecto a z y calcular las raı́ces:
dK
=0
dz
Si la raı́z obtenida proporciona un valor positivo de K, dicho valor es un
punto de dispersión y/o de confluencia.
9. Las ramas del lugar de las raı́ces que salen de polos complejos o entran en
ceros complejos, lo hacen con un ángulo que se calcula por aplicación del
criterio del argumento, tal y como se muestra en la figura 4.30.
Im Im
Plano z Plano z
p As z All
1 2 1 1 2 1
Re Re
3 2
Figura 4.30:
Problemas Resueltos
Problemas Resueltos 119
Problema 4.1
• |an | < a0
Cumple por ser 0.02 < 1
• p(1) > 0
p(1) = 1 − 0.5 − 0.04 + 0.02 = 0.48
Cumple por ser 0.48 > 0
z0 z1 z2 z3
fila 1 0.02 −0.04 −0.5 1
fila 2 1 −0.5 −0.04 0.02
fila 3 b2 b1 b0 0
0.02
1
b2 =
= −0.99
1 0.02
0.02 −0.04
b0 =
= 0.03
1 −0.5
Problema 4.2
• |an | < a0
Cumple por ser 0.05 < 1
• p(1) > 0
p(1) = 1 + 0.6 + 0.05 = 1.55
Cumple por ser 1.55 > 0
Problema 4.3
p(z) = 1 + KG(z) = 0
z
p(z) = 1 + K =0
z 3 + 0.4z 2 + 0.47z + 0.13
• |an | < a0
Cumple por ser 0.13 < 1
• p(1) > 0
p(1) = 1 + 0.4 + (0.47 + K) + 0.13 > 0
2 + K > 0 luego K > −2 para estable.
z0 z1 z2 z3
fila 1 0.13 0.47 + K 0.4 1
fila 2 1 0.4 0.47 + K 0.13
fila 3 b2 b1 b0 0
0.13
1
b2 =
= −0.9831
1 0.13
0.13 0.47 + K
b0 =
= −0.418 − K
1 0.4
122 Problemas Resueltos
z z −2
M (z) = =
z 3 + 0.47z 2 + 1.47z + 0.13 1 + 0.4z −1 + 1.47z −2 + 0.13z −3
La ecuación en diferencias se podrá expresar a partir de la anterior de la
forma:
w3 0.48 3.6
w2 2.48 1.44
w a 0
w0 b
Problema 4.5
Y (z) 1
=
X(z) 1 − 0.6z − 0.81z + 0.67z −3 − 0.12z −4
−1 −2
Y (z) z4
= 4
X(z) z − 0.6z 3 − 0.81z 2 + 0.67z − 0.12
La condición previa a0 > 0 se cumple por ser a0 = 1
124 Problemas Resueltos
• |an | < a0
Cumple por ser 0.12 < 1
• p(1) > 0
p(1) = 1 − 0.6 − 0.81 + 0.67 − 0.12 = 0.14
Cumple por ser 0.14 > 0
• p(−1) > 0 por n=4 par
p(−1) = 1 + 0.6 − 0.81 − 0.67 − 0.12 = 0
Al ser 0, existe al menos una raı́z en z = −1
• Tabla de Jury con 2n − 3 = 5 filas:
z0 z1 z2 z3 z4
fila 1 −0.12 0.67 −0.81 −0.6 1
fila 2 1 −0.6 −0.81 0.67 −0.12
fila 3 b3 b2 b1 b0 −
fila 4 b0 b1 b2 b3 −
fila 5 c2 c1 c0 − −
−0.12
1
b3 =
= −0.9856
1 −0.12
−0.12 −0.6
b2 =
= 0.5196
1 0.67
−0.12 −0.81
b1 =
= 0.9072
1 −0.81
−0.12 0.67
b0 =
= −0.598
1 −0.6
b
b0
c2 =
3
= 0.6138
b0 b3
b
b2
c0 =
3
= −0.5834
b0 b1
Problemas Resueltos 125
Problema 4.6
1
M1 (z) =
z3 − 4.1z 2
+ 1.1z + 6.2
1
M2 (w) =
w − 0.5w − 3.34w2 + 3.08w − 0.48
4 3
w4 1 −3.34 −0.48
w3 −0.5 3.08 0
w2 a b 0
w c 0 −
w0 d 0 −
126 Problemas Resueltos
1
1 −3.34
a=−
= 2.82
(−0.5)
−0.5 3.08
1
1 −0.48
b=−
= −0.48
(−0.5)
−0.5 0
1
−0.5 3.08
c=−
= 2.99
a
a b
1
a b
d=−
= −0.48
c
c 0
Hay tres cambios de signo en la primera columna, luego el sistema tiene tres
polos inestables en las posiciones w1 = 1.5, w2 = 0.8, w3 = 0.2.
Problema 4.7
1
G1 (z) =
z + 0.5
1
G2 (z) =
1 − 0.5z −1
1
G3 (z) =
1 + 0.5z −1
1
1. Para G1 (z) = z+0.5
yk = −0.5yk−1 + xk−1
Resolviendo cuando la entrada es el escalón xk = {1, 1, 1, . . . 1},
y0 = −0.5y−1 + x−1 = 0
y1 = −0.5y0 + x0 = 1
Problemas Resueltos 127
y2 = −0.5y1 + x1 = 0.5
y3 = −0.5y2 + x2 = 0.75
y4 = −0.5y3 + x3 = 0.625
..
.
1 1
y1 (∞) = lim (1 − z −1 ) · G1 (z)X(z) = lim (1 − z −1 ) = 0.66
z→1 z→1 (z + 0.5) (1 − z −1 )
1 1
y1 (0) = lim G1 (z)X(z) = lim =0
z→∞ z→∞ (z + 0.5) (1 − z −1 )
1
2. Para G2 (z) = 1−0.5z −1
yk = 0.5yk−1 + xk
Por aplicación del teorema del valor inicial y final se obtiene,
1 1
y2 (0) = lim G2 (z)X(z) = lim =1
z→∞ z→∞ (1 − 0.5z −1 ) (1 − z −1 )
1 1
y2 (∞) = lim (1 − z −1 ) · G2 (z)X(z) = lim (1 − z −1 ) =2
z→1 z→1 (1 − 0.5z ) (1 − z −1 )
−1
1
3. Para G3 (z) = 1+0.5z −1
yk = −0.5yk−1 + xk
Por aplicación del teorema del valor inicial y final se obtiene,
Problemas Resueltos 129
1 1
y3 (0) = lim G3 (z)X(z) = lim =1
z→∞ z→∞ (1 + 0.5z ) (1 − z −1 )
−1
1 1
y3 (∞) = lim (1 − z −1 ) · G3 (z)X(z) = lim (1 − z −1 ) = 0.66
z→1 z→1 (1 + 0.5z ) (1 − z −1 )
−1
Problema 4.8
np2 = 2
ns2 = 5
La secuencia y3k :
1.5 − 1.09
Mp3 = = 37%
1.09
np3 = 1
ns3 = 5
1. Ordenando las caracterı́sticas Mp , np y ns :
Problema 4.9
Obtener:
0.74z + 0.52 1
y(∞) = lim (1 − z −1 ) · M (z)X(z) = lim (1 − z −1 ) =1
z→1 z→1 z2 − 0.63z + 0.89 (1 − z −1 )
siendo Kp ,
0.74z + 0.52
Kp = lim B0 G(z) = lim =∞
z→1 z→1 z 2 − 1.37z + 0.37
Por lo tanto,
1
e∞ = =0
1+∞
k = 2 → yk = 1.73, luego np = 2
5. Muestreando con T = 2,
−1 4z −1 2 2
B0 G(z) = (1 − z ) − + =
(1 − z −1 )2 (1 − z −1 ) (1 − 0.135z −1 )
Problema 4.10
X(z) + Y(z)
1 1
K z-0.5 z-0.7
-
1
z+0.1
Obtener:
K(z + 0.1)
M (z) =
(z − 0.5)(z − 0.7)(z + 0.1) + K
con ecuación caracterı́stica,
• |an | < a0
(0.035 + K) < 1 luego K < 0.965 para estable.
Problemas Resueltos 135
• p(1) > 0
p(1) = 1 − 1.1 + 0.23 + (0.035 + K) > 0
0.165 + K > 0 luego K > −0.165 para estable.
z0 z1 z2 z3
fila 1 (0.035 + K) 0.23 −1.1 1
fila 2 1 −1.1 0.23 (0.035 + K)
fila 3 b2 b1 b0 0
0.035 + K
1
b2 =
= (0.035 + K)2 − 1
1 0.035 + K
0.035 + K
0.23
b0 =
= −1.1(0.035 + K) − 0.23
1 −1.1
Calculando,
K1 < 0.45 cumple en
(0.035 + K)2 − 1
> |−1.1(0.035 + K) − 0.23|
K1 > 0.45 no cumple en
(0.035 + K)2 − 1
> |−1.1(0.035 + K) − 0.23|
Con todas las restricciones positivas de K se concluye que para K < 0.45 el
sistema será siempre estable.
1 0.1
e∞ = lim (1−z −1 )E(z) = lim (1−z −1 ) 1−
z→1 z→1 (1 − z −1 ) (z − 0.5)(z − 0.7)(z + 0.1) + 0.1
0.1
e∞ = lim 1 − = 0.623
z→1 (z − 0.5)(z − 0.7)(z + 0.1) + 0.1
Problema 4.11
- -
1
s
G(s)
4
s T = 0.7
H(s)
Obtener:
Problemas Resueltos 137
−1 s 1 1 − z −1
B0 Gdirecta (z) = (1 − z )Z = (1 − z −1 ) =
s(s + 1) 1 − e−T z −1 1 − 0.5z −1
4
H(z) =
1 − z −1
B0 Gdirecta (z) 1 − z −1
M (z) = =
1 + B0 Gdirecta (z)H(z) 5 − 0.5z −1
1 − z −1
y(0) = lim M (z)X(z) = lim = 0.2
z→∞ z→∞ (5 − 0.5z −1 )
1 − z −1
y(∞) = lim (1 − z −1 ) · M (z)X(z) = lim (1 − z −1 ) =0
z→1 z→1 (5 − 0.5z −1 )
Problema 4.12
-
H(s)
5
s
Obtener:
0.43z −1 5
p(z) = 1 + KB0 G(z)H(z) = 1 + K
(1 − 0.14z ) (1 − z −1 )
−1
z + 0.4
p(z) = 1 + K = z 2 + (K − 1.5)z + (0.4K + 0.5) = 0
z 2 − 1.5z + 0.5
2 2
T T T T
1+ w + (K − 1.5) 1 + w 1 − w + (0.4K + 0.5) 1 − w =0
2 2 2 2
(T w)2
(3 − 0.6K) + T w(0.5 − 0.4K) + 1.4K = 0
4
Al ser una ecuación de orden 2 bastará con imponer la condición de que todos
los términos en w sean mayor que cero, de manera que:
Problema 4.14
Y (z) z −2 1
M (z) = = −1 −2
= 2
X(z) 1 − z + 0.5z z − z + 0.5
p(z) = z 2 − z + 0.5 = 0
El sistema resulta ser estable con polos en las posiciones z1 = 0.5 + 0.5i y
z2 = 0.5 − 0.5i, ambos dentro del cı́rculo unidad.
z −2 1
y(0) = lim M (z)X(z) = lim =0
z→∞ z→∞ (1 − z −1 + 0.5z ) (1 − z −1 )
−2
z −2 1
y(∞) = lim (1−z −1 )·M (z)X(z) = lim (1−z −1 ) =2
z→1 z→1 (1 − z −1 + 0.5z −2 ) (1 − z −1 )
Problema 4.15
En el sistema de la figura,
Obtener:
2
R(z) =
1 − z −1
1 − e−T s
B0 (s) =
s
1
G(s) =
s+1
R(z)B0 G(z)
M (z) =
1 + R(z)B0 G(z)
siendo B0 G(z),
1 0.39z −1
B0 G(z) = (1 − z −1 )Z =
s(s + 1) 1 − 0.61z −1
2 · 0.39z −1
M (z) =
(1 − z −1 )(1 − 0.61z −1 ) + 2 · 0.39z −1
0.78z −1 0.78z
M (z) = −1 −2
= 2
1 − 0.83z + 0.61z z − 0.83z + 0.61
−1 ) 2 0.39z −1
(1 − z −1 )R(z)B0 G(z) (1 − z (1−z −1 ) (1−0.61z −1 )
Kv = lim = lim =4
z→1 T z→1 0.5
Problema 4.16
Se pide:
0.5z −2 0.5
G1 (z) = −1 −2
= 2
1 + 0.6z + 0.05z z + 0.6z + 0.05
z −2 + 0.5z −3 z + 0.5
G2 (z) = −1 −2
= 3
1 − 2.8z + 1.6z z − 2.8z 2 + 1.6z
Con los polos en z1 = −0.1 y z2 = −0.5, ambos dentro del cı́rculo unidad.
0.5 1
y1 (∞) = lim (1−z −1 )·G1 (z)X(z) = lim (1−z −1 ) = 0.303
z→1 z→1 (z 2 + 0.6z + 0.05) (1 − z −1 )
El segundo sistema es inestable con un polo fuera del cı́rculo unidad, por lo
tanto y2 (∞) = ∞
y1k = {0, 0, 0.5, 0.2, 0.35, 0.27, 0.32, 0.29, 0.31, 0.30, . . . 0.30}
y2k = {0, 0, 1, 4.3, 11.9, 28.1, 60.9, 127, 260, 526, 1059, . . . ∞}
146 Problemas Resueltos
Problema 4.17
1
B0(s)
s T=1
H(s)
Obtener:
B0 G(z)
M (z) =
1 + B0 G(z)B0 H(z)
0.4325z −2
B0 G(z)B0 H(z) =
1 − 1.135z −1 + 0.1353z − 2
3. Los polos conjugados dentro del semicı́rculo positivo dan una respuesta
subamortiguada, tal y como se muestra en la figura.
Problema 4.18
G1(s) G2(s)
−1 2z −1 2 2
B0 Gdirecta (z) = (1 − z ) −1 2
− + =
(1 − z ) (1 − z ) (1 − 0.37z −1 )
−1
(0.74z −1 + 0.52z −2 )
Kp = lim B0 Gdirecta (z) = lim =∞
z→1 z→1 (1 − 1.37z −1 + 0.37z −2 )
(0.74z + 0.52) 1
y(∞) = lim (1 − z −1 ) =1
z→1 (z 2− 0.63z + 0.89) (1 − z −1 )
e∞ = x(∞) − y(∞) = 1 − 1 = 0
150 Problemas Resueltos
(0.74z −1 + 0.52z −2 )
B0 Gdirecta (z) =
(1 − z −1 )(1 − 0.37z −1 )
(1 − z −1 )(0.74z −1 + 0.52z −2 )
Kv = lim =2
z→1 (1 − z −1 )(1 − 0.37z −1 )
Problema 4.19
Se desea:
X(z) + Y1(z)
0.1 1
z-0.2 (z-0.1)(z-0.2)
-
0.1
RG(z) (z−0.2)(z−0.1)(z−0.2) 0.1
M1 (z) = = 0.1 = 2
1 + RG(z) 1 + (z−0.2)(z−0.1)(z−0.2) (z − 0.2) (z − 0.1) + 0.1
N(z) + Y2(z)
1
(z-0.1)(z-0.2)
-
R(z)
0.1
z-0.2
1
G(z) (z−0.1)(z−0.2) (z − 0.2)
M2 (z) = = 0.1 =
1 + RG(z) 1+ (z−0.2)(z−0.1)(z−0.2)
(z − 0.2)2 (z − 0.1) + 0.1
Luego,
0.1 1 (z − 0.2) z −3
Y (z) = · + ·
(z − 0.2)2 (z − 0.1) + 0.1 (1 − z −1 ) (z − 0.2)2 (z − 0.1) + 0.1 (1 − z −1 )
Problema 4.20
0.2
z-0.7
+
H(z)
1
-2
z
N(z)
G(z)
Y1(z)
X(z) + z-0.2
z(z-0.5)
-
0.2
z-0.7
H(z)
(z−0.2)
G(z) z(z−0.5) (z − 0.2)(z − 0.7)
M1 (z) = = 0.2(z−0.2)
=
1 + GH(z) 1+ z(z − 0.5)(z − 0.7) + 0.2(z − 0.2)
z(z−0.5)(z−0.7)
0.2
z-0.7
+
H(z)
1
z-2
N(z)
1 z −2 1 1
N (z) · = · −2 =
z −2 −1
(1 − z ) z (1 − z −1 )
Aplicando superposición,
2 1 1
Y (z) = M1 (z) −1
+ M2 (z) −1
= M1 (z)
(1 − z ) (1 − z ) (1 − z −1 )
(z − 0.2)(z − 0.7) 1
y(∞) = lim (1 − z −1 ) = 0.77
z→1 z(z − 0.5)(z − 0.7) + 0.2(z − 0.2) (1 − z −1 )
−1 1 0.2 1
e2(∞) = − lim (1 − z ) + M2 (z) =
z→1 (1 − z −1 ) (z − 0.7) (1 − z −1 )
Problemas Resueltos 155
0.2 (z − 0.2)(z − 0.7)
= − lim 1 − · = −0.484
z→1 (z − 0.7) z(z − 0.5)(z − 0.7) + 0.2(z − 0.2)
Problema 4.21
G1(s) G2(s)
x(t) y(t)
1 1
B0(s)
T=1 s+1 s
Y1 (z) B0 G1 G2 (z)
M1 (z) = =
X(z) 1 + B0 G1 G2 (z)
156 Problemas Resueltos
−1 1 0.37z + 0.26
B0 G1 G2 (z) = (1 − z )Z 2 = 2
s (s + 1) z − 1.37z + 0.37
Luego,
0.37z + 0.26
M1 (z) =
z 2 − z + 0.63
El valor de la salida en régimen permanente debido sólo a la entrada x(t)
valdrá:
G1(s) G2(s)
y2(t)
1 1
B0(s)
T=1 s+1 s
N G2 (z)
Y2 (z) =
1 + B0 G1 G2 (z)
0.2 0.2z −1 0.2z
N G2 (z) = Z 2 = −1 2
=
s (1 − z ) (z − 1)2
Luego,
E2 (z) = 0 − Y2 (z)
Aplicando el teorema del valor final,
Problema 4.22
N(z)
+ Y(z)
X(z) +
1 + z-0.1
z z-0.5
-
0.1
z-0.6
X(z) + Y1(z)
1 z-0.1
z z-0.5
-
0.1
z-0.6
+ Y2(z)
1 - z-0.1
z z-0.5
0.1
z-0.6
Problemas Resueltos 159
Aplicando superposición,
y(∞) = lim (1 − z −1 ) [M1 (z)X(z) + M2 (z)N (z)] = lim [M1 (z) + M2 (z)]
z→1 z→1
(z − 0.1)(z − 0.6) z(z − 0.1)(z − 0.6)
y(∞) = lim +
z→1 z(z − 0.5)(z − 0.6) + 0.1(z − 0.1) z(z − 0.5)(z − 0.6) + 0.1(z − 0.1)
y(∞) = 2.48
Por lo tanto,
1
e1(∞) = = 0.69
1 + 0.45
0.1 0.1
E2 (z) = 0 − · Y2 (z) = − · M2 (z)N (z)
(z − 0.6) (z − 0.6)
−1 0.1 z(z − 0.1)(z − 0.6) 1
e2(∞) = − lim (1 − z ) ·
z→1 (z − 0.6) z(z − 0.5)(z − 0.6) + 0.1(z − 0.1) (1 − z −1 )
0.1 z(z − 0.1)(z − 0.6)
e2(∞) = − lim · = −0.31
z→1 (z − 0.6) z(z − 0.5)(z − 0.6) + 0.1(z − 0.1)
Problema 4.23
Y(z)
X(z) + 1 z-0.1
z z-0.5
-
+
0.1
z-0.6
+
1
z-3
N(z)
Y1(z)
X(z) + 1 z-0.1
z z-0.5
-
0.1
z-0.6
0.1
z-0.6
+
1
-3
z
N(z)
1 z −3 1 1
N (z) · = · =
z −3 (1 − z −1 ) z −3 (1 − z −1 )
Aplicando superposición,
1 1
Y (z) = Y1 (z) + Y2 (z) = M1 (z) −1
+ M2 (z)
(1 − z ) (1 − z −1 )
(z − 0.1)(z − 0.6) −0.1(z − 0.1)
y(∞) = lim +
z→1 z(z − 0.5)(z − 0.6) + 0.1(z − 0.1) z(z − 0.5)(z − 0.6) + 0.1(z − 0.1)
y(∞) = 0.93
Por lo tanto,
1
e1(∞) = = 0.69
1 + 0.45
0.1 0.1
E2 (z) = 0 − [X(z) + Y2 (z)] = − X(z) [1 + M2 (z)]
(z − 0.6) (z − 0.6)
−1 0.1 −0.1(z − 0.1)
e2(∞) = lim (1 − z ) − X(z) 1 +
z→1 (z − 0.6) z(z − 0.5)(z − 0.6) + 0.1(z − 0.1)
0.1 −0.1(z − 0.1)
e2(∞) = lim − 1+ = −0.17
z→1 (z − 0.6) z(z − 0.5)(z − 0.6) + 0.1(z − 0.1)
Problema 4.24
• Polos y ceros en lazo abierto: tres polos (n = 3) en las posiciones del plano
z1 = 0.2, z2 = 0.7, z3 = 0.9, y un cero (m = 1) en la posición z1 = 0.
0.2+0.7+0.9−0
• Centroide: σ = 2 = 0.9
164 Problemas Resueltos
Problemas Resueltos 165
Problema 4.25
Se desea:
• Polos y ceros en lazo abierto: dos polos (n = 2) en las posiciones del plano
z1 = 0.37, z2 = 1, y un cero (m = 1) en la posición z1 = −0.7.
(z − 1)(z − 0.37)
K=−
0.37(z + 0.7)
0.37(z + 0.7)
1+ = 0 ⇒ z 2 − z + 0.63 = 0
(z − 1)(z − 0.37)
Cuya solución da dos polos complejos conjugados en las posiciones z1,2 =
0.5 ± 0.62i, y que al quedar dentro del cı́rculo unidad (|z1,2 | = 0.8) el sistema será
estable.
π π
np = = = 3.5 Se da entre el instante 3 y 4.
θ 0.89
π π
ns = = = 14.28 Se da entre el instante 14 y 15.
σT 0.22 · 1
(0.37z + 0.26) 1
y(∞) = lim (1 − z −1 ) =1
z→1 (z 2 − z + 0.63) (1 − z −1 )
e∞ = x( ∞) − y( ∞) = 1 − 1 = 0
Problema 4.26
G(z)
Y(z)
X(z) + z+0.5
z2+0.25
-
1
z
H(z)
• Polos y ceros en lazo abierto: tres polos (n = 3) en las posiciones del plano
z1 = 0, z2,3 = ±0.5i, y un cero (m = 1) en la posición z1 = −0.5.
• Angulo de salida del polo complejo: calculando previamente los ángulos que
forman el polo complejo con el resto de polos y ceros se obtiene,
0.5 X
X
-0.5
-0.5 X
(z + 0.5)
1+ = 0 ⇒ z 3 + 1.25z + 0.5 = 0
z(z 2 + 0.25)
Cuya solución da tres polos en las posiciones z1,2 = 0.18 ± 1.16i, z3 = −0.36,
de los cuales la pareja de complejos conjugados queda fuera del cı́rculo unidad,
por lo tanto el sistema es inestable.
Problema 4.27
R(z) G(z)
Y(z)
X(z) +
1 1
z z2+0.3z-0.4
-
Se desea:
3. Con el valor anterior, determinar las nuevas posiciones de los polos en lazo
cerrado.
• Polos y ceros en lazo abierto: tres polos (n = 3) en las posiciones del plano
z1 = 0, z2 = 0.5, z3 = −0.8, y ningún cero (m = 0).
0+0.5−0.8
• Centroide: σ = 3 = −0.1
Resolviendo se obtienen los dos puntos del plano z1 = 0.28, z2 = −0.48 que
proporcionan un valor de K:
Cuya solución da los tres polos en las posiciones z1,2 = 0.35 ± 0.42i, z3 = −1.
Problema 4.28
G(z)
Y(z)
X(z) + z+0.5
2
z -z+0.5
-
• Polos y ceros en lazo abierto: dos polos (n = 2) en las posiciones del plano
z1,2 = 0.5 ± 0.5i, y un cero (m = 1) en la posición z1 = −0.5.
(z 2 − z + 0.5)
K=−
(z + 0.5)
Resolviendo se obtienen los dos puntos del plano z1 = 0.6, z2 = −1.6 que
proporcionan unos valores de K:
(z12 − z1 + 0.5)
K|z=z1 = − = −0.24
(z1 + 0.5)
(z22 − z2 + 0.5)
K|z=z2 = − = 4.2
(z2 + 0.5)
Problemas Resueltos 175
• Angulo de salida del polo complejo: calculando previamente los ángulos que
forman el polo complejo con el resto de polos y ceros se obtiene,
0.5 X
0.5
-0.5
-0.5 X
(z + 0.5)
1+ = z2 + 1 = 0
(z 2 − z + 0.5)
Cuya solución da los polos en las posiciones z1,2 = ±i. Por lo tanto el sistema
es crı́ticamente estable por tener los polos sobre el cı́rculo unidad.
Problema 4.29
G(z)
Y(z)
X(z) + 1
2 2
(z -0.25)(z +0.25)
-
0.5−0.5+0.5i−0.5i
• Centroide: σ = 4 =0
(z 2 − 0.25)(z 2 + 0.25)
K=− = −(z 4 − 0.0625)
1
dK
= −4z 3 = 0
dz
• Angulo de salida del polo complejo: calculando previamente los ángulos que
forman el polo complejo con el resto de polos y ceros se obtiene,
0.5 X
X X
-0.5 0.5
-0.5 X
b X Kcrítica
0.5 X d4
d1
d3
45º
X X
-0.5 0.5 a
d2
-0.5 X
La aplicación del critero del modulo en ese punto del lugar dará el valor de K
buscado:
Kcrítica = d1 · d2 · d3 · d4
Por trigonometrı́a:
a = b = 1 · sen(45o ) = 0.71
d1 = (a − 0.5)2 + b2 = 0.74
d4 = a2 + (b − 0.5)2 = 0.74
d2 = a2 + (0.5 + b)2 = 1.4
d3 = (0.5 + a)2 + b2 = 1.4
Luego,
Problema 4.30
G(z)
Y(z)
X(z) + 1
z2-z+0.74
-
Se desea:
3. Con el valor anterior, determinar las nuevas posiciones de los polos en lazo
cerrado.
• Polos y ceros en lazo abierto: dos polos (n = 2) en las posiciones del plano
z1,2 = 0.5 ± 0.7i y ningún cero (m = 0).
0.5+0.7i+0.5−0.7i
• Centroide: σ = 2 = 0.5
Problemas Resueltos 181
Kcrítica
a X
d1
0.7 X
d2
0.5
-0.7 X
La aplicación del critero del módulo en ese punto del lugar dará el valor de K
buscado:
Kcrítica = d1 · d2
Por trigonometrı́a:
a= 12 − 0.52 = 0.867
d1 = a − 0.7 = 0.167
d2 = a + 0.7 = 1.567
Luego,
3. Los polos en lazo cerrado están en las posiciones z1,2 = 0.5 ± 0.867i
Aplicando la transformación z = eT s ,
σ = ζωn = 0.872
ωd = ωn 1 − ζ 2 = 9.937
0.872
ζ=
ωn
0.8722
ωd = ωn2 1− = 9.9372 ⇒ ωn2 − 0.76 = 98.74
ωn2
π π
np = = = 3.16 Se da aproximadamente en el instante 3.
θ 0.9937
π π
ns = = = 36 Se da en el instante 36.
σT 0.872 · 0.1
Problema 4.31
G(z)
Y(z)
X(z) + (z-0.7)
2
(z-0.9)(z -1.4z+0.98)
-
Se desea:
• Polos y ceros en lazo abierto: tres polos (n = 3) en las posiciones del plano
z1 = 0.9, z2,3 = 0.7 ± 0.7i, y un cero (m = 1) en z1 = 0.7.
0.9+0.7+0.7i+0.7−0.7i−0.7
• Centroide: σ = 2 = 0.8
• Angulo de salida del polo complejo: calculando los ángulos a polos y ceros
tal y como aparecen en el plano z de la figura se obtiene,
0.7 X
X
0.7 0.9
-0.7 X
(z − 0.7)
1+ = z 3 − 2.3z 2 + 3.24z − 1.582 = 0
(z − 0.9)(z 2 − 1.4z + 0.98)
Problema 4.32
Se desea:
dK
1. Dibujar el lugar de las raı́ces. (Nota: dz = 0, tiene una raı́z en z=0.176)
2. Determinar el valor de k para obtener una respuesta con coeficiente de
amortiguamiento (ζ) igual a uno, obteniendo las nuevas posiciones de los
polos en lazo cerrado.
3. calcular el error de posición en régimen permanente suponiendo una entrada
de valor 4.
4. ¿pertenece el punto p = 0.861 + j0.5141 al lugar de las raı́ces?.
5. calcular K en el punto p = 0.861 + j0.5141.
• Polos y ceros en lazo abierto: dos polos (n = 2) en las posiciones del plano
z1 = 0.0, z1 = 0.5, z1 = 1.1, y un cero (m = 1) en la posición z1 = −0.5.
Problemas Resueltos 187
(z 3 − 1.6z 2 + 0.55z)
K=−
(z + 0.5)
(z 3 − 1.6z 2 + 0.55z)
K=−
(z + 0.5)
(z 3 − 1.6z 2 + 0.55z)
K|z=z1 = − = −0.0779
(z + 0.5)
(z 3 − 1.6z 2 + 0.55z)
K|z=z2 = − = 0.0557
(z + 0.5)
(z22 − z2 + 0.5)
K|z=z3 = − = −6.2752
(z2 + 0.5)
5RRW/RFXV
/XJDUGHODVUDLFHV
,PDJLQDU\$[LV
5HDO$[LV
multiples y por tanto es criticamente amortiguado. Esto quiere decir que el valor
de k es igual a 0.557.
3.
(z + 0.5)
Kp = lim G(z)H(z) = lim =∞
z→1 z→1 z(z − 0.5)(z − 1.1)
que representa la diferencia entre los ángulos a los polos y los ceros desde el
punto del plano considerado, tal y como se muestra en la figura 4.31.
θ1 + θ2 + θ3 − α1 = ±180(2k + 1) k = 0, 1, 2, · · ·
0.5141 0.5141 0.5141 0.5141
(180 − arctan( )) + arctan( ) + arctan( ) − arctan( ) ± 180
0.259 0.361 0.861 1.361
(4.27)
n
j=1 |(z − pj )|
|K| = m
i=1 |(z − zi )|
190 Problemas Resueltos
Figura 4.31:
Obtener:
3. Error de posición.
Sampling time: 1
0.5
−0.5
−1
−1.5
−2.5 −2 −1.5 −1 −0.5 0 0.5 1
x(t) + 2 1-e
-Ts
1 y(t)
-1
T = 0.5 1-z s s+1
-
Ejercicios con MATLAB 193
Polos en lazo cerrado dentro del cı́rculo unidad, luego el sistema es estable.
Ejercicio 4.3 Para el sistema realimentado de la figura, en el que todos los blo-
ques son discretos, la señal de entrada X(z) es el escalón unitario y la señal
perturbación N(z) es también un escalón unitario pero retrasado tres instantes,
obtener el valor de la salida y el error en régimen permanente:
N(z)
R(z) G(z)
+ Y(z)
X(z) + 0.1 + 1
z-0.2 (z-0.1)(z-0.2)
-
>> ef = 1-yf
ef =
-0.3333
Ejercicios con MATLAB 195
Y(z)
X(z) + 1 z-0.1
z z-0.5
-
+
0.1
z-0.6
+
1
z-3
N(z)
1
z
H(z)
1.5
0.5
−0.5
−1
−1.5
−2
−2.5
−1 −0.5 0 0.5 1
Dos polos en lazo cerrado están fuera del cı́rculo unidad, luego el sistema es
inestable.
199
Capı́tulo 5
Diseño de Reguladores
Discretos
5.1 Introducción
El comportamiento de los sistemas reales no es ideal, están sometidos a pertur-
baciones, presentan errores y en muchos casos caracterı́sticas no lineales en su
funcionamiento. Para hacer frente a estos aspectos el presente capı́tulo tratará
el análisis y diseño de los reguladores discretos.
En la teorı́a de los sistemas continuos, el bloque denominado controlador o
regulador ( ver figura 5.1), y que en el esquema realimentado se sitúa entre la
señal de error y el actuador, es el encargado de conseguir que la respuesta del
sistema cumpla los requerimientos impuestos por el diseñador, entre los que cabe
citar la minimización del error del sistema, limitación del pico de sobreoscilación,
del tiempo de pico, tiempo de subida o tiempo de establecimiento para conseguir,
por ejemplo, respuestas más rápidas o más lentas.
En los sistemas continuos el regulador por excelencia y que mejores resultados
ofrece es el PID. Este regulador implementa tres acciones de control: la acción
proporcional, acción integral y acción derivativa. La ecuación que describe el
funcionamiento del PID estándar es la siguiente:
t
1 de(t)
u(t) = K e(t) + e(t)dt + Td (5.1)
Ti 0 dt
siendo,
200 Capı́tulo 5: Diseño de Reguladores Discretos
Figura 5.1:
Figura 5.2:
Figura 5.3:
−1 R(s)
Z[R(s)] = R(z) = Z[B0 (s)R(s)] = (1 − z )Z
s
Figura 5.4:
Ts + 1
Z[R(s)] = R(z) = Z[B1 (s)R(s)] = (1 − z −1 )2 Z R(s)
T s2
de e(k) − e(k − 1)
=
dt T
Aplicando la transformada de Laplace al término de la izquierda de la ex-
presión anterior:
de
L(
) = sE(s) (5.2)
dt
Aplicando la transformada Z al término de la derecha:
e(k) − e(k − 1) 1 − z −1
Z( )= E(z) (5.3)
T T
Igualando (5.2) y (5.3) se obtiene:
1 − z −1
s=
T
1−z −1
De manera que, haciendo el cambio s = T , se puede obtener una expresión
aproximada para el regulador discreto.
Este método no pone restricción alguna a la relación entre el orden del nu-
merador y el orden del denominador del regulador continuo R(s).
GE = lim R(s)
s→0
GE = lim R(z)
z→1
Figura 5.5:
KI
U (z) = Kp + + Kd (1 − z −1 ) E(z)
1 − z −1
Y su función de transferencia:
U (z) KI
= R(z) = Kp + + Kd (1 − z −1 ) (5.4)
E(z) 1 − z −1
siendo,
KT
Kp = K − ganancia proporcional
2Ti
KT
KI = ganancia integral
Ti
KTd
KD = ganancia derivativa (5.5)
T
206 Capı́tulo 5: Diseño de Reguladores Discretos
Figura 5.6:
Kp (1 − z −1 ) + KI + KD (1 − z −1 )2
R(z) =
1 − z −1
(Kp + KI + KD ) − (Kp + 2KD )z −1 + KD z −2
=
1 − z −1
−1
q 0 − q1 z + q2 z −2
= (5.6)
1 − z −1
con,
q0 = Kp + KI + KD
q1 = Kp + 2KD
q 2 = KD (5.7)
5.3 Reguladores PID Discretos 207
U (z) = K.E(z)
U (z)
R(z) = =K (5.8)
E(z)
Aporta una ganancia K al resto del sistema controlado. Siempre que se diseñe
un regulador para controlar un sistema se empezará con un tipo P, por ser su
implementación muy sencilla y sólo implica ajustar una ganancia. En muchos
casos la simple variación de la ganancia va a conseguir modificar el régimen
transitorio y permanente de la respuesta de la salida cumpliendo los criterios
impuestos por el diseñador.
Td
U (z) = K 1 + (1 − z −1 ) E(z)
T
T + Td Td −1
R(z) = K − z
T T
T +T
T
d
z − Td
T
= K (5.9)
z
Este regulador aporta un polo en el origen y un cero ajustable al resto del
sistema controlado. La posición del cero en el plano z será:
Td
z=
T + Td
Este cero modificará el comportamiento transitorio de la respuesta del sis-
tema.
T 1 + z −1
U (z) = K 1 + ) E(z)
2Ti 1 − z −1
(2Ti + T ) + (T − 2Ti )z −1
R(z) = K
2Ti (1 − z −1 )
208 Capı́tulo 5: Diseño de Reguladores Discretos
(2Ti + T )z + (T − 2Ti )
= K
2Ti (z − 1)
(5.10)
Este regulador aporta un polo en z=1 y un cero ajustable al resto del sistema
controlado. La posición del cero en el plano será:
2Ti − T
z=
2Ti + T
La acción integral, con la adición de un polo en z=1, hace que aumente el tipo
del sistema mejorarando la precisión del mismo; esto permite corregir el error que
presenta el sistema en régimen permanente. El cero, por el contrario va a afectar
al régimen transitorio, pudiendo complicar la estabilidad relativa del sistema.
q0 − q1 z −1 + q2 z −2
R(z) =
1 − z −1
2
q0 z − q1 z + q2
=
z(z − 1)
(5.11)
Este regulador aporta al resto del sistema un polo en z=1, otro en el origen,
y dos ceros ajustables. El correcto ajuste de los ceros permitirá al diseñador
modificar a voluntad el régimen transitorio y permanente de la respuesta del
sistema.
2. Localizar estos polos sobre el lugar de las raı́ces del sistema, y determinar
si éste pasa por los polos dominantes, gráficamente o mediante el criterio
del argumento. En caso afirmativo bastará con ajustar el parámetro K, es
decir, incluir un regulador tipo P. Si se cumplen las condiciones estáticas
(error), el diseño del regulador finaliza con la obtención de un regulador de
tipo P.
4. En el caso de que el lugar de las raı́ces no pase por los polos dominantes,
se diseñará un regulador tipo PD (polo en el origen y un cero ajustable).
La posición del cero será tal que cumpla el criterio del argumento. A con-
tinuación habrá que volver a ajustar el valor de K mediante el criterio del
módulo.
Figura 5.7:
Problemas Resueltos 211
Problemas Resueltos
212 Problemas Resueltos
Problema 5.1
1. Aproximación derivada.
2. Método Tustin.
R(s) G(s)
x(t) + y(t)
2.84 (s+0.4) 2
s (s+0.5)(s+2)
-
−1 2 0.0092(z + 0.92)
Z[B0 G(s)] = B0 G(z) = (1 − z )Z =
s(s + 0.5)(s + 2) (z − 0.95)(z − 0.82)
1 − z −1
s= = 10(1 − z −1 )
T
Sustituyendo en la expresión de R(s) se obtiene:
2 (1 − z −1 ) (1 − z −1 )
s= · = 20 ·
T (1 + z −1 ) (1 + z −1 )
Sustituyendo en la expresión de R(s) se obtiene:
Respuesta al escalón
1.4
1.2
Derivada
0.8
Sistema continuo
Tustin
0.6
0.4
0.2
0
0 1 2 3 4 5 6 7
214 Problemas Resueltos
Problema 5.2
R(s) G(s)
x(t) + y(t)
1.5 (s+2.25) 2
(s+3.5) s (s+1)
-
−1 2 0.0097(z + 0.97)
Z[B0 G(s)] = B0 G(z) = (1 − z )Z 2 =
s (s + 1) (z − 1)(z − 0.90)
(z − 0.798)
R(z) = 1.41 ·
(z − 705)
Respuesta al escalón
1.4
Sistema discreto equivalente
1.2
0.8
Sistema continuo
0.6
0.4
0.2
0
0 2 4 6 8
Problema 5.3
G(z)
X(z) + Y(z)
1
R(z) (z-0.5)(z-0.7)
-
216 Problemas Resueltos
3
ts = = 1.83 s ⇒ σ = 1.64 rad/s
σ
1.64
σ = ζ · ωn ⇒ ωn = = 8.2 rad/s
0.2
ωd = ωn 1 − ζ 2 = 8.2 1 − 0.22 = 8.03 rad/s
θ = T ωd = 0.8 rad
El lugar de las raı́ces del sistema dado muestra como el polo dominante es un
punto del lugar.
Sólo será necesario ajustar una K mediante el criterio del módulo para que el
sistema sitúe sus polos en lazo cerrado en la posición del polo dominante.
Pd
0.6
d2 d1
X X
0.5 0.7
-0.6
Problema 5.4
Diseñar el regulador discreto R(z) más sencillo que cumpla las siguientes
caracterı́sticas dinámicas para el sistema mostrado en la figura:
G(z)
X(z) + Y(z)
1
R(z) (z-0.5)(z-0.7)
-
0.37 1
y(∞) = lim (1 − z −1 ) = 0.71
z→1 (z − 0.5)(z − 0.7) + 0.37 (1 − z −1 )
La condición del error no se cumple, luego habrá que añadir una acción in-
tegral mediante un regulador PI, el cual introduce un polo en z = 1 y un cero
situado a 1/3 de la distancia entre el polo dominante y el cı́rculo unidad.
Problemas Resueltos 219
di
Pd
-T
e
-1 1
Pd
0.6
d1
d3
d2 d4
X X oX
0.5 0.7
c
-0.6
d1 = (1 − 0.6)2 + 0.62 = 0.72
d2 = d3 = 0.12 + 0.62 = 0.61
d4 = (0.94 − 0.6)2 + 0.62 = 0.69
(z − 0.94)
R(z) = 0.39 ·
(z − 1)
Problema 5.5
G(z)
X(z) + Y(z)
(z+0.2)
R(z) (z-0.5)(z-0.7)
-
3
ts = = 0.38 s ⇒ σ = 7.89 rad/s
σ
7.89
σ = ζ · ωn ⇒ ωn = = 13.15 rad/s
0.6
ωd = ωn 1 − ζ 2 = 13.16 1 − 0.62 = 10.53 rad/s
θ = T ωd = 1.05 rad
El lugar de las raı́ces del sistema dado muestra como el polo dominante no es
un punto del lugar de las raı́ces.
222 Problemas Resueltos
Hay que incorporar un regulador tipo PD para que las nuevas ramas del lugar
pasen por el polo dominante. Añadiendo un polo en el origen y un cero ajustable
se tendrá un nuevo mapa de polos y ceros como el de la figura.
Problemas Resueltos 223
Pd
0.39
d
d4 d3 d
d5
o X o X X
-0.2 0.22 0.5 0.7
-0.39
a
θ1 + θ2 + θ3 − α1 − α2 = 180o ⇒ α2 = 104.28o
d1 = (0.7 − 0.22)2 + 0.392 = 0.62
d2 = (0.5 − 0.22)2 + 0.392 = 0.48
d3 = 0.392 + 0.222 = 0.45
d4 = (0.2 + 0.22)2 + 0.392 = 0.57
d5 = (a − 0.22)2 + 0.392 = 0.40
Con el criterio del módulo, el valor global que debe tomar K en ese punto del
lugar será:
(z − 0.32)
R(z) = 0.59 ·
z
d6 = (1 − 0.22)2 + 0.392 = 0.87
d7 = (0.82 − 0.22)2 + 0.392 = 0.72
(z − 0.32)(z − 0.82)
R(z) = 0.71 ·
z(z − 1)
La figura adjunta muestra las respuestas del sistema sin regulador y con el
regulador PID diseñado
226 Problemas Resueltos
con R(z)
1.2
0.8
0.6
sin R(z)
0.4
0.2
0
0 0.5 1 1.5 2 2.5 3
Problema 5.6
ts |±2% ≈ 2 s ζ = 0.5
4
ts = = 2 s ⇒ σ = 2 rad/s
σ
Problemas Resueltos 227
2
σ = ζ · ωn ⇒ ωn = = 4 rad/s
0.5
ωd = ωn 1 − ζ 2 = 4 1 − 0.52 = 3.46 rad/s
θ = T ωd = 0.69 rad
228 Problemas Resueltos
Pd
0.43
o X X
-0.88 0.52 0.67
-0.43
0.43
θ1 = 180 − arctan = 138.17o
(1 − 0.52)
0.43
θ2 = 180 − arctan = 108.95o
(0.67 − 0.52)
0.43
α = arctan = 17.07o
(0.88 + 0.52)
θ1 + θ2 − α = ±180o
por tanto,
Luego el polo dominante no pertenece al lugar de las raı́ces. Hay que incor-
porar un regulador tipo PD para que las nuevas ramas del lugar pasen por el
polo dominante. Añadiendo un polo en el origen y un cero ajustable se tendrá
un nuevo mapa de polos y ceros como el de la figura.
Problemas Resueltos 229
Pd
0.43 d
d5 d3 d4 d
o X o X X
-0.88 0.52 0.67
a
-0.43
θ1 + θ2 + θ3 − α1 − α2 = ±180o ⇒ α2 = 89.64o
De manera que la posición que ocupa el cero sobre el eje real será:
0.43
tan α2 = tan 89.64 = 159.15 = ⇒ a = 0.52
(0.52 − a)
d3 = 0.522 + 0.432 = 0.67
d4 = (0.52 − a)2 + 0.432 = 0.43
230 Problemas Resueltos
d5 = (0.88 + 0.52)2 + 0.432 = 1.46
Con el criterio del módulo, el valor global que debe tomar K en ese punto del
lugar será,
(z − 0.52)
R(z) = 17.45 ·
z
Problema 5.7
G(z)
Y(z)
X(z) +
z+0.5
R(z) z (z-1)
-
4
ts = = 1.03 s ⇒ σ = 3.88 rad/s
σ
Problemas Resueltos 231
3.88
σ = ζ · ωn ⇒ ωn = = 16.167 rad/s
0.24
ωd = ωn 1 − ζ 2 = 16.167 1 − 0.242 = 15.69 rad/s
θ = T ωd = 1.569 rad
Pd
0.678
d3 d1
d2
o X X
-0.5
-0.678
0.678
θ1 = 180 − arctan = 145.86o
1
θ2 = 90o
0.678
α = arctan = 53.59o
0.5
θ1 + θ2 − α = ±180o
por tanto,
d1 = 12 + 0.6782 = 1.208
Problemas Resueltos 233
d2 = 0.678
d3 = 0.52 + 0.6782 = 0.842
d1 · d 2 1.208 · 0.678
K= = = 0.973
d3 0.842
Ası́ pues,
1
e∞ = =0
1+∞
Problema 5.8
ts |±2% ≈ 5 s ζ = 0.5
G(z)
X(z) + Y(z)
(z+0.5)
R(z) (z-0.1)(z-0.8)
-
234 Problemas Resueltos
4
ts = = 5 s ⇒ σ = 0.8 rad/s
σ
0.8
σ = ζ · ωn ⇒ ωn = = 1.6 rad/s
0.5
ωd = ωn 1 − ζ 2 = 1.6 1 − 0.52 = 1.39 rad/s
θ = T ωd = 1.39 rad
Problemas Resueltos 235
θ1 + θ2 − α = 180o
Pd
0.49
o X X
-0.5 0.1 0.8
-0.49
por tanto,
Pd
0.49
d3 d
d4
d d5
o XX o X
-0.5 0.1 0.8
θ1 + θ2 + θ3 − α1 − α2 = 180o ⇒ α2 = 96.45o
De manera que la posición que ocupa el cero sobre el eje real será:
180 − α2 = 83.55o
0.49
tan 83.55 = ⇒ a = 0.145
(a − 0.09)
d3 = 0.092 + 0.492 = 0.5
Problemas Resueltos 237
d4 = (0.09 + 0.5)2 + 0.492 = 0.77
d5 = (a − 0.09)2 + 0.492 = 0.49
Con el criterio del módulo, el valor global que debe tomar K en ese punto del
lugar será,
(z − 0.145)
R(z) = 0.56 ·
z
Problema 5.9
G(z)
Y(z)
X(z) +
1
R(z)
z(z-0.4)
-
3
ts = = 5.2 s ⇒ σ = 0.58 rad/s
σ
0.58
σ = ζ · ωn ⇒ ωn = = 1.35 rad/s
0.43
238 Problemas Resueltos
ωd = ωn 1 − ζ 2 = 1.35 1 − 0.432 = 1.22 rad/s
θ = T ωd = 1.22 rad
Pd
0.53
d1
d2
X X
0.4
-0.53
Problema 5.10
Determinar un regulador discreto R(z) de tipo PI para que los polos del
sistema en lazo cerrado ocupen las posiciones z1 = 0.8 y z2 = 0.05
R(z) G(z)
Y(z)
X(z) +
0.1
PI z-0.1
-
240 Problemas Resueltos
Para situar los polos en lazo cerrado en las posiciones indicadas, la ecuación
caracterı́stica del conjunto será:
(z − ci ) 0.1
p(z) = 1 + R(z) · G(z) = 1 + K · =0
(z − 1) (z − 0.1)
(z − 0.24)
R(z) = 2.5 ·
(z − 1)
Pd = 0.05 Pd = 0.8
X o X
0.1
ci
El valor de la K en cada uno de los polos en lazo cerrado del sistema será
el mismo. Aplicando el criterio del módulo en ambas posiciones e igualando se
obtiene:
K|0.05 = K|0.8
K = K|R(z) · K|G(z)
242 Problemas Resueltos
0.25
K|R(z) = = 2.5
0.1
Obteniéndose la misma función de transferencia del regulador:
(z − 0.24)
R(z) = 2.5 ·
(z − 1)
Problema 5.11
0.53
1
d2 d1
X X
0.8 0.9
-0.53
Aplicando el criterio del módulo desde este punto se obtiene el valor global
de K,
d1 = d 2 = (0.85 − 0.8)2 + 0.532 = 0.53
Sampling time: 1
Respuesta al escalón
1
0.9
0.8
0.7
0.6
Amplitud
0.5
0.4
0.3
0.2
0.1
0
0 2 4 6 8 10 12
Indice
0.9
0.8
0.7
0.6
Amplitud
0.5
0.4
0.3
0.2
0.1
0
0 2 4 6 8 10 12
Indice
246 Problemas Resueltos
R(z)
x(t) y(t) {yk}
17.45 (z-0.52) 1
B0(s)
T=0.2 z s (s+2) T=0.2
Respuesta al escalón
1.4
1.2
1
Amplitud
0.8
0.6
0.4
0.2
0
0 1 2 3 4 5
Indice
0.8
0.6
Amplitud
0.4
0.2
−0.2
−0.4
0 1 2 3 4 5
Indice
248 Problemas Resueltos
Ejercicio 5.3 Dado el sistema de la figura, comparar las respuestas que se ob-
tendrán si se implementa un regulador discreto R(z) como resultado de discretizar
el continuo R(s) = 2.84(s + 0.4)/s empleando bloqueador de orden 0, la aproxi-
mación de Tustin y la equivalencia polo-cero. La entrada es el escalón unitario:
G(s)
x(t) + y(t)
2
R(z) B0(s)
T=0.1 (s+0.5)(s+2)
-
1.2
0.4
0.2
0
0 1 2 3 4 5 6 7
Indice
249
Capı́tulo 6
Espacio de Estado.
Introducción y Análisis
6.1 Introducción
Hasta ahora, se han llevado a cabo una serie de hipótesis sobre los sistemas
estudiados que en ocasiones no son muy adecuadas con la realidad de los mismos.
En concreto, con objeto de simplificar el análisis y el diseño de los controladores, a
estos sistemas se les ha considerado lineales, monovariables - esto es, con una única
entrada y una salida - , invariantes en el tiempo y su funcionamiento consideraba
unas condiciones iniciales nulas. Sin embargo, por desgracia no son pocas las
ocasiones en los que los sistemas reales son más complejos, y de hecho suelen
presentar caracterı́sticas no lineales, además de ser multivariables (más de una
entrada y una salida). Para poder afrontar el estudio de estos sistemas más
acordes con la realidad surge la teorı́a del control en el espacio de estados. Esta
teorı́a permite el estudio de los sistemas con condiciones iniciales no nulas e incluso
de aquéllos que son variantes en el tiempo (sus parámetros no son constantes a
lo largo del tiempo).
La representación de los sistemas en el espacio de estado se traduce en un
conjunto de “n” ecuaciones en diferencias de primer orden, siendo “n” el orden
del sistema. A este tipo de representación se le conoce también con el nombre
de representación interna, frente a la representación externa de un sistema que
es su relación entrada-salida. La denominación de representación interna deriva
del concepto de estado que se define a continuación.
250 Capı́tulo 5: Espacio de Estado. Introducción y Análisis
Figura 6.1:
Y (z) b0 + b1 z −1 + · · · + bn z −n
=
U (z) 1 + a1 z −1 + · · · + an z −n
Y (z) b0 z n + b1 z n−1 + · · · + bn
= (6.8)
U (z) z n + a1 z n−1 + · · · + an
6.4 Representación en el Espacio de Estado 253
⎡ ⎤
x1 (k)
⎢
⎢ x2 (k) ⎥
⎥
y(k) = (bn − an b0 ) (bn−1 − an−1 b0 ) . . . (b1 − a1 b0 ) . ⎢
⎢ .. ⎥ + b0 u(k)
⎥
⎣ . ⎦
xn (k)
(6.10)
2. Forma canónica observable
La representación en el espacio de estado del sistema discreto dado por
las ecuaciones (6.7) y (6.8)
⎡ ⎤ ⎡ ⎤⎡ ⎤ ⎡ ⎤
x1 (k + 1) 0 0 ... 0 −an x1 (k) (bn − an b0 )
⎢ x2 (k + 1) ⎥ ⎢ 1 0 . . . 0 −an−1 ⎥⎢ x2 (k) ⎥ ⎢ (bn−1 − an−1 b0 ) ⎥
⎢ ⎥ ⎢ ⎥⎢ ⎥ ⎢ ⎥
⎢ .. ⎥ ⎢ .. .. .. .. .. ⎥⎢ .. ⎥ ⎢ .. ⎥
⎢ . ⎥=⎢ . . . . . ⎥⎢ . ⎥+⎢ . ⎥ u(k)
⎢ ⎥ ⎢ ⎥⎢ ⎥ ⎢ ⎥
⎣ xn−1 (k + 1) ⎦ ⎣ 0 0 ... 0 −a2 ⎦ ⎣ xn−1 (k) ⎦ ⎣ (b2 − a2 b0 ) ⎦
xn (k + 1) 0 0 ... 1 −a1 xn (k) (b1 − a1 b0 )
(6.11)
254 Capı́tulo 5: Espacio de Estado. Introducción y Análisis
⎡ ⎤
x1 (k)
⎢
⎢ x2 (k) ⎥
⎥
y(k) = 0 0 ... 1 .⎢
⎢ .. ⎥ + b0 u(k)
⎥
⎣ . ⎦
xn (k)
(6.12)
Y (z) c1 c2 cn
= + + ··· + + b0
U (z) z − p 1 z − p2 z − pn
La representación en el espacio de estado del sistema discreto dado por
las ecuaciones (6.7) y (6.8)
⎡ ⎤
x1 (k + 1) ⎡ ⎤⎡ ⎤ ⎡ ⎤
⎢ ⎥ p1 0 . . . 0 x1 (k) 1
⎢ x2 (k + 1) ⎥ ⎢ ⎥⎢ ⎥ ⎢ ⎥
⎢ .. ⎥ ⎢ 0 p2 . . . 0 ⎥⎢ x2 (k) ⎥ ⎢ 1 ⎥
⎢ ⎥=⎢ .. .. .. .. ⎥⎢ .. ⎥+⎢ .. ⎥ u(k)
⎢ . ⎥ ⎢ ⎥⎢ ⎥ ⎢ ⎥
⎢ ⎥ ⎣ . . . . ⎦⎣ . ⎦ ⎣ . ⎦
⎣ xn−1 (k + 1) ⎦
0 0 . . . pn xn (k) 1
xn (k + 1)
(6.13)
⎡ ⎤
x1 (k)
⎢⎢ x2 (k) ⎥
⎥
y(k) = c1 c2 . . . cn . ⎢
⎢ .. ⎥ + b0 u(k)
⎥
⎣ . ⎦
xn (k)
(6.14)
⎡ ⎤
x1 (k)
⎢⎢ x2 (k) ⎥
⎥
y(k) = c1 c2 . . . cn . ⎢
⎢ .. ⎥ + b0 u(k)
⎥
⎣ . ⎦
xn (k)
(6.16)
256 Capı́tulo 5: Espacio de Estado. Introducción y Análisis
P −1 GP = matriz diagonal
k−1
x(k) = Gk x(0) + Gk−j−1 Hu(j), k = 1, 2, 3, . . . (6.26)
j=0
k−1
y(k) = CGk x(0) + C Gk−j−1 Hu(j) + Du(k) (6.27)
j=0
k−1
Gk−j−1 Hu(j) = Z −1 [(zI − G)−1 HU (z)] (6.33)
j=0
6.7 Matriz de Transición de Estado 259
x(k + 1) = Gx(k)
(6.35)
como sigue:
x(k) = Ψ(k)x(0)
(6.36)
Ψ(k) = Gk (6.38)
k−1
x(k) = Ψ(k)x(0) + Ψ(k − j − 1)Hu(j) (6.39)
j=0
k−1
y(k) = CΨ(k)x(0) + C Ψ(k − j − 1)Hu(j) + Du(k) (6.40)
j=0
Y (z)
M (z) = = C(zI − G)−1 H + D
U (z)
M(z) se conoce como matriz de función de transferencia pulso que es una
matriz de dimensión m × r.
Al sustituir el valor de (zI − G)−1 ,
C · adj(zI − G)T H
M (z) = +D
|zI − G|
Se observa que los polos de M(z) se obtienen al resolver
p(z) = |zI − G| = 0
Ya que |zI − G| representa la ecuación caracterı́stica del sistema. Dicho
de otra manera, los autovalores o valores propios de la matriz G son los
polos del sistema.
262 Problemas Resueltos
Problemas Resueltos
Problemas Resueltos 263
Problema 6.1
a2 = 0.4, a1 = 1.3
b2 = 1, b1 = 1, b0 = 0
La representación mediante la forma canónica controlable es:
x1 (k + 1) 0 1 x1 (k) 0
= + u(k)
x2 (k + 1) −0.4 −1.3 x2 (k) 1
x1 (k)
y(k) = 1 1)
x2 (k)
x1 (k + 1) = x2 (k)
x2 (k + 1) = −0.4x1 (k) − 1.3x2 (k) + u(k)
y(k) = x1 (k) + x2 (k)
264 Problemas Resueltos
x1 (k + 1) −0.8 0 x1 (k) 1
= + u(k)
x2 (k + 1) 0 −0.5 x2 (k) 1
x1 (k)
y(k) = −0.67 1.67
x2 (k)
k = 0
x1 (1) = x2 (0) = 0
x2 (1) = −0.4x1 (0) − 1.3x2 (0) + u(0) = 0 + 0 + 1 = 1
y(0) = x1 (0) + x2 (0) = 0
k = 1
x1 (2) = x2 (1) = 1
x2 (2) = −0.4x1 (1) − 1.3x2 (1) + u(1) = 0 − 1.3 + 1 = −0.3
y(1) = x1 (1) + x2 (1) = 0 + 1 = 1
k = 2
x1 (3) = x2 (2) = −0.3
x2 (3) = −0.4x1 (2) − 1.3x2 (2) + u(2) = −0.4 + (1.3 × 0.3) + 1 = 0.99
y(2) = x1 (2) + x2 (2) = 1 − 0.3 = 0.7
Problemas Resueltos 265
k = 3
x1 (4) = x2 (3) = 0.99
x2 (4) = −0.4x1 (3) − 1.3x2 (3) + u(3) = (−0.4 ∗ −0.3) − (1.3 × 0.99) + 1
= −0.17
y(3) = x1 (3) + x2 (3) = −0.3 + 0.99 = 0.69
k = 4
x1 (5) = x2 (4) = −0.17
x2 (5) = −0.4x1 (4) − 1.3x2 (4) + u(4) = (−0.4 ∗ 0.99) + (1.3 × 0.17) + 1
= 0.825
y(4) = x1 (4) + x2 (4) = 0.99 − 0.17 = 0.82
z 0 0 1 z −1
[zI − G] = − =
0 z −0.4 −1.3 0.4 z + 1.3
z 0.4
adj
adj [zI − G] T
−1 z + 1.3
[zI − G]−1 = =
|zI − G| z 2 + 1.3z + 0.4
1 z + 1.3 1
= 2
z + 1.3z + 0.4 −0.4 z
Luego,
1 1 z
= 2
z + 1.3z + 0.4 z z−1
⎡ ⎤
z
(z 2 +1.3z+0.4)(z−1)
= ⎣ z2
⎦
(z 2 +1.3z+0.4)(z−1)
z(z + 1)
Y (z) = CX(z) = 1 1 X(z) =
(z 2 + 1.3z + 0.4)(z − 1)
Problema 6.2
Se pide:
a2 = 3, a1 = −4
b2 = −1, b1 = 2, b0 = 0
Problemas Resueltos 267
x1 (k + 1) 0 1 x1 (k) 0
= + u(k)
x2 (k + 1) −3 4 x2 (k) 1
x1 (k)
y(k) = −1 2)
x2 (k)
x1 (k + 1) = x2 (k)
x2 (k + 1) = −3x1 (k) + 4x2 (k) + u(k)
y(k) = −x1 (k) + 2x2 (k)
x1 (k + 1) 3 0 x1 (k) 1
= + u(k)
x2 (k + 1) 0 1 x2 (k) 1
x1 (k)
y(k) = 2.5 −0.5)
x2 (k)
x1 (k + 1) 0 −3 x1 (k) −1
= + u(k)
x2 (k + 1) 1 4 x2 (k) 2
x1 (k)
y(k) = 0 1)
x2 (k)
Problema 6.3
Figura 6.2:
Se pide:
Y (z) 0.5
= 2
U (z) z − 0.5z
a2 = 0, a1 = −0.5
b2 = 0.5, b1 = 0, b0 = 0
x1 (k + 1) 0 1 x1 (k) 0
= + u(k)
x2 (k + 1) 0 0.5 x2 (k) 1
x1 (k)
y(k) = 0.5 0)
x2 (k)
x1 (k + 1) = x2 (k)
x2 (k + 1) = 0.5x2 (k) + u(k)
y(k) = 0.5x1 (k)
Y (z) −1 1
= +
U (z) z z − 0.5
polos : p1 = 0, p2 = 0.5
coeficientes : c1 = −1, c2 = 1, b0 = 0
x1 (k + 1) 0 0 x1 (k) 1
= + u(k)
x2 (k + 1) 0 0.5 x2 (k) 1
270 Problemas Resueltos
x1 (k)
y(k) = −1 1)
x2 (k)
x1 (k + 1) = u(k)
x2 (k + 1) = 0.5x2 (k) + u(k)
y(k) = −x1 (k) + x2 (k)
Problema 6.4
x1 (k + 1) 0 1 x1 (k) 2
= + u(k)
x2 (k + 1) −0.21 −1 x2 (k) −2
x1 (k)
y(k) = 1 0)
x2 (k)
−1
−1 z −1
(zI − G) =
0.21 z + 1
z+1 1
(z+0.3)(z+0.7) (z+0.3)(z+0.7)
= −0.21 z
(z+0.3)(z+0.7) (z+0.3)(z+0.7)
1.75 0.75 2.5 2.5
(z+0.3)
− (z+0.7) (z+0.3)
− (z+0.7)
= −0.525 0.525 −0.75 1.75
(z+0.3)
+ (z+0.7) (z+0.3)
+ (z+0.7)
2z
Sabiendo que U (z) = z−1
2z 2z 2
2z
zx(0) + HU (z) = + z−1
2z = z−1
−2z 2 +4z
−2z z−1 z−1
⎡ ⎤
(2z 2 +4)z
(z+0.3)(z+0.7)(z−1) ⎦
= ⎣ (−2z 2 +3.58z)z
(z+0.3)(z+0.7)(z−1)
−8.04z 7.32z 2.71z
(z+0.3)
+ (z+0.7) + (z−1)
= 2.41z 5.13z 0.71z
(z+0.3)
− (z+0.7) + (z−1)
−1 −8.04(−0.3)k + 7.32(−0.7)k + 2.71
x(k) = Z [X(z)] =
2.41(−0.3)k − 5.13(−0.7)k + 0.71
Problema 6.5
x1 (k + 1) 0 1 x1 (k) 0
= + u(k)
x2 (k + 1) −0.4 −1.3 x2 (k) 1
x1 (k)
y(k) = 1 0)
x2 (k)
2 0
P =
−1 2
(k + 1)
x = P −1 GP x
(k) + P −1 Hu(k)
2 0 0 1 2 0 1 (k)
x 2 0 0
= + u(k)
1 2 −0.4 −1.3 −1 2 2 (k)
x 1 2 1
−2 4 1 (k)
x 0
= + u(k)
0 −3.2 2 (k)
x 2
ẋ = Ax + Bu
x˙1 0 1 x1 0
= + u
x˙2 −0.45 −1.4 x2 1
y = Cx + Du
x
1
y= 1 0
x2
Obtener:
G =
0.9553 0.3529
-0.1588 0.4612
H =
0.0994
0.3529
H =
0.1814
0.7256
Cd =
1 0
Dd =
0
Respuesta al escalón
2.5
1.5
Amplitud
0.5
0
0 5 10 15 20 25
Indice (sec)
x1 (k + 1) 0.9553 0.3529 x1 (k) 0.0994
= + u
x1 (k + 1) −0.1588 0.4612 x2 (k) 0.3529
x1 (k)
y(k) = 1 0)
x2 (k)
Obtener:
1. La representación en el espacio de estado del sistema continuo equiv-
alente para un tiempo de muestreo de T=0.5 segundos.
B =
0
1.0000
C =
1 0
D =
0
Respuesta al escalón
2.5
1.5
Amplitud
0.5
0
0 2 4 6 8 10 12
Tiempo (sec)
x1 (k + 1) 0 1 x1 (k) 0
= + u(k)
x2 (k + 1) −0.4 −1.3 x2 (k) 1
x1 (k)
y(k) = 1 1)
x2 (k)
x1 (k + 1) 0 1 x1 (k) 0
= + u(k)
x2 (k + 1) −0.4 −1.3 x2 (k) 1
x1 (k)
y(k) = 1 1)
x2 (k)
x=
0 0
0 1.0000
1.0000 -0.3000
-0.3000 0.9900
0.9900 -0.1670
-0.1670 0.8211
x1 (k + 1) 0 1 x1 (k) 0
= + u(k)
x2 (k + 1) −3 4 x2 (k) 1
x1 (k)
y(k) = −1 2)
x2 (k)
Ejercicios con MATLAB 281
x=
0 0
0 1
1 5
5 18
18 58
58 179
x1 (k + 1) 3 0 x1 (k) 1
= + u(k)
x2 (k + 1) 0 1 x2 (k) 1
x1 (k)
y(k) = 2.5 −0.5)
x2 (k)
282 Problemas Resueltos
x=
0 0
1 1
4 2
13 3
40 4
121 5
Ejercicio 6.7 Convertir el siguiente sistema discreto representado por el
espacio de estado a una función de transferencia:
x1 (k + 1) 0.9553 0.3529 x1 (k) 0.0994
= + u
x1 (k + 1) −0.1588 0.4612 x2 (k) 0.3529
x1 (k)
y(k) = 1 0)
x2 (k)
Ejercicios con MATLAB 283
0 0.0994 0.0787
den =
0.0994z+0.0787
------------------------
z∧ 2- 1.4164 z + 0.4966
Sampling time: 1
Ejercicio 6.8 Convertir el siguiente sistema discreto dado por una función
de transferencia pulso a la representación en el espacio de estado:
Transfer function :
1
------------------------
z∧ 2+ z + 0.21
Sampling time: 1
284 Problemas Resueltos
-1.0000 -0.2100
1.0000 0
H =
1
0
C =
0 1
D =
0
x1 (k + 1) 0.9553 0.3529 x1 (k) 0.0994
= + u(k)
x2 (k + 1) −0.1588 0.4612 x2 (k) 0.3529
x1 (k)
y(k) = 1 0)
x2 (k)
Ejercicios con MATLAB 285
p =
0.7788
0.6376
k =
0.0994
>> fdt=zpk([-0.7919],[0.7788;0.6376],[0.0994],1)
Zero/pole/gain:
0.0994(z+0.7919)
------------------------
(z-0.7788)(z-0.6376)
Sampling time: 1
Índice alfabético
287
288 Índice alfabético