感谢国家天文台LAMOST项目之“宇宙驿站”提供网络空间和数据库资源! 感谢国家天文台崔辰州博士等人的多方努力和技术支持！
科学空间欢迎您转载本站文章，但在转载本站原创文章时，希望您能够尊重版权，注明来自科学空间，谢谢！（本文已被刊登在2012年1月的《天文爱好者》上，于笔者而言这是一份很棒的新年礼物！）
在去年第七期《天爱》上，我们看到了N体问题所呈现出来的一些对称、漂亮的周期轨道，这体现了N体问题和谐有序的一面。但是这仅仅是N体问题的冰山一角，笔者也提到过N体问题的本质是混沌、无序的，通俗来讲就是非常乱，无法用数学方程来精确描述。这看起来是一种不完美。但试想，探索当初伽利略将望远镜对准月球后，看到的是如想象中光滑的月面，那么他还会惊叹宇宙的神奇吗？
本文就让我们来更深入地了解一下N体问题的研究历史。
观测&拟合时代
由于人类的自我优越感以及日月星辰东升西落的经验，让我们长期都认为地球是宇宙的中心。第一个比较系统提出地心说的人当属天文学家欧多克斯（Eudoxus，死于公元前347年左右），但他的地心说是非常粗糙的，以至于无法解释很多基本现象，如无法准确预言日食和解释行星逆行等。但亚里士多德接受了地心说，并且由于他在政治和科学上的权威，使地心说免去了夭折的命运。后来托勒密通过他的本轮，完善了地心说，使之延续到了16世纪。
转载自：http://www.matrix67.com/blog/archives/3979
源文件来自这里，提供了多体问题中颇具代表性的 47 个解的数据（本文附带下载），是gnuplot 格式。Matrix67选择了其中 30 个，用 Mathematica 读出数据，生成了 30 个直观的 gif 动画。大家将会看到，在引力的作用下，多颗星体可能会形成的一些极其诡异的轨道。后面的解越来越不平凡，可见多体问题之难。
The New Calculation Of Lagrangian Point 4,5
上一回我们已经求出了拉格朗日点L1,L2,L3，并且希望能够求出L4,L5两个点。由于L4,L5与“地球-太阳”连线已经不共线了，所以前边的方法貌似不能够用了。为了得到一个通用的定义，我们可以采用以下方法来描述拉格朗日点：位于拉格朗日点的天体，它与太阳的连线以及地球与太阳的连线所组成的角的大小是恒定的。（这里为了方便，采用了地日系的拉格朗日点来描述，对于一般的三体问题是一样的）
对于L4,L5来说，我们或许可以设置一个新的向量来描述这两点的向径（如$\vec{R}$）。当我们这样做后，很快就会发现这样会令我们的计算走向死胡同。因为我们发现：已知两个向量的夹角和其中一个向量，我们很难把另一个向量用已知向量的式子表达出来。不能做到这一点，就不能找出$\vec{R}$与$\vec{r}$的关系，就无法联立方程求解。难道，我们这一条路走到尽头了吗？一开始BoJone也冥思苦想不得头绪，但是...
The New Calculation Of Lagrangian Point 1,2,3
关于n体问题，选择质心或其他定点为参考点，我们可以列出下面的运动方程：
$\ddot{\vec{r}}_k=\sum_{i=1,i != k}^{n} Gm_i\frac{\vec{r}_i-\vec{r}_k}{|\vec{r}_i-\vec{r}_k|^3}$————(19)
现在我们只考虑三体问题。天文学家一直希望能够找到三体问题的简洁解，可是很遗憾，庞加莱已经证明了三体问题的解是混沌的，也就是说任何微小的扰动都有可能造成不可预料的后果（可以形象的比喻为：巴西的一只蝴蝶翅膀的扇动，有可能因此美国的一场龙卷风）。
这两个程序都是为了天体力学服务的...
最近评论