Control Automatico 2013 - High
Control Automatico 2013 - High
El libro es especialmente til para estudiantes y profesores en las carreras de Ingeniera Elctrica,
Electrnica, Mecnica, Mecatrnica e incluso Robtica.
Victor Manuel Hernndez Guzmn
Ramn Silva Ortigoza
Roberto Valentn Carrillo Serrano
COLECCIN CIDETEC
Roberto Valentn Carrillo Serrano. Naci
Victor Manuel Hernndez Guzmn. Naci
en Quertaro, Qro., Mxico en 1976. Recibi
en Quertaro, Qro., Mxico. Recibi el ttulo
el ttulo de Ingeniero en Instrumentacin y
de Ingeniero Industrial en Elctrica por parte
Control de Procesos por la Universidad
del Instituto Tecnolgico de Quertaro en
Autnoma de Quertaro, donde tambin
1988, el grado de Maestro en Ciencias en
obtuvo los grados de Maestro en Ciencias en
Ingeniera Elctrica (Control) por parte del
Instrumentacin y Control Automtico, y de
Instituto Tecnolgico de la Laguna, en
Doctor en Ingeniera en 2000, 2008 y 2011
Torren, Coah., en 1991 y el grado de Doctor
respectivamente. Trabaj en Kellogg de
en Ciencias en Ingeniera Elctrica
Mxico de 1999 a 2006. Las reas de inters del
(Mecatrnica) por parte del CINVESTAV-
Dr. Carrillo Serrano son los robots
IPN, en Mxico, D.F., en 2003. Actualmente
manipuladores, control de mquinas elctricas
es Profesor en los programas de Licenciatura y
y control de sistemas mecatrnicos.
Maestra en Instrumentacin y Control de la
Actualmente es Profesor en el programa de
Universidad Autnoma de Quertaro. Su
Licenciatura en Ingeniera en Automatizacin
trabajo de investigacin trata sobre el control
y de la Maestra en Instrumentacin y Control
de robots manipuladores y sistemas electro-
Automtico de la Universidad Autnoma de
mecnicos. Tiene particular inters en la
Quertaro, adems de ser miembro del
construccin de prototipos didcticos para la
Sistema Nacional de Investigadores, Mxico.
enseanza de tcnicas de control clsicas y
En cuanto a publicaciones se refiere, estas se
modernas (no lineales).
encuentran en estndares internacionales que
pertenecen a la base de datos Journal Citation
Reports (JCR) e ISI-Thomson.
CONTROL AUTOMATICO
Autores:
Victor Manuel Hernandez Guzm an
Ramon Silva Ortigoza
Roberto Valentn Carrillo Serrano
CIDETEC
Av. Juan de Dios Batiz S.N. esq. Miguel Oth
on de Mendiz
abal
Col. Nueva Industrial Vallejo
Del. Gustavo A. Madero C.P. 07700
ISBN: 978-607-414-362-1
Hecho en Mexico / Made in Mexico
V
Enero de 2013.
Indice general
1. Introducci on . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.1. Sistema de control de un ca non antiaereo . . . . . . . . . . . . . . . . . . . 2
1.2. Historia del Control Automatico . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.3. Prototipos didacticos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.4. Resumen del captulo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.5. Preguntas de repaso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
Referencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
Referencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
Referencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 175
Referencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 219
Indice general XV
Referencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 289
Referencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 391
Referencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 445
Indice general XVII
Referencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 497
Referencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 529
Referencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 573
Referencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 601
Referencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 643
Referencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 667
Referencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 701
Indice alfab
etico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 703
1
Introducci
on
can
on
radar antiaereo
d
motor
v engranes
d
controlador amplificador
de
potencia
d
d
Figura 1.2. El ca
non antiaereo siempre debe moverse hacia la orientaci
on deseada.
v = kp (d ) (1.1)
amplificador
d
+ controlador de v motor
potencia can
on
kp
1
0
0
tiempo
T(t)
l
g
m
Una vez que se ha descrito a groso modo que es lo que se persigue al dise nar
un sistema de control as como algunos conceptos u tiles para conseguirlo, a
continuacion se presenta una breve historia del Control Automatico. El obje-
tivo es que el lector se de cuenta de como han sido formulados estos conceptos
y, ademas, pueda apreciar cuales han sido las necesidades ingenieriles que han
motivado estas ideas. Tambien se indican las partes del presente libro en don-
de se abordan los principales conceptos y tecnicas de los sistemas de control.
El contenido de esta seccion esta basado en la informacion reportada en [1].
El Control Automatico se ha usado desde hace mas de 2000 a nos. Se tiene
conocimiento de la existencia de relojes de agua, construidos por Ktesibios,
hacia el ano 270 antes de Cristo, as como una variedad de mecanismos in-
geniosos construidos en Alejandra y descritos por Heron. Sin embargo, des-
de el punto de vista de la ingeniera, el avance mas importante en Control
Automatico se debio a James Watt en 1789 al introducir un regulador de ve-
locidad para su maquina de vapor. Sin embargo, a pesar de su importancia, el
regulador de velocidad de Watt tena varios problemas. Se observaba que en
ocasiones la velocidad variaba de manera oscilatoria en lugar de permanecer
constante o creca sin lmite. Tratando de determinar bajo que condiciones se
poda asegurar un funcionamiento estable, entre 1826 y 1851 J.V. Poncelet
y G.B. Airy mostraron que era posible utilizar ecuaciones diferenciales para
representar el funcionamiento completo de la maquina de vapor junto con el
regulador de velocidad (vease el captulo 2 para el problema del modelado de
sistemas fsicos).
En esas fechas los matematicos saban que la estabilidad de una ecuacion
diferencial estaba determinada por la ubicacion de las races de la ecuacion
caracterstica correspondiente y se saba que haba inestabilidad si la parte
real de una raz era positiva (vease el captulo 3, seccion 3.8). Sin embargo, no
era sencillo calcular el valor numerico de dichas races y en ocasiones ni siquie-
ra era posible. As, en 1868 J.C. Maxwell mostro la manera de determinar la
estabilidad de las maquinas de vapor con el regulador de Watt simplemente
examinando los coeficientes de la ecuacion diferencial que la representa. Sin
embargo, este resultado solo era u til para ecuaciones diferenciales de segundo,
1.2 Historia del Control Autom
atico 9
tercero y cuarto orden. Mas tarde, entre 1877 y 1895 y de manera indepen-
diente, E.J. Routh y A. Hurwitz dedujeron un metodo para determinar la
estabilidad de sistemas de cualquier orden, resolviendo el problema que Max-
well haba dejado abierto. Este metodo ahora es conocido como el Criterio de
Routh o de Routh-Hurwitz (vease el captulo 4, seccion 4.3).
Durante gran parte del siglo XIX se desarrollaron muchas aplicaciones
relacionadas con el Control Automatico entre las cuales estaban el control
de temperatura, de presion, de nivel de lquidos y la velocidad de maquinas
rotativas. Por otro lado, se empezo a usar el vapor para mover grandes ca nones
y para actuar sobre el sistema de direccion de barcos cada vez mas grandes.
Incluso, por esa epoca, en Francia se introdujeron los terminos de servomotor
y servomecanismo para describir un movimiento generado por un servidor o
esclavo. Sin embargo, la mayora de los controladores de ese entonces eran
del tipo encendido-apagado y fueron personas como E. Sperry y M.E. Leeds
quienes se dieron cuenta de que los mejores operadores humanos empleaban
el sentido de anticipacion disminuyendo la potencia conforme la variable a
controlar se acercaba a su valor deseado. Fue N. Minorsky quien en 1922
presento un analisis claro de los sistemas de control de posicion y formulo lo
que hoy en da se conoce como controlador PID (vease el captulo 5, seccion
5.2.5). Este controlador fue deducido observando la manera en que el piloto
humano de un barco controla su direccion.
Por otro lado, desde 1920 la amplificacion haba producido muchos proble-
mas a las compa nas telefonicas ya que se distorsionaba fuertemente la se nal
de audio. Fue en esa epoca que H.S. Black encontro que si una peque na can-
tidad de la senal obtenida a la salida de un amplificador se utilizaba para ser
realimentada a la entrada del mismo se poda reducir la distorsion producida
por el amplificador. Durante el desarrollo de este trabajo, Black fue ayudado
por H. Nyquist quien, a partir de esas experiencias, publico en 1932 un trabajo
titulado Regeneration Theory en donde establecio las bases de lo que ahora
es conocido como el Analisis de Nyquist (vease el captulo 6, seccion 6.5).
Durante el perodo 1935-1940, las compa nas telefonicas deseaban ampliar
el ancho de banda de sus sistemas de comunicacion para aumentar el n umero
de sus usuarios. Para esto necesitaban que sus lneas telefonicas presentaran
una buena caracterstica de respuesta en frecuencia (vease el captulo 6): una
ganancia constante sobre un amplio rango de frecuencias con un peque no
angulo de atraso y una aguda pendiente de atenuacion a partir de una deter-
minada frecuencia de corte. Motivado por este problema, H. Bode estudio la
relacion existente entre una caracterstica de atenuacion dada y el mnimo
corrimiento de fase que se le puede asociar. Como resultado introdujo los con-
ceptos de margen de ganancia y margen de fase (vease el captulo 6, seccion
6.6) y empezo a manejar el punto (1, 0) del plano complejo como un punto
crtico en lugar del punto (+1, 0) manejado por Nyquist. Detalles completos
del trabajo de Bode aparecieron en 1945 en su libro Network Analysis and
Feedback Amplifier Design.
10 1 Introducci
on
Pero fue la Segunda Guerra Mundial la que hizo que el trabajo en sistemas
de control se concentrara en unos pocos problemas especficos. El mas impor-
tante de estos fue el relacionado con el direccionamiento de ca nones antiaereos.
El trabajo en este problema motivo el desarrollo de nuevas ideas en el control
de servomecanismos. G. S. Brown del Instituto Tecnologico de Massachusetts
mostro que muchos sistemas electricos y mecanicos pueden ser representados
y manipulados usando diagramas de bloques (vease el captulo 4, seccion 4.1)
y A. C. Hall mostro en 1943 que manejando los bloques como funciones de
transferencia (vease el captulo 3) poda obtenerse la funcion de transferencia
del sistema completo para finalmente usar el criterio de estabilidad de Nyquist
y determinar los margenes de ganancia y de fase.
Los investigadores en el Instituto Tecnologico de Massachusetts usaron
circuitos de adelanto (vease el captulo 10, seccion 10.2.2) en el trayecto directo
para modificar el desempe no del sistema de control, mientras que en el Reino
Unido se usaron varios lazos internos para modificar la respuesta del sistema
de control.
Hacia el final de la Segunda Guerra Mundial las tecnicas de respuesta en
frecuencia basadas en el metodo de Nyquist y las graficas de Bode ya estaban
bien establecidas, describiendo el desempe no del sistema de control en termi-
nos de ancho de banda, frecuencia de resonancia, margen de fase y margen
de ganancia (vease el captulo 6). El enfoque alternativo a estas tecnicas se
basaba en la solucion de las ecuaciones diferenciales usando la Transforma-
da de Laplace y describan el desempe no del sistema de control en terminos
del tiempo de subida, sobre paso, error en estado estacionario y el amorti-
guamiento (vease el captulo 3, seccion 3.3). Muchos ingenieros preferan este
metodo porque los resultados estaban expresados en terminos reales. Pero
este enfoque tena la desventaja de que no exista una manera sencilla que
permitiera al disenador relacionar los cambios en los parametros con cambios
en la manera en que responda el sistema. Fue precisamente el metodo del
Lugar de las Races (vease el captulo 5, secciones 5.1 y 5.2) introducido en
1948 y 1950 por W. Evans el que permitio librar estos obstaculos. As, hacia
esas fechas las ahora llamadas tecnicas de control clasico estaban bien estable-
cidas y estaban orientadas a sistemas que podan ser descritos por ecuaciones
diferenciales lineales con coeficientes constantes y con una sola entrada.
Entonces vino la era de los vuelos supersonicos y espaciales. Era necesario
utilizar modelos fsicos detallados representables en ecuaciones diferenciales
que podan ser lineales o no lineales. Los ingenieros que trabajaban en las
industrias aeroespaciales encontraron que siguiendo las ideas de Poincare era
posible formular ecuaciones diferenciales generales en terminos de un conjunto
de ecuaciones diferenciales de primer orden y as empezo a nacer lo que hoy
se conoce como la tecnica de las variables de estado (vease el captulo 7).
Un gran impulsor de esta tecnica fue R. Kalman quien, alrededor de 1960,
presento los conceptos de controlabilidad y observabilidad (vease la seccion
7.6).
1.3 Prototipos did
acticos 11
Finalmente, se debe decir que a partir de esas fechas se han detectado nue-
vos problemas en los sistemas de control que han motivado la introduccion
de diversas tecnicas de control que a un hoy en da se siguen desarrollando.
Por ejemplo, las no linealidades encontradas en los servomecanismos ha mo-
tivado el desarrollo de las tecnicas de control no lineal. El control de aviones
supersonicos, que deben operar bajo amplios rangos de temperatura, presion,
velocidad, etc., motivo el desarrollo de control adaptable. La introduccion
de sistemas de control basados en radar motivo el desarrollo de tecnicas de
control para sistemas en tiempo discreto, etc.
Darse cuenta de que el modelo matem atico de los sistemas fsicos est
a cons-
tituido por ecuaciones diferenciales.
Aprender a obtener el modelo matem atico de sistemas mec anicos, electri-
cos y electromecanicos.
Identificar a los sistemas fsicos como la interconexion de procesadores de
energa.
Identificar las analogas entre sistemas de diferente naturaleza, a partir de
la manera en que procesan la energa.
Tradicionalmente, cuando se aborda el modelado de sistemas fsicos se re-
curre a los metodos especficos desarrollados para cada tipo de sistema. Es
decir, los sistemas electricos son modelados usando los metodos desarrollados
en la teora de circuitos electricos y cuando se trata de sistemas mecanicos
se usan los metodos desarrollados en ingeniera mecanica. Sin embargo, pro-
ceder de esta manera tiene varias desventajas: i) la persona no especializada
en el tema en cuestion debe simplemente aceptar (sin muchas explicaciones)
los metodos de modelado desarrollados en esa rama del conocimiento, ii) el
estudiante no se da cuenta de que las leyes fsicas utilizadas para el modelado
de sistemas de diferente naturaleza estan relacionadas por analogas; aunque
en todo curso de control automatico siempre se tratan de resaltar las ana-
logas, estas tienen que ser explicadas a partir de las ecuaciones resultantes
del modelado y no se puede ver que es la esencia del fenomeno la que tiene su
contraparte en sistemas de diferente naturaleza.
Una manera de estudiar el modelado de sistemas fsicos de diferente natu-
raleza bajo una perspectiva unificada es utilizando un concepto que es com un
a varios ambitos de la ingeniera: la energa del sistema. Por esta misma razon
en esta obra el modelado de sistemas se limita a sistemas electricos y mecani-
cos los cuales pueden estudiarse usando conceptos generales de energa que
son comunes a todos estos sistemas. La idea fundamental del uso de la energa
para el modelado es que los componentes de cada sistema bajo estudio suminis-
tran, almacenan o disipan energa y cuando se unen para formar un sistema
complejo es cuestion de determinar como esa energa es transmitida de un
componente a otro.
Este enfoque de usar la energa para el modelado (y, por tanto, el presente
captulo) esta basado en las ideas introducidas en [1].
f
Sistema Sistema
1 e 2
puerto
w = ef
Para cada uno de los componentes del sistema se debe definir una convencion
de signos para las variables generalizadas de esfuerzo y de flujo (vease la figura
2.2). De este modo, si la potencia w tiene signo positivo para el componente i
entonces este componente recibe energa en ese instante. Si la potencia w tiene
signo negativo para el elemento i entonces el componente entrega energa a
los otros componentes del sistema en ese instante.
La energa intercambiada (E) en el intervalo de tiempo [0, t] esta dada
como la integral de la potencia instantanea:
18 2 Modelado matem
atico de sistemas fsicos
Componente
e
i
Figura 2.2. Sentidos definidos como positivos para las variables de esfuerzo y flujo.
Z t
E= ef dt (2.1)
0
e dt = dea , f dt = dfa
f = (ea )
Por otro lado, para poder calcular la integral en (2.3) es necesario que el
esfuerzo e se pueda escribir como funcion de la acumulacion de flujo fa , es
decir que se pueda escribir:
e = (fa )
Almacenamiento de flujo
F = ma (2.4)
p = mv12 (2.5)
es muy conveniente porque se puede integrar (2.5) y (2.4) para concluir que el
momentum es la acumulacion de flujo (p = fa ), es decir que se puede escribir:
Z t
p= F dt + p(0)
0
m F
2
Almacenamiento de esfuerzo
v 12
1 2
F F
v1 v2
Z t
x12 = v12 dt + x12 (0) (2.6)
0
22 2 Modelado matem
atico de sistemas fsicos
Almacenamiento de flujo
T ! 12
2
Eje de rotacin
T = I (2.8)
h = I12 (2.9)
Almacenamiento de esfuerzo
Z t
12 = 12 dt + 12 (0) (2.10)
0
T = k 12 , (f = (ea )) (2.11)
2.2.3. Sistemas el
ectricos
Almacenamiento de esfuerzo
En la figura 2.7 se muestra un inductor a traves del cual circula una corrien-
te electrica i (i = f , flujo) bajo el efecto de un voltaje v12 (v12 = e, esfuerzo)
aplicado entre sus extremos. Se supone que no existen efectos parasitos, es
decir, el inductor no tiene resistencia electrica interna ni existe capacitancia
entre sus espiras. Los sentidos mostrados para la corriente electrica y el voltaje
son los sentidos que se definen como positivos para este tipo de componente
de sistema.
Faraday fue el primero en darse cuenta de que un inductor (y en general
cualquier circuito electrico suficientemente largo, dado que un inductor es un
conductor electrico arrollado sobre un n ucleo) tiene propiedades analogas a
las que tiene el momentum en sistemas mecanicos. Lo que Faraday llamo el
! 12
T !1 !2 T
1 2
i
1 2
+
v 12
= Li (2.12)
Almacenamiento de flujo
i
1 2
+
v 12
e = (f )
w = ef (2.17)
v12 = (F )
donde v12 (e, esfuerzo) es la velocidad relativa entre los cuerpos y F (f , flujo)
es la fuerza aplicada. Un caso particular muy importante de friccion es la
friccion viscosa, la cual esta representada por una funcion constitutiva lineal
de la forma:
1
v12 = F, friccion viscosa (2.18)
b
donde b es una constante positiva conocida como el coeficiente de friccion
viscosa. Este tipo de friccion ocurre, por ejemplo, cuando una placa plana
y llena de orificios se desplaza dentro de una camara cerrada llena de aire,
como se muestra en la figura 2.9. Debido a que el aire debe fluir a traves de
los orificios conforme la placa se mueve a velocidad v12 , es necesario aplicar
una fuerza F para mantener la velocidad del movimiento. Ademas, v12 =
v1 v2 donde v1 y v2 son las velocidades de los extremos 1 y 2 del dispositivo,
28 2 Modelado matem
atico de sistemas fsicos
v 12
v1 F
F
1
2
v2
F F F
Fs Fc
b
v 12 v 12 v 12
Fs Fc
estacionario diferente de cero. Por otro lado, dado que las fricciones estatica
y de Coulomb tienen funciones constitutivas no lineales y discontinuas, no
pueden ser manejadas usando tecnicas de control por sistemas lineales ni por
las tecnicas de control para sistemas no lineales suaves. Esta es la razon
por las que con frecuencia no son consideradas en el modelado de sistemas.
Sin embargo, hay que tener presente que en cualquier situacion experimental
se observaran los efectos de estas dos fricciones. Se sugiere consultar [4] si se
desea conocer algunos metodos para medir los parametros de distintos tipos
de friccion.
2.3.3. Sistemas el
ectricos
T !1 !2 T
! 12
y esta definida por la Ley de Ohm [2], pag. 530, en el caso de que el disipador
sea lineal:
v12 = iR (2.22)
v 12
i
1 2
+
f
e
e1
+
e1 Recibe Entrega
potencia potencia
f
(a)
f
f1
+
f1 Recibe Entrega
e potencia potencia
e
(b)
Una fuente de fuerza es una fuente de flujo, por lo que la fuente entrega
una fuerza cuyo valor es independiente de la velocidad (esfuerzo) del compo-
nente del sistema sobre el cual se esta aplicando dicha fuerza. Una fuente de
velocidad es una fuente de esfuerzo, por lo que la velocidad que produce es
independiente de la fuerza (flujo) ejercida sobre el componente del sistema
correspondiente. Claramente es mucho mas facil visualizar lo que es una fuen-
te de fuerza que una fuente de velocidad. Sin embargo, algunos fenomenos
pueden ser modelados de manera conveniente como fuentes de velocidad.
2.4.3. Sistemas el
ectricos
Una fuente de corriente electrica es una fuente de flujo, por lo que la fuente
entrega una corriente electrica cuyo valor es independiente del voltaje (esfuer-
zo) entre sus terminales. Una fuente de voltaje es una fuente de esfuerzo, por
lo que el voltaje entre sus terminales es independiente de la corriente electrica
a traves de ella. Tal como sucede en los sistemas mecanicos, es mas facil visua-
lizar una clase de fuentes: el concepto de fuente de voltaje es mas claro que el
de una fuente de corriente electrica. Sin embargo, el uso de fuentes de corrien-
te es un artificio muy u til para modelar algunos efectos que se presentan, por
ejemplo, en circuitos electronicos.
32 2 Modelado matem
atico de sistemas fsicos
2.5. Convertidores
Los convertidores son componentes que simplemente permiten que la
energa fluya a traves de ellos sin almacenarla ni disiparla. Estos componentes
estan conectados a traves de dos puertos (vese la figura 2.14). Cada puerto
esta definido por una variable generalizada de esfuerzo y una de flujo. El pa-
pel de los convertidores es modificar las variables de esfuerzo y de flujo en
el segundo puerto respecto de las variables de esfuerzo y de flujo que existen
en el primer puerto. Este proceso es realizado de manera que la energa en el
segundo puerto es igual a la energa en el primer puerto. Esos componentes
se pueden clasificar en dos tipos: i) los transformadores, en los cuales las va-
riables de esfuerzo y de flujo en ambos puertos son de la misma naturaleza e
ii) los adaptadores, en lo cuales las variables de esfuerzo y flujo en un puerto
son de naturaleza diferente a las del otro puerto.
f1 f2
2.5.1. Transformadores
Sistemas el
ectricos
1 = L1 i1 + M12 i2 , (2.23)
2 = M21 i1 + L2 i2 (2.24)
i2
+
v2
n2
i1
+ n1
v1
Por otro lado, dividiendo (2.24) entre (2.23), con M12 = M21 = M , y usando
(2.26):
r
2 M i1 + L2 i2 L2
= =
1 L1 i1 + M i2 L1
A partir de esto se puede escribir:
i2 L1 L2 M
= q
i1 L2 L L1 M
2
r
L1 ( L1 L2 M )
L2 L1
= = (2.27)
L1 L2 M L2
y de (2.27):
n1
i2 = i1
n2
Lo cual significa que la potencia en ambos puertos es igual porque:
n2 n1
w2 = v2 i2 = v1 i1 = v1 i1 = w1
n1 n2
Sistemas mec
anicos traslacionales
x1 = a sin(), x2 = b sin()
F1 x1 = F2 x2
x1 = a sin(), x2 = b sin()
a
F2 = F1
b
Notese que la potencia en ambos puertos es igual:
b a
w2 = v2 F2 = v1 F1 = F1 v1 = w1
a b
F2
b
v2
v1
F1
Sistemas mec
anicos rotativos
en cada rueda. El hecho de que las ruedas sean dentadas evita que las ruedas
puedan patinar. Esto asegura que la longitud de arco de crculo s generado
por un angulo 1 descrito en el eje 1 sea igual a la longitud de arco generado
por un angulo 2 en el eje 2. Existe una relacion geometrica que relaciona al
radio de un crculo, el angulo descrito por el mismo y la longitud del arco de
crculo generado. Aplicando esta relacion en ambos ejes se tiene:
s = 1 r1 , s = 2 r2 (2.29)
1 r1 = 2 r2 (2.30)
T1 T2
1 s 2
!1 r1 r2 !2
n1
n2
F
Por otro lado, se sabe que el paso de los engranes (ancho de cada diente)
es igual en ambos ejes, pues esto es necesario para que dichos dientes se
adapten correctamente al girar. Esto significa que:
p1 = n1 , p2 = n2
p1 = 2r1 , p2 = 2r2
T1 = F r1 , T2 = F r2
2.5.2. Adaptadores
Sistemas mec
anicos
T1 = rF2 (2.34)
v2 = r1 (2.35)
Sistemas electromec
anicos
r
T1
F2 v2
Figura 2.18. Pi
non-cremallera.
N S
v+
R
L
i
+
E
(a) Vista superior.
F
B
N S
F
(b) Vista frontal.
F = il
uB (2.36)
F = ilB
= BA cos()
v = BA sin() = BA sin()
Notese que, de acuerdo a la Ley de Lenz [3], pag. 327, el voltaje inducido v
tiene una polaridad tal (vease la figura 2.19(a)) que se opone a la causa que
lo produce la cual en este caso es, a fin de cuentas, la corriente que circula i.
Ahora suponga que se colocan varias espiras con diferentes inclinaciones, de
manera que un conmutador mecanico automaticamente conecte y desconecte
las espiras de tal modo que siempre este conectada solamente la espira para
la cual = 90 . Entonces, el voltaje en las terminales del conjunto esta dado
como:
v = BA
40 2 Modelado matem
atico de sistemas fsicos
v = BA (2.37)
1
i= T
BA
Si se considera que en el puerto 1 se aplica un par T a una velocidad ,
produciendo en el puerto 2 un voltaje v y una corriente electrica i, entonces
se tiene un generador electrico de CD y ambas expresiones en (2.37) a un son
validas. Notese que en cualquiera de estos casos, la potencia en ambos puertos
es igual:
1
w1 = vi = T BA = T = w2
BA
En situaciones mas generales se acostumbra escribir las expresiones en (2.37)
del siguiente modo:
v = ke (2.38)
T = km i (2.39)
Sistemas electromec
anicos: fuerza debida a un campo magn
etico
i
+
F
v
dEM = F dx (2.40)
d dL(x)
v= =i
dt dt
Por tanto:
Por otro lado, como el cambio en la energa magnetica solo se debe al cambio
en la posicion x, entonces se puede escribir:
Em
dEm = dx
x
Entonces, de acuerdo a esto, (2.40), (2.45) y (2.41) se concluye que:
Em 1 L(x)
F = = i2 (2.46)
x 2 x
El desarrollo aqu presentado se basa en las ideas reportadas en [5], cap. 5.
2.6. Ejemplos
Una vez que se han establecido las funciones constitutivas de cada ele-
mento de sistema se debe abordar el problema de conectar varios de estos
elementos para obtener el modelo matematico de sistemas mas o menos com-
plejos. La clave para interconectar los elementos de un sistema mecanico es
considerar que todos ellos se conectan mediante uniones rgidas. Esto significa
que la velocidad a ambos lados de una union es igual. Para ser mas claros, a
continuacion se presentan varios ejemplos.
Ejemplo 2.2 Considere el sistema masa-resorte-amortiguador mostrado en
la figura 2.21(a) donde se aplica una fuerza externa F (t) sobre el cuerpo de
masa m. Este sistema quiz a sea el mas sencillo desde el punto de vista del
modelado, pero es muy importante desde el punto de vista de control porque
representa el comportamiento b asico de los sistemas de control de posici
on.
En la figura 2.21(b) se muestra el diagrama de cuerpo libre correspondien-
te. Notese que en cada union entre componentes existen dos fuerzas iguales
2.6 Ejemplos 43
K F (t) b
m
(a)
vK vm vb
F (t)
FK Fb
m
FK Fb
(b) Diagrama de cuerpo libre.
que en las figuras 2.4 y 2.9 tienen una velocidad y una fuerza aplicadas en el
mismo sentido. N otese tambien que se ha definido:
dxK
vK = (2.50)
dt
Considerando que los diferentes elementos de sistema se conectan mediante
uniones rgidas, de acuerdo a la figura 2.21(b) se puede escribir:
vK = vm = vb (2.51)
De acuerdo a esto se concluye que se puede escribir:
xK = xm + c1 = xb + c2
donde c1 y c2 son dos constantes. Si se usa la convenci
on de definir xK = 0,
xm = 0 y xb = 0 como aquellos puntos ocupados por el extremo m ovil del
resorte, el centro del cuerpo de masa m y el extremo movil del amortiguador,
respectivamente, cuando el sistema esta en reposo con F (t) = 0, entonces se
puede decir que c1 = c2 = 0 y, por tanto:
xK = xm = xb
Usando estas relaciones, tambien se puede escribir (2.47), (2.48), (2.49) co-
mo:
d2 xm dxm
m +b + Kxm = F (t) (2.52)
dt2 dt
Esta ultima expresi on es una ecuacion diferencial ordinaria, de segundo or-
den, lineal y de coeficientes constantes, que representa el modelo matem atico
del sistema masa-resorte-amortiguador. Esto significa que si se conocen los
ametros m, b y K, as como las condiciones iniciales xm (0), dxdtm (0) y la
par
fuerza externa aplicada F (t) como una funci on del tiempo, entonces se puede
resolver dicha ecuaci on diferencial para obtener la posici
on xm (t) y la velo-
cidad dxdtm (t) del cuerpo como funciones del tiempo. Dicho de otro modo, se
puede conocer la posici on y la velocidad del cuerpo en cualquier instante de
tiempo presente o en el futuro, es decir, para todo t 0. Finalmente, cuando
se estudia la solucion de las ecuaciones diferenciales es conveniente expresar
una ecuaci on diferencial de manera que el coeficiente de la potencia mayor
sea igual a la unidad, es decir, es mas conveniente escribir (2.52) como:
b K 1
x
m + x m + xm = F (t) (2.53)
m m m
d2 xm dxm
donde x
m = dt2 y x m = dt .
v = dx
dt , donde x es la posici
on del cuerpo.
v1 = dxdtK = dx
dt
b
, donde xK = xb representan la posici
on del extremo m
ovil
del resorte y del amortiguador, respectivamente.
De acuerdo a la figura 2.22(b) y a (2.4), (2.7), (2.18), se obtienen las siguien-
tes expresiones:
dv
masa: m= F (t) + FK + Fb (2.54)
dt
resorte: KxK = FK
amortiguador: bv1 = Fb
m F(t)
(a)
v1 v
FK
FK
Fb m F(t)
Fb
v1
(b) Diagrama de cuerpo libre.
v = v1 (2.55)
xK = xb = x (2.56)
46 2 Modelado matem
atico de sistemas fsicos
b K 1
x
+ x + x = F (t)
m m m
2
donde x = ddt2x y x = dx
dt . Este modelo matem
atico es identico al mostrado en
(2.53) si se considera que x = xm .
Ejemplo 2.4 En la figura 2.23(a) se muestran dos cuerpos unidos por un
resorte cuando se aplica una fuerza externa F (t) sobre el cuerpo de la izquier-
da. Esta situacion representa el caso practico en el que el cuerpo 1 transmite
movimiento al cuerpo 2 pero ambos cuerpos est an conectados mediante una
pieza mec anica que presenta flexibilidad. Esta flexibilidad no es introducida
de manera intencional sino que es un problema no deseado que aparece debido
a la rigidez finita del material con el que est a hecha la pieza mec anica que
une a los cuerpos 1 y 2. Con frecuencia es muy importante estudiar cual es
el efecto de esta flexibilidad en el desempe
no del sistema de control y por eso
debe ser tomada en consideraci on durante el modelado. Los amortiguadores
se introducen para considerar el efecto de la friccion de cada cuerpo con sus
soportes (o el suelo).
F (t)
K
m1 m2
b1 b2
(a)
v m1 vK v m2 v b2
F (t) F K1 K F b2
F b1
m1 m2
v b1 F b1 F K1 F K2 F K2 F b2
v1 v2
(b) Diagrama de cuerpo libre.
Las fuerzas que tienen el mismo sentido que vm1 y vm2 , y que act uan sobre
cada uno de los cuerpos, aparecen afectadas con un signo positivo y con un
signo negativo en caso contrario. N
otese tambien que se ha definido:
dxK
vK = = v1 v2
dt
Considerando que todos los elementos del sistema se conectan mediante unio-
nes rgidas entonces, de acuerdo a la figura 2.23(b) se tiene que:
dvm1
masa 1: m1 = F (t) + b1 vb1 KxK
dt
dvm2
masa 2: m2 = KxK b2 vb2
dt
48 2 Modelado matem
atico de sistemas fsicos
d2 xm1 dxm1
masa 1: m1 2
= F (t) b1 KxK
dt dt
d2 xm2 dxm2
masa 2: m2 = KxK b2
dt2 dt
Por otro lado, de acuerdo al p arrafo que sigue a (2.6) se tiene que xK =
x1 x2 , por lo que se puede escribir:
d2 xm1 dxm1
masa 1: m1 = F (t) b1 K(x1 x2 )
dt2 dt
d2 xm2 dxm2
masa 2: m2 = K(x1 x2 ) b2
dt2 dt
Finalmente, usando de nuevo (2.60) se obtiene que el modelo matem atico
correspondiente est
a dado por las dos siguientes ecuaciones diferenciales, las
cuales deben ser resueltas simult
aneamente:
d2 xm1 b1 dxm1 K 1
masa 1: 2
+ + (xm1 xm2 ) = F (t)
dt m1 dt m1 m1
d2 xm2 b2 dxm2 K
masa 2: 2
+ (xm1 xm2 ) = 0
dt m2 dt m2
Estas dos ecuaciones diferenciales no pueden ser combinadas para obtener
una sola ecuacion diferencial, ya que las variables xm1 y xm2 son linealmente
independientes, es decir, no se cuenta con ninguna expresi on algebraica que
relacione a estas variables.
Ejemplo 2.5 En la figura 2.24(a) se muestran dos cuerpos conectados a tres
resortes y un amortiguador. En este caso, el amortiguador representa la fric-
ci
on existente entre los dos cuerpos y, por esto, no puede ser sustituido por la
posible fricci
on existente entre cada cuerpo y el suelo (si no estuvieran colo-
cadas las ruedas debajo de cada cuerpo). El diagrama de cuerpo libre corres-
pondiente se muestra en la figura 2.24(b). La nomenclatura utilizada es la
siguiente:
x 1 = dxdt , donde x1 es la posici
1
on del cuerpo 1.
dx2
x 2 = dt , donde x2 es la posici on del cuerpo 2.
vb1 = dxdtb1 , donde xb1 es la posici on del extremo del amortiguador que
est a conectado al cuerpo 1.
vb2 = dxdtb2 , donde xb2 es la posici on del extremo del amortiguador que
est a conectado al cuerpo 2.
v1 = dx dt , donde x es la posici
on del extremo del resorte unido al cuerpo 1.
dy
v2 = dt , donde y es la posici on del extremo del resorte unido al cuerpo 2.
vK1 = dxdtK1 , donde xK1 es la posici on del extremo m ovil del resorte co-
nectado s olo al cuerpo 1.
2.6 Ejemplos 49
F(t) x1 K2 x2
K1 K3
m1 m2
b
(a)
F K2 F K2
v K1 x1 x2
F K2 F K3
F K1 F K2 F K3
m1 v1 v2 m2
Fb
F K1 F(t) Fb v K3
Fb Fb
v b1 v b2
(b) Diagrama de cuerpo libre.
masa 1: m1 x
1 = F (t) + FK1 FK2 Fb (2.61)
masa 2: m2 x
2 = FK2 FK3 + Fb (2.62)
resorte entre cuerpos: K2 xK2 = FK2
resorte izquierdo: K1 xK1 = FK1
resorte derecho: K3 xK3 = FK3
amortiguador: b(vb1 vb2 ) = Fb
Considerando que todos los elementos del sistema se conectan mediante unio-
nes rgidas entonces, de acuerdo a la figura 2.24(b) se tiene que:
masa 1: m1 x
1 = F (t) + K1 xK1 K2 xK2 b(vb1 vb2 )
masa 2: m2 x
2 = K2 xK2 K3 xK3 + b(vb1 vb2 )
masa 1: m1 x
1 = F (t) K1 x1 K2 xK2 b(x 1 x 2 )
masa 2: m2 x
2 = K2 xK2 K3 x2 + b(x 1 x 2 )
Por otro lado, de acuerdo al p arrafo que sigue a (2.6) se tiene que xK2 =
x1 x2 , por lo que se puede escribir:
masa 1: m1 x
1 = F (t) K1 x1 K2 (x1 x2 ) b(x 1 x 2 )
masa 2: m2 x
2 = K2 (x1 x2 ) K3 x2 + b(x 1 x 2 )
Finalmente, se obtiene que modelo matem atico correspondiente esta dado por
las dos siguientes ecuaciones diferenciales, las cuales deben ser resueltas si-
multaneamente:
b K1 K2 1
masa 1: x
1 + (x 1 x 2 ) + x1 + (x1 x2 ) = F (t)
m1 m1 m1 m1
b K3 K2
masa 2: x
2 (x 1 x 2 ) + x2 (x1 x2 ) = 0
m2 m2 m2
Tal como ocurre en el ejemplo anterior, estas dos ecuaciones diferenciales no
pueden ser combinadas para obtener una sola ecuaci on diferencial, ya que las
variables x1 y x2 son linealmente independientes, es decir, no se cuenta con
ninguna expresi on algebraica que relacione a estas variables. N otese tambien
que los terminos correspondientes al resorte central y al amortiguador aparecen
con signo contrario en cada una de estas ecuaciones diferenciales. Esto es
debido a que la fuerza que produce el resorte central o el amortiguador sobre
el cuerpo 1 es en sentido contrario a la fuerza que ejerce sobre el cuerpo 2.
Ejemplo 2.6 En la figura 2.25(a) se muestra un cuerpo rotativo unido a dos
muros a traves de un resorte y un amortiguador. Sobre este cuerpo se aplica
un par externo T (t). En este caso, el amortiguador representa la fricci on
existente entre el cuerpo y los rodamientos sobre los que se apoya. El resorte
representa la flexibilidad del eje o flecha del cuerpo giratorio. El diagrama de
cuerpo libre correspondiente se muestra en la figura 2.25(b). La nomenclatura
utilizada es la siguiente:
2.6 Ejemplos 51
= d
dt , donde es la posici on angular del cuerpo.
b = d dt
b
, donde b es la posici
on angular del extremo movil del amorti-
guador.
K = ddtK , donde K es la posici on angular del extremo m
ovil del resorte.
Usando la figura 2.25(b), as como (2.8), (2.11), (2.21), se obtienen las si-
guientes expresiones:
b K
I
T(t)
(a)
Tb Tb TK TK K
I
b
!b T(t)
(b) Diagrama de cuerpo libre.
Notese que los pares que se aplican sobre el cuerpo I y tienen el mismo
aparecen afectados por un signo positivo en
sentido que la velocidad angular ,
(2.65). En caso contrario aparecen afectados por un signo negativo, como es
el caso de TK . Considerando que todos los elementos del sistema se conectan
mediante uniones rgidas entonces, de acuerdo a la figura 2.25(b) se tiene que:
= b = K (2.66)
= b = K (2.67)
I = T (t) + bb KK
= T (t) b K
Motor
T(t) I1 I2
Cremallera m
(a)
T2
!1
r r
T1
F1 F2 F2
F1 w2
m
v1 vm v2
Fb
vb
Fb
(b) Diagrama de cuerpo libre.
b1 T b1 T1 !1
I1
! b1 T b1 T(t)
T2 !2 T b2
b2
I2
T b2 ! b2
(c) Diagrama de cuerpo libre (cont.).
Usando las figuras 2.26(b) y 2.26(c), as como (2.4), (2.8), (2.18), (2.21),
(2.34), (2.35) se obtienen las siguientes expresiones:
d1
inercia 1: I1 = T (t) + T1 Tb1 (2.69)
dt
d2
inercia 2: I2 = T2 Tb2 (2.70)
dt
54 2 Modelado matem
atico de sistemas fsicos
dvm
masa: m = F1 F2 + Fb (2.71)
dt
amortiguador izquierdo: b1 b1 = Tb1
amortiguador derecho: b2 b2 = Tb2
amortiguador central: bvb = Fb
pi
non-cremallera izquierda: T1 = rF1 , v1 = r1 (2.72)
pi
non-cremallera derecha: T2 = rF2 , v2 = r2 (2.73)
Notese que las fuerzas que tienen el mismo sentido que 1 , 2 y vm aparecen
afectadas con un signo positivo en (2.69), (2.70), (2.71) y con un signo nega-
tivo en caso contrario. Considerando que todos los elementos del sistema se
conectan mediante uniones rgidas entonces, de acuerdo a las figuras 2.26(b)
y 2.26(c) se tiene que:
vm = v1 = v2 = vb (2.74)
1 = b1 , 2 = b2 (2.75)
xm = xb (2.76)
1 = b1 , 2 = b2
Por tanto, las expresiones en (2.69), (2.70) y (2.71) se pueden escribir como:
d1
inercia 1: I1 = T (t) + rF1 b1 b1 (2.77)
dt
d2
inercia 2: I2 = rF2 b2 b2 (2.78)
dt
dvm
masa: m = F1 F2 + bvb (2.79)
dt
Recordando que v1 = v2 , se pueden usar (2.72) y (2.73) para encontrar
que 1 = 2 . Entonces se puede usar este resultado as como (2.75) para,
restando (2.78) a (2.77), encontrar:
d1
(I1 + I2 ) = T (t) + r(F1 F2 ) (b1 + b2 )1 (2.80)
dt
Por otro lado, usando vm = v1 = r1 = vb se puede escribir (2.79) como
d1
mr = F1 F2 bvm
dt
2.6 Ejemplos 55
b1 K b2
I1 I2
T(t)
(a)
!1 T1 !2 T3 T3 !4 T4 T6 T6
I1 I2
T1 T(t) !3 T4 !5 !6
(b) Diagrama de cuerpo libre.
2 = 1 = 3 , 5 = 6 = 4 (2.83)
2 = 1 = 3 , 5 = 6 = 4 (2.84)
Estas dos ecuaciones diferenciales no pueden ser agrupadas en una sola porque
las variables 2 y 5 son linealmente independientes, es decir, no se cuenta
con una expresion algebraica que relacione a estas dos variables.
2.6 Ejemplos 57
b 1 T(t) n1
I1 b2
I2
n2
(a)
T b1 T b1 ! T n1 ! n1
I1 T n2 T3 !4 T b2 T b2
! b1 T(t) I2
! n2 ! b2
(b) Diagrama de cuerpo libre.
= b1 = n1 , 4 = b2 = n2 (2.88)
= b1 = n1 , 4 = b2 = n2 (2.89)
Por tanto, usando (2.88) y (2.89) las expresiones en (2.85) y (2.86) se pueden
escribir como:
n1
cuerpo 1: I1 = T (t) b1 T3 (2.90)
n2
cuerpo 2: I2 4 = T3 b2 4 (2.91)
Por otro lado, de (2.87) y (2.88) se tiene = nn21 4 . Usando esto, despejan-
do T3 de (2.91), sustituyendo en (2.90) y despues de acomodar terminos se
encuentra:
2 ! 2 !
n2 n2 n2
I2 + I1 4 + b2 + b1 4 = T (t) (2.92)
n1 n1 n1
Notese que las dos ecuaciones diferenciales en (2.90) y (2.91) se han podi-
do combinar en una sola ecuaci on gracias a que las variables y 4 son
linealmente dependientes pues est an relacionadas a traves de = nn12 4 . Esta
ecuacion diferencial representa el modelo matematico buscado: si se conoce el
par externo aplicado T (t) entonces se puede conocer la posicion 4 del cuerpo
2 al resolver la ecuaci
on diferencial anterior.
Ejemplo 2.10 En la figura 2.29(a) se muestra el mismo mecanismo de la
figura 2.26(a) (estudiado en el ejemplo 2.7) pero ahora la cremallera es fle-
xible. Esta flexibilidad no se introduce de manera intencional, sino que es un
2.6 Ejemplos 59
problema no deseado que aparece debido a la rigidez finita del material con
el que est
a hecha la cremallera. Para encontrar el modelo matem atico en este
caso se puede aprovechar el desarrollo realizado en el ejemplo 2.7. Es decir, se
puede partir de las ecuaciones (2.77), (2.78) y (2.79), las cuales se reescriben
a continuacion para facilitar la referencia:
d1
inercia 1: I1 = T (t) + rF1 b1 b1
dt
d2
inercia 2: I2 = rF2 b2 b2
dt
dvm
masa: m = F1 F2 + bvb
dt
considerando ahora que, de acuerdo a la figura 2.29(b) y las convenciones en
la figura 2.4 y (2.7):
F1 = K1 (xA xB ), F2 = K2 (xC xD )
con:
vA = dx dxB
dt y vB = dt , donde xA y xB son las posiciones de los extremos
A
xm
I1 K1 K2 I2
m
A B C D
(a)
v AB
F1 F1
vA K1 vB
v CD
F2 F2
vC K2 vD
(b) Flexibilidad en la cremallera.
d1
inercia 1: I1 = T (t) + rK1 (xA xB ) b1 1 (2.93)
dt
d2
inercia 2: I2 = rK2 (xC xD ) b2 2 (2.94)
dt
masa: m
xm = K1 (xA xB ) K2 (xC xD ) bx m (2.95)
Por otro lado, de acuerdo a (2.72), (2.73), (2.74), (2.76) y considerando que
todos los elementos del sistema se conectan mediante uniones rgidas se en-
cuentra que:
xB = xC = xm , xA = r1 , xD = r2 (2.96)
donde 1 y 2 son las posiciones angulares de las ruedas dentadas 1 y 2 defi-
nidas de acuerdo a la figura 2.26(b). Sustituyendo (2.96) en (2.93), (2.94) y
(2.95):
inercia 1: I1 1 + b1 1 rK1 (r1 xm ) = T (t)
inercia 2: I2 2 + b2 2 rK2 (xm r2 ) = 0
masa: m
xm + bx m K1 (r1 xm ) + K2 (xm r2 ) = 0
El modelo matem atico est
a dado por estas tres ecuaciones diferenciales que
deben ser resueltas simult aneamente. En este caso no se pueden combinar
estas tres ecuaciones para obtener una sola ecuaci on diferencial equivalente
debido a que las variables 1 , 2 y xm son linealmente independientes pues no
se cuenta con una expresi on algebraica que relacione a estas tres variables.
Notese que desde el punto de vista del modelado este es el efecto que introdu-
cen los dos resortes considerados en este ejemplo. Comp arese con el modelo
matem atico obtenido en el ejemplo 2.7.
b1 n1
I1 K b2
I2
T(t)
n2
(a)
T n2 !B
K
!A T n2
(b) Flexibilidad en
los engranes.
izquierdo recibe un par externo T (t). Por tanto, este ejemplo representa el
caso de un motor (cuerpo izquierdo) que transmite movimiento a una carga
(cuerpo derecho) a traves de una caja de engranes cuyos dientes son flexibles.
Esta flexibilidad es un problema no deseado que, sin embargo, aparece en la
pr
actica debido a la rigidez finita del material con el que se hacen los engranes
(acero). Para resolver este problema se puede aprovechar el desarrollo presen-
tado en el ejemplo 2.9 para modelar el sistema mostrado en la figura 2.28(a).
Por tanto, se puede regresar a las expresiones en (2.90) y (2.91), las cuales
se reescriben a continuacion para facilitar la referencia:
n1
cuerpo 1: I1 = T (t) b1 Tn2 (2.97)
n2
cuerpo 2: I2 4 = Tn2 b2 4 (2.98)
Tn2 = K(A B )
Estas dos ecuaciones diferenciales deben ser resueltas simult aneamente y re-
presentan el modelo matem atico buscado. Notese que estas dos ecuaciones
diferenciales no pueden ser combinadas en una sola debido a que las variables
y 4 son linealmente independientes, es decir no se cuenta con una expresion
algebraica que relacione a estas variables. Esta es la principal diferencia con
el problema resuelto en el ejemplo 2.9 y es debida a la introducci
on del resorte
en el presente ejemplo.
n1
+ Pion
u Motor
b1
r
n2
Cremallera
b2
R L
+ T n1
+ b1
u i eb
Im r
IP
n2 b2
M
b3
(a)
b1 T b1 T T1 !1 n 1
Im T2 !2
T b1 ! b1 m
n2
(b)
!P !P
T b2 T b2
r
IP
! b2 TP
T2 v
F
M x
F b3
v b3
F b3
(c)
Notese que, en las expresiones (2.99), (2.100) y (2.101), las fuerzas que tienen
el mismo sentido que m , P y x aparecen afectadas con un signo positivo y
con un signo negativo en caso contrario. Considerando que todos los elementos
del sistema se conectan mediante uniones rgidas entonces, de acuerdo a la
figura 2.32 se tiene que:
v = x = vb3 (2.104)
m = b1 = 1 , 2 = b2 = P (2.105)
x = xb3 (2.106)
m = b1 (2.107)
Por otro lado, usando la Ley de voltajes de Kirchhoff (vease el ejemplo 2.14)
al circuito de la figura 2.32(a) se obtiene:
di
L + Ri + eb = u (2.113)
dt
Usando (2.38), (2.39) y m = nn12r x se encuentra que el voltaje inducido y el
par generado est
an dados como:
n2
eb = ke m = ke x,
T = km i
n1 r
Entonces, usando este resultado, (2.113) y (2.112) se encuentra, finalmente,
que el modelo matem atico buscado esta representado por las dos ecuaciones
diferenciales siguientes, las cuales deben ser resueltas simult
aneamente:
di n2
+ Ri + ke L x = u (2.114)
dt n1 r
2 ! 2 !
n2 1 n2 1 n2
Im + IP 2 + M x + b1 + b2 2 + b3 x = km i
n1 r r n1 r r n1 r
(2.115)
Otra manera de expresar este modelo matem atico es mediante una sola ecua-
ci
on diferencial obtenida al combinar (2.114) y (2.115). Para esto, dervese
66 2 Modelado matem
atico de sistemas fsicos
di
una vez respecto al tiempo a (2.115). Esto produce la aparici on de dt en el
miembro derecho de la ecuaci on resultante. A continuaci on se sustituye el va-
di
lor de dt obtenido al despejar de (2.114). Con esto se obtiene una ecuaci on
d3 x
diferencial en terminos de de dt3 , x, x e i. Para eliminar i, se despeja esta
variable de (2.115) y se sustituye. As, se obtiene una ecuacion diferencial de
tercer orden en la variable x con el voltaje u como excitaci on. De este modo
el modelo matem atico queda representado por una sola ecuaci on diferencial a
partir de la cual se puede conocer la posici on x de la masa M si se conoce el
voltaje u aplicado al motor. Se deja como ejercicio para el lector el realizar el
procedimiento descrito.
I = ml2
Notese que los pares que se aplican sobre el cuerpo I y tienen el mismo sentido
aparecen afectados por un signo positivo en (2.117).
que la velocidad angular ,
En caso contrario aparecen afectados por un signo negativo, como es el caso
de Tg . Considerando que todos los elementos del sistema se conectan mediante
uniones rgidas entonces, de acuerdo a la figura 2.33(b) se tiene que:
= b (2.118)
2.6 Ejemplos 67
T(t)
l
Tg
g Tb Tb
m
I
b
d !b T(t)
(a) (b) Diagrama de cuerpo libre.
= b (2.119)
Notese que esta es una ecuaci on diferencial de segundo orden pero no lineal,
caracterstica que es debida al hecho de que aparece la funci on sin(). En
la secci
on 7.2 del captulo 7 se explica la manera de manipular este tipo de
ecuaciones diferenciales. M as aun, en los captulos 11, 13 y 14 se muestra
como dise nar controladores para este tipo de ecuaciones diferenciales.
Usando las convenciones de signo de las figuras 2.7, 2.8, 2.12 se obtiene
el diagrama mostrado en la figura 2.34(b). Entonces, a partir de esta figura y
68 2 Modelado matem
atico de sistemas fsicos
C L
+
vi R
(a)
+ vC + vL
C L
+ +
vi i R vR
(b)
donde L{vC (t)} = VC (s), L{vL (t)} = VL (s) y L{i(t)} = I(s). Por tanto, aho-
ra se tiene una relacion algebraica entre las transformadas de Laplace de los
voltajes y las corrientes en un inductor y un capacitor. Al factor que relaciona
a estas variables se le denomina impedancia y se representa como:
V (s)
Z(s) = (2.126)
I(s)
VL (s)
ZL (s) = = sL (2.127)
I(s)
VC (s) 1
ZC (s) = = (2.128)
I(s) sC
ii L C R
(a)
+
ii L C R v
iL iC iR
(b)
Usando las convenciones de signo de las figuras 2.7, 2.8, 2.12 se obtiene
el diagrama mostrado en la figura 2.35(b). Entonces, a partir de esta figura y
2.6 Ejemplos 71
Z1 (s)Z2 (s)
ZT (s) = (2.140)
Z1 (s) + Z2 (s)
2.6 Ejemplos 73
I(s) R C
+ +
V i (s) C R V o (s)
C C C
+ +
V 1(s) I1 R I2 R I3 R V 2(s)
1
I1 (s) + (I1 (s) I2 (s))R = V1 (s) (2.144)
sC
1
(I2 (s) I1 (s))R + I2 (s) + (I2 (s) I3 (s))R = 0
sC
1
(I3 (s) I2 (s))R + I3 (s) + V2 (s) = 0, V2 (s) = RI3 (s)
sC
De la tercera ecuaci
on en (2.144) se obtiene:
1 1
I2 (s)R = R+ V2 (s) + V2 (s) (2.145)
R sC
ecuaci
on mostrada en (2.144) se obtiene V1 (s) como funci
on de V2 (s) so-
lamente. Finalmente, acomodando terminos convenientemente se obtiene la
expresi
on buscada:
V2 (s) R3 C 3 s3
= 3 3 3 = F (s) (2.146)
V1 (s) R C s + 6R2 C 2 s2 + 5RCs + 1
L D
i
E Q v R
E v C R
u
(b) Convertidor con interruptor ideal.
L
i a
E I v C R
E I v C R
di
L = E (1 uav ) (2.155)
dt
d
C = (1 uav )i (2.156)
dt R
donde i y representan, respectivamente, la corriente y el voltaje promedio en
el inductor y el capacitor. Mientras que la entrada uav ahora denota la entrada
de control promedio, la cual puede tomar valores en el intervalo continuo [0,1),
siendo esta la funcion de posicion promedio del interruptor.
Como consecuencia de imponer que la entrada del convertidor sea uav =
di
U = constante, se encuentra que en estado estacionario (cuando dt = 0 y
d
dt = 0) el sistema (2.155)-(2.156) genera que:
0 (1 U ) i E
= (2.157)
(1 U ) R1 0
1 E
i= (2.158)
R (1 U )2
E
= (2.159)
(1 U )
De esta manera, a partir de las relaciones en estado estacionario, determi-
nadas por (2.158)-(2.159), es claro que se tiene la siguiente relaci
on est
atica
(en estado estacionario) promedio:
1
H(U ) = = (2.160)
E (1 U )
De (2.160) se puede observar que el cociente E es mayor o igual a uno,
pues U [0, 1), de ah que a este convertidor se le conozca como elevador de
voltaje. Finalmente, en la figura 2.40 se puede observar la curva caracterstica
est
atica del convertidor Boost, donde claramente se ve cumplida la desigualdad
E.
10
6
H(U)
0
0 0.2 0.4 0.6 0.8 1
U
por:
d1 1
C1 = i2 , (2.163)
dt R1
d2 2
C2 = . (2.164)
dt RL
b) Para el caso en el que Q1 est a encendido y Q2 apagado, i. e., u1 = 1 y
u2 = 0, a partir de la figura 2.41, el circuito equivalente es el que muestra en
la figura 2.43(b). Aplicando la LVK a las mallas I y II de la figura 2.43(b)
se obtienen las ecuaciones para i1 e i2 respectivamente:
di1
L1 = E, (2.165)
dt
di2
L2 = 1 2 . (2.166)
dt
De igual forma, aplicando la LCK en los nodos A y B de la figura 2.43(b), se
obtienen las ecuaciones que gobiernan a los voltajes 1 y 2 :
d1 1
C1 = i2 , (2.167)
dt R1
d2 2
C2 = i2 (t) . (2.168)
dt RL
c) Cuando Q1 est a apagado y Q2 encendido, i. e., u1 = 0 y u2 = 1, el circuito
resultante es el que muestra en la figura 2.43(c) y aplicando la LVK en las
mallas I y II, se logra el planteamiento de las ecuaciones que describen el
comportamiento de las corrientes i1 e i2 respectivamente:
di1
L1 = 1 + E, (2.169)
dt
di2
L2 = 1 . (2.170)
dt
Asimismo, valiendose de la LCK para los nodos A y B de la figura 2.43(c),
se tiene lo siguiente:
d1 1
C1 = i1 i2 , (2.171)
dt R1
d2 2
C2 = . (2.172)
dt RL
d) Finalmente, cuando u1 = 0 y u2 = 0, el cual corresponde al caso en el que
los transistores se encuentran apagados, el circuito equivalente es el que se
muestra en la figura 2.43(d). Para obtener las ecuaciones que gobiernan a los
voltajes y corrientes ilustrados en la figura 2.43(d), de igual forma, se aplica
la LVK para las mallas I y II obteniendo lo siguiente:
2.6 Ejemplos 81
di1
L1 = 1 + E, (2.173)
dt
di2
L2 = 1 2 . (2.174)
dt
Mientras que, aplicando la LCK para los nodos A y B, se tiene:
d1 1
C1 = i1 i2 , (2.175)
dt R1
d2 2
C2 = i2 . (2.176)
dt RL
As, mediante inspecci
on, al comparar las ecuaciones (2.161)-(2.162),
(2.165)-(2.166), (2.169)-(2.170) y (2.173)-(2.174), se producen como resul-
tado las ecuaciones que gobiernan el comportamiento de i1 e i2 , determinadas
respectivamente por (2.177) y (2.179). Realizando un proceso similar para las
variables 1 y 2 , determinadas por las ecuaciones (2.163)-(2.164), (2.167)-
(2.168), (2.171)-(2.172) y (2.175)-(2.176), se obtienen las ecuaciones (2.178)
y (2.180) respectivamente. De esta manera, el modelo matem atico del converti-
dor Boost-Boost queda determinado por el sistema de ecuaciones diferenciales
de cuarto orden no lineal, siguiente:
di1
L1 = (1 u1 )1 + E, (2.177)
dt
d1 1
C1 = (1 u1 )i1 i2 , (2.178)
dt R1
di2
L2 = 1 (1 u2 )2 , (2.179)
dt
d2 2
C2 = (1 u2 )i2 . (2.180)
dt RL
As, las ecuaciones (2.177)-(2.180) constituyen el modelo matem atico del con-
vertidor de potencia Boost-Boost. Si se desea corroborar que estas ecuaciones
reproducen los cuatro modos de operaci on del convertidor (o circuitos equiva-
lentes de la figura 2.43), basta con asignarle los valores correspondientes a las
entradas u1 y u2 .
Con la finalidad de conocer las propiedades del convertidor de potencia
Boost-Boost en estado estacionario se parte de su modelo promedio asociado,
el cual se obtiene suponiendo que ambos transistores conmutan a alta frecuen-
cia y esta determinado como:
di1
L1 = (1 u1av ) 1 + E, (2.181)
dt
d
1 1
C1 = (1 u1av )i1 i2 , (2.182)
dt R1
di2
L2 = 1 (1 u2av )
2 , (2.183)
dt
82 2 Modelado matem
atico de sistemas fsicos
d
2 2
C2 = (1 u2av )i2 , (2.184)
dt RL
donde i1 , 1 , i2 y 2 representan las corrientes y los voltajes promedio asocia-
dos a L1 , C1 , L2 y C2 respectivamente. Las entradas u1av y u2av representan
las entradas de control promedio, las cuales pueden tomar valores en el in-
tervalo continuo [0, 1), siendo estas las funciones de posici on promedio de los
interruptores.
De esta manera, partiendo del modelo promedio (2.181)-(2.184), se procede
a obtener los valores en estado estacionario del sistema que aparecen cuando
las entradas del convertidor u 1av y u
2av son constantes en alg un valor dentro
di2
del intervalo [0, 1). Esto se consigue suponiendo que ddti1 = 0, d 1
dt = 0, dt = 0
2.7 Caso de estudio 83
d
2
y dt= 0, para encontrar:
0 (1 u1av ) 0 0 1 E
(1 u 1 1 0
1av ) R1 1 0 = ,
0 1 0 (1 u2av ) 2 0
0 0 (1 u
2av ) R1L 2 0
De esta manera, (2.185) permite conocer de manera directa los valores que
toman las variables promedio i1 , 1 , i2 y 2 , en estado estacionario. Asi-
mismo, a partir de la u ltima relaci
on en (2.185) se puede observar que la
salida satisface la desigualdad 2 E, pues recordemos que u1av [0, 1) y
u1av [0, 1).
i C
v -
+
L
+
E C0 R v0
(a) Serie
i i0
L0
+
L
+
v C0 R
E v0
- -
(b) Paralelo
L C
Q1 Q4 i +v
+ D1 D4 +
E
R v0
Q3
C0
Q2
D2 D3
L C
v-
+
i
+
E(t)
- Co
v0
una fuente de voltaje alterno y cuadrado E(t) que toma valores E. Notese
que, debido a la presencia de los diodos rectificadores, el funcionamiento del
circuito puede ser separado en dos casos: 1) cuando i > 0, figura 2.47(a), y
2) cuando i < 0, figura 2.47(b). A partir de la figura 2.47(a) se obtienen las
ecuaciones para i > 0. Aplicando la Ley de Kirchhoff de voltajes al trayecto
cerrado de la izquierda se obtiene:
di
E(t) = L
+ v + v0 , i > 0
dt
Es importante resaltar que se ha tomado en consideracion el sentido indicado
para la corriente i, as como la polaridad indicada para v0 . Por otro lado, de
acuerdo a (2.120) la relacion entre el voltaje y la corriente en el capacitor
resonante esta dada como:
dv
C = i, i > 0
dt
Finalmente, aplicando la Ley de Kirchhoff de corrientes al nodo com un a
ambos capacitores se obtiene:
86 2 Modelado matem
atico de sistemas fsicos
dv0 v0
i = C0 + , i>0
dt R
L C
+ -
+
i v
v0
+
E(t) - C0 R
-
(a) i > 0
L C
i +
v- -
+
E(t) - C0 R v0
+
(b) i < 0
Por otro lado, a partir de la figura 2.47(b) se obtienen las ecuaciones para
i < 0. Aplicando la Ley de Kirchhoff de voltajes al trayecto cerrado de la
izquierda se obtiene:
di
E(t) = L + v v0 , i < 0
dt
Notese que v0 esta afectado por un signo porque la polaridad indica-
da para este voltaje esta en sentido opuesto al definido por el sentido de la
corriente electrica i. Sin embargo, la relacion entre el voltaje y la corriente en
el capacitor resonante no sufre cambios:
dv
= i, i < 0
C
dt
Finalmente, aplicando la Ley de Kirchhoff de corrientes al nodo com
un a
ambos capacitores se obtiene:
v0 dv0
i= C0 , i<0
R dt
Notese que el sentido correcto de la corriente en el capacitor C0 y en la resis-
tencia es de la terminal marcada con + a la terminal marcada con , es
decir, es contrario al sentido indicado para i.
2.8 Resumen del captulo 87
R R
+
v1 v v2
+
C
+
+ vo
vi R
Q
L
i
E D v C R
i L
u
E v C R
u
(b) Representaci
on ideal del convertidor.
d
C = i (2.190)
dt R
Las variables i, representan la corriente presente en el inductor L y el
voltaje entre las terminales del capacitor C, respectivamente. Mientras
que E representa el voltaje de entrada al sistema y R es la resistencia de
salida del sistema. Finalmente, la variable u denota la funcion de posicion
del interruptor o conmutador, la cual act ua como variable de control y
solo toma el valor de 0 o de 1.
9. El sistema Convertidor de potencia Buck-Motor de CD, mostrado en la
figura 2.51, representa una forma de conveniente de alimentar a un motor
de CD. El objetivo de control de este sistema consiste en lograr que la
velocidad angular de la flecha del motor , se estabilice a un valor de
velocidad deseada a traves del voltaje , obtenido del convertidor de
potencia Buck, mediante el dise no adecuado de u.
Partiendo del hecho de que el modelo matematico de un motor de CD
esta determinado por:
dia
La = Ra ia ke ,
dt
d
J = km ia B.
dt
Demuestre que el modelo matematico del sistema Convertidor de potencia
Buck-Motor de CD, representado idealmente en la figura 2.51(b), esta de-
terminado por el siguiente sistema de ecuaciones diferenciales lineales:
di
L = + Eu, (2.191)
dt
d
C = i ia ,
dt R
92 2 Modelado matem
atico de sistemas fsicos
dia
La = Ra ia ke ,
dt
d
J = km ia B.
dt
donde:
i es la corriente presente en el inductor del convertidor Buck.
es el voltaje de salida del convertidor y a su vez es el voltaje aplicado
en las terminales de armadura del motor.
u representa la funcion de posicion del interruptor o conmutador, la
cual act
ua como variable de control y solo toma el valor de 0 o de 1.
L denota la inductancia del convertidor Buck.
C representa la capacitancia del convertidor Buck.
E es el voltaje de alimentacion del sistema convertidor Buck-Motor de
CD.
R denota la resistencia de salida del convertidor.
ia es la corriente electrica de armadura.
es la velocidad angular del rotor del motor.
La representa la inductancia de armadura.
Ra es la resistencia de armadura.
D !
(a) Construcci
on del sistema convertidor Buck-Motor de CD
mediante transistor-diodo.
E(t)
50
[V] 0
-50
0 0.5 1 1.5 2
t [s] x 10
-4
30
P(t)
20
[Watt]
10
0
0 2 4 6 8
t [s] x 10
-4
1.5
i(t)
0.5
[A]
-0.5
v(t)
-1.5
-500 -250 0 250 500
[V]
Otro resultado importante es el teorema del valor final [1], pag. 25, [3], pag.
304:
3.1. Ecuaci
on diferencial de primer orden
Considere la siguiente ecuacion diferencial ordinaria, lineal, de coeficientes
constantes de primer orden:
y + a y = k u, y(0) = y0 (3.4)
k 1
Y (s) = U (s) + y0 (3.6)
s+a s+a
Para continuar es necesario especificar la funcion del tiempo u. Suponga que:
A
u = A, U (s) = (3.7)
s
100 3 Base matem
atica: ecuaciones diferenciales
donde A es una constante. Por tanto, ahora se puede escribir (3.6) como:
kA 1
Y (s) = + y0 (3.8)
s(s + a) s + a
De acuerdo al metodo de expansion en fracciones parciales [2], Cap. 4, [3],
pag. 319, si a 6= 0 entonces se debe escribir:
kA B C
= + (3.9)
s(s + a) s+a s
donde B y C son dos constantes. El valor de B se obtiene multiplicando ambos
lados de (3.9) por el factor (s + a) y evaluando las expresiones resultantes en
s = a:
kA B C
(s + a) = (s + a) + (s + a)
s(s + a) s=a s+a s=a s s=a
es decir:
kA kA
B= = (3.10)
s s=a a
es decir:
kA kA
C= = (3.11)
s + a s=0 a
y(t)
= kAeat ay0 eat
y se sustituye y(t)
e y(t) en (3.4) para obtener:
kA at kA
+ ay(t) = kAeat ay0 eat + a
y(t) e + + y0 eat = kA = ku
a a
porque en (3.7) se ha definido u = A. Esto significa que la expresi on para
y(t) presentada en (3.14) satisface a (3.4) y, por tanto, se ha comprobado que
(3.14) es la soluci
on de (3.4). Este procedimiento debe seguirse para comprobar
la soluci
on de cualquier ecuaci on diferencial.
La solucion dada en (3.14) se puede descomponer en dos partes:
Estudio gr
afico de la soluci
on
kA
y( ) = 0.632 (3.20)
a
Finalmente, es importante se nalar que y(t) crece sin lmite si a < 0. Si a = 0
entonces se tiene el caso estudiado en la seccion 3.2, es decir y(t) = kAt crece
sin lmite al aumentar el tiempo.
y(t)
kA
a
0:632 kA
a
tiempo
Funci
on de transferencia
o tambien:
Y (s) k
= G(s) = (3.21)
U (s) s+a
Ejemplo 3.2 Suponga que se tiene un tanque que almacena un lquido incom-
presible como el agua (vease la figura 3.3). El tanque es de paredes paralelas,
es decir de secci
on constante cuya a rea es representada por C. El lquido es
3.1 Ecuaci
on diferencial de primer orden 105
Im (s)
a 0 a Re (s)
a>0 a<0
qi
h
R qo
dh h
C = qi
dt R
Entonces, el modelo matem
atico del sistema de nivel de la figura 3.3 queda
expresado como:
dh 1
C + h = qi
dt R
Dividiendo entre C, se puede expresar este modelo como la siguiente ecuaci on
diferencial ordinaria, lineal, de coeficientes constantes, de primer orden:
dh
+ ah = kqi
dt
1 1
a= , k= (3.26)
RC C
La razon por la cual esta ecuaci on diferencial es el modelo de este sistema
de nivel de lquido es porque mediante su soluci on se puede conocer como
evoluciona el nivel h(t) en el tiempo, si se conocen el flujo del lquido que se
introduce qi (t), el area de la secci
on del tanque C, si se tiene informacion a
cerca de que tan grande es la abertura de la v avula de salida, es decir si se
conoce R, y si se conoce el nivel inicial.
3.1 Ecuaci
on diferencial de primer orden 107
1
tanque decrece (n
otese que a = RC > 0) hasta que el tanque se vaca
lmt h(t) = 0. Esto ocurre m as r
apido si el producto RC es menor
(cuando a es mayor) y m as lentamente cuando RC es mayor (cuando a
es menor).
h0
h(t)
R 1C 1 R 2C 2
tiempo
Figura 3.4. Comportamiento del nivel cuando h0 > 0 y qi (t) = A = 0 (R1 C1 <
R2 C2 ).
3.2. Un integrador
Considere la siguiente ecuacion diferencial:
y = ku (3.29)
donde k es una constante real diferente de cero. Notese que esta ecuacion se
obtiene a partir de (3.4) cuando a = 0. Integrando directamente:
Z y(t) Z t
dy = ku dt
y(0) 0
Z t
y(t) = k u(t) dt + y(0) (3.30)
0
3.2 Un integrador 109
En este caso yn (t) permanece constante cuando el tiempo crece. Notese que si
u = A es constante, entonces la respuesta forzada crece sin lmite yf (t) = kA t.
En la figura 3.5 se muestra graficada la solucion y(t) = y(0) + kA t, para los
casos a) y(0) 6= 0 y u = A =constante positiva, y b) y(0) 6= 0 y u = A = 0. A
los sistemas fsicos representados por este tipo de ecuacion diferencial se les
llama integradores porque, si y(0) = 0 la solucion y(t) (salida) es la integral
de la entrada u.
Ejemplo 3.4 Considere el tanque con agua mostrado en la figura 3.3. Su-
ponga que la v
alvula de salida est
a totalmente cerrada por lo que R .
Entonces, el modelo en (3.26) se reduce a:
dh
= kqi
dt
1
porque a = RC 0 cuando R . Entonces, la evoluci
on del nivel en el
tiempo esta descrita por (3.30), es decir:
Z t
h(t) = h(0) + k qi (t) dt
0
Se puede hacer uso de la figura 3.5, haciendo h(t) = y(t) y qi (t) = u, para
apreciar gr aficamente la soluci
on en (3.31) en los casos i) h(0) > 0 y qi =
A =constante positiva, e ii) h(0) > 0 y qi = A = 0. N otese que el nivel
permanece constante si no se introduce lquido (cuando A = 0) y que el nivel
crece sin lmite cuando A > 0 es constante. El lector puede usar su experiencia
previa para comprobar que estas dos situaciones verdaderamente suceden en
un sistema de nivel de lquido real.
Ejemplo 3.5 La ecuaci on diferencial en (3.29) tambien puede resolverse
usando expansion en fracciones parciales, es decir, usando el metodo de la
secci
on 3.1. Usando (3.1) en (3.29) se obtiene:
sY (s) y0 = k U (s)
k 1
Y (s) = U (s) + y0 (3.32)
s s
110 3 Base matem
atica: ecuaciones diferenciales
y(t)
kA
y0
0
0
tiempo
(a) y0 6= 0 y u = A =constante positiva
y(t)
y0
tiempo
(b) y0 6= 0 y u = A = 0
kA B C
2
= + 2 (3.34)
s s s
donde B y C son dos contantes. El valor de C se obtiene multiplicando ambos
lados de (3.34) por s2 y evaluando en s = 0, es decir:
2 kA B 2 C 2
s 2 = s + 2s
s s=0 s s=0 s s=0
donde yn (t) y yf (t) reciben los nombres de respuesta natural y respuesta for-
zada, respectivamente.
Tal como se explic o en la seccion 3.1, la respuesta natural yn (t) recibe
dicho nombre porque est a determinada u nicamente por la estructura (o natu-
raleza) de la ecuacion diferencial, es decir, por la transformada de Laplace de
y presente en el lado izquierdo de (3.29). Esto da origen al termino 1s y0 en
(3.35), lo cual permite tener yn (t) = y0 en (3.38).
La respuesta forzada yf (t) es debida a la excitaci on (entrada) u = A.
Sin embargo, en este caso la respuesta forzada tambien recibe el efecto de la
transformada de Laplace del termino y presente en el lado izquierdo de (3.29),
pues ambos generan el termino kA s2 , del cual se obtiene yf (t) = kAt.
112 3 Base matem
atica: ecuaciones diferenciales
s
olo tiene un polo en s = 0. N otese que la entrada U (s) = As tambien tiene
un polo en s = 0. La combinaci on de estos dos polos repetidos en (3.33)
produce una respuesta forzada de la forma yf (t) = kAt (un polinomio del
tiempo de primer grado), cuando la entrada es u = A (un polinomio del
tiempo de grado cero). N otese que la respuesta forzada correspondiente a la
ecuaci on diferencial en (3.4) tiene la forma yf (t) = kA a (un polinomio del
tiempo de grado cero) cuando la entrada es u = A (un polinomio del tiempo
de grado cero). A partir de estos dos ejemplos se llega a la siguiente conclusi on
importante:
La respuesta forzada de una ecuaci on diferencial de primer orden excita-
da con una entrada constante, tambien ser a constante si el polinomio carac-
terstico no tiene races en s = 0.
En la figura 3.2 se muestra la ubicacion en el plano s del polo de la funci
on
de transferencia en (3.39), es decir, el polo en s = 0. El lector debe aprender a
relacionar la ubicaci on en el plano s de los polos de la funci
on de transferencia
correspondiente con la forma en el tiempo de la soluci on y(t) correspondiente
(veanse los casos listados antes del ejemplo 3.2).
3.3. Ecuaci
on diferencial de segundo orden
Considere la siguiente ecuacion diferencial ordinaria, lineal, con coeficientes
constantes, de segundo orden:
para obtener:
n2 A y0 (s + 2n ) + y 0
Y (s) = k + 2 (3.41)
(s2 + 2n s + n2 )s s + 2n s + n2
Suponga primero que las condiciones iniciales son cero y 0 = 0, y0 = 0, enton-
ces:
A n2
Y (s) = k
(s2 + 2n s + n2 )s
3.3 Ecuaci
on diferencial de segundo orden 113
s2 + 2n s + n2 = (s a)(s a
) (3.42)
a = + jd , a = jd , j = 1
p
2
d = n 1 > 0, = n 0
s2 + 2n s + n2 = (s [ + jd ])(s [ jd ]) = (s )2 + d2
para obtener:
An2
k = B( + jd ) + C
+ jd
Multiplicando el lado izquierdo por ( jd )/( jd ):
( jd )
kAn2 = B + C + jBd
2 + d2
114 3 Base matem
atica: ecuaciones diferenciales
(s )2 + d2 (s )2 + d2
kAn2 (s ) kAn2 1
= +
( + d ) [(s ) + d ] ( + d ) [(s )2 + d2 ]
2 2 2 2 2 2
se obtiene:
( )
Bs + C kAn2 t kAn2
L1 = e cos( d t) + et sin(d t)
(s )2 + d2 2
2 + d d ( 2 + d2 )
kAn2 t
= e [ cos(d t) + sin(d t)]
2 + d2 d
s
2
kAn2 t
= e [sin cos(d t) + cos sin(d t)] +1
2 + d2 d
Notese que en el ultimo paso se han usado algunas relaciones de la figura
3.6(a). Por otro lado, se puede seguir reduciendo:
1
sin cos(d t) + cos sin(d t) = [sin( d t) + sin( + d t)] +
2
1
+ [sin(d t ) + sin(d t + )]
2
= sin(d t + )
3.3 Ecuaci
on diferencial de segundo orden 115
=! d
CA
1 ( )
2 + 1
r
CO
!d ( )
(a) (b)
sin(d t + ) = sin (d t + + )
= sin (d t + )
p
1 2
= arctan
por lo que ahora se puede reescribir:
r 2
( )
+1
Bs + C d
L1 = kAn2 et sin (d t + )
(s )2 + d2 2 + d2
p
2 + d2 t
= kAn2 e sin (d t + )
d ( 2 + d2 )
1
= kAn2 p et sin (d t + )
d 2 + d2
116 3 Base matem
atica: ecuaciones diferenciales
( )
Bs + C 1
L1
= kA p et sin (d t + ) (3.47)
(s )2 + d2 1 2
p
donde se ha usado = n y d = n 1 2 . Por otro lado, el valor de
D en (3.43) se obtiene facilmente multiplicando ambos miembros de dicha
expresion por el factor s y evaluando en s = 0:
kA n2
D= 2 = kA (3.48)
s + 2n s + n s=0
2
donde:
1 y0 (s + 2n ) + y 0
p(t) = L (3.50)
s2 + 2n s + n2
ya que U (s) = A/s. Notese que de acuerdo al parrafo que sigue a (3.50),
cuando las condiciones iniciales son diferentes de cero estas solo afectan a la
respuesta natural en el sentido de que su coeficiente se vera afectado por las
condiciones iniciales y0 y y 0 , pero la forma de la funcion del tiempo siempre
sera en t sin (d t + ). De hecho, p(t) forma parte de la respuesta natural
cuando las condiciones iniciales son diferentes de cero.
El lector puede revisar el procedimiento presentado despues de (3.44) para
comprobar que en el caso en que 1 < < 0 con k y n positivos, la solucion
en (3.51) o, equivalentemente, en (3.44) se transforma en:
cuando las condiciones iniciales son cero, es decir, yn (t) es una funcion
oscilatoria cuya amplitud no aumenta pero tampoco disminuye. Esto sig-
nifica que aunque y(t) no crece sin lmite, sin embargo no alcanzara de
manera permanente a la respuesta forzada yf (t).
4. El comportamiento obtenido cuando esta fuera del rango 1 < < 1 es
estudiado en las secciones que siguen bajo los casos en los que las races
del polinomio caracterstico son reales y repetidas o reales y diferentes.
Finalmente, notese que de acuerdo a lo anteriormente expuesto y a (3.49),
(3.50), la solucion de una ecuacion diferencial de segundo orden puede pre-
sentar oscilaciones sostenidas ( = 0) a un cuando la entrada o excitacion sea
cero (u = 0) si las condiciones iniciales son diferentes de cero. Esto explica el
118 3 Base matem
atica: ecuaciones diferenciales
Estudio gr
afico de la soluci
on
A continuacion se estudia de manera grafica la solucion presentada en
(3.51) para la ecuacion diferencial en (3.40). Recuerdese que se supone que
0 < 1, que n , k son constantes positivas y que las condiciones iniciales
son cero y0 = 0, y 0 = 0. Notese que la respuesta natural tiende a cero conforme
el tiempo crece si 0 < < 1:
!
en t
lm yn (t) = lm kA p sin (d t + ) = 0 (3.56)
t t 1 2
porque n > 0. Recuerdese tambien que en el caso en que las condiciones
iniciales no son cero, la respuesta natural tambien tiende a cero al crecer
el tiempo. Por tanto, para tiempos muy grandes la solucion de la ecuacion
diferencial es igual a la respuesta forzada u
nicamente:
lm y(t) = yf (t) = kA (3.57)
t
Esto significa que mientras mas rapido desaparezca la respuesta natural yn (t)
(lo cual ocurre conforme el producto n > 0 es mayor) mas rapido la solucion
completa y(t) alcanzara a la respuesta forzada yf (t). Tal como se menciono pa-
ra las ecuaciones diferenciales de primer orden, la respuesta natural permite
que la respuesta total y(t) transite de manera continua desde el valor dictado
por las condiciones iniciales (iguales a cero en este caso) hasta la respuesta
forzada yf (t).
Dos caractersticas importantes de la respuesta en (3.51) son el tiempo de
subida tr y el sobre paso Mp , las cuales se indican en la figura 3.7. La manera
de medir el sobre paso es mediante la expresion:
yp yf
Mp ( %) = 100 (3.58)
yf
donde yf = kA representa la respuesta forzada. Recuerdese que se supone que
y0 = y 0 = 0. A partir del analisis detallado de la expresion en (3.51) se puede
demostrar que [1], pag. 232:
" p !#
1 1 2
tr = arctan (3.59)
d
Mp ( %) = 100 e 1 2
y(t)
yp
yf
tr
tiempo
Funci
on de transferencia
kn2
Y (s) = U (s) (3.60)
s2 + 2n s + n2
o tambien:
Y (s) kn2
= G(s) = 2 (3.61)
U (s) s + 2n s + n2
y(t)
2kA
=0
0:1
0:2
0:3
0:4
0:5
0:6
kA
0:7
1:0
2:0
tiempo
(races del polinomio caracterstico). Por tanto, G(s) es una funcion de trans-
ferencia de segundo orden. De acuerdo a (3.42) estos polos estan ubicados en
s = a y en s = a , donde:
a = + jd , a
= jd
Es importante subrayar que estos polos son complejos conjugados solo bajo la
suposicion de que 1 < < 1. El lector puede verificar que los polos son reales
si no satisface la condicion anterior. Ese caso es estudiado en las secciones
que siguen.
De acuerdo al estudio realizado en las secciones previas se pueden hacer
las siguientes afirmaciones en lo que se refiere a la salida y(t) de la funcion de
transferencia G(s):
Si los polos tienen parte real negativa = n < 0, es decir > 0,
entonces al crecer el tiempo la respuesta natural yn (t) desaparece y la
respuesta total y(t) alcanza a la respuesta forzada yf (t). Cuando la entrada
es una constante u(t) = A, entonces la respuesta forzada es yf (t) = kA.
La rapidez con la que yn (t) desaparece es mayor conforme el producto n
es mayor. En la figura 3.10 se muestra como la respuesta del sistema de
segundo orden permanece envuelta por dos funciones exponenciales cuya
rapidez de decaimiento esta determinada por el producto n .
3.3 Ecuaci
on diferencial de segundo orden 121
y(t)
!n 3 !n 2 !n 1
kA
tiempo
y(t)
e ! nt
kA 1 + p 1 2
kA
e ! nt
p
kA 1 1 2
tiempo
kn2 A kn2
lm y(t) = lm sY (s) = lm s = A = A(3.62)
t s0 s0 s2 + 2n s + n2 s n2
Im (s)
j! d
!n
! n Re (s)
Figura 3.11. Uno de los polos de G(s) en (3.61) cuando 0 < < 1.
de transferencia en (3.61).
Si 1 < < 0 entonces la parte real de los polos = n es positiva por
lo que la respuesta total y(t) crece sin lmite al crecer el tiempo. Notese
que esto implica que los polos complejos conjugados se encuentran en el
lado derecho del plano s.
El coeficiente recibe el nombre de factor de amortiguamiento porque es
el que determina que tan oscilatoria es la respuesta y(t) (vease (3.59), para
el sobre paso, y la figura 3.8).
3.3 Ecuaci
on diferencial de segundo orden 123
m
x = Kx bx + f (3.63)
b K 1
x
+ x + x = f
m m m
lo cual, a su vez, se puede escribir como la expresi
on en (3.40):
+ 2n x + n2 x = kn2 f
x (3.64)
b K 1
2n = , n2 = , k =
m m K
donde x juega el papel de la salida y mientras que f juega el papel de la entrada
u. Por tanto, x(t) tiene la forma mostrada para y(t) en (3.51) si 0 < 1
cuando las condiciones iniciales son cero. N otese que todas las constantes
en (3.63) son positivas: la masa m, el coeficiente de fricci on viscosa b y la
constante de rigidez del resorte K s olo pueden ser positivas. De acuerdo a
(3.64), el factor de amortiguamiento est a dado como:
b
= (3.65)
2 mK
Por tanto, si se aplica una fuerza constante f = A se pueden usar los resul-
tados en las figuras 3.7, 3.8 y 3.9 para concluir lo siguiente.
Si no hay friccion, es decir si b = 0, entonces = 0 y el carrito osci-
1
lara permanentemente alrededor de la posici on xf = kA = K A. En el
caso en que la fuerza externa aplicada sea cero f = 0 pero alguna de las
condiciones iniciales x 0 o x0 sea diferente de cero, entonces el carrito tam-
bien oscilar
a permanentemente pero alrededor de la posici on x = 0, como
lo describe la funci
on p(t) que aparece en (3.49).
124 3 Base matem
atica: ecuaciones diferenciales
Si la fricci
on b > 0 se incrementa poco a poco, entonces > 0 tambien
se incrementa y el carrito oscilar a cada vez menos veces porque el sistema
masa-resorte-amortiguador estar a cada vez mas amortiguado (cuando =
1 el sistema tiene dos races reales repetidas y ya no presenta oscilaciones,
como se ve en las secciones siguientes). En todos estos casos, cuando el
1
carrito deje de oscilar se detendr a en la posici
on x = xf = kA = K A. En
el caso en que la fuerza externa aplicada sea cero f = 0 pero alguna de
las condiciones iniciales x 0 o x0 sea diferente de cero, entonces el carrito
tambien oscilara cada vez menos (conforme > 0 se aproxima a 1) y
se detendr a en la posicion x = 0. De nuevo, esto est a descrito por la
funcion p(t) que aparece en (3.49). N otese que de acuerdo a (3.65) el
amortiguamiento tambien aumenta si decrece la masa m o la constante
de rigidez del resorte K. Para entender la raz on de esto considerese el
caso contrario: valores mayores de m y K producen fuerzas m as grandes
(mx para la masa y Kx para el resorte) que, por tanto, pueden vencer con
mayor facilidad a la fuerza de fricci on (bx)
y por ello el movimiento puede
permanecer durante m as tiempo.
Debido a que el tiempo de subida qtr depende inversamente de la frecuencia
natural no amortiguada n = K m , entonces un resorte mas rgido (con
una K mayor) o una masa m as peque
na producen respuestas m as r
api-
das (tr menores). Por otro lado, el hechoq que la frecuencia de oscilaci on
p K
2
este dada por d = n 1 con n = m es el principio mediante el
cual se ajusta la cadencia o ritmo de un reloj mec anico mediante el uso de
un sistema masa-resorte-amortiguador oscilatorio: si el reloj se atrasa
hay que aumentar su frecuencia de oscilaci on (ritmo o cadencia), lo cual
se consigue tensando m as el resorte para aumentar K. Se hace lo contrario
si el reloj se adelanta.
El caso cuando 1 < < 0 no es posible en este ejemplo porque todas las
constantes b, m y K son positivas. Sin embargo, la situaci
on en la que 1 <
< 0 ocurre con frecuencia en sistemas de control realimentados, debido a la
interconexi
on de diversos componentes.
K
m f
b
k P (s)
Y (s) = U (s) n
sn + an1 sn1 + + a1 s + a0 s + an1 sn1 + + a1 s + a0
k A P (s)
Y (s) =
N (s) s N (s)
N (s) = sn + an1 sn1 + + a1 s + a0
Suponga por lo pronto que todas las condiciones iniciales son cero. Esto sig-
nifica que P (s) = 0 y se puede escribir:
k A
Y (s) =
N (s) s
Suponga que N (s) solo tiene races reales, distintas entre s y, ademas, distintas
de s = 0, es decir:
N (s) = (s p1 )(s p2 ) (s pn )
k A c1 c2 cn f
Y (s) = = + + + + (3.67)
N (s) s s p1 s p2 s pn s
Adem
as, ambas races son negativas:
p1 < 0, p2 < 0
y(t)
0<<1
kA
=1
>1
tiempo
(a) Respuesta en el tiempo
Im (s)
p2 ! n p1 Re (s)
(b) Ubicaci
on de los polos cuando > 1
Suponga por lo pronto que todas las condiciones iniciales son cero, es decir
P (s) = 0, entonces:
k A
Y (s) =
N (s) s
130 3 Base matem
atica: ecuaciones diferenciales
Suponga que N (s) solo tiene una raz la cual esta repetida n veces y que es
diferente de cero, es decir:
N (s) = (s p)n
k A cn cn1 cn2 c2 c1 f
= + + + + + +
N (s) s (s p)n (s p)n1 (s p)n2 (s p)2 sp s
(3.72)
kA
cn =
p
Por otro lado, las constantes ci , i = 1, . . . , n 1, se pueden calcular multipli-
cando ambos miembros de (3.72) por el factor (s p)n , derivando n i veces
respecto a s y evaluando en s = p para obtener:
dni kA
ci = , i = 1, . . . , n 1
dsni s s=p
se obtiene finalmente:
tn1 tn2 tn3
y(t) = cn ept + cn1 ept + cn2 ept +
(n 1)! (n 2)! (n 3)!
+ c2 t ept + c1 ept + f
3.5 Races reales repetidas 131
p = n
la cual est
a repetida dos veces. La posici
on x(t) del carrito evoluciona exac-
tamente de la misma forma en que lo hace y(t) en la figura 3.13. Este tipo
de respuesta puede ser de mucha utilidad si se desea que el movimiento del
carrito sea r
apido y sin oscilaciones.
Ejemplo 3.12 Suponga que en la figura 3.12 no existe ning un resorte ni
amortiguador, es decir, suponiendo que b = 0 y K = 0 la ecuaci
on (3.63) se
reduce a:
m
x=f (3.76)
x(t)
"0
m t1
"0 2
2m
t
t1 t [s]
Figura 3.14. Posicion del carrito de la figura 3.12 cuando es perturbado y no existe
ning
un resorte ni amortiguador.
p1 = p2 = n > 0
k A P (s)
Y (s) =
N (s) s N (s)
N (s) = sn + an1 sn1 + + a1 s + a0
Suponga por lo pronto que todas las condiciones iniciales son cero, es decir
P (s) = 0, entonces:
k A
Y (s) =
N (s) s
Suponga que N (s) tiene n/2 races complejas conjugadas diferentes, es decir:
i=n/2
Y
N (s) = [(s i )2 + di
2
]
i=1
donde f es una constante que se puede calcular como en las secciones previas,
es decir:
kA
f=
N (s) s=0
mientras que Fi y Ci son constantes que pueden calcularse del mismo modo
en que se calculan B y C en la seccion 3.3. Por tanto, usando (3.47) se puede
escribir:
donde:
1 P (s)
q(t) = L
N (s)
Es sencillo verificar que q(t) esta dado por funciones de la misma forma que
las funciones que componen a yn (t). De hecho q(t) forma parte de la respuesta
natural cuando las condiciones iniciales son diferentes de cero.
Es importante senalar que si 1 o 1 entonces se obtienen races
reales y se recupera uno de los casos de las secciones 3.4 o 3.5.
Es conveniente aclarar que se asume que las races complejas siempre apa-
recen con sus parejas conjugadas debido a que esa es la u nica manera en
que el producto (s a)(s a ), con a un numero complejo y a su conjugado,
resulte en un polinomio cuyos coeficientes son todos n umeros reales. Enton-
ces, cualquier polinomio que contenga races complejas y tambien sus parejas
conjugadas solo tendra coeficientes reales. Este es un aspecto importante en
ingeniera de control donde las ecuaciones diferenciales representan sistemas
fsicos que, por tanto, solo manejan variables y parametros reales. Por ejemplo,
el polinomio caracterstico del sistema masa-resorte-amortiguador presentado
b
en (3.64) esta dado como s2 + m s+ Km donde los coeficientes representan la
masa, el coeficiente de friccion viscosa y la constante de rigidez del resorte.
Es claro que ninguno de estos parametros pueden tener parte imaginaria pues
representan propiedades cuyos efectos pueden ser apreciados experimental-
mente.
Ejemplo 3.14 En el ejemplo 2.5 del captulo 2 se ha obtenido el modelo
matem
atico del sistema compuesto por dos cuerpos y tres resortes mostrado
3.6 Races complejas conjugadas diferentes 137
F(t) x1 K2 x2
K1 K3
m1 m2
Con el fn de simplificar el a
lgebra correspondiente, supongamos que b = 0,
entonces:
1 2 K2 +K3
m1 s + m2
X1 (s) = 2
F (s)
s2 + mK2
s2 + K2 +K3 K2
1 m2 m1 m2
o, equivalentemente:
138 3 Base matem
atica: ecuaciones diferenciales
1 K2 +K3
m1 s2 + m2
X1 (s) = F (s) (3.80)
K2 K2 +K3 K2 K3
s4 + m1 + m2 s2 + m1 m2
Por otro lado, notese que si se multiplican dos polinomios de segundo orden
con coeficientes reales a, d, c, e, se obtiene:
es posible si y s
olo si:
a+d=0 y cd + ea = 0 (3.82)
K2 K2 + K3
2e a2 = +
m m2
r1
K2 K3
e=
m1 m2
Combinando ambas expresiones se obtiene:
r
2 K2 K3 K2 K3 K2
a =2 + <0 (3.83)
m1 m2 m1 m2 m2
ejemplo). Entonces, el u
nico caso posible es i) d = a = 0, el cual, cuando es
aplicado al polinomio en el denominador de (3.80), implica:
K2 K2 + K3
c+e = + (3.84)
m1 m2
K2 K3
ce = (3.85)
m1 m2
K2
Definiendo f = m 1
+ K2m+K2
3
yg= m K2 K3
1 m2
se tiene, de la segunda expresi
on,
c = g/e y, sustituyendo en la primera expresi on, g/e + e = f . Finalmente,
multiplicando este resultado por e se obtiene:
e2 ef + g = 0
c2 cf + g = 0
K2 +K3 K2 K2 +K3
m2 m1 m2 = 2 K2m+K
2
3
> 0 esto no asegura que e2 = c2 > 0
4K 2
por el termino adicional m1 m2 2 dentro del radical. Este problema se resuelve
recurriendo a (3.85): los valores de e y c deben tener el mismo signo porque
el miembro derecho de (3.85) es positivo. Entonces, ambos e y c deben ser
positivos porque de acuerdo a los u nicos casos posibles a) y b) se debe usar
e = e1 > 0 (lo cual fuerza que c = c2 > 0) o c = c1 > 0 (lo cual fuerza
que e = e2 > 0). Finalmente, el razonamiento al principio de este p arrafo
tambien muestra que, a excepci on de algunos posibles valores muy especiales
para m1 , m2 , K2 , K3 , ambos e1 = c1 y e2 = c2 , son diferentes de K2m+K
2
3
. Esto
asegura que no existen cancelaciones entre los polinomios del numerador y
del denominador en (3.86) y el hecho de que e y c sean positivos y diferentes
asegura que el polinomio caracterstico de (3.86) tiene dos pares de races
complejas (imaginarias puras) conjugadas diferentes.
Suponga por lo pronto que todas las condiciones iniciales son cero, es decir
P (s) = 0, entonces:
k A
Y (s) =
N (s) s
Suponga que el polinomio caracterstico N (s) tiene dos races complejas con-
jugadas repetidas n/2 veces (para tener un total de n races):
k A F1 s + C1 F2 s + C2
= 2
+ +
N (s) s 2
[(s ) + d ] n/2 [(s )2 + d2 ]n/21
Fn/21 s + Cn/21 Fn/2 s + Cn/2 f
+ 2 2 2
+ 2 2 +
[(s ) + d ] (s ) + d s
En el caso de races reales se encontro que cuando estas son sencillas intro-
ducen terminos de la forma ept donde p es la raz en cuestion. Cuando dicha
raz es repetida j veces, se encontro que los terminos introducidos se transfor-
man en la forma tj1 ept . En el caso de races complejas conjugadas repetidas
ocurre algo similar. En este caso no se hara una demostracion formal de lo
que sigue debido a la complejidad del tema. Solo se presenta la idea intuitiva
de la solucion para comprender el porque de la forma obtenida.
Usando la transformacion inversa obtenida en (3.47) se encontro que una
raz compleja conjugada no repetida (sencilla) introduce en el tiempo un termi-
no de la forma:
( )
Bs + C
L1
= et sin(d t + )
(s )2 + d2
Notese que en este caso la respuesta natural, yn (t), crece sin lmite confor-
me el tiempo crece en el caso en que la raz se repita al menos dos veces.
Pero si la raz es sencilla entonces yn (t) esta formada por una funcion os-
cilatoria cuya amplitud no crece ni disminuye, es decir, aunque y(t) nunca
alcanzara de manera permanente a yf (t) la solucion y(t) permanece sin
crecer demasiado.
Finalmente, es sencillo verificar que si las condiciones iniciales no son cero
entonces:
donde:
1 P (s)
q(t) = L
N (s)
La funcion q(t) esta dada por las mismas funciones que componen a yn (t) y
las tres posibilidades enumeradas previamente tambien siguen vigentes en este
caso. De hecho q(t) forma parte de la respuesta natural cuando las condiciones
iniciales son diferentes de cero.
Ejemplo 3.15 (Tomado de [2], cap. 4)Considere el sistema masa-resorte
mostrado en la figura 3.12. Suponga que no existe ning un amortiguador, es
decir que b =q0, y que se aplica una fuerza externa con la forma f = F0 sin(t)
K
donde = m . Se desea describir la posici
on del carrito x(t) si x(0) = 0
y x(0)
= 0. No es necesario calcular el valor numerico de las constantes que
aparecen en x(t).
Soluci
on. La ecuaci
on diferencial que describe la situaci
on en la figura 3.12
cuando b = 0 es:
r
K
m x + Kx = f, f = F0 sin(t), = (3.89)
m
Aplicando la transformada de Laplace a (3.89):
K 1
s2 X(s) + X(s) = F (s)
m m
N
otese que esta expresi
on tiene la forma:
F0
F (s) =
s2 + 2
Entonces:
3 F0
X(s) =
(s2 + 2 )2
Expandiendo en fracciones parciales:
As + B Cs + D
X(s) = + 2
(s2 + 2 )2 s + 2
Por tanto, de acuerdo a lo expuesto en la presente secci
on, se puede escribir:
x(t) = 1 t sin(t + 1 ) + 2 sin(t + 2 )
Donde 1 , 2 , 1 , 2 son algunas constantes
q reales. En la figura 3.16 se mues-
tra una grafica de x(t) cuando = K m = 10[rad/s], m = 1[Kg] y F0 = 1[N].
A este fen omeno se le conoce como resonancia y significa que aunque la
entrada es acotada la salida puede alcanzar valores muy grandes a pesar de
que el polinomio caracterstico no tiene races con parte real positiva ni repe-
tidas en s = 0, por lo que este es un fen
omeno diferente al que se presenta en
el ejemplo 3.12. N otese que el problema de la resonancia aparece cuando un
sistema (ecuaci on diferencial) est
a mal amortiguado ( 0) y se excita con
una senal oscilatoria cuya frecuencia es igual (o muy cercana) a la frecuencia
natural del sistema. Debido a que la posici on del carrito crece sin lmite la
resonancia es un fen omeno que se considera peligroso.
x [m]
0,8
0,6
0,4
0,2
-0,2
-0,4
-0,6
-0,8
0 2 4 6 8 10 12 14 16 18 20
t [s]
Y (s) B(s)
= G(s) = (3.91)
U (s) N (s)
y del factor 1s introducido por U (s) = A/s. La variable y(t) se conoce como
la salida y u(t) como la entrada de la funcion de transferencia G(s). El orden
de G(s) se define como el grado del polinomio caracterstico N (s). Como un
polinomio de grado n tiene n races, entonces una funcion de transferencia de
orden n tiene n polos (vease el parrafo que sigue a (3.21)). Por otro lado, las
races del polinomio B(s) se conocen como los ceros de la funcion de transfe-
rencia G(s). Si B(s) es de grado m entonces G(s) tiene m ceros. Recuerdese
la condicion n m impuesta al principio de esta seccion. En las secciones
previas se ha visto que siempre se puede escribir y(t) = yn (t) + yf (t) y que
el efecto de las condiciones iniciales siempre se pueden considerar como parte
de yn (t). Tambien se ha visto que yn (t) depende de la naturaleza de las races
del polinomio caracterstico N (s), es decir, si son reales, sencillas o repetidas,
si son complejas conjugadas, sencillas o repetidas y si tienen parte real nega-
tiva, positiva o cero. Si yn (t) tiende a cero conforme el tiempo crece se dice
que la ecuacion diferencial en (3.90) o, equivalentemente, que la funcion de
transferencia G(s) es estable. Resumiendo todos los casos estudiados en las
secciones anteriores ahora se puede afirmar lo siguiente:
qi = k1 u (3.93)
u = kp (hd h) (3.94)
dh
+ b1 h = b2 h d (3.95)
dt
b1 = a + kk1 kp > 0, b2 = kk1 kp > 0, b1 > b2
b2 hd b1 t b2 hd
h(t) = e + + h0 eb1 t , h0 = h(0)
b1 b1
N
otese que, debido a que b1 > b2 > 0:
b2 h d
lm h(t) = < hd (3.96)
t b1
es el valor que alcanza en estado estacionario el nivel de agua en el tanque y es
menor (diferente) que el nivel deseado hd . Esto puede explicarse usando (3.92)
del siguiente modo. La funci on de transferencia de la ecuacion diferencial en
(3.95) es:
H(s) b2
G(s) = =
Hd (s) s + b1
Entonces:
b2 h d b2
lm h(t) = lm sH(s) = lm s = h d < hd
t s0 s0 s + b1 s b1
Ahora bien, este resultado tambien puede explicarse usando la experiencia
cotidiana del siguiente modo. Supongase por un momento que lmt h(t) =
hd . Entonces, de acuerdo a (3.93) y (3.94):
qi = k1 kp (hd h) = 0
el flujo de agua que entra al tanque es cero porque h = hd . Como el agua con-
tinua fluyendo a traves de la v
alvula de salida, qo 6= 0, lo anterior resultar
a en
h < hd de nuevo. Esto significa que no es posible mantener h = hd de manera
permanente, lo que justifica el resultado en (3.96).
Ahora suponga que se usa (3.93) junto con:
Z t
u = kp (hd h) + ki (hd h(r))dr (3.97)
0
lo cual es completamente cierto si a + kk1 kp > 0 y kk1 ki > 0 (se deja como
ejercicio consultar la secci on 4.2.1, del captulo 4, para verificar que estas
condiciones aseguran que las dos races del polinomio caracterstico s2 + (a +
kk1 kp )s + kk1 ki tienen parte real negativa). De nuevo, este resultado puede
ser explicado usando la experiencia cotidiana. Dado que h = hd en estado
estacionario, entonces es de esperarse que la se nal de error e = hd h tenga un
comportamiento
Rt como el mostrado en la figura 3.17. Recuerdese que la integral
0
(h d h(r))dr es el area debajo de la curva mostrada en la figura 3.17 la cual
es positiva. Esto significa que cuando h = hd la integral permanece constante
(el integrando vale cero) en un valor positivo. Entonces, de acuerdo a (3.93)
y (3.97) contin ua entrando agua al tanque con un flujo qi que permanece
constante en este valor positivo multiplicado por k1 ki > 0. Este flujo resulta
ser exactamente igual al flujo del agua que sale qo de manera que el nivel de
agua h = hd permanece igual al nivel deseado, tal como lo predice (3.98).
m
x + bx = f
f = kp (xd x) (3.99)
3.8 Una ecuaci
on diferencial general 149
e [m]
t [s]
Rt
Figura 3.17. La integral 0 (hd h(r))dr es el
area sombreada debajo de la curva
definida por e = hd h. Se supone que hd > h(0).
m
x + bx = kp (xd x)
y agrupando:
b kp kp
x
+ x + x = xd
m m m
Usando la transformada de Laplace bajo la suposicion de condiciones iniciales
iguales a cero se encuentra la siguiente funci
on de transferencia:
kp
X(s) m
G(s) = = kp
Xd (s) s2 + b
ms + m
es igual a la posici
on deseada xd . La raz on de este resultado puede explicar-
se usando, de nuevo, la experiencia: cuando x = xd la fuerza producida de
acuerdo a (3.99) es f = 0 por lo que el carrito puede detenerse y permanecer
en dicha posici
on. M as a
un, si x < xd entonces f > 0 y el carrito incrementa
su posici
on x acercandose a xd (vease la figura 3.12 para recordar la direcci
on
en se aumenta x y en la que f es positiva). Si x > xd entonces f < 0 y el
carrito decrementa su posicion x acerc andose, de nuevo, a xd .
150 3 Base matem
atica: ecuaciones diferenciales
3.9.1. Cancelaci
on polo-cero y modelos reducidos
k(s d)
Y (s) = U (s) (3.100)
(s p1 )(s p2 )
k(s d) A B C D
Y (s) = = + + (3.101)
(s p1 )(s p2 ) s s p1 s p2 s
d A
B k 1 (3.113)
a RC
dA
C k 2 D (3.114)
a
Por tanto, usando (3.109) y despreciando C se obtiene:
d A 1 t d A
y(t) k 1 e
RC + k
1 a (3.115)
a RC RC
es la solucion de:
kd
Y (s) = 1 U (s) (3.117)
a(s + RC )
3.9 Polos y ceros en sistemas de orden superior 153
con U (s) = A/s. Por tanto, se puede utilizar (3.117) en lugar de (3.108). La
1
condicion a RC se interpreta diciendo que el polo en s = a es muy
1 1
rapido comparado con el polo en s = RC . Al polo en s = RC se le conoce
como el polo dominante porque su efecto es el que sobresale en la respuesta del
sistema. Por tanto, se concluye. Si un polo es mucho mas rapido que los otros,
entonces se puede obtener una funcion de transferencia de orden reducido si
se desprecia el polo rapido y se conservan los polos dominantes (lentos). Un
criterio aceptado es que los polos rapidos (los que se pueden despreciar) tienen
una parte real de al menos 5 veces la parte real de los polos dominantes (los que
se conservan). Notese, sin embargo, que no es cuestion de simplemente hacer
desaparecer el factor s + a en el denominador de la funcion de transferencia:
la constante a aun aparece dividiendo en (3.117) porque es necesaria para
conservar la ganancia en estado estacionario de la funcion de transferencia en
(3.108).
Una manera sencilla de obtener (3.117) a partir de (3.108) es la siguiente:
k d k d
Y (s) = 1 s + a U (s) = 1 U (s) (3.118)
s + RC s + RC a( a1 s + 1)
1
Si a es muy grande se puede suponer que as 1 y, por tanto:
kd
Y (s) 1 U (s) (3.119)
a(s + RC )
dh
+ ah = kqi
dt
1 1
a= , k=
RC C
Usando la transformada de Laplace (3.1) y considerando todas las condiciones
iniciales iguales a cero, no es difcil comprobar que las ecuaciones correspon-
dientes adquieren la forma:
1/La
I(s) = (V (s) ke (s)) (3.121)
s+ R a
La
km /J
(s) = I(s) (3.122)
s+ BJ
1/C
H(s) = 1 Qi (s) (3.123)
s + RC
km /J 1/La
(s) = (V (s) ke (s))
s+ B Ra
J s + La
k /J 1/La
(s) = B
Jm (V (s) ke (s))
Bs + 1
Ra La
J La Ra s + 1
1 km
(s) = (V (s) ke (s))
s+ B
J
JR a
Para continuar, se pueden reagrupar los terminos que contienen (s) para
reescribir esta expresi
on del siguiente modo:
km
JR V
(s) = (s)
B km ke
s+ J + JRa
o bien:
1 km
C JR V
H(s) = 1 (s)
s+ B km ke
RC J + JRa
Z 1(s) Z 5(s)
Z 3(s)
(a)
Z 3(s)
(b)
I 3(s)
I a (s)
Z 3(s)
(c)
V1 (s) V2 (s)
I1 (s) = , I2 (s) = (3.126)
Z1 (s) Z5 (s)
Z 1(s)
+
Z 3(s)
+
V 1(s) Z 2(s) Z 4(s)
+ V (s)
2
(a)
Z 1(s) Z 3(s)
(b)
Z 3(s)
(c)
Z a (s) Z 3(s) I(s)
V 3(s) + + V 2(s)
(d)
Z3 (s)
Vz3 (s) = (V3 (s) V2 (s))
Za (s) + Z3 (s)
Z3 (s)
= (Za (s)I1 (s) V2 (s))
Za (s) + Z3 (s)
Z3 (s) Za (s)
= V1 (s) V2 (s)
Za (s) + Z3 (s) Z1 (s)
di
L = v v0 sign(i) + E(t) (3.128)
dt
dv
C =i (3.129)
dt
dv0 v0
C0 = abs(i) (3.130)
dt R
donde:
+1, i > 0
sign(i) = (3.131)
1, i < 0
Simplificando se encuentra:
z1 = z2 z3 sign(z1 ) + u (3.133)
z2 = z1 (3.134)
z3
z3 = abs(z1 ) (3.135)
Q
donde el punto representa lapderivada respecto al tiempo normalizado
mientras que = C0 /C, Q = R C/L y u es una variable que solo toma los
valores de +1 (cuando E(t) = +E) y 1 (cuando E(t) = E). Normalmente,
el valor del capacitor de salida C0 se selecciona muy grande comparado con
el valor del capacitor resonante C porque esto asegura que z3 , es decir v0 , se
mantendra aproximadamente constante durante varios ciclos de la corriente
162 3 Base matem
atica: ecuaciones diferenciales
Esto significa que si la solucion de (3.138) se dibuja sobre un plano donde el eje
horizontal es z4 y el eje vertical es z4 , entonces se obtiene un crculo centrado
en el origen, (z4 , z4 ) = (0, 0), que tiene
p un radio (constante) determinado
por las condiciones iniciales, es decir z40 2 + z 2 . Mas a
un, de acuerdo a los
40
cambios de variable z4 = z2 ve y z4 = z1 , se concluye que si la solucion
de (3.136), (3.137), se dibuja sobre un plano donde el eje horizontal es z2
y el eje vertical es z1 , entonces se obtiene un crculo centradop en el punto
(z2 , z1 ) = (ve , 0), que tiene un radio (constante) igual a (z20 ve )2 + z10 2
Z1
Z2
( 1 + Z3; 0) (1 Z3; 0)
z1
s=0
z 1 = z 2
z 1 > z 2
z1 > 0
Q1; Q3
z 1 < z 2
z1 > 0
D2; D4
1 z3 1 + z3 1 z3 1 + z3 z2
D1; D3 Q2; Q4
z 1 > z 2 z 1 < z 2
z1 < 0 z1 < 0
Notara que se empezara a obtener una trayectoria con forma de espiral, como
la lnea punteada mostrada en la figura 3.21, que finalmente terminara en una
trayectoria cerrada. El hecho de que dicha trayectoria cerrada sea alcanzada a
partir de cualquier punto inicial sobre el plano z2 z1 es una muestra de que
tal trayectoria cerrada es estable (o atractiva). El lector puede proceder de la
misma manera usando un valor z3 > 1 para verificar que en este caso no se
consigue alcanzar ninguna trayectoria cerrada. Esto es un indicativo de que
en un convertidor resonante serie no es posible obtener un voltaje de salida v0
que sea mayor que el voltaje de suministro E, es decir no es posible obtener
un valor de z3 mayor que la unidad (recuerdese que v0 = Ez3 ). Si se desea
un voltaje de salida tal que v0 > E se puede utilizar el convertidor resonante
paralelo mostrado en la figura 2.44(b).
Por otro lado, en la figura 3.21, los arcos de crculo con centro en
(z2 , z1 ) = (1 z3 , 0) y (z2 , z1 ) = (1 z3 , 0) se dibujan de manera conse-
cutiva. Entonces, si ambos estuvieran centrados en el origen (z2 , z1 ) = (0, 0),
juntos completaran un angulo de 180 grados. Sin embargo, los centros de estos
arcos de crculo estan colocados en el cuadrante donde z2 tiene signo contrario
al cuadrante donde esta ubicado dicho arco de crculo. Esta observacion junto
con el uso de geometra basica permite concluir que los arcos de crculo con
centro en (z2 , z1 ) = (1 z3 , 0) y (z2 , z1 ) = (1 z3 , 0) describen juntos un
angulo menor que 180 grados a pesar de que juntos completan un semiciclo
completo en ambas variables resonantes z2 , z1 . Por tanto, los cuatro arcos de
crculo que componen a la trayectoria cerrada completa en la figura 3.21 des-
criben juntos un angulo menor que 360 grados a pesar de que describen un
ciclo completo en ambas variables resonantes. Notese que la frecuencia con la
que se recorren los arcos de crculo mencionados sigue siendo n = 1 pero, de
acuerdo a lo recien explicado, ahora se requiere menos tiempo para realizar
una oscilacion completa porque los cuatro arcos de crculo describen un angu-
lo menor a la misma velocidad angular. Por tanto, la frecuencia de operacion
o de las variables resonantes es mayor que n = 1, es decir o > n . Usando
(3.132) se puede comprobar que:
n 1 o
nt = = , ot = > nt
LC LC LC
donde nt y ot representan, respectivamente, la frecuencia de resonancia del
circuito y la frecuencia de operacion del circuito, expresadas respecto a la
base de tiempo real t. Por tanto, se concluye que el circuito debe operar a
frecuencias mayores que la frecuencia de resonancia para entregar voltajes de
salida menores que el voltaje de suministro v0 < E (es decir z3 < 1). Por esta
razon, cuando los transistores se activan de acuerdo a la regla u = sign(s),
el circuito en las figuras 2.44(a) y 2.45 se denomina convertidor electronico
de potencia de CD a CD tipo resonante serie de alta frecuencia. Para mayor
informacion sobre convertidores de potencia de CD a CD tipo resonantes
se recomienda consultar [8], donde se dise nan, construyen, controlan y se
prueban experimentalmente (ambos tipos: serie y paralelo), y [7], donde se
3.12 Resumen del captulo 167
sea una constante, esto esta determinado por los terminos independientes
de los polinomios de la funcion de transferencia.
3. Respuesta transitoria. Este aspecto se refiere a darle la forma adecuada a
la respuesta natural y depende de la seleccion adecuada de los polos de la
funcion de transferencia en lazo cerrado.
10. Dibuje el plano complejo s y sobre el ubique la zona donde deben colocarse
i) los polos reales que corresponden a respuestas rapidas, ii) los polos
complejos conjugados que corresponden a respuestas poco oscilatorias, iii)
los polos complejos conjugados que corresponden a respuestas rapidas, iv)
los polos complejos conjugados que corresponden a respuestas que oscilan
de manera permanente y v) los polos inestables.
h
[m]
t [s]
x
[m]
t [s]
donde U (s) = A/s, con A una constante. Realice simulaciones en las que
use valores de a que aumentan desde valores peque
nos hasta valores muy
3.14 Ejercicios propuestos 171
y + 4y + y = 2
y + y + 10y = 3t
y + 7y = 2 sin(t)
y + 3y + 2y = 2et
y + 2y + 2y = 2 cos(2t) 4 sin(2t)
y 4y = 8t2 4
y (4) + 11y (3) + 28
y = 10
8. Para las siguientes ecuaciones diferenciales diga cuanto vale lmt y(t).
Notese que el mismo resultado se obtiene sin importar el valor de las
condiciones iniciales.
y (3) + y 2y + y = 8
y + y 10y = cos(5t)
y + ay + by = kA, a > 0, b > 0, a, b, A, k son constantes
3y + 3y + y6 = 2A
y + 4y + 10y = 9A
A es una constante.
10. Considere la solucion de una ecuacion diferencial de segundo orden ante
una entrada constante o tipo escalon, considerando condiciones iniciales
cero, presentada en (3.44). Obtenga los maximos de dicha respuesta y
resolviendo para el primero de dichos maximos demuestre que el sobre
paso se puede calcular como:
Mp ( %) = 100 e 1 2
!d + I q + vq
PI 1= M PI !
vd (dq) 1 PMSM
PI
Id dq
I d = 0
+ Iq
Los sistemas de control en lazo cerrado pueden ser complejos. En ellos, los
componentes del sistema de control se conectan entre s dando origen a los
diagramas de bloques. La manera de trabajar estos problemas es simplificando
los diagramas de bloques para obtener una u nica funcion de transferencia de
lazo cerrado que representa al sistema completo. Una vez conseguido esto, se
puede determinar la estabilidad del sistema completo estudiando la ubicacion
de los polos de la funcion de transferencia de lazo cerrado. Este problema
puede ser complejo incluso para sistemas de segundo orden. Por esta razon,
se propone un metodo que determina la estabilidad simplemente analizando
los signos de los coeficientes del polinomio caracterstico. Para sistemas de
orden superior se utiliza el criterio de Routh. Por otro lado, el estudio del
error en estado estacionario permite determinar que tipo de controlador se
debe utilizar para conseguir que la respuesta forzada sea igual o muy parecida
a la entrada del sistema en lazo cerrado (o valor deseado a la salida).
178 4 Criterios de estabilidad y error
para obtener:
Ejemplo 4.2 Suponga que se tienen dos sistemas conectados en paralelo co-
mo se muestra en la figura 4.2. Entonces se puede escribir:
para obtener:
Y1 (s) + Y2 (s) = (G1 (s) + G2 (s))U (s) G(s) = G1 (s) + G2 (s)
Y1 (s) + Y2 (s) = G(s)U (s)
Y 1(s)
G 1(s)
U(s)
+
Y 1(s) + Y 2(s)
+
G 2(s)
Y 2(s)
U(s)
) G(s) Y 1(s) + Y 2(s)
donde:
C(s) G(s)
M (s) = = (4.1)
R(s) 1 + G(s)H(s)
se conoce como la funci
on de transferencia de lazo cerrado.
nke s
(a)
T p ( s)
I (s) + (s)
+ 1 I(s) + 1
KAp nkm
Ls+R Js 2 +bs
nke
KA p
s
(b)
T p ( s)
I (s) + KA p I(s) + 1 (s)
nkm
Ls+R+KA p Js 2 +bs
nke
KA p
s
(c)
G(s) (s)
M (s) = = = G1 (s)
1 + G(s)H(s) I (s)
con:
KAp nkm nke s
G(s) = , H(s) =
(Ls + R + KAp )(Js2 + bs) KAp
para obtener:
nkm
G1 (s) = h i (4.3)
sL+R n2 km ke
KAp + 1 (sJ + b) + KAp s
4.1 Diagramas de bloques 183
I(s) + (s)
KA p 1
Ls+R+KA p
nkm Js 2+bs
nk e
KA p
s
(a) Diagrama de bloques cuando Tp (s) = 0.
T p(s) 1
(s)
Js 2+bs
nk mKA p nk e
Ls+R+KA p KA p
s
(b) Diagrama de bloques cuando I (s) = 0.
I(s) G1(s)
+
(s)
+
T p(s) G2(s)
(c) Diagrama de bloques equivalente a cualquiera de los diagra-
mas mostrados en la figura 4.4.
Por tanto, cualquiera de los diagramas de bloques de las figuras 4.4(a), 4.4(b)
o 4.4(c) se puede representar como el diagrama de bloques de la figura 4.5(c)
con G1 (s) y G2 (s) dadas en (4.3) y (4.4).
Ejemplo 4.5 Considere el sistema en lazo cerrado mostrado en la figura
4.6(a). Este ejemplo representa un esquema de control conocido como contro-
lador con dos grados de libertad. Este es un sistema con dos entradas y una
salida. Por tanto, se puede proceder como en el ejemplo previo para, usando
el principio de superposici
on, escribir:
para obtener:
(Gc1 (s) + Gc2 (s))Gp (s)
C(s) = F (s) (4.6)
1 + (Gc1 (s) + Gc2 (s))Gp (s)
4.1 Diagramas de bloques 185
D(s)
R(s) + + C(s)
G c1 ( s)
+ + G p ( s)
G c2 ( s)
G c2 ( s)
G c1 ( s)
+
+
G c2 ( s)
(c) Caso cuando R(s) = 0
G c2(s)
R(s) +
V(s) C(s)
G c1(s) + G c2(s) G p ( s)
+
(b)
G c2(s)
G c1(s)+G c2(s)
R(s)
F(s) V(s) C(s)
G c1(s) + G c2(s) G p ( s)
+
+
1
(c)
se obtiene:
Gp (s)
C(s) = D(s)
1 + (Gc1 (s) + Gc2 (s))Gp (s)
es decir:
Gp (s)
G2 (s) = (4.8)
1 + (Gc1 (s) + Gc2 (s))Gp (s)
D(s) C(s)
+ G p ( s)
G c1(s) + G c2(s)
R(s) G 1(s)
+ C(s)
+
D(s) G 2(s)
p(s) = s2 + cs + d
188 4 Criterios de estabilidad y error
c2 > c2 4d > 0
De acuerdo a (4.9), esto significa que ambas races son reales, distintas y
negativas en este caso:
c + c2 4d
s1 = <0
2
c c2 4d
s2 = <0
2
Caso iii): c2 4d = 0.
En este caso ambas races son reales, iguales y negativas:
c
s1,2 = (4.10)
2
Ejemplo 4.6 N otese que el polinomio s2 + s + 1 satisface el caso i) por lo que
sus races son complejas conjugadas con parte real negativa. De hecho, no es
difcil verificar que estas races son 0.5 + j0.866 y 0.5 j0.866. Por otro
lado, el polinomio s2 + 2.5s + 1 satisface el caso ii) por lo que sus dos races
son reales, distintas y negativas. No es difcil verificar que dichas races son
2 y 0.5. Finalmente, el polinomio s2 + 2s + 1 satisface el caso iii) por lo
que sus dos races son reales, iguales y negativas. No es difcil verificar que
dichas races son 1 y 1.
4.2 Regla de los signos 189
p(s) = s2 + cs + d
c2 < c2 4d
Esto significa que ambas races son reales y distintas con una de ellas
positiva:
c + c2 4d
s1 = <0
2
c c2 4d
s2 = >0
2
Si d > 0 con c < 0 entonces:
c2 > c2 4d
190 4 Criterios de estabilidad y error
Como c < 0, ya que d > 0, entonces hay dos races reales positivas:
c + c2 4d
s1 = >0
2
c c2 4d
s2 = >0
2
Caso iii): c2 4d = 0.
En este caso d = c2 /4 > 0, por lo que c < 0 (vease 1)). Esto implica que
ambas races son reales, iguales y positivas:
c
s1,2 = >0 (4.12)
2
Ejemplo 4.7 A continuaci on se presentan varios polinomios de segundo or-
den y sus correspondientes races. Se deja como ejercicio para el lector que
verifique estos resultados, que identifique a que caso de los anteriores corres-
ponde cada uno y que compruebe que se cumplen las predicciones hechas en el
an
alisis arriba presentado.
s2 + 2s 1, races: 2.4142 y 0.4142.
s2 2.5s + 1, races: 2 y 0.5.
s2 s 1, races: 1.618 y 0.618.
s2 s + 1, races: 0.5 + j0.866 y 0.5 j0.866.
s2 2s + 1, races: 1 y 1.
En este caso es muy facil demostrar que se cumple algo similar a los poli-
nomios de segundo grado:
Criterio 4.3 Si los dos coeficientes tienen el mismo signo entonces la u nica
raz es real y negativa. Si un coeficiente tiene signo contrario al otro coeficiente
entonces la u nica raz es real y positiva.
Criterio 4.4 Si al menos un coeficiente tiene signo contrario a los otros coefi-
cientes entonces es seguro que existe al menos una raz con parte real positiva.
4.3 Criterio de estabilidad de Routh 191
Tabla 4.1. Tabla que se debe construir para aplicar el criterio de Routh.
sn an an2 an4 an6 ... 0
n1
s an1 an3 an5 an7 ... 0
sn2 b1 b2 b3 b4 ... 0
sn3 c1 c2 c3 c4 ... 0
sn4 d1 d2 d3 d4 ... 0
.. .. .. .. ..
. . . . . ... 0
s2 p1 p2 0
s1 q1 0
s0 p2
s2 + as + b
con a y b dos constantes reales. Dado que este polinomio es de segundo grado
(tiene dos races) entonces se puede usar la regla de los signos vista en la
secci
on 4.2 para concluir que las dos races tiene parte real negativa si
a > 0, b>0
a > 0, b>0
deben cumplirse simultaneamente para asegurar que ambas races tienen parte
real negativa. De este modo, el criterio de Routh y la regla de los signos de la
secci
on 4.2 dan el mismo resultado, como se esperaba.
s3 + 5s2 + 2s 8
Dado que existe un coeficiente de signo contrario al de los otros (8) se puede
usar la regla de los signos vista en la secci
on 4.2 para concluir que existe al
menos una raz con parte real positiva. Ahora se har
a uso del criterio de Routh
194 4 Criterios de estabilidad y error
para verificar este resultado. De acuerdo a las reglas del criterio de Routh se
construye la tabla 4.3. El criterio de Routh establece que el n
umero de cambios
de signo al recorrer los elementos de la primera columna es igual al n umero
de races con parte real positiva. Por tanto, se concluye que el polinomio s3 +
5s2 + 2s 8 tiene una raz con parte real positiva. De hecho, usando software
especializado se pueden usar metodos numericos para encontrar que:
s = 2, s = 1, s = 4
son las races del polinomio en cuesti
on.
s3 + as2 + bs + c
De acuerdo a las reglas del criterio de Routh se construye la tabla 4.6. Sin
embargo, la construcci
on de la tabla se detiene en el rengl
on correspondiente a
196 4 Criterios de estabilidad y error
s3 + 3s2 + s + 3
De acuerdo a las reglas del criterio de Routh se construye la tabla 4.8. Sin
embargo, la construccion de la tabla se detiene en el rengl on correspondien-
te a s1 porque todo ese renglon est
a formado exclusivamente por ceros. Este
es un caso especial que indica la existencia de races: 1) simetricas y reales,
2) simetricas e imaginarias, o 3) 4 races colocadas sobre los vertices de un
rect
angulo centrado en el origen. Para continuar con el metodo se debe for-
mar un polinomio con los elementos del rengl on inmediato superior al renglon
formado exclusivamente por ceros:
P (s) = 3s2 + 3,
se obtiene su derivada:
dP (s)
= 6s
ds
y se sustituyen estos coeficientes en el rengl
on formado exclusivamente por
ceros, como se muestra en la tabla 4.9, para continuar construyendo la tabla.
Esto significa que las races del polinomio 3s2 + 3 tambien son races del
polinomio s3 + 3s2 + s + 3. Estas races se pueden obtener como:
3s2 + 3 = 0, s = j
Por tanto, el polinomio s3 + 3s2 + s + 3 tiene una raz en s = j y otra en
s = j. De hecho, se puede usar software especializado para encontrar, usando
metodos numericos, que las races del polinomio s3 + 3s2 + s + 3 son:
s = 3, s = j, s = j
Ejemplo 4.15 (Caso especial. Un rengl on esta formado por ceros
exclusivamente o un rengl on est a formado por un s olo elemento el
cual, adem
as, es cero) Sea el siguiente polinomio:
s4 + s2 + 1
De acuerdo a las reglas del criterio de Routh se construye la tabla 4.10. Como
hay un renglon formado exclusivamente por ceros se procede como en el ejem-
plo anterior. El polinomio formado con los elementos del rengl on inmediato
superior (igual al polinomio original en este caso) es derivado respecto a s:
dP (s)
P (s) = s4 + s2 + 1, = 4s3 + 2s
ds
como se muestra en la tabla 4.11. Como hay dos cambios de signo al recorrer
la primera columna de la tabla 4.11 se concluye que existen dos races con
parte real positiva. Esto significa que es posible uno de los siguientes casos:
hay races simetricas y reales, o hay 4 races colocadas sobre los vertices de
un rectangulo centrado en el origen. De hecho, se puede hacer uso de softwa-
re especializado para encontrar, usando metodos numericos que las races del
polinomio s4 + s2 + 1 son:
Notese que hay dos races repetidas sobre el eje imaginario. De acuerdo a las
reglas del criterio de Routh se construye la tabla 4.12. Como hay un rengl on
formado exclusivamente por ceros se obtiene la derivada del siguiente polino-
mio:
dP (s)
P (s) = s4 + 2s2 + 1, = 4s3 + 4s
ds
y se contin
ua con el metodo como se muestra en la tabla 4.13.
el metodo indica, correctamente, que no hay races con parte real positiva.
Adem as, las races del polinomio P1 (s) = s2 + 1 son s = j, por lo que se
concluye estabilidad marginal. N otese que esto es incorrecto porque el poli-
nomio s4 + 2s2 + 1 tiene dos races imaginarias repetidas dos veces, lo cual
implica inestabilidad. As que debe tenerse cuidado con estos casos.
decir, cuando el tiempo es muy grande. Por esta razon, en lo que sigue se
estudia el comportamiento del error en estado estacionario.
Considere el sistema en lazo cerrado con realimentacion unitaria mostrado
en la figura 4.10. Como se muestra en esta seccion, un concepto fundamental
para la determinacion del error en estado estacionario es el tipo del sistema,
el cual se define a continuacion. Una funcion de transferencia de lazo abierto
de orden n siempre se puede escribir del siguiente modo:
k(s z1 )(s z2 ) (s zm )
G(s) =
sN (s p1 )(s p2 ) (s pnN )
donde n es el numero de polos de lazo abierto y m es el n umero de ceros de
lazo abierto con n > m, zj 6= 0, j = 1, . . . , m y pi 6= 0, i = 1, . . . , n N .
De este modo, N es el n umero de polos que la funcion de transferencia de
lazo abierto tiene en el origen (una vez cancelados los ceros que la funcion de
transferencia de lazo abierto pueda tener en el origen). Por definicion, el tipo
del sistema es igual a N , es decir, es el n
umero de integradores que el sistema
tiene en lazo abierto.
Por otro lado, la se
nal de error esta dada como la diferencia entre la salida
del sistema realimentado y la referencia o salida deseada:
r(t)
0 t
Figura 4.11. Entrada de prueba escal
on.
r(t)
0 t
Figura 4.12. Entrada de prueba rampa.
Entrada de prueba par abola. Suponga que ahora el avion pasa enfrente
del canon con aceleracion constante A (tratando de escapar), dirigiendose
hacia otro punto, y que inicialmente el ca non esta apuntando hacia el avion
(en cero). Entonces la referencia o la direccion en la que se desea apunte
el canon tiene la forma mostrada en la figura 4.13. Esta se nal se conoce
como parabola e indica que se desea que la direccion en la que apunta el
canon cambie con una aceleracion constante A (la aceleracion del avion).
Si r(t) es una parabola:
0, t<0
r(t) = 1 2
2 At , t0
2
donde A = d dtr(t)
2 es una constante que representa la aceleracion con que
r(t) cambia, entonces R(s) = L{r(t)} = sA3 es una funcion suficientemente
sencilla para ser manejada matematicamente.
Como cualquiera de las tres situaciones anteriores puede aparecer durante
el funcionamiento del sistema de control de posicion del ca non antiaereo, el
dise
no del controlador debe ser tal que el error de estado estacionario sea cero
o muy peque no ante cualquiera de las tres senales de prueba. A continuacion
se determinan las condiciones que permiten conseguir estas especificaciones
para cada senal de prueba.
204 4 Criterios de estabilidad y error
r(t)
1
2
At2
0 t
Figura 4.13. Entrada de prueba par
abola.
s A
ess = lm
s0 1 + G(s) s
A
= , kp = lm G(s)
1 + kp s0
por lo que:
por lo que:
4.4 Error en estado estacionario 205
s A
ess = lm
s0 1 + G(s) s2
A
= , kv = lm sG(s)
kv s0
por lo que:
Esto significa que el error en estado estacionario crece sin lmite (tiende a
infinito):
A
ess =
kv
Sistemas tipo igual a 1 (N = 1). En este caso la funcion de transferencia
de lazo abierto tiene la forma:
k(s z1 )(s z2 ) (s zm )
G(s) =
s(s p1 )(s p2 ) (s pn1 )
por lo que:
Suponga que R(s) = A/s3 es una parabola con segunda derivada igual a
A. Usando la expresion en (4.17) se obtiene:
s A
ess = lm
s0 1 + G(s) s3
A
= , ka = lm s2 G(s)
ka s0
r(t)
e ss6=0
e ss = 0 c(t)
e ss ! 1
0 t
Figura 4.14. Comportamiento de la salida c(t) cuando la referencia r(t) es una
rampa. N otese que ess = 0 cuando el tipo es mayor o igual a 2, ess 6= 0 cuando el
tipo es igual a 1 y ess cuando el tipo es igual a 0.
k
por la funci
on de transferencia Gm (s) = s+a , con k > 0 y a > 0, mientras que
la constante positiva kp representa la funci
on de transferencia de un controla-
dor proporcional. Esto significa que la corriente usada como se nal de control
se calcula como:
i = kp (d )
d ( s) + I (s) k
(s)
kp s ( s+a)
d ( s) + I (s) (s)
+ k
kp s ( s+a)
kv s
d ( s) + I (s) k
(s)
s+d
s+c s ( s+a)
d(s) + kp
s 2+k s+k i
k I (s) (s)
d d k
kd s s(s+a)
Xd (s) s+ b + + k X(s)
Ax s+ c
s ( s+ a )
s2
+ 1 2
k s A
Por otro lado, considere el sistema mostrado en la figura 4.22. Se trata del
sistema de control de posici on de un sistema ball and beam (vease el captulo
12). N otese que la funci on de transferencia de lazo abierto cuenta con dos
polos en s = 0. En este caso, estos dos polos de lazo abierto en el origen son
parte del modelo del sistema ball and beam a controlar, es decir, la planta posee
de manera natural dichos polos y no es necesario introducir un controlador
de tipo integral para generarlos. Por tanto, el error en estado estacionario
sera igual a cero si la posici
on deseada xd es un escalon o una rampa, pero si
la referencia es una par abola entonces el error en estado estacionario ser
a una
constante diferente de cero.
Ejemplo 4.21 Suponga que tiene un motor de CD con funci on de transferen-
(s) k
cia Gm (s) = I (s) = s(s+a) y que se desea que el error en estado estacionario
sea igual a cero cuando la referencia de posicion d sea una par abola. En este
caso se necesita que el sistema sea de tipo 3 y como la planta tiene un polo
4.5 Resumen del captulo 213
J + f + K = T
G1 (s)G2 (s)/m1
Xm2 (s) = F (s)
1 G1 (s)G2 (s)K/m1
F (t)
K
m1 m2
b1 b2
R(s) + C(s)
G(s)
H(s)
Suponiendo que los polos y los ceros de la planta se conocen, el metodo del
lugar de las races es una herramienta grafica que es util para determinar cuales
seran los polos de lazo cerrado si se usa un controlador que contiene ciertos
polos y ceros que se proponen o que se deben calcular. Por otro lado, tambien
es posible utilizar el lugar de las races para conseguir que uno o varios ceros de
lazo cerrado se cancelen con uno o varios polos de lazo cerrado. De acuerdo a la
seccion 3.9.1 esto permite reducir el orden del sistema y eliminar la presencia
de ceros de manera que la respuesta transitoria del sistema en lazo cerrado
puede ser aproximada por uno de los casos sencillos estudiados en las secciones
3.1 y 3.3. Estas caractersticas facilitan la tarea de dise nar el controlador de
manera que la funcion de transferencia de lazo cerrado M (s) tenga los polos
y los ceros que aseguren las caractersticas deseadas de respuesta transitoria.
El lugar de las races esta representado por varias curvas cuyos puntos cons-
tituyen los polos del sistema de lazo cerrado obtenidos conforme un parametro,
k, es variado desde k = 0 hasta k = +. De acuerdo a (5.1), todo punto s
que sea un polo de lazo cerrado (y que, por tanto, pertenezca al lugar de las
races) debe cumplir 1 + G(s)H(s) = 0. El lugar de las races se construye
proponiendo puntos s en el plano complejo s, los cuales deben ser verificados
si satisfacen la condicion para ser polos de lazo cerrado: 1 + G(s)H(s) = 0.
Para verificar esta condicion de una manera sencilla se propone un conjunto
de reglas que se listan a continuacion. Estas reglas dan lugar a la construccion
grafica del lugar de las races el cual se dibuja a partir de los polos y los ceros
de la funcion de transferencia de lazo abierto G(s)H(s). Esto significa que los
224 5 Dise
no usando la respuesta en el tiempo
polos de lazo cerrado seran determinados a partir de los polos y los ceros de
lazo abierto.
La funcion de transferencia de lazo abierto se puede escribir como:
Qm
j=1 (s zj )
G(s)H(s) = k Qn , n>m (5.2)
i=1 (s pi )
s zj = lj j
s pi = li i
Im (s)
s
s
zj
s
lj li
pi
i
j
zj pi Re (s)
Qm Qm
Xm Xn
j=1 l j j j=1 l j
G(s)H(s) = k Qn = k Qn j i (5.3)
i=1 li i i=1 li j=1 i=1
Por otro lado, los valores de s que satisfacen lo siguiente son los polos de lazo
cerrado:
1 + G(s)H(s) = 0 (5.4)
Esto significa que los polos de lazo cerrado son los valores de s que satisfacen:
si k = 0. Lo que significa que los polos de lazo cerrado (los que satisfacen
1 + G(s)H(s) = 0) son identicos a los polos de lazo abierto (s = pi ,
i = 1, . . . , n).
2. El lugar de las races termina (k ) en los ceros de lazo
abierto.
Partiendo de (5.8):
n
Y m
Y
1 + G(s)H(s) = (s pi ) + k (s zj ) = 0
i=1 j=1
m
Y
k (s zj ) = 0
j=1
Qn Qm
porque i=1 (s pi ) k j=1 (s zj ) si k . Lo que significa que los
polos de lazo cerrado (los que satisfacen 1 + G(s)H(s) = 0) son identicos
a los ceros de lazo abierto (s = zj , j = 1, . . . , m).
3. Cuando k hay nm ramas del lugar de las races que tienden
hacia alg un punto en el infinito del plano s. Esto significa que
la funcion de transferencia de lazo abierto G(s)H(s) tiene n m
ceros en el infinito.
Estas ramas pueden ser identificadas mediante la inclinacion de la linea
226 5 Dise
no usando la respuesta en el tiempo
asntota
Im (s) s
polos y ceros Re (s)
de lazo abierto
Figura 5.3. Polos y ceros de lazo abierto para determinar las asntotas de la regla
3.
punto s que se encuentra muy lejos del origen sobre una asntota, todos
los polos y todos los ceros de G(s)H(s), o de lazo abierto, se observan
todos reunidos en un punto del plano s, como se muestra en la figura 5.3.
Por tanto, el angulo con que contribuye cada polo y cada cero de lazo
abierto a la condicion en (5.7) (vease la figura 5.3) es identico al de los
demas, es decir, representando dicho angulo por la condicion de angulo
(5.7) se puede escribir como:
m
X n
X
j i = (m n) = (2q + 1)180 , q = 0, 1, 2, . . .
j=1 i=1
Im (s)
i 180 180
s
0 0
Re (s)
0i 0j
para alg
un q = 0, 1, 2, . . ., es decir (No. cerosDNo. polosD) es impar,
entonces tambien se cumple:
m
X n
X
j j = N (+180 ) = (2q + 1)180
j=1 i=1
(5.10)
s Im (s)
pv
pv
i j
pi zj
Re (s)
punto s es un punto que pertenece al lugar de las races y que esta infini-
tesimalmente cercano al polo en el punto pv . El angulo pv medido como
el angulo formado con respecto al eje horizontal positivo por el vector
trazado desde el polo en pv hacia el punto s es igual al angulo de salida
del polo complejo en pv . La condicion de angulo (5.7) establece que:
m
X n
X
j i =
j=1 i=1
5.1 Dise
no con el lugar de las races 229
= (z1 + z2 + + zm ) (p1 + p2 + + pv + + pn )
= (2q + 1)180 , q = 0, 1, 2, . . .
donde zj y pi representan los angulos debidos a los ceros y a los polos de
lazo abierto, respectivamente (puede considerarse que s = pv ). Despejando
pv :
pv = (2q + 1)180 + (z1 + z2 + + zm ) (p1 + p2 + + pn )
para q = 0, 1, 2, . . ., que equivale a la formula (5.10).
8. El angulo de llegada a un cero complejo se calcula como:
angulo de llegadaPa un cero complejo=
(2q + 1)180 [angulos de vectores desde los otros ceros hacia el cero
enP cuestion]
+ [angulos de vectores desde los polos hacia el cero en cuestion]
(5.11)
La explicacion de esta formula es como sigue. Considere la figura 5.6. El
s Im (s)
zv
zv
i j
pi zj
Re (s)
punto s es un punto que pertenece al lugar de las races y que esta infini-
tesimalmente cercano al cero en el punto zv . El angulo zv medido como
el angulo formado con respecto al eje horizontal positivo por el vector
trazado desde el cero en zv hacia el punto s es igual al angulo de llegada
al cero en zv . La condicion de angulo (5.7) establece que:
m
X n
X
j i =
j=1 i=1
= (z1 + z2 + + zv + + zm ) (p1 + p2 + + pn )
= (2q + 1)180
donde q = 0, 1, 2, . . ., mientras que zj y pi representan los angulos de-
bidos a los ceros y a los polos de lazo abierto, respectivamente (puede
considerarse que s = zv ). Despejando zv :
230 5 Dise
no usando la respuesta en el tiempo
10. Los puntos de cruce con el eje imaginario se pueden obtener usando el
criterio de Routh (vease la seccion 4.3).
(p1 + p2 ) = 180
s Im (s)
p2 p1
p2 p1 Re (s)
Figura 5.7. Un punto perteneciente al lugar de las races de un sistema con dos
polos en lazo abierto.
Im (s) s
p1
p2 p3
p2 p3 p1 Re (s)
(p1 + p2 + p3 ) = 180
Esto significa que la suma p1 + p2 debe ser menor en la figura 5.8 que
en la figura 5.7 lo cual se consigue si el punto s en la figura 5.7 se recorre
hacia la derecha como en la figura 5.8, es decir, si el lugar de las races es
doblado hacia la derecha en la figura 5.8. Es facil ver que esta desviacion
es mayor (el lugar de las races es empujado cada vez mas hacia la zona
de inestabilidad) conforme p3 es mayor, es decir, conforme p3 es mas
cercano al origen. Esto es una muestra de que un controlador integral
tiende a inestabilizar a un sistema en lazo cerrado.
12. Un cero de lazo abierto (en el semiplano izquierdo) agregado al
lugar de las races tiende a aumentar la estabilidad del sistema
en lazo cerrado. Este efecto estabilizante es mayor conforme el
cero se coloque cada vez m as cerca del origen.
Considere de nuevo la figura 5.7. De acuerdo a la condicion de angulo
(5.7):
(p1 + p2 ) = 180
z1 (p1 + p2 ) = 180
Esto significa que la suma p1 + p2 debe ser mayor en la figura 5.9 que
en la figura 5.7 lo cual se consigue si el punto s en la figura 5.7 se recorre
hacia la izquierda como en la figura 5.9, es decir, si el lugar de las races es
doblado hacia la izquierda en la figura 5.9. Es facil ver que esta desviacion
es mayor (el sistema es halado cada vez mas hacia la zona de mayor
estabilidad) conforme z1 es mayor, es decir, conforme z1 es mas cercano
232 5 Dise
no usando la respuesta en el tiempo
s Im (s)
p2 z1 p1
p2 z1 p1 Re (s)
Figura 5.9. El lugar de las races es halado hacia la izquierda al introducir un cero
adicional de lazo abierto.
5.2. Ejemplos
5.2.1. Control proporcional de posici
on
d(s) + I(s) k
(s)
kp s(s+a)
variada por el metodo para que tome valores desde 0 hasta +. Primero se
reescribe G(s)H(s) en la siguiente forma:
kp k
G(s)H(s) = (1 + 2 )
l1 l2
donde se han definido los vectores s 0 = l1 1 , s (a) = l2 2 (vease
la figura 5.11). Las dos condiciones fundamentales que definen el lugar de las
races son las condiciones de angulo (5.7) y de magnitud (5.6) las cuales se
expresan, respectivamente, como:
(1 + 2 ) = (2q + 1)180 , q = 0, 1, 2, . . .
kp k
=1
l 1 l2
La regla 5 indica que sobre el eje real solo existe lugar de las races entre los
puntos s = 0 y s = a. Ademas, de acuerdo a la regla 1, el lugar de las races
inicia (kp = 0) en s = 0 y s = a. Por otro lado, de acuerdo a las reglas 2 y
3, cuando kp tiende a + las ramas que inician en los puntos s = 0 y s = a
deben terminar en alg un cero de lazo abierto (no existe ninguno en este caso)
o en algun punto en el infinito del plano s. Por tanto, el lugar de las races
debe alejarse de los puntos s = 0 y s = a, sobre el eje real, conforme kp
crece de manera que debe existir un punto de separacion en alg un lugar entre
dichos puntos para luego dirigirse hacia el infinito sobre el plano s.
Por otro lado, de acuerdo a la condicion de angulo, (1 + 2 ) = 180 =
( + 1 ) en la figura 5.11, se concluye que = 2 y, por lo que los triangulos
t1 y t2 deben ser identicos para cualquier polo de lazo cerrado s. Entonces,
las dos ramas mencionadas anteriormente y que son mostradas en la figura
5.11 deben ser paralelas al eje imaginario. Esto significa que el punto donde
las dos ramas se separan del eje real esta ubicado en el punto medio entre los
polos ubicados en s = 0, s = a, es decir en (0 a)/2 = a/2. Esto tambien
234 5 Dise
no usando la respuesta en el tiempo
Im (s)
s
l2 l1
t1 t2
2 1
a a2 0 Re (s)
kp k
Figura 5.11. Lugar de las races para G(s)H(s) = s(s+a)
.
l1 l2 a2 a
kp = = , con l1 = l2 = ,
k 4k 2
los dos polos de lazo cerrado son reales, repetidos, negativos y ubicados en
s = a2 ; por lo que se alcanza la mayor rapidez sin que haya oscilaciones, iii)
a2
conforme kp > 4k se incrementa los dos polos se alejan del eje real (uno hacia
arriba y el otro hacia abajo) sobre la linea vertical que pasa por s = a2 ;
esto significa que el sistema es mas rapido (porque n , la distancia de los
polos al origen, aumenta; vease la seccion 3.3) y la respuesta del sistema
5.2 Ejemplos 235
en lazo cerrado (de segundo orden) es cada vez mas oscilatoria (porque el
angulo 90 disminuye y el amortiguamiento, dado como = sin(90 ),
disminuye; vease la seccion 3.3).
Todo lo anterior significa que no es posible conseguir simultaneamente una
respuesta tan rapida y tan amortiguada como se desee. Esto es una consecuen-
cia directa de que los polos de lazo cerrado no pueden ser ubicados en cualquier
lugar del plano s y solo pueden ser colocados sobre la linea gruesa mostrada
en la figura 5.11. En el siguiente ejemplo se muestra que la introduccion de
un cero en la funcion de transferencia de lazo abierto permite colocar los dos
polos de lazo cerrado en cualquier punto del plano s.
kd k l3
G(s)H(s) = 3 (1 + 2 ) (5.16)
l1 l2
donde se han definido los vectores s 0 = l1 1 , s (a) = l2 2 , s (c) =
l3 3 . Las condiciones de angulo y de magnitud se expresan, respectivamente,
como:
3 (1 + 2 ) = (2q + 1)180 , q = 0, 1, 2, . . .
kd k l3
=1
l1 l2
El lugar de las races correspondiente a este caso se puede obtener a par-
tir del obtenido para la funcion de transferencia en (5.13) considerando que
simplemente se ha adicionado un cero en s = c.
De acuerdo a la regla 12 el cero sera la causa para que las dos ramas del
lugar de las races mostrado en la figura 5.11 se doblen hacia la izquierda.
Esto puede comprobarse usando la regla 3 para encontrar que ahora existe una
sola rama del lugar de las races que tiende hacia una asntota que forma un
angulo de 1800 con el eje real positivo. Como se muestra en la figura 5.13 el
Im(s)
c a Re(s)
(a) c > a
Im(s)
Re(s)
a c
(b) c < a
kd k(s+c)
Figura 5.13. Lugar de las races para G(s)H(s) = s(s+a)
.
lugar de las races puede tener dos formas diferentes dependiendo de en donde
5.2 Ejemplos 237
Esto significa que usando kp y kd se pueden colocar los polos de lazo cerrado en
cualquier punto del semiplano complejo izquierdo s porque se puede asignar
cualquier valor a y a n . Por a otro lado, es com un ubicar los polos de
lazo cerrado (complejos conjugados) de manera que se consiga el tiempo de
subida tr y el sobre paso Mp ( %) deseados. Para esto, normalmente se usan
las expresiones en (3.59) para encontrar:
abs(ln(Mp ( %)/100))
= q (5.20)
2 + ln2 (Mp ( %)/100)
p !!
1 1 2
d = atan
tr
d
n = p (5.21)
1 2