通常来说，选取惯性系为参考系，列出的三体问题方程为
$\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}$

历史上出现过很多不同形式的变换，使得三体问题的运动方程有了各样的形式，如Lagrange形式、Jacobi形式、Hamilton形式等。这些变换形式都各有特点，都能够在一定程度上化简三体问题。BoJone在研究摆弄等质量型三体问题的运动方程时，也发现了一种很有趣的变换，在此贴出与大家分享。

设$\vec{R}_1=\vec{r}_1-\vec{r}_2,\vec{R}_2=\vec{r}_2-\vec{r}_3,\vec{R}_3=\vec{r}_3-\vec{r}_1$，则三体问题的运动方程变为

阅读剩余部分...

在上一些关于限制性三体问题的探讨中，我们得出了在平面上的方程：
$\ddot{R}+2i\omega \dot{R}=\omega^2 R-GM\frac{R-l_1}{|R-l_1|^3}-Gm\frac{R-l_2}{|R-l_2|^3}$————(32)

能量积分为：
$\frac{1}{2}|\dot{R}|^2=\frac{1}{2} \omega^2 |R|^2+\frac{GM}{|R-l_1|}+\frac{Gm}{|R-l_2|}-C$————(33)

下面就以这两个方程为基础，再说说限制性三体问题的那些事儿...

阅读剩余部分...

平面圆形限制性三体问题运动方程及能量积分
plane circular restricted three-body problem
02.04有重要修正！！

寒假一个很大的目标就是能够在三体问题的周期轨道上有点突破，于是就出动了“向量”、“复分析”、“微分方程”等理论“核武”，遗憾的是，“有心栽花花不开”，到今天还是没有多少进展。不过俗语也说“无心插柳柳成荫”，也不错。今天回看《天体力学引论》中的“圆形限制性三体问题”，经过一番思考，利用这些天的思考方法重新推导出了其运动方程和能量积分，也算是“意外收获”在此作为春节礼物与大家分享。

所谓“圆形限制性三体问题”，就是指两个大质量天体（质点）在它们相互引力作用下做圆周运动，假设第三天体（质量趋于0）只受到这两个天体的引力作用而不影响两个天体运行的一种运动情况。由于普通三体问题无法积分，而这个“限制性模型”能够把问题化简不少（不过还是不能积分出来的），因此也得到了一定应用。它的应用条件是：第三体质量小（如当前航天器与地球、太阳）、短程。注意短程也是相当重要的条件之一，注意短程也是相当重要的条件之一，质量越小应用范围越大。要是质量大的话，就不能计算太长的路程。

阅读剩余部分...

The Three Body Problem and its Classical Integration

很多天文爱好者都已经接触到了“二体问题”（我们在高中学习到的“开普勒三定律”就是内容之一），由于在太阳系中行星质量相对较小而且距离相对较远，应用“二体问题”的解对天体进行计算、预报等能够满足一定的近似需求。不过，如果需要更高精度的计算，就不能把其他行星的引力给忽略掉了，于是就产生了所谓N体问题（N-Body Problem），即N个质点尽在它们各自引力的相互作用下的运动规律问题。最简单的二体已经被彻底解决，而三体或更多体的问题则与二体大相径庭，因为庞加莱证明了，三体问题不能严格求解，而且这是一个混沌系统，任何微小的扰动都会造成不可预期的效果。

根据牛顿力学，选择惯性参考系，设三个质点分别为$M_1,M_2,M_3$，向径分别为$\vec{r_1},\vec{r_2},\vec{r_3}$，可以列出运动方程（以下的导数都默认是对时间t求导）

阅读剩余部分...

2010已经成为历史了，在2011的第一天，BoJone祝大家新年快乐，生活、学习、工作都更上一层楼！我愿一直与大家探讨科学，分享科学！

一直想好好地总结一下过去的一年内的事情，无奈事情太多，一拖再拖。其实在2010年里，最值得纪念的当然就是完完整整地经历了一次天文竞赛。从3月的预选，到五月的宁夏固原决赛，接着是7月的北京集训，最后是9月下旬的北京IOAA。一步步走来的足迹，浮现在脑海，历历在目。

阅读剩余部分...

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)

现在我们只考虑三体问题。天文学家一直希望能够找到三体问题的简洁解，可是很遗憾，庞加莱已经证明了三体问题的解是混沌的，也就是说任何微小的扰动都有可能造成不可预料的后果（可以形象的比喻为：巴西的一只蝴蝶翅膀的扇动，有可能因此美国的一场龙卷风）。

阅读剩余部分...

不知道各位读者还记得BoJone在《方程与宇宙》这一章中写了整整三篇文章来学习天体力学中的二体问题吗？虽然对二体问题基本上做了一个描述，但是依旧是冰山一角。而在最近写的几篇文章中，BoJone又强调了“向量”的巨大作用。那么，当天体力学与向量碰头后，会发生什么大事呢？难道，火星撞上了地球？

阅读剩余部分...