Tema 3.

Ejercicios adicionales (completos)

Sistemas Electrónicos de Control - 0910a

Tema 3: Ejercicios adicionales (completos)

Para recuperar entregas o para subir nota

Escoger sólo 1 de cada apartado
Fecha límite de entrega: Día del Examen Final

Contenido
1. Especificaciones ....................................................................................................................... 1 
2. Métodos empíricos ................................................................................................................... 4 
3. Métodos gráficos ...................................................................................................................... 6 
4. Métodos de estado .................................................................................................................. 12 
5. Control robusto....................................................................................................................... 15 
6. Dos grados de libertad ............................................................................................................ 18 
7. Control óptimo ....................................................................................................................... 26 

1. Especificaciones

Ejercicio 1. 1. Construcción de M(s). Guías. Se trata de construir/sintetizar una función

de transferencia M ( s) cuya respuesta indicial presente aproximadamente las siguientes
características (especificaciones):
(a) De estabilidad: R pt  5% .
(b) De precisión: error permanente nulo a una entrada en escalón y, además, una k v  100
(c) De velocidad de respuesta: t r  0.1 .
(d) De realizabilidad: exceso de polos sobre ceros igual a 2.
Se pide: Construir/sintetizar dicha M ( s) a partir del correspondiente bloque dominante de
segundo orden. (Recuérdese que los ceros también influyen en la forma del transitorio a través
de los residuos.)

Ejercicio 1. 2. Especificaciones. Si Sr = 1 (caso del sistema óptimo LQR)), estimar MG y

MF.

Ejercicio 1. 3. Loopshaping. Traducción de las especificaciones del servo a la forma de

L(j). En el marco del control, el enfoque directo (selección de polos y ceros de M(s)) va
dirigido a construir la función deseable del servo M(s) para que este satisfaga unas determinadas
especificaciones, en general temporales. Aquí vamos a cambiar el enfoque en dos aspectos: 1)
En vez de M(s) construiremos L(s), 2) En vez de hablar de polos y ceros lo haremos de
ganancias del lazo L ( j ) . Este es el enfoque clásico de Bode y Horowitz que hoy, con el
renovado interés por los métodos frecuenciales, retorna con fuerza en los sistemas MIMO y en

ETSETB. Sistemas Electrónicos de Control 0910a 1

 Tema 3. Ejercicios adicionales (completos)

control robusto a través del análisis vía valores singulares. En concreto, se trata de determinar
las características del lazo L ( j ) para que un determinado servo satisfaga las siguientes
especificaciones:
(a) Precisión: - error permanente nulo a una entrada en escalón.
- error permanente inferior a 0.002 a una entrada en rampa.
- atenuación máx. de 0.02 para las frecuencias situadas en la banda de 0 a 1rad.
(b) Velocidad:  b  2 .
(c) Estabilidad: MF positivo.
(d) Rechazo al ruido: Atenuación mínima de 0.01 para las frecuencias superiores a 1000 rad/s
(e) Planta: Presenta un exceso de P  Z  2 .
Estudio Previo:
1) Dibujar en el diagrama de Bode las zonas de bajas y altas frecuencias por donde no puede
pasar L ( j ) .
2) Dibujar posibles asíntotas de bajas frecuencias (precisión estática), altas frecuencias (filtrado
del ruido de los sensores) e intermedias (recordando que la fase asociada no debe alcanzar los -
180).
3) Unir estas asíntotas mediante otras asíntotas de pendientes adecuadas.
4) Formular la expresión de L(s) así determinada.
5) Comprobar k p , kv , MF (dibujar la curva de fase) y  b (con ayuda del ábaco de Nichols).

Simulación:
1) Hallar M ( j ) y comprobar  b , la atenuación en la gama de 0 a 1 y de   1000 .
2) Hallar la respuesta indicial y a la rampa y comprobar la precisión del resultado.

Ejercicio 1. 4. Sistema de orden n. Con los mismos datos que el ¡Error! No se encuentra el
origen de la referencia., se pide:
1) Sin tener en cuenta la planta, construir la M(s) por los siguientes métodos:
1.1) Usando una forma prototipo de orden n (polinomios ITAE) de las tablas.
1.2) Minimizando el ISE. En este caso, al no dar tablas, se trata de obtener numéricamente
1
(MATLAB) el valor óptimo de los coeficientes de M 2 ( s)  2 ,
s  a1 s  1
1
M 3 ( s)  3 . Comprobarlo teóricamente (Parseval).
s  a 2 s  a1 s  1
2

2) Teniendo en cuenta la planta, construir una M(s) óptima (criterio ISE con penalización de
10
control) con ayuda del lugar cuadrático/simétrico de las raíces suponiendo G p ( s )  ,
s ( s  10)
20
  1 y Gp ( s)  ,   0.1 .
( s  2 )( s  6)
3) Comprobar por simulación los resultados anteriores y, guiados interactivamente por los
resultados, acabar de ajustar los parámetros de cada uno de ellos para alcanzar las
especificaciones exigidas. 

Ejercicio 1. 5. Polos y ceros adicionales. Se desea ajustar aproximadamente una M(s) para que
presente una respuesta caracterizada por las siguientes especificaciones: (a) R pt  15% , (b)
 b  40 , (c) M ( 0)  1 y (d) k v  60 .
Se pide:
k
1) Tratar de lograrlo con la forma M ( s)  2 .
s  2 n s   2n
(Nota: Obsérvese la dificultad de ajustar las especificaciones con sólo tres parámetros)

ETSETB. Sistemas Electrónicos de Control 0910a 2

 Tema 3. Ejercicios adicionales (completos)

2) En caso de que k v resulte insuficiente, mejorarla

2.1) Añadiendo un prefiltro con un cero.
2.2) Añadiendo un prefiltro con un dipolo.
3) Comparar ambas síntesis/formulaciones con la de segundo orden estricta. 

Ejercicio 1. 6. Waterbed effect. Si un servo presenta un lazo con un exceso de polos sobre
ceros (P-Z  2) demostrar que S   1 . (Nota: Demostrar previamente que

0
ln 1  L ( j ) d  0 a partir de los teoremas de Bode). 

ETSETB. Sistemas Electrónicos de Control 0910a 3

 Tema 3. Ejercicios adicionales (completos)

2. Métodos empíricos

Ejercicio 2. 1. Compensación anti-reset windup. Se trata de comprobar el efecto de la

saturación de la válvula asociada a un regulador con acción integral (reset windup) y su
compensación con un lazo auxiliar no lineal (anti-reset windup).
Datos:
KP  1
KI  1

Figura 1

0) Compensación antireset-windup. Reflexionar sobre el problema de la saturación de la válvula

(con PI) y trazar la evolución de uc ( t ) .
1) Compensación antireset-windup. Hallar la respuesta indicial c ( t ) a un escalón de r ( t ) y
representar, además, la evolución de las señales de control u ( t ) y uc ( t ) . Seguir los siguientes
pasos:
1.1) Sistema simple sin saturación ni compensación: Tomar KI  4 . Hacer s  10 y A  0 (esto
equivale a eliminar la aparición de la saturación y de la necesidad del compensador).
1.2) Introducir el efecto de saturación tomando s  11. . Observar el comportamiento de la salida
y las señales u y uc .
1.3) Añadir el control con zona muerta poniendo A  1 y s  11 . . Observar la salida y las señales
u y uc . ¿Qué diferencias hay respecto al caso anterior?

