100% found this document useful (1 vote)
450 views34 pages

EDO2

1. Se explica que cosh(x) y senh(x) son soluciones de la ecuación diferencial lineal homogénea porque pueden expresarse como combinaciones lineales de las soluciones ex y e-x, las cuales forman un conjunto de funciones linealmente independientes. 2. Se usa la técnica de reducción de orden para encontrar una segunda solución xe2x de la ecuación diferencial y'' - 4y' + 4y = 0, cuya primera solución conocida es e2x. 3. Se aplica la misma técnica para hall

Uploaded by

Jesús Villegas
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
100% found this document useful (1 vote)
450 views34 pages

EDO2

1. Se explica que cosh(x) y senh(x) son soluciones de la ecuación diferencial lineal homogénea porque pueden expresarse como combinaciones lineales de las soluciones ex y e-x, las cuales forman un conjunto de funciones linealmente independientes. 2. Se usa la técnica de reducción de orden para encontrar una segunda solución xe2x de la ecuación diferencial y'' - 4y' + 4y = 0, cuya primera solución conocida es e2x. 3. Se aplica la misma técnica para hall

Uploaded by

Jesús Villegas
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 34

1.

Ejercicio 1
Suponga que y1 = ex y y2 = e−x son dos soluciones de una ecuación
diferencial lineal homogénea. Explique por qué y3 = cosh(x) y y4 = senh(x)
son también soluciones de la ecuación.

y1 , y2 satisfacen la EDO lineal homogénea si y solo si son linealmente indepen-


dientes entre si, esto implica que su Wronskiano debe ser distinto de 0

y1 = ex , y10 = ex , y100 = ex , y1000 = ex

y2 = e−x , y20 = −e−x , y200 = e−x , y2000 = −e−x



y1 y2
W (y1 , y2 ) = 0 6= 0
y1 y20
x −x

e e
⇒ W (y1 , y2 ) = x 6= 0
e −e−x
Como evidencia y por inspección ya que
ex 6= ∀x ∈ (−∞, +∞) podrı́amos asimilar que W (y1 , y2 ) 6= 0, y1 , y2 son lineal-
mente independientes formando un C.F.I.
Para demostrarlo desarrollando los productos del determinante
x
−x
e e
⇒ W (y1 , y2 ) = x 6= 0 = ex ∗ (−e−x ) − ex ∗ e−x = −1 − 1 = −2
e −e−x
Ya que -2 es una cantidad real siempre será distinta de cero por lo tanto y1 y
y2 forman un C.F.I
y(x) = c1 ex + c2 e−x
Si y3 = cosh(x) y y4 = senh(x) bien podrı́amos probar que forman parte
del C.F.I y desarrollar un Wronskiano 4 * 4 bastante largo en cálculos, pero
ex + e−x ex e−x
por inspección sabiendo que cosh(x) = = + y senh(x) =
2 2 2
ex − e−x ex e−x
= − ya que 12 y − 12 son constantes si tomamos y(x) con
2 2 2
c1 = c2 = 21 obtenemos

1 1
y(x) = ex + e−x = cosh(x)
2 2
1
si tomamos y(x) con c1 = 2
y c2 = − 12 obtenemos

1 1
y(x) = ex − e−x = senh(x)
2 2
Por lo tanto y3 = cosh(x) y y4 = senh(x) son soluciones de la ecuación
diferencial.

1
2. Ejercicio 2
Use la reducción de orden para encontrar una segunda solución:

(1) y 00 − 4y 0 + 4y = 0; y1 = e2x
y10 (x)
= 2e2x y100 (x) = 4e2x
Considerando y2 (x) = u(x)y1 (x)

y20 (x) = u0 (x)y1 (x) + u(x)y10 (x)


= e2x u0 (x) + 2e2x u(x)
= e2x [u0 (x) + 2u(x)]
y200 (x) = u00 (x)y1 (x) + 2u0 (x)y10 (x) + u(x)y100 (x)
= e2x u00 (x) + 4e2x u0 (x) + 4e2x u(x)
= e2x [u00 (x) + 4u0 (x) + 4u(x)]
Sustituyendo en (1)

e2x [u00 (x) + 4u0 (x) + 4u(x)] − 4e2x [u0 (x) + 2u(x)] + 4e2x u(x) = 0

Factor común e2x

e2x [u00 (x) + 4u0 (x) + 4u(x) − 4 (u0 (x) + 2u(x)) + 4u(x)] = 0

e2x [u00 (x) + 4u0 (x) + 4u(x) − 4u0 (x) − 8u(x) + 4u(x)] = 0
e2x [u00 (x))] = 0
Por teorema del factor nulo e2 x = 0 o u00 (x) ya que e2x 6= 0 ∀x ∈ (−∞, +∞)

u00 (x) = 0

Considerando w(x) = u0 (x) ⇒ w0 (x) = u00 (x)

w0 (x) = 0
dw
=0
dx
Z Z
dw = 0dx + c

w(x) = c
u0 (x) = w(x)
du
= w(x)
dx
Z Z
du = c dx + c1

2
u(x) = cx + c1
Considerando c = 1 y c1 = 0

u(x) = x

y2 (x) = xe2x
∵ y(x) = c1 e2x + c2 xe2x Forman un C.F.S

(2) x2 y 00 + 2xy 0 − 6y = 0; y1 = x2

Considerando
y2 (x) = u(x)y1 (x) = x2 u(x)
y10 (x) = 2x , y100 (x) = 2
y20 (x) = u0 (x)y1 (x) + u(x)y10 (x) = x2 u0 (x) + 2xu(x)
y20 (x) = x [xu0 (x) + 2u(x)]
y200 (x) = u00 (x)y1 (x) + 2u0 (x)y10 (x) + u(x)y100 (x) = x2 u00 (x) + 4xu0 (x) + 2u(x)
Sustituyendo en (2)

x2 x2 u00 (x) + 4xu0 (x) + 2u(x) + 2x2 [xu0 (x) + 2u(x)] − 6x2 u(x) = 0
 

 

x2 x2 u00 (x) + 4xu0 (x) + 2xu0 (x) + 2u(x) + 4u(x) − 6u(x) = 0
| {z }
0

x4 u00 (x) + 6x3 u0 (x) = 0


