- 1、本文档共12页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
卫星轨道与位置
水星的轨道和位置
摘要
本文主要在已知水星的远日点和绕日运行的线速度的条件下,通过建立微分方程模型,使用解析法和数值方法求解水星的轨道方程与位置。解析法的求解的过程中,结合了开普勒三大定律,准确的给出了微分方程的精确解,求得水星到太阳的最近距离,水星绕太阳运行的周期约为88天。数值计算求解水星自远日点运行50天后的位置时,本文分别采用了Simpson求积法,基于压缩映射的求根方法以及经典的四阶龙格—库塔法,使用matlab数学软件编程,得到了较为合理的行星运行模型的近似解,三种方法所得结果对应分,,,及,。
关键词 行星轨道 微分方程 Simpson法 四阶龙格—库塔法 matlab
问题重述
水星到太阳的最远距离为m,此时水星绕太阳运行的线速度为
m/s。试求
问题一 水星到太阳的最近距离
问题二 水星绕太阳运行的周期
问题三 从远日点开始的第50天(地球天)结束时水星的位置并画出轨道曲线
问题分析
求水星到太阳的最近距离以及水星绕太阳运行的周期等,需要先将水星轨道方程求出,因此可以根据Newton第二定律及万有引力定律,建立微分方程模型,将原问题转化为求解带有初值条件的微分方程问题,进而采用解析法或数值方法求解远日点和周期。
模型假设
1.水星运行的轨道是以太阳为一个焦点的椭圆
2.从太阳指向水星的线段在单位时间内扫过的面积相等
3.水星运行周期的平方与其运行轨道椭圆长轴的立方之比为常量
符号系统
1. 水星在远日点的线速度
2. 太阳的质量
3. 水星的质量
4. 水星在远日点的距离
5. 周期
建立模型与求解
模型一 水星的轨迹方程
设太阳中心所在的位置为复平面的原点O,在时刻t,水星位于
所表示的点P。这里均为t的函数,分别表示的模和辐角。于是水星的速度为,加速度为(1.1),而太阳对行星的引力依万有引力定律,大小为,方向由行星位置P指向太阳的中心O,故为,其中为太阳的质量,m为水星的质量,为万有引力常数。
依Newton定律,我们得到 (1.2),将(1.1)代入(1.2),然后比较实部与虚部,就有
这是两个未知函数的二阶微分方程组。在确定某一行星轨道时,需要加上定解条件。假设当t=0时,行星正处于远日点,而远日点位于正实轴上,距原点O为,行星的速度为。那么就有初值条件:
因此问题转化为求解带初值问题的微分方程组
又将两边同乘以r ,即得,从而(1.3),其中,这样有向线段在时间内扫过的面积等于,这个正是Kepler的第二定律,从太阳指向水星的线段在单位时间内扫过的面积相等。
将(1.3)代入得,于是我们可以得到水星运行的较为简单形式的数学模型:
为了求得行星的轨迹方程,要消去变量t,令,那么可以改写为从而将上式代入,化简后为 (1.4),其中,引进,立即可以求出,这里A和是待定的常数。记,上式可以写为
这个就是水星的轨道方程,是一条平面二次曲线。由于水星绕太阳运行,故必有。由于r在t=0时取道最大值(远日点),这个就意味着此时函数取道最大值1.于是就有 ,从而轨迹方程为 。对于水星而言,,又水星的近日点到太阳的距离。依据已知数据,可知,,,从而计算水星到太阳的最近距离为
模型二 水星的运行周期
设水星的周期为T,那么利用Kepler第二定律,我们有 (1.4)
上式左端为水星轨迹椭圆所围的面积,记为S,由于椭圆的半长轴,半短轴,从而有将上式代入式(1.4),解得 (1.5)将有关数据代入,易得
模型三 水星的位置
由于水星的运行满足Kepler第二定律,则该式可改写为,从而可得
如果我们要求时相应的和,则意味着首先要解方程 , ,其中
在求出了时的后,立即可以由得到相应的r。
下面用数值方法求解水星的位置
1. Simpson法
由被积函数的恒正性可知单调,从而方程的根必存在且唯一。取,记。若,那么位于与之间,在h适当小时,可取。
计算可采用不同的数值积分法,本文采用Simpson法,取步长h=0.001,具体求解过程见附录一,最后结果为,
2. 基于压缩映像的求根方法
我们引入水星轨道椭圆的参数方程,由于椭圆的半长轴,半短轴,从而中心到焦点的距离为。因左焦点为原点,故椭圆中心位于(ae,0),于是得到参数方程
它们与的关系为
此式可改写成
当时解方程
记,,那么上式即,就是说要去求函数的不动点,求解方程不动点可以采用简单迭代法,对于水星,我们已计算出,由于e很小,因此迭代收敛理论上可以很快,当时间从远日点开始的第50天结束
文档评论(0)