Examen Final EDP
Universidad de Guanajuato, División de Ciencias e Ingenierías
Ecuaciones Diferenciales Parciales
Joseph Leonel Manzo Bustos
14 de diciembre de 2022
1. Resolver la ecuación de calor siguiente:
∂u ∂2u
= + cos (πx)
∂t ∂x2
0 < x < 1, 0<t<∞
u(x = 0, t) = 1
u(x = 1, t) = −1
u(x, t = 0) = 0
La ecuación es no homogenea, por lo tanto, estará compuesta de una solución
particular y una solución asociada a la ecuación homegénea.
u(x, t) = u0 (x, t) + upart (x)
La inhomogeneidad de nuestra ecuación depende de x, por lo tanto, tenemos
que encontrar ua solución particular que solamente dependa de x.
∂ 2 upart
+ cos (πx) = 0
∂x2
Para obtener la solución particular, integramos dos veces respecto a x.
∂upart 1
= − sin (πx) − A
∂x π
1
Ecuaciones Diferenciales Parciales
1
upart (x) = cos (πx) − Ax + B
π2
Teniendo la solución, aplicamos las condiciones de frontera de u(x, t)
u(x = 0, t) = 1 → upart (0) = 1
u(x = 1, t) = −1 → upart (1) = −1
usando la primera condición frontera
1
upart (0) = +B =1
π2
π 2 −1
Por lo tanto B tiene que ser π2
. Aplicando la segunda condición de frontera
1 π2 − 1
upart (1) = − − A + = −1
π2 π2
2
Entonces A tiene que cumplir en ser 2(1−π
π2
)
.
De esta forma, la solución particular de nuestra ecuación en derivadas par-
ciales es
1 2(1 − π 2 ) π2 − 1
upart (x) = 2 cos (πx) − x + (1)
π π2 π2
A partir de nuestra solución particular, buscaremos las condiciones de fron-
tera para nuestra solución homogénea
u0 (x = 0, t) = u(x = 0, t) − upart (0) → 1−1=0
u0 (x = 1, t) = u(x = 1, t) − upart (1) → −1 + 1 = 0
1 2(1 − π 2 ) π2 − 1
u0 (x, 0) = u(x, t = 0) − upart (x) → cos−(πx) + x −
π2 π2 π2
Teniendo las condiciones iniciales y de frontera, proponemos una solución
separable para la homegénea u0 (x, t) = X(x)T (t)
∂u0 ∂ 2 u0
=
∂t ∂x2
Separando variables
1 ∂T 1 ∂2X
=
T ∂t X ∂x2
La unica forma en que esta ecuación se cumpla, es si las dos partes de la
ecuación son iguales a una constante, por lo tanto
1 ∂T 1 ∂2X
= = α2
T ∂t X ∂x2
2
Ecuaciones Diferenciales Parciales
1 ∂T 1 ∂2X
= α2 = α2
T ∂t X ∂x2
Tomando todos los posibles casos donde α2 = 0, α2 < 0 y α2 > 0.
El caso que cumple con las ecuaciones de frontera es α2 < 0, por lo tanto, la
solución para X(x) y T(t) será
X(x) = A sin (αx) + B cos (αx)
2t
T (t) = Dn e−α
usando las condiciones de frontera para X(x)
X(0) = B = 0 X(1) = A sin (α) = 0
Para que esta relación se cumpla α = nπ ∴
X(x) = An sin (nπx)
Tomando la propuesta inicial para la separación de variables, u0 (x, t) =
X(x)T (t), An Dn = Bn
n=1
2 π2 t
X
u0 (x, t) = [Bn sin (nπx)] e−n (2)
∞
Usando la condición inicial
n=1
1 2(1 − π 2 ) π2 − 1 X
u0 (x, t = 0) = − cos (πx) − x − = [Bn sin (nπx)]
π2 π2 π2 ∞
De (2) podemos ver que las constantes Bn son los coeficientes de Fourier para
la función u0 (x, t). Para calcular estos coeficientes hay que recordar que
< en |f >
cn =
||en ||2
donde, en este caso, cn = Bn y en = sin (nπx), ahora calculando los coefi-
cientes
Z 1
2(1 − π 2 ) π2 − 1
1
Bn = 2 − 2 cos (πx) + x− sin (nπx)
0 π π2 π2
Realizando el calculo de la integral
((4π − 3π 3 )n3 + (3π 3 − 3π)n)(−1)n + π 3 n3 + (π − π 3 )n
Bn = −2
π 4 n2 (n2 − 1)
3
Ecuaciones Diferenciales Parciales
Ahora, sustituyendo el coeficiente en la solución de la homogénea
n=1
X ((4π − 3π 3 )n3 + (3π 3 − 3π)n)(−1)n + π 3 n3 + (π − π 3 )n
−n2 π 2 t
u0 (x, t) = −2 sin (nπx)e
∞
π 4 n2 (n2 − 1)
Como último paso sumamos las soluciones para obtener la solución general.
1 2(1 − π 2 ) π2 − 1
u(x, t) = cos (πx) − x +
π2 π2 π2
n=1
X ((4π − 3π 3 )n3 + (3π 3 − 3π)n)(−1)n + π 3 n3 + (π − π 3 )n
−n2 π 2 t
+ −2 sin (nπx)e
∞
π 4 n2 (n2 − 1)
(3)
