“天枢Ⅱ号”X射线脉冲星导航动态模拟系统及实验验证_盛立志

《“天枢Ⅱ号”X射线脉冲星导航动态模拟系统及实验验证_盛立志》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

航空学报Feb.152023Vol.44No.3ActaAeronauticaetAstronauticaSinicaISSN1000-6893CN11-1929/Vhttp://hkxb.buaa.edu.cnhkxb@buaa.edu.cn引用格式：盛立志，郑伟，苏桐，等.“天枢Ⅱ号”X射线脉冲星导航动态模拟系统及实验验证［J］.航空学报，2023，44（3）：526656.SHENGLZ，ZHENGW，SUT，etal.GroundtestbenchforX-raypulsarnavigationdynamicsimulation［J］.ActaAeronauticaetAstronauticaSinica，2023，44（3）：526656（inChinese）.doi：10.7527/S1000-6893.2022.26656“天枢Ⅱ号”X射线脉冲星导航动态模拟系统及实验验证盛立志1，*，郑伟2，苏桐1，张大鹏2，王奕迪2，杨向辉1，徐能1，李治泽21．中国科学院西安光学精密机械研究所，西安7101192．国防科技大学空天科学学院，长沙410073摘要：X射线脉冲星导航真实运行场景的模拟验证手段尚不够完善，限制了导航和探测理论的发展及工程化应用的实现。提出了基于脉冲星辐射特性和航天器轨道模型的X射线脉冲星动态信号模拟方法，设计了能实现多种实验场景模拟的X射线脉冲星模拟源，并根据脉冲星辐射特性构建了“天枢Ⅱ号”X射线脉冲星导航地面实验系统，可以高质量复现脉冲星导航空间运行场景。设计了静态及动态不同类型的实验对实验系统的性能进行了验证测试，针对PSRB0531+21、PSRB1937+21两颗脉冲星，开展了静态模拟实验，获得的脉冲轮廓相似度分别为99.5%、99.1%；开展了200km轨道高度的动态近地圆轨模拟实验，PSRB0531+21、PSRB1937+21的周期测试值和理论值的偏差分别为38451、350ps，还原到SSB（SolarSystemBarycenter）处的轮廓相似度分别为99.86%、99.99%。实现了椭圆轨道的超实时仿真实验，仿真时长可压缩50%，轮廓相似度为99.89%。实现了基于霍曼变轨模型的轨道机动模拟，周期变化的测试值与理论值的偏差标准差为637ps。该地面实验系统性能稳定，可以满足不同类型的X射线脉冲星导航模拟实验需求。关键词：X射线脉冲星导航；动态实验；X射线脉冲源；X射线探测器；半实物验证中图分类号：V11文献标识码：A文章编号：1000-6893（2023）03-526656-11X射线脉冲星导航（X-rayPulsarNavigation，positionExplorer，NICER），NICER搭载了56个XPNAV）是以高稳定的脉冲星辐射信号的时空口径为10cm的X射线聚焦镜，开展了X射线脉信息为参考，从而实现航天器的位置、速度信息确冲星计时研究及脉冲星导航验证研究，并通过对定的新型导航方式，导航系统运行过程中不需要毫秒脉冲星的观测得出了导航精度可达10km的［7］地面观测站的校准，具有完全自主导航的能结论。中国于2016年11月发射了XPNAV-1，［1-6］力。目前，国内外已围绕X射线脉冲星导航这是中国首颗脉冲星导航试验卫星。XPNAV-1XPNAV进行了大量的理论及初步实验研究。美中搭载了准直型和聚焦型两种X射线探测器，其国国家航空航天局（NASA）于2017年发射了中中，聚焦型探测器有效探测能段为0.5~10keV，2［8］子星内部结构探测器（Neutron-starInteriorCom⁃有效面积为30cm。卫星对PSRB0531+21脉收稿日期：2021-11-15；退修日期：2021-12-10；录用日期：2022-01-27；网络出版时间：2022-03-1909：27网络出版地址：https：//hkxb.buaa.edu.cn/CN/Y2023/V44/I3/526656基金项目：国家自然科学基金（61901470）；中国科学院“西部之光”人才计划（XAB2020YN13）；强脉冲辐射环境模拟与效应国家重点实验室开放课题（SKLIPR2021）*通信作者．E-mail：shenglizhi80@163.com526656-1

