Sistemas Control
Sistemas Control
REGULACIÓN AUTOMÁTICA
pa FT
w
FC TC
Fv
TT
q T
Condensado
Francisco Vázquez
Jorge Jiménez
Fernando Morilla
1.1 INTRODUCCIÓN
El objetivo de un Ingeniero de Control es, por tanto, diseñar el equipo necesario para que
actúe sobre las entradas de un sistema de manera que las salidas se mantengan en un valor
predecible y modificable. El término sistema se refiere a toda combinación de componentes que
actúan conjuntamente y cumplen un determinado objetivo, no estando el concepto limitado a
objetos físicos. Dicho concepto puede aplicarse a fenómenos dinámicos abstractos, como los que
se encuentran en economía. Por tanto, el término sistema hay que interpretarlo como referido a
entes físicos, biológicos, económicos y otros. Debido a esta variedad de posibles sistemas a
controlar, el ingeniero de control deberá trabajar en colaboración con otros especialistas o utilizar
información bibliográfica adicional que le ayuden a conocer con más profundidad el sistema a
controlar, de forma que le permita establecer la estrategia de control a seguir y los medios a
utilizar.
Así, pueden encontrarse muchos ejemplos habituales como podrían ser los siguientes:
Esta edad de oro de la automática, que se denomina ya como Automática Clásica tuvo
muchos más protagonistas, como W.R Evans que en 1948 publicó un trabajo que permitía
conocer de forma gráfica la estabilidad de un sistema (en concreto estuvo relacionado con
sistemas de vuelo) permitiendo analizar sistemas con dinámicas estables o condicionalmente
estables. Este método, denominado técnica del lugar de las raíces, permite representar
gráficamente la evolución de las raíces de la ecuación característica cuando algún parámetro
del sistema variaba.
Durante la década de los años 50 varios autores, entre ellos Bellman, Kalman y
Pontryagin, vuelven a considerar las ecuaciones diferenciales ordinarias como la herramienta
matemática básica para el estudio de los sistemas de control, dando lugar a nuevas técnicas que
han venido a conformar lo que se denomina Automática Moderna. Ello estuvo motivado por un
nuevo campo del control (el control de satélites artificiales) y por la aparición de las computadoras.
Como la descripción de sistemas de control de las plantas modernas (con muchas entradas y
salidas) requiere una gran cantidad de ecuaciones, la teoría clásica que trata sistemas de entrada y
salida única, se vuelve absolutamente impotente. Esto ha provocado un rápido desarrollo de la
teoría de control moderna a partir de los años 60.
deterministas como estocásticos), que permiten obtener un control que minimiza alguna función de
coste, en el estudio de sistemas de control adaptativo (que permiten la autosintonía de los
controladores ante un cambio en el sistema que controlan) o en la incorporación de técnicas de
sistemas expertos a los sistemas de control. El desarrollo actual está principalmente encauzado
hacia el control robusto, con cuya teoría se obtienen controladores que consiguen sus objetivos
para una familia de plantas (no sólo para una planta nominal) obtenidas al utilizar modelos de las
perturbaciones que afectan a la planta nominal en el modelado de éstas o hacia el denominado
control predictivo, que permite tener cierto carácter anticipativo a las acciones de control.
Se ha visto que un sistema se define como una combinación de componentes que actúan
conjuntamente y cumplen un determinado objetivo. Más concretamente, en un sistema de control
el objetivo principal es conseguir que las salidas del sistema evolucionen de una forma adecuada
para unas entradas determinadas. La figura 1.2 representa un sistema de control cualquiera.
Perturbaciones
Señal de Señal de PV o
Referencia OP Salida
Error Control
+ Controlador Actuador Proceso
-
Salida
medida
Sensor
Al introducirnos en el mundo del Control Automático se tienen que tener claro cuáles son
sus componentes básicos, las señales y las variables de todo Sistema de Control. Los componentes
básicos principales son:
El proceso: Se define como una operación artificial o voluntaria que consiste en una serie
de acciones controladas o movimientos dirigidos sistemáticamente hacia un determinado resultado
o fin. En definitiva se denomina proceso a cualquier operación que se pueda controlar. En algunos
casos, en lugar de proceso se habla de planta, que es un término muy utilizado en la industria,
como un objeto físico que ha de ser controlado (por ejemplo un motor, un intercambiador de calor,
un satélite artificial, etc...). El objetivo de control vendrá impuesto sobre una o varias de las salidas
de este proceso y el ingeniero de control tendrá que decidir sobre qué variable de entrada debe
actuarse sin modificar la estructura física del proceso.
El actuador: Para actuar sobre las entradas del proceso se debe convertir la señal
generada por el controlador (normalmente de tipo eléctrico de baja potencia) en otra variable
física. Por ejemplo, una válvula controlada por un motor paso a paso convierte la tensión de
entrada en forma de impulsos en una apertura de válvula y, consiguientemente, en un cambio de
caudal másico. A veces el controlador asume esta función. Sin embargo, cuando se modela un
proceso, el efecto de este elemento se incorpora en el modelo de la planta y en el diseño o la
elección del actuador también debe encargarse el ingeniero de control.
acumulada (historia pasada del error) y su derivada (es decir, cómo se modifique su velocidad de
cambio con respecto al tiempo).
La señal de control: Es la señal que genera el controlador para gobernar el elemento
final de control (actuador). El controlador podrá ser analógico, en cuyo caso las entradas y las
salidas son magnitudes eléctricas continuas o digital, en cuyo caso deben existir conversores A/D y
D/A para adaptar las señales digitales a analógicas.
Ejemplo 1.1
En la figura 1.3 se puede ver un intercambiador de calor. El líquido que fluye a través de la
tubería (en forma de serpentín), se calienta mediante vapor de agua que se inyecta en el
intercambiador de calor. Por tanto, la misión del intercambiador de calor es calentar el líquido
desde una temperatura de entrada Ti(t) hasta una temperatura de salida To(t).
En condiciones ideales, el intercambiador y las tuberías están bien aisladas y, por lo tanto, la
energía suministrada al líquido es igual a la energía perdida en la condensación del vapor de
agua. En este proceso hay muchas variables que pueden cambiar, por ejemplo, la temperatura
Ti(t) o el caudal de líquido q(t) de entrada, provocando que la salida del proceso se desvíe del
valor deseado. El objetivo del sistema de control será mantener la temperatura del líquido de
salida a un valor determinado a pesar de las variaciones en otras variables. Una forma de
conseguir este objetivo es:
3. Basado en esta comparación decidir lo que se debe hacer para corregir cualquier
desviación.
Una posibilidad de acción consistiría en variar la cantidad de vapor de agua que llega al
intercambiador. La medida de To(t) se realiza mediante un sensor (termopar, termorresistencia,
etc.). Este sensor se conecta físicamente a un transmisor que convierte la salida del sensor en
una señal eléctrica (tensión o corriente, de suficiente nivel para evitar ruidos e interferencias) y
la transmite al controlador. El controlador recibe la señal eléctrica (que está directamente
relacionada con la temperatura del líquido de salida), la compara con el valor de temperatura
deseado y, en base a esta comparación, envía una señal de control al actuador (que en este caso
es una válvula) que gobierna la cantidad de vapor de agua que llega al intercambiador.
En contraposición a los sistemas de control en lazo cerrado, están los sistemas de control
en lazo abierto, en los que la salida del proceso no tiene ningún efecto sobre la acción de control.
Un ejemplo práctico de control en lazo abierto es una lavadora automática: el remojo, lavado y
enjuague de la ropa se cumplen sobre una base de tiempos determinada, accionados por un
programa de lavado previamente seleccionado. Una perfecta lavadora automática debería tener
medios para medir la limpieza de la ropa a lo largo del lavado y ajustar los tiempos, por sí sola, a la
limpieza deseada. Queda así perfectamente claro que una lavadora es un sistema de control en
lazo abierto porque no mide la salida del sistema, es decir, la limpieza de la ropa. Otro ejemplo de
control en lazo abierto es el control de tráfico mediante semáforos. Sin embargo, si continuamente
se mide la cantidad de automóviles en espera en cada señal de tráfico y esa información se tiene
en cuenta para modificar el tiempo de cambio en los semáforos, el sistema se convierte en un
sistema de control en lazo cerrado. Esta práctica es cada vez más habitual en áreas urbanas
congestionadas por el tráfico. La teoría de autómatas finitos se encarga de describir
matemáticamente este tipo de dispositivos de control.
La principal ventaja de los sistemas de control en lazo cerrado sobre los de lazo abierto es
que el uso de la realimentación hace al sistema relativamente insensible a perturbaciones externas
y a posibles variaciones internas de los parámetros del sistema. Así, los objetivos básicos de
control, como ha quedado de manifiesto en los párrafos y ejemplo anterior pueden resumirse en
dos:
En otros sistemas, por el contrario, los cambios más importantes se dan en la señal de
referencia (que es por tanto una función del tiempo) y el control se diseña de manera que
la salida del proceso siga de la mejor forma posible a la señal de referencia. Este tipo de
control se le conoce como servomecánica y al conjunto del sistema como un
servomecanismo. El objetivo de control es el de seguimiento de trayectorias.
De acuerdo con esta clasificación, no del todo ‘pura’ ya que la mayoría de los sistemas de
control tienen ambos objetivos ponderados de alguna manera, se pueden poner algunos ejemplos
que aclaren ambos conceptos. De esta forma, el motor que mueve una escalera mecánica de unos
grandes almacenes o de una estación de trenes o metro, tiene un objetivo de control claro y es el
de rechazo de perturbaciones, vistas éstas como el aumento o descenso brusco de carga
soportada, o sea, el número de usuarios. En este caso la referencia, es decir, la velocidad de salida
del sistema (la velocidad de subida o bajada del usuario) debe ser lo menos variable posible (no se
cambian las referencias de velocidad en distintas horas del día ni en función de la carga) y debe de
mantenerse lo más constante posible ante un aumento de la carga en un momento de ocupación
máxima de la escalera. Por otro lado, el control de posicionamiento del lector láser del CD de un
ordenador tiene continuamente cambios de consignas, cambios en la referencia de posición
teniendo que hacer frente a este objetivo mediante exactos movimientos que lo sitúen en las
distintas coordenadas polares del CD.
Según el tipo de sistema de control: Esta clasificación puede hacerse mucho más
completa, pero a grosso modo podemos hablar de sistemas continuos, que son aquéllos
que en instantes de tiempo finitos pueden cambiar infinitas veces o de sistemas discretos,
que en instantes finitos pueden cambiar un número de veces finito. Este tipo de sistemas
está relacionado con la tecnología digital, y el uso de ordenadores o microprocesadores en
el lazo de control, de forma que éste necesita muestrear la variable salida, pasar un
proceso de acondicionamiento de estas señales (utilizando conversores A/D y D/A),
tratarlas mediante algún algoritmo discreto de control y, posteriormente, adaptar sus
niveles de tensión para poder atacar al “bloque” del actuador. Todo este proceso dura un
cierto tiempo, denominado período de muestreo, tras el cual puede comenzar la lectura
de otra muestra. De esta forma, el control discreto sólo da un nuevo valor de control en
instantes determinados, separados todos un espacio de T segundos, siendo T el período de
muestreo. Los sistemas de control continuo, por su parte, utilizan sistemas analógicos,
aunque la mayoría de estos están siendo sustituidos por los anteriores, gracias a su
capacidad de cálculo.
d 2 y (t ) dy (t )
2
P(t ) Q (t ) y (t ) R(t ) donde P(t), Q(t) y R(t) son polinomios en t. La
dt dt
propiedad más importante de los sistemas lineales es que se les puede aplicar el principio
de superposición. El principio de superposición establece que la respuesta de un sistema
producida por la aplicación simultánea de varias fuentes excitadoras es la suma de las
respuestas individuales. Los sistemas dinámicos que son lineales y están constituidos por
componentes invariables en el tiempo, se pueden describir por ecuaciones diferenciales
lineales con coeficientes constantes. Este tipo de sistemas reciben el nombre de sistemas
lineales invariantes en el tiempo (o lineales de coeficientes constantes). Los sistemas
representados por ecuaciones diferenciales cuyos coeficientes son funciones del tiempo,
reciben el nombre de sistemas lineales variantes en el tiempo. Un ejemplo de sistemas de
este tipo es el sistema de control de un vehículo espacial, dado que la masa del vehículo
cambia debido al consumo de combustible y la fuerza de la gravedad, que se modifican a
medida que el vehículo se aleja de la Tierra. El análisis y diseño de sistemas lineales
variantes en el tiempo es tan complejo que la mayoría de los casos se supondrá que los
parámetros son constantes o lentamente variables. Sistemas no lineales son aquellos
sistemas que se representan mediante ecuaciones diferenciales no lineales. Ejemplos
posibles de modelos de este tipo de sistemas son:
2
d 2 y (t ) dy (t )
y (t ) R (t )
dt 2 dt
d 2 y (t ) dy (t )
2
P( y) y (t ) 0
dt dt
Hablando estrictamente, todos los sistemas son no lineales pues aunque muchos
fenómenos físicos se describen mediante ecuaciones lineales, un estudio cuidado de los
sistemas físicos indica que incluso los denominados sistemas lineales lo son realmente en
un restringido rango de operación. Ejemplos típicos de no linealidades son: a) Las
saturaciones que se pueden producir en la salida de algunos componentes del sistema de
control para valores elevados de la señal de entrada; y b) las zonas muertas que afectan a
las señales de entrada pequeñas (una zona muerta es un rango en la entrada a la que no
responde el componente).
Ejemplo 1.2
Motor de corriente continua. Este tipo de motor se suele utilizar con bastante frecuencia en
aquellos procesos donde se requiere fijar una velocidad o una posición. Por ejemplo, en los
procesos de laminación en frío, la materia prima es un rollo de chapa de un grosor
determinado y el producto final es una chapa de espesor uniforme, menor que el inicial. Ello
se consigue haciendo pasar la chapa inicial entre unos rodillos fijos al mismo tiempo que se le
somete a una tensión constante. El espesor de la chapa resultante viene determinado por la
posición de los rodillos, pero ese espesor sólo será uniforme si la tensión en cualquier punto
del tren de laminado ha sido siempre la misma y la única forma de conseguir esto es mantener
una velocidad constante en el arrastre de la chapa, para lo que se utilizan motores de corriente
continua controlados en velocidad. En una máquina de escribir o una impresora con cabezal
de tipo margarita, un motor de corriente continua es el encargado de posicionar el carácter
correspondiente de la margarita antes de que se produzca la impresión en el papel.
Ejemplo 1.3
Depósito de líquido. En determinados procesos químicos es necesario mantener constante el
caudal de un líquido. Actualmente el control de caudal de un líquido se suele reducir al control
de nivel de líquido en un depósito, ya que un pequeño orificio en un depósito producirá un
caudal de salida constante si el nivel de líquido se mantiene constante y la ventaja de hacerlo
así es que dimensionando adecuadamente el depósito el caudal de salida puede ser
prácticamente insensible a pequeñas variaciones en el nivel. El mecanismo de control de nivel
más simple, pero aún ampliamente utilizado, es la válvula de flotador (típica boya de las
cisternas) que combina en el mismo dispositivo el sensor y el actuador.
La figura 1.5 nos muestra un depósito de líquido con un orificio de salida variable de forma
manual, según la posición del grifo y una tubería que le suministra un caudal de líquido
también variable (en este caso de forma electromecánica) mediante una válvula de control.
Dada la posición del grifo de salida se desea tener control del nivel del líquido.
Una posibilidad para conseguir este objetivo se indica en la figura, consistente en medir el
nivel (mediante el sensor de nivel, por ejemplo un electrodo medidor de capacidad),
compararlo con el nivel de referencia y realizar la acción de control abriendo o cerrando la
válvula que gobierna el caudal de entrada (mediante el controlador).
Ejemplo 1.4
Satélite artificial. Los satélites artificiales han supuesto un gran avance en el campo de las
telecomunicaciones terrestres y en el conocimiento del espacio. Un satélite artificial consta,
normalmente, de un cuerpo principal (que contiene toda la electrónica de comunicación, los
propulsores, etc..) y de dos brazos que contienen los paneles solares, las antenas y demás
sensores. El éxito de las misiones encomendadas a un satélite de comunicaciones depende de
que éste mantenga en todo momento una orientación adecuada de sus antenas y sensores
respecto a la tierra. Con este fin se incorpora en el satélite un sistema de control de la
orientación (“attitude” en la literatura inglesa) respecto a un sistema de referencia inercial.
La figura 1.6 nos presenta de forma esquemática una sección del satélite y su orientación
respecto a uno de los ejes del sistema inercial.
En primer lugar se deben dar las especificaciones de diseño. Éstas podrán venir
expresadas en el dominio del tiempo (ver capítulo 3), con condiciones como pueden ser
que el sistema tenga un tiempo de asentamiento determinado, una velocidad de respuesta
dada, una máxima sobreelongación ante un cambio escalón en la referencia determinado,
etc... También pueden venir expresadas en el dominio de la frecuencia, como pueden ser
ancho de banda, márgenes de ganancia y fase, etc... (ver capítulo dedicado a la respuesta
en frecuencia de los sistemas). La mayoría de estas especificaciones deben ser en muchas
ocasiones traducidas desde el lenguaje coloquial a alguna de las anteriores.
Una vez que se tienen claras las especificaciones el siguiente paso, y uno de los más
comprometidos, es establecer un modelo matemático del sistema. La complejidad de
éste tiene que estar en consonancia con los objetivos propuestos. Por ejemplo, puede no
hacer falta establecer el modelo de comportamiento térmico de un motor para controlar su
velocidad. En otras palabras, deben utilizarse el mayor número de hipótesis que
simplifiquen el modelo sin que se pierda información necesaria en la posterior fase de
diseño.
A partir del modelo linealizado, por ejemplo una función de transferencia obtenida
mediante transformada de Laplace de las ecuaciones diferenciales, comienza la fase de
diseño del controlador propiamente dicha. Generalmente se utilizan herramientas de
simulación para esta fase, ya que la mayoría del software de este tipo permite la inclusión
de modelos con este formato. Las ayudas gráficas de estos paquetes informáticos permiten
realizar sucesivos rediseños hasta que se cumplan las especificaciones dadas en la primera
fase del proceso.
Por supuesto, no en todo diseño es preciso pasar por todas las etapas anteriores, ya que
se pueden emplear experiencias previas para evitar etapas complicadas, como puede ser la de
modelado.
Especificaciones
de diseño
Obtención del
modelo
NO
Validación
Linealización
NO
Validación
Diseño del
controlador con
el modelo lineal
¿Cumple las NO
especificaciones
dadas?
Validación con el
modelo no lineal
¿Cumple las NO
especificaciones
dadas?
Validación con el
proceso real
Implementación
física
NO
¿FUNCIONA?
Como se ha visto en el último apartado del tema anterior, uno de los puntos críticos a la
hora de diseñar un controlador es obtener un buen modelo del sistema que se va a controlar. En
este punto, debe existir un compromiso entre la precisión del modelo y la utilidad de éste. En
ocasiones, puede ser más útil crear un modelo en un punto concreto de trabajo, haciendo
múltiples suposiciones que lo simplifiquen antes que obtener un modelo cualitativo del sistema
completo.
En principio, existen dos fuentes de conocimiento de las propiedades de los sistemas, que
aclaran un poco el compromiso expresado en el párrafo anterior. Una de estas fuentes es la
recopilación de la experiencia de los expertos y de la literatura del área en cuestión. Si un
ingeniero de control tiene que diseñar un controlador que satisfaga ciertos requerimientos de
calidad en el producto de una columna de destilación, y necesita hacer un modelo de ésta,
necesitará de información previa de este tipo de plantas, haciendo uso de la bibliografía existente.
Como el tipo de sistemas sobre el que se puede establecer criterios de control es prácticamente
infinito, no se puede pretender que el diseñador sea experto en cada uno de ellos. Dentro de esta
fuente se incluyen las leyes de la naturaleza en las que han estado trabajando generaciones de
científicos. El arco iris de Karplus, en la figura 2.1, nos da una idea de los modelos que se pueden
modelar con más exactitud utilizando este tipo de técnicas (y en general esto sería aplicable a
cualquier técnica).
Sistemas
Sistemas Sistemas
Biológicos
Económicos Químicos
Sistemas
Sistemas Mecánicos
Sociales Análisis
Predicción Control
Especulación Diseño
EDO’s
Sistemas EDif’s Circuitos
Psicológicos EDP’s
Eléctricos
EDO’s
EA’s
La otra fuente es el sistema por sí mismo. Observaciones del sistema y experimentos con
él son la base de todas las descripciones de sus propiedades. Se pueden excitar las entradas de los
sistemas, para después de analizar las salidas producidas, obtener modelos que produzcan la
misma respuesta ante los mismos estímulos. En este caso es mucho más importante hacer un
estudio del rango de validez de los modelos obtenidos.
Por tanto existen dos principios básicos y bastante diferentes para la construcción de
modelos, que pueden esquematizarse en la siguiente figura.
SISTEMA
F=ma
U=RI
MODELADO
IDENTIFICACIÓN
FÍSICO
MODELO
Identificación: el otro principio básico consiste en usar las observaciones de los sistemas
para conseguir fijar las propiedades del modelo a las del sistema. En general se usa como
complemento del método anterior. Entre sus características generales se pueden incluir las
siguientes:
Debe reflejar adecuadamente aquellas características del sistema que son de nuestro
interés.
Debe ser suficientemente sencillo como para resultar manejable.
A lo largo de este capítulo se van a exponer algunos ejemplos propios de la primera de las
metodologías comentadas anteriormente. Las técnicas de análisis y diseño de los sistemas
dinámicos caracterizados por ecuaciones diferenciales, requieren la descripción matemática del
sistema en cuestión. Para obtener las ecuaciones diferenciales que rigen el comportamiento del
sistema nos basaremos sencillamente en las leyes físicas de Newton para sistemas mecánicos, en
las leyes de Kirchhoff para sistemas eléctricos, etc.
Los sistemas que pueden representarse por los mismos modelos matemáticos, pero que
son físicamente diferentes, se denominan sistemas análogos. Si un sistema mecánico está
descrito por medio de las mismas ecuaciones diferenciales que uno eléctrico, aunque no tienen
nada que ver el uno con el otro desde un punto de vista físico, diremos que ambos son análogos.
A la derecha del arco iris de la figura 2.1 se encuentran los modelos eléctricos y
mecánicos, es decir, los modelos obtenidos a partir de las leyes físicas reproducen con bastante
exactitud las dinámicas de los sistemas que representan. Como representación se exponen algunos
ejemplos de este tipo de sistemas, comenzando por los mecánicos y continuando con los
eléctricos. Posteriormente se mostrará algún ejemplo de que ambos modelos representan a
sistemas análogos.
Ejemplo 2.1
m
X1
K C
X0
k = coeficiente de rigidez.
C = coeficiente de amortiguamiento.
Aplicando la segunda ley de Newton al sistema presente se obtiene la ecuación del modelo
considerado.
Ejemplo 2.2
m1 = masa de la rueda
m2
X2
Ks C
m1
X1
Kr
X0
x1 C x1 x2 K s x1 x2 K r x1 x0 m1 g 0
m1
x2 C x2 x1 K s x2 x1 m2 g 0
m2
Ejemplo 2.3
Supongamos que al carro se le aplica una fuerza u y que el centro de gravedad de la barra del
péndulo se encuentra en su punto medio. La masa del péndulo es m, y su longitud 2l, y la
masa del carro es M. Sabiendo que el péndulo se desplaza un ángulo de la línea vertical, y
que las coordenadas de enganche de la barra y el carro son (x, y), la posición del centro de
gravedad de la barra queda perfectamente definida:
xG x lsen
yG lcos
Si se considera el diagrama del cuerpo libre, el movimiento rotacional de la barra del péndulo
alrededor de su centro de gravedad se describe mediante:
d2
maxG FxG m x lsen H
dt 2
d2
ma yG FyG m 2 lcos V mg
dt
d 2x
MaxC FxC M 2 u H
dt
d2y
Ma yC FyC M 2 0
dt
I Vl Hl
m
x l H
0 V mg
Mx u H
M m x ml u
Estas dos ecuaciones son el modelo matemático lineal de nuestro sistema, ya que definen
perfectamente el comportamiento del mismo.
Junto con los anteriores, son los sistemas cuyos modelos reflejan con mayor precisión sus
dinámicas. Utilizando las conocidas leyes de Kirchoff, que indican que la suma algebraica de las
intensidades que entran y salen de un nudo es cero y que las sumas algebraicas de las tensiones
en una malla cerrada es nula, junto con las expresiones que relacionan tensión e intensidad en
cada uno de los componentes pasivos (resistencias, condensadores y bobinas), no resulta difícil
obtener expresiones como las mostradas a continuación.
Ejemplo 2.4
Considérese el circuito eléctrico que se muestra en la figura 2.6. Aplicando las leyes de
Kirchhoff a cada una de las mallas que configuran el circuito se obtiene las siguientes
ecuaciones diferenciales:
1
i1 t i2 t dt R1i1 t ei t
C1
1 1
i2 t i1 t dt R2i2 t i2 t dt 0
C1 C2
1
i2 t dt e0 t
C2
R1 R2
ei i1 C1 i2 C2 e0
Evidentemente son aquellos que utilizan elementos que convierten la energía mecánica en
eléctrica o viceversa. Los componentes más utilizados son los motores y generadores de corriente
continua, así como los motores bifásicos de corriente alterna. Por lo general, el motor más usado
es el de corriente continua, ya que en su eje se pueden generar potencias adecuadas para
multitud de aplicaciones. Dentro de estos motores existe una gran variedad de tipos, siendo el más
común el motor con excitación independiente. La figura 2.7 muestra el esquema eléctrico de este
tipo de motor, donde el inducido se ha modelado como un circuito eléctrico con una resistencia Ra
en serie con una inductancia La, y una fuente de tensión eb que representa la fuerza
contraelectromotriz generada en el inducido por el giro del rotor. La excitación se ha representado
con una resistencia Rf en serie con una inductancia Lf:
Ejemplo 2.5
Se trata de determinar un modelo matemático que represente la dinámica del comportamiento
del motor de corriente continua controlado mediante excitación independiente por inducido.
El esquema se mostró en la figura 2.7.
Según la figura, podemos controlar el motor de dos formas distintas: mediante la tensión ea
(por inducido), o mediante la tensión ef (por campo).
t k f i f t
También sabemos que el par desarrollado por el motor es proporcional al flujo magnético y a
la corriente en la armadura:
Tm km t ia t
En el control por inducido, ea es variable, pero ef se mantiene constante, lo que da lugar a que
if y, por tanto, también sean constantes. Ello da lugar a que el par motor solo dependa de ia:
Tm k1ia t
Siendo k1 km k f i f
eb t kb t
dia t
La Ra ia t eb t ea t
dt
d 2 t d t
Jm Bm Tl t Tm t
dt 2 dt
Para llegar al sistema de ecuaciones que definen el comportamiento del sistema sustituimos los
valores de eb y Tm determinados anteriormente:
dia t
La kb t Ra ia t ea t
dt
d 2 t d t
Jm 2
Bm k1ia t Tl t
dt dt
Ejemplo 2.6
Considérese un sistema tal como el de la figura 2.8. donde se presenta un tanque de agua.
Determine las ecuaciones diferenciales que rigen el comportamiento del sistema sabiendo que
la energía térmica se conserva en todo momento.
Qr Qa Qc Q f Q p
dTt
Qa c
dt
donde c es la capacidad térmica del agua del tanque, cuyo valor podría ponerse en función de
la masa m y del peso específico v.
- Flujos de calor existente en el agua del tanque que sale y en el agua fría que entra en el
tanque:
Qc VvTt
Q f VvT f
- Flujo de calor que pasa a través de las paredes del tanque, ya que se supone que éstas no son
adiabáticas y permiten el paso de calor a través de las mismas.
Tt Te
Qp
R
Una vez que se conocen todos los términos de la ecuación de partida, se sustituyen y se
obtiene el siguiente resultado:
dTt T T
Qr c Vv Tt T f t e
dt R
dTu 1
Qr c Vv Tu
dt R
Antes de comenzar con el ejemplo, sería conveniente repasar una serie de cuestiones
referentes a mecánica de fluidos. Un flujo en un conducto se considera laminar cuando su número
de Reynolds (este número adimensional relaciona la velocidad del flujo, las dimensiones del
conducto y la viscosidad del fluido) es menor de 2200. Cuando el número de Reynolds es mayor de
2200 se considera turbulento.
Laminar Q KH
Turbulento Q K H
Ejemplo 2.7
Se trata de hallar la ecuación diferencial que gobierna la altura de líquido en el depósito de área
A y volumen V (figura 2.9), suponiendo que la descarga es turbulenta, y que en el instante
inicial el depósito se encuentra vacío.
Qe
H
Válvula
Qs
A
dV d AH t dH
A Qe t Qs t
dt dt dt
Combinando ambas ecuaciones y teniendo en cuenta las condiciones iniciales nos queda:
dH t
A Qe t K H t
dt
H t 0 0
Por otro lado, la transformada permite también convertir funciones comunes en otras más
simples en el dominio de s, tales como exponenciales, senoides amortiguadas, etc... En definitiva,
funciones que aparecen con relativa frecuencia en el análisis de las ecuaciones diferenciales que
representan los sistemas dinámicos.
Una ventaja adicional que supone el empleo de la Transformada de Laplace es que nos
proporciona tanto el componente transitorio como estacionario de la respuesta temporal (solución
de la ecuación diferencial).
G s Gx jG y
G s Gx jG y
G s Gx2 G y2
Gy
atan
Gx
Se dice que una función compleja G(s) es analítica en una región si G(s) y todas sus
derivadas existen en dicha región. La derivada de una función analítica G(s), se determina
mediante:
dG s G s s G s G
lim lim
dt s 0 s s 0 s
Debido a que s = + j, el límite anterior puede tender a cero por infinidad de
caminos diferentes. Es demostrable que si se calcula la derivada por dos caminos distintos y ésta
coincide en ambos casos, entonces la derivada existe y, además, es única.
dG s Gx G y Gx G y
lim j j
ds 0
dG s Gx G y G G y
lim j j x
ds j 0
j j
Para que la derivada por los dos caminos seguidos sea única es necesario que se cumpla:
Gx G G G
j y j x y
Gx G y G y G
y x
Aquellos puntos del plano s donde la función G(s) sea analítica se denominan puntos
ordinarios. Los puntos donde G(s) no sea analítica se llaman puntos singulares y aquellos
puntos singulares donde la función o sus derivadas tiendan a infinito se les llama polos. Los
puntos en los que G(s) se anule se llaman ceros.
G s s p
n
con n 1, 2,3,... tiene un valor finito distinto de cero en el mismo punto,
entonces se dice que el polo es de orden n. Por ejemplo, dada la función de variable compleja:
K s 2 s 10
G s
s s 1 s 5 s 15
2
Hay que decir que si se incluyen los puntos en el infinito, toda función de variable compleja
tiene el mismo número de ceros que de polos. Por ejemplo, la función:
K
G s
s3
Teorema de Euler. Este teorema es una gran importancia puesto que relaciona las
funciones trigonométricas con las exponenciales, permitiendo resolver integrales en las que
aparecen estas funciones trigonométricas de forma más sencilla. Utilizando las expansiones en
serie de potencias:
2 4 6
cos 1 ...
2! 4! 6!
3 5 7
sen ...
3! 5! 7!
entonces:
j j j
2 3 4
x 2 x3
ex 1 x ...
2! 3!
se observa que:
cos jsen e j
El complejo conjugado del obtenido mediante el Teorema de Euler se puede escribir como:
cos jsen e j
e j e j
cos
e j cos jsen 2
j j j
e cos jsen sen e e
2j
L f t F s f t e st dt
0
1 c j
L1 F s f t F s e st ds, t 0
j 2 c j
donde c es una constante real que se elige mayor que las partes reales para todos los
puntos singulares de F(s). En consecuencia, la trayectoria de integración es paralela al eje
imaginario desplazada una cantidad c de él y, además, se encuentra a la derecha de todos los
puntos singulares de F(s).
Así pues la transformada de Laplace L actúa sobre cualquier función f(t) [f(t) = 0 t < 0]
para la que esa integral exista, y produce su transformada de Laplace L[f(t)] = F(s), una función
de la variable compleja s. Recordamos al lector que la integral impropia que aparece en la
definición está definida por el límite:
b
f t e st dt lim f t e st dt
b
0 0
y existe sólo cuando tal límite existe también, en cuyo caso se dice que la integral impropia
es convergente.
Linealidad: L f1 t f 2 t L f1 t L f 2 t F1 s F2 s
df t
Diferenciación real.: L sF s f 0
dt
d 2 f t 2
L 2 s F s sf 0 f 0
dt
......
d n f t n n2 n 1
L n s F s s n 1
f 0 s n2
f 0 ... s f 0 f 0
dt
F s f t dt t 0
L f t dt
s s
t F s
Caso de integrales definidas: L f t dt
0 s
t
Cambio de escala: L f aF as
a
dF s
L tf t
ds
d 2F s
L t 2 f t
ds
......
d nF s
L t n f t 1
n
, con n=1,2,3,...
ds
Teorema del valor final: Sólo es aplicable si el límite temporal cuando t existe:
lim f t lim sF s
t s 0
Teorema del valor inicial: a diferencia del anterior, este teorema es válido en todos los
casos:
f 0 lim sF s
s
t
L f1 t * f 2 t L f1 t f 2 d F1 s F2 s
0
f t 0 para t 0
f t A para t 0
f t 0 para t 0
f t At para t 0
Su transformada es:
A
L At Ate st dt
0
s2
f t 0 para t 0
f t Asen t para t 0
e j t e j t
Sabemos que sent se puede escribir como sen t
2j
Por lo tanto su transformada es:
A A 1 1 A
L Asen t
2j 0
e j t e j t e st dt 2
2 j s j s j s 2
De forma similar, la transformada de Laplace de Acos(t) queda:
As
L Acos t
s 2
2
f t 0 para t 0
f t Ae
t
para t 0
0 0
s
f t 0 para t 0, t0 t
A
f t t para 0 t t0
0
A A
Podemos poner f(t) como f t 1 t 1 t t0
t0 t0
La transformada nos queda:
A A A A st0 A
L f t L 1 t L 1 t t0 e
1 e st0
t0 t0 t0 s t0 s t0 s
Función impulso. La función impulso es un caso especial limitativo de la función pulso. Sea la
función impulso siguiente:
f t 0 para t 0, t0 t
A
f t lim
t0 0 t
para 0 t t0
0
A
L f t lim
t0 0 t s
1 e st0 A
0
El área cubierta bajo el impulso es igual a A. La función impulso cuya área es igual a la
unidad, recibe el nombre de función impulso unitario o delta de Dirac. La función impulso unitario
que se produce en el tiempo t = t0 se designa generalmente por (t-t0) y satisface las siguientes
condiciones:
t t0 0, t t0
t t 0 , t t 0
t t0 dt 1
-
d
t t0 1 t t0
dt
Tabla de transformadas.
En la siguiente tabla se pueden observar las transformadas de Laplace más comunes
utilizadas durante el resto de capítulos:
f(t) F(s)
t 1
1
1 t
s
1
t
s2
t n 1 1
n 1, 2,3,...
n 1! sn
n!
tn n 1, 2,3,...
s n 1
1
e at
sa
1
te at
s a
2
1 1
t n 1e at n 1, 2,3,...
n 1 ! s a
n
n!
t n e at
s a
n 1
sen t
s 2
2
s
cos t
s 2
2
senh t
s 2
2
s
cosh t
s 2
2
Tabla 2.1.
1 c j
L1 F s f t F s e st ds, t 0
j 2 c j
donde c, la abscisa de convergencia, es una constante real y elegida mayor que las partes
reales de todos los puntos singulares de F(s). Evaluar esta integral de inversión es complicado
pero, afortunadamente, existen procedimientos más simples que realizar esta integración
directamente, como es la descomposición en fracciones simples.
B s
F s
A s
Donde tanto A(s) como B(s) son polinomios en la variable s. Partiendo de esta situación,
para poder aplicar el procedimiento de expansión en fracciones parciales es necesario que la
potencia mayor de A(s) sea superior a la potencia mayor de B(s). En caso contrario, se precisará
dividir B(s) entre A(s) para obtener un polinomio en s además de un residuo al que se le podrá
aplicar el procedimiento (ya que en este residuo sí se cumplirá la condición). Además, también hay
que calcular las raíces del polinomio del denominador A(s).
F s F1 s F2 s ... Fn s
Las funciones temporales obtenidas por este procedimiento son únicas salvo en aquellos
puntos en los que las funciones temporales sean discontinuas. Si dichas funciones temporales son
continuas entonces la correspondencia entre la función temporal y su transformada es uno a uno.
Por otro lado, una ventaja adicional que presenta el procedimiento es que, tras aplicarlo,
las funciones de variable compleja s que aparecen son muy sencillas y su transformada inversa es
fácil de memorizar.
Esta descomposición en fracciones simples se realiza de manera distinta en función del tipo
de polos que presente la función de variable compleja sobre la que se va a trabajar. Los casos
posibles son:
B s K s z1 s z2 ... s zm
F s mn
A s s p1 s p2 ... s pn
B s a a a
F s 1 2 ... n
A s s p1 s p2 s pn
siendo ai constantes denominadas residuos y que se pueden obtener, para este caso, de la
siguiente manera:
ak s pk F s s p
k
a a a a
1 s pk 2 s pk ... k s pk ... n s pk ak
s p1 s p2 s pk s pn s pk
Si dos de los polos que aparecen en la descomposición son complejos conjugados, sus
residuos correspondientes también lo serán, con el consiguiente ahorro de cálculo que supone.
Una vez calculados todos los residuos obtendríamos para cada una de las fracciones
simples:
a
L1 k ak e pk t
s pk
Y, finalmente:
Ejemplo 2.8
s
F s
s 2 s 3
Dado que sólo aparecen polos simples y distintos podemos hacer la descomposición en
fracciones parciales descrita con anterioridad:
s a1 a
F s 2
s 2 s 3 s2 s3
a1 F s s 2 s 2 2
a2 F s s 3 s 3 3
En consecuencia:
2 1 3
f t L1 L 2e 2t 3e 3t t0
s 2 s 3
Ejemplo 2.9
s7
F s
s 1 s 2 2s 5
s 2 2 s 5 s 1 j 2 s 1 j 2
Se observa que las raíces de este polinomio son complejas conjugadas y distintas. En este caso,
a pesar de poder aplicarse el procedimiento general ya visto, es más recomendable realizar otro
tipo de expansión tendente a la búsqueda de funciones temporales que se corresponden con
senoides o cosenoides amortiguadas por exponenciales. Este tipo de descomposición se realiza
de la forma siguiente:
a bs c
F s 2
s 1 s 2s 5
3
a F s s 1 s 1
2
Para el resto de residuos no podemos aplicar las expresiones sino que es necesario resolver el
sistema de ecuaciones que se plantea al operar e igualar los numeradores de las dos funciones
F(s):
a s 2 2 s 5 bs c s 1 a b s 2 2a b c s 5a c s 7
3
a2
ab 0
3
2a b c 1 b
5a c 7 2
1
c 2
3 1 1 3s 1
F s 2
2 s 1 2 s 2s 5
La transformada inversa del primer término es sencilla y se corresponde, según las tablas de
transformadas, con una exponencial. En cuanto al segundo término, es necesario operar un
poco con él para reconocer la transformada de una función temporal conocida. En concreto,
vamos a tratar de llegar a senoides y cosenoides amortiguadas. Como ya es sabido:
L e at sen t
s a 2
2
sa
L e at cos t
s a 2
2
3s 1 3s 1 s 1 s 1 2 s 1 2
3 3 3 3 3 3
s 2 2 s 5 s 12 4 s 1 4 s 1 4 s 1 4 s 1 4
2 2 2 2
s 1 1 2
3 3e t cos 2t e t sen 2t
s 1 4 s 1 4
2 2
3
De esta manera hemos conseguido obtener las funciones senoide y cosenoide amortiguadas
por exponenciales. Si a esta función obtenida le añadimos la transformada inversa del primer
término de la descomposición en fracciones parciales, obtenemos:
3 1 1
f t e t 3e t cos 2t e t sen 2t e t 3 3cos 2t sen 2t t 0
2 2 2
B s B s a a2 an
F s 1 ...
A s s p n
s p s p 2
s p
n
an s p F s
n
s p1
d
an 1 s p F s
n
ds s p1
...
1 dj
j s p F s
n
an j
j ! ds s p1
...
1 d n 1
n 1 s p F s
n
a1
n 1! ds s p1
Ejemplo 2.10
5 s 2
Calcular la transformada inversa de F s
s 2 s 1 s 3
Dado que existe un polo múltiple en s = 0, así como dos polos simples en s = -1 y en s = -3,
luego la expansión en fracciones simples correspondiente es:
a1 a b b
F s 2 1 22
s 1 s 3 s s
5
a1 s 1 F s s 1
2
5
a2 s 3 F s s 3
18
d 25
b1 s 2 F s
ds s 0 9
10
b2 s 2 F s
s 0 3
Por lo tanto:
25 1 10 1 5 1 5 1 25 10 5 t 5 3t
F s f t t e e t0
9 s 3 s 2
2 s 1 18 s 3 9 3 2 18
Ejemplo 2.11
x 7 x 3 x 3,
Hallar la solución x(t) de la ecuación diferencial 2 x 0 3, x 0 3
L x t X s
L x t sX s x 0
x t s 2 X s sx 0 x 0
L
2 s 2 X s sx 0 x 0 7 sX s x 0 3 X s 3
s
2s 2 X s 6s 6 7 sX s 21 3 X s 3
s
X s s 2s 7 s 3 3 s 27 6s
2
6 s 2 27 s 3 6 s 2 27 s 3
X s
s 2 s 2 7 s 3 s s 3 s 1
2
a1 a a
X s 2 3
s s3 s 1
2
a1 sX s s 0 1
8
a2 s 3 X s s 3
5
1 18
a2 s X s
2 s 1 5
2
1 8 1 18 1 8 18 1 t
X s x t 1 e3t e 2 , t 0
s 5 s3 5 s 1 5 5
2
Ejemplo 2.12
A
sX s x 0 aX s
s 2 2
A bs 2 b 2 A
s a X s 2 2 b
s s2 2
bs 2 b 2 A bs 2 b 2 A
X s 2
s 2 s a s j s j s a
a1 b s c1
X s 12
s a s 2
b a 2 2 A
a1 s a X s s a
a2 2
a1 s 2 2 b1s c1 s a bs 2 b 2 A
a1 b1 b
ab1 c1 0
a ac b 2 A
2
1 1
En este sistema tan sólo dos ecuaciones son necesarias, puesto que el residuo a1 es conocido.
En consecuencia, eliminamos la 3ª ecuación que, a priori, parece la más larga de operar. Por lo
tanto, empleando las dos primeras ecuaciones:
b a 2 2 A A
b1 b a1 b
a
2 2
a 2
2
Aa
c1 ab1 2
a 2
b a 2 2 A 1 A sa
X s 2
a
2 2
s a a s 2
2 2
x t f1 t f 2 t
b a 2 2 A
f1 t e at
a2 2
A s a A s a
F2 s 2 2 2
2 2
a s 2
s a s s 2 2
2 2 2 2
A a
f2 t 2 2
sen t cos t
a
Por lo tanto:
b a 2 2 A A a
x t e at sen t cos t
a
2 2
a 2
2
Sea el sistema lineal invariante en el tiempo definido por las siguiente ecuación diferencial:
2.4.2. Ejemplos
Ejemplo 2.13
En el ejemplo 2.1 llegamos a la siguiente ecuación diferencial para el sistema simplificado de
una suspensión de automóvil:
Hay que recordar que, para calcular la función de transferencia entre la salida (en nuestro caso,
x1) y la excitación (x0) que estamos considerando, se suponen el resto de entradas nulas, así
como las condiciones iniciales.
Luego en este caso hay que anular la excitación exterior debida al peso, con lo que la ecuación
diferencial queda:
ms 2
Cs k X 1 s Cs k X 0 s 0
X1 s Cs k
2
X 0 s ms Cs k
Ejemplo 2.14
x1 C x1 x2 k s x1 x2 kn x1 x0 m1 g 0
m1
x2 C x2 x1 ks x2 x1 m2 g 0
m2
x1 C x1 x2 k s x1 x2 kn x1 x0 0
m1
x2 C x2 x1 ks x2 x1 0
m2
m1s 2 Cs k s kn X 1 s Cs k s X 2 s kn X 0 s 0
m2 s Cs k s X 2 s Cs k s X 1 s 0
2
m1s 2 Cs k s kn Cs k s X 1 s kn X 0 s
Cs k s m2 s 2 Cs ks X 2 s 0
m1s 2 Cs k s kn kn X 0 s
Cs k s 0
X 2 s
m1s 2 Cs k s kn Cs k s
Cs ks m2 s 2 Cs k s
kn Cs k s X 0 s
X 2 s
m s
1
2
Cs k s kn m2 s 2 Cs ks Cs k s
2
X2 s kn Cs ks
X 0 s m1s Cs k s kn m2 s 2 Cs k s Cs k s 2
2
Ejemplo 2.15
Para el circuito eléctrico visto en el ejemplo 2.3, las ecuaciones que se obtenían eran:
1
i1 t i2 t dt R1i1 t ei t
C1
1 1
i2 t i1 t dt R2i2 t i2 t dt 0
C1 C2
1
i2 t dt e0
C2
1
I1 s I 2 s R1 I1 s Ei s
C1 s
1 1
I 2 s I1 s R2 I 2 s I2 s 0
C1s C2 s
1
I 2 s E0 s
C2 s
Eliminando I1(s) e I2(s), se encuentra que la función de transferencia entre E0(s) y Ei(s), es:
E0 s 1
Ei s R1 R2C1C2 s R1C1s R2C2 R1C2 s 1
2
Ejemplo 2.16
Existen robots industriales que tienen una gran flexibilidad en los brazos aún con una carga
pesada en las pinzas. La figura muestra un modelo de sistema flexible con dos masas.
Encuéntrese la función de transferencia Y(s)/F(s).
x y
f
F(t)
M m
x y
k(y - x) k(y - x)
F(t)
M m
Tomando como sentido positivo el de la fuerza aplicada, se establecen las dos ecuaciones de
movimiento:
F f y x k y x Mx
f y x k y x mx
a las que, suponiendo condiciones iniciales nulas, se les aplica la transformada de Laplace para
obtener:^
F s fs Y s X s k Y s X s Ms 2 X s
fs Y s X s k Y s X s ms Y s
2
ms 2 fs k
F s fs k Y s
fs k
Ms 2 fs k Y s
Y s fs k
2
F s s Mms M m fs k
2
Ejemplo 2.17
Obtengamos la función de transferencia para el sistema del péndulo invertido, de forma que se
muestre la relación entre el ángulo de desplazamiento del péndulo y la fuerza externa u
aplicada al sistema. Aplicando la transformada de Laplace a las dos ecuaciones que resultaron,
se obtiene el siguiente resultado:^
En este sistema se consideró que U(s) era la entrada al sistema, (s) el ángulo que forma el
péndulo con la vertical (salida) y X(s) el desplazamiento del carro (salida), estando las tres
señales en el dominio complejo de Laplace. Al existir dos salidas y una entrada, se obtendrán
dos funciones de transferencia. Para obtener la relación entre el ángulo y la fuerza externa, se
despeja X(s) de una de las ecuaciones y se sustituye en la otra, de modo que quede eliminada.
Se despejará de la segunda ecuación:
U s mls 2 s
X s
M m s2
Sustituyendo en la primera de las ecuaciones del sistema y despejando hasta dejar la relación
(s)/U(s) se llegará a la siguiente función de transferencia:
U s mls 2 s
I ml s s mls M m s 2 mgl s
2 2 2
s 1
G s
U s M m I ml 2 m 2l 2
M m g s2
ml
Los circuitos electrónicos con amplificadores operacionales también pueden ser tratados
como sistemas eléctricos para hallar su función de transferencia. En el siguiente ejemplo muestra
cómo aplicando las características de los AA.OO. se puede obtener la función de transferencia del
sistema donde se encuentra integrado. En electrónica tiene gran aplicación el hecho de determinar
la función de transferencia de un circuito, ya que de ella se puede conocer para qué tipo de
señales se puede aplicar el circuito en función de la frecuencia mediante la representación de
Bode. Una vez representado el sistema se puede ver claramente la evolución de la amplitud y la
fase de la respuesta senoidal del circuito en un rango de frecuencias ante una entrada senoidal
(esto sólo es válido en el tipo de sistemas que nos ocupa, es decir, lineales e invariantes en el
tiempo).
Ejemplo 2.18
Se considera que la corriente que circula hacia el interior del A.O. es prácticamente nula, por
lo que las intensidades que recorren el circuito son:
i1 i2 i3
ei e d e e0 e0 e
C
R1 dt R2
A continuación, se aplica la conexión virtual del punto de tensión e’ a tierra, de forma que e’
= 0. Si se sustituye este valor en la expresión anterior se obtendrá el siguiente resultado:
ei de e
C 0 0
R1 dt R2
Ei s R2Cs 1
E0 s
R1 R2
E0 s R2 1
Ei s R1 R2Cs 1
Ejemplo 2.19
Mostrar que los sistemas de las figuras son análogos, probando que sus funciones de
transferencia son del mismo tipo.
E0 s X 0 s
Ei s X i s
A B
+ Z1 + 1 R1
Z1 s
1 1 R1C1s
C1s
I(s) R1
Z2
Ei(s) Eo(s) 1 1 R2C2 s
Z 2 s R2
C2 s C2 s
C
- -
La caída de tensión entre B y C viene dada por Eo(s) = I(s) Z2(s) y la caída de tensión entre A
y C por Ei(s) = I(s) ( Z1(s) + Z2(s)) , luego:
xi
f2 k2
a x0
f1
k1
Cuando hay un desplazamiento vertical del sistema, en el punto a concurren tres fuerzas: la
debida al rozamiento f2, la debida al muelle k2 y la debida al rozamiento f1. Y en el punto b,
dos fuerzas: la debida al rozamiento f1 y la debida al muelle k1.Tomando como sentido
positivo del desplazamiento el dibujado en la figura (hacia abajo), un desplazamiento (xo-xi) <
0 provoca los siguientes equilibrios de fuerzas:
f1 f2
1 s 1 s
X0 s k1 k2
Xi s f f f f f
s2 1 2 s 1 2 1 1
k1k2 k1 k2 k2
Donde se observa una perfecta analogía con la función de transferencia del sistema anterior,
con la relación entre variables:
R f
1
C
k
Ya se vio en el primer tema que los sistemas podían ser lineales o no lineales, dependiendo
del tipo de ecuaciones diferenciales que rigen su dinámica (modelo matemático). En Ingeniería de
Control es mucho más habitual encontrar métodos de diseño que traten con modelos lineales pero
incluso los sistemas más simples presentan no linealidades como pueden ser zonas muertas,
histéresis, saturaciones, etc. Incluso su misma dinámica puede estar representada por una
ecuación en sí misma no lineal, como se pudo ver en los modelos del depósito anteriores. En estos
casos resulta conveniente aplicar alguna linealización de nuestro modelo con el objeto de poder
aplicar dichos procedimientos de diseño. Para realizar esto, y como veremos a continuación, es
necesario elegir el punto de trabajo o de operación en el que se va a linealizar. Posteriormente se
comprobará, por ejemplo mediante simulación, el rango de validez de la linealización, puesto que
no siempre la linealización es práctica, bien por el modelo o por el punto de trabajo elegido (por
ejemplo, elegir un punto de trabajo próximo a una zona de saturación de un actuador conduce a
desavenencias entre el modelo linealizado y el modelo no lineal, o más en concreto desavenencias
con el sistema real al que representa).
Para poder aproximar el sistema no lineal, que podría estar representado por una función
como la de la figura 2.12, es necesario fijar el sistema en un punto de trabajo, es decir, suponer
que las variables se desvían muy poco dentro de una condición de operación. Supóngase un
sistema con una entrada x(t) y una salida y(t). La relación que une la entrada y la salida es una
función f. Supóngase ahora que el punto de operación fijado es (x,y)=(a,f(a)) =(a,b). Se puede
llegar a la una relación aproximada entre la entrada y la salida del sistema mediante un desarrollo
en serie de Taylor:
df 1 d2 f
y f x f a x a x a
2
...
dx xa 2! dx 2 xa
Como se dijo anteriormente, hay muy poca desviación de las variables dentro de la
condición de operación. Ello hace que los factores x – a se hagan despreciables a partir de los
términos de segundo orden. El resultado es evidente:
df
y f x f a b K x a
dx xa
y b K x a
Como puede apreciarse en esta última expresión, existe una clara relación entre y–b y x–a,
resultando así una relación lineal a través de K. Podemos decir que hemos transformado el modelo
no lineal en un modelo lineal.
Este procedimiento es válido para modelos de más de una variable. En el caso de una
función de dos variables la relación viene dada por la función z = f(x,y). Del mismo modo que para
el caso de una variable, se fijan ahora las dos variables en torno a una condición normal de trabajo
(x,y,z) = (a,b,f(a,b,)) = (a,b,c) y se desarrolla esta función en serie de Taylor en torno a dicho
punto. La función se transforma en:
f f
z f x, y f a , b x a y b
x x a y xa
y b y b
1 2 f 2 f 2 f
x a y b x a y b ...
2 2
2! x 2 xa y 2 xa xy x a
y b y b y b
Por el mismo motivo que en el caso de una variable despreciamos el término de segundo
orden y superiores. El resultado simplificado es la siguiente expresión:
f f
z f x, y f a , b x a y b c K1 x a K 2 y b
x xa y xa
y b y b
z c K1 x a K 2 y b
La linealización de los sistemas sólo es válida para un punto de operación en el que apenas
haya desviaciones. Cuando se produzcan grandes desviaciones, el modelo linealizado dejará de
tener validez, o tendremos que relinealizar el modelo para las nuevas condiciones de operación.
Ejemplo 2.20
En el ejemplo 2. 7, las ecuaciones que se obtuvieron para el depósito fueron:
dH t
Qs K H t A Qe t Qs t
dt
Vemos que se trata de una ecuación no lineal (la variable H aparece afectada del exponente ½)
y, por lo tanto, es necesario linealizarla para poder aplicar la transformada de Laplace. Vamos
a suponer un punto de equilibrio en torno al cual las magnitudes varían muy poco:
Qe t Qe qe t
Qs t Qs qs t
H t H h t
Aplicando Taylor:
dQs dQs 1 1 1 1
Qs Qs
dH
H H dH
2
K
H
2
K
H H H H H H H
1 1 K K
Qs Qs
2
K H H qs h K1 qs K1h
H 2 H 2 H
dH t
A Qe t Qs t
dt
d H h
A Q qe Q qs
dt
dh
A qe qs
dt
dh
Usando la relación linealizada para qs en esta última ecuación nos queda: A qe K1h .
dt
L h t H s
L qe t Qe s
AsH s K1 H s Qe s
H s 1
Qe s As K1
Ejemplo 2.21
a) Con los agujeros A y C abiertos y B tapado, escribir las ecuaciones del sistema para h1 y h2
en función del caudal de entrada.
c) Repetir los apartados (a) y (b) con el agujero A cerrado y los agujeros B y C abiertos.
a).
Al tapar el agujero B, los depósitos están comunicados por el agujero A (véase figura). En
estas condiciones el caudal entre el primer y segundo depósito QA(t) viene dado por:
QA t K g H 1 t H 3 t
Qc t K g H 2 t
dH1 t dH 2 t
D Q t QA t D QA t QC t
dt dt
donde se ha supuesto que los depósitos son uniformes y que tienen el mismo área D.
Sustituyendo las expresiones de QA(t) y QC(t) se obtienen las ecuaciones para H1(t) y H2(t) en
función del caudal de entrada Q(t) y de las características de los depósitos (Kg, H3 y D).
c)
Al tapar el agujero A los depósitos están comunicados por el agujero B (véase figura). En estas
condiciones el caudal entre el primer y segundo depósito QB(t) viene dado por:
QB t K g H1 t H 2 t
donde se ha supuesto que el orificio B tiene también el mismo tamaño que el A. El caudal de
salida QC(t) está dado por la misma expresión que en el apartado (a) y, por el teorema de
conservación de masas, se tiene:
dH1 t dH 2 t
D Q t QB t D QB t QC t
dt dt
La relación entre caudales y niveles de líquido en los depósitos es una relación no lineal pero
que se puede linealizar suponiendo que las variaciones de caudales y de niveles (q(t), qB(t),
q t Q t Q, qB t QB t QB , qC t QC t QC
h1 t H1 t H1 t h2 t H 2 t H 2 t
El estado de equilibrio se alcanza cuando los caudales de entrada y de salida son iguales para
cada depósito. Entonces los niveles no varían y están relacionados con el caudal de entrada
como sigue:
Q2
Q QB QC H1 2 H 2 H 2
K g2
QB t Q t
qB t h1 t B h2 t
H1 t H1 t H1 H 2 t H1 t H1
H 2 t H 2 H 2 t H 2
qB t K B h1 t h2 t
siendo:
Kg Kg K g2
KB
2 H1 H 2 2 H2 2Q
qC t K C h2 t
siendo:
QC t Kg
KC KB
H 2 t H 2 H2
2 t H 2
dh1 t
D q t qB t q t K B h1 t h2 t
dt
dh t
D 2 qB t qC t K B h1 t h2 t K C h2 t
dt
a las que se les puede aplicar la transformada de Laplace, con el siguiente resultado:
Q s KB H2 s
H1 s
Ds K B
H2 s KB
Q s D s D 2K B KC s K B KC
2 2
K B KC K B K g2 2K B KC 3
n 1,5
D D 2 DQ 2 n D 2
Este tipo de representaciones son muy útiles en Ingeniería de Control, donde se supone
que cada sistema tiene una dinámica, representada por su función de transferencia, que no afecta,
o dicho de otra manera, que no carga, a las de los sistemas a los que está conectado. Así, las
flechas de estas representaciones no indican un flujo energético sino de señales. Esto es posible
porque la mayoría de los sistemas se conectan con los demás utilizando señales eléctricas de
forma que a su salida presentan una impedancia muy pequeña y a su entrada una impedancia muy
alta (infinita para consideraciones prácticas). Esto se consigue con el uso de amplificadores
operacionales o dispositivos similares.
Dos conceptos muy usados en la representación con diagramas de bloques son los
siguientes:
Punto de suma. En relación con la figura 2.14, un círculo con una cruz constituye el
símbolo que indica la operación suma. El signo más o menos indica si la señal ha de sumarse o
restarse. Es imprescindible que las cantidades a sumar o restar tengan las mismas dimensiones y
unidades:
a a-b
+
-
Punto de
bifurcación
R(s) E(s) C(s)
+ G(s)
-
B(s)
H(s)
B( s)
Función de transferencia en lazo abierto : G ( s ) H ( s ) FTLA
E ( s)
C (s)
Función de transferencia directa: G ( s ) FTTD
E (s)
Para determinar la función de transferencia en lazo cerrado es necesario tan sólo tener
presente que:
C s G s R s H s C s
C s G s
FTLC
R s 1 G s H s
Una relación entre las tres expresiones anteriores de FTLC, FTLA y FTTD es la siguiente:
FTTD
FTLC
1 FTLA
En la figura 2.16, podemos ver un sistema en lazo cerrado sometido a una perturbación.
Cuando dos entradas (la señal de referencia y la perturbación) están presentes en un sistema
lineal, aplicando el principio de superposición, cada entrada puede tratarse independientemente de
la otra, siendo la salida total la suma de las salidas individuales.
N(s)
R(s) + C(s)
+ G1(s) + G2(s)
-
H(s)
CN s G2 s
N s 1 G1 s G2 s H s
CR s G1 s G2 s
R s 1 G1 s G2 s H s
Si se hace |G1(s)| muy grande y |H(s)| aproximadamente la unidad para las frecuencias de
interés, se pueden obtener aproximaciones como las dos siguientes. La primera muestra que el
efecto de la realimentación con las premisas anteriores consiste en atenuar los efectos de las
perturbaciones, más conocido como rechazo de perturbaciones:
CN s G2 s G2 s 1
0
N s 1 G1 s G2 s H s G1 s G2 s H s G1 s H s
CR s G1 s G2 s G1 s G2 s
1
R s 1 G1 s G2 s H s G1 s G2 s H s
En la figura 2.17 aparecen las reglas más importantes para simplificar diagramas de
bloques, más o menos complicados. Las reglas mostradas en esta figura, se pueden resumir como:
figura 2.9.
Capítulo 2 Regulación Automática
Junto con las funciones de transferencia de los componentes, permite representar las
relaciones causa-efecto a través del sistema.
B C C B
(1)
C
B B C
(2)
AG1 AG1G2 A AG1 AG1G2
G1 G2 G2 G1
(3)
A AG1 AG1G2 A AG1G2
G1 G2 G1 G2
(4)
A AG1 AG1+AG2 A AG1+AG2
G1 + G1+G2
+
AG2
G2
(5)
B
A AG A AG-B
AG-B A G
G + + G
- -
B
G 1 B
B G
(6)
(7)
A AG A AG
G G
AG
AG
G
(8)
A AG A AG
G G
A
1 A
AG
G
(9)
A A-B - A-B
+ +
-
A-B A A-B
B +
-
(10)
(11)
A B A AG1+AG2
G1 1
+ + G2 G1
- G2 -
G2
(12)
A B
+ G1 A G1 B
-
1 G1G2
G2
(13)
2.5.3 EJEMPLOS.
Veamos algunos ejemplos de simplificación de diagramas de bloques.
Ejemplo 2.22
Simplifique el diagrama de bloques que aparece en la figura 2.18, y obtenga la función de
transferencia de lazo cerrado C(s)/R(s).
H1
G1
R(s) + + C(s)
+ G2 + +
-
+ G3
-
H2
H1
R(s) + C(s)
+ G1+G2 +
-
G3-H2
H1
G1 G2
R(s) + C(s)
+ G1+G2 +
-
G3-H2
R(s) G1 G2 H1 C(s)
1 G1 G2 G3 H 2 1
G1 G2
R(s) G1 G2 H1 C(s)
1 G1 G2 G3 H 2
Ejemplo 2.23
Simplifiquemos el sistema que aparece en la figura a). Si se mueve el punto del lazo de
realimentación negativa que contiene H2 hacia fuera del lazo de realimentación positiva
contiene H1, obtenemos la figura. b).
H2
R(s) C(s)
-
+ + G1 + G2 G3
- +
H1
Figura a)
H2
G1
R(s) C(s)
-
+ + + G1 G2 G3
- +
H1
Figura b)
H2
G1
R(s) C(s)
- G1G2
+ + G3
- 1 G1G2 H1
Figura c)
R(s) C(s)
G1G2G3
+
- 1 G1G2 H1 G2G3 H 2
Figura d)
R(s) C(s)
G1G2G3
+
- 1 G1G2 H1 G2G3 H 2 G1G2G3
Figura e)
Ejemplo 2.24
El sistema representado en la figura 2.20 es un intercambiador de calor controlado por una
válvula que regula la entrada de fluido caliente según una señal obtenida en el regulador, la cual
es función del error existente entre la temperatura del fluido calentado y una de referencia.
Se trata de representar el diagrama de bloques del sistema así como la función de transferencia
entre la temperatura de referencia Tr y la temperatura de salida del fluido calentado T.
Ti Td
r q T
+
Regulador Válvula Intercambiador
Tr
-
Tt
Termómetro
e(t): error
Intercambiador:
qs(t): flujo calorífico que sale del cambiador con el fluido calentado.
qp: calor perdido por conducción a través de las paredes del intercambiador.
Con estas definiciones se pueden plantear las siguientes ecuaciones para el intercambiador:
q qs qe qM q p
qs m2CeT
qe m2CeTi
dT
qM M
dt
q p R T Td
La primera ecuación indica el balance de flujo calorífico en él. Las restantes son aplicaciones
de las leyes elementales de Termotecnia.
La obtención de la función de transferencia se puede hacer operando con las ecuaciones. Así
que, sustituyendo todas las ecuaciones en la primera de ellas, tenemos:
dT
q m2Ce(T Ti ) M R(T Td )
dt
1 m2 Ce R
T q T Td
Ms m2 Ce R Ms m2 Ce R Ms m 2 Ce R
Es decir, el sistema tiene, como ya se ha dicho, tres variables que actúan sobre T. Se
considerarán Ti y Td constantes, uniendo el intercambiador al resto del bucle de regulación por
medio de la variable q. Según la función de transferencia se podría representar el sistema por el
diagrama de bloques de la figura.
Td
+ T
q 1
+ Ms m2Ce R
+
m2Ce
Ti
Tt ( s ) k
T ( s) s k
r( s )
Regulador: Se supone que es de acción proporcional, es decir: kr
E( s )
x( s ) A
Válvula: Su función de transferencia será
r( s ) ms fs C
2
Con todas las expresiones obtenidas, se puede ya trazar el diagrama de bloques de la totalidad
del sistema de regulación.
Td
R
Tr + E r x m1 q + T
A 1
Kr Kp Kq
ms 2 fs C + Ms m2Ce R
- +
m2Ce
Ti
K
sK
Ejemplo 2.25
Encontrar la función de transferencia C(s)/R(s) para el diagrama de bloques de la figura.
G1(s)
+ + C(s)
R(s)
G2(s) G3(s)
+
- -
H1(s)
H2(s)
Partiendo del diagrama original, en la figura (a) se ha descompuesto el punto suma, en (b) se
ha desplazado a la derecha el punto diferencia, en (c) se han combinado los bloques G1 y G2 y
se ha reducido el lazo de realimentación H1, en (d) se ha reducido el lazo de realimentación G2
H2.
Por último, combinando los dos bloques en serie de la figura (d), se obtiene la función de
transferencia del sistema.
C s G1 s G2 s G3 s
R s 1 G3 s H1 s G2 s G3 s H 2 s
G1(s)
+
R(s) + C(s)
+ +
G2(s) G3(s)
- -
H1(s)
H2(s)
Figura (a)
G1(s)
+
R(s) C(s)
+ + +
G2(s) G3(s)
- -
H1(s)
G2(s)H2(s)
Figura (b)
R(s) C(s)
+ G3 s
G1(s)+G2(s)
1 G3 s H1 s
-
G2(s)H2(s)
Figura (c)
R(s) G3 s C(s)
G1(s)+G2(s)
1 G3 s H 1 s G2 s G3 s H 2 s
Figura (d)
sistemas condicionalmente estables. Por otro lado, el error estacionario, es una medida de la
exactitud del sistema de control, es decir, es el error que aparece cuando la salida del sistema en
estado estacionario no coincide exactamente con la entrada. En este capítulo se analizará para los
distintos tipos de sistemas.
C s K
R s s 1
1
Al tratarse de una función escalón su transformada de Laplace viene dada por R( s)
s
K 1
La transformada de la respuesta queda C ( s )
s 1 s
1
Al expandir C(s) en fracciones simples C ( s) K
s s 1
t
Tomando la transformada inversa queda c (t ) K 1 e
(t 0)
Si el escalón de entrada no tuviera amplitud unitaria, sino que tuviera una amplitud n
obtendríamos:
K n n n
C s . K
s 1 s s s 1
t
c t Kn 1 e
dc(t ) K t K
e
dt t 0 t 0
En t 4 se alcanza el 98.2% del valor final, por lo que este valor puede considerarse
como el tiempo de asentamiento, es decir, aquél en el que la salida alcanza la
referencia. En realidad, matemáticamente sólo se alcanzará el valor estable cuando
t pero se suele considerar como criterio para la estabilización cuando toda la salida
se encuentra dentro de una franja del 2% alrededor de dicho valor estable.
c(t)
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0 1 2 3 4 5 6
En vista de los resultados anteriores, cuanto menor sea , es decir, cuanto más a la
izquierda esté situado el polo en el plano complejo, la salida alcanzará antes a la entrada o, lo
que es lo mismo, tendrá un comportamiento más rápido o agresivo. Por lo tanto, la constante de
tiempo de un sistema de primer orden es una medida de la velocidad de respuesta del mismo.
1
La transformada de Laplace de la entrada viene dada por R ( s )
s2
K 1
La transformada de la respuesta queda C ( s ) 2
s 1 s
1 2
Al expandir C(s) en fracciones simples tenemos C ( s ) K 2
s s s 1
t
Y tomando la transformada inversa, c (t ) K t e (t 0)
K n n n n 2
C ( s) 2 K 2
s 1 s s s s 1
t
c t Kn t e
De esta expresión puede comprobarse que el error estacionario (cuando t) si K=1
equivale a . Si K1 entonces el error estacionario valdrá + (0<K<1) ó - (K>1).
Evidentemente, para K=1, cuanto menor sea menos error se cometerá en la salida. En la
figura 3.2 aparecen representadas la entrada r(t), la salida c(t) y el error e(t) en estado
estacionario, cuando K=1:
10
8
7 r(t)
6
c(t)
5
4
pte = K
3
0
0 1 2 3 4 5 6 7 8 9 10
t/
C (s) n2
R( s) s 2 2 n s n2
donde:
n = frecuencia natural no amortiguada.
= factor de amortiguamiento del sistema.
n = = factor de atenuación.
Ejemplo 3.1
Un ejemplo típico de sistema de segundo orden es el servomecanismo de posición angular
de la carga mecánica colocada en el eje del motor de corriente continua controlado por
inducido, que ya vimos, con determinadas simplificaciones, en temas anteriores. Este proceso
es de segundo orden, pero inestable, debido al polo de origen. Con la incorporación, en el lazo
de control, del controlador proporcional (que genera una señal de control proporcional a la
señal de error) de función de transferencia Kc, se pretende estabilizar el sistema y poder llevar
la carga mecánica a cualquier posición angular. Su esquema se representa en la figura 3.3
r(s) + (s)
K
Kp
s s 1
-
KK c
Función de transferencia:
s s KK c
2
Para hallar los polos lo único que hay que hacer es resolver la ecuación característica del
sistema:
s 2 2 n s n2 0
Los polos son, por lo tanto:
s1,2 n j n 1 2 j d
j
Raíz
n
cos n 1 2
d
n
En ésta gráfica se aprecian las relaciones entre los distintos parámetros, entre las que
destaca que = cos .
Vamos a realizar una clasificación general de los sistemas de segundo orden, atendiendo al
valor que toma :
j
s1 n j n 1 2
s2 n j n 1 2
j
s1 = s2 = n
En sistemas en los que s1 s2, se toma un único polo doble haciendo que el sistema se
comporte con amortiguamiento crítico.
j
s1 n n 2 1
s2 n n 2 1
j
s1 j n
s2 j n
5) Sistemas inestables:
j
s1 n j n 1 2
s2 n j n 1 2
5 b) ( = 1):
Polo doble real y positivo:
j
s1 = s2 = n
5 c) ( < 1):
Polos reales distintos y positivos:
j
s1 n n 2 1
s2 n n 2 1
n2 n2 n2
C ( s)
( s 2 2 n s n2 ) s s n j d s n j d s ( s n ) 2 d2 s
1 s n n
C ( s)
s ( s n ) d ( s n ) 2 d2
2 2
s n nt
L1 2
e cos d t
( s n ) d
2
n n 1 d
L1 2
L 2
e nt sen d t
( s n ) 2
d d ( s n ) 2
d 1 2
c(t ) 1 e nt cos d t sen d t o también:
1 2
e n t 1 2
c (t ) 1 sen d t siendo arctg (t 0)
1 2
e n t
e (t ) r (t ) c (t ) sen d t
1 2
Esta señal presenta una frecuencia de oscilación igual también a d, y en estado
estacionario:
lim e(t ) 0
t
0.25
0.2
c(t)
0.15
0.1
0.05
0
0 1 2 3 4 5 6 t
= 0
1.8
1.4
= 0.4
1.2
= 0.6
1
0.8 = 1
0.6
= 0.8
0.4
0.2
0
0 5 10 t
n2 1 1 n
La transformada de la salida queda C ( s)
( s n ) s s s n s n 2
2
c(t ) 1 e nt (1 n t )
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0 5 10 15
t
Figura 3.6. Respuesta temporal de un sistema de 2º orden con amortiguamiento crítico a un salto
escalón
En este caso se obtenían polos complejos puros. Para obtener la respuesta temporal tan
sólo es necesario hacer = 0 en la ecuación general de C(s) y se obtiene:
n2 1 s
C ( s) 2 2 y c(t ) 1 cos n t (t 0)
( s n ) s s s n2
2
1.8
1.6
1.4
1.2
0.8
0.6
0.4
0.2
0
0 2 4 6 8 10 12
t
En este caso los polos de la ecuación característica en lazo cerrado son reales negativos y
diferentes. Ante una entrada escalón la transformada de la salida C(s) es:
n2
C ( s)
( s n n 2 1)( s n n 2 1) s
1 2 1) n t 1 2 1) n t
c(t) 1 e ( e ( (t0)
2 2 1( 2 1) 2 2 1( 2 1)
n e s1t e s2t
c(t ) 1
2 2 1 s1 s2
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0 5 10 15 20 25
t
Cuando es bastante mayor que la unidad quiere decir que s1 s2 y, por lo tanto,
podemos despreciar s1 en la solución ya que el término exponencial e-s1t tiende mucho antes a cero
que el término exponencial e-s2t y quedarnos sólo con s2, convirtiendo la respuesta en una
aproximada a la de un sistema de primer orden:
C ( s ) n2 1
R( s ) s1 s s2
n2
y la respuesta temporal queda c (t )
s1s2
1 e s2t
En este caso, las raíces se encuentran en el semiplano derecho, es decir, la parte real de
las raíces es positiva, lo que representa términos exponenciales crecientes en la respuesta del
sistema, como puede verse en las figuras incluidas en la tabla y, por tanto, la respuesta crece
indefinidamente.
CLASIFICACIÓN SEGÚN
Im(s)
Re(S)
0 1
e n t 1 2
c (t ) 1 sen d t siendo arctg
1 2
SISTEMA CON
s1,2 n
AMORTIGUAMIENTO CRITICO
Im(s)
Re(S)
1
1 1 1 1 n
L1 n
L 2
1 e nt 1 nt
n s s n s n
2
s s
Im(s)
Re(S)
1
c(t ) L1
n 2
s s s1 s s 2
1
e n t
2 2 1
2 1 e n 2 1 t
2 1 e n 2 1 t
CLASIFICACIÓN SEGÚN
SISTEMA SIN
s1,2 j n
AMORTIGUAMIENTO
Im(s)
=0
Re(s)
n2 1
1 cos n t
1
L 2
s n s
2
Im(s)
-1<<0
Re(s)
Im(s)
<-1
C(t)
Tolerancia
1
Mp
0.9
td
0.5
0.1
0 tr
tp t
ts
Tiempo de retardo, td: Es el tiempo que tarda la respuesta en alcanzar por primera vez
el 50% de su valor final.
Tiempo de subida, tr: Se define como el tiempo necesario para que la respuesta del
sistema pase del 10 al 90% de su valor final. En algunos casos, concretamente en sistemas
de segundo orden subamortiguados, se define como el tiempo entre 0 y 100%.
Tiempo de pico, tp: Es el tiempo requerido por la respuesta para alcanzar la máxima
sobreelongación.
Máxima sobreelongación, Mp: Se define como la máxima desviación de la respuesta
transitoria respecto a su valor estacionario. Se define siempre en tanto por ciento cuando
el valor estacionario es distinto de la unidad.
c(t p ) c()
M p (%) 100
c ( )
Tiempo de asentamiento, ts: Se define como el tiempo necesario para que la
respuesta alcance y se mantenga dentro de un determinado rango (habitualmente el
5% o el 2% de su valor final). En los sistemas de primer orden ya se observó qué se
entendía por este concepto y equivalía aproximadamente a cuatro o cinco veces la
constante de tiempo.
e nt
c(t ) 1 sen d t
1 2
e n td
0.5 sen d t d
1 2
ntd 1 0.7
luego:
1 0.7
td (0 1)
n
e n t r
c (t r ) 1 1 sen d t r
1 2
d t r n ( n 0,1, 2...)
Como tr debe ser positivo y corresponder al primer instante en que c(t)=1, se tiene:
tr
d
0.8 2.5
tr (0 1)
n
Para tp, la función c(t) tiene un máximo. Luego derivando c(t) nos queda:
n e nt
d e nt
sen d t cos d t 0
1 2 1 2
1 2
tg( d t ) tg
Por tanto, dt = n (n = 0,1,2...) son las infinitas soluciones, es decir los instantes en los
que se producen los máximos y mínimos, como se aprecia en la figura 3.10.
R e s p u e s ta te m p o ra l
C (t)
M a x s o b re e lo n g a c ió n
1 .0
M in .
0
2 3 4 n t
1 2
1 2
1 2
1 2
tp (0 1)
d n 1 2
Máxima sobreelongación, Mp
1 2
e
c (t p ) 1 sen
1 2
Como sen( ) sen 1 , operando nos queda c ( t p ) 1 e 1 2
2
1 2
M p (% ) e 100%
sobreelongación”. Para valores de 0<<0.6 se puede aproximar como M p 100 1 , que
0.6
también se representa en la misma gráfica.
Mp 100
90
Código MATLAB
80
70
clf;
clear;
60 clc;
50
indice=1;
% Dibuja Mp frente a delta
40 delta=[0:0.01:1];
30 Mp=exp(-delta*pi./sqrt(1-
delta.^2));
20
plot(delta,100*Mp,'b');
10 % Dibuja Mp frente a delta
aprox.
0
0 0 .1 0 .2 0 .3 0 .4 0 .5 0 .6 0 .7 0 .8 0 .9 1 delta l=[0:0.2:0.6];
Tiempo de asentamiento, ts
La expresión exacta es muy difícil de calcular. Sin embargo, se puede obtener una
expresión aproximada, haciendo uso en el caso subamortiguado (0 < < 1) de las curvas 1exp(-
2
nt/ 1 ). Estas curvas son las envolventes de la sinusoide amortiguada correspondiente a la
R e s p u e s ta te m p o r a l
1 .8
1 .6
e δ ωn t
1
1 .4 1 δ2
1 .2
Amplitude
0 .8
0 .6
e δ ωn t
0 .4 1
1 δ2
A p r o x . ts
0 .2 ts
0
0 10 20 30 40 50 60
T im e ( s e c .)
e n ts
1 1.05 y, despejando ts:
1 2
1
ts ln(0.05 1 2 )
n
Ejemplo 3.2
Considere el sistema de la figura. Obtengamos los valores de y n y los tiempos de subida tr,
el tiempo de pico tp, el sobrepaso máximo Mp y el tiempo de asentamiento ts, cuando el sistema
esta sujeto a una entrada escalón unitario.
25
s 6s 25
2
d n 1 2 4
A partir de los valores dados de y n obtenemos :
n 3
β 3.14
Tiempo de Subida tr: t r
d ω 4
4
donde se obtiene mediante : tan 1 tan 1 0.93 rad
3
3.14 0.93
Por tanto, el tiempo de subida tr es : tr 0.55 seg
4
3.14
Tiempo de pico tp : t p 0.785 seg
d 4
π
d
Sobrepaso máximo Mp: M p e ꞏ100 9,5%
Tiempo de establecimiento ts: Tanto para el criterio del 2% como para el de 5% tomaremos la
consideración t s 1
Ejemplo 3.3
A) Analizar el siguiente sistema en respuesta a una entrada escalón para los siguientes valores
de la ganancia K: 1) K= 0.02 2) K= 2.5 . Comprobar para los dos casos el tp, ts y Mp.
B) ¿ Cuál debe ser el valor de K para que se cumplan las condiciones de los siguientes casos?
R(s) + 10 C(s)
K
s s 1
-
Solución:
10 K
s(s 1) 10K
FTLC 2
10 K s s 10 K
1
s(s 1)
1
n 0.2 0.447 1.118 d 0.223 0.5
2 0.447
s1,2 0.5 0.223 0.273 , 0.723
Dado que > 1 no hay tp ni Mp. Tampoco nos sirve la expresión obtenida anteriormente de ts
puesto que sólo es para sistemas subamortiguados.
Respuesta temporal
-0.723 -0.273
A.2) K=2.5
1
25 5 0.1 0.5 d 4.975 s1,2 0.5 4.975 j
n 10
tp 0.63 , t s 2π , M p 100 e d 73%
d
-0.5+j4.975
-0.5-j4.975
B.1)
tp d 1 rad/seg n 1 2 1
d
0.5
t s 2 0.5 n 0.5 , n
0.5
1 2 1 0.447 n 1.118
10 K n2 K 0.125
B.2)
tp d 1 rad/seg
d
1 2
M p 20.8 % e 0.208 1.5702
1 2
2 0.1998 0.447
n 1 2 1 n 1.118 10 K n2 K 0.125
Ejemplo 3.4
Cuando el sistema de la figura esta sujeto a una entrada escalón unitario la salida del sistema
responde como se aprecia en la gráfica temporal. Determine los valores de K y T sabiendo que
el sobrepaso máximo es de 25,4 %, para un tiempo de 3 seg.
R(s) + K C(s)
s Ts 1
-
Solución.
tp 3
d n 1 0.42
C s K
A partir del diagrama de bloques tenemos que 2
R s Ts s K
de la cual:
K 1
n 2 n
T T
T = 1.09 K = 1.42
Ejemplo 3.5
Determine los valores de K y Kh del sistema en lazo cerrado de la figura para que el sobrepaso
máximo de la respuesta escalón unitario sea del 25 % y el tiempo de pico sea de 2 seg.
Suponga que J= 1 kg m2 y B=0.
+ + K 1
Js b s
- -
Kh
K
FTLC J
B KK K
s2 h
s
J J
Observe que
K KK h B
n 2 n
J J
1 2
M p e 0.25
que se especifica como del 25 %. Por tanto, despejando , su valor será = 0.404.
tp 2 , d 1.57
d
d
n 1.72
1 2
K n 2 2.95 Nm
Por tanto, obtenemos: 2 n
Kh 0.471 seg
K
a n ai
C (s)
s i 1 s pi
donde ai es el residuo en el polo s = -pi.
Lo habitual es que parte de los polos sean complejos, con lo cual podemos reescribir la
ecuación anterior como:
a q aj r
bk ( s k k ) ck k 1 k2 nq
C ( s) r
s j 1 s p j k 1 s 2 2 k k s k2 2
La respuesta temporal será pues:
q r r
c(t ) a a j e bk e kk t cos k 1 k2 t ck e kk t sen k 1 k2 t
p jt
(t 0)
j 1 k 1 k 1
Por lo tanto, la respuesta de un sistema de orden superior se puede escribir como suma de
funciones simples halladas en la respuesta de sistemas de primer y segundo orden. El tipo de
respuesta transitoria esta determinado por los polos en lazo cerrado, mientras que la forma
concreta de dicha respuesta depende principalmente de los ceros en lazo cerrado. Los polos de la
entrada R(s) producen los términos de la respuesta en estado estacionario en la solución, c() =
a, mientras que los polos de C(s)/R(s) entran en los términos exponenciales de la respuesta
transitoria, y/o en los términos senoidales amortiguados de la respuesta transitoria. Los ceros de
C(s)/R(s) no afectan los exponentes de los términos exponenciales pero sí las magnitudes y signos
de los residuos.
Ejemplo 3.6
Calcular el tiempo de asentamiento del sistema ante una entrada escalón.
+ 0.05 1
s 0.4 s3
-
1
s2
0.05 s 2 0.05 s 2
FTLC
s 3 s 2
2.4 s 0.85 s 3 s 0.43 s 1.97
Podemos considerar que, prácticamente, se puede cancelar (s+2) con (s+1.97) y que como la
parte real del polo s= -3 es menor que 6 veces la parte real del polo dominante lo podemos
descartar de nuestra función, dado que su influencia será mínima. Hay que hacer notar que se
cancelan los polos pero hay que incluir la contribución de ambos a la ganancia total en estado
estacionario.
Respuesta gráfica del sistema con los polos cancelados y sin cancelar:
Step Response
0.04
Con cancelaciones
0.035
0.03
0.025
Amplitud
0.015
0.01
0.005
0
0 5 10 15
Tiempo
Ejemplo 3.7
Atendiendo a la función de trasferencia considera los siguientes casos para C y analiza el
sistema.
5
G(s)
s C s s 4.25
2
Solución.
Obtenemos las raíces del polinomio de segundo grado: s1,2= -0.5 j 0.2.
5 1
FTLC 2
5 s s 4.25 s s 4.25
2
d
0.5 Mp e 100 45%
d 2
tp 1.57 seg ts 6.28seg
d
Respuesta gráfica :
Step Response
0.35
0.3
0.25
0.2
Amplitude
0.15
0.1
0.05
0
0 2 4 6 8 10 12
Tim e (sec.)
Respuesta gráfica:
S te p R e s p o n s e
0 .2 5
0 .2
0 .1 5
Amplitude
0 .1
0 .0 5
0
0 1 2 3 4 5 6 7 8
T im e ( s e c . )
C = 0.1. No se pueden despreciar los polos conjugados pues p<6 pero el sistema se parece
más a uno de primer orden. Si se eliminase, el sistema quedaría:
5 5 1.1765
G(s)
s 0.1 ( s 2
s 4.25) ( s 0.1)4.25 s 0.1
1
T 10 t 30 seg .
0.1 s
Comparemos las respuestas temporales con y sin cancelación ante un salto escalón en la
entrada.
Step Response
12
10
8
Amplitude
0
0 10 20 30 40 50 60
Time (sec.)
Ejemplo 3.8
El control de grandes procesos químicos es una parte importante del campo de la ingeniería
de control. En la figura se muestra el diagrama de bloques para el circuito de control de
temperatura en un proceso químico de Xyleno. Las funciones de transferencia para los
componentes son:
0.03 2
Proceso: G p s Sensor: Gs s
s 0.25s 0.01
2
3s 1
K s 2 s 0.1 2
Regulador: Gc s Válvula: Gv s
s 2s 1
b) Obtener una expresión aproximada para la respuesta del sistema a un cambio escalón
unitario en la referencia, cuando K=0.2.
D(s)
+
R(s) + C(s)
Gc(s) Gv(s) Gp(s)
+
-
Gs(s)
Solución.
Gc s Gv s G p s Gp s
C s R s D s
1 Gc s Gv s G p s Gs s 1 Gc s Gv s G p s Gs s
C s 3s 1 K s 2 s 0.1ꞏ2ꞏ0.03
R s 3s 1 s 2 s 1 s 2 0.25s 0.01 K s 2 s 0.1ꞏ2ꞏ0.03ꞏ2
-0.4602±0.0705j
-0.1308
-0.0160±0.1177j
todos ellos en el semiplano s izquierdo, por lo que el sistema es estable y donde se aprecia una
pareja de polos complejos muy próximos al eje imaginario que, por tanto, son los polos
dominantes del sistema.
Se trata ahora de aproximar la función utilizando únicamente el tercer término del desarrollo,
que va asociado a los polos dominantes. De esta forma se aproxima el régimen transitorio
pero falta por conseguir el mismo régimen estacionario. Para ello es preciso hacer un ajuste de
ganancias como sigue:
siendo:
C s 0.0002
kcc funcion completa lim 0.5
s 0 R s 0.0004
0.0044s 0.0077 0.0077
kcc termino del desarrollo lim 0.5461
s 0 s 0.0321s 0.0141 0.0141
2
El resultado es:
Para determinar la transformada inversa de esta expresión se puede hacer uso de los pares de
transformadas, con la correspondencia:
c t 0.2837ꞏ0.1187ꞏ1.0093ꞏe 0.0161t sen 0.1177t 0.5 1 1.0093ꞏe0.0161t sen 0.1177t 1.4352
En la figura están representadas las respuestas del sistema completo y del sistema aproximado,
para un intervalo de 100 seg. La repuesta aproximada presenta sobreelongaciones ligeramente
menores.
0.8 completa
aproximada
0.6
c(t)
0.4
0.2
0
0 20 40 60 80 100
t (sg.)
B(s) H(s)
Los coeficientes de error estáticos definidos a continuación son cifras de mérito de los
sistemas de control. En adelante, se denominará posición a la salida, velocidad al ritmo de
variación de la salida, etc.
El error estacionario del sistema, para una entrada de escalón unitario, es:
s 1
ess lim
s 0 1 G ( s) H ( s) s
La constante Kp de error estático de posición se define como:
K p lim G ( s ) H ( s )
s 0
El error estacionario del sistema, para una entrada de rampa unitaria, es:
s 1
ess lim
s 0 1 G ( s) H ( s) s 2
La constante Kv de error estático de velocidad se define como:
K p lim sG ( s ) H ( s )
s 0
1
ess Tipo 0
Kv
1 1
ess Tipo 1
Kv K
1
ess 0 Tipo 2 o mayores
Kv
El error estacionario del sistema, para una entrada parabólica unitaria (entrada
aceleración), definida por
t2
r (t ) para t 0
2
0 para t 0
está dado por:
s 1
ess lim
s 0 1 G ( s ) H ( s ) s 3
K p lim s 2G ( s ) H ( s )
s 0
1
ess
Ka
Para un sistema tipo 0:
s 2 K ( a s 1)( b s 1)
K a lim 0
s 0 ( 1s 1)( 2 s 1)
Para un sistema tipo 1:
s 2 K ( a s 1)( b s 1)
K a lim 0
s 0 s ( 1s 1)( 2 s 1)
Para un sistema tipo 2:
s 2 K ( a s 1)( b s 1)
K a lim 2 K
s 0 s ( s 1)( s 1)
1 2
s 2 K ( a s 1)( b s 1)
K a lim N ( N 3)
s 0 s ( s 1)( s 1)
1 2
Luego, para una entrada escalón unitario, el error estacionario ess queda:
1
ess Tipo 0 y tipo 1
Ka
1 1
ess Tipo 2
Ka K
1
ess 0 Tipo 3 o mayores
Ka
A continuación, se muestra la tabla 3.1, con todos los tipos de errores:
Tipo 1 0 1
K
Tipo 2 0 0 1
K
Ejemplo 3.9
En la siguiente figura se muestra el diagrama de bloques del control automático de freno de un
tren de alta velocidad, donde:
Y donde la respuesta del freno a una entrada escalón unitario, despreciando las fuerzas de
fricción, obedece a la ecuación f(t) = 100 ( 1 - e-10t ).
c) Si la velocidad del tren se quiere mantener de forma estable a 60 m/s, ¿cuál debe ser el
valor de Vr?
Vr + e eb f(t) v(t)
Amplificador
Freno Tren
K
-
et
Tacómetro
Kt
Solución.
a) El bloque tren tiene una salida, que es la velocidad v en m/s, y una entrada, f(t), que es la
fuerza ejercida por el freno. La ecuación que relaciona ambas variables es f(t) = M dv(t)/dt,
siendo M la masa del tren. De la que, aplicando transformada de Laplace suponiendo
condiciones iniciales nulas, se obtiene la función de transferencia para el tren:
V s 1 1
F s Ms 25000 s
1 1 1000
F s 100
s s 10 s s 10
y, teniendo en cuenta que está función es resultado de un escalón unitario en la entrada del
freno eb(t), la función de transferencia que describe el freno es:
F s 1000
Eb s s s 10
1000 1
V s 100
s 10 25000s 4
2
Vr s 1 100 1000 1
ꞏ0.15 s 10s 0.6
s 10 25000s
c) El sistema tiene sus dos polos, -0.0604 y -9.9396, en el semiplano s izquierdo, luego es
estable. Tiene una ganancia en estado estacionario de:
V s 4
kcc lim
s 0 Vr s 0.6
Por tanto, si se quiere mantener una velocidad del tren de 60 m/s, se necesita una tensión Vr
de 60/kcc = 9 V.
Ejemplo 3.10
Un objetivo importante en el proceso de manufactura del papel es mantener una consistencia
uniforme en la pulpa de salida a medida que ésta pasa al secado y enrollamiento. En la figura
(a) se muestra un esquema del sistema de control de consistencia por disolución en agua. En
la figura (b) se muestra el diagrama de bloques. Sean H(s)=1 y
K 1
Gc s G s
10s 1 2s 1
(b)
Solución.
C s Gc s G s K
R s 1 Gc s G s H s 20s 12s 1 K
2
1 1
ess
1 lim Gc s G s H s 1 K
s 0
si se quiere que este error sea inferior al 1%, la ganancia K debe ser mayor de 99. Para que esto
sea cierto nos queda por comprobar que el sistema es estable para K>99. Como se trata de un
sistema de segundo orden, con ecuación característica 20s2+ 12s + 1 + K = 0, la condición de
estabilidad es K+1>0 que, efectivamente, se cumple para K>99 y será posible conseguir un
error inferior al 1%.
Ejemplo 3.11
El sistema de control de la figura se ha diseñado para que cumpla las siguientes
especificaciones: (1) error estacionario menor del 10% para entrada en rampa, (2) máxima
sobreelongación menor del 5%, (3) tiempo de asentamiento (para un 2%) menor de 3 seg.
+ C(s)
R(s) K
- s s 0.4
1 K2s
Solución.
K
C s s s 2 K1
R s 1 K1
1 K 2 s s 2 2 1 K1K 2 s K1
s s 2 2
1 2
ess
K1 K1
lim s 1 K 2 s
s 0 s s 2
n2 K1 n K1
K1 K 2 2 K1 K 2
n 1
2 2 K1
Cada especificación impone alguna condición sobre los polos del sistema, concretamente:
2
(1) ess 0.1 n 20
K1
1 2
(2) M p 100e 5 0.69
4 4
(3) ts 2% 3 n
n 3
pero como 0.69 20 4 / 3 , cumpliendo las dos primeras especificaciones siempre se cumple
la tercera. La conjunción de ambas condiciones delimita la región del plano-s donde pueden
estar los polos del sistema para que se cumplan las especificaciones, concretamente la zona
rallada de la figura:
2
ess 0.1 K1 20
K1
2 K1 K 2
0.69 K 2
2 0.69 K1 1
2 K1 K1
El hecho de que todos los polos del sistema queden en el semiplano izquierdo del plano
complejo, no garantiza características satisfactorias de la respuesta transitoria. Por ejemplo si hay
polos dominantes complejos conjugados cercanos al eje imaginario, la respuesta transitoria puede
representar excesivas oscilaciones o puede ser muy lenta. Por tanto, para garantizar características
de respuesta transitorias rápidas y, sin embargo, bien amortiguadas, es necesario que los polos del
sistema queden fuera de una zona determinada del plano s, como puede ser la región rayada de la
figura 4.1.
R(s) + C(s)
G(s)
-
H(s)
1. La primera y segunda fila se forman con los coeficientes ai con i par del polinomio y
con los coeficientes ai con i impar, respectivamente.
2. Las siguientes filas se van formando con las dos filas inmediatamente superiores como
un producto cruzado de términos (véase la tabla 4.1 para un polinomio de sexto
orden).
3. El número total de filas es igual al número de coeficientes, y es útil indicar las
potencias de s asociadas (en orden decreciente) a cada fila de la tabla, como se ve en
la tabla 4.1.
s6 a0 a2 a4 a6
5
s a1 a3 a5 0
a1a2 a0 a3 a1a4 a0 a5 a1a6 a0 0
s4 A B a6 0
a1 a1 a1
Aa3 a1 B Aa5 a1a6 A 0 a1 0
s3 C D 0 0
A A A
BC AD Ca6 A 0 C 0 A 0
s2 E a6 0 0
C C C
ED Ca6
s1 F 0 0 0
E
Fa6 E 0
s0 a6 0 0 0
F
Una vez obtenida la tabla, el criterio de estabilidad de Routh se completa con el examen de
los signos de los coeficientes en la primera columna de la tabla. El criterio establece que:
Las raíces del polinomio están todas en el semiplano izquierdo si todos los elementos
de la primera columna de la tabla de Routh tienen el mismo signo.
El número de cambios de signos en esta columna determina el número de raíces del
polinomio con parte real positiva, es decir, el número de raíces en el semiplano
derecho.
Ejemplo 4.1
Sea el polinomio genérico de tercer grado, con todos los coeficientes positivos. Determinar la
condición de estabilidad en función de los coeficientes.
Solución.
a0 s 3 a1s 2 a2 s a3 0
s3 a0 a2
s2 a1 a3
a1a2 a0 a3
s1 0
a1
s0 a3 0
Como a0, a1, a2 y a3 son positivos, la condición de estabilidad es que a1a2>a0a3, ya que de esta
forma todos los elementos de la primera columna son positivos.
Ejemplo 4.2
Dada la siguiente ecuación característica, determinar si el sistema es estable o inestable:
s 4 2s 3 3s 2 4s 5 0
Solución.
s4 1 3 5 1 3 5
s3 2 4 0 1 2 0
2
s 1 5 1 5
s1 6 1
0
s 5 5
Queda en este ejemplo claro que el resultado de aplicar el criterio de Routh es invariable si se
multiplican o dividen los coeficientes de una misma fila por un número positivo, para
simplificar los cálculos. Esto es lo que se ha realizado en algunas de las filas de la tabla de
Routh anterior.
Ejemplo 4.3
Dada la ecuación: s 3 s 2 0 determinar cuantas raíces están localizadas en el semiplano
derecho del plano s.
Solución.
s3 1 1
s 2
0 2
2
s1 1
0
s 2
2
En la tabla de Routh que se ha formado, como es muy pequeño y positivo, ( 1 ) es
negativo. Hay por lo tanto dos cambios de signo, lo que indica la existencia de dos raíces con
parte real positiva. Las raíces de s 3 s 2 0 son 0.5000 1.3229i y -1.0000 de las que,
efectivamente, dos tienen parte real positiva.
Segundo: Todos los elementos obtenidos en una fila de la tabla de Routh son nulos.
Esto indica que existen raíces radicalmente opuestas respecto al origen en el plano s, es
decir:
Ejemplo 4.4
Dada la siguiente ecuación característica, determinar la localización de las raíces:
s 6 s 5 2 s 4 3s 3 7 s 2 4s 4 0
Solución.
Como los coeficientes son de distintos signos, evidentemente el sistema va a ser inestable.
s6 1 2 7 4
s 5
1 3 4
s4 1 3 4 P( s ) s 4 3s 2 4
dP( s )
s3 0 4 0 6 0 4s3 6s
ds
s2 1.5 4
s1 16.7
s0 4
1 3 1 3
Las raíces de s 6 s 5 2 s 4 3s 3 7 s 2 4s 4 0 son: 2 ; -2 ; j ; j ; j ; -j
2 2 2 2
En la fila s3 se ha presentado el 2º caso especial (todos los elementos de la fila iguales a cero),
entonces se ha formado el polinomio auxiliar con los elementos de la fila s4 y se ha
reemplazado la fila de ceros con los coeficientes del polinomio derivado de P(s). En la primera
columna de la tabla obtenida, sólo existe un cambio de signo, por lo tanto, sólo una raíz de la
ecuación principal está situada en el semiplano derecho (parte real positiva).
Ejemplo 4.5
Dado el sistema de control en lazo cerrado de la siguiente figura, determinar el rango de
valores de la ganancia K del controlador para los que el sistema es estable.
R(s) + C(s)
1
s s 2 s 2 s 1
K
-
Solución.
s4 1 3 K
s3 3 2
7
s2 K
3
9
s1 2 K
7
s0 K
Para que el sistema sea estable deben ser positivos todos los elementos de la primera columna.
Por tanto, el rango de estabilidad para este sistema en función de la ganancia del controlador,
viene dado por:
14 14
0K . Cuando K , el sistema se hace oscilatorio con amplitud de oscilación
9 9
constante, para una entrada escalón. Para K>14/9 , el sistema es inestable.
Ejemplo 4.6
Últimamente han adquirido gran importancia las grandes antenas para microondas en
radioastronomía y en el rastreo de satélites. Este tipo de antenas está expuesto a momentos de
torsión muy grandes debidos a las ráfagas de viento. Concretamente para una antena de 20 m
de diámetro, los experimentos muestran que un viento de 56 km/h ejerce una perturbación
máxima de 2 V en la entrada, Td(s) de la amplidina. Otro de los problemas del manejo de
antenas grandes es su resonancia estructural, que se pone de manifiesto en el siguiente modelo
para el conjunto (antena, motor de mando y amplidina):
100
G s
s 8s 100
2
ka
G1 s
0.2s 1
b) Determínese el error del sistema en estado estacionario, operando en lazo abierto (ks=0)
con R(s)=0, para una ráfaga de viento de 56 km/h.
Td(s)
Sensor
ks=1
Solución.
100ka
1 0.2 s 3 2.6s 2 28s 100 ka 1 0
0.2s 1 s 8s 100
2
se le pueden aplicar las condiciones de estabilidad determinadas en el ejemplo 4.1. Ello nos
lleva a que:
0.2ꞏ100 ka 1
100 ka 1 0 ka 1 28 ka 2.64
2.6
Por tanto, el rango de estabilidad del sistema en función de la ganancia del amplificador es:
1 ka 2.64
b) Con el sistema en lazo abierto, ks=0 y R(s)=0, la función de transferencia que relaciona a la
posición con la perturbación es la función del conjunto (antena, motor de mando y
amplidina):
s 100
Td s s 8s 100
2
100
ess 2 lim 2 rad 115º
s 0 s 8s 100
2
c) Con el sistema en lazo cerrado, ks=1 y R(s)=0, la función de transferencia que relaciona a la
posición con la perturbación es:
el error que producirá la ráfaga de viento en la salida del sistema es ahora función de ka:
100 0.2s 1 2
ess 2 lim
s 0
0.2s 1 s 2
8s 100 100ka 1 ka
si se quiere que este error sea inferior a 2o, por tanto inferior a 0.0349 radianes, ka debe ser
mayor que 56.31. Para este valor de ka, según el apartado (a), el sistema no es estable por tanto
nunca se podrá alcanzar el estacionario. El mínimo error posible es de:
2 2
0.55 rad 31.5º
1 ka maxima 1 2.64
Los polos del sistema en lazo cerrado son las raíces de la ecuación característica del
sistema. El procedimiento de descomposición en factores para hallar estas raíces, se hace bastante
laborioso si el grado del polinomio es tres o mayor, y no es conveniente porque es necesario
repetir los cálculos cada vez que uno de los parámetros del sistema (ya sea del controlador, del
proceso o del bloque convertidor) se modifica. Una solución para esto es la utilización de un
ordenador. Por ejemplo, MATLAB aporta una solución simple a este problema.
Lugar de las raíces (LR): Es el lugar geométrico del plano complejo s formado por las
raíces que se obtienen cuando el parámetro toma valores positivos, es decir, LR
cuando 0≤ K<∞.
Lugar de las raíces inverso (LRI): Es el lugar de las raíces cuando K toma valores
negativos.
Contorno de las raíces: Es el diagrama del lugar geométrico de las raíces, para un
sistema cuando tienes parámetros múltiples, es decir cuando varia más de un
parámetro.
R(s) + C(s)
G(s)
-
H(s)
La ecuación característica para este sistema en lazo cerrado se obtiene haciendo que el
denominador del segundo miembro sea cero. Es decir,
1+G(s)H(s)=0 o bien G(s)H(s)= -1
En muchos casos, G(s)H(s) contiene un parámetro de ganancia K. La ecuación
característica se podría expresar ahora:
1+KG(s)H(s)=0 KG(s)H(s)= -1
Como G(s)H(s) es una magnitud compleja y K es una variable real positiva, la ecuación
característica se podría descomponer en dos expresiones que determinan las condiciones de
magnitud y ángulo del LR.
Condición de magnitud o módulo:
1
G ( s) H ( s)
K
Condiciones de ángulo o argumento:
arg G ( s ) H ( s ) 2n 1 (n = 0,1,2,3,...)
De acuerdo con todo lo anterior, se pueden dar dos definiciones alternativas del lugar de
las raíces, desde el punto de vista de sistemas de control:
Ubicando los puntos y las asíntotas específicas y calculando los ángulos de salida de los
polos complejos y los ángulos de llegada a los ceros complejos, podemos construir la forma
general de los lugares geométricos de las raíces sin dificultad.
Para localizar los polos en el plano complejo se suele utilizar el símbolo ‘x’ y para los ceros,
la ‘o’. Como los polos complejos sólo ocurren en pares conjugados, los lugares geométricos de las
raíces son simétricos con respecto al eje real del plano s, lo que facilitará su trazado.
Los puntos del lugar geométrico que corresponden a K=0 son los polos en lazo abierto:
1 ( s z1 )( s z2 )...( s zm )
lim
K 0 K ( s p1 )( s p2 )...( s pn ) s p
i
Conforme K disminuye, el valor de s debe tender a uno de los polos en lazo abierto. Cada
rama del lugar geométrico de las raíces se origina en un polo de la función de transferencia en lazo
abierto G(s)H(s). Conforme K tiende a cada lugar geométrico tiende al cero de la función de
transferencia en lazo abierto o al infinito del plano complejo:
1 ( s z1 )( s z2 )...( s zm )
lim 0
K K ( s p1 )( s p2 )...( s pn ) s z
i
El valor de s debe aproximarse a uno de los ceros finitos en lazo abierto o a un cero en
lazo abierto en infinito. Así, los lugares geométricos de las raíces empiezan en los polos de
G(s)H(s) y terminan en los ceros de G(s)H(s) conforme K aumenta de cero a infinito (los polos y
los ceros incluyen tanto aquellos finitos como infinitos en el plano complejo).
Otro punto importante a tener en cuenta en el trazado consiste en que la gráfica del lugar
geométrico de las raíces tendrá tantas ramificaciones como raíces tenga la ecuación característica,
ya que estas raíces se moverán a lo largo de su tramo respectivo para los distintos valores de K. El
número de raíces o el de tramos equivale al número de polos del sistema en lazo abierto (o al de
ceros, el que sea mayor de los dos números). Así:
Regla 2: Determinar los lugares geométricos de las raíces sobre el eje real.
Se determinan a partir de los polos y los ceros en lazo abierto que se encuentran sobre él.
Los polos y los ceros complejos conjugados de la función de transferencia en lazo abierto no
afectan a la ubicación de los lugares geométricos de las raíces sobre el eje real porque la
contribución del ángulo de un par de polos o ceros complejos conjugados es 360º.
Para determinar si un punto si cualquiera sobre eje real pertenece al LR, la cantidad total
de polos y ceros reales a la derecha de este punto debe ser impar, para que dicho punto verifique
la condición del argumento.
arg G ( si ) H ( si ) 2n 1 (n = 0,1,2,3,...)
s Imag
s3 s4
Real
s2
s3 s4 s3
r3 r4 s4
r4 r3
Real Real
r2 r2
s2
s2
Aquí h=0 corresponde a las asíntotas con el ángulo más pequeño con respecto al eje real.
Aunque h puede tomar infinitos valores, el ángulo se repite, siendo el número final de asíntotas
diferentes, n - m. Todas las asíntotas cortan al eje real en un punto (centroide) que se obtiene del
modo siguiente:
n m
p j zi
j 1 i 1
nm
siendo pj y zi , los polos y los ceros finitos de la función G(s)H(s).
Debido a la propiedad de simetría del LR, las asíntotas también son simétricas respecto al
eje real:
Si n - m=1 la asíntota forma un ángulo de 180º
Si n - m=2 las asíntotas forman ángulos de 90º y –90º
Si n - m=3 las asíntotas forman ángulos de 60º,180º y –60º
Si n - m=4 las asíntotas forman ángulos de 45º, 135º, -135º y –45º.
Además las asíntotas forman ángulos entre sí de 360º/(n-m).
En caso de que haya polos y ceros complejos, éstos aparecen como pares complejos
conjugados. Al aplicar la expresión de , las partes imaginarias de cada par de polos o de ceros
complejos conjugados se cancelan entre sí, quedando la ecuación:
n m
Re( p j ) Re( zi )
j 1 i 1
nm
valores de K en los que se encuentran polos múltiples en lazo cerrado, tal y como se muestra en la
figura 4.4. Al ir aumentando el valor de K sobre el LR se obtienen los polos en lazo cerrado
marcados con en número 1. Al seguir aumentado los polos en lazo cerrado serían los marcados
con el 2. Si K sigue aumentado llega un momento en que los dos puntos se superponen (en el
punto p), correspondiendo a un valor máximo (o mínimo, si el tramo está entre dos ceros en vez
de entre dos polos) por lo que estos puntos se obtienen escribiendo la ecuación característica
como A(s) +K B(s)=0, despejando la K y buscando las soluciones reales de dK/ds. De las posibles
soluciones de esta ecuación, sólo serán puntos de dispersión aquellos que pertenezcan al LR
(puede dar más soluciones).
Esta derivada puede obtenerse de forma analítica o numérica por tanteo, como se verá en
los ejemplos de páginas posteriores, ya que en determinados casos la derivada no es inmediata y
es preferible obtener un valor de K aproximado.
Regla 5: Determinar el ángulo de partida (ángulo de llegada) del lugar geométrico de las
raíces, desde los polos complejos (ceros complejos).
Debemos encontrar las direcciones de los lugares geométricos de las raíces cercanas a los
polos y ceros complejos. Si se selecciona un punto de prueba y se mueve en la cercanía del polo
complejo (o del cero complejo), se considera que no cambia la suma de las contribuciones
angulares de los otros polos y ceros.
s1 Imag
1
r1
r4
s3 r3 4
3 2
s4
Real
r2
s2
En la figura, suponiendo que en los alrededores de s1 existe un tramo del LR, se quiere
obtener el ángulo de salida. Como un punto muy próximo al s1 debe cumplir la condición del
argumento se tiene que:
arg(G(si)H(si)) = arg(r4)- arg(r1)- arg(r2)- arg(r3) = 4 - 1 - 2- 3 = 180º
Como suponemos que el punto está muy próximo, prácticamente si s1, los valores de
los ángulos se toman como los de los radiovectores ri formados como los vectores que van desde
los polos y ceros hasta el polo s1.
donde rceros es el radiovector desde un cero hasta el polo sp y rpolos desde un polo hasta sp.
En posteriores ejemplos se aclarará este punto.
Regla 6: Encontrar los puntos en los que los lugares geométricos de las raíces cruzan el
eje imaginario.
Estos puntos se pueden encontrar con facilidad mediante:
El criterio de estabilidad de Routh, tratando de hacer cero filas completas de la tabla.
Suponiendo que s=j en la ecuación característica, igualando con cero la parte real y
la parte imaginaria y despejando y K. Los valores encontrados de representan las
frecuencias en las cuales los lugares geométricos de las raíces cruzan el eje
imaginario. El valor de K corresponde a la ganancia para cada frecuencia de cruce.
Regla 7: Tomando una serie de puntos de prueba en la vecindad amplia del origen del
plano s, trazar los lugares geométricos.
La parte más importante de los lugares geométricos de las raíces no está sobre el eje real ni en las
asíntotas, sino en la vecindad amplia del eje j y el origen, que puede hacer falta obtener con más
precisión. Este apartado puede no ser necesario para obtener una representación clara del LR.
Este valor puede calcularse en forma gráfica o analítica, como se verá en los ejemplos del
siguiente apartado.
La condición del módulo se mantiene pero la del argumento se modifica, de forma que
un punto pertenece al LRI si cumple que:
arg(G(s)H(s)) = 2 n, con n =0,1,2,...
Las ramas comienzan en los ceros para K=- y terminan en los polos para K=0.
Los puntos del eje real que pertenecen al LRI deben dejar a su derecha un número par
de polos y ceros. De esta manera, un tramo cualquiera del eje real pertenece al LR o
al LRI.
Los ángulos de las asíntotas se obtienen en el LRI mediante la ecuación:
360º h
Angulos de las asíntotas = (h=0,1,2,...)
nm
Ejemplo 4.7
El diagrama de bloques de un sistema de control es el mostrado en la figura:
R(s) + C(s)
K s 1
s s 2 s 4
-
Solución.
K ( s 1)
G ( s) H ( s)
s ( s 2)( s 4)
1 G s H s 0
En nuestro caso:
s s 2 s 4 K s 1 0 s 3 6s 2 8 K s K 0
3.- Puntos de partida y de llegada: el lugar arranca en los polos (K=0) y termina en -1 (una de
las ramas) y en el infinito (las otras dos ramas).
4.- Asíntotas: n – m = 3 – 1 = 2.
180º
0 90º
(2i 1) 2
i
nm 3 180º
1 270º
2
n m
Re( p ) Re( z )
j 1
j
i 1
i
0 2 4 1
2.5
nm 2
7.- Puntos de dispersión y de confluencia en el eje real: para determinar el punto en el que el
lugar abandona el eje real (valores de K máximo), se obtiene una expresión de K y se deriva:
( s 2 8)( s 4) dK 2s 3 9s 2 12s 8
K 0
s 1 s 1
2
ds
cuyo resultado es s=-2.91 y que está asociado a un valor de K de 1.5, obtenido en la expresión
anterior al sustituir el valor de s dado.
8.- Puntos de corte con el eje imaginario: para determinarlos aplicamos el Criterio de Routh a
la ecuación característica: (Regla 6).
s 3 6s 2 8 K s K
s3 1 8 K
s2 6 K
s1 b1 0
s0 K
1 8 K
6 K
b1
6
Ejemplo 4.8
Dibujar el LR para el siguiente sistema:
K ( s 3)
G(s) H (s)
( s 2) 2 4
Solución.
3.- Puntos de partida y de llegada: el lugar arranca en los polos -2-2j y -2+2j (K=0) y termina en
–3 y en el infinito (K= )
4.- Asíntotas: n – m = 2 – 1 = 1.
6.- Puntos de dispersión y de confluencia en el eje real: para determinar el punto en el que el
lugar abandona el eje real (valores de K máximo), se obtiene una expresión de K y se deriva:
s 2 4s 8 dK s 2 6s 4
K 0
s3 s 3
2
ds
cuyo resultado es s=-5.23, que pertenece al lugar de las raíces y s=-0.76, que no pertenece al
LR. El primer valor está asociado a un valor de K de 6.47.
1 2 3 2i 1ꞏ180º
2
atan 2 90º 180º 2 153.4º
3 2
8.- Puntos de corte con el eje imaginario: aplicamos el Criterio de Routh a la ecuación
característica:
s 2 4 K s 8 3K 0
s2 1 8 3K
s1 4 K
s0 b1
1 8 3K
4 K 0
b1
4 K
Con los puntos, ángulos y líneas determinados se obtiene un trazado del LR, tal como
representa la figura:
Ejemplo 4.9
Dada la siguiente función de transferencia:
K ( s 3)( s 2)
G(s) H (s)
( s 1)( s 2 2s 2)
Solución.
K ( s 3)( s 2)
G(s)
( s 1 j )( s 1 j )( s 1)
Asíntotas: n - m = 4 – 3 = 1.
Direcciones asintóticas:
(1 2 ) 180º 1 180º
0 180º
nm 1
i n j m
p z
i 1
i
j 1
j
a 2
nm
1 2 3 4 5 180º
3 180º 90º 1 2 90º 71.57º
Siendo 1 y 2 los ángulos correspondientes a los ceros en lazo abierto y 3, 4 y 5 los
correspondientes a los polos en lazo abierto.
( s 1)( s 2 2 s 2)
Puntos de dispersión en el eje real: Sabiendo que K y
( s 3)( s 2)
dK
haciendo 0 encontramos para s -5.75 una K 10.85
ds
s 1 s 2 2s 2 K s 2 s 3 s 3 3 K s 2 4 5K s 2 6 K 0
s3 1 4 5K
s2 3 K 2 6K
5K 2 13K 12
s1 0
K 3
s0 2 6K 0
La tercera fila no puede ser nula porque la ecuación de segundo grado no tiene solución real y
es positiva para K>-3. La cuarta nos indica un corte con los ejes para K=-1/3 que pertenecería
al LRI.
Ejemplo 4.10
El diagrama de bloques de la figura muestra un sistema realimentado de control:
R(s) + K C(s)
s 2
s
s 1
- 2600 26
1
1 0.04s
Determinar:
a) El LR del sistema.
Solución.
K
G ( s) H (s)
s 2
s
s 1 (1 0.04 s )
2600 26
65000 K
G ( s) H (s)
s ( s 50 10 j )( s 50 10 j )( s 25)
Asíntotas: n - m = 4 – 0 = 4.
1 180º
0 45º
4
3 180º
1 135º
(1 2 ) 180º 4
nm 5 180º
2 225º
4
7 180º
3 315º
4
in j m
pi z j
i 1 j 1 50 50 25
a 31.25
nm 4
1
K s ( s 50) 2 100 ( s 25)
65000
Derivando con respecto a s e igualando a cero se obtiene s=-45.96 que no pertenece al LR,
s=-38.64, que tampoco pertenece y s=-9.15 para el que se obtiene K=3.946
s4 1 5100 65000 K
3
s 125 65000
s 2 b1 b2
1
s c1 0
s 0 d1
1 51000
125 65000
b1 4580 ; b2 65000 K
125
125 65000
4580 65000 K
c1 65000 1774 K
125
d1 65000 K
Para que:
Los puntos de intersección del LR con el eje imaginario se obtienen resolviendo la ecuación
auxiliar, una vez sustituido el parámetro K por su valor:
%--- Gráfica del lugar geométrico de las raíces. Ejemplo 4.10 ---
sys= zpk([],[0 –25 –50+10i –50-10i],65000);
rlocus (sys);
grid;
b) Para hallar los polos dominantes para un factor de amortiguamiento = 0.5, obtenemos el
punto de corte del lugar de la raíces con la recta que partiendo del origen forma un ángulo
con el eje real negativo, siendo:
Ejemplo 4.11
Dado un sistema con realimentación unitaria cuya función de transferencia en bucle abierto
viene dada por:
K ( s 2 2 s 5)( s 5)
G ( s)
s ( s 2)( s 2)( s 4)
determinar:
b) Valor de K para que los polos dominantes del sistema tengan un factor de amortiguamiento
= 0.1 y una frecuencia natural n = 2.4 rad/s.
Solución.
K ( s 1 2 j )( s 1 2 j )( s 5)
G(s)
s( s 2)( s 2)( s 4)
Asíntotas: n - m = 4 – 3 = 1.
Direcciones asintóticas:
(1 2 ) 180º 1180º
0 180º
nm 1
in j m
pi z j
i 1 j 1
a 3
nm
s 4 4 K s 3 7 K 4 s 2 15 K 16 s 25 K
Obteniéndose que, para K = 3.938, el sistema en bucle cerrado presenta dos polos en:
s = 2.33j
b) Para que los polos dominantes presenten un factor de amortiguamiento = 0.1 y una
frecuencia natural n = 2.4 rad/s estarán situados en:
Comprobamos si los polos dominantes pertenecen al LR para algún valor de K. Para ello,
aplicamos el Criterio del Argumento a los Pd:
j m i n
err = 1.46º
p i Pd
2.39 x3.27 x 2.96 x 4.45
Kd K i 1
j m
5.32 x0.849 x 4.44
z
j 1
j Pd
resultando: K = 5.13
Ejemplo 4.12
Para el servosistema de gran antena, analizado en el ejemplo 4.6, represéntese el LR en
función de la ganancia del amplificador.
Solución.
100ka
G s H s
0.2s 1 s 2 8s 100
que no tiene ceros y tiene los polos en -4±9.1652j y -5. Por tanto sólo habrá LR en el tramo
del eje real (-∞ -5]. El número de asíntotas es nº de polos - nº de ceros = 3, que formarán
ángulos de 60o, 180o y 300o respecto al eje real, y tendrán su punto de corte en:
4 4 5
4.33
3
Habrá corte con el eje imaginario, puesto que en el ejemplo 4.6 se determinó que para ka>2.64
el sistema era inestable. Para determinar las raíces en el eje imaginario se puede resolver la
ecuación característica del sistema 0.2 s3 + 2.6 s2 + 28 s + 100 (ka + 1) = 0 sustituyendo
ka=2.64, o bien rescatar la fila s2 de la tabla de Routh del sistema de tercer orden analizado en
el problema 4.11.
a1s 2 a3 0
son los puntos de corte con el eje imaginario. Sin embargo, no habrá puntos de ruptura. Ya se
tiene información suficiente para trazar el LR representado en la figura.
Ejemplo 4.13
El siguiente diagrama de bloques corresponde a un servomecanismo de posición con
realimentación tacométrica (debido al bloque KT) y control proporcional (debido al bloque K).
b) Indicar sobre el lugar (a) las raíces correspondientes a K=4, y determinar para ellas las
características asociadas de respuesta transitoria: tr (tiempo de subida), Mp (máxima
sobreelongación) y ts (tiempo de asentamiento).
d) Sobre el lugar (c) determinar KT de forma que Mp=0.05 (que equivale a =0.707), y estimar
el tr y ts asociados a estas raíces.
e) Calcular el error estacionario en velocidad para el sistema con los parámetros del apartado
(d).
R(s) + + 1 1 C(s)
K
s1 s
- -
KT
Solución.
K
1 0
s s 1
Respecto al LR, la función de transferencia en lazo abierto a considerar tiene solamente dos
polos ( 0 y -1), y la representación es sobradamente conocida, véase figura:
b) Las raíces para K=4 se pueden obtener resolviendo la ecuación característica para este valor
de K o aplicando la condición de magnitud del LR. Si suponemos que para ese valor las raíces
son complejas conjugadas, según el LR deben ser del tipo -0.5±jx. Cualquiera de estas raíces
equidista del polo en el origen y del polo en -1, una distancia:
d 0.25 x 2
K 4
1 1 x 3.75 1.94
d1d 2 0.25 x 2
En este apartado se pide también determinar las características tr, Mp y ts asociadas a dichas
raíces. Para ello es necesario calcular en primer lugar el coeficiente de amortiguamiento y la
frecuencia natural de las raíces:
al ser pequeño se pueden aplicar las expresiones exactas o las aproximadas para el cálculo de
tr, Mp y ts:
1 2
atan
0.8 2.5
tr 0.94 tr 0.7125
n 1 2 n
1 2
M p 100e 44.43% M p 100 1 58.33%
0.6
ts 2%
n
1
ln 0.02 1 2 7.89 ts 2%
4
n
4
0.5
8
R(s) + C(s)
1
4
s s 1 KT
-
4
1 0
s s 1 KT
KT s
1 0
s s4
2
d 1
0 s s 4 2 s 1 0 s 4 0 s1,2 2
2 2
ds G s H s
como +2 no pertenece al LR, el único punto de ruptura es el -2. Con la información obtenida
se traza el LR de la figura:
d) Sobre el lugar anterior se pide determinar el valor de KT para que las raíces tengan
3º) Si el dibujo del LR es bastante exacto, como es el representado en la figura, basta con
obtener gráficamente el punto de intersección.
1er procedimiento: Dada la raíz -x+jx, los ángulos y módulos de los radiovectores que la unen
con el cero y los dos polos son respectivamente, véase figura:
2
Imag
1.94
r2
x
r1 1
x
r3 Real
3
-1.94
1 135º r1 2 x 2 x 2
x 1.94
0.5 x x 1.94
2 2
2 atan r2
0.5 x
x 1.94
0.5 x x 1.94
2 2
3 atan r3
0.5 x
tan 2 tan 3
tan 2 3
1 tan 2tan 3
x 1.94 x 1.94
0.5 x 0.5 x 2 x 2 0.52 1.942 0 x 2
x 1.94 x 1.94
1 ꞏ
0.5 x 0.5 x
x2 x2 x 2 2 x 2
r2 r3 1.05ꞏ3.48
KT 1.83
r1 2
También se piden en este apartado los valores de tr y ts asociados a las raíces. Como el
e) Con K=4 y KT=1.83, el sistema es estable, la función de transferencia en lazo abierto del
sistema equivalente es:
4
G s H s
s s 1.83 1
2
1 2.83
ess 71%
4 4
lim s
s 0 s s 2.83
En muchos problemas de diseño es necesario estudiar los efectos que sobre el sistema
tienen los cambios en varios parámetros diferentes a la ganancia K sobre los polos en lazo cerrado.
Esto último es lo que permite el lugar geométrico de las raíces. Cuando varían dos o más
parámetros los lugares geométricos de las raíces se denominan contornos de las raíces.
Veamos con un ejemplo como se pueden aplicar las mismas reglas del LR para construir el
contorno de las raíces.
Ejemplo 4.14
Estudiar los efectos que sobre las raíces del sistema en lazo cerrado (con realimentación
unitaria y control proporcional), tienen la ganancia K y el parámetro a. La función de
transferencia en lazo cerrado es:
C ( s) K
2
R( s ) s as K
Solución.
as
Que se puede expresar como 1 0
s K
2
K
que se puede expresar como 1 0
s2
La gráfica del lugar geométrico de las raíces de la ecuación aparece en la figura 4.5a.
Para construir los contornos de las raíces, supongamos que K es una constante (K=4). Los
polos en lazo abierto están en s=j2. El cero finito en lazo abierto está en el origen. La gráfica
del lugar geométrico que corresponde a este valor de K aparece en la figura 4.5b. Para
valores deferentes de K se producen lugares geométricos similares.
Figura 4.5: (a) Gráfica del lugar geométrico de las raíces para la ecuación (2). (b) Gráfica del lugar
geométrico de las raíces (0 a , K = 4). (c) Gráfica de los contornos de las raíces.
Los contornos de las raíces, diagrama que muestra los lugares geométricos de las raíces que
corresponden a 0 K se dibujan como en la figura 4.5c. Los contornos de las raíces
empiezan en los polos de la función de transferencia y terminan en sus ceros. Las flechas
indican la dirección del incremento en el valor de a.
Los contornos de las raíces muestran los efectos de las variaciones de los parámetros de
sistema sobre los polos en lazo cerrado. A partir de la gráfica de los contornos de las raíces de
la figura 4.5c observamos que, para 0 K , 0 a , los polos en lazo cerrado se
encuentran en el semiplano izquierdo del plano s y el sistema es estable.
5.1 INTRODUCCIÓN
Si el sistema es estable, la salida obtenida para una entrada senoidal x(t)=X sent será
y(t)=Y sen(t + )
en donde:
Y X G j
y
Im G ( j )
arg G j atan
Re G ( j )
es decir, el sistema presentará una salida senoidal de la misma frecuencia que la entrada
pero de amplitud y fase de salida diferentes de las de la entrada. De hecho, la amplitud de la
salida se obtiene del producto X|G(j)| en tanto que el ángulo de fase difiere del de la entrada en
una cantidad =arg G(j).
Por ello, la característica de respuesta de un sistema para una entrada senoidal se obtiene
directamente de Y(j) / X(j)= G(j). La función de transferencia senoidal G(j) es, por tanto,
una cantidad compleja y como tal admite varios tipos de representaciones:
Las características que presenta el trazado del diagrama de Bode pueden resumirse en las
siguientes:
Al calcular en decibelios la amplitud 20 log |G(j)| del sistema, los productos y las
divisiones (contribuciones de los ceros y de los polos) aparecen como sumas y
diferencias.
La fase de cada una de las contribuciones también puede sumarse algebraicamente a
la expresión total del diagrama de argumentos.
La mayoría de los diagramas se puede aproximar mediante tramos lineales, en lo que
se denomina aproximaciones asintóticas, que posteriormente pueden ampliarse
mediante una representación más exacta, mediante algunas correcciones.
Permite una fácil ampliación del rango de frecuencia, aunque no permite la extensión
del diagrama hasta frecuencias nulas, por el propio escalado logarítmico.
5.2.2.Trazado asintótico
Como ya se ha dicho la ventaja principal de usar una traza logarítmica es la facilidad
relativa de obtener las curvas de la respuesta en frecuencia. Los factores básicos que suelen
aparecer en una función de transferencia arbitraria G(j)H(j) son los obtenidos al representarla
de forma factorizada:
(1 s z ) i
G ( s) H ( s) K i 1
n
(1 s
j 1
pj )
y el vector de argumentos:
La ganancia K.
j
1
Polos y ceros en el origen
Veamos estos términos simples uno a uno y las contribuciones de cada uno a los dos
diagramas:
La ganancia K.
Un número mayor que la unidad tiene un valor positivo en decibelios, en tanto que un
número menor que la unidad tiene un valor negativo. La curva de magnitud logarítmica para una
ganancia constante K es una recta horizontal cuya magnitud es de 20 log K decibelios.
Es fácil observar que las diferencias en las respuestas en frecuencia de los factores 1/j y
j estriban en los signos de las pendientes de las curvas de magnitud logarítmica y en los signos
de los ángulos de fase. Ambas magnitudes logarítmicas se vuelven iguales a 0 dB en =1.
1
20 log n 20 log j 20n log dB
j
n
o bien:
Por tanto, las pendientes de las curvas de magnitud logarítmica para los factores (1/j)n y
(j)n son –20n dB/década y 20n dB/década, respectivamente, así como un valor para los ángulos
de fase de –n·90º y n·90º en todo el rango de frecuencia.
1
20 log 20 log 1 2 / p 2 dB
1 j / p
Para frecuencias bajas, tales que <<p, la magnitud logarítmica se aproxima mediante:
De esta forma, para el rango de altas frecuencias, la curva de magnitud logarítmica es una
línea recta con una pendiente de –20 dB/década.
El ángulo de fase exacto del factor 1/(1+j/p) es =-arctg /p. En una frecuencia cero,
p
el ángulo de fase es 0º. En la frecuencia de corte este ángulo es arctg arctg1 45º
p
mientras que en el infinito, el ángulo de fase se convierte en –90º.
Se puede calcular el error cometido al aproximar la curva de magnitud real por las dos
asíntotas. El error máximo ocurre en la frecuencia de corte y es aproximadamente igual a –3 dB
dado que
Una ventaja de las trazas de Bode es que, para factores recíprocos, por ejemplo, el factor
1+j/z, las curvas de magnitud logarítmica y de ángulo de fase sólo necesitan cambiar de signo.
Dado que:
1
20 log 1 j / z 20 log
1 j / z
y
1
arg 1 j / z arctg / z arg
1 j / z
la frecuencia de corte es igual para ambos casos. La pendiente de la asíntota de alta
frecuencia de 1+j/z es 20 dB/década, y el ángulo de fase varía de 0º a 90º conforme la
frecuencia se incrementa de cero a infinito.
La curva de magnitud logarítmica, junto con las asíntotas y la curva del ángulo de fase
para el factor 1+jT se representa en la figura 5.4:
Para el caso en que una función de transferencia determinada contenga términos como
(1+j/z)n, se hace una construcción asintótica similar. La frecuencia de corte esta todavía en =z
y las asíntotas son rectas. La asíntota de frecuencia baja es una recta horizontal en 0 dB, en tanto
que la asíntota de frecuencia alta tiene la pendiente de –20n dB/década o 20n dB/década.
Asimismo el error implícito en las ecuaciones asintóticas es n veces el que existe para (1+j/z)1 y
el ángulo de fase es n veces el de (1+j/z)1 en cada punto de frecuencia.
n2
s 2 2 n s n2
que para s=j se normaliza como
1
1 2 ( j ) ( j )2
n n
Si >1, este factor cuadrático se expresa como un producto de dos factores de primer
orden con polos reales. Si 0<<1, este factor cuadrático es el producto de dos factores complejos
conjugados. Las aproximaciones asintóticas para las curvas de respuesta en frecuencia no son
precisas para un factor con valores bajos de , ya que la magnitud y la fase del factor cuadrático
dependen de la frecuencia de corte y del factor de amortiguamiento relativo .
La curva asintótica de respuesta en frecuencia se obtiene del modo siguiente. Dado que:
2 2
1 2
20 log 20 log 1 2 2
n n
2
1 2 j j
n n
Para frecuencias bajas tales que <<n, la magnitud logarítmica se convierte en -20 log 1
= 0 dB. Por tanto, la asíntota de frecuencia baja es una recta horizontal en 0 dB. Para frecuencias
altas tales que >>n, la magnitud logarítmica se vuelve:
2
20 log 40 log dB
n 2
n
es decir, la ecuación para la asíntota de alta frecuencia es una recta con pendiente de –
10
40dB por década, dado que 40 log 40 40 log
n n
La asíntota de alta frecuencia intersecta a la de baja frecuencia en =n, dado que en esta
n
frecuencia 40 log 40 log1 0 dB . Esta frecuencia, n, es la frecuencia de corte para el
n
factor cuadrático considerado.
Cerca de la frecuencia de corte =n, ocurre un pico de resonancia. El factor de
amortiguamiento relativo determina la magnitud de este pico de resonancia. Es obvio que la
aproximación mediante las asíntotas genera errores, dependiendo la magnitud del error del valor
de . Para valores pequeños de este, el error es grande.
La figura 5.5 muestra las curvas exactas de magnitud logarítmica y las curvas exactas del
ángulo de fase para el factor cuadrático visto con varios valores de . En este caso el valor de va
desde 0.1 hasta 1.0 correspondiéndose respectivamente con las curvas de mayor sobrepaso a las
de menor sobrepaso.
Figura 5.5. Diagrama de Bode de un sistema de segundo orden para =0.01, 0.1, 0.2, 0.4, 0.6 y 1
2
1 n
arg 2
arctg 2
1 2 j
j 1
n n n
es decir, es una función de y de . En =0, el ángulo de fase es igual a 0º. En la
frecuencia de corte =n, el ángulo de fase es –90º, y en =, el ángulo de fase se convierte en
–180º. La curva del ángulo de fase tiene una pendiente simétrica respecto al punto de inflexión,
punto en el que =-90º.
Las curvas de respuesta en frecuencia para el factor 1 2 j j pueden
n n
obtenerse si simplemente se invierte el signo de la magnitud logarítmica y el del ángulo de fase del
1
factor 2
1 2 j j
n n
La magnitud de:
1
G ( j ) 2
1 2 j j
n n
es:
1
G j
2 2
2
1 2 2
n n
2 2
2
g 1 2 2
n n
sea mínima. Dado que g() se puede escribir como:
2 n 2 1 2 2
2
g 4 2 1 2
n
2
el valor mínimo de g() ocurre en =n(1-22)1/2. Por tanto, la frecuencia de resonancia r
es:
r n 1 2 2
Para >0.707, Mr = 1. Conforme tiende a cero, Mr tiende a infinito. Esto significa que, si
el sistema no amortiguado se excita en su frecuencia natural, la magnitud de G(j) se vuelve
infinita.
1 2 2
arg( j r ) arctg 90º arcsen
1 2
Ejemplo 5.1
El uso de las trazas de Bode con aproximaciones asintóticas toma mucho menos tiempo que
otros métodos utilizados para calcular la respuesta en frecuencia de una función de
transferencia. La facilidad de dibujar las curvas de respuesta en frecuencia para una función de
transferencia determinada y la facilidad para modificar la curva de respuesta en frecuencia
conforme se agrega una compensación, son las principales razones por las cuales las trazas de
Bode se usan tanto en la práctica.
10 s 3
G j
s s 2 s 2 s 2
Para evitar posibles errores al trazar la curva de magnitud logarítmica, es conveniente escribir
G(j) en la siguiente forma normalizada, en la que las asíntotas de baja frecuencia para los
factores de primer orden y el factor de segundo orden son la línea 0 dB.
j
7.5 1
3
G j
j j 2
j
j 1 1
2 2 2
j
2
a) 7.5 b) 1 j c) 1/ j d) 1 1 j e) 1 1 j
3 2 2 2
Las frecuencias de corte del segundo, cuarto y quinto términos son =3, =2 y =(2)1/2,
respectivamente. Obsérvese que el último término tiene el factor de amortiguamiento relativo
de 0.3536.
Para dibujar las trazas de Bode, es necesario representar previamente las curvas asintóticas
para cada uno de los factores por separado. A continuación se obtiene la curva compuesta
Hay que tener en cuenta, que cuando se agregan las curvas asintóticas individuales a cada
frecuencia, la pendiente de la curva compuesta es acumulativa. Debajo de =21/2, la gráfica
tiene la pendiente de –20dB/década. En la primera frecuencia de corte =21/2, la pendiente
cambia a –60dB/década y continúa a la siguiente corte de frecuencia en =2, en donde la
pendiente se convierte en –80dB/década. En la última frecuencia de corte =3, la pendiente
cambia a –60dB/década.
Una vez dibujada una curva aproximada de magnitud logarítmica, la curva real se obtiene
agregando correcciones a todas las frecuencias de corte y a las frecuencias una década por
encima y por debajo de las frecuencias de corte.
Empezamos por dibujar la parte de la recta de frecuencia más baja (es decir, la recta con la
pendiente de –20dB/década para < 21/2). Conforme la frecuencia aumenta, obtenemos el
efecto de los polos complejos conjugados (el término cuadrático) en la frecuencia de corte
=21/2. Los polos complejos conjugados provocan que las pendientes de la curva de magnitud
cambien de –20 a –60dB/década. En la siguiente corte de frecuencia, =2, el efecto del polo
es cambiar la pendiente a –80dB/década. Por último, en la frecuencia de corte =3, el efecto
del cero es cambiar la pendiente de –80 a –60 dB/década.
Para dibujar la curva de ángulo de fase completa, deben trazarse las curvas de ángulo de fase
de todos los factores. La suma algebraica de todas las curvas de ángulo de fase proporciona la
curva completa de ángulo de fase.
Ejemplo 5.2
Dibujar las trazas de Bode para la siguiente función de transferencia:
9 s 2 0.2s 1
G ( s)
s s 2 2.2s 7
num=[0 9 1.8 9] ;
Las funciones de transferencia que no tienen polos ni ceros en el semiplano derecho del
plano s son funciones de transferencia de fase mínima y los sistemas que las contienen se
denominan sistemas de fase mínima. Por el contrario aquellas funciones de transferencia que
tienen polos y/o ceros en el semiplano derecho del plano s son funciones de transferencia de fase
no mínima y los sistemas que las contienen son sistemas de fase no mínima.
Ejemplo 5.3
Considere la siguiente función de transferencia:
10( s 5)
G ( s)
( s 1)( s 7)
Sustituimos s por j, y ponemos la función en forma que podamos distinguir cada término
independientemente, tenemos:
j
7.143( 1)
G ( j ) 5
j
( j 1)( 1)
7
a) K = 7.143 0º
G ( s ) H ( s) F ( s )e t0 s
donde t0 es el retardo en el tiempo. La respuesta en frecuencia del sistema se puede
obtener a partir de las expresiones del módulo y del argumento. Así:
G ( j ) H ( j ) F ( j ) e j t0 F ( j )
arg G ( j ) H ( j ) arg F ( j ) arg e j t0 arg F ( j ) t0
De esta forma se aprecia que el retardo en el tiempo no afecta a la amplitud sino que
modifica la fase introduciendo un desfase proporcional al tiempo de retardo y a la frecuencia de la
señal de entrada.
Si se desea trazar el diagrama de Bode de un sistema con retardo, habrá que trazar el de
F(s) y añadir al diagrama de argumentos un desfase de to en cada frecuencia.
El retardo de transporte tiene un comportamiento de fase no mínima y tiene un atraso de
fase excesivo sin atenuación en frecuencias altas. Estos retardos de transporte ocurren, por lo
general en los sistemas térmicos, hidráulicos y neumáticos.
Ejemplo 5.4
Considere la siguiente función de transferencia:
1
G(s) e 0.1s
( s 1)
Sustituimos s por j, y ponemos la función en forma que podamos distinguir cada término
independientemente, tenemos:
1
G ( j ) e 0.1 j
j 1
K=1 0º
Como puede apreciarse, el desfase aumenta con la frecuencia y de forma proporcional a ella,
no convergiendo a ningún argumento constante, como ocurre con los sistemas vistos
anteriormente.
lim G ( j ) K p
0
dB
-20dB/década
20 log Kp
-40 dB/década
0
en la escala
logarítmica
-20dB/década
dB
20 log Kv
2 1 en la escala
logarítmica
La intersección del segmento inicial de –20 dB/década con la línea =1 tiene la magnitud
de 20 log Kv. Esto se observa del modo siguiente: en un sistema de tipo 1, G(j)=Kv/j, para
Kv
<<1. Por tanto, 20 log 20 log K v
j 1
La intersección del segmento inicial –20dB/década con la línea 0 dB tiene una frecuencia
cuyo valor numérico es igual a Kv. Para entender esto, es necesario definir la frecuencia en esta
Kv
intersección como 1; así, 1 o bien Kv = 1.
j1
Como ejemplo, considérese el sistema de tipo 1 con realimentación unitaria cuya función
de transferencia en lazo abierto es:
K
G s
s Js F
Si definimos la frecuencia de corte como 2 y la frecuencia de la intersección del segmento
K 1 3
Dado que 1 K v se deduce que 1 2 3 ,
2
o bien
F 3 2
En el diagrama de Bode, log 1-log 3 = log 3-log 2 Por tanto, el punto 3 está justo en
la mitad entre los puntos 2 y 1. Entonces, el factor de amortiguamiento relativo del sistema es
F
2
2 KJ 23
-40 dB/década
dB
-60 dB/década
20 log Ka
-20 dB/década
a=(Ka)1/2 en la escala
logarítmica
=1
La intersección del segmento inicial –40dB/década (o su extensión) con la línea =1 tiene
Ka
una magnitud de 20 log Ka. Dado que, a bajas frecuencias, G j para <<1 se
j
2
Ka
deduce que 20 log 20 log K a
j
2
1
Ka
20 log 20 log1 0
j
2
lo cual produce a Ka
Una ventaja de usar una traza polar es que representa, en una sola gráfica, las
características de la respuesta en frecuencia de un sistema en el rango de frecuencia completo.
Una desventaja es que la traza no indica en forma clara la contribución de todos los factores
individuales de la función de transferencia en lazo abierto, como sí ocurre en el diagrama de Bode.
1 1
G j arg(45º )
G j 0 1arg(0º ) T 2
y
Si tiende a infinito, la magnitud de G(j) tiende a cero y el ángulo de fase tiende a –90º.
1
1 2T 2
Im
1
T 0.5 Re
1 2T 2
Im
=0
Re
2
El punto de corte con los ejes se obtiene cuando 1 0 , es decir cuando =n.
n2
= =0
=
= = =
=r
=
Para el caso sobreamortiguado, conforme aumenta mucho más allá de la unidad, el lugar
geométrico G(j) tiende a ser un círculo. Esto se observa pues, para un sistema muy amortiguado,
las raíces características son reales y una es mucho más pequeña que la otra. Dado que para un
suficientemente grande el efecto de la raíz mayor sobre la respuesta se vuelve muy pequeño, el
sistema se comporta como uno de primer orden.
2 2
G j 1 2 j j 1 2 j
n n n n
Dado que la parte imaginaria de G(j) es positiva para >0 y aumenta en forma
monótona, además de que la parte real de G(j) se decrementa en forma monótona a partir de la
unidad, la forma general de la traza polar de G(j) es la que aparece en la figura 5.17. El ángulo
de fase está entre 0º y 180º.
=0.6 =0.8
=1
=0.4
=0.2
=0
Ejemplo 5.5
Consideremos la siguiente función de transferencia de segundo orden:
1
G s
s Ts 1
Vamos a obtener la traza polar de esta función de transferencia. Dado que la función de
transferencia senoidal se puede escribir como
1 T 1
G j j
j 1 jT 1 T
2 2
1 2T 2
lim G j 0 j 0 0 arg(180º )
Im
-T
Re
La traza de G(j) es asintótica para la línea vertical que pasa por el punto (-T,0). Dado que
esta función de transferencia contiene un integrador (1/s), la forma general de la traza polar
difiere sustancialmente de las funciones de transferencia se segundo orden que no poseen un
integrador.
en donde n>m, o el grado del polinomio de denominador es mayor que el del numerador,
tendrá las formas generales siguientes:
Para =0 o sistemas de tipo 0: el punto inicial de la traza polar (que corresponde a =0) es
finito y está sobre el eje real positivo. La tangente para la traza polar en =0 es perpendicular al
eje real. El punto terminal, que corresponde a =, está en el origen y la curva es tangente a uno
de los ejes.
Para =1 sistemas de tipo 1: El término j del denominador contribuye –90º al ángulo de fase
total de G(j) para 0 . En =0, la magnitud de G(j) es infinita y el ángulo de fase se
convierte en –90º. En frecuencias bajas, la traza polar es asintótica para una línea paralela al eje
imaginario negativo. En =, la magnitud se vuelve cero y la curva converge hacia el origen y es
tangente a uno de los ejes. Esto es lo que sucede con sistemas de tipo 0 compensados mediante
un PI, es decir, al introducirle un integrador puro en el lazo directo.
Para =2 o sistemas de tipo 2: el término (j)2 del denominador contribuye –180º al ángulo de
fase total de G(j) para 0 . En =0, la magnitud de G(j) es infinita y el ángulo de fase es
igual a –180º. En frecuencias bajas, la traza polar es asintótica para una línea paralela al eje real
negativo. En =, la magnitud se vuelve cero y la curva es tangente a uno de los ejes.
Las formas generales de las partes de frecuencia baja de las trazas polares de los sistemas
de tipo 0, tipo 1 y tipo2 aparecen en la siguiente figura.
Im
Sistema de tipo 1
Re
Sistema de tipo 2
Sistema de tipo 0
Se observa que, si el grado del polinomio del denominador de G(j) es mayor que el del
denominador, entonces los lugares geométricos G(j) convergen al origen en el sentido de las
manecillas del reloj. En =, los lugares geométricos son tangentes a uno u otro de los ejes.
Las forma de las curvas de las trazas polares serán más o menos complicadas en función
de la dinámica del numerador, es decir, de las constantes de tiempo del numerador de la función
de transferencia. La figura 5.20 muestra dos ejemplos de trazas polares de funciones de
transferencia con una dinámica distinta del numerador. Ante todo, en el análisis de los sistemas de
control, debe determinarse con precisión la traza polar de G(j) en el rango de frecuencias que
interese, dado que es en esta zona donde se centrará el estudio para, por ejemplo, un diseño
posterior.
Im = Im =
Re Re
1. La técnica del lugar de las raíces, que permite estudiar no sólo la estabilidad absoluta sino
también la relativa, ya que permite ver la situación exacta de los polos en lazo cerrado.
2. El criterio de Routh, que permitía determinar la estabilidad absoluta pues indicaba cuántos
polos en lazo cerrado estaban situados en el semiplano derecho
En este apartado se presenta uno de los métodos más utilizados para analizar la
estabilidad de sistemas que presenta las siguientes características:
Es un método gráfico.
Da información sobre la estabilidad absoluta y relativa, cuantificando ésta mediante
algunos parámetros.
Indica cómo mejorar esos índices de estabilidad, es decir orienta en la síntesis del
esquema de control.
Permite incorporar en la gráfica los retardos de los sistemas.
Al ser un método en el dominio de la frecuencia, da una idea de la respuesta en
frecuencia de los sistemas.
Teorema de mapeo: Dada una función F(s) racional y univaluada, que es analítica en todo s
excepto en un número finito de puntos, si se elige un camino cerrado arbitrario en el plano s tal
que la función F(s) sea analítica en todos los puntos de este camino , entonces el correspondiente
camino imagen del anterior F() es también cerrado y rodea al origen tantas veces como
diferencia entre el número de ceros y polos finitos de F(s) que son rodeados por el camino en el
plano s.
Si nz > np entonces N > 0 y por tanto el camino F() rodea N veces al origen en el mismo
sentido que . En caso de que N < 0 entonces el camino F() rodea N veces al origen en sentido
contrario que .
Imag
Tramo III
Tramo I
radio
0
Tramo VIII
Tramo VII
Tramo VI
-j1
Tramo V
El trazado del diagrama de Nyquist se obtiene calculando la imagen de cada uno de los
tramos de la figura 5.21. Los pequeños círculos del camino se han hecho rodeando a posibles
polos (en el origen o complejos puros) que tuviese la función F(s) , ya que el camino no debe
pasar por ninguna de las singularidades de F(s). En un caso normal los tramos II, VI y VIII no
existirán.
Una vez especificado el camino de Nyquist hay que obtener la función F(s). La estabilidad
del sistema de control en lazo cerrado viene dada por la ecuación característica, y por ello se elige
F(s)=1+G(s)H(s). A esta función se le puede aplicar el teorema de mapeo anterior. Sin embargo,
como lo que normalmente se conoce del sistema es la función de transferencia en lazo abierto
G(s)H(s), el criterio de Nyquist se puede establecer en términos de esta función respecto al punto
–1 + j 0 (simplemente desplazando los ejes una unidad). Así el teorema de Nyquist se puede
establecer como:
Teorema de Nyquist: Si la función de transferencia en lazo abierto G(s)H(s) tiene np polos finitos
en el semiplano derecho del plano s, el sistema de control en lazo cerrado es estable si y sólo si los
valores de G(s)H(s), obtenidos al recorrer el camino de Nyquist (es decir, la imagen de este camino
con F(s)=G(s)H(s) ) rodean np veces al punto crítico –1 + j0 del plano imagen y en sentido
contrario al elegido para el camino de Nyquist (en la figura 5.21 este camino es en el sentido de
las agujas del reloj).
Justificación: Si se aplica el teorema de mapeo a F(s), se tiene que N=nz – np, siendo en
este caso np y nz el número de polos y ceros de F(s)=1+G(s)H(s) en el semiplano derecho del
plano s (interior del camino de Nyquist ) y N el número de veces que la imagen de F() rodea al
origen del plano s.
Como los polos de F(s) son los mismos que los de G(s)H(s) y, como el origen del plano
1+G(s) H(s) corresponde al punto –1+j0, si se conoce np (los polos de G(s)H(s) en el semiplano
derecho de s que son los mismo que los polos de 1+G(s)H(s) ) y si se determina gráficamente el
valor de N (rodeos en torno al punto crítico –1+j0) nz será cero si y sólo si N=-np, es decir np
rodeos en sentido contrario al del camino de Nyquist , en el caso actual, sentido antihorario.
Ejemplo 5.6
Trazar el diagrama de Nyquist de:
1
G s
( s 1) 2
Se elige el camino de Nyquist de la figura 5.21. En este caso existen sólo tres tramos: el
primero el eje imaginario positivo, el segundo el semicírculo de radio infinito y el tercero el
semieje imaginario negativo, que es simétrico del primer tramo.
1 1 1 2 2 j
G j
( j 1) 2 (1 2 2 j ) (1 2 ) 2 4 2
lim G ( j ) 1 j 0 lim G ( j ) 0 j 0
0
Cálculo de cortes con los ejes: Parte real nula = 1. Tomamos solo = +1 ya que el
trazado es semieje positivo
G ( j ) = -0.5j
1
1
G j j
(Re 1) 2
1 2 j
lim G ( j ) lim e Si va desde 90º a -90º con el sentido 90º, 45º, 0º, -45º, -90º
R R R
entonces el límite anterior es un arco de circunferencia de radio nulo que va desde 180º hasta
180º y que recorre el sentido 180º, 90º y 180º
Tercer tramo: Sustituimos la s por -j, obteniéndose el simétrico del primer camino.
1
Figura 5.22. Diagrama de Nyquist de
( s 1) 2
Se ha añadido un círculo de radio unidad para mejor visualización de las características del
diagrama. Como se observa, el número de rodeos al punto crítico –1+j0 es cero. Por tanto el
número de ceros de la ecuación característica en el semiplano derecho, obtenido al despejar nz
de la ecuación N = nz - np con N=0 y np=0 es cero, y por tanto el número de polos del sistema
en lazo cerrado dentro del semiplano derecho también es nulo y el sistema es estable.
Ejemplo 5.7
Diagrama de Nyquist del sistema
1
G s
s ( s 1) 2
Se procede como en el problema anterior, teniendo en cuenta que el camino de Nyquist tiene
que rodear al polo en el origen, es decir, tendrá un nuevo tramo.
1 j (1 j ) 2 2 1 2
G j j
j ( j 1) 2 2 (1 2 ) 2 (1 2 ) 2 (1 2 ) 2
lim G ( j ) 2 j
0
lim G ( j ) 0 j 0
Cálculo de cortes con los ejes: Haciendo la parte imaginaria nula = 1. Tomamos sólo
=+1 ya que el trazado es semieje positivo. Se obtiene:
G ( j ) 0.5
1
1
G j
Re (Re j 1) 2
j
1 3 j
lim G ( j ) lim e , que representa a un arco de circunferencia con radio nulo y ángulo
R R R3
3 desde -270º a 270º siguiendo la trayectoria -270º, -135, 0º, 135º, 270º.
Tercer tramo: Este tramo es el simétrico del primer tramo respecto del eje real.
Cuarto tramo: Sustituimos la s por e j , un arco de radio 0 y desde –90º hasta 90º.
1 1
lim lim e j , arco de radio infinito y ángulo con recorrido 90º, 45º, 0º, -45º, –90º
0 e j 0
1
Figura 5.23. Diagrama de Nyquist de G s
s ( s 1) 2
El diagrama de Nyquist permite analizar la estabilidad del sistema en lazo cerrado para
cualquier valor de la ganancia K, tal y como posibilitan el diagrama de las raíces o el criterio de
Routh, métodos con los que obviamente se llega a las mismas conclusiones que con el presente.
Para realizar este análisis es necesario sustituir el punto crítico. Si para la ecuación
característica 1+G(s)H(s)=0 el punto crítico es el –1+j0, para la que incluye a la ganancia variable
K , es decir para 1+KG(s)H(s) =0 este punto crítico sería –1/K+j0. Y es con este punto con el que
hay que interpretar el diagrama de Nyquist.
El punto –1/K puede estar a la derecha del punto –0.5+j0 o a su izquierda. Veamos estos dos
casos:
a) Si –1/K está a la izquierda de –0.5+j0, es decir, si –1/K < -0.5, o puesto de otra forma, si
K<2, entonces N=0 y np=0, por lo que nz =0 y el sistema es estable.
b) Si –1/K está a la derecha de –0.5+j0, es decir si –1/K > -0.5, o K>2, entonces N=2 y np=0,
por lo que nz = 2 y el sistema es inestable.
Por tanto el rango de estabilidad es K<2, como muestran otros métodos de análisis de la
estabilidad.
Ejemplo 5.8
s 1
Trazar el diagrama de Nyquist de G s K
s
s ( 1)
10
El camino de Nyquist será el mismo que el del caso anterior (tendrá el tramo adicional que
rodea el origen, debido a la presencia del polo).
110
lim G ( j ) j y lim G ( j ) 0 j 0
0 100
10(Re j 1)
G j
Re j (Re j 10)
lim G ( j ) 0e j , que representa un arco de radio cero y ángulo de -90º a 90º pasando por
R
e j 1 1
lim lim e j ( 180)
0 j 0
e ( j
e 1)
10
que es una circunferencia de radio y arco desde -90º a –270º con trayectoria -135º, -180º, -
225º.
s 1
Figura 5.24. Diagrama de Nyquist de G s
s
s ( 1)
10
Se ha calculado que el cruce con el eje real se produce en la frecuencia 10 y en el punto
G ( j ) 1 . Con esta información y con la suministrada por la gráfica del diagrama de
10
Nyquist es posible analizar su estabilidad para cualquier valor de K. Existen dos posibilidades:
a) -1/K < -1. En este caso, es decir, si K < 1, el punto crítico, -1/K queda dentro de la
semicircunferencia de radio infinito de la figura 5.24, estando rodeado una vez en sentido
horario por la imagen del camino de Nyquist. Como np =1, ya que existe un polo en lazo
abierto dentro del semiplano derecho, el sistema es inestable (nz=2).
b) –1/K > -1. Aquí, con K > 1, el punto crítico –1/K está rodeado por una vez por la imagen
del camino de Nyquist, pero en sentido antihorario, por lo que N=-1. Como np=1, el sistema
es estable y nz=0 (no tendrá ningún polo en lazo cerrado dentro del semiplano derecho).
Ejemplo 5.9
Trazar el diagrama de Nyquist del sistema
s 1
G s
( s 2 4)
En este caso el camino de Nyquist debe evitar los dos puntos singulares que el sistema tiene
en el eje imaginario, en concreto j2. Se muestran entonces los diferentes tramos:
j 1
G j
2 4
1
=0 lim G ( j )
0 4
=2 lim G ( j ) j
2
G ( j ) arg G j atan 63º
1
= -2 lim G ( j ) j
2
G ( j ) arg G j atan 63º
1
2 j e j 1
G j
(2 j e j ) 2 4
lim G ( j ) e j atan 2 90º , que corresponde a un arco de circunferencia de radio 0 y ángulo
0
= lim G ( j ) 0 j 0 G ( j )
Re j 1 1
lim lim e j
R R 2 e 2 j 1 R R
s 1
Figura 5.25. Diagrama de Nyquist de G s
( s 2 4)
Se observa que el sistema es estable ya que el número de vueltas es cero, por lo que el número
de polos en lazo cerrado ubicados en el semiplano derecho del plano complejo es nulo.
Ejemplo 5.10
Considérese un sistema en lazo cerrado, cuya función de transferencia en lazo abierto, se
obtiene mediante:
2
G s
(3s 1)(2 s 1)
2
Figura 5.26. Diagrama de Nyquist de G s
(3s 1)(2s 1)
Vemos que G(s)H(s) no tiene polos en el semiplano derecho del plano s y el punto -1 + j0 no
está encerrado por el lugar geométrico G(s)H(s). Por tanto, este sistema es estable.
Ejemplo 5.11
Sea el sistema cuya función de transferencia en lazo abierto es:
3
G s
s (2 s 1)(3s 1)
num = 3
den = [ 6 5 1 0 ]
sys=tf(num,den)
nyquist (sys)
3
Figura 5.27. Diagrama de Nyquist de G s
s (3s 1)(2s 1)
Vemos en la gráfica que el punto -1 + j0 queda encerrado dos veces en sentido de las agujas
del reloj. El primer encierro se ve claramente en el gráfico y el segundo viene dado por el
contorno que une rama de arriba con la de abajo mediante un arco de radio infinito que
engloba todo el semiplano derecho, característica de la gráfica de Nyquist.
Por tanto como el contorno encierra dos veces el punto - 1 + j0, el sistema tiene dos polos en
el semiplano derecho del plano s, y el sistema es inestable.
Ejemplo 5.12
Veamos un sistema similar, pero en lugar de darle una ganancia grande, numerador de G(s)=3,
le daremos una mas pequeña, por ejemplo 0.5. En este caso:
0.5
G s
s (2s 1)(3s 1)
num = 0'5;
den = conv([2 1 0],[3 1 ]);
sys=tf(num,den);
Nyquist (sys)
axis([-4 4 -4 4]);
0.5
Figura 5.28. Diagrama de Nyquist de G s
s (3s 1)(2s 1)
En la que vemos que el contorno no encierra el punto - 1 + j0, y puesto que el número de
polos de G(s)H(s) en el semiplano derecho es 0; el sistema es estable, como se vio en el
problema anterior. Vemos pues que, para lograr una buena precisión, el valor de la ganancia
debe de ser grande. Sin embargo, una ganancia grande provoca estabilidad deficiente o incluso
inestabilidad; por tanto, es necesario incluir una red de compensación en el sistema (esto se
verá más adelante).
Ejemplo 5.13
Determine la estabilidad del sistema cuya función de transferencia en lazo abierto es:
1
G s H s
s (3 / 4 s 1)
1
Figura 5.29. Diagrama de Nyquist de G s
3
s s 1
4
Con MATLAB:
num = 1;
den = [3/4 -1 0];
sys=tf(num,den);
nyquist (sys)
Aunque en la gráfica no se aprecia, el contorno queda cerrado por un arco de radio infinito
que une la rama de abajo con la de arriba por la izquierda en sentido horario. Por tanto, el
contorno encierra una vez el punto - 1 + j0 en el sentido de las agujas del reloj. Por tanto
N=1. Dado que Z = N + P y P = 1 (s = 4/3), entonces Z = 2; el sistema tiene dos polos en
lazo cerrado en el semiplano derecho del plano s y es inestable.
Este es uno de los inconvenientes de utilizar la resolución numérica (como hace MATLAB)
para dibujar el diagrama de Nyquist, ya que no permite la visualización de los círculos de radio
infinito, y por tanto, analizar si el diagrama rodea o no al punto crítico.
Ejemplo 5.14
2( s 3)
Vea si el sistema cuya función de transferencia en lazo abierto G s H s es
s ( s 1)
estable o no.
num = [ 2 6 ]
den = [1 -1 0]
sys=tf(num,den);
nyquist (sys);
2 s 3
Figura 5.30. Diagrama de Nyquist de G s
s s 1
En la figura anterior, vemos que el diagrama de Nyquist encierra una vez el punto - 1 + j0 en
sentido inverso a las agujas del reloj; o sea N = -1, luego como Z = N + P, tenemos Z = 0 que
indica que no hay un cero de 1+G(s)H(s) en el semiplano derecho del plano s, luego el sistema
en lazo cerrado es estable. Este es un caso en el que un sistema inestable en lazo abierto se
vuelve estable cuando se cierra el lazo.
Ejemplo 5.15
Sistemas con retardo: Como ya se vio en apartados anteriores, si se desea trazar el diagrama
de Bode de un sistema con retardo, habrá que trazar el de F(s) y al diagrama de fases añadir un
desfase de to en cada frecuencia. Para analizar la estabilidad en lazo cerrado, el criterio de
Routh no permitía tratar este tipo de sistemas. Sin embargo, el criterio de Nyquist sí. Para ello,
a cada frecuencia hay que añadir el desfase comentado anteriormente, pero sin alterar el
módulo. Por tanto, un punto correspondiente a una frecuencia sufre un giro hacia la izquierda,
como se aprecia en la figura 5.31 donde los puntos de las curvas con las distintas
circunferencias corresponderían a puntos con la misma frecuencia. En esta figura se muestra la
representación de los sistemas dados por la función de transferencia siguiente:
1
G ( s) e st0
s ( s 1)
y con retardos t0 =0, 0.5 , 1 y 5. Se aprecia que el efecto de un retardo puede desestabilizar a
un sistema, como ocurre con el último de ellos. Lo dibujamos en MATLAB mediante las
instrucciones:
El ancho de banda (BW del inglés band width) se define como la frecuencia a la que el
módulo de la función de transferencia en lazo cerrado es el 70.7% del módulo a baja frecuencia.
En el diagrama de Bode, esto corresponde a una caída de 3 dB.
El ancho de banda es una medida del rechazo al ruido de los sistemas y a un ancho de
banda mayor, más sensible será el sistema ante estas señales con componentes de altas
frecuencias. Por la misma razón, existe una relación directa entre el ancho de banda y la velocidad
de respuesta de los sistemas y las medidas que caracterizan a ésta, como son el tiempo de pico y
el de subida. Por tanto, a la hora de diseñar un sistema de control existirá un compromiso entre
conseguir una respuesta lo suficientemente alta y un rechazo de perturbaciones aceptable.
El margen de ganancia (MG o GM del inglés Gain Margin) es uno de los criterios más
empleados para medir la estabilidad relativa de sistemas de control. En el dominio de la frecuencia,
el margen de ganancia se emplea para indicar la cercanía de la intersección del eje real negativo
con el diagrama de Nyquist de G(j) al punto crítico (-1,j0). Se define como
1
Margen de ganancia = MG = 20 log 20 log G ( j cmg ) en dB
G ( j cmg )
donde cmg es la frecuencia de cruce de fase que es donde arg G ( j cmg ) 180º .
Esta frecuencia es aquélla en la que el diagrama de Nyquist corta el eje real negativo. Por tanto
podemos definir el margen de ganancia como: La cantidad de ganancia en (dB) que se
pueden añadir al lazo antes de que el sistema en lazo cerrado se vuelva inestable.
Este hecho puede apreciarse en la figura 5.33
G( j cmg )
cmg
Cuando el diagrama de Nyquist no corta al eje real negativo, lo que ocurre en sistemas de
primer o segundo orden puros, el margen de ganancia es infinito en dB. Esto indica, que en estos
sistemas se puede aumentar sin límite su ganancia sin que se haga inestable.
Cuando el diagrama de Nyquist pasa por el punto –1+j0, el margen de ganancia es 0 dB,
lo que significa que la ganancia de lazo no puede crecer más, ya que el sistema está en el margen
de inestabilidad. El sistema es críticamente estable.
Cuando el cruce de fase está a la izquierda del punto -1+j0, el margen de ganancia es
negativo en dB, y la ganancia de lazo se debe reducir en el margen de ganancia para alcanzar la
estabilidad. El sistema es inestable.
5.5.4.Margen de fase
El margen de ganancia por sí solo puede resultar insuficiente para indicar la estabilidad
relativa de un sistema. El margen de fase (MF o PM del ingles Phase Margin) se define como el
ángulo en grados que la traza G(j) debe rotar alrededor del origen, para que el cruce de ganancia
pase por el punto –1+j0. El cruce de ganancia es un punto sobre el diagrama polar de G(j) en el
cual su magnitud es igual a 1. Esto puede apreciarse en la figura 5.34.
MF
cmf
Al igual que ocurre con el margen de ganancia, el margen de fase por sí solo no aporta
información suficiente sobre la estabilidad relativa, debiéndose complementar con el margen de
ganancia.
Para un sistema de fase mínima, los márgenes de ganancia y de fase deben ser positivos
para que sea el sistema estable. Los márgenes negativos indican inestabilidad. Para obtener un
sistema satisfactorio, el margen de fase debe estar entre 30 y 60, y el margen de ganancia debe
ser mayor de 6 dB. Con estos valores un sistema de fase mínima tiene una estabilidad garantizada.
Para sistemas de fase no mínima el criterio de estabilidad en relación a los márgenes de fase y de
ganancia es el opuesto, es decir, estos han de ser negativos para estabilidad.
Ejemplo 5.16
Veamos como ejemplo los márgenes de ganancia y de fase del sistema cuya función de
transferencia de lazo es:
1
G s
s ( s 1)( s 3)
MG: La frecuencia cmg a la que se produce el margen de ganancia será aquella que verifique:
90º atan atan 180º
3
1
El módulo para esta frecuencia es G ( j cmg ) 0.0833
cmg 1 2
cmg 32 cmg
2
MF: El margen de fase se mide a la frecuencia cmf a la que el módulo de G(j) vale la unidad.
Así:
1
G ( j cmf ) 1
cmf 1 2
cmf 32 cmf
2
Resolviendo esta ecuación resulta cmf =0.316. De esta forma, el margen de fase será:
cmf
MF 90 arctg cmf arctg 180º 66.4º
3
La figura 5.36 ilustra cómo obtener los márgenes de fase y ganancia a partir del diagrama de
Bode.
Las gráficas han sido obtenidas a partir de las siguientes funciones, donde además se incluye la
función margins, que devuelve el margen de ganancia (no en dB) el de fase, y las frecuencias a la
que se producen.
Ejemplo 5.17
Este ejemplo muestra la influencia del cambio de la ganancia en un sistema de tercer orden en
su estabilidad relativa, analizando ésta mediante el margen de fase, obtenido a partir del
diagrama de Nyquist.
Kp
GH
s s s 1
2
GH K
3 2 p
1 GH s s s K p
La traza de Nyquist corta al eje real negativo en un punto que está muy lejano a la derecha del
punto (-1,j0). En concreto el margen de ganancia es de 3.3 (unos 10 dB) y el de fase de 70.7º.
La respuesta al escalón está bien amortiguada. Ambas gráficas se muestran a continuación:
La intersección se acerca al punto (-1+j0); el sistema aún es estable, ya que el punto (-1+j0) no
está encerrado, pero la respuesta al escalón tiene un sobrepaso máximo grande.
Los márgenes de ganancia y fase son ahora de 2 (6 dB) y 50.3º respectivamente. Se muestran
las mismas gráficas del anterior para compararlas.
En este caso el diagrama de Nyquist pasa a través del punto (-1+j0) y el sistema es
críticamente estable, estando ahora en el límite de estabilidad. La respuesta al escalón se
vuelve oscilatoria con amplitud constante, como se observa en las siguientes gráficas.
El margen de fase es ahora de 0º, típico de un sistema con estas características y el margen de
ganancia de 1 (0 dB).
A continuación se muestra el código MATLAB utilizado para generar las ocho gráficas
anteriores:
close all
clear all
for kp=[0.3 0.5 1 1.5];
a=tf(kp,[1 1 1 0]);
% Muestra los márgenes de estabilidad
[Gm,Pm,wcg,wcp] = margin(a)
w=logspace(-2,2,1000);
% Dibuja el diagrama de Nyquist
[re,im]=nyquist(a,w);
figure
plot(re(1,:),im(1,:),'r');
axis([-2 2 -2 2])
circle(1,0,0,200,'g')
xlabel('Real')
ylabel('Imag')
title('Diagrama de Nyquist')
figure
% Dibuja la respuesta temporal en lazo cerrado
b=feedback(a,1);
[y,t]=step(b,25);
plot(t,y,'r');
ylabel('Amplitud')
xlabel('Tiempo')
title('Respuesta temporal ante salto escalón')
end
Ejemplo 5.18
Un aeroplano de despegue vertical (VTOL) es un vehículo inherentemente inestable y requiere
un sistema automático de estabilización. En la figura aparece el diseño, en forma de diagrama de
bloques, para un sistema de estabilización de un VTOL. A 40 nudos, la dinámica del vehículo se
representa aproximadamente por la función de transferencia
10
G(s)
s 0.25
2
K ( s 7)
Gc ( s ) H s s
s2
d) Determínese la amplitud máxima de la respuesta en frecuencia del sistema en lazo cerrado, así
como la frecuencia a la que se presenta.
+
r(s) + (s)
+
Gc G
-
s
840 1 s
60( s 7) s 7
Gc ( s )G ( s ) H ( s )
( s 2)( s 0.25) s s 2
2
1 1
2 0.5
expresados adecuadamente para el trazado del diagrama de Bode. En este diagrama hay que
considerar las siguientes contribuciones:
- Una constante 840 que contribuye en el diagrama de amplitud con 20 log 840 58.5 db y en
el diagrama de fase con 0o.
- Dos polos imaginarios puros a la frecuencia 0.5, que contribuye en el diagrama de amplitud
2
con 20 log 1 2 , contribución que se hace infinito a la frecuencia 0.5 y en el diagrama
0.5
de fase con -180o a partir de la frecuencia 0.5, por lo que éste presentará una discontinuidad
fuerte a esa frecuencia.
2
20 log 1 2
y en el diagrama de fase con arctg
2 2
2
20 log 1 2
y en el diagrama de fase con arctg
7 7
b) En el diagrama de fase se observa que nose produce el corte con la línea de -180o, luego el
margen de ganancia del sistema es MGdB = . En el diagrama de amplitud se observa que el
corte con la línea de 0 dB se produce a la frecuencia 60 rad/s , para la que el diagrama de
fase presenta un valor de aproximadamente -95o, luego el margen de fase es aproximadamente de
180o-95o=85o.
G s 10 s 2 10 s 2
3 2
1 G s H s Gc s s 2 s 0.25 10 Ks s 7 s s 2 10 K s 0.25 70 K 0.5
2
que para K=6 corresponde a un sistema estable, puesto que anteriormente se han determinado
margen de ganancia infinito y margen de fase positivo. La ganancia en estado estacionario de esta
función de transferencia es igual a 40, independiente del valor de K. Luego un escalón unitario en
Td(t) provocan 40 rad de desviación en la posición (t) respecto a la referencia r(t).
d) Respecto a la entrada R(s) y con K=6 la función de transferencia del sistema en lazo cerrado es:
Gc s G s 60 s 7
3
1 Gc s G s H s s 62s 2 420.5s 0.5
que tiene todos sus ceros (-7) y polos (-0.0012, -7.7446, -54.2542) situados en el eje real del
semiplano-s izquierdo, correspondiendo la menor frecuencia a uno de los polos. Sin representar
el diagrama de Bode, por la naturaleza de los ceros y polos y por la situación relativa de éstos se
sabe que la máxima ganancia del sistema se presenta a baja frecuencia, con valor
Ejemplo 5.19
La función de transferencia en lazo abierto de un sistema de control con realimentación unitaria
es:
K
G s
s 0.1s 1 s 1
a) Determinar el valor de K para que el pico de resonancia Mr del sistema sea igual a 1.4.
b) Determinar el valor de K para que el margen de ganancia del sistema sea 20 dB.
c) Determinar el valor de K para que el margen de fase del sistema sea 60o.
G s K
M s
1 G s s 1 0.1s s 1 K
1
Mr
2 1 2
para el cálculo de pico de resonancia Mr. El problema se debería resolver de forma gráfica
haciendo uso de la carta de Nichols. Esta carta nos permite obtener la respuesta en frecuencia del
sistema en lazo cerrado a partir de la representación magnitud-fase del sistema en lazo abierto.
No obstante, también se puede resolver de forma analítica si suponemos que la ecuación
característica del sistema tiene dos raíces complejas dominantes, y se exige que estas raíces tengan
un coeficiente de amortiguamiento =0.4 para que Mr=1.4. Concretamente, expresando la
ecuación característica de la forma:
1 0.1 n2 0.8 n p
K 0.1 n2 p K 1.24
la parte real de las raíces complejas, n=0.44, es mucho menor que la raíz real, p=10.1, luego
efectivamente con K=1.24 se puede conseguir el objetivo (Mr=1.4).
b) Este apartado se puede resolver gráficamente, representando el diagrama de Bode de G(s) con
K=1, obteniendo su margen de ganancia, y determinando cuánto hay que desplazar el diagrama
de amplitud para que el margen de ganancia sea 20 dB. Este último paso equivale a aplicar la
expresión
Hay que recordar que, normalmente, un aumento de K disminuye el margen de ganancia del
sistema. Utilizar la siguiente figura, que representa el diagrama de Bode exacto para G(s)/K, para
calcular el valor de K según el procedimiento descrito.
arg G j 180º 90º atan atan 180º a b 90º
10
2
tan a b 1 cg 10
10
G j 10 dB
20 log 10 20 log 1 10 20 log 1
10
100
20.82 MG 20.82 dB
Luego para disminuir el MG, la ganancia K aumentar tal que aporte 20logK=0.82dB => K=1.1.
Compárese con el resultado obtenido gráficamente.
c) De manera similar al apartado (b), éste se puede resolver gráficamente sobre el diagrama de
Bode de G(s) con K=1, determinando cuánto hay que desplazar el diagrama de amplitud para que
el margen de fase sea 60o. Este último paso equivale a aplicar la expresión
Utilizar el diagrama de Bode de G(s)/K para calcular el nuevo valor de K. Otra forma, analítica,
difiere de la anterior en que la amplitud a la que se presenta el margen de fase deseado se calcula
analíticamente. Así:
MF 60º arg G j 120º 90º atan atan 120º a b 30º
10
10 tan 30º 0.58 0.058 2 1.1 0.58 0 0.51
2 cf
1
10
0.512
G j 0.51 dB 20 log 0.51 20 log 1 0.512 20 log 1 20 log1.74
100
como la amplitud es mayor de 0 dB, la ganancia K debe disminuir tal que 1.74K=1 => K=0.57.
Compárese con el resultado obtenido gráficamente.
Elegir una frecuencia de diseño c para la red superior a la frecuencia cf (frecuencia a la
que el diagrama de amplitud toma el valor de 0 dB) del sistema sin compensar. De esta
forma la amplitud del sistema sin compensar a la frecuencia c es menor que la unidad y
Calcular la fase que debe aportar la red de adelanto para que a la frecuencia c se
Determinar la frecuencia del cero so y del polo sp de la red tal que a la frecuencia c la
ganancia de la red compense a la ganancia del proceso. Equivale a utilizar las expresiones
sen sen
so c sp c
1
cos cos G j c
G j c
intento. En el caso de que el margen de ganancia del sistema sin compensar sea finito, se recomienda
que la frecuencia de diseño c de la red se elija entre cf y cg.
Ejemplo 5.20
En la figura se muestra el sistema de control para el vehículo del módulo de excursión lunar del
Apolo 11. La amortiguación del vehículo es despreciable y la posición se controla por chorros de
Vehículo
Compensador Actuador
r(s) + (s)
V T 1
Gc(s) K2
-
Js 2
K1
Englobando la dinámica del actuador, del vehículo y del sensor, el proceso a compensar tiene la
siguiente función de transferencia:
K K1 K 2
G s siendo K
s2 J
En este caso concreto, la cf del sistema sin compensar se obtiene de la siguiente expresión:
K
20 log K 40 log 0 1 cf K
2
luego si se elige c=3 > cf , se está imponiendo la siguiente condición a la ganancia K<9.
K
arg G j 180º G j c 0.1111
9
50º 180º (180º ) 50º
0.766 0.766
so 3ꞏ 0.275 ; sp 3ꞏ 4.3224
9 0.6428 0.6428 0.1111
Los diagramas de amplitud y de fase de la red diseñada se determinan a partir de las expresiones
siguientes:
s
1 15.72 s 0.275
Gc s 0.275
s
1 s 4.3224
4.3224
2 2
Gc j 20 log 1 20 log 1
0.2752 4.32242
arg Gc j atan atan
0.275 4.3224
wsp
20 log 24dB
wso
=50º
c=3
La siguiente figura representa los diagramas de amplitud (parte superior) y de fase (parte inferior)
del sistema sin compensar (gráfica inferior) y del sistema compensado (gráfica superior). Los de
este último se obtienen sumando las contribuciones en amplitud y en fase del sistema sin
compensar y las de la red de adelanto. Los márgenes, MG=infinito y MF=50o, están dentro de
las especificaciones.
MF=50
c
Nos queda por comprobar si la otra especificación, tiempo de asentamiento menor de 2 seg, se
cumple. Para ello se calcula la función de transferencia del sistema en lazo cerrado:
Gc s G s 15.72s 4.32
3
1 Gc s G s s 4.32 s 2 15.72 s 4.32
cuyos polos son -0.3 y -2±3.24j, en el que domina la constante de tiempo =1/0.3=3.33 asociada
al polo real y al que corresponde un tiempo de asentamiento aproximado de 4=13.33 seg, muy
superior al especificado. Por tanto el diseño realizado no es totalmente adecuado. Habría que
probar otros valores de c y K (superiores, para que el ancho de banda lo sea también, y por tanto
disminuya el tiempo de asentamiento).
b) La dificultad de encontrar solución al problema planteado en (a) con una red de adelanto tiene
fácil explicación mediante el lugar de las raíces. Concretamente el lugar de las raíces de la planta
coincide con el eje imaginario, por lo que el sistema con control proporcional estaría siempre en
el límite de la estabilidad. La incorporación de una red de adelanto, introduce un cero y un polo,
estando el cero siempre a la derecha del polo, por tanto el lugar de las raíces tendría una forma
similar a la representada en la figura (que corresponde al sistema con la red compensadora del
apartado a).
que buscan las asíntotas y con parte real n > 2, necesaria para que el tiempo de asentamiento
sea inferior a 2 seg, pero siempre existirá una tercera raíz real que es dominante respecto a ellas
(es lo que ha ocurrido en el apartado a). Al probar con otros valores de c y K superiores se
consigue situar el cero de la red mucho más a la izquierda, de tal manera que la raíz real, aún
siendo dominante, tiene una constante de tiempo asociada inferior a 0.5 seg. Otra solución
posible es que el cero del compensador cancele uno de los polos de la planta, y el lugar de las
raíces tenga la forma representada en la siguiente figura, sobre el que sí es más fácil conseguir las
especificaciones.
Figura 5.50. Lugar de las raíces del sistema con cancelación de un polo y un cero
s K
1 0 s 2 ps K 0
s p s2
luego de las especificaciones, se tiene que el sistema se debe comportar como de segundo orden
con =0.5 y n>4, que se traducen a:
p 2 n 4 y K n2 16
Esta expresión equivale a suponer que la contribución de la red de retardo a esa frecuencia
va a ser = -5o. El margen de fase deseado debe ser mayor o igual que el especificado.
Situar la frecuencia so asociada al cero de la red de retardo una década antes de c.:
so 0.1 c
Determinar la frecuencia sp asociada al polo de la red tal que la ganancia de la red a la
so
sp
G j c
Observación: Este procedimiento garantiza que se consigue un margen de fase algo menor que el
deseado. Si además existe especificación de margen de ganancia u otra especificación, puede ocurrir
que con la primera elección de la frecuencia so no se consiga la especificación y sea necesario
Ejemplo 5.21
Un sistema de control con realimentación unitaria tiene una planta con la siguiente función de
transferencia:
K
G s
s s
s 10 1
2 6
Se desea tener un coeficiente de error de velocidad Kv=20, un margen de fase de 45o y un ancho
de banda igual o mayor que 2 rad/s. Diseñar una red de retardo de fase para cumplir dichas
especificaciones.
K
K v lim sG s
s 0 10
luego para que Kv sea igual a 20, la ganancia debe ser de K=200.
Para trazar el diagrama de Bode de la planta sin compensar pero con K=200 hay que tener en
cuenta las siguientes contribuciones: una constante de 20, un polo a la frecuencia 0, un polo real
negativo a la frecuencia 6 y un polo real negativo a la frecuencia 200.
En la figura, que se incluye más adelante, se puede ver dicho diagrama junto con el del sistema
compensado.
De las gráficas superiores, que corresponden al sistema sin compensar, se observa que los
márgenes de ganancia y de fase son muy pequeños y, por tanto, están muy por debajo de las
especificaciones.
wsp
20 log 14dB
wso
-5º
wc
Se pide diseñar una red de retardo. Como hay que conseguir Kv=20, se recomienda considerar la
ganancia obtenida anteriormente como parte del sistema a compensar y aplicar el procedimiento
de diseño anterior. En este caso concreto, si se elige MF=45o se obtienen los siguientes
resultados:
s
1
0.35 s 0.35
Gc s
1 5ꞏ s 0.07
s
0.07
La siguiente gráfica representa los diagramas de amplitud y de fase de los sistemas sin compensar
(gráficas superiores) y compensado (gráficas inferiores), el nuevo margen MF=45o está dentro de
las especificaciones.
Nos queda por comprobar si la otra especificación, el ancho de banda, es igual o mayor que
2rad/s. Para ello se calcula la función de transferencia del sistema en lazo cerrado:
Gc s G s 498.065s 175.5417
4
1 Gc s G s s 26.0731s 121.9017 s 2 506.842s 175.5417
3
que tiene un cero en -0.35, debido a la red, y polos en -0.3780, -2.1983±j4.1447, -21.4786. La
proximidad de uno de los polos al cero hace que el ancho de banda del sistema venga
prácticamente fijado por el polo complejo conjugado. Concretamente, se sabe que WB es
superior a n , siendo n=4.69 en este caso, por tanto, WB es mayor que 2 rad/s y se cumple la
especificación. En la siguiente figura se puede ver el diagrama de Bode del sistema en lazo
cerrado, donde se pone de manifiesto la afirmación anterior. Gráficamente se puede determinar
que el punto de corte con la recta de -3 dB y por tanto el ancho de banda se produce a la
frecuencia 6 rad/s.
condición particular:
De esta forma se garantiza que el controlador elegido va a poder compensar la fase del
sistema a la frecuencia c.
Calcular la fase que debe aportar el controlador para que a la frecuencia c se consiga
PD hacer KI=0.
KI sen
K D c
c G j c
Ejemplo 5.22
Una mesa estabilizada de precisión usa un tacómetro y un motor de corriente continua como se
muestra en la figura. Se desea mantener una alta precisión en estado estacionario para el control
de velocidad, para lo que se elige un controlador proporcional integral. Selecciónense los
parámetros del controlador tal que el sistema tenga una sobreelongación de aproximadamente el
10% (margen de fase aproximado de 60o) y un tiempo de asentamiento entre 0.5 y 1 s.
Motor
R(s) + C(s)
0.1
25 Gc(s)
s 0.1
-
25
0.1s 1
Tacómetro
La función de transferencia del sistema en lazo abierto sin compensar se obtiene del producto de
las funciones de transferencia del motor y carga y del tacómetro:
2.5 25
FTLA
s 0.1 0.1s 1 s 1 s 1
0.1 10
Para trazar el diagrama de Bode de la planta sistema sin compensar, que se incluye en una figura
más adelante junto con el diagrama del sistema compensado, hay que tener en cuenta las
siguientes contribuciones: una constante de 25, un polo real negativo a la frecuencia 0.1 y un polo
a a
MF 60º arg G j a 120º atan atan 120º a b 120º
0.1 10
a a
0.1 10 tan 120º 1.73 1.73 2 10.1 1.73 0 6
a a a
1 a a
0.1 10
arg G j c 115.4º
1 0.0802
Kp 2.22 K I 5ꞏ 0.89
0.45 0.45
0.89
Gc s 2.22
s
Se muestra a continuación la figura del compensador. La siguiente figura representa los diagramas
de amplitud del sistema sin compensar (gráfica inferior) y del sistema compensado (gráfica
superior), y los diagramas de fase sin compensar (gráfica superior) y compensado (gráfica
inferior). En los que se observa la reducción del margen de fase.
20logKp7dB
-4.6º
c
Nos queda por comprobar si la otra especificación, tiempo de asentamiento entre 0.5 y 1 s, se
cumple. Para ello se calcula la ecuación característica del sistema en lazo cerrado:
0.89 0.1 25
1 Gc s G s Gt s 1 2.22 0 s 3 10.1s 2 56.5s 22.25 0
s s 0.1 0.1s 1
cuyas raíces son -0.43 y -4.84±5.38j, en las que domina la raíz real con constante de tiempo
=1/0.43=2.33 asociada al polo real y a la que corresponde un tiempo de asentamiento
aproximado de 4=9.32 sg, muy superior al especificado. Por tanto el diseño realizado no es
totalmente adecuado. A continuación se hace una discusión sobre el lugar de las raíces que aclara
la dificultad de encontrar una solución al diseño PI.
El lugar de las raíces del sistema compensado con un bloque proporcional tiene la forma típica de
un sistema de segundo orden, en el que es posible conseguir raíces complejas con parte real igual
a (0.1+10)/2=5.05, y por tanto con tiempo de asentamiento ts=4/5.05=0.79 seg, que está
comprendido entre 0.5 y 1 como se especifica en el diseño. Bastaría imponer la especificación de
10% de sobreelongación (=0.6) para obtener el valor de la ganancia proporcional necesaria. El
único inconveniente de este diseño es que, debido al error en estado estacionario, no se alcanzaría
Ejemplo 5.23
Repítase el problema anterior usando como compensador una red de adelanto y compárense los
resultados.
En el enunciado de este problema se plantea el diseño de una red de adelanto. La solución más
inmediata que evite la aparición de un polo real, más próximo al eje imaginario que las raíces
complejas, es la compensación mediante cancelación del polo en -0.1. La ecuación característica
debe ser:
s 0.1 0.1 25
1 0 s 2 10 p s 10 p 25 0
s p s 0.1 0.1s 1
Por otra parte, según las especificaciones, el sistema se debe comportar como de segundo orden
con Mp=10% (=0.6) y:
4
0.5 ts 1 6.67 n 13.33
0.6 n
8 10 p 16 y 44.48 10 p 25 177.69
Así pues, p=4 puede ser una solución. La siguiente figura recoge el lugar de las raíces en este caso.
El diseño realizado sigue teniendo el inconveniente del diseño con controlador proporcional
respecto al diseño con controlador PI, ya que no elimina el error en estado estacionario.
Los distintos tipos de acciones básicas de control que podemos encontrarnos se clasifican en
las siguientes:
Control
Automático
+ C(s)
R(s)
Amplificador Actuador Planta
-
Sensor
En este tipo de controlador el actuador sólo presenta dos posiciones predeterminadas que, en
la mayor parte de los casos, son solamente de encendido y apagado. Su principal ventaja es su
simplicidad y bajo coste, que lo hace ideal para su empleo en muchos sistemas de control, tanto
industriales como domésticos.
u (t ) K p e(t )
U(s)
KP
E(s)
+
Kp
-
En el control proporcional la señal u(t) siempre sigue al error de forma que, a mayor error,
con más intensidad se produce la acción de control.
du(t)
Kie(t) u(t) Kie(t)dt
dt
U(s) K i
E(s) s
+ Ki
-
s
La primera consecuencia de la aplicación del control integral es que, para aquellas plantas que
no posean un integrador en el origen en la función de transferencia en lazo abierto (sistemas tipo 0),
el error pasa a ser cero ante la entrada escalón. Ya se demostró con anterioridad que en los sistemas
tipo 0 el error ante entrada escalón era constante. Por tanto, se consigue una importante mejora en la
respuesta estacionaria del sistema.
Otra característica del control integral es que la señal de control es, en todo momento, el área
encerrada bajo la curva temporal del error. Por ello, cuando el error sea cero la salida del controlador
integral será constante y, así mismo, cuando el error se mantenga constante, la señal de control será
una rampa, debido a la acumulación de área.
t
Kp
u(t) Kp e(t) e(t) dt
Ti 0
U (s) 1
kp 1
E(s) Ti s
+ K p (1 Ti s)
- Ti s
+
K p (1 Td s )
-
La acción de control derivativa presenta cierto carácter de anticipación puesto que responde a
la velocidad de cambio del error (su pendiente instantánea), con lo que se consigue que el valor que
adquiere dicho error no sea demasiado grande. Por lo tanto, se emprenden con antelación acciones
correctoras que, en definitiva, aumentan la estabilidad del sistema. Por otro lado, la acción derivativa
reduce la sobreoscilación de la salida añadiendo amortiguamiento, con lo que puede incrementarse el
valor de la ganancia del sistema para reducir el error en estado estacionario.
Por otra parte, el control derivativo tiene la desventaja de que aumenta el ancho de banda del
sistema global, es decir, amplifica las señales de ruido (según su diagrama de Bode, este controlador
amplifica las señales de alta frecuencia), lo que puede traer como consecuencia la saturación de los
actuadores.
Por último, comentar que la acción de control derivativa nunca se emplea sola, sino junto con
un control proporcional, de modo que en la señal de control no intervenga solo la pendiente del error
sino también el error mismo.
Todo sistema de control se diseña para conseguir unas ciertas especificaciones que,
generalmente, hacen referencia a la respuesta transitoria (velocidad de respuesta, sobreoscilaciones,
etc), estacionaria (precisión o error estacionario), estabilidad absoluta y estabilidad relativa.
En sistemas complejos o que requieren para el diseño del sistema de control de un análisis
más elaborado, las especificaciones anteriores se suelen indicar tanto de forma cualitativa como
cuantitativa. En problemas de diseño sencillos, por el contrario, se proporcionan de manera general
sólo cuantitativamente como valores numéricos precisos. Por supuesto, en los primeros sistemas
comentados puede ocurrir que no sea posible que se cumplan todos los requerimientos
simultáneamente, debiendo llegar en este caso a una solución de compromiso dando una mayor
prioridad a aquellas especificaciones de relevancia superior.
Para cumplir con estas especificaciones, en algunos casos tan sólo será necesario introducir
una ganancia en el sistema (control proporcional), pero puede ser que con esto se disminuya la
estabilidad relativa o que incluso se llegue a la inestabilidad, por no decir las modificaciones que
puede sufrir la respuesta transitoria (todo ello observable a través del diagrama lugar de las raíces).
En estos casos, se hace imprescindible añadir un controlador más complejo, denominado
compensador, de manera que se puedan conseguir las especificaciones requeridas.
Emplearemos aquí el enfoque del lugar geométrico de las raíces para, gráficamente, con la
introducción de los ceros y polos en lazo abierto del compensador, determinar sus ubicaciones y
establecer los polos en lazo cerrado ajustando también la ganancia para cumplir las especificaciones
de desempeño. Esto implica utilizar un método de ensayo y error variando la posición de los polos y
ceros del compensador hasta que consigamos nuestro objetivo.
Por lo tanto, es preciso construir de nuevo el diagrama lugar de las raíces para cada variación
del compensador, por lo que el empleo del ordenador para el diseño es primordial por el ahorro de
tiempo que supone. El objetivo final será que los polos dominantes en lazo cerrado cumplan
efectivamente las especificaciones de desempeño impuestas al principio.
El efecto de la adición de polos se observa en las dos figuras siguientes, en las que se tiene el
lugar de las raíces de un sistema (a la izquierda) y el correspondiente al del sistema total después de
añadir un polo (a la derecha).
1 2
0.8
1.5
0.6
1
0.4
0.5
0.2
Imag Axis
Imag Axis
0 0
-0.2
-0.5
-0.4
-1
-0.6
-1.5
-0.8
-1 -2
-2 -1.5 -1 -0.5 0 0.5 1 -2 -1.5 -1 -0.5 0 0.5 1
Real Axis Real Axis
Incluso si se añade otro polo más, el lugar de las raíces se transforma en esto:
1.5
0.5
Imag Axis
0
-0.5
-1
-1.5
-2 -1.5 -1 -0.5 0 0.5 1
Real Axis
Este efecto se aprecia partiendo del lugar de las raíces de un sistema, que es el que se
muestra a continuación:
1.5
0.5
Imag Axis
-0.5
-1
-1.5
-2 -1.5 -1 -0.5 0 0.5 1
Real Axis
2 2
1.5 1.5
1 1
0.5 0.5
Imag Axis
Imag Axis
0 0
-0.5 -0.5
-1 -1
-1.5 -1.5
-2 -2
-2 -1.5 -1 -0.5 0 0.5 1 -2 -1.5 -1 -0.5 0 0.5 1
Real Axis Real Axis
2.5
1.5
0.5
Imag Axis
-0.5
-1
-1.5
-2
-2.5
-2 -1.5 -1 -0.5 0 0.5 1
Real Axis
Control PID:
Con este controlador conseguimos combinar las tres acciones vistas anteriormente, de forma
que la señal de control viene dada por la expresión:
t
kp de(t )
u (t ) k p e(t )
Ti e(t ) dt
0
k p Td
dt
en forma de función de transferencia:
U ( s) 1
K p (1 Td s )
E ( s) Ti s
El diagrama de bloques representativo de este tipo de controladores es el siguiente:
+ K p (1 Ti s TiTd s 2 )
- Ti s
La acción proporcional varía cada vez que se modifica el error y alcanza un valor estacionario
cuando lo alcanza éste. La acción integral es acumulativa (es el área bajo la curva del error, como se
comentó con anterioridad) y, por tanto, tiene en cuenta la historia pasada del error. Solo alcanza un
estacionario cuando el error es cero. En cuanto a la acción derivativa, es anticipativa (proporcional a
las variaciones instantáneas del error) y se anula cuando el error se hace constante.
Proporcional (P): Td = 0, Ti
Proporcional-integral (PI): Td = 0
Proporcional-derivativo (PD): Ti
1
s
G ( s) aT
1
s
T
Se observa que está constituida por un cero y un polo. La ubicación relativa entre los dos
determina el tipo de red:
Para a > 1 tenemos una red de adelanto, en la que el cero se encuentra a la derecha del
polo en el plano complejo.
Para a < 1 tenemos una red de retraso y, en este caso, el cero se encuentra a la
izquierda del polo.
Si se analiza el diagrama de Bode para los dos tipos de redes encontramos que la red de
adelanto representa un filtro paso-alto y que, además, introduce un desfase positivo entre la entrada
y salida sinusoidales. Es por esta razón por la que a este tipo de red se le denomina así. En cuanto a
la red de retraso, representa un filtro paso-bajo y el desfase que introduce es negativo entre la
entrada y la salida.
A la vista de la función de transferencia se extrae que los parámetros de diseño para estos
controladores son:
T: nos da la posición del polo.
a: nos da la posición relativa entre el cero y el polo.
Indicar, por último, que también existen redes de adelanto-atraso, que combinan los dos tipos
de redes anteriores introduciendo una fase positiva en un rango determinado de frecuencias y una
fase negativa en otro rango distinto.
e1 -
e2 +
e0
e0 K e1 e2
Como primer circuito sencillo que se empleará como complemento, como se verá después, de
otros más complejos, se analizará el amplificador operacional en su configuración inversora.
R2
i2
R1
i1
-
e'
+
ei e0
ei e ' e ' e0
i1 i2
R1 R2
Suponiendo que sólo una pequeña parte de la corriente i1 fluye por los terminales de entrada
del amplificador, tal y como se supuso con anterioridad, las corrientes i1 e i2 deben ser
aproximadamente iguales. En consecuencia:
ei e ' e ' e0
R1 R2
ei e
0
R1 R2
O, lo que es lo mismo:
R2 E ( s) R
e0 ei 0 2
R1 Ei ( s) R1
Se observa de esta expresión que la salida se invierte en signo respecto a la entrada, además
de amplificarse o atenuarse en el factor R2 / R1. En el caso en que R1 = R2 se tendría un circuito
inversor de signo.
Por estas razones, y para una configuración de circuito específica con amplificadores
operacionales en la que se basan todos los controladores que hemos visto con anterioridad, se puede
plantear el tratar de llegar a un procedimiento sistemático que nos permita determinar la función de
transferencia de estos circuitos. Concretamente, la configuración que se va a emplear es la siguiente:
I(s)
Z2(s)
I(s)
Z1(s) -
+
Ei(s) E0(s)
Ei ( s) Z1 ( s ) I ( s ) E0 ( s ) Z 2 ( s ) I ( s )
E0 ( s ) Z ( s)
2
Ei ( s ) Z1 ( s )
Por lo tanto, la función de transferencia siempre puede obtenerse para este tipo de circuitos
como el cociente cambiado de signo de la transformada de Laplace de las impedancias Z2 y Z1
indicadas en la figura anterior. De esta manera, de forma sistemática se puede calcular la función de
transferencia sin necesidad de tener que trabajar con las tensiones, corrientes y mallas del circuito.
Como ejemplo se puede plantear calcular la función de transferencia del siguiente circuito:
C
i3
i2
R1 R2
i1
-
+
ei e0
Se puede llegar con facilidad a las siguientes expresiones para las transformadas de Laplace
de las impedancias Z1 y Z2 si comparamos el circuito propuesto con el de la figura anterior:
1 1 1 1 R2
Z1 ( s ) R1 Cs Z 2 (s)
Z 2 ( s) 1 R2 R2 1 R2Cs
Cs
E0 ( s ) R 1
2
Ei ( s ) R1 1 R2Cs
R2
R4
R1
R3
-
-
+
ei +
e0
E0 ( s ) R4 R2
Ei ( s ) R3 R1
Siendo Ei(s) la señal de error y E0(s) la señal de control, según todo lazo cerrado de control.
Lógicamente, si de la segunda etapa inversora lo único que se desea es el cambio de signo se
dimensionará la red de polarización del segundo amplificador de modo que R3 = R4. Se obtendría
entonces como ganancia del controlador P el cociente R2 / R1.
Como se vio con anterioridad, el controlador integral proporciona una señal de control
proporcional a la integral del error y, básicamente, tiene como función de transferencia un polo en el
origen. El circuito con amplificadores operacionales que permite realizar esto es el siguiente (con su
correspondiente etapa inversora adicional por las razones de cambio de signo expuestas en el
apartado anterior):
C2
R4
R1
R3
-
-
+
ei +
e0
1 E (s) 1
Z1 ( s ) R1 Z 2 ( s)
C2 s Ei ( s ) R1C2 s
E0 ( s ) R4 1 R4 1
Ei ( s ) R3 R1C2 s R1 R3C2 s
Se observa que el circuito introduce un polo en el origen, es decir, un integrador puro, por lo
que efectivamente se trata de un controlador integral o controlador I.
C1 R2
R4
R1
R3
-
-
+
ei +
e0
1 1 E ( s) R
C1s Z 2 ( s ) R2 2 1 R1C1s
Z1 ( s ) R1 Ei ( s ) R1
E0 ( s ) R4 R2 RRC 1
1 R1C1s 2 4 1 s
Ei ( s ) R3 R1 R3 R1C1
Se observa que el controlador introduce un cero además de una ganancia, por lo que, en
efecto, el circuito corresponde a un controlador PD.
El controlador proporcional-integral (PI) tiene como salida una señal que es proporcional al
error y a la integral del error. Como tal introduce un cero y un polo en el origen. El circuito con
amplificadores operacionales que responde a este comportamiento es el siguiente:
R2 C2
R4
R1
R3
-
-
+
ei +
e0
La función de transferencia resultante muestra un cero con un polo en el origen, por lo que el
circuito representa un controlador PI.
Este tipo de controlador combina las tres acciones básicas de control, introduciendo dos ceros
y un polo en el origen. El circuito que implementa este controlador es la combinación de los que se
vieron antes:
C1 R2 C2
R4
R1
R3
-
-
+
ei +
e0
Comparando términos:
1 K p 1 Ti s Td Ti s 2
K p 1 Td s
Ti s Ti s
R R C R2C2
Kp 4 1 1
R1 R3C2
Ti R1C1 R2C2
R1 R2C1C2
Td
R1C1 R2C2
C2
C1
R4
R2
R1
R3
-
-
+
ei +
e0
1
s
E0 ( s ) R4C1 R1C1
Ei ( s ) R 3 C2 s 1
R2C2
1
s
T R2C2 R4C1
Kc T R1C1 Kc
1 R1C1 R3C2
s
T
Dado que el tipo de red está determinado por el valor del parámetro , éstas serán las
condiciones que se habrán de cumplir para cada tipo de red:
Veremos en este epígrafe una serie de ejemplos de aplicación de cada uno de los
controladores vistos con anterioridad.
6.2.6.1 Control P
Ejemplo 1:
Supongamos un proceso genérico que tiene la siguiente función de transferencia:
K
G p ( s)
Ts 1
C ( s) K
R( s ) Ts 1 K
Si trazamos el diagrama lugar de las raíces para este sistema encontramos la siguiente
representación:
0.8
0.6
0.4
0.2
Imag Axis
-0.2
-0.4
-0.6
-0.8
-1
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
Real Axis
1 1 1
ess
1 K p 1 lim G p ( s ) 1 K
s 0
C (s) KK p
R( s ) Ts 1 KK p
Con lo que el polo en lazo cerrado se desplaza más a la izquierda (disminución del tiempo de
KK p
asentamiento) y, además, la ganancia del sistema nos queda , expresión que tiende a la
1 KK p
unidad cuanto mayor sea el producto KKp.
Ejemplo 2:
1
G p ( s)
s2
1.5
0.5
Imag Axis
-0.5
-1
-1.5
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
Real Axis
Ejemplo 3:
K
G p ( s)
s (Ts 1)
El diagrama lugar de las raíces, para unos valores concretos de K y T, es el siguiente:
1.5
0.5
Imag Axis
-0.5
-1
-1.5
-2
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
Real Axis
Al tratarse de un sistema tipo 1, sabemos que el error de posición va a ser cero, por lo que el
control proporcional en este caso no va a mejorar aún más ese aspecto, pero sí en lo que respecta al
error de velocidad (que es constante en este caso), de modo que un valor mayor de la ganancia da
como resultado un error ante entrada rampa menor.
Por otro lado, la variación de la ganancia que introduce el control proporcional puede hacer
que la respuesta transitoria del sistema sea sustancialmente distinta, de forma que se puede
conseguir que los polos en lazo cerrado sean reales (respuesta sobreamortiguada) o que sean
Conclusiones:
De los ejemplos anteriores se pueden extraer las siguientes conclusiones más o menos
generales aunque, claro está, todo depende de cada caso concreto:
El control proporcional puede ser adecuado para procesos que incorporen un integrador
puro (polo en el origen, sistemas de tipo 1).
Dado un sistema de tipo n, con control P la salida del sistema nunca alcanzará la
referencia ante una entrada de orden n.
6.2.6.2 Control PI
Por lo que introduce a la función de transferencia en lazo abierto del sistema un polo en el
1
origen y un cero en , incrementando en 1 el orden del sistema global y teniendo, en
Ti
consecuencia, una influencia evidente sobre el error estacionario: convierte el error de posición
constante de un sistema tipo 0 a cero (siempre que sea estable), aparte de la incidencia en sistemas
de órdenes superiores.
Sin embargo, la introducción del polo en el origen puede inestabilizar un sistema, por lo que
los parámetos Kp y Ti tienen que elegirse adecuadamente (e incluso puede ser que la introducción del
PI sea totalmente inapropiada). El valor de Ti permite fijar la posición del cero modificando, en
Kp
consecuencia, el lugar de las raíces del sistema original, mientras que el cociente determinará la
Ti
posición de los polos en lazo cerrado en el nuevo lugar de las raíces (Kp será el parámetro variable).
Ejemplo:
Volvamos al caso anterior del sistema que tenía esta función de transferencia:
K
G p (s)
Ts 1
KK p (Ti s 1)
G p (s)
Ti s (Ts 1)
El sistema ha pasado a ser de tipo de 1, por lo que se ha anulado el error de posición, pero
ahora, en función del valor de Ti el diagrama lugar de las raíces variará sustancialmente. Empleando
la herramienta “rltool” de MATLAB y editando el compensador podemos encontrar dos formas muy
diferentes del lugar de las raíces modificado:
0.8
0.6
0.4
0.2
Imag Axis
-0.2
-0.4
-0.6
-0.8
-1
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
Real Axis
Caso de Ti > T
0.8
0.6
0.4
0.2
Imag Axis
-0.2
-0.4
-0.6
-0.8
-1
-2 -1.5 -1 -0.5 0 0.5 1
Real Axis
Caso de Ti < T
Observamos que en el caso de Ti > T los dos polos en lazo cerrado son reales para cualquier
valor de la ganancia, existiendo uno de ellos como dominante. Por lo tanto, la respuesta alcanzará la
referencia en la forma de un sistema de segundo orden sobreamortiguado.
En el caso de Ti < T, existen valores de la ganancia para los que los polos en lazo cerrado son
complejos conjugados, dando lugar a respuestas oscilatorias ante entrada escalón.
En conclusión, en función de cuáles sean las especificaciones que se deseen, habrá que
ajustar el valor de las constantes Ti y Kp.
6.2.6.3 Control PD
U ( s)
K p (1 Td s )
E ( s)
1
Este controlador introduce un cero en en la función de transferencia en lazo abierto del
Td
sistema, por lo que no tiene influencia sobre el error estacionario. Sin embargo, esta modificación,
dependiendo del caso concreto, puede hacer el sistema más estable.
El valor de la constante Td fija la posición del cero en el lugar de las raíces, modificando el del
sistema original. El parámetro en este nuevo lugar geométrico corresponde a Kp, que permitirá tratar
de imponer las especificaciones de respuesta transitoria y estacionaria.
Ejemplo:
1
G p (s)
s2
Este sistema presenta el problema de que no estable, puesto que los polos en lazo cerrado
siempre se encuentran sobre el eje imaginario.
0.8
0.6
0.4
0.2
Imag Axis
-0.2
-0.4
-0.6
-0.8
-1
-2 -1.5 -1 -0.5 0 0.5 1
Real Axis
Con lo que se ha conseguido estabilizar el sistema llevando los dos polos en lazo cerrado
resultantes al semiplano izquierdo del plano complejo para cualquier valor de la ganancia. Además se
verifica que la introducción de un cero desplaza hacia la izquierda el lugar de las raíces original.
Además, dada la forma del lugar de las raíces, se observa que se puede conseguir que la
respuesta sea oscilatoria (sistema subamortiguado) o parecida a la de un sistema de primer orden
(muy amortiguada), en función de la elección adecuada de Kp y Td.
U ( s) 1 K (1 Ti s TT 2
i ds )
K p (1 Td s ) p
E ( s) Ti s Ti s
De aquí se desprende que el controlador PID introduce dos ceros y un polo en lazo abierto
(tener presente que los dos ceros aparecerán también en el lazo cerrado, con las implicaciones en la
respuesta que ello conlleva) y que ahora es necesario el ajuste de tres parámetros.
Los dos ceros mencionados siempre se encontrarán en el semiplano izquierdo del plano
complejo mientras Ti y Td sean positivas, pero serán de tipos distintos en función de la relación entre
ellas:
Ti > 4Td ceros reales y distintos.
Ti = 4Td ceros reales múltiples.
Ti < 4Td ceros complejos conjugados.
En algún caso y para reducir la simplicidad, suele suponerse que Ti = 4Td, con lo que se evita
el ajuste de uno de los parámetros y se obtiene un cero doble. Una vez que se tenga una forma
concreta del lugar de las raíces modificado se tomará como parámetro Kp para tratar de cumplir las
especificaciones de respuesta transitoria y estacionaria.
1
s
G ( s) aT , con a > 1
1
s
T
1 1
Esta red introduce un cero (en ) y un polo (en ) en la función de transferencia en
aT T
lazo abierto estando siempre el primero a la derecha del segundo.
Normalmente, la red de adelanto se emplea junto con un bloque de control proporcional que
servirá de parámetro para la modificación de la respuesta transitoria, una vez fijado el lugar de las
raíces del sistema global. En concreto, se persiguen varios objetivos básicos: estabilización del sistema
en el caso de que originalmente no sea estable y, en caso contrario, la disminución del tiempo de
asentamiento y del sobrepaso máximo, mejorando la estabilidad relativa. Para ello, los polos
dominantes en lazo cerrado deben cumplir con las especificaciones y, como en cualquier otro caso,
ser lo suficientemente dominantes como para garantizar los requerimientos deseados.
Por otro lado, hay que tener en cuenta también los ceros que pudieran aparecer en lazo
cerrado, ya que si son dominantes (están muy cercanos al eje imaginario) variarán sustancialmente la
respuesta en comparación con la inicialmente prevista con la simple ubicación de los polos en lazo
cerrado.
Ejemplo:
K
G p ( s)
s (Ts 1)
2.5
1.5
0.5
Imag Axis
-0.5
-1
-1.5
-2
-2.5
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2
Real Axis
La consecuencia más inmediata es que, con el valor de ganancia apropiado, se pueden seguir
cumpliendo las especificaciones en términos de sobrepaso máximo, pero además puede reducirse el
tiempo de asentamiento al doblar hacia la izquierda las ramas del lugar de las raíces. Concretamente,
cuanto más desplacemos a la izquierda el polo de la red más a la izquierda se doblarán las ramas y
mejoraremos la estabilidad relativa y el tiempo de asentamiento. Así mismo, cuanto más a la derecha
desplacemos el cero de la red se producirá el mismo efecto, pero quizá con el inconveniente de que
los dos polos más cercanos al eje imaginario sean menos dominantes y se modifique la respuesta
deseada. En este ejemplo, se producen estos efectos con la red por la modificación de las
coordenadas del centroide de las asíntotas.
1
s
G ( s) K c aT , con 0 < a < 1
1
s
T
Al igual que en el caso anterior, también se introducen un cero y un polo pero con la
diferencia de que ahora el polo se encuentra a la derecha del cero. En la mayoría de las ocasiones,
esto suele tender a empeorar la estabilidad del sistema, por lo que normalmente se ubican el polo y el
cero muy próximos entre sí y cercanos al eje imaginario.
Si se coloca una red de retraso en serie con la planta con el polo y el cero muy cercanos y se
persigue que el transitorio no se modifique (es decir, que la ubicación de los polos en lazo cerrado sea
la original), es preciso que la ganancia que aporte la red sea aproximadamente la unidad.
Se pretende que el cero y el polo de la red se coloquen muy cerca uno del otro y, a su vez,
muy cercanos al eje imaginario. Para ello, es preciso que T adquiera un valor elevado y, además,
como el cero y el polo han de estar cercanos y con valor pequeño (cerca del eje imaginario), se puede
hacer el valor de a bastante pequeño y aún 1/aT y 1/T tendrán un valor similar. No obstante, tampoco
se puede elegir el valor de T arbitrariamente muy grande, puesto que eso puede provocar que la red
no sea realizable físicamente.
Como se desea que el estacionario no se modifique es preciso que la ganancia del lugar de las
raíces permanezca casi inalterada (para que los polos en lazo cerrado se ubiquen en casi las mismas
posiciones iniciales) y es por eso por lo que se hace la ganancia Kc de la red aproximadamente la
unidad.
Por otro lado, se puede demostrar que las constantes de error estático se incrementan en el
factor Kc/a, por lo que si a se ha elegido pequeña, se mejora ostensiblemente el error en estado
estacionario.
Todo lo anterior razona la necesidad de que el polo y el cero se encuentren muy próximos y
cercanos al eje imaginario.
Ejemplo:
K
G p ( s)
s (Ts 1)
Si introducimos la red de retardo como controlador el lugar de las raíces nos queda de esta
forma:
0.6
0.4
0.2
Imag Axis
-0.2
-0.4
-0.6
El lugar de las raíces apenas se ha modificado pero ahora, para ubicar los polos en lazo
cerrado en una misma posición del lugar geométrico original, se necesita un valor de ganancia
superior, con la repercusión que esto tiene en el error de velocidad, cumpliendo así como el
requerimiento que deseado.
En este epígrafe se analizarán las modificaciones que producen los controladores vistos con
anterioridad en la respuesta en frecuencia del sistema. El estudio permitirá reafirmar algunas de las
conclusiones a las que se llegaron con la compensación en el lugar de las raíces pero
fundamentalmente permitirá establecer técnicas alternativas de diseño en el dominio de la frecuencia.
Para ello emplearemos los diagrama de Bode y polar.
La zona de frecuencias medias, que se encuentran situadas en las cercanías del punto crítico
(en la traza de Nyquist el punto -1+j0 y en el diagrama de Bode en la amplitud 0 dB y fase –180º),
nos ofrecen información acerca de la estabilidad relativa del sistema, ya que en estas zonas es donde
se mide el margen de ganancia y el margen de fase. Como se vio en el tema anterior, cuanto mayor
sea el valor de estos dos parámetros y positivos, mayor será la estabilidad del sistema.
Finalmente, la zona de frecuencias altas nos indica la medida del rechazo al ruido del sistema.
Si se desea un sistema poco sensible a las señales de alta frecuencia la ganancia en esta zona debe
ser pequeña.
Se analizará aquí la influencia de cada uno de los tipos de controladores vistos con
anterioridad en la respuesta en frecuencia del sistema global después de la introducción en el lazo de
cada uno de ellos.
6.3.1.1 Control P
6.3.1.2 Control PI
Bode Diagrams
From: U(1)
60
40
Phase (deg); Magnitude (dB)
20
-20
-20
-40
To: Y(1)
-60
-80
-100
10-2 10-1 100 101
1
Frequency (rad/sec)
Ti
Se observa que, a bajas frecuencias, la amplitud es infinita, lo que se traduce en que el error
de posición es nulo (ver tema anterior para determinar en Bode la constante de error de posición Kp).
1
Para frecuencias próximas a (que es uno de los parámetros del controlador) se observa
Ti
que se producirá una disminución de la fase del sistema global (FTLA), ya que la curva
correspondiente de la función de transferencia del controlador se va aproximando cada vez más a
cero. Esto puede provocar una reducción importante del MG y del MF, reduciendo la estabilidad
relativa del sistema o incluso llegando a inestabilizarlo. El riesgo de hacer inestable el sistema, en
general, es inversamente proporcional a Ti (mayor cuanto menor sea Ti).
6.3.1.3 Control PD
Bode Diagrams
From: U(1)
30
25
Phase (deg); Magnitude (dB)
20
15
10
80
60
To: Y(1)
40
20
0
10-1 100 101
1
Frequency (rad/sec)
Td
1
Para frecuencias próximas o superiores a este controlador introduce un incremento en la
Td
fase de la función de transferencia global que puede producir, en general y dependiendo del caso
concreto, un aumento del MG y del MF, con lo que se mejora la estabilidad relativa del sistema de
control. Como consejo práctico para el diseño del PD en frecuencia es conveniente tratar de que la
1
que frecuencia en la que se mide el MF esté próxima a la frecuencia .
Td
El mayor inconveniente que presenta este tipo de controlador es que a frecuencias altas la
amplificación de las señales de entrada es muy grande (véase la curva de magnitud), con lo que se
obtiene una alta sensibilidad al ruido, siendo éste un inconveniente importante en el trabajo con
sistemas reales.
Bode Diagrams
From: U(1)
0
-10
Phase (deg); Magnitude (dB)
-20
-30
-40
80
60
To: Y(1)
40
20
0
10-2 10-1 100 101 102
1 Frequency (rad/sec) 1
aT T
La principal función de una red de adelanto es incrementar la fase del sistema global en las
proximidades de la frecuencia en la que se mide el MF, con el objetivo de mejorar éste.
Generalmente, una mejora del MF suele traer consigo una mejora del MG.
Bode Diagrams
From: U(1)
40
30
Phase (deg); Magnitude (dB)
20
10
-20
To: Y(1)
-40
-60
-80
10-2 10-1 100 101 102
1 Frequency (rad/sec) 1
T aT
Como en el caso anterior, el nombre de red de retardo está plenamente justificado cuando se
observa la curva de fase (siempre es negativo, por lo que la salida está retrasada respecto de la
entrada senoidal). La curva de fase presenta un mínimo m a la frecuencia m. Las expresiones
anteriores para la red de adelanto también son válidas para la red de retardo, salvo que ahora
obtendremos un mínimo de fase en lugar de un máximo.
Esta red también introduce una amplificación a bajas frecuencias que depende de a,
concretamente, 20 log a , que se traduce en una disminución del error estacionario, siendo ésta la