当前位置:[北京同好会]>[编程]>[大行星历表计算]

 

大行星历表计算

MatLab源程序

以下球面天文学(岁差、章动)和天体力学(行星星历表)的原始计算公式、算法和程序主要选自:http://www.moshier.net/并改写为以.m为扩展名的Matlab的源程序,可以使用任何文本浏览器查看和编辑,并在Matlab 5.3.0下调试通过。

Matlab的编程语法与标准的数学表示方式非常相似,很容易阅读理解。从源程序中可以看出计算所依据的物理学、天文学的原理和定律。

 

通用程序

sign1.m 数字符号
d2r.m 度、分、秒转换为弧度
r2d.m 弧度转换为度、分、秒

r2h.m

弧度转换为小时、分、秒
header.m 通用常数
hel2geo.m 日心坐标转为地心坐标
getdate.m 输入日期

时间换算

cal2julian.m

公历/儒略历转换为儒略日
julian2cal.m 儒略日转换为公历/儒略历

deltat.m

地球力学时与世界时之差

球面天文学

epsiln.m

黄赤交角

ec2eq.m

黄道坐标与赤道坐标的转换
pol2rec.m 球坐标转换为直角坐标
rec2pol.m 直角坐标转换为球坐标
precess.m 黄道坐标岁差改正
nutate.m 章动系数
annuab.m 周年光行差

星历表计算

mer404,ven404,ear404 大行星摄动参数表
mar404,jup404,sat404
ura404
,nep404,plu404
大行星摄动参数表
moonlr.m, moonlat.m 月亮地心坐标摄动参数表
g3plan.m 地球日心黄道坐标三角函数级数解(J2000)
gplan.m 大行星日心黄道坐标三角函数级数解
g1plan,g2plan 月亮地心黄道坐标三角函数级数解
pearth.m 地球日心黄道坐标(J2000)
pmoon.m 月亮地心黄道坐标(J2000)
mean_elements.m 大行星平黄经

应用程序

appsun.m 太阳位置计算
appplanet.m 大行星位置计算
astro.zip 以上所有程序的压缩文件(98K)