1航空学报冲星进行了长时间的观测并获得了有效的数据，倍增管PMT（PhotomultiplierTubes），通过减实现了对X射线脉冲星“看得见”的目标。此外，光片将光子流量控制在单光子探测范畴，实现［14-15］中国于2016年9月发射的天宫二号上的伽马暴偏脉冲星信号地面模拟实验研究。中国科学［9］振探测仪（POLAR）和2017年6月发射的用于院西安光学精密机械研究所对X射线模拟源及空间天文研究的硬X射线调制望远镜（HardX-探测器等关键技术进行了大量研究，提出了［10］rayModulationTelescope，HXMT）均对X射一种用于产生脉冲X射线的栅控X射线管，可线脉冲星进行了观测，并开展了导航试验空间验以对任意脉冲星的轮廓进行模拟，开展了基于［16-18］证。虽然上述空间观测和验证已经获得不少成X射线脉冲星导航的静态模拟实验。目前果，但XPNAV工程化应用仍有一定的距离，完善文献报道的基于X射线波段的半实物实验系统的脉冲星数据库的建立、脉冲星信号降噪及高精开展的大多局限于静态实验，并没有包含航天度脉冲到达时间TOA（TimingofArrival）估计、器轨道运动特性。轻型化高效率X射线探测系统研制等关键技术仍本文面向脉冲星导航应用研究的需求，采用需要深入研究。栅控X射线管的方法开展基于脉冲星辐射特性X射线无法穿透大气层，导航工程化之前必和航天器轨道模型的X射线脉冲星动态信号模须要开展X射线脉冲星地面模拟实验研究，为拟，构建地面实验系统，对X射线脉冲星静态辐导航相关理论提供实验验证平台。NASA的戈射特性及包含各种动态效应的动态脉冲信号进达德空间飞行中心GSFC开发了X射线导航实行高精度模拟和效果评估。验平台（X-rayNavigationLaboratoryTestbed，1“天枢Ⅱ号”X射线脉冲星导航动态模GXLT），X射线模拟源采用调制X射线源拟系统（ModulatedX-RaySource，MXS），探测器为硅漂移探测器SDD（SiciconDriftDetector）。“天枢Ⅱ号”X射线脉冲星导航地面实验系统MXS主要是通过紫外LED调制光电阴极，产生的组成如图1所示，第1个模块主要实现对X射的光电子经过电子倍增器件后在阳极高压作用线脉冲星及背景噪声的模拟。PC软件控制终下轰击阳极靶材，从而产生调制的X射线。端、动态信号发生器、栅控信号调制器组成了脉GXLT可以进行3个层次的仿真：纯软件层面、冲星模拟控制器，用该控制器驱动栅控X射线管不模拟光子到达过程的仿真，纯软件模拟光子发射所要模拟的X射线脉冲信号，模拟不同的工到达过程的软件层面闭环模拟实验，以及全物作状态，利用GPS（GlobalPositionSystem）驯服［11］理闭环仿真实验。中国也开展了脉冲星导航的原子钟来给出精确时标；直流X射线管和二次地面模拟实验验证技术的相关研究，并取得了靶转台构成了单能X射线模拟器，在上位机软件一定的进展。国防科技大学于2011年研制了国控制终端及驱动电源的作用下可以实现对背景X内首套X射线脉冲星导航半实物仿真系统，可射线噪声的模拟。第2个模块为真空实验平台，以利用探测器实时探测X射线光子信号，将光主要由真空计、真空管道、分子泵、机械干泵等构子到达时间通过数据采集传入导航计算机并进成，模拟太空中的真空环境。X射线聚焦光学行处理，实现TOA估计和导航滤波算法，从而（NestedX⁃rayFocusingOptics，NXFO）SDD探［12］实现X射线脉冲星导航的半实物仿真。中国测器、制冷及电压控制模块、SDD电子学等构成航天科技集团五院建立了脉冲星导航探测器地了聚焦型探测系统，用于接收光子信息；接收到面标定装置，用于探测器脉冲轮廓还原能力和的光子信息在信号处理计算机中进行处理，完成［13］探测效率、能量分辨率等指标的标定。西安相关实验验证工作。系统实物如图2所示。电子科技大学基于可见光搭建了脉冲星地面验X射线脉冲星模拟源由脉冲星模拟控制终端证系统，模拟源采用线性较好的蓝色发光二极和脉冲X射线发射源组成，实现对X射线脉冲星管LED（LightEmittingDiode），探测器为光电的模拟。脉冲星模拟控制终端用于实验参数的526656-2

2航空学报GPS驯服的原子钟同步信号PC软动态信号栅控信号聚信件发生器调制器栅控X焦SDD号控真空管道SDD电子处射线管光制栅控X射线管电源学理学器分子泵机械干泵PC制冷、电压控制X射线脉冲星信号模拟真空平台探测及数据处理系统图1X射线脉冲星导航地面实验系统组成Fig.1CompositionofX-raypulsarnavigationgroundexperimentsystem图2“天枢Ⅱ号”X射线脉冲星导航动态模拟系统Fig.2“TIANSHUⅡ”XPNAVgroundexperimentsystem配置，包含航天器及星体的轨道信息和脉冲星辐射特征等信息，并根据设定的条件实时计算脉冲轮廓数据，用于脉冲X射线发射源的控制，实现导航的静态及动态实验仿真。脉冲X射线发射源的核心器件为栅控X射线管，由灯丝、阴极罩、栅极、聚焦极、阳极靶、玻璃外壳组成，如图3所示。当灯丝加上电流，阳极加上高压，灯丝发射的电子就会在高压电场下加速运动并轰击阳极靶材产生X射线。如果栅极加Vg−栅极电压；Va−阳极电压上调制信号，轰击阳极靶的电子数目就随调制信图3　栅控X射线管结构图号的强度的变化而变化，发射X射线的强度也随Fig.3StructurediagramofgridcontrolledX-raytube之变化。聚焦极的主要作用是改变球管内部电场分布实现对从阴极发射的电子汇聚，缩小电子定的实验条件实时计算脉冲轮廓幅度数据，然后束斑尺寸，从而降低电子束的弥散时间，提高频经过数字信号传输接口发送给数据转换电子学，率响应能力。数据转换电子学将轮廓幅度信号转换为栅极电在模拟实验中，脉冲星模拟控制终端根据设压信号，在栅极电压的作用下栅控X射线管发射526656-3