Reescribiendo en su forma estándar
6
u00 + u0 = 0
x
Considerando w = u0 ⇒ w0 = u00
6
w0 + w = 0
x
6
w0 = − w
x
dw 6
=− w
dx x
Z Z
dw dx
= −6
w x

ln |w| = −6ln |x| + c0

3
eln|w| = e−ln|x |+c0
6

c0
w= 6
x
Ya que u0 = w ⇒ u = wdx
R

Z
dx c0
u = c0 6
= − x−5 + c1 Si c0 = −5 y c1 = 0 u = x−5
x 5

y2 (x) = u(x)y1 (x) = x−5 x2 = x−3


∴ y(x) = c1 x2 + c2 x−3 Forman un C.F.S

(3) xy 00 + y 0 = 0 ; y1 (x) = ln(x)


1 1
y10 (x) = , y100 (x) = − 2
x x
u 2 1
y2 = uy1 = uln(x) , y10 = u0 ln(x) + , y200 = u00 ln(x) + u0 − 2 u
x x x
Sustituyendo en (3)
 
2 0 1 u
x u ln(x) + u − 2 u + u0 ln(x) + = 0
00
x x x
1 u
xln(x)u00 + 2u0 − u + u0 ln(x) + = 0
x x
00 0 0
xln(x)u + 2u + u ln(x) = 0
xln(x)u00 + (2 + ln(x)) u0 = 0
Reescribiendo la ED en su forma estándar
 
00 2 + ln(x)
u + u0 = 0
xln(x)

Considerando w = u0 ⇒ w0 = u00
 
0 2 + ln(x)
⇒w + w=0
xln(x)
 
dw 2 1
=− + w
dx xln(x) x
Z Z Z
dw dx dx
= −2 −
w xln(x) x
Z
dx
ln |w| = −2 −ln |x|
xln(x)
| {z }
I0

4
Integrando I0 por partes
dx
u = ln(x) ⇒ du =
x
Z
du
⇒ I0 = = ln |u| = ln |ln(x)|
u
1


x
⇒ ln |w| = −2I0 − ln |x| = −2ln |ln(x)| − ln |x| + c0 = ln 2 + c0
ln (x)


1
ln |w| = ln + c0
xln2 (x)

1
ln

ln|w| 2 +c0
e =e xln (x)

c0
w=
xln2 (x)
Ya que u0 = w ⇒ u =
R
wdx
Z
dx
u = c0
xln2 (x)
dx
Integrando por partes u = ln(x) ⇒ u2 = ln2 (x) ⇒ du = x
Z
du co c0
u = c0 2
= − + c1 = − + c1
u u ln(x)

Si c0 = 1 y c1 = 0
1
u=−
ln(x)
ln(x)
y2 (x) = uy1 = − = −1
ln(x)
∴ y(x) = c1 ln(x) − c2 Es un C.F.S para la ED

3. Ejercicio 3
3.1. A)
A)Hallar la solución de las EDO de Orden superior. Use el método de
Coeficientes Indeterminados para determinar la solución
d2 y dy
a) dx 2 − 3 dx + 2y = x
00 0
y − 3y + 2y = x

Resolviendo la parte homogénea

5
00 0
y − 3y + 2y = 0
⇒ λ2 − 3λ + 2 = 0
(λ − 2)(λ − 1) = 0
λ1,2 = 2, 1
∴ yH (x) = C1 e2x + C2 ex
Parte particular

yp (x) = Ax + B
0
yp (x) = A
00
yp (x) = 0
Sustituyendo en la ED

0 − 3A + 2(Ax + B) = x
2Ax + 2B − 3A = x
1
⇒ 2A = 1 ⇒ A =
2
3
⇒ 2B − 3A = 0 ⇒ 2B = 3A ⇒ B =
4

1 3
∴ yp (x) = x +
2 4
Entonces:
1 3
y(x) = C1 e2x + C2 ex + x +
2 4

d3 x 2 000 00
b) dt3
− 2 ddt2x + x = t2 − 3 = x − 2x + x = t2 − 3
Haciendo uso de la ecuación auxiliar para resolver la parte homogénea
000 00
x − 2x + x = 0
⇒ λ3 − 2λ2 + 1 = 0
La raiz esp?a es: λ = 1
Entonces:
2
(λ − 1)(λ − λ − 1)

1± 5
∴ λ1 = 1, λ2,3 =
2
 √   √ 
1+ 5 1− 5
t t t
∴ xH (t) = C1 e + C2 e 2
+ C3 e 2

6
Parte particular

xp (t) = At2 + Bt + C
ẋp (t) = 2At + B
ẍp (t) = 2A
...
x p (t) = 0

Sustituyendo en la ED

0 − 2(2A) + At2 + Bt + C = t2 − 3
At2 + Bt − 4A + C = t2 − 3
⇒ A = 1, B = 0
⇒ −4A + C = −3 ⇒ C = −3 + 4 = 1

∴ xp (t) = t2 − 1

Entonces:
 √   √ 
1+ 5 1− 5
t t t
x(t) = C1 e + C2 e 2
+ C3 e 2
+ t2 − 1

( 00 0
y = −2y + 3y + e−x
c) 0
y(0) = 1, y (0) = 2

Reescribiendo la ecuaci?n

00 0
y + 2y − 3y = e−x

Resolviendo la parte homog?nea

00 0
y + 2y − 3y = 0
⇒ λ2 + 2λ − 3 = 0
(λ + 3)(λ − 1) = 0
λ1,2 = −3, 1

∴ yH (x) = C1 e−3x + C2 ex

Parte particular

7
P m(x)eαx = e−x
α = −1, λ1,2 = −3, 1
yp (x) = Ae−x
0
yp (x) = −Ae−x
00
yp (x) = Ae−x

Sustituyendo en la ED

Ae−x − 2(Ae−x ) − 3(Ae−x ) = e−x


−4Ae−x = e−x
1
⇒ −4A = 1 ⇒ A = −
4

1
∴ yp (x) = − e−x
4
Entonces:

