欢迎光临散文网 会员登陆 & 注册

用微分方程来推导开普勒三定律

2022-12-25 19:43 作者:求导宗师的线性空间  | 我要投稿

hello,大家好!

早在高中物理的必修课程中,我们就学过了开普勒三定律,它是行星运动的三条定律:

①椭圆定律:行星运动的轨迹是椭圆,且恒星处于椭圆的一个焦点上;

②面积定律:行星与恒星的连线在相同时间内扫过相同的面积;

③调和定律:行星轨迹的半长轴立方与行星公转周期平方的比为定值。

今天,我们就来用微分方程来推导这三条定律吧。

在此之前先明确几点:

①认为恒星静止不动,不受行星引力的影响

②只研究一个行星的运动,且该行星仅受恒星的引力影响

③不考虑相对论,只使用经典力学的观点

ok!开始吧!

首先,以恒星为原点建立平面直角坐标系。设 t 时刻行星的坐标为 (rcos%5Ctheta%20%2Crsin%5Ctheta%20) ,其中 r 与 %5Ctheta%20 都是关于时间 t 的函数;引力常数为 G,恒星质量为 M,行星质量为 m 。

将引力分解到 x 方向与 y 方向,根据牛顿第二定律有:

x 方向:-%5Cfrac%7BGMm%7D%7Br%5E2%7Dcos%5Ctheta%20%3Dma_%7Bx%7D

y 方向:-%5Cfrac%7BGMm%7D%7Br%5E2%7Dsin%5Ctheta%20%3Dma_%7By%7D

(注意引力的方向,要添个负号)

