Reportes científicos de la FACEN, Enero - Junio 2023, 14(1): 91–97

999-999

http://doi.org/10.18004/rcfacen.2023.14.1.91
Artículo Original
Resolviendo la ecuación del movimiento armónico amortiguado mediante algunos métodos
numéricos

Solving the equation of dampening harmonic motion using some numerical methods

Salustiano Vega1,*, , Jorge Jara2, , Osvaldo Vega3, , Gustavo González3, & Crispín Vargas2,
1
Universidad Nacional de Asunción, Facultad Politécnica, Departamento de Ciencias Básicas, San Lorenzo, Paraguay.
2
Universidad Nacional de Asunción, Facultad de Ciencias Exactas y Naturales, Departamento de Física, San Lorenzo, Paraguay.
3
Universidad Nacional de Asunción, Facultad de Ciencias Exactas y Naturales, Departamento de Matemática, San Lorenzo, Paraguay.
*Autor correspondiente: [email protected].

Resumen: El modelo del movimiento armónico amortiguado es uno de los temas abordados en el área de la
matemática y de física, su representación mediante ecuaciones diferenciales es el motivo principal de su estudio.
Usualmente se analiza el caso más simple para obtener una solución analítica real, en el que las técnicas usuales
de resolución impartidas para describir el movimiento armónico amortiguado consideran pequeñas intensidades
de amortiguamiento. Sin embargo, el análisis mediante ciertos métodos numéricos que resuelven de manera
aproximada ecuaciones diferenciales ordinarias permitirá resolver numéricamente este modelo y proporcionar
diferentes técnicas de resolución. En este trabajo se resolverán numéricamente la ecuación diferencial que describe
el movimiento armónico amortiguado. Se propondrán algunas condiciones iniciales a la ecuación diferencial y se
implementarán las funciones ODE23, ODE45 y ODE 113 de Matlab y algunos algoritmos numéricos clásicos,
tales como los métodos, de Euler Back de primer orden, de Runge Kutta de segundo orden, de Adams Moulton
de tercer orden y de Runge Kutta de cuarto orden. Se realizará finalmente, una comparación de los algoritmos
numéricos y de las funciones de Matlab implementadas con la solución analítica exacta de la ecuación diferencial,
así como con la solución exacta para la velocidad vertical de la masa.
Palabras clave: Ecuaciones diferenciales, Métodos numéricos, Algoritmo.

Abstract: The model of damped harmonic motion is one of the topics addressed in the area of ​​mathematics and
physics, its representation through differential equations is the main reason for its study. The simplest case is
usually analyzed to obtain a real analytical solution, in which the usual resolution techniques taught to describe
damped harmonic motion consider small damping intensities. However, analysis using certain numerical methods
that approximate ordinary differential equations will allow us to solve this model numerically and provide differ-
ent resolution techniques. In this work, the differential equation that describes the damped harmonic motion will
be solved numerically. Some initial conditions for the differential equation will be proposed and the functions
ODE23, ODE45 and ODE 113 of Matlab and some classical numerical algorithms will be implemented, such as
the first-order Euler Back methods, the second-order Runge Kutta method, of Adams Moulton of the third order
and de Runge Kutta of the fourth order. Finally, a comparison of the numerical algorithms and the implemented
Matlab functions will be made with the exact analytical solution of the differential equation, as well as with the
exact solution for the vertical velocity of the mass.
Keywords: Differential equations, Numerical methods, Algorithm.

Introducción tiguado se desarrolla por medio de las ecuaciones

La apariencia del movimiento armónico amortigua- diferenciales que, en la mayoría de los casos, por
do generalmente está relacionada con un cuerpo de medio de métodos usuales de resolución, se obtiene
cierta masa unido a un extremo de una fuente elás- una solución analítica para pequeñas intensidades
tica, permaneciendo fijo el otro extremo de dicha de amortiguamiento.
fuente. De hecho, es un tema estudiado en modelos Los trabajos realizados hasta el momento presen-
matemáticos y físicos. tan diferentes enfoques y aplicaciones, por ejemplo,
En términos matemáticos según Costa et al., en (Pascual et al.,2008) se realiza una simulación
(2011), el estudio del movimiento armónico amor- en Maple de la modelación del movimiento tanto
Recibido: 10/02/2023 Aceptado: 14/04/2023
2078-399X/2023 Facultad de Ciencias Exactas y Naturales - Universidad Nacional de Asunción, San Lorenzo, Paraguay. Este
es un artículo de acceso abierto bajo la licencia CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/deed.es).
92 Reportes científicos de la FACEN, Enero - Junio 2023, 14(1): 999-999
91–97

armónico simple como amortiguado, en (Serrano

et al.,2016) se presenta un caso del modelado y
solución de una ecuación diferencial ordinaria de
segundo orden considerando un sistema masa re-
sorte con movimiento libre amortiguado, resuelta
analíticamente y simulada con el apoyo del soft-
ware GeoGebra, en González (2009) se estudian
las ecuaciones diferenciales ordinarias lineales
de segundo orden de coeficientes constantes que
describe este problema con algunas aplicaciones
físicas, en Rodríguez (2019). se construye un mo-
delo matemático para explicar o describir la acción
de los terremotos sobre edificios de varios pisos,
mediante ecuaciones diferenciales de segundo
orden de sistema masa resorte.
La motivación de este trabajo es encontrar
soluciones numéricas adecuadas para la ecuación
diferencial que describe el movimiento armónico
amortiguado, implementando algunas funciones
ODE de Matlab, así como métodos numéricos de
diferentes órdenes, que permiten resolver ecua- Figura 1. Sistema de movimiento amortiguado.
ciones diferenciales ordinarias con condiciones
estará sometida a una fuerza de origen elástica dada
iniciales establecidas.
por la ecuación F = – k x, siendo k la constante
El objetivo de este artículo es resolver numé-
elástica del resorte que depende de la geometría y
ricamente la ecuación diferencial del movimiento
del material del cual está hecho y x es la elongación
armónico amortiguado. Se implementarán las
del resorte a partir de su posición de equilibrio.
funciones ODE23, ODE45 y ODE113 de Matlab
El efecto del peso no se considera porque se
y algunos métodos numéricos de orden 1, de orden
compensa con la deformación inicial que sufre el
2, orden 3 y orden 4 con condiciones iniciales a
resorte al colgarle la masa m y de esta manera, se
la ecuación diferencial del problema estudiado
dice que el sistema está en la posición de equilibrio,
y posteriormente, se compararán con la solución
se toma como positivo el desplazamiento hacia
analítica obtenida para pequeñas intensidades de
abajo, por eso los signos de la fuerza elástica y la
amortiguamiento.
retardadora se consideran negativos, pues actúan
Según (Sears et al., 2013) y Molina (2016), un
hacia arriba.
sistema amortiguado es un esquema sencillo que
Para obtener la ecuación de movimiento se
se puede observar en la Fig. 1.
aplica la segunda ley de Newton al sistema masa
Se puede observar un cuerpo de masa m sus-
resorte y se obtiene el siguiente resultado:
pendido de un resorte de constante k, que oscila,
tiene movimiento armónico amortiguado, cuando la dx
b kx ma
fricción del medio sobre la masa fr es directamente dt
dx
proporcional a la velocidad, esto es f r , b d 2x
dt donde a  es la aceleración de la masa, de
dx dt 2
siendo b la constante de amortiguamiento y es
dt esta forma, el movimiento armónico amortiguado
la velocidad de la masa. de la Fig. 1, quedará descrito por la ecuación di-
Como la masa está sujeta a un resorte, entonces
Salustiano Vega et al.
 Reportes científicos de la FACEN, Enero - Junio 2023, 14(1): 999-999
91–97 93

ferencial: b
d 2x dx x0 v0
m
dt 2
b
dt
kx 0............(1) arctg 2m
wx0
De esta manera, mediante la aplicación de los
métodos usuales de resolución, la solución analítica
de la ecuación (1) para pequeñas intensidades de Métodos
Métodos numéricos para ecuaciones diferenciales
amortiguamiento b k , está dada por:
2m m En la mayoría de los casos, es muy difícil encontrar
b soluciones exactas de una ecuación diferencial por
t medio de métodos usuales, en este contexto, la apli-
x t Ae 2m cos wt ; con cación de métodos numéricos es indispensable para
obtener soluciones aproximadas. Generalmente
k b2 según (Prima et al.,2000), León (2015) y (Atkinson
w ............(2)
m 4m 2 et al., 2011), los métodos numéricos son utilizados
En la ecuación (2), las constantes de integración, añadiendo condiciones iniciales a la ecuación dife-
A, y ϕ (amplitud y desfase) se calculan de mediante rencial a estudiar.
las condiciones iniciales: Considere el problema de valor inicial de primer
orden:
2
b y ' f x , y ; y x0 y0 ............(3)
x 0 x0 x0 v0
A x0 1 2m ; Sea y = ϕ(x) una solución exacta única de la
x' 0 v0 wx0 ecuación (3) en algún intervalo centrado en x = x0,
con el supuesto de que f y fy son continuas en al-

