《方程与宇宙》:抛物线与双曲线轨道(三)
By 苏剑林 | 2010-04-03 | 48562位读者 |经过上两回的讨论,我们已经基本摸清了二体问题的运动情况。我们已经找到了二体问题在轨道为椭圆的时候的所有积分,给出了“活力公式”等常用公式的证明,并且留下了一些没有解答的问题。那就是在轨道为抛物线和双曲线时的最后一个积分还没有找出来,现在我们解决这两个问题。其中的关键积分依旧是
$\dot{r}^2={2\mu}/r-{\mu a(1-e^2)}/r^2-\frac{\mu}{a}$——(12)
一、抛物线轨道
对于抛物线轨道,我们有$e=1,a->\infty$,但是近日点的关系式依然没有变,也就是说近日距依旧是$a(1-e)=p$,而由于a趋向无穷大,那么$\frac{\mu}{a}$可以忽略掉,于是(12)可以简化成:
$$\dot{r}^2=\frac{2\mu}{r}-\frac{2p\mu}{r^2}$$
不难推出
$$\sqrt{2\mu}dt=\frac{rdr}{\sqrt{r-p}}$$
并且令$R^2=r-p$,则上式变为
$$\sqrt{2\mu}dt=\frac{(R^2+p)d(R^2+p)}{R}=(2R^2+2p)dR$$
两端积分
$$\sqrt{2\mu}t=2/3 R^3+2pR+K_3$$
取t=0为近日点的时刻,那么可以得到$K_3=0$,于是抛物线的轨道问题就解决了:
$\sqrt{\frac{\mu}{2}}t=1/3 R^3+pR$,$$r=R^2+p\tag{16}$$
(附:这个不同于《天体力学引论》等教程的解法,这是直接求出向径r的长度,然后再求角度,过程刚好与原著相反。个人认为这种方法更好理解,而且更符合我们的思维。另外这是一道关于R的三次方程,我们也不必去寻找它的求根公式,只要使用方程的迭代解法针对某一数据就可以求出相应的R)
二、双曲线轨道
对于双曲线,有$e>1,a<0$,所以我们令$A=-a$,避免负数开平方的情况。(12)变为
$$\dot{r}^2={2\mu}/r-{\mu A(e^2-1)}/r^2+\frac{\mu}{A}\tag{17}$$从而推出
$$\sqrt{\frac{\mu}{A}}dt=\frac{rdr}{\sqrt{(A+r)^2-A^2 e^2}}$$
这个与椭圆轨道的方程类似,方法也是用换元法,但是不是用三角函数,而是用“双曲函数”。令$A+r=Ae cosh(E)$,则可以变换成
$$\sqrt{\frac{\mu}{A^3}}dt=(e \cos h E -1)dE$$
两端积分
$$\sqrt{\frac{\mu}{A^3}}t=e \sin h E -E+K_4$$
同样取t=0为近日点的时刻,那么可以得到$K_4=0$,于是可以得到双曲线的开普勒方程
$$\sqrt{\frac{\mu}{A^3}}t=e \sin h E -E\tag{18}$$这与椭圆轨道的开普勒方程有点类似,同样是超越方程。
附:简介一下双曲函数
根据欧拉的表示法,我们可以把余弦表示成
$$\cos \theta=\frac{e^{i\theta}+e^{-i\theta}}{2}$$
这里的$e=2.71828....,i^2=-1$,要是我们把里面的“i”去掉,可以定义一种新的函数,称为“双曲余弦函数”,记为\cos h,也就是说
$$\cos h \theta=\frac{e^{\theta}+e^{-\theta}}{2}$$
类似的:
$$\begin{aligned}\sin h \theta=\frac{e^{\theta}-e^{-\theta}}{2} \\ \tan h \theta=\frac{\sin h \theta}{\cos h \theta}\end{aligned}$$
对于求导数
$$\begin{aligned}(\sin h \theta)'=\cos h \theta \\ (\cos h \theta)'=\sin h \theta\end{aligned}$$
更多内容可以参考:
http://zh.wikipedia.org/wiki/%E5%8F%8C%E6%9B%B2%E5%87%BD%E6%95%B0
最后,二体问题部分就讲到这里了,还有其他一些诸如怎样求解开普勒方程的小问题,大家其实不用想的那么复杂。对于一道难以直接求解的方程的求解,数学上有通用的数值方法,可以用来求任意精确度的特解,对于天文应用来说来说,这已经足够了。另外《方程与宇宙》这个主题可能会有很长时间不更新了,因为要准备天文奥赛,以及想钻研一下三体问题的摄动解法,所以希望大家等待了。当然,blog还会坚持更新的。
转载到请包括本文地址:https://kexue.fm/archives/582
更详细的转载事宜请参考:《科学空间FAQ》
如果您还有什么疑惑或建议,欢迎在下方评论区继续讨论。
如果您觉得本文还不错,欢迎分享/打赏本文。打赏并非要从中获得收益,而是希望知道科学空间获得了多少读者的真心关注。当然,如果你无视它,也不会影响你的阅读。再次表示欢迎和感谢!
如果您需要引用本文,请参考:
苏剑林. (Apr. 03, 2010). 《《方程与宇宙》:抛物线与双曲线轨道(三) 》[Blog post]. Retrieved from https://kexue.fm/archives/582
@online{kexuefm-582,
title={《方程与宇宙》:抛物线与双曲线轨道(三)},
author={苏剑林},
year={2010},
month={Apr},
url={\url{https://kexue.fm/archives/582}},
}
May 23rd, 2010
望站主多多更新有关这方面的内容
这需要时间,还要我自己去研究。因为网站上的内容基本是我研究过后的心得与笔记。
当然我会努力去更新相关内容的。现在奥赛考完了,轻松了一点,可以抽多一点时间来建设
August 1st, 2010
我来挑骨头喽~~~e=2.718281828……上面打错了。
August 7th, 2010
呵呵,我继续纠错:
1、抛物线轨道的偏心率=1;
2、双曲线的开普勒方程中好像应该把A^3改为A;
3、(纯属挑骨头行为)最后一段话中应该把“人以”改为“任意”。
=================
向站长请教一下:那个抛物线轨道上含有R的一元三次方程尽管我采用“天文单位、平太阳日、太阳质量”单位制,取0为初值,但还是迭代得天昏地暗,有没有什么办法(莫非需要采用您的“终极算法”?)能加快收敛速度呢?
好吧,我会写一篇文章专门用来求解三次方程的。
我感觉你迭代出错了。你可以用年为时间单位,以AU为长度单位,然后$\mu=4\pi^2$,采用牛顿法,应该很快的。
A^3那里没有错误
不用了。我采用牛顿法,终于解决了这个关于迭代的问题。
关于A^3那里,是我自己积分代换时出现了问题,确实没有错误,抱歉……