4 - Fundamentos Matematicos para El Control Automatico de Procesos PDF
4 - Fundamentos Matematicos para El Control Automatico de Procesos PDF
4 - Fundamentos Matematicos para El Control Automatico de Procesos PDF
Fundamentos Matemáticos
para el Control Automático de
Procesos
Conceptos básicos de la
transformada de Laplace
Discusiones Preliminares
+ b
L t = e − st
tdt = lim dt
te − st
b →
0 0
1
L t = 2 para s 0
s
Señales de entrada
comunes en el
estudio de respuestas
de sistemas de
control.
(a) Función escalón
unitario: u(t).
(b) Pulso.
(c) Función impulso
Figura # 1: Las cuatro señales
unitario: presentadas en la figura se aplican
(t ) comúnmente como entradas de pro
(d) Onda sinusoidal, cesos e instrumentos para estudiar
su respuesta dinámica. Ahora se
2 usará la definición de la
sen t = transformada de Laplace para
T
deducir sus transformadas.
a) Función escalón unitario
• Éste es un cambio súbito de magnitud
unitaria, como se ilustra gráficamente en
la figura # 1 acápite (a). Su representación
algebraica es:
0 t0
u (t ) = ( 4 − 2)
1 t0
• Sustituyendo en la definición de la
transformada:
1 − st 1
L u ( t ) = u ( t ) e dt = e
− st
0 = − ( 0 − 1)
0 S S
1
L u ( t ) = ( 4 − 3)
S
b) Pulso de magnitud H y duración T
• El pulso se muestra gráficamente en la figura
# 1 y se representa como:
0 t 0, t T
f (t ) = ( 4 − 4)
H 0 1 T
• Sustituyendo en la ecuación se obtiene:
H − st T
L f ( t ) = 0 f ( t )e dt = 0 He dt = e 0
− st − st
s
1 − e − st
( t ) = − ( e − 1) = H
H − st
L f ( 4 − 5)
S s
c) Función de impulso unitario
• Esta función, también conocida como función
delta Dirac y representada por el símbolo δ(t),
se grafica en la figura # 1 con el acápite c.
• Es un pulso ideal de duración cero y área
unitaria.
• Toda su área se concentra en el tiempo cero.
Como su función es cero en todo momento
excepto en cero, y el término e-st en la
ecuación es igual a la unidad en t=0, la
transformada de Laplace es
c) Función de impulso unitario
L ( t ) = ( t )e dt = 1
− st
( 4 − 6)
0
L af ( t ) = aL f ( t ) = aF ( s ) ( 4 − 9)
• La propiedad distributiva de la adición también
se sigue de la propiedad de la linealidad:
L af ( t ) + bg ( t ) = aL f ( t ) + bL g ( t )
= aF ( s ) + bG ( s ) ( 4 − 10 )
• donde a y b son constantes.
• Es sencillo deducir ambas fórmulas
mediante la aplicación del la definición de la
transformada de Laplace.
Teorema de diferenciación real
= 0 − f ( 0 ) + s f ( t ) e − st dt
0
= sF ( s ) − f ( 0 )
• La extensión a derivadas de orden superior
es directa
d 2 f (t ) d df ( t ) df ( t ) sf
L 2 = L = sL − t =0
dt
dt dt dt dt
df df
= s sF ( s ) − f ( 0 ) − t = 0 = s F ( s ) − sf ( 0 ) −
2
t =0
dt dt
• En general
d n f (t ) d n −1
f
= s F ( s ) − s f ( 0 ) − ... − n −1 ( 4.12 )
n n −1
L n t =0
dt dt
0 s
• La demostración de este teorema se hace
integrando por partes la definición de la
transformada de Laplace. La transformada
de Laplace de la integral n-ésima de una
función es la transformada de la función
dividida entre s.
Teorema de traslación real
• Este teorema se refiere a la traslación de
una función en el eje del tiempo.
• Cómo se ilustra en la figura 2. La función
trasladada es la función original con retraso
en el tiempo.
• Como se verá mas adelante, los retrasos de
tiempo son causados por retrasos de
transporte, un fenómeno también conocido
como tiempo muerto. El teorema establece
que
L f ( t − t0 ) = e − st0
F (s) ( 4 − 15)
• Como la transformada de Laplace no
contiene información acerca de la función
original para tiempo negativo, la función
retrasada debe ser cero para todos los
tiempos menores al tiempo de retraso
(véase la figura 2).
• Esta condición se satisface si las variables
del proceso se expresan como desviaciones
de la condición inicial de estado
estacionario.
Demostración
• Por la definición de la transformada de Laplace
Sea = t −t0 ( 0 t = t−0st +t + ) y sustituyendo:
L f ( t − t0 ) = f ( ) e (0 )
d ( t0 − )
r =− t0
= f ( ) e − st0 − st
e dt = e − st0
0 f ( ) e dt
− sr
r =0
= e − st0 F ( s ) q.e. d
En esta demostración se utilizó el hecho de que :
f ( ) = 0 para 0 ( t t0 ) .
Figura # 2: La función retrasada
es cero para todos los tiempos
menores que el tiempo de
retardo t0.
Teorema del valor final
• Este teorema permite calcular el valor final
o de estado estacionario de una función a
partir de su transformada.
• También es útil para verificar la validez de
la transformada que se deduce.
• Si el límite de f(t) cuando t→∞ existe, se
puede calcular a partir de la transformada
de Laplace como sigue:
lim f ( t ) = lim sF ( s ) ( 4 − 16 )
t → s →
Teorema de la diferenciación compleja
d
L tf ( t ) = − F ( s ) ( 4 − 17 )
ds
Teorema de la traslación compleja
L e f ( t ) = F ( s − a )
at
( 4 − 18)
Teorema del valor inicial
• Este teorema permite calcular el valor
inicial de una función a partir de su
transformada.
• Constituiría otra demostración de la validez
de las transformadas deducidas si no fuera
por el hecho de que, en el análisis de
dinámica de procesos, las condiciones inicial
es de las variables por lo general son cero.
El teorema establece que
lim f ( t ) = lim sF ( s ) ( 4 − 19 )
t → s →
Ejemplo # 3
• Obtener la transformada de Laplace de la
siguiente ecuación diferencial
d y (t )
2
dy ( t )
9 2
+6 + y (t ) = 2x (t )
dt dt
• Con condiciones iniciales cero en estado
estacionario, es decir,
dy
y (0) = 0 =0
dt
Solución
• Aplicando la propiedad de linealidad, se toma
la transformada de Laplace de cada término:
d 2 y (t ) dy ( t )
9L 2 + 6L + L y ( t ) = 2L x ( t )
dt dt
9s Y ( s ) + 6sY ( s ) = 2 X ( s )
2
• Finalmente, despejando Y(s), se tiene
2
Y (s) = 2 X (s)
9s + 6s + 1
• El ejemplo muestra cómo la transformada de
Laplace convierte la ecuación diferencial
original en una ecuación algebraica, la cual
puede reordenarse para despejar la variable
dependiente Y(s).
Ejemplo # 4
• Obtener la transformada de Laplace en la
siguiente función:
e ( t ) = u ( t − 3) 1 − e −( t −3) 4
• Nota: el término u(t-3) en esta expresión
indica que la función es cero para t<3.
• Recuérdese de la figura # 1 , que u(t-3) es
un cambio de cero a uno en t=3, lo cual
significa que la expresión entre corchetes e
stá multiplicada por cero hasta que t=3 y
que después está multiplicada por la unidad.
Solución
f ( t ) = u ( t ) 1 − e −1 4
= u ( t ) − u ( t ) e t 4
• Aplicar, la propiedad de linealidad, y utilizar
la tabla 4-1, con a=1/4
1 1 2
F (s) = − =
s s + 1 s ( 4s + 1)
4
• A continuación, aplicar el teorema de
traslación real
−3 s
e
C ( s ) = L f ( t − 3) = e F ( s ) =
3s
s ( 4s + 1)
Solución de ecuaciones
diferenciales mediante la
Transformada de Laplace
• En este tema se presenta el uso de la
transformada de Laplace para resolver las
ecuaciones diferenciales que representan la
dinámica de los procesos y sus sistemas de
control.
• Debido a que el objetivo es encontrar la forma
en que las señales de salida responden a las
funciones de forzamiento de entrada,
asumimos siempre que las condiciones iniciales
están en estado estacionario (derivadas con
respecto al tiempo cero).
• También definimos todas las variables como
desviaciones de sus valores iniciales.
• Esto hace que todas las variables de
desviación tengan un valor inicial igual a
cero.
Procedimiento de solución por la
Transformada de Laplace
• El procedimiento para resolver una ecuación
diferencial utilizando la transformada de
Laplace consta de tres pasos:
1. Transformar la ecuación diferencial en una
ecuación algebraica en función de la
variable s de la transformada de Laplace
2. Despejar La transformada de la variable de
salida (o dependiente)
3. Invertir la transformada para obtener la
respuesta de la variable de salida en
función del tiempo, t.
• Considérese la siguiente ecuación diferencial
de segundo orden:
d 2 y (t ) dy ( t )
a2 2
+ a1 + a0 y ( t ) = bx ( t ) ( 4 − 19 )
dt dt
dy ( t )
L = sY ( s ) − y ( 0 )
dt
• A continuación se sustituyen estos términos
en la ecuación 4-20 y se reordenan para
obtener
( a2 s + a1s + a0 ) Y ( s ) − ( a2 s + a1 ) y ( 0) − a2 dt = bX ( s )
2 dy
t =0
A1 = lim ( s − r1 ) Y ( s )
2
s → r1
1 d
A2 = lim ( s − r1 ) Y ( s ) ( 4 − 30 )
2
s → r1 1! ds
• Nuevamente, la inversión de la ecuación
4-29 se lleva a cabo hallando la
correspondencia entre sus términos y los de
la tabla 4-1.
• El primer término concuerda con la sexta
entrada de la tabla para a=-r1 obteniéndose
la inversa.
s → r1
1 d k −1
Ak = lim ( s − r1 ) Y ( s ) ( 4 − 32 )
m
s → r1 ( k − 1) ! ds k −1
• Para k = 2,………m. Por lo tanto, la
función inversa es:
A1t m −1 A2t m − 2 r1t
y (t ) + + ... + Am e + ... ( 4 − 33)
( m − 1) ! ( m − 2 )!
• El siguiente ejemplo está diseñado para
ilustrar numéricamente el procedimiento de
expansión en fracciones parciales y el
proceso completo de inversión.
• Se consideran tres casos: raíces reales no
repetidas. raíces repetidas y raíces
conjugadas complejas.
Ejemplo 2
• Dada la ecuación diferencial cuadrática
considerada en la discusión anterior, ecuación
4-19, con condiciones iniciales en estado
estacionario iguales a cero, se obtendrá la
respuesta escalón unitario de la variable de
salida y(t) para tres conjuntos diferentes de
parámetros.
• a) Raíces reales no repetidas. Sean: a2 = 9,
a1 = 10, a0 = 1 y b = 2, en la ecuación 4-19.
Entonces, la respuesta escalón unitario
obtenida con la ecuación 4-24 es:
• Entonces, la respuesta escalón unitario
obtenida con la ecuación 4-24 es:
2 1
y (t ) = L 2
−1
9 x + 6s + 1 s
• Las raíces, obtenidas con la fórmula
cuadrática. son r1=-1/9, r2=-1. El
denominador se factoriza como sigue:
2 1 A1 A2 A3
Y (s) = = + +
1 s s + 1 s +1 s
9 s + ( s + 1)
9 9
• Los coeficientes se calculan utilizando la
ecuación 4-28: A = lim 2 = 2.25
s →1 9 9 ( s + 1) s
2
A2 = lim = 0.25
s →−1 1
9 s + s
9
2
A3 = lim =2
s →0 1
9 s + ( s + 1)
9
2 A1 A2 A3
Y (s) = 2
= 2
+ +
1 1 1 s
9 s + s s + s+
3 3 3
• Los coeficientes que se obtienen. con la
ecuación 4-32, son:
2 2
A1 = lim =−
s →1 3 9 s 3
1 d 2 2
A2 = lim = lim − 2 = −2
s →1 3 1! ds 9 s
s →1 3 9s
• Y A3 = 2, como antes. La respuesta
escalón se obtiene por correspondencia con
las entradas de tabla 4-1.
2 −t 3
y ( t ) = − t − 2 e + 2u ( t )
3
Par de raíces conjugada complejas
• Sea a1=3 y sean los otros parámetros como
antes. Las raíces obtenidas con la fórmula
para las ecuaciones cuadráticas.
• son :
−0,167 i 0, 289; donde i = −1
es la unidad de los numeros inmaginarios.
• La transformada de la salida es
2
Y (s) =
9 ( s + 0.167 − i 0.289 )( s + 0.167 + i 0.289 ) s
A1 A2 A3
Y (s) = + +
s + 0.167 − i 0.289 s + 0.167 + i 0.289 s
• Una vez más los coeficientes se calculan
utilizando la ecuación 4-28:
2
A1 = lim = −1 + i 0577
s →0.167 + i 0.289 9 ( s + 0.167 + i 0.289 ) s
2
A2 = lim = −1 − i 0577
s →0.167 −i 0.289 9 ( s + 0.167 − i 0.289 ) s