CN107144274A - 在轨x射线脉冲星计时模型构建方法 - Google Patents

在轨x射线脉冲星计时模型构建方法 Download PDF

Info

Publication number
CN107144274A
CN107144274A CN201710501797.7A CN201710501797A CN107144274A CN 107144274 A CN107144274 A CN 107144274A CN 201710501797 A CN201710501797 A CN 201710501797A CN 107144274 A CN107144274 A CN 107144274A
Authority
CN
China
Prior art keywords
mrow
msub
mtd
msup
mover
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710501797.7A
Other languages
English (en)
Other versions
CN107144274B (zh
Inventor
孙海峰
方海燕
孙雄
包为民
李小平
刘彦明
沈利荣
苏剑宇
李铁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201710501797.7A priority Critical patent/CN107144274B/zh
Publication of CN107144274A publication Critical patent/CN107144274A/zh
Application granted granted Critical
Publication of CN107144274B publication Critical patent/CN107144274B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/02Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Astronomy & Astrophysics (AREA)
  • Automation & Control Theory (AREA)
  • General Physics & Mathematics (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明提供一种在轨X射线脉冲星计时模型构建方法。所述在轨X射线脉冲星计时模型构建方法通过X射线脉冲星光子到达时间预处理、高精度的脉冲星自转频率计算、计时模型的初值区间估计、观测脉冲轮廓重构、高精度的脉冲到达时间估计和鲁棒的计时模型估计几种计算方法有机结合,有效解决了X射线信号强度弱、观测平台动态强所引入的数据问题,可提供X射线脉冲星导航所需的高精度脉冲星位置和时间等星历信息,实现精准导航。本发明适用于地球或其它行星的轨道卫星或星座的自主运控。

Description

在轨X射线脉冲星计时模型构建方法
技术领域
本发明涉及天文导航领域,尤其涉及一种在轨X射线脉冲星计时模型构建方法。
背景技术
X射线脉冲星导航是一种新型自主导航方式,可以在单一仪器上自主完成航天器位置、速度和时间等信息测量,能够为近地、深空以及行星际飞行航天器提供有效的导航与授时解决方案,被认为是未来深空自主导航领域最具潜力的导航与授时方式之一。
构建X射线计时模型是脉冲星导航系统的基础性工作。对于脉冲星自主导航而言,随着传输距离的增加以及星际遮挡等原因,无法依靠地面测控系统实时、连续地传输射电计时模型参数,此外,当导航脉冲星存在周期跃变等自转周期不规则现象时,需要利用X射线观测数据,重新评估原模型参数,并建立新的脉冲星计时模型。因此,利用X射线观测数据独立自主地建立计时模型对航天器自主导航具有明显的支撑作用。
构建X射线计时模型是获得脉冲星时的关键技术。脉冲星辐射的信号的重复周期均匀性好、频率间隔精度高、并且能连续稳定运行,符合公认的钟的原理与特征,在建立脉冲星时过程中,需要描述脉冲星的自转频率、脉冲到达时间与原子时的数学关系,因此,精确的脉冲星计时模型是分析和建立脉冲星时的关键环节。
构建X射线计时模型还是脉冲星辐射机理研究的重要补充手段。对于有些类型的脉冲星,比如PSR J0205+6449,先测到X射线接着才测到射电波段信号,随着X射线巡天的深入开展,未来会发现更多此类脉冲星,在未知其射电计时模型时,建立其X射线计时模型研究其辐射机理就显得至关重要。同时,建立X射线计时模型还可以与射电计时模型相互验证,为射电计时模型中存在的色散等误差源提供参考。
2016年10月8日,我国的“脉冲星试验卫星”成功发射,目前也已将构建X射线计时模型列为此次观测试验的一项科学和技术目标。
由于脉冲星在射电和X射线电磁波段均有辐射,所以脉冲星计时模型又可以分为射电计时模型和X射线计时模型,但与射电计时数据处理相比,X射线观测主要存在以下两个主要问题:
(1)X射线观测平台的动态性高。射电观测平台是地面的射电望远镜,而X射线观测平台则是星载X射线探测器。由于脉冲星观测可以认为是一种相对观测,地球自转周期24小时,而近地轨道卫星的自转周期仅100分钟,射电观测站的角加速度要远远低于卫星平台,两种平台的动态性会给信号处理带来较大差异,从而影响时间校正方法。
(2)X射线信号强度弱,射电信号通过大型天线也能转换成电信号,而X射线是随机分立的光子信息,这种弱信号的特性给X射线数据处理方法带来了很大的不同,包括到达时间校正方法、频率搜索方法、脉冲TOA估计方法、以及在初始相位不精确条件下,计时模型的代价函数设计等,这些方法均需要重新改正与设计。
现有的射电计时模型处理方法较为成熟,但由于X射线数据形式、平台动态、信号强度等与射电观测有很大的差异,尤其是以上两个问题的存在,需要改进和设计X射线计时模型构建方法。
发明内容
本发明的目的在于提供一种利用X射线观测数据计算脉冲星计时模型的方法框架。
本发明提供一种在轨X射线脉冲星计时模型构建方法,步骤一,X射线脉冲星光子到达时间预处理:将近地轨道航天器探测的光子到达时间的固有时转换到太阳系质心坐标系下的坐标时,对于航天器处测量的X射线光子到达时间假设用τobs表示,那么相应的SSB处的光子到达时间为
tSSB=τobs+△I+△C+△P+△R+△E+△S
公式中,△I为仪器响应延迟,用于修正X射线探测器的硬件响应时间延迟;
C是时钟校正,该项校正由于记录X射线到达时间的原子钟长期稳定度降低而引入的时间误差;
P是周年视差校正,校正地球绕太阳周年运动所产生的视差;
R是Roemer延迟,在真空中传播的同一光脉冲从航天器位置到SSB处的时间差;
E为太阳系Einstein延迟,该项是校正由于地球的运动,以及地球的引力势在卫星位置和SSB处的不同而导致的时延;
S Shapiro延迟,该项是校正光在经过太阳系内大质量天体时所引入的时间延迟;
步骤二,高精度的脉冲星自转频率计算:假设第i个光子到达时间分别用ti表示,对于假定的自转频率f,那么每个光子相位可以表示为
φi=2πmod(fti,1)
mod为取余运算符,那么一种高精度的方法的代价函数表示为
其中的分别为
最优的参数M的快速计算方法为
式中:为允许的积分平均误差,参数M可以通过给定的直接由该式计算;
步骤三,计时模型的初值区间估计:在建模X射线脉冲星计时模型时,需要计算并估计频率参数的初值区间;
步骤四,观测脉冲轮廓重构:在获得步骤二中的各数据包的自转频率后,利用历元折叠的方法获取观测脉冲轮廓;
步骤五,高精度的脉冲到达时间估计:根据观测脉冲轮廓模型估计初始相位的似然函数,计算脉冲峰值的达到时间。
与相关技术相比,本发明提供的在轨X射线脉冲星计时模型构建方法避免了传统的卡方量估计方法精度低带来的问题,同时解决了解决传统小波变换对脉冲轮廓消噪处理导致的信号急剧变化部分产生的震荡现象,同时保留脉冲轮廓的线性相位,提高计算精度。
附图说明
图1为本发明提供的在轨X射线脉冲星计时模型构建方法流程图;
图2本发明利用方法处理卫星实测数据得到的代价函数图;
图3为对7组卫星数据搜索的频率拟合的结果示意图;
图4为初始相位误差为0.2时,传统频率搜索代价函数得到的结果示意图;
图5为初始相位误差为0.2时,本发明提出的频率搜索代价函数得到的结果示意图。
具体实施方式
以下将参考附图并结合实施例来详细说明本发明。
如图1所示,所述在轨X射线脉冲星计时模型构建方法包括如下步骤:
1、X射线脉冲星光子到达时间预处理:将近地轨道航天器探测的光子到达时间的固有时转换到太阳系质心坐标系下的坐标时,对于航天器处测量的X射线光子到达时间假设用τobs表示,那么相应的SSB处的光子到达时间为
tSSB=τobs+△I+△C+△P+△R+△E+△S
公式中,△I为仪器响应延迟,用于修正X射线探测器的硬件响应时间延迟;
C是时钟校正,该项校正由于记录X射线到达时间的原子钟长期稳定度降低而引入的时间误差;
P是周年视差校正,校正地球绕太阳周年运动所产生的视差;
R是Roemer延迟,在真空中传播的同一光脉冲从航天器位置到SSB处的时间差;
E为太阳系Einstein延迟,该项是校正由于地球的运动,以及地球的引力势在卫星位置和SSB处的不同而导致的时延;
S Shapiro延迟,该项是校正光在经过太阳系内大质量天体时所引入的时间延迟。
该步骤与射电计时观测数据处理的差别:1)射电需要进行大气色散校正、星际色散校正,而X射线计时观测数据处理并不需要;2)射电仅需要对某个脉冲到达时间的固有时校正,而X射线需要对所有光子到达时间的固有时进行校正。
2、高精度的脉冲星自转频率计算:假设第i个光子到达时间分别用ti表示,对于假定的自转频率f,那么每个光子相位可以表示为
φi=2πmod(fti,1)
mod为取余运算符,那么一种高精度的方法的代价函数表示为
其中的分别为
最优的参数M的快速计算方法为
式中:为允许的积分平均误差,参数M可以通过给定的直接由该式计算。
利用方法得到的代价函数如图2所示,取代价函数最大值对应的频率即为最优频率,不过,为了提高频率估计的精度通常在最大频率附近采用插值拟合的方法,计算最优频率。
该步骤与现有频率计算方法的差别在于,这里提出一种的高精度频率估计方法,而且其中的参数m可以通过给定的允许的积分平均误差快速计算,避免了传统的卡方量估计方法精度低带来的问题。
3、计时模型的初值区间估计:在进行X射线脉冲星计时模型建模时,需要利用频率参数的初值区间。这里给出了一种基于最小二乘法的频率参数初值估计方法。
(1)将步骤二估计的m个数据包的脉冲星自转频率参数记为{(△t1,f1),(△t2,f2),…,(△tm,fm)},其中△ti=ti-t0,即观测时间与参考时间t0的时间间隔,那么测量数据可以表示为Y=[f1,f2,…,fm]T,要估计的参数向量表示为β=[f(L)/L!,f(L-1)/(L-1)!,…,f(0)/0!]T,根据最小二乘法,有以下关系成立
Y=Xβ+ε
其中的X用矩阵的形式展开为
(2)待估计的参数向量可以表示为
其中的R、Q分别是矩阵X进行QR分解的上、下三角矩阵;
(3)的最小二乘估计协方差矩阵为
(4)对于中的每个参数相应的参数估计区间为
其中,βi=f(i)/i!,tα/2(m-L-1)为学生氏t分布,置信水平为α,Cii为矩阵(XTX)-1第i行、第i列的元素,方差的无偏估计量
如图3所示,利用最小二乘法确定的频率范围(29.48,29.96),频率一阶导数范围(-5.543×10-10,1.279×10-9),频率二阶导数范围(-1.564×10-18,1.731×10-19)。经过该步骤,就可以为计时模型中的各频率参数提供初始估计区间。
4、观测脉冲轮廓重构:在获得步骤二中的各数据包的自转频率后,利用历元折叠的方法获取观测脉冲轮廓。
(1)选择各数据包的中点时刻tmid作为历元折叠的参考时间。
(2)计算光子到达时间相对于中点参考时间的相位,对于第i个数据包,其相应的光子到达相位为φj=fi(tj-tmid)。
(3)借助历元折叠方法计算观测脉冲轮廓
式中:N0为初始相位Φ(t0)引入的相位整数,Nb是历元折叠的bin块数量,NP为积分时间内包含的整周期,Tobs≈NP/f,Ci,k是第i个整相位第k个bin块内的光子数。
(4)借助具有平移不变性质的小波变换对观测脉冲轮廓进行消噪处理,此处使用平移不变小波消噪,以解决传统小波变换对脉冲轮廓消噪处理导致的信号急剧变化部分产生的震荡现象,同时保留脉冲轮廓的线性相位,以便于后续的脉冲到达时间估计处理。
5、高精度的脉冲到达时间估计:利用每个bin内的光子计数值的异方差泊松模型,估计初始相位φ0的似然函数为:
式中:λ(φ)为标准速率函数,观测脉冲轮廓的联合概率密度函数为:
其相对于φ0的最大值即可得到φ0的最大似然估计值的最大似然函数为:
其中:
脉冲到达时间差D可以表示为
脉冲峰值的到达时间为
tpk=tmid+D。
该步骤提出一种快速最大似然脉冲到达时间估计方法,计算脉冲峰值的到达时间,该步骤中使用快速最大似然方法,具更高的计算精度。
6、鲁棒的计时模型估计:该方法对初始参考频率参数敏感,当初始参考频率参数区间的误差较大时,由于光子到达时间相位取整过程发生变化,导致估计精度降低,甚至产生错误估计。这里提出了一种鲁棒的频率估计方法,频率搜索的代价函数设计如下:
(1)根据步骤三的计算结果,设置中各参数的搜索区间。
(2)设置中各参数的搜索步长,步长越精细,搜索的精度越高,但相应的搜索时间也会增加。
(3)计算在各假定频率参数下的频率搜索的代价函数χ2(β)。
(4)将多维问题转换为一维的量,绘制代价函数χ2(β)曲线
(5)计算代价函数χ2(β)曲线的最小值,计算相应的频率参数值
该方法克服了传统方法在频率区间误差较大时,计算的计时模型出错的问题,而且该方法还可以在初始相位位置的条件下进行计算,计算方法的鲁棒性得到大大加强。
如图4、图5所示,当初始相位误差为0.2时,利用传统的计时模型构建方法得到的代价函数是发散的,代价函数的整体误差较大,在中心频率处也没有有效地降低,利用传统代价函数得到的参数估计结果为:频率估计误差为﹣3.5×10-4,频率一阶导数估计误差为4.69×10-13,频率二阶导数估计误差为1.97×10-20。而利用本发明的方法,代价函数在中心频率处的代价函数值迅速减小,频率估计误差为﹣5.9067×10-9,频率一阶导数估计误差为9×10-15,频率二阶导数估计误差为9.7×10-21
从参数估计结果的比较来看,本发明提出的方法能够有效地降低初始相位误差对估计精度的影响,具有较好的鲁棒性。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其它相关的技术领域,均同理包括在本发明的专利保护范围内。

