基于Bayesian估计的X射线脉冲星微弱信号检测

累积轮廓、流量和周期是X射线脉冲星辐射信号的三个重要特征,将其应用于X射线脉冲星信号检测中,提出了一种基于Bayesian估计的X射线脉冲星周期辐射信号时域检测方法.该方法以非脉冲区噪声观测为先验知识,利用X射线脉冲星辐射信号的泊松分布模型推导了信号概率密度分布函数,以该函数的累积分布函数为判据,对X射线脉冲星微弱信号进行检测,并提取位相偏移量.利用仿真数据和RXTE卫星的实测数据进行实验验证,结果表明:本文方法性能优于同类的基于高斯分布模型的检测方法,在检测信号的同时能在一定精度下给出信号位相偏移值. 关键词： 脉冲星 Bayesian估计 位相测量 时域检测     

非等间隔计时数据的X射线脉冲星周期快速搜索算法

利用X射线脉冲星信号的计时信息实现航天器深空导航具有重要意义,精确的时间模型(周期及周期导数)是获得高精度自主导航的基础.针对X射线光子并非等间隔地到达探测器的特点,在对比分析两种非等间隔数据处理方法(x2量评估方法和Lomb算法)的基础上,提出了用Lomb算法估计脉冲星周期初值,进一步使用x2量评估方法精化周期值的方法.同时运用快速傅里叶变换算法的思想对Lomb算法进行改进,较大地提高了程序的运算效率.最后利用X射线源仿真系统的计时数据对算法进行试验分析,搜索出准确的脉冲星周期和折叠出正确的脉冲轮廓.     

一种新的X射线脉冲星信号模拟方法

由于无法在地面直接观测到X射线脉冲星信号, 且空间飞行探测耗时长, 成本高, 模拟具有真实物理特征的X射线脉冲星信号对X射线脉冲星信号处理方法及导航方案的验证具有重要意义. 本文提出了一种利用太阳系质心处脉冲星信号模型和航天器轨道信息, 建立航天器处实时光子到达速率函数, 再利用尺度变换法产生光子到达时间序列的脉冲星信号模拟新方法. 该方法真实还原了脉冲星信号的频率缓变特性, 考虑了动态探测环境下的相对论效应, 且避免了现有模拟航天器处光子序列方法中的迭代过程, 产生非齐次泊松光子到达时间序列的运算量也低于常用的齐次泊松过程筛选法和反函数递推法. 数值仿真从频率特性、流量特性、轮廓相似度及与实测数据的接近程度四个方面验证了该模拟方法的正确性.     

一种基于最优频段的X射线脉冲星累积轮廓时延估计方法

为提高X射线脉冲星导航中累积脉冲轮廓的时间延迟估计精度, 分析了X射线脉冲星累积脉冲轮廓的频谱特性和现有Taylor快速傅立叶变换时延估计算法的缺陷, 提出了一种基于最优频段的累积轮廓时延估计算法, 并通过建立不同信噪比下时延估计误差与所采用频段之间的关系以确定最优频段. 数值及实测数据实验结果表明: 在短时观测或光子流量较小时, 该算法优于常用的近似最大似然 (FAML)、相关 (CC)、最小二乘 (NLS) 及加权最小二乘 (WNLS) 方法; 在观测时间较长或光子流量较大时, 该算法的估计精度与CC及NLS方法相当, 但其运算量低于NLS, FAML 及WNLS方法. 本文所提算法适用于短时观测脉冲轮廓或低流量脉冲星的高精度时延估计.     

基于X射线脉冲星导航的模拟调制仿真源研究

为实现X射线脉冲星导航系统在地面的演示验证, 提出并设计了一种高精度X射线脉冲星仿真源. 该仿真源由模拟调制信号发生器和栅控X射线球管组成. 模拟调制信号发生器根据脉冲星标准脉冲轮廓数据和栅控球管特性曲线, 利用直接数字频率合成技术产生模拟调制信号, 加载在X射线球管控制栅极, 从而控制轰击阳极靶的高速电子数目来实现X射线的调制, 产生与脉冲星标准脉冲轮廓高吻合度的X射线光子统计分布. 在X射线脉冲星地面仿真系统中对Crab脉冲星仿真源的性能进行测试, 测试结果为: 采集脉冲轮廓与标准脉冲轮廓时域相关度和频率相关度分别达到0.9774和0.9853, 辐射流量为1. 90 ph·cm^-2·s^-1, 脉冲辐射流量与总辐射流量之比为76.15%, 脉冲半宽度为1.879 ms. 结果表明: 该仿真源具有符合度高, 成本低, 操作简单灵活, 对X 射线脉冲星导航关键技术的攻关具有重要的意义. 关键词： X射线脉冲星导航 仿真源 栅控X射线球管 直接数字频率合成     

