simulacion y studio de n osciladores acoplados

simulation and study of n coupled oscillators

H. Torres 1
1
Estudiante del programa de Física, Universidad del Atlántico, Barranquilla-Colombia

Resumen

se analizó la simulación de un sistema de n osciladores acoplados, para observar el cambio de este sistema en el tiempo,
para esto se examinó el sistema desde 3 osciladores hasta 100 de estos acoplados, para el primer problema se obtuvieron las
gráficas de posición, velocidad y cambio del sistema en función del tiempo, además se verifico la conservación de la energía
en este problema, con los 100 osciladores acoplados se realizó un análisis estadístico de cómo se distribuían las energías y
amplitudes y que tan rápido era esa cambio.

Palabras clave: osciladores acoplados, energía , distribucion.

Abstract
analyzed the simulation of a system of n coupled oscillators, to observe the change of this system in time, for this the
system was examined from 3 oscillators to 100 of these coupled, for the first problem were obtained the graphs of position,
speed and change of the system as a function of time, in addition to the conservation of energy in this problem, with the 100
coupled oscillators, a statistical analysis of how the energies and amplitudes were distributed and how fast that change was
made.

Keywords: coupled oscillators, energy, distribution.

1. INTRODUCCIÓN da una de las masa conformantes del oscilador, a de-

más puede ser caracterizado los grados de participa-
El oscilador armónico simple es uno de los siste- ción de los osciladores acoplados con sus amplitudes.
mas de mayor estudio en la física, ya que la mayoría En el modelado computacional de un sistema de n
de fenómenos que nos rodean y despiertan nuestra de osciladores acoplados realizado en el programa de
curiosidad, oscilan alrededor de un punto de equili- Matlab se puede visualizar como es el comportamien-
brio estable, todos estos sistemas pueden ser modela- to de las masas y el cambio de sus posiciones en el
dos y estudiados en una primera aproximación como tiempo, sus velocidades y encontrar las amplitudes de
un oscilador armónico simple; este tipo de sistemas se cada una de las n masas, y con esta información rea-
caracterizan por su periodo que es el tiempo necesa- lizar un modelado de como es la interacción de cada
rio para un ciclo completo del movimiento, una masa uno de los osciladores acoplados, y como se distribu-
sujeta al extremo de un hilo (péndulo) o a un resor- ye la energía de este sistema.
te(sistema masa resorte), la carga eléctrica almacena-
da en un condensador, las cuerdas de un instrumento
musical, y las moléculas de una red cristalina son al-
2. MARCO TEÓRICO.
gunos de los ejemplos físicos que a menudo realizan
2.1. OSCILADOR ARMÓNICO SIMPLE.
un movimiento oscilatorio.
Otra característica principal de este sistema es la El movimiento armónico simple es un movimien-
conservación de la energía, y como podemos estudiar to periódico de vaivén, en el que un cuerpo oscila de
el aporte de las energías, cinéticas y potenciales de ca- un lado a otro de su posición de equilibrio, en una di-

1
rección determinada y en intervalos iguales de tiem-
Figura 1: oscilador acoplado de dos masas
po. Este sistema se caracteriza por su frecuencia la cual
se define como el recorrido de una posición extrema y
volver a la primera, pasando dos veces por la posición
de equilibrio, se denomina ciclo, el periodo de una
oscilación es el tiempo transcurrido entre dos puntos
equivalentes de la oscilación, y una fuerza restaurado-
ra que tiende a devolverlo al punto de equilibrio esta-
ble, con una intensidad proporcional a la separación
respecto de dicho punto:
→
−
F = −K (→
−
x −−
→)
x 0 (1) d 2 x1
m1 = −K 1 x 1 + K 2 (x 2 − x 1 ) (2)
dt2
donde k es la constante de recuperación, y x 0 es la
posición de equilibrio, que sin pérdida de generalidad d 2 x2
podemos tomar x 0 = 0. Aplicando la segunda ley de m 2 = −K 3 x 2 − K 2 (x 2 − x 1 ) (3)
dt2
Newton, el movimiento armónico simple se define en-
El sistema de ecuaciones que generan 2 y 3 ,aplican-
tonces en una dimensión, mediante la ecuación dife-
do un poco de algebra básica y reacomodando se ge-
rencial:
nera una matriz de frecuencias con un vector de ace-
d 2→
−
x →
−
m = −K x leraciones:
dt2
d2
· ¸ · ¸
La solución de la ecuación diferencial puede escri- x1 x1
=M c
birse en la forma x(t ) = A cos(w t + φ). Lafuerza ejer- dt2 x2 x2
cida por un resorte ideal es conservativa y las fuerzas
verticales no efectúan trabajo, así que se conserva la donde M
ces la matriz de frecuencias.
energía mecánica total del sistema.  
−(k 1 + k 2 ) k2
M