1
y(x) = C1 e−3x + C2 ex − e−x
4
Del PVI:
1
y(x) = C1 e−3x + C2 ex − e−x
4
0 1
y (x) = −3C1 e−3x + C2 ex + e−x
4
1
y(0) = C1 + C2 − = 1 ........1
4
0 1
y (0) = −3C1 + C2 + = 2 ........2
4
De 3×1+ 2
3
3C1 + 3C2 − = 3
4
+
1
−3C1 + C2 + =2
4
11 11
⇒ 4C2 = ⇒ C2 =
2 8
De 1
11 1 1
C1 + − = 1 ⇒ C1 = −
8 4 8
1 11 1
∴ y(x) = − e−3x + ex − e−x
8 8 4

8
d6 x d4 x
d) − 4 = 1 = x(6) − x(4) = 1
dt6 dt
Resolviendo la parte homog?nea

x(6) − x(4) = 0
⇒ λ6 − λ4 = 0
λ4 (λ2 − 1) = 0
λ1,2,3,4 = 0, λ5,6 = ±1

∴ xH (t) = C1 + C2 t + C3 t2 + C4 t3 + C5 et + C6 e−t

Parte particular

xp (t) = At4
0
xp (t) = 4At3
00
xp (t) = 12At2
000
xp (t) = 24At
x(4)
p (t) = 24A

x(5)
p (t) = 0

x(6)
p (t) = 0

Sustituyendo en la ED

1
0 − 24A = 1 ⇒ A = −
24
1 4
∴ xp (t) = − t
24
Entonces:

1 4
x(t) = C1 + C2 t + C3 t2 + C4 t3 + C5 et + C6 e−t − t
24
Es la Familia de soluciones para la ED

2 00 0
e) ddt2x − 4 dx
dt
+ 4x = et + e2t + 1 = x − 4x + 4x = et + e2t + 1

Resolviendo la parte homogénea

9
00 0
x − 4x + 4x = 0
⇒ λ2 − 4λ2 + 4 = 0
(λ − 2)2 = 0
λ1,2 = 2
∴ xH (t) = C1 e + C2 te2t
2t

Parte particular 1

xp1 (t) = Aet


0
xp1 (t) = Aet
00
xp1 (t) = Aet

Sustituyendo en la ED

Aet − 4(Aet ) + 4Aet = et


Aet = et ⇒ A = 1
∴ xp1 (t) = et

Parte particular 2

λ1,2 = α = 2
xp2 (t) = t2 Ae2t

0
xp2 (t) = 2tAe2t + 2t2 Ae2t = Ae2t 2t + 2t2

00
xp2 (t) = (2 + 4t) Ae2t + Ae2t 4t + 4t2 = Ae2t 4t2 + 8t + 2
 

Sustituyendo en la ED

Ae2t 4t2 + 8t + 2 − 4Ae2t 2t + 2t2 + 4t2 Ae2t = e2t


  

Ae2t 4t2 + 8t + 2 − 8t − 8t2 + 4t2 = e2t




1
2Ae2t = e2t ⇒ A =
2

1
∴ xp2 (t) = t2 e2t
2
Parte particular 3

xp3 (t) = A
0
xp3 (t) = 0
00
xp3 (t) = 0

10
Sustituyendo en la ED

0 + 4(0) + 4A = 1
1
4A = 1 ⇒ A =
4
1
∴ xp3 (t) =
4
Entonces:

1 1
x(t) = C1 e2t + C2 te2t + et + t2 e2t +
2 4

f)ẍ + x = sen(t) − cos(2t)


Resolviendo la parte homogénea

ẍ + x = 0
λ2 + 1 = 0
λ1,2 = ±i, α = 0, β = 1

∴ xH (t) = C1 cos(x) + C2 sen(x)

Parte particular 1

ẍ + x = (0)cos(t) + sen(t)
α = 0, β = 1
xp1 (t) = t (Acos(t) + Bsen(t))
0
xp1 (t) = Acos(t) + Bsen(t) + t(Bcos(t) − Asen(t))
00
xp1 (t) = Bcos(t) − Asen(t) + Bcos(t) − Asen(t) − t(Bsen(t) + Acos(t))
= −2Asen(t) + 2Bcos(t) − t(Acos(t) + Bsen(t))

Sustituyendo en la ED

−2Asen(t) + 2Bcos(t) − t(Acos(t) + Bsen(t)) + t(Acos(t) + Bsen(t)) = (0)cos(t) + sen(t)


−2Asen(t) + 2Bcos(t) = (0)cos(t) + sen(t)
1
−2A = 1 ⇒ A = − , B = 0
2
1
∴ xp1 (t) = − t cos(t)
2

11
Parte particular 2

ẍ + x = −cos(2t) + (0)sen(2t)
α = 0, β = 2
xp2 (t) = Acos(2t) + Bsen(2t)
0
xp2 (t) = −2Asen(2t) + 2Bcos(2t)
00
xp2 (t) = −4Acos(2t) − 4Bsen(2t)

Sustituyendo en la ED

−4Acos(2t) − 4Bsen(2t) + Acos(2t) + Bsen(2t) = −cos(2t) + (0)sen(2t)


−3Acos(2t) − 3Bsen(2t) = −cos(2t) + (0)sen(2t)
1
−3A = −1 ⇒ A = , B = 0
3
1
∴ xp2 (t) = cos(2t)
3
Entonces:

1 1
x(t) = C1 cos(x) + C2 sen(x) − t cos(t) + cos(2t)
2 3

3.2. B)
B)Use el método de Variación de Parámetros para determinar la solución no
homogénea:
d2 y dy 00 0
a) 2 + 4 + 5y = cos(x) = y + 4y + 5y = cos(x)
dx dx
Resolviendo la parte homogénea

00 0
y + 4y + 5y = 0
⇒ λ2 + 4λ + 5 = 0
λ1,2 = −2 ± i
α = −2, β = 1
−2x
∴ yH (x) = e (C1 cos(x) + C2 sen(x))

y1 (x) = e−2x cos(x) , y2 (x) = e−2x sen(x), f (x) = cos(x)

12
y10 (x) = −e−2x (2cos(x) + sen(x)) , y20 (x) = e−2x (cos(x) − 2sen(x))
−2x −2x

e cos(x) e sen(x)
W = −2x −2x

−e (2cos(x) + sen(x)) e (cos(x) − 2sen(x))
W = e−4x

0 e−2x sen(x)
W1 =
−2x

