matlab 转换 地平赤道坐标系,赤道坐标系变为地平坐标系

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

高度角(地平纬度)= arcsin[sin(纬度)sin(天体赤纬) + cos(纬度)cos(天体赤纬)cos(经度 - 天体光线直射点经度)]

如果是计算视高度角，还需要考虑大气折射的影响，实际看到天体的位置比这里计算出来的稍稍高一点。

其中天体光线直射点经度可以利用天体的赤经以及恒星时计算出来。

恒星时计算公式为

恒星时 = {0.671262+1.0027379094*TJD}*24

其中{}表示取小数部分，TJD为简化儒略日，它等于儒略日减去2440000.5，其计算方法可以自己百度。

天体光线直射点经度 = (赤经 - 恒星时)*15

如果算出来的天体光线直射点经度不在[-180，180]范围内，则对它加上或减去360°，使它落到这个区间内。

方位角(地平经度)= 180 ± arccos{[sin(高度角)sin(纬度) - sin(天体赤纬)] / [cos(高度角)cos(纬度)]}

判断天体是否已经过中天。过中天之前“±”号取-，中天之后取+。

这里计算的方位角是真方位角，即正北方向为0°，从正北开始顺时针旋转一周得到360°。例如西偏南70°用真方位角表示就是200°。

看以下例子：

计算2015年12月7日16:30(北京时间)时，北京(116°23'50" E，39°54'27" N)的太阳高度角和方位角。已知当时的太阳赤纬角为Dec = -22°34'52.6"，赤经为RA = 16h54m26.3s

北京时间2015年12月7日16:30时的TJD为1