地球轨道计算方法,开普勒三大定律的证明,补充-解析几何之椭圆

地球轨道计算方法
http://math.sjtu.edu.cn/jidi/sxsy/

开普勒三大定律的证明
http://zh.wikipedia.org/wiki/%E5%BC%80%E6%99%AE%E5%8B%92%E5%AE%9A%E5%BE%8B#.E5.BC.80.E6.99.AE.E5.8B.92.E7.AC.AC.E4.B8.89.E5.AE.9A.E5.BE.8B

开普勒定律

维基百科,自由的百科全书

跳转到: 导航, 搜索
图示遵守开普勒行星运动定律的两个行星轨道与。 (1) 第一个行星的椭圆轨道焦点是 f1\,\!f2\,\! ,第二个行星的椭圆轨道焦点是 f1\,\!f3\,\! 。太阳的位置是在点 f1\,\! 。 (2) A1 与 A2 是两个面积相等的阴影区域。太陽与第一个行星的連線,扫过这两个阴影区域,所需的時間相等。(3) 各个行星绕太阳公转周期的比率为a1^{3/2}:a2^{3/2}\,\! ;这里,a1\,\!a2\,\! 分别为第一个行星与第二个行星的半长轴长度。

开普勒定律开普勒发现的关于行星运动的定律。他於1609年在他出版的《新天文学》上发表了关于行星运动的两条定律,又於1618年,发现了第三条定律。

开普勒很幸运地能够得到,著名的丹麦天文学家第谷·布拉赫所观察与收集的,非常精确的天文資料。大约于 1605 年,根据布拉赫的行星位置資料,开普勒发现行星的移动遵守三条相当简单的定律。