基于TOA信息的脉冲星信号周期估计方法

针对脉冲星信号周期估计需逐步搜索、计算量大等问题,提出脉冲星信号周期直接估计算法.根据周期不确定误差对时域上脉冲到达时间(Pulse Time of Arrival,TOA)估计方法的影响,推导了基于TOA信息的脉冲星信号周期估计方法的数学模型.该方法将一组观测到的脉冲光子到达时间序列(Photon Time of Arrivals,TOAs)进行等时间间隔分段,在时域上对每段TOAs进行相位估计获得相应的TOA信息,并根据建立的TOA信息与周期误差的关系采用最小二乘原理估计脉冲星信号周期.对周期准确性的评价不同于传统的依赖折叠轮廓的好坏,而是以TOA与时间图像的斜率作为依据,这一依据更直观、易于理解.理论分析与物理、计算机仿真数据结果表明,本文所提的周期估计方法能够根据短时间脉冲光子到达时间序列获得高精度和高分辨率的脉冲信号周期,有利于X射线脉冲星导航的工程应用.     

空间X射线观测确定脉冲星星历表参数精度分析

脉冲星星历表维持着脉冲星导航所需的时空基准, 其精度直接影响着航天器导航定位结果, 是脉冲星导航系统的基本要素. 本文分析了脉冲星空间观测精度的估计方法, 探索性地研究了基于空间X射线观测获取星历表参数的可行性. 通过建立星历表参数拟合模型, 采用大样本重复事件仿真分析了空间X射线观测精度, 研究了星历表参数确定精度与观测精度、观测时间及观测频次的关系. 研究结果表明, 空间X射线观测可以确定脉冲星星历表参数, 但受限于脉冲星信号特征及探测器技术水平, 当前高精度导航用脉冲星星历表难以通过空间X射线观测手段获得, 可通过地面射电观测技术较好地建立与维持, 提出了推进我国大口径射电望远镜建设的建议.     

X射线脉冲星导航半物理仿真实验系统研究

由于费用巨大,X射线脉冲星导航初步研究阶段不可能进行空间搭载实验.为此,文章设计了一种X射线脉冲星导航半物理仿真实验系统.该系统由脉冲星信号模拟部分和导航参数解算部分组成.模拟部分采用非齐次泊松过程对光子到达太阳系质心的时间建模,时间转换后,模拟X射线探测器观测脉冲星时输出的脉冲信号.导航参数解算部分利用Delta-Correction方法解算模拟信号中蕴含的导航信息.该系统可同时模拟4颗脉冲星信号,成本低,精度高,可对脉冲星导航的信号处理和参数解算过程进行光子级仿真研究,并为后续空间搭载实验中原理样机的设计提供一定的参考. 关键词： X射线脉冲星 导航 仿真实验系统 非齐次泊松过程     

基于TOA信息的脉冲星信号周期估计方法（英文）

针对脉冲星信号周期估计需逐步搜索、计算量大等问题,提出脉冲星信号周期直接估计算法.根据周期不确定误差对时域上脉冲到达时间(Pulse Time of Arrival,TOA)估计方法的影响,推导了基于TOA信息的脉冲星信号周期估计方法的数学模型.该方法将一组观测到的脉冲光子到达时间序列(Photon Time of Arrivals,TOAs)进行等时间间隔分段,在时域上对每段TOAs进行相位估计获得相应的TOA信息,并根据建立的TOA信息与周期误差的关系采用最小二乘原理估计脉冲星信号周期.对周期准确性的评价不同于传统的依赖折叠轮廓的好坏,而是以TOA与时间图像的斜率作为依据,这一依据更直观、易于理解.理论分析与物理、计算机仿真数据结果表明,本文所提的周期估计方法能够根据短时间脉冲光子到达时间序列获得高精度和高分辨率的脉冲信号周期,有利于X射线脉冲星导航的工程应用.     

提升小波的X射线脉冲星信号降噪EI北大核心CSCD

对于短时间内X射线脉冲星信号形成的脉冲轮廓的信噪比低,提出基于提升小波的X射线脉冲星信号的降噪方法。分析了基于提升小波的信号降噪理论模型。在周期叠加的基础上,利用探测设备及空间环境的先验信息建立预降噪模型,并对X射线脉冲星的脉冲轮廓做降噪预处理。通过选择小波和自适应的阈值函数,利用提升小波对脉冲轮廓降噪。通过对多组罗西X射线计时探测器(RXTE)卫星探测的实测X射线脉冲星数据的仿真分析,结果表明,降噪后的脉冲轮廓信噪比得到提高,且算法有更快的处理速度。     

具有多物理特性的X射线脉冲星导航地面验证系统

