Search This Blog

Thursday, July 22, 2010

Oscilador Armónico

Robert Hooke fue un enemigo de Newton. Según la nota anterior de Wikipedia que pongo en este párrafo, no desarrollaba sus ideas, aunque tenía muchas. Se pelearon por la autoría de la Ley Universal de la Gravitación que será tema de otra nota. Ahora me interesa la llamada ``Ley de Hooke.''

Lo bueno de empezar una área del conocimiento es que el nombre de uno aparece en cuestiones fundamentales. En este caso, de la segunda liga de Wikipedia del párrafo anterior podemos leer:

``La forma más común de representar matemáticamente la Ley de Hooke es mediante la ecuación del muelle o resorte, donde se relaciona la fuerza F ejercida sobre el resorte con la elongación o alargamiento [;\delta;] producido:
F = -k\delta \, ''
Aquí doy la cinemática, y más adelante en dinámica, veremos ¿qué representa la fuerza, y porqué esta relación se satisface en muchos casos?

Cinemáticamente, es decir matemáticamente, nos interesa conocer la posición de un cuerpo, cuya aceleración es conocida, no en función del tiempo, sino de la posición misma. Sucede que este problema se puede resolver exactamente con diferencias finitas, y por eso lo pongo aquí.

Del problema de aceleración constante presentado antes, vimos que si la velocidad varía linealmente se puede sustituir esa variación por la velocidad promedio.

[;\bar v_n = \frac{v_{n+1} + v_n}{2};]

Poniendo atención en ésto; vemos que es la matemática de la variación la que nos lleva a una situación en que la sustitución de una constante en lugar de una variable no produce errores de cálculo. En pocas palabras, lo que corre de menos Aquiles al principio de la carrera, exactamente lo corre de más al final.

Según Hooke, la aceleración varía linealmente con la posición. Con esta observación empírica, obtenemos el siguiente problema en una dimensión.

En lugar de:

[;a_n = -k x_n;]

se puede obtener una solución exacta con:

[;a_n = -k \bar x_n;]

Donde:

[;\bar x_n = \frac{x_{n+1} + x_n}{2};]

Con la definición cinemática de aceleración, y la ecuación anterior, podemos plantear y resolver el problema del oscilador armónico exactamente en diferencias finitas.

[;\bar v_n = \frac{x_{n+1} - x_n}{t_{n+1} - t_n};]

Queda un sistema lineal de dos ecuaciones y dos incógnitas recursivo. Usaremos matrices de dos dimensiones para despejar la velocidad y la posición en el tiempo [;(n+1)\epsilon;], en función de la velocidad y la posición en [;n\epsilon;], en la siguiente nota.

Para no complicar la presentación tomo el caso de [;k = 1;], y [;m = 1;]. Así tenemos:

[;\frac{v_{n+1} + v_n}{2}\times \epsilon = x_{n+1} - x_n;]

[;v_{n+1} - v_n = - \frac{x_{n+1} + x_n}{2}\times \epsilon;]

No comments:

Post a Comment