开普勒的定律给予亚里士多德派与托勒密派在天文学与物理学上极大的挑战。他主张地球是不斷地移动的;行星轨道不是周转圆 (epicycle 的,而是椭圆形的;行星公转的速度不等恒。这些论点,大大地动摇了当时的天文学与物理学。经过了几乎一世纪披星戴月,废寝忘食的研究,物理学家终于能够用物理理论解释其中的道理。牛顿利用他的第二定律万有引力定律,在数学上严格地証明开普勒定律,也让人们了解其中的物理意义。

目录

[隐藏]

[编辑] 开普勒定律

开普勒的三条行星运动定律改变了整个天文学,彻底摧毁了托勒密复杂的宇宙体系,完善并简化了哥白尼日心说

[编辑] 开普勒第一定律

开普勒第一定律

开普勒第一定律,也称椭圆定律:每一个行星都沿各自的椭圆轨道环绕太阳,而太阳则处在椭圆的一个焦点中。

[编辑] 开普勒第二定律

开普勒第二定律

开普勒第二定律,也称面积定律:在相等时间内,太阳运动着的行星的连线所扫过的面积都是相等的。

这一定律实际揭示了行星绕太阳公转的角动量守恒。用公式表示为

S_{AB}=S_{CD}=S_{EK}\,\!

[编辑] 开普勒第三定律

开普勒第三定律

开普勒第三定律,也称调和定律:各个行星太阳公转周期平方和它们的椭圆轨道半长轴立方正比

由这一定律不难导出:行星与太阳之间的引力与半径的平方成反比。这是牛顿万有引力定律的一个重要基础。

用公式表示为

\frac{a^3}{\tau^2}=K\,\!

这里,a\,\! 是行星公转轨道半长轴,\tau\,\! 是行星公转周期,K\,\! 是常数。

[编辑] 数学导引

开普勒定律是关于行星环绕太阳的运动,而牛顿定律更广义的是关于几个粒子因万有引力相互吸引而产生的运动。在只有两个粒子,其中一个粒子超轻于另外 一个粒子,这些特别状况下,轻的粒子会环绕重的粒子移动,就好似行星根据开普勒定律环绕太阳的移动。然而牛顿定律还容许其它解答,行星轨道可以呈抛物线运动或双曲线运动。这是开普勒定律无法预测到的。在一个粒子并不超轻于另外一个粒子的状况下,依照广义二体问题的解答,每一个粒子环绕它们的共同质心移动。这也是开普勒定律无法预测到的。

开普勒定律,或者是用几何语言,或者是用方程式,将行星的坐标及时间跟轨道参数相连结。牛顿第二定律是一个微分方程式。开普勒定律的导引涉及解微分方程式的艺术。我们会先导引开普勒第二定律,因为开普勒第一定律的导引必须建立于开普勒第二定律。

[编辑] 开普勒第二定律导引

牛顿万有引力定律阐明:任意两个粒子由通过连线方向的力相互吸引。该引力的的大小与它们的质量乘积成正比,与它们距离的平方成反比。由于太阳超重于行星,我们可以假设太阳是固定的。用方程式表示,

 \mathbf{F}= - G \frac{mM}{r^2}\ \hat{\mathbf{r}}\,\!

这里,\mathbf{F}\,\! 是太阳作用於行星的万有引力、m\,\! 是行星的质量、M\,\! 是太阳的质量、\mathbf{r}\,\! 是行星相对于太阳的位移向量、\hat{\mathbf{r}}\,\!\mathbf{r}\,\!单位向量

牛顿第二定律声明:物體受力後所产生的加速度 \ddot{\mathbf{r}}\,\! ,和其所受的淨力 \mathbf{F}\,\! 成正比,和其質量 m\,\! 成反比。用方程式表示,

\mathbf{F}=m\ddot{\mathbf{r}}\,\!

合并这两个方程式,

\ddot{\mathbf{r}}= - G \frac{M}{r^2}\ \hat{\mathbf{r}}\,\!(1)

思考位置向量 \mathbf{r}= r \hat{\mathbf{r}}\,\! ,随时间 t\,\! 微分一次可得到速度向量,再微分一次则可得到加速度向量:

\dot\mathbf{r} =\dot{r}\hat\mathbf{r} + r \dot{\hat\mathbf{r}}
=\dot{r}\hat\mathbf{r} +r\dot{\theta}\hat{\boldsymbol{\theta}}\,\!
\ddot\mathbf{r} 
= (\ddot{r}\hat{\mathbf{r}} +\dot{r}\dot\hat{\mathbf{r}} )
+ (\dot{r}\dot{\theta}\hat{\boldsymbol{\theta}} + r\ddot{\theta}\hat{\boldsymbol{\theta}}
+ r\dot{\theta} \dot\hat{\boldsymbol{\theta}})
= (\ddot{r} - r\dot{\theta}^2) \hat{\mathbf{r}} + (r\ddot{\theta} + 2\dot{r}\dot{\theta}) \hat{\boldsymbol{\theta}}\,\!(2)

在这里,我们用到了单位向量微分方程式:

 \dot\hat{\mathbf{r}} = \dot{\theta}  \hat{\boldsymbol{\theta}}\,\!
 \dot\hat{\boldsymbol{\theta}} = - \dot{\theta} \hat{\mathbf{r}}\,\!

合并方程式 (1) 与 (2) ,可以得到向量运动方程式:

 (\ddot{r} - r\dot{\theta}^2) \hat{\mathbf{r}} + (r\ddot{\theta} + 2\dot{r} \dot{\theta}) \hat{\boldsymbol{\theta}}=  - \frac{GM}{r^2}\hat{\mathbf{r}}\,\!

取各个分量,我们得到两个常微分方程式,一个是关于径向加速度,另一个是关于切向加速度:

\ddot{r} - r\dot{\theta}^2 = - \frac{GM}{r^2}\,\!(3)
r\ddot{\theta} + 2\dot{r}\dot{\theta} = 0\,\!(4)

导引开普勒第二定律只需切向加速度方程式。试想行星的角动量 \ell=m r^2\dot{\theta}\,\! 。由于行星的质量是常数,角动量随时间的导数为

\dot{\ell}=mr (2\dot{r}\dot{\theta}+r\ddot{\theta})=0\,\!

角动量 \ell\,\! 也是一个运动常数,即使距离 r\,\! 与角速度 \dot{\theta}\,\! 都可能会随时间变化。

从时间 t_1\,\! 到时间 t_2\,\! 扫过的区域 \Delta A\,\!

\Delta A=\int_{t_1}^{t_2}\frac{1}{2}\cdot r\cdot r\dot \theta\cdot dt=\int_{t_1}^{t_2}\frac{\ell}{2m}dt=\frac{\ell}{2m}\cdot(t_2-t_1) \,\!

行星太阳连线扫过的区域面积相依于间隔时间 t_2 - t_1\,\! 。所以,开普勒第二定律是正确的。

[编辑] 开普勒第一定律导引

设定 u =\frac{1}{r}\,\! 。这样,角速度是

\dot{\theta} =\frac{\ell}{mr^2}=\frac{\ell u^2}{m}\,\!

随时间微分与随角度微分的关系为

\frac{d}{dt}=\dot{\theta}\frac{d}{d\theta}=\frac{\ell u^2}{m} \frac{d}{d\theta}\,\!

随时间微分徑向距離 r=\frac{1}{u}\,\!

\dot{r}=\frac{\ell u^2}{m} \frac{d}{d\theta}\frac{1}{u} = - \frac{\ell u^2}{m}\frac{1}{u^2}\frac{du}{d\theta}= - \frac{\ell }{m}\frac{du}{d\theta} \,\!

再微分一次:

\ddot{r} =\frac{\ell u^2}{m} \frac{d\dot{r}}{d\theta}=\frac{\ell u^2}{m}  \frac{d}{d\theta} ( - \frac{\ell}{m}\frac{du}{d\theta})=  - \frac{\ell^2 u^2}{m^2}\frac{d^2 u}{d\theta^2}\,\!

代入径向运动方程式 (3) ,\ddot{r} - r\dot{\theta}^2 = - \frac{GM}{r^2}\,\!

 - \frac{\ell^2 u^2}{m^2}\frac{d^2 u}{d\theta^2} - \frac{\ell^2 u^3}{m^2} = - GMu^2\,\!

将此方程式除以  - \frac{\ell^2 u^2}{m^2}\,\! ,则可得到一个简单的常係数非齐次线性全微分方程式来描述行星轨道:

\frac{d^2u}{d\theta^2} + u = \frac{GMm^2}{\ell^2}\,\!

特征方程式为

u=\frac{GMm^2}{\ell^2}\,\!

求解剩馀的常係数齐次线性全微分方程式

\frac{d^2u}{d\theta^2} + u = 0\,\!

其特解方程式为

u=\epsilon\ \cos(\theta - \theta_0) \,\!

这里,\epsilon\,\!\theta_0\,\! 都是任意积分常数。综合特征方程式与特解方程式,

u=\frac{GMm^2}{\ell^2}+\epsilon\ \cos(\theta-\theta_0)\,\!

选择坐标轴,让 \theta_0=0\,\! 。代回 u=\frac{1}{r}\,\!

\frac{1}{r}=\frac{GMm^2}{\ell^2}+ \epsilon\ \cos{\theta}\,\!

假若 \epsilon\ <\ 1\,\! ,则 r\,\! 所描述的是椭圆轨道。所以,开普勒第一定律是正确的。

[编辑] 开普勒第三定律导引

在建立牛顿万有引力定律的概念与数学架构上,开普勒第三定律是牛顿依据的重要线索之一。假若我们接受牛顿运动定律。试想一个虚拟行星环绕着太阳公转,行星的移动轨道恰巧呈圆形,轨道半径为 r\,\! 。那末,太阳作用于行星的万有引力为 F=\frac{mv^2}{r}\,\! 。行星移动速度为 v=\frac{2\pi r}{\tau}\,\! 。依照开普勒第三定律,这速度 v\,\! 与半径的平方根 \sqrt{r}\,\! 成反比。所以,万有引力 F\propto \frac{1}{r^2}\,\! 。猜想这大概是牛顿发现万有引力定律的思路,虽然我们并不能完全确定,因为我们无法在他的计算本裡,找到任何关于这方面的证据。

行星环绕太阳(焦点 F1 )的椭圆轨道。

开普勒第一定律阐明,行星环绕太阳的轨道是椭圆形的。椭圆的面积是 \pi ab\,\! ;这里,a\,\!b\,\! 分别为椭圆的半長軸半短軸。在开普勒第二定律导引里,行星-太阳连线扫过区域速度 \frac{dA}{dt}\,\!

\frac{dA}{dt}=\frac{\ell}{2m}\,\!

所以,行星公转周期 \tau\,\!

\tau=\frac{2m\pi ab}{\ell}\,\!(5)

关于此行星环绕太阳,椭圆的半長軸 a\,\! ,半短軸 b\,\!近拱距 r_A\,\!近拱点 A 与引力中心之间的距离),远拱距 r_B\,\!远拱点 B 与引力中心之间的距离)的关系分别为