c=  m1 m1 

2.2. OSCILADORES ACOPLADOS  k2 −(k 2 + k 3 ) 
m2 m2
Un sistema de n osciladores acoplados es aquel que
Como se dijo anteriormente se pueden seguir agre-
consta de muchos osciladores individuales interco-
gando masas y resortes al sistema de estudio, para en-
nectados entre sí. Este planteamiento de problema de
contrar como se modela este sistema para n oscilado-
estudio se puede aplicar tanto a sistemas mecánicos
res se aplica el método de inducción matemática, se
como a modelos atómicos de sólidos, además este sis-
observa de qué manera queda el sistema de ecuacio-
tema tiene asociado un conjunto de modos de oscila-
nes para 3 masas como se muestra en la figura 2.
ción con frecuencias características bien definidas, el
modelo de osciladores acoplados mas simple es el for-
mado por dos masas y tres resortes: el primer resorte Figura 2: 3 masas acopladas
con un extremo fijo a una pare y el otro extremo a la
primera masa, a continuación el segundo resorte a la
primera masa y la segunda masa y por último el tercer
resorte unido a la segunda masa por un extremo y a la
pared por el otro extremo, como lo muestra la figura 1
Este tipo de problemas sencillos como el de la figura
1 es sencillo de solucionar, se buscan las dos ecuacio- El sistema de ecuaciones diferenciales en forma ma-
nes diferenciales que modelan el sistema, para realizar tricial es generado de la siguiente manera:
esto tenemos dos caminos uno mediante una análisis
dinámico y otro con la energía potencial de los osci- x x1
   
d2  1 
ladores ambas formas de solucionarlo nos llevaran al x2 =M
c x 2 
mismo sistema de ecuaciones: dt2
x3 x3

2
 La matriz de frecuencias M c que generan estas tres 2.4. CARACTERIZACIÓN.
masas es una matriz de tamaño 3x3, y se puede encon-
Para realizar un mayor análisis del sistema de osci-
trar un patrón de como cada una de las componentes
ladores acoplados y de como se comporta este mien-
va ha estar ubicada; en la diagonal principal encontra-
tras las masas que lo conforman esta en movimiento
remos la suma de las contantes k 1 , k 2 yk 3 .
y como de una en una empiezan a repartir su energía
y el cambio de sus amplitudes, esto se realiza median-
 
−(k 1 + k 2 ) k2 te los grados de participación (R) que están definidos
0

 m1 m1 
 para las amplitudes y energías de cada uno de los osci-
k 2 −(k 2 + k 3 ) k 3 ladores, están definidos por las siguientes ecuaciones:
M
 
c=  

 m2 m2 m2 

 k3 −(k 3 + k 4 )  ( A 2 )2
P
0 R = (6)
m3 m3 P 4
A
Con esto se realiza una generalización del conjun- Donde A es la amplitud del movimiento.
to de estudio de n osciladores acoplados, los cuales se
( E T2 )2
P
describirán como:
R1 = P 4 (7)
    ET
x1 x1
 x   x  Con la energía también podemos encontrar los gra-
 2  2
d2 
 
 x3