cos(x) e (cos(x) − 2sen(x))
W1 = −e−2x sen(x)cos(x)

e−2x cos(x) 0
W2 = −2x

−e (2cos(x) + sen(x)) cos(x)
W2 = e−2x cos2 (x)

W1 e−2x sen(x)cos(x) 1
C10 (x) = =− −4x
= −e2x cos(x)sen(x) ya que cos(x)sen(x) = sen(2x)
W e 2
Z
1
C1 (x) = − e2x sen(2x)dx
2
| {z }
I0

Integrando I0 por partes


1
u = sen(2x) ⇒ du = 2cos(2x)dx , dv = e2x dx ⇒ v = e2x
2
Z
1
I0 = e2x sen(2x) − e2x cos(2x)dx
2
| {z }
I1

Resolviendo para I1 por partes


1
u = cos(2x) ⇒ du = −2sen(2x)dx , dv = e2x dx ⇒ v = e2x
2
Z
1
I1 = e2x cos(2x) + e2x sen(2x)dx
2
| {z }
I0

Sustituyendo I1 en I0
1
I0 = e2x sen(2x) − I1
2
1 1
I0 = e2x sen(2x) − e2x cos(2x) + I0
2 2
1 2x
I0 = e (sen(2x) − cos(2x))
4
1
C1 (x) = e2x (cos(2x) − sen(2x))
8

13
W2 e−2x cos2 (x) 1
C20 (x) = = = e2x
cos 2
(x) ya que cos 2
(x) = (1 + cos(2x))
W e−4x 2
1
C20 (x) = e2x (1 + cos(2x))
2
Z
C2 (x) = e2x cos2 (2x)dx
| {z }
I0
Z Z
1 2x 1
I0 = e (1 + cos(2x)) dx = e2x + e2x cos(2x)dx
2 2
 
Z Z
1 
I0 =  e dx + e2x cos(2x)dx
 2x
2 
| {z }
I1
 
1 1 2x
I0 = e dx + I1
2 2
1
Integrando por partes I1 u = cos(2x) ⇒ du = −2sen(2x)dx , dv = e2x dx ⇒ v = e2x
2
Z
1
I1 = e2x cos(2x) + e2x sen(2x)dx
2
| {z }
I2

Integrando I2 por partes


1
u = sen(2x) ⇒ du = 2cos(2x)dx , dv = e2x dx ⇒ v = e2x
2
Z
1
I2 = e2x sen(2x) − e2x cos(2x)dx
2
| {z }
I1

Sustituyendo I2 en I1
1
I1 = e2x cos(2x) + I2
2
1 1
I1 = e2x cos(2x) + e2x sen(2x) − I1
2 2
1
2I1 = e2x (cos(2x) + sen(2x))
2
1
I1 = e2x (cos(2x) + sen(2x))
4
 
1 1 2x
I0 = e dx + I1
2 2

14
 
1 1 2x 1 2x
I0 = e dx + e (cos(2x) + sen(2x))
2 2 4
1
I0 = e2x (sen(2x) + cos(2x) + 2)
8
1
∴ C2 (x) = e2x (cos(2x) + sen(2x) + 2)
8

yp (x) = C1 (x)y1 (x) + C2 (x)y2 (x)


1 1
yp (x) = e2x e−2x cos(x) (cos(2x) − sen(2x))+ e2x e−2x sen(x) (cos(2x) + sen(2x) + 2)
 
8 8
1
yp (x) = {cos(x) [(cos(2x) − sen(2x))] + sen(x) [(cos(2x) + sen(2x) + 2)]}
8
1
yp (x) = [cos(2x)cos(x) − cos(x)sen(2x) + sen(x)cos(2x) + sen(2x)sen(x) + 2sen(x)]
8
 
1
yp (x) = cos(2x)cos(x) + sen(2x)sen(x) + sen(x)cos(2x) − cos(x)sen(2x) +2sen(x)

8 | {z } | {z }
cos(2x−x) sen(x−2x

1
∴ yp (x) = [cos(x) + sen(x)]
8
d3 y 2
d y dy
b) dx3
− 3 dx2 + 3 dx − y = x − 2

000 00 0
y − 3y + 3y − y = x − 2

Resolviendo la parte homogénea

000 00 0
y − 3y + 3y − y = 0
⇒ λ3 − 3λ2 + 3λ − 1 = 0
(λ − 1)3 = 0
λ1,2,3 = 1
∴ yH (x) = C1 e + C2 xe + C3 x2 ex
x x

y1 = ex , y2 = xex , y3 = x2 ex

x x 2 x

e xe x e
x x
w = e e (x + 1) xe (x + 2) = 2e3x
x
ex ex (x + 2) ex (x2 + 4x + 2)

15

0
xex x2 ex

x x
w1 = 0
e (x + 1) xe (x + 2)
x − 2 ex (x + 2) ex (x2 + 4x + 2)
w1 = x1 e2x (x − 2)

x
e
x 0 x2 e x

w2 = e 0 xex (x + 2)
ex x − 2 ex (x2 + 4x + 2)
=x2 e2x (x − 2) − xe2x (x2 − 4)
=e2x (x3 − 2x2 ) − e2x (x3 − 4x)
=e2x (x3 − x3 − 2x2 + 4x)
=2xe2x (−x + 2)

x x

e xe 0
x x
w3 = e e (x + 1) 0
ex ex (x + 2) x − 2
=e2x (x + 1)(x − 2) − xe2x (x − 2)
=e2x (x2 − x + 2) − e2x (x2 − 2x)
=e2x (x2 − x2 − x + 2x − 2)
=e2x (x − 2)

Entonces:

x2 e2x (x − 2) x2 (x − 2)
C10 (x) = =
2e3x 2ex
2xe2x (2 − x) x(2 − x)
C20 (x) = 3x
=
2e ex
2x
e (x − 2) (x − 2)
C30 (x) = 3x
=
2e 2ex

x2 (x − 2)
Z Z
1
C1 (x) = x
dx = (x3 − 2x2 )(e−x )dx
2e 2
Integrando por partes

16
u = x3 − 2x2 du = 3x2 − 4x
dv = e−x v = −e−x

Z  Z 
1 3 2 −x 1 3 2 −x −x 2
(x − 2x )(e )dx = (x − 2x )(−e ) − (−e )(3x − 4x)dx
2 2

