法律状态公告日
法律状态信息
法律状态
2022-06-28
公开
发明专利申请公布
技术领域
本发明属于航天导航技术领域，具体涉及一种卫星对地观测过顶时刻快速高精度计算方法。
背景技术
在进行卫星对地观测任务规划时，需要计算卫星对地面目标的可见弧段，从而决定卫星在哪一弧段内可以对特定地面目标进行观测，而在可见弧段中卫星对地面目标的过顶时刻，卫星与目标距离最近，观测效果最好。在通常的多卫星多地面目标观测任务规划时，对任一卫星和任一地面目标的均需计算过顶时刻，而卫星和地面目标的组合数常常在数千组量级，这导致过顶时刻计算的效率至关重要；另一方面，因卫星运动速度很快，对某一地面目标的最佳观测时机稍纵即逝，这导致过顶时刻的计算精度同样重要。现有卫星对地观测过顶时刻计算方法要么采用高精度轨道外推模型，从大步长不断减小到小步长进行搜索，精度高但计算速度很慢；要么采用解析模型进行计算，但与卫星的真实轨道偏差过大，精度不足以用于真实的观测任务规划，其结果仅能在卫星轨道初始设计阶段作为参考。
发明内容
本发明的目的在于提供一种卫星对地观测过顶时刻快速高精度计算方法，用于克服现有技术中计算效率和计算精度难以兼顾的问题，在快速完成过顶时刻计算的同时保证高精度。
本发明所采用的技术方案是，一种卫星对地观测过顶时刻快速高精度计算方法，具体按以下步骤实施：
步骤1，记计算时间区间为[t
步骤2，选取k＝2,3,L,n-1，若L
步骤3，选取i＝1,2,L,m，以
步骤4，选取i＝1,2,L,m，若
步骤5，按时间先后次序对
本发明的特点还在于：
其中Δt范围为0＜Δt≤P/2，P为卫星轨道周期，对于高精确动力学模型，优选的，为减少后续局部优化时的轨道外推耗时，优选的，Δt取为1分钟，对于简化动力学模型，轨道外推耗时很少，为减少时间点个数，优选的，Δt取为P/2；
其中步骤1中计算并记录t
计算各时间点处卫星轨道根数时采用高精度动力学模型，如式(1)所示：
式中，r为卫星矢径，μ为地球引力常数，a
卫星星下点与目标点的地面距离计算方法如式(6)所示：
式中，L为卫星星下点与目标的地面距离，R
其中步骤1中计算各时间点处卫星轨道根数时，若精度要求适中的环境，为获得更快的计算速度，采用仅考虑J
式中，a为轨道半长轴，e为偏心率，incl为轨道倾角，Ω为升交点赤经，ω为近地点幅角，M为平近点角，t为时刻，t
其中步骤3中的一维搜索算法为黄金分割法，卫星轨道外推采用的动力学模型与所述步骤1相同；
其中步骤4中计算卫星对地观测的等效半幅宽的计算过程如下：
式中，W
本发明的有益效果是：
本发明的一种卫星对地观测过顶时刻快速高精度计算方法，适用于快速且高精度计算卫星对地观测过顶时刻的问题。本方法首先计算等间隔时间点处卫星星下点与目标点的地面距离，然后通过对相邻时间点所对应星下点与目标点地面距离的对比确定可能包含距离极小值的一系列较短时间区间，最后分别在各时间区间内通过一维搜索算法求解过顶时刻。本方法首先查找过顶时刻的可能区间，而后利用一维搜索算法在较小可能区间内精确获得过顶时刻，避免了对过顶时刻的大范围重复搜索计算，且本方法对卫星运动模型无限制，可采用于不同精度的动力学模型进行轨道外推计算，具备计算速度快、计算精度高和适应度广的优点，解决了传统卫星对地观测过顶时刻算法计算效率和计算精度难以兼顾、适应度差的问题。
附图说明
图1为本发明的一种卫星对地观测过顶时刻快速高精度计算方法的流程图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明的一种卫星对地观测过顶时刻快速高精度计算方法，如图1所示，具体计算过程包括以下步骤：
步骤1，记计算时间区间为[t
计算各时间点处卫星轨道根数时采用高精度动力学模型，如式(1)所示：
式中，r为卫星矢径，μ为地球引力常数，a
对于精度要求适中的场合，为获得更快的计算速度，也可采用仅考虑J
式中，a为轨道半长轴，e为偏心率，incl为轨道倾角，Ω为升交点赤经，ω为近地点幅角，M为平近点角，t为时刻，t
时间间隔Δt可以在(0,P/2]区间内任意选取，但根据卫星轨道外推动力学模型的复杂程度合理选取可以减少总的计算耗时，优选的，对于高精确动力学模型，为减少后续局部优化时的轨道外推耗时，建议Δt取为1分钟，对于简化动力学模型，轨道外推耗时很少，为减少时间点个数，建议Δt取为P/2；
卫星星下点与目标点的地面距离计算方法如式(6)所示：
式中，L为卫星星下点与目标的地面距离，R
步骤2，对k＝2,3,L,n-1，如果L
步骤3，对i＝1,2,L,m，以
步骤4，对i＝1,2,L,m，如果
卫星对地观测的等效半幅宽的计算方法如式(7)所示：
式中，W
步骤5，按时间先后次序对
