Medp U2 A3 Jucl

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 11

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
Actividad 3 Resolución de Problemas de Cauchy

Julio Cueva Lozano


Es 1821000315
Resuelva los siguientes problemas de Cauchy

1.- 𝑈𝑥𝑥 + 𝑢𝑥 + (2 − sin(𝑥) − cos(𝑥))𝑢𝑦 − (3 + cos2 (𝑥)) 𝑢𝑦𝑦 = 0


𝑥
𝑆𝑖 𝑢(𝑥, cos(𝑥)) = 0 𝑦 𝑢𝑦 (𝑥, cos(𝑥)) = 𝑒 −2 cos(𝑥)

Hay que ver qué tipo de ecuación diferencial parcial tenemos

Primero hay que obtener el discriminante

𝑎 = 1 2𝑏 = 0 𝑦 𝑐 = −(3 + cos 2 𝑥)
𝑏 − 4𝑎𝑐 = 02 − (−(3 + cos 2 𝑥) = 3 + cos2 𝑥
2

Y como el coseno va de -1 a 1, es obvio que el discriminante será positivo por el


+3, por lo tanto, será una ecuación hiperbólica

La ecuación tiene la forma, las letras a b c d e l y f depende de x y y, son variables

𝑎𝑢𝑥𝑥 + 2𝑏𝑢𝑥𝑦 + 𝑐𝑢𝑦𝑦 + 𝑑𝑢𝑥 + 𝑒𝑢𝑦 + 𝑙𝑢 = 𝑓


Para ello tenemos que
𝑎 = 1, 2𝑏 = 0, 𝑐 = −3 − cos2 𝑥 , 𝑑 = 1, 𝑒 = (2 − 𝑠𝑒𝑛 𝑥 − cos 𝑥), 𝑙 = 0𝑓 =0

La ecuación característica asociada es la siguiente


𝑎𝑑𝑦 2 − 2𝑏𝑑𝑥𝑑𝑦 + 𝑐𝑑𝑥 2 = 0

Dividiendo todo entre 𝑑𝑥 2 tenemos

𝑎𝑑𝑦 2 2𝑏𝑑𝑦
− − +𝑐 = 0
𝑑𝑥 2 𝑑𝑥

Haciendo una pequeña sustitución donde


𝑑𝑦
𝜆=
𝑑𝑥
𝑎𝜆2 − 2𝑏𝜆 + 𝑐 = 0
Por teoría de la bibliografía anexa, sabemos que al ser hiperbólicas significa que
𝜆1 ≠ 𝜆2
Y se puede factorizar
(𝑑𝑦 − 𝜆1 𝑑𝑥)(𝑑𝑦 − 𝜆2 𝑑𝑥) = 0
Donde las soluciones son curvas características de la forma
𝜙(𝑥, 𝑦) = 𝑐1 𝑦 𝜑 = 𝑐2

Mediante la siguiente transformación, podremos resolverla


𝜙(𝑥, 𝑦) = 𝜉 𝑦 𝜑(𝑥, 𝑦) = 𝜂
Donde
𝜕(𝜉, 𝜂)
≠0
𝜕(𝑥, 𝑦)

Haremos la transformación de 𝑢(𝑥, 𝑦) 𝑎 𝑈( 𝜉, 𝜂)

𝐴(𝜉, 𝜂)𝑈𝜉𝜉 + 2𝐵(𝜉, 𝜂)𝑈𝜉,𝜂 + 𝐶(𝜉. 𝜂)𝑈𝜂𝜂 + 𝐷(𝜉, 𝜂)𝑈𝜉 + 𝐸(𝜉, 𝜂)𝑈𝜂 + 𝐿(𝜉, 𝜂)𝑈 = 𝐹(𝜉, 𝜂)

Con
𝐴(𝜉, 𝜂) = 𝑎𝜙𝑥2 + 2𝑏𝜙𝑥 𝜙𝑦 + 𝑐𝜙𝑦2
𝐵(𝜉, 𝜂) = 𝑎𝜙𝑥 𝜑𝑥 + 𝑏(𝜙𝑥 𝜑𝑦 + 𝜙𝑦 𝜑𝑥 ) + 𝑐𝜙𝑦 𝜑𝑦
𝐶(𝜉, 𝜂) = 𝑎𝜑𝑥2 + 2𝑏𝜑𝑥 𝜑𝑦 + 𝑐𝜑𝑦2
𝐷(𝜉, 𝜂) = 𝑎𝜙𝑥𝑥 + 2𝑏𝜙𝑥𝑦 + 𝑐𝜙𝑦𝑦 + 𝑑𝜙𝑥 + 𝑒𝜙𝑦
𝐸(𝜉, 𝜂) = 𝑎𝜑𝑥𝑥 + 2𝑏𝜑𝑥𝑦 + 𝑐𝜑𝑦𝑦 + 𝑑𝜑𝑥 + 𝑒𝜑𝑦
𝐿(𝜉, 𝜂) = 𝑙
𝐹(𝜉, 𝜂) = 𝑓

Con la ecuación característica asociada

𝑎(𝑑𝜂)2 − 2𝑏(𝑑𝜉)(𝑑𝜂) + 𝑐(𝑑𝜂)2 = 0

Donde la relación de la original y la transformada son


2
2
𝜕(𝜉, 𝜂)
𝐵 − 𝐴𝐶 = (𝑏 2 − 𝑎𝑐) [ ]
𝜕(𝑥, 𝑦)

Entonces diferenciamos para encontrar las características

𝜙𝑥 𝑑𝑥 + 𝜙𝑦 𝑑𝑦 = 0
𝜑𝑥 𝑑𝑥 + 𝜑𝑦 𝑑𝑦 = 0
En otras palabras
𝜙𝑥 𝜑𝑥
𝑑𝑦 + 𝑑𝑥 = 0 𝑦 𝑑𝑦 + 𝑑𝑥 = 0
𝜙𝑦 𝜑𝑦
Donde tenemos los valores
𝜙𝑥 𝜑𝑥
𝜆1 = −𝑦 𝜆2 = −
𝜙𝑦 𝜑𝑦
Usando la siguiente función para encontrar las lambdas, tenemos entonces

𝑏 ± √𝑏 2 − 𝑎𝑐
𝜆= 𝑎 = 1 2𝑏 = 0 𝑐 = −3 − cos2 𝑥
𝑎

0 − √(02 − 1(−3 − cos 2 𝑥) 1


𝜆1 = = −(3 + cos 2 𝑥)2
1
Análogamente para
1
𝜆2 = (3 + cos 2 𝑥)2

Como
𝑑𝑦
𝜆=
𝑑𝑥
𝑑𝑦 1
= −(3 + cos 2 𝑥)2
𝑑𝑥
Pasando dx a la derecha
1
𝑑𝑦 = −(3 + cos 2 𝑥)2 𝑑𝑥
Integrando
1
∫ 𝑑𝑦 = ∫ − (3 + cos2 𝑥)2 𝑑𝑥
La integral
1
= ∫ − (3 + cos 2 𝑥)2 𝑑𝑥

No la puedo resolver de manera analítica porque no hay las condiciones


suficientes para unicidad y existencia de una solución.

2.- 𝑢𝑥𝑥 + 2 cos 𝑥 𝑢𝑥𝑦 − sin2 𝑥 𝑢𝑦𝑦 − 𝑠𝑖𝑛𝑥 𝑢𝑦 = 0


𝑆𝑖 𝑢(𝑥, sin 𝑥) = 𝑥 + cos 𝑥 𝑢𝑦 = (𝑥, sin 𝑥) = sin 𝑥

Como el anterior, primero debemos de saber qué tipo edp es

𝑎 = 1 3𝑏 = 2 cos 𝑥 𝑦 𝑐 = −𝑠𝑒𝑛2 𝑥

𝑏 2 − 𝑎𝑐 = cos2 𝑥 − 1(−𝑠𝑒𝑛2 𝑥) = cos 2 𝑥 + 𝑠𝑒𝑛2 𝑥 = 1


Y como siempre dará 1, entonces siempre será mayor que 0, por lo tanto, es una
ecuación diferencial parcial hiperbólica

Es la misma teoría del ejercicio anterior, asi que iré directamente primero a
resolver las lambas

𝑏 ± √𝑏 2 − 𝑎𝑐
𝜆= 𝑎 = 1 2𝑏 = 2 cos 𝑥 𝑐 = − sin2 𝑥
𝑎

− cos 𝑥 + √(cos2 𝑥 − (−𝑠𝑒𝑛2 𝑥))


𝜆1 = = cos 𝑥 + 1
1
De manera análoga para 𝜆2
− cos 𝑥 − √(cos 2 𝑥 − (−𝑠𝑒𝑛2 𝑥))
𝜆2 = = cos 𝑥 − 1
1
𝑑𝑦
Para 𝜆1 = 𝑑𝑥

𝑑𝑦
cos 𝑥 + 1 =
𝑑𝑥
Pasando dx a la izquierda, e integrando ambos lados tenemos

∫ (cos 𝑥 + 1)𝑑𝑥 = ∫ 𝑑𝑦

𝑠𝑒𝑛 𝑥 + 𝑥 + 𝑐1 = 𝑦

𝑐1 = 𝑦 − 𝑠𝑒𝑛 𝑥 − 𝑥

Para 𝜆2
𝑑𝑦
cos 𝑥 − 1 =
𝑑𝑥

∫ (cos 𝑥 − 1)𝑑𝑥 = ∫ 𝑑𝑦

𝑠𝑒𝑛 𝑥 − 𝑥 + 𝑐1 = 𝑦

𝑐2 = 𝑦 − 𝑠𝑒𝑛 𝑥 + 𝑥

Tenemos ahora que sustituir los valores en las soluciones con el cambio de
coordenadas
𝜉 = 𝑐1 𝑦 𝜂 = 𝑐2

Tenemos ahora que


𝜉 = 𝑦 − 𝑠𝑒𝑛 𝑥 − 𝑥 𝜂 = 𝑦 − 𝑠𝑒𝑛𝑥 + 𝑥

Calculamos los valores derivando con respecto a 𝜉 𝑦 𝜂 𝑑𝑒 𝜉


𝜙𝑥 = − cos 𝑥 − 1
𝜙𝑦 = 1
𝜙(𝑥, 𝑦) = 𝜉 = 𝑦 − 𝑠𝑒𝑛 𝑥 − 𝑥 𝜙𝑥𝑥 = 𝑠𝑒𝑛 𝑥
𝜙𝑦𝑦 = 0
{ 𝜙𝑥𝑦 = 0

Ahora calculamos los valores derivando con respecto a 𝜉 𝑦 𝜂 𝑑𝑒 𝜂


𝜑𝑥 = − cos 𝑥 + 1
𝜑𝑦 = 1
𝜑(𝑥, 𝑦) = 𝜂 = 𝑦 − 𝑠𝑒𝑛 𝑥 + 𝑥𝑥 𝜑𝑥𝑥 = 𝑠𝑒𝑛 𝑥
𝜑𝑦𝑦 = 0
{ 𝜑𝑥𝑦 = 0

Sustituimos los valore en


𝐴(𝜉, 𝜂) = 𝑎𝜙𝑥2 + 2𝑏𝜙𝑥 𝜙𝑦 + 𝑐𝜙𝑦2
𝐵(𝜉, 𝜂) = 𝑎𝜙𝑥 𝜑𝑥 + 𝑏(𝜙𝑥 𝜑𝑦 + 𝜙𝑦 𝜑𝑥 ) + 𝑐𝜙𝑦 𝜑𝑦
𝐶(𝜉, 𝜂) = 𝑎𝜑𝑥2 + 2𝑏𝜑𝑥 𝜑𝑦 + 𝑐𝜑𝑦2
𝐷(𝜉, 𝜂) = 𝑎𝜙𝑥𝑥 + 2𝑏𝜙𝑥𝑦 + 𝑐𝜙𝑦𝑦 + 𝑑𝜙𝑥 + 𝑒𝜙𝑦
𝐸(𝜉, 𝜂) = 𝑎𝜑𝑥𝑥 + 2𝑏𝜑𝑥𝑦 + 𝑐𝜑𝑦𝑦 + 𝑑𝜑𝑥 + 𝑒𝜑𝑦
𝐿(𝜉, 𝜂) = 𝑙
𝐹(𝜉, 𝜂) = 𝑓

Del original tenemos


𝑎 = 1 2𝑏 = 2 cos 𝑥 𝑐 = −𝑠𝑒𝑛2 𝑥 𝑑=0 𝑒 = −𝑠𝑒𝑛 𝑥 𝑙 = 0 𝑓 = 0

𝐴(𝜉, 𝜂) = (1)(− cos 𝑥 − 1)2 + 2(cos 𝑥)(− cos 𝑥 − 1) + (−𝑠𝑒𝑛2 𝑥)(12 )


= cos2 𝑥 + 2 cos 𝑥 + 1 − 2 cos 2 𝑥 − 2 cos 𝑥 − 𝑠𝑒𝑛2 𝑥 = − cos 2 𝑥 − 𝑠𝑒𝑛2 𝑥 + 1 = −1 +
1=0

𝐵(𝜉, 𝜂) = (1)(− cos 𝑥 − 1)(− cos 𝑥 + 1) + cos 𝑥(− cos 𝑥 − 1 + 1)(− cos 𝑥 + 1) −
𝑠𝑒𝑛2 𝑥(1)(1)
= (− cos 𝑥 − 1)(− cos 𝑥 + 1) + cos 𝑥 (− cos 𝑥 − 1) + (− cos 𝑥 + 1) − 𝑠𝑒𝑛2 𝑥
= cos2 𝑥 − 1 − 2 cos 2 𝑥 − 𝑠𝑒𝑛2 𝑥 = − cos2 𝑥 − 1 − 𝑠𝑒𝑛2 𝑥 = −1 − 1 = −2

𝐶(𝜉, 𝜂) = (1)(− cos 𝑥 + 1)2 + 2(cos 𝑥)(− cos 𝑥 + 1)(1) + (−𝑠𝑒𝑛2 𝑥)(12 )
= cos2 𝑥 − 2 cos 2 𝑥 − 2 cos 𝑥 + 2 cos 𝑥 − 𝑠𝑒𝑛2 𝑥 + 1 = − cos 2 𝑥 − 𝑠𝑒𝑛2 𝑥 + 1 = 0

𝐷(𝜉, 𝜂) = (1)(𝑠𝑒𝑛 𝑥) + 2 cos 𝑥 (0) + (−𝑠𝑒𝑛2 𝑥)(0) + (0)(− cos 𝑥 − 1) + (−𝑠𝑒𝑛 𝑥)(1)
= 𝑠𝑒𝑛 𝑥 − 𝑠𝑒𝑛 𝑥 = 0

𝐸(𝜉, 𝜂) = (1)(𝑠𝑒𝑛 𝑥) + 2 cos 𝑥 (0) + (−𝑠𝑒𝑛2 𝑥)(0) + (0)(− cos 𝑥 − 1) + (−𝑠𝑒𝑛 𝑥)(1)
= 𝑠𝑒𝑛 𝑥 − 𝑠𝑒𝑛 𝑥 = 0

Tenemos entonces ahora

𝐴(𝜉, 𝜂)𝑈𝜉𝜉 + 2𝐵(𝜉, 𝜂)𝑈𝜉,𝜂 + 𝐶(𝜉. 𝜂)𝑈𝜂𝜂 + 𝐷(𝜉, 𝜂)𝑈𝜉 + 𝐸(𝜉, 𝜂)𝑈𝜂 + 𝐿(𝜉, 𝜂)𝑈 = 𝐹(𝜉, 𝜂)

Sustituyendo:

0𝑈𝜉𝜉 + 2(−2)(𝜉, 𝜂)𝑈𝜉,𝜂 + 0𝑈𝜂𝜂 + 0𝑈𝜉 0𝑈𝜂 + 0𝑈 = 0

−4𝑈𝜉𝜂 = 0
𝑈𝜉𝜂 = 0
Para obtener la solución general

𝜉 = 𝑦 − 𝑠𝑒𝑛 𝑥 − 𝑥 𝜂 = 𝑦 − 𝑠𝑒𝑛𝑥 + 𝑥

𝑢(𝑥, 𝑦) = ℎ1 𝜙(𝑥, 𝑦) + ℎ2 𝜑(𝑥, 𝑦)


𝑢(𝑥, 𝑦) = ℎ1 (𝑦 − 𝑠𝑒𝑛 𝑥 − 𝑥) + ℎ2 (𝑦 − 𝑠𝑒𝑛 𝑥 + 𝑥)
Es hora de tomar las condiciones iniciales

𝑢(𝑥. sin 𝑥 ) = 𝑥 + cos 𝑥 𝑢𝑦 = (𝑥, sin 𝑥) = 𝑠𝑒𝑛 𝑥


Sustituyendo
𝑢(𝑠, 𝑠𝑒𝑛 𝑥) = ℎ1 (𝑠𝑒𝑛 𝑥 − 𝑠𝑒𝑛 𝑥 − 𝑥) + ℎ2 (𝑠𝑒𝑛 𝑥 − 𝑠𝑒𝑛 𝑥 + 𝑥)
𝑢(𝑥, 𝑠𝑒𝑛 𝑥) = ℎ1 (−𝑥) + ℎ2 (𝑥) = 𝑥 + cos 𝑥

Si derivamos tenemos
1 − 𝑠𝑒𝑛 𝑥 = ℎ1 (−𝑥) + ℎ2 (𝑥)
Sustituyendo
𝑢𝑦 = 𝜙𝑦 ℎ1 (𝑦 − 𝑠𝑒𝑛 𝑥 − 𝑥) + 𝜑𝑦 ℎ2 (𝑦 − 𝑠𝑒𝑛 𝑥 + 𝑥)
𝑢𝑦 = ℎ1 (1)(𝑦 − 𝑠𝑒𝑛 𝑥 − 𝑥) + ℎ2 (𝑦 − 𝑠𝑒𝑛 𝑥 + 𝑥)
𝑢𝑦 (𝑥, 𝑠𝑒𝑛 𝑥) = ℎ1 (−𝑥) + ℎ2 (𝑥)
𝑢𝑦 (𝑥, 0) = ℎ1 (−𝑥) + ℎ2 (𝑥) = 𝑠𝑒𝑛 𝑥
Integrando tenemos
𝑥
∫ ℎ1 (−𝑥) + ∫ ℎ2 (𝑥) = ∫ 𝑠𝑒𝑛 𝑥 𝑑𝑥
0
−ℎ1 (−𝑥) + ℎ2 (𝑥) = − cos 𝑥
Sumando
ℎ1 (−𝑥) + ℎ2 (𝑥) = 𝑥 + cos 𝑥
−ℎ1 (−𝑥) + ℎ2 (𝑥) = − cos 𝑥
Tenemos que
2ℎ2 (𝑥) = 𝑥
𝑥
ℎ2 (𝑥) =
2

Sustituyendo
ℎ1 (−𝑥) + ℎ2 (𝑥) = 𝑥 + cos 𝑥
𝑥
ℎ1 (−𝑥) + = 𝑥 + cos 𝑥
2

Para la fórmula general


𝑢(𝑥, 𝑦) = ℎ1 (𝑦 − 𝑠𝑒𝑛 𝑥 − 𝑥) + ℎ2 (𝑦 − 𝑠𝑒𝑛 𝑥 + 𝑥)
Tenemos entonces
1 1
𝑢(𝑥, 𝑦) = (−𝑦 + 𝑠𝑒𝑛 𝑥 + 𝑥) + cos(𝑦 − 𝑠𝑒𝑛 𝑥 − 𝑥) + (𝑦 − 𝑠𝑒𝑛𝑥 + 𝑥)
2 2
𝑦 𝑠𝑒𝑛 𝑥 𝑥 𝑦 𝑠𝑒𝑛 𝑥 𝑥
𝑢(𝑥, 𝑦) = − + + + cos(−𝑦 + 𝑠𝑒𝑛 𝑥 + 𝑥) + − +
2 2 2 2 𝑥 2

Y finalmente la solución particular es

𝒖(𝒙, 𝒚) = 𝒙 + 𝐜𝐨𝐬(−𝒚 + 𝒔𝒆𝒏 𝒙 + 𝒙)


2𝑦
3.- 4𝑦 2 𝑢𝑥𝑥 − 𝑢𝑦𝑦 + 2(−𝑦 2 + 1)𝑢𝑥𝑦 − 1+𝑦 2 (2𝑢𝑥 + 𝑢𝑦 ) = 0
𝑠𝑖(𝑥, 0) = 𝑓(𝑥) 𝑦 𝑢𝑦 (𝑥, 0) = 𝑔(𝑥)
¿Qué tipo de edp tenemos?

𝑎 = 4𝑦 2 2𝑏 = 2(−𝑦 2 + 1) 𝑐 = −1
𝑏 2 − 𝑎𝑐 = (−𝑦 2 + 1)2 − (4𝑦 2 )(−1) = 𝑦 4 − 2𝑦 2 + 1 + 4𝑦 2
= (𝑦 2 + 1)2
El +1 hará que la edp siempre será >0, entonces es hiperbólica

Es de la forma
𝑎𝑢𝑥𝑥 + 2𝑏𝑢𝑥𝑦 + 𝑐𝑢𝑦𝑦 + 𝑑𝑢𝑥 + 𝑒𝑢𝑦 + 𝑙𝑢 = 𝑓

2𝑦 2𝑦
𝑎 = 4𝑦 2 2𝑏 = 2(−𝑦 2 + 1) 𝑐 = −1 𝑑 = −2 ( ) 𝑒 = − ,𝑙 = 0 𝑓 = 0
1 + 𝑦2 1 + 𝑦2

De nuevo directamente con las lambdas

𝑏 ± √𝑏 2 − 𝑎𝑐
𝜆= 𝑎 = 4𝑦 2 2𝑏 = 2(−𝑦 2 + 1) 𝑐 = −1
𝑎

−𝑦 2 + 1 + √{((−𝑦 2 + 1)2 −} −𝑦 2 + 1 + √(𝑦 2 + 1)2


𝜆1 = =
4𝑦 2 4𝑦 2
−𝑦 2 + 1 + 𝑦 2 + 1 −𝑦 2 + 1 + 𝑦 2 + 1 2 1
= 2
= 2
= 2= 2
4𝑦 4𝑦 4𝑦 2𝑦

−𝑦 2 + 1 − √{((−𝑦 2 + 1)2 −} −𝑦 2 + 1 − √(𝑦 2 + 1)2


𝜆2 = =
4𝑦 2 4𝑦 2
−𝑦 + 1 − 𝑦 − 1 −𝑦 + 1 − 𝑦 − 1 −2𝑦 2
2 2 2 2
1
= 2
= 2
= 2
=−
4𝑦 4𝑦 4𝑦 2

𝑃𝑎𝑟𝑎 𝜆1
1 𝑑𝑦
=
2𝑦 2 𝑑𝑥
𝑑𝑥 = 2𝑦 2 𝑑𝑦
Integramos
∫ 𝑑𝑥 = ∫ 2𝑦 2 𝑑𝑦
2
𝑥 = 𝑦 3 + 𝑐1
3
2
𝑐1 = 𝑦 3 − 𝑥
3
Para 𝜆2
1 𝑑𝑦
− =
2 𝑑𝑥
∫ 𝑑𝑥 = ∫ − 2𝑑𝑦
𝑐2 = 2𝑦 + 𝑥

Calculamos los valores derivando con respecto a 𝜉 𝑦 𝜂 𝑑𝑒 𝜉


𝜙𝑥 = 1
𝜙𝑦 = −2𝑦 2
2 3
𝜙(𝑥, 𝑦) = 𝜉 = 𝑥 − 𝑦 𝜙𝑥𝑥 = 0
3 𝜙𝑦𝑦 = −4𝑦
{ 𝜙𝑥𝑦 = 0

Ahora calculamos los valores derivando con respecto a 𝜉 𝑦 𝜂 𝑑𝑒 𝜂


𝜑𝑥 = 1
𝜑𝑦 = 2
𝜑(𝑥, 𝑦) = 𝜂 = 2𝑦 + 𝑥 𝜑𝑥𝑥 = 0
𝜑𝑦𝑦 = 0
{𝜑𝑥𝑦 = 0

Hagamos las sustituciones para la transformación

𝐴(𝜉, 𝜂) = (4𝑦 2 )(12 ) + 2(−𝑦 2 + 1)(1)(−2𝑦 2 ) + (−1)(−2𝑦 2 )2


= 4𝑦 2 + 4𝑦 4 − 4𝑦 2 − 4𝑦 4 = 0

𝐵(𝜉, 𝜂) = (4𝑦 2 )(1)(1) + (−𝑦 2 + 1)(2 + (−2𝑦 2 ) + (−1)(−2𝑦 2 )(2)


= 4𝑦 2 + (−𝑦 2 + 1)(2 − 2𝑦 2 ) + 4𝑦 2 = 4𝑦 2 + 2 + 2𝑦 4 − 4𝑦 2 = 4𝑦 2 + 2 + 2𝑦 4

𝐶(𝜉, 𝜂) = (4𝑦 2 )(12 ) + 2(−𝑦 2 + 1)(1)(2) + (−1)(22 ) = 4𝑦 2 + 4(−𝑦 2 + 1) − 4 +


= 4𝑦 2 − 4𝑦 2 + 4 − 4 = 0

2𝑦 2𝑦
𝐷(𝜉, 𝜂) = (4𝑦 2 )(0) + 2(−𝑦 2 + 1)(0) + (−1)(−4𝑦) − 2 ( 2
) (1) − (−2𝑦 2 )
1+𝑦 1 + 𝑦2
4𝑦 4𝑦 3 4𝑦+4𝑦 3 1+𝑦 2
= 4𝑦 − 1+𝑦 2 + 1+𝑦 2 = 4𝑦 − = 4𝑦 − 4𝑦 (1+𝑦 2 ) = 4𝑦 − 4𝑦(1) = 0
1+𝑦 2

2𝑦 2𝑦
𝐸(𝜉, 𝜂) = (4𝑦 2 )(0) + 2(−𝑦 2 + 1)(0) + (−1)(0) − 2 (1+𝑦 2 ) (1) + 1+𝑦 2 (2) = 0

Tenemos entonces ahora

𝐴(𝜉, 𝜂)𝑈𝜉𝜉 + 2𝐵(𝜉, 𝜂)𝑈𝜉,𝜂 + 𝐶(𝜉. 𝜂)𝑈𝜂𝜂 + 𝐷(𝜉, 𝜂)𝑈𝜉 + 𝐸(𝜉, 𝜂)𝑈𝜂 + 𝐿(𝜉, 𝜂)𝑈 = 𝐹(𝜉, 𝜂)

Sustituyendo:

0𝑈𝜉𝜉 + 2(4𝑦 2 + 2 + 2𝑦 4 )𝑈𝜉,𝜂 + 0𝑈𝜂𝜂 + 0𝑈𝜉 + 0𝑈𝜂 + 0𝑈 = 0


2(4𝑦 2 + 2 + 2𝑦 4 )𝑈𝜉,𝜂 = 0
(8𝑦 2 + 4 + 4𝑦 4 )𝑈𝜉,𝜂 = 0
𝑈𝜉,𝜂 = 0

Para obtener la solución general

2
𝜉 = 𝑥 − 𝑦3 𝜂 = 2𝑦 + 𝑥
3

𝑢(𝑥, 𝑦) = ℎ1 𝜙(𝑥, 𝑦) + ℎ2 𝜑(𝑥, 𝑦)

2
𝑢(𝑥, 𝑦) = ℎ1 (𝑥 − 𝑦 3 ) + ℎ2 (2𝑦 + 𝑥)
3

(𝑥, 0) = 𝑓(𝑥) 𝑦 𝑢𝑦 (𝑥, 0) = 𝑔(𝑥)


Sustituyendo
2
(𝑥, 0) = ℎ1 (𝑥 − (0)3 ) + ℎ2 (2(0) + 𝑥)
3

𝑢(𝑥, 0) = ℎ1 (𝑥) + ℎ2 (𝑥) = 𝑓(𝑥)

ℎ1 (𝑥) + ℎ2 (𝑥) = 𝑓(𝑥)

Derivando tenemos que


𝑓´(𝑥) = ℎ1′ (𝑥) + ℎ2′ (𝑥)
Sustituyendo 𝑢𝑦 (𝑥, 0) = 𝑔(𝑥) 𝑡𝑒𝑛𝑒𝑚𝑜𝑠 𝑞𝑢𝑒 𝑥 = 𝑥 𝑦 = 0
2 2
𝑢𝑦 = 𝜙𝑦 ℎ1 (𝑥 − 𝑦 3 ) + 𝜑𝑦 ℎ2 (2𝑦 + 𝑥) = ℎ1 (−2𝑦 2 ) (𝑥 − 𝑦 3 ) + ℎ2 (2)(2𝑦 + 𝑥)
3 3

Teniendo que x=x y y=0


2
𝑢𝑦 (𝑥, 0) = ℎ1 (−2(0)2 ) (𝑥 − (0)3 ) + ℎ2 (2)(2(0) + 𝑥)
3
𝑢𝑦 (𝑥, 0) = ℎ1 (0) + ℎ2 (2𝑥)
𝑢𝑦 (𝑥, 0) = 2ℎ2 (𝑥) = 𝑔(𝑥)
Haciendo el correcto despeje
2ℎ2 (𝑥) = 𝑔(𝑥)
𝑔(𝑥)
ℎ2 (𝑥) =
2
𝑥
𝑔(𝑠)
ℎ2 (𝑥) = ∫ 𝑑𝑠
𝑜 2

𝑓(𝑥) = ℎ1 (𝑥) + ℎ2 (𝑥)

ℎ1 (𝑥) = 𝑓(𝑥) − ℎ2 (𝑥)


𝑥
𝑔(𝑠)
ℎ1 (𝑥) = 𝑓(𝑥) − ∫ 𝑑𝑠
𝑜 2

Y al sustituir
2
𝑢(𝑥, 𝑦) = ℎ1 (𝑥 − 𝑦 3 ) + ℎ2 (2𝑦 + 𝑥)
3

Tenemos que la solución particular 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.

También podría gustarte