一种基于分区搜索的脉冲信号到达时间快速计算方法
【技术领域】
[0001] 本发明涉及航天器导航领域，特别地，涉及一种基于分区搜索的处理脉冲星信号 的脉冲到达时间的快速计算方法，适用于不同飞行任务的X射线脉冲星导航任务。
【背景技术】
[0002] X射线脉冲星导航是一种新兴的天文导航方法。同传统的天文导航方法相比， X射线脉冲星导航可提供良好的外部参考时间频率基准。同以GPS (Global Positioning System)为例的卫星导航系统相比，X射线脉冲星导航不局限于近地空间，可服务于深空探 测器。同其他波段相比，X射线脉冲星在X射线频段的特征明显，可避免空间各种信号的干 扰。在满足流量和时空分辨率的前提下可降低探测器的体积和功耗。此外，X射线脉冲星 的数目较少，信号结构简单，降低了信号识别的难度。
[0003] 脉冲星的脉冲信号到达时间（Time of Arrival，T0A)是脉冲星导航的基本测量 信息。由于脉冲星信号微弱且不连续，需要长期积累光子到达时间信息才能获得理想精度 的脉冲信号到达时间。假设固定了脉冲星信号的整周模糊度，脉冲信号的到达时间即在一 个周期内变化。目前主流的脉冲到达时间求取方法为历元折合法（Epoch Folding，EF)和 直接利用光子到达时间法Oirect Use of TOAs，DUT)。由于利用了所有光子的信息，DUT 可获得的T0A理论精度高于EF。然而，DUT需要利用每个光子的信息建立似然函数。在相 同观测时间的前提下，DUT的计算复杂度和计算时间都远超EF。目前对DUT的通常解法是 在0到1之间的相位区间内进行网格搜索，即为网格搜索法。假设如下场景：观测流量为 7800ph/m2/s的Crab脉冲星1小时，DUT的搜索步长为0. 0001个相位。在VC6. 0的编译环 境下，利用主频为2. 0GHz的台式机进行500次蒙特卡洛打靶。统计分析表明，DUT的计算 时间约为10小时。显然，该计算时间难以被要求强调实时性的航天器导航接受。
【发明内容】
[0004] 本发明目的在于提供一种基于分区搜索的脉冲信号到达时间快速计算方法，以解 决脉冲信号到达时间计算过程过于冗长，以致影响到航天器导航反应速度的技术问题。
[0005] 为实现上述目的，本发明提供了一种基于分区搜索的脉冲信号到达时间快速计算 方法，包括以下步骤：
[0006] A、以tvigin个相位的步长建立似然函数LF
[0007] B、搜索的单峰区间，所述单峰区间内有且仅有一个最值点；
[0008] C、利用hFurthOT= 0. 5h QHgin个相位的步长建立似然函数LF FurthOT;
[0009] D、搜索LFFurthOT的单峰区间；
[0010] E、判断步骤D的分区结果和步骤B是否相同，如果相同进入步骤F，否则令 =hFurthOT，返回步骤A ;
[0011] F、在步骤E得到的所有单峰区间内，以似然函数为指标函数，搜索每个单峰区间 内的似然函数最值ML』及相应的相位伞 ma!U，j = 1、2、3……N，其中，N为单峰区间的个数；
[0012] G、比较N个ML」的大小，最大的最值ML」对应的相位巾即为脉冲星脉冲信号到 达时间对应的相位；
[0013] H、计算脉冲到达时间TOA=(i>ma!U><P，其中P为脉冲星信号的周期。
[0014] 优选的，tvigin为0.02个相位。
[0015] 本发明具有以下有益效果：
[0016] (1)同网格搜索法相比，本发明通过分单峰区间，在各单峰区间求似然函数最值的 策略，有效地降低了计算量，节省计算成本，提高航天器导航反应速度。
[0017] DUT求解方法的计算复杂度取决于似然函数的计算次数和计算复杂度。考虑到似 然函数的计算复杂度取决于航天器接收到的光子个数，可以在给定航天器接收光子数的前 提下，通过比较不同求解方法计算似然函数的次数来概略比较方法的计算复杂度。
[0018]若要求搜索相位的精度为0.0001个相位，通过网格搜索法求解，需要计算10000次似然函数。利用本发明的方法，令hMgin为0.02个相位，步骤A-B最多需要计算100次似 然函数，步骤C-D最多需要计算200次似然函数。
[0019]1000次的蒙特卡洛数值分析表明：可以保证单峰区间的搜索结果可靠，且单峰区 间个数不超过20个；在每个单峰区间利用一维搜索法，最多30步可以搜到最值。一维搜索 法运行一次需要计算1次似然函数。因此，本发明获得最值点最多需要计算900次似然函 数。同网格搜索法相比，本发明少计算了 9100次似然函数。因此，同网格搜索法相比，本发 明可大幅降低计算量，节省了计算成本。
[0020] ⑵本发明提出了在0到1的相位区间内，先搜索确定单峰区间，再求似然函数最 值的步骤，可有效保证最终获得的结果为全局最优解，避免了直接利用一维搜索算法容易 导致局部最优解的情况。
[0021] DUT的求解问题可以归结为一维搜索问题。然而，现有的一维搜索方法（如黄金分 割法、抛物线法、牛顿迭代法等）的适用区间为单峰区间。然而，DUT对应的似然函数即一 维搜索的目标函数在〇到1的相位区间内为多峰。若不进行单峰区间的划分，直接利用一 维搜索方法，会导致一维搜索方法失效，或仅能找到局部最优解。先划分单峰区间，在各个 单峰区间进行一维搜索，能够获得每个最值点。在最值点集合中，再选择能够获得似然函数 最值的点即可得到全局最优解。
[0022] (3)本发明提出的方法可操作性强，易推广和使用，提高了脉冲星得脉冲信号到达 时间计算的可靠性和实时性，为脉冲星导航的推广提供了条件。
[0023] 除了上面所描述的目的、特征和优点之外，本发明还有其它的目的、特征和优点。 下面将参照图，对本发明作进一步详细的说明。
【附图说明】
[0024] 构成本申请的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实 施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
[0025] 图1是本发明具体应用实例的流程示意图。
[0026] 图2是本发明具体应用实例中使用的Crab脉冲星轮廓函数。
【具体实施方式】
[0027] 以下结合附图对本发明的实施例进行详细说明，但是本发明可以根据权利要求限 定和覆盖的多种不同方式实施。
[0028] 参见图1、图2,如图1所示，本发明的具体步骤如下：
[0029] (1)以hQrigin个相位的步长建立似然函数LF Qrigin〇
[0030] 光子到达航天器服从非齐次泊松分布，在[t山]之间，有I个光子到达的概率为
【主权项】
1. 一种基于分区捜索的脉冲信号到达时间快速计算方法，其特征在于，包括w下步 骤： A、 W Vw。个相位的步长建立似然函数LF hw。； B、 捜索UVw。的单峰区间，所述单峰区间内有且仅有一个最值点； C、 利用hFurthet= 0.化。tw。个相位的步长建立似然函数LFrurthet; D、 捜索LFpurthJ勺单峰区间； E、 判断步骤D的分区结果和步骤B是否相同，如果相同进入步骤F，否则令 hFurther，返回步骤A ; F、 在步骤E得到的所有单峰区间内，W似然函数为指标函数，捜索每个单峰区间内的 似然函数最值MLj及相应的相位d) j = 1、2、3……N，其中，N为单峰区间的个数； G、 比较N个ML斯大小，最大的最值ML对应的相位4 ma。即为脉冲星脉冲信号到达时 间对应的相位； H、 计算脉冲到达时间TOA = (K^jXP，其中P为脉冲星信号的周期。
2. 根据权利要求1所述的一种基于分区捜索的脉冲信号到达时间快速计算方法，其特 征在于，h〇rw。为0. 02个相位。
【专利摘要】本发明提供了一种基于分区搜索的脉冲信号到达时间快速计算方法，包括步骤：先以初始相位的步长建立似然函数后搜索单峰区间，再以0.5个初始相位的步长建立似然函数后搜索单峰区间，若两次搜索的结果相同，则搜索每个单峰区间内的似然函数最值及相应的相位，最大的最值对应的相位即为脉冲星脉冲信号到达时间对应的相位，最后以该相位与脉冲星信号的周期的积得到脉冲到达时间。同网格搜索法相比，本发明通过分单峰区间，在各单峰区间求似然函数最值的策略，有效地降低了计算量，节省计算成本，提高航天器导航反应速度。
【IPC分类】G01C21-20
【公开号】CN104535067
【申请号】CN201510017525
【发明人】郑伟, 王奕迪, 汤国建, 张大鹏
【申请人】中国人民解放军国防科学技术大学
【公开日】2015年4月22日
【申请日】2015年1月14日