Ejercicio 2. 2. Control de procesos con retardos. El predictor de Smith. Se trata de

estudiar el efecto de los retardos puros, asociados al transporte de materia o al carácter
distribuido del sistema, sobre la estabilidad del sistema y como se pueden extrapolar las técnicas
de configuración clásicas (sin retardos) al caso en cuestión
Datos:
1
Planta: GP ( s)  y retardo puro: T0  2
s 1
 1 
Regulador PI: Gc ( s)  0.4 1  
 0.4 s 
(Calculado para la planta hipotética GP ( s) )

ETSETB. Sistemas Electrónicos de Control 0910a 4

 Tema 3. Ejercicios adicionales (completos)

Figura 2

Se pide estimar la respuesta a una señal de mando r ( t )  u ( t ) y a una perturbación

L ( t )  u ( t  20 ) , en el supuesto de:
1) T0  1 y sin predictor
2) T0  2 y sin predictor
3) T0  2 con predictor
(Nota: Aquí u( t ) simboliza la función escalón unitario)
4) Hallar las respuestas especificadas

Ejercicio 2. 3. Efecto de los algoritmos P, PI, PD y PID sobre el transitorio y permanente

de un servo. Se trata de elegir un Gc ( s ) adecuado para insertarlo en serie con la planta
4500
G ( s)  según el esquema de la Figura 3.
s( s  300)

Figura 3
Para ello se pretende estudiar sus efectos sobre las características de la respuesta valiéndonos
del ordenador. 
Simulación:
Estudiar su comportamiento con los siguientes algoritmos de control:
1) Control P con kc  182 . Determinar tp, Rpt y el error permanente a una entrada en rampa
er().
2) Control PD (para estabilizarlo más) con Gc ( s )  182 (1  0. 002 s) . Determinar tp, Rpt y er().
 0.11
3) Control PI (para reducir a cero el error de velocidad) con Gc ( s)  182 0.075  .
 s 
Determinar tp, Rpt y er().
 4.5 
4) Control PID (para mejorar el transitorio anterior) con Gc ( s)  182 0.31   0.0006s .
 s 
Determinar tp, Rpt y er().
5) Comprobar que el efecto del polo adicional (realizabilidad) de los reguladores prácticos no
afectan esencialmente el transitorio si n10. 

ETSETB. Sistemas Electrónicos de Control 0910a 5

 Tema 3. Ejercicios adicionales (completos)

3. Métodos gráficos

Ejercicio 3. 1. Compensación tacométrica (en paralelo). Lazo auxiliar. Se trata de

calcular los parámetros del compensador tacométrico (compensación en paralelo o control con
lazo auxiliar) de la Figura 4.
Las especificaciones son k v  10 y m  35 . La corrección se realiza en paralelo con un
tacómetro y un filtro quasi derivador que evita las altas frecuencias/ruido de forma que la
Ts
transmitancia total del compensador es H ( s)  s kT .
1  Ts

Figura 4 Figura 5

1) Fijar k p para satisfacer la precisión.

2) Obtener, con la ayuda de un diagrama de Bode, unos valores adecuados para k T y T ( RC ) .
Simulación:
1) Representar en un diagrama de Bode la respuesta frecuencial del lazo principal.
2) Representar la respuesta indicial del sistema.

Ejercicio 3. 2. Servomecanismo de posición con lazo interior tacométrico. El modelo

de un determinado servomecanismo es el siguiente:

G1(s)

H2(s)

H3(s)

Figura 6
donde
k A  9  10 4 k 1
( s  5)( s  30)
G1 ( s) 
( s  20)( s  285)
s( s  50)( s  26)( s 2  15s  10 6 )
G p ( s)  12
( s  300)( s 2  176s  300 2 )
( s  26)
H 2 ( s) 
( s  01. )( s  20)( s  285)
k T  12 k 2
0.07
H 3 ( s) 
s2

ETSETB. Sistemas Electrónicos de Control 0910a 6

 Tema 3. Ejercicios adicionales (completos)

Se trata de estudiar su comportamiento y elegir adecuadamente k T y k A para un correcto

funcionamiento. Para ello se analiza y ajusta el lazo interior (vía Bode, Nichols, Evans y
Nyquist) y luego el lazo exterior, comparando los efectos de los distintos diseños. Se pide:
1) Lazo interior:
1.1) Trazar el lugar de las raíces correspondiente a la evolución de los polos del lazo interior al
variar k2. Determinar el valor límite y la frecuencia de oscilación.
1.2) Trazar los diagramas de Bode del lazo interior L1(j) con k 2  1 y estudiar la estabilidad.
1.3) Trazar el diagrama polar de L1(j) con k 2  1 y comprobar la estabilidad con el criterio de
Nyquist.
1.4) Hallar, con ayuda del ábaco de Nichols, la respuesta en frecuencia del lazo menor
G p ( s)
M 1 ( s)  y representarla en coordenadas de Bode (log-log) con k 2  0.1 .
1  k T G p ( s) H 1 ( s)
2) Lazo externo o principal:
2.1) Trazar el diagrama de Bode del lazo principal (usando las curvas de M1(j) halladas) para
los casos k 1  1 y k 1  0.1 . Estudiar la estabilidad.
2.2) Con ayuda de Nichols hallar el diagrama de Bode de M(j) para k 1  0.1 .
Simulación:
1) Lazo interno:
1.1) Representar el lugar geométrico de las raíces de Evans (para k2 variable). Hallar k1 , crítico
y .
1.2) Representar el diagrama de Bode para k 2  1 .
1.3) Representar el diagrama en polares para k 2  1 .
1.4) Representar el diagrama de Bode de M1(j) con k 2  0.1 .
1.5) Hallar y representar los polos y ceros de M1(s).
2) Lazo externo:
2.1) Representar el lugar geométrico de las raíces de Evans (para k1 variable). Hallar k1, crítico y
.
2.2) Trazar el diagrama de Bode del lazo principal (para k1 = 0.05; 0.063; 0.1; 0.15). Para
k1=0.15 pasarlo a A- (amplitud-fase) y determinar r y el valor de la resonancia.
2.3) Trazar el diagrama de Bode de M(j) principal (para k1 = 0.05; 0.063; 0.1; 0.15).
2.4) SIMULINK: Representar la respuesta indicial para los valores de k1 indicados. 

Ejercicio 3. 3. Pseudoderivada. Debido a que cualquier circuito derivador (o su equivalente

numérico) es muy sensible al ruido (inevitable), en la práctica no se usan los derivadores puros
Ts
Gc ( s)  Td s sino una versión más realista Gc ( s )  d llamada pseudoderivada. Se trata de
1  Tp s
valorar las características de esta aproximación comparando los circuitos de las Figura 7 y
Figura 8.

Figura 7 Figura 8

Estudio Previo:

ETSETB. Sistemas Electrónicos de Control 0910a 7

 Tema 3. Ejercicios adicionales (completos)

