Métodos II, curso 2021–2022 Grupos D y E.

SOLUCIONES EXAMEN FINAL, 24 de junio de 2022

Resolver dos problemas entre I, II, III y tres problemas entre IV, V, VI, VII.

I Sea la ecuación: xux + (x + 2y)uy = 2x + 2y.

(1) Calcule la solución general.
(2) Calcule una solución particular que verifique: u(1, y) = 0.
(3) Utilizando el Teorema de existencia y unicidad para este tipo de ecuaciones, discuta
la unicidad de la solución hallada en el apartado b).

Resolución.
(1) Las ecuaciones de Lagrange-Charpit son
dx dy du
= = .
x x + 2y 2y + 2x
Las caracterı́sticas son dadas por:
dx dy dy 2y
= , = +1 .
x x + 2y dx x
La ecuación homogénea tiene solución yhom = cx2 . Utilizando el método de variación de
constantes, o bien la fórmula Z
2 dx
yp (x) = x
x2
obtenemos la solución particular yp (x) = −x, ası́ que y = cx2 − x. Entonces, un invariante
es dado por
y+x
ξ= .
x2
Por tanto, la solución general es dada por las relaciones:
dx du dx du du
= ⇐⇒ = 2
⇐⇒ 2xdx = ⇐⇒ u = ξx2 + K(ξ)
x 2y + 2x x 2ξx ξ
es decir:  
y+x
u(x, y) = y + x + K
x2
(2) La solución particular es:
u(1, y) = 1 + y + K (1 + y) = 0, K(1 + y) = −(1 + y), K(w) = −w
 
x+y 1
u(x, y) = x + y − 2
= (x + y) 1 − 2
x x
(3) De acuerdo con el teorema de existencia y unicidad, como los coeficientes de la
ecuación:
a(x, y)ux + b(x, y)uy = c(x, y)u + d(x, y)
1
 2

es decir a(x, y) = x, b(x, y) = x + 2y, c(x, y) = 0, d(x, y) = 2y + 2x, son suficiente-

mente regulares (analı́ticos reales), se trata de ver si la curva de datos es tangente a las
caracterı́sticas. Como la curva de datos es γ(s) = (1, s), la condición de transversalidad
es:
(1, 2s + 1) · (1, 0) = 1 6= 0
luego la solución es única en un entorno de la curva de datos.

II Sea la ecuación en derivadas parciales (EDP)

4uxx + 10uxy + 4uyy + (y − 2x)2 (ux + 2uy ) = 0.
(1) Clasifique la EDP.
(2) Reduzca la EDP a forma normal mediante un cambio de variable.
(3) Resolviendo la ecuación reducida, halle la solución general de la EDP.

