Matemática Aplicada: Clase 17
Alex Antequeda
Universidad de La Frontera
A. Antequeda IME278 IME278
Segundo semestre 2023 1 / 12
Resumen
Aplicación: Resolución de EDPs
Producto de convolución
Delta de Dirac
A. Antequeda IME278 IME278
Segundo semestre 2023 2 / 12
Aplicación: Resolución de EDPs
Ejemplo
Problema “físico”: Distribución de temperatura T = T (x, t) de una barra
suficientemente aislada, de longitud infinita
2
Condición inicial: T (x, 0) = e−x
2
∂T = ∂ T ,
t > 0, x∈R
PVI: ∂t ∂x2
T (x, 0) = e−x2 ,
x∈R
Para cada t > 0 fijo pero arbitrario, podemos pensar T (x, t) como una función de
x, por lo que, aplicando la transformada de Fourier y la propiedad de la derivada
se obtiene lo siguiente:
" #
∂T ∂2T
Fx (λ) = Fx = −λ2 Fx [T ](λ)
∂t ∂x2
donde Fx denota la transformada de Fourier sólo en la variable x
A. Antequeda IME278 IME278
Segundo semestre 2023 3 / 12
Ejemplo
Pero Z ∞
∂T 1 ∂T
Fx (λ) = √ (x, t)e−iλx dx
∂t 2π −∞ ∂t
y si intercambiamos la derivada con la integral se tiene:
Z ∞
∂T ∂ 1 ∂
−iλx
Fx (λ) = √ T (x, t)e dx = Fx [T ](λ)
∂t ∂t 2π −∞ ∂t
∂
=⇒ Fx [T ](λ) = −λ2 Fx [T ](λ)
∂t
Se sigue que Fx [T ](λ) = Fx [T (x, t)](λ) , como función de t, es solución de la
dy
ecuación = −λ2 y(t). Por tanto,
dt
2t
Fx [T ](λ) = α(λ)e−λ
para alguna constante α(λ).
A. Antequeda IME278 IME278
Segundo semestre 2023 4 / 12
Ejemplo
Evaluando en t = 0 y usando la condición inicial:
2
Fx [T (x, 0)](λ) = Fx [e−x ](λ) = α(λ)
Ejercicio propuesto (Ver nota en el Campus Virtual)
2 1 λ2
Para a > 0, F[e−ax ] = √ e− 4a
2a
Luego
2 1 2 1+4t
Fx [T (x, t)](λ) = α(λ)e−λ t = √ e−λ ( 4 )
2
y aplicando nuevamente la transformada de Fourier,
1 h 2 1+4t
i 1 λ2
T (−λ, t) = √ Fx e−x ( 4 ) (λ) = √ e− 1+4t
2 1 + 4t
1 x2
Si cambiamos −λ por x, la solución nos queda: T (x, t) = √ e− 1+4t
1 + 4t
A. Antequeda IME278 IME278
Segundo semestre 2023 5 / 12
Producto de convolución
Convolución
Definición
Sean f, g : R → R dos funciones que tienden rápidamente a cero. La convolución
de f y g se define como
Z ∞
(f ∗ g)(x) = f (u)g(x − u) du
−∞
Observaciones.
• La integral converge: como g es continua y lim |g(x)| = 0, entonces g es
|x|→∞
acotada (¿por qué?)
Si |g(x)| ≤ M para algún M > 0 y todo x ∈ R, se tiene
Z ∞ Z ∞
f (u)g(x − u) ≤ M |f (u)|du < ∞
−∞ −∞
• Es posible, pero tedioso, demostrar que si f y g tienden rápidamente a cero,
entonces f ∗ g tiende rápidamente a cero.
A. Antequeda IME278 IME278
Segundo semestre 2023 6 / 12
Propiedades de la convolución
Teorema
Bajo las hipótesis de la definición anterior, tenemos que
• La convolución es conmutativa: f ∗ g = g ∗ f
• La convolución es asociativa: (f ∗ g) ∗ h = f ∗ (g ∗ h)
• La convolución es distributiva: f ∗ (g + h) = f ∗ g + f ∗ h
• La convolución respeta shifts en el tiempo: si f (t) = g(t) ∗ h(t), entonces
f (t − t0 ) = g(t − t0 ) ∗ h(t)
Observación: Todas estas propiedades se deducen fácilmente. Por ejemplo
Z ∞ Z −∞
(f ∗ g)(x) = f (u)g(x − u) du = f (x − v)g(v)(−dv) = (g ∗ f )(x)
−∞ | {z } ∞
v
A. Antequeda IME278 IME278
Segundo semestre 2023 7 / 12
Transformada de Fourier de la convolución
Teorema
Sean f, g : R → R dos funciones que tienden rápidamente a cero. Entonces:
√
F[f ∗ g] = 2πF[f ]F[g]
Demostración imprecisa: Saque la cuenta
1
Z∞
F[f ∗ g](λ) = √ (f ∗ g)(u)e−iλu du
2π −∞
Z ∞ Z ∞
1
=√ f (t)g(u − t) dt e−iλu du
2π −∞ −∞
Z ∞ Z ∞
1
=√ f (t)g(u − t)e−iλu dt du
2π −∞ −∞
Pero e−iλu = e−i(u−t)λ e−iλt
A. Antequeda IME278 IME278
Segundo semestre 2023 8 / 12
Transformada de Fourier de la convolución
v = u − t
Hacemos el cambio de variables que tiene Jacobiano
w = t
! !
∂v ∂v
∂u ∂t 1 −1
det ∂w ∂w = det =1
∂u ∂t 0 1
Entonces tenemos
1 ∞Z ∞ Z
F[f ∗ g](λ) = √ f (w)g(v)e−ivλ e−iwλ dv dw
2π −∞ −∞
Z ∞ Z ∞
1
=√ f (w)e−iwλ dw g(v)e−ivλ dv
2π | −∞ −∞
√ √
{z }| {z }
2πF [f ] 2πF [g]
√
= 2π F[f ] F[g]
A. Antequeda IME278 IME278
Segundo semestre 2023 9 / 12
Delta de Dirac
Función de Heaviside
Definición
La función escalón de Heaviside se define como
0, x<a
µ(x − a) =
1, x≥a
1.0
0.8
0.6
0.4
0.2
1 2 3 4
A. Antequeda IME278 IME278
Segundo semestre 2023 10 / 12
Función cajón
Dado a > 0 y r > 0, defina una función cajón apropiada
r
0, |x − a| > 2
hr (x − a) =
1, r
r |x − a| < 2
En términos de la función de Heaviside, podemos describir esta función como
1 r r
hr (x − a) = µ x−a+ −µ x−a−
r 2 2
Observación. Para cualquier r > 0 se verifica lo siguiente:
Z ∞ Z a+ r
2 1
hr (x − a)dx = dx = 1
−∞ a− r2 r
Por lo tanto, Z ∞
lim hr (x − a)dx = 1
r→0 −∞
A. Antequeda IME278 IME278
Segundo semestre 2023 11 / 12
Idea intuitiva de la delta de Dirac
Lo anterior sugiere definir el límite δ(x − a) := lim hr (x − a), llamada delta de
r→0
Dirac que conserve las dos propiedades:
Z ∞
δ(x − a) = 0 si x ̸= a y δ(x − a)dx = 1
−∞
OJO: Esto es puramente moral. En realidad ninguna de esas cosas existe
Más precisamente: las propiedades listadas contradicen matemática básica, pues si
f (x) = 0 para todo x ̸= a, entonces
Z ∞
f (x)dx = 0
−∞
La verdad: La delta de Dirac es una distribución
Sugerencia. Para la siguiente clase repase álgebra lineal. Por ejemplo, la definición
de espacio vectorial y de transformación lineal.
A. Antequeda IME278 IME278
Segundo semestre 2023 12 / 12
