Ecuaciones Diferenciales Parciales: Guía Completa
Ecuaciones Diferenciales Parciales: Guía Completa
Unidad I
Conceptos Básicos.
Solución de una E.D.P. mediante las E.D.O.
Ecuaciones en derivadas parciales lineales de primer orden.
Ecuaciones en derivadas parciales lineales de segundo orden.
Método de solución de una E.D.P. de segundo orden.
Unidad II
Series de Fourier conceptos básicos.
Teorema de convergencia de una serie de Fourier.
Serie compleja de Fourier.
Integral de Fourier.
Integral compleja de Fourier.
Unidad III
Transformada de Fourier.
Propiedades de la transformada de Fourier.
Ecuación de la onda.
Ecuación de la onda. Modelamiento de la ecuación de la cuerda vibrante.
Ecuación de la onda. Modelamiento de la ecuación del calor.
Ecuación de Laplace.
1
1. Ecuaciones Diferenciales Parciales
1.1. Definición
Se denomina ecuación diferencial en derivadas parciales (E.D.P.) a una ecuación don-
de interviene una función incógnita de dos o mas variables independientes y derivadas
parciales de orden uno o mayor de la función incógnita.
Sea µ = µ(x1 , x2 , ..., xn ) una función en n variables independiente. Luego una E.D.P.
se puede representar por:
F (x1 , x2 , ..., xn , µ, µx1 , µx2 , ..., µxn , µx1 x2 , ...) = 0 (1)
Luego (1) es una E.D.P. en la función µ, ademas F es una función en estos argumentos.
1.1.1. Ejemplo
Sea µ = µ(x, y) una función en dos variables independientes, algunas E.D.P. son:
1. µx + µy + µ = 0, esto es F (µ, µx , µy ).
2. µxx − 3µxy + 2µyy − µ = 0, esto es F (x, y, µ, µxx , µxy , µyy ) = 0
∂µ
3. µxx − 5µyy + µx − 3µy − 8µ = 0, esto es F (x, y, µ, µx , µy , µxx , µyy ) donde µx = ∂x
;
2
µy = ∂µ
∂y
; µxx = ∂∂xµ2 ; µxy = ∂y
∂ ∂µ
( ∂x )
1.2. Definición
Se denomina orden de una E.D.P., al mayor orden superior de las derivadas parciales
de la función incógnita.
1.3.1. Ejemplo
Sea la E.D.P. de segundo orden µyy (x, y)−µx (x, y)−µ(x, y) = 0. Verificar si la función
µ(x, y) = e−2x cos(y) es una solución de la E.D.P.
solución
De lo estudiado se sabe que la función µ(x, y) = e−2x cos(y) es una función diferencia-
ble en todo R2 y ademas µ es de clase C ∞ en R2 .
2
Derivando se obtiene
µx (x, y) = −2e2x cos(y); µyy (x, y) = −e−2x cos(y)
Sustituyendo en la ecuación
−e−2x cos(y) − (−2e−2x cos(y)) − e−2x cos(y) = 0
luego µ(x, y) = e2x cos(y), si es una solución de la ecuación dada.
1.3.2. Observación
1.- Ahora surge una pregunta, ¿como se puede resolver una E.D.P. de una función
incógnita µ?.
1.3.3. Ejemplo
∂µ(x,y)
Halle la solución de la E.D.P. ∂x =0
Solución
∂µ(x,y)
Si ∂x = 0, ello significa que la función µ(x, y) no depende de la variable x, solo
depende de la variable y.
Luego la función µ(x, y) va a ser una función cualquiera de la forma ϕ(y), donde ϕ es una
función real arbitraria pero derivable. Luego la solución de la E.D.P. es µ(x, y) = ϕ(y).
1.3.4. Ejemplo
Halle la solucion de la E.D.P.:
∂ 2 µ(x, y)
=0 (i)
∂x∂y
Solución
∂ 2 µ(x, y)
Consideremos la ecuación = 0 de la forma:
∂x∂y
∂ ∂µ(x, y)
( )=0 (ii)
∂x ∂y
∂µ(x, y)
para que la E.D.P. (ii) sea cero, significa que la función no depende de la variable
∂y
x, sino de la variable y entonces hacemos
∂µ(x, y)
v(x, y) = (iii)
∂y
3
∂v(x, y)
sustituyendo en (ii) se obtiene = 0, luego de forma similar, del ejemplo anterior,
∂x
se obtiene v(x, y) = ϕ(y), donde ϕ es una función real arbitraria.
Luego reemplazando esto en (iii) se obtiene
∂µ(x, y)
= ϕ(y) (iv)
∂y
Z Z
integrando (iv) respecto a y (fijado x) se obtiene ∂µ(x, y) = ϕ(y) dy, luego se obtiene
Z
µ(x, y) = ϕ(y) dy + h(x) (v)
R
donde h es una función arbitraria que depende de x, hacemos ψ(y) = ϕ(y)dy y reem-
plazamos en (v) y se obtiene µ(x, y) = ψ(y) + h(x), esto representa la solución general de
la E.D.P.(i).
1.3.5. Observación
1.- En una E.D.O. de orden n se sabe que la solución general de la ecuación, involucra
n constantes reales Estas constantes cuando la solución requerida, satisface ciertas
condiciones iniciales o de frontera.
2.- En el caso de las E.D.P., la solución general involucra funciones arbitrarias en vez
de constantes reales. El numero de estas funciones arbitrarias es generalmente igual
al orden de la E.D.P. Ademas, si la función incógnita tiene m variables independien-
tes (µ = µ(x1 , x2 , ..., xm )), entonces las funciones arbitrarias dependen de (m − 1)
independientes.
1.3.6. Ejemplo
Halle la solución general de la E.D.P.:
µxy = cos(x) + y (i)
Solución
∂ ∂µ
Se sabe µxy (x, y) = ∂y ( ∂x ), luego la E.D.P. es de la forma
∂ ∂µ
( ) = cos(x) + y (ii)
∂y ∂x
∂µ
supondremos que ∂x solo depende de la variable y (fijando x), entonces hacemos
∂µ
v(y) = (iii)
∂x
reemplazando en (ii) se obtiene
dv(y)
= cos(x) + y (iv)
dy
integrando (iv) respecto a la variable y se obtiene
4
Z Z
y2
dv(y) = (cos(x) + y) dy, v(y) = y cos(x) + 2
+ ψ(x)
∂µ y2
= y cos(x) + + ψ(x) (v)
∂x 2
donde ψ es una función real arbitraria.
Ahora integramos la ecuación (v) respecto a x (fijando y)
y2
Z
(y cos(x) + + ψ(x)) dx
2 Z
2
µ(x, y) = y sen(x) + x y2 + ψ(x) dx + h(y)
2
de aquı́ se obtiene µ(x, y) = y sen(x) + x y2 + ϕ(x) + h(y) donde ϕ y h son funciones reales
arbitrarias.
1.3.7. Observación
1.- En ciertas E.D.P. especiales, solamente aparece la derivada parcial de una de las
variables.
2.- En estos casos, para resolver dicha E.D.P., se puede tratar como una E.D.O., donde
se fija la variable independiente cuya derivada parcial no aparece.
1.3.8. Ejemplo
Halle la solución general de la E.D.P.
Solución
Se observa que en la E.D.P. (i) no aparece la derivada respecto a la variable indepen-
diente y. Luego la E.D.P.(i) se puede considerar como una E.D.O. respecto a la variable
independiente x, fijando la variable y.
Luego se tiene
5
1.3.9. Ejemplo
Resolver la E.D.P.
Solución
∂ ∂µ
Sabemos µyx = ( ), expresando la E.D.P. de la forma (i)
∂x ∂y
∂ ∂µ ∂µ
x2 ( ) + 4x = y2 (ii)
∂x ∂y ∂y
∂µ
supongamos que la función ∂y
solo dependa de x,luego
∂µ
= v(x) (iii)
∂y
sustituyendo (ii) en (iii) se obtiene
d
x2 v(x) + 4xv(x) = y 2 (iv)
dx
la ecuación (iv), tiene la forma de una E.D.O. lineal de primer orden en la variable x, de
(iv) se obtiene
4 y2
v 0 (x) + v(x) = 2 ; x 6= 0 (v)
x x
resolviendo la ecuación (v) se obtiene
Z
4
− dx Z R 4 y 2
v(x) = e x e x
dx
dx + φ(y)
x2
2
Z
−4 ln x 4 ln x y
v(x) = e e dx + φ(y)
x2
Z 4 2
−4 xy
v(x) = x dx + φ(y)
x2
Z
−4 2 2
v(x) = x x y dx + φ(y)
x3 2
v(x) = x−4 ( y + φ(y))
3
y2
v(x) = + x−4 φ(y)
3x
y2 φ(y)
v(x) = + 4 (vi)
3x x
6
Sustituyendo (iii) en (vi) se obtiene
∂µ(x, y) y2 φ(y)
= + 4 (vii)
∂y 3x x
esto es
y3
µ(x, y) = + x−4 h(y) + ψ(x) (viii)
9x
Z
donde h(y) = φ(y) dy y ψ son funciones arbitrarias.
1.3.10. Observación
1.- En ciertas E.D.P., consiste en obtener en algún dominio de las variables indepen-
dientes, una solución de la E.D.P. que satisface ciertas condiciones, las cuales son
llamadas condiciones auxiliares.
1.3.11. Ejemplo
Dada la E.D.P.:
∂ 2 µ(x, y)
= x + y3 (i)
∂x∂y
cuya solución satisface la condición siguiente
7
Solución
∂ ∂µ(x, y)
De la E.D.P. (i) se obtiene = x + y3
∂x ∂y
∂µ(x, y)
supongamos que la función solo depende de la variable x (fijo y) entonces hace-
∂y
∂µ(x, y)
mos v(x) = , sustituyendo esto la E.D.P. queda expresado por:
∂y
dv(x)
= x + y3 (ii)
dx
integrando respecto a x (fijo y) se obtiene
x2
Z
v(x) = (x + y 3 ) dx = + xy 3 + φ(y)
2
Luego retornando a la expresión inicial, se obtiene
∂µ(x, y) x2
= + xy 3 + φ(y) (iii)
∂y 2
φ es una función real arbitraria. Integrando (iii) respecto a la variable y (fijo x) se obtiene
yx2 xy 4
Z
µ(x, y) = + + φ(y) dy + ψ(x)
2 4
yx2 xy 4
µ(x, y) = + + h(y) + ψ(x) (iv)
2 4
la expresión (iv) representa la solución general de la E.D.P.(i). Ahora emplearemos las
condiciones de problema, de (iv) y la condición dada se obtiene
y y4
µ(1, y) = 2y 2 − 4y = + + h(y) + ψ(1)
2 4
despejando se obtiene
y y4
h(y) = 2y 2 − 4y − − − ψ(1) esto es
2 4
9y y 4
h(y) = 2y 2 − − − ψ(1) (v)
2 4
sustituyendo (v) en (iv) se obtiene
x2 y4 9y y 4
µ(x, y) = y + x + 2y 2 − − + ψ(x) (vi)
2 4 2 4
De la otra condición del problema se obtiene
µ(x, −2) = −x2 + 4x + 8 + 9 − 4 + ψ(x) = x + 8
µ(x, −2) = −x2 + 4x + 13 + ψ(x) = x + 8
8
Despejando ψ(x) se obtiene
ψ(x) = x + 8 + x2 − 4x − 13+ = x2 − 3x − 5 (vii)
reemplazando (vii) en (vi) se obtiene
x2 y4 9y y 4
µ(x, y) = y + x + 2y 2 − − + x2 − 3x − 5
2 4 2 4
el cual representa la única solución de la E.D.P.
1.3.12. Ejemplo
Dada la E.D.P.
2µxx (x, y) − µxy (x, y) − µyy (x, y) = 0 (i)
Demostrar que la función µ(x, y) = f (x − 2y) + g(x + y) es una solución general de la
E.D.P. (i); donde f y g son funciones arbitrarias de clase C 2 en R.
Solución
vx = 1 wx = 1
Luego: Luego expresamos
vy = −2 wy = 1
µ(x, y) = f (v(x, y)) + g(w(x, y)) (ii)
9
También
∂ ∂µ ∂ 0 ∂v ∂w
µxy = ( )= (f (v) + g 0 (w)) = f 00 (v) + g 00 (w)
∂y ∂x ∂y ∂y ∂y
00 00
µxy = −2f (v) + g (w) (v)
F (x,y, µ, µx , µy ) = 0 (i)
considerando que µ(x, y) = z y f (x, y) = p(x, y) + φ(x) + ϕ(y), entonces (ii) queda
expresado por:
z = f (x, y)
La solución de la E.D.P. (i) se interpreta como una superficie en R3 . Luego para funciones
arbitrarias φ(x) y ϕ(y), se tiene una familia de superficies. Sin embargo, si elegimos
una función particular φ(x) y otra función particularϕ(y), entonces se obtiene una única
solución particular.
1.5.1. Ejempo
1.- Dada la función
10
el cual es solución de una E.D.P. donde ϕ(·) es una función derivable en x. Halle
una E.D.P. de primer orden de tal forma que µ sea la solución general.
Solución
Sea la función
µ(x, y) = −3y 3 ϕ(x) − 6x + y (ii)
deirvando respecto a y
Solución
11
d xµx − µ
derivando (iii) respecto a y, se obtiene = 0; luego
dy y
(xµxy − µy )y − (xµx − µ)
= 0;
y2
yxµxy − yµy − xµx + µ = 0
luego se obtiene la E.D.P.: xyµxy − yµy − xµx + µ = 0
1.6.1. Definición
Sean V y W espacios vectoriales de funciones diferenciables hasta el orden k; esto es
V = {µ : D ⊂ Rn → R; µ ∈ C k en D}
W = {v : D ⊂ Rn → R; v ∈ C k en D}
Una aplicación lineal L : V → W es llamado operador de funciones, si ello transforma
cada función µ ∈ V en otra función Lµ ∈ W de otra clase de funciones.
1.6.2. Ejemplo
Sean los espacios V = {µ : R2 → R; µ es de clase C 2 },
W = {v : R2 → R; v es continua en R2 }, se define el operador L : V → W
∂ ∂2
L= −c 2
∂x ∂y
Luego si µ ∈ V
∂µ ∂ 2µ
L(µ) = −c 2,
∂x ∂y
luego se verifica que L(µ) ∈ W; pues:
∂µ ∂ 2µ
µ ∈ V es de clase C 2 ; esto es es continua, también es continua; luego L(µ) es
∂x ∂y 2
continua. Por lo tanto L es un operador de funciones.
1.6.3. Observación
1.- En el ejemplo el operador L transforma cada función µ que admite derivada parcial
∂µ ∂ 2µ
hasta segundo orden con continuidad, en una nueva función L(µ) = − c 2 , el
∂x ∂y
cual es continua.
2.- Este operador L, es un operador diferencial por que involucra básicamente las deri-
vadas parciales.
12
1.6.4. Definición
Un operador de funciones L : V → W es llamado operador lineal si se cumple:
∀u, v ∈ V y ∀α, β ∈ R
L(αu + βv) = αL(u) + βL(v)
1.6.5. Ejemplo
∂ ∂2
El operador diferencial L = − c 2 del ejemplo (1.6.2) es lineal.
∂x ∂y
Solución
Sea u, v ∈ V, α, β ∈ R, luego
∂2
∂
L(αu + βv) = − c 2 (αu + βv)
∂x ∂y
∂ ∂2 ∂ ∂2
= (αu + βv) − c 2 (αu + βv) esto pues y son funciones.
∂x ∂y ∂x ∂y 2
Luego,
2
∂ (αu) ∂ 2 (βv)
∂(αu) ∂(βv)
L(αu + βv) = + −c +
∂x ∂x ∂y 2 ∂y 2
∂u ∂ 2u ∂v ∂ 2v
=α − αc 2 + β − βc 2
∂x ∂y ∂v ∂y
2
∂ 2v
∂u ∂ u ∂v
=α −c 2 +β −c 2
∂x ∂y ∂x ∂y
= αL(u) + βL(v) Luego L es un operador lineal
1.6.6. Observación
1.- La definición de linealidad de un operador se puede extender para cualquier numero
finito de funciones.
Si u1 , u2 , ..., un son n funciones y c1 , c2 , ..., cn son n constantes arbitrarias, entonces
n
! n
X X
L ci ui = ci L(ui )
i=1 i=1
1.6.7. Definición
Sean L1 y L2 dos operadores lineales definidos en un mismo espacio común V. Se define
el operador suma por:
(L1 + L2 )(u) = L1 (u) + L2 (u) ∀u ∈ V
Ademas el operador suma es lineal.
13
1.6.8. Definición(Operador Producto)
Sean L1 y L2 dos operadores lineales definidos en un mismo espacio de funciones V,
se define el operador producto por:
De tal forma que L2 (u) y L1 (L2 (u)) estén bien definidas y pertenezcan a un cierto espacio
de funciones W.
Si L1 y L2 son operadores lineales se verifica que el operador producto (L1 · L2 ) es
lineal.
En efecto:
∀u, v ∈ V, ∀α, β ∈ R se tiene:
def
(L1 · L2 ) (αu + βv) = L1 (L2 (αu + βv))
L2 lineal
= L1 (αL2 (u) + βL2 (v))
L1 lineal
= L1 (αL2 (u)) + L1 (L2 (v))
= αL1 (L2 (u)) + βL1 (L2 (v))
= α(L1 · L2 )(u) + β(L1 · L2)(v).
1.6.9. Definición
Sea L un operador diferencial lineal, luego la ecuación de la forma: L(u) = f . . . (i)
donde u ∈ V y f una función dada, es llamado una E.D.P. lineal; luego:
1.6.11. Observación
1.- El principio de superposición es empleado en forma extensa para resolución de la
E.D.P. lineales.
Lvi = fi ; ∀i = 1, 2, ..., n
14
n
X
Donde f1 , f2 , ..., fn son funciones dadas, entonces las función V = vi satisface
i=1
n
X
la siguiente E.D.P.: L(V ) = fi .
i=1
3.- Dada la E.D.P. lineal no homogénea L(u) = f . . . (i). Supongamos que v es una
solución de la E.D.P. homogénea asociada a (i). Si consideramos a la función w una
solución particular de la E.D.P.(i), entonces la función z = v + w es solución de la
E.D.P.(i).
En efecto:
L(z) = L(v + w) = L(v) + L(w) = 0 + L(v) = f .Luego se verifica que z = v + w es
solución de la E.D.P.(i).
donde a, b y c son constantes reales tal que |a| + |b| > 0 y la función incógnita es u.
Ademas g(·) es una función conocida y esta definido en un cierto dominio Ω ⊂ R2
Luego fijando la variable x, se observa que la E.D.P.(ii) tiene la forma de una E.D.O. lineal
no homogénea en la variable y. Resolviendo de la E.D.P.(ii) asociada con una E.D.O. de
primer orden, se obtiene
Z Z
c c
− dy Z dy g(x, y)
u(x, y) = e b e b dy + φ(x)
b
15
donde φ es una función arbitraria en la variable x.
Nuestro objetivo en este caso, es anular una de las derivadas parciales, para ello se intro-
duce una nueva variable
ε = bx − ay (iii)
∂ ∂w ∂x ∂w ∂ε
ux = w(x, ε) = + = 0 + wε b
∂x ∂x ∂x ∂ε ∂x
Esto es pues ε = bx − ay;
εx = b
εy = −a
ε + ay
Despejando x = ;
b
∂x
=0
∂x
a
xy =
b
También
∂w(x, ε) ∂w ∂x ∂w ∂ε a
uy (x, y) = = + = wx − awε
∂y ∂x ∂y ∂ε ∂y b
16
Si fijamos la nueva variable ε, la E.D.P.(v), toma la forma de una E.D.O. lineal de primer
orden en la variable x esto es
c g(x, ε)
wx (x, ε)+ w(x, ε) =
a a
2.2.1. Observación
1.- De la E.D.P.(i), si g(x, y) = 0, es decir la E.D.P. es homogénea, entonces la solución
general de la E.D.P. es
c
u(x, y) = e− a x φ(bx − ay)
2.2.2. Ejemplos
1.- Dada la E.D.P.: Aux (x, y) + Buy (x, y) + Cu(x, y) = D . . . (i)
Halle la solución general. A, B, C, D son constantes reales.
Solución
C D
wx (x, ε) + w(x, ε) = (ii)
A A
fijando la variable ε, la E.D.P.(ii) toma la forma de la E.D.O. de primer orden,
donde su solución general es
Z
C
−C x x
w(x, y) = e A e A D
A
dx + φ(ε)
17
esto es
C C C
w(x, ε) = e− A x DA CA e A + e− A φ(ε)
C
w(x, ε) = DC + φ(ε)e− A x ; donde φ es una función arbitraria
Luego retornando a la variable original, se obtiene:
C
u(x, y) = D
C
+ φ(Bx − Ay)e− A x
2.- Dado la E.D.P. 3ux (x, y)−6uy (x, y)+8u(x, y) = 2x . . . (i). Halle la solución general.
Solución
18
Luego retornando a la variable original se tiene
x 9 8
u(x, y) = − + φ(6x + +3y)e− 3 x
4 64
donde φ es una función arbitraria.
Donde a(·), b(·) y c(·) son funciones que dependen de las variables independientes x, y
y de la función incógnita u(·). Es decir, son de la forma, a(x, y, u(x, y)) = h(x, y)u(x, y),
para obtener la solución de la E.D.P.(i) emplearemos el método de de las curvas integrales
y sus direcciones. Esto es, haremos uso de la interpretación geométrica y las derivadas,
para hallar la solución del problema.
De acuerdo con la interpretación geométrica, se tiene que el gráfico (superficie) de la solu-
ción u(·) del problema(i), tiene la propiedad que el vector normal a la superficie (gráfico
de u) es siempre ortogonal al vector (a, b, c); esto es, si la función u(x, y) = z es solución
dela ecuación(i), entonces considerando la función G(x, y, z) = u(x, y) − z, su gradiente es
→
−
∇G(x, y, z) = ( ∂u(x,y)∂x
, ∂u(x,y)
∂y
, −1). Luego el vector N = (ux , uy , −1) es el vector normal
a la superficie (gráfico de u) en un punto (x0 , y0 ) ∈ Dom(u); por lo tanto, como u es
→
−
solución de la ecuación(i), entonces el vector N es ortogonal al vector P0 = (a, b, c), esto
es (a, b, c)(ux , uy , −1) = 0 esto es aux + buy − c = 0, es decir, aux + buy = c.
De otro lado, si C es una curva en R3 , cuya ecuación parametrica es R(t) = (x(t), y(t), z(t));
entonces para que esta curva C pertenezca al gráfico de u (solución (i)) se debe cumplir
lo siguiente: el vector tangente a la curva C ( en cada punto), debe ser paralela al vector
direccional (a, b, c); esto es:
d
R(t) = (x0 (t), y 0 (t), z 0 (t)) k (a, b, c) o simplemente
dt
dx dy dz
( , , ) = (a, b, c) (ii)
dt dt dt
19
2.3.1. Observación
1- De la ecuación(ii) se obtiene las ecuaciones:
dx dy dz
dt
= a(x, y, z), dt
= b(x, y, z), dt
= c(x, y, z)
2.- Una de las propiedades importantes que poseen las curvas integrales es: si un punto
de la curva C se encuentra sobre la superficie (gráfico u), entonces toda la curva C
debe pertenecer a la superficie.
2.3.2. Teorema
Dada la E.D.P. semilineal de primer orden
x = f (s)
y = g(s)
z = h(s) donde ”s” es un parametro
Ahora el problema consiste en encontrar una solución de la ecuación(iv) de tal forma que
se verifique la relación h(s) = u(f (x), g(s)) . . . (v)
Nuestro objetivo es obtener una solución local próximo a los valores de x0 = f (s0 ),
y0 = g(s0 ), en una vecindad del punto (x0 , y0 ) (V ((x0 , y0 ), ρ)).
De otro lado el problema de Cauchy también puede ser expresado como un problema
de valor inicial, esto es, en vez que cumpla la relación (v), encontrar una solución de la
E.D.P.(iv), satisfaciendo la condición u(x, 0) = h(x) . . . (vi).
20
2.3.4. Teorema(Existencia y Unicidad)
Dada la E.D.P. semilineal de primer orden
x = f (s)
y = g(s)
z = h(s) y satisface la condición f 0 b − g 0 a 6= 0,
entonces existe una única función z = u(x, y) localmente definida, que es solución de la
ecuación (vii) y satisface la relación h(s) = u(f (s), g(s)).
2.3.5. Observación
1.- Si la E.D.P. de primer orden es de la forma:
2.3.6. Ejemplos
1.- Resolver el problema de Cauchy:
Solución
Primero se debe obtener la curva parametrica C , sobre la cual esta prescrita la condición
inicial.(α2 ).
Esto es:
x = f (s) = s
y = g(s) = 0 (α3 )
z = h(s)
21
se introduce la variable s y t por:
s = S(x, y)
t = T (x, y)
2.3.7. Observación
1.- Del ejemplo 1. si la condición inicial hubiera sido u(x, 0) = h(x) = 2x2 −3, entonces
la única solución del problema de Cauchy hubiera sido:
y 2
u(x, y) = 2 x − a −3
b
2.-Resolver el problema de Cauchy:
Satisfaciendo:
u(x, 0) = 6 − x2 (α2 )
22
Solución
x = f (s) = s
y = g(s) = 0
z = h(s) = 6 − 3s2 ;
y = t; x = 3t + s; s = x − 3t
satisfaciendo
x = s2
y=s
z = 3s
23
Solución
x = f (s) = s2
y = g(s) = s
z = h(s) = 3s;
dividiendo se obtiene
x3 = e6t s6 ; y 2 = e6t s2
luego
x3 e6t s6
2
= 6t 2
= s4 ;
y e s
Luego
s
4 x3
s=
y2
24
Ademas,
y y 1 y
= e3t ; ln( ) = 3t; ln q = t,
s s 3 4 x3
y2
esto es
√ 3/2 1/3 1/2
1 yy y y
t = ln √ = ln 3/4
= ln .
3 4
x 3 x x1/4
Reemplazando en (α4 )
1/4 !
x3 x3/4 y 1/2
1/2
ln( y 1/4 )
u(x, y) = h e x =h
y2 y 1/2 x1/4
x3/4 y 1/2
esto es u(x, y) = 3 = 3x1/2 .
y 1/2 x 1/4
Solución
x = f (s) = s
y = g(s) = 0
z = h(s) = 3s
25
esto es
X(s, t) = −3ts + s
Y (s, t) = 2t
Z(s, t) = −3s despejando s y t se obtiene
X X
X = s(1 − 3t); = s; =s
1 − 3t 1 − 3y
2
t = Y /2
2x
Luego t = y/2; s = ; luego la solución es
2 − 3y
26
3. Ecuaciones Diferenciales Parciales De Segundo
Orden
3.1. Definición
Sea µ : D ⊂ R2 → R una función de clase C 2 , es decir µ ∈ C 2 (D), donde D es un
conjunto abierto y acotado.
Una E.D.P. de segundo orden lineal no homogenea, es de la forma
A(x, y)µxx (x, y) + B(x, y)µxy (x, y) + c(x, y)µyy (x, y)+
+D(x, y)µx (x, y) + E(x, y)µy (x, y)+F (x, y)µ(x, y) = G(x, y)
ε = αx + βy, η = γx + ρy (ii)
Donde α, β, γ y ρ son constantes, las cuales seran determinados sujeto a la condición del
Jacobiano de la transformación.
Donde el jacobiano es diferente de cero; esto es
∂(ε, η) εx εy α β
= = = αρ − γβ 6= 0
∂(x, y) ηx ηy γ ρ
εx = α; εy = β; ηx = γ; ηy = ρ, tambien
28
∂ 2 µ(x, y) ∂ ∂ ∂
2
= α · wε + γ · wη = α wε + γ wη , luego
∂x ∂x ∂x ∂x
∂ 2 µ(x, y)
= α (w εε · ε x + w εη · η x ) + γ w ηε · ε x + w ηη · η x
∂x2
2
∂ µ(x, y)
= α2 wεε + αγwεη + γαwηε + γ 2 wηη
∂x2
Como se esta considerando w ∈ C 2 (D) entonces wεη = wηε . Luego se obtiene
∂ 2 µ(x, y)
= α2 wεε + 2αγwεη + γ 2 wηη (α2 )
∂x2
∂µ ∂w(ε, η)
También = = wε · εy + wη · ηy = βwε + ρwη , luego
∂y ∂y
∂ 2 µ(x, y) ∂ ∂ ∂
= β · wε + ρ · wη = β · wε + ρ · w
∂y 2 ∂y ∂y ∂y η
= β wεε · εy + wεη · ηy + ρ wηε · εy + wηη · ηy
= β 2 wεε + βρwεη + βρwηε + ρ2 wηη
∂ 2 µ(x, y)
= β 2 wεε + 2βρwερ + ρ2 wηη (α3 )
∂y 2
donde:
a = Aα2 + 2αβB + Cβ 2
b = αγA + B(αρ + βγ) + βρC (iv)
c = γ 2 A + 2γρB + ρ2 C
ademas se verifica:
De otro lado, puesto que las constantes α, β, γ y ρ estan a nuestra disposición, se puede
elegir de tal forma que uno o dos de los términos de segundo orden de (iii) se puede anular.
29
3.2.1. Formas Canónicas De La E.D.P. De Segundo Orden
Considerando a B 2 − AC de (v) como el discriminante de la E.D.P. (i), donde A, B
y C son los coeficientes de las derivadas de segundo orden de la ecuación(i), se tienen los
siguientes tipos de E.D.P.
a) E.D.P. De Tipo Hiperbólico
3.2.2. Ejemplo
(1) Dado la E.D.P. µxx (x, y) − 4µxy (x, y) = 0, averiguar que tipo de E.D.P. es
Solución
Solución
Solución
30
3.2.3. Observación
1.- Si los coeficientes de la E.D.P. A, B y C fueran funciones en las variables x e y y
estubieran definidas en cierto dominio D ⊂ R2 , entonce se tiene lo siguiente:
3.2.4. Ejemplo
Dado la E.D.P.: xµxx − 2xµxy + yµyy = 0
Solución
Se observa
Luego
31
3.3. Determinación De Las Formas Canónicas
De la sección anterior, la E.D.P. de segundo orden a estudiar es
ε = αx + βy, η = γx + ρy (vi)
Se desea que a = c = 0 y b 6= 0; veamos algunos casos que se puede dar para los coefi-
cientes de la E.D.P. (i).
Primer Caso: Si A = C = 0
Segundo Caso: Si A 6= 0
32
α
puesto que a y c deben ser iguales a cero, se elige α, β, γ y ρ de tal forma que y
β
γ
sean las raı́ces de la ecuación cuadrática:
ρ
Como el discriminante de esta ecuación es:4 = B 2 − AC > 0, entonces las raı́ces reales
de (vii) son números reales y diferentes.
√ √
−B + B 2 − AC −B − B 2 − AC
Sean las raı́ces: m1 = , m2 =
A A
Luego igualando con las raices ya asumidas se tiene:
√ √
α −B + B 2 − AC γ B − B 2 − AC
= ; =
β A ρ A
√ √
α = −B + B 2 − AC; γ = −B − B 2 − AC
Luego se puede elegir:
β = A; ρ=A
Con esta eleccion de los coeficientes de (iv) se garantiza que: a = 0, c = 0 y el Jacobiano
αρ − βγ 6= 0.
Luego para este segundo caso, las nuevas variables de transformación queda especificado:
√
ε = (−B + B 2 − AC)x + Ay
√
η = (−B − B 2 − AC)x + Ay (x)
wεη + · · · = 0
Tercer Caso: Si A = 0, C 6= 0
0 = 2Bαβ + Cβ 2
0 = 2Bγρ + Cρ2 (xii)
33
Ademas de (v) se sabe: (B 2 − AC)(αρ − βγ)2 = b2 − ac.
Luego como A = 0 y ”C” debe ser cero entonces se obtiene B 2 (αρ − βγ)2 = b2 . Lue-
go α y γ no pueden ser ceros en forma simultanea, esto es; α 6= 0 y γ 6= 0. Dividiendo la
primera ecuación de (xii) por α2 y la segunda ecuación por γ 2 se obtiene:
2 2
β β ρ ρ
2B +C = 0 ∧ 2B +C =0
α α γ γ
β ρ
Consideremos las raı́ces de estas ecuaciones cuadráticas por y , esto es, son
α γ
raı́ces de la ecuación:
−B + |B| −B − |B|
m1 = , m2 =
C C
Luego igualando con las raı́ces supuestas se obtiene
β −B + |B| ρ −B − |B|
= , =
α C γ C
Luego, eligiendo los coeficientes:
β = −B + |B|; ρ = −B − |B|
α = C; γ = C
ε = Cx + (−B + |B|)y
η = Cx + (−B − |B|)y
wεη + · · · = 0
34
Primer Caso: Si A 6= 0
De la condición de tipo parabólico, para obtener su forma canónica se debe cumplir que
b = 0 y c = 0; esto es:
α = 1, γ = B
β = 0, ρ = −A
ε=x
η = Bx − Ay
wεε + · · · = 0
Segundo Caso: Si C 6= 0
Ejercicio.
[Link] Ejemplo Dada la E.D.P.: µxx (x, y) + 6µxy (x, y) − 16µyy (x, y) = 0 (i)
Halle la forma canónica de la E.D.P.
Solución
35
E.D.P. es de tipo hiperbólico. Ademas como A 6= 0, empleamos el segundo caso de la
E.D.P. de tipo hiperbólico. Luego las variables de transformación son:
√
ε = −B + B 2 − AC x + Ay
√
η = −B − B 2 − AC x + Ay
ε = 2x + y
reemplazando, se obtiene ; con estas nuevas variables se debe transformar
η = −8x + y
la E.D.P.(i).
εx = 2; εy = 1
Consideremos: µ(x, y) = w(ε, η), ademas
ηx = −8, ηy = 1
Luego derivando y aplicando la regla de la cadena se obtiene
∂
µx = w = wε · εx + wη · ηx = 2wε − 8wη
∂x
∂ ∂ ∂
µxx = 2wε − 8wη = 2 wε − 8 wη
∂x ∂x ∂x
= 2 wεε · εx + wεη · ηx − 8 wηε · εx + wηη · ηx
µxx = 4wεε + 64wηη − 32wεη (α1 )
∂ ∂
µxy = (µx ) = 2wε − 8wη
∂y ∂y
= 2 wεε · εy + wεη · ηy − 8 wηε · εy · ηy
µxy = 2wεε − 8wηη − 6wηε (α2 )
∂
µy = w = w ε · εy + w η · η y = w ε + w η
∂y
∂
= wε + wη = wεε · εy + wεη · ηy + wηε · εy + wηη · ηy
∂y
µyy = wεε + 2wεη + wηη (α3 )
[Link] Ejemplo Dado la E.D.P.: µxx (x, y) + 4µxy (x, y) + 4µyy (x, y) = 0 (i)
Halle la forma canónica de la E.D.P.
Solución:
ε = x; η = Bx − Ay es decir, ε = x, η = 2x − y
36
con esta variable de transformación, se va reducir la E.D.P.(i) a su forma canónica. con-
ε = 1 εy = 0
sideremos µ(x, y) = w(ε, η), donde x .
ηx = 2 ηy = −1
derivando y aplicando la regla de la cadena, se obtiene:
∂
µx = w = wε · εx + wη · ηx = wε + 2wη
∂x
∂
µxx = wε + 2wη = wεε · εx + wεη · ηx + 2 wηε · εx + wηη · ηx
∂x
µxx = wεε + 4wηη + 4wεη (α1 )
∂
µy = w = wε · εy + wη · ηy = −wη
∂y
∂ ∂u ∂
µxy = = wε + 2wη
∂y ∂x ∂y
µxy = wεε · εy + wεη · ηy + 2 wηε · εy + wηη · ηy
µxy = −wεη − 2wηη (α2 )
∂
µyy = −wη = − wηε · εy + wηη · ηy = wηη (α3 )
∂y
reemplazando (α1 ), (α2 ) y (α3 ) en la E.D.P.(i) se obtiene:
wεε + 4wεη + wηη + 4(−wεη − 2wηη ) + 4wηη = 0
Luego se obtiene wεε = 0, el cual es la forma canónica.
[Link] Ejemplo Dada la E.D.P.: µxx (x, y) − 2µxy (x, y) + 5µyy (x, y) = 0 (i)
reducir a su forma canónica.
Solución
37
con estas nuevas variables de transformación, la E.D.P.(i) se reducirá a su forma canónica.
ε = 1, εy = 1
Consideremos µ(x, y) = w(ε, η), donde x derivando y empleando la regla
ηx = 2, ηy = 0
de la cadena, se obtiene:
∂
µx = w = wε · εx + wη · ηx = wε + 2wη
∂x
∂
µxx = wε + 2wη = wεε + 4wεη + 4wηη
∂x
∂
µy = w = wε · εy + wη · ηy = wεε
∂y
∂
µyy = (w ) = wεε · εy + wεη · ηy = wεε
∂y ε
∂
µxy = wε + 2wη = wεε · εy + wεη · ηy + 2 wηε · εy + wηη · ηy
∂y
µxy = wεε + 2wεη
L(µ) = 0
∂ ∂ ∂2 ∂2
donde Dx = , Dy = , Dx Dy = , Dx2 =
∂x ∂y ∂x∂y ∂x2
El operador L es un operador diferenciable de orden dos. Algunos operadores de orden
dos se pueden expresar de la forma L = L1 · L2 , donde L1 y L2 son operadores lineales y
38
diferenciables de orden uno.
Se verifica que el producto de operadores lineales de primer orden es conmutativo, esto
es, L1 · L2 = L2 · L1 . En efecto:
3.4.2. Ejemplo
1.- Sea el operador L = Dx2 + 2Dx − 3I
averiguar si L es factorable.
Solución
L = Li · L2 ; L es factorable.
39
3.4.3. Observación
1. Dado una E.D.P. lineal de segundo orden
L(µ) = 0
L1 (µ) = 0 (α2 )
En efecto:
Sea z = v1 + v2 P.D. L(z) = 0
∴ L(z) = 0
3.4.4. Ejemplo
Dado una E.D.P. lineal de segundo orden y homogéneo
aµxx (x, y) + bµxy (x, y) + cµyy (x, y) + dµx (x, y) + eµy (x, y) + f µ(x, y) = 0 (i)
40
a = a1 a2 , b = a1 b 2 + b 1 a2
ademas, .
c = b1 b2 , d = a1 c2 + c1 a2 , e = b1 c2 + b2 c1
Estos operadores lineales L1 y L2 generan dos E.D.P.s homogéneas de primer orden,
esto es
3.4.5. Observación
1. si a = a1 · a2 = 0 donde a1 = 0 o a2 = 0, entonces la solución general de la E.D.P.(ii)
contiene los términos
c1 c2
− x − x
a1 a2
f1 (x) = e o f2 (x) = e donde f1 y f2 son funciones arbitrarias.
3.4.6. Ejemplo
Halle la solución de la E.D.P.
41
Solución
donde L1 = Dx + Dy y L2 = Dx + Dy − I
Sea L(µ) = g(x), donde g(x) = sen x, puesto que g solo depende de la variable x, en-
tonces emplearemos el método de coeficientes indeterminados que se estudia en las E.D.O.
De aquı́ se obtiene: A − B = 1; −A − B = 0
donde el operador L es
L1 (L1 µ) = 0 (ii)
consideremos
L1 (µ) = v (iii)
L1 (v) = 0 (iv)
a1 vx +b1 vy + c1 v = 0 (v)
c1
− x
a1
v(x, y) = g1 (b1 x − a1 y)e (vi)
c1
x
−
L1 (µ) = g1 (b1 x − a1 y)e a1
esto es la E.D.P.
c1
− x
a1
a1 µx + b1 µy + c1 µ = g1 (b1 x − a1 y)e (vii)
43
para resolver (vii) consideremos la nueva variable de transformación ε = b1 x − a1 y,
luego, considerando u(x, y) = w(x, ε), la E.D.P.(vii) se reduce a la siguiente ecuación
c1
x
−
a1 wx (x, ε) + c1 w(x, ε) = g1 e a1
esto es
c1
− x
c1 g (ε) a1
wx (x, ε) + w(x, ε) = 1 e (viii)
a1 a1
resolviendo la ecuación (viii) como una E.D.O. se obtiene
Z
c1
c1 Z x c 1
− a1 −
a e a 1
w(x, ε) = e 1 g1 (ε)e dx + g2 ε
a1
operando se obtiene
c1
x
−
a1 x
w(x, ε) = e g (ε) + g2 (ε)
a1 1
retornando a la variable original se obtiene
c1
−
a1 x x
u(x, y) = e g (b x − a1 y) + g2 (b1 x − a1 y)
a1 1 1
3.4.8. Ejemplo
Resolver la E.D.P.: µxx (x, y) − 2µxy (x, y) + µyy (x, y) = 0 . . . (i)
Solución
Luego L = L1 · L2 ; donde L1 = Dx − Dy
44
esto es L1 (L1 (µ)) = 0 . . . (ii)
reemplazando se obtiene
3.4.9. Observación
1. Consideremos una E.D.P. de segundo orden lineal y no homogénea: L(µ) = G . . . (i)
donde L es un operador factorable: L = L1 · L2 esto es L1 6= L2 o L1 = L2 .
Para obtener la solución particular de (i) se produce en forma similar que se realizo
para resolver la E.D.P. L21 = 0, donde L = L1 · L1 .
L1 (v) = G (iv)
45
3.4.10. Ejemplo
Resolver la E.D.P.:
µxx (x, y) − 2µxy (x, y) − 3µyy (x, y) − µx (x, y) − µy (x, y) = x · y (i)
solución
factorizando se obtiene
L = (Dx − 3Dy − I)(Dx − Dy ) = L1 · L2
donde L1 = Dx − 3Dy − I y L2 = Dx + Dy
vx − 3vy − v = x · y (v)
Para resolver (v) introducimos la variable ε = −3x − y, luego considerando
v(x, y) = w(x, ε) , la E.D.P.(v) se reduce a
wx (x, ε) − w(x, ε) = x(−3x − ε) = −3x2 − εx (vi)
resolviendo esta E.D.P. como una E.D.O. se obtiene
Z Z
dx Z − dx
w(x, ε) = e · e
(−3x2 − εx)dx
; luego
Z
x −x 2
w(x, ε) = e e (−3x − εx)dx ; operando se obtiene
46
w(x, ε) = 3x2 + (6 + ε)(x + 1), esto es la solución de (vi).
Retornando a la variable original se obtiene
v(x, y) = 3x2 + (6 − 3x − y)(x − 1) (vii)
reemplazando (vii) en (iv) se obtiene
x3 x2
Luego w(x, ε) = + (8 − ε) + x(ε − 6).
3 2
Retornando a la variable original se obtiene
x3 x2
µ(x, y) = + (8 − x + y) + x(x − y − 6)
3 2
esto representa la solución particular de la E.D.P.(i).
47
donde X e Y son funciones que dependen de las variables x e y respectivamente.
Reemplazando (ii) en la E.D.P.(i) se obtiene una proporción formada por dos razones to-
talmente independientes. Esta razón, adjuntado la constante de proporcionalidad genera
dos E.D.O., una en la variable x y la otra en la variable y.
Resolviendo independientemente cada E.D.O. obtenida, se obtienen sus soluciones respec-
tivas que reemplazados en la expresión (ii), nos permite obtener una solución general de
la E.D.P.(i). Este método de se conoce como Método de Fourier o también el de Variables
Separables.
3.4.12. Ejemplo
Resolver la E.D.P. empleando el método de Fourier
µx (x, y) + 3µy (x, y) = 0 (i)
Solución
48
3.4.13. Observación
1. Si al problema anterior se le añade una condición de frontera, esto es
µ(0, y) = 4e−2y − 3e−6y , entonces debemos hallar la solución particular del problema
µx + 3µy = 0 satisfaciendo
(vii) −2y −6y
µ(0, y) = 4e − 3e
3.4.14. Ejemplo
∂µ(x, t) ∂ 2 µ(x, t)
Resolver el problema: −2 = 0 · · · (i)
∂t ∂x2
satisfaciendo la condición de frontera:
µ(0, t) = 0 3πx
∀t ≥ 0 · · · (ii) y µ(x, 0) = 50 sen + 20 sen πx − 5 sen 4πx · · · (iii)
µ(10, t) = 0 2
Solución
49
µ(x, t) = X(x) · T (t) · · · (iv); derivando se obtiene: X(x) · T (t) − 2X 00 (x) · T (t) = 0 de
aquı́ se obtiene
T 0 (t) X 00 (x)
= =c
2T (t) X(x)
si se resuelve estas E.D.O., se obtiene solución trivial, esto es X(x) ≡ 0; por lo tanto si
reemplazamos en (iii) se obtiene µ(x, t) = 0, esto es la solución trivial, y esta solución no
nos interesa puesto que la E.D.P. es de segundo orden, entonces para no tener solución
trivial, se debe elegir a la constante de proporcionalidad por c = −λ2 ; luego reemplazando
en (v) se obtiene las E.D.O.
50
resolviendo la otra E.D.O.(vii) se obtiene:
T (t) = Be−2λt , donde B es una constante
Luego la solución de la E.D.P.(i) satisfaciendo (ii) es de la forma
nπ !2
−2 t nπx
u(x, t) = Be 10 · sen (xi)
10
Ahora, emplearemos la condición de frontera (iii) en (xi)
nπx nπx
µ(x, 0) = Be0 sen = B sen
10 10
3πx
= 50 sen + 20 sen(2πx) − 5 sen(4πx) (xii)
2
sin embargo, no podemos encontrar un solo par de constantes n y B que satisfagan la
igualdad (xii) puesto que las constantes B y n son arbitrarias, consideremos tres soluciones
de la forma
n1 π !2
n2 π 2
−2 t
−2 t n πx n πx
10
µ1 (x, t) = B1 e 10 · sen 1
, µ2 (x, t) = B2 e · sen 2
10 10
n3 π 2
−2 t
n πx
10
µ3 (x, t) = B3 e · sen 3
10
aplicando el principio de superposición se obtiene:
n21 π 2 t n22 π 2 t
n πx n πx
µ(x, y) =B1 e 50 · sen 1 + B2 e 50 · sen 2 +
10 10
n23 π 2 t
n πx
+ B3 e 50 · sen 3 (xiii)
10
el cual también es solución.
51
4. Series De Fourier
4.1. Introducción
Las series de Fourier son series con términos de senos y cosenos, surgen por la impor-
tancia de representar funciones periódicas generales mediante una serie infinita.
Las series de Fourier constituyen una herramienta importante para la resolución de ecua-
ciones diferenciales. Las series de Fourier, en cierto sentido son mas generales que las
series de Taylor, puesto que muchas funciones periódicas son discontinuas y sin embargo
pueden desarrollarse en series de Fourier, lo cual no es posible desarrollarlo en series de
Taylor.
Joseph Fourier (1768-1830) utilizo las series de Fourier en su obra principal (Theorie
Analytique del calor 1822) en la que desarrolla la teorı́a de la conducción del calor.
Prueba
Prueba
T
Sea M 6= 0, entonces f M x + = f (M x + T ) = f (M x), pues M x ∈ R y
M
52
T
f es periódica de periodo T ; luego f M x + = f (M x); luego f (M x) es
M
T
periódica de periodo .
M
4.2.3. Ejemplos
(1) Si f (x) = cos 3x, halle el periodo de f
Solución
Sea g(x) = cos x, se sabe que su periodo es T = 2π, luego si g(3x) = cos 3x,
entonces el periodo es
T 2π
=
3 3
x
(2) Si f (x) = sen halle el periodo de f
4
Solución
x x
Sea g(x) = sen x, su periodo es T = 2π, luego g = sen , su periodo es
4 4
2π
= 8π
1
4
x x
(3) Sea la función f (x) = sen + cos , halle el periodo de f .
3 4
Solución
x 2π
Sea g(x) = sen , entonces su periodo es = 6π, luego
3 1/3
x x
g(x) = sen = sen + 6π (i)
3 3
x
sea h(x) = cos , entonces su periodo es 8π, luego
4
x x
h(x) = cos = cos + 8π (ii)
4 4
Luego, buscaremos un periodo único para (i) y (ii)
x x
sabemos g(x) = sen + 6π = sen + n(6π) , n∈N
x3 x3
h(x) = cos + 8π = cos + m(8π) , m∈N
4 4
53
Luego si n = 4 y m = 3 se obtiene
x x x x
g(x) = sen = sen + 24π y h(x) = cos = cos + 24π (iii)
3 3 4 4
Luego reemplazando (iii) en la función original
x x x x
f (x) = sen + cos = sen + 24π + cos + 24π
3 4 3 4
Luego f es periódica, de periodo T = 24π.
x x
(4) Si f (x) = sen x + sen + cos , halle el periodo de f .
2 3
Operando se obtiene que el periodo de f es T = 12π
4.2.4. Observación
1. Si f es periódica, de periodo T , ademas f es integrable en [0, T ], entonces se cumple
Z T Z T +C
f (x) dx = f (x) dx donde C es una constante real
0 C
4.2.6. Observación
1. Si y = f (x) es una función par, entonces su gráfico es simétrico respecto al eje y.
Z a Z a
2. Si f es una función par e integrable en [−a, a] entonces f (x) dx = 2 f (x) dx.
−a 0
Z a
3. Si f es una función impar e integrable en [−a, a] entonces f (x) dx = 0
−a
4.2.7. Ejemplos
Z π
(1) Sea f (x) = sen x; entonces sen x dx = 0, f es impar.
−π
Z π Z π
2 2
(2) Sea f (x) = cos x; entonces f (x) dx = f (x) dx, f es par.
−π
2
0
Z 1
5
(3) Sea f (x) = x , entonces x5 dx = 0, f es impar.
−1
54
4.2.8. Observación
1. Si f es par y g es impar entonces f · g es impar.
4.2.10. Observación
1. Toda función continua en [a, b] es seccionalmente continua en [a, b].
4.2.11. Ejemplos
(1) Sea f (x) = x2 , averiguar si f es seccionalmente continua en todo intervalo acotado.
Solución
Solución
55
Grafico de f
y
4.3.2. Definición
Sea S = {ϕ0 , ϕ1 , ϕ2 , · · · } un conjunto de funciones reales definidos en [a, b]. Diremos
que el conjunto S es ortogonal en [a, b] si:
Z b
ϕi · ϕj dx = 0, ∀ i 6= j
a
4.3.3. Ejemplo
Dado el conjunto de funciones, definidos en −L ≤ x ≤ L;
nπx nπx
{1; sen , cos ; n = 1, 2, 3, . . . }. Probar que dicho conjunto es un sistema ortogonal
L L
en [−L, L]
Solución
Se debe probar que para cada par de funciones diferentes del conjunto, la integral del
producto de ambas es cero en [−L, L], veamos los siguientes casos:
56
Z L
nπx
Caso 1: 1 · sen
dx
−L L
nπx
Se sabe que la función sen es una función impar ∀n ∈ N.
L
Luego
Z L
nπx
sen dx = 0
−L L
Z L
nπx
Caso 2: 1 · cos
dx
−L L
nπx
Se sabe que la función cos es una función par, luego
L
Z L Z L L
nπx nπx L nπx
cos dx = 2 cos dx = 2 sen =0
−L L 0 L nπx L 0
Z L nπx nπx
Caso 3: sen · cos dx
−L L L
Sabemos que el producto de una función par con una función impar es función impar,
luego.
nπx nπx
sen · cos es función impar luego, de propiedad se obtiene
L L
Z L nπx nπx
sen · cos dx = 0.
−L L L
Z L nπx mπx
Caso 4: cos · cos dx, n 6= m
−L L L
nπx mπx 1 πx πx
Sabemos: cos · cos = cos (m + n) + cos (m − n)
L L 2 L L
Luego integrando se obtiene:
Z L
1 L
nπx mπx Z
πx πx
cos · cos dx = cos (m + n) + cos (m − n) dx
−L L L 2 −L L L
πx πx
L sen(m + n) sen(n − m)
= L + L =0
2π m+n n−m
Z L nπx mπx
Caso 5: sen · sen dx
−L L L
(Ejercicio)
57
4.4. Series de Fourier
4.4.1. Definición
Sea f una función periódica de periodo T = 2L y ademas f es seccionalmente continua
en [−L, L]. La serie de Fourier asociada a la función f esta definido por
∞
a X nπx nπx
f (x) = 0 + an cos + bn sen
2 n=1
L L
4.4.2. Lema
Sea f una función periódica de periodo T = 2L y ademas f es seccionalmente continua
en [−L, L]; si la serie de Fourier de f es:
∞
a0 X nπx nπx
f (x) = + an cos + bn sen (i)
2 n=1
L L
1 L
Z
(I) a0 = f (x) dx.
L −L
1 L
Z nπx
(II) an = f (x) · cos dx ∀n = 1, 2, . . . .
L −L L
1 L
Z nπx
(III) bn = f (x) · sen dx ∀n = 1, 2, . . .
L −L L
Prueba
58
nπx nπx
empleando el hecho de que la función sen es impar, y que la función cos
L L
es par, se obtiene de (ii)
Z L ∞
X Z L nπx ∞
X Z L nπx
f (x) dx = a0 L + 2an cos dx + bn sen dx
−L n=1 −L L −L L
|n=1 {z }
0
nπx L
Z L ∞ sen
L
X
f (x) dx = a0 L + 2an
nπ
+0
−L n=1
L 0
Z L ∞
X
Luego f (x) dx = a0 L + 2an · (0) + 0, esto es
−L n=1
Z L
1
a0 = f (x) dx
L −L
mπx
Prueba de II: Multiplicando la serie (i) por cos , e integrando desde −L hasta
L
L se obtiene:
Z L mπx Z L mπx X∞ Z L nπx mπx
f (x) · cos dx = a0 · cos dx + an cos cos dx
−L L −L L n=1 −L L L
X∞ Z L nπx mπx
+ bn sen cos dx · · · (iii)
n=1 −L L L
nπx mπx
Ademas, sabemos que la función: sen cos es una función impar, entonces
L L
Z L nπx mπx Z L mπx
sen cos dx = 2a0 cos dx
−L L L 0 L
∞
X Z L nπx mπx
+ bn sen cos dx
n=1 −L L L
Z L mπx
ademas se sabe que dx = 0 · · · (iv)
cos
L 0
Luego de la igualdad de (iv) queda de la forma siguiente
Z L mπx ∞
X Z L nπx mπx
f (x) cos dx = 2an cos cos dx · · · (v)
−L L n=1 −L L L
59
Caso (1). Si m 6= n
integrando se obtiene:
πx πx L
Z L mπx ∞ sen(m + n) sen(n − m)
L + L
X
f (x) cos dx = an · π π
−L L n=1 (m + n) (n − m)
L L 0
Z L mπx
De aquı́ se obtiene: f (x) cos dx = 0
−L L
Luego no se obtiene nada.
Caso 2: m = n
integrando
Z L mπx Z L
2πx
f (x) cos dx = 2am cos + 1 dx
−L L −L L
L
2mπx
sen
L
= am
+ x = am · L
2mπ
L
0
Z L
1 mπx
Luego am = f (x) cos dx ∀n = 1, 2, 3, . . .
L −L L
4.4.3. Ejemplos
(1) Halle la serie de Fourier de la función
0; −π < x < 0
f (x) = ; f (x + 2π) = f (x)
π − x; 0≤x<π
60
Solución
de otro lado
Z π Z π Z π
1 1
an = f (x) · cos(nx) dx = 0 · cos(nx) dx + (π − x) · cos(nx) dx
π −π π −π −π
Z π
a0 = (π − x) · cos(nx) dx; integrando por partes se obtiene
−π
π π π
1 sen(nx) sen(nx) cos(nx)
an = π· −x +
π n 0 n 0 n2 0
(−1)n − 1
1 cos(nπ) − 1 1
a0 = 0−0+ = (cos(nπ) − 1) =
π n2 πn2 πn2
(−1)n − 1
Luego an = . . . (iii)
πn2
1 π
Z Z 0 Z π
1
En forma análoga bn = f (x) sen(nx) dx = 0 dx + (π − x) sen(nx) dx
Z π π −π π −π 0
1
Luego bn = (π − x) sen(nx) dx = . . . (iv)
0 n
Finalmente reemplazando (ii), (iii) y (iv) en (i) se obtiene
∞
π X ((−1)n − 1)
sen(nx)
f (x) = + cos(nx) +
4 n=1 n2 · π n
La serie de Fourier de f .
61
Solución
1 π 1 π+7π 1 8π
Z Z Z
a0 = f (x) dx = f (x) dx = f (x) dx
π −π π −π+7π π 6π
También
Z π Z 8π
1 1
an = f (x) · cos(nx) dx = f (x) · cos(nx) dx
π −π π 6π
Z 7π Z 8π
1
= (x − 6π) cos(x) dx + (8π − x) cos(nx) dx
π 6π 7π
(cos(nx) − 1)
=2
πn2
Ademas
Z 8π
1
bn = f (x) · sen(nx) dx
π 6π
Z 7π Z 8π
1
= (x − 6π) sen(x) dx + (8π − x) sen(nx) dx
π 6π 7π
operando se obtiene bn = 0 ∀n = 1, 2, . . .
62
4.4.4. Operaciones con Series de Fourier
Sean f y g funciones periódicas de periodo T = 2L que admite serie de Fourier, donde:
∞
a0 X nπx nπx
f (x) = + an · cos + bn · sen
2 n=1
L L
∞
ab X nπx nπx
g(x) = 0 + abn · cos + bbn · sen
2 n=1
L L
entonces se verifica.
(i) (f ± g) admite serie de Fourier en [−L, L], donde sus coeficientes de Fourier son:
(ii) La función (c · f ) admite serie de fourier en [−L, L] donde sus coeficientes de Fourier
son:
Si f una función suave por tramos en h−L, Li, entonces para cada x ∈ h−L, Li, la
serie de Fourier converge a f (x); excepto en los puntos de discontinuidad esto es, si
x0 ∈ h−L, Li es un punto de discontinuidad, entonces la serie de fourier de f converge a
1
f (x+ ) + f (x 0 ) +
2 0
63
4.5.3. Observación
1. si f es periódica de periodo de T = 2L y continua en [−L, L] entonces la serie de
Fourier converge a f (x) ∀x ∈ [−L, L]; esto es;
∞
a0 X nπx nπx
+ an · cos + bn · sen → f (x)
2 n=1
L L
3. f (x+
0
) = lı́m+ f (x), f (x−
0
) = lı́m− f (x)
x→x0 x→x0
4.5.4. Ejemplos
0; −π ≤ x < 0
1.- Sea f (x) = ; f (x + 2π) = f (x)
π − x; 0≤x≤π
Halle la serie de Fourier de f y halle la función donde la serie converge.
Solución
Por lo tanto f es una función suave por tramos puesto que f es discontinua en x = 0
entonces
f (0+ ) + f (0− ) 0+π π
= =
2 2 2
64
Luego la serie de Fourier de f converge a
∞ 0, −π < x < 0
π X ((−1)n − 1) sen(nx) π
+ · cos(nx) + = , x=0
4 n=1 n2 π n 2
π − x, 0<x<π
Sin necesidad de hallar la serie de Fourier de f , halle la función donde la serie converge.
Solución
Derivando se obtiene.
π
5 · cos x; −2π < x < −
2
π
f 0 (x) = 2x; − <x<2
2
−7 · sen x; 2<x<π
4; π < x < 2π
π
Luego f 0 es discontinua en γ1 = − , γ2 = 2, γ3 = π; sin embargo, en cada uno de los
2
puntos de discontinuidad existe los limites laterales de f 0 . Luego f 0 es continua por tramos
en h−π, πi; por lo tanto f es suave por tramos.
65
f (2+ ) + f (2− ) 4 + 7 · cos(2)
=
2 2
f (π + ) + f (π − ) 4π − 7
=
2 2
Luego la serie de Fourier f converge a la función siguiente
π
5 · sen x; −2π < x < −
2
π 2 + 20 π
x=−
;
2 2
π
x2 ; −
<x<2
2
∞ nπx
a0 X nπx
+ an · cos + bn · sen = 4 + 7 · cos(2)
2 n=1
2π 2π ; x=2
2
7 · cos x; 2<x<π
4π − 7
; x=π
2
4x; π < x < 2π
4.5.5. Observación
1. Del ejemplo (1), la serie de Fourier de f nos puede ser útil para hallar ciertas series
números reales.
Sabemos
∞ 0, −π < x < 0
π X ((−1)n − 1) sen(nx) π
+ · cos(nx) + = , x=0
4 n=1 n2 π n 2
π − x, 0<x<π
Ahora, si x = 0 se obtiene
∞
π X ((−1)n − 1) sen(n · 0) π
+ 2
· cos(n · 0) + =
4 n=1 nπ n 2
∞
π X ((−1)n − 1) π
+ =
4 n=1 n2 π 2
∞
X ((−1)n − 1) π
2
=
n=1
nπ 4
66
luego el valor de la serie de números reales es:
∞
X ((−1)n − 1) π
2
=
n=1
nπ 4
4.5.7. Ejemplo
Sea la función f definida por:
−x
e ; −2 ≤ x < 1
f (x) = ; f (x + 4) = f (x)
−2x2 ; 1≤x≤2
Sin necesidad de hallar la serie de Fourier, halle la función donde converge la serie
∀x ∈ [−2, 2]
Solución
Se observa que f es periódica de periodo T = 4 y continua por tramos en [−2, 2], donde
f es discontinua en x = 1. Veamos si f es suave por tramos
−e−x ;
0 −2 < x < 1
f (x) =
−4x; 1<x<2
Luego f 0 es continua por tramos en h−2, 2i y discontinua en x = 1, donde existe los limites
laterales de f 0 . Luego f es suave por tramos.
Ahora hallaremos los limites laterales en el punto discontinuidad x = 1 y en los extremos
x = ±2.
Se tiene
f (1+ ) + f (1− ) e−1 − 2 0
= ; f (−2+ ) = −e2 ; f 0 (2− ) = −8
2 2
67
Ademas
f (−2+ ) + f (2− ) e2 − 8
=
2 2
Luego empleando el segundo teorema de convergencia, la serie de Fourier: converge a la
función
−x
e ; −2 < x < 1
−1
∞
e − 2
a0 X nπx nπx ; x=1
+ an · cos + bn · sen = 22
2 L L −2x ; 1<x<2
n=1 2
e − 8
; x = ±2
2
4.5.8. Observación
1.- Consideremos la función periódica de periodo T = 2π
−1; −π ≤ x < 0
f (x) =
1; 0≤x≤π
∞
X 2
Esta función admite una serie de Fourier y esta dado por [1 − (−1)n ] sen(nx),
n=1
nπ
ademas esta serie converge a
−1; −π < x < 0
∞
X 2
0; x=0
[1 − (−1)n ] sen(nx) = (α1 )
nπ 1; 0 <x<π
n=1
0; x = ±π
68
4.6. Derivación e Integración de Series de Fourier
4.6.1. Teorema.(Derivación)
Sea f una función periódica de periodo T = 2L, continua en [−L, L] cuya serie de
Furier esta dado por
∞
a0 X nπx nπx
f (x) = + an · cos + bn · sen
2 n=1
L L
Si f 0 es continua por tramos en [−L, L], ademas existe f 00 (x) ∀x ∈ h−L, Li y se verifica
que f (−L) = f (L), entonces
∞
X nπ nπx nπ nπx
f 0 (x) = bn · · cos − an · · sen
n=1
L L L L
4.6.2. Teorema.(Integración)
Sea f ua función periódica de periodo T = 2L y seccionalmente continua en [−L, L];
entonces para x ∈ [−L, L], la integral de la serie de Fourier esta dado por
Z x Z x ∞ Z x
a0 X nπt nπt
f (t) dt = dt + an · sen + bn · sen dt
−L −L 2 n=1 −L
L L
esto es
Z x ∞
a0 LX1n nπx h nπx io
f (t) dt = (x + L) + an · sen − bn · cos − (−1)n
−L 2 π n=1 n L L
4.6.3. Ejemplos
1. Sea f (x) = |x|, −π ≤ x ≤ π; f (x + 2π) = f (x)
Solución
−x; −π ≤ x < 0
f (x) = ; f es continua en todo R; luego f es continua en
x; 0≤x≤π
[−π, π]; L = π.
69
Luego la serie de Fourier de f es
∞
π X ((−1)n − 1)
f (x) = + 2 · cos(nx), ∀x ∈ [−π, π]
2 n=1 πn2
Solución
• Si x ∈ h−π, 0i:
Z x x
x
r2
Z
1
−r dr = − = − x2 − π 2
f (r) dr =
−π −π 2 −π 2
70
• Si x ∈ [0, π]:
Z x Z 0 Z x Z 0 Z x
f (r) dr = f (r) dr + f (r) dr = −r dr + r dr
−π −π 0 −π 0
Luego
x
0 x
r2 r2 π 2 x2
Z
f (r) dr = − + = +
−π 2 −π 2 0 2 2
Fácilmente se obtiene
2 2
Z x π −x ;
−π ≤ x < 0
f (r) dr = 2
2 2 (ii)
−π
π + x
; 0≤x≤π
2
Luego de (i) y (ii) se obtiene
2
∞
π − x2
n ; −π ≤ x < 0
π 2 X [(−1) − 1]
2
(x − π) + · sen(nx) = 2 2 (iii)
2 π n=1
n3 π +x ;
0≤x≤π
2
4.6.4. Observación
π
1.- Si reemplazamos x = en la expresión (iii), se obtiene
2
∞
π π 2 X [(−1)n − 1] π π 2 + π/4
( − π) + · sen(n · ) =
2 2 π n=1 n3 2 2
∞
3π 2 2 X [(−1)n − 1] π 5π 2
+ · sen(n · ) =
4 π n=1 n3 2 8
∞
2 X [(−1)n − 1] π 5π 2 6π 2 π2
· sen(n · ) = − = −
π n=1 n3 2 8 8 8
∞
X [(−1)n − 1] π π3
luego · sen(n · )=
n=1
n3 2 16
71
Sin embargo en muchas aplicaciones existe la necesidad de expresar mediante series de
Fourier definidas en ciertos intervalos de la forma h0, Li, esto es f : h0, Li → R, donde f
no es periódica ni definida en todo R.
Si deseamos representar mediante serie de Fourier a este tipo de funciones, primeramente
de debe extender a la función f , con una extension periódica, donde la función extendida
sea una función par o impar.
Luego, la serie de Fourier de cosenos para esta función esta dado por:
∞
a0 X nπx
+ an · cos = fp (x) = f (x); 0 < x < L
2 n=1
L
Z L Z L
2 2 nπx
donde a0 = f (x) dx; an = f (x) · cos dx
L 0 L 0 L
Luego, la serie de Fourier de senos para esta función esta dado por:
∞
X nπx
bn · sen = fi (x) = f (x); 0 < x < L
n=1
L
Z L
2 nπx
donde bn = f (x) · sen dx
L 0 L
4.7.3. Ejemplos
(1) Sea f (x) = x2 ; 0 < x < π
72
Solución
Solución
73
a) Halle la extension par de f en [−L, L] y grafique
-L -L/2 L/2 L
74
4k h nπ
n
i
an = · 2 · cos − (−1) − 1 , luego se obtiene
n2 π 2 L
∞
X 4k h nπ
n
i nπx
f (x) = · 2 · cos − (−1) − 1 · cos
n=1
n2 π 2 L L
c) La extensión impar periódica de f es
2k
− (L + x); −L < x < −L/2
L
2k
L x; −L/2 < x < 0
fi (x) = ; fi (x + 2L) = fi (x)
2kx
; 0 < x < L/2
L
2k
(L − x); L/2 < x < L
L
graficando fi se obtiene
-L -L/2 L/2 L
-k
n=1
n2 π 2 L
75
4.8. Serie de Fourier de forma Compleja
4.8.1. Teorema.(Serie de Fourier Compleja)
Sea f una función compleja de periodo T = 2L y secionalmente continua en [−L, L],
si la serie de Fourier de f es
∞
a0 X nπx nπx
f (x) = + an · cos + bn · sen
2 n=1
L L
La serie de Fourier compleja esta dado por:
nπx nπx
+∞ i 1
Z L −i
cn · e L , donde cn = L
X
f (x) = f (x) · e
n=−∞
2L −L
Prueba
76
ordenando adecuadamente se obtiene
nπx nπx
a
∞
(an − ibn ) i (a + ibn ) −i
e L + n L
X
f (x) = 0 + e
2 n=1
2 2
Luego
nπx nπx
a0 X
∞ i −i
f (x) = + cn e L + k n e L (vi)
2 n=1
nπx
1
Z L −i 1
Z L
Ademas cn = f (x) · e L dx, entonces c0 = f (x) dx
2L −L 2L −L
1 L
Z
a
Luego 2c0 = f (x) dx = a0 , esto es c0 = 0
L −L 2
reemplazando en (vii) se obtiene
nπx
+∞ i
cn · e L ;
X
f (x) = c0 + (viii)
n=−∞
n6=0
nπx
1
Z L −i
donde cn = f (x) · e L dx
2L −L
77
4.8.2. Observación
1.- A partir de la serie de Fourier compleja, se puede hallar la serie de Fourier real de f .
an = cn + kn (x)
luego de (x), (xi) y (xii) se obtienen los coeficientes de Fourier de la serie de Fourier real
de f .
4.8.3. Ejemplo
Sea la función definida por:
−1; −π < x < 0
f (x) = ; f (x + 2π) = f (x)
1; 0<x<π
Solución
78
luego esta serie compleja converge a
1; 0<x<π
+∞
X i
−1; −π < x < 0
[(−1)n − 1] einx =
nπ 0; x=0
n=−∞
0; x = ±π
4.8.4. Ejemplo
Dado la función f (x) = x2 ; −π < x < π; f (x + 2π) = f (x)
• Halle la serie de Fourier complejo de f .
• A partir de la serie de Fourier complejo, halla la serie de Fourier real.
Solución
79
Luego se ha obtenido los coeficientes de Fourier de la serie de Fourier real de f ; esto es:
∞
2π 2 X 4(−1)n
+ , 2
· cos(nx) = f (x); −π ≤ x ≤ π
3 n=1
n
4.8.5. Observación
1.- Si x = 0, de (iii) se obtiene:
∞ ∞
2π 2 X 4(−1)n 2π 2 X 4(−1)n
+ , · cos(0) = 0; + , =0
3 n=1
n2 3 n=1
n2
∞ ∞
X 4(−1)n−1 2π 2 X 2(−1)n−1 π2
luego: = ; =
n=1
n2 3 n=1 n2 3
2.- Si x = 1: de (iii) se obtiene:
∞
2π 2 X 4(−1)n
+ , · cos(n) = 1, Luego se obtiene
3 n=1
n2
∞
X 2(−1)n−1 2π 2 − 3
, · cos(n) =
n=1
n2 6
80
5. La integral y Transformada de Fourier
5.1. La integral de Fourier
Sabemos de la sección anterior que si una función f es periódica entonces dicha función
se puede expresar mediante serie de Fourier.
Supongamos que una función f : R → R es no periódica, deseamos estudiar que condi-
ciones debe satisfacer dicha función para poder expresarlo mediante una serie o mediante
una integral. Sea fL : h−L, Li → R una función periódica de periodo T = 2L y continua
por tramos en h−L, Li; luego su representación en serie de Fourier es
∞
a0 X nπx nπx
fL (x) = + an · cos + bn · sen
2 n=1
L L
nπ
Si consideramos wn = se tiene la serie
L
∞
a X
fL (x) = 0 + an · cos wn x + bn · sen wn x (i)
2 n=1
1 L 1 L 1 L
Z Z Z
donde a0 = f (v) dv; an = f (v) cos(wn v) y bn = f (v) sen(wn v)
L −L L L −L L L −L L
reemplazando en (i) se obtiene:
Z L ∞
1 L
Z
1 X
fL (x) = f (v) dv+ f (v) cos wn v dv cos wn x
2L −L L n=1
L −L L
Z L
1
+ f (v) sen wn v dv sen(wn x)
L −L L
esto es
Z L ∞ Z L
1 1X
fL (x) = fL (v) dv + cos wn x · fL (v) cos wn v dv
2L −L L n=1 −L
Z L
+ sen wn x · fL (v) sen wn v dv (ii)
−L
(n + 1)π nπ
Ahora, hacemos 4w = wn+1 − wn , esto es 4w = −
L L
π 4w 1
Luego 4w = y = . . . (iii)
L π L
Reemplazamos (iii) en (ii)
Z L ∞ Z L
1 4w X
fL (x) = f (v) dv + cos wn x · fL (v) cos wn v dv
2L −L L π n=1 −L
Z L
+ sen wn x · fL (v) sen wn v dv
−L
81
Luego
Z L ∞ Z L
1 1X
fL (x) = fL (v) dv + cos wn x · 4w · fL (v) cos wn v dv
2L −L π n=1 −L
Z L
+ sen wn x · 4w · fL (v) sen wn v dv (iv)
−L
La representación de (iv) es valida para cualquier L fijo arbitrariamente grande pero finito.
lı́m fL (x) = f (x) donde f : R → R, es una función real no periodica definida en todo R.
L→∞
π
Ademas, desde que 4w = → 0, entonces parece plausible que la serie infinita (v) se
L
convierta en una integral desde 0 hasta +∞, esto es
1 +∞
Z Z +∞ Z +∞
f (x) = cos wx · f (v) cos wv dv + sen wx · f (v) · sen wv dv dw
π 0 −∞ −∞
Luego
Z +∞ Z +∞
1
f (x) = cos wx · f (v) cos wv dv
0 π −∞
Z +∞
1
+ sen wx · f (v) · sen wv dv dw (vi)
π −∞
82
Ahora se introduce la notación
1 +∞ 1 +∞
Z Z
A(w) = f (v) · cos wv dv; B(w) = f (v) · sen wv dv
π −∞ π −∞
Luego se concluye, que la función no periódica se puede expresar mediante una integral,
denominada la integral de Fourier de f .
5.1.1. Teorema
Sea f : R → R una función real; si f es continua por tramos en todo intervalo acotado
de R, f admite derivadas laterales en cada punto x ∈ R y si ademas f es absolutamente
integrable en todo R, entonces f se puede representar por las integral de Fourier (vii),
ademas la serie de Fourier converge a:
f (x); Si f es continua en x
Z +∞
[cos wx · A(w) + sen wx · B(w)] dw = + −
0 f (x0 ) + f (x0 ) ;
Si f es discontinua en x0
2
5.1.2. Ejemplo
1.- Dado la función
1; |x| < 1
f (x) = Halle la integral de Fourier
0; |x| ≥ 1
Solución
83
Ahora se introduce la notación
1 +∞ 1 +∞
Z Z
A(w) = f (v) · cos wv dv; B(w) = f (v) · sen wv dv
π −∞ π −∞
Luego se concluye, que la función no periódica se puede expresar mediante una integral,
denominada la integral de Fourier de f .
5.1.1. Teorema
Sea f : R → R una función real; si f es continua por tramos en todo intervalo acotado
de R, f admite derivadas laterales en cada punto x ∈ R y si ademas f es absolutamente
integrable en todo R, entonces f se puede representar por las integral de Fourier (vii),
ademas la serie de Fourier converge a:
f (x); Si f es continua en x
Z +∞
[cos wx · A(w) + sen wx · B(w)] dw = + −
0 f (x0 ) + f (x0 ) ;
Si f es discontinua en x0
2
5.1.2. Ejemplo
1.- Dado la función
1; |x| < 1
f (x) = Halle la integral de Fourier
0; |x| ≥ 1
Solución
Redefiniendo f se tiene
1; −1 < x < 1
f (x) = 0; 1<x<∞ ; es claro que f es continua
0; −∞ < x < −1
83
Ahora hallaremos A(w) y b(w)
Z −1
1 +∞
Z Z 1 Z +∞
1
A(w) = f (v) cos wv dv = 0 · cos wv dv + cos wv dv + 0 · cos wv dv
π −∞ π −∞ −1 1
1 1
Z
2 sen x
A(w) = cos wv dv = ; también
π −1 πw
1 +∞ 1 1
Z Z
B(w) = f (v) · cos wv dv = sen wv dv = 0
π −∞ π −1
Luego la integral de Fourier de f esta dado por:
Z +∞
2 · sen x · cos wx
f (x) = dw; luego la integral converge a
0 πw
∀x ∈ R − ±1
(
Z +∞
2 · sen w · cos wx f (x) :
dw = 1
0 πw x = ±1
2
1;
−1 < x < 1
= 0; x > 1 ∨ x < −1 (i)
1
; x = ±1
2
Luego de (i) se obtiene
π
; |x| < 1
2
+∞
sen w · cos wx
Z
= 0; |x| > 1
0 w
π;
x = ±1
4
5.1.3. Observación
1. De la integral de Fourier de f obtenida en el ejemplo anterior, si x = 0, se obtiene
Z ∞
sen · cos w(0) π
dw = esto es
0 w 2
Z ∞
sen w π
dw =
0 w 2
∞
sen w · cos w
Z
π
2. Si x = 1; se obtiene dw =
0 w 4
5.1.4. Ejemplo
e−x ;
x>0
Dado la función f (x) =
ex ; x<0
84
Halle la integral de Fourier de f
Solución
5.1.5. Observación
1. Del ejemplo anterior
Z ∞ Z ∞
2 cos 3w cos 5w π
si x = 5, entonces dw = e−5 . Luego dw = e−5
π 0 1 + w2 0 1+w 2 2
2 R ∞ cos w(−1)
2. Si x = −1, entonces dw = e−1
π 0 1 + w2
Z ∞
cos wx π
Luego 2
dw = e−1
0 1+w 2
3. La integral de Fourier de una función f , también se emplea para calcular ciertas
integrales impropias de funciones en R.
85
entonces la integral de Fourier de f es
Z ∞
f (x) = [A(w) cos wx + B(w) sen wx] dw (i)
0
donde
1 +∞
Z
A(w) = f (v) · cos wv dv (ii)
π −∞
1 +∞
Z
B(w) = f (v) · sen wv dv (iii)
π −∞
1 +∞
Z Z a
1
A(w) = f (v) · cos wv dv = lı́m f (v) cos wv dv = 0
π −∞ π a→∞ −a
1 +∞
Z Z a
1
B(w) = f (v) · sen wv dv = lı́m f (v) sen wv dv = 0
π −∞ π a→∞ −a
5.1.7. Ejemplo
Sea la función f (x) = e−x , x > 0
86
Solución
reemplazando se obtiene
+∞ +∞
2w · sen wx w · sen wx
Z Z
2
f (x) = 2
dw = dw
0 π(1 + w ) π 0 1 + w2
87
Z +∞ Z +∞
1 1
donde A(w) = f (v) · cos wv dv; B(w) = f (v) · sen wv dv
π −∞ π −∞
luego
1 +∞
Z Z +∞
f (x) = f (v) {cos wv · cos wx + sen wv · sen wx} dv dw
π 0 −∞
1 +∞
Z Z +∞
f (x) = f (v) · {cos(wv − wx)} dv dw (ii)
π 0 −∞
R +∞
consideremos G(w) = −∞ f (v) {cos(w + wx)} dv
se verifica que G es una función par, esto es, G(w) = G(−w), luego reemplazando en
(ii) se obtiene
1 +∞
Z
f (x) = G(w) dw (iii)
π 0
Ademas de la propiedad de función par se sabe
Z +∞ Z +∞
1 +∞
Z Z +∞
G(w) dw = 2 G(w) dw; esto es; G(w) dw = G(w) dw
−∞ 0 2 −∞ 0
En efecto
Z +∞
Sea H(w) = f (v) · sen(wx − wv) dv probaremos que H(w) es una función impar
−∞
Z +∞
H(−w) = f (v) · sen(−wx + wv) dv
−∞
Z +∞
=− f (v) · sen(wx − wv) dv = −H(w)
−∞
88
Luego
Z +∞ Z +∞
1
f (v) · sen(wx − wv) dv dw = 0
2π −∞ −∞
también
Z +∞ Z +∞
i
f (v) · sen(wx − wv) dv dw = 0 (vi)
2π −∞ −∞
Luego se obtiene
Z +∞ Z +∞
1
f (x) = f (v) {cos(wx − wv) + i sen(wx − wv)} dv dw
2π −∞ −∞
Luego se obtiene
Z +∞ Z +∞
1 iw(x−v)
f (x) = f (v) · e dv dw (viii)
2π −∞ −∞
5.2.2. Ejemplo
1.- Sea f (x) = e−a|x| , ∀x ∈ R, halle la integral compleja de Fourier de f
Solución
−ax
e ; x≥0
f (x) = ax Se observa que f es continua en todo R y es absolutamen-
e ; x<0
te convergente.
Luego f admite la integral compleja de Fourier
Z +∞
1
f (x) = cw eiwx dw donde
2π −∞
Z +∞ Z 0 Z +∞
−iwv −iwv
cw = f (v) · e dv = av
e ·e dv + e−av · e−iwv dv
−∞ −∞ 0
Z 0 Z +∞
2
cw = ev(a−iv) + e−v(a−iw) dv = 2
−∞ 0 a + w2
89
Luego la integral compleja de Fourier es
Z +∞
1 2
f (x) = eiwx dw
2π −∞ a + w2
2
5.3.1. Definición.(Existencia)
Sea f : R → R una función definida en R. Diremos que f admite transformada de
Fourier si cumple lo siguiente:
5.3.2. Observación
1.- Otra forma de denotar a la transformada de una función f es fˆ(w) o F (w).
90
5.3.3. Ejemplo
−x
e ; x>0
Sea f (x) = ; Halle F (f )
0; x<0
Solución
1
Luego F (f ) = fˆ(w) =
1 + iw
Ademas la transformada inversa de Fourier es
−x
e ;
x>0
F (F (w)) = f (x) =
−1 0; x<0
1;
x=0
2
5.3.4. Ejemplo
Si f (x) = e−|x| , x ∈ R, Halle F (f )
Solución
ex ; x<0
Se sabe f (x) = ; f es continua en todo R y es absolutamente integrable;
e−x ; x≥0
luego
Z +∞ Z 0 Z +∞
F (f ) = F (w) = f (x)e = −iwx x
e · e−iwx dx + e−x · e−iwx dx
−∞ −∞ 0
Z 0 Z +∞
1 1 2
= ex(1−iw) dx + e−x(1+iw) dx = + =
−∞ 0 1 − iw 1 + iw 1 + w2
2
Luego F (f ) =
1 + w2
Ademas como f es continua en todo R, entonces la transformada inversa de Fourier
es:
f (x) = e−|x| ∀x ∈ R.
5.3.5. Ejemplo
2
Sea f (x) = e−x halle: F (w)
91
Solución
x2 + y 2 = r 2
x = r · cos θ
; 0 ≤ θ ≤ ; J(r, θ) = r; Jacobiano de la transformada
y = r · sen θ
r≥0
esto es
Z 2π Z ∞ Z +∞ Z 2π
2 −r2 −r2
(I) = e · r dr dθ = e · r dθ dr
0 00 0
Z +∞
1 2 2 +∞
= 2π − e−r d(−r2 ) = −πe−r =π
2 0 0
√
Luego (I)2 = π, entonces I = π; reemplazando en (ii)
w2 √
F (w) = e 4 · π
5.3.6. Ejercicio
w2 √
−at2 e− 4a · π
Si f (t) = e , probar F (w) = √ ; a 6= 0
a
92
5.4. Transformada de Fourier de Seno y Coseno
5.4.1. Definición.(Transformad de Fourier de Coseno)
Si f es una función real par, continua por tramos y absolutamente integrable entonces,
la integral de Fourier de Coseno es
Z ∞
f (x) = A(w) cos wx dw (i)
0
Z +∞
2
donde A(w) = f (v) · cos wv dv (ii)
π 0
√
2
Denotemos por A(w) = √ · fˆc (w) (iii)
π
Luego igualando (ii) y (iii) se obtiene
√
2 +∞
Z
2 ˆ
√ · fc (w) = f (v) · cos wv dv; Luego
π π 0
√ Z +∞
2
fˆc (w) = √ f (v) · cos wv dv o
π 0
√ Z +∞
ˆ 2
fc (w) = √ f (x) · cos wx dx (iv)
π 0
2 R +∞
donde B(w) = f (v) · sen wv dv (vii)
π 0
r
2 ˆ
Denotemos por: B(w) = · fs (w) (viii)
π
93
Luego reemplazando (vii) en (viii) se obtiene
r
2 +∞
Z
2 ˆ
f (v) · sen wv dv = · fs (w)
π 0 π
r Z −∞ √ Z +∞
2 2
fˆs (w) = f (v) · sen wv dv o fˆs (w) = √ · f (x) · sen wx dx (ix)
π 0 π 0
5.4.3. Ejemplos
k; 0<x<a
1.- Sea f (x) = , halle: fˆc (w) y fˆs (w)
0; x>a
Solución
r Z ∞ r Z a Z +∞
2 2
fˆc (w) = f (x) · cos wx dx = k · cos wx dx + 0 · cos wx dx
π 0 π 0 a
r Z a r
ˆ 2 k 2
fc (w) = k · cos wx dx + 0 = · sen aw
π 0 w π
√
k 2
esto es fˆs (w) = √ · sen aw
w π
También
r Z ∞ r Z a Z +∞
2 2
fˆs (w) = f (x) · sen wx dx = k · sen wx dx + 0 · sen wx dx
π 0 π 0 a
r Z a √
ˆ 2 2 1 − cos wa
fs (w) = k · sen wx dx = √ k
π 0 π w
esto es
√
ˆ 2 1 − cos wa
fs (w) = √ k
π w
94
2.- Si f (x) = e−x ; hallar fˆc (w)
Solución
√ Z +∞ r Z +∞
ˆ 2 2
fc (w) = √ f (x) · cos wx dx = · e−x · cos wx dx
π 0 π 0
integrando por partes se obtiene
r
2
π
fˆc (w) =
1 + w2
95
5.4.5. Observación
1.- Si f 0 es continua y absolutamente integrable en R+ y ademas f 00 es continua por
tramos en todo intervalo acotado, en donde lı́m f 0 (x) = 0, entonces:
x→∞
r
2 0
(I) Fc (f 00 ) = −w2 · Fc (f ) − · f (0)
π
r
2
(II) Fc (f ) = −w · Fs (f ) +
00 2
· w · f (0)
π
Demostración
5.4.6. Ejemplo
Empleando la propiedad de la transformada de la derivada hallar Fc (f ) si f (x) = e−ax ,
a>0
Solución
96
Sin embargo en muchas aplicaciones existe la necesidad de expresar mediante series de
Fourier definidas en ciertos intervalos de la forma h0, Li, esto es f : h0, Li → R, donde f
no es periódica ni definida en todo R.
Si deseamos representar mediante serie de Fourier a este tipo de funciones, primeramente
de debe extender a la función f , con una extension periódica, donde la función extendida
sea una función par o impar.
Luego, la serie de Fourier de cosenos para esta función esta dado por:
∞
a0 X nπx
+ an · cos = fp (x) = f (x); 0 < x < L
2 n=1
L
Z L Z L
2 2 nπx
donde a0 = f (x) dx; an = f (x) · cos dx
L 0 L 0 L
Luego, la serie de Fourier de senos para esta función esta dado por:
∞
X nπx
bn · sen = fi (x) = f (x); 0 < x < L
n=1
L
Z L
2 nπx
donde bn = f (x) · sen dx
L 0 L
4.7.3. Ejemplos
(1) Sea f (x) = x2 ; 0 < x < π
72
Solución
Solución
73
a) Halle la extension par de f en [−L, L] y grafique
-L -L/2 L/2 L
74
4k h nπ
n
i
an = · 2 · cos − (−1) − 1 , luego se obtiene
n2 π 2 L
∞
X 4k h nπ
n
i nπx
f (x) = · 2 · cos − (−1) − 1 · cos
n=1
n2 π 2 L L
c) La extensión impar periódica de f es
2k
− (L + x); −L < x < −L/2
L
2k
L x; −L/2 < x < 0
fi (x) = ; fi (x + 2L) = fi (x)
2kx
; 0 < x < L/2
L
2k
(L − x); L/2 < x < L
L
graficando fi se obtiene
-L -L/2 L/2 L
-k
n=1
n2 π 2 L
75
4.8. Serie de Fourier de forma Compleja
4.8.1. Teorema.(Serie de Fourier Compleja)
Sea f una función compleja de periodo T = 2L y secionalmente continua en [−L, L], si la
serie de Fourier de f es
∞
a0 X nπx nπx
f (x) = + an · cos + bn · sen
2 n=1
L L
La serie de Fourier compleja esta dado por:
nπx nπx
+∞ i 1
Z L −i
cn · e L , donde cn = L
X
f (x) = f (x) · e
n=−∞
2L −L
Prueba
76
ordenando adecuadamente se obtiene
nπx nπx
a
∞
(an − ibn ) i (a + ibn ) −i
e L + n L
X
f (x) = 0 + e
2 n=1
2 2
Luego
nπx nπx
a0 X
∞ i −i
f (x) = + cn e L + k n e L (vi)
2 n=1
nπx
1
Z L −i 1
Z L
Ademas cn = f (x) · e L dx, entonces c0 = f (x) dx
2L −L 2L −L
1 L
Z
a
Luego 2c0 = f (x) dx = a0 , esto es c0 = 0
L −L 2
reemplazando en (vii) se obtiene
nπx
+∞ i
cn · e L ;
X
f (x) = c0 + (viii)
n=−∞
n6=0
nπx
1
Z L −i
donde cn = f (x) · e L dx
2L −L
77
4.8.2. Observación
1.- A partir de la serie de Fourier compleja, se puede hallar la serie de Fourier real de f .
an = cn + kn (x)
luego de (x), (xi) y (xii) se obtienen los coeficientes de Fourier de la serie de Fourier real
de f .
4.8.3. Ejemplo
Sea la función definida por:
−1; −π < x < 0
f (x) = ; f (x + 2π) = f (x)
1; 0<x<π
Solución
78
luego esta serie compleja converge a
1; 0<x<π
+∞
X i
−1; −π < x < 0
[(−1)n − 1] einx =
nπ 0; x=0
n=−∞
0; x = ±π
4.8.4. Ejemplo
Dado la función f (x) = x2 ; −π < x < π; f (x + 2π) = f (x)
• Halle la serie de Fourier complejo de f .
• A partir de la serie de Fourier complejo, halla la serie de Fourier real.
Solución
79
Luego se ha obtenido los coeficientes de Fourier de la serie de Fourier real de f ; esto es:
∞
2π 2 X 4(−1)n
+ , 2
· cos(nx) = f (x); −π ≤ x ≤ π
3 n=1
n
4.8.5. Observación
1.- Si x = 0, de (iii) se obtiene:
∞ ∞
2π 2 X 4(−1)n 2π 2 X 4(−1)n
+ , · cos(0) = 0; + , =0
3 n=1
n2 3 n=1
n2
∞ ∞
X 4(−1)n−1 2π 2 X 2(−1)n−1 π2
luego: = ; =
n=1
n2 3 n=1 n2 3
2.- Si x = 1: de (iii) se obtiene:
∞
2π 2 X 4(−1)n
+ , · cos(n) = 1, Luego se obtiene
3 n=1
n2
∞
X 2(−1)n−1 2π 2 − 3
, · cos(n) =
n=1
n2 6
80
5. La integral y Transformada de Fourier
5.1. La integral de Fourier
Sabemos de la sección anterior que si una función f es periódica entonces dicha función
se puede expresar mediante serie de Fourier.
Supongamos que una función f : R → R es no periódica, deseamos estudiar que condi-
ciones debe satisfacer dicha función para poder expresarlo mediante una serie o mediante
una integral. Sea fL : h−L, Li → R una función periódica de periodo T = 2L y continua
por tramos en h−L, Li; luego su representación en serie de Fourier es
∞
a0 X nπx nπx
fL (x) = + an · cos + bn · sen
2 n=1
L L
nπ
Si consideramos wn = se tiene la serie
L
∞
a X
fL (x) = 0 + an · cos wn x + bn · sen wn x (i)
2 n=1
1 L 1 L 1 L
Z Z Z
donde a0 = f (v) dv; an = f (v) cos(wn v) y bn = f (v) sen(wn v)
L −L L L −L L L −L L
reemplazando en (i) se obtiene:
Z L ∞
1 L
Z
1 X
fL (x) = f (v) dv+ f (v) cos wn v dv cos wn x
2L −L L n=1
L −L L
Z L
1
+ f (v) sen wn v dv sen(wn x)
L −L L
esto es
Z L ∞ Z L
1 1X
fL (x) = fL (v) dv + cos wn x · fL (v) cos wn v dv
2L −L L n=1 −L
Z L
+ sen wn x · fL (v) sen wn v dv (ii)
−L
(n + 1)π nπ
Ahora, hacemos 4w = wn+1 − wn , esto es 4w = −
L L
π 4w 1
Luego 4w = y = . . . (iii)
L π L
Reemplazamos (iii) en (ii)
Z L ∞ Z L
1 4w X
fL (x) = f (v) dv + cos wn x · fL (v) cos wn v dv
2L −L L π n=1 −L
Z L
+ sen wn x · fL (v) sen wn v dv
−L
81
Luego
Z L ∞ Z L
1 1X
fL (x) = fL (v) dv + cos wn x · 4w · fL (v) cos wn v dv
2L −L π n=1 −L
Z L
+ sen wn x · 4w · fL (v) sen wn v dv (iv)
−L
La representación de (iv) es valida para cualquier L fijo arbitrariamente grande pero finito.
lı́m fL (x) = f (x) donde f : R → R, es una función real no periodica definida en todo R.
L→∞
π
Ademas, desde que 4w = → 0, entonces parece plausible que la serie infinita (v) se
L
convierta en una integral desde 0 hasta +∞, esto es
1 +∞
Z Z +∞ Z +∞
f (x) = cos wx · f (v) cos wv dv + sen wx · f (v) · sen wv dv dw
π 0 −∞ −∞
Luego
Z +∞ Z +∞
1
f (x) = cos wx · f (v) cos wv dv
0 π −∞
Z +∞
1
+ sen wx · f (v) · sen wv dv dw (vi)
π −∞
82
Ahora se introduce la notación
1 +∞ 1 +∞
Z Z
A(w) = f (v) · cos wv dv; B(w) = f (v) · sen wv dv
π −∞ π −∞
Luego se concluye, que la función no periódica se puede expresar mediante una integral,
denominada la integral de Fourier de f .
5.1.1. Teorema
Sea f : R → R una función real; si f es continua por tramos en todo intervalo acotado
de R, f admite derivadas laterales en cada punto x ∈ R y si ademas f es absolutamente
integrable en todo R, entonces f se puede representar por las integral de Fourier (vii),
ademas la integral de Fourier converge a:
f (x); Si f es continua en x
Z +∞
[cos wx · A(w) + sen wx · B(w)] dw = + −
0 f (x0 ) + f (x0 ) ;
Si f es discontinua en x0
2
5.1.2. Ejemplo
1.- Dado la función
1; |x| < 1
f (x) = Halle la integral de Fourier
0; |x| ≥ 1
Solución
Redefiniendo f se tiene
1; −1 < x < 1
f (x) = 0; 1<x<∞ ; es claro que f es continua
0; −∞ < x < −1
83
Ahora hallaremos A(w) y B(w)
Z −1
1 +∞
Z Z 1 Z +∞
1
A(w) = f (v) cos wv dv = 0 · cos wv dv + cos wv dv + 0 · cos wv dv
π −∞ π −∞ −1 1
1 1
Z
2 sen x
A(w) = cos wv dv = ;
π −1 πw
1 +∞ 1 1
Z Z
B(w) = f (v) · cos wv dv = sen wv dv = 0
π −∞ π −1
Luego la integral de Fourier de f esta dado por:
Z +∞
2 · sen x · cos wx
f (x) = dw; luego la integral converge a
0 πw
∀x ∈ R − ±1
(
Z +∞
2 · sen w · cos wx f (x) :
dw = 1
0 πw x = ±1
2
1;
−1 < x < 1
= 0; x > 1 ∨ x < −1 (i)
1
; x = ±1
2
Luego de (i) se obtiene
π
; |x| < 1
2
+∞
sen w · cos wx
Z
= 0; |x| > 1
0 w
π;
x = ±1
4
5.1.3. Observación
1. De la integral de Fourier de f obtenida en el ejemplo anterior, si x = 0, se obtiene
Z ∞
sen · cos w(0) π
dw = esto es
0 w 2
Z ∞
sen w π
dw =
0 w 2
∞
sen w · cos w
Z
π
2. Si x = 1; se obtiene dw =
0 w 4
5.1.4. Ejemplo
e−x ;
x>0
Dado la función f (x) =
ex ; x<0
84
Halle la integral de Fourier de f
Solución
5.1.5. Observación
1. Del ejemplo anterior
Z ∞ Z ∞
2 cos 5w cos 5w π −5
si x = 5, entonces dw = e−5 . Luego dw = e
π 0 1 + w2 0 1 + w2 2
2 R ∞ cos w(−1)
2. Si x = −1, entonces dw = e−1
π 0 1 + w2
Z ∞
cos w π
Luego 2
dw = e−1
0 1+w 2
3. La integral de Fourier de una función f , también se emplea para calcular ciertas
integrales impropias de funciones en R.
85
entonces la integral de Fourier de f es
Z ∞
f (x) = [A(w) cos wx + B(w) sen wx] dw (i)
0
donde
1 +∞
Z
A(w) = f (v) · cos wv dv (ii)
π −∞
1 +∞
Z
B(w) = f (v) · sen wv dv (iii)
π −∞
1 +∞
Z Z a
1
A(w) = f (v) · cos wv dv = lı́m f (v) cos wv dv = 0
π −∞ π a→∞ −a
1 +∞
Z Z a
1
B(w) = f (v) · sen wv dv = lı́m f (v) sen wv dv = 0
π −∞ π a→∞ −a
5.1.7. Ejemplo
Sea la función f (x) = e−x , x > 0
86
Solución
reemplazando se obtiene
+∞ +∞
2w · sen wx w · sen wx
Z Z
2
f (x) = 2
dw = dw
0 π(1 + w ) π 0 1 + w2
87
Z +∞ Z +∞
1 1
donde A(w) = f (v) · cos wv dv; B(w) = f (v) · sen wv dv
π −∞ π −∞
luego
1 +∞
Z Z +∞
f (x) = f (v) {cos wv · cos wx + sen wv · sen wx} dv dw
π 0 −∞
1 +∞
Z Z +∞
f (x) = f (v) · {cos(wv − wx)} dv dw (ii)
π 0 −∞
R +∞
consideremos G(w) = −∞ f (v) {cos(w + wx)} dv
se verifica que G es una función par, esto es, G(w) = G(−w), luego reemplazando en
(ii) se obtiene
1 +∞
Z
f (x) = G(w) dw (iii)
π 0
Ademas de la propiedad de función par se sabe
Z +∞ Z +∞
1 +∞
Z Z +∞
G(w) dw = 2 G(w) dw; esto es; G(w) dw = G(w) dw
−∞ 0 2 −∞ 0
En efecto
Z +∞
Sea H(w) = f (v) · sen(wx − wv) dv probaremos que H(w) es una función impar
−∞
Z +∞
H(−w) = f (v) · sen(−wx + wv) dv
−∞
Z +∞
=− f (v) · sen(wx − wv) dv = −H(w)
−∞
88
Luego
Z +∞ Z +∞
1
f (v) · sen(wx − wv) dv dw = 0
2π −∞ −∞
también
Z +∞ Z +∞
i
f (v) · sen(wx − wv) dv dw = 0 (vi)
2π −∞ −∞
Luego se obtiene
Z +∞ Z +∞
1
f (x) = f (v) {cos(wx − wv) + i sen(wx − wv)} dv dw
2π −∞ −∞
Luego se obtiene
Z +∞ Z +∞
1 iw(x−v)
f (x) = f (v) · e dv dw (viii)
2π −∞ −∞
5.2.2. Ejemplo
1.- Sea f (x) = e−a|x| , ∀x ∈ R, halle la integral compleja de Fourier de f
Solución
−ax
e ; x≥0
f (x) = ax Se observa que f es continua en todo R y es absolutamen-
e ; x<0
te convergente.
Luego f admite la integral compleja de Fourier
Z +∞
1
f (x) = cw eiwx dw donde
2π −∞
Z +∞ Z 0 Z +∞
−iwv −iwv
cw = f (v) · e dv = av
e ·e dv + e−av · e−iwv dv
−∞ −∞ 0
Z 0 Z +∞
2
cw = ev(a−iv) + e−v(a−iw) dv = 2
−∞ 0 a + w2
89
Luego la integral compleja de Fourier es
Z +∞
1 2
f (x) = eiwx dw
2π −∞ a + w2
2
5.3.1. Definición.(Existencia)
Sea f : R → R una función definida en R. Diremos que f admite transformada de Fourier
si cumple lo siguiente:
5.3.2. Observación
1.- Otra forma de denotar a la transformada de una función f es fˆ(w) o F (w).
90
5.3.3. Ejemplo
−x
e ; x>0
Sea f (x) = ; Halle F(f )
0; x<0
Solución
1
Luego F(f ) = fˆ(w) =
1 + iw
Ademas la transformada inversa de Fourier es
−x
e ;
x>0
−1
F (F (w)) = f (x) = 0; x<0
1;
x=0
2
5.3.4. Ejemplo
Si f (x) = e−|x| , x ∈ R, Halle F(f )
Solución
ex ; x<0
Se sabe f (x) = ; f es continua en todo R y es absolutamente integrable;
e−x ; x≥0
luego
Z +∞ Z 0 Z +∞
−iwx −iwx
F(f ) = F (w) = f (x)e = e ·e x
dx + e−x · e−iwx dx
−∞ −∞ 0
Z 0 Z +∞
1 1 2
= e x(1−iw)
dx + e−x(1+iw) dx = + =
−∞ 0 1 − iw 1 + iw 1 + w2
2
Luego F(f ) =
1 + w2
Ademas como f es continua en todo R, entonces la transformada inversa de Fourier
es:
f (x) = e−|x| ∀x ∈ R.
5.3.5. Ejemplo
2
Sea f (x) = e−x halle: F (w)
91
Solución
x2 + y 2 = r 2
x = r · cos θ
; 0 ≤ θ ≤ 2π ; J(r, θ) = r; Jacobiano de la transformada
y = r · sen θ
r≥0
esto es
Z 2π Z ∞ Z +∞ Z 2π
2 −r2 −r2
(I) = e · r dr dθ = e · r dθ dr
0 0 0 0
Z +∞
1 2 2 +∞
= 2π − e−r d(−r2 ) = −πe−r =π
2 0 0
√
Luego (I)2 = π, entonces I = π; reemplazando en (ii)
w2 √
F (w) = e 4 · π
5.3.6. Ejercicio
w2 √
−at2 e− 4a · π
Si f (t) = e , probar F (w) = √ ; a 6= 0
a
92
5.4. Transformada de Fourier de Seno y Coseno
5.4.1. Definición.(Transformad de Fourier de Coseno)
Si f es una función real par, continua por tramos y absolutamente integrable entonces, la
integral de Fourier de Coseno es
Z ∞
f (x) = A(w) cos wx dw (i)
0
Z +∞
2
donde A(w) = f (v) · cos wv dv (ii)
π 0
√
2
Denotemos por A(w) = √ · fˆc (w) (iii)
π
Luego igualando (ii) y (iii) se obtiene
√
2 +∞
Z
2 ˆ
√ · fc (w) = f (v) · cos wv dv; Luego
π π 0
√ Z +∞
2
fˆc (w) = √ f (v) · cos wv dv o
π 0
√ Z +∞
ˆ 2
fc (w) = √ f (x) · cos wx dx (iv)
π 0
2 R +∞
donde B(w) = f (v) · sen wv dv (vii)
π 0
r
2 ˆ
Denotemos por: B(w) = · fs (w) (viii)
π
93
Luego reemplazando (vii) en (viii) se obtiene
r
2 +∞
Z
2 ˆ
f (v) · sen wv dv = · fs (w)
π 0 π
r Z −∞ √ Z +∞
2 2
fˆs (w) = f (v) · sen wv dv o fˆs (w) = √ · f (x) · sen wx dx (ix)
π 0 π 0
5.4.3. Ejemplos
k; 0<x<a
1.- Sea f (x) = , halle: fˆc (w) y fˆs (w)
0; x>a
Solución
r Z ∞ r Z a Z +∞
2 2
fˆc (w) = f (x) · cos wx dx = k · cos wx dx + 0 · cos wx dx
π 0 π 0 a
r Z a r
ˆ 2 k 2
fc (w) = k · cos wx dx + 0 = · sen aw
π 0 w π
√
k 2
esto es fˆc (w) = √ · sen aw
w π
También
r Z ∞ r Z a Z +∞
2 2
fˆs (w) = f (x) · sen wx dx = k · sen wx dx + 0 · sen wx dx
π 0 π 0 a
r Z a √
ˆ 2 2 1 − cos wa
fs (w) = k · sen wx dx = √ k
π 0 π w
esto es
√
ˆ 2 1 − cos wa
fs (w) = √ k
π w
94
2.- Si f (x) = e−x ; hallar fˆc (w)
Solución
√ Z +∞ r Z +∞
ˆ 2 2
fc (w) = √ f (x) · cos wx dx = · e−x · cos wx dx
π 0 π 0
integrando por partes se obtiene
r
2
π
fˆc (w) =
1 + w2
95
5.4.5. Observación
1.- Si f 0 es continua y absolutamente integrable en R+ y ademas f 00 es continua por tramos
en todo intervalo acotado, en donde lı́m f 0 (x) = 0, entonces:
x→∞
r
2 0
(I) Fc (f 00 ) = −w2 · Fc (f ) − · f (0)
π
r
00 2 2
(II) Fs (f ) = −w · Fs (f ) + · w · f (0)
π
Demostración
5.4.6. Ejemplo
Empleando la propiedad de la transformada de la derivada hallar Fc (f ) si f (x) = e−ax ,
a>0
Solución
96
5.4.7. Ejercicio
Si f (x) = 4e−3x+1 ; halle: Fc (f ) y Fs (f )
5.5.1. Propiedad.(Linealidad)
Si F(f ) = F (w), F(g) = G(w), entonces
Prueba.
Z +∞ Z +∞ Z +∞
−iwx −iwx
F(αf ± βg) = (αf (x) ± βg(x)) e dx = αf (x)e dx ± βg(x)e−iwx dx
−∞ −∞ −∞
= αF(f ) ± βF(g) = αF (w) ± βG(w)
5.5.2. Propiedad(Dilatación)
1 w
Si F(f ) = F (w), entonces F (f (ax)) = F , a 6= 0
|a| a
Z +∞
du
Prueba. F(f (ax)) = f (ax)e−iwx dx: hacemos ax = u, dx = .
−∞ a
Si a > 0: si x → +∞, entonces u → +∞; luego
1 +∞
Z Z +∞
−iw( u ) 1 w 1 w
F(f (ax)) = f (u)e a du = f (u)e−i( a )u du = F (i)
a −∞ a ∞ a a
1 +∞
Z w 1 w
w
F (f (ax)) = f (u)e−iw( a ) d = F (ii)
a −∞ a a a
97
5.5.3. Propiedad(Transformada del producto con escalar)
Se sabe que la función escalón unidad esta definido por
1; si t ≥ 0 1; si t ≥ a
u(t − 0) = u(t) = u(t − a) =
0; si t < 0 0; si t < a
1
Si f (x) = e−x u(x), entonces F(f ) =
1 + iw
Prueba.
Z +∞ Z +∞ Z 0 Z +∞
−iwx −x −iwx −x −iwx
F(f ) = f (x)e dx = e u(x)e dx = e 0e dx + e−x e−iwx dx
−∞ −∞ −∞ 0
+∞ −x[1+iw] +∞
Z
e 1
= e−x[1+iw] dx = − =
0 1 + iw 0 1 + iw
5.5.4. Observación
e−a(1+iw)
1. Si f (x) = e−x u(x − a); entonces F(f ) =
1 + iw
1
2. Si f (x) = e−ax u(t); entonces F(f ) = ;a>0
a + iw
5.5.5. Ejemplos
1. Si f (t) = e−|at| , hallar F(f )
Solución. Sabemos del ejemplo anterior si:
2
f (t) = e−|t| , entonces F(f ) = = F (w)
1 + w2
ahora aplicando la propiedad de dilatación.
1 w
Si f (at) = e−|at| ; entonces F(f (at)) = F (i)
|a| a
2 w 2
Sabemos F (w) = 2
; entonces F = w 2 (ii)
1+w a
1+
a
de (ii) en (i) se obtiene
1 2 2|a|2 2|a|
F(f (at)) = w 2 = = 2
|a|
2 2
|a| (a + w ) a + w2
1+
a
2|a|
∴ F e−|at| =
a2 + w2
98
2. Si f (t) = e−3t u(t − 4); hallar F(f )
Solución.
Z +∞ Z 4 Z +∞
−3t −iwt −3t −iwt
F(f (t)) = e u(t − 4)e dt = e 0e dt + e−3t u(t − 4)e−iwt dt
−∞ −∞ 4
+∞ +∞ −t(3+iw) +∞
Z Z
−3t −iwt −t(3+iw) e
= e e dt = e dt = −
4 4 3 + iw 4
−4(3+iw)
e
=
3 + iw
e−b(a+iw)
3. Si f (t) = e−at u(t − b), entonces F(f (t)) =
a + iw
u = t − a, du = dt y t = u + a
reemplazando se obtiene
Z +∞ Z +∞
−iw(u+a)
F(f (t − a)) = e du = f (u)e−iwu e−iwa du
−∞ −∞
Z +∞
= e−iwa f (u)e−iwu du = e−iwa F (w)
−∞
99
5.5.8. Ejemplos
2
1. Si f (t) = e−t eit halle F(f (t))
2
Solución. Sea h(t) = e−t , entonces f (t) = h(t)eit (i)
de propiedad (5.5.6)
2 w2 √
pero se sabe de lo estudiado F(h(t)) = F(e−t ) = e− 4 π, reemplazando en (ii) se obtiene
(w−1)2 √
F(f (t)) = H(w − 1) = e− 4 π
2. Sea f (t) = g(t) sen at; donde g(t) admite transformada de Fourier; Halle F(f (t))
Solución. Se sabe que eait = cos at+i sen at, eait = cos at−i sen at restando ello se obtiene
eita − e−iat
sen at =
2i
eita − e−iat
Luego, f (t) = g(t)
2i
aplicando transformada
1 1 −ait 1 1
ait
= F g(t)eiat − F g(t)e−iat
F(f (t)) = F g(t)e − g(t)e
2i 2i 2i 2i
e − e−i3t
i3t
h(t) i3t h(t) −i3t
F(f (t)) = F (h(t) cos 3t) = F h(t) =F e + e
2 2 2
1 1
= F h(t)ei3t + F h(t)e−i3t
2 2
100
de propiedad (5.5.6)
1 1
F(f (t)) = H(w − 3) + H(w + 3) (iii)
2 2
1 1
pero sabemos H(w + 3) = ; H(w − 3) = , reemplazando en (iii)
1 + i(w + 3) 1 + i(w − 3)
se obtiene
1 1 1 1
F(f (t)) = +
2 1 + i(w − 3) 2 1 + i(w + 3)
Luego
F(f 0 (t)) = 0 + 0 + iwF(f (t))
F(f 0 (t)) = iwF(f (t))
5.5.10. Observación
1. S f 0 : R → R es continua y f 00 es absolutamente integrable en R y ademas
lı́m f 0 (t) = 0, entonces
|t|→∞
F(f 00 (t)) = (iw)2 F(f (t))
2. Aplicando la transformada para derivadas superiores se obtiene:
F(f 000 (t)) = (iw)3 F(f (t))
101
5.5.11. Ejemplos
1. Sea f (t) = u(t)e−3t ; hallar F(f (t))
1
Solució[Link] sabe de (5.5.4) F(f (t)) = F (u(t)e−3t ) = .
3 + iw
2. Sea g(t) = u(t − 5)e−3(t−5) ; hallar F(g(t))
2 2 2 g 0 (t) 2
Solución. Sea f (t) = te−t y sea g(t) = e−t ; g 0 (t) = −2te−t , luego − = te−t ,
2
g 0 (t)
esto es − = f (t) aplicando transformada
2
0
g (t) 1
F(f (t)) = F − = − F(g 0 (t)); aplicando (5.5.9)
2 2
1 1 iw
F(f (t)) = − F(g 0 (t)) = − [(iw)F(g(t))] = − F(g(t)) (i)
2 2 2
Ahora de lo estudiado en la sección anterior se sabe:
2 w√
F(g(t)) = F(e−t ) = e− 4 π (ii)
de (ii) en (i)
iw − w2 √
F(f (t)) = − e 4 π
2
dF
Luego F(tf (t)) = (i) (w)
dw
5.5.13. Observación
1. Si n = 2 Z +∞
2
F(t f (t)) = t2 f (t)e−iwt dt (iv)
−∞
de otro lado Z +∞
F (w) = f (t)e−iwt dt (v)
−∞
103
luego en (i),Z cambiando el signo de ”t”
+∞
2πf (−t) = F (w)e−iwt dw; luego hacemos w = t, entonces
−∞
Z +∞
2πf (−w) = F (t)e−iwt dt = F(F (t)); luego
−∞
5.5.15. Ejemplos
2
1. Si f (t) = , halle F(f (t))
1 + t2
2
Solución. Si g(t) = e−|t| , entonces sabemos F(g(t)) = G(w) =
1 + w2
2 2
G(w) = 2
; luego G(t) = = f (t)
1+w 1 + t2
2
deseamos hallar F(G(t)) = F
1 + t2
se propiedad (5.5.14) se obtiene
2
F = F(G(t)) = 2πg(−w) (i)
1 + t2
t2
2. Sea f (t) = , halle F(f (t))
4 + t2
Solución. Sea F(f (t)) = F (w), de la propiedad de la dilatación se sabe
1 w
F(f (at)) = F
|a| a
104
1 4 1 4 1
2
, luego F(f (2t)) = 2
, esto es F(f (2t)) = , luego F f (2t) =
4+w 4+w 4 4 + w2 4
1 1
2
; sea h(t) = f (2t) entonces
4+w 4
1
F(h(t)) = = H(w) (ii)
4 + w2
de la propiedad de simetrı́a se sabe:
Si F(h(t)) = H(w) entonces F(H(t)) = 2πh(−w) (iii)
aplicando (iii) en (ii)
1 1
F 2
= 2π f (−2w) (iv)
4+t 4
de (i) se obtiene f (−2w) = e−|−2w| = e−|2w| reemplazando en (iv)
1 1 −2|w| π −2|w|
F = 2π e = e
4 + t2 4 2
1 π
esto es F 2
= e−2|w| (v)
4+t 2
De otro lado, se sabe de propiedad
2
t 2 1
F 2
=F t 2
= F(t2 H(t))
4+t 4+t
1
donde H(t) = ; esto es
4 + t2
2 2
2 F 2d 2d F 1
F(t H(t)) = (−i) (H(w)) = (−i) (vi)
dw2 dw2 4 + w2
reemplazando (vi) en (v)
t2 d2 π −2|w|
2
F(t H(t)) = F = (−i2 ) e
4 + t2 dw2 2
t2 2
2π d
Luego F = (−i) e−2|w| (vii)
4 − t2 2 dw2
d −2|w| w
sabemos e = −2 e−2|w| , luego, operando se obtiene
dw |w|
d2 −2|w|
e = 4e−2|w| (viii)
dw2
de (viii) en (vii) se obtiene
2
t π
F 2
= (−i)2 (4)e−2|w| = 2(−i)2 πe−2|w|
4+t 2
105
2
3. Si f (t) = t2 e−5t ; halle F(f (t))
2
Solución. De lo estudiado se sabe si g(t) = e−5t , entonces
r
−5t2 π − w2
F(g(t)) = F(e )= e 20 = G(w) (i)
5
d2 G
F(f (t)) = F(t2 g(t)) = (i)2 (w) (ii)
dw2
de (ii) en (i) √
d2
2 π − w2
F(f (t)) = (−i) √ e 20 (iii)
dw2 5
de otro lado
w2
d2 − w2 e− 20 w2
e 20 = − 1 (esto es derivando dos veces) (iv)
dw2 10 10
t
4. Si f (t) = ; hallar F(f (t))
25 + t2
2
Prueba. Sabemos si h(t) = e−|t| entonces H(w) = , ademas si h(5t) = e−5|t| ,
1 + w2
entonces
1 t 1 w 1 2 25 5
F(h(5t)) = F(h( )) == H( ) = w 2 = =
|5| 5 5 5 5 2
5 [25 + w ] 25 + w2
1+
5
esto es si h(5t) = g(t), entonces
5
F(g(t)) = F(h(5t)) = = G(w)
25 + w2
5
aplicando la propiedad de simetrı́a F(G(t)) = 2πg(−w), esto es F = 2πe−5|−w| ;
25 + t2
esto es
5 −5|w| 1 2π −5|w|
F = 2πe esto es F = e
25 + t2 25 + t2 5
106
1 2π −5|w|
esto es F = e (i)
25 − t2 5
de otro lado
1 dF 1 d 2π −5|w|
F t 2
= (i) 2
= (i) e
25 + t dw 25 + t dw 5
luego
t 2π d −5|w|
F 2
= (i) e
25 + t 5 dw
t w
F 2
= −2πi e−5|w|
25 + t |w|
5.6.2. Observación
1. La convolución o producto de dos funciones no siempre existe; ejemplo:
Si f (t) = g(t) = 1 ∀t ∈ R; entonces
Z +∞ Z +∞
(f ∗ g)(t) = f (y)g(t − y) dy = 1 dy = ∞
−∞ −∞
107
5.6.3. Ejemplo
1; |t| ≤ 1
Sea f (t) = ; g(t) = 1 ∀t ∈ R, halle (f ∗ g)(t)
0; |t| > 1
Z +∞
Solución. Se observa que |f (t)| dt < ∞ y |g(t)| ≤ 1; ∀t ∈ R, entonces
−∞
Z +∞ Z −1 Z 1 Z +∞
f (y)g(t − y) dy = 0 · 1 dy + 1 · 1 dy + 0 · 1 dy
−∞ −∞ −1 1
Z +∞
(f ∗ g)(t) = f (y)g(t − y) dy = 2
−∞
P1 ) (f ∗ g)(t) = (g ∗ f )(t) ∀t ∈ R
108
t−y =z t=y+z
hacemos un cambio de variable reemplazando se obtiene
dt = dz
Z +∞ Z +∞
−iw(y+z)
= e f (y)g(z) dz dy; esto es
−∞ −∞
Z +∞ Z +∞ Z +∞ Z +∞
= e−iw(y+z) f (y)g(z) dz dy = e−iwy e−iwz f (y)g(z) dz dy
−∞ −∞ −∞ −∞
Z +∞ Z +∞
= f (y)e−iwy dy g(z)e−iwz dz esto es
−∞ −∞
Z +∞ Z +∞
−iwy
= f (y)e dy g(z)e−iwz dz
−∞ −∞
5.6.6. Observación
1. Del teorema se sabe F((f ∗ g)(t)) = F(f (t))F(g(t)) (i)
aplicando la transformada inversa en (i) se obtiene
109
6. Ecuación de la onda y del calor
Consideremos una E.D.P. lineal de segundo orden no homogénea de tipo hiperbólico
∂ 2 u(x, t) 2
2 ∂ u(x, t)
− c = F (x, t) (i)
∂t2 ∂x2
Donde x e t son las variables independientes y ”c” es una constante real
x : es la variable espacial
t : es la variable temporal (t ≥ 0)
110
Cuando es establece una ecuación diferencial correspondiente a un problema fı́sico da-
do, por lo general es necesario hacer supuestos de simplificación, para asegurar de que la
ecuación no sea muy complicado.
2. La tensión causada por el estiramiento de la cuerda es tan grande, que puede des-
preciarse la fuerza gravitacional sobre la cuerda; ası́ el peso de la cuerda puede ser
despreciada.
Bajo todos estos supuestos, se puede esperar que la solución u(x, t) de la ecuación dife-
rencial que se va a obtener, describirá de una forma razonable las pequeñas vibraciones
de la cuerda fı́sica.
111
6.1.3. Modelamiento
De la figura se tiene el segmento AB arbitrario de la cuerda, de longitud ∆S. Sean x,
x + ∆x dos puntos de la cuerda en reposo (eje X) tal que las coordenadas de los puntos
A y B son A(x, u(x, t)); B(x + ∆x, u(x + ∆x, t)) respectivamente.
Si T1 y T2 son las fuerzas de tensión actuando en los puntos A y B respectivamente; y
estas fuerzas son tangenciales al perfil de la cuerda en los puntos A yB. Puesto que el
movimiento es vertical, la suma de los componentes horizontales de la fuerza de tension
debe ser igual a cero, esto es:
Donde α1 y α2 son los ángulos agudos; formados por las tangentes en los puntos A y B
con el eje X.
1 1
Se sabe sec2 α1 = 1+tan2 α1 ; luego cos2 α1 = 2
, esto es; cos α1 = ;
1 + tan α1 [1 + u2x (x, t)]1/2
1
en forma similar cos α2 = , donde tan α1 = ux (x, t),
[1 + u2x (x + ∆x, t)]1/2
tan α2 = ux (x + ∆x, t)
De los supuestos realizados, se sabe u2x → 0; luego se obtiene
1 ∼ 1 ∼
cos α1 = 1/2
=1 ∧ cos α2 = =1 (ii)
[1 + u2x (x, t)] [1 + u2x (x + ∆x, t)]1/2
Esto significa que la tensión de la cuerda es constante. De otro lado la fuerza vertical
actuante (componente vertical) sobre el segmento AB es: |T | sen α2 , con dirección hacia
arriba y |T1 | sen α1 con dirección hacia abajo. Luego la fuerza resultante es
s 2
u(x + ∆x, t) − u(x, t)
esto es ∆S = ∆x 1 +
∆x
112
luego si ∆x → 0 entonces se tiene
s 2
u(x + ∆x, t) − u(x, t)
∆S = ∆x 1 + lı́m
∆x→0 ∆x
p
∆S = ∆x 1 + u2x (x, t) (v)
∆S ≈ ∆x (vi)
De otro lado, si ρ es la densidad de la cuerda AB, entonces la masa ”m” del segmento
AB es:
F ≈ ρ · ∆x · a (viii)
∂ 2 u(ε, t)
puesto que a = , donde x < ε < x + ∆x
∂t2
Luego en (viii) se obtiene
∂ 2 u(ε, t)
F ≈ ρ · ∆x · (ix)
∂t2
luego igualando (ix) y (ix) se obtiene
∂ 2 u(ε, t)
R = ρ · ∆x · = |T | [ux (x + ∆x, t) − ux (x, t)] , (x)
∂t2
donde x < ε < x + ∆x
dividiendo (x) por ∆x y tomando limite, cuando ∆x → 0 se obtiene
∂ 2 u(x, t)
u(x + ∆x, t) − u(x, t)
ρ = |T | lı́m
∂t2 ∆x→0 ∆x
2 2
∂ u(x, t) ∂ u(x, t)
ρ 2
= |T | ; luego se tiene
∂t ∂t2
∂ 2 u(x, t) 2
2 ∂ u(x, t)
− c =0 (xi)
∂t2 ∂t2
1/2
|T |
donde c =
ρ
113
Luego, la ecuación (x), es llamado la ecuación unidimensional de la onda, el cual gobierna
el problema planteado. Esto es una E.D.P. de tipo hiperbólico homogéneo con coeficientes
constantes. La constante c tiene la dimensión de una velocidad; esto se verifica sustitu-
yendo apropiadamente las unidades de |T | y ρ. Esta constante c representa la velocidad
máxima de proporción de una perturbación en la cuerda.
Ahora del problema planteado, si existe una fuerza vertical externa, supongamos f (x, t)
por unidades de longitud, actuando sobre la cuerda, entonces el termino f (x, t) · ∆x, se
debe añadir en el lado derecho de la ecuación (x), esto es
∂ 2 u(ε, t)
∂u ∂u
ρ · ∆x · = |T | (x + ∆x, t) − (x, t) + f (x, t) · ∆x (xii)
∂t2 ∂x ∂x
satisfaciendo
u(x, t) = f (x) (ii)
, x ∈ R, t ≥ 0
ut (x, t) = g(t) (iii)
114
De lo estudiado, la ecuación (i) es de tipo hiperbólico y a la solución general de la ecuación
(i) es:
u(x, t) = ϕ(x + ct) + ψ(x − ct) (iv)
utilizando la condición inicial (ii) en (iv) se obtiene
ϕ(x) + ψ(x) = f (x) (v)
Además, derivando (iv) se tiene
esto es
Z x−ct x+ct x+ct
f (x + ct) + f (x − ct)
Z Z
1
u(x, t) = + g(r) dr + g(r) dr − g(r) dr
2 2c 0 x−ct 0
luego
x+ct
f (x + ct) + f (x − ct)
Z
1
u(x, t) = + g(r) dr (ix)
2 2c x−ct
luego la expresión (ix) representa la única solución del problema propuesto (i) satisfaciendo
(i) y (iii). Esto solución (ix) se denomina la formula de D’alembert del problema planteado.
115
6.2.1. Ejemplos
1. Resolver el problema:
Solución
2. Resolver el problema
Solución
∂ 2 u(x, t) 2
2 ∂ u(x, t)
− c =0 (1)
∂t2 ∂x2
x ∈ [0, L], t > 0
116
satisfaciendo la ecuación de frontera u(0, t) = 0, u(L, t) = 0 (2)
y la condición inicial
u(x, 0) = f (x) (3)
ut (x, 0) = g(x) (4)
Las vibraciones de una cuerda elástica, están gobernadas por la ecuación unidimensional
de la onda (1). La cuerda elástica es de longitud ”L” y esta fijado en sus extremos x = 0,
x = L. La condición (2) significa que en los extremos de la cuerda el desplazamiento o la
deflexión de la cuerda es cero. Anotar que u(x, t) es la deflexión de la cuerda en el punto
x de la cuerda y en el instante t.
Además (3) significa la deflexión inicial de la cuerda en el instante t = 0; esto es la forma
del movimiento de la cuerda se comporta como la función f (x). También (4) denota la
velocidad inicial de un punto de la cuerda en el instante t = 0; y esta velocidad se com-
porta como la función g(x).
El problema consiste en encontrar una solución de la ecuación (1) satisfaciendo las con-
diciones de frontera y las condiciones (2), (3) y (4).
Para resolver este problema, emplearemos el método de Fourier (variables separables).
Consiste en suponer que la solución u(x, t) se puede expresar como un producto de dos
funciones independientes, uno que solo dependa de ”t” y el otro solo dependa de x, luego
de reemplazar en la ecuación (1) y derivarlo adecuadamente se obtendrá dos E.D.O’s; de
los cuales se obtendrá sus respectivas soluciones, donde uno de ellos satisface la condición
de frontera. Finalmente, usando series de Fourier, se realiza la composición de estas solu-
ciones a fin de llegar a una solución de la E.D.P.(1) satisfaciendo la condición de frontera
y la condición inicial. Supongamos que la solución de la ecuación (1) se puede expresar
por
u(x, t) = F (x)G(t) (5)
Donde F y G son funciones derivables hasta el orden dos, derivando (5) dos veces respecto
a la variable x y respecto a t se obtiene:
uxx (x, t) = F 00 (x)G(t) y utt (x, t) = F (x)G00 (t)
luego reemplazando en la ecuación (1) se obtiene
F 00 (x) G00 (t)
c2 F 00 (x)G(t) = F (x)G00 (t), esto es = 2
F (x) c G(t)
puesto que F solo depende de x y G solo depende de t, entonces ambas razones son total-
mente independientes, luego ambas razones son constantes, luego se tiene la proporción
G00 (t) F 00 (x)
= =k (6)
c2 G(t) F (x)
donde k es la constante de proporcionalidad, y de (6) se obtiene dos E.D.O.’s
F 00 (x) − kF (x) = 0 (7)
G00 (t) − kc2 G(t) = 0 (8)
117
De otro lado, la condición de frontera (2), debe relacionarse con algunas de las E.D.O.s
obtenidas.
Sabemos de (5) que la solución de la ecuación (1) es de la forma:
u(x, t) = F (x)G(t), empleando la condición (2) se obtiene u(0, t) = F (0)G(t) = 0 y
u(L, t) = F (L)G(t) = 0 de estas dos ecuaciones se obtiene:
caso contrario se obtendrı́a, G(t) ≡ 0; y por lo tanto la solución (5) del problema seria la
solución trivial (u ≡ 0), que es lo que no deseamos.
Luego de (7) y (9) se tiene una E.D.O. con condición frontera
F 00 (x) − kF (x) = 0
(10)
F (0) = 0, F (L) = 0
F 00 (x) + γ 2 F (x) = 0
(11)
F (0) = 0, F (L) = 0
118
luego de (12) y (13) se obtienen:
nπx
un (x, t) = Fn (x)Gn (t) = sen [Bn cos λn t + cn sen λn t] ∀n = 1, 2, . . . (14)
L
los cuales son soluciones de la E.D.P.(1), además satisface la condición de frontera (2).
Luego se tiene una familia de soluciones de (1) y (2). Del principio de superposición, se
sabe que toda combinación lineal finita de soluciones, también es solución de la E.D.P. Sin
embargo nuestra familia de soluciones es infinito y deseamos que (14) satisfaga la condición
inicial (3) y (4), para ello extenderemos el principio de superposición para una familia
infinita de soluciones, es decir la combinación lineal de infinitas soluciones lo denotamos
por
∞ ∞ h
X X nπx nπx i
u(x, t) = Fn (x)Gn (t) = Bn sen cos λn t + cn sen sen λn t (15)
n=1 n=1
L L
esta claro que esta serie infinita satisface (1) y (2). Ahora veamos que satisfaga (3) y (4).
Evaluando (15) en t = 0 se tiene
∞
X nπx
u(x, 0) = sen [Bn cos λn (0) + cn sen λn (0)] = f (x)
n=1
L
esto es
∞
X nπx
Bn sen = f (x) (16)
n=1
L
Ahora, debemos obtener los coeficientes Bn de la serie para que dicho serie converja a
f (x).
De lo estudiado, se observa que la serie que aparece en (16) es la serie de Fourier de seno
de la función f (x), luego se sabe que sus coeficientes son de la forma
Z L
2 nπx
Bn = f (x) sen dx, n = 1, 2, 3, . . . (17)
L 0 L
evaluando en t = 0
∞
∂u(x, t) X nπx
= cn λn sen = g(x) (18)
∂t n=1
L
119
luego se debe obtener los coeficientes (cn λn ) de la serie para que dicha serie converja a
g(x). Pero, la serie que aparece en (18), es la serie de Fourier de seno de la función g(x),
luego sabemos que
2 L
Z nπx
cn λn = g(x) sen dx n = 1, 2, . . .
L 0 L
Z L
2 nπx
cn = g(x) sen dx
λn L 0 L
Z L
2 nπx
cn = g(x) sen dx (19)
cnπ 0 L
Luego reemplazando (17) y (19) en (15) se obtiene
∞
2 L nπx
X Z
u(x, t) = f (x) sen dx cos λn t+
n=1
L 0 L
Z L nπx
2 nπx
+ g(x) sen dx sen λn t sen (20)
cnπ 0 L L
Luego la expresión (20) representa la solución del problema planteado, la solución esta
expresad mediante la serie de Fourier; además, para que esta serie converja, f (·) debe ser
derivable hasta el orden dos y g(·) debe ser derivable hasta el orden uno.
6.3.1. Observación
1. Del problema planteado (1)-(4) se presentan casos particulares. Dado el problema
2. Dado el problema:
utt (x, t) − c2 utt (x, t) = 0 0 < x < L; t > 0
u(0, t) = u(L, t) = 0 t>0
u(x, 0) = 0
x ∈ [0, L]
ut (x, 0) = g(x)
120
3. De (14) se obtuvo
nπx
un (x, t) = sen [Bn cos λn t + cn sen λn t]
L
Las funciones un son llamadas funciones propias (o funciones caracterı́sticas) y los
cnπ
valores λn = son llamados valores propios o valores caracterı́sticos de la cuerda
L
vibratoria. Además el conjuntos{λ1 ; λ2 ; λ3 ; . . .} es llamado espectro.
Cada solución un (x, t), representa un movimiento armónico que tiene la frecuencia
λn cn
= ciclos por unidad de tiempo.
2π 2L
Este movimiento armónico, es llamado el n-ésimo modo normal de la onda.
6.3.2. Ejemplos
1. Resolver el siguiente problema de V.I.F.
satisfaciendo
u(0, t) = u(x, π) = 0
5x
u(x, 0) = x cos 0≤x≤π
2
ut (x, 0) = 0
Solución
5x
En este caso f (x) = x cos ; g(x) = 0, luego la forma de la solución es
2
∞ Z π nπx
X 2 nπx
u(x, t) = f (x) sen dx cos λn t sen
n=1
π 0 π π
cnπ 2nπ
sabemos λn = = = 2n
L π
2 π
Z
Bn = f (x) sen(nx) dx ∀n = 1, 2, . . .
π 0
2 π
Z
5x
operando se obtiene Bn = x cos sen(nx) dx
π 0 2
8n(−1)n+1
Bn = , luego la solución es
(5 + 2n)2 (5 − 2n)2
∞
8 X n(−1)n+1
u(x, t) = sen(nx) cos(2nt)
π n=1 (25 + 4n)2
121
2. Resolver el problema
utt (x, t) − 16uxx (x, t) = 0, 0 < x < π, t > 0
u(x, 0) = 0; ut (x, 0) = sen 2x, x ∈ [0, 2π]
u(0, t) = u(π, t) = 0, t≤0
Solución
cnπ
En este caso, f (x) = 0, g(x) = sen 2x, c = 4, L = π, λn = = cn, la solución
π
es de la forma
∞ Z π
X 2 nπx nπx
u(x, t) = g(x) sen dx sen λn t sen
n=1
4nπ 0 π π
π
1 sen(2 − n)π sen(2 + n)π
Z
1
cn = sen 2x sen(nx) dx = −
2nπ 0 4nπ 2−n 2+n
reemplazando se obtiene
∞
X 1 sen(2 − n)π sen(2 + n)π
u(x, t) = − sen(4nt) sen(nx)
n=1
4π 2−n 2+n
3. Resolver el problema
utt (x, t) − 9uxx (x, y) = 0
Satisfaciendo
u(0, t) = u(4, t) = 0
u(x, 0) = 2 sen(πx)
ut (x, 0) = 0
Solución
L = 4, c = 3, f (x) = 2 sen(πx), g(x) = 0
Z 4
2 nπx 3nπ
Luego Bn = f (x) sen dx, λn =
4 0 4 4
la solución es ∞
X nπx 3nπ
u(x, t) = Bn sen cos t
n=1
4 4
desarrollando la serie, se obtiene
π π 3π
B1 sen x + B2 sen x + B3 sen x + B4 sen πx + · · · = 2 sen(πx)
4 2 4
igualando coeficientes se obtiene: B4 = 2, Bn = 0, n 6= 4
Luego Bn = 2; la solución es
122
6.4. Conducción o difusión de calor(Ecuación del calor)
6.4.1. Introducción
Supongamos que se tiene una barra delgada de metal de longitud L que se coloca en
el eje 0X, esto es
Figura 6.4.2
123
consideremos el elemento de volumen de la barra incluido entre los dos planos vecinos
paralelos a A y que se ha anotado por B y C, a distancias x y x + ∆x respectivamente,
de A. Denotemos la temperatura en el plano B en el instante t por u(x, t). Entonces en
el plano C en el instante t su temperatura sera u(x + ∆x, t).
Para poder continuar en la obtención de la formulación matemática necesitamos tener en
cuenta las dos leyes fı́sicas correspondientes a la transferencia del calor.
Segunda ley La cantidad de calor que fluye a través de un área (B o C) por unidad
de tiempo, es proporcional a la tasa de cambio de la temperatura con respecto a la
distancia perpendicular del área.
Si tomamos como positiva la dirección de izquierda a derecha en la Figura 6.4.2, se
puede escribir
∂u
Q = −KA∆t (i)
∂x
siendo
El signo menos en (i) muestra que Q es positivo (esto es, el flujo es a la derecha) cuando
∂u
es negativo(esto es, cuando la temperatura está decreciendo a medida que se va hacia
∂x
∂u
la derecha). De forma similar Q es negativo cuando es positivo. Esto va en sintonia
∂x
con los hechos fı́sicos. Usando (i) se puede decir que la cantidad de calor que fluye de
izquierda a derecha a través del plano B en la Figura 6.4.2 es
∂u
−KA∆t
∂x x
En forma similar, la cantidad de calor que fluye de izquierda a derecha a través del plano
C, en la Figura 6.4.2 es
∂u
−KA∆t
∂x x+∆x
De aquı́ se obtiene que la cantidad neta que se acumula en el volumen entre C y B, es la
cantidad que entra por B menos la cantidad que sale por C, esto es
∂u ∂u ∂u ∂u
−KA∆t − −KA∆t = KA∆t − (ii)
∂x x ∂x x+∆x
∂x x+∆x ∂x x
124
Esta cantidad de calor acumulado eleva o baja la temperatura del elemento de volumen
si (ii) es positivo o negativo. Luego de la primera ley se tiene
∂u ∂u
KA∆t − = m · S · ∆u = ρ · A · ∆x · S · ∆u (iii)
∂x x+∆x ∂x x
pues, la masa del elemento de volumen es ρ · A · ∆x.
Se debe mencionar que (ii) es solo aproximadamente cierto, siendo el grado de aproxima-
ción mayor, cuando más pequeño sean los valores ∆x, ∆u y ∆t.
Dividiendo en (iii) por A · ∆x · ∆t se obtiene
K ∂u ∂u ρ · S · ∆u
(x + ∆x, t) − ((x, t)) = (iv)
∆x ∂x ∂x ∆t
aplicando limite en (iv) cuando ∆x → 0 se obtiene
∂u ∂u
(x + ∆x, t) − (x, t)
∂x ∂x ρ · S · ∆u ρ · S · ∆u
K lı́m = lı́m =
∆x→0 ∆x ∆x→0 ∆t ∆t
esto es
∂ 2u ρ · S · ∆u
K2
(x, t) = (v)
∂x ∆t
Ahora aplicando limite en (v) cuando ∆t → 0
∂ 2u
∆u u(x, t + ∆t) − u(x, t)
K 2 (x, t) = ρ · S lı́m = ρ · S lı́m
∂ x ∆t→0 ∆t ∆t→0 ∆t
∂ 2 u(x, t) ∂u(x, t)
luego K = ρ · S esto es
∂x2 ∂t
∂u(x, t) K ∂ 2 u(x, t) K
− 2
= 0; tomando c =
∂t ρ · S ∂x ρ·S
se obtiene
∂u(x, t) ∂ 2 u(x, t)
−c =0 (vi)
∂t ∂x2
la ecuación (vi), es llamado la ecuación del flujo de calor o de conducción de calor en una
dimensión. También se le llama la ecuación del calor unidimensional. Además c representa
el coeficiente de difusividad.
6.4.3. Observación
1. Si la superficie de la barra no estuviera aislada se tendrı́a que considerar un termino
adicional en la ecuación (iii); el cual serı́a la cantidad de calor que escapa (o fluye
dentro) del elemento, y cuya ecuación en este caso será:
∂u(x, t) ∂ 2 u(x, t)
=c − r(u − u0 ) (vii)
∂t ∂x2
siendo ”r” una constante y u0 la temperatura de los alrededores.
125
6.5. Problema de Cauchy Para la conductividad térmica (Barra
infinita)
consideremos la ecuación de conductividad térmica
∂u(x, t) ∂ 2 u(x, t)
− a2 =0 (viii)
∂t ∂x2
k
donde a2 = es positivo.
ρS
El problema de Cauchy se plantea de la siguiente forma:
Hallar una función u(x, t) que satisfaga la ecuación
ut (x, t) − a2 uxx (x, t) = 0, t > 0, x ∈ h−∞, +∞i
y satisface la condición (ix)
u(x, 0) = ϕ(x), x ∈ h−∞, +∞i
El sentido fı́sico del problema consiste en determinar la temperatura de una varilla ho-
mogénea ilimitado en cualquier tiempo t > 0, según su temperatura conocida ϕ(x) en el
instante t = 0. Considerando que la superficie lateral de la varilla está termoaislada de
tal manera que a través de ello el calor no se evacua.
La expresión Z +∞
1 (x−λ)2
u(x, t) = √ ϕ(λ)e− 4a2 t dλ, t > 0 (x)
2a πt −∞
nos proporciona, una solución del problema de valor inicial (ix) y se denomina la integral
de Poisson.
Luego si la función ϕ(x) es continua y acotada, entonces la función u(x, t) de (x) admite
derivada de cualquier orden respecto a x y t, cuando t > 0, satisface la ecuación (viii) y
la condición inicial dada
x−λ √
de (x) consideramos √ = v, λ = x − 2a tv
√ 2a t
Luego dλ = −2a t dv, luego reemplazando
1
Z +∞ √ 2 √
u(x, t) = √ ϕ(x,2a tv)e−v (−2a t dv)
2a πt −∞
1
Z +∞ √ 2
=√ ϕ(x − 2a tv)e−v dv
π −∞
luego hacemos t → 0+
1 +∞
Z √ −v 2 1
Z +∞
2
u(x, t) = lı́m+ ϕ(x − 2a tv)e dv = √ ϕ(x)e−v dv
π −∞ t→0 π −∞
Z +∞ √
ϕ(x) 2 ϕ(x) π
u(x, 0) = √ e−v dv = √ = ϕ(x)
π −∞ π
En la clase de funciones acotadas
Cu(x,t) = {u(x, t)/ |u(x, t)| < M, x ∈ h−∞, +∞i, t > 0}
La solución del problema de Cauchy (ix) es única y depende continuamente de los dos
datos iniciales.
126
6.6. Propagación del calor en una varilla finita
6.6.1. Planteamiento del problema
Si la varilla tiene longitud finita L y ocupa el segmento 0 ≤ x ≤ L del eje 0X, entonces
para plantear el problema de propagación del calor en tal varilla, además de la ecuación
2
∂u(x, t) 2 ∂ u(x, t)
−a =0 (i)
∂t ∂x2
con condición inicial
u(x, t)|t=0 = ϕ(x) (ii)
hace falta fijar también el régimen de temperatura en los extremos de la varilla, es decir
en x = 0 x = L, esto es, plantear las condiciones de frontera. Las condiciones de frontera
pueden ser diferentes en función del régimen de temperaturas en los extremos de la varilla.
Se examinará tres tipos principales de las condiciones de frontera
Tipo (I) La varilla mantiene sus extremos a una temperatura constante igual a cero. En
este caso
u(0, t) = 0; u(0, t) = 0 ∀t ≥ 0
Tipo (II) La varilla tiene extremos aislados. En este caso consideraremos que la varilla
mantiene sus extremos aislados, de modo que no hay flujo de calor que entre ni
salga. Luego
ux (0, t) = 0; ux (x, t) = 0 ∀t ≥ 0
Tipo (III) En los extremos de la varilla vienen dadas las relaciones lineales entre la
función u y su derivada
∂u ∂u
= λ [u(0, t) − Θ(t)] ; = −λ [u(L, t) − Θ(t)]
∂x x=0 ∂x x=L
∂u(x, t) ∂ 2 u(x, t)
− c2 =0 (i)
∂t ∂x2
x ∈ h0, Li, t > 0
127
satisfaciendo la condición de frontera
u(0, t) = 0, u(L, t) = 0 t ≥ 0 (ii)
y la condición inicial
u(x, 0) = f (x) x ∈ [0, L] (iii)
Donde u(x, t) representa la temperatura de un cuerpo (alambre) de material homogéneo.
K
Aquı́ c2 representa la difusión térmica (c2 = ), G es el calor especifico y ρ es la densidad
Gρ
del material del cuerpo. La ecuación (i), es llamado la ecuación conductiva térmica o la
ecuación del calor.
Como una importante aplicación de esta ecuación, consideremos la temperatura de una
barra o alambre delgado de una sección transversal constante y material homogéneo,
el cual se orienta a lo largo del eje X. Dicho alambre es de longitud L y está aislado
perfectamente en su superficie lateral, de tal forma que el calor solo fluye en la dirección
del eje X. Luego u depende únicamente de la variable x y del tiempo t, por lo tanto, la
ecuación del calor queda como la ecuación unidimensional del calor.
La condición (ii), condición de frontera significa que en los extremos del alambre en
x = 0 ∧ x = L la temperatura del objeto es cero.
La condición (iii), condición inicial del alambre, esta representado por f (x) en el instante
t = 0, ∀x ≥ 0.
Nuestro objetivo es determinar una solución de la ecuación (i), satisfaga la condición (ii) y
(iii) respectivamente para ello emplearemos el método de Fourier(Separación de variables).
Sea
u(x, t) = F (x)G(t) (iv)
derivando (iv) respecto a x dos veces, respecto a t una vez y luego reemplazando en (i)
se obtiene
G0 (t) F 00 (x)
= (v)
c2 G(t) F (x)
se observa que ambas razones de (v) son independientes, pues la primera razón depende
de t y la segunda razón sólo depende de x; luego se tiene
G0 (t) F 00 (x)
= =k (vi)
c2 G(t) F (x)
donde k es la constante de proporcionalidad de la proporción.
Para que la solución de la ecuación (i) no sea la trivial, consideraremos a esta constante
por k = −ρ2 .
luego reemplazando esto en (vi) y despejando se obtiene dos E.D.O.s
F 00 (x) + ρ2 F (x) = 0 (vii)
G0 (t) + c2 ρ2 G(t) = 0 (viii)
De otro lado de la condición de frontera (ii) y de (iv) se obtiene.
u(0, t) = F (0)G(t) = 0
u(L, t) = F (L)G(t) = 0
128
puesto que se desea que la solución (iv) no sea la trivial, entonces se concluye que
luego de (vii) y (ix) se tiene una E.D.O. de segundo orden con condición de frontera
F 00 (x) + ρ2 F (x) = 0
(x)
F (0) = F (L) = 0
Es claro que la expresión (xiv) satisface la ecuación (i) y la condición (ii). Veamos ahora
que (xiv) satisfaga la condición (iii).
Empleando (xiv) y la condición inicial (iii) se obtiene
X∞ nπx
u(x, 0) = Bn sen e0 = f (x), esto es
n=1
L
∞
X nπx
Bn sen = f (x) n = 1, 2, 3, . . . (xv)
n=1
L
Luego para que la serie obtenida converja hacia la función f (x), los coeficientes Bn deben
ser los coeficientes de la serie de Fourier de seno de la función f (x).
De lo estudiado se sabe que los coeficientes Bn se obtiene de
2 L
Z nπx
Bn = f (x) sen dx n = 1, 2, 3, . . . (xvi)
L 0 L
La expresión (xvii) representa la solución del problema planteado satisfaciendo (i), (ii) y
(ii).
Además, la serie infinita converge a la solución del problema planteado si f (·) es una
función continua por tramos en [0, L] y tiene derivadas laterales en cada punto interior
del intervalo [0, L].
6.6.3. Ejemplos
(1). Halle la solución del problema de Cauchy
∂u(x, y) ∂ 2 u(x, t)
− = 0, t > 0, x ∈ h−∞, +∞i
∂t ∂x2
x2
u(x, t)|t=0 = e− 2 , x ∈ h−∞, +∞i
130
Solución
La ecuación mostrada, es la ecuación del calor o de la conductividad térmica, donde
la varilla es una varilla infinita. Para resolver este problema, emplearemos la formula
de Poisson.
esto es Z +∞
1 (x−λ)2
u(x, t) = √ f (λ)e 4a2 t dλ, t > 0
2a πt −∞
λ2
donde a = 1, f (λ) = e− 2 ; luego
Z +∞
1 λ2 (x−λ)2
u(x, t) = √ e− 2 e− 4t dλ (α1 )
2 πt −∞
2
Z +∞ 2
x 1 (1+2t) x
=e − 2(1+2t)
e− 2 2t [λ− 1+2t ] dλ (α2 )
−∞
r
x 1 + 2t 1 + 2t
hacemos λ − = α, r2 = , r = , dλ = dα; reemplazando en
1 + 2t 2t 2t
(α2 ) se obtiene
−x2
Z +∞ Z +∞
−x2 (rα)2
− 2 e 2(1+2t) (rα)2
=e 2(1+2t) e dλ = e− 2 d(rα) (α3 )
−∞ r −∞
| {z }
(?)
Z +∞
u2
√
Z +∞
u2 √
De otro lado, se sabe; e du = πy e− 2 du = 2π, reemplazando en
−∞ −∞
(α3 ) se obtiene:
−x2 √
1 e 2(1+2t) √ 1 2t √ −x2
u(x, t) = √ r 2π = √ √ 2πe 2(1−2t)
2 πt 1 + 2t 2 πt 1 + 2t
2t
−x 2 −x2
e 2(1+2t) e 2(1−2t)
u(x, t) = √ ⇒ u(x, t) = √ , t>0
1 + 2t 1 + 2t
esto representa la único solución del problema.
131
(2). Problemas
Supongamos que se tiene una varilla homogénea infinita. Demostrar que si la tem-
2 2
peratura inicial es: u(x, 0) = f (x); donde f (x) = u0 e−% x , x ∈ h−∞, +∞i, (u0 >
0, % > 0 constantes) en un momento cualquiera t > 0, entonces la temperatura es
% 2 x2
−
u0 e 1+4a2 %2 t ut − a2 uxx = 0
u(x, t) = p
1 + 4a2 %2 t u(x, 0) = f (x)
(3). Una barra metálica de 100 cm de longitud tiene sus extremos x = 0, x = 100
mantenidos a 0 ºC inicialmente, la mitad de la barra esta a 60 ºC, mientras que la
otra mitad esta a 40 ºC. Asumiendo un coeficiente de difusión de 0,16 unidades e.g.s
y un entorno aislado, encontrar la temperatura en cualquier posición de la barra en
cualquier tiempo.
Solución
La ecuación de conducción de calor es
∂u(x, t) ∂ 2 u(x, t)
− 0,16 =0 (α1 )
∂t ∂x2
satisfaciendo
F 00 (x) + λ2 F (x) = 0
; resolviendo este problema se obtiene
F (0) = 0, F (100) = 0
nπx
Fn (x) = sen , n = 1, 2, . . .
100
Resolviendo la otra ecuación se obtiene
2
Gn (t) = Be−λn t , n = 1, 2, 3, . . .
132
0, 16n2 π 2
donde λ2n =
(100)2
nπx 2
Luego un (x, t) sen Bn e−λn t , es una familia de soluciones que satisface (i) y
100
(ii), consideremos la combinación lineal infinita de soluciones por
∞ ∞ nπx 2
X X
u(x, t) = un (x, t) = Bn sen e−λn t (α5 )
n=1 n=1
100
luego
Z 100 Z 50 Z 100
2 nπx 1 nπx nπx
Bn = f (x) sen dx = 60 sen dx + 40 sen dx
100 0 100 50 0 100 50 100
6 50 4 100
Z nπx Z nπx
Bn = sen dx + sen dx
5 0 100 5 50 100
6
100 nπx 50 4 100 nπx 100
= − cos + − cos
5 nπ 100 0 5 nπ 100 50
40 nπ 80
Bn = 3 − cos − (−1)n
nπ 2 nπ
reemplazando en (α4 ) se obtiene
∞
X 40 nπ 80 nπx 2
u(x, t) = 3 − cos − n
(−1) sen e−λn t (α6 )
n=1
nπ 2 nπ 100
∂u(x, t) ∂ 2 u(x, t)
−4 = 0, 0 < x < 2, t > 0
∂t ∂x2
u(0, t) = 0; u(2, t) = 0 t≥0
2
u(x, 0) = x (2 − x); 0≤x≤2
Solución
La ecuación representa la ecuación del calor, donde la longitud de la varilla es
L = 2, c = 2 y f (x) = x2 (2 − x)
Aplicando el método de Fourier se obtiene
X nπx 2
u(x, t) = Bn sen e−λn t
2
133
4n2 π 2
donde λ2n = = n2 π 2
2Z2
2 2 2 nπx
Luego Bn = x (2 − x) sen dx, operando se obtiene
2 0 2
−32
Bn = [1 + 2(−1)n ]
n3 π 3
luego la solución del problemas es
X 32 nπx 2 2
u(x, t) = − 3 3
(1 + 2(−1)n
) sen e−n π t
nπ 2
∂u(x, t) ∂ 2 u(x, t)
− c2 = 0 0 < x < π, t > 0
∂t ∂x2
u(0, t) = u(π, t) = 0, t≥0
u(x, 0) = sen x; 0≤x≤π
desarrollando la serie
134
6.8. Ecuación del calor para una barra de extremos aislados
Consideremos la conducción la conducción del calor en una barra con extremos aisla-
dos; aquı́ no hay perdida de energı́a por los extremos. Si la temperatura inicial es f (x),
la función temperatura esta modelado por el problema con valor de frontera
∂u(x, t) ∂ 2 u(x, t)
− c2 = 0; 0 < x < L, t > 0 (i)
∂t ∂x2
ux (x, t) = 0, ux (L, t) = 0; t>0 (ii)
u(x, 0) = f (x), 0≤x≤L (iii)
reemplazando en (i)
G0 (t) F 00 (x)
F (x)G0 (t) = c2 F 00 (x)G(t) luego, = =k
c2 G(t) F (x)
luego, considerando la constante de proporcionalidad por
k = −ρ2
de otro lado, de la condición de frontera (ii) se tiene ux (x, t) = F 0 (x)G(t), luego evaluando
en los extremos de la barra
ux (0, t) = F 0 (0)G(t) = 0
(vii)
ux (L, t) = F 0 (L)G(t) = 0
135
utilizando la condición de frontera
6.8.1. Observación
1. Luego los un (x, t) son las funciones propias correspondientes a los valores propios
cnπ
λn = .
L
2. En caso que el valor propio sea λ0 = 0, la función propia u0 =cte, el cual es la
solución del problema si la temperatura inicial es constante.
esto es ∞
X nπx 2
u(x, t) = An cos e−λn t (xiii)
n=0
L
Ahora, de la condición inicial se tiene
∞
X nπx
u(x, 0) = An cos e0 = f (x) (xiv)
n=0
L
136
la expresión (xiv) es la serie de Fourier de cosenos de f . Luego de lo estudiado se tiene
Z L Z L
2 nπx 1
An = f (x) cos dx; A0 = f (x) dx
L 0 L L 0
n = 1, 2, 3 . . .
Luego la solución del problema (i)-(iii) esta dado por
∞ Z L
X 2 nπx nπx
u(x, t) = f (x) cos dx cos e−λn t (xv)
L n=0 0 L L
cnπ
donde λn =
L
6.8.2. Ejemplo
Obtener la temperatura de una barra delgada aislada lateralmente y también por los
extremos, cuya longitud es L, sabiendo que su temperatura inicial está dado por
x; 0 < x < L/2
f (x) =
L − x; L/2 < x < L
Solución
En este caso se tiene la ecuación del calor
ut − c2 uxx = 0 t > 0, 0 < x < L
ux (0, t) = ux (L, t) = 0 t≥0
u(x, 0) = f (x) 0≤x≤L
la solución es
∞
L X 2L h nπ i nπx 2
u(x, t) = + 2 2
2 cos − (1 − (−1)n
) cos e−λn t
8 n=1 n π 2 L
137
6.8.3. Observación
L
1. Se observa que cuando t crece la solución u(x, t) decrece y que se aproxima a esto
8
L
u(x, t) →
es 8 .
t→∞
2. Esto es plausible por que no es posible que escape calor de esta barra totalmente
aislada.
3. En contraste, con el enfriamiento de los extremos, esto lleva a una perdida de calor
y u(x, t) → 0 por la temperatura que se conserva en los extremos de la barra.
138
b) busquemos la solución del problema (5)-(7) en la forma de la serie
∞
X nπx
v(x, t) = Tn (t) sen (11)
n=1
L
(nπ)
según las funciones propias sen x del problema de contorno
L
X 00 (x) + λX(x) = 0
X(0) = X(L) = 0
donde Z L
2 nπ
fn (t) = f (ε, t) sen ε dε, n = 1, 2, . . . (14)
L 0 L
comparando los desarrollos de (12) y (13) de la forma f (x, t) en la serie de Fourier,
obtendremos
∞ ∞
n 2 π 2 a2
X nπ X nπ
0
Tn (t) + Tn (t) sen x= fn (t) sen x
n=1
L2 L n=1
L
n 2 π 2 a2
nπ nπ
esto es Tn0 (t) + Tn (t) sen x = f n (t) sen x
L2 L L
luego
n2 π 2 a2
+ Tn0 (t) Tn (t) = fn (t) ∀x 6= 0, L (15)
L2
aprovechando la condición para v(x, t) en (11) se tiene
∞
X nπ
v(x, 0) = Tn (0) sen x=0 0≤x≤L
n=1
L
139
de aquı́ resulta que
Tn (0) = 0 n = 1, 2, . . . (16)
luego de (15) y (16) se tiene el problema
2 2 2
n π a
Tn0 (t) +
2
Tn (t) = fn (t) (17)
L
satisfaciendo Tn (0) = 0
o
uxx (x, y, z) + uyy (x, y, z) + uzz (x, y, z) = 0 (2)
con ciertas condiciones apropiadas de frontera.
Las ecuaciones de la forma (1) o (2) son llamados las ecuaciones de Laplace. La ecuación
(1) es una E.D.P. de tipo elı́ptico con dos variables independientes (la E.D.P.(2) en forma
similar en 3 variables independientes).
Posteriormente se estudiará los problemas de Dirichlet y Neumann que son las más cono-
cidas en los problemas de E.D.P., tanto por la riqueza de las propiedades de sus soluciones,
como también por la cantidad de aplicaciones conocidas; por ejemplo: Teorı́a del potencial
gravitatorio (o potencial newtoniano), potencial elástico, potencial armónico o bi-armóni-
co en elasticidad y muchos otros campos de la fı́sica matemática.
Una de las propiedades fundamentales de las soluciones de problemas de tipo elı́ptico; es
su grado de regularidad. Se muestra que toda solución continua de la ecuación de Laplace
es analı́tica, esto es, se puede representar en series de potencias convergentes, lo cual no
acontece con la ecuación de la onda.
140
6.10.2. La ecuación de Laplace. Campos conservativos, potencial electrosta-
tico
Estamos interesados en determinar funciones u = u(x, t) satisfaciendo la ecuación
uxx (x, y) + uyy (x, y) = 0
o más general, determinar u = u(x, y, z) tal que satisfaga la ecuación
uxx + uyy + uzz = 0
Estas ecuaciones son llamadas las ecuaciones de Laplace
141
Sea g : [0, 1] → Ω ⊂ R3 una parametrización diferenciable de la curva Γ, entonces por la
ley de Newton se sabe
→
−
F (x(t), y(t), z(t)) = m · g 00 (t) (5)
d2
donde g 00 (t) = 2 g(t), es la aceleración y g(t) = (x(t), y(t), z(t))
→
− dt
Luego si F es un campo conservativo, esto es
→
−
F = −(ux , uy , uz ) (6)
luego
md 0 2
mg 00 (t)g 0 (t) = |g (t)| (9)
2 dt
además
d ∂u ∂x ∂u ∂y ∂u ∂z
u(g(t)) = + +
dt ∂x ∂t ∂y ∂t ∂z ∂t
∂u ∂u ∂u ∂x ∂y ∂z
= , , · , ,
∂x ∂y ∂z ∂t ∂t ∂t
= grad u(g(t))g 0 (t)
esto es
∂
− u(g(t)) = −grad u(g(t))g 0 (t) (10)
∂t
luego de (10), (9) en (8) se obtiene
d md 0 2
− u(g(t)) = |g (t)| esto es
dt 2 dt
d m 0 2 d
|g (t)| + u(g(t)) = 0
dt 2 dt
d hm 0 2 i
|g (t)| + u(g(t)) = 0 (11)
dt 2
luego de (11) se obtiene que la función
m 0 2
|g (t)| + u(g(t)) es constante
2
m 0 2
La función C(t) = |g (t)| es denominado la energı́a cinética de la partı́cula en el
2
m 0 2
instante t, y la suma: E(t) = |g (t)| + u(g(t)) es llamado energı́a total de la curva Γ.
2
142
Luego la función u(x, y, z) es una función posición en Ω ⊂ R3 y es denominada función
→
−
potencial del campo F .
Ahora, nuestro objetivo es saber como encontrar campos que sean conservativos. Para
ello recordamos un resultado conocido en el calculo vectorial. Sea Ω ⊂ R3 una región
simplemente conexa (esto es toda región tal que al unir dos puntos cualesquiera de Ω
→
−
mediante un segmento, este debe estar en Ω). Sea el campo vectorial F : Ω ⊂ R3 → R3
→
−
continuamente diferenciable en Ω. Si la matriz jacobiana de F (DF ) es simétrica en Ω,
→
−
entonces F es un campo conservativo. Este teorema es valido en cualquier dimensión.
6.10.5. Ejemplo
Supongamos que en cierto punto P0 (x0 , y0 , z0 ) ∈ R3 se encuentra una masa ”m”, luego
por la ley gravitacional de Newton, se sabe que, sobre una masa ”m1 ” situado en un punto
p(x, y, z) actúa una fuerza de atracción
m · m1 x − x0 y − y0 z − z0
F (x, y, z) = −k , , (12)
d2 d d d
donde k es una constante gravitacional y ”d” es la distancia de P0 (x0 , y0 , z0 ) a P (x, y, z).
k · m · m1
Ası́ mismo la intensidad de F esta dado por: |F | =
d2
En particular, si escogemos un sistema de coordenadas tal que k = 1, m1 = 1 y
P0 (x0 , y0 , z0 ) = (0, 0, 0), se obtiene de (12) que
!
x y z
F (x, y, z) = −m ; ; (13)
(x2 + y 2 + z 2 )3/2 (x2 + y 2 + z 2 )3/2 (x2 + y 2 + z 2 )3/2
suponiendo que P (x, y, z) esta variando en una región simplemente conexa que no con-
tiene al origen, entonces se verifica que la matriz jacobiana de F en (13) es simétrica,
además de ello F es continuamente diferenciable, luego, del teorema del calculo mencio-
nado anteriormente se sigue que existe una función u(x, y, z) tal que
∆u = 0 (1)
Aunque la ecuación (1) aparecı́a por primera vez en un articulo de Euler sobre hidrodi-
namica en 1752, la ecuación quedó con el nombre de Laplace (1749-1827) que, a partir
de 1782, estudio exhaustivamente sus soluciones mientras investigaba la atracción gravi-
tacional entre cuerpos en el espacio.
143
El problema de Dirichlet consiste en lo siguiente:
Sea Ω ⊂ R2 el interior de un rectángulo en R2 , denotamos la frontera de Ω por ∂Ω.
Sea f : ∂Ω → R una función continua bien definida, entonces el problema de Dirichlet
para la ecuación de Laplace es, determinar una función u : Ω̄ → R tal que
∆u(x, y) = 0 en Ω
(2)
u(x, y) = 0 en ∂Ω
2
X ∂2
donde Ω̄ = Ω ∪ ∂Ω y ∆ = 2
es el operador laplaciano.
i=1
∂x i
Resolvemos el problema Dirichlet (2) mediante el método de Fourier (variables separables).
Comenzaremos resolviendo el problema en el caso que Ω = (0, a) × (0, b) y la función f
se anula en y = 0, y = b, y x = 0, esto es
uxx (x, y) + uyy (x, y) = 0 ∀(x, y) ∈ (0, a) × (0, b)
u(x, 0) = 0 = u(x, b), x ∈ [0, a] (3)
u(0, y) = 0, u(a, y) = f (y) y ∈ [0, b]
144
de la condición de frontera se obtiene:
u(0, y) = X(0)Y (y) = 0 ⇒ X(0) = 0
(6)
u(x, 0) = X(x)Y (0) = 0 ⇒ Y (0) = 0
luego de (4), (5) y (6) se obtiene dos E.D.O. con condición de frontera
X 00 (x) − λX(x) = 0 en (0, a)
(7)
X(0) = 0
y
145
imponiendo la condición de frontera no homogénea
+∞
X kπa kπy
f (y) = Dk senh sen = u(a, y), y ∈ [0, b]
k=1
b b
luego
∞
X kπa kπy
Dk senh sen = f (y)
k=1
b b
tomando
kπa
bk = Dk senh (12)
b
y usando la serie de Fourier de senos de f se tiene
2 b
Z
kπy
bk = f (t) sen dt, k∈N
b 0 b
Luego la expresión (14) es un candidato a solución del problema (3). Si es posible derivar
cada termino de la serie (14), entonces de hecho que u(x, y) es solución de la ecuación de
Laplace.
satisfaciendo
146
donde A es una constante positiva. El término Ax en la ecuación de la onda representa
una fuerza externa que en x tiene una magnitud Ax. En la ecuación consideramos c = 1.
Aplicando el método de Fourier para resolver el problema (1)-(4) se tiene
El procedimiento es elegir la función ψ(x) de una forma adecuada para reducir el problema
original.
Sustituyendo (6) en la E.D.P. (1) se obtiene
∂2 ∂2
(v(x, t) + ψ(x)) = (v(x, t) + ψ(x)) + Ax
∂t2 ∂x2
∂2 ∂2
v(x, t) = v(x, y) + ψ 00 (x) + Ax (7)
∂t2 ∂x2
la ecuación (7) se reducirá si se elige la función ψ de tal forma que
ψ 00 (x) + Ax = 0 (8)
Ax2
ψ 0 (x) = − + c luego una vez más
2
Ax3
ψ(x) = − + cx + D (9)
6
donde c y D son constantes de integración, y se puede elegir convenientemente.
Puesto que estamos asumiendo que la solución es u(x, y) = v(x, y) + ψ(x), entonces la
condición de frontera es de la forma
de otro lado, de (9) asumiremos que D = 0, luego ψ(0) = 0, reemplazando en (10) se tiene
147
de otro lado de la condición de frontera
AL3
pero de (9) se tiene ψ(L) = − + cL, reemplazando ello en (12) se tiene
6
AL3
u(L, t) = v(L, t) − + cL = 0 (13)
6
debemos reducir (13) a u(L, t) = v(L, t) = 0
AL3
para ello, debemos elegir una constante adecuada ”c” de tal forma que − + cL = 0
6
AL3 AL2
esto es = cL y c = .
6 6
Luego reemplazando en (9), la función ψ(x) debe ser de la forma
Ax3 AL2 x
ψ(x) = − + esto es
6 6
Ax 2
L − x2
ψ(x) = (14)
6
con esta elección de la función ψ en (14) se obtiene que ψ(L) = 0, luego (12) es de la
forma
u(L, t) = v(L, t) = 0 (15)
Ahora veamos que sucede con las condiciones iniciales de (6) se tiene
u(x, 0) = v(x, 0) + ψ(x) = 0; entonces de aquı́ se tiene v(x, 0) = −ψ(x) es decir de (14) se
tiene
Ax 2
x − L2
v(x, 0) = (16)
6
además
∂ ∂ ∂
u(x, 0) = v(x, 0) + 0 = 1, esto es v(x, 0) = 1 (17)
∂t ∂t ∂t
finalmente de (7), (8), (11), (15) y (17) se tiene el siguiente problema
∂ 2 v(x, t) ∂ 2 v(x, t)
− = 0; 0 < x < L; t > 0 (18)
∂t2 ∂x2
satisfaciendo
Luego el problema (18)-(21) es una ecuación de la onda homogénea con condición mixta
Ax 2
(inicial y frontera) donde f (x) = (x − L2 ) representa el desplazamiento inicial y
6
g(x) = 1, representa la velocidad inicial de un punto x de la cuerda en el instante t = 0,
148
respectivamente.
Ahora si, usando el método de Fourier se puede resolver el problema (18)-(21).
De lo estudiado anteriormente se sabe que la solución del problema reducido es de la
forma X nπx
v(x, t) = (Bn cos λn t + Cn sen λn t) sen (22)
L
donde
2 L
Z nπx cnπ
Bn = f (x) sen dx, λn = , c=1
L 0 L L
Z L
2 nπx nπ
Cn = g(x) sen dx λn =
λn c 0 L L
reemplazando f (x), g(x) e integrando se obtiene
2 L Ax 2
Z nπx
x − L2 sen
Bn = dx
L 0 6 L
2AL3 (−1)n
Bn = , también
π3 n3
Z L
2 nπx
Cn = (1) sen dx
nπ 0 L
2L
Cn = [1 − (−1)n ]
(nπ)2
satisfaciendo
149
(1). Problema
Resolver el problema mixto
satisfaciendo
150