V0 ( s)
1) Calcular la función de transferencia G ( s)  del circuito de la Figura 7, comprobar que
V1 ( s)
es un derivador y representar su diagrama de Bode de ganancia.
V ( s)
2) Repetir el cálculo de Ga ( s)  0 para el circuito de la Figura Figura 8, comprobar que es
V1 ( s)
un pseudoderivador y representar su diagrama de Bode de ganancia.
3) Comentar en qué condiciones podemos suponer que en primera aproximación, Ga(s) se
comporta como un derivador y en qué condiciones no se ve afectado por el ruido.
Simulación:
1) En el caso derivador fijar los valores R  C  1 y excitar el circuito con las siguientes señales:
escalón, rampa, sinusoide y ruido blanco. Comentar los resultados observados.
2) En el caso pseudoderivador, elegir unos valores de R  C  R p  1 y excitarlo con las
siguientes señales: escalón, rampa, senoide y ruido blanco. Comentar los resultados. 

En los próximos Ejercicios se pretende comprobar la mejora en la dinámica de la planta

Y ( s) 1
G P ( s)   al ser gobernada, en lazo cerrado, mediante un
U ( s) ( s  1)( s  4 )( s  16)
controlador Gc ( s ) en serie (ver Figura 9).
P

Figura 9

Ejercicio 3. 4. Bode. Resolver el mismo problema de los Ejercicios anteriores, pero esta vez
con la ayuda de los diagramas de Bode y mediante un controlador serie de avance
1  ns
Gc ( s)  kc , n 1, de manera que se satisfagan las siguientes especificaciones:
1  s
- Precisión: kv  10
- Estabilidad relativa: R pt  25% 

Estudio Previo:
1) Si Gc  kc , calcular k c para que el servo tenga una k p  20 .
2) Si Gc  kc , con ayuda de los diagramas de Bode de la función de transferencia del lazo
L ( j ) , determinar kcu y Tu .
3) Convertir la especificación R pt  25% en otra equivalente en forma de m .
4) Calcular, por el método desarrollado en clase, los parámetros kc , n y  .
Simulación:
1) Mediante compensación estática o ajuste de la ganancia ( Gc  kc ), determinar  c , mG y m
correspondientes al valor de k c calculado en el apartado 1) del Estudio Previo.
2) Mediante compensación por avance  Gc ( s)  kc 1  ns ; n  1 , hallar los nuevos valores de
 1  s 
 c , mG y m si ahora, además de ajustar k c al valor anterior, añadimos un filtro corrector con
1  0. 2 s
transmitancia .
1  0. 02 s

ETSETB. Sistemas Electrónicos de Control 0910a 8

 Tema 3. Ejercicios adicionales (completos)

3) Mediante compensación por retardo  Gc ( s)  k c 1  s ; n  1 hallar los nuevos valores de

 1  ns 
 c , mG y m si ahora, además de ajustar k c al valor del apartado 1) añadimos un filtro
1  2s
corrector con función de transferencia .
1  20 s
4) Mediante compensación por retardo y avance  Gc ( s)  kc 1  1s 1  n22 s ; n1 , n2  1 , hallar los
 1  n11s 1  2 s 
nuevos valores de  c , mG y m si ahora, además de ajustar k c , al valor del apartado 1),
1  2 s 1  0. 2 s
añadimos un filtro corrector con función de transferencia .
1  20 s 1  0. 02 s
5) Comparar las respuestas a perturbaciones en escalón de los sistemas de los apartados 1), 2),
3) y 4).
6) Tanteando, hallar un nuevo corrector de retardo que consiga el mismo m que el de avance
del apartado 2). Hallar la respuesta indicial (de la consigna) y comprobar su R pt . (En caso de
no satisfacer los requisitos, retocar los parámetros)
7) Tanteando, hallar un nuevo corrector de retardo-avance que alcanzando el mismo m que el
de avance (apartado 2)) retenga la  c obtenida con el sólo ajuste de la ganancia (apartado
1)). 

Ejercicio 3. 5. Ábaco de Nichols. Resolver el mismo problema del ¡Error! No se encuentra el

origen de la referencia., pero esta vez con la ayuda del ábaco de Nichols.

Estudio Previo:
1) Estimar, por el método desarrollado en clase, los parámetros k c , n y  (¿Qué inconveniente
tiene este método con relación al de Bode?)
2) Dibujar en el plano fase-ganancia la respuesta frecuencial de L ( s ) con la k c calculada en el
apartado 1) del estudio previo del ¡Error! No se encuentra el origen de la referencia..
3) Estimar  b y tr .
4) En el caso de que pudiéramos disponer de un corrector ideal de avance (es decir, sin que
modificase la ganancia), determinar con ayuda del ábaco de Nichols el avance de fase
necesario para producir una resonancia de Mr  2 dB .
Simulación:
1) Dibujar en el plano fase-ganancia la respuesta frecuencial del lazo L ( j ) , antes y después
de la compensación, determinando la variación de m , mG ,  c ,  b ,  r y M r . Estimar la n
del servo.
2) Hallar la respuesta indicial y determinar R pt y tr . Hallar el offset, t p y t s . Comparar el valor
de t r con el previsto por la expresión  btr  3 .
3) Comprobar el producto  btr  3 . 

Ejercicio 3. 6. Evans. Se trata de resolver el mismo problema del ¡Error! No se encuentra el

origen de la referencia., manteniendo la especificación de estabilidad pero cambiando la de
precisión por la de velocidad (concretamente  n  11 rad / s ) y utilizando el lugar de las raíces
de Evans. 
Estudio Previo:
1) Convertir la especificación de R pt  25% en el correspondiente coeficiente de
amortiguamiento  y, junto con n , fijar la situación de los polos dominantes. Estimar el
valor previsible de m ( m  100  ) .

ETSETB. Sistemas Electrónicos de Control 0910a 9

 Tema 3. Ejercicios adicionales (completos)

1  ns sz
2) Reformular el controlador Gc ( s)  kc en la forma Gc ( s)  kc' (¿Por qué?).
1  s s p
3) Trazar aproximadamente el lugar de las raíces de 1  L ( s ) correspondiente a la variación de
Gc  kc .
4) Determinar gráficamente la situación de un cero adicional que hiciera m  45 . Calcular, por
el método desarrollado en clase, k c' , z y p y, a partir de ellos, k c , n y  .
5) Calcular la k v que se obtiene de esta manera.
6) Fijando k c para el valor calculado en los ejercicios anteriores, determinar los polos
dominantes (comprobar si lo son), su coeficiente de amortiguamiento  y hacer una
estimación de t r .
Simulación:
1) Comprobar el lugar geométrico del estudio previo y determinar las raíces para la k c indicada.
1  0. 2 s
2) Trazar el lugar de Evans cuando al lazo se le inserta el filtro , observar el cambio que
1  0. 02 s
ha experimentado y determinar las nuevas raíces.
1  2s
3) Ídem para el filtro .
1  20 s
1  2 s 1  0. 2 s
4) Ídem para el filtro .
1  20 s 1  0. 02 s
5) Hallar las respuestas indiciales (de consigna) y comprobar R pt y  n (a partir de  d ). 

Ejercicio 3. 7. Compensación en paralelo. Nyquist inverso. Aquí, con la misma planta que
Y ( s) 1
venimos estudiando, G P ( s)   , en vez de la compensación en
U ( s) ( s  1)( s  4 )( s  16)
serie, se supone la configuración de la Figura 10 que consiste en una compensación tacométrica
(en paralelo) además del amplificador k c (en serie).

