5-11-2015

Ecuaciones diferenciales Parciales.

Unidad 2. El problema de Cauchy

para EDP de segundo orden
Actividad 2. Discusión de dos problemas

[Link]

Joel Alberto Montalvo Hernández

AL12523631
ECUACIONES DIFERENCIALES PARCIALES
[Link] DEL CARMEN LOZANO ARIZMENDI

Joel Alberto Montalvo Hernández 05 Noviembre 2015

Dra. María del Carmen Lozano Arizmendi
Ecuaciones diferenciales parciales
Unidad 1. El problema de Cauchy para EDP de segundo orden.
Actividad 2. Discusión de dos problemas.
1. I) Para una cuerda infinita (i.e. no hay que preocuparse por las condiciones iniciales), ¿qué
condiciones iniciales darían lugar a una onda que se mueve sólo hacia adelante? Exprese su
respuesta en términos del desplazamiento inicial 𝒖 ( 𝒙, 𝟎 ) = 𝒇 (𝒙) y velocidad inicial 𝒖𝒕 (𝒙, 𝟎 ) = 𝒈( 𝒙)
y sus derivadas 𝒇 ′(𝒙 ), 𝒈 ‘( 𝒙 ) . Interpretar el resultado intuitivamente.
La solución de la ecuación de la cuerda nos la proporciona la fórmula de D’Alambert
𝑓(𝑥 + 𝑐𝑡) + 𝑓(𝑥 − 𝑐𝑡) 1 𝑥+𝑐𝑡
𝑢(𝑥, 𝑡) = 𝐹(𝑥 − 𝑐𝑡) + 𝐺(𝑥 + 𝑐𝑡) = + ∫ 𝑔(𝑠) 𝑑𝑠,
2 2𝑐 𝑥−𝑐𝑡
𝑐𝑜𝑛 𝑐 𝑙𝑎 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑑𝑒 𝑝𝑟𝑜𝑝𝑜𝑝𝑎𝑔𝑎𝑐𝑖ó𝑛
𝐹(𝑥 − 𝑐𝑡)𝑟𝑒𝑝𝑟𝑒𝑠𝑒𝑛𝑡𝑎 𝑢𝑛𝑎 𝑜𝑛𝑑𝑎 𝑞𝑢𝑒 𝑠𝑒 𝑚𝑢𝑒𝑣𝑒 𝑎 𝑙𝑎 𝑑𝑒𝑟𝑒𝑐ℎ𝑎 𝑦 𝐺(𝑥 + 𝑐𝑡) 𝑢𝑛𝑎 𝑞𝑢𝑒 𝑠𝑒 𝑚𝑢𝑒𝑣𝑒 𝑎 𝑙𝑎 𝑖𝑧𝑞𝑢𝑖𝑒𝑟𝑑𝑎
𝑃𝑎𝑟𝑎 𝑜𝑏𝑡𝑒𝑛𝑒𝑟 𝑢𝑛𝑎 𝑜𝑛𝑑𝑎 𝑞𝑢𝑒 𝑠𝑒 𝑚𝑢𝑒𝑣𝑒 𝑠𝑜𝑙𝑜 𝑎 𝑙𝑎 𝑑𝑒𝑟𝑒𝑐ℎ𝑎 𝐺(𝑥 + 𝑐𝑡) = 𝐺 ′ (𝑥 + 𝑐𝑡) = 0
0, 𝑥 < 0
𝑙𝑎𝑠 𝑐𝑜𝑛𝑑𝑖𝑐𝑖𝑜𝑛𝑒𝑠 𝑖𝑛𝑖𝑐𝑖𝑎𝑙𝑒𝑠 𝑠𝑒𝑟á𝑛 𝑢(𝑥, 0) = { 𝑛𝑜 𝑡𝑒𝑛𝑒𝑚𝑜𝑠 𝑢𝑛 𝑝𝑒𝑟𝑓𝑖𝑙 𝑞𝑢𝑒 𝑠𝑒 𝑝𝑟𝑜𝑝𝑎𝑔𝑢𝑒 𝑎 𝑙𝑎 𝑑𝑒𝑟𝑒𝑐ℎ𝑎
𝑓(𝑥), 𝑥 > 0
0, 𝑥 < 0
𝑢𝑡 (𝑥, 0) = 𝑔(𝑥, 0) = { 𝑛𝑜 𝑡𝑒𝑛𝑒𝑚𝑜𝑠 𝑢𝑛 𝑖𝑚𝑝𝑢𝑙𝑠𝑜 ℎ𝑎𝑐𝑖𝑎 𝑙𝑎 𝑑𝑒𝑟𝑒𝑐ℎ𝑎
𝑔(𝑥), 𝑥 > 0
Físicamente tendríamos una cuerda infinita, sujeta por el punto medio (x=0) a la cual se la da un perfil y un
impulso inicial para un x>0, con esto se genera una onda que se moverá siempre hacia adelante (no hay impulso
que se proporcione hacia la izquierda) como la cuerda es infinita no tendremos una onda reflejada.
El equivalente sería una cuerda semi-infinita sujeta por un extremo.

ii) Una vez más para una cuerda infinita, supongamos que 𝒖 ( 𝒙, 𝟎 ) = 𝒇 (𝒙) y 𝒖𝒕 (𝒙, 𝟎 ) = 𝒈( 𝒙 ) son
cero para |𝒙 | > 𝒂, por algún número real 𝒂 > 𝟎. Probar que si 𝒕 + 𝒙 > 𝒂 y 𝒕 − 𝒙 > 𝒂, entonces el
desplazamiento 𝒖 ( 𝒙 , 𝒕) de la cadena es constante. Relacionar esta constante de 𝒈 ( 𝒙 ) .

Como u(x,0) f(x)=0 fuera del intervalo [-a,a] tenemos para t=0 una onda cuyos extremos son los puntos son
(-a,0) y (a,0), no importa el perfil de la onda puede ser triangular, cuadrada, semicircular, senoidal, etc. lo que
importa es que su longitud de onda es 2a, como g(x) que es la velocidad inicial, el impulso inicial de la onda
(que no es lo mismo que la velocidad de propagación c), solo se presenta en el intervalo (-a,a).
0, 𝑥 + 𝑡 > 𝑎
𝑆𝑒 𝑡𝑖𝑒𝑛𝑒 𝑔(𝑥) = {𝑔(𝑥), −𝑎 ≤ 𝑥 ≤ 𝑎
0, 𝑥 − 𝑡 < −𝑎
De acuerdo a la ecuación de D’Alambert la solución de u(x,t) es
𝑓(𝑥 + 𝑐𝑡) + 𝑓(𝑥 − 𝑐𝑡) 1 𝑥+𝑐𝑡
𝑢(𝑥, 𝑡) = 𝐹(𝑥 − 𝑐𝑡) + 𝐺(𝑥 + 𝑐𝑡) = + ∫ 𝑔(𝑠) 𝑑𝑠 =
2 2𝑐 𝑥−𝑐𝑡
−𝑎 𝑎 𝑥+𝑐𝑡
𝑓(𝑥 + 𝑐𝑡) + 𝑓(𝑥 − 𝑐𝑡) 1
= + [∫ 𝑔(𝑠) 𝑑𝑠 + ∫ 𝑔(𝑠)𝑑𝑠 + ∫ 𝑔(𝑠)𝑑𝑠] =
2 2𝑐 𝑥−𝑐𝑡 −𝑎 𝑎
−𝑎 𝑎 𝑥+𝑐𝑡
𝑓(𝑥 + 𝑐𝑡) + 𝑓(𝑥 − 𝑐𝑡) 1
= + [∫ (0) 𝑑𝑠 + ∫ 𝑔(𝑠)𝑑𝑠 + ∫ (0)𝑑𝑠] =
2 2𝑐 𝑥−𝑐𝑡 −𝑎 𝑎

Joel Alberto Montalvo Hernández 05 Noviembre 2015

Dra. María del Carmen Lozano Arizmendi.
Ecuaciones diferenciales parciales
Unidad 1. El problema de Cauchy para EDP de segundo orden.
Actividad 2. Discusión de dos problemas.
𝑎
𝑓(𝑥 + 𝑐𝑡) + 𝑓(𝑥 − 𝑐𝑡) 1 𝑓(𝑥 + 𝑐𝑡) + 𝑓(𝑥 − 𝑐𝑡) 1
= + [∫ 𝑔(𝑠)𝑑𝑠] = + (𝑘), 𝑘 𝑒𝑠 𝑒𝑙 𝑣𝑎𝑙𝑜𝑟 𝑑𝑒 𝑙𝑎 𝑖𝑛𝑡𝑒𝑔𝑟𝑎𝑙
2 2𝑐 −𝑎 2 2𝑐
𝑓(𝑥 + 𝑐𝑡) + 𝑓(𝑥 − 𝑐𝑡) 1
𝑠𝑒 𝑡𝑖𝑒𝑛𝑒 𝑢(𝑥, 𝑡) = + 𝑙𝑜 𝑐𝑢𝑎𝑙 𝑐𝑜𝑟𝑟𝑒𝑠𝑝𝑜𝑛𝑑𝑒 𝑎 𝑢𝑛𝑎 𝑜𝑛𝑑𝑎 𝑞𝑢𝑒 𝑠𝑒 𝑑𝑒𝑠𝑝𝑙𝑎𝑧𝑎 𝑐𝑜𝑛 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑
2 2𝑐
𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 𝑐

2. Considere el problema de Cauchy

𝒖𝒕𝒕 − 𝒖𝒙𝒙 = 𝟎, − ∞ < 𝒙 < ∞, 𝒕>𝟎

𝟎, − ∞ < 𝒙 < −𝟏
𝒙 + 𝟏, − 𝟏 ≤ 𝒙 ≤ 𝟎
𝒖(𝒙, 𝟎) = 𝒇(𝒙) = {
𝟏 − 𝒙, 𝟎≤𝒙≤𝟏
𝟎, 𝟏<𝒙<∞

𝟎, − ∞ < 𝒙 < −𝟏
𝒖𝒕 (𝒙, 𝟎) = 𝒈(𝒙) = { 𝟏, −𝟏 ≤ 𝒙 ≤ 𝟎
𝟎, 𝟏<𝒙<∞

𝐹𝑎𝑙𝑡𝑎 𝑒𝑙 𝑣𝑎𝑙𝑜𝑟 𝑑𝑒 𝑔(𝑥) 𝑝𝑎𝑟𝑎 0 < 𝑥 ≤ 1, 𝑝𝑜𝑟 𝑠𝑖𝑚𝑒𝑡𝑟í𝑎 𝑦 𝑐𝑜𝑛𝑠𝑖𝑑𝑒𝑟𝑎𝑛𝑑𝑜 𝑢(𝑥, 0), 𝑢𝑡 (𝑥, 0) = 𝑔(𝑥) = −1, 𝑡𝑒𝑛𝑒𝑚𝑜𝑠

0, − ∞ < 𝑥 < −1 0, − ∞ < 𝑥 < −1

1, − 1 ≤ 𝑥 ≤ 0 1, − 1 ≤ 𝑥 ≤ 0
𝑢𝑡 (𝑥, 0) = 𝑔(𝑥) = { 𝑜 𝑏𝑖𝑒𝑛 𝑢𝑡 (𝑥, 0) = 𝑔(𝑥) = {
−1, 0 < 𝑥 ≤ 1 1, 0 < 𝑥 ≤ 1
0, 1<𝑥<∞ 0, 1<𝑥<∞

Graficamos con Geogebra

Joel Alberto Montalvo Hernández 05 Noviembre 2015

Dra. María del Carmen Lozano Arizmendi.
Ecuaciones diferenciales parciales
Unidad 1. El problema de Cauchy para EDP de segundo orden.
Actividad 2. Discusión de dos problemas.

𝟏
a) Evalúa 𝒖 en el punto (𝟏, 𝟐). Sugerencia utilice la fórmula de D’Alembert.

𝑓(𝑥 + 𝑐𝑡) + 𝑓(𝑥 − 𝑐𝑡) 1 𝑥+𝑐𝑡

𝑢(𝑥, 𝑡) = + ∫ 𝑔(𝑠) 𝑑𝑠 , 𝑐 2 = 1 𝑐 = 1
2 2𝑐 𝑥−𝑐𝑡
1 1
𝑓 (1 + (1) (2)) + 𝑓 (1 − (1) (2)) 1
1+(1)( )
1 1 1 2
𝑒𝑛 (1, ) 𝑢 (1, ) = + ∫ 𝑔(𝑠) 𝑑𝑠 =
2 2 2 2(1) 1−(1)( )
1
2
3 1 3 3 1 3
𝑓 (2) + 𝑓 (2) 1 2 𝑓 (2) + 𝑓 (2) 1 1 1 2
= + ∫ 𝑔(𝑠) 𝑑𝑠 = + ∫ 𝑔(𝑠) 𝑑𝑠 + ∫ 𝑔(𝑠) 𝑑𝑠 =
2 2 1 2 2 1 2 1
2 2
1 3 1 1
0 + (1 − 2) 1 1 1 2 1 1 1 1 1 1 1
= + ∫ −1 𝑑𝑠 + ∫ 0 𝑑𝑠 = 2 + |−𝑥|1 + 0 = − (1 − ) = − ( ) =
2 2 1 2 1 2 2 4 2 2 4 2 2
2 2
1 1
= − =0
4 4

b) ¿Qué puedes decir acerca de las singularidades de la solución?

Joel Alberto Montalvo Hernández 05 Noviembre 2015

Dra. María del Carmen Lozano Arizmendi.
Ecuaciones diferenciales parciales
Unidad 1. El problema de Cauchy para EDP de segundo orden.
Actividad 2. Discusión de dos problemas.
0, − ∞ < 𝑥 < −1
1, − 1 ≤ 𝑥 ≤ 0
Consideremos 𝑔(𝑥) = {
1, 0 < 𝑥 ≤ 1
0, 1<𝑥<∞

𝑓(𝑥 + 𝑡) + 𝑓(𝑥 − 𝑡) 1 𝑥+𝑡

𝐷𝑎𝑑𝑎 𝑢(𝑥, 𝑡) = + ∫ 𝑔(𝑠) 𝑑𝑠 , 𝑐 2 = 1 𝑐 = 1
2 2𝑐 𝑥−𝑡
𝑖) 𝐹𝑖𝑗𝑎𝑚𝑜𝑠 𝜀 ∈ 𝑅 𝑦 𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑚𝑜𝑠 lim 𝑓( 𝜀 + 𝑡) = 0, lim 𝑓( 𝜀 − 𝑡) = 0,
𝑡→∞ 𝑡→∞
𝑥+𝑡 1
1 𝟏
lim ∫ 𝑔(𝑠) 𝑑𝑠 = ∫ 1 𝑑𝑠 = 𝑥 | = 1 − (−1) = 1 + 1 = 2 ∴ 𝐥𝐢𝐦 𝒖( 𝒙, 𝒕) = 𝟎 + 𝟎 + (𝟐)( ) = 𝟏
𝑡→∞ 𝑥−𝑡 −1 −1 𝒕→∞ 𝟐

𝒑𝒐𝒓 𝒍𝒐 𝒕𝒂𝒏𝒕𝒐 𝐥𝐢𝐦 𝒖( 𝒙, 𝒕) = 𝟎 + 𝟎 + 𝟏 = 𝟏

𝒕→∞

Calculemos el max(u(x,t))
max[𝑓(𝑥 + 𝑡)] = 1 𝑝𝑎𝑟𝑎 − 1 ≤ 𝑥 + 𝑡 ≤ 1, max[𝑔(𝑥 − 𝑡)] = 1 𝑝𝑎𝑟𝑎 − 1 ≤ 𝑥 − 𝑡 ≤ 1
𝑥+𝑡 1 1
max [∫ 𝑔(𝑠)𝑑𝑠] = ∫ 𝑔(𝑠) 𝑑𝑠 = ∫ 1 𝑑𝑠 = 1 − (−1) = 2
𝑥−𝑡 −1 −1

𝑒𝑙 𝑐𝑢𝑎𝑙 𝑠𝑒 𝑎𝑙𝑐𝑎𝑛𝑧𝑎 𝑒𝑛 𝑙𝑎 𝑖𝑛𝑡𝑒𝑟𝑠𝑒𝑐𝑐𝑖ó𝑛 𝑑𝑒 𝑙𝑜𝑠 𝑠𝑒𝑚𝑖𝑝𝑙𝑎𝑛𝑜𝑠

𝑥+𝑡 =1
𝑥 + 𝑡 ≥ 1 𝑦 𝑥 − 𝑡 ≤ 1 𝑠𝑖 𝑖𝑛𝑡𝑒𝑟𝑠𝑒𝑐𝑡𝑎𝑚𝑜𝑠 𝑎𝑚𝑏𝑜𝑠 𝑐𝑜𝑛𝑗𝑢𝑛𝑡𝑜𝑠 𝑠𝑒 𝑑𝑒𝑏𝑒𝑛 𝑠𝑎𝑡𝑖𝑓𝑎𝑐𝑒𝑟 𝑠𝑖𝑚𝑢𝑙𝑡á𝑛𝑒𝑎𝑚𝑒𝑛𝑡𝑒 →
𝑥 − 𝑡 = −1
→ 𝑥 = 0, 𝑡 = 1
𝑓(0 + 1) + 𝑓(0 − 1) 1 1 1+1 1 1 1
𝑒𝑣𝑎𝑙𝑢𝑎𝑛𝑑𝑜 𝑢(0,1) = + ∫ 𝑔(𝑠) 𝑑𝑠 = + ∫ 1 𝑑𝑠 = 1 + (2) = 2
2 2 −1 2 2 −1 2
max[𝑢(𝑥, 𝑡)] = 𝑢(0,1) = 2
Consideremos el otro caso para g(x)
0, − ∞ < 𝑥 < −1
1, − 1 ≤ 𝑥 ≤ 0
𝑔(𝑥) = {
−1, 0 < 𝑥 ≤ 1
0, 1<𝑥<∞
𝑖) 𝐹𝑖𝑗𝑎𝑚𝑜𝑠 𝜀 ∈ 𝑅 𝑦 𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑚𝑜𝑠 lim 𝑓( 𝜀 + 𝑡) = 0, lim 𝑓( 𝜀 − 𝑡) = 0,
𝑡→∞ 𝑡→∞
𝑥+𝑡 0 1
0 1
lim ∫ 𝑔(𝑠) 𝑑𝑠 = ∫ −1 𝑑𝑠 + ∫ 1 𝑑𝑠 = −𝑥 | + 𝑥 | = −(0 − (−1)) + (1 − 0) = −1 + 1 = 0 𝑙𝑜 𝑐𝑢𝑎𝑙 𝑐𝑜𝑖𝑛𝑐𝑖𝑑𝑒
𝑡→∞ 𝑥−𝑡 −1 0 −1 0
𝑐𝑜𝑛 𝑒𝑙 ℎ𝑒𝑐ℎ𝑜 𝑑𝑒 𝑞𝑢𝑒 𝑔(𝑥)𝑒𝑠 𝑖𝑚𝑝𝑎𝑟 𝑒𝑛 𝑒𝑙 𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑜 [−1,1],
𝒑𝒐𝒓 𝒍𝒐 𝒕𝒂𝒏𝒕𝒐 𝐥𝐢𝐦 𝒖( 𝒙, 𝒕) = 𝟎
𝒕→∞

Calculemos el max(u(x,t))
max[𝑓(𝑥 + 𝑡)] = 1 𝑝𝑎𝑟𝑎 − 1 ≤ 𝑥 + 𝑡 ≤ 1, max[𝑔(𝑥 − 𝑡)] = 1 𝑝𝑎𝑟𝑎 − 1 ≤ 𝑥 − 𝑡 ≤ 1

Joel Alberto Montalvo Hernández 05 Noviembre 2015

Dra. María del Carmen Lozano Arizmendi.
Ecuaciones diferenciales parciales
Unidad 1. El problema de Cauchy para EDP de segundo orden.
Actividad 2. Discusión de dos problemas.
𝑥+𝑡 0 1
max [∫ 𝑔(𝑠)𝑑𝑠] = max [ ∫ 𝑔(𝑠) 𝑑𝑠 + ∫ 𝑔(𝑠)] = 0
𝑥−𝑡 −1 0

𝑒𝑙 𝑐𝑢𝑎𝑙 𝑠𝑒 𝑎𝑙𝑐𝑎𝑛𝑧𝑎 𝑒𝑛 𝑙𝑎 𝑖𝑛𝑡𝑒𝑟𝑠𝑒𝑐𝑐𝑖ó𝑛 𝑑𝑒 𝑙𝑜𝑠 𝑠𝑒𝑚𝑖𝑝𝑙𝑎𝑛𝑜𝑠

𝑥+𝑡 =1
𝑥 + 𝑡 ≥ 1 𝑦 𝑥 − 𝑡 ≤ 1 𝑠𝑖 𝑖𝑛𝑡𝑒𝑟𝑠𝑒𝑐𝑡𝑎𝑚𝑜𝑠 𝑎𝑚𝑏𝑜𝑠 𝑐𝑜𝑛𝑗𝑢𝑛𝑡𝑜𝑠 𝑠𝑒 𝑑𝑒𝑏𝑒𝑛 𝑠𝑎𝑡𝑖𝑓𝑎𝑐𝑒𝑟 𝑠𝑖𝑚𝑢𝑙𝑡á𝑛𝑒𝑎𝑚𝑒𝑛𝑡𝑒 →
𝑥 − 𝑡 = −1
0 1
𝑓(0 + 1) + 𝑓(0 − 1) 1 1+1
→ 𝑥 = 0, 𝑡 = 1 𝑒𝑣𝑎𝑙𝑢𝑎𝑛𝑑𝑜 𝑢(0,1) = + (∫ 𝑔(𝑠) 𝑑𝑠 + ∫ 𝑔(𝑠) 𝑑𝑠) = +0=1
2 2 −1 0 2
𝑴𝒂𝒙[𝒖(𝒙, 𝒕)] = 𝒖(𝟎, 𝟏) = 𝟏

c) Discute la continuidad de la solución 𝒖

La condición inicial u(x,0)=f(x) es continua, la condición inicial ut (x,0)=g(x,0) es continua excepto en los
puntos x=±1, por lo cual la solución u(x,t) es continua pero es C1 a trozos.

Referencias
 [Link]
df.
 Tijonov, A., Samarsky, A., Ecuaciones de la física matemática, 2ª ed., 1980, Moscú, Ed. Mir.
 Petrovski, I., Lecciones sobre ecuaciones en derivadas parciales, 1969, Cuba, Instituto del
[Link] 75.

Joel Alberto Montalvo Hernández 05 Noviembre 2015

Dra. María del Carmen Lozano Arizmendi.