加速度是位移的二阶导数,因此有:a_%7Bx%7D%3D(rcos%5Ctheta%20)''%3D(r''-r(%5Ctheta%20')%5E2)cos%5Ctheta%20-(2r'%5Ctheta'%2Br%5Ctheta'')sin%5Ctheta%20

a_%7By%7D%3D(rsin%5Ctheta%20)''%3D(r''-r(%5Ctheta%20')%5E2)sin%5Ctheta%20%2B(2r'%5Ctheta'%2Br%5Ctheta'')cos%5Ctheta%20

(注意复合函数求导)

可将括号内的式子看作整体

记 A%3Dr''-r(%5Ctheta')%5E2B%3D2r'%5Ctheta'%2Br%5Ctheta''%20%20

可得:

Acos%5Ctheta-Bsin%5Ctheta%2B%5Cfrac%7BGM%7D%7Br%5E2%7Dcos%5Ctheta%3D0%20

Asin%5Ctheta%2BBcos%5Ctheta%2B%5Cfrac%7BGM%7D%7Br%5E2%7Dsin%5Ctheta%3D0

解得:

A%3D-%5Cfrac%7BGM%7D%7Br%5E2%7D ①,B%3D0 ②

先研究 ② 式,分离变量,得到:

%5Cfrac%7B2r'%7D%7Br%7D%3D-%5Cfrac%7B%5Ctheta''%7D%7B%5Ctheta'%7D 即 2%5Cfrac%7Bdr%7D%7Br%7D%3D-%5Cfrac%7Bd%5Ctheta'%7D%7B%5Ctheta'%7D

两边积分,得到:

r%5E2%5Ctheta'%3DC_%7B1%7D

这个式子很重要,我们记为(*)式。将(*)式代入 ① 式得到:

r''-%5Cfrac%7BC_%7B1%7D%5E2%7D%7Br%5E3%7D%2B%5Cfrac%7BGM%7D%7Br%5E2%7D%3D0

这是一个不含一阶导数的式子,我们有一种常用的降次方法:两边同时乘 2r'。得到:

2r'r''-2%5Cfrac%7BC_%7B1%7D%5E2r'%7D%7Br%5E3%7D%2B2%5Cfrac%7BGMr'%7D%7Br%5E2%7D%3D0

(r'%5E2%2B%5Cfrac%7BC_%7B1%7D%5E2%7D%7Br%5E2%7D-2%5Cfrac%7BGM%7D%7Br%7D)'%3D0

积分得:r'%5E2%2B%5Cfrac%7BC_%7B1%7D%5E2%7D%7Br%5E2%7D-2%5Cfrac%7BGM%7D%7Br%7D%3DC_%7B2%7D

将 r' 分离出来,再结合(*)式替换掉 t 得到:

r'%3D%5Cpm%20%5Csqrt%7BC_%7B2%7D%2B2%5Cfrac%7BGM%7D%7Br%7D-%5Cfrac%7BC_%7B1%7D%5E2%7D%7Br%5E2%7D%7D

r'%3D%5Cfrac%7Bdr%7D%7Bdt%7D%3D%5Cfrac%7Bdr%7D%7Bd%5Ctheta%7D%5Cfrac%7Bd%5Ctheta%7D%7Bdt%7D%3D%5Cfrac%7Bdr%7D%7Bd%5Ctheta%7D%5Cfrac%7BC_%7B1%7D%7D%7Br%5E2%7D

我们终于得到了关于 r 和 %5Ctheta%20 的微分方程:

d%5Ctheta%3D%5Cpm%5Cfrac%7Bdr%7D%7Br%5E2%5Csqrt%7B%5Cfrac%7BC_%7B2%7D%7D%7BC_%7B1%7D%5E2%7D%2B2%5Cfrac%7BGM%7D%7BC_%7B1%7D%5E2r%7D-%5Cfrac%7B1%7D%7Br%5E2%7D%7D%7D

解之得:

%5Cpm%20%5Ctheta%2BC_%7B3%7D%3Darccos(%5Cfrac%7B%5Cfrac%7B1%7D%7Br%7D-%5Cfrac%7BGM%7D%7BC_%7B1%7D%5E2%7D%7D%7B%5Csqrt%7B%5Cfrac%7BC_%7B2%7D%7D%7BC_%7B1%7D%5E2%7D%2B%5Cfrac%7BG%5E2M%5E2%7D%7BC_%7B1%7D%5E4%7D%7D%7D)

不妨取 C_%7B3%7D%3D0,整理得:

r%3D%5Cfrac%7B%5Cfrac%7BC_%7B1%7D%5E2%7D%7BGM%7D%7D%7B1%2Bcos%5Ctheta%20%5Csqrt%7B%5Cfrac%7BC_%7B2%7DC_%7B1%7D%5E2%7D%7BG%5E2M%5E2%7D%2B1%7D%7D

%5Cfrac%7BC_%7B1%7D%5E2%7D%7BGM%7D%3Dpe%5Csqrt%7B%5Cfrac%7BC_%7B2%7DC_%7B1%7D%5E2%7D%7BG%5E2M%5E2%7D%2B1%7D%3De

得到:r%3D%5Cfrac%7Bpe%7D%7B1%2Becos%5Ctheta%20%7D

这是圆锥曲线的极坐标方程,另外行星运动的轨迹应是封闭的,于是行星的运动轨迹为椭圆,且恒星位于坐标原点,也就是椭圆的一个焦点,这便是开普勒第一定律。

回到(*)式,由极坐标面积公式,任意 %5CDelta%20t 时间内,行星与恒星的连线扫过的面积

S%3D%5Cint_%7Bt_%7B0%7D%7D%5E%7Bt_%7B0%7D%2B%5CDelta%20t%7D%20%5Cfrac%7B1%7D%7B2%7Dr%5E2d%5Ctheta%20%3D%5Cint_%7Bt_%7B0%7D%7D%5E%7Bt_%7B0%7D%2B%5CDelta%20t%7D%20%5Cfrac%7B1%7D%7B2%7Dr%5E2%5Cfrac%7Bd%5Ctheta%7D%7Bdt%7Ddt

    %3D%5Cint_%7Bt_%7B0%7D%7D%5E%7Bt_%7B0%7D%2B%5CDelta%20t%7D%20%5Cfrac%7B1%7D%7B2%7DC_%7B1%7Ddt%3D%5Cfrac%7B1%7D%7B2%7DC_%7B1%7D%5CDelta%20t

这表明任意 %5CDelta%20t 时间内行星与恒星连线扫过的面积为定值,这便是开普勒第二定律。

而一个周期内,连线扫过的面积是整个椭圆的面积 %5Cpi%20ab,再结合轨迹方程有:

%5Cpi%20ab%3D%5Cfrac%7BC_%7B1%7DT%7D%7B2%7Dpe%3D%5Cfrac%7Bb%5E2%7D%7Ba%7D%3D%5Cfrac%7BC_%7B1%7D%5E2%7D%7BGM%7D

可得 %5Cfrac%7Ba%5E3%7D%7BT%5E2%7D%3D%5Cfrac%7BGM%7D%7B4%5Cpi%20%5E2%7D

即半长轴立方与周期平方之比为定值,这便是开普勒第三定律。

至此我们就证明了开普勒三定律!

进一步观察行星运动轨迹的方程,它包含了所有二次曲线。也就是说,只要某个天体的速度合适,它的离心率 e 就足够大,只要大于等于 1,它就能脱离恒星引力的束缚,逃逸出去,这也为航空航天奠定了理论基础。


好了,就玩这么多,感谢大家收看!

拜拜~~

用微分方程来推导开普勒三定律的评论 (共 条)

分享到微博请遵守国家法律