技术特征：
1.一种基于脉冲星辐射信号先验信息的脉冲星自转频率估计方法，其特征在于，包括如下步骤：(1)初始化参数：初始化在观测时间[t
a
,t
b
]内航天器所探测的i个x射线脉冲星光子到达自身的toa序列为其中t
i
为第i个x射线脉冲星光子到达航天器的toa时刻，i＞10；(2)获取x射线脉冲星辐射信号的折叠轮廓：以h为步长对脉冲星自转频率估计范围[f
l
,f
h
]进行遍历，得到p个预设频率值通过用每个预设频率值f
p
对x射线脉冲星光子toa序列进行bin块数为m的周期折叠，得到对应的x射线脉冲星辐射信号的折叠轮廓的集合其中f
l
和f
h
分别为脉冲星自转频率估计范围的下限和上限，f
p
＝f
l
+(p-1)h，p＝1,2,
…
,n+1，p＝n+1，y
p
为f
p
对应的包含m个bin块的x射线脉冲星辐射信号的折叠轮廓，m＝2
s
,s∈n
+
，n
+
为正整数集；(3)获取x射线脉冲星辐射信号标准轮廓与折叠轮廓的时域互相关函数的频谱：以x射线脉冲星辐射信号的标准轮廓x为先验信息，对x、每个预设频率值f
p
所对应的x射线脉冲星辐射信号的折叠轮廓y
p
分别进行m点快速傅里叶变换，得到x的频谱f
x
、折叠轮廓的频谱的集合将f
x
和每个f
yp
的共轭结果相乘，得到x射线脉冲星辐射信号标准轮廓x与折叠轮廓的时域互相关函数的频谱的集合其中(4)获取每个预设频率对应的代价函数值：根据设定的x射线脉冲星辐射信号标准轮廓的频谱的截断点q，对每个时域互相关函数的频谱f
rp
中0点至q之间的频谱分量进行快速傅里叶逆变换，得到频谱截断的时域互相关函数的集合并将每个的最大值作为每个预设频率值f
p
对应的代价函数值，其中q∈n
+
,q∈[15,45]，τ∈n
+
,τ∈[1,m]；(5)获取x射线脉冲星自转频率的估计结果：对以p个预设频率值为自变量、以p个预设频率值所对应的p个代价函数值为因变量所组成的代价函数进行单峰值拟合，并将代价函数最大值所对应的频率值作为脉冲星自转频率的估计值应的频率值作为脉冲星自转频率的估计值2.根据权利要求1中所述方法，步骤(2)中用预设频率值f
p
对x射线脉冲星光子toa序列
进行周期折叠，得到一个对应于预设频率值f
p
的x射线脉冲星辐射信号的折叠轮廓y
p
，实现步骤为：(2a)将x射线脉冲星光子toa序列按照预设周期值t
p
分为k个子序列，其中t
p
＝1/f
p
，k＝(t
b-t
a
)/t
p
；(2b)对第k个子序列进行m等分，统计第k个子序列中第m个bin块内的x射线脉冲星光子数c
k
(t
m
)，其中m∈[1,m]，k∈[1,k]；(2c)将k个子序列中第m个bin块内的x射线脉冲星光子数c
k
(t
m
)相加，得到一个对应于预设频率值f
p
的x射线脉冲星辐射信号的折叠轮廓y
p
：
技术总结
本发明提出了一种基于脉冲星辐射信号先验信息的脉冲星自转频率估计方法，实现步骤为：(1)初始化参数；(2)获取X射线脉冲星辐射信号的折叠轮廓；(3)获取X射线脉冲星辐射信号标准轮廓与折叠轮廓的时域互相关函数的频谱；(4)获取每个预设频率对应的代价函数值；(5)获取X射线脉冲星自转频率的估计结果。本发明由于选用了脉冲星辐射信号的标准轮廓作为先验信息，以周期折叠轮廓和标准轮廓的时域互相关函数的最大值作为代价函数值，不同于现有技术中仅利用周期折叠轮廓计算统计量的方式，弥补了在短时观测条件下观测数据信噪比低对统计量的计算所造成的负面影响，从而提高了在短时观测条件下脉冲星自转频率的估计精度。观测条件下脉冲星自转频率的估计精度。观测条件下脉冲星自转频率的估计精度。
技术研发人员：方海燕 姚云峰 张魁 孙泽阳 王正品 冯龄慰 张同
受保护的技术使用者：西安电子科技大学
技术研发日：2022.03.28
技术公布日：2022/7/1