感谢国家天文台LAMOST项目之“宇宙驿站”提供网络空间和数据库资源! 感谢国家天文台崔辰州博士等人的多方努力和技术支持！
科学空间欢迎您转载本站文章，但在转载本站原创文章时，希望您能够尊重版权，注明来自科学空间，谢谢！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)
现在我们只考虑三体问题。天文学家一直希望能够找到三体问题的简洁解，可是很遗憾，庞加莱已经证明了三体问题的解是混沌的，也就是说任何微小的扰动都有可能造成不可预料的后果（可以形象的比喻为：巴西的一只蝴蝶翅膀的扇动，有可能因此美国的一场龙卷风）。
对于解方程，代数学家希望能够从理论上证明解的存在性以及解的求法，所以就有了1到4次方程的求根公式、5次及以上的代数方程没有根式可解等重要理论；然而，通常的学者（如物理学家、天文学家）都不需要这些内容，他们只关心如何尽可能快地求出指定方程的根（尤其是实数根），所以他们通常关注的是方程的数值算法，当然，如果能有一个相对简单的求根公式，也是他们所希望的。而接下来所要介绍的内容，则是满足了这一需要的三次方程的求根公式，其中用到的相当一部分的理论，是与三角函数相关的。
知识储备：
$\frac{2}{tan 2A}=\frac{1}{tan A}-tanA$————(1)
$\frac{2}{sin 2A}=\frac{1}{tan A}+tanA$————(2)
$cos(3A)=4cos^3 A-3cosA$————(3)
不知道各位读者还记得BoJone在《方程与宇宙》这一章中写了整整三篇文章来学习天体力学中的二体问题吗？虽然对二体问题基本上做了一个描述，但是依旧是冰山一角。而在最近写的几篇文章中，BoJone又强调了“向量”的巨大作用。那么，当天体力学与向量碰头后，会发生什么大事呢？难道，火星撞上了地球？
呵呵，做了一回标题党，可能说得夸张了一点。说是“终极算法”，主要是因为它可以任意提高精度、而且几乎可以应付任何非线性方程（至少理论上是这样），提高精度是已知的迭代式上添加一些项，而不是完全改变迭代式的形式，当然在提高精度的同时，计算量也会随之增大。其理论基础依旧是泰勒级数。
我们考虑方程$x=f(y)$，已知y求x是很容易的，但是已知x求y并不容易。我们考虑把y在$(x_0,y_0)$处展开成x的的泰勒级数。关键是求出y的n阶导数$\frac{d^n y}{dx^n}$。我们记$f^{(n)}(y)=\frac{d^n x}{dy^n}$，并且有
$\frac{dy}{dx}=\frac{1}{(\frac{dx}{dy})}=f'(y)^{-1}$
经过上两回的讨论，我们已经基本摸清了二体问题的运动情况。我们已经找到了二体问题在轨道为椭圆的时候的所有积分，给出了“活力公式”等常用公式的证明，并且留下了一些没有解答的问题。那就是在轨道为抛物线和双曲线时的最后一个积分还没有找出来，现在我们解决这两个问题。其中的关键积分依旧是
$\dot{r}^2={2\mu}/r-{\mu a(1-e^2)}/r^2-\frac{\mu}{a}$——(12)
在上一回的讨论中，我们已经解决了大部分的问题，并且表达了找到r或者$\theta$关于时间t的函数的希望。在最后的内容中，我们做了以下工作：
由(7)得到$\dot{\theta}=h/r^2$，代入(6)得到：
$\ddot{r} -h^2/r^3=-\frac{\mu}{r^2}$——(10)
这是一个二阶微分方程，它的解很容易找出，但是这个积分太复杂：
$\dot{r}\frac{d\dot{r}}{dr}=h^2/r^3-\frac{\mu}{r^2}$
$\dot{r}d\dot{r}=(h^2/r^3-\frac{\mu}{r^2})dr$，两端积分
$\dot{r}^2={2\mu}/r-h^2/r^2+K_1$——(11)
$=>{dt}/{dr}=\frac{r}{\sqrt{K_1 r^2+2\mu r-h^2}}$
$t=\int \frac{r}{\sqrt{K_1 r^2+2\mu r-h^2}}dr$
最近评论