a=(r_A+r_B) / 2\,\!(6)
b=\sqrt{r_A\ r_B}\,\!(7)

如果想要知道半長軸与半短軸,必须先求得近拱距与远拱距。依据能量守恒定律

E=\frac{1}{2}m\dot{r}^2+\frac{1}{2}mr^2\dot{\theta}^2 - G\frac{mM}{r}\,\!

在近拱点 A 与远拱点 B,径向速度都等于零:

\dot{r}=0\,\!

所以,

E=\frac{1}{2}mr^2\dot{\theta}^2 - G\frac{mM}{r}= \frac{\ell^2}{2mr^2} - G\frac{mM}{r}\,\!

稍为加以编排,可以得到 r\,\! 的一元二次方程式:

r^2+\frac{GmM}{E}r -  \frac{\ell^2}{2mE}=0\,\!

其兩個根分别为椭圆轨道的近拱距 r_A\,\! 与远拱距 r_B\,\!

r_A=\left( - \frac{GmM}{E} - \sqrt{\left(\frac{GmM}{E}\right)^2+\frac{2\ell^2}{mE}}\right) / 2\,\!
r_B=\left( - \frac{GmM}{E} + \sqrt{\left(\frac{GmM}{E}\right)^2+\frac{2\ell^2}{mE}}\right) / 2\,\!