c x 3
 =M
 
 dos de participación de sistema de osciladores.
dt2 
 ..

 ..
 
(E T )2 i
 
 .  .
P
 
xn xn R1 = P 2 (8)
ET
La matriz de frecuencia M c generalizada de rango Esta es una medida para este conjunto de datos, ya
nxm se formaría de la siguiente manera. que cada uno de ellos tiene una importancia relativa,
donde i es la posición de cada uno de los osciladores,
 −(k +k )
1 2 k2
 nos permite visualizar la tendencia del sistema.
m1 m1 0 0 ···
k2 −(k 2 +k 3 ) k3
(E T i )2
 
0 ··· 
 P
m2 m2 m2
M =
c 
k3 −(k 3 +k 4 ) k n+1
 R1 = P 2 (9)
0 ···  ET

 m 3 m 3 m n
..
 
k n+1 −(k n +k n+1 )
. 0 m n+1 m n+1
Desplazamiento cuadrático medio con esto se pue-
de dar una idea que tan rápido se da el cambio de la
Al generalizar el problema para n osciladores en- tendencia de nuestro sistema.
contramos m numero de ecuaciones diferenciales de
segundo orden:
3. MODELADO DEL SISTEMA DE OSCILADORES
ACOPLADOS
D 2 |X c〉 = M c|X c〉 (4)

donde X c es un vector columna de n componentes, M c En el software de Matlab se realizo el programa de

2
es una matriz de nxn componentes y D es el opera- simulación para n osciladores acoplados, como se dis-
dor derivada de segundo orden en el tiempo. cutió anteriormente este problema es modelado por
un sistema de ecuaciones diferenciales de segundo or-
den, este fue el primer problema ya que hasta el mo-
2.3. ENERGÍAS DEL SISTEMA.
mento estos programas y sus algoritmos internos no
El conjunto de osciladores acoplados cumple con el resuelven este tipos de ecuaciones, para solucionar es-
principio de conservación de energía, que está dada to se realiza un cambio de variable en las derivadas y
por el cambio de energía cinética a potencial o de for- disminuimos el grado de la ecuación.
ma contraria y que la suma de estas energías es una Para generar la matriz de frecuencia se utilizó la fun-
constante. ción (r and (n, 1)0 ) y ones para generar un vector de
constantes elásticas aleatorias en los intervalos desea-
ET = k + V (5) dos, además, utilizando d i ag se realizo el cambio de

3
las diagonales de la matriz de frecuencias M
c, Por ulti- función del tiempo de 3 osciladores acoplados, tam-
mo se genero un matriz M’ de tamaño (2nx2m) la cual bién es notoria que la máxima amplitud del sistema es
esta compuesta por la matriz de frecuencias, una ma- ±1 y que tiene cierto grado de desorden.
triz identidad y dos matrices de ceros.

0 0 1 0
 Figura 4: v vs t
 0 0 0 1 
M 10 =  k1+k2

k2
 − m1 0 0 

m1
k2
m2 − k1+k2
m2 0 0

en este ejemplo se tiene un sistema de dos oscila-

dores acoplados, y se puede observar que su matriz de
frecuencias M c de tamaño 2x2 y es una componente
de la M 1 que es el doble del tamaño de M
0 c, para ter-
minar las condiciones iniciales del sistema y las reglas
de evolución de este, se produce un arreglo de tiem-
po desde un (t 0 , t f ) y la posición de la masa que será
excitada (y(n) = 1), estas son las condiciones iniciales
del sistema y tienen como objeto ser la semilla para
el algoritmo de integración utilizada para resolver el
problema, en este estudio se empleara un runge kutta
de orden 4, con un paso adaptativo en Matlab este al-
goritmo es llamado od e45 este utilizara la matriz M 10 en la figura 4 se están graficando las velocidades de
una función y las condiciones iniciales y como resul- 3 osciladores acoplados contra el tiempo, esta es muy
tado de este proceso entregara un arreglo de n com- parecida a la gráfica 3, tiene una amplitud de oscila-
ponentes las cuales son tiempo, posición y velocidad ción mayor, pero a diferencia de la anterior es más ar-
y = [t , (x, v)] este arreglo contiene la información de mónica.
la simulación, con este arreglo se puede encontrar la
energía cinética, potencial y total del sistema, además Figura 5: ENERGIA DEL SISTEMA
nos permite graficar las posiciones y velocidades de
cada uno de los osciladores en función del tiempo, y
podemos generar la grafica del espacio de fase.

