Análisis de Fourier y EDP – Grado en Matemáticas: Prueba de Evaluación Continua (PEC)

 
27 de Abril 2024

Justifique todas sus respuestas. Las respuestas sin justificar recibirán una calificación de 0 puntos.

* Problema 1 (3 puntos): Considere la ecuación

xux + yuy + u = 1.

Para cada uno de los datos iniciales

1 ii.) u(x, 2x) = 3x, iii.) u(1, y) = m(y),
i.) u(x, x) = 1 + ,
x

comprobar la existencia y unicidad de soluciones del correspondiente problema de Cauchy. Si tiene

solución, encontrarla. En caso de tener multiples soluciones, explicar por qué y proporcionar dos
soluciones distintas. En caso de no tener solución, explicar por qué.
Solución: Las curvas caracterı́sticas vienen dadas por el sistema
dx


 = x,


 dτ
dy

= y,

 dτ
 du = 1 − u,



dτ
o, de manera equivalente,
dx dy du
= = .
x y 1−u
De la primera igualdad obtenemos
dy dx
= ⇒ y = ξx,
y x
mientras que de
du dx h(ξ)
= ⇒ u=1− ,
1−u x x
para una función arbitraria h ∈ C 1 (R) Ası́, la solución general de la EDP es

h xy


u(x, y) = 1 − . (1)
x
Estudiamos ahora el problema de Cauchy correspondiente a cada uno de los datos iniciales dados

i.) La curva dato Γ(s) = (γ1 (s), γ2 (s), γ3 (s)) se parametriza como
 
1
Γ(s) = s, s, 1 +
s

y la matriz asociada, digamos

 
a(γ1 (s), γ2 (s), γ3 (s)) b(γ1 (s), γ2 (s), γ3 (s)) c(γ1 (s), γ2 (s), γ3 (s))
A(s) =
γ1′ (s) γ2′ (s) γ3′ (s)
 viene entonces dada por
− 1s
!
s s
A(s) =
1 1 − s12
La matriz A(s) tiene rango 1 y, por tanto, el problema de Cauchy tiene infinitas soluciones,
una para cada h ∈ C 1 (R) con h(1) = −1, ya que el dato se proporciona sobre una curva
caracterı́stica. En efecto, de la solución general (1), tenemos

h(1) 1
u(x, x) = 1 − =1+ ⇒ h(1) = −1.
x x
Por ejemplo, tomando h(t) = −1, obtenemos la solución
1
u(x, y) = 1 + ,
x
y, para h(t) = −t2 , obtenemos
y2
u(x, y) = 1 + .
x3
ii.) En este caso la curva dato y la matriz asociada son
!
s 2s 1 − 3s
Γ(s) = (s, 2s, 3s) y A(s) =
1 2 3

Ahora bien,

s 2s s 1 − 3s 2s 1 − 3s
∆= =0 pero ∆1 = ̸= 0 y ∆2 = ̸= 0,
1 2 1 3 2 3

luego no se cumple la condición de transversalidad y el problema de Cauchy no tiene solución.

De hecho, de la solución general (1), tenemos

h(2)
u(x, 2x) = 1 − = 3x ⇒ h(2) = x − 3x2
x
lo cual es imposible pues h(2) es constante.
iii.) En este caso la curva dato y la matriz asociada son
!
1 s 1 − m(s)
Γ(s) = (1, s, m(s)) y A(s) =
0 1 m′ (s)

Ahora bien,
1 s
∆= ̸= 0,
0 1
por lo que se cumple la condición de transversalidad y el problema de Cauchy tiene solución
única. De la solución general (1), obtenemos

u(1, y) = 1 − h(y) = m(y) ⇒ h(·) = 1 − m(·),

y, por tanto, la solución del problema de Cauchy es

y

1−m x
u(x, y) = 1 − .
x

* Problema 2 (2+1 puntos): Calcule los autovalores y autofunciones del problema



 −y ′′ = λy, x ∈ (0, 1) ,
y(0) + y ′ (0) = 0, (Eλ )

y(1) = 0.

 Analice para qué valores de λ ∈ R el problema


 −y ′′ = λy + 1 − x, x ∈ (0, 1) ,
