Matemáticas V
UNIDAD 3 TRANSFORMADA DE LAPLACE
-FUNCIÓN DELTA DIRAC
FUNCION DELTA DIRAC. Impulso unitario con frecuencia, sobre los sistemas mecánicos actúan fuerzas externas (o fem sobre los circuitos eléctricos) de gran magnitud solo durante un lapso muy breve por ejemplo, en un ala de aeroplano que se encuentre oscilando, puede caer un rayo, se puede dar un golpe brusco a una masa en un resorte con un martillo de bola, o una bola de beisbol (golf o tenis), podría mandarse volando golpeándola violentamente con algún tipo de garrote, como un bate, palo de golf o una raqueta. La función
a(t – t0)=
Cuando a> 0, t 0>0 se ven en la figura a), y podrían servir como modelo matemático de este tipo de fuerzas. Para los valores pequeños de es, esencialmente, una función
∫
constante de gran magnitud que se encuentra “encendida” sólo durante un lapso muy
pequeño alrededor de figura b). Esta función integración,
. El comportamiento de cuando a 0 se muestra en la , se llama impulso unitario porque tiene la prioridad de la .
1
1
Zill, Dennis G. (1997). Ecuaciones diferenciales con aplicaciones de modelado. México: International Thompson Editores.
Función delta Dirac. En la practica conviene trabajar con otro tipo de impulso unitario, con una función que aproxima a(t – t0), definida con el limite.
(t – t0) =
Esta ultima expresión, que por ningún motivo es una función, se puede caracterizar mediante las dos propiedades siguientes:
, a
(t – t0)
Impulso unitario
∫ –
(ii)
dt = 1
se denomina función delta Dirac.
Es posible obtener la transformada de Laplace de la función delta de Dirac con la hipótesis formal de que
} }
Teorema: Transformada de la función delta de Dirac .
Para
}
( )( )
DEMOSTRACIÓN: DEMOSTRACIÓN: comenzaremos expresando escalón unitario, de acuerdo con las ecuaciones.
en terminos de la función
Según la linealidad, la transformada de Laplace de esta expresión es
}
Como esta ecuación tiene la forma indeterminada 0/0 cuando a 0, aplicamos la regla de L’Hôpital:
}} }
Cuando que
parece lógico suponer de acuerdo con la ecuación descrita en el teorema,
}
Este resultado subraya el hecho de que no es el tipo normal de función que tenemos manejando porque de acuerdo con el teorema, esperaríamos que cuando
∫
}
Definición de función delta:
es un objeto matematico conocido como la función delta. Es un ejemplo de una distribución o función generalizada. generalizada . Tiene las siguientes propiedades: ii) ii)
iii)
iv) si g(t) es una función continua en
entonses
Es posible construir una definición lógicamente rigurosa pero no se hará intuitivamente, se puede pensar en como en una aproximación de un impulso físico de magnitud 1 en el tiempo t=0 por ejemplo podría ser la transferencia rápida de una unidad de carga en el tiempo de carga en el tiempo cero. Se puede demostrar que si a es una constante entonces
v)
∫
vi)
vii) si g(t) es una función continua en (
entonces
2
De estas formulas se obtiene que.
Entonces, formalmente
,
Y se puede considerar, en cierto sentido, como la derivada de la función Heaviside. Una propiedad agradable de la transformada de Laplace es que se aplica a manera casi c asi tan sencilla a impulsos y distribuciones como a las funciones ordinarias.
No es necesario decir que en los problemas reales que involucran impulsos, debe tenerse cuidado de asegurar que las ecuaciones que se están usando sigan proporcionando modelos exactos en la presencia de valores grandes, pero breves dados por el impulso.
2
Stephen L. Campbell, Richard Haberman
Introducción a las ecuaciones diferenciales con problemas de valor de frontera, 1ra. Ed.
Función Delta Dirac El producto de la convolución puede usarse para encontrar la respuesta de un sistema dinámico a una fuerza súbita de amplitud larga y duración corta (como golpear una bola de beisbol). Esto puede modelarse mediante la función delta de Dirac. Empezaremos en un plano hipotético, “definiendo” la función delta por una propiedad que se requiere que
posea.
La función delta de Dirac y sus propiedades Comencemos por ampliar los dominios de las funciones f en E a toda la recta real al tomar de valor de f(t) como cero para toda t negativa. Este nuevo conjunto de funciones se denota con . Ahora ya podemos dar la definición de un objeto muy extraño al que Dirac llamó función generalizada.
Función delta de Dirac. Supóngase que existe un elemento δ en toda y toda función que es continua en ,
∫
tal que para
(1)
El elemento δ se denomina función delta de Dirac. Dirac.
En realidad no puede haber mas de una función delta de Dirac. Para verlo, digamos que hay dos funciones que satisfacen la definición (1). Entonces, para cualquier punto es continua,
3
3
Paul Adrien Maurice Dirac (1902-1084), físico teórico británico, fue premiado (junto con Erwin Schrödinger) con el premio Nobel de física por su trabajo en mecánica cuántica. Ente otros logros, Dirac describió el movimiento movimiento de un electrón mediante cuatro ecuaciones diferenciales simultaneas. Con e hallazgo de estados negativos de energía que predijo con este modelo matemático, Dirac planteo la hipótesis de la existencia de positrones, o antipartículas electrónicas, hipótesis que se co nfirmó mas tarde. 2
Robert L. Borreli, Courthey S. Coleman. Ecuaciones diferenciales. Ed. OXFORD UNIVERSITY PRESS. Pág. 397
Por tanto,
porque
en todos los valores t donde δ y µ son continuas.
Ahora supongamos que existe una función delta de Dirac. Es posible que el lector sienta incomodidad con el calculo de una función que podría no existir; no obstante, proseguiremos por este camino.
TEOREMA: Propiedades de la función delta de Dirac. Se tiene
Para comprobar la primera propiedad, sea f(t)=escalón(t), que es la función de Heaviside. De (1) puede establecerse que u=t-u y demostrarse que
∫ ∫ ∫
Por tanto,
(3)
. Y se demuestra la primera desigualdad de
(2).
Para demostrar la segunda propiedad se emplea el teorema de la convolución.
(4)
Por otro lado de tiene que
(5)
Y que (después de declarar
)
(6)
Donde se han utilizado identidad (3). Por tanto (4),(5) y (6) se observa que
y la
.
La función delta se utiliza a menudo para resolver PVI. Supóngase que para las constantes reales a y b. Entonces, el siguiente resultado da una formula de solución para un PVI:
Teorema: la función delta y problemas de valor inicial. La solución del PVI.
Donde f está en
, está dada por
Donde para cada u la función
es la solución del PVI
Para ver por qué la formula (8) da la solución del PVI (7) se razona como sigue: si
*+ ∫ ∫ ∫ , entonces la solución del PVI(7) está dada por
=
Donde se ha aprovechado que otra parte, se ve que la transformada transformada G(t,u) G(t,u) del PVI está dada por
En consecuencia, por (11) y la fórmula del teorema del corrimiento se observa que.
Por ultimo de (10) y (12) se tiene la fórmula buscada (8):
. Por
Por consiguiente, se ha planteado otra forma de encontrar la solución única de un PVI. ¿Existe la solución delta de Dirac?
Por el momento, supóngase que es una función en y veamos a donde nos lleva esto. Digamos que es un punto de continuidad para δ y que . Entonces hay un intervalo en el que es positiva y continua. Consideremos la función Entonces, como que por (1) se tiene que
Donde también se ha utilizado el hecho de que cambio de variables se observa por la elección de
al hacer el
que
Puesto que hemos supuesto que es positiva para esta contradicción indica que no puede ser positiva en ningún punto de la continuidad. De manera similar, tampoco puede ser negativa y, por tanto, δ desaparece en todo punto donde es continua. Esto significa que el miembro izquierdo de (1) desaparece para toda t, sin importar qué función Pero esto es absurdo de modo que cualquiera que sea , no es una función en . Teníamos la sospecha de que a partir de la segunda fórmula de (2) ya que δ no puede ser una función, pues cualquier función en debe tener una transformada que tiende a cero cuando . Desde los tiempos de Dirac, las “funciones” como δ han sido muy importantes en las
aplicaciones. En los tratamientos avanzados de las modernas matemáticas aplicadas se construye una teoría lógicamente rigurosa que incluye objetos, conocidos como distribuciones, funciones generalizadas o función simbólica. Que se comportan como la “función delta”.
Ejercicio 1: Función de forzamiento impulsiva Resuelva Solución:
Tomando la transformada de Laplace en ambos lados se tiene:
* + * +
Despejando Y(s),
De manera que:
(1)
Físicamente, este ejercicio se puede ver como el circuito RC lineal simple de la figura 3.12.1, donde y es la carga en el capacitor del tiempo, y existe una carga inicial de 1 en el capacitor. Para , el voltaje e es cero, y el capacitor se esta descargando. En el tiempo t=1, hay un impulso de voltaje, es decir se aplica un voltaje muy alto durante un periodo breve, lo que recarga el capacitor. Después el voltaje es 0 otra vez y el capacitor continúa descargándose. La grafica de (1) se puede ver en la figura 3.12.2. Esta grafica debe interpretarse como que, en una problema real, y(t) estaría dada por una función como la de la figura f igura 3.12.3.
Fig. 3.12.1 4
4
Fig. 3.12.2 Gráfica de (1)
Fig. 3.12.3
Stephen L. Campbell, Richard Haberman . Introduccion a las ecuaciones diferenciales como problemas de valor de frontera. Pág. 288
Ejercicio 2: Una masa unida a un resorte se libera desde el reposo, a un metro por bajo de la posición de equilibrio para el sistema masa-resorte y comienza a vibrar. Después de π segundos, la masa es golpeada por un martillo que ejerce un impulso sobre la masa. El sistema queda descrito por el problema simbólico con valores iniciales.
(I)
Donde x(t) denota el desplazamiento con respecto del equilibrio en el instante t. Determinar x(t).
Solución:
} } } } }
Sea
Como
Al calcular la transformada de de Laplace de ambos lados de (I) y despejar X(s) tenemos:
(s)
Usamos la propiedad de translación para determinar la transformada inversa de Laplace de X(s) Translación en t. Teorema: Suponga que positiva, entonces:
} } } } existe para
Si α es una constante
,
Recíprocamente, una una transformada inversa inversa de Laplace de
esta dada por:
sí tenemos que:
, { √
La grafica de x(t) aparece en color en la figura 3.12.4. Como comparación, la curva punteada exhibe el desplazamiento de un resorte vibrante sin perturbaciones. Observe que el impulso suma 3 unidades al momento en el instante t= π.
Desplazamiento de un resorte vibrante golpeado por un martillo t= π
Grafica 3.12.4
5
5
R. Kent Nagle, Edward B. Saff, Arthur David Snider . Ecuaciones diferenciales con problemas con valores de
frontera .Pág. 433
Ejercicio 3: Suponga que tratamos de resolver:
[ ]
Con condiciones iniciales . Supongamos que buscamos una solución continua. Tomando transformadas de Laplace con como es usual, obtenemos
Para esta inversa necesitamos el segundo teorema de translación.
Segundo teorema de translación Sea f(t) una función que tiene una una transformada transformada de Laplace y sea sea positiva. Entonces
Forma 1: Forma 2: Se establece el resultado que es:
α u
, ,
La grafica de esta solución se muestra en la figura 3.12.5
Fig. 3.12.5
6
6
Daniel A. Marcus . Ecuaciones Diferenciales. Pág. 448
Bibliografía
Zill, Dennis G. (1997). Ecuaciones diferenciales con aplicaciones de modelado. modelado. México: International Thompson Editores. Stephen L. Campbell, Richard Haberman Introducción a las ecuaciones diferenciales con problemas de valor de frontera, 1ra. Ed.
McGRAW-HILL, México, 1998
Robert L. Borrelli, Courtney S. Coleman Ecuaciones Diferenciales, una perspectiva de modelación, 1a. ed.
Oxford university press, México, 2002.
Daniel A. Marcus Ecuaciones Diferenciales , 1a. ed.
Compañía editorial continental, México, 1993
R. Kent Nagle, Edward B. Saff, Arthur David Snider. Ecuaciones diferenciales diferenciales y problemas problemas con valores de frontera, 3a. ed.
Pearson Educación, México, 2001
Stephen L. Campbell, Richard Haberman Introducción a las ecuaciones diferenciales con problemas de valor de frontera ,
1ra. Ed. McGRAW-HILL, México, 1998