Figura 10 
Estudio Previo:
1) Calcular M ( s ) . A partir de ella indicar el efecto que cada uno de los parámetros k c y kT
tiene sobre los polos y ceros de M ( s ) . Fijar k c para alcanzar la precisión exigida.
2) Tomando k c  1280 y, con ayuda del lugar de Evans, determinar el valor de kT para obtener
  0. 7 ( R pt  25% ). (Nota: Para resolver este problema hay que hacer dos lugares de raíces:
N ( s)
uno, previo, para determinar las raíces de D(s) en la expresión 1  k T  0 , y el segundo,
D ( s)
de diseño, para fijar el valor de kT .).
3) Comprobar si los polos complejos del servo son o no dominantes.
kcGp
4) Con k c  1280 , y con ayuda del diagrama polar inverso de M1 ( s)  y del ábaco
1  kcGp kT s
de Hall inverso, determinar el valor de kT para   0. 7 . (Dibujar en polares la respuesta
frecuencial de la transmitancia inversa del lazo interno ( M 11 ( s )) y, comparándola con la
correspondiente circunferencia de Mr , determinar la kT necesaria).

ETSETB. Sistemas Electrónicos de Control 0910a 10

 Tema 3. Ejercicios adicionales (completos)

Simulación:
1) Calcular las raíces del polinomio D(s) del estudio previo y compararlos con el resultado allí
obtenido.
2) Trazar el lugar de los polos del servo M ( s ) con k c  1280 y kT variable.
3) Trazar el diagrama polar inverso de la transmitancia M1 ( s ) del lazo interno ( k c  1280 y kT
variable) y, por tanteo, variando kT , determinar el valor de kT para que M r  1. 2 . (Dibujar
en polares M 11 ( j ) y determinar kT añadiendo verticalmente kT  ( kT  0) hasta conseguir
la tangencia en la correspondiente circunferencia del ábaco de Hall). Determinar  r y
comprobar que k T  r es el desplazamiento que ha experimentado la curva sin compensar en
la  r 0 .
4) Comparar ambos valores y hallar las respuestas indiciales. Verificar R pt .
5) Dibujar la respuesta a una entrada en escalón sin el tacómetro y con el tacómetro. 

ETSETB. Sistemas Electrónicos de Control 0910a 11

 Tema 3. Ejercicios adicionales (completos)

4. Métodos de estado

Ejercicio 4. 1. Realizaciones de una H(s). Formas canónicas. Dada la función de

transferencia

s2 s3
H ( s)  H ( s) 
s  4s  3
2
s  2s  4
2

Se pide:
1) Obtener su realización en Forma Companion Controlable (FCC) (versión Matlab). (Matlab:
tf2ss)
2) Obtener su realización en Forma Companion Controlable (FCC) (versión teoría). (Matlab:
rot90, flipud, fliplr)
3) Obtener su realización en Forma Companion Controlable (FCC) a partir de la trasformación
de similaridad x  Tx , T  M c W . (Matlab: ctrb, poly, inv)
4) Obtener su realización en Forma Companion Observable (FCO) a partir de la trasformación
de similaridad x  Tx , T  ( WMo ) 1 . (Matlab: obsv, poly, inv)

Ejercicio 4. 2. Solución de ecuaciones de estado. ZS, ZI y matriz de transición. Dada

la función de transferencia

s2 s3
H ( s)  H ( s) 
s  4s  3
2
s  2s  4
2

se trata de obtener su respuesta temporal.

Se pide:
1) Calcular la respuesta zero-input tomando como condiciones iniciales (CI) a x1 ( 0)  1 y
x2 ( 0)  1. (Matlab: tf2ss, initial, ss, plot, title, xlabel)
2) Calcular la respuesta zero-state a una excitación escalón unitario. (Matlab: step, ss, plot,
title, xlabel)
3) Calcular y representar la respuesta total (zero-state más zero-input). (Matlab: lsim, ss, plot
title, xlabel)
4) Presentar en el plano de estado ( x2 ( x1 ) ) la respuesta del sistema a un conjunto de condiciones
iniciales (variando entre –1 y 1). (Matlab: for..end, initial, ss, plot, axis, hold on hold off, title,
xlabel, ylabel)

Ejercicio 4. 3. Formas canónicas. Dado el sistema con:

5  2  4   1  1
     3  1  2
A  2 1  2  B3 C
0 
; 1 ;
1  1 0   1  1  1 0

Se pide:

ETSETB. Sistemas Electrónicos de Control 0910a 12

 Tema 3. Ejercicios adicionales (completos)

1) Convertir las matrices A, B y C a la Forma Canónica Diagonal (Nota: Calcularla a

través de x  Mz).
2) Convertir las matrices A, B y C a la Forma Canónica Diagonal (Nota: Calcularla a
través de la función canon).
3) Obtener la Forma Canónica Controlable de Kalman (FCCK), separando el
subsistema controlable ( A c ) del no controlable ( A c ).
A c A12 

Nota: La transformación es x  Tx , T  T1 T2  A     , siendo T1 una
 0 Ac 
submatriz con n1  n columnas linealmente independientes de M c y el resto (T2)
columnas arbitrarias que no hagan singular la matriz T. La función ctrbf  de
 Ac 0
matlab da como resultado A   .
A 21 A c 

Ejercicio 4. 4. Análisis de estabilidad. Dada la función de transferencia

s2
H ( s) 
s  4s  3
2

se trata de estudiar su estabilidad: 1) obteniendo sus autovalores (por tratarse de un

SLI), y 2) aplicando el segundo método de Lyapunov (de gran utilidad en los sistemas
no lineales).

Se pide:
1) Calcular los autovalores de dos de sus realizaciones y comprobar que coinciden con
los polos de H ( s) . (Matlab: roots, tf2ss, canon, ssdata, eig)
2) A partir de la matriz del sistema correspondiente a la Forma Companion Controlable,
resolver la ecuación de Lyapuov ( AT P  PA   I ) y determinar el signo de la forma
cuadrática cuya matriz de coeficientes es P, mediante el criterio de Sylvester. (Matlab:
lyap, eye, det)

Ejercicio 4. 5. Controlabilidad y observabilidad. Estabilizabilidad y detectabilidad.

Dado el sistema definido por las siguientes matrices:

2 8 0 0  1
0   
, c t  0 1 1 0
3 6 0 1
A  , b 
0 0 4 0  0
   
0 0 0  5  0

se trata de estudiar su controlabilidad (estabilizabilidad) y su observabilidad

(detectabilidad).
Se pide:
 
1) Calcular la matriz de controlabilidad M c  b Ab  A n 1b , determinar por
inspección su rango e interpretar el resultado. (Matlab: ctrb,rank)

ETSETB. Sistemas Electrónicos de Control 0910a 13

 Tema 3. Ejercicios adicionales (completos)

 ct 
 t 
2) Calcular la matriz de observabilidad c A , determinar por inspección su
Mo  
  
 t n 1 
c A 
rango e interpretar el resultado. (Matlab: obsv, rank)
3) Hallar la Forma Canónica Diagonal e identificar los modos observables y/o
controlables. (Matlab: canon)
4) A la vista del resultado anterior indicar si el sistema es estabilizable y/o detectable.

ETSETB. Sistemas Electrónicos de Control 0910a 14

 Tema 3. Ejercicios adicionales (completos)

5. Control robusto

Ejercicio 5. 1. Control Robusto. Principio del modelo interno. El hecho de que el