y(0) + y ′ (0) = 0, (Pλ )

y(1) = 0,


tiene solución y cuándo la solución es única. Indicación: Utilice el Teorema de la alternativa de Fredholm.
√
Solución: El polinomio caracterı́stico de (Eλ ) es p(r) = r2 + λ, cuyas raı́ces son r± = ± −λ.
Distinguimos los siguientes casos
λ = 0: La solución general de y ′′ = 0 es y(x) = c1 + c2 x. Teniendo en cuenta las condiciones
de frontera obtenemos la condición c1 + c2 = 0. Ası́, λ = 0 es un autovalor asociado a la
autofunción ϕ0 (x) = 1 − x.
λ < 0: Tomemos1 λ = −µ2 . La solución general es en este caso y(x) = c1 eµx + c2 e−µx . Las
condiciones de frontera nos dan el sistema de ecuaciones
(
c1 (1 + µ) + c2 (1 − µ) = 0,
c1 + c2 = 0,
cuya única solución (nótese que µ ̸= 0) es la solución trivial c1 = c2 = 0, lo que produce la
solución trivial y(x) = 0. Luego concluimos que no hay2 autovalores λ < 0.
λ > 0: Tomando λ = µ2 , el polinomio caracterı́stico tiene en este caso raı́ces complejas
r± = ±iµ. La solución general es en este caso y(x) = c1 sen(µx) + c2 cos(µx). Las condiciones
de frontera nos dan el sistema de ecuaciones
(
c2 + c1 µ = 0,
(2)
c1 sen(µ) + c2 cos(µ) = 0,
Despejando c2 de la primera ecuación y sustituyendo en la segunda llegamos a la ecuación
tan(µ) = µ (3)
La ecuación anterior tiene un numero infinito de soluciones µn , n = 1, 2, . . . las cuales se
aproximan conforme crece n a los múltiplos impares de π2 , esto es µn ≈ (2n+1) π2 , n = 1, 2, . . ..
Más aún, dado que la función tan(·) es una función impar, esto es, la ecuación (3) permanece
invariante bajo el cambio µ → −µ, las únicas soluciones no triviales de la ecuación (3) son
los valores µn ≈ ±(2n + 1) π2 , n = 1, 2, . . .. Por tanto, el problema (Eλ ), tiene una sucesión
infinita de autovalores
π2
λn = µ2n ≈ (2n + 1)2 , n = 1, 2, . . . (4)
4
Teniendo en cuenta que, de la segunda ecuación de (2), tenemos c2 = −c1 sen(µ)
cos(µ) , obtenemos

yn (x) = c1 sen(µn x) + c2 cos(µn x)

sen(µ)
= c1 sen(µn x) − c1 cos(µn x)
cos(µ)
c1
= (sen(µn x) cos(µn ) − sen(µ) cos(µn x))
cos(µn )
c1
= sen(µn (1 − x))
cos(µn )
√
Ası́, a cada autovalor λn le corresponde la autofunción ϕn (x) = sen( λn (1 − x)).
Concluimos entonces que el problema (Eλ ) tiene como autofunciones
p
ϕ0 (x) = 1 − x, ϕn (x) = sen( λn (1 − x)),
asociadas a los autovalores λ0 = 0 y λn dados por (4), respectivamente.
La existencia y unicidad de soluciones de (Pλ ) es una aplicación directa del Teorema de la
alternativa de Fredholm. En particular, notando que el término no homogéneo f (x) = 1 − x
coincide con la autofunción ϕ0 (x), tenemos:
1
Por el Teorema 3.30, 2ª ed. texto base (4.31 1ª ed.) sabemos que λ ≥ 0, luego no harı́a falta considerar el caso λ < 0.
2
Como ya habı́amos anunciado.
 a) Si λ ̸= 0 y λ ̸= λn , (esto es, λ no es autovalor del problema (Eλ )) entonces el problema (Pλ )
tiene solución única.
b) Si λ = λn con λ ̸= 0: por el Teorema 3.18, 2ª ed. del texto base (Teorema 4.19, 1ª ed.) la
autofunciones de un problema de Sturm-Liouville son ortogonales y dado que f (x) = 1 − x,
luego f (x) = ϕ0 (x), se satisface la condición de ortogonalidad del Teorema de la alternativa
de Fredholm y, por tanto, el problema (Pλ ) tiene infinitas soluciones de la forma
1 p
u(x) = (1 − x) + c sen( λn (1 − x)), con c ∈ R.
λn
c) Si λ = 0, entonces el problema no tiene solución pues no se satisface la condición de
ortogonalidad del Teorema de la alternativa de Fredholm: f (x) no es ortogonal a ϕ0 (x),
de hecho
1
⟨f, ϕ0 ⟩L2 (0,1) = ⟨ϕ0 , ϕ0 ⟩L2 (0,1) = ∥ϕ0 ∥2L2 (0,1) = .
3
En este caso, podemos llegar a la misma conclusión sin invocar el Teorema de la alternativa
de Fredholm. En concreto, para λ = 0, tenemos


 −y ′′ = 1 − x, x ∈ (0, 1) ,