Claims (5)

1.一种在轨X射线脉冲星计时模型构建方法,其特征在于,包括:
步骤一,X射线脉冲星光子到达时间预处理:将近地轨道航天器探测的光子到达时间的固有时转换到太阳系质心坐标系下的坐标时,对于航天器处测量的X射线光子到达时间假设用τobs表示,那么相应的SSB处的光子到达时间为
tSSB=τobs+△I+△C+△P+△R+△E+△S
公式中,△I为仪器响应延迟,用于修正X射线探测器的硬件响应时间延迟;
C是时钟校正,该项校正由于记录X射线到达时间的原子钟长期稳定度降低而引入的时间误差;
P是周年视差校正,校正地球绕太阳周年运动所产生的视差;
R是Roemer延迟,在真空中传播的同一光脉冲从航天器位置到SSB处的时间差;
E为太阳系Einstein延迟,该项是校正由于地球的运动,以及地球的引力势在卫星位置和SSB处的不同而导致的时延;
S Shapiro延迟,该项是校正光在经过太阳系内大质量天体时所引入的时间延迟;
步骤二,高精度的脉冲星自转频率计算:假设第i个光子到达时间分别用ti表示,对于假定的自转频率f,那么每个光子相位可以表示为
φi=2πmod(fti,1)
mod为取余运算符,那么一种高精度的方法的代价函数表示为
<mrow> <msubsup> <mi>Z</mi> <mi>m</mi> <mn>2</mn> </msubsup> <mo>=</mo> <mn>2</mn> <mi>n</mi> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <mrow> <mo>(</mo> <msubsup> <mover> <mi>&amp;alpha;</mi> <mo>^</mo> </mover> <mi>k</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mover> <mi>&amp;beta;</mi> <mo>^</mo> </mover> <mi>k</mi> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> </mrow>
其中的分别为
<mrow> <msub> <mover> <mi>&amp;alpha;</mi> <mo>^</mo> </mover> <mi>k</mi> </msub> <mo>=</mo> <mn>1</mn> <mo>/</mo> <mi>n</mi> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <msub> <mi>k&amp;phi;</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mover> <mi>&amp;beta;</mi> <mo>^</mo> </mover> <mi>k</mi> </msub> <mo>=</mo> <mn>1</mn> <mo>/</mo> <mi>n</mi> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>k&amp;phi;</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow>
最优的参数M的快速计算方法为
<mrow> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <mi>&amp;pi;</mi> </mrow> </mfrac> <munderover> <mo>&amp;Integral;</mo> <mn>0</mn> <mrow> <mn>2</mn> <mi>&amp;pi;</mi> </mrow> </munderover> <msup> <mrow> <mo>(</mo> <mo>&amp;lsqb;</mo> <mn>1</mn> <mo>+</mo> <mn>2</mn> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msub> <mover> <mi>&amp;alpha;</mi> <mo>^</mo> </mover> <mi>k</mi> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mo>(</mo> <mrow> <mi>k</mi> <mi>&amp;phi;</mi> </mrow> <mo>)</mo> <mo>+</mo> <msub> <mi>&amp;beta;</mi> <mi>k</mi> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mo>(</mo> <mrow> <mi>k</mi> <mi>&amp;phi;</mi> </mrow> <mo>)</mo> <mo>&amp;rsqb;</mo> <mo>-</mo> <mi>h</mi> <mo>(</mo> <mi>&amp;phi;</mi> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mi>d</mi> <mi>&amp;phi;</mi> <mo>&amp;le;</mo> <msup> <mover> <mi>&amp;epsiv;</mi> <mo>&amp;OverBar;</mo> </mover> <mn>2</mn> </msup> </mrow>
式中:为允许的积分平均误差,参数M可以通过给定的直接由该式计算;
步骤三,计时模型的初值区间估计:在建模X射线脉冲星计时模型时,需要计算并估计频率参数的初值区间;
步骤四,观测脉冲轮廓重构:在获得步骤二中的各数据包的自转频率后,利用历元折叠的方法获取观测脉冲轮廓;
步骤五,高精度的脉冲到达时间估计:根据观测脉冲轮廓模型估计初始相位的似然函数,计算脉冲峰值的达到时间。
2.根据权利要求1所述的在轨X射线脉冲星计时模型构建方法,其特征在于,在步骤三中,计时模型的初值区间估计给出了一种基于最小二乘法的频率参数初值估计方法,具体为:
1)将步骤二估计的m个数据包的脉冲星自转频率参数记为{(△t1,f1),(△t2,f2),…,(△tm,fm)},其中△ti=ti-t0,即观测时间与参考时间t0的时间间隔,那么测量数据可以表示为Y=[f1,f2,…,fm]T,要估计的参数向量表示为β=[f(L)/L!,f(L-1)/(L-1)!,…,f(0)/0!]T,根据最小二乘法,有以下关系成立
Y=Xβ+ε
其中的X用矩阵的形式展开为
<mrow> <mi>X</mi> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msup> <msub> <mi>&amp;Delta;t</mi> <mn>1</mn> </msub> <mi>L</mi> </msup> </mrow> </mtd> <mtd> <mrow> <msup> <msub> <mi>&amp;Delta;t</mi> <mn>1</mn> </msub> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </msup> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <msub> <mi>&amp;Delta;t</mi> <mn>2</mn> </msub> <mi>L</mi> </msup> </mrow> </mtd> <mtd> <mrow> <msup> <msub> <mi>&amp;Delta;t</mi> <mn>2</mn> </msub> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </msup> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <msub> <mi>&amp;Delta;t</mi> <mi>m</mi> </msub> <mi>L</mi> </msup> </mrow> </mtd> <mtd> <mrow> <msup> <msub> <mi>&amp;Delta;t</mi> <mi>m</mi> </msub> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </msup> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>
2)待估计的参数向量可以表示为
<mrow> <mover> <mi>&amp;beta;</mi> <mo>^</mo> </mover> <mo>=</mo> <msup> <mi>R</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mi>Q</mi> <mi>T</mi> </msup> <mi>Y</mi> </mrow>
其中的R、Q分别是矩阵X进行QR分解的上、下三角矩阵;
3)的最小二乘估计协方差矩阵为
<mrow> <mi>C</mi> <mi>o</mi> <mi>v</mi> <mrow> <mo>(</mo> <mover> <mi>&amp;beta;</mi> <mo>^</mo> </mover> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mi>&amp;sigma;</mi> <mn>2</mn> </msup> <msup> <mrow> <mo>(</mo> <msup> <mi>X</mi> <mi>T</mi> </msup> <mi>X</mi> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> </mrow>
4)对于中的每个参数相应的参数估计区间为
<mrow> <mi>C</mi> <mi>I</mi> <mo>=</mo> <mo>&amp;lsqb;</mo> <msub> <mover> <mi>&amp;beta;</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>&amp;PlusMinus;</mo> <msub> <mi>t</mi> <mrow> <mi>&amp;alpha;</mi> <mo>/</mo> <mn>2</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <mi>L</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mover> <mi>&amp;sigma;</mi> <mo>^</mo> </mover> <msqrt> <msub> <mi>C</mi> <mrow> <mi>i</mi> <mi>i</mi> </mrow> </msub> </msqrt> <mo>&amp;rsqb;</mo> <mi>i</mi> <mo>!</mo> </mrow>
其中,βi=f(i)/i!,tα/2(m-L-1)为学生氏t分布,置信水平为α,Cii为矩阵(XTX)-1第i行、第i列的元素,方差的无偏估计量
3.根据权利要求2所述的在轨X射线脉冲星计时模型构建方法,其特征在于,获取观测脉冲轮廓的方法为:
1)选择各数据包的中点时刻tmid作为历元折叠的参考时间;
2)计算光子到达时间相对于中点参考时间的相位,对于第i个数据包,其相应的光子到达相位为φj=fi(tj-tmid);
3)借助历元折叠方法计算观测脉冲轮廓
<mrow> <mover> <mi>&amp;lambda;</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>N</mi> <mi>b</mi> </msub> <mi>f</mi> </mrow> <msub> <mi>N</mi> <mi>P</mi> </msub> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <msub> <mi>N</mi> <mn>0</mn> </msub> <mo>+</mo> <mn>1</mn> </mrow> <mrow> <msub> <mi>N</mi> <mn>0</mn> </msub> <mo>+</mo> <msub> <mi>N</mi> <mi>P</mi> </msub> </mrow> </munderover> <msub> <mi>C</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>,</mo> <mi>k</mi> <mo>&amp;Element;</mo> <mo>(</mo> <mn>0</mn> <mo>,</mo> <msub> <mi>N</mi> <mi>b</mi> </msub> <mo>&amp;rsqb;</mo> </mrow>
式中:N0为初始相位Φ(t0)引入的相位整数,Nb是历元折叠的bin块数量,NP为积分时间内包含的整周期,Tobs≈NP/f,Ci,k是第i个整相位第k个bin块内的光子数;
4)借助具有平移不变性质的小波变换对观测脉冲轮廓进行消噪处理,获得高信噪比观测脉冲轮廓,以进行后续时延估计处理。
4.根据权利要求3所述的在轨X射线脉冲星计时模型构建方法,其特征在于,利用每个bin内的光子计数值的异方差泊松模型,估计初始相位φ0的似然函数为:
<mrow> <mi>P</mi> <mrow> <mo>(</mo> <msub> <mi>c</mi> <mi>i</mi> </msub> <mo>|</mo> <msub> <mi>&amp;phi;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msup> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>(</mo> <msub> <mi>&amp;phi;</mi> <mi>i</mi> </msub> <mo>)</mo> <msub> <mi>T</mi> <mrow> <mi>o</mi> <mi>b</mi> <mi>s</mi> </mrow> </msub> <mo>/</mo> <msub> <mi>N</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>c</mi> <mi>i</mi> </msub> </msup> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&amp;lambda;</mi> <mo>(</mo> <msub> <mi>&amp;phi;</mi> <mi>i</mi> </msub> <mo>)</mo> <msub> <mi>T</mi> <mrow> <mi>o</mi> <mi>b</mi> <mi>s</mi> </mrow> </msub> <mo>/</mo> <msub> <mi>N</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>c</mi> <mi>i</mi> </msub> <mo>!</mo> </mrow> </mfrac> </mrow> 2
式中:λ(φ)为标准速率函数,观测脉冲轮廓的联合概率密度函数为:
<mrow> <mi>P</mi> <mrow> <mo>(</mo> <msubsup> <mrow> <mo>{</mo> <msub> <mi>c</mi> <mi>i</mi> </msub> <mo>}</mo> </mrow> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>b</mi> </msub> </msubsup> <mo>|</mo> <msub> <mi>&amp;phi;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Pi;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>b</mi> </msub> </munderover> <mi>P</mi> <mrow> <mo>(</mo> <msub> <mi>c</mi> <mi>i</mi> </msub> <mo>|</mo> <msub> <mi>&amp;phi;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Pi;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>b</mi> </msub> </munderover> <mfrac> <mrow> <msup> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>(</mo> <msub> <mi>&amp;phi;</mi> <mi>i</mi> </msub> <mo>)</mo> <msub> <mi>T</mi> <mrow> <mi>o</mi> <mi>b</mi> <mi>s</mi> </mrow> </msub> <mo>/</mo> <msub> <mi>N</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>c</mi> <mi>i</mi> </msub> </msup> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&amp;lambda;</mi> <mo>(</mo> <msub> <mi>&amp;phi;</mi> <mi>i</mi> </msub> <mo>)</mo> <msub> <mi>T</mi> <mrow> <mi>o</mi> <mi>b</mi> <mi>s</mi> </mrow> </msub> <mo>/</mo> <msub> <mi>N</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>c</mi> <mi>i</mi> </msub> <mo>!</mo> </mrow> </mfrac> </mrow>
其相对于φ0的最大值即可得到φ0的最大似然估计值的最大似然函数为:
<mrow> <mi>L</mi> <mi>L</mi> <mi>F</mi> <mrow> <mo>(</mo> <msub> <mi>&amp;phi;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>A</mi> <mo>-</mo> <mi>B</mi> <mo>+</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>b</mi> </msub> </munderover> <msub> <mi>c</mi> <mi>i</mi> </msub> <mi>l</mi> <mi>n</mi> <mo>&amp;lsqb;</mo> <mi>&amp;lambda;</mi> <mrow> <mo>(</mo> <msub> <mi>&amp;phi;</mi> <mn>0</mn> </msub> <mo>+</mo> <mi>i</mi> <mo>/</mo> <msub> <mi>N</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow>
其中:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>A</mi> <mo>=</mo> <mi>l</mi> <mi>n</mi> <munderover> <mo>&amp;Pi;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mn>5</mn> </msub> </munderover> <mfrac> <msup> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mrow> <mi>o</mi> <mi>b</mi> <mi>s</mi> </mrow> </msub> <mo>/</mo> <msub> <mi>N</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>c</mi> <mi>i</mi> </msub> </msup> <mrow> <msub> <mi>c</mi> <mi>i</mi> </msub> <mo>!</mo> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>B</mi> <mo>=</mo> <mfrac> <msub> <mi>T</mi> <mrow> <mi>o</mi> <mi>b</mi> <mi>s</mi> </mrow> </msub> <msub> <mi>N</mi> <mi>b</mi> </msub> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mn>5</mn> </msub> </munderover> <mi>&amp;lambda;</mi> <mrow> <mo>(</mo> <msub> <mi>&amp;phi;</mi> <mn>0</mn> </msub> <mo>+</mo> <mi>i</mi> <mo>/</mo> <msub> <mi>N</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
脉冲到达时间差D可以表示为:
<mrow> <mi>D</mi> <mo>=</mo> <msub> <mover> <mi>&amp;phi;</mi> <mo>^</mo> </mover> <mn>0</mn> </msub> <mi>T</mi> <mo>=</mo> <mi>arg</mi> <munder> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mrow> <msub> <mi>&amp;phi;</mi> <mn>0</mn> </msub> <mo>&amp;Element;</mo> <mo>&amp;lsqb;</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>)</mo> </mrow> </munder> <mi>T</mi> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>b</mi> </msub> </munderover> <msub> <mi>c</mi> <mi>i</mi> </msub> <mi>l</mi> <mi>n</mi> <mo>&amp;lsqb;</mo> <mi>&amp;lambda;</mi> <mrow> <mo>(</mo> <msub> <mi>&amp;phi;</mi> <mn>0</mn> </msub> <mo>+</mo> <mi>i</mi> <mo>/</mo> <msub> <mi>N</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow>
脉冲峰值的到达时间为
tpk=tmid+D。
5.根据权利要求4所述的在轨X射线脉冲星计时模型构建方法,其特征在于,还包括步骤六鲁棒的计时模型估计,该步骤为频率搜索的代价函数,公式如下:
<mrow> <msup> <mi>&amp;chi;</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mi>&amp;beta;</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mfrac> <msup> <mrow> <mo>(</mo> <msub> <mover> <mi>&amp;phi;</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>(</mo> <mrow> <mi>&amp;beta;</mi> <mo>,</mo> <msub> <mi>t</mi> <mi>i</mi> </msub> </mrow> <mo>)</mo> <mo>-</mo> <msub> <mover> <mi>&amp;phi;</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mo>(</mo> <mrow> <mi>&amp;beta;</mi> <mo>,</mo> <msub> <mi>t</mi> <mi>i</mi> </msub> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mi>N</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> <msubsup> <mi>&amp;sigma;</mi> <mi>i</mi> <mn>2</mn> </msubsup> </mrow> </mfrac> </mrow>
1)根据步骤三的计算结果,设置中各参数的搜索区间;
2)设置中各参数的搜索步长,步长越精细,搜索的精度越高,但相应的搜索时间也会增加;
3)计算在各假定频率参数下的频率搜索的代价函数χ2(β);
4)将多维问题转换为一维的量,绘制代价函数χ2(β)曲线;
5)计算代价函数χ2(β)曲线的最小值,计算相应的频率参数值。
CN201710501797.7A 2017-06-27 2017-06-27 在轨x射线脉冲星计时模型构建方法 Active CN107144274B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710501797.7A CN107144274B (zh) 2017-06-27 2017-06-27 在轨x射线脉冲星计时模型构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710501797.7A CN107144274B (zh) 2017-06-27 2017-06-27 在轨x射线脉冲星计时模型构建方法

Publications (2)

Publication Number Publication Date
CN107144274A true CN107144274A (zh) 2017-09-08
CN107144274B CN107144274B (zh) 2019-06-21

Family

ID=59785071

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710501797.7A Active CN107144274B (zh) 2017-06-27 2017-06-27 在轨x射线脉冲星计时模型构建方法

Country Status (1)

Country Link
CN (1) CN107144274B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110986922A (zh) * 2019-12-30 2020-04-10 西安电子科技大学 X射线脉冲星短时观测高信噪比轮廓获取方法
CN113219815A (zh) * 2021-05-06 2021-08-06 中国科学院国家授时中心 一种基于x射线脉冲星的深空授时方法
CN113375697A (zh) * 2021-06-23 2021-09-10 西安电子科技大学 基于最大后验估计的x射线脉冲星信号时延估计方法
CN114739414A (zh) * 2022-04-22 2022-07-12 西安电子科技大学 基于频率预报的航天器处脉冲星导航观测量获取方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101038169A (zh) * 2007-02-13 2007-09-19 北京空间飞行器总体设计部 基于x射线脉冲星的导航卫星自主导航系统与方法
CN103776454A (zh) * 2014-01-21 2014-05-07 西安电子科技大学 基于x射线脉冲星的最大似然相位估计方法
CN104296755A (zh) * 2014-10-23 2015-01-21 中国空间技术研究院 一种x射线脉冲星导航脉冲toa的确定方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101038169A (zh) * 2007-02-13 2007-09-19 北京空间飞行器总体设计部 基于x射线脉冲星的导航卫星自主导航系统与方法
CN103776454A (zh) * 2014-01-21 2014-05-07 西安电子科技大学 基于x射线脉冲星的最大似然相位估计方法
CN104296755A (zh) * 2014-10-23 2015-01-21 中国空间技术研究院 一种x射线脉冲星导航脉冲toa的确定方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
周庆勇: "脉冲星计时模型及应用", 《四川兵工学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110986922A (zh) * 2019-12-30 2020-04-10 西安电子科技大学 X射线脉冲星短时观测高信噪比轮廓获取方法
CN110986922B (zh) * 2019-12-30 2022-09-06 西安电子科技大学 X射线脉冲星短时观测高信噪比轮廓获取方法
CN113219815A (zh) * 2021-05-06 2021-08-06 中国科学院国家授时中心 一种基于x射线脉冲星的深空授时方法
CN113375697A (zh) * 2021-06-23 2021-09-10 西安电子科技大学 基于最大后验估计的x射线脉冲星信号时延估计方法
CN113375697B (zh) * 2021-06-23 2022-09-23 西安电子科技大学 基于最大后验估计的x射线脉冲星信号时延估计方法
CN114739414A (zh) * 2022-04-22 2022-07-12 西安电子科技大学 基于频率预报的航天器处脉冲星导航观测量获取方法
CN114739414B (zh) * 2022-04-22 2024-08-06 西安电子科技大学 基于频率预报的航天器处脉冲星导航观测量获取方法

Also Published As

Publication number Publication date
CN107144274B (zh) 2019-06-21

Similar Documents

Publication Publication Date Title
US7197381B2 (en) Navigational system and method utilizing sources of pulsed celestial radiation
CN107144274A (zh) 在轨x射线脉冲星计时模型构建方法
Liu et al. Novel algorithm for X-ray pulsar navigation against Doppler effects
Iess et al. Gravity, geodesy and fundamental physics with BepiColombo’s MORE investigation
CN102997922B (zh) 一种利用光学导航信息的脉冲到达时间差确定方法
CN105203101A (zh) 一种基于目标天体星历修正的深空探测器捕获段天文导航方法
US6859690B2 (en) Method for using GPS and crosslink signals to correct ionospheric errors in space navigation solutions
Tartaglia et al. A null frame for spacetime positioning by means of pulsating sources
Wang et al. Absolute navigation for Mars final approach using relative measurements of X-ray pulsars and Mars orbiter
Feng et al. Autonomous orbit determination and its error analysis for deep space using X-ray pulsar
CN104316048A (zh) 一种普适性的脉冲星自主导航测量模型构建方法
Crassidis et al. Efficient and optimal attitude determination using recursive global positioning system signal operations
Huang et al. Navigation using binary pulsars
Kim et al. Flight performance analysis of GRACE K-band ranging instrument with simulation data
Mengna et al. Frequency steering of spaceborne clocks based on XPNAV-1 observations
Syndergaard A new algorithm for retrieving GPS radio occultation total electron content
Butala et al. Dynamic tomographic imaging of the solar corona
Moeller Nanosatellites: The next big chapter in atmospheric tomography
Liang et al. New observation strategy for X‐ray pulsar navigation using a single detector
Sheikh et al. Deep space navigation augmentation using variable celestial x-ray sources
Wandelt Primordial Non-Gaussianity
Zhou et al. Autonomous navigation for lunar satellite using X‐ray pulsars with measurement faults
Chen et al. Spacecraft autonomous GPS navigation based on polytopic linear differential inclusion
Huang et al. On time transfer in X-ray pulsar navigation
LARESE DE PASQUA Space Navigation with Optical Pulsars

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant