Universidad Abierta y a Distancia de México

Ecuaciones Diferenciales Parciales

Dr. Luis Fernando Orozco Cortés

Unidad 2. El problema de Cauchy para EDP de

segundo orden
Evidencia de Apredizaje. Análisis al planteamiento de
problemas

Julio Cueva Lozano

Es 1821000315
Resuelva los siguientes problemas de Cauchy

1.- Resuelva el siguiente problema de dimensión uno

𝑢𝑡𝑡 (𝑥, 𝑡) − 𝑢𝑥𝑥 (𝑥, 𝑡) = 𝑥𝑒 −𝑡−1 𝑡 > 0
𝑠𝑖 𝑢(𝑥, 0) = cos2 𝑥 𝑦 𝑢𝑡 (𝑥, 0) = 𝑒 −𝑥

El problema anterior de cauchy para la ecuación no homogénea esta dada por

𝑢𝑡𝑡 = 𝑎2 𝑢𝑥𝑥 + 𝑓(𝑥, 𝑡)

Sujeta a las condiciones
𝑢(𝑥, 0) = 𝜙(𝑥) 𝑢𝑡 (𝑥, 0) = 𝜑(𝑥)

Resolviendo por el método de las características entonces tenemos que

𝜉 = 𝑥 + 𝑎𝑡 𝜂 = 𝑥 − 𝑎𝑡
Lo cual nos lleva a la forma
𝜉+𝜂 𝜉−𝜂 1 𝜉+𝜂 𝜉−𝜂
𝑢𝜉𝜂 ( , ) = − 2𝑓( , )
2 2𝑎 4𝑎 2 2𝑎
𝜉
Y al integrar con respecto a 𝜂 tenemos
𝜉+𝜂 𝜉−𝜂 1 𝜉+𝜂 𝜉−𝜂
𝑢𝜉𝜂 ( , ) = 𝑢𝜂 (𝜂, 0) − 2 ∫ 𝑓 ( , ) 𝑑𝜉
2 2𝑎 4𝑎 2 2𝑎
𝜂
Y ahora, caso contrario al integra con respecto a
𝜉
1 1
𝑢𝜂 (𝜂, 0) = 𝑢𝑥 (𝜂, 0) − 𝑢 (𝜂, 0)
2 2𝑎 𝑡
Lo que entonces obtenemos