代入方程式 (6) 与 (7) ,

a= - \frac{GmM}{2E}\,\!
b=\frac{\ell}{\sqrt{ - 2mE}}=\frac{\ell}{m}\frac{\sqrt{a}}{\sqrt{GM}}\,\!

代入方程式 (5) ,周期的方程式为

\tau=\frac{2\pi a^{3/2}}{\sqrt{GM}}\,\!

Q.E.D.

补充-解析几何之椭圆

椭圆

维基百科,自由的百科全书

跳转到: 导航, 搜索
椭圆和它的某些数学性质。

数学中,椭圆是平面上到两个固定点的距离之和是常数的轨迹。这两个固定点叫做焦点

经由这个定义,这样画出一个椭圆:先准备一条线,将这条线的两端各绑在一点上(这两个点就当作是椭圆的两个焦点);取一支笔,将线绷紧,这时候两个点和笔就形成了一个三角形;然后拉着线开始作图,持续的使线绷紧,最后就可以完成一个椭圆的图形了。

目录

[隐藏]

[编辑] 概述

椭圆是一种圆锥曲线: 如果一个平面切截一个圆锥面,且不与它的底面相交,也不与它的底面平行,则圆锥和平面交截线是个椭圆。

在代数上说,椭圆是在笛卡尔平面上如下形式的方程所定义的曲线

A x^2 + 2B xy + C y^2 + D x + E y + F = 0 \,

使得 B^2 < AC \,,这里的系数都是实数,并存在定义在椭圆上的点对 (x, y) 的多于一个的解。

穿过两焦点并终止于椭圆上的线段 AB 叫做长轴。长轴是通过连接椭圆上的两个点所能获得的最长线段。穿过中心(两焦点的中点)垂直于长轴并且终止于椭圆的线段 CD 叫做短轴半长轴(图中指示为 a)是长轴的一半: 从中心通过一个焦点到椭圆的边缘的线段。类似的,半短轴(图中指示为 b)是短轴的一半。

如果两个焦点重合,则这个椭圆是;换句话说,圆是离心率为零的椭圆的特殊情况。