Resolución.
(1) Como a = 4, b = 5, c = 4, tenemos que δ = b2 − ac = 9 > 0, y por tanto la ecuación
es hiperbólica.
(2) El cambio de variables que reduce la ecuación a forma canónica es
(
ξ = y − x2
η = y − 2x
Obtenemos
1
ux = − uξ − 2uη , uy = uξ + uη
2

1
uxx = 4 uξξ + 2uξη + 4uηη ,

uyy = uξξ + 2uξη + uηη ,
uxy = − 12 uξξ − 25 uξη − 2uηη



La EDP original se convierte en la ecuación

6uξη − η 2 uξ = 0 .
(3) Deducimos
η3
uξ = g(ξ, η), 6gη = η 2 g, g(ξ, η) = K(ξ)e 18
e integrando obtenemos
3 /18
u(ξ, η) = A(ξ)eη + B(η).
La solución general de la ecuación de partida es
3 /18
u(x, y) = A(2y − x)e(y−2x) + B(y − 2x).
 3

III Sea la ecuación

(E) : x2 y 00 + x (x − 3) y 0 + 4y = 0.
(1) Calcule ai , i = 1, 2, 3, los coeficientes de los primeros términos del desarrollo en
serie:
y1 (x) = xr (1 + a1 x + a2 x2 + a3 x3 + · · · )
y el exponente r de una solución de la ecuación (E).
(2) Demuestre que todas las soluciones de la ecuación (E) están acotadas en el origen.

Resolución.
(1) x = 0 es un punto singular:
x−3 4
p(x) =, q(x) = 2
x x
p(x) no es analı́tica en x = 0. Pero, xp(x) = x − 3, x2 q(x) = 4, son analı́ticas en x = 0
que es singular regular. Entonces:
p0 = lim xp(x) = −3, q0 = lim x2 q(x) = 4 .
x→0 x→0
La ecuación indicial es:
r(r − 1) − 3r + 4 = 0, r2 − 4r + 4 = 0, r1 = 2, doble
Según el teorema de Frobenius, existe una solución:
∞
X
2
an xn = x2 a0 + a1 x + a2 x2 + a3 x3 + · · ·


y(x) = x
n=0
de esta forma, r = 2 y a0 = 1. Sustituyendo en la ecuación:
∞
X ∞
X ∞
X
x2 (n + 2)(n + 1)an xn + x(x − 3) (n + 2)an xn+1 + 4 an xn+2 = 0
n=0 n=0 n=0
∞
X ∞
X ∞
X ∞
X
(n + 2)(n + 1)an xn+2 + (n + 2)an xn+3 − 3(n + 2)an xn+2 + 4an xn+2 = 0
n=0 n=0 n=0 n=0
∞
X ∞
X
[(n + 2)(n + 1) − 3(n + 2) + 4]an xn+2 + (n + 2)an xn+3 = 0
n=0 n=0
∞
X ∞
X ∞
X ∞
X
n2 an xn+2 + (n + 2)an xn+3 = 0, n2 an xn+2 + (n + 1)an−1 xn+2 = 0
n=0 n=0 n=0 n=1
∞
X  2
n an + (n + 1)an−1 xn+2 = 0

n=2
n+1
a0 = 1, n2 an + (n + 1)an−1 = 0, an = − an−1 , n = 1, 2, 3, . . .
n2
 4

3 3 2
a1 = −2a0 = −2, a2 = − a1 = , a3 = −
4 2 3
y la solución es:
 
3 2
y1 (x) = x2 1 − 2x + x2 − x3 + · · ·
2 3
(2) Según el teorema de Frobenius la solución general es
∞
!
X
y(x) = c1 y1 (x) + c2 y1 (x) log x + x2 bn xn
n=1

que está acotada en x = 0 pues

∞
2

X
y1 (x) log x = x ln x an xn , lim x2 ln x = 0
x→0
n=0

IV Calcule la solución del problema de Cauchy:


9utt + 12utx + 4uxx = 0,
 x ∈ R, t > 0,
u(x, 0) = 9x2 ,

ut (x, 0) = −12x,


utilizando transformadas de Fourier en x.

Resolución.

La transformada de Fourier en x de la ecuación nos da:

9ûtt + 12ikût − 4k 2 û = 0 .
El polinomio caracterı́stico y sus raı́ces son:
√
2 2 −144k 2 + 144k 2
−12ik ± 2ik
9λ + 12ikλ − 4k = 0, λ= =− , doble
18 3
Por tanto, la solución general de la ecuación diferencial es:
 
û(k, t) = Â(k) + t B(k) e−2ikt/3 .

Aplicando la transformada inversa obtenemos:

   
2t 2t
u(x, t) = A x − + tB x −
3 3
Imponiendo las condiciones iniciales, obtenemos:
2 2
u(x, 0) = A(x) = f (x), ut (x, 0) = − A0 (x) + B(x) = g(x), B(x) = g(x) + f 0 (x)
3 3
 5

y la solución general es:

      
2t 2t 2 2t
u(x, t) = f x− +t g x− + f0 x −
3 3 3 3
Si f (x) = 9x2 , g(x) = −12x:
2t 2
      
2t 2t
u(x, t) = 9 x − + t −12 x − + 12 x − = (3x − 2t)2 .
3 3 3

V Sea el problema de contorno

u00 + λu = 0, u(0) + u0 (0) = 0, u(2π) + u0 (2π) = 0, x ∈ (0, 2π)
(1) Calcule los autovalores y las autofunciones del problema.
(2) Estudie la existencia y unicidad de soluciones del problema no homogéneo u00 −u = e2x .

Resolución.
(1) Distinguimos tres casos.
(i) λ > 0 √ √
u(x) = A cos λ x + Bsen λ x
Imponiendo las condiciones de contorno:
√
u(0) + u0 (0) = A + B λ = 0
√ √ √ √ √ √
u(2π) + u0 (2π) = A cos 2 λ π + Bsen2 λ π − λ Asen2 λ π + λ B cos 2 λ π
 √ √ √   √ √ √ 
= A cos 2 λ π − λ sen2 λ π + B sen2 λ π + λ cos 2 λ π
 √ 
√ 1√ √ √ √ λ √
det
cos 2 λ π − λ sen2 λ π sen2 λ π + λ cos 2 λ π
√ √ √ √  √ √ √ 
= sen2 λ π + λ cos 2 λ π − λ cos 2 λ π − λ sen2 λ π
√ n2
(1 + λ)sen2 λ π = 0, λ= , n = 1, 2, . . .
4
Las autofunciones son:
n n n nx nx
A+B = 0, A = − B, un (x) = cos − sen
2 2 2 2 2
(ii) λ = 0:
u(x) = Ax + B
Imponiendo las condiciones de contorno:
u(0) + u0 (0) = B + A = 0, u(2π) + u0 (2π) = 2πA + B + A = 0
 6

luego λ = 0 no es
√ autovalor.
(iii) λ < 0, µ = −λ > 0
u(x) = A cosh µ x + B sinh µ x
Imponiendo las condiciones de contorno:
u(0) + u0 (0) = A + B µ = 0

u(2π) + u0 (2π) = A cosh 2µ π + B sinh 2µ π + µ A sinh 2µ π + µ B cosh 2µ π

= A (cosh 2µ π + µ sinh 2µ π) + B (sinh 2µ π + µ cosh 2µ π)
 
1 µ
det
cosh 2µ π + µ sinh 2µ π sinh 2µ π + µ cosh 2µ π
= sinh 2µ π + µ cosh 2µ π − µ (cosh 2µ π + µ sinh 2µ π)

(1 − µ2 ) sinh 2µ π = 0, µ = 1, λ = −1
La autofunción es:
A + B = 0, A = −B, u−1 (x) = cosh x − sinh x = e−x
que verifica la condición u0 + u = 0 en todo punto (la ecuación de segundo orden se
factoriza).
(2) Escribiendo el problema en la forma Sturm-Liouville Ly = µy + F (x), es decir −y 00 =
−y+e2x , observamos que, como λ = −1 es autovalor del problema homogéneo, el problema
no homogéneo admite infinitas soluciones o ninguna. Considerando que
Z 2π
e−x e2x dx = e2π − 1 6= 0
0
entonces el problema no admite ninguna solución.

VI Mediante separación de variables, resuelva el problema mixto siguiente para la ecuación

de ondas: 

 utt − uxx = 0, x ∈ (0, π), t > 0, k ∈ R\{0}

u(x, 0) = x,

 ut (x, 0) = k cos x,

ux (0, t) = ux (π, t) = 0.


Resolución.
Separando variables en la forma u(x, t) = X(x)T (t), obtenemos para la función X(x)
el problema de Sturm-Liouville homogéneo, regular y simétrico
(
X 00 + λX = 0,
X 0 (0) = X 0 (π) = 0,
 7

con condiciones de tipo Neumann. Se comprueba fácilmente que su solución es

λ n = n2 , Xn (x) = cos nx, n = 0, 1, 2, . . .
La solución que buscamos tendrá por tanto la forma general
∞
X
u(x, t) = Tn (t) cos nx.
n=0

Sustituyendo esta expresión para u(x, t) en la ecuación de ondas, tenemos

∞
X
(Tn00 (t) + n2 Tn (t)) cos nx = 0
n=0

La condición u(x, 0) = x nos da

∞ ∞
X a0 X
Tn (0) cos nx = x = + an cos nx
2
n=0 n=1
con Z π Z π
2 2 4
a0 = xdx = π, an = x cos nxdx = −
π 0 π 0 π(2n − 1)2
π 4 P∞ cos(2n−1)x
Por tanto x = 2 − π n=1 (2n−1)2 , ası́ que
∞ ∞
X π 4 X cos(2n − 1)x
Tn (0) cos nx = −
2 π (2n − 1)2
n=0 n=1

La condición ut (x, 0) = k cos x se escribe como

∞
X
Tn0 (0) cos nx = k cos x .
n=0

Compaginando los resultados anteriores, obtenemos los problemas de Cauchy siguientes:


00
T0 (t) = 0,

π
T0 (0) = π2 ⇐⇒ T0 (t) =
 0
 2
T0 (0) = 0,

00
T1 (t) + T1 (t) = 0,

4
T1 (0) = − π4 ⇐⇒ T1 (t) = − cos t + k sin t
 0
 π
T1 (0) = k,

00 2
T2n (t) + 4n T2n (t) = 0,

T2n (0) = 0 ⇐⇒ T2n (t) = 0 n = 1, 2, . . .
 0

T2n (0) = 0,
 8


00 2
T2n−1 (t) + (2n − 1) T2n−1 (t) = 0,

4 cos(2n − 1)t
4
T2n−1 (0) = − π(2n−1) 2 ⇐⇒ T2n−1 (t) = − n = 2, 3, . . .

 0 π(2n − 1)2
T2n−1 (0) = 0,
En definitiva, tenemos
  ∞
π 4 4 X cos[(2n − 1)t] cos[(2n − 1)x]
u(x, t) = + − cos t + k sin t cos x −
2 π π (2n − 1)2
n=2

VII Calcule la solución del siguiente problema en la corona circular de radios 1 y 2:

∆u = 0, u(1, ϕ) = 1, ur (2, ϕ) = senϕ, ϕ ∈ [0, 2π]

Resolución.
En polares:
1 1
urr + ur + 2 uϕϕ = 0
r r
Si u(r, ϕ) = R(r)Φ(ϕ):
r2 R00 + rR0 Φ00
=− =λ
R Φ
La ecuación en el ángulo ϕ es (condiciones periódicas):
Φ00 + λΦ = 0, Φ(0) = Φ(2π), Φ0 (0) = Φ0 (2π)
λ = 0 es un autovalor, con autofunción 1 y los autovalores positivos son λn = n2 con
autofunciones cos nx, sennx.
∞
X
u(r, ϕ) = u0 (r) + (un (r) cos nϕ + vn (r)sennϕ)
n=1

Sustituyendo en la ecuación:
∞ 
n2 n2
   
00 1 0 X 00 1 0 00 1 0
u0 + u0 + un + un − 2 un cos nϕ + vn + vn − 2 vn sennϕ = 0
r r r r r
n=1

Tenemos las ecuaciones:

1
u000 + u00 = 0
r
1 n2 1 n2
u00n + u0n − 2 un = 0, vn00 + vn0 − 2 vn = 0, n = 1, 2, . . .
r r r r
con soluciones:
Bn Dn
u0 (r) = A0 + B0 log r, un (r) = An rn + , vn (r) = Cn rn + , n = 1, 2, . . .
rn rn
 9

Las condiciones de contorno son:

∞
X
u(1, ϕ) = u0 (1) + (un (1) cos nϕ + vn (1)sennϕ) = 1
n=1
∞
X
ur (2, ϕ) = u00 (2) + u0n (2) cos nϕ + vn (2)0 sennϕ = senϕ


n=1
es decir:
u0 (1) = 1, un (1) = 0, vn (1) = 0
u00 (2) = 0, u0n (2) = 0, v10 (2) = 1, vn0 (2) = 0, n 6= 1
Cuando las condiciones son nulas, las soluciones son 0. Las soluciones no triviales son:
u0 (r) =A0 + B0 log r, u0 (1) = 1, u00 (2) = 0
D1
v1 (r) =C1 r + , v1 (1) = 0, v10 (2) = 1
r
y las soluciones:
 
B0 D1 4 1
A0 = 1, = 0, u0 (r) = 1; C1 + D1 = 0, C1 − = 1, v1 (r) = r−
2 4 5 r
De esta forma:  
4 1
u(r, ϕ) = 1 + r− senϕ
5 r