Figura 3: x vs t

en la figura 5 fue graficada la energía cinética poten-

cial y total, para verificar que el sistema de 3 oscilado-
res conservara la energía total y de esta manera eva-
luar si lo hecho hasta el momento se encontraba de la
forma correcta, como se puede observar la energía ci-
en la gráfica 3 se observa el cambio de la posición en nética y potencial, oscilan como una onda y la suma

4
de estas dos energía total es una línea recta constan- por el ode45 se realiza la gráfica de los grados de par-
te en 1, asegurando el principio de conservación de la ticipación de las amplitudes y las energías del sistema,
energía en el sistema de 3 osciladores acoplados. como se muestra en la figura 7, donde para un t = 0
la amplitud del sistema empieza en 1 lo cual muestra
Figura 6: intercambio de energia de cada OSCILADOR un oscilador excitado y la línea de R1 para este tiempo
tiene una posición a la altura de dos, y se nota como
a medida que le tiempo transcurre esta dos líneas co-
mienza con una pendiente positiva ruidosa, hasta lle-
gar a un máximo este punto nos dice que todos los os-
ciladores se encuentra en movimiento y acaban de re-
botar con las paredes.

Figura 8: M1 VS t

El eje x de la gráfica 6 están las posiciones de las 3

masas y se observa que para un instante t = 0 se tiene
m1 y m3 con posiciones iguales a cero es decir fijas y
m2 se ha dotado con una posición inicial diferente y
posee una energía potencial inicial, como lo muestra
la gráfica 6 ya que en el centro en t = 0 se observa que
tiene un color amarillo, luego transcurre un instante y la figura 8 se observa el cambio de tendencia del sis-
esta comienza a transferir energía a las otras dos masas tema de los 3 osciladores acoplados y como a partir
hasta que estas se colocan totalmente amarillas y m2 de un t = 0 la segunda masa trasmite a las otras ma-
azul en este instante ha transferido toda su energía . sas su energía, hasta llegar a un mínimo donde estas
empiezan a rebotar y transferir de nuevo la energía la
Figura 7: R,R1 VS t masa 2 hasta alcanzar de nuevo un máximo y empe-
zar a decaer de nuevo, esto mostrando una tendencia
de nuestro sistema. la gráfica 9 ilustra el cambio de las

Figura 9: M2 vs t

con las ecuaciones 6 y 7 y con los datos obtenidos

5
partículas de una manera aleatoria y como para un so- lo para sus constantes de (0,9, 1,1), con estas constan-
la repetición de este sistema el número de experimen- tes aleatorias se repitió la simulación 100 generando
tos no es suficiente para mostrar cómo se distribuyen la información suficiente, para encontrar la tendencia
las masas en el tiempo, para realizar un análisis más de distribución de nuestro sistema con respeto a sus
profundo se debe realizar un numero de experimen- amplitudes, Se observa también que las figuras tienen
tos altos y promediarlos de esta manera encontrando todas las mismas tendencia creciente y con pequeñas
la distribución de los osciladores. variaciones o saltos.

Figura 10: R vs t Figura 11: R1 vs t

(a) constante (0.7,1.3) (a) constante (0.7,1.3)

(b) constante (0.8,1.2) (b) constante (0.8,1.2)

(c) constante (0.9,1.1) (c) constante (0.9,1.1)

