- 中华农历论坛 (http://bbs.nongli.net/index.asp)
-- 历法知识 (http://bbs.nongli.net/list.asp?boardid=2)
---- 大行星星历计算程序源代码 (http://bbs.nongli.net/dispbbs.asp?boardid=2&id=18614)
-- 发布时间：2009/4/21 8:44:00
-- 大行星星历计算程序源代码
MatLab源程序
以下球面天文学（岁差、章动）和天体力学（行星星历表）的原始计算公式、算法和程序主要选自：http://www.moshier.net/并改写为以.m为扩展名的Matlab的源程序，可以使用任何文本浏览器查看和编辑，并在Matlab 5.3.0下调试通过。
Matlab的编程语法与标准的数学表示方式非常相似，很容易阅读理解。从源程序中可以看出计算所依据的物理学、天文学的原理和定律。
|
通用程序 |
||
| sign1.m | 数字符号 | |
| d2r.m | 度、分、秒转换为弧度 | |
| r2d.m | 弧度转换为度、分、秒 | |
| 弧度转换为小时、分、秒 | ||
| header.m | 通用常数 | |
| hel2geo.m | 日心坐标转为地心坐标 | |
| getdate.m | 输入日期 | |
|
时间换算 |
||
| 公历/儒略历转换为儒略日 | ||
| julian2cal.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 |
-- 发布时间：2009/4/21 8:45:00
--
上面是别人的 和许老师的相比，不知那个精度高
同时也希望许老师加上 紫气 月渤 计都 罗睺
这样才算是完整的农历
-- 发布时间：2017/11/5 21:41:00
--
请问楼主有用MATLAB计算行星速度的代码吗？
-- 发布时间：2018/3/29 12:18:00
--
谢谢楼主分享！！！