[发明专利]基于超越不等式的卫星覆盖预报方法

摘要

基于超越不等式的卫星覆盖预报方法，（1）构造函数f1(x1,x2)、f2(x1,x2)；（2）构造一个空的序列vlist，用于记录0/1；构造一个空的有序序列tlist，用于记录变量值，当预报时间范围时，此变量值代表时间，当预报侧摆范围时，此变量值代表侧摆角度，有序序列tlist按照变量值从小到大排序；tlist与vlist始终长度一致且值一一对应；设定计数器counter，初始值为1；（3）将搜索区间[A,B]分成n等分，每一等分用区间[a,b]表示，counter值加上n-1；对每一等分[a,b]进行覆盖预报计算，覆盖预报计算大体分为两个部分，第一部分是递归细分，构建0/1序列，也就是F(a,b)符号的序列。第二部分是根据0/1序列找寻不等式的解。

主权项

基于超越不等式的卫星覆盖预报方法，其特征在于步骤如下：（1）构造函数f1(x1,x2)、f2(x1,x2)；（2）构造一个空的序列vlist，用于记录0/1；构造一个空的有序序列tlist，用于记录变量值，当预报时间范围时，此变量值代表时间，当预报侧摆范围时，此变量值代表侧摆角度，有序序列tlist按照变量值从小到大排序；tlist与vlist始终长度一致且值一一对应；设定计数器counter，初始值为1；（3）将搜索区间[A,B]分成n等分，每一等分用区间[a,b]表示，counter值加上n‑1；对每一等分[a,b]同时进行如下处理：（3.1）求解函数F(a,b)，当预报时间范围时，令F(a,b)=f1(a,b)，当预报侧摆范围时，F(a,b)=f2(a,b)；（3.2）判断F(a,b)值的符号，若F(a,b)>0，转步骤（3.3）；否则，在有序序列tlist中按照a的值二分搜索到位置p；往序列vlist中的位置p插入0，往序列tlist的位置p插入a，counter值减1，转步骤（4）；（3.3）判断等分[a,b]的区间长度b‑a是否大于等于阈值，若是，则转步骤（3），将当前等分[a,b]作为搜索区间[A,B]，重新进行等分及后续处理；否则在有序序列tlist中按照a的值二分搜索到位置p；往序列vlist中的位置p插入1，往序列tlist的位置p插入a，counter值减1，转步骤（4）；（4）判断counter值是否等于0，若等于0，则转步骤（5）；否则，等待；（5）初始化一个二元组列表ulist，用于记录不等式的解，并令搜索起始位置为q，q=1；（6）在序列vlist中从位置q开始，搜索由0到1的跳变，找到0的位置p；从tlist中找到位置p对应的变量值，记录到临时变量s；（7）在序列vlist中从位置p开始，搜索1到0的跳变，找到0对应的位置p，从tlist中找到位置p对应的变量值e，与之前记录的变量s，组成一个二元组[s,e]作为不等式的一个解记入列表ulist；（8）判断是否搜索若到序列末尾，若是，则停止，否则令搜索起始位置q=p，转步骤（6）重复进行，二元组列表ulist中的解[s,e]为所有满足预报要求的时间范围或者侧摆角度范围。