CN103196451A - 一种基于Fisher信息矩阵的导航脉冲星选择方法 - Google Patents
一种基于Fisher信息矩阵的导航脉冲星选择方法 Download PDFInfo
- Publication number
- CN103196451A CN103196451A CN2013101171031A CN201310117103A CN103196451A CN 103196451 A CN103196451 A CN 103196451A CN 2013101171031 A CN2013101171031 A CN 2013101171031A CN 201310117103 A CN201310117103 A CN 201310117103A CN 103196451 A CN103196451 A CN 103196451A
- Authority
- CN
- China
- Prior art keywords
- pulsar
- centerdot
- navigation
- sigma
- information matrix
- 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
Links
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种基于Fisher信息矩阵的导航脉冲星选择方法,属于深空探测技术领域。本发明方法建立了脉冲星导航的测量简化模型,基于信息论相关概念构建了导航状态的Fisher信息矩阵,根据使得Fisher信息矩阵的行列式最大的原则选择导航脉冲星,提高导航状态可观测性;提高脉冲星导航的可观测性及导航精度。
Description
技术领域
本发明涉及一种基于Fisher信息矩阵的导航脉冲星选择方法,属于深空探测技术领域。
背景技术
脉冲星是一类旋转的中子星,可以发射独特、稳定、可预测的X射线信号。脉冲星的特性为空间自主导航提供了新的发展方向。通过对脉冲星发射的X射线的测量可以实现不依赖地面深空网的自主导航。近年来美国航天局及欧空局都在考虑将脉冲星作为潜在的导航信息源,国内外学者也针对脉冲星导航的基础理论及导航精度进行了大量的研究。
脉冲星作为导航信息源,其位置及空间分布必将影响到导航性能,但已有的研究主要侧重脉冲星导航可行性的验证,对如何合理地选择导航脉冲星缺少***的研究。所以需要结合信息论的思想,对导航脉冲星进行优化选择以提高导航***状态的可观测性。
发明内容
本发明的目的是为了提高脉冲星导航的可观测性及导航精度,提出一种基于Fisher信息矩阵的导航脉冲星选择方法,结合信息论中Fisher信息矩阵的概念,通过构建导航状态的Fisher信息矩阵并对其行列式进行分析,优化选择导航脉冲星使得导航方案可观测性最强。
基于Fisher信息矩阵的导航脉冲星选择方法流程如下:
步骤1:建立脉冲星导航测量模型
通过接收脉冲星发射的X射线,并与通过地面长期观测得到的标准波形进行比对,得到X射线到达探测器与到达太阳系质心SSB的时间差
i=1,2,…,m
式中c为光速,ni为太阳系质心SSB到第i颗脉冲星的单位矢量,b是SSB在日心惯性坐标系下的位置矢量,rb为探测器相对SSB的位置矢量,满足:
rs=b+rb (2)
式中rs=[rx,ry,rz]T为探测器位置状态矢量。另外D0i为第i颗脉冲星到日心的距离,m为所测量的脉冲星数量。
由于太阳质量占太阳系质量的99%以上,并且脉冲星通常距离太阳系非常遥远,将(1)式简化为:
式中εi为测量噪声,服从零均值、σ标准差的高斯分布。由于估计探测器状态需要同时利用至少三颗脉冲星的测量信息,所以取m≥3,并且得到脉冲星导航测量模型为:
y=[y1,y2,....,ym]T (4)
yi=Δti=hi(rs)+εi, i=1,2,…,m
步骤2:构建Fisher信息矩阵
由于εi为高斯白噪声,将所测量的m颗脉冲星的联合似然函数表示为
取(5)式的自然对数,并只保留与状态rs相关项,得到对数似然函数
从而Fisher信息矩阵F由下式计算:
步骤3:导航脉冲星选择
根据Cramér-Rao定理,Fisher信息矩阵奇异性越强导航状态可观测性越差。而Fisher信息矩阵的行列式可以反映其奇异性,所以应选择导航脉冲星使得Fisher信息矩阵的行列式最大。Fisher信息矩阵的行列式只与矩阵的行列式相关,且
有益效果
(1)建立了脉冲星导航的测量简化模型。
(2)基于信息论相关概念构建了导航状态的Fisher信息矩阵。
(3)通过使得Fisher信息矩阵的行列式最大选择导航脉冲星,提高导航状态可观测性。
附图说明
图1为本发明方法的流程图;
图2为具体实施方式中所选导航脉冲星的空间位置图;
图3为具体实施方式中导航脉冲星构建火星最终接近段自主导航误差曲线,其中图a)显示x轴位置误差;图b)显示y轴位置误差;图c)显示z轴位置误差。
具体实施方式
本实例针对火星最终接近段基于脉冲星测量的自主导航方案,构建导航状态的Fisher信息矩阵并对其行列式进行分析,优化选择导航脉冲星使得导航方案可观测性最强,并利用数学仿真验证了方法的可行性。本实例的具体实施方法如下:
步骤1:脉冲星导航测量模型建立
通过脉冲星发射的X射线的接收,并与通过地面长期观测得到的标准波形的比对,可以得到X射线到达探测器与到达太阳系质心SSB的时间差
i=1,2,…,m
式中c为光速,ni为太阳系质心SSB到第i颗脉冲星的单位矢量,b是SSB在日心惯性坐标系下的位置矢量,rb为探测器相对SSB的位置矢量,满足:
rs=b+rb (2)
式中rs=[rx,ry,rz]T为探测器位置状态矢量。另外D0i为第i颗脉冲星到日心的距离。这里考虑利用三颗脉冲星的测量信息,所以取m=3。
由于太阳质量占太阳系质量的99%以上,并且脉冲星通常距离太阳系非常遥远,所以(1)式可以简化为:
式中εi为测量噪声,服从零均值、σ标准差的高斯分布。进一步得到脉冲星导航测量模型为:
y=[y1,y2,y3]T (4)
yi=Δti=hi(rs)+εi, i=1,2,3
步骤2:Fisher信息矩阵构建
由于εi为高斯白噪声,将联合似然函数表示为
将(5)式取自然对数,只保留与状态rs相关项,得到对数似然函数
进一步Fisher信息矩阵F可由下式计算:
步骤3:导航脉冲星选择
根据Cramér-Rao定理Fisher信息矩阵奇异性越强导航状态可观测性越差。而Fisher信息矩阵的行列式可以反映其奇异性,所以应选择导航脉冲星使得Fisher信息矩阵的行列式最大。Fisher信息矩阵的行列式只与矩阵的行列式相关,且
所以从备选导航脉冲星中选择三颗导航脉冲星使得(n1·(n2×n3))2最大。实例中所涉及的备选导航脉冲星通如表1所示。
最终选择J0030+0451、B0633+17和B0540-69作为导航脉冲星,其空间位置如图2所示。
采用选择的导航脉冲星构建火星最终接近段自主导航方案,针对火星大气进入前12小时的自主导航性能进行仿真分析。进入前6小时地面深空网测量数据切断,随后探测器利用对所选三颗导航脉冲星测量信息进行自主导航,导航误差结果如图3所示。其中图a)、b)、c)分别显示x、y、z轴位置误差。
通过对选择的导航脉冲星进行测量,火星最终接近段自主导航***的收敛速度更快、精度高、可观测性强,证明基于Fisher信息矩阵的导航脉冲星选择方法的可行性。
表1备选导航脉冲星参数
Claims (2)
1.一种基于Fisher信息矩阵的导航脉冲星选择方法,其特征在于:具体包括如下步骤:
步骤1:建立脉冲星导航测量模型;
通过接收脉冲星发射的X射线,并与通过地面长期观测得到的标准波形进行比对,得到X射线到达探测器与到达太阳系质心SSB的时间差
i=1,2,…,m
式中c为光速,ni为太阳系质心SSB到第i颗脉冲星的单位矢量,b是SSB在日心惯性坐标系下的位置矢量,rb为探测器相对SSB的位置矢量,满足:
rs=b+rb (2)
式中rs=[rx,ry,rz]T为探测器位置状态矢量;D0i为第i颗脉冲星到日心的距离,m为所测量的脉冲星数量;
将(1)式简化为:
式中εi为测量噪声,服从零均值、σ标准差的高斯分布;得到脉冲星导航测量模型为:
y=[y1,y2,....,ym]T (4)
yi=Δti=hi(rs)+εi,i=1,2,…,m
步骤2:构建Fisher信息矩阵;
将所测量的m颗脉冲星的联合似然函数表示为
取(5)式的自然对数,并只保留与状态rs相关项,得到对数似然函数
Fisher信息矩阵F由下式计算:
步骤3:选择导航脉冲星;
2.根据权利要求1所述的一种基于Fisher信息矩阵的导航脉冲星选择方法,其特征在于:m≥3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310117103.1A CN103196451B (zh) | 2013-04-07 | 2013-04-07 | 一种基于Fisher信息矩阵的导航脉冲星选择方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310117103.1A CN103196451B (zh) | 2013-04-07 | 2013-04-07 | 一种基于Fisher信息矩阵的导航脉冲星选择方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103196451A true CN103196451A (zh) | 2013-07-10 |
CN103196451B CN103196451B (zh) | 2015-10-14 |
Family
ID=48719152
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310117103.1A Expired - Fee Related CN103196451B (zh) | 2013-04-07 | 2013-04-07 | 一种基于Fisher信息矩阵的导航脉冲星选择方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103196451B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103913173A (zh) * | 2014-04-18 | 2014-07-09 | 中国人民解放军国防科学技术大学 | 单x射线脉冲星导航观星序列选择方法 |
CN103954287A (zh) * | 2014-03-31 | 2014-07-30 | 北京控制工程研究所 | 深空探测自主导航的一种路标规划方法 |
CN106871909A (zh) * | 2017-02-20 | 2017-06-20 | 中国人民解放军国防科学技术大学 | 一种多航天器***下基于Fisher信息矩阵的脉冲星选星方法 |
CN115355915A (zh) * | 2022-10-19 | 2022-11-18 | 中国人民解放军国防科技大学 | 一种卫星方向矢量导航的观星序列规划方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005043189A (ja) * | 2003-07-28 | 2005-02-17 | National Institute Of Information & Communication Technology | 宇宙航行システムおよび方法 |
WO2007008195A1 (en) * | 2004-10-28 | 2007-01-18 | University Of Maryland | Navigation system and method using modulated celestial radiation sources |
CN102243311A (zh) * | 2011-04-15 | 2011-11-16 | 中国人民解放军国防科学技术大学 | 一种x射线脉冲星导航使用的选星方法 |
-
2013
- 2013-04-07 CN CN201310117103.1A patent/CN103196451B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005043189A (ja) * | 2003-07-28 | 2005-02-17 | National Institute Of Information & Communication Technology | 宇宙航行システムおよび方法 |
WO2007008195A1 (en) * | 2004-10-28 | 2007-01-18 | University Of Maryland | Navigation system and method using modulated celestial radiation sources |
CN102243311A (zh) * | 2011-04-15 | 2011-11-16 | 中国人民解放军国防科学技术大学 | 一种x射线脉冲星导航使用的选星方法 |
Non-Patent Citations (3)
Title |
---|
SUNEEL I.SHEIKH,DARRYJ.PINES,PAUL S.RAY,KENT S.WOOD: "Spacecraft Navigation Using X-Ray Pulsars", 《JOUUNAL OF GUIDANCE,CONTROL,AND DYNAMICS》, vol. 29, no. 1, 28 February 2006 (2006-02-28), pages 49 - 63 * |
褚永辉,王大轶,黄翔宇: "脉冲星导航中最优脉冲星组合选取方法", 《中国空间科学技术》, no. 5, 31 October 2011 (2011-10-31), pages 64 - 69 * |
谢振华,许录平,倪广仁: "基于最大似然的X射线脉冲星空间定位研究", 《宇航学报》, vol. 28, no. 6, 30 November 2007 (2007-11-30), pages 1605 - 1608 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103954287A (zh) * | 2014-03-31 | 2014-07-30 | 北京控制工程研究所 | 深空探测自主导航的一种路标规划方法 |
CN103954287B (zh) * | 2014-03-31 | 2016-08-17 | 北京控制工程研究所 | 深空探测自主导航的一种路标规划方法 |
CN103913173A (zh) * | 2014-04-18 | 2014-07-09 | 中国人民解放军国防科学技术大学 | 单x射线脉冲星导航观星序列选择方法 |
CN103913173B (zh) * | 2014-04-18 | 2017-03-01 | 中国人民解放军国防科学技术大学 | 单x射线脉冲星导航观星序列选择方法 |
CN106871909A (zh) * | 2017-02-20 | 2017-06-20 | 中国人民解放军国防科学技术大学 | 一种多航天器***下基于Fisher信息矩阵的脉冲星选星方法 |
CN115355915A (zh) * | 2022-10-19 | 2022-11-18 | 中国人民解放军国防科技大学 | 一种卫星方向矢量导航的观星序列规划方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103196451B (zh) | 2015-10-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102591343B (zh) | 基于两行根数的卫星轨道维持控制方法 | |
CN102565814B (zh) | 卫星导航***的信号精度和定位服务可用性的评估方法 | |
CN102175241B (zh) | 一种火星探测器巡航段自主天文导航方法 | |
CN104240541B (zh) | 一种4d航迹生成方法 | |
CN105066994A (zh) | 一种嵌入式大气数据***与惯性导航***的数据融合方法 | |
CN104236546A (zh) | 一种卫星星光折射导航误差确定与补偿方法 | |
Neubauer et al. | Titan's near magnetotail from magnetic field and electron plasma observations and modeling: Cassini flybys TA, TB, and T3 | |
CN110132261A (zh) | 一种基于数值拟合的高精度星上轨道预报方法 | |
CN104848862B (zh) | 一种环火探测器精密同步定位守时方法及*** | |
CN101968542A (zh) | 一种利用地球站对月球探测器进行跟踪的方法 | |
CN105116444B (zh) | 一种地面微地震监测各向异性速度模型 | |
CN103196445A (zh) | 基于匹配技术的地磁辅助惯性的载体姿态测量方法 | |
CN103196451A (zh) | 一种基于Fisher信息矩阵的导航脉冲星选择方法 | |
CN103234538A (zh) | 一种行星最终接近段自主导航方法 | |
CN103105624A (zh) | 基于数据库技术的纵横波时差定位方法 | |
CN105509750A (zh) | 一种天文测速与地面无线电组合的火星捕获段导航方法 | |
CN102944238B (zh) | 一种行星探测器接近目标过程中相对位置确定方法 | |
CN103744052A (zh) | 一种应用于空中目标定位的双星测时差测向方法及装置 | |
CN105487405B (zh) | 低低跟踪重力测量卫星半物理仿真*** | |
CN103217161A (zh) | 一种脉冲星导航位置和速度联合估计方法 | |
CN103759726A (zh) | 一类循环平稳泊松信号快速模拟方法及其硬件*** | |
CN108020866A (zh) | 一种星体重力场反演的方法和***、以及处理器 | |
CN103852812B (zh) | 一种前向模拟掩星方法 | |
Weng et al. | Advanced data assimilation for cloud-resolving hurricane initialization and prediction | |
Zheng et al. | Improved forecast skill through the assimilation of dropsonde observations from the Atmospheric River Reconnaissance program |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151014 Termination date: 20180407 |
|
CF01 | Termination of patent right due to non-payment of annual fee |