中心位于原点的椭圆  A x^2 + 2B xy + C y^2 = 1 \, 可以被看作单位圆在关联于对称矩阵 A^\prime =\begin{bmatrix}A & B\\B & C\end{bmatrix} = PDP^T \,线性映射下的图像,这里的 D 是带有 A^\prime特征值对角矩阵,二者沿着主对角线都是正实数的,而 P 是拥有 A^\prime特征向量作为纵列的实数的酉矩阵。椭圆的长短轴分别沿着 A^\prime 的两个特征向量的方向,而两个与之对应的特征值分别是半长轴半短轴的长度的平方的倒数。

椭圆可以通过对一个圆的所有点的 x 坐标乘以一个常数而不改变 y 坐标来生成。

[编辑] 离心率

椭圆的形状可以用叫做椭圆的离心率的一个数来表达,习惯上指示为 \varepsilon \,。离心率是小于 1 大于等于 0 的正数。离心率 0 表示着两个焦点重合而这个椭圆是

对于有半长轴 a 和半短轴 b 的椭圆,离心率是

\varepsilon = \sqrt{1 - \frac{b^2}{a^2}}.

离心率越大,ab比率就越大,因此椭圆被更加拉长。

如果 c 等于从中心到任一焦点的距离,则

\varepsilon = \frac{c}{a}.

距离 c 叫做椭圆的线性离心率。在两个焦点间的距离是 2aε。

[编辑] 方程

中心位于点 (h,k) 的主轴平行于 x 轴的椭圆由如下方程指定

\frac{(x-h)^{2}}{a^{2}} + \frac{(y-k)^{2}}{b^{2}} = 1 .

这个椭圆可以参数化表达为

x = h+a\,\cos t,\,\!
y = k+b\,\sin t\,\!

这里的 t 可以限制于区间 -\pi\leq t \leq \pi\,\!

如果 h = 0k = 0 (就是说,如果中心是原点(0,0)),则

椭圆方程 \frac{x^2}{a^2}+\frac{y^2}{b^2}=1(a>b>0) \frac{y^2}{a^2}+\frac{x^2}{b^2}=1(a>b>0)
图像
范围 -a\le x\le a,  -b\le y\le b -a\le y\le a,  -b\le x\le b
椭圆(用红色绘制)可以表达为内旋轮线在 R=2r 时的特殊情况。

用极坐标可表达为

r = \frac{ab}{\sqrt{a^2 \sin^2 \theta + b^2 \cos^2 \theta}}=\frac{b}{\sqrt{1-\varepsilon^2 \cos^2 \theta}}

这里的 \varepsilon 是椭圆的离心率。

有一个焦点在原点的椭圆的极坐标方程是

r = \frac{ a\cdot(1-\varepsilon^{2})}{1 + \varepsilon\cdot\cos\theta} .

[编辑] 半正焦弦和极坐标

椭圆的半正焦弦,通常指示为 l\,\!),是从椭圆的一个焦点到椭圆自身,沿着垂直主轴的直线测量的距离。它有关于 a\,\!b\,\! (椭圆的半轴),通过公式 al=b^2\,\! 或者如果使用离心率的话 l=a\cdot(1-\varepsilon^2)\,\!

椭圆,使用半正焦弦展示

极坐标中,一个焦点在原点而另一个焦点在负 x 轴上的椭圆给出自方程

r\cdot(1 + \varepsilon\cdot \cos \theta) = l \,\!

椭圆可以被看作是圆的投影: 在与水平面有角度 φ 的平面上的圆垂直投影到水平面上给出离心率 sin φ 的椭圆,假定 φ 不是 90°。

[编辑] 面积和周长

椭圆所包围的面积是 \pi ab \,,这里的 'a' 和 'b' 是半长轴和半短轴。在圆的情况下 a = b,表达式简化为 πa2

椭圆的周长是 4 a E(\varepsilon),这里的函数 E 是第二类完全椭圆积分

精确的无穷级数为:

C = 2\pi a \left[{1 - \left({1\over 2}\right)^2\varepsilon^2 - \left({1\cdot 3\over 2\cdot 4}\right)^2{\varepsilon^4\over 3} - \left({1\cdot 3\cdot 5\over 2\cdot 4\cdot 6}\right)^2{\varepsilon^6\over5} - \dots}\right]\!\,

