Universidad de Antioquia

Facultad de Ingeniería

Tarea #2

Paula Valentina Pinzón Gómez CC. 1007423535

Bioingenieria

Adrian Montoya Lince

Ecuaciones Diferenciales

Medellín
2020
a) Determine la solución del PVI sabiendo que y1 = ex
xy’’ - (2x + 1)y’ + (x + 1)y = 𝑒 𝑥 , con y(1) = 0; y’(1) = 0
Normalizando la ED dividiendo entre “x”:
1 1 𝑒𝑥
y’’ - (2 + 𝑥)y’ + (1 + 𝑥)y = 𝑥

Hallamos el factor integrante

H
a
1
l =
∅ 𝑒 2𝑥 ∗ 𝑒 ∫ −(2+𝑥)𝑑𝑥
l
∅
a= 𝑒 2𝑥 ∗ 𝑒 −2𝑥−ln(𝑥)
m
∅= 𝑒 2𝑥 ∗ 𝑒 −2𝑥 ∗ 𝑒 −ln(𝑥)
o
s =
∅ 𝑥 −1

Luego
e hallamos y2 mediante 𝑦2 = 𝑦1 ∫ ∅−1𝑑𝑥:
l
𝑦2 = 𝑒 𝑥 ∫(𝑥 −1 )−1 𝑑𝑥
f
a 1
𝑦2 = 𝑒 𝑥 ∗ ( 𝑥 2 )
c 2
t 1 2 𝑥
𝑦
o2 = 𝑥 𝑒
2
r

i
Por lo tanto, la solución homogénea será de a forma:
n
t𝑦 = 𝐶1 𝑒 𝑥 + 𝐶2 (1 𝑥 2 𝑒 𝑥 )
2
e
g
Hallamos la solución particular:
r
𝑟(𝑥)
a
𝑦𝑝 = 𝑦1 ∗ ∫ ∅−1 [∫ ∅ ∗ 𝑑𝑥]𝑑𝑥
n 𝑦1
t 𝑒 𝑥 ∗ 𝑥 −1
e𝑝 = 𝑒 𝑥 ∗ ∫(𝑥 −1 )−1 [∫ 𝑥 −1
𝑦 𝑑𝑥]𝑑𝑥
𝑒𝑥
,

d
o
n
d
e
𝑦𝑝 = 𝑒 −𝑥 ∗ ∫(𝑥 −1 )−1 [∫ 𝑥 −2 𝑑𝑥]𝑑𝑥

−1
𝑦𝑝 = 𝑒 −𝑥 ∗ ∫(𝑥 −1 )−1 [ ] 𝑑𝑥
𝑥
−1
𝑦𝑝 = 𝑒 −𝑥 ∗ ∫ 𝑥 𝑑𝑥
𝑥
𝑦𝑝 = 𝑒 𝑥 ∗ −𝑥

La solución general será:

𝑦𝑠𝑠 = 𝐶1 𝑒 𝑥 + 𝐶2 (𝑥 2 𝑒 𝑥 ) + 𝑦𝑝

𝑦𝑠𝑠 = 𝐶1 𝑒 𝑥 + 𝐶2 (𝑥 2 𝑒 𝑥 ) − 𝑥𝑒 𝑥

Mediante las condiciones iniciales y(1) = 0; y’(1) = 0 encontramos 𝐶1 y 𝐶2 :

𝑦𝑠𝑠 = 𝐶1 𝑒 𝑥 + 𝐶2 (𝑥 2 𝑒 𝑥 ) − 𝑥𝑒 𝑥

𝑦′𝑠𝑠 = 𝐶1 𝑒 𝑥 + 𝐶2 (𝑒 𝑥 𝑥 2 + 2𝑒 𝑥 𝑥) − 𝑥𝑒 𝑥 − 𝑒 𝑥
Entonces, reemplazando e igualando:

0 = 𝐶1 𝑒1 + 𝐶2 (12 𝑒 𝑥 ) − 1𝑒 1

0 = 𝐶1 𝑒1 + 𝐶2 (𝑒 1 12 + 2𝑒 1 1) − 1𝑒 1 − 𝑒 1
Nuestras ecuaciones quedarían de la siguiente forma:
0 = 𝐶1 𝑒 + 𝐶2 𝑒 − 𝑒
0 = 𝐶1 𝑒 + 𝐶2 (𝑒 + 2𝑒) − 2𝑒
Resolviendo el sistema de ecuaciones, encontramos que:
1 1
𝐶1 = 2 𝐶2 = 2

Así, la solución al P.V.I:

1 𝑥 1 2 𝑥
𝑦𝑠𝑠 = 𝑒 + (𝑥 𝑒 ) − 𝑥𝑒 𝑥
2 2
Y la gráfica será de la siguiente manera:
b) Con ayuda de software de cálculo numérico resuelva el PVI y grafique el intervalo
indicado. Compare con la solución obtenida en el numeral anterior.

Podemos usar Matlab para buscar la solución directamente introduciendo el P.V.I:

Su grafica será:

ANALISIS: Podemos ver que las gráficas nos dieron idénticas, pues la obtención
del PVI manual y por software nos dieron resultados exactos.
a. Escriba las EDO y determine los CFS para los desplazamientos x1 y x2. ¿Qué tipos
de movimientos experimentaran las masas?