Integrando por partes

u = 3x2 − 4x du = 6x − 4
dv = −e−x v = e−x

 Z 
1 3 2 −x −x 2
(x − 2x )(−e ) − (−e )(3x − 4x)dx =
2
  Z 
1 3 2 −x 2 −x −x
(x − 2x )(−e ) − (3x − 4x)(e ) − (e )(6x − 4)dx
2

Integrando por partes

u = 6x − 4 du = 6
dv = e−x v = −e−x

  Z 
1 3 2 −x 2 −x −x
(x − 2x )(−e ) − (3x − 4x)(e ) − (e )(6x − 4)dx
2
   Z 
1 3 2 −x 2 −x −x −x
= (x − 2x )(−e ) − (3x − 4x)(e ) − (6x − 4)(−e ) − 6(−e )dx
2
1 3
(x − 2x2 )(−e−x ) − (3x2 − 4x)(e−x ) + (6x − 4)(e−x ) + 6e−x
 
=
2
1 3
(x − 2x2 )(−e−x ) − (3x2 − 4x)(e−x ) − (6x − 4)(e−x ) − 6e−x

=
2
1 1 1 1
= − e−x (x3 − 2x2 ) − e−x (3x2 − 4x) − e−x (6x − 4) − e−x (6)
2 2 2 2
1 −x 3
= − e (x − 2x2 + 3x2 − 4x + 6x − 4 + 6)
2
1
= − e−x (x3 + x2 + 2x + 2) = C1 (x)
2

x(2 − x)
Z Z
C2 (x) = dx = (2x − x2 )(e−x )dx
ex
Integrando por partes

17
u = 2x − x2 du = 2 − 2x
dv = e−x v = −e−x

Z
(2x − x2 )(e−x )dx
Z
−x
= (2x − x )(−e ) − 2
(2 − 2x)(−e−x )dx

Integrando por partes

u = 2 − 2x du = −2
dv = −e−x v = e−x

Z
−x
2
(2x − x )(−e ) − (2 − 2x)(−e−x )dx
 Z 
2 −x −x −x
= (2x − x )(−e ) − (2 − 2x)(e ) − −2e dx

= (2x − x2 )(−e−x ) − (2 − 2x)(e−x ) − 2e−x


 

= (2x − x2 )(−e−x ) − (2 − 2x)(e−x ) + 2e−x


= e−x (−2x + x2 − 2 + 2x + 2)
= e−x (x2 ) = C2 (x)

(x − 2)
Z Z
1
C3 (x) = x
dx = (x − 2)(e−x )dx
2e 2
Integrando por partes

u = x − 2 du = 1
dv = e−x v = −e−x

Z
1
(x − 2)(e−x )dx
2
 Z 
1 −x −x
= (x − 2)(−e ) − (−e )dx
2
1 1
= (x − 2)(−e−x ) − (e−x )
2 2
1 −x
= e (−x + 2 − 1)
2
1
= − e−x (x − 1) = C3 (x)
2

18
Entonces:
1 1
yp (x) = − (x3 + x2 + 2x + 2) + x3 − x2 (x − 1)
2 2
1 3
= − (x + x2 + 2x + 2 − 2x3 + x3 − x2 )
2
1
= − (2x + 2) = −(x + 1)
2

∴ y(x) = C1 ex + C2 xex + C3 x2 ex − (x + 1)
( 00 0
y − 6y + 10y = 0
c)
y(0) = 10, y(0) = 5
Resolviendo la parte homogénea
00 0
y − 6y + 10y = 0
⇒ λ2 − 6λ + 10 = 0
λ1,2 = 3 ± i
α = 3, β = 1

∴yH (x) = e3x (C1 cos(x) + C2 sen(x))

y1 = e3x cos(x), y2 = e3x sen(x)

3x 3x

e cos(x) e sen(x)
w = 3x 3x 3x 3x

−e sen(x) + 3e cos(x) e cos(x) + 3e sen(x)

= e3x cos(x)(e3x cos(x) + 3e3x sen(x)) − e3x sen(x)(−e3x sen(x) + 3e3x cos(x))
= e6x cos2 (x) + 3e6x cos(x)sen(x) + e6x sen2 (x) − 3e6x cos(x)sen(x)
= e6x (cos2 (x) + sen2 (x)) = e6x
3x

0 e sen(x)
w1 = 3x 3x
= −100e3x sen(x)
100 e cos(x) + 3e sen(x)
3x

e cos(x) 0
= 100e3x cos(x)
w2 = 3x 3x
−e sen(x) + 3e cos(x) 100

19
Entonces:
−100e3x sen(x)
C10 (x) = = −100e−3x sen(x)
e6x
100e3x cos(x)
C20 (x) = = 100e−3x cos(x)
e6x

Z Z
−3x
C1 (x) = (−100e sen(x))dx = 100 (−e−3x sen(x))dx

Integrando por partes

u = sen(x) du = cos(x)
1
dv = −e−3x v = e−3x
3
Z
100 (−e−3x sen(x))dx
 Z 
1 −3x 1 −3x
= 100 e sen(x) − e cos(x)dx
3 3

Integrando por partes

u = cos(x) du = −sen(x)
1
dv = −e−3x v = − e−3x
3
 Z 
1 1
100 e−3x sen(x) − e−3x cos(x)dx
3 3
  Z 
1 −3x 1 1 −3x 1 −3x
= 100 e sen(x) − − e cos(x) + −(e )sen(x)dx
3 3 3 3
Z
100 −3x 100 −3x 100
= e sen(x) + e cos(x) − −(e−3x )sen(x)dx
3 9 9

Z Z
−3x 100 −3x 100 −3x 100
⇒ 100 (−e sen(x))dx = e sen(x) + e cos(x) − −(e−3x )sen(x)dx
3 9 9
Z
10 1 1
⇒ (−e−3x sen(x))dx = e−3x sen(x) + e−3x cos(x)
9 3 9
Z
3 1
⇒ −(e−3x )sen(x)dx = e−3x sen(x) + e−3x cos(x)
10 10

Z
∴ 100 (−e−3x sen(x))dx = 30e−3x sen(x) + 10e−3x cos(x) = C1 (x)