Figura 2. Solución exacta y puntos aproximados.

Resolviendo la ecuación del movimiento armónico amortiguado mediante algunos métodos numéricos
94 Reportes científicos de la FACEN, Enero - Junio 2023, 14(1): 999-999
91–97

gún rectángulo del plano xy que contenga al punto Método de Euler backward de orden 1.
(x0, y0).
Para encontrar soluciones aproximadas, en mé- Datos: x0 , y0 ,h , f x , y , N .
todos numéricos se emplean algoritmos iterativos Algoritmo:
de resolución que calculan los valores aproximados
xn xn h.
y0, y1, y2, …, yn, … de la solución y = ϕ(x) de la 1
ecuación (3) en un conjunto de puntos x0, x1, x2, …, yn 1 yn h f xn 1 , yn 1 .
xn, …, como se muestra en la Fig. 2. n 0 ,1, 2 ,3,… , N 1.
En la Fig. 2 se tiene que y(x 0) = y por el
problema de valor inicial, pero ϕ(xn) ≠ yn para Método de Runge Kutta de orden 2.
n ≥ 1, además, ' x0 f x0 , y0 y0 ' pero
Datos: x0 , y0 ,h , f x , y , N .
' xn f xn , xn yn ' p a r a n ≥ 1 .
Algoritmo:
También sobre el eje x se realizará un espaciamiento xn 1 xn h.
o tamaño de paso uniforme h, así entonces
1 1
x1 x0 h; x2 x0 2h;… ; xn x0 nh. yn 1 yn h f xn h , yn h f xn , yn .
2 2
Los métodos numéricos se dividen en dos n 0 ,1, 2 ,3,…
… , N 1.
grupos, métodos de un paso donde se utiliza la
Método de Adams Moulton de orden 3.
solución aproximada yn en xn para obtener la so-
lución aproximada yn+1 en xn+1 y en los métodos Datos: x0 ,x1 , y0 , y1 , h , f x , y , N .
multipaso la solución yn+1 se obtiene utilizando el Algoritmo:
valor o los valores de solución calculados en los xn 1 xn h.
puntos anteriores. 1
yn 1 yn
h 5 f xn 1 , yn 1 8 f xn , yn f xn 1 , yn 1 .
El orden de un método numérico es un entero 12
n 1, 2 ,3,… , N 1.
no negativo p de modo que xn yn c h p para
alguna constante c. Método de Runge Kutta de orden 4.
En este artículo, se utilizarán métodos numéricos Datos: x0 , y0 ,h , f x , y , N .
de un paso y multipaso conjuntamente, además se
Algoritmo:
implementarán y se compararán según el orden de
xn xn h.
los mismos. 1
1 h h
yn 1 yn h f xn , yn 2 f xn , yn f xn , yn
6 2 2
Métodos numéricos implementados en la reso-
1 h h h h
lución 6
h 2 f xn
2
, yn
2
f xn
2
, yn
2
f xn , yn

Los siguientes métodos numéricos implementados 1 h h h h

h f xn h , yn hf xn , yn f xn , yn f xn , yn .
en este artículo, resuelven numéricamente el proble- 6 2 2 2 2
ma de valor inicial (3), en un intervalo cerrado [a,b], n 0 ,1, 2 ,3,… , N 1.
donde x0 = a, xn+1 = b y en N pasos. Primeramente
se divide el intervalo [a,b] en N partes iguales de Las funciones ODE de Matlab.
b a Las funciones ODE de Matlab proporcionan una
longitud h y posteriormente se aplica el
N solución numérica de ecuaciones diferenciales or-
algoritmo correspondiente a cada método. dinarias con problemas de valores iniciales, cuando
es difícil encontrar la solución analítica de la men-

Salustiano Vega et al.

 Reportes científicos de la FACEN, Enero - Junio 2023, 14(1): 999-999
