以文本方式查看主题

-  中华农历论坛  (http://bbs.nongli.net/index.asp)
--  历法知识  (http://bbs.nongli.net/list.asp?boardid=2)
----  大行星星历计算程序源代码  (http://bbs.nongli.net/dispbbs.asp?boardid=2&id=18614)

--  作者:ymy111
--  发布时间:2009/4/21 8:44:00
--  大行星星历计算程序源代码

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.rar

 

 


--  作者:ymy111
--  发布时间:2009/4/21 8:45:00
--  

上面是别人的 和许老师的相比,不知那个精度高

同时也希望许老师加上 紫气 月渤 计都 罗睺

这样才算是完整的农历


--  作者:wyf96
--  发布时间:2017/11/5 21:41:00
--  
请问楼主有用MATLAB计算行星速度的代码吗?
--  作者:CPU
--  发布时间:2018/3/29 12:18:00
--  
谢谢楼主分享!!!