20
Z Z  
−3x 1 −3x 1 −3x
De 100 (−e sen(x))dx = 100 e sen(x) − e cos(x)dx
3 3
Z Z
−3x 1 −3x 1
⇒ e sen(x) = e sen(x) − e−3x cos(x)dx
3 3
Z Z
⇒ 3 e sen(x) = e sen(x) − e−3x cos(x)dx
−3x −3x

Z  
−3x −3x 3 −3x 1 −3x
e cos(x)dx = e sen(x) − 3 e sen(x) + e cos(x)
10 10
9 3
= e−3x sen(x) − e−3x sen(x) − e−3x cos(x)
10 10
1 −3x 3 −3x
= e sen(x) − e cos(x)
10 10

Z Z
C2 (x) = 100e cos(x)dx = 100 e−3x cos(x)dx
−3x

 
1 −3x 3 −3x
= 100 e sen(x) − e cos(x)
10 10
−3x −3x
= 10e sen(x) − 30e cos(x)
Entonces:

yp (x) = e3x cos(x) 30e−3x sen(x) + 10e−3x cos(x) + e3x sen(x) 10e−3x sen(x) − 30e−3x cos(x)
   

= 30sen(x)cos(x) + 10cos2 (x) + 10sen2 (x) − 30sen(x)cos(x)


= 10(cos2 (x) + sen2 (x)) = 10

∴ y(x) = e3x (C1 cos(x) + C2 sen(x)) + 10

Del PVI:

y(0) =C1 cos(0) + C2 sen(0) + 10 = 10


⇒ C1 , C2 = 0
∴ y(x) = 10
y(0) =C1 cos(0) + C2 sen(0) + 10 = 5
⇒ C1 + 10 = 5
C1 = −5, C2 = 0
∴ y(x) = −5e3x cos(x) + 10

21
2
d) ddt2x − 4 dx
dt
+ 4x = et + e2t + 1

00 0
x − 4x + 4x = et + e2t + 1
Resolviendo la parte homogénea
00 0
x − 4x + 4x = 0
⇒ λ2 − 4λ2 + 4 = 0
(λ − 2)2 = 0
λ1,2 = 2

∴ xH (t) = C1 e2t + C2 te2t

x1 = e2t , x2 = te2t

2t 2t

e te
w = 2t 2t

2t
2e 2e + te
2t 2t

e te
= 2t 2t
2e e (2t + 1)

= e4t (2t + 1) − 2te4t


= e4t (2t + 1 − 2t) = e4t


0 te2t
= −te3t − te4t − te2t
w1 = t

e + e2t + 1 e2t (2t + 1)
2t
e 0
w2 = 2t t 2t
= e3t + e4t + e2t
2e e + e + 1

Entonces:
−te3t − te4t − te2t
C10 (t) = 4t
= −te−t − t − te−2t
e
e3t + e4t + e2t
C20 (t) = = e−t + 1 + e−2t
e4t

Z Z Z Z
−t −2t −t
−te−2t dt

C1 (t) = −te − t − te dt = −te dt − tdt +

22
Integrando por partes

u = t du = 1
dv = −e−t v = e−t
m = t dm = 1
1
dn = −e−2t n = e−2t
2

Z Z Z
−te dt − tdt + −te−2t dt
−t

t2 1 −2t 1
Z Z
−t −t
= te − e dt − + te − e−2t dt
2 2 2
t2 1 1
= te−t + e−t − + te−2t + e−2t
2 2 4

Z Z Z Z
−t −2t −t
e−2t dt

C2 (t) = e +1+e dt = e dt + dt +
1
= −e−t + t − e−2t
2
Entonces:

t2 1 −2t 1 −2t
   
2t −t −t 2t −t 1 −2t
xp (t) = e te + e − + te + e + te −e + t − e
2 2 4 2
2 2t
te 1 1 1
= tet + et − + t + − tet + t2 e2t − t
2 2 4 2
2 2t
te 1
= et + +
2 4

t2 e2t 1
∴ y(x) = C1 e2t + C2 te2t + et + +
2 4

e) ẍ + x = sent − cos(2t) Tomando la ecuación auxiliar para encontrar una


solución particular
ẍ + x = 0
Usando la ecuación caracterı́stica.

λ2 + 1 = 0 ⇒ λ1,2 = ±i

⇒ xH (t) = c1 cos(t) + c2 sen(t)


x1 (t) = cos(t) , x2 (t) = sen(t)

23
x01 (t) = −sen(t) , x02 (t) = cos(t)

cos(t) sen(t)
W = =1
−sen(t) cos(t)

0 sen(t)
W1 =
= cos(2t)sen(t) − sen2 (t)
sen(t) − cos(2t) cos(t)
1
Simplificando W1 usando las identidades sen(α)cos(β) = sen(α − β) +
2
2 1
sen(α + β) y sen (α) = (1 − cos(2α))
2
1
W1 = (sen(3t) + cos(2t) − sen(t) − 1)
2

cos(t) 0
W2 = = sen(t)cos(t) − cos(2t)cos(t)
−sen(t) sen(t) − cos(2t)
1
Simplificando W2 usando la identidad cos(α)cos(β) = cos(α−β)+cos(α+β)
2
1
W2 = sen(t)cos(t) − (cos(3t) + cos(t))
2
Z Z Z Z 
1
C1 = sen(3t)dt + cos(2t)dt − sen(t)dt − dt
2
 
1 1 1
C1 = − cos(3t) + sen(2t) + cos(t) − t
2 3 2
1 1 1 1
C1 = − cos(3t) + sen(2t) + cos(t) − t
6 4 2 2
Z Z Z 
1
C2 = sen(t)cos(t)dt − cos(3t)dt + cos(t)dt
2
| {z }
I0

Integrando I0 por partes

u = cos(t) ⇒ −du = sen(t)dt


Z
1 1
I0 = − udu = − u2 = − cos2 (t)
2 2
Z Z 
1
C2 = I0 − cos(3t)dt + cos(t)dt
2
 
1 2 1 1
C2 = − cos (t) − sen(3t) + sen(t)
2 2 3
1 1 1
C2 = − cos2 (t) − sen(3t) − sen(t)
2 6 2

24
xp (x) = C1 x1 + C2 x2
 
1 1 1 1
xp (x) = − cos(3t) + sen(2t) + cos(t) − t cos(t)
6 4 2 2
 
1 1 1
+ − cos2 (t) − sen(3t) − sen(t) sen(t)
2 6 2
1 1 1 1
xp (x) = − cos(3t)cos(t) + sen(2t)cos(t) + cos2 (t) − tcos(t)
6 4 2 2
1 1 1
− cos2 (t)sen(t) − sen(3t)sen(t) − sen2 (t)
2 6 2

   
1 2 2  1  1
xp = cos (t) − sen (t) − cos(3t)cos(t) − sen(3t)sen(t) − tcos(t)
2 | {z } 6 | {z } 2
cos(2t) cos(2t)

 
1 1 1
− (1 + cos(2t)) sen(t) + sen(2t)cos(t)
2 2 4
 
1 1 1 1 1
cos(2t) − cos(2t) − tcos(t) − sen(t) + sen(2t)cos(t) − cos(2t)sen(t)

2 6 2 4 4 | {z }
sen(t)

1 1 1 1
cos(2t) − tcos(t) − sen(t) + sen(t)
3 2 4 4
1 1
∴ xp (t) = cos(2t) − tcos(t) Es la solución particular de la ED
3 2

4. Ejercicio 4
Usando el método de Euler, resuelva la Ecuación de Cauchy-Euler:

a)x2 y 00 − xy 0 + y = ln(x) Haciendo uso de la ecuación auxiliar

x2 y 00 − xy 0 + y = 0

y = xm , y 0 = mxm−1 , y 00 = m(m − 1)xm−2


x2 m(m − 1)xm−2 − xmxm−1 + xm = 0
xm (m2 − m) − xm (m) + xm = 0
xm (m2 − m − m + 1) = 0
xm (m2 − 2m + 1) = 0
Ya que xm 6= 0 ∀x ∈ (0, +∞)
m2 − 2m + 1 = 0 (m − 1)2 = 0

25
⇒ m1,2 = 1
yH (x) = C1 x + C2 xln(x)
Reescribiendo a) en su forma estándar para hacer uso del método de variación
de parámetros y encontrar una solución particular

1 1 ln(x)
y 00 − y 0 + 2 y =
x x x2
y1 = x , y2 = xln(x)

x xln(x)
W = =x
1 ln(x) + 1
2

= − ln (x)
0 xln(x)
W1 = lnx
x2
ln(x) + 1 x

x 0 ln(x)
W2 = ln(x) =
1 x2 x
W1 ln2 (x)
C10 = =−
W x2
ln2 (x)
Z
C1 = − dx
x2
| {z }
I0

Integrando por partes:

2ln(x)
u = ln2 (x) ⇒ du = dx
x
dx 1
dv = 2
⇒v=−
x x
2 Z
ln (x) ln(x)
I0 = − +2 dx
x x2
| {z }
I1

Integrando I1 por partes


1
u = ln(x) ⇒ du = dx
x
dx 1
dv = 2
⇒v=−
x x
Z
ln(x) 1 ln(x) 1
I1 = − + 2
dx = − −
x x x x
ln2 (x) ln2 (x)
 
ln(x) 1
I0 = − + 2I1 = − −2 +
x x x x

26
ln2 (x) 2ln(x) 2
C1 = −I0 = + +
x x x
2
ln x + 2ln(x) + 2
C1 =
x
W2 ln(x)
C20 = =
W x2
Z
ln(x)
C2 = dx
x2
Integrando por partes:
dx)
u = ln(x) ⇒ du =
x
dx 1
dv = ⇒ v = −
x2 x
Z Z
ln(x) ln(x) 1 ln(x) 1
dx = − − − dx = − −
x2 x x2 x x
ln(x) + 1
C2 = −
x
yP (x) = c1 (x)y1 (x) + c2 (x)y2 (x)
 2   
ln x + 2ln(x) + 2 ln(x) + 1
yP (x) = x− xln(x)
x x
yP (x) = ln2 (x) + 2ln(x) + 2 + c2 (x)y2 (x) − ln2 (x) − ln(x)
yP (x) = ln(x) + 2
∴ y(x) = C1 x + C2 xln(x) + ln(x) + 2
Es la familia de soluciones para la ED
b)x2 y 00 − 4xy 0 + 6y = 2

x2 y 00 − 4xy 0 + 6y = 0
y = xm , y 0 = mxm−1 , y 00 = m(m − 1)xm−2
x2 m(m − 1)xm−2 − 4xmxm−1 + 6xm = 0
xm (m2 − m) − 4xm m + 6xm = 0
xm (m2 − m − 4m + 6) = 0
xm (m2 − 5m + 6) = 0
Ya que
xm 6= 0 ∀x ∈ (0, +∞)
m2 − 5m + 6 = 0
m1 = 2 , m2 = 3

27
Yh (x) = C1 x2 + C2 x3
Usando Variación de parámetros para encontrar una solución particular

y1 = x2 , y2 = x3

Reescribiendo la ecuación en su forma estándar:


4 6 2
y 00 − y 0 + 2 y = 2
x x x
2
x x3 4
W = x
2x 3x2

0 x3
W1 = 2 2 = −2x

x 2 3x
2
x 0
W2 = 2 = 2

2x x2
W1 2x 2
C10 = =− 4 =− 3
W x x
Z
dx
C1 = −2
x3
 
1
C1 = −2 − 2
2x
1
C1 = 2
x
W 2 2
C20 = = 4
W x
Z
dx
C2 = 2
x4
 
1
C2 = 2 − 3
3x
2
C2 = − 3
3x
   
1 2 2
yp (x) = x − x3
x2 3x3
3 2
yp (x) = −
3 3
1
yp (x) =
3
1
y(x) = C1 x2 + C2 x3 +
3
Es la familia de soluciones para la ED

28
5. Ejercicio 5
A)Un punto material se hunde lentamente en un lı́quido. Hallar la ley de
movimiento, considerando que durante una inmersión lenta la resistencia del
lı́quido es proporcional a la velocidad de inmersión.

Figura 1: Objeto puntual hundiéndose en un liquido

Considerando a Rl como la resistencia del lı́quido y Vi como la velocidad