regulador capaz de eliminar el offset a una entrada o a una perturbación en escalón sea el tipo I
(integrador) y que éste sea justamente el modelo matemático del escalón, no es una casualidad.
La idea puede extenderse al seguimiento asintótico de señales no transitorias (escalones,
oscilaciones, rampas, exponenciales positivas, etc.) poniendo como controlador un
Gc ( s )  N c ( s ) / Dc ( s ) siendo Gc ( s ) la transformada de Laplace de la señal que se quiera seguir (o
filtrar). A dicha función de transferencia se le pueden añadir ceros adicionales para conformar el
transitorio.
Para ilustrar lo antedicho vamos a controlar dos plantas:
1
Planta A) GP ( s )  , de manera que no presente error permanente a entradas en escalón y,
s3
además, que su polinomio característico sea  c ( s)  s2  5s  25 .
1
Planta B) GP ( s)  , de manera que siga asintóticamente a una sinusoide con   1 y
s( s  1)
cuyo polinomio característico sea  c ( s)  s2  5s  25 .

Figura 11 Figura 12 

Planta A) Dada la estructura de la Figura 11, en la que k2 (¡retroacción de estado!) está para
k
modificar los polos y 1 (corrector serie) para añadir precisión (ceros), se pide:
s
1) Calcular M ( s) y determinar por inspección k1 y k2 .
2) Calcular E y comprobar que, efectivamente, presenta un cero que bloquea la señal de
R
continua.
3) Con objeto de acelerar la respuesta se añade un término anticipativo Nr, de manera que:
u  k1  edt  k2 y  Nr
3.1) Dibujar el nuevo esquema de bloques con la nueva señal Nr.
3.2) Repetir el dibujo en caso de que r no sea accesible pero sí lo sea e y, obviamente, y.
3.3) Calcular M ( s) y comprobar que la señal Nr ha introducido un cero, ajustable a través del
parámetro N.
Simulación: Dibujar la respuesta indicial en los siguientes casos
1) N  0
2) N  5
3) N  0 , y el polo (-3) de la planta perturbada al valor (-1)
4) Ídem perturbado a (-6)

Planta B) Considerar la estructura de la Figura 12.

4) Calcular k3 para modificar los polos de la planta a  c ( s)  s2  2.5s  1.

ETSETB. Sistemas Electrónicos de Control 0910a 15

 Tema 3. Ejercicios adicionales (completos)

E
5) Calcular M ( s) y y comprobar la existencia de ceros que pueden bloquear las señales
R
sinusoidales.
Simulación: Dibujar la respuesta a r  4 sen( t ) en los siguientes casos:
5) k1  k2  1
6) k1  k2  5
7) k1  1 ; k2  5
8) k1  k2  1 ; si el polo (-1) ha sido perturbado y su nuevo valor es (-0.5)
9) Ídem, perturbado a (-2).

Ejercicio 5. 2. Control robusto: Ceros de anclaje y prefiltro.

Enfoque del diseño:
1) Se elige Gc de modo que se satisfaga una adecuada sensibilidad en el intervalo de frecuencias
de interés.
2) Se elige GF de manera que se asegure una adecuada transmisión (M).
Datos:
2500
- G p ( s)  k
s( s  25)
- k tiene un valor nominal de 1 pero puede variar entre 0.5 y 20.
- Polos deseados: -12.4  j9.6

Se trata de:
1  10s
1) Comprobar que un diseño convencional con un compensador de retardo Gc ( s) 
1  100s
presenta una fuerte sensibilidad (frecuencial, de los polos y temporal) a las variaciones de k.
2) Comprobar que esto puede mejorarse con un controlador robusto en serie con ceros cercanos
s 2  26s  269
a los polos deseados G c ( s)  y con un prefiltro (ver Figura 13) para
269
269
cancelar los ceros introducidos GF ( s)  2 .
s  26s  269

D
+
+ + Y
GF(s) Gc(s) Gp(s)
-

Figura 13

Para ello se pide:
1) Con corrector convencional (retardo):
1.1) Representar el LGR de Evans. Indicar los polos para k = 0.5, 1, 2.
1.2) Representar la respuesta indicial: Comparar los tr (0.1% a 0.9%) para los valores
indicados de k.
2) Con corrector robusto y prefiltro:
2.1) Representar el LGR de Evans. Indicar los polos para k = 0.5, 1, 2.
2.2) Representar la respuesta indicial: Comparar los tr (0.1% a 0.9%) para los valores
indicados de k.
3) Efecto de las perturbaciones:

ETSETB. Sistemas Electrónicos de Control 0910a 16

 Tema 3. Ejercicios adicionales (completos)

Y ( s)
3.1) Calcular M D ( s)  (Nota: Obsérvese que en este caso coincide con la función de
D( s )
sensibilidad).
3.2) Representar su módulo (Bode) en las siguientes condiciones:
- Servo sin compensar.
- Servo con compensador de retardo.
- Servo con compensador robusto y prefiltro.
(¿Por qué aquí no se respeta la condición de la integral de la sensibilidad?)

ETSETB. Sistemas Electrónicos de Control 0910a 17

 Tema 3. Ejercicios adicionales (completos)

6. Dos grados de libertad

 Ejercicio 6. 1. Feedforward. Tracking perfecto de una r(t) cuya forma de onda es

conocida (continua y derivable). El sistema ahora considerado es el de la Figura 14

GFF(s)

+
r + 1 3 y
T(s) +
R(s) s(s-1)
-

S(s)
Figura 14
y se trata de que y(t) siga una señal cuyo perfil de evolución temporal es r (t )  2  t 2 . El
controlador elegido es T ( s)  ( s  2)( s 2  2.4s  2.44) , S ( s)  4.2 s  163
. , R ( s)  s  5.4 ,
GFF ( s)  0.333s 2  0.33s  4.2 .
Se pide:
1) Calcular M(0).
2) Calcular L(s)
1
3) Comprobar que G FF ( s)  .
G p ( s)
4) Hallar la respuesta y(t) y compararla con r(t).

 Ejercicio 6. 2. Controlador con dos grados de libertad (2 DOF). Ajuste

independiente de seguimiento y robustez mediante el procesado independiente de y(t) y
r(t). Dado el servosistema de la Figura 15 se trata de conseguir que sus características de
2.44
transmisión (seguimiento) correspondan a las del modelo M m ( s)  elegido
s  2.4s  2.44
2

de forma que presente poco rebase ( 3%) y un tr  2.5s. Y, por otra parte, que se puedan
modificar las características de sensibilidad (lazo) sin modificar por ello la respuesta del servo.

+ 1 3
T(s) R(s)
- s(s-1)

S(s)
Figura 15
Para ello se eligen tres controladores:
1) T1 ( s)  0.81( s  0.2) ; R1 ( s)  s  3.6 ; S1 ( s)  2.17s  0.16
2) T2 ( s)  0.81( s  2) ; R2 ( s)  s  5.4 ; S2 ( s)  4.21s  163
.
3) T3 ( s)  0.81( s  20) ; R3 ( s)  s  3.23.4 ; S3 ( s)  4.6s  16.3 
Se pide:
1) Estimar a priori tr y Rpt de la respuesta indicial del modelo Mm(s).

ETSETB. Sistemas Electrónicos de Control 0910a 18

 Tema 3. Ejercicios adicionales (completos)

2) Calcular y representar a escala la respuesta indicial del modelo Mm(s).

