UNAD Ecuaciones Diferenciales Actividad Colaborativa 2
UNAD Ecuaciones Diferenciales Actividad Colaborativa 2
UNAD Ecuaciones Diferenciales Actividad Colaborativa 2
UNAD
Ecuaciones Diferenciales
Grupo: 100412_263
Unidad 2
TUTOR:
Yenifer Elizabeth Galindo
ESTUDIANTE:
Erika Rocha Código 1071868434
Fabián Ardila Código 1095922175
Johan Rueda Código 1073712808
Diego Julián Andrade 1073506346
Gustavo Sánchez
Bogotá
Octubre de 2017
Introducción
𝑑 2 𝑦(𝑥) 𝑑𝑦(𝑥)
= 2
=2 − 3𝑦(𝑥)
𝑑𝑥 𝑑𝑥
𝑑2 ℷ𝑥 𝑑
2
(𝑒 ) − 2 (𝑒 ℷ𝑥 ) + 3𝑒 ℷ𝑥 = 0
𝑑𝑥 𝑑𝑥
Después de hacer esta sustitución, debemos hacer solo la otra parte para
𝑑2
despejar y desarrollar el ejercicio, en esta ocasión vamos a sustituir 𝑑𝑥 2 (𝑒 ℷ𝑥 ) por
𝑑
ℷ2 𝑒 ℷ𝑥 y 𝑑𝑥 (𝑒 ℷ𝑥 ) por ℷ𝑒 ℷ𝑥
ℷ2 𝑒 ℷ𝑥 − 2ℷ𝑒 ℷ𝑥 + 3ℷ𝑥 = 0
(𝑒 ℷ𝑥 ) ℷ2 − 2ℷ + 3 = 0
Después de este proceso y ya que 𝑒 ℷ𝑥 ≠ 𝑑𝑒 0 para cualquier finito ℷ los ceros
deben venir del polinomio: ℷ2 − 2ℷ + 3 = 0 vamos a solucionar: ℷ
ℷ = 1 + 𝒾√2 𝑜 ℷ = 1 − 𝒾√2
Al sacar la raíz nos da como resultado:
𝑦1(𝑥) = 𝑐1 𝑒 (1+𝒾√2)𝑥
𝑦2(𝑥) = 𝑐2 𝑒 (1−𝒾√2)𝑥
Como debemos sacar la solución general que incluya las ecuaciones pasadas
vamos a hacer la suma: 𝑦(𝑥) = 𝑦1(𝑥) + 𝑦2(𝑥)
𝑐1 𝑒 (1+𝒾√2)𝑥 + 𝑐2 𝑒 (1−𝒾√2)𝑥
Al ver el resultado con las posibles respuestas dadas en el ejercicio no está
completo. Por esto vamos a aplicar la identidad de Euler:
𝑑𝑛 𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
𝑎𝑛 (𝑥) 𝑛 + 𝑎𝑛−1 (𝑥) 𝑛−1 + ⋯ 𝑎1 (𝑥) + 𝑎0 (𝑥)𝑦 = 𝑔(𝑥)
𝑑𝑥 𝑑𝑥 𝑑𝑥
cuya solución general se escribe como la suma de las soluciones de una ecuación
homogénea y una particular.
𝑦 = 𝑦𝑐 + 𝑦𝑝
𝑦𝑐 se determina haciendo 𝑔(𝑥) = 0 para convertir la ecuación a una homogénea con
coeficientes constantes. Esta es la llamada solución asociada 𝑦𝑐 y se encuentra una solución
particular de la ecuación no homogénea. Esta es la llamada solución particular 𝑦𝑝 . Dicha
solución depende de la forma de la función g(𝑥). De acuerdo a lo mencionado
anteriormente la solución de la ecuación diferencial no homogénea 4𝑦 ´´ + 36𝑦 = csc 3𝑥 es:
1 1 1
A. 𝑦 = 𝐶1 𝑐𝑜𝑠 3𝑥 + 𝐶2 𝑠𝑖𝑛 3𝑥 − 12 𝑥 𝑠𝑖𝑛 3𝑥 + 36 (𝑐𝑜𝑠 3𝑥) 36 𝑙𝑛|𝑠𝑒𝑐 3𝑥|
1 1 1
B. 𝑦 = 𝐶1 𝑐𝑜𝑠 3𝑥 + 𝐶2 𝑠𝑖𝑛 3𝑥 + 12 𝑥 𝑠𝑖𝑛 3𝑥 − 36 (𝑐𝑜𝑠 3𝑥) 36 𝑙𝑛|𝑠𝑒𝑐 3𝑥|
1 1 1
C. 𝑦 = 𝐶1 𝑐𝑜𝑠 3𝑥 + 𝐶2 𝑠𝑖𝑛 3𝑥 + 12 𝑥 𝑐𝑜𝑠 3𝑥 − 36 (𝑠𝑖𝑛 3𝑥) 36 𝑙𝑛|𝑠𝑒𝑐 3𝑥|
1 1 1
D. 𝑦 = 𝐶1 𝑐𝑜𝑠 3𝑥 + 𝐶2 𝑠𝑖𝑛 3𝑥 − 12 𝑥 𝑐𝑜𝑠 3𝑥 + 36 (𝑠𝑖𝑛 3𝑥) 36 𝑙𝑛|𝑠𝑒𝑐 3𝑥|
Procedimiento:
𝟒𝒚´´ + 𝟑𝟔𝒚 = 𝒄𝒔𝒄 𝟑𝒙 Ecuación Diferencial
𝟏 Simplificamos y usamos el
𝒚′′ + 𝟗𝒚) 𝐜𝐬𝐜 𝟑𝒙 divisor 4
𝟒
𝟏 La función complementaria
𝒇(𝒙) = 𝐜𝐬𝐜 𝟑𝒙 para la ecuación determinado
𝟒
por 𝑦𝑐 = 𝑐1 cos 3𝑥 + 𝑐2 sin 3𝑥;
usando para 𝑦1 = cos 3𝑥 y
𝑦2 = sin 3𝑥
𝐜𝐨𝐬 𝟑𝒙 𝐬𝐢𝐧 𝟑𝒙 De la expresión anterior
𝑾(𝐜𝐨𝐬 𝟑𝒙 , 𝐬𝐢𝐧 𝟑𝒙) = =𝟑
−𝟑 𝐬𝐢𝐧 𝒙 𝟑 𝐜𝐨𝐬 𝟑𝒙 determinamos lo siguiente:
𝟎 𝐬𝐢𝐧 𝟑𝒙 𝟏
𝑾𝟏 = 𝟏 =−
𝐜𝐬𝐜 𝟑𝒙 𝟑 𝐜𝐨𝐬 𝟑𝒙 𝟒
𝟒
𝐜𝐨𝐬 𝟑𝒙
𝟎 𝟏 𝐜𝐨𝐬 𝟑𝒙
𝑾𝟐 = 𝟏 =
−𝟑 𝐬𝐢𝐧 𝟑𝒙 𝐜𝐬𝐜 𝟑𝒙 𝟒 𝐬𝐢𝐧 𝟑𝒙
𝟒
𝑊1
Ahora para: 𝑢1 ′ = y 𝑢2 ′ =
𝟏 𝟏 𝐜𝐨𝐬 𝟑𝒙 𝑊
− 𝟏 𝟏 𝐜𝐨𝐬 𝟑𝒙 𝑊2
𝒖′𝟏 = 𝟒
= − 𝟏𝟐 y 𝒖′𝟐 = 𝟒 𝐬𝐢𝐧 𝟑𝒙
= 𝟏𝟐 𝒔𝒊𝒏 𝟑𝒙 obtenemos:
𝟑 𝟑 𝑊
𝟏 𝒙 Integramos 𝑢1 ′ y 𝑢2 ′
𝒖𝟏 = ∫
𝒅𝒙 =
𝟏𝟐 𝟏𝟐
𝟏 𝐜𝐨𝐬 𝟑𝒙
𝒖𝟐 = ∫
𝟏𝟐 𝐬𝐢𝐧 𝟑𝒙
𝒙 𝟏 Reemplazamos en la ecuación
𝒚𝒑 = − 𝐜𝐨𝐬 𝟑𝒙 + 𝐥𝐨𝐠|𝒔𝒊𝒏 𝟑𝒙| 𝒔𝒆𝒄 𝟑𝒙
𝟏𝟐 𝟑𝟔 𝑦𝑝 = 𝑢1 𝑦1 + 𝑢2 𝑦2
𝒚𝒄 = 𝒄𝟏 𝐜𝐨𝐬 𝟑𝒙 + 𝒄𝟐 𝐬𝐢𝐧 𝟑𝒙 Por lo tanto, obtenemos:
𝒙 𝟏
− 𝐜𝐨𝐬 𝟑𝒙 + 𝐥𝐨𝐠|𝒔𝒊𝒏 𝟑𝒙| 𝐬𝐞𝐜 𝟑𝒙
𝟏𝟐 𝟑𝟔
4. Una ecuación diferencial de segundo orden homogéneo tiene dos soluciones
independientes. Para el caso 2 al resolver la ecuación característica las soluciones deben ser
reales repetidas 𝑚1=𝑚2 y su solución general es de la forma: 𝑦 = 𝐶1 𝑒 𝑚1 𝑥 + 𝐶2 𝑥𝑒 𝑚2 𝑥 .
Teniendo en cuenta la información anterior la solución general de la ecuación diferencial
𝑦 ′′ − 14𝑦 ′ + 49𝑦 = 0 corresponde a:
a. 𝑦 = 𝐶1 𝑒 −7𝑥 + 𝐶2 𝑥𝑒 −7𝑥
b. 𝑦 = 𝐶1 𝑒 7𝑥 + 𝐶2 𝑒 2𝑥
c. 𝑦 = 𝐶1 𝑒 7𝑥 + 𝐶2 𝑥𝑒 7𝑥
d. 𝑦 = 𝐶1 𝑒 −7𝑥 + 𝐶2 𝑒 −2𝑥
Procedimiento:
𝒀’’ − 𝟏𝟒𝒚’ + 𝟒𝟗 = 𝟎 Ecuación diferencial
𝒅𝒚(𝒙) 𝒅𝟐 𝒚(𝒙) Tenemos los pasos intermedios
−𝟏𝟒 + + 𝟒𝟗𝒚(𝒙)
𝒅𝒙 𝒅𝒙𝟐
=𝟎
𝟐
𝒅 𝒅 𝝀𝒙 Asumiendo que la solución será proporcional a:
𝟐
(𝒆𝝀𝒙 ) − 𝟏𝟒 (𝒆 ) + 𝟒𝟗𝒆𝝀𝒙 𝑒 𝜆𝑥 para una constante 𝜆.
𝒅𝒙 𝒅𝒙
=𝟎 Sustituimos y(x) por 𝑒 𝜆𝑥 En la ecuación
diferencial.
𝝀𝟐 𝒆𝝀𝒙 − 𝟏𝟒𝝀𝒆𝝀𝒙 + 𝟒𝟗𝒆𝝀𝒙 = 𝟎 𝑑2 𝑑
Sustituimos 𝑑𝑥 2 (𝑒 𝜆𝑥 ) = 𝜆2 𝑒 𝜆𝑥 y 𝑑𝑥 (𝑒 𝜆𝑥 ) =
𝜆𝑒 𝜆𝑥
(𝝀𝟐 − 𝟏𝟒𝝀 + 𝟒𝟗)𝝀𝒆𝝀𝒙 Tenemos que: 𝜆𝑒 𝜆𝑥 es el factor común.
𝝀𝟐 − 𝟏𝟒𝝀 + 𝟒𝟗 = 𝟎 Dado que 𝑒 𝜆𝑥 ≠ 0 para cualquier número finito
𝜆, los ceros deben provenir del polinomio
(𝝀 − 𝟕)𝟐 = 𝟎 Obteniendo el factor:
𝝀=𝟕o𝝀=𝟕 Despejando 𝜆 podemos decir:
La multiplicidad de la raíz 𝜆 = 7 es 2 que da
𝑦1 (𝑥) = 𝑐1 𝑒 7𝑥 , 𝑦2 (𝑥) = 𝑐2 𝑒 7𝑥 como soluciones,
𝒚(𝒙) = 𝒚𝟏 (𝒙) + 𝒚𝟐 (𝒙) donde C1 y C2 son constantes arbitrarias.
=
𝒄𝟏 𝒆𝟕𝒙 + 𝒄𝟐 𝒙𝒆𝟕𝒙 La solución general es la suma de las soluciones
anteriores obteniendo:
𝑑𝑛 𝑦
5. Una ecuación diferencial de orden superior es de la forma 𝑎𝑛 (𝑥) 𝑑𝑥 𝑛 +
𝑑𝑛−1 𝑦 𝑑𝑦
𝑎𝑛−1 (𝑥) 𝑑𝑥 𝑛−1 + ⋯ 𝑎1 (𝑥) 𝑑𝑥 + 𝑎0 (𝑥)𝑦 = 𝑔(𝑥) y puede ser solucionada por
diferentes métodos. La ecuación diferencial: 𝑦 ´´ − 𝑦 ´ + 𝑦 = 2 sin 3𝑥, puede ser
solucionada por los siguientes métodos y tiene como solución general:
Solución:
2𝑠𝑒𝑛(3𝑥)
𝑑𝑦(𝑥) 𝑑 2 𝑦(𝑥)
+ + 𝑦(𝑥) = 0
𝑑𝑥 𝑑𝑥 2
𝑑2 𝑑
2 ℷ𝑥
− − 𝑒 ℷ𝑥 = 0
𝑑𝑥 𝑒 𝑑𝑥𝑒 ℷ𝑥
𝑑2
Al tener esta primera sustitución, vamos a anunciar que 𝑑𝑥 2 𝑒 ℷ𝑥 = ℷ2 𝑒 ℷ𝑥 y
𝑑
= ℷ𝑒 ℷ𝑥 :
𝑑𝑥𝑒 ℷ𝑥
ℷ2 𝑒 ℷ𝑥 − ℷ𝑒 ℷ𝑥 − 𝑒 ℷ𝑥 = 0
En este paso debemos sacar factor de 𝑒 ℷ𝑥
(𝑒 ℷ𝑥 ) − ℷ2 − ℷ − 1 = 0
1 𝒾√3 1 𝒾√3
ℷ=− + 𝑜 ℷ=− −
2 2 2 2
Ahora hallamos ℷ
1 𝒾√3
(− + )𝑥
2 2
𝑌1(𝑥) = 𝐶1 𝑒
1 𝒾√3
(− − )𝑥
2 2
𝑌2(𝑥) = 𝐶2 𝑒
Ahora realizamos la suma de y1 y y2
1 𝒾√3 1 𝒾√3
(− + )𝑥 (− − )𝑥
2 2 2 2
𝑦(𝑥) = 𝐶1 𝑒 + 𝐶2 𝑒
Asociamos 𝐶1 + 𝐶2 = 𝐶1 𝑦 𝐶1 − 𝐶2 = 𝐶2
√𝟑(𝒙) √𝟑(𝒙)
𝒄𝒐𝒔
𝒚(𝒙) = 𝑪𝟏 𝟐 + 𝑪 (𝒔𝒆𝒏 𝟐
𝝅 𝟐 𝝅
𝒆𝟐 𝒆𝟐
𝑑2 𝑦(𝑥) 𝑑𝑦(𝑥)
La Solución Particular la vamos a resolver de la forma: − =
𝑑𝑥 2 −𝑦(𝑥) 𝑑𝑥
−2𝑠𝑒𝑛(3𝑥) por el método indeterminado de coeficientes:
𝑑𝑦𝑝 (𝑥) 𝑑
= + 𝑎2 𝑠𝑒𝑛(3𝑥 =
𝑑𝑥 𝑑𝑥𝑎1 cos(3𝑥)
𝟏 √𝟑 √𝟑 𝟏𝟔 𝟔
𝒚 = 𝒆−𝟐𝒙 (𝑪𝟏 𝒄𝒐𝒔 𝒙 + 𝑪𝟐 𝒔𝒆𝒏 𝒙) + 𝐜𝐨𝐬 𝟑𝒙 + − 𝐬𝐞𝐧 𝟑𝒙
𝟐 𝟐 𝟕𝟑 𝟕𝟑
Solución
Y”+y=0 en c1 cos x + c2
G(x)=4x +10 sin x
𝑦𝑝 = 𝐴𝑥 + 𝐵 + 𝐶𝑥 cos 𝑥 + 𝐸 sin 𝑥
𝑦𝑝 = 𝐴𝑥 + 𝐵 + 𝐶𝑥 cos 𝑥 + 𝐸 sin 𝑥
𝑦"𝑝 + 𝑦𝑝 = 𝐴𝑥 + 𝐵 − 2𝐶 sin 𝑥 + 2𝐸 cos 𝑥 = 4𝑥 + 10 𝑠𝑖𝑛
Entonces:
Se obtiene
𝑦𝑝 = 4𝑥 − 5𝑥 𝑐𝑜𝑠𝑥
Ecuación general:
an Dn y+an−1Dn−1 y+…+a1Dy+a0 y = g ( x )
𝑑𝑘𝑦
Donde 𝐷𝑘 𝑦 = , 𝑘 = 0,1,2,3, … … , 𝑛 Cuando se cumple la ecuación anterior
𝑑𝑥 𝑘
también se escribe como L (y) =g (x), donde L denota el operador diferencial o polinomial,
lineal de n-ésimo orden
an Dn+an−1 Dn−1+…+a1D+a0
La notación de operador no sólo es una abreviatura útil, sino que en un nivel muy práctico la
aplicación de operadores diferenciales permite justificar las reglas para determinar la forma
de la solución particular yp. Ésta se deduce casi de manera automática una vez se encuentra
un operador diferencial lineal adecuado que anula a g (x). Por lo anterior de la ecuación
diferencial 𝑦 ′′ − 3𝑦 ′ = 8𝑒 3𝑥 + 4𝑠𝑒𝑛(𝑥), se puede afirmar que:
Solución:
𝑦 ′′ − 3𝑦 ′ = 8𝑒 3𝑥 + 4𝑠𝑒𝑛(𝑥)
Se define la solución particular yc (Solución homogénea)
𝑦 ′′ − 3𝑦 ′ = 𝑂
𝐷2 − 3𝐷 = 0
Factorizamos: 𝐷(𝐷 − 3) = 0
Entonces tenemos que: 𝑀1 = 0 𝑦 𝑀2 = 3
Solución particular yc : 𝑦𝑐 = 𝐶1 𝑒 0𝑥 + 𝐶2 𝑒 3𝑥 = 𝐶1 + 𝐶2 𝑒 3𝑥
indeterminados
𝑦𝑝 = 8𝑒 3𝑥 + 4 sen 𝑥
8𝑒 3𝑥 = Axe3𝑥
4 sen 𝑥 = 𝐵𝑠𝑖𝑛(𝑥) + 𝐶𝑠𝑖𝑛(𝑥)
Nuestra solución particular seria:
𝑦 ′′ − 3𝑦 ′ = 8𝑒 3𝑥 + 4𝑠𝑒𝑛(𝑥)
Teniendo en cuenta que las raíces de 𝑦 ′′ − 3𝑦 ′ corresponden a 𝐷2 − 3𝐷,
proseguimos a determinar el operador diferencial que anula según la tabla para las
funciones:
𝑎=0
𝛽=1
𝑛=1
Teniendo en cuenta lo anterior, podemos calcular el operador diferencial que anula,
así:
(𝐷 − 3)8𝑒 3𝑥 = 0
(𝐷2 + 1)4𝑠𝑒𝑛(𝑥) = 0
Factorizamos (𝑚 − 5)(𝑚 − 5) = 0
Tenemos que:
𝑚−5=0 𝑚−5=0
𝑚1 = 5 𝑚2 = 5
𝒚 = 𝑪𝟏 𝒆 + 𝑪𝟐 𝒙𝒆𝟓𝒙
𝟓𝒙
𝐹𝑔 = −𝑚𝑔
−70𝑘𝑔 ∗ 9.8𝑚
𝐹𝑔 =
𝑠2
𝐹𝑔 = −686.7 𝑁
Al tener esta información debemos sacar la aceleración derivada segunda de y y
despejamos la formula apropiada
𝐹𝑒 + 𝐹𝑔 = 𝑚𝑎
−350𝑦 − 686.7 = 75 ∗ 𝑦 𝑖
75𝑦" + 350𝑦 = −686.7
350𝑦 −686.7
𝑦" + =
75 75
𝑦" + 4.6𝑦 = −9.156
𝑘 2 + 4.6 = 0
𝑘 2 = 4.6
𝑘 = ±√4.6𝑖
𝑘 = ±2.6𝑖
Ahora damos la solución general a la ecuación diferencial homogénea
𝑦𝑝(𝑡) = 𝐴
0 + 4.6𝐴 = −9.156
−9.156
𝐴=
4.6
𝐴 = −1.962
𝑦(0) = −8
𝐶1 − 1.962 = −8
𝐶1 = −6.038
𝑦(0) = 30
2.1602469𝐶2 = 30
𝐶2 = 13.8873015
Un sistema vibratorio que consiste en una masa unida a un resorte como se muestra en la
figura
1 1
Se suelta desde el reposo a unidades debajo de la posición de equilibrio. La masa es de 𝐾𝑔
2 5
𝑁
y la constante elástica es 𝑘 = 2 𝑚. El movimiento es amortiguado (𝛽 = 1,2) y está siendo
𝜋
impulsado por una fuerza periódica externa (𝑇 = 2 𝑠), comenzando en 𝑡 = 0. Dicha fuerza
está definida como 𝑓(𝑡) = 5 cos 4𝑡. Para esta situación, procedemos a encontrar la ecuación
diferencial que describe el movimiento
En los sistemas físicos acelerados la sumatorio de fuerzas se expresa de acuerdo a la
formulación de la segunda ley de Newton: ∑ 𝐹 = 𝑚𝑎
De acuerdo al problema planteado se tiene un Movimiento forzado con amortiguamiento. En
concordancia con la ley anterior:
𝑑2𝑥 𝑑𝑥
𝑚 2 = −𝑘𝑥 − 𝛽 + 𝑓(𝑡)
𝑑𝑡 𝑑𝑡
𝑑2 𝑥 𝑑𝑥
Donde la aceleración y la velocidad están dadas por 𝑎 = 𝑑𝑡 2 y 𝑣 = 𝑑𝑡
Transponiendo términos en la ecuación:
𝑑2 𝑥 𝑑𝑥
𝑚 2 +𝛽 + 𝑘𝑥 = 𝑓(𝑡)
𝑑𝑡 𝑑𝑡
Y reemplazando los valores dados en esta se tiene:
1 𝑑2 𝑥 𝑑𝑥 1
+ 1,2 + 2𝑥 = 5 cos 4𝑡 𝑥(0) = 𝑥 ´ (0) = 0
5 𝑑𝑡 2 𝑑𝑡 2
𝑑2 𝑥 𝑑𝑥
Equivalente a: + 4 𝑑𝑡 + 5𝑥 = 25 cos 4𝑡
𝑑𝑡 2
−16𝐴 cos 4𝑡 − 16𝐵 sin 4𝑡 + 4(−4𝐴 sin 4𝑡 + 4𝐵 cos 4 𝑡) + 5(𝐴 cos 4𝑡 + 𝐵 sin 4𝑡)
= 25 cos 4𝑡
−16𝐴 cos 4𝑡 − 16𝐵 sin 4𝑡 + 6(−4𝐴 sin 4𝑡 + 4𝐵 cos 4 𝑡) + 10(𝐴 cos 4𝑡 + 𝐵 sin 4𝑡)
= 25 cos 4𝑡
−16𝐴 cos 4𝑡 − 16𝐵 sin 4𝑡 − 24𝐴 sin 4𝑡 + 24𝐵 cos 4 𝑡 + 10𝐴 cos 4𝑡 + 10𝐵 sin 4𝑡
= 25 cos 4𝑡
−6𝐴 cos 4𝑡 − 6𝐵 sin 4𝑡 − 24𝐴 sin 4𝑡 + 24𝐵 cos 4 𝑡 = 25 cos 4𝑡
24𝐵 − 6𝐴 = 25
−6𝐵 − 24𝐴 = 0
𝐵 = 4𝐴
24(4𝐴) − 6𝐴 = 25
96𝐴 − 6𝐴 = 25
25
𝐴=
90
5
𝐴=
18
5
𝐵 =4∗
18
10
𝐵=
9
Operando:
−16 𝐴cos 4𝑡 − 16𝐵 sin 4𝑡 − 16𝐴 sin 4𝑡 + 16𝐵 cos 4𝑡 + 5𝐴 cos 4𝑡 + 5𝐵 sin 4𝑡
= 25 cos 4𝑡
Reuniendo términos semejantes:
−11𝐴 cos 4𝑡 − 11𝐵 sin 4𝑡 − 16𝐴 sin 4𝑡 + 16𝐵 cos 4𝑡 = 25 cos 4𝑡
Factorizando:
(−11𝐴 + 16𝐵) cos 4𝑡 + (−16𝐴 − 11𝐵) sin 4𝑡 = 25 cos 4𝑡
El sistema de ecuaciones resultante:
−11𝐴 + 16𝐵 = 25
−16𝐴 − 11𝐵 = 0
25 50
Se cumple que: 𝐴=− 𝑦 𝐵=
102 51
Reescribiendo: 𝑦𝑝 = 𝐴 cos 4𝑡 + 𝐵 sin 4𝑡
25 50
𝑦𝑝 = − cos 4𝑡 + sin 4𝑡
102 51
5 10
𝑦𝑝 = cos 4𝑡 + sin 4𝑡
18 9
La solución sería: 𝑦 = 𝑦𝑐 + 𝑦𝑝
25 50
𝑦 = 𝑒 −2𝑡 (𝐶1 cos 𝑡 + 𝐶2 sin 𝑡) − cos 4𝑡 + sin 4𝑡
102 51
5 10
𝑦 = 𝐶1 𝑒 −3𝑡 cos 𝑡 + 𝐶2 𝑒 −3𝑡 sin 𝑡 + cos 4𝑡 + sin 4𝑡
18 9
Haciendo 𝑡 = 0
1
CONDICION INICIAL 𝑦(0) = 2
1 5 10
= 𝐶1 𝑒 −3(0) cos 𝑡 + 𝐶2 𝑒 −3(0) sin 0 + cos 4(0) + sin 4(0)
2 18 9
1 5
= 𝐶1 +
2 18
5 1
𝐶1 = − +
18 2
2
𝐶1 =
9
25 50
𝑦(0) = 𝑒 −2(0) [𝐶1 cos(0) + 𝐶2 sin(0)] − cos 4(0) + sin 4(0)
102 51
1 25 50
= 𝑒 −2(0) [𝐶1 cos(0) + 𝐶2 sin(0)] − cos 4(0) + sin 4(0)
2 102 51
1 25
𝐶1 = +
2 102
38
𝐶1 =
51
Como parte del reposo y la velocidad inicial es cero, y se sabe que la velocidad es la derivada
de la posición, entonces se tiene que:
𝑦′(0) = 0
5 10
𝑦 = 𝐶1 𝑒 −3𝑡 cos 𝑡 + 𝐶2 𝑒 −3𝑡 sin 𝑡 + cos 4𝑡 + sin 4𝑡
18 9
10
𝑦 ′ = −3𝐶1 𝑒 −3𝑡 cos 𝑡 − 𝐶1 𝑒 −3𝑡 sen 𝑡 − 3𝐶2 𝑒 −3𝑡 sin 𝑡 + 𝐶2 𝑒 −3𝑡 cos 𝑡 − sin 4𝑡
9
40
+ cos 4𝑡
9
0 = −3𝐶1 𝑒 −3(0) cos( 0) − 𝐶1 𝑒 −3(0) sen(0) − 3𝐶2 𝑒 −3(0) sin(0) + 𝐶2 𝑒 −3(0) cos(0)
10 40
− sin 4(0) + cos 4(0)
9 9
40
0 = −3𝐶1 + 𝐶2 +
9
Como:
2
𝐶1 =
9
34
𝐶2 = −
9
Derivando la expresión 𝑡 = 0
86
𝐶2 = −
51
La ecuación en movimiento es :
38 86 25 50
𝑦 = 𝑒 −2𝑡 ( 𝑐𝑜𝑠 𝑡 − 𝑠𝑖𝑛 𝑡) − 𝑐𝑜𝑠 4𝑡 + 𝑠𝑖𝑛 4𝑡
51 51 102 51
2 −3𝑡 34 5 10
𝑦(𝑡) = 𝑒 cos 𝑡 − 𝑒 −3𝑡 sin 𝑡 + cos 4𝑡 + sin 4𝑡
9 9 18 9
Bibliografía