91–97 95

T devuelve el tiempo del vector de columna.

Y devuelve el vector de la columna de so-
lución correspondiente a T.

Resultados y discusión
Para la implementación de los algoritmos mencio-
nados anteriormente, se consideraron la ecuación
d 2x dx
diferencial m b kx
0 que describe
dt 2 dt dx
Figura 3. Soluciones numéricas para x(t) en el intervalo [0,30], el movimiento armónico amortiguado, donde
con h = 0.0075, N = 400, m = 1, b = 0.2 y k = 9.81. dt
representa la velocidad del movimiento vertical de
cionada ecuación diferencial. Se implementaron la masa, con las siguientes condiciones iniciales,
debido a la ineficiencia del método de Euler en las x 0 1 y x' 0 0 .
aplicaciones. [9] Se presentarán a continuación, los resultados
En este trabajo se implementarán las funcio- obtenidos para la solución analítica x(t) y para la
nes ODE23, ODE45 y ODE113 para resolver la dx
ecuación (1) con algunas condiciones iniciales velocidad del movimiento vertical de la masa .
dt
establecidas. A) Resultados de la implementación de las
El uso de las funciones ODE en código de Mat- funciones ODE de Matlab.
lab se describe a continuación.
En la Fig. 3 se muestra la comparación de
T ,Y ODE odefun ,tspan ,Y0 las funciones ODE23, ODE45 y ODE113
con la solución exacta x(t), mediante una
odefun es un identificador de función, que
experimentación numérica en el intervalo
puede ser un nombre de archivo de
[0,30].
función, un identificador de función anó-
En la Fig. 4 se muestra la compara-
nimo o un nombre de función en línea.
ción de las funciones ODE23, ODE45 y
tspan es el intervalo x0 , xn 1 o una serie ODE113 con la velocidad del movimiento
de puntos dispersos x0 , x1 ,… , xn 1 .
dx
vertical de la masa , mediante una la
dt
Y0 es el vector de valor inicial. experimentación numérica en el intervalo
[0,30].

dx
Figura 4. Soluciones numéricas para en el intervalo
dt Figura 5. Soluciones numéricas para x(t) en el intervalo
[0,30], con h = 0.0075, N = 400, m = 1, b = 0.2 y k = 9.81. [0,30], con h = 0.0075, N = 400, m = 1, b = 0.2 y k = 9.81.
Resolviendo la ecuación del movimiento armónico amortiguado mediante algunos métodos numéricos
96 Reportes científicos de la FACEN, Enero - Junio 2023, 14(1): 999-999
91–97

B) Resultados de la implementación de los solución analítica exacta del modelo estudiado y

métodos numéricos de orden 1, de orden con la velocidad del movimiento vertical de la masa.
2, de orden 3 y de orden 4 descritos en A continuación, se presentan las conclusiones
este trabajo. sobre la comparación de las distintas funciones
En la Fig. 5 se muestra la comparación ODE de Matlab con las soluciones exacta y de
de las soluciones numéricas de los métodos velocidad vertical.
eulerbackward de orden 1, Runge Kutta de • Las funciones ODE23, ODE45 y ODE113
orden 2, Adams Moulton de orden 3 y Run- incorporadas en Matlab permiten resolver de
ge Kutta de orden 4 con la solución exacta manera eficiente ecuaciones diferenciales ordi-
x(t), mediante una la experimentación nu- narias con problema de valor inicial, sobre la
mérica en el intervalo [0,30]. implementación de estas funciones al modelo
En la Fig. 6 se muestra la comparación estudiado se puede concluir que, la función
de las soluciones numéricas de los métodos ODE113 presenta más aproximación a las solu-
Euler backward de orden 1, Runge Kutta ciones analíticas exacta y de velocidad vertical,
de orden 2, Adams Moulton de orden 3 y debido a su número de orden.
Runge Kutta de orden 4 con la velocidad
dx • En general, no hay mucha diferencia entre las
del movimiento vertical de la masa , funciones ODE de Matlab y se pueden ser
dt consideradas como las soluciones numéricas
mediante una la experimentación numérica
más exactas para ecuaciones diferenciales or-
en el intervalo [0,30].
dinarias con problema de valor inicial.
Con relación a los métodos numéricos aplicados,
el método de orden 1 Eulerback, es el que presenta
menor aproximación a las soluciones analíticas
exactas y de velocidad vertical, el resultado del
método numérico de orden 2, Runge Kutta, pre-
senta una notable mejoría respecto al resultado del
método de orden 1. Sin embargo, los métodos de
orden 3 y de orden 4 son los que se aproximan de
manera eficiente a las soluciones analíticas exactas
y de velocidad vertical, siendo el método Runge
dx Kutta de orden 4 la mejor.
Figura 6. Soluciones numéricas para en el intervalo
dt Para la resolución numérica de este tipo de
[0,30], con h = 0.0075, N = 400, m = 1, b = 0.2 y k = 9.81. problema de valor inicial se recomienda la imple-
mentación de cualesquiera de las funciones ODE
Conclusiones de Matlab y de los métodos numéricos, el método
En este trabajo se han realizado soluciones numéri- de Runge Kutta de orden 4 es el indicado.
cas para pequeñas intensidades de amortiguamiento
del problema que describe el movimiento armónico Contribución de los autores
amortiguado representado mediante ecuaciones Todos los autores contribuyeron de manera equita-
diferenciales, la cual presenta solución analítica tiva en la elaboración de este artículo.
real para la condición dada, posteriormente, se han
implementados algunos métodos numéricos, así Conflictos de interés
como las funciones ODE23, ODE45 y ODE113 Los autores declaran no tener conflictos de interés.
de Matlab, para realizar una comparación con la

