Semana 5 Teoría Alumnos

Ecuaciones
diferenciales
sistemas no homogéneas
y modelos con ED de
orden superior
Semana 05: Teoría
Profesores del curso:

Orlando Galarza Gerónimo
Patricia Reynoso Quispe
Profesores: Utec-Ciencias
Índice
1 Sistemas no homogéneos de
ED lineales
2 Sistemas masa-resorte libre
no amortiguado
3 Sistemas masa-resorte libre
amortiguado
4 Sistemas masa-resorte
amortiguado-forzado
Objetivos
Resolver sistemas de ED de primer orden no homogéneos.
Identificar los modelos lineales que describen el comportamiento de los
sistemas resorte-masa, en movimientos libres no amortiguado y amortiguado.
Resolver modelos lineales que describen el comportamiento de los sistemas
resorte-masa, en movimientos libres no amortiguado y amortiguado.
Identificar modelos lineales que describen el comportamiento de sistemas
masa-resorte en movimiento amortiguado-forzado.
Analizar y resolver modelos lineales que describen el comportamiento de
sistemas masa-resorte en movimiento amortiguado-forzado.
Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 2 / 67

Sistemas no
homogéneos de ED
lineales 1
Logros
?
? Resuelve sistemas ED de primer orden no homogéneo.

Sistemas de EDL no homogéneo
Recordemos que, en general los sistemas de ecuaciones diferenciales lineales
tienen la forma:
dx1
= a11 (t)x1 + a12 (t)x2 + · · · + a1n (t)xn + f1 (t)
dt
dx2
= a21 (t)x1 + a22 (t)x2 + · · · + a2n (t)xn + f2 (t)
dt
.. ..
. .
dxn
= an1 (t)x1 + an2 (t)x2 + · · · + ann (t)xn + fn (t)
dt
Los cuales se pueden escribir matricialmente como X ′ = AX + F (t), donde
x1 a11 (t) ··· a1n (t) f1 (t)
     
 x2  a21 (t) ··· a2n (t) f2 (t)
X =  . , A= . ..  , F (t) =  . 
  
..
  
 ..   .. . .   .. 
xn an1 (t) ... ann (t) fn (t)

Anteriormente hemos estudiado dichos sistemas cuando F = 0. Ahora analizare-
mos el caso cuando F ̸= 0. Tales sistemas de EDL son llamados no homogéneos.
En ese contexto, dado un sistema de EDL lineal no homogéneo:
X ′ = AX + F (t) (1)
se le asocia el sistema de EDL homogéneo:
X ′ = AX . (2)
Sean Xp una solución particular de (1) y sea
Xh = c1 X1 + c2 X2 + · · · + cn Xn ,
una solución de (2). Entonces una solución general de (1) será:
X = Xh + Xp

Resolver un sistema no homogéneo
Primero: Resolver el sistema homogéneo asociado utilizando las técnicas aprendi-
das anteriormente.
Segundo: Encontrar la solución particular mediante alguno de los siguientes
métodos:
Coeficientes indeterminados
Variación de parámetros

Método de coeficientes indeterminados
El método de coeficientes indeterminados puede ser utilizado para encontrar una
solución particular para el sistema de EDL X ′ = AX + F (t), solo cuando se cumplan:
A es una matriz de tamaño n × n con entradas constantes.

F (t) es un vector cuyas entradas son constantes, polinomios, exponenciales,
senos, cosenos, o alguna combinación finita de estas (mediante suma o
producto).
Así, este método consiste en proponer una solución particular Xp motivado por la
clase de funciónes que aparezcan en el vector F (t).

Ejemplo 1
Resuelva el sistema:
x ′ = −x + 2y − 8
y′ = −x + y + 3
Solución:
Lo representamos en la forma matricial

′ −1 2 −8 x
X = X+ donde X =
−1 1 3 y
Primero:
Resolvemos el sistema homogéneo asociado,

′ −1 2
X = X
−1 1
para tal fin, calculemos los eigenvalores de A

−1 − λ 2
P(λ) = det(A − λI) = det =0
−1 1−λ
λ2 + 1 = 0
Luego, los eigenvalores son λ1 = i y λ2 = −i (notemos que son conjugados)

Ahora debemos encontrar los eigenvectores,
Para λ1 = i : Debemos analizar el sistema
(
−1 − i 2 k1 0 (−1 − i)k1 + 2k2 = 0
= →
−1 1−i k2 0 −k1 + (1 − i)k2 = 0
→ k1 = (1 − i)k2

k1 1−i 1−i
Luego = k2 y así K1 = es un eigenvector asociado a
k2 1 1
λ1 = 0 + i.
Para λ2 = −i : Sabemos que para calcular el eigenvector asociado a λ2 basta

1−i 1+i
calcular K1 , así = es el eigenvector requerido.
1 1

Finalmente la solución (con coordenas complejas) es:

1 − i it 1 + i −it
X = c1 e + c2 e
1 1
| {z } | {z }
K1 K1
Además, usando lo definido anteriormente: B1 = 21 (K1 + K1 ), B2 = 2i (−K1 + K1 )

1 −1
B1 = B2 =
1 0
Entonces la solución general sin componentes complejas será: X = c1 X1 + c2 X2

donde:
X1 = B1 cos(βt) − B2 sen(βt) e αt

X2 = B2 cos(βt) + B1 sen(βt) e αt


1 −1 α −1 1
X = c1 cos βt − t
sin βt e + c2 cos βt + sin βt e αt
1 0 0 1
Reemplazando α = 0 y β = 1, tenemos

1 −1 −1 1
X = c1 cos t − sin t + c2 cos t + sin t
1 0 0 1

cos(t) + sin(t) sin(t) − cos(t)
= c1 + c2
cos(t) sin(t)

Segundo:
Para determinar la solución particular Xp , vemos que:

−8
F (t) =
3
es un vector constante, luego la solución particular se supone también constante

a
Xp =
b
entonces, derivando

0
Xp′ =
0

reemplazando en el sistema original

′ −1 2 −8 a
Xp = Xp + donde X =
−1 1 3 b
se tiene:
0 = −a + 2b − 8
0 = −a + b + 3
resolviendo el sistema se obtiene a = 14 y b = 11 de esta manera la solución

particular es

14
Xp =
11

Finalmente, la solución general del sistema de ED es X = Xh + Xp

cos(t) + sin(t) sin(t) − cos(t) 14
X = c1 + c2 +
cos(t) sin(t) 11

Para el alumno

6 1
′ 6t
X = X+ en < −∞, +∞ >
4 3 −10t + 4
Solución:
Primero:
Resolvemos el sistema homogéneo asociado.

′ 6 1
X = X
4 3

Usando el método de los valores y vectores propios, calculamos los eigenvalores
de A:

6−λ 1
det =0
4 3−λ
(λ − 2)(λ − 7) = 0
Luego, los eigenvalores son λ1 = 2 y λ2 = 7 Determinando sus vectores propios

asociados respectivamente

1 1
K1 = K2 =
−4 1

Por lo tanto la solución homogénea es

1 1 7t
Xc = c1 e 2t + c2 e
−4 1
Segundo:

6t 6 0
La forma de F (t) = = t+ nos hace suponer que Xp debe
−10t + 4 −10 4
tener la misma forma

a2 a
Xp = t+ 1
b2 b1
Sustituyendo esta forma de Xp en el sistema original dado

6 1 6t
Xp′ = X +
4 3 p −10t + 4

se obtiene:
   
a2 (6a2 + b2 + 6)t + 6a1 + b1 − a2
 = 
b2 (4a2 + 3b2 − 10)t + 4a1 + 3b1 − b2 + 4
y este sistema genera a su vez cuatro ecuaciones con cuatro incógnitas
6a2 + b2 + 6 = 0 6a1 + b1 − a2 =0
y
4a2 + 3b2 − 10 = 0 4a1 + 3b1 − b2 + 4 = 0
resolviendo el sistema,a partir de las dos primeras ecuaciones del lado izquierdo se
tiene a2 = −2 y b2 = 6.
Sustituyendo estos resultados en las siguientes dos ecuaciones, se obtiene
a1 = − 47 y b1 = 107 .

De esta manera se tiene una vector solución particular
   4
−2 −7
Xp =   t +  
6 10
7
Por lo tanto, la solución general X = Xh + Xp es:

       4
1 1 −2 −7
X = c1   e 2t + c2   e 7t +   t +  
−4 1 6 10
7

Para el alumno
x′ = 5x + 3y − 2e −t + 1
y′ = −x + y + e −t − 5t + 7
Sugerencia
Reconocer y representar a F (t) en términos matriciales
−2e −t + 1

−2 −t 0 1
F (t) = −t = e + t+
e − 5t + 7 −1 −5 7
Suponer que la solución Xp es de la forma

a a a
Xp = 3 e −t + 2 t + 1 .
b3 b2 b1
Método de variación de parámetros
El método de variación de parámetros, es una alternativa de resolver cuando
hay dificultades para conseguir la solución particular Xp con el método de los
coeficientes indeterminados.
Sean X1 , X2 , ..., Xn un conjunto fundamental de soluciones del sistema homogéneo

X ′ = AX , entonces la solución general
X = c1 X1 + c2 X2 + ... + cn Xn
se puede representar de la siguiente manera:
       
x11 x12 x1n c1 x11 + c2 x12 + · · · + cn x1n
x21  x22  x2n  c1 x21 + c2 x22 + · · · + cn x2n 
X = c1  .  + c2  .  + · · · + cn  .  = 
       
..
 ..   ..   ..  

. 
xn1 xn2 xnn xn1 + c2 xn2 + · · · + cn xnn
| {z }
Φ(t)C

es decir, la solución homogénea se puede escribir como
X = Φ(t)C (3)
donde C es un vector (de constantes) de tamaño n × 1 y Φ es una matriz de

tamaño n × n cuyas columnas consisten en los vectores X1 , X2 , . . . , Xn :
 
x11 x12 · · · x1n
x21 x22 · · · z2n 
Φ(t) =  .
 
. .. 
 . . 
xn1 xn2 · · · xnn
Llamaremos a Φ matriz fundamental.

Propiedades de la Matriz fundamental
Se requiere dos propiedades de la matriz fundamental
Una matriz fundamental Φ tiene inversa (det Φ(t) ̸= 0, para todo t)
Φ es una matriz fundamental del sistema X ′ = AX entonces Φ′ (t) = AΦ(t)
De forma análoga a variación de parámetros para una ecuación diferencial, se
plantea una solución particular para el sistema no homogéneo:
X ′ = AX + F (t)
como 
u1 (t)
u2 (t)
Xp = Φ(t)U (t) donde U (t) =  .  (4)
 
.
 . 
un (t)

Reemplazando Xp = Φ(t)U (t) en X ′ = AX + F (t) se obtiene:
Z
U (t) = Φ−1 (t)F (t)dt (5)
volviendo a la ecuación (4):

Z
Xp = Φ(t) Φ−1 (t)F (t)dt (6)
De esta forma se tiene la solución general del sistema de ecuaciones diferenciales
Z
X = Φ(t)C + Φ(t) Φ−1 (t)F (t)dt

Para el alumno

′ −3 1 3t
X = X + −t
2 −4 e
Solución:
Primero:
Resolvemos el sistema homogéneo asociado.

′ −3 1
X = X
2 −4

Usando el método de los valores y vectores propios, calculamos los eigenvalores
de A:

−3 − λ 1
det =0
2 −4 − λ
(λ + 2)(λ + 5) = 0
Luego, los eigenvalores son λ1 = −2 y λ2 = −5. Determinando sus vectores

propios asociados respectivamente

1 1
K1 = K2 =
1 −2

Por lo tanto la solución homogénea es

1 −2t 1
Xc = c1 e + c2 e −5t
1 −2
Damos la forma para reconocer los términos

−2t −5t
e e
Xc = c1 −2t +c2
e −2e −5t
| {z } | {z }
X1 X2
Y formamos la matriz fundamental Φ(t)

−2t
e −5t

e
Φ(t) = −2t
e −2e −5t

Segundo: Usamos variación de parámetros para determinar a Xp .
Determinamos la inversa de Φ(t)

2 −2t 1 2t

3e 3e
Φ−1 (t) =  
1 5t
3e − 13 e 5t
Utilizando la fórmula (6)
 −2t 2 1 2t
e −5t −2t
 
Z e Z 3e 3e

3t
Xp = Φ(t) Φ−1 (t)F (t)dt =     −t dt
e
e −2t −2e −5t 1 5t
3e
1 5t
−3e
| {z }
2te 2t + 13 e t
 −2t
e −5t
  
e Z
=     dt
e −2t −2e −5t 5t 1 4t
te − 3 e
La integral en el vector se realiza para cada componente, obteniéndose
te 2t − 12 e 2t + 13 e t
 −2t  6 27
e −5t

e 5t − 50
Z
Xp = Φ(t) Φ−1 (t)F (t)dt =   =
e −2t −2e −5t 1 5t 1 5t
5 te − 25 e − 12 e
1 4t 3
5t −
21
50
Por lo tanto, la solución general X = Xh + Xp es:
6 27
+ 14 e −t

5t −

1 −2t 1 50
X = c1 e + c2 e −5t +  
1 −2 3 21 1 −t
5t − 50 + 2e
    6  27   1 
1 1 5 50 4
⇒X = c1   e −2t + c2   e −5t +   t −   +   e −t
1 −2 3 21 1
5 50 2

2
SISTEMAS
MASA-RESORTE
LIBRE NO
AMORTIGUADO
Logros
?
?
Identifica los modelos lineales que describen el comportamiento de los sis-
temas resorte-masa, en movimientos libres no amortiguado y amortiguado.
Resuelve modelos lineales que describen el comportamiento de los sistemas
resorte-masa, en movimientos libres no amortiguado y amortiguado.

Sistema masa-resorte: Movimiento libre
Suponga que un resorte se suspende verticalemtne de un soperte rígido y luego se
le fija una masa m a su extremo libre.
El resorte se va estirar, la longitud de elon-

gación depende de la masa del objeto que
se coloque.
Ley de Hooke
El resorte ofrece una fuerza restau-
radora opuesta y proporcional a la elon-
gación s:
F = ks
donde k es la constante del resorte.

En el equilibrio mg = ks.
Si la masa se desplaza una distancia x de su posición de equilibrio la fruerza
restauradora del resorte es k (x + s).
Suponiendo que no hay fuerzas restauradoras que actúen sobre el sistema y
suponiendo que la masa vibra libre de otras fuerzas externas “Movimiento libre” se
puede aplicar la segunda ley de Newton:
d2x d2x
m = −k (s + x ) + mg ⇒ m = −ks − kx + mg
dt 2 dt 2
d2x k d2x
= − x ⇒ + ω2 x = 0 (7)
dt 2 m dt 2
donde ω2 = k /m. La ecuación auxiliar de la EDO (7) está dado por
r 2 + ω2 = 0 ⇒ r1,2 = ±i ω
Por lo tanto, la solución está dado por
2π
x (t ) = c1 cos(ωt ) + c2 sin(ωt ), periodo T = (8)
ω
Ejemplos
5
1Una masa de 15 kg alarga 32 m un resorte. En t = 0 se libera la masa desde un
2
punto que está a 3 m abajo de la posición de equilibrio con una velocidad ascendente
de 43 m/s. Determine la ecuación del movimiento. g = 10 m/s2 .
Solución:
Por la ley de Hooke:

5 N
150 = k ⇒ k = 960
32 m
k 960
⇒ ω2 = = = 64
m 15
Por lo tanto, la ecuación que debemos
resolver es
′′
x + 64x = 0
x(0) = 32 , x ′ (0) = − 43

La ecuación caracteristica es:
r 2 + 64 = 0, r1,2 = ±8i
Por lo tanto, la solución es
x(t) = c1 cos(8t) + c2 sin(8t) (9)
derivando
x ′ (t) = −8c1 sin(8t) + 8c2 cos(8t). (10)
Utilizando las condiciones iniciales en las ecuaciones (9) y (10), se obtiene
2 1
c1 = , c2 = −
3 6
Por lo tanto, la solución final es
2 1
x(t) = cos(8t) − sin(8t)
3 6
2 Una masa de 20 g alarga 5 cm un resorte. En t = 0 se libera la masa desde un punto
que está a 15 cm abajo de la posición de equilibrio con una velocidad descendente de
2 m/s. Determine la ecuación del movimiento. g = 10 m/s2 .
Solución:
Por la ley de Hooke (Convertimos a
unidades SI):
N
(0.02)(10) = k(0.05) ⇒ k=4
m
k 4
⇒ ω2 = = = 200
m 0.02
Por lo tanto, la ecuación que debemos
resolver es
′′
x + 200x = 0
x(0) = 0.15, x ′ (0) = 2

La ecuación caracteristica es:
√
r 2 + 200 = 0, r1,2 = ±10 2i
Por lo tanto, la solución es
√ √
x(t) = c1 cos 10 2t + c2 sin 10 2t (11)
derivando
√ √ √ √
x ′ (t) = −10 2c1 sin 10 2t + 10 2c2 cos 10 2t . (12)
Utilizando las condiciones iniciales en las ecuaciones (11) y (12), se obtiene

1
c1 = 0.15, c2 = √
5 2
Por lo tanto, la solución final es
√ 1 √
x(t) = 0.15 cos 10 2t + √ sin 10 2t
5 2
Para el alumno
Se tiene un resorte de longitud inicial L0 = 10 m unido a un bloque de masa
m = 10 kg. El sistema se libera desde el reposo a una distancia de 2 m a la
derecha de la posición de equilibrio. Halle la ecuación que describe la posición del
objeto si k = 360 N /m. Considere que el origen del sistema coordenado es el
punto de equilibrio.
Respuesta:
x(t) = 2 cos (6t)
3
SISTEMAS
MASA-RESORTE
LIBRE
AMORTIGUADO
Sistema masa-resorte: Movimiento libre
amortiguado
El sistema libre no amortiguado o movimiento armónico libre es un caso irreal, a
menos que la masa se suspenda en un vacío perfecto. Ya que por lo menos habrá
una fuerza de resistencia debido al medio circundante (Fricción).
En el estudio de la mecánica, se considera que las fuerzas de amortiguamiento son
proporcional a la velocidad instantánea.
d2x dx
m 2
+γ + kx = 0, (13)
dt dt
donde γ es la constante de amortiguamiento.

Modelamiento
Plantee la ecuación de un sistema masa resorte suspendido verticalmente si el aire
ofrece una resistencia proporcional a la velocidad.
De la segunda ley de Newton

X
m ẍ = F.
Por lo tanto
m ẍ = mg − γ ẋ − k(s + x)
m ẍ + γ ẋ + kx = mg − ks
⇒ m ẍ + γ ẋ + kx = 0

Análisis de la Ecuación Diferencial
Del modelamiento usando la segunda Ley de Newton y Ley de Hooke, considerando
además el amortiguamiento proporcional a la velocidad, se tiene la siguiente ED:
d2x dx
m +γ + kx = 0, (14)
dt 2 dt
Esta EDO tiene coeficientes constantes, así que podemos usar la ecuación auxiliar:
p
2 γ γ 2 − 4mk p
mr + γ r + k = 0 ⇒ r1,2 = − ± = −λ ± λ2 − ω2 (15)
2m 2m
γ k
donde λ = 2m y ω2 = m . Dependiendo del valor del discriminante △ = λ2 −
2
ω , existen tres posibles formas de movimiento (Soluciones reales y diferentes,
soluciones iguales o repetidas y soluciones complejas conjugadas).

Caso I: ∆ = λ2 − ω2 > 0
√ √
−λt λ2 −ω2 t − λ2 −ω2 t
x(t) = e c1 e + c2 e
Es un movimiento uniforme y no oscilatorio.El sistema está sobreamortiguado:

porque el coeficiente de amortiguamiento γ es grande comparado con la constante
del resorte k.
Caso II: ∆ = λ2 − ω2 = 0
x(t) = e −λt (c1 + c2 t)
El sistema está criticamente amortiguado: porque cualquier disminución en la

fuerza de estiramiento daria como resultado un movimiento oscilatorio.
Caso III: ∆ = λ2 − ω2 < 0
p p
x(t) = e −λt c1 cos w2 − λ2 t + c2 sen w2 − λ2 t
El sistema está subamortiguado porque el coeficiente de amortiguamiento es

pequeño comparado con la constante del resorte.

Ejemplos
1 Una masa de 0.5 kg se une a un resorte que tiene una longitud de 5 m de
largo. En equilibrio mide 6 m. Si al inicio la masa se libera desde el reposo
en un punto 2 m arriba de la posición de equilibrio. Encuentre la posición x (t )
en todo momento si además se sabe que el medio circundante ofrece una
resistencia numéricamente igual a la velocidad instantáneas. g = 10 m /s2 .
Solución:
De la ley de Hooke:
N
0.5(10) = k (6 − 5) ⇒ k = 5 .
m
Por lo tanto
d2x dx 1 d2x dx
m 2 = −kx − γ ⇒ 2
+ + 5x = 0.
dt dt 2 dt dt
La solución de esta ecuación sujeta a las condiciones: x (0) = −2, x ′ (0) = 0, está
dado por:
−t 2
x (t ) = e −2 cos(3t ) − sin(3t )
3
2 Grafique la solución de la ecuación mx ′′ + γ x ′ + kx = 0, sujeta a x(0) = −2, x ′ (0) = 0 para
diferentes valores de m, γ , k.
x ′′ + 2x ′ + 10x = 0
5
x ′′ + x ′ + 10x = 0
2
x ′′ + 3x ′ + 10x = 0
x ′′ + 4x ′ + 10x = 0
x ′′ + 6x ′ + 10x = 0

en todo momento si además se sabe que el medio circundante ofrece una
resistencia igual a 3.5 veces la velocidad instantánea. g = 10 m /s2 .
Solución:
De la ley de Hooke:
N
0.5(10) = k (6 − 5) ⇒ k=5 .
m
Por lo tanto
d2x dx 1 d2x 7 dx
m 2
= −kx − γ ⇒ 2
+ + 5x = 0.
dt dt 2 dt 2 dt
dado por:
10 4
x (t ) = − e −2t + e −5t
3 3
Gráfica de las ecuaciones dadas sujetas a x(0) = −2, x ′ (0) = 0:
x ′′ + 6x ′ + 10x = 0
x ′′ + 7x ′ + 10x = 0
x ′′ + 20x ′ + 10x = 0
La curva roja es la solución del problema anterior.

en todo momento √ si además se sabe que el medio circundante ofrece una
resistencia igual a 10 veces la velocidad instantánea. g = 10 m /s2 .
Solución:
De la ley de Hooke:
N
0.5(10) = k (6 − 5) ⇒ k=5 .
m
Por lo tanto
d2x dx 1 d 2 x √ dx
m 2
= −kx − γ ⇒ + 10 + 5x = 0.
dt dt 2 dt 2 dt
dado por: √ √ √
x (t ) = −2e − 10t − 2 10te − 10t
5 Una masa de 0.25 kg alarga 58 m a un resorte, suponiendo que una fuerza amor-
tiguadora es igual a dos veces la velocidad instantánea. Determine la ecuación del
movimiento si la masa inicial se libera desde la posición de equilibrio con una velocidad
ascendente de 3 m/s. Considere g = 10 m/s2 .
Solución:
Resultado: x(t) = −3te −4t

Plano inclinado
Plantee la ecuación diferencial y sus condiciones iniciales, el resorte se estira y
luego se suelta.

Para el alumno
En un plano inclinado de 30◦ se encuentra un bloque de piedra de 100 kg en reposo
sujetado por un resorte cuya constante elástica es de 2500 N /m. Suponiendo que
no existe rozamiento determinar la elongación del muelle. Considere g = 10 m/s2 .
Solución:
Respuesta:
s = 0.2 m

Para el alumno
Se tiene el sistema del ejemplo anterior, donde la longitud inicial del resorte es
L0 = 4 m. Se coloca una objeto de masa m = 8 kg en el extremo inferior del
resorte y este se estira hasta que dicho extremo está en la posición x = 6 del
sistema coordenado donde el vértice está en el punto más alto del plano inclinado
y el eje x está en la dirección de elongación del resorte. Luego se libera desde el
reposo el sistema (masa-resorte) en la posición x = 5. Si el aire ejerce una fuerza
de resistencia numéricamente igual a tres veces la velocidad instantánea, plantee el
problema de valores iniciales que describe al objeto en el sistema de coordenadas
dado. Considere g = 10 m/s2 y el ángulo del plano (liso) inclinado es θ = 30◦ .
Respuesta (
8x ′′ + 3x ′ + 20x = 120
x(0) = 5, x ′ (0) = 0

4
SISTEMAS
MASA-RESORTE
AMORTIGUADO-
FORZADO
Logros
?
?
Identifica modelos lineales que describen el comportamiento de sistemas
masa-resorte en movimiento amortiguado-forzado.
Analiza y resuelve modelos lineales que describen el comportamiento de
sistemas masa-resorte en movimiento amortiguado-forzado.

Sistema masa-resorte forzado
En esta sección investigaremos qué sucede cuando una fuerza externa f (t), la cual
depende sólo del tiempo, es aplica al sistema masa-resorte.
d2x dx
m 2
+γ + kx = f (t), (16)
dt dt
donde γ es la constante de amortiguamiento.
La discusión está restringida a funciones de forzamiento de tipo sinusoidales, es
decir:
f (t) = F cos(ωt) (17)
¿Por qué es importante considerar esta fuerza del tipo sinusoidal?
Este tipo de fuerzas suelen ser periódicas (vibraciones)

Corrientes alternas en los circuitos.
El estudio de la respuesta forzada puede ser dividido en dos casos, dependiendo
de si la fricción está presente, o no.
Caso 1: La fricción no está presente (γ = 0). La ecuación resulta ser:
d2x
m + kx = F cos(ωt)
dt 2
cuya solución está dado por:
k
⇒ x(t) = xp + C1 cos(ω0 t) + C2 sin(ω0 t), siendo ω0 2 = .
m
El valor ω0 se le conoce como frecuencia natural. Para cuando ω0 ̸= ω, entonces
F
xp = cos(ωt)
k − m ω2
Para cuando ω0 = ω se dice que el sistema está en resonancia pura.

Caso 2: La fricción está presente (γ ̸= 0). La ecuación resulta ser:
d2x dx
m +γ + kx = F cos(ωt)
dt 2 dt
donde la solución resulta ser
x(t) = xH (t) + xp (t)
Para determinar la solución particular podemos utilizar el método de variación de

parámetros o el método de coeficientes indeterminados.

Ejercicios
1 Una masa de 15 kg alarga 5/32 m un resorte. En t = 0 se libera la masa desde un
punto que está 2 m arriba de la posición de equilibrio con una velocidad ascendente
de 4 m/s. Esta masa tiene asociado una fuerza externa de f (t) = −945 cos(t).
Determine la ecuación de movimiento. Considere g = 10 m/s2 .
Solución:
Por Hooke:
5 N
150 = k ⇒ k = 960 .
32 m
De la segunda ley de Newton:
d2x
m = −kx + f (t)
dt 2
(
x ′′ + 64x = −63 cos(t)
x(0) = −2, x ′ (0) = −4
La solución de PVI anterior es
1
x(t) = − cos(8t) − sin(8t) − cos(t)
2

Podemos comparar el resultado anterior con el resultado de la ecuación
(
x ′′ + 64x = −63 cos(8t)
x(0) = −2, x ′ (0) = −4
la cual está en resonancia pura (¿Por qué?). En este caso el resultado es:
1 63
x(t) = −2 cos (8t) − sin (8t) − t sin (8t)
2 16

Conclusiones
?
?
1 Los métodos de coeficientes indeterminados y variación de parámetros
utilizados para determinar la solución particular de EDOs lineles no
homogéneas, se puede adaptar para determinar la solución particular
para un sistemas de ecuaciones lineales no homogéneas.
2 La dirección del eje x positivo puede estar en cualquier dirección.
3 En oscilaciones masa-resorte la fuerza del peso no participa de la
ecuación de movimiento.

Gracias

Semana 5 Teoría Alumnos

Cargado por

Copyright:

Formatos disponibles

Semana 5 Teoría Alumnos

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Semana 5 Teoría Alumnos

Cargado por

Copyright:

Formatos disponibles

Ecuaciones

Profesores del curso:

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 2 / 67

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 4 / 67

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 5 / 67

se le asocia el sistema de EDL homogéneo:

Sean Xp una solución particular de (1) y sea

una solución de (2). Entonces una solución general de (1) será:

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 6 / 67

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 7 / 67

A es una matriz de tamaño n × n con entradas constantes.

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 8 / 67

Lo representamos en la forma matricial

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 10 / 67

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 11 / 67

Además, usando lo definido anteriormente: B1 = 21 (K1 + K1 ), B2 = 2i (−K1 + K1 )

Entonces la solución general sin componentes complejas será: X = c1 X1 + c2 X2

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 12 / 67

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 13 / 67

Para determinar la solución particular Xp , vemos que:

es un vector constante, luego la solución particular se supone también constante

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 14 / 67

resolviendo el sistema se obtiene a = 14 y b = 11 de esta manera la solución

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 15 / 67

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 16 / 67

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 17 / 67

Luego, los eigenvalores son λ1 = 2 y λ2 = 7 Determinando sus vectores propios

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 18 / 67

tener la misma forma    

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 19 / 67

y este sistema genera a su vez cuatro ecuaciones con cuatro incógnitas

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 20 / 67

Por lo tanto, la solución general X = Xh + Xp es:

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 21 / 67

Sean X1 , X2 , ..., Xn un conjunto fundamental de soluciones del sistema homogéneo

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 23 / 67

donde C es un vector (de constantes) de tamaño n × 1 y Φ es una matriz de

Llamaremos a Φ matriz fundamental.

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 24 / 67

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 25 / 67

volviendo a la ecuación (4):

De esta forma se tiene la solución general del sistema de ecuaciones diferenciales

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 26 / 67

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 27 / 67

Luego, los eigenvalores son λ1 = −2 y λ2 = −5. Determinando sus vectores

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 28 / 67

Damos la forma para reconocer los términos

Y formamos la matriz fundamental Φ(t)

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 29 / 67

Determinamos la inversa de Φ(t)

Por lo tanto, la solución general X = Xh + Xp es:

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 31 / 67

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 33 / 67

El resorte se va estirar, la longitud de elon-

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 34 / 67

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales February 6, 2024 36 / 67

Por lo tanto, la solución es

x(t) = c1 cos(8t) + c2 sin(8t) (9)

tener la misma forma