技术特征：
1.一种静止轨道卫星在轨自主地影计算方法，其特征在于，包括：步骤s1：根据太阳圆盘张角和地球圆盘张角计算地影区圆盘的半径；步骤s2：根据卫星定点位置星下点经度计算卫星午夜时刻相对于格林尼治时间的累计秒数；步骤s3：根据卫星午夜时刻相对于格林尼治时间的累计秒数，计算午夜时刻太阳高度角；步骤s4：根据地影区圆盘半径和午夜时刻的太阳高度角，计算地影期间太阳划过的角度；步骤s5：根据地影期间太阳划过的角度，计算卫星地影持续时间、卫星进出影时刻并进行地影区判断。2.根据权利要求1所述的静止轨道卫星在轨自主地影计算方法，其特征在于，在所述步骤s1中：太阳圆盘张角为由卫星定点位置看向太阳的圆盘张角，由于地心-日心平均距离远大于卫星-地心平均距离，在计算过程中将该张角简化为地心位置处的太阳圆盘张角，依据太阳平均半径r
sun
及地-日平均距离r
se
计算当前卫星位置处的太阳半张角度；地球圆盘张角为由卫星定点位置看向地球的圆盘张角，由于地球经圈半径与赤道半径相差极小忽略不计，在计算过程中将地球简化为半径为6371.23km的球体；依据地球平均半径r
e
及星-地平均距离r
ce
计算当前卫星位置处的地球半张角度；地影区圆盘的半径a
se
＝a
sun
+a
eth
＝8.8591度。3.根据权利要求1所述的静止轨道卫星在轨自主地影计算方法，其特征在于，在所述步骤s2中：依据卫星定点位置星下点经度λ
sat
，单位：度，计算卫星相对于格林尼治时间的时间差δt＝240λ
sat
秒；在卫星定点位置午夜时刻，地心-太阳矢量在地球赤道面内的投影、地心-卫星矢量在地球赤道面内的投影之间的夹角为180
°
，即卫星星下点当地地方时为00:00:00，此时的格林尼治时间为：t
g
＝-δt；计算卫星午夜时刻相对于格林尼治时间2000年1月1日12:00:00时刻的累计秒数t0＝t
g0
+t
g
＝t
g0-δt；式中，t
g0
为格林尼治时间当日零点相对于格林尼治时间2000年1月1日12:00:00的累积秒计数。4.根据权利要求1所述的静止轨道卫星在轨自主地影计算方法，其特征在于，在所述步骤s3中：太阳高度角描述为地心-太阳矢量与地球赤道面内之间的夹角，当太阳位于北半球上方时该角度为正，太阳位于南半球上方时该角度为负；春/秋分时刻太阳高度角变化最快，整个地影持续时间内太阳高度角变化值<0.03
°
，若以午夜时刻太阳高度角作为整个地影期间太阳高度角，最大角度偏差<0.015
°
，对地影持续时间计算结果的影响忽略；
午夜时刻t0的太阳高度角：β＝arcsin[sin(23.4366)cos(1.14079552584275
×
10-5
t
0-169.699038446726)]。5.根据权利要求1所述的静止轨道卫星在轨自主地影计算方法，其特征在于，在所述步骤s4中：以卫星作为中心参考点，地影过程近似描述为太阳圆盘水平划过地影区圆盘，计算地影期间太阳划过的角度6.根据权利要求1所述的静止轨道卫星在轨自主地影计算方法，其特征在于，在所述步骤s5中：计算太阳圆盘划过地影区圆盘的角速度ω0＝360/86400＝0.004167
°
/s；地影持续时间秒。7.根据权利要求1所述的静止轨道卫星在轨自主地影计算方法，其特征在于，在所述步骤s5中：卫星进地影时刻卫星出地影时刻8.根据权利要求1所述的静止轨道卫星在轨自主地影计算方法，其特征在于，在所述步骤s5中：当t
in
≤t＜t
out
时，卫星处于地影区，否则卫星处于光照区；其中地影区包含了半影区和全影区；其中，t表示当前时刻。9.一种静止轨道卫星在轨自主地影计算系统，其特征在于，包括：模块m1：根据太阳圆盘张角和地球圆盘张角计算地影区圆盘的半径；模块m2：根据卫星定点位置星下点经度计算卫星午夜时刻相对于格林尼治时间的累计秒数；模块m3：根据卫星午夜时刻相对于格林尼治时间的累计秒数，计算午夜时刻太阳高度角；模块m4：根据地影区圆盘半径和午夜时刻的太阳高度角，计算地影期间太阳划过的角度；模块m5：根据地影期间太阳划过的角度，计算卫星地影持续时间、卫星进出影时刻并进行地影区判断。10.根据权利要求1所述的静止轨道卫星在轨自主地影计算方法，其特征在于，在所述模块m1中：太阳圆盘张角为由卫星定点位置看向太阳的圆盘张角，由于地心-日心平均距离远大于卫星-地心平均距离，在计算过程中将该张角简化为地心位置处的太阳圆盘张角，依据太阳平均半径r
sun
及地-日平均距离r
se
计算当前卫星位置处的太阳半张角
度；地球圆盘张角为由卫星定点位置看向地球的圆盘张角，由于地球经圈半径与赤道半径相差极小忽略不计，在计算过程中将地球简化为半径为6371.23km的球体；依据地球平均半径r
e
及星-地平均距离r
ce
计算当前卫星位置处的地球半张角度；地影区圆盘的半径a
se
＝a
sun
+a
eth
＝8.8591度；在所述模块m2中：依据卫星定点位置星下点经度λ
sat
，单位：度，计算卫星相对于格林尼治时间的时间差δt＝240λ
sat
秒；在卫星定点位置午夜时刻，地心-太阳矢量在地球赤道面内的投影、地心-卫星矢量在地球赤道面内的投影之间的夹角为180
°
，即卫星星下点当地地方时为00:00:00，此时的格林尼治时间为：t
g
＝-δt；计算卫星午夜时刻相对于格林尼治时间2000年1月1日12:00:00时刻的累计秒数t0＝t
g0
+t
g
＝t
g0-δt；式中，t
g0
为格林尼治时间当日零点相对于格林尼治时间2000年1月1日12:00:00的累积秒计数；在所述模块m3中：太阳高度角描述为地心-太阳矢量与地球赤道面内之间的夹角，当太阳位于北半球上方时该角度为正，太阳位于南半球上方时该角度为负；春/秋分时刻太阳高度角变化最快，整个地影持续时间内太阳高度角变化值<0.03
°
，若以午夜时刻太阳高度角作为整个地影期间太阳高度角，最大角度偏差<0.015
°
，对地影持续时间计算结果的影响忽略；午夜时刻t0的太阳高度角：β＝arcsin[sin(23.4366)cos(1.14079552584275
×
10-5
t
0-169.699038446726)]；在所述模块m4中：以卫星作为中心参考点，地影过程近似描述为太阳圆盘水平划过地影区圆盘，计算地影期间太阳划过的角度在所述模块m5中：计算太阳圆盘划过地影区圆盘的角速度ω0＝360/86400＝0.004167
°
/s；地影持续时间秒；在所述模块m5中：卫星进地影时刻卫星出地影时刻在所述模块m5中：当t
in
≤t＜t
out
时，卫星处于地影区，否则卫星处于光照区；其中地影区包含了半影区和
全影区；其中，t表示当前时刻。
技术总结
本发明提供了一种静止轨道卫星在轨自主地影计算方法及系统，包括：步骤S1：根据太阳圆盘张角和地球圆盘张角计算地影区圆盘的半径；步骤S2：根据卫星定点位置星下点经度计算卫星午夜时刻相对于格林尼治时间的累计秒数；步骤S3：根据卫星午夜时刻相对于格林尼治时间的累计秒数，计算午夜时刻太阳高度角；步骤S4：根据地影区圆盘半径和午夜时刻的太阳高度角，计算地影期间太阳划过的角度；步骤S5：根据地影期间太阳划过的角度，计算卫星地影持续时间、卫星进出影时刻并进行地影区判断。本发明不依赖地面操控处理，在轨自主完成地影判断，避免由于地面系统漏报造成的卫星能源危机和安全问题。题。题。
技术研发人员：王皓 许海玉 李绿萍 边志强 曾擎 徐凯
受保护的技术使用者：上海卫星工程研究所
技术研发日：2022.03.07
技术公布日：2022/6/28