𝒅𝟐 𝒙𝟏
𝒎𝟏 = −𝒌𝟏 𝒙𝟏 + 𝒌𝟑 (𝒙𝟐 − 𝒙𝟏 ) + 𝒇(𝒕)
𝒅𝒕𝟐
𝒅𝟐 𝒙𝟐
{ 𝒎𝟐 = −𝒌𝟐 𝒙𝟐 − 𝒌𝟑 (𝒙𝟐 − 𝒙𝟏 )
𝒅𝒕𝟐

Resolvemos

𝒅𝟐 𝒙𝟏
𝒎𝟏 = −𝒌𝟏 𝒙𝟏 + 𝒌𝟑 𝒙𝟐 − 𝒌𝟑 𝒙𝟏 + 𝒇(𝒕)
𝒅𝒕𝟐
𝒅𝟐 𝒙𝟐
{ 𝒎𝟐 = −𝒌𝟐 𝒙𝟐 − 𝒌𝟑 𝒙𝟐 + 𝒌𝟑 𝒙𝟏
𝒅𝒕𝟐

Normalizamos y pasamos los respectivos términos al otro lado del igual

𝒅𝟐 𝒙𝟏 𝒌𝟏 𝒌𝟑 𝒌𝟑 𝒇(𝒕)
𝟐
+ 𝒙𝟏 − 𝒙𝟐 + 𝒙𝟏 =
𝒅𝒕 𝒎𝟏 𝒎𝟏 𝒎𝟏 𝒎𝟏
𝒅𝟐 𝒙𝟐 𝒌𝟐 𝒌𝟑 𝒌𝟑
+ 𝒙 + 𝒙 − 𝒙 =𝟎
{ 𝒅𝒕𝟐 𝒎𝟏 𝟐 𝒎𝟐 𝟐 𝒎𝟐 𝟏
 𝒅𝟐 𝒙𝟏 𝒌𝟏 𝒌𝟑 𝒌𝟑 𝒇(𝒕)
𝟐
+( + )𝒙𝟏 − 𝒙𝟐 =
𝒅𝒕 𝒎𝟏 𝒎𝟏 𝒎𝟏 𝒎𝟏
𝟐
𝒅 𝒙𝟐 𝒌𝟐 𝒌𝟑 𝒌𝟑
+( + )𝒙𝟐 − 𝒙 =𝟎
{ 𝒅𝒕
𝟐 𝒎𝟐 𝒎𝟐 𝒎𝟐 𝟏

𝒌𝟏 𝒌𝟑 𝒌𝟑 𝒇(𝒕)
(𝑫𝟐 + ( + ))𝒙𝟏 − 𝒙𝟐 =
𝒎𝟏 𝒎𝟏 𝒎𝟏 𝒎𝟏
𝒌𝟑 𝒌𝟐 𝒌𝟑
− 𝒙 + (𝑫𝟐 + ( + ))𝒙 = 𝟎
{ 𝒎𝟐 𝟏 𝒎𝟐 𝒎𝟐 𝟐

Y se procede a montar el sistema matricial:

𝑘1 𝑘3 𝑘3
𝐷2 + ( + ) − 𝑓(𝑡)
𝑚1 𝑚1 𝑚1 𝑥1
[𝑥 ] = [ 𝑚1 ]
𝑘3 𝑘2 𝑘3 2
− 𝐷2 + ( + ) 0
[ 𝑚2 𝑚2 𝑚2 ]

Reemplazamos los valores conocidos de 𝑘1 = 1.044, 𝑘2 = 2.395, 𝑘3 = 8.004, 𝑚1 =

10.552 𝑦 𝑚2 = 1.623

𝑓(𝑡)
2
[𝐷 + 0.857 −0.758 ] [𝑥1 ] = [
−4.956 𝐷 2 + 6.407 𝑥2 10.552]
0

Mediante regla de Cramer, hallamos las determinantes:

2
∆= det ([𝐷 + 0.857 −0.758 ]) = 𝐷 4 + 7.264𝐷2 + 1.734151
−4.956 𝐷 2 + 6.407

𝑓(𝑡) (𝐷 2 + 6.407)𝑓(𝑡)
∆𝑥1 = det ([10.552 −0.758 ]) =
10.552
0 𝐷 2 + 6.407
2 𝑓(𝑡)
∆𝑥2 = det ([𝐷 + 0.857
1.239𝑓(𝑡)
10.552]) = 2.638
−4.956 0

De esta manera, las ecuaciones resultan de la forma:

(𝐷 2 + 6.407)𝑓(𝑡)
(𝐷 4 + 7.264𝐷2 + 1.734151)𝑥1 =
10.552

1.239𝑓(𝑡)
(𝐷 4 + 7.264𝐷2 + 1.734151)𝑥2 =
2.638
Así, de las ecuaciones obtenemos el polinomio característico asociado al
movimiento:

𝐿(𝜆) = 𝜆4 + 7.264𝜆2 + 1.734151

Ahora para los CFS del polinomio

𝑃(𝜆) = 𝜆4 + 7.264𝜆2 + 1.734151

Las raíces de este polinomio serán:
𝜆1 = 0.4970𝑖 𝜆2 = −0,4970𝑖 𝜆3 = 2.648𝑖 𝜆4 = −2.648𝑖
𝐶𝐹𝑆 = {cos(0.4970𝑡) , sin(0.4970𝑡) , cos(2.648𝑡) , sin(2.648𝑡)}
𝑦ℎ = 𝑐1 cos(0.4970𝑡) + 𝑐2 sin(0.4970𝑡) + 𝑐3 cos(2.648𝑡) + 𝑐3 sin(2.648𝑡)

Por lo tanto, podemos afirmar que el movimiento es sub-amortiguado.

La solución homogénea es la misma para x1 y x2. Ahora, sabiendo que F0 = 2N

Hallamos Xss1
(𝐷 2 +6.407)𝑓(𝑡)
(𝐷 4 + 7.264𝐷2 + 1.734151)𝑥1 =
10.552

(𝐷 4 + 7.264𝐷2 + 1.734151)𝑥1 = 0.0947𝐷2 + 0.5730𝑓(𝑡)

0.0947𝐷 2 + 0.5730
𝑥1 = (2)
(𝐷 4 + 7.264𝐷2 + 1.734151)

0.0947(0)2 + 0.5730
𝑥1 = (2)
((0)4 + 7.264(0)2 + 1.734151)

𝑥1 = 0.6608

Por lo tanto, la solución general para x1 será:

 𝑥 = 𝑐1 cos(0.4970𝑡) + 𝑐2 sin(0.4970𝑡) + 𝑐3 cos(2.648𝑡) + 𝑐3 sin(2.648𝑡) + 0.6608

Se sabe que el sistema parte del reposo y se estudia a partir de t>0, las condiciones
iniciales son:
𝑓 2
= = 0.1895
𝑚1 10,552

Ahora se calcula la solución al PVI para 𝑥1 :

𝑥 = −0.66cos(0.4970𝑡) − 3.6e − 3 cos(2.648𝑡) + 0.6608

La solución homogénea es la misma para x1 y x2. Ahora, sabiendo que F0 = 2N

Hallamos Xss2

1.239𝑓(𝑡)
(𝐷 4 + 7.264𝐷 2 + 1.734151)𝑥2 =
2.638

(𝐷 4 + 7.264𝐷2 + 1.734151)𝑥2 = 0.4696𝑓(𝑡)

0.4696
𝑥2 = (2)
(𝐷 4 + 7.264𝐷2 + 1.734151)

0.4696
𝑥2 = (2)
((0)4 + 7.264(0)2 + 1.734151)

𝑥2 = 0.541590

Por lo tanto, la solución general para x2 será:

𝑥2 = 𝑐1 cos(0.4970𝑡) + 𝑐2 sin(0.4970𝑡) + 𝑐3 cos(2.648𝑡) + 𝑐3 sin(2.648𝑡) + 0.541590

Se sabe que el sistema parte del reposo y se estudia a partir de t>0, además se
sabe que en este cuerpo no se presenta fuerza neta, para x2 las condiciones iniciales
son:
 Ahora se calcula la solución al PVI para 𝑥1 :
𝑥 = −0.56cos(0.4970𝑡) + 0.02 cos(2.648𝑡) + 0.6608

GRAFICAS
1) Graficas de las soluciones obtenidas en el tiempo para x1 , x2
2) Diagrama de fase para x1

La solución homogénea es la misma para x1 y x2 ahora, sabiendo que

F0 = sin(0.517t) N
Hallamos Xss2

1.239𝑓(𝑡)
(𝐷 4 + 7.264𝐷2 + 1.734151)𝑥2 =
2.638
(𝐷 4 + 7.264𝐷2 + 1.734151)𝑥2 = 0.4696𝑓(𝑡)

0.4696
𝑥2 = (sin(0.517𝑡))
(𝐷 4 + 7.264𝐷2 + 1.734151)

0.4696
𝑥2 = (sin(0.517𝑡))
((0.517)4 + 7.264(0.517)2 + 1.734151)

𝑥2 = 0.125320 (sin(0.517𝑡))

Por lo tanto, la solución general para x2 será:

𝑥2 = 𝑐1 cos(0.4970𝑡) + 𝑐2 sin(0.4970𝑡) + 𝑐3 cos(2.648𝑡) + 𝑐3 sin(2.648𝑡)
+ 0.125320 (sin(0.517𝑡))
Se sabe que el sistema parte del reposo y se estudia a partir de t>0, además se
sabe que en este cuerpo no se presenta fuerza neta, para x2 las condiciones iniciales
son:
 𝑥 = −0.1299 sin(0.4970𝑡) − 0.000076 sin(2.648𝑡) + 0.125320 (sin(0.517𝑡))
GRAFICAS
1) Diagrama de fase:

La solución homogénea es la misma para x1 y x2 ahora, sabiendo que

F0 = sin(2.531t) N
 Hallamos Xss1

(𝐷 2 + 6.407)𝑓(𝑡)
(𝐷 4 + 7.264𝐷2 + 1.734151)𝑥1 =
10.552

(𝐷 4 + 7.264𝐷2 + 1.734151)𝑥1 = 0.0947𝐷2 + 0.5730𝑓(𝑡)

0.0947𝐷 2 + 0.5730
𝑥1 = (sin(2.531t) )
(𝐷 4 + 7.264𝐷2 + 1.734151)

0.0947(2.531)2 + 0.5730
𝑥1 = (sin(2.531t) )
((2.531)4 + 7.264(2.531)2 + 1.734151)

𝑥1 = 0.01320 ∗ sin(2.531t)

Por lo tanto, la solución general para x1 será:

𝑥 = 𝑐1 cos(0.4970𝑡) + 𝑐2 sin(0.4970𝑡) + 𝑐3 cos(2.648𝑡) + 𝑐3 sin(2.648𝑡) + 0.01320

∗ sin(2.531t)

Se sabe que el sistema parte del reposo y se estudia a partir de t>0, las condiciones
iniciales son:
𝑓 sin(2.531t)
= = 0.09476876 ∗ sin(2.531t)
𝑚1 10,552

ANALISIS
Tenido en cuenta de que lo único que cambia respecto a los ejercicios anteriormente
desarrollados es nuestro F0 = sin(2.531t), se puede concluir que lo único que cambia es la
solución particular, los otros aspectos procedimentales se mantienen.