导航地面验证是X射线脉冲星导航研究必不可少的环节.针对导航算法验证需要真实连续的脉冲星信号的需求,同时避免X射线调制及探测难度大、成本高的问题,提出了一种基于可见光源的X射线脉冲星导航地面验证系统.该系统利用太阳系质心处脉冲星信号模型和航天器轨道信息,建立航天器处实时光子到达速率函数,再通过硬件系统转换成电压信号,利用该电压控制线性光源输出,最后经衰减、探测及甄别后获得航天器处的实时光子到达时间序列.该时间序列不仅具有导航脉冲星的轮廓特性、自转特性,还包括空间传播时间效应及宇宙X射线背景.本系统利用半物理装置对可见光进行调制及衰减,实时判断轨道各位置处导航脉冲星的可见性,实现X射线脉冲星信号传播过程的模拟.该系统提供四路可控输出信号,支持多种导航模式的验证.仿真系统的性能分析和功能验证结果表明,该系统具有良好的性能,可提供真实便捷的地面验证环境.     

基于S变换的脉冲星辐射脉冲信号检测

为了提高脉冲星辐射脉冲信号的检测速度和在低信噪比下的检测效果, 提出了一种基于S变换的脉冲星辐射脉冲信号检测算法. 文中证明了高斯白噪声S变换域功率谱服从自由度为2的卡方分布, 基于此对累积信号S变换域功率谱进行阈值处理,累加阈值处理后的时频功率谱作为统计量进行检测. 此外阈值处理后的功率谱也可用来测量脉冲星信号的时间延迟. 仿真实验验证了本文算法的有效性,其检测性能优于同类的基于高斯分布模型的检测算法, 同时还可以在一定精度下给出脉冲星信号的时间延迟值. 关键词： 脉冲星 卡方分布 S变换域检测 时延测量     

基于X射线脉冲星导航的地面模拟系统研究

为进行X射线脉冲星导航的关键技术研究,搭建了基于X射线脉冲星导航的地面模拟系统.地面模拟系统由模拟X射线脉冲源、基于微通道板的高灵敏X射线光子探测器、电荷灵敏前放和主放电路、时间测量单元、X射线脉冲轮廓构造及X射线脉冲到达时间测量系统组成.该模拟系统可在地面模拟X射线脉冲星导航的星源的强度、周期及脉冲轮廓,实现对X射线脉冲星单光子到达时间的记录,构造X射线脉冲星脉冲轮廓,计算X射线脉冲到达时间.描述了基于X射线脉冲星导航的地面模拟系统的组成和工作原理,报道了基于X射线脉冲星导航的地面模拟系统的初步结果. 关键词： X射线脉冲星导航 微通道板光子探测器 脉冲轮廓     

X射线脉冲星导航系统模拟光源的研究

介绍了X射线脉冲星导航地面模拟光源研究的必要性及非伺服的机械调制方法所存在的问题和缺陷, 提出了基于栅控X射线球管的X射线脉冲星辐射脉冲模拟方法, 通过电子光学设计计算, 对栅控X射线管的电极结构进行设计优化, 研制了栅控X射线管和脉冲星模拟光源装置.实验测试了栅控球管的性能, 测试结果与理论计算结果基本相符, 实现了对X射线的调制; 通过基于FPGA的直接数字频率合成方法, 产生脉冲星的任意形状脉冲轮廓电压信号, 加载至球管控制栅极, 并对其出射脉冲轮廓进行测试, 结果表明产生的X射线脉冲轮廓逼真程度在95%以上, 模拟源频率稳定度约为2×10^-11. 关键词： 脉冲星导航 X射线球管 栅极控制     

一种脉冲星信号模拟新方法

在脉冲星导航地面实验中,为解决地面无法接收脉冲星辐射的X射线信号的问题,必须在实验系统中模拟脉冲星辐射信号. 本文提出两种全新脉冲星信号模拟方法. 第一种方法基于统计物理模型,该方法克服了常用泊松模型模拟方法中仿真适应性窄、仿真速度慢和时间分辨率低等不足,实现了对任意随机过程的脉冲信号进行快速高时间分辨率的仿真;第二种方法基于优化统计模型,该方法针对非齐次泊松过程信号的特性进行优化,仿真速度比常用泊松模型模拟方法提高至少30 倍,对于低流量脉冲星,甚至提高4个数量级,可实现纳秒级时间分辨率仿真. 关键词： 脉冲星导航 非齐次泊松过程 X射线光子信号     

高稳定度X射线脉冲星信号模拟

X射线脉冲星导航是一种完全自主的导航方式, 在深空乃至行星际空间具有潜在的工程应用价值.由于空间飞行试验系统复杂, 成本巨大, 在实验室环境下高精度地模拟X射线脉冲星信号对数据处理方法和导航方案的验证具有重要意义. 针对当前机械转盘式模拟系统中时间稳定度和轮廓精度的不足, 提出了一种通过产生的轮廓电压信号直接控制可见光光源, 再利用衰减获得光子流, 最后经单光子探测和处理电路输出光子到达时间序列的模拟新方法.该方法实现成本低, 支持任意X射线脉冲星信号的模拟, 且具有高时间稳定度和轮廓精度. 详细地讨论了该方法的原理和涉及的关键技术, 搭建了X射线脉冲星信号模拟系统, 并进行了实验. 实验结果表明: 该系统大幅提高了X射线脉冲星信号的模拟效果, 将模拟脉冲星自转周期的稳定度从现有的10^-4提高到10^-9; 当探测器面积为1 m², 探测能谱范围为2–10 keV, 积分时间为1200 s时, 模拟的PSR B1509-58 观测脉冲轮廓与标准脉冲轮廓的相关系数达到了0.993. 关键词： X 射线光子信号 时间稳定度 观测脉冲轮廓 导航算法     

基于循环平稳信号相干统计量的脉冲星周期估计新方法

讨论了脉冲星的周期估计在脉冲星搜索中的重要性.将脉冲星观测信号建模为二阶循环平稳模型,仿真验证了其合理性.在此基础上,提出了一种基于双谱相干统计量的周期估计新方法,并给出了稳健有效的周期搜索策略以消除其对数据量的敏感性.分别对单脉冲脉冲星(PSR J0437-4715)的实测信号和双脉冲脉冲星(PSR B1821-24)的仿真信号进行了周期估计,实验表明,相比较于常用的傅里叶频谱法,该时域方法直观、有效,在低信噪比情况下仍具有很好的性能,且适用于非连续观测数据.尽管方法运算量较大,但仍然可为微弱脉冲星的周 关键词： 脉冲星 脉冲星搜索 周期估计 循环平稳信号     

一种脉冲星周期的时域估计新方法

脉冲星信号的周期估计在脉冲星搜索、观测以及定时研究中具有重要作用.在分析脉冲星周期频域识别方法的基础上,基于脉冲星观测信号的二阶循环平稳模型,提出了一种周期估计的时域新方法-----最大相关方差搜索法.通过对观测信号按假定周期(可变)分组,计算分组后不同时刻的相关系数,进一步得到系数的方差,然后寻找不同周期下的方差最大值,即可得到脉冲星的辐射周期估计.分别对单峰脉冲星(PSRJ0437-4715)的实测信号、双峰脉冲星(PSRB1821-24)的仿真信号以及ATNFEPN数据库中的多颗脉冲星数据进行了周期估计,均得到了令人满意的结果.对众多脉冲星的计算结果表明,与傅里叶频谱法相比,该时域方法对噪声不敏感,尤其在低信噪比情况下仍表现出较优越的性能;与双谱相干统计量方法^[6]相比,运算复杂度更低;尽管其在处理非连续观测数据时能力不足,但仍然可为微弱脉冲星的实时周期估计提供有力参考.     

X射线脉冲星累积轮廓建模及信号辨识

为减少X射线脉冲星信号辨识所需的观测时间, 提出一种基于轮廓光子分布统计量的辨识算法. 用该方法分析了观测数据按不同周期累积所获轮廓的光子分布差异, 按累积周期正确与否, 对轮廓进行分类建模. 基于两类轮廓模型, 研究了光子分布统计量的性质差异, 以此进行信号辨识. 利用罗西X射线时变探测器数据进行了辨识实验, 结果表明该算法有效地减少了辨识所需的观测时间, 且不存在误检问题.     

多谐波脉冲星信号时延估计方法

针对脉冲星导航技术中延时估计这一关键问题, 提出了频域上直接使用脉冲星信号测量到达时间集合进行时延估计的方法——多谐波脉冲星信号时延估计(MHSPE)方法. 该方法建立在频域上相位时延的极大似然估计的基础上, 通过高次谐波对脉冲星观测信号提取出各谐波相位的极大似然估计, 然后取频谱上各谐波的幅值进行归一化作为各谐波相位的权值, 最后取各谐波相位的加权平均作为该时刻的相位估计. 理论上证得MHSPE算法对相位的估计是无偏、一致的, 相比于频域上一次谐波的极大似然估计, MHSPE方法的信噪比随谐波数m的增加而增加, 当各谐波幅值相同时, 信噪比可提高m_1/2倍; 与脉冲星信号时延的克拉美罗界比较, 脉冲星信号时域的导数在频域上的反映就是各谐波分量的数量, 因此随着谐波次数的增加脉冲星信号时延估计可极大趋近克拉美罗界. 采用RXTE航天器对Crab脉冲星的实测数据检验MHSPE方法的性能, 实验结果表明, 针对低信噪比的脉冲星信号, MHSPE可获得高精度的相位估计, 随观测时间增加, 估计精度快速收敛于克拉美罗界.