3航空学报与模拟脉冲星轮廓一致的X射线脉冲信号。等时间间隔存储有限的数据点。固有时间序列τ0，τ1，…，τk是等间隔排列的，而其对应的SSB坐2X射线脉冲星导航动态模拟信号标时时间序列tSSB0，tSSB1，…，tSSBk并非等间隔的。在地面模拟系统的脉冲信号发生器中，信号采样2.1X射线脉冲星动态信号产生原理的频率是固定值，因而是等间隔的，因此需要通脉冲星的标准脉冲轮廓是将长时间观测的过插值的方法从等间隔的点轮廓数据点中恢复光子到达时间序列转换到太阳系质心坐标SSB出任意时刻的轮廓数据。（SolarSystemBarycenter）处的太阳系质心坐标航天器所在位置的TCB（BarycentricCoor⁃时，然后进行周期折叠得出的累积轮廓。在常规dinateTime）时间与SSB处TCB时间关系为的静态模拟实验中，脉冲星模拟控制器依照标准n·r12(t+■·r-r2■+脉冲轮廓模型计算X射线脉冲星轮廓数据。而SSB-t)=(n)c2cD0■■在动态实验中，由于航天器及星体的运动会产生2μsr+n·r-8各种动态效应，航天器记录的光子到达时间序列ln+1+O(10)（2）c3(b-n·b)周期折叠所得到的脉冲轮廓不再是标准脉冲轮式中：tSSB为SSB处的TCB时；t为航天器处TCB廓，动态信号发生器依据航天器的运动状态和标时；n为太阳系到脉冲星视线方向的单位向量；r准轮廓模型计算出动态脉冲轮廓数据，流程如为SSB指向航天器的距离矢量；D0为太阳质心到图4所示。脉冲星的距离；b为SSB到太阳质心的距离矢量；μs=GMs为太阳引力常数，G为引力常数，Ms为太阳质量。同样，对于选定的航天器轨道，r可以确定，其他量都为常量。对于航天器所在位置TCB时间tk，就能转化成SSB处TCB时间tSSBk。图4　动态脉冲轮廓数据计算过程将式（2）表达为tSSBk=g(tk)，式（1）表达为tk=Fig.4Calculationprocessofdynamicpulseprofiledatah(τk)，则由于脉冲星模拟控制器为数字系统，模拟脉-1-1τk=h[g(tSSBk)]（3）冲的相位并非连续的，因此只需根据计算周期转航天器接收的脉冲相位与所选坐标系无关，该相换等间隔的有限时间点。将航天器固有时τ按照位值对应着SSB处tSSBk时刻的相位值，即等时间间隔Δτ0排列为τ0，τ1，…，τk，与之对应的坐Øτ=Øt（4）SC(k)SSB(SSBk)标时为t0，t1，…，tk，有式中：ØSC(τk)为航天器处接收的脉冲相位；1tk-t0=k∆τ+2(rEAk·VEk+2rEAk·VEAk+ØSSB(tSSBk)为SSB处的脉冲相位。c1根据SSB处相位预测模型有2rEk·VEk)-2(rEA0·VE0+2rEA0·VEA0+cØSSB(tSSBk)=ØSSB(tSSB0)+f(tSSBk-tSSB0)+2rE0·VE0)（1）f2f3(tSSBk-tSSB0)+(tSSBk-tSSB0)（5）式中：k∆τ=τk-τ0；k为固有时τ的等分数；rEAk为26地心到航天器的距离矢量；rE为SSB到地心的距式中：f为脉冲星频率；f、f分别为其一阶导数和二离矢量；VE为地球相对于SSB的速度矢量；VEAk阶导数，本文模拟的Crab脉冲星和PSR为航天器相对地心的速度矢量。B1937+21周期稳定性较好，导数项可以略去，根据式（1）就可将等间隔的固有时间序列转脉冲周期为P，有换为坐标系时间序列，坐标系时间序列并非等间tSSBk-tSSB0ØSSB(tSSBk)=ØSSB(tSSB0)+（6）隔。由于脉冲轮廓没有具体函数表达式，通常是P526656-4

4航空学报极高压设置为15kV。X射线脉冲信号经过真空2.2　生成脉冲星导航动态模拟轮廓数据［19］管道后被聚焦型探测系统接收，经过信号处理等间隔固有时间序列τ0，τ1，…，τk对应的SSB后提取光子到达时间信息和能量信息，并存储在时间序列tSSB0，tSSB1，…，tSSBk并不是等间隔的，因信号处理计算机中。为了消除脉冲星模拟发射源此ØSC(τ1)，ØSC(τ2)，…，ØSC(τk)也不是等间隔的，和探测系统的电子学之间的时钟误差造成的对导不与标准轮廓的已知相位点重合。为得到这些相航模拟精度的影响，两部分电子学均采用经过位点对应的脉冲强度，需要对数据序列进行插值。GPS驯服的铷原子钟作为基准时钟，其长时间稳-12由于模拟源信号发生器的采样间隔为10μs，远小定性和短时间稳定性都可高达10s/s。同时为于毫秒脉冲星的周期，因而采用辛格函数SINC了降低模拟源和探测器的同步误差，在每组实验插值方法可以获得较为平滑的数据序列。起始时刻栅控信号调制器都会向SDD信号处理设标准轮廓已知相位点的时间间隔为T0，时电子学发送同步触发信号，使时间记录清零。间点tSSBk可表示为tSSBk=mT0+δk，m为整数，δk实验系统对静态脉冲信号的模拟与真实情表示与相邻已知相位点时间偏移量，则tSSBk对应况的符合程度主要由轮廓模拟相似度衡量。若脉冲强度为实验中模拟的脉冲星为固定周期的，则对采集到S(tSSBk)=S(mT0+δk)=的光子到达时间序列以真实周期进行周期折叠+∞得到累积轮廓，比较累积轮廓与真实轮廓即标准∑S(nT0)·sinc(mT0+δk-nT0)（7）轮廓的相似度及脉冲TOA差；若模拟的脉冲星n=-∞为精确计算tSSBk处强度值，实际上不可能取遍无周期是变化的，则通过对采集到的光子到达时间限个点。由于sinc函数在原点具有最大值，而随序列进行周期搜索比较测试周期值与理论周期［20-22］离开原点距离的增加其函数值衰减很快，所以值的误差。式（7）中只需取n=m附近较少的点就可得到比分别模拟了在SSB处接收PSRB0531+21、较精确的脉冲强度值。PSRB1937+21两颗脉冲星的X射线脉冲信号标由式（4）有S(τk)=S(tSSBk)，于是时间序列准轮廓数据源自欧洲EPN（TheEuropeanPulsar［23］Network）脉冲星数据库，不计脉冲星相位预测τ0，τ1，…，τk对应的脉冲强度就可以确定了，它对模型中频率高阶导数项，即脉冲周期固定不变，分应的是航天器在轨飞行时所接收到的脉冲星信别为33730、1560μs。每组实验时长45min，探测号轮廓。器接收的光子计数率控制在300count/s。动态信号脉冲发生器采用模块化设计，主要对采集的光子到达时间序列进行周期折叠，包括顶层控制、指令通信、高速数据传输、航天器设采集时间为N个周期，每个周期分为Nb个bin轨道积分、地球星历查询、相位计算、轮廓插值等区间，每个bin区间的大小为10μs，得出的累积脉模块。控制模块对其他各功能模块的运行状态冲轮廓如图5所示。累积脉冲轮廓与标准脉冲轮进行管理，各功能模块之间通过控制模块实现数据的交换，独立并行运行。3X射线脉冲星导航地面模拟实验3.1　脉冲星导航静态模拟实验在X射线脉冲星导航地面静态模拟实验中，脉冲信号发生器根据脉冲星轮廓数据及相位预测模型计算脉冲轮廓数据，脉冲轮廓数据经过栅控信号调制器后作用到栅控X射线管，栅控X射线管发射所需模拟的X射线脉冲，栅控X射线管阳526656-5

5航空学报3.2.1X射线脉冲星脉冲周期变化及轮廓还原在近地圆轨道中航天器做匀速圆周运动，接收的X射线脉冲周期是变化的。模拟了航天器在近地轨道飞行一周的情形，初始状态假定为T［x0，y0，z0，vx0，vy0，vz0］=［514306.42m，-1380476.65m，6723879.58m，T-6342.58m/s，3984.04m/s，1298.94m/s］，6即航天器轨道半径6.88337×10m，速度7602m/s，轨道周期5700s。初始时间为图5XPVAV地面静态模拟实验结果TMJD57723，对应地球位置状态为［Ex0，Ey0，Ez0］=Fig.5XPNAVgroundsimulationresultsinstatic［53400782398m，126870439795m，experimentT54972975401m］。分别对航天器在轨飞行时探廓之间的相似度用Pearson相关系数表示，两颗测器接收PSRB0531+21、PSRB1937+21脉冲脉冲星累积脉冲轮廓与标准脉冲轮廓之间的相星的脉冲轮廓情形进行了模拟，每组实验时长似度分别为99.5%、99.1%，均高于99%。1.5h。对每组实验采集到的光子到达时间序列进3.2　脉冲星导航动态模拟实验行周期搜索，得到的脉冲周期变化趋势如图6所在动态模拟实验中，动态信号发生器根据轨示。PSRB0531+21的周期搜索每150s搜索道力学模型可以由一个初始轨道状态值推导出后续一定有效时间内的轨道状态，因此只需提供初始时刻的航天器位置和速度坐标。实验开始后，动态信号发生器实时计算包含各种动态效应的变化周期的动态脉冲轮廓幅值数据，再通过栅控信号调制器控制栅控X射线管输出变化周期的动态X射线脉冲信号。探测系统接收到X射线光子信息后存储在信号处理计算机中，通过对采集到的光子到达时间序列进行周期搜索，并消除动态效应后进行脉冲轮廓还原，然后与理论值进行对比从而对实验的质量和模拟相似性进行验证。实验系统对动态脉冲信号的模拟与真实情况的符合程度主要由脉冲周期测试值与理论值误差、动态效应消除后累积轮廓相似度、脉冲TOA估计值与理论值的偏差衡量。对采集到的光子到达时间序列进行周期搜索，比较周期测试值与理论值的误差；此外，对采集到的光子到达时间序列，先将其转换到SSB处的TCB时间序列，然后根据真实周期进行周期折叠得到累积轮廓，比较累积轮廓与真实轮廓即标准轮廓的相似图6　动态模拟实验周期变化趋势测试度和脉冲TOA差。Fig.6Pulseperiodchangesindynamicexperiment526656-6

6航空学报10一次周期，bin的长度为100μs，周期测试值与理论索，测试得出的周期值分别为3.373×10、值偏差的标准差为38451ps。PSRB1937+21每91.560×10ps，与理论值偏差可以忽略。将新的30s搜索一次周期，bin的长度为10μs，周期测试时间序列进行周期折叠，得出的累积脉冲轮廓如值与理论值偏差的标准差为350ps。分析两颗脉图8所示，累积轮廓与标准轮廓的相似度分别为冲星地面模拟的周期测试值与理论值得偏差是99.86%、99.99%。由于信号模拟的栅控信号调制器的产生精度、X射线管长时间工作出射光子的强度稳定度等综合因素引起。在所有引起周期变化的动态效应中，由航天器的运动产生的多普勒效应占主导因素，因此，两组数据的周期变化趋势近似为正弦曲线，表征了航天器绕地球飞行一周时速度在脉冲星方向向量上投影的变化趋势。由于地球的运动同样会产生多普勒效应，会使中心周期（最大与最小周期的平均值）相对于静态周期有一定的偏移，两组数据的偏移量分别为844、-41ns。对每组实验数据，分别对观测时间为0~300、0~1800、0~3600、0~5400s的数据进行周期搜索，然后使用搜索到的周期进行周期折叠得到的累积轮廓如图7所示。随着观测时间的增加，脉冲展宽越明显，轮廓畸变越严重。图8　动态实验推算到SSB处的还原轮廓Fig.8PulseprofileattheSSBindynamicexperiment3.2.2　椭圆轨道超实时仿真实验超实时仿真的主要目的是在保证仿真信息图7　动态模拟实验不同时间段轮廓还原量不变的情况下缩短仿真时间，以满足长时间实Fig.7Pulseprofileatdifferenttimeperiodindynamic验的需要。在本次实验中，模拟航天器在椭圆轨experiment道飞行一周，假定轨道周期P为6h，根据开普勒第三定律，轨道半长轴a的表达式为探测系统采集到的光子到达时间为航天器32固有时间，根据实验设定的轨道信息及脉冲星相Pa=μ()（8）2π关参数，先将固有时间转换为TCB时间，再将其532转换到SSB处的TCB时间，即可消除动态效应式中：µ为开普勒常数，其值为3.9861×10km/s。7的影响。时间转换的步长为10μs，即每10μs计则半长轴a为1.6763523×10m，为了使实验现算一次轨道状态，并将该时间间隔内包含的光子象明显，将轨道偏心率设定较大为0.6，则半短轴7事件进行时间补偿，得到其对应的SSB处的TCBb为1.3314309×10m，近地点为6.578000×6时间值。对新的光子到达时间序列进行周期搜10m。椭圆轨道航天器速度表达式为526656-7

7航空学报21v=μ(-)（9）ra式中：r为航天器到地球的距离。根据式（9）计算得近地点速度为9869.98m/s，若将仿真起始点设置为近地点，则轨道初始状态设置为［x0，y0，Tz0，vx0，vy0，vz0］=［0m，6578000m，0m，T9869.98m/s，0m/s，0m/s］。初始时间设定为MJD57723历元，对应地球位置状态为［Ex0，TEy0，Ez0］=［53400782398m，126870439795m，T图10　椭圆轨道消除动态后PSRB1937+21还原轮廓54972975401m］。采用双倍速率的超实时仿Fig.10PSRB1937+21pulseprofileofellipticalorbit真，即实际的实验时长为3h。脉冲星选取PSRwithdynamiceffecteliminatedB1937+21，周期1.56ms，探测器光子计数率为100counts/s。可以有效对航天器在椭圆轨道飞行的情形进行在超实时仿真模式下，探测系统时标会压缩模拟，超实时仿真对脉冲周期变化的模拟及轮廓到与脉冲X射线源相同的尺度。对采集到的光模拟相似度与正常速率仿真模式处于同一水平，子到达时间序列进行周期搜索，bin的长度设置为仍能保持较高的实验质量。10μs，每120s搜索一次，得出的脉冲周期变化趋3.2.3　轨道机动模拟势如图9所示。周期测试值与理论值的偏差的标准差为592ps，中心周期与静态脉冲周期的偏移在航天器的飞行任务中，轨道调整与机动是量主要由地球公转产生的多普勒效应引起，偏移经常执行的环节。实验模拟了霍曼变轨的过程，量为-40.4ns，与图6（b）一致。初始轨道为高度是200km的圆轨道，经过转移轨道将初始轨道高度增加1200km，到达高度为1400km的目标轨道，轨道机动示意图如图11所示。起始时刻，航天器在初始轨道的A点处，飞行半周后到达转移轨道的近地点B点处，经过轨道机动后沿着转移轨道飞行到远地点C点，再经过轨道机动后进入目标轨道，最后在目标轨道上飞行半周后到达D点。轨道机动模拟主要通过设置不同轨道状态的轨道根数，从而改变轨道积分模型的输入状态。A点处航天器状态为［x0，Ty0，z0，vx0，vy0，vz0］=［0m，-6578000m，0m，T-7784.43m/s，0m/s，0m/s］，机动前B点状态图9　航天器为椭圆轨道条件下PSRB1937+21脉冲周为［0m，6578000m，0m，7784.43m/s，0m/s，期变化趋势T0m/s］，机动后B点的状态为［0m，6578000m，Fig.9PSRB1937+21pulseperiodchangeswithellip⁃Tticalorbit0m，8103.25m/s，0m/s，0m/s］，其中椭圆轨6道半长轴为7.178000×10m，半短轴为根据设定的轨道信息及脉冲星初始参数，将7.152879×106m，机动前C点状态为［0m，采集光子到达时间序列转换到SSB处的时间序-7778000m，0m，-6853.07m/s，0m/s，列，转换步长为10μs，得到的累积脉冲轮廓如T0m/s］，机动后C点的状态为［0m，-7778000m，图10所示，bin的大小为10μs，与标准轮廓相似T0m，-7158.80m/s，0m/s，0m/s］。脉冲星度为99.89%。由实验结果可知，地面实验系统选取PSRB1937+21，探测器光子计数率为526656-8

8航空学报效应模拟，其中，对动态效应中的短时间缓慢变化过程及长时间缓慢变化过程都能实现高相似度的模拟。因此，基于该地面实验系统，可以实现接近太空实验场景的模拟实验，开展更丰富的动态脉冲数据处理方法的研究以及系统性的导航滤波算法验证研究。基于该地面实验系统，可以为动态信号处理中各种影响因素的研究提供有力的实验数据支撑。4　结论1）提出了基于脉冲星辐射特性和航天器轨图11　霍曼变轨示意图Fig.11Hohmannorbitchangediagram道模型的X射线脉冲星动态信号模拟方法，设计了能实现多种实验场景模拟的X射线脉冲星模100counts/s，A点处起始时刻为MJD57723历拟源，并根据脉冲星辐射特性构建了X射线脉冲元，实验时长2.53h。星导航地面实验系统。对采集到的光子到达时间序列进行周期搜2）针对PSRB0531+21、PSRB1937+21索，bin的长度设置为10μs，每60s搜索一次，得出两颗脉冲星，开展了静态模拟实验，获得的脉冲的脉冲周期变化趋势如图12所示，其中，周期测轮廓相似度分别为99.5%、0.99.1%。试值与理论值的偏差的标准差为637ps。变化曲3）开展了200km轨道高度的动态近地圆轨9线中心周期为1.559959232×10ps，曲线与中心模拟实验，周期测试值和理论值的偏差分别为周期的交点分别在第0、2655、5681、9094s，初始38451、350ps，还原到SSB处的轮廓相似度分别轨道、转移轨道、目标轨道的半周期分别为为99.86%、99.99%。2655、3026、3413s，因此，中心周期的取值分别4）实现了椭圆轨道的超实时仿真实验，仿在轨道的A、B、C、D点，与轨道模型一致。真时长可压缩50%，轮廓相似度为99.89%。5）实现了基于霍曼变轨模型的轨道机动模拟，周期变化的测试值与理论值的偏差标准差为637ps。基于“天枢Ⅱ号”X射线脉冲星导航动态模拟系统，可以开展各种不同类型的脉冲星导航模拟实验，实现对探测器性能评估、脉冲信号处理方法及导航算法的验证，从而实现对影响导航性能因素的综合性分析，以对X射线脉冲星导航这样一项复杂的工程，提供有力的理论及实践支撑。参考文献图12　霍曼变轨过程中PSRB1937+21脉冲周期变化［1］SHEIKHSI.TheuseofvariablecelestialX-raysources趋势forspacecraftnavigation［D］.Maryland：UniversityofFig.12PulseperiodvariationofPSRB1937+21inMaryland，2005.theHohmannorbittransfer［2］HANSONJE.PrinciplesofX-raynavigation［D］.Stan⁃ford：StanfordUniversity，1996.根据以上动态模拟实验结果可知，“天枢Ⅱ［3］帅平，李明，陈绍龙，等.X射线脉冲星导航系统原理与方号”地面实验系统可以实现对不同类型的航天器法［M］.北京：中国宇航出版社，2009：11-22.轨道模拟以及航天器在轨运动产生的各种动态SHUAIP，LIM，CHENSL，etal.X-raypulsarnaviga⁃526656-9

9航空学报tionsystemprincipleandmethods［M］.Beijing：ChinaAs⁃SUNSM.StudyonautonomousnavigationmethodoftronauticPublishingHouse，2009：11-22（inChinese）.spacecraftbasedonX-raypulsars［D］.Changsha：Na⁃［4］GENDREAUKC，ARZOUMANIANZ，OKAJIMAtionalUniversityofDefenseTechnology，2011：116-124T.TheNeutronstarInteriorCompositionExploreR（inChinese）.（NICER）：anExplorermissionofopportunityforsoft［13］HUHJ，SONGJ，SHAOSP，etal.AnewmethodofX-raytimingspectroscopy［C］∥SpaceTelescopesandIn⁃evaluationofX-raypulsardetectorsensitivity［C］∥Fifthstrumentation2012：UltraviolettoGammaRay.Belling⁃SymposiumonNovelOptoelectronicDetectionTechnol⁃ham：SPIE，2012：322-329.ogyandApplication.Bellingham：SPIE，2019，11023：［5］WOODKS，DETERMANJR，RAYPS，etal.Using452-460.theunconventionalstellaraspect（USA）experimenton［14］方海燕，丛少鹏，孙海峰，等.具有多物理特性的X射线ARGOStodetermineatmosphericparametersbyX-ray脉冲星导航地面验证系统［J］.物理学报，2019，68（8）：occultation［C］∥SPIE，OpticalSpectroscopicTech⁃089701.niques，RemoteSensing，andInstrumentationforAtmo⁃FANGHY，CONGSP，SUNHF，etal.Groundveri⁃sphericandSpaceResearchIV.Bellingham：SPIE，ficationsystemofX-raypulsarnavigationwithmulti-2002，4485：258-265.physicalproperties［J］.ActaPhysicaSinica，2019，68［6］郑伟，王奕迪，汤国建.X射线脉冲星导航理论与应用（8）：089701（inChinese）.［M］.北京：科学出版社，2015：20-26.［15］苏哲，许录平，王婷.X射线脉冲星导航半物理仿真实验ZHENGW，WANGYD，TANGGJ.X-raypulsar-系统研究［J］.物理学报，2011，60（11）：119701.basednavigation：theoryandapplications［M］.Beijing：SUZ，XULP，WANGT.X-raypulsar-basednavigationSciencePress，2015：20-26（inChinese）.semi-physicalsimulationexperimentsystem［J］.Acta［7］WINTERNITZLM，MITCHELLJW，HASSOUNEHPhysicaSinica，2011，60（11）：119701（inChinese）.MA，etal.SEXTANTX-raypulsarnavigationdemonstra⁃［16］盛立志，赵宝升，吴建军，等.X射线脉冲星导航系统模tion：Additionalon-orbitresults［C］∥2018SpaceOpsCon⁃拟光源的研究［J］.物理学报，2013，62（12）：574-579.ference.Reston：AIAA，2018：2538.SHENGLZ，ZHAOBS，WUJJ，etal.Researchof［8］黄良伟，帅平，张新源，等.脉冲星导航试验卫星时间数X-raypulsarnavigationsimulationsource［J］.Acta据分析与脉冲轮廓恢复［J］.中国空间科学技术，2017，37PhysicaSinica，2013，62（12）：574-579（inChinese）.（3）：1-10.［17］徐能，盛立志，张大鹏，等.X射线脉冲星导航动态模拟HUANGLW，SHUAIP，ZHANGXY，etal.实验系统研制与性能测试［J］.物理学报，2017，66（5）：XPNAV-1Satellitetimingdataanalysisandpulseprofile334-340.recovery［J］.ChineseSpaceScienceandTechnology，XUN，SHENGLZ，ZHANGDP，etal.Development2017，37（3）：1-10（inChinese）.andperformancetestofdynamicsimulationsystemfor［9］郑世界，葛明玉，韩大炜，等.基于天宫二号POLAR的X-raypulsarnavigation［J］.ActaPhysicaSinica，2017，脉冲星导航实验［J］.中国科学：物理学力学天文学，66（5）：334-340（inChinese）.2017，47（9）：099505.［18］徐能.用于导航的X射线脉冲星地面模拟关键技术［D］.ZHENGSJ，GEMY，HANDW，etal.Testofpulsar西安：西安交通大学，2019：74-96.navigationwithPOLARonTG-2spacestation［J］.Sci⁃XUN.CoreTechnologiesofX-rayPulsarGroundSimu⁃entiaSinica（Physica，Mechanica＆Astronomica），lationforNavigationApplication［D］.Xi’an：Xi’anJiao⁃2017，47（9）：099505（inChinese）.tongUniversity，2019：74-96（inChinese）.［10］ZHENGSJ，ZHANGSN，LUFJ，etal.In-orbitdem⁃［19］李连升，梅志武，吕政欣，等.掠入射聚焦型X射线脉冲onstrationofX-raypulsarnavigationwiththeInsight-星望远镜及在轨数据分析［J］.兵器装备工程学报，2017，HXMTsatellite［J］.TheAstrophysicalJournalSupple⁃38（12）：175-179.mentSeries，2019，244（1），doi：10.3847/1538-4365/LILS，MEIZW，LYUZX，etal.Grazingincidencefo⁃ab371.cusingX-raypulsartelescopeandanalysisofIn-orbitob⁃［11］WINTERNITZLM，HASSOUNEHMA，servationdata［J］.JournalofOrdnanceEquipmentEngi⁃MITCHELLJW.X-raypulsarnavigationalgorithmsneering，2017，38（12）：175-179（inChinese）.andtestbedforSEXTANT［C］∥2015IEEEAerospace［20］周庆勇，姬剑锋，任红飞.非等间隔计时数据的X射线脉Conference.Piscataway：IEEEPress，2015：1-14.冲星周期快速搜索算法［J］.物理学报，2013，62（1）：［12］孙守明.基于X射线脉冲星的航天器自主导航方法研究019701.［D］.长沙：国防科技大学，2011：116-124.ZHOUQY，JIJF，RENHF.Quicksearchalgorithm526656-10

10航空学报ofX-raypulsarperiodbasedonunevenlyspacedtimingchanica＆Astronomica），2011，41（8）：1015-1020（indata［J］.ActaPhysicaSinica，2013，62（1）：019701（inChinese）.Chinese）.［22］WANGYD，ZHENGW.PulsarphaseandDopplerfre⁃［21］胡慧君，赵宝升，盛立志，等.一种基于泊松分布的提高quencyestimationforXNAVusingon-orbitepochfoldingX射线脉冲星脉冲轮廓信噪比的方法［J］.中国科学：物理［J］.IEEETransactionsonAerospaceandElectronic学力学天文学，2011，41（8）：1015-1020.Systems，2016，52（5）：2210-2219.HUHJ，ZHAOBS，SHENGLZ，etal.Amethodfor［23］EPNpulsardataarchive［EB/OL］.［2022-01-31］.improvingtheSNRofX-raypulsarpulseprofilebasedonhttp：∥www.jb.man.ac.Uk/research/pulsar/resources/Poissondistribution［J］.ScientiaSinica（Physica，Me⁃epn/browser.Html.（责任编辑：李丹）GroundtestbenchforX-raypulsarnavigationdynamicsimulationSHENGLizhi1，*，ZHENGWei2，SUTong1，ZHANGDapeng2，WANGYidi2，YANGXianghui1，12XUNeng，LIZhize1.Xi’anInstituteofOpticsandPrecisionMechanics，ChineseAcademyofSciences，Xi’an710119，China2.CollegeofAerospaceScienceandEngineering，NationalUniversityofDefenseTechnology，Changsha410073，ChinaAbstract：TheoreticalandpreliminaryexperimentalstudiesofX-raypulsar-basednavigationhavebeencarriedoutathomeandabroad，buttherearestillsomedefectssuchaslackofcredibleexperimentalverificationandimperfecttheo⁃reticalmodel.Inthispaper，amethodforsimulationofX-raypulsardynamicsignalsisproposed，whichcanbeusedtogeneratethepulsardynamicsignalbasedonX-raypulsarpropertiesandthespacecraftorbitmodel.AnX-raypulsarsimulationsourceandthewholegroundtestbencharedeveloped.Staticanddynamicgroundsimulationtestsareper⁃formed.BasedontheparametersofthepulsarsPSRB0531+21andPSRB1937+21，thepulseprofilesimilarityob⁃tainedwiththestaticsimulationtestis99.5%and99.1%，respectively.Dynamicsimulationtestsofthetwopulsarsatthecircularorbitheightof200kmareperformed.Deviationofthetestresultsandthetheoryresultsofthepulseperiodis38451psand350psandthepulseprofilesimilarityis99.8%and99.9%whenthetimingcoordinatetransfertoSSB，forPSRB0531+21andPSRB1937+21，respectively.Thetestbenchsystemcanrealizeorbitmaneuversimu⁃lationbasedontheHofmannmodel.Thegroundexperimentsystemhasstableperformanceandcanmeettheneedsofdifferenttypesofsimulationexperiments.Keywords：X-raypulsar-basednavigation；dynamicexperiment；X-raysimulationsource；X-raydetector；hardware-in-the-loopsimulationReceived：2021-11-15；Revised：2021-12-10；Accepted：2022-01-27；Publishedonline：2022-03-1909：27URL：https：//hkxb.buaa.edu.cn/CN/Y2023/V44/I3/526656Foundationitems：NationalNaturalScienceFoundationofChina（61901470）；TheCAS“LightofWestChina”Program（XAB2020YN13）；StateKeyLaboratoryFundofIntensePulsedRadiationSimulationandEffect（SKLIPR2021）*Correspondingauthor.E-mail：shenglizhi80@163.com526656-11