[发明专利]一种脉冲星周期跃变的检测方法

权利要求书

1.一种脉冲星周期跃变的检测方法，其特征在于，该检测方法基于KL距离标准差获得，具体包括如下步骤：

第一步，将探测数据分成M段折叠轮廓，M取大于等于1的自然数；

第二步，计算第N组折叠轮廓与第一组折叠轮廓的KL距离，N为大于等于2且小于等于M的自然数；

第三步，计算KL距离的标准差，并判断是否超过门限范围，若超过门限范围，则判定为脉冲星周期跃变，否则，则判定为脉冲星周期未发生跃变；其中，采用下列公式计算KL距离的标准差σ_KL：

其中：KL_i为第i组折叠轮廓与第一组折叠轮廓的KL距离。

2.根据权利要求1所述的脉冲星周期跃变的检测方法，其特征在于，第一步中利用历元折叠法将探测数据分段折叠轮廓。

3.根据权利要求2所述的脉冲星周期跃变的检测方法，其特征在于，将探测数据分段折叠轮廓，具体为：

假设观测周期为T_obs，含有N_p个脉冲星周期，即N_p＝floor(T_obs/P)；脉冲星周期P进一步划分为N_b个封装段，每个封装段的长度为T_b；

历元折合包括以下步骤：

步骤1.1、将记录在后续周期内的光子TOA序列折合到第一个周期；

步骤1.2、计算每个封装段中的光子数c_i，其中：i＝1,2,…,N_b；

步骤1.3、通过归一化光子数，得到表达式1)的折叠轮廓

其中，c_j(T_i)为第j个周期中折合入第i个封装段的光子数，c_i是第i个封装段内的光子总数，T_i为第i个封装段的中间时刻。

4.根据权利要求1所述的脉冲星周期跃变的检测方法，其特征在于，所述第二步中计算第N组折叠轮廓与第一组折叠轮廓的KL距离KL_N具体为：

设p(x)和q(x)是两个随机变量X₁与X₂的概率密度函数，随机变量X₁与X₂之间的KL距离定义为表达式2)：

对于离散型分布X₁与X₂有表达式3)：

5.根据权利要求1所述的脉冲星周期跃变的检测方法，其特征在于，根据拉依达准则，如果第n+1个KL距离导致其标准差增大至前n个KL距离标准差的三倍，n为大于等于1的自然数，则判定第n+1个KL距离对应的折叠轮廓出现异常，即发生周期跃变；否则，则判定第n+1个KL距离对应的折叠轮廓未出现异常，即未发生周期跃变。