首先,设定一个天体互绕系统:
$$ d_1 = \frac{m_2 a}{m_1 + m_2}, d_2 = \frac{m_1 a}{m_1 + m_2} $$
$$ \omega = \sqrt{\frac{G(m_1 + m_2)}{a^3}} $$
然后,求解拉格朗日点的位置。
我们以这个旋转参考系为参照物。那么检验质点的势包括三部分:
三者可以线性叠加。那么当 $m$ 位于 $(x,y)$ 的时候,其总势能:
$$ E(x,y) = -\dfrac{G m_1}{\sqrt{(x+d_1)^2+y^2}}-\dfrac{Gm_2}{\sqrt{(x-d_2)^2+y^2}}-\dfrac{G(m_1+m_2)}{2a^3}(x^2+y^2) $$
平衡条件为势能在各个方向的一阶偏导数为 0,即:
$$ \begin{cases} \dfrac{\partial E}{\partial x}=\dfrac{Gm_1\left(x+\dfrac{m_2 a}{m_1+m_2}\right)}{\left(\left(x+\dfrac{m_2 a}{m_1+m_2}\right)^2+y^2\right)^{3/2}}+\dfrac{Gm_2\left(x-\dfrac{m_1 a}{m_1 + m_2}\right)}{\left(\left(x-\dfrac{m_1 a}{m_1+m_2}\right)^2+y^2\right)^{3/2}}-\dfrac{G(m_1+m_2)}{a^3}x = 0 \newline \dfrac{\partial E}{\partial y}=\left[\dfrac{Gm_1}{\left(\left(x+\dfrac{m_2 a}{m_1+m_2}\right)^2+y^2\right)^{3/2}}+\dfrac{Gm_2}{\left(\left(x-\dfrac{m_1 a}{m_1+m_2}\right)^2+y^2\right)^{3/2}}-\dfrac{G(m_1+m_2)}{a^3}\right]y=0 \end{cases} $$
根据第二条,满足 $y$ 方向为 0 有两种可能: