1.一种x射线脉冲星短时观测高信噪比轮廓获取方法,其特征在于,包括如下:
(1)对光子到达时间序列进行折叠,获取时域波形:
1a)选取一段观测时间t的光子到达时间ti,i=1,2,3...,根据脉冲星星历的自转频率参数,通过脉冲星计时模型计算光子到达时间ti所对应的相位φ(ti);
1b)选用(211,∞)个bin块对光子到达时间ti进行周期折叠或对ti对应的相位φ(ti)进行历元折叠,获得时域波形;
(2)对时域波形进行傅里叶变化获得频谱;
(3)对频谱进行信号与噪声的分离:
3a)选用代价函数,计算频谱的截断点;
3b)根据信号的频域主要集中在低频,而噪声分布在整个频域且幅值较小的不同频域特性,利用截断点对频谱进行截断,选取0点至截断点之间的频谱,使得信号与噪声分离;
(4)恢复高信噪比时域观测轮廓:
4a)对截取的频谱去掉直流分量,即去掉零次频率分量的值;
4b)对去掉零次频率分量的频谱进行傅里叶反变换得到不包含直流分量的时域复信号,对不包含直流分量的时域复信号取实部得到不包含直流分量的时域观测轮廓;
4c)对零次频率分量的值进行傅里叶逆变换得到直流分量;
4d)不包含直流分量的时域观测轮廓加上直流分量得到高信噪比时域观测轮廓。
2.根据权利要求1所述的方法,其中1a)中通过脉冲星计时模型计算光子到达时间ti所对应的相位,公式如下:
其中:φ(t0)是参考时间t0对应的初始相位,f(l)分别是自转频率的l阶偏导数,l取值为2。
3.根据权利要求1所述的方法,其中步骤1b)中选用(211,∞)个bin块对光子到达时间ti进行周期折叠,公式如下:
式中,ts为脉冲星周期,p=t/ts,m为bin块数,其值范围(211,∞),tb=ts/m,cj(ti)为第j个周期内第i个bin块内的光子数,
4.根据权利要求1所述的方法,其中步骤1b)中选用(211,∞)个bin块对光子到达时间ti对应的相位φ(ti)进行历元折叠,公式如下:
式中:m是bin块数,其值范围(211,∞);f为脉冲星的频率,q为脉冲星相位周期数,cj,i为第j个整数相位第i个bin块内的光子数,
5.根据权利要求1所述的方法,其中2)中对时域波形
式中:m是bin块数,其值范围(211,∞);
6.根据权利要求1所述的方法,其中3a)中计算频谱的截断点,通过下式计算:
式中,sk为标准信号傅里叶变换,
7.根据权利要求1所述的方法,其中3b)中的信号与噪声分离,其公式如下:
式中:m是bin块数,其值范围(211,∞);
8.根据权利要求1所述的方法,其中4b)中对去掉零次频率分量的频谱进行傅里叶反变换得到不包含直流分量的时域复信号,其公式如下:
式中:m是bin块数,其值范围(211,∞);n为时域采样点数,q为频域截断点,
9.根据权利要求1所述的方法,其中4b)中得到的不包含直流分量的时域观测轮廓,其公式如下:
式中:n为时域采样点数,
10.根据权利要求1所述的方法,其中4c)中得到的直流分量,其公式如下:
式中:m是bin块数,其值范围(211,∞);
11.根据权利要求1所述的方法,其中4d)中得到的高信噪比时域观测轮廓,其公式如下:
式中:n为时域采样点数,s0为直流分量,