或:

C = 2\pi a \sum_{n=0}^\infty {\left\lbrace - \left[\prod_{m=1}^n \left({ 2m-1 \over 2m}\right)\right]^2 {\varepsilon^{2n}\over 2n - 1}\right\rbrace}

一个好的近似是拉马努金的:

C \approx \pi \left[3(a+b) - \sqrt{(3a+b)(a+3b)}\right]\!\,

它还可以写为:

C \approx \pi a \left[ 3 (1+\sqrt{1-\varepsilon^2}) - \sqrt{(3+ \sqrt{1-\varepsilon^2})(1+3 \sqrt{1-\varepsilon^2})} \right] \!\,

[编辑] 标准方程的推导

  • 如果在一个平面内一个动点到两个定点距离等于定长,那么这个动点的轨迹叫做椭圆。

假设(注意所有假设只是为了导出椭圆方程时比较简便)动点为P(x,y) \,,两个定点为F_1(-c,0) \,F_2(c,0) \,,则根据定义,动点P的轨迹方程满足(定义式):

|PF_1|+|PF_2|=2a (a>0) \,,其中2a \,为定长。

用两点的距离公式可得:|PF_1|=\sqrt{(x+c)^2+y^2} \,|PF_2|=\sqrt{(x-c)^2+y^2} \,,代入定义式中,得:

\sqrt{(x+c)^2+y^2}=2a-\sqrt{(x-c)^2+y^2} \,

整理上式,并化简,得:

(a^2-c^2)x^2+a^2y^2=a^2(a^2-c^2) \,

a>c \,时,并设a^2-c^2=b^2 \,,则①式可以进一步化简:

b^2x^2+a^2y^2=a^2b^2 \,

因为a^2b^2>0 \,,将②式两边同除以a^2b^2 \,,可得:

x^2/a^2+y^2/b^2=1 \,

则该方程即动点P的轨迹方程,即椭圆的方程。这个形式也是椭圆的标准方程

  • 椭圆的图像如果在直角坐标系中表示,那么上述定义中两个定点被定义在了x轴。若将两个定点改在y轴,可以用相同方法求出另一个椭圆的标准方程
\frac{y^2}{a^2}+\frac{x^2}{b^2}=1(a>b>0) \,
  • 在方程中,所设的2a \,称为长轴长,2b \,称为短轴长,而所设的定点称为焦点,那么2c \,称为焦距。在假设的过程中,假设了a>c \,,如果不这样假设,会发现得不到椭圆。当a=c \,时,这个动点的轨迹是一个线段;当a<c \,时,根本得不到实际存在的轨迹,而这时,其轨迹称为虚椭圆。另外还要注意,在假设中,还有一处:a^2-c^2=b^2 \,
  • 通常认为是椭圆的一种特殊情况。

[编辑] 椭圆的旋转和平移

对于平面上任意椭圆  A x^2 + 2B xy + C y^2 + D x + E y + F = 0 \,,我们总可以将之转化为

A(x-u)^2 + 2B(x-u)(y-v) + C(y-v)^2 + f = 0 \,

的形式。具体步骤为,将后式的各乘积乘方项展开,根据与前式对应项系数相等的法则便可求得u,v,f的值。其中,(u,v) \,便是该椭圆的中心。

若将

x=x^\prime - u
y=y^\prime - v

带入式中便可得到平移前的椭圆。

B\ne 0,则表示椭圆的长短轴与坐标系的坐标轴并不平行或垂直,即发生了旋转。设旋转的角度为\displaystyle \phi, 则有

\displaystyle  tan(2 \phi)=2B/(A-C) 当A-C=0,则说明\phi=\pm \frac{\pi}{4}

若将

x=x^\prime cos \phi + y^\prime sin \phi
y=y^\prime cos \phi - x^\prime sin \phi

带入式中便可得到旋转前的椭圆。

[编辑] 参见

[编辑] 外部链接


posted on 2009-02-15 00:35  cutepig  阅读(31149)  评论(1编辑  收藏  举报

导航