为什么行星的运行轨道都是椭圆的?
2024-11-19 阅读 114
更新于 2024年11月21日
似乎没有计算的，我来抛砖引玉一下吧
先说结论：经典框架下，若恒星质量足够大，则行星的运行轨迹为圆锥曲线。
一、运动方程不妨设恒星位于原点，
行星运动所对应的轨迹向量为 \vec s(t)=(x(t),y(t))\vec s(t)=(x(t),y(t)) ,
求出其速度向量为 \vec{v}(t)=(x'(t),y'(t))\vec{v}(t)=(x'(t),y'(t)) ,
而加速度向量显然为 \vec{a}=(x'',y'')\vec{a}=(x'',y'')
由牛顿定理 F=ma,F=\frac{\color{blue}{GM}m}{\color{green}{r^2}}F=ma,F=\frac{\color{blue}{GM}m}{\color{green}{r^2}} 不难得到矢量等式
\vec{a}=\frac{\color{blue}k}{\color{green}{\left| s^2 \right|}}(-\frac{\vec s}{\left| s \right|})\\\vec{a}=\frac{\color{blue}k}{\color{green}{\left| s^2 \right|}}(-\frac{\vec s}{\left| s \right|})\\ 两个分量分别相等可以得到
x''=-k\frac{x}{(x^2+y^2)^{\frac{3}{2}}}\\ y''=-k\frac{y}{(x^2+y^2)^{\frac{3}{2}}}\\x''=-k\frac{x}{(x^2+y^2)^{\frac{3}{2}}}\\ y''=-k\frac{y}{(x^2+y^2)^{\frac{3}{2}}}\\ 两式相加、相除，整理一下就是
\begin{align} x''^2+y''^2&=k^2\frac{1}{(x^2+y^2)^2}&(1)\\ x''y&=xy''&(2)\\ \end{align}\\\begin{align} x''^2+y''^2&=k^2\frac{1}{(x^2+y^2)^2}&(1)\\ x''y&=xy''&(2)\\ \end{align}\\ 为了方便解这个微分方程，进行换元（这里人傻了，应该一开始就用极坐标系求解的）
x(t)=r(t)\cos(\theta(t))\\ y(t)=r(t)\sin(\theta(t))\\x(t)=r(t)\cos(\theta(t))\\ y(t)=r(t)\sin(\theta(t))\\ 对两式求二阶导
x''=r''\cos\theta-2r'\sin\theta \cdot\theta'-r\cos\theta \cdot\theta'^2-r\sin\theta\cdot\theta''\\ y''=r''\sin\theta+2r'\cos\theta \cdot\theta'-r\sin\theta \cdot\theta'^2+r\cos\theta\cdot\theta''\\x''=r''\cos\theta-2r'\sin\theta \cdot\theta'-r\cos\theta \cdot\theta'^2-r\sin\theta\cdot\theta''\\ y''=r''\sin\theta+2r'\cos\theta \cdot\theta'-r\sin\theta \cdot\theta'^2+r\cos\theta\cdot\theta''\\ 代入 (2)(2) 化简得到
2r'\theta'+r\theta''=0\\2r'\theta'+r\theta''=0\\ 利用这个式子化简微分
x''=(r''-r\theta'^2)\cos \theta\\ y''=(r''-r\theta'^2)\sin \theta\\x''=(r''-r\theta'^2)\cos \theta\\ y''=(r''-r\theta'^2)\sin \theta\\ 代入 (1)(1)
(r''-r\theta'^2)^2=\frac{k^2}{r^4}\\(r''-r\theta'^2)^2=\frac{k^2}{r^4}\\ 最后得到了一个微分方程组
\begin{cases} 2r'\theta'+r\theta''=0\\ r''-r\theta'^2=\frac{k}{r^2}\\ \end{cases}\\\begin{cases} 2r'\theta'+r\theta''=0\\ r''-r\theta'^2=\frac{k}{r^2}\\ \end{cases}\\
二、方程求解 第一个解出两者的关系
\begin{align} 2r'r\theta'+r^2\theta''&=0\\ \frac{d(r^2\theta')}{dt}&=0\\ r^2\theta'&=C \end{align}\\\begin{align} 2r'r\theta'+r^2\theta''&=0\\ \frac{d(r^2\theta')}{dt}&=0\\ r^2\theta'&=C \end{align}\\ 这个其实就是著名的开普勒第二定律，可以写作
\frac{1}{2}r^2d\theta=Cdt\\\frac{1}{2}r^2d\theta=Cdt\\
即相等时间间隔，行星与太阳的连线扫过的面积相等
solve problem NO.1,now to problem NO.2
有了关系一之后，可以将原表达式转换为不含 tt 的形式
视 rr 为 \theta\theta 的函数，并令 u=\frac{1}{r}u=\frac{1}{r}
\begin{align} r'&=\frac{dr}{d\theta}\frac{d\theta}{dt}=\frac{dr}{d\theta}\frac{C}{r^2}\\ &=-C\frac{du}{d\theta}\\ r''&=-C\frac{d^2u}{d\theta^2}\frac{d\theta}{dt}=-C^2\frac{d^2u}{d\theta^2}u^2 \end{align}\\\begin{align} r'&=\frac{dr}{d\theta}\frac{d\theta}{dt}=\frac{dr}{d\theta}\frac{C}{r^2}\\ &=-C\frac{du}{d\theta}\\ r''&=-C\frac{d^2u}{d\theta^2}\frac{d\theta}{dt}=-C^2\frac{d^2u}{d\theta^2}u^2 \end{align}\\ 代入第二个等式，
\begin{align} -C^2\frac{d^2u}{d\theta^2}u^2-C^2u^3&=ku^2\\ \frac{d^2u}{d\theta^2}+u&=A \end{align}\\\begin{align} -C^2\frac{d^2u}{d\theta^2}u^2-C^2u^3&=ku^2\\ \frac{d^2u}{d\theta^2}+u&=A \end{align}\\ 得到一个非常清爽的二阶常系数常微分方程
通解显然为
u=A+B\cos(\theta-\theta_0)\\u=A+B\cos(\theta-\theta_0)\\ 即
r=\frac{p}{1-e\cos(\theta-\theta_0)}\\r=\frac{p}{1-e\cos(\theta-\theta_0)}\\ 这是极坐标下的圆锥曲线一般方程,但是仍需要一些办法来求解未知的参数
三、初值求解根据能量守恒，引力势能和行星动能守恒，有
\begin{align} \frac{1}{2}mv_2^2-\frac{1}{2}mv_1^2=-\int_{r_1}^{r_2}\frac{GMm}{r^2}\\ v_2^2-v_1^2=2GM(\frac{1}{r_2}-\frac{1}{r_1}) \end{align}\\\begin{align} \frac{1}{2}mv_2^2-\frac{1}{2}mv_1^2=-\int_{r_1}^{r_2}\frac{GMm}{r^2}\\ v_2^2-v_1^2=2GM(\frac{1}{r_2}-\frac{1}{r_1}) \end{align}\\
另外，还有行星的曲率半径关系，
F=\frac{GMm}{r^2}=m\frac{v^2}{R}\\F=\frac{GMm}{r^2}=m\frac{v^2}{R}\\ 以下记 r'=\frac{dr}{d\theta}r'=\frac{dr}{d\theta} ,那么曲率半径就是
R=\frac{(r^2+r'^2)^{\frac{3}{2}}}{\left| r^2+2r'^2-rr'' \right|}=\frac{v^2R^2}{GM}\\R=\frac{(r^2+r'^2)^{\frac{3}{2}}}{\left| r^2+2r'^2-rr'' \right|}=\frac{v^2R^2}{GM}\\ 利用这些关系式，再加上 t_0t_0 时的初值 (r_0,\theta_0),v_0,GM(r_0,\theta_0),v_0,GM ，可以解出相关的参数
而如果行星轨道为椭圆，有
v=\sqrt{GM(\frac{2}{r}-\frac{1}{a})}\\v=\sqrt{GM(\frac{2}{r}-\frac{1}{a})}\\ 其中， aa 是半长轴，利用这个式子确定初值会方便一些。
附录1.圆锥曲线极坐标
简单添加几句关于圆锥曲线极坐标统一方程的问题
圆锥曲线的统一定义：到定点 FF 的距离与到此点外定直线 ll 的距离之比为 ee 的点的集合
其中， 01e>1 为双曲线。
在极坐标下，若 FF 位于原点，若设点 FF 到直线 ll 的距离为 pp ,则圆锥曲线的方程可以统一写作
\rho=\frac{ep}{1-e\cos \theta}\\\rho=\frac{ep}{1-e\cos \theta}\\ 2.水星进动等问题
本文主要讲的是经典物理体系下，大质量二体运动模型的解，即只可能是圆锥曲线（带退化），而行星轨道为什么都是椭圆的，主要是由于恒星质量往往非常大，且不同轨道之间的行星相互影响很小，所以可以近似认为是二体运动，而行星的轨道要求封闭，而圆锥曲线中只有椭圆是封闭轨道，因此行星的运动轨道可以认为是椭圆。
然而，很早以前，人们就观测到了水星的进动现象，
可以简单认为是水星运动分解为绕太阳的椭圆运动与椭圆轴绕太阳旋转的叠加，这得等到广义相对论出来之后才能解释了。当年爱因斯坦在发表狭义相对论后不久就曾研究过这个问题，答主曾经读过那篇大胆的论文，然而如果我没有搞错的话，他的第一篇广义相对论的论文（即解释水星进动的论文）是有问题的，不过这都是后话了，扯远了哈，总之，本篇回答主要还是集中在经典领域下的。
行星的运行轨道可以是圆形，椭圆形，抛物线形，以及双曲线形，但是太阳系演化46亿年后的今天，太阳系内所有行星的轨道都是椭圆的。
任何在引力作用下运行的行星都会沿着圆锥截面运动，如下图，这些圆锥截面曲线包括圆、椭圆、抛物线和双曲线。同时，轨道的焦点/圆点是恒星。
圆锥截面曲线这四种轨道的不同之处在于它们的偏心率。圆轨道的离心率是0，一个完美的圆轨道；椭圆轨道离心率在0和1之间，也是绝大多数行星的轨道形状；抛物线轨道离心率为1；双曲线轨道离心率大于1，这两类轨道对应的天体从太阳系外经过不会返回。
四种轨道有了这四种行星轨道的选项，我们看看为什么行星的运行轨道绝大多数都是椭圆。
首先，太阳系有 46 亿年的历史，任何具有抛物线或双曲线轨道的行星都早已不复存在了，这些轨道上的行星或者彗星，由于轨道不闭合，对于太阳系来讲都是过客；其次， 圆形轨道需要实现正好为零的偏心率，这很难，但椭圆轨道的偏心率可以在 0 到 1 之间的任何位置，这就很简单。所以，我们看到的绝大多数的行星轨道都是椭圆。
但实际上，太阳系的很多行星的轨道离心率都很小，基本上也可以近似为圆形，但严格讲还是椭圆。比如：地球离心率是0.0167，近似为圆形；离心率最小的是金星，只有0.0067，最大的是水星，为0.2056，所以水星的轨道是太阳系中最椭圆的。
牛顿发现主导行星运动的力是一个只与距离平方反比相关的定律。当时牛顿的证明是很不容易的。即使到今天，用现代微积分的做法推导出来也是不容易的，但标准的理论力学书一般都有推导。
1.我们可以先来证明，如果行星的运动轨道是椭圆，那么力是平方反比的。这里还要用到开普勒第二定律或者角动量守恒，也即
r^2\dot{\phi}=Cr^2\dot{\phi}=C （1）
另外还要用到机械能守恒，也即
\frac{1}{2}m(\dot{r}^2+r^2\dot{\phi}^2)+V(r)=E\frac{1}{2}m(\dot{r}^2+r^2\dot{\phi}^2)+V(r)=E （2）