de inmersión, en base al enunciado ”durante una inmersión lenta la resistencia
del lı́quido es proporcional a la velocidad de inmersión”tenemos la relación
proporcional:
Rl ∝ Vi
Si agregamos una constante de proporcionalidad ”k” la relación se vuelve di-
rectamente proporcional:
Rl = kVi
Haciendo un análisis de fuerzas considerando el sentido positivo hacia abajo
tenemos:

Figura 2: Fuerzas externas actuando sobre el material

29
X
Fy = mg − Rl = ma
mg − kVi = ma
Ya que el análisis es para cualquier instante dentro de la trayectoria de
inmersión se considerará a la velocidad y a la aceleración como instantáneas:
dx d2 x
mg − k =m 2
dt dt
d2 x dx
m 2 +k = mg
dt dt
mẍ(t) + k ẋ(t) = mg
Reescribiendo la ecuación en su forma estándar:
k
ẍ(t) + ẍ = g
m
Resolviendo la parte homogénea:
k
ẍ(t) + ẋ(t) = 0
m
k
λ2 + λ = 0
m
 
k
λ λ+ =0
m
λ1 = 0
k
λ2 = −
m
k
− t
xH (t) = C1 + C2 e m
Resolviendo la parte particular mediante el método de coeficientes constantes
xP (t) = At
ẋP (t) = A
ẍP (t) = 0
Sustituyendo en la ED:
k
(0) + A=g
m
mg
A=
k
mg
xP (t) = t
k
k

t mg
∴ x(t) = C1 + C2 e m + t
k

30
B)Encuentre la carga en el capacitor de un circuito en serie LRC en t =
0.01s cuando L = 0.05h, R = 2ohms, C = 0.01f, el voltaje de entrada E(t) =
tV , q(0) = 5C, e i(0) = 0A. Determine la primera vez en que la carga del
capacitor es igual a cero.

Figura 3: Circuito LRC

ε(t) = VL + VR + VC
dq
i=
dt
dq d2 q
VL = L ⇒ VL = L 2
dt dt
dq
VR = Ri ⇒ VR = R
dt
q
VC =
c
Sustituyendo los valores dados:

1 d2 q dq
2
+ 2 + 100q = t
20 dt dt
Reescribiendo la ecuación en su forma estándar

q 00 + 40q 0 + 2000q = 20t

Resolviendo la parte homogénea:

q 00 + 40q 0 + 2000q = 0

λ2 + 40λ + 2000 = 0

−40 ± 2 1600 − 8000
λ=
2
λ1,2 = −20 ± 40i

31
qH (t) = e−20t (C1 cos(40t) + C2 sen(40t))
Planteando la solución particular:

qP (t) = At + B

qP0 (t) = A
qP00 (t) = 0
Entonces:
40A + 2000(At + B) = 20t
2000A = 20 ⇒ A = 0,01
40A + 2000B = 0 ⇒ B = 0,0002
Sustituyendo los valores:

qP (t) = 0,01t + 0,0002

∴ q(t) = e−20t (C1 cos(40t) + C2 sen(40t)) + 0,01t + 0,0002


Usando
 el PVI
q(0) = 5C
i(0) = 0A
q(0) = e0 (C1 cos(0) + C2 sen(0)) 0,0002 = 5
Ya que ,0002 es una cantidad muy pequeña es despreciable.

C1 = 5

q(t) = e−20t (5cos(40t) + C2 sen(40t)) + 0,01t + 0,0002


i(t) = −20e−20t (5cos(40t) + C2 sen(40t))+e−20t (−250sen(40t) + 40C2 cos(40t))+,01V
Ya que el voltaje mide la variación de potencial, sin amperaje no hay tal por
lo cual ,01v es despreciable para el análisis en i(0)

i(0) = −100 + 40C2 = 0


5
C2 =
2
 
−20t 5
⇒ q(t) = e 5cos(40t) + sen(40t) + 0,01t + 0,0002
2
Dado que la corriente eléctrica no varia a lo largo del circuito la carga en el
capacitor será la misma que en todo el circuito

q(,01) ' 4,108

32
La carga del capacitor en t=.01s =4.108C
Mediante un sofware y un método numérico se estimo que la primera ves que
la carga del capacitor sea 0 será en t ' 0,509s

Figura 4: Gráfica del comportamiento de la carga en el primer cuarto de se-


gundo

C) Calcule la carga en el capacitor y la corriente en un circuito LC cuando


L = 0.1h, C = 0.1f, el voltaje de entrada E(t) = 100sen(2t)V , q(0) = 0C e
i(0) = 0A.

Figura 5: Circuito LC

ε(t) = VL + VC
dq
i=
dt

33
dq d2 q
VL = L ⇒ VL = L 2
dt dt
q
VC =
c
Sustituyendo los valores dados:
d2 q
0,1 + 10q = 100sen(2t)
dt2
q 00 + 100q = 1000sen(2t)
Resolviendo la parte homogénea :
q 00 + 100q = 0
λ2 + 100 = 0
λ1,2 = ±10i
qH (t) = C1 cos(10t) + C2 sen(10t)
qP (t) = Acos(2t) + Bsen(2t)
qP0 (t) = −2Asen(2t) + 2Bcos(2t)
qP00 (t) = −4Acos(2t) − 4Bsen(2t)
Entonces:
−4Acos(2t) − 4Bsen(2t) + 100(Acos(2t) + Bsen(2t)) = 1000sen(2t)
(−4A + 100A)cos(2t) + (−4B + 100B)sen(2t) = 1000sen(2t)
96A = 0 ⇒ A = 0
125
96B = 1000 ⇒ B =
12
Sustituyendo los valores de A y B:
125
qP (t) = sen(2t)
12
125
∴ q(t) = C1 cos(10t) + C2 sen(10t) + sen(2t)
12
Tomando en cuenta la consideración inicial q(0) = 0C
q(0) = C1 = 0
125
q 0 (t) = −10C1 sen(10t)10C2 cos(10t) + cos(2t)
6
Y como q 0 (0) = i(0) = 0A
125 25
i(0) = 10C2 + = 0 ⇒ C2 = −
6 12
Sustituyendo los valores de C1 y C2 en q(t) la solución del problema es:
25 125
q(t) = sen(10t) − sen(2t)
12 12

34

You might also like