𝜉+𝜂 𝜉−𝜂
𝑢(𝜉, 0) − 𝑢 ( , )=
2 2𝑎
𝜉
1 1 1 𝜉 𝜉 𝜉 ̅ + 𝜂̅ 𝜉 ̅ − 𝜂̅
∫ ( 𝑢𝑥 (𝜂̅ , 0) − 𝑢𝑡 (𝜂̅ , 0)𝑑𝜂̅ − 2 ∫ ∫ 𝑓( , ) 𝑑𝜂̅
𝜂 2 2𝑎 4𝑎 𝜂 𝜂 2 2𝑎

Donde tenemos que

𝜉 𝜉 𝜉
1 1
∫ ( 𝑢𝑥 (𝜂̅ , 0) − 𝑢𝑡 (𝜂̅ , 0)) 𝑑𝜂̅ = ∫ 𝑢𝑥 (𝜂̅ , 0)𝑑𝜂̅ = ∫ 𝜙 ′ (𝜂̅ )𝑑𝜂̅ = 𝜙(𝜉) − 𝜙(𝜂)
𝜂 2 2𝑎 𝜂 𝜂

Resolviendo la segunda integral volveremos a las variables transformando y

obteniendo como solución
1 1 𝑥+𝑎𝑡 1 𝑡 𝑥+𝑎(𝑡−𝜏)
𝑢(𝑥, 𝑡) = {𝜙(𝑥 − 𝑎𝑡) − 𝜙(𝑥 + 𝑎𝑡)} + ∫ 𝜑(𝜉)𝑑𝜉 + ∫ ∫ 𝑓(𝜉, 𝜏)𝑑𝜉𝑑𝜏
2 2𝑎 𝑥−𝑎𝑡 2𝑎 0 𝑥−𝑎(𝑡−𝜏)

Que no es más que la fórmula de D’Alembert donde tenemos que

𝑢𝑡𝑡 (𝑥, 𝑡) − 𝑢𝑥𝑥 = 𝑥𝑒 −𝑡−1
𝑢(𝑥, 0) = cos 2 𝑥
𝑢𝑡 (𝑥, 0) = 𝑒 −𝑥

Entonces tenemos que

𝜙(𝑥, 0) = cos2 𝑥
𝜑(𝑥, 0) = 𝑒 −𝑥
𝑓(𝑥, 𝑡) = 𝑥𝑒 −𝑡−1
𝑎=1

Ahora usando la fórmula

𝑥+𝑡 𝑡 𝑥+(𝑡−𝑠)
1 2 (𝑥 2 (𝑥 −𝑦
𝑢(𝑥, 𝑡) = (cos + 𝑡) + cos − 𝑡) + ∫ 𝑒 𝑑𝑦 + ∫ ∫ 𝑦𝑒 −𝑠−1 𝑑𝑦𝑑𝑠)
2 𝑥−𝑡 0 𝑥−(𝑡−𝑠)

1
𝑢(𝑥, 𝑡) = (cos 2 (𝑥 + 𝑡) + cos 2 (𝑥 − 𝑡) + 𝑒 𝑡−𝑥 − 𝑒 −𝑡−𝑥
2
𝑡 (𝑥
+ 𝑡 − 𝑠)2 − (𝑥 − 𝑡 + 𝑠)2 −𝑠−1
+ ∫ 𝑒 𝑑𝑠)
0 2

1
𝑢(𝑥, 𝑡) = (cos2 (𝑥 + 𝑡) + cos2 (𝑥 − 𝑡) + 𝑒 𝑡−𝑥 − 𝑒 −𝑡−𝑥 + 2𝑥(𝑒 −𝑡−1 + 𝑒 −1 (𝑡 − 1))
2

La solución única viene dada por

1 1
𝑢(𝑥, 𝑡) = (cos 2 (𝑥 + 𝑡) + cos2 (𝑥 − 𝑡)) + (𝑒 𝑡−𝑥 − 𝑒 −𝑡−𝑥 ) + 𝑥(𝑒 −𝑡−1 + 𝑒 −1 (𝑡 − 1))
2 2

2.- Una onda esférica es la solución de la ecuación de onda en tres dimensiones

1
𝑢𝑡𝑡 = 𝑐 2 ∇2 𝑢 de la forma 𝑢(𝑥, 𝑦, 𝑧, 𝑡) = 𝑈((𝑥 2 + 𝑦 2 + 𝑧 2 )2 , 𝑡). Encuentre la EDP para
𝑢(𝑟, 𝑡) y su solución general

Tenemos que usar coordenadas esféricas

𝜕2 𝜕2 𝜕2 1 𝜕 𝑟 2𝜕 1 𝜕 1 𝜕 𝜕2
+ + = ∙ ( )+ 2 (𝑠𝑒𝑛 𝜃 ) + 2
𝜕𝑥 2 𝜕𝑦 2 𝑑𝑧 2 𝑟 2 𝜕𝑟 𝜕𝑟 𝑟 𝑠𝑒𝑛𝜃𝜕𝜃 𝜕𝜃 𝑟 𝑠𝑒𝑛2 𝜃 𝜕𝜑 2

Ahora tomando en cuenta que la ecuación de onda se transforma mediante

𝑈(𝑟, 𝜃, 𝜑) = 𝑢(𝑟𝑠𝑒𝑛𝜃𝑐𝑜𝑠𝜑, 𝑟𝑠𝑒𝑛𝜃𝑠𝑒𝑛𝜑, 𝑟𝑐𝑜𝑠(𝜃, 𝑡))

Sabiendo que
 𝑥 = 𝑟𝑠𝑒𝑛𝜃𝑐𝑜𝑠𝜑
2 1
𝑈𝑡𝑡 = 𝑐 2 (𝑈𝑟𝑟 + 𝑈𝑟 ) = 𝑐 2 (𝑟𝑈)𝑟𝑟
𝑟 𝑟
Sabiendo que
𝑥 = 𝑟𝑠𝑒𝑛𝜃𝑐𝑜𝑠𝜑
𝑦 = 𝑟𝑠𝑒𝑛𝜃𝑠𝑒𝑛𝜑
𝑧 = 𝑟𝑐𝑜𝑠𝜃
1
(𝑥 2 + 𝑦 2 + 𝑧 2 )2 = √𝑥 2 + 𝑦 2 + 𝑧 2
1
(𝑥 2 + 𝑦 2 + 𝑧 2 )2 = √(𝑟𝑠𝑒𝑛𝜃𝑐𝑜𝑠𝜑)2 + (𝑟𝑠𝑒𝑛𝜃𝑠𝑒𝑛𝜑)2 + (𝑟𝑐𝑜𝑠𝜃)2
= √𝑟 2 𝑠𝑒𝑛2 𝜃 cos2 𝜑 + 𝑟 2 𝑠𝑒𝑛2 𝜃𝑠𝑒𝑛2 𝜑 + 𝑟 2 cos 2 𝜃
= √𝑟 2 𝑠𝑒𝑛2 𝜃(cos2 𝜑 + 𝑠𝑒𝑛2 𝜑) + cos 2 𝜃
= √𝑟 2 (𝑠𝑒𝑛2 𝜃 + cos 2 𝜃) = 𝑟

Entonces tenemos
1
𝑢(𝑥, 𝑦, 𝑧, 𝑡) = 𝑈 ((𝑥 2 + 𝑦 2 + 𝑧 2 )2 , 𝑡)
= 𝑢(𝑟𝑠𝑒𝑛𝜃𝑐𝑜𝑠𝜑, 𝑟𝑠𝑒𝑛𝜃𝑠𝑒𝑛𝜑, 𝑟𝑐𝑜𝑠(𝜃, 𝑡))
Entonces determinamos el laplaciano en coordenadas esféricas como:

1 𝜕 𝑟 2𝜕 1 𝜕 1 𝜕 𝜕2
∙ ( ) + (𝑠𝑒𝑛 𝜃 ) +
𝑟 2 𝜕𝑟 𝜕𝑟 𝑟 2 𝑠𝑒𝑛𝜃𝜕𝜃 𝜕𝜃 𝑟 2 𝑠𝑒𝑛2 𝜃 𝜕𝜑 2

Donde
𝑈𝑛 = 𝑐 2 Δ𝑢
2
2
1 𝜕 𝑟 𝜕 1 𝜕 1 𝜕 𝜕2
𝑈𝑛 = 𝑐 ( 2 ∙ ( )+ 2 (𝑠𝑒𝑛 𝜃 ) + 2 ) 𝑈(𝑟, 𝑡)
𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝑠𝑒𝑛𝜃𝜕𝜃 𝜕𝜃 𝑟 𝑠𝑒𝑛2 𝜃 𝜕𝜑 2

1 𝜕𝑈(𝑟, 𝑡) 𝑟 2 𝜕 1 𝜕𝑈(𝑟, 𝑡) 1 𝜕 𝜕 2 𝑈(𝑟, 𝑡)

𝑈𝑛 = 𝑐 2 ( 2 ∙ ( )+ 2 (𝑠𝑒𝑛 𝜃 )+ 2 )
𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝑠𝑒𝑛𝜃𝜕𝜃 𝜕𝜃 𝑟 𝑠𝑒𝑛2 𝜃 𝜕𝜑 2

1 𝜕𝑈(𝑟, 𝑡) 𝑟 2 𝜕
2
1 1 𝜕
𝑈𝑛 = 𝑐 ( 2 ∙ ( )+ 2 (𝑠𝑒𝑛 𝜃(0)) + 2 (0))
𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝑠𝑒𝑛𝜃𝜕𝜃 𝑟 𝑠𝑒𝑛2 𝜃

1 𝜕𝑈(𝑟, 𝑡) 𝑟 2 𝜕 1 𝜕 1
𝑈𝑛 = 𝑐 2 ( 2
∙ ( )) = 𝑐 2 ( 2 (𝑟 2 𝑈𝑟 )) = 𝑐 2 ( 2 (2𝑟𝑈𝑟 + 𝑟 2 𝑈𝑟𝑟 ))
𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝑟 𝑟

2 𝑐2 𝑐2
𝑈𝑛 = 𝑐 2 ( (𝑈𝑟 + 𝑈𝑟𝑟 )) = (2𝑈𝑟 + 𝑟𝑈𝑟𝑟 ) = (𝑈𝑟 + 𝑈𝑟 + 𝑟𝑈𝑟𝑟 )
𝑟 𝑟 𝑟
 𝑐2 𝜕 𝜕 𝑐2 𝜕 𝑐2 𝜕 𝜕
𝑈𝑛 = ( 𝑈 + 𝑟𝑈𝑟 ) = ( (𝑈 + 𝑟𝑈𝑟 )) = ( ( (𝑟𝑈)))
𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝑟 𝑟 𝜕𝑟 𝜕𝑟

𝑐2𝜕 2 1 2
𝑈𝑛 = (𝑟𝑈) 𝑟𝑟 = 𝑐 (𝑟𝑈)𝑟𝑟
𝑟𝜕𝑟 2 𝑟

Donde sabemos que la ecuación de onda se transforma en

1
𝑈𝑛 = 𝑐 2 (𝑟𝑈)𝑟𝑟
𝑟
Ahora a partir de aquí introducimos la función 𝑉(𝑟, 𝑡) = 𝑟𝑈(𝑟, 𝑡) para llegar a la
ecuación de onda unidimensional dado que

1 2
𝑈𝑡𝑡 =𝑐 (𝑟𝑈)𝑟𝑟
𝑟
Multiplicando en ambos lados por r de la ecuación
𝑟1
𝑟𝑈𝑡𝑡 = 𝑐 2 (𝑟𝑈)𝑟𝑟 𝑖𝑚𝑝𝑙𝑖𝑐𝑎 𝑟𝑈𝑡𝑡 = 𝑐 2 (𝑟𝑈)𝑟𝑟
𝑟

Dado que r no es función de t, entonces 𝑟𝑈𝑡𝑡 = (𝑟𝑈)𝑡𝑡 donde (𝑟𝑈)𝑡𝑡 = 𝑐 2 (𝑟𝑈)𝑟𝑟

Si sabemos que 𝑉(𝑟, 𝑡) = 𝑟𝑈(𝑟, 𝑡) entonces reducimos la ecuación a

𝑉𝑡𝑡 = 𝑐 2 𝑉𝑟𝑟

Que representa la ecuación de onda unidimensional con velocidad de

programación c, hagamos un cambio de escala y sea:
𝑟
𝑦=
𝑐

Entonces
𝜕2 𝜕 𝜕𝑉𝜕𝑦 𝜕 1 1𝜕𝑉𝑦 𝜕𝑦 1
(𝑉) = ( ) = ( 𝑉𝑦 ) = = 𝑉
𝜕𝑟 2 𝜕𝑟 𝜕𝑦𝜕𝑟 𝜕𝑟 𝑦 𝑐𝜕𝑦𝜕𝑟 𝑐 2 𝑦𝑦

Sustituyendo en
𝑉𝑡𝑡 = 𝑐 2 𝑉𝑟𝑟

Donde tenemos que 𝑉𝑡𝑡 = 𝑉𝑦𝑦 que es la ecuación unidimensional, a partir de aquí
vamos a la solución general 𝑉(𝑟, 𝑡) dado que

1
𝑈(𝑟, 𝑡) = 𝑉(𝑟, 𝑡)
𝑟
Donde se tiene que 𝑉𝑡𝑡 − 𝑉𝑦𝑦 = 0
Ahora tenemos que ver qué de tipo de ecuación es:

𝑎 = 1 2𝑏 = 0 𝑐 = −1
 𝑏 2 − 𝑎𝑐 = 02 − 1(−1) = 1

La ecuación es hiperbólica, ahora nos toca calcular el valor de 𝜆 donde

𝑎 = 1 𝑦 𝑐 = −1 𝑑𝑜𝑛𝑑𝑒

𝑏 − √𝑏 2 − 𝑎𝑐
𝜆=
𝑎
0 + √02 − 1(−1)
𝜆1 = =1
1

0 − √02 − 1(−1)
𝜆2 = = −1
1

𝑃𝑎𝑟𝑎 𝜆1
𝑑𝑦
1=
𝑑𝑡

∫ 𝑑𝑡 = ∫ 𝑑𝑦

𝑡 + 𝑐1 = 𝑦

Para 𝜆2
𝑑𝑦
−1 =
𝑑𝑡

∫ −𝑑𝑡 = ∫ 𝑑𝑦
−𝑡 + 𝑐2 = 𝑦

Sustituyendo los valores en las soluciones con el cambio de coordenadas

𝛼 = 𝑐1 𝛽 = 𝑐2 𝑡𝑒𝑛𝑒𝑚𝑜𝑠 𝑞𝑢𝑒 𝛼 = 𝑦 − 𝑡 𝛽 = 𝑦 + 𝑡

Calculamos los valores derivando con respecto a y y a t de 𝛼

𝛼 =𝑦−𝑡
𝛼𝑦 = 1 𝛼𝑡 = −1 𝛼𝑦𝑦 = 𝛼𝑡𝑡 = 𝛼𝑡𝑦 = 0

Calculamos los valores derivando con respecto a 𝛼 𝑦 𝑎 𝑡 𝑑𝑒 𝛽

𝛽 =𝑦+𝑡

𝛽𝑦 = 𝛽𝑡 = 1 𝛽𝑦𝑦 = 𝛽𝑡𝑡 = 𝛽𝑦𝑡 = 0

 𝑉𝑦𝑦 = 𝑉𝛼𝛼 𝛼𝑦2 + 2𝑉𝛼𝛽 𝛼𝑦 𝛽𝑦 + 𝑉𝛽𝛽 𝛽𝑦2 + 𝑉𝛼 𝛼𝑦𝑦 + 𝑉𝛽 𝛽𝑦𝑦

𝑉𝑡𝑡 = 𝑉𝛼𝛼 𝛼𝑡2 + 2𝑉𝛼𝛽 𝛼𝑡 𝛽𝑡 + 𝑉𝛽𝛽 𝛽𝑡2 + 𝑉𝛼 𝛼𝑡𝑡 + 𝑣𝑉𝛽 𝛽𝑡𝑡

Donde
𝑉𝑦𝑦 = 𝑉𝛼𝛼 (1)2 + 2𝑉𝛼𝛽 (1)(1) + 𝑉𝛽𝛽 (1)2 + 𝑉𝛼(0) + 𝑉𝛽 (0)
𝑉𝑦𝑦 = 𝑉𝛼𝛼 + 2𝑉𝛼𝛽 + 𝑉𝛽𝛽

𝑉𝑡𝑡 = 𝑉𝛼𝛼 (−1)2 + 2𝑉𝛼𝛽 (−1)(1) + 𝑉𝛽𝛽 (1)2 + 𝑉𝛼 (0) + 𝑉𝛽 (0)

𝑉𝑡𝑡 = 𝑉𝛼𝛼 − 2𝑉𝛼𝛽 + 𝑉𝛽𝛽
Sustituyendo en
𝑉𝑡𝑡 − 𝑉𝑦𝑦 = 0
(𝑉𝛼𝛼 − 2𝑉𝛼𝛽 + 𝑉𝛽𝛽 ) − (𝑉𝛼𝛼 − 2𝑉𝛼𝛽 + 𝑉𝛽𝛽 ) = 0
Y nos queda
−4𝑉𝛼𝛽 = 0

Y nos queda finalmente

𝑉𝛼𝛽 = 0

Integrando en la solución general

𝑉(𝛼, 𝛽) = 𝐹(𝛼) + 𝐺(𝛽)

Haciendo el cambio de variables tenemos que

𝛼 =𝑦−𝑡 𝛽 = 𝑦 + 𝑡 𝑑𝑜𝑛𝑑𝑒
𝑟 𝑟
𝑉(𝑦, 𝑡) = 𝐹(𝑦 − 𝑡) + 𝐺(𝑦 + 𝑡) = 𝐹 ( − 𝑡) + 𝐺 ( + 𝑡) = 𝐹(𝑟 − 𝑐𝑡) + 𝐺(𝑟 + 𝑐𝑡)
𝑐 𝑐

Como 𝑉(𝑟, 𝑡) = 𝑟𝑈(𝑟, 𝑡) entonces tenemos que

𝑟𝑈(𝑟, 𝑡) = 𝐹(𝑟 − 𝑐𝑡) + 𝐺(𝑟 + 𝑐𝑡)

Sustituyendo en
1
𝑈(𝑟, 𝑡) = 𝑉(𝑟, 𝑡)
𝑟
1
𝑈(𝑟, 𝑡) = [𝐹(𝑟 − 𝑐𝑡) + 𝐺(𝑟 + 𝑐𝑡)]
𝑟
Ahora bien, si tomamos la misma ecuación
 𝑉𝑡𝑡 − 𝑉𝑦𝑦 = 0
Resolviendo por el método de medidas esféricas debido a Poisson donde a cada
𝑔 ∈ 𝐶 0 (𝑅 𝑁 ) asociando a su medida G(x,y) sobre la superficie de la bola de centro
x y radio r

1
𝐺(𝑥; 𝑟) = ∫ 𝑔(𝑦)𝑑𝑠(𝑦)
𝑤𝑁 𝑟 𝑛−1 𝜕𝐵(𝑥;𝑟)

Mediante esta expresión podemos suponer que G está definida en r positivos

menores o iguales a cero, aplicando el teorema de derivación respecto de
parámetros en 𝑅 𝑁 𝑥𝑅 y deducir
𝑔 ∈ 𝐶 𝑘 (𝑅 𝑁 )(𝑅 𝑁 𝑥𝑅) para 𝑔 ∈ (𝑅 𝑁 )(𝑅{𝑁 𝑥 [0,∞) ) entonces k>2

Ahora multiplicando por 𝑟 1−𝑛 y derivando con respecto a r

𝜕𝐺 𝑟
lim (𝑥; 𝑟) = lim ( ∫ Δ𝑔(𝑥 + 𝑟𝜉)𝑑𝜉 ) = 0
𝑟→0 𝜕𝑟 𝑟→0 𝑤𝑁 𝐵(0;1)
En donde vemos la ecuación de Derbaux puesto que G es una función contínua
en𝑥 ∈ 𝑅 𝑁

lim 𝐺(𝑥; 𝑟) = 𝑔(𝑥)

𝑟→0

Probando que 𝑔 ∈ 𝑐 𝑘 (𝑅 𝑁 )(𝑅𝑁 𝑥[0, ∞)) ahora supongamos que u es la solución

clásica al problema de Cauchy para los datos iniciales de 𝑢𝑜 ∈ 𝐶 2 (𝑅 𝑁 ) 𝑦
𝑢1 ∈ 𝐶 1 (𝑅 𝑁 ) 𝑐𝑜𝑛 𝑁 ≥ 2 entonces definimos para cada (𝑥, 𝑡)𝑥[0, ∞) las cantidades
Y de ahí podemos decir que la solución clásica del problema de Cauchy es

Bibliografía

UNAD. Ecuaciones Diferenciales Parciales. Unidad 2, Preliminares. Fecha de

recuperación: 14/Abril/2022.
https://campus.unadmexico.mx/contenidos/DCEIT/BLOQUE2/MT/06/MEDP/U1/descargabl
es/MEDP_U1_contenido.pdf

Claudia Maercela Giordano. (2017). Ecuaciones diferenciales Parciales.

Buenos Aires Argentina: Editorial de la Universidad la Plata.

Gabriel López Garcia y Hugo Martinez Ortiz. (2013). Ecuaciones

Diferenciales Parciales. México D:F. Libro libre.
http://sgpwe.izt.uam.mx/files/users/uami/hect/Ecuaciones_en_Derivadas_Parciales/libro_G
abrile.pdf

Mauch, S. (nero 2004). Introduction to Methods of Applied Mathematics or Advanced

Mathematical Methods for Scientists and Engineers.