Moonpos用书中的方法计算月球的位置
天文算法，简·米乌斯，第45章，第307 - 313页。的
所使用的算法在经度上精确到10弧秒，在经度上精确到4弧秒
在纬度上根据书，虽然时间范围不是
提及。此算法中使用的周期项是基于
Chapront ELP-2000/82月球理论，包括后来的改进
Chapront。这一理论的结果被称为平均春分
但月球将这些数据转换为明显的地心位置
月之分指的是真正的春分日。

语法

B = moonpos(月亮)
B = moonpos(A,S1)
B = moonpos(A,S1,S2，…)

描述

A =包含日期值的Nx1双向量。默认情况下，数据
的格式，但这可以是
改为儒略日，修改儒略日，简化儒略日，或
用可选字符串截断的Julian Dat(见下文)。

= Nx6双日期矩阵，其中每一行都是一个完整的日期向量
按此顺序指定年、月、日、小时、分、秒。
日期必须与MATLAB的数据枚举函数兼容。

= Nx3双日期矩阵，其中每一行都是部分日期向量
按此顺序指定年、月和日。日期必须是
兼容MATLAB的数据枚举函数。

的日期兼容的NxM字符数组
datenum函数。

对于最后三种情况，datenum函数被自动调用
将日期转换为SDN格式。

B =包含月球位置数据的NxM双向量。默认为
Nx3矢量返回三列为赤经(rad)的位置，
赤纬(rad)和距离(km)，但这可以改变
可选字符串(见下文)。输出是来自地球的真实日期
中心。

Sn =可选string参数，用于更改输入和输出设置。你
可以指定多个可选字符串。它们被左至处理
和不区分大小写。如果指定互斥
选项(如'jd'和'rjd')，以最右边的选项为准。

"朱利安" "朱利安日"或" jd "
=输入向量A为Julian Day格式。也就是时间的间隔
从公元前4713年1月1日格林尼治中午开始

'缩减朱利安日'或'缩减'或'rjd'
=输入向量A为简化朱利安日格式，RJD = JD - 2400000。

'modified Julian day'或者'modified'或者'mjd'
=输入向量A为Modified Julian Day格式，MJD = JD - 2400000.5

'截断的朱利安日'或'截断'或'tjd'
=输入向量A为截断朱利安日格式，TJD = JD - 2440000.5

“radec”
=输出Nx3矩阵的赤纬、赤纬、距离格式
单位为弧度、弧度、公里(这是默认值)。
但请注意，你可以改变角度或距离的单位
(见下文)，你可以加上x, y, z笛卡尔输出在
同时使用'xyz'选项，得到Nx6矩阵结果。

“xyz”
=输出是x, y, z笛卡尔坐标的Nx3矩阵。请注意
如果你愿意，这个可以在'radec'之外指定
包含两个数据集的Nx6矩阵输出。的顺序
哪个“radec”和“xyz”选项出现在参数列表中
确定输出在结果中出现的顺序。

'feet'或'ft'
=输出距离将以英尺代替公里。

'meter'或者'meters'或者'm'
=输出距离将以米为单位，而不是千米。

'degree'或'degrees'或'deg'
=输出角度将以角度代替弧度。

“号”
=输出赤经将在hh.mmsss…格式。

dms的
=产量下降将在dd.mmsss…格式。

“nutationlow”
=使用底部稍低精度的章动公式
而不是完整的表格21。一个术语。

例子

包含的文件moonpos_examples.m中给出了几个示例