Ejercicio 1: Taller 1. Ecuaciones de Primer Y Segundo Orden
Ejercicio 1: Taller 1. Ecuaciones de Primer Y Segundo Orden
de Física y Electrónica.
TALLER 1.
ECUACIONES DE PRIMER Y SEGUNDO ORDEN.
Luis Favian Cordero Cantero1 , María Camila Martinez Betin2 , and Yeimi Paola Orozco3
1
Programa de Física, Departamento de Física y Electrónica, Universidad de Córdoba, Cra. 6 #77-305,
Montería - Córdoba, Colombia, Teléfono: +57 304 378 6036,
Septiembre 24 de 2024.
1 Ejercicio 1
uxx + 2uxy + ux = 0
Solución
Sean los valores de a, b, c:
a = 1, b = 1, c = 0
Para encontrar el cambio de variable adecuado y poder transformar a la forma canónica, primero debemos
ver qué tipo de ecuación tenemos. Para eso evaluamos el determinante:
1 1
det(M ) = = 0 − 1 = −1 < 0
1 0
Como det(M ) < 0, la EDP es de tipo hiperbólica.
Para llevarla a la forma canónica, las curvas características satisfacen la EDO:
2
dy dy
a − 2b +c=0
dx dx
Remplazando los valores de los coeficientes:
2
dy dy
−2 =0
dx dx
Para resolver la ecuación cuadrática asociada, fíjese que podemos factorizar:
dy dy
−2 =0
dx dx
Nos queda entonces que:
dy dy
=0 o =2
dx dx
De esta forma:
1
Universidad de Córdoba EDP
dy = 0 · dx o dy = 2 · dx
Integrando, obtenemos los dos cambios de variable:
y=C y y = 2x + C
Por lo tanto:
η=y y ξ = y − 2x
Ahora, procedemos a calcular uxx , uxy , ux mediante la regla de la cadena:
ux = uξ · ξx = uξ (−2) = −2uξ
uxx = 4uξξ
−4uξη − 2uξ = 0
La forma canónica queda entonces:
1
uξη + uξ = 0
2
Resolviendo la EDP:
1
(uξ )η + uξ = 0
2
Sea entonces:
w = uξ
Nos queda una EDO:
1
wη + w = 0
2
Resolviendo por separación de variables:
∂w 1
=− w
∂η 2
∂w 1
= − ∂η
w 2
1
ln |w| = − η + f (w)
2
Física 2 TALLER 1.
Universidad de Córdoba EDP
1
w = e− 2 η+f (w)
Recuperando variables:
1
uξ = f (ξ)e− 2 η
Para encontrar u, integramos respecto a ξ:
Z
1
u(ξ, η) = f (ξ)e− 2 η dξ
1
u(ξ, η) = F (ξ)e− 2 η + f (η)
Recordando el primer cambio de variable ξ = y − 2x y η = y, pasamos a expresar u en términos de x e y:
y
u(x, y) = F (y − 2x)e− 2 + f (y)
Esta es la solución de nuestra EDP.
2 Ejercicio 2
1 −2
det(M ) = = (1)(4) − (−2)(−2) = 0
−2 4
Como det(M ) = 0, la EDP es de tipo parabólica.
Para llevar a la forma canónica, las curvas características satisfacen la EDO:
dy b dy −2
= ⇒ =
dx a dx 1
dy = −2 dx ⇒ y = −2x + C
De aquí sacamos nuestro primer cambio de variable, el cual es:
ξ = y + 2x
η=x
Física 3 TALLER 1.
Universidad de Córdoba EDP
3 Ejercicio 3
A = x, B = 2x2 , C=0
Calculamos el discriminante:
Física 4 TALLER 1.
Universidad de Córdoba EDP
Como ∆ > 0, podemos decir que la EDP es hiperbólica, la cual satisface la siguiente EDO:
√ √
dy B± ∆ 2x2 ± 4x4 2x2 ± 2x2
= = =
dx 2A 2x 2x
Nos queda:
dy 4x2 dy
= = 2x ó =0
dx 2x dx
Por lo tanto, integrando:
dy = 2x dx ⇒ y = x2 + C ó dy = 0 dx ⇒ y=C
ξ = y − x2 y η=y
Física 5 TALLER 1.
Universidad de Córdoba EDP
1
u(ξ, η) = √ + F (ξ) + g(η)
2 ξ−η
Recuperando las variables originales x e y, con ξ = y − x2 y η = y, nos queda la solución de la EDP:
1
u(x, y) = p + F (y − x2 ) + g(y)
2 y − x2
4 Ejercicio 4
A = x2 , B = −2xy, C = y2
Calculamos el discriminante:
Como ∆ = 0, podemos decir que la EDP es parabólica, la cual satisface la siguiente EDO:
dy B −2xy y
= = 2
=−
dx 2A 2x x
Resolviendo la EDO:
dy dx
=−
y x
Z Z
dy dx
=−
y x
ln |y| = − ln |x| + C
ln |y| + ln |x| = C
ln |xy| = C
xy = eC = C1
Por lo tanto, nuestra primera sustitución será:
ξ = xy
Dado que estamos en el caso parabólico, podemos tomar una segunda variable η = x. Aplicamos la regla de
la cadena.
Para ux :
∂ξ ∂η
ux = uξ · + uη · = uξ · y + uη · 1
∂x ∂x
ux = yuξ + uη
Para uxx :
∂ ∂
uxx = (ux ) = (yuξ + uη )
∂x ∂x
∂ ∂
uxx = y · (uξ ) + (uη )
∂x ∂x
∂ξ ∂η ∂ξ ∂η
= y · (uξξ · + uξη · ) + (uηξ · + uηη · )
∂x ∂x ∂x ∂x
Física 6 TALLER 1.
Universidad de Córdoba EDP
x2 y 2 uξξ + 2x2 yuξη + x2 uηη − 2xyuξ − 2x2 y 2 uξξ + y 2 x2 uξξ + xyuξ + xuη + xyuξ = 0
Cancelamos términos similares:
x2 y 2 uξξ + 2x2 yuξη + x2 uηη − 2x2 y 2 uξξ + y 2 x2 uξξ − 2xyuξ + 2xyuξ + xuη = 0
Nos queda:
x2 uηη + xuη = 0
Dividimos por x2 :
uη
uηη +
=0
x
Ahora, recordemos que η = x, por lo tanto, tenemos:
uη
uηη = −
η
Sea uη = w, tenemos:
dw w
=−
dη η
Resolvemos por separación de variables:
dw dη
=−
w η
ln |w| = − ln |η| + C
Física 7 TALLER 1.
Universidad de Córdoba EDP
C
w=
η
Por lo tanto, tenemos:
f (ξ)
uη =
η
Finalmente, integramos:
u(ξ, η) = f (ξ) ln |η| + g(ξ)
Volviendo a las variables originales x e y, tenemos la solución final de la EDP:
5 Ejercicio 5
uxx − xuyy = 0
Solución
Sea:
A = 1, B = 0, C = −x
Nótese que si ∆ > 0, podemos decir que la EDP es hiperbólica (i). Nótese que si ∆ < 0, podemos decir que
la EDP es elíptica (ii).
Entonces:
i) si ∆ > 0, es hiperbólica y satisface la EDO:
√ √ √
dy B± ∆ 0 ± 4x 4x √
= = =± =± x
dx 2A 2(1) 2
C = y − 23 x3/2 y C = y + 23 x3/2 .
2 2
η = y − x3/2 ξ = y + x3/2
3 3
Aplicando la regla de la cadena para las nuevas variables:
Para ux : √ √
ux = uξ · ξx + uη · ηx ⇒ ux = uξ (− x) + uη ( x)
√ √
ux = − xuξ + xuη
Para uxx :
∂ √ √ 1 √ 1 √
uxx = (− xuξ + xuη ) ⇒ uxx = − √ uξ − x(uξξ ξx + uξη ηx ) + √ uη + x(uηξ ξx + uηη ηx )
∂x 2 x 2 x
1 √ √ √ 1 √ √ √
uxx = − √ uξ − x(uξξ (− x) + uξη x) + √ uη + x(uηξ (− x) + uηη x)
2 x 2 x
Física 8 TALLER 1.
Universidad de Córdoba EDP
1 1
uxx = − √ uξ + xuξξ − xuξη + √ uη + xuηη
2 x 2 x
Para uy :
uy = uξ · ξy + uη · ηy ⇒ uy = uξ (1) + uη (1)
uy = uξ + uη
Para uyy :
uyy = (uξ + uη )y ⇒ uyy = uξξ + 2uξη + uηη
Reemplazando en la EDP original:
1 1
− √ uξ + xuξξ − xuξη + √ uη + xuηη − x(uξξ + 2uξη + uηη ) = 0
2 x 2 x
1 1
− √ uξ + xuξξ − xuξη + √ uη + xuηη − xuξξ − 2xuξη − xuηη = 0
2 x 2 x
1 1
− √ uξ + √ uη − 3xuξη = 0
2 x 2 x
1
uξη = − √ (uξ − uη )
6x x
Esta es la forma canónica, sin embargo está expresada en términos de η, x e y. Para expresarlo solamente en
razón de las nuevas variables, consideremos:
Si:
2 2
η − ξ = y + x3/2 − y + x3/2
3 3
4
η − ξ = x3/2
3
Multiplicando por 6 la expresión:
6(η − ξ) = 8x3/2
Esta expresión, si la reemplazamos en la forma canónica antes calculada, nos da la EDP en su función de ξ
y η:
1
uξη = − (uξ − uη )
6(η − ξ)
ii) si ∆ < 0, es elíptica y satisface la EDO:
√ √
dy B±i ∆ 0 ± i 4x √
= = = ±i x
dx 2A 2(1)
2i 3/2 2i 3/2
C =y+ x y C =y− x
3 3
.
2i 3/2 2i 3/2
ℑ=y+ x τ =y− x
3 3
Física 9 TALLER 1.
Universidad de Córdoba EDP
Ahora, estas sustituciones son complejas. Para encontrar soluciones reales, podemos tomar un ξ, η de la
siguiente forma:
ℑ+τ ℑ−τ 2
ξ= =y η= = x3/2
2 2i 3
Aplicando la regla de la cadena:
Para ux : √
ux = uξ · ξx + uη · ηx ⇒ ux = uη · x
Para uxx :
1 √
uxx = √ uη + xuηη
2 x
Para uy :
uy = uξ
Para uyy :
uyy = uξξ
Reemplazando en la EDP original:
1 √
√ uη + xuηη − xuξξ = 0
2 x
Organizando:
1
uηη − uξξ +uη = 0
2x
Finalmente, expresamos todo en términos de las nuevas variables:
1
uηη − uξξ + uη = 0
3η
6 Ejercicio 6
uxx + x2 uyy = 0
Solución
Sea:
A = 1, B = 0, C = x2
Como ∆ < 0, podemos decir que la EDP es elíptica, la cual satisface la siguiente EDO:
p √
dy B ± i |∆| 0 ± −4x2
= = = ±ix
dx 2A 2(1)
Nos queda:
dy dy
= ix o = −ix
dx dx
Operando las EDOs resultantes:
i)
dy
= −ix ⇒ dy = −ixdx
dx
i
y = − x2 + C
2
Nuestra primera sustitución es entonces:
i
ℑ = y + x2
2
ii)
dy
= ix ⇒ dy = ixdx
dx
Física 10 TALLER 1.
Universidad de Córdoba EDP
i
y = x2 + C
2
Nuestra segunda sustitución es entonces:
i
τ = y − x2
2
Estas sustituciones encontradas son de carácter imaginario. Para encontrar sustituciones reales, tomemos los
siguientes ξ y η:
ℑ+τ
ξ= =y
2
ℑ−τ 1
η= = x2
2i 2
Aplicamos regla de la cadena para las sustituciones encontradas:
Para ux :
∂ξ ∂η
ux = uξ · + uη · = uξ · (0) + uη · (x)
∂x ∂x
ux = xuη
Para uxx :
∂
uxx = (xuη ) = uη + xuηη
∂x
Para uy :
∂ξ ∂η
uy = uξ · + uη · = uξ
∂y ∂y
Para uyy :
∂
uyy = (uξ ) = uξξ
∂y
Reemplazamos en la EDP original:
uη + xuηη + x2 uξξ = 0
Dividiendo por x2 y recordando que η = 12 x2 , tenemos la EDP en su forma canónica:
1
uηη + uξξ + uη = 0
2ξ
7 Ejercicio 7
y 2 uxx − x2 uyy = 0
Solución
Sean:
A = y2
B=0
∆ = B − 4(A)(C) = 4x2 y 2 > 0
2
C = −x2
Como ∆ > 0, podemos decir que la EDP es hiperbólica, la cual satisface la siguiente EDO:
√
dy B± ∆ dy x
= ⇒ =±
dx 2A dx y
Reorganizando e integrando, nos queda entonces:
Z Z
1 1
y dy = x dx ⇒ y 2 = x2 + C
2 2
Física 11 TALLER 1.
Universidad de Córdoba EDP
C = y 2 − x2 ⇒ ξ = y 2 − x2
Haciendo el mismo procedimiento con la otra solución, obtenemos que la segunda sustitución es:
C = y 2 + x2 ⇒ η = y 2 + x2
De esta forma, aplicaremos la regla de la cadena asociada a nuestras dos nuevas sustituciones para obtener
la forma canónica asociada.
Para ux :
∂ξ ∂η
ux = uξ · + uη · ⇒ ux = uξ (−2x) + uη (2x)
∂x ∂x
ux = −2xuξ + 2xuη
Para uxx :
∂ ∂
uxx = (ux ) ⇒ uxx = (−2xuξ + 2xuη )
∂x ∂x
∂ξ ∂η ∂ξ ∂η
uxx = −2(uξ + x(uξξ · + uξη · )) + 2(uη + x(uηξ · + uηη · ))
∂x ∂x ∂x ∂x
uxx = −2(uξ + x(uξξ (−2x) + uξη (2x))) + 2(uη + x(uηξ (−2x) + uηη (2x)))
uxx = −2uξ + 4x2 uξξ − 4x2 uξη + 2uη − 4x2 uηξ + 4x2 uηη
Para uy :
∂ξ ∂η
uy = uξ · + uη · ⇒ uy = uξ (2y) + uη (2y)
∂y ∂y
uy = 2yuξ + 2yuη
Para uyy :
∂ ∂
uyy = (uy ) ⇒ uyy = (2yuξ + 2yuη )
∂y ∂y
∂ξ ∂η ∂ξ ∂η
uyy = 2(uξ + y(uξξ · + uξη · )) + 2(uη + y(uηξ · + uηη · ))
∂y ∂y ∂y ∂y
y 2 (−2uξ + 4x2 uξξ − 4x2 uξη + 2uη − 4x2 uηξ + 4x2 uηη ) − x2 (2uξ + 4y 2 uξξ + 4y 2 uξη + 2uη + 4y 2 uηξ + 4y 2 uηη ) = 0
Reduciendo términos:
−2y 2 uξ + 4x2 y 2 uξξ − 4x2 y 2 uξη + 2y 2 uη − 4x2 y 2 uηξ + 4x2 y 2 uηη − ...
2x2 uξ − 4x2 y 2 uξξ − 4x2 y 2 uξη − 2x2 uη − 4x2 y 2 uηξ − 4x2 y 2 uηη = 0
Nos queda:
−16y 2 x2 uξη − 2x2 uξ − 2x2 uη = 0
Física 12 TALLER 1.
Universidad de Córdoba EDP
8 Ejercicio 8
ξ = y − ln | sin(x)|
Ahora, por ser este un caso parabólico, se puede escoger una variable η = x, la cual para efectos prácticos
será nuestra segunda sustitución:
η=x
Aplicando regla de la cadena para ξ y η:
ηx = 1, ξx = − cot(x)
Física 13 TALLER 1.
Universidad de Córdoba EDP
ηy = 0, ξy = 1
Para ux :
ux = uξ · ξx + uη · ηx = uξ (− cot(x)) + uη (1)
ux = − cot(x)uξ + uη
Para uxx :
∂
uxx = (− cot(x)uξ + uη ) = csc2 (x)uξ − cot(x) (uξξ · ξx + uξη · ηx ) + uηξ · ξx + uηη · ηx
∂x
uxx = csc2 (x)uξ − cot(x) (uξξ (− cot(x)) + uξη · (1)) + uηξ (− cot(x)) + uηη (1)
uxx = csc2 (x)uξ + cot2 (x)uξξ − cot(x)uξη − cot(x)uηξ + uηη
Para uy :
uy = uξ · ξy + uη · ηy = uξ
Para uyy :
∂
uyy = (uξ ) = uξξ
∂y
Para uyx :
∂ ∂
uyx = (uy ) = (uξ ) = uξξ · ξx + uξη · ηx
∂x ∂x
uyx = uξξ (− cot(x)) + uξη
Reemplazando en la EDP original:
sin2 (x) csc2 (x)uξ + cot2 (x)uξξ − cot(x)uξη − cot(x)uηξ + uηη + sin(2x) (uξξ (− cot(x)) + uξη ) + cos2 (x)uξξ = 0
Simplificando:
uξ + cos2 (x)uξξ − 2 sin(x) cos(x)uξη + sin2 (x)uηη = 0
Esta es finalmente la forma canónica. Sin embargo, aún podemos expresarla en términos solo de las sustitu-
ciones ξ y η, recordando que η = x, nos queda entonces:
sin2 (η)uηη + uξ = 0
uξ
uηη = −
sin2 η
9 Ejercicio 9 (23)
(t − x)
ux + ut = 0
t+x
tenemos por sustitución
dx x−t x
= ⇒ =w ⇒ x = wt
dt t+x t
x
dx t −1 dx
= x ⇒ = wt t + w
dt t +1 dt
Reemplazamos, obtenemos:
Física 14 TALLER 1.
Universidad de Córdoba EDP
dw w−1
+w =
dt w+1
dw w−1 w − 1 − w(w + 1)
⇒ = −w →
dt w+1 w+1
dw w − 1 − w(w + 1) dw −(1 + w2 )
⇒ t= ⇒ t=
dt w+1 dt w+1
(w + 1) dt w 1 dt
2
dw = − ⇒ ( 2
+ 2
)dw = −
1+w t 1+w 1+w t
Integrando:
Z Z
w 1 1
2
+ dw = − dt
1+w 1 + w2 t
1
⇒ ln(1 + w2 ) + tan−1 (w) = − ln |t| + C
2
Reemplazamos w:
x2
1 x
ln 1 + 2 + tan−1 = − ln |t| + C
2 t t
1 1 x
⇒ ln |x2 + t2 | − ln |t2 | + tan−1 = − ln |t| + C
2 2 t
1 x
⇒ ln |x2 + t2 | + tan−1 +C =0
2 t
tenemos:
1 x
u(x, t) = ln |x2 + t2 | + tan−1 +C
2 t
Finalmente, reemplazando u(x, 0) = f (x):
1 x
u(x, 0) = ln |x2 + 02 | + tan−1
2 0
x
Valor indefinido
0
10 Ejercicio 10 (5)
a = 1, b = −2, c=1
= 1 − 4 = −3 < 0
Por lo tanto, la EDP es de tipo hiperbólica.
Física 15 TALLER 1.
Universidad de Córdoba EDP
√ √
2y = −4 ± 12x + c y ξ = −4 + 12x + 2y
√ √ √
η =4− 12x + 2y y ξx = 4 + 12, ηx = 4 − 12
ξxx = 0, ηxx = 0 y ξy = 2, ηy = −2
Física 16 TALLER 1.