Salustiano Vega et al.

 Reportes científicos de la FACEN, Enero - Junio 2023, 14(1): 999-999
91–97 97

Literatura citada Exposición. I Jornada de Investigación y

Atkinson, K., Han, W. & Stewart, D. E. (2011). Nu- Transferencia de la Facultad de Ingeniería.
merical solution of ordinary differential La Plata: Universidad Nacional de La Plata.
equations. Hoboken: John Wiley & Sons. 8 pp.
xii + 252. Pascual, F.E.B., Mariño, M.L. & Sada, A.E.G.
Di Prima, B. & Boyce, W. (2000). Ecuaciones (2008). Análisis del Movimiento Oscilato-
diferenciales y problemas con valores en rio. Reporte Técnico. Veracruz: Instituto
la frontera. México, D.F.: Limusa / Grupo Tecnológico y de Estudios Superiores de
Noriega Editores. 758 pp. Monterrey. 16 pp.
González, H. (2009). Ecuaciones diferenciales Rodríguez Chuquimango, S.P. (2019). Modela-
ordinarias lineales de segundo orden de miento matemático de la acción de los terre-
coeficientes constantes y algunas aplica- motos sobre edificios de varios pisos usando
ciones físicas. Repositorio Institucional, sistemas de ecuaciones diferenciales de
Universidad de Carabobo. [Consulted: segundo orden. Informe Final de Proyecto de
18.vi.2022]. <http://riuc.bc.uc.edu.ve/ Investigación. Callao: Universidad Nacional
handle/123456789/756>. del Callao. 79 pp.
León Camejo, P.D. (2015). Métodos numéricos Sears, F. W., Zemansky, M. W., Young, H. D. &
para ecuaciones diferenciales ordinarias: Freedman, R. A. (2013). Física Universitaria:
métodos Runge-Kutta explícitos. Trabajo de Volumen 1. (13ª Ed.). Naucalpan de Juárez:
Grado. San Cristóbal de La Laguna: Univer- Pearson Educación de México. xxvii + 686
sidad de La Laguna. 43 pp. + 25 pp.
Molina M.C. (2016). Métodos Matemáticos de la Serrano, M.D.C.C., Oliver, E.B.V., Resendiz,
Física II: Ecuaciones Diferenciales y Fun- S.M. & Galván, W.G.A. (2016). Sistema
ciones Especiales. Texto de Curso. Granada: masa resorte con movimiento libre amor-
Universidad de Granada. 152 pp. tiguado, casos: sobreamortiguado, críti-
Costa, V.A., Bordogna, C.M. & Torroba, P.L. camente amortiguado y subamortiguado,
(2011). El movimiento armónico simple y su modelado y solución, con el apoyo de
las ecuaciones diferenciales de segundo or- GeoGebra. Pistas Educativas, 38(121):
den lineal: su articulación en la enseñanza. 61–80.

Resolviendo la ecuación del movimiento armónico amortiguado mediante algunos métodos numéricos