tag:blogger.com,1999:blog-11464109.post-1143566198555621082006-03-28T20:15:00.000+02:002006-03-28T20:18:04.306+02:00El problema de los tres cuerpos<p>Ojeando la última revista del Círculo de Lectores, me llamó la atención un libro de título <q><a href="http://www.circulo.es/Contenido/Libros/Libro.asp?Codigo=44222">La incógnita Newton</a></q>, y leí el resumen para comprobar si se trataba de otro libro más que se suma a la moda iniciada por El Codigo DaVinci. Decía así:</p><blockquote cite="http://www.circulo.es/Contenido/Libros/Libro.asp?Codigo=44222"><p>Cambridge, año 1888. El mundo de la joven maestra Vanessa Duncan sufre una sacudida cuando un profesor de matemáticas de la prestigiosa universidad es hallado muerto de un golpe en la nuca. El profesor Akers se encontraba trabajando en la resolución del <q>problema de los tres cuerpos</q>, un enigma matemático que sir Isaac Newton fue el primero en plantear</p></blockquote><p>Y resulta que el problema de los tres cuerpos no es ningún <q>enigma</q>. Es un problema sin solución analítica, que es algo muy distinto. ¿El qué? Veamos primero en qué consiste el famoso problema.</p><p> En el cole nos enseñaron la famosa <a href="http://es.wikipedia.org/wiki/Gravedad">Ley de Gravitación Universal</a> de Newton, cuya fórmula nos indicaba la fuerza gravitatoria existente entre dos cuerpos: <span style="font-weight: bold;">F = G·M·m/r<sup>2</sup></span>, donde <span style="font-weight: bold;">F</span> es la fuerza, <span style="font-weight: bold;">M</span> la masa de uno de los cuerpos, <span style="font-weight: bold;">m</span> la masa del otro cuerpo, <span style="font-weight: bold;">r</span> la distancia que los separa, y <span style="font-weight: bold;">G</span> la constante de gravitación universal. También nos enseñaron que la fuerza aplicada sobre un cuerpo, producía en este una aceleración, dada por la fórmula <span style="font-weight: bold;">a=F/m</span> (<span style="font-weight: bold;">a</span> es la aceleración, <span style="font-weight: bold;">F</span> la fuerza y <span style="font-weight: bold;">m</span> la masa). Con estas dos fórmulas, pero expresadas en forma vectorial (la fuerza y la aceleración tienen dirección), es bastante sencillo calcular la trayectoria de un objeto en órbita, conociendo su posición y velocidad en un instante dado.</p><p>Pero para ello siempre consideramos que uno de los objetos tiene una masa muy pequeña con respecto al otro. Esto nos vale para satélites artificiales y vehículos espaciales, pero no para sistemas en los que ambos cuerpos tienen masas no demasiado distintas, como una <a href="http://es.wikipedia.org/wiki/Estrella_binaria">estrella binaria</a> (dos estrellas orbitando una alrededor de la otra). En este último caso, la masa del segundo cuerpo afectará a la trayectoria del primero. Es el llamado <a href="http://en.wikipedia.org/wiki/Two-body_problem">problema de los dos cuerpos</a>, y aunque es más complicado, se pueden obtener dos funciones que nos indiquen la posición de cada cuerpo, en función del tiempo.</p><p>¿Y qué pasa si en vez de dos cuerpos, tenemos tres? Pues como todos imaginaréis, estamos ante el <a href="http://en.wikipedia.org/wiki/Three-body_problem">problema de los tres cuerpos</a>. En este caso, no existe solución analítica, es decir, no podemos encontrar una función para cada cuerpo que nos de su posición en función del tiempo. Pero el que no tenga solución analítica no quiere decir que sea un enigma. Mediante análisis numérico se pueden realizar cálculos suficientemente aproximados sobre la trayectoria de los cuerpos.</p><p>¿Solución analítica? ¿Análisis numérico? ¿Qué es todo esto? Veamos un ejemplo sencillo:</p><p> Supongamos que queremos saber la posición de un cuerpo en movimiento rectilíneo y uniformemente acelerado, esto es, con aceleración constante. Conocemos la distancia inicial <span style="font-weight: bold;">S<sub>0</sub></span> a nuestro origen de coordenadas, la velocidad inicial <span style="font-weight: bold;">V<sub>0</sub></span>, y la aceleración <span style="font-weight: bold;">a</span>. En el colegio nos enseñaron que la distancia del cuerpo al origen de coordenadas en función del tiempo es <span style="font-weight: bold;">S(t) = S<sub>0</sub> + V<sub>0</sub> · t + (1/2) · a · t<sup>2</sup></span>. Esto es una solución analítica. Tenemos una función que nos permite calcular la distancia <span style="font-weight: bold;">S</span> para cualquier valor de <span style="font-weight: bold;">t</span> (tanto positivo como negativo, es decir, tanto en el futuro como en el pasado).<br /></p><p>Imaginemos ahora que desconocemos dicha fórmula, y no sabemos calcularla. ¿Qué haríamos? Bueno, sabemos que la distancia inicial es <span style="font-weight: bold;">S<sub>0</sub></span>, por lo que en nuestros cálculos deberemos partir de ahí. Sabemos también que la distancia recorrida por un cuerpo a velocidad constante es <span style="font-weight: bold;">v·t</span>. Sabemos incluso que la velocidad varía de la forma <span style="font-weight: bold;">v=a·t</span>. Pero tenemos el problema de que no sabemos calcular la distancia recorrida en esas circunstancias. Bien, supongamos que calculamos la distancia recorrida en un intervalo de tiempo muy pequeño, de forma que la velocidad apenas varía. Entonces, podemos calcular esa pequeña distancia como <span style="font-weight: bold;">v·t</span>. Ya, pero yo quiero calcularla para un intervalo de tiempo más grande. Pues dividimos ese intervalo en otros más pequeños, calculamos la distancia recorrida en cada intervalo (cada uno con velocidad distinta) y lo sumamos todo. El resultado será aproximado, y esa aproximación será mejor cuanto más pequeño sea el intervalo en el que suponemos que la velocidad es constante. Pues bien, a grandes rasgos (y que me perdonen los matemáticos por la simplificación), eso sería el <a href="http://es.wikipedia.org/wiki/An%C3%A1lisis_num%C3%A9rico">análisis numérico</a>.</p><p>El problema de los tres cuerpos no tiene solución analítica, es decir, no podemos obtener una función <span style="font-weight: bold;">S(t)</span> que nos diga la posición del objeto en cualquier instante de tiempo. Pero no es un problema irresoluble, ya que podemos atacarlo mediante análisis numérico. Y para eso contamos en la actualidad con una herramienta extraordinaria: el ordenador. Efectivamente, para realizar este tipo de cálculos tenemos que realizar muchas veces la misma operación, pero con distintos valores. Y esto es precisamente lo que mejor se le da a un ordenador (de hecho, podríamos decir que es una de las pocas cosas que sabe hacer).</p><p>Así que el problema de los tres cuerpos dista mucho de ser un <q>enigma matemático</q>. Cualquier PC actual tiene potencia más que suficiente para realizar los cálculos necesarios en milésimas de segundo, y pintar una animación en tiempo real, o incluso acelerada, de las trayectorias de los cuerpos. Un ejemplo de ello lo tenemos en <a href="http://faculty.ifmo.ru/butikov/Projects/Collection.html">esta web</a>, con varios <a style="font-style: italic;" href="http://es.wikipedia.org/wiki/Java_applet">applets</a> que nos pintan diversas trayectorias (hay que tener instalado el <a style="font-style: italic;" href="http://es.wikipedia.org/wiki/Plugin">plugin</a> de <a href="http://es.wikipedia.org/wiki/JRE">Java</a> en el navegador).</p><p>Para no cargar toda al culpa sobre el Círculo de Lectores, hay que decir que el mismo error está presente en otros sitios, como en <a href="http://www.casadellibro.com/fichas/fichabiblio/0,1094,2900001062870,00.html?codigo=2900001062870&titulo=LA+INCOGNITA+NEWTON">La Casa del Libro</a>. Supongo que el error original estará en la propia editorial del libro.<br /></p><div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/11464109-114356619855562108?l=www.malaciencia.info'/></div>Alfhttp://www.blogger.com/profile/00534038511753640060alf@malaciencia.info20