la sub graficas de la figura 10 se observa el grado de la sub graficas de la figura 11 se observa el grado de
participación de 100 osciladores acoplados con cons- participación de 100 osciladores acoplados con cons-
tantes diferentes en cada uno de sus resortes de una tantes diferentes en cada uno de sus resortes de una
manera aleatoria para la gráfica (a) las constantes de manera aleatoria para la gráfica (a) las constantes de
sus resortes se encontraban en un intervalo (0,7, 1,3), sus resortes se encontraban en un intervalo (0,7, 1,3),
para el grafico de la figura (b) se utilizó el intervalo de para el grafico de la figura (b) se utilizó el intervalo de
(0,8, 1,2) y por último en la gráfica de (c) un interva- (0,8, 1,2) y por último en la gráfica de (c) un interva-

6
lo para sus constantes de (0,9, 1,1), con estas constan- de sus resortes de una manera aleatoria para la gráfi-
tes aleatorias se repitió la simulación 100 generando ca (a) las constantes de sus resortes se encontraban en
la información suficiente, para encontrar la tenden- un intervalo (0,7, 1,3), para el grafico de la figura (b) se
cia de distribución de nuestro sistema con respeto a utilizó el intervalo de (0,8, 1,2) y por último en la gráfi-
sus energías, y como estas se distribuyen en el siste- ca de (c) un intervalo para sus constantes de (0,9, 1,1),
ma de estudio, Se observa también que las figuras tie- con estas constantes aleatorias se repitió la simulación
nen todas las mismas tendencia creciente hasta llegar 100 para realizar un promedio del cambio de este sis-
al máximo donde se estabilizan y después de un tiem- tema pero como se observa en las sub figuras este pro-
po en las primeras dos graficas con más aleatoriedad blema no tiene una tendencia, cada una de las gráficas
comienzan a decaer, también tienen pequeñas varia- difiere y el único punto en común es en t = 0 donde se
ciones o saltos. inicia la perturbación del sistema.

Figura 12: M1 vs t Figura 13: M2 vs t

(a) constante (0.7,1.3) (a) constante (0.7,1.3)

(b) constante (0.8,1.2) (b) constante (0.8,1.2)

(c) constante (0.9,1.1) (c) constante (0.9,1.1)

la sub graficas de la figura 12 se observa como es el la sub graficas de la figura 13 se observa que tan rá-
cambio medio en la distribución del los 100 oscilado- pido y aleatorio es el cambio medio en la distribución
res acoplados con constantes diferentes en cada uno del los 100 osciladores acoplados con constantes dife-

7
rentes en cada uno de sus resortes de una manera alea-
toria para la gráfica (a) las constantes de sus resortes se
encontraban en un intervalo (0,7, 1,3), para el grafico
de la figura (b) se utilizó el intervalo de (0,8, 1,2) y por
último en la gráfica de (c) un intervalo para sus cons-
tantes de (0,9, 1,1), con estas constantes aleatorias se
repitió la simulación 100, se puede observar que aun-
que la tendencia de cambio de los sistemas son dife-
rentes la velocidad y distribución en la que cambian
son muy parecida ya que en las tres sub graficas de la
figura 13 tienen tendencias crecientes con pequeños
saltos o variaciones.

4. CONCLUSIÓN
Mediante la herramienta computacional de matlab
se realizó la simulación de un sistema de n oscilado-
res acoplados, con los cuales se generó una caracte-
rización de su amplitud y se verifico la conversación
de sus energías, además se encontró la forma cambio
en la tendencia de distribución de este sistema, esto
mediante a la repetición de la simulación 100 veces,
obteniendo el desplazamiento cuadrático medio y la
tendencia de cambio con diferentes y aleatorias cons-
tantes elásticas en unos intervalos definidos, también
se examinó como se repartía la energía en estos siste-
mas es mediante el grado de participación, con lo cual
se llega a la conclusión de que distribuyen su anergia
de una manera parecida sin importar el cambio en sus
constantes.

5. Referencias
[ 1 ] fisica conceptual 10 edicion.págs. 22,112-114,

[ 2 ] Physical Review A 91, 033815 (2015).