y(0) + y ′ (0) = 0, (P0 )

y(1) = 0,

3 x2
Integrando la ecuación −y ′′ = 1 − x, obtenemos y(x) = x6 − 2 + c1 x + c2 . Las condiciones
de frontera producen entonces el sistema
(
c2 + c1 = 0,
1 1
6 − 2 + c1 + c2 = 0,
1 1
que implica 6 − 2 = 0 lo cual, claramente, carece de sentido.

* Problem 3 (2+2 puntos):

1
a) Calcule la serie de Fourier de f (x) = x + x2 para −π < x < π. Utilice el resultado para
4
probar, justificadamente, que
∞
X 1 π2
2
= .
n 6
n=1
b) Considere la función 
0 si − π < x < 0,
f (x) = √
x si 0 ≤ x < π.
Pruebe que
∞ Z π π
5π 3
X  Z 
f (x) cos(nx)dx f (x) sin(nx)dx ≤ .
−π −π 36
n=1

Solución:
a) El cálculo directo de los coeficientes de la serie de Fourier nos da
∞
1 π 2 X (−1)n
x + x2 ∼ + (cos(nx) − 2n sen(nx)), π < x < π.
4 12 n2
n=1

La función f (x) = x + 14 x2 y la aproximación dada por los 15 primeros términos de su serie de Fourier.
 Dado que f (x) = x + 14 x2 es una función C 1 (−π, π), por el Teorema 4.51, 2ª ed. del texto
base (Teorema 5.48, 1ª ed.) tenemos

2 ∞ n  f (x)
 si − π < x < π,
π X (−1)
+ (cos(nx) − 2n sen(nx)) = f (π − ) + f (−π + )
12 n2  si x = ±π,
n=1 
2
por tanto, para x = π, tenemos
∞
π 2 X (−1)n f (π − ) + f (−π + )
+ (cos(nπ) − 2n sen(nπ)) = ,
12 n2 2
n=1

esto es
∞
π + 14 π 2 + −π + 41 π 2



π 2 X (−1)n n
+ ((−1) − 0) = ,
12 n2 2
n=1
de donde
∞
X 1 π2 π2 π2
= − = .
n2 4 12 6
n=1

b) Sea
∞ Z
X π  Z π 
S := f (x) cos(nx)dx f (x) sin(nx)dx .
n=1 −π −π

El resultado se sigue usando la desigualdad de Young

a2n b2n
an bn ≤ + ,
2 2
junto con la identidad de Parceval
π
a20 X 2
Z
1
+ (an + b2n ) = f 2 (x)dx,
2 π −π
n∈N

donde
Z π Z π
1 1
an = f (x) cos(nx)dx, n≥0 y bn = f (x)sen(nx)dx, n ≥ 1.
π −π π −π

De este modo
π2 X 2 π2 1 π 2 a20 π2 π 2 5π 3
X  Z   
2 2
S=π an bn ≤ (an + bn ) = f (x)dx − = − π = .
2 2 π −π 2 2 2 9 36
n∈N n∈N