3) A partir de la gráfica determinar tr y Rpt. Comparar con 1).
4) Hallar la expresión matemática de u(t) a partir de la de y(t). calcular M(0) aplicando el
Teorema del Valor Inicial.
5) Calcular y representar (Bode de amplitud) L(j) para el primer caso.
1
6) Ídem para la función de sensibilidad S ( j )  .
1  L ( j )
7) Determinar el valor de S   S ( j ) max en cada caso.
8) ¿Asegura una estabilidad robusta frente a cambios paramétricos?
Simulación:
1) Representar la y(t) en los tres casos. Comprobar que es la misma.
2) Representar la u(t) en los tres casos. Ídem.
3) Representar L ( j ) en una misma gráfica. Comparar.
4) Representar S ( j ) en una misma gráfica. Comparar.
5) Comprobar la robustez del diseño:
5.1) (MG) Comprobar que permanece estable al variar el parámetro de 3 a 5.
5.2) (MF) Ídem si la planta presenta un retardo de e-0.15

 Ejercicio 6. 3. Métodos Polinómicos. Controlador RST. Se trata de estabilizar la planta

b( s) 1
G ( s)   2 mediante los dos controladores de la configuración mostrada en la
a ( s) s  1
Figura 16.

Figura 16
Se pide:
Y ( s)
1) Hallar la expresión de M ( s)  .
R ( s)
2) Hallar la solución (polinomios p ( s) y q ( s) ) de la ecuación diofántica obtenida al fijar el
polinomio característico del servo a DM ( s)   C ( s)  s2  s  0.5 y siendo (por elección)
d ( s)   O ( s)  s 2  2 s  2 .
db b
3) Comprobar que M  se reduce a M  .
pa  bq c
4) Determinar la respuesta indicial y comprobar que se ha estabilizado. 

Ejercicio 6. 4. Controladores con dos grados de libertad. Dados los controladores de la

Figura 17 y la Figura 18.

ETSETB. Sistemas Electrónicos de Control 0910a 19

 Tema 3. Ejercicios adicionales (completos)

C3(s)

+
+ + N1(s) + N2(s)
C1(s) C2(s) G(s) G(s)
- D1(s) D2(s)
- -

C4(s) C5(s)

Figura 17 Figura 18

se pide:
1) En el caso de la Figura 17, calcular las funciones de transferencia L(s), S(s) y M(s).
2) En el caso de la Figura 18 calcular los polinomios R, S y T definidos en la Figura 19.
+ 1
T G(s)
R
-

Figura 19

Ejercicio 6. 5. Diseño con un grado de libertad (1 DOF). Mantenimiento de las

características del lazo y degradación de la respuesta. El sistema ahora considerado es el de
la Figura 20

+ S(s) 3
R(s) s(s-1)
-

Figura 20
tomando el controlador S/R los tres pares de valores del .
Estudio Previo:
1) Calcular la expresión de L(s) y comprobar que es la misma
2) Comprobar que M(0) es en los tres casos la misma
3) Estimar gráficamente, a partir de la gráfica de y(t) anterior, la forma de las respuestas y(t) en
los tres casos (Nota: Al tener un prefiltro con un cero (s + z), y z ( t )  y ( t )  zy ( t ) )
4) A partir de las gráficas de y(t) estimar las de u(t). Comentar el posible problema de
saturación.
Simulación:
1) Representar y(t) en los tres casos. Comparar con lo previsto.
2) Representar u(t). Comparar. 

Ejercicio 6. 6. Regulador de mínima varianza. En este caso, al ser un problema de

regulación (que no de seguimiento) la solución se reduce a la simple asignación de los polos,
objetivo que puede alcanzarse con sólo la parte retroactiva  S / R . No obstante, resulta
ilustrativo constatar que un problema originalmente tratado por variables de estado (LQG)
acepta también una formulación polinómica cuya solución depende de una (¡otra!) ecuación
diofántica.

ETSETB. Sistemas Electrónicos de Control 0910a 20

 Tema 3. Ejercicios adicionales (completos)

Problema: Dada una planta descrita por el modelo ARMAX :

y ( n  1)  0.8 y ( n )  u ( n )  ( n  1)  0. 7 ( n )
siendo ( n ) una secuencia de ruido blanco gaussiano de media nula y varianza 2 ( n) , se trata
de generar una señal de control u ( n ) que reduzca lo más posible el efecto del ruido. Efecto que,
sobre la salida, cuantificaremos mediante el criterio J ( n )  E y 2 ( n ) , es decir, mediante la
varianza. 
Estudio Previo:
1) Dibujar un esquema de bloques que relacione u ( n ), ( n ) y y ( n ) . Identificar la H ( z ) del
filtro generador.
2) Razonar el término "modelo ARMAX".
3) Si la solución en este caso es (con referencia a la Figura 21) T  0 , R  FB , S  q 1G , donde
F y G se obtienen de C ( q 1 )  A( q 1 ) F ( q 1 )  q  ( d 1)G ( q 1 ) , resolver la ecuación
U ( z)
diofántica (por identificación de coeficientes) y determinar la Hc ( z ) del corrector  .
Y ( z)

Figura 21
Nota: Estructura de los polinomios:
A( z 1 )  1  a1z 1 an z  n
B( z 1 )  b1z 1 bn z  n
C ( z 1 )  1  c1 z 1  cn z  n
F ( z 1 )  1  f 1 z 1  f d z  d
G ( z 1 )  g 0 z 1  g n 1 z  ( n 1)
Simulación:
1) Obtener la respuesta, en lazo abierto, a una entrada aleatoria.
2) Ídem en lazo cerrado y comparar los resultados. 

Ejercicio 6. 7. Control anticipativo (feedforward) de las perturbaciones y señal de

mando. Se trata de estudiar el efecto de la acción feedforward (FF), tanto de las perturbaciones
como de la señal de mando.

Figura 22 Figura 23 
Estudio Previo:
1) Dibujar aproximadamente la forma de y ( t ) (sin FF en la Figura 22) cuando p ( t ) es un
impulso unitario.

ETSETB. Sistemas Electrónicos de Control 0910a 21

 Tema 3. Ejercicios adicionales (completos)

2) Si la sensibilidad del sensor es k s  1, calcular la G FF ( s ) necesaria para reducir el efecto de

la perturbación p ( t ) .
3) Dibujar aproximadamente las señales u( t ) y y ( t ) (sin y con FF en la Figura 23) cuando r ( t )
es un escalón unitario.
Simulación:
1) Representar la respuesta a un impulso de p ( t ) (sistema de la Figura 22) sin y con FF.
2) Representar la respuesta a un escalón de r ( t ) (sistema de la Figura 23) sin y con FF. 

Ejercicio 6. 8. Retroacción condicional. Dado el sistema de la Figura 24, con

100
G p ( s)  ,
( s  1)( s  5)
w
+
R + + Y
Gp(s)
+
-
Hc(s)
+

Gc(s)

Figura 24
se pide:
Y
1) Hallar si Gc(s) = Gp(s).
R
2) En el supuesto de diseñar una buena regulación frente a una perturbación w de tipo escalón,
seleccionar la Hc(s) más adecuada de entre las siguientes opciones:
10 10 10( s  1) k
H 1 ( s)  , H 2 ( s)  , H3 ( s)  , H 4 ( s)  2 . 
s( s  1) ( s  1)( s  2) ( s  2) s

Ejercicio 6. 9. Compensación con prefiltro (forward compensation) y anticipativa

