transformadas de laplace
Post on 02-Jun-2015
25.230 Views
Preview:
DESCRIPTION
TRANSCRIPT
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
1I N G E N I E R Í A E N S I S T E M A S - UNCAus
MATEMATICA SUPERIOR – RESUMEN UNIDAD Nº3
Alumno: Nuñez, Juan
Docentes:
Prof. Marina BloeckLic. Nori Cheein de Auat
INGENIERIA EN SISTEMAS
AÑO 2011AÑO 2011
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
Transformada de Laplace
Para comenzar con la exposición primero vamos a demostrar de donde viene la transformada de Laplace, para ello vamos a tener en cuenta la siguiente fórmula matemática, que es la trasformada de Laplace:
F ( s )=∫0
∞
f ( t ) e−st dt(1)
La función F de la variable independiente s se la llama transformada de Laplace de la función original f de la variable independiente t. Para demostrar (1) tenemos que tener en cuenta SERIES DE POTENCIAS pues la transformada de Laplace es un caso especial de SERIES DE POTENCIAS, esta serie tiene la forma:
∑n=0
∞
an xn(2)
El resultado de la sumatoria de (2), si converge, converge a una función de x:
∑n=0
∞
an xn=A (x )
anes una función de la variable independiente n, entonces podemos expresar esto de la siguiente manera funcional:
f :→ (3)
Donde nosotros tomamos un elemento del conjunto de los naturales y atravez de una función llegamos a:
f :❑→f f (n )=a(n)
Ahora vamos a expresar (2) de manera funcional
∑n=0
∞
a(n) xn=A (x )
Esta sucesión es una sucesión discreta ¿Por qué es discreta?, pues simple porque el dominio de lo que estamos siguiendo hasta ahora va desde los naturales a los reales (3), esto significa que entre los valores 1 y 2 no hay ningún valor.
Análogamente trabajaremos la SERIE DE POTENCIAS pero en un dominio continuo y en este caso entre los valores 1 y 2 tenemos muchos valores, redefinimos la manera funcional de (2) cambiando la variable n (discreta) por una variable t (continua) y t está definida en un intervalo:
0< t<∞
2I N G E N I E R Í A E N S I S T E M A S - UNCAus
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
Como ahora la serie de potencias está definida en dominio continuo en lugar de trabajar una sumatoria vamos a trabajar con integrales llamando a (n ) como f (t)
∫0
∞
f ( t ) x tdt=F (x) (4)
(4) es lo que conocemos como integrales impropias. En esta integral hay ciertas cosas que tenemos que cambiar, y mirando (1) nos damos una idea de los cambios que tenemos que realizar, claramente tenemos que trabajar x t
En la integral impropia (4) “x” es una constante entonces:
x=e ln x
e y ln xson funciones inversas entonces al operar regresamos a la x. Pero como dijimos nosotros queremos sustituir x t entonces:
x t=[e ln x ]t
Aparte de esta transformación hay otras cosas que tenemos que tener en cuenta, si x es un valor que es mayor que 1 cuando evaluemos la integral vamos a notar que x t va a crecer
demasiado rápido con respecto a f ( t ) porque tiene un comportamiento exponencial. La única chance que tenemos para que (4) converja es que el valor de x sea menor a 1 y mayor que 0:
0<x<1
Con este intervalo definido para x tenemos que el ln x :
ln 0=−∞ ; ln 1=0
Entonces:
ln x<0
Como la variable x del logaritmo natural es continua vamos a definir una nueva variable que sea igual al ln x:
−S=ln x con S>0
Para hacer ver que el resultado es negativo escribimos –S. Ya tenemos listo todo, ahora replanteamos (4) con las nuevas correspondencias:
∫0
∞
f (t ) x tdt=∫0
∞
f ( t ) [e ln x ]tdt=∫0
∞
f (t ) [ e−S ]tdt=∫0
∞
f ( t ) e−St dt
3I N G E N I E R Í A E N S I S T E M A S - UNCAus
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
∫0
∞
f ( t ) e−St dt=F (s)
Y llegamos al resultado que esperábamos. Observemos algo MUY IMPORTANTE de las transformadas de Laplace, nosotros dentro de la integral tenemos f ( t ) (una función de la variable independiente t) y nos devuelve F ( s )(una función de la variable independiente s).
Nosotros sabemos que:
Con todo esto visto estamos en condiciones de de presentar la definición formal de la transformada de Laplace.
Sea f de la variable independiente t, para toda t ≥0. Multiplicando dicha función por e−St e integrando con respecto a t desde cero hasta el infinito; entonces, si la integral resultante existe, es una función de s, por ejemplo F(s):
F ( s )=∫0
∞
f ( t ) e−st dt
La función F de la variable independiente s, se llama transformada de Laplace de la función f de la variable independiente t y se denota por:
F ( s )=£ ( f )=∫0
∞
f ( t ) e−st dt
4I N G E N I E R Í A E N S I S T E M A S - UNCAus
f(x) OPERADOR
h(x)
Transformada f(t) F(s)
DERIVADA f(x) 𝑓′(𝑥)
La transformada de Laplace es un método para resolver ecuaciones diferenciales. El proceso de resolución consta de traspasos principales:
Paso 1: El problema “duro” dado se transforma en una ecuación “sencilla” que la llamaremos ecuación subsidiaria
Paso 2: La ecuación subsidiaria se resuelve a través de manipulaciones algebraicas.
Paso 3: La solución de la ecuación subsidiaria se transforma en sentido inverso, a fin
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
Algunas propiedades
La transformada de Laplace, como trabaja con integrales hereda algunas propiedades de ella:
1) Suma de transformadas de Laplace (s1)
F ( s )+G ( s)=£ (f +g )=£ ( f )+£ (g )=∫0
∞
f ( t ) e−st dt+∫0
∞
g (t ) e−st dt
2) Trasformada de Laplace de una constante por una función (s2)
cF (s )=£ (cf (t))=c £ ( f (t ))=c∫0
∞
f (t )e− st dt
Linealidad de la transformada de Laplace
De (s1) y (s2) desprende el teorema de linealidad de transformadas de Laplace, lo vemos:
£ (af (t )+bg( t))=∫0
∞
e−st [af (t )+bg (t)] dt=
∫0
∞
e−st [af ( t ) ] dt+∫0
∞
e−st [bg (t ) ]dt=a∫0
∞
e−st [ f ( t ) ]dt+b∫0
∞
e−st [ g ( t ) ]dt
Entonces:
£ (af ( t )+bg( t))=a∫0
∞
e−st [ f ( t ) ]dt+b∫0
∞
e−st [ g ( t ) ]dt
Existencia de la transformada de Laplace
Ahora veremos que no es necesario que la función sea continua. Este hecho tiene importancia practica, ya que son precisamente las entradas (fuerzas impulsoras) discontinuas para las que las que la transformación de Laplace se vuele
5I N G E N I E R Í A E N S I S T E M A S - UNCAus
La transformada de Laplace es un método para resolver ecuaciones diferenciales. El proceso de resolución consta de traspasos principales:
Paso 1: El problema “duro” dado se transforma en una ecuación “sencilla” que la llamaremos ecuación subsidiaria
Paso 2: La ecuación subsidiaria se resuelve a través de manipulaciones algebraicas.
Paso 3: La solución de la ecuación subsidiaria se transforma en sentido inverso, a fin
t
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
particularmente útil. Basta con requerir que la función de la variable independiente t sea seccionalmente continua sobre todo un intervalo finito para t ≥0 , entonces las únicas discontinuidades que puede tener una función seccionalmente continua son los saltos finitos; estas se conocen como discontinuidades ordinarias:
LA DEMOSTRACION DE ESTE TEOREMA LA VEREMOS MAS ADELANTE, YA QUE NECESITAMOS SABER EL PRIMER TEOREMA DE TRASLACION Y ALGUNAS PROPIEDADES MAS.
Primer teorema de traslación
Primero vamos a calcular la transformada de Laplace de 1:
£ ( f (1))=∫0
∞
1e−St dt=limb→∞
∫0
b
1e−St dt=¿ limb→∞ [−1
Se−St ]
0
b
¿=limb→∞
1
S(−e−Sb+1)
Entonces:
£ (1 )= 1S (5)
Para que esto converja condicionamos a “S” con: S>0
Podemos generalizar (5) utilizando (s2):
£ (cf (1))=c∫0
∞
1e−St dt= limb→∞
c∫0
b
1e−St dt=¿ limb→∞
c [−1Se−St ]
0
b
¿=limb→∞
c1S(−e− Sb+1)
6I N G E N I E R Í A E N S I S T E M A S - UNCAus
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
Entonces:
£ (cf (1))= cS (6)
Para que esto converja condicionamos a “S” con: S>0
Prestar especial atención aquí porque voy a calcular la transformada de Laplace de:
f (t)eat
£ ( f (t)eat )=∫0
∞
f ( t ) eat e−St dt=¿∫0
∞
f ( t ) e−( s−a) tdt=¿∫0
∞
f ( t ) e−( s−a) tdt ¿¿
Hacemos un cambio de variable y replanteando:
w=s−a
∫0
∞
f ( t ) e−wt dt(7)
Si miramos con buen ojo a (7), veremos que es la transformada de Laplace:
∫0
∞
f ( t ) e−wt dt ≡∫0
∞
f ( t ) e−st dt
Solo que:
∫0
∞
f ( t ) e−wt dt=F(w) y ∫0
∞
f (t ) e−st dt=F (s)
Pero como w=s−a:
∫0
∞
f (t ) e−(s−a)tdt=F (s−a)
Y por ultimo concluimos que:
£ ( f (t)eat )=F(s−a)
Es lo que se conoce como teorema de translación, y nos dice que cuando se tiene una función de la variable independiente t y multiplicado esta por eat provocamos un corrimiento de la función en (s-a). Para entender esto imaginamos una función de la variable independiente t igual al tiempo (el dominio de esa función es el tiempo), la multiplicamos por eat y esto va a provocar un corrimiento en la línea del tiempo hacia “a” (entonces tenemos un corrimiento provocado poreat ).
7I N G E N I E R Í A E N S I S T E M A S - UNCAus
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
Ahora vamos a calcular la transformada de Laplace de eat solo que la función de la variable independiente igual a 1:
£ ( f (1)eat )
Sabemos que:
£ ( f (1))=1S
Aplicando el teorema de traslación (porque a la función la multiplicamos por eat y como sabemos esto provoca el corrimiento de la función hacia a):
£ ( f (1)eat )= 1s−a
El cual sería el resultado de la transformada de Laplace de eat
Teorema de la existencia para las transformadas de Laplace
Sea la función f de la variable independiente t seccionalmente continua sobre todo el intervalo del semieje t>0 y que satisface:
|f (t)|≤M eat
Para toda t ≥0 y para algunas constantes a y M; entonces la transformada de Laplace de f (t) existe para toda s>a. Como la función f de la variable independiente t seccionalmente continua; e−St f (t) es integrable sobre cualquier intervalo finito del eje t. Suponiendo que s>a se obtiene:
8I N G E N I E R Í A E N S I S T E M A S - UNCAus
Teoremas y propiedades que necesitaremos para demostrar el teorema de existencia.
£ ( f (1)) 1S
£ ( f (1)eat ) 1s−a
£ (cf (1)) cS
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
|£ (f ( t ))|=|∫0
∞
f (t)e−St dt|≤∫0
∞
|f (t )|e−St dt ≤∫0
∞
Meat e−St dt
Vemos que a la integral ∫0
∞
M eat e−St dt la podemos expresar como:
£ (M f (1)eat )=∫0
∞
M eat e−St dt
Aplicando las transformadas y teoremas que hemos visto anteriormente (transformada de una constante multiplicado por f(1) y el primer teorema de traslación) tenemos que:
£ (M f (1)eat )=∫0
∞
M eat e−St dt=M∫0
∞
eat e−St dt=¿=M∫0
∞
e−(S−a)tdt=¿ MS−a
¿¿
Entonces:
|£ (f ( t ))|=|∫0
∞
f (t)e−St dt|≤∫0
∞
|f (t )|e−St dt ≤∫0
∞
Meat e−St dt= MS−a
Que era lo que queríamos demostrar. En donde se necesito la condición de que S>a
Transformada de Laplace del seno y coseno
Vamos a ver una manera diferente de calcular la transformada de Laplace del cos(at) y el sen(at).No vamos a utilizar la definición de la transformada de Laplace para calcular eso, para ello vamos a recurrir a la formula de euler:
cos (at )= eiat+e−iat
2=1
2(e iat+e−iat)
£ (cos (at))=12£ (e iat+e−iat )=1
2¿
Por el teorema de traslación sabemos cuál es la transformada de eat, entonces:
£ (cos (at))=12 [ 1s−ia
+ 1s+ia ]
Calculando el común denominador (multiplicación del complejo por su conjugado):
9I N G E N I E R Í A E N S I S T E M A S - UNCAus
CALCULO AUXILIAR
u dv
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
£ (cos (at))=12 [ (s+ia )+(s−ia)
(s−ia )(s+ia) ]=12 [ 2 ss2+a2 ]=¿
Finalmente:
£ (cos (at))=¿
De manera análoga por formulas de Euler procedemos a calcular el sen(at):
£ (sin (at))= 12i [ 1
s−ia− 1s+ia ]
Calculando:
£ (sin (at))= 12i [ ( s+ia )−(s−ia)
( s−ia )(s+ia) ]= 12 i [ 2ia
s2+a2 ]=¿
Finalmente:
£ (sin (at))=¿
Transformada de Laplace de t n
Todavía no hemos trabajado transformadas de Laplace de polinomios, por ello vamos a empezar calculando t n haciéndolo por definición:
£ (t n )=∫0
∞
t n e−St dt
Para resolver esta integral vamos a utilizar el método de integración por partes, realizamos cálculos auxiliares:
10I N G E N I E R Í A E N S I S T E M A S - UNCAus
CALCULO AUXILIAR
u=tn ;du=n tn−1
∫ d v=∫e−S t dt
v=−1se−St
P(t)
Q(t)
Indeterminación de tipo:
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
Y recordando el método de integración por partes:
£ (t n )= −t n e−St
S |0
b
b→∞
lim ¿−∫0
∞−1se−St nt n−1¿
dt = −t ne−St
S |0
b
b→∞
lim ¿(i1)+nS∫0∞
e−St tn−1¿dt (i2)
Vamos a analizar (i1):
Para el extreme inferior = 0
−0n e−S0
S=−0S
=0
Para el extreme superior = ∞
−t n e−St
St →∞
lim ¿=∞ne−S ∞
S=∞0S
¿
Pero como ∞ 0 es indeterminado vamos a trabajar un poco más:
Los posibles casos son:
Resultado Condición∞ Si el grado de P(t) es mayor que el grado de Q(t)
11I N G E N I E R Í A E N S I S T E M A S - UNCAus
Regla mnemotécnica: "Un Día Vi Una Vaca sin rabo (menos integral) Vestida
De Uniforme".
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
0 Si el grado de Q(t) es mayor que el grado de P(t)ab
Si el grado de P(t) es igual que el grado de Q(t), siendo a y b los coeficientes de los términos de mayor grado de P(t) y Q(t)
En nuestro caso el grado de Q(t) es mayor que el grado de P(t), porque sabemos que la exponencial (Q(t)) crece más rápido que t n (P(t)), entonces:
−t n
S eStt→∞
lim ¿=0¿
Con esto queda claro que (i1) se hace cero cuando lo evaluamos en infinito y en cero, ahora solo nos queda trabajar con (i2):
nS∫0
∞
e−St t n−1dt
Esto no es más que la transformada de Laplace de t n−1; nS £ (t n−1 )
Replanteando tenemos que:
£ (t n )= nS£ ( tn−1 )
nS£ (t n−1 )=n(n−1)
s2 £ (t n−2 )
n(n−1)s2 £ (t n−2 )=n(n−1)(n−2)
s3 £ (t n−3 )=⋯=n!sn£ (t 0 )
Pero como t 0=1 :
n !
sn£ (1 )
Y sabemos que:
£ (1 )= 1S
A partir de las transformadas demostradas formamos la siguiente tabla:
12I N G E N I E R Í A E N S I S T E M A S - UNCAus
£ (t n )=n!sn£ (1 )= n !
sn+1
1
1
f(t)
t
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
f(t) F(s)£ (1 ) 1
S£ (c ) c
S£ (cos (at)) ¿£ (sin (at)) ¿
£ (t n ) n !
sn+1
Obtención de la transformada de Laplace a partir de una grafica de función
La función que vemos en la grafica está definida por trozos, en el intervalo para t de 0 a 1 tenemos una función (todavía no sabemos cuál) y cuando va de 1 al infinito dicha función vale 0
Para calcular cual es la función en el intervalo para t de 0 a 1 suponemos que:
f ( t )= y ; x=t
Para encontrar la ecuación de la recta hacemos:
x+ y=1
Para: x=0
y=1
Para: y=0
x=1→x+ y=1verifica
13I N G E N I E R Í A E N S I S T E M A S - UNCAus
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
Entonces:
f ( t )=1−t
Ahora trabajamos con la definición de la transformada de Laplace y como en nuestra grafica tenemos una función dividida en trozos tenemos que:
∫0
1
f ( t ) e−St+¿∫1
∞
f ( t ) e−St=¿∫0
1
(1−t)e−St+¿∫1
∞
0e−St ¿¿¿
Solo nos resta trabajar con: ∫0
1
(1−t )e−St ; integramos por partes
£ (1−t )=(1−t)e−St
S |0
1
−∫0
1−1se−St−dt=¿=
(1−t)e−St
S |0
1−1s∫0
1
e−St dt=1S−1s∫0
1
e−St dt= 1S+ 1S2 e
−St|0
1
=1S+¿ 1S2 [e−S−1 ]=1
S+ e
−S
S2 − 1S2 ¿¿
Entonces:
F (1−t )=1S+ e
−S
S2 − 1S2
Segundo teorema de traslación
En síntesis, recordando el primer teorema de traslación nos decía que el corrimiento de la función era sobre la variable s. Lo que nos dice este segundo teorema es que al multiplicar la función por eat vamos a generar un corrimiento, pero ese corrimiento es especial porque la función de la variable independiente t solo está definida para valores de t positivos (t>0), entonces el corrimiento tiene que ser solo para valores positivos. Para demostrar este corrimiento necesitamos hablar de una función que se llama ESCALON UNITARIA:
14I N G E N I E R Í A E N S I S T E M A S - UNCAus
CALCULO AUXILIAR
u=1−t ;du=−dt
∫ d v=∫e−St dt
v=−1se−St
1
a
f(t)
t
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
U (t−a )={0 parat<a1 parat ≥a
La función escalón unitaria (U ( t−a )) está definida por tramos, y esto nos indica que para cualquier valor t<a de t la función esta con el valor 0, cuando esta t ≥ a dicha función toma el valor 1.
Una explicación un poco mas practica seria decir que la función para valores de t<a esta apagada (con valor 0, color azul en la grafica), cuando los valores t ≥ a la función se enciende (toma el valor 1, color rojo en la grafica).
Con la función escalón unitario podemos escribir como se comportan funciones definidas por pedazos, por ejemplo:
f (t )={ 0; para t<0g ( t ); para0≤t<ah ( t ) ; para a≤t<b
0 ; para t ≥b
Esto nos indica que nosotros necesitamos que g(t) funcione en el valor 0 de la variable independiente t (es decir que g(t) tome el valor 1) para esto hacemos que el escalón unitario se levante a 1 en 0 de la siguiente manera:
f ( t )=U ( t )g (t)
Ahora cuando a≤ t<b necesitamos que se active la h(t) y g(t) se desactive, entonces hacemos:
f ( t )=U ( t )g (t)+U ( t−a ) [h (t )−g(t )]
Y cuando t ≥ b necesitamos que h(t) se desactive:
f ( t )=U ( t )g (t)+U ( t−a ) [h ( t )−g(t )]−U (t−b )h(t )
15I N G E N I E R Í A E N S I S T E M A S - UNCAus
a
f(t)
0
g(t) h(t)1
t
t
a
t
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
Entonces ya tenemos definida f(t) como una función de escalones unitarios, veamos todo esto gráficamente:
Ahora vamos el segundo teorema de traslación:
A continuación vemos la grafica de f ( t )=t 3 con t ≥0
Veamos lo que pasa si nosotros escribimos f ( t−a ):
Lo que sucede es cuando le restamos una constante al dominio, corremos la función hacia a (en lugar de estar el origen en “0” está en “a”), y como consecuencia de esto empezamos a ver una parte de la función que antes no veíamos por la condición de t ≥0 (en el grafico es la parta roja). Para calcular la transformada de Laplace de la
16I N G E N I E R Í A E N S I S T E M A S - UNCAus
b
a
t
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
función corrida en “a” tendremos que quitarle esa porción de curva que apareció, lo que hacemos entonces:
f ( t−a )U (t−a)
Multiplicamos la función corrida en a por el escalón unitario de menos “a”, y esto produce:
Que es lo que nosotros queríamos y tenemos un corrimiento puro de f ( t−a ), ahora si calculamos la transformada de Laplace de f (t−a )U (t−a):
£ ( f ( t−a )U (t−a))=∫0
∞
f (t−a )U ( t−a)e−St dt
En el intervalo de 0< t<∞ U (t−a) vale “0” entonces:
∫0
a
f ( t−a )U (t−a)e−St dt=0
Ahora trabajamos con
∫a
∞
f (t−a ) e−St dt
∫0
∞
f (u ) e−S (u+a)du=∫0
∞
f (u ) e−Su e−Sadu=e−Sa∫0
∞
f (u )e−Su du
17I N G E N I E R Í A E N S I S T E M A S - UNCAus
CALCULO AUXILIAR
u=t−a;du=dt;t=u+a
Limites de integración
u=a−a=;u=∞−a=∞
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
Podemos observar aquí que:
∫0
∞
f (u ) e−Sudu=∫0
∞
f ( t ) e−St dt
Entonces:
e−Sa∫0
∞
f (u ) e−Su du=e−Sa £ (f (u))
O bien
e−Sa∫0
∞
f ( t ) e−St dt=e−Sa£ ( f (t))
Por lo tanto:
£ ( f ( t−a )U (t−a))=e−SaF (s)
Se puede ver muy bien que el segundo teorema de traslación es muy parecido al primer teorema de traslación, pero lo que está pasando es que al multiplicar por una exponencial a la transformada de Laplace provoca una traslación en nuestro dominio t, y dicha traslación es pura porque solo nos queda la parte positiva de f de la variable independiente t, lo que es distinto al primer teorema que nos provocaba un corrimiento de s hacia ”a”
Transformada de una derivada
La derivación de una función f de la variable independiente t corresponde simplemente a la multiplicación de la transformada F por s. Esto permite reemplazar las operaciones del cálculo por simples operaciones algebraica sobre transformadas
£ { f (t)}=∫0
∞
f ( t)e−St dt
£ { f ' (t)}=∫0
∞
f ' (t )e−St dt
Vamos a trabajar con: £ { f ' (t)} y como es una integral impropia tenemos que cambiar los límites de integración
18I N G E N I E R Í A E N S I S T E M A S - UNCAus
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
£ { f ' (t)}=∫0
∞
f ' (t )e−St dt= limb→∞
∫0
b
f ' (t)e−St dt
Tenemos f '( t)dt y eso se siente como que es el diferencial de la función. Y como hay un diferencial vamos a hacer una integral por partes
Entonces:
£ { f ' (t)}=∫0
∞
f ( t)e−St dt=limb→∞
∫0
b
f ( t)e−St dt=limb→∞
e−St f ( t ) ¿0b−∫
0
b
f ( t )−Se−St dt=¿ limb→∞
e−St f ( t )¿0b+S∫
0
b
f (t )e−St dt ¿
Miramos:
limb→∞
e−St f (t ) ¿0b= lim
b→∞e−Sb f (b )−e−S0 f (0 )
En el caso de: limb→∞
e−Sb f (b ), el signo menos del exponencial hace que este se convierta
en denominador y esto ya lo hemos visto anteriormente. Sea Q (t )=e−Sb y P ( t )=f (b ) y P(t )Q(t)
, como en el denominador se encuentra el exponencial (todo esto para que la
función converja y exista la transformada de Laplace), y este crece mucho mas rápido que el denominador cuando b tiene al infinito el resultado de ese cociente va a ser “0”. Nos queda:
£ { f ' (t)}=−e−S 0 f ( 0 )+S∫0
∞
f (t ) e−St dt
Y como e0=1:
£ { f ' (t)}=−f (0 )+S∫0
∞
f ( t ) e−St dt
Y esto es equivalente a decir:
£ { f ' (t)}=S £ {f ( t)}−f (0 )
19I N G E N I E R Í A E N S I S T E M A S - UNCAus
CALCULO AUXILIAR
u=e−St ;du=−s e−St dt
v=f ( t ); dv=f ' ( t )dt
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
Que es lo que queríamos mostrar, y en base a esto ya podemos ir subiendo el orden de derivación:
£ { f ' ' (t)}=S £ {f ' (t)}−f ' (0 )
Y £ { f ' (t)}=S £ {f ( t)}− f (0 ) y entonces:
£ { f ' ' (t ) }=S [S£ { f (t ) }−f (0 ) ]−f ' (0 )=S2£ {f ( t ) }−Sf (0 )−f ' (0 )
Y así podemos generalizar para derivadas de orden “n”:
£ { f (n )(t ) }=Sn£ {f (t ) }−Sn−1 f (0 )−Sn−2 f (0 )−⋯−f (n−1) (0 )
Transformada de Laplace: Aplicaciones
Un sistema dinámico puede definirse conceptualmente como un ente que recibe unas acciones externas o variables de entrada, y cuya respuesta a estas acciones externas son las denominadas variables de salida. Las acciones externas al sistema se dividen en dos grupos, variables de control, que se pueden manipular, y perturbaciones sobre las que no es posible ningún tipo de control. La Figura 3 ilustra de un modo conceptual el funcionamiento de un sistema de control.
El campo de aplicación de los sistemas de control es muy amplio y una herramienta que se utiliza en el diseño de control clásico es precisamente la TRANSFORMADA DE LAPLACE.
¿Por qué Transformada de Laplace?
En el estudio de los procesos es necesario considerar modelos dinámicos, es decir, modelos de comportamiento variable respecto al tiempo, esto trae como consecuencia el uso de ecuaciones diferenciales respecto al tiempo para representar matemáticamente el comportamiento de un proceso. El problema es que las ecuaciones diferenciales que representan el comportamiento dinámico de los procesos naturales son muy complejas y es ahí donde entra en juego la transformada de Laplace que permite resolver ecuaciones diferenciales lineales mediante la transformación en ecuaciones algebraicas con lo cual se facilita su estudio
Ejemplo: Suspensión de automóviles
Para poder diseñar un sistema de control automático, se requiere:
20I N G E N I E R Í A E N S I S T E M A S - UNCAus
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
Conocer el proceso que se desea controlar, es decir, conocer la ecuación diferencial que describe su comportamiento, utilizando las leyes físicas, químicas y/o eléctricas.
A esta ecuación diferencial se le llama modelo del proceso. Una vez que se tiene el modelo, se puede diseñar el controlador.
Después:
21I N G E N I E R Í A E N S I S T E M A S - UNCAus
f(t)z(t)
k b
m
Fuerza de entrada
Desplazamiento, salida del sistema
∑ F=ma
f ( t )−kz ( t )−bdz ( t )dt
=md2 z ( t )dt 2
f ( t )−kz ( t )−bdz( t )dt
=md2 z( t )dt 2
Aplicando la transformada de Laplace a cada término(considerando condiciones iniciales igual a cero)F ( s )−kZ ( s )−bsZ ( s )=ms2 Z (s )F ( s )=Z( s ) [ms2+bs+k ]Z ( s )F ( s )
=1
ms2+bs+k
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
Función de transferencia
Una función de transferencia es un modelo matemático que a través de un cociente relaciona la respuesta de un sistema (modelada) a una señal de entrada o excitación (también modelada). Entonces resumiendo este enunciado tenemos que la función de transferencia:
Representa el comportamiento dinámico del proceso Nos indica cómo cambia la salida de un proceso ante un cambio en la entrada
Para el ejemplo de la suspensión del automóvil tenemos que:
Convolucion
22I N G E N I E R Í A E N S I S T E M A S - UNCAus
Proceso
Entrada del proceso
(función forzante o
estímulo)
Salida del proceso
(respuesta al
estímulo)
Entrada
(Bache)
Salida
(Desplazamiento del coche)kbsms 2
1
0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000-10
-8
-6
-4
-2
0
2
4
6
8
10
0 0.5 1 1.5 2 2.5 3
x 104
-4
-3
-2
-1
0
1
2
3x 10
-3
Y ( s )X ( s )
=Cambio en la salida del procesoCambio en la entrada del proceso
Y ( s )X ( s )
=Respuesta del procesoFunción forzante
MATEMATICA SUPERIOR – INGENIERIA EN SISTEMAS
Si f y g de la variable independiente t son las transformadas inversas de F(s) y G(s), respectivamente, y satisfacen las hipótesis del teorema de existencia, entonces la transformada inversa h(t) del producto H(s) F(s)G(s) es la convolucion de f (t) y g(t ) denotando ( f∗g )( t) y se define por:
h (t )=( f∗g ) ( t )=∫0
t
f (∂ )g ( t−∂ ) . d ∂
Las propiedades de la convolucion de transformadas son:
f∗g=g∗f (ley conmutativa)
f∗(g1+g2 )=f∗g1+ f g2(ley distributiva)
( f∗g )∗v=f∗(g∗v)(ley asociativa)
23I N G E N I E R Í A E N S I S T E M A S - UNCAus
top related