(feedforward). Aunque presenta varios parámetros de ajuste la compensación serie sólo tiene un
grado de libertad y así los objetivos quedan limitados no pudiendo, a veces, superar la
contradicción entre estabilidad y sensibilidad (que requiere ganancia alta). Una solución es
recurrir a estrategias con dos grados de libertad, así
1) El prefiltro (forward) en cascada con el lazo (con su Gc serie).
2) El anticipativo (feedforward) en paralelo con Gc.
En ambos casos el corrector no afecta los polos del lazo. Se elige de manera que cancele polos y
ceros. Vamos a ilustrarlo con un ejemplo: control con un grado de libertad (Gc) ampliado a dos
grados por dos métodos: a) GF y b) GFF.

GFF(s)
+
GF(s) + Gc(s) Gp(s) + Gc(s) + Gp(s)
- -

ETSETB. Sistemas Electrónicos de Control 0910a 22

 Tema 3. Ejercicios adicionales (completos)

(a) (b)
Figura 25
2500(1  10s)
Nota: Tomar G c ( s) G p ( s)  . 
s( s  25)(1  100s)
Estudio Previo:
1) Calcular GFF equivalente a GF ( s)  1  0.05s
Simulación:
1) Hallar Rpt y tr del sistema básico (control con Gc).
2) Para hacerlo más rápido sin degradar excesivamente la estabilidad se añade
GF ( s)  1  0.05s . Hallar de nuevo Rpt y tr.
3) Con la GFF calculada anteriormente, hallar Rpt y tr. 

Ejercicio 6. 10. Método de Truxal-Guillemin. Se trata de familiarizarse con el método

directo de Guillemin-Truxal y, en particular, con la construcción (síntesis) de M ( s) a partir de
las especificaciones de la respuesta temporal (Análisis Inverso). En concreto se trata de diseñar
C ( s) 1
un sistema de control para una planta con G p ( s)   cuya dinámica se quiere
U ( s) s( s  1)
modificar de manera que su salida c( t ) responda a las señales de mando r ( t ) tipo escalón según
las siguientes especificaciones:
Estabilidad: R pt  5% (  0. 7 , si el sistema es de segundo orden)
Precisión: e(  )  0
Velocidad: t p  0.1
Realizabilidad: Exceso P  Z  2
Adicional: k v  50 , en caso de entrada en rampa 
Estudio Previo:
1) Construcción de M ( s) .
1.1) Partiendo del concepto de polos dominantes: Estimar  n y . Añadir un polo adicional
(para la realizabilidad), un cero adicional (para la precisión) y ajustar la k. (Nota:
Recuérdese que el cero adicional disminuye el t p y aumenta el R pt )
1.2) Partiendo de un denominador con polinomios óptimos: Elegir un ITAE del orden
adecuado (3º en este caso) con la  n indicada por la velocidad y determinar el
numerador necesario para que el tipo del sistema sea 1 y la k v la especificada (ver
Apéndices).
2) Estimar (en caso necesario) los polos y ceros y trazar aproximadamente la respuesta indicial.
Simulación:
1) Polos dominantes
1.1) A partir de la M(s) obtenida anteriormente, calcular el corrector serie necesario para
realizarla.
1.2) Calcular k v a partir de L ( s) y a partir de M ( s) .
1.3) Hallar la respuesta indicial del servo y comprobar los valores de  y t p .
2) Polinomio ITAE
2.1) Igual que en 1.1).
2.2) Igual que en 1.2).
2.3) Igual que en 1.3). 

Ejercicio 6. 11. Asignación de polos por métodos polinómicos (problema regulador).

Controlador con dos grados de libertad y observador implícito. Ecuación diofántica.

ETSETB. Sistemas Electrónicos de Control 0910a 23

 Tema 3. Ejercicios adicionales (completos)

N ( s) s5
Problema: Dada la planta G ( s )   2 se trata de diseñar un controlador con
D ( s ) s  3s  2
observador implícito, o su equivalente serie/paralelo (ver figuras), con objeto de modificar sus
polos.
Solución:
- Comentario previo: Esta configuración viene sugerida por los métodos desarrollados en la
Teoría de Control Moderno, y su solución (condicionada a que N y D sean polinomios primos
entre sí y a que los polinomios N1 y N 2 sean de orden inferior al de D0 ) se lleva a cabo por
métodos polinómicos (ecuación diofántica)

Figura 26 Figura 27

- Datos: Polinomio característico deseado: DM  s  6s  21  c (en los métodos de estado)

2

Observador, elegido con una dinámica: D0  s2  15 s  50   0 

Estudio Previo:
' D D .
1) Calcular DM M 0

2) Resolver la ecuación diofántica D N1  D0  NN 2  DM  

' . (Nota: Hacerlo por
identificación de coeficientes, habiendo tenido en cuenta las restricciones, en cuanto a su grado,
indicadas en el comentario previo).
Simulación:
1) Comprobar la solución de la ecuación diofántica.
2) Comparar las respuestas indiciales del sistema antes y después de ser controlado. 

Ejercicio 6. 12. Fijación de polos y ceros (problema servo). Seguimiento de un modelo.

Una determinada planta (previamente discretizada) viene descrita por
B( z) z  0. 8
G(z)   , y se desea controlarla de manera que la transmisión desde la señal de
A( z ) ( z  1)( z  0. 9 )
Bm ( z ) 1. 08 z
mando (problema servo) venga descrita por el modelo Mm ( z )   . Se pretende
Am ( z ) z2  0.56 z  0. 64
resolver este problema mediante la configuración (con dos grados de libertad) de la Figura 28,
donde R, S, T, A y B son polinomios en z, y donde se supone un observador (implícito) con
polinomio característico D0  1 (  A0 a partir de ahora siguiendo la nomenclatura de Åström).

Figura 28 
Estudio Previo:

ETSETB. Sistemas Electrónicos de Control 0910a 24

 Tema 3. Ejercicios adicionales (completos)

1) ¿Por qué, dentro de la limitación de su grado, hemos podido elegir arbitrariamente Am y A0 ?

Y
2) Calcular M  .
Yr
3) Suponiendo que conseguimos M  M m , y teniendo en cuenta que D M'  Am A0 B  (*),
escribir las relaciones entre los numeradores y denominadores de M y M m .
4) Resolver la ecuación diofántica AR  BS  Am A0 B  (que tiene solución (**) por cuanto
A y B son primos entre sí).
T S
5) Calcular los controladores y y verificar que se alcanzan las especificaciones de
R R
transmisión. ¿Qué puede decirse sobre la sensibilidad a variaciones lentas de los parámetros
de la planta?
(*) Nota: Esta ecuación se demostrará en clase, y es una consecuencia del teorema de
separación que aparece en el diseño del regulador por retroacción del estado, cuando se usa una
estimación del mismo obtenida mediante un observador, en cuyo caso resulta D M'   0 c . Al
ser extendido al caso servo aparece el factor B  proveniente de los factores estables de
B ( B  B B ).
(**)Nota: Para que, además de existir, la solución sea única hay que elegir, en este caso, S de
primer grado y R de grado cero.
Simulación:
1) Comprobar la solución de la ecuación diofántica.
2) Hallar la respuesta del sistema (y, u).
3) Valorar la sensibilidad del sistema. 

ETSETB. Sistemas Electrónicos de Control 0910a 25

 Tema 3. Ejercicios adicionales (completos)

7. Control óptimo

Ejercicio 7. 1. Solución paramétrica. Se trata de controlar la planta del ¡Error! No se

encuentra el origen de la referencia. optimizando, ahora, su dinámica. Para ello se parte de la
siguiente configuración (k>0):

Figura 29
Estudio Previo:
1) Razónese la elección del algoritmo de control (el polo, los ceros y el signo de k).
2) Bosquejar el lugar de Evans correspondiente al parámetro k.
3) Calcular los valores de k que estabilizan el servo.
4) Estimar el valor de k para obtener unos polos dominantes con   0. 7 .
5) Calcular el valor del criterio ISE ( J 3 ( k ) ) a una entrada en escalón (ver Apéndices).
6) Valorar la dificultad de hallar analíticamente el valor de k que minimiza el criterio J 3 .
Simulación:
1) Trazar el lugar de Evans y verificar los resultados obtenidos en los apartados 2), 3) y 4) del
estudio previo.
2) Hallar el valor óptimo de k por medio de los siguientes procedimientos:
2.1) Representando J 3 ( k ) y estimando el mínimo gráficamente.
 J3
2.2) Resolviendo la ecuación  0.
k
2.3) Usando directamente la función de MATLAB fmin (Notar la eficacia de esta función en la
detección de mínimos relativos).
3) Con el valor de k obtenido en el apartado 4) del estudio previo trazar la respuesta indicial y,
sobre la gráfica, calcular aproximadamente el valor del ISE. Comparar el resultado con el
obtenido al calcular J 3 ( k ) . 

Ejercicio 7. 2. Solución funcional. Dada la planta del ¡Error! No se encuentra el origen de la

referencia., se trata de hallar la función de transferencia óptima del servo, correspondiente a una



entrada en escalón y siendo J  ( e 2 ( t )  u 2 ( t )) dt el criterio a optimizar, mediante la
0

aplicación de la factorización espectral de Wiener. (Nota: La inclusión de u 2 ( t ) introduce

realismo en la optimización al restringirla indirectamente a valores "pequeños" de u( t ) -no
saturación). 
Estudio Previo:
1) Realizar la factorización espectral solución de NN   2 DD (ver Apéndices).
2) Calcular M O ( s) .
3) Calcular el controlador serie (con H=1) necesario para obtener dicha M O ( s) .
Simulación:
1) Obtener la factorización espectral calculando las raíces de NN  2 DD  0 .
2) Trazar la respuesta impulsional y comprobar que el sistema se ha estabilizado.
3) Calcular el valor del criterio ISE a un escalón unitario. 

ETSETB. Sistemas Electrónicos de Control 0910a 26

 Tema 3. Ejercicios adicionales (completos)

Ejercicio 7. 3. Solución Paramétrica (estructura fija). Se trata de estudiar el efecto del ruido
de entrada sobre el error del servo (medido en términos de e2 ) y estudiar cómo depende de los
parámetros del sistema.
N E 1 C
+
 s 
s  2 
 n 

Figura 30
En este caso se trata del servo de la Figura 30, siendo  el parámetro ajustable y n( t ) una señal
aleatoria cuyo espectro de potencia es:
S0
Sn ( )  2 

1  
 0 
Estudio Previo:
E ( s)
1) Calcular la función de transferencia .
N ( s)
2) Calcular la densidad espectral del error Se (  ) .
3) Calcular e 2 ( ) , utilizando las tablas de los Apéndices. (Para simplificar la expresión,
n
utilícese la frecuencia normalizada u  )
0
4) Hallar el valor óptimo de  y particularizarlo para el caso  0  10 n .
Simulación:
1) Generar 100 muestras de ruido n( t ) blanco gaussiano de media cero y varianza unidad.
2) Simular la respuesta del sistema a n( t ) en lazo cerrado para  n  1 y    0 .
3) Obsérvese la respuesta e( t ) para los valores anteriores y calcular e 2 cuando    0 .
4) Comprobar la reducción del ruido al variar  y acercarse a su valor óptimo. 

Ejercicio 7. 4. Regulador óptimo determinista (LQ). Solución polinómica. Se trata de

comprobar, una vez más, la versatilidad del controlador con dos grados de libertad, aplicado
aquí a implementar un regulador óptimo LQ. La solución de este problema tiene, como paso
previo, la solución de una ecuación diofántica (muy sencilla) que permite obtener un predictor
de d pasos.
0. 01z  0. 01
Problema: Se requiere regular una planta (lineal) H ( z )  z 2 2 de manera que
z  2. 002 z  1. 002
los estados perturbados se recuperen de forma que se minimice el criterio (cuadrático)

 . un  un 1  .
J (n)  y 2 (n  d )  01
2

Solución: La teoría de control óptimo demuestra que la solución viene dada ¡también! por la
configuración de la Figura 28, satisfaciendo aquí los polinomios R ( z 1 ), S ( z 1 ), T ( z 1 ) las
siguientes expresiones:
r
R  PBF  0 Ru ; S  PG ; T  Q
b0
siendo:

ETSETB. Sistemas Electrónicos de Control 0910a 27

 Tema 3. Ejercicios adicionales (completos)

P, Q, Ru los polinomios de ponderación de J ( n)

que, en este caso, valen P  1 , Q  0 , y Ru  0.1  0.1z 1.
B ( z 1 ) el numerador de la planta.
r0 , b0 los términos constantes de los polinomios Ru ( z 1 ) y B ( z 1 ) respectivamente.
F, G los polinomios solución de la ecuación diofántica 1  A( z 1 ) F ( z 1 )  z  d G ( z 1 ) . 
Estudio Previo:
1) Comprobar que la planta es inestable.
1
2) Resolver la ecuación diofántica realizando la división: , de donde resulta que F es el
A ( z 1 )
cociente y G el resto.
3) Calcular R , S , T .
4) Hallar los polos del sistema y comprobar su estabilidad.
Comentario sobre las notaciones: Obsérvese que la G ( z ) de la ecuación diofántica no es la
función de transmitancia de la planta, aquí representada por H ( z ) . El cambio en la notación es
debido a la generalización de su uso, extendido por Åström.
Simulación:
1) Calcular M ( z ) a partir de la estructura y componentes del sistema.
2) Determinar la respuesta a una entrada de r ( n ) en escalón.
3) Determinar la respuesta a una perturbación impulsional, p( n) , a la entrada de la planta. 

Ejercicio 7. 5. Comparación de métodos de cálculo. Se trata de controlar una planta

1
G p ( s)  , con una misma estructura de control (H = 1 y corrector serie), para satisfacer
( s  1)
especificaciones diversas y por diferentes métodos de cálculo:
1  ns
1) Bode. Corrector de avance Gc ( s)  k . Especificaciones: MF = 45, kv = 10.
1  s
sz
2) Evans. Corrector de avance G c ( s)  k . Especificaciones (dinámica del servo):  =
s p
0.5; n = 2.
3) Truxal. Dinámica del servo:  = 0.5; n = 2. Precisión: M(0) = 1.

4) Optimización paramétrica (frente a entradas en escalón): Gc = k, J   0
( e 2  u 2 ) dt

5) Optimización funcional (frente a entradas en escalón): J  0
( e 2  u 2 ) dt
3
6) Optimización estocástica. Señales de entrada: Referencia: S r ( s)  ; Ruido aditivo no
1 s2
correlado Sn = 1. Criterio: min e 2   

ETSETB. Sistemas Electrónicos de Control 0910a 28