CN103940424A - 一种x射线脉冲星导航信号整周模糊度检测与估计方法 - Google Patents
一种x射线脉冲星导航信号整周模糊度检测与估计方法 Download PDFInfo
- Publication number
- CN103940424A CN103940424A CN201410147927.8A CN201410147927A CN103940424A CN 103940424 A CN103940424 A CN 103940424A CN 201410147927 A CN201410147927 A CN 201410147927A CN 103940424 A CN103940424 A CN 103940424A
- Authority
- CN
- China
- Prior art keywords
- integer ambiguity
- pulsar
- tdoa
- moment
- xnav
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/02—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供一种X射线脉冲星导航信号整周模糊度检测与估计方法。技术方案是一个递推解算过程,在任意的测量数据获取时刻k,k>0,首先修正XNAV信号整周模糊度,然后估计航天器轨道和接收机钟差,利用TDOA测量残差采用稀疏寻优方法检测XNAV信号的整周模糊度,再采用标准Kalman滤波方法估计XNAV信号的整周模糊度的浮点解,最后采用四舍五入的取整运算计算k时刻的整周模糊度。本发明过程实现简单,实时性强,并可有效降低检测的漏警率和虚警率,提高整周模糊度检测的稳健性和整周模糊度估计的精度。
Description
技术领域
本发明涉及导航信息获取与处理技术领域,特别涉及一种X射线脉冲星导航信号的整周模糊度检测与估计方法。
背景技术
XNAV(X-ray pulsar-based navigation,X射线脉冲星导航)是近年来发展迅速的一种天文导航技术,它基于X射线脉冲星辐射的高稳定脉冲信号和已知的脉冲星空间位置信息,为观测者提供航天器的位置、速度、姿态乃至时间等运动参数信息,相比地面雷达和GNSS(Global Navigation Satellite System,全球导航卫星系统)等传统导航技术,具有完全自主、稳定、可靠和高精度等诸多优势。因此,将XNAV系统应用于航天器自主导航,特别是深空探测中的行星际航天器自主导航,有重要的理论和实际意义。
X射线脉冲星导航信号是脉冲星辐射信号到航天器和空间参考点(如太阳系质心、地球质心等)的TDOA(time difference of arrival,到达时间差)测量,它是导航接收机获取的脉冲星辐射信号经过降噪和相位估计等处理后获得的,其整周模糊度的快速检测和准确估计是提高XANV导航性能的关键因素之一。TDOA测量的整周模糊度解算是一个三维空间解模糊问题,传统的空间搜索方法计算复杂度很高。为提高算法实时性,有学者提出直接利用TDOA测量或航天器轨道预测值估计整周模糊度的方法,然而由于TDOA测量误差较大,以及系统工作初期或发生故障时轨道预测精度较低,容易出现整周模糊度的漏检和虚警,可靠性较低,这就要求我们提出高效稳健的整周模糊度检测与估计方法。
发明内容
本发明要解决的技术问题是,提供一种高效稳健的XNAV信号整周模糊度检测和估计方法。本方法利用XNAV信号整周模糊度的相位整周特性和稀疏变化特性,将整周模糊度的检测和估计看成系统误差的估计问题,采用稀疏寻优方法,快速高效的完成整周模糊度的检测与估计。
本发明技术方案的基本思路是,由于XNAV信号整周模糊度的出现和周跳变化通常只在少数几个脉冲星测量通道内发生,因此通过将整周模糊度等效为一种系统测量误差,在每个新的TDOA测量时刻,首先利用之前估计得到的整周模糊度修正TDOA测量,然后估计航天器的轨道和接收机钟差,再利用TDOA测量残差进行稀疏寻优,检测每个测量通道内整周模糊度的存在性,最后利用TDOA测量残差估计得到新的XNAV信号整周模糊度。
本发明的技术方案是,已知XNAV系统具有N颗脉冲星的测量通道,各测量通道的信号积累周期为ΔTi,i=1,…,i,N,第i颗脉冲星TDOA测量的信号周期为Ti,i=1,…,i,N,第i颗脉冲星TDOA测量的噪声标准差为σi,i=1,…,i,N,第i颗脉冲星TDOA测量的整周模糊度估计值的初值在任意的测量数据获取时刻k,k>0,航天器接收到N颗X射线脉冲星的TDOA测量值ti,k,i=1,2,...,N,进行下述处理步骤:
步骤一:修正XNAV信号整周模糊度。
将TDOA测量值ti,k减去k-1时刻的整周模糊度估计值得到TDOA测量修正量t'i,k。
步骤二:估计航天器轨道和接收机钟差。
根据航天器轨道动力学特征和接收机钟差的时间演化特性,采用典型的脉冲星非线性导航滤波方法(具体方法参照《关于X射线脉冲星导航的轨道力学问题》,帅平等,中国科学E辑,2009.3,第556页到第561页),设计航天器状态滤波器,以TDOA测量修正量t'i,k,i=1,2,...,N为滤波器输入,估计航天器状态Xk=[rk,vk,τk]T,其中rk表示k时刻航天器的轨道位置,vk表示k时刻航天器的轨道速度,τk表示k时刻的导航接收机钟差,并记滤波输出的第i颗脉冲星TDOA测量残差为Δti,k,i=1,…,N。
步骤三:检测XNAV信号的整周模糊度。
计算第i颗脉冲星的归一化TDOA测量残差yi,k:
yi,k=Δti,k/σi,i=1,…,N (公式一)
对yk=[y1,k,…,yi,k…,yN,k]采用基于匹配追踪的稀疏寻优方法(具体方法参照《FastBayesian Matching Pursuit》,Philip Schniter,Information Theory and ApplicationsWorkshop,2008,第326页到第333页),获得yk中幅度大于α||yk||的分量,相应通道的脉冲星组成k时刻可能存在整周模糊度的待选脉冲星序列。其中,稀疏度百分比α(0<α<1)作为稀疏度标准,根据实际情况确定。当k<5时,各通道TDOA测量的整周模糊度估计值保持为初值0,结束k时刻的数据处理步骤,转入k+1时刻的数据处理步骤;当k≥5时,若某脉冲星在最近5个数据获取时刻有不少于3次位于待选脉冲星序列内,则判定该脉冲星在k时刻的TDOA测量存在整周模糊度。
步骤四:估计XNAV信号的整周模糊度。
将k时刻XNAV信号的整周模糊度Sk=[s1,k,…si,k,…,sN,k]作为TDOA测量的系统误差,建立如下的整周模糊度状态方程:
Sk=CkSk-1+Uk (公式二)
其中Ck为N×N的对角矩阵,步骤三判定在k时刻存在整周模糊度的脉冲星序号对应的Ck对角元素取值为1,否则取值为0,Sk-1为k-1时刻XNAV信号的整周模糊度,Uk=[u1,k,…ui,k,…,uN,k]为零均值高斯白噪声向量,其第i个元素ui,k的标准差为βσi,其中比例系数β可根据实际情况确定。
建立k时刻的整周模糊度测量方程:
Δti,k=si,k+εk i=1,2,...,N (公式三)
其中εk是标准差为σi的零均值高斯白噪声向量。
利用公式二和公式三采用标准Kalman滤波方法,估计得到XNAV信号的整周模糊度浮点解最后计算k时刻的第i颗脉冲星整周模糊度如下:
其中round表示四舍五入的取整函数。
本发明的有益效果是:(1)采用基于匹配追踪的稀疏寻优方法进行整周模糊度检测,过程实现简单,计算量远小于传统的空间搜索法,实时性强;(2)将航天器状态滤波器的TDOA测量残差作为整周模糊度检测与估计的输入,可减小原始TDOA测量中随机误差和粗差的影响,有效降低检测的漏警率和虚警率,提高整周模糊度检测的稳健性和整周模糊度估计的精度。
附图说明
图1是本发明的原理流程示意图;
图2是本发明实施例输入的各脉冲星测量通道的相关参数;
图3是本发明实施例输入的航天器轨道根数;
图4是整周模糊度检测和估计结果;
图5是航天器轨道位置估计误差;
图6是航天器轨道速度估计误差;
图7是航天器接收机钟差估计误差。
具体实施方式
图1是本发明的原理流程示意图。技术方案是一个递推解算过程,在任意的测量数据获取时刻k,k>0,首先修正XNAV信号整周模糊度,然后估计航天器轨道和接收机钟差,利用TDOA测量残差采用稀疏寻优方法检测XNAV信号的整周模糊度,再采用标准Kalman滤波方法估计XNAV信号的整周模糊度的浮点解,最后采用四舍五入的取整运算计算k时刻的整周模糊度。
图2至图7是利用本发明进行仿真实验的结果。
仿真实验中,航天器XNAV系统的脉冲星测量通道数N为6,航天器轨道根数见图2所示,6颗脉冲星TDOA测量通道的信号周期等相关参数见图3所示。导航接收机钟差τk为1微秒的常值,各通道信号积累周期ΔTi均为500s,各脉冲星TDOA测量的噪声标准差σi均为1.267微秒,仿真时长为48小时。第1颗脉冲星的TDOA测量在80000s至100000s时间段内出现1个周期的整周模糊度;第2颗脉冲星的TDOA测量在全仿真时段内有4个周期的整周模糊度;第3颗脉冲星在20000s至120000s时间段内出现1个周期的整周模糊度,在120000s至仿真结束出现3个周期的整周模糊度;第4、5、6颗脉冲星均未出现整周模糊度。步骤三中参数α取0.1,步骤四中参数β取1000。
图4为6颗脉冲星TDOA测量的整周模糊度检测与估计结果。其中横坐标为时间(105s),纵坐标为整周模糊度对应的周期数(以下简称为模糊数),空心圆圈形状的曲线为模糊数真值,实心圆点形状的曲线为模糊数估计值,箭头指向的小块曲线是模糊数曲线跳变处的局部放大图。第1、3、5、6颗脉冲星均获得了理想的整周模糊度检测与估计结果;第2颗脉冲星在经历了一段收敛时间后也得到理想的估计结果;第4颗脉冲星仅在初始阶段出现了短暂的-1个周期的整周模糊度虚警。总体上,6颗脉冲星均可快速高精度的获得令人满意的整周模糊度检测与估计结果。
图5、图6和图7分别为航天器XNAV的轨道位置、轨道速度和接收机钟差的估计误差曲线。其中横坐标为时间(s),纵坐标为各项估计误差,实线为测量数据中无整周模糊度时的估计误差,虚线为测量数据中有整周模糊度时并采用本方法后的估计误差。从虚线中可以看出,导航滤波器在每次整周模糊度变化后快速收敛,各项误差与无整周模糊度时相当,取得了令人满意的导航精度和稳健性。
Claims (2)
1.一种X射线脉冲星导航信号整周模糊度检测与估计方法,X射线脉冲星导航记为XNAV,已知XNAV系统具有N颗脉冲星的测量通道,各测量通道的信号积累周期为ΔTi,i=1,…,i,N,第i个测量通道即第i颗脉冲星的TDOA测量的信号周期为Ti,,其中TDOA表示到达时间差,第i颗脉冲星TDOA测量的噪声标准差为σi,第i颗脉冲星TDOA测量的整周模糊度估计值的初值在任意的测量数据获取时刻k,k>0,航天器接收到N颗X射线脉冲星的TDOA测量值ti,k,其特征在于,进行下述处理步骤:
步骤一:修正XNAV信号整周模糊度:
将TDOA测量值ti,k减去k-1时刻的整周模糊度估计值得到TDOA测量修正量t'i,k;
步骤二:估计航天器轨道和接收机钟差:
根据航天器轨道动力学特征和接收机钟差的时间演化特性,采用典型的脉冲星非线性导航滤波方法设计航天器状态滤波器,以TDOA测量修正量t'i,k为滤波器输入,估计航天器状态Xk=[rk,vk,τk]T,其中rk表示k时刻航天器的轨道位置,vk表示k时刻航天器的轨道速度,τk表示k时刻的导航接收机钟差,并记滤波输出的第i颗脉冲星TDOA测量残差为Δti,k;
步骤三:检测XNAV信号的整周模糊度:
计算第i颗脉冲星的归一化TDOA测量残差yi,k:
yi,k=Δti,k/σi (公式一)
对yk=[y1,k,…,yi,k…,yN,k]采用基于匹配追踪的稀疏寻优方法获得yk中幅度大于α||yk||的分量,相应通道的脉冲星组成k时刻可能存在整周模糊度的待选脉冲星序列;其中,稀疏度百分比α作为稀疏度标准,根据实际情况确定且0<α<1;当k<5时,各通道TDOA测量的整周模糊度估计值保持为初值0,结束k时刻的数据处理步骤,转入k+1时刻的数据处理步骤;当k≥5时,若某脉冲星在最近5个数据获取时刻有不少于3次位于待选脉冲星序列内,则判定该脉冲星在k时刻的TDOA测量存在整周模糊度;
步骤四:估计XNAV信号的整周模糊度:
将k时刻XNAV信号的整周模糊度Sk=[s1,k,…si,k,…,sN,k]作为TDOA测量的系统误差,建立如下的整周模糊度状态方程:
Sk=CkSk-1+Uk (公式二)
其中Ck为N×N的对角矩阵,步骤三判定在k时刻存在整周模糊度的脉冲星序号对应的Ck对角元素取值为1,否则取值为0,Sk-1为k-1时刻XNAV信号的整周模糊度,Uk=[u1,k,…ui,k,…,uN,k]为零均值高斯白噪声向量,其第i个元素ui,k的标准差为βσi,其中比例系数β可根据实际情况确定;
建立k时刻的整周模糊度测量方程:
Δti,k=si,k+εk i=1,2,...,N (公式三)
其中εk是标准差为σi的零均值高斯白噪声向量;
利用公式二和公式三采用标准Kalman滤波方法,估计得到XNAV信号的整周模糊度浮点解最后计算k时刻的第i颗脉冲星整周模糊度如下:
其中round表示四舍五入的取整函数。
2.根据权利要求1所述的X射线脉冲星导航信号整周模糊度检测与估计方法,其特征在于,步骤三中参数α的取值是0.1,步骤四中参数β的取值是1000。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410147927.8A CN103940424B (zh) | 2014-04-14 | 2014-04-14 | 一种x射线脉冲星导航信号整周模糊度检测与估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410147927.8A CN103940424B (zh) | 2014-04-14 | 2014-04-14 | 一种x射线脉冲星导航信号整周模糊度检测与估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103940424A true CN103940424A (zh) | 2014-07-23 |
CN103940424B CN103940424B (zh) | 2016-09-14 |
Family
ID=51188193
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410147927.8A Active CN103940424B (zh) | 2014-04-14 | 2014-04-14 | 一种x射线脉冲星导航信号整周模糊度检测与估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103940424B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104864875A (zh) * | 2015-04-03 | 2015-08-26 | 北京控制工程研究所 | 一种基于非线性h∞滤波的航天器自主定位方法 |
CN106871909A (zh) * | 2017-02-20 | 2017-06-20 | 中国人民解放军国防科学技术大学 | 一种多航天器系统下基于Fisher信息矩阵的脉冲星选星方法 |
CN107421533A (zh) * | 2017-06-22 | 2017-12-01 | 北京航空航天大学 | 一种深空探测器x射线脉冲星toa/dtoa组合导航方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5596511A (en) * | 1993-01-27 | 1997-01-21 | Fuji Xerox Co., Ltd. | Computing method and apparatus for a many-body problem |
CN101178312A (zh) * | 2007-12-12 | 2008-05-14 | 南京航空航天大学 | 基于多信息融合的航天器组合导航方法 |
CN103674022A (zh) * | 2013-12-19 | 2014-03-26 | 中国空间技术研究院 | 一种快速脉冲星导航整周模糊度解算方法 |
-
2014
- 2014-04-14 CN CN201410147927.8A patent/CN103940424B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5596511A (en) * | 1993-01-27 | 1997-01-21 | Fuji Xerox Co., Ltd. | Computing method and apparatus for a many-body problem |
CN101178312A (zh) * | 2007-12-12 | 2008-05-14 | 南京航空航天大学 | 基于多信息融合的航天器组合导航方法 |
CN103674022A (zh) * | 2013-12-19 | 2014-03-26 | 中国空间技术研究院 | 一种快速脉冲星导航整周模糊度解算方法 |
Non-Patent Citations (3)
Title |
---|
孙景荣等: "一种基于遗传算法的XPNAV整周模糊度求解方法", 《第三届中国卫星导航学术年会电子文集》 * |
谢强等: "脉冲星导航解周期模糊匹配搜索算法", 《系统工程与电子技术》 * |
黄震等: "脉冲星导航的整周模糊度解算方法研究", 《空间控制技术与应用》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104864875A (zh) * | 2015-04-03 | 2015-08-26 | 北京控制工程研究所 | 一种基于非线性h∞滤波的航天器自主定位方法 |
CN104864875B (zh) * | 2015-04-03 | 2018-01-05 | 北京控制工程研究所 | 一种基于非线性h∞滤波的航天器自主定位方法 |
CN106871909A (zh) * | 2017-02-20 | 2017-06-20 | 中国人民解放军国防科学技术大学 | 一种多航天器系统下基于Fisher信息矩阵的脉冲星选星方法 |
CN107421533A (zh) * | 2017-06-22 | 2017-12-01 | 北京航空航天大学 | 一种深空探测器x射线脉冲星toa/dtoa组合导航方法 |
CN107421533B (zh) * | 2017-06-22 | 2019-07-30 | 北京航空航天大学 | 一种深空探测器x射线脉冲星toa/dtoa组合导航方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103940424B (zh) | 2016-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109738917B (zh) | 一种北斗变形监测中的多路径误差削弱方法及装置 | |
CN104296755B (zh) | 一种x射线脉冲星导航脉冲toa的确定方法 | |
CN101403790B (zh) | 单频gps接收机的精密单点定位方法 | |
CN108919321B (zh) | 一种基于尝试法的gnss定位粗差探测方法 | |
EP2570823A1 (en) | Method and apparatus for differential global positioning system (DGPS) - based real time attitude determination (RTAD) | |
JP5352422B2 (ja) | 測位装置及びプログラム | |
CN103217162B (zh) | 采用稀疏表示的脉冲星累积脉冲轮廓时间延迟测量方法 | |
CN109459778A (zh) | 基于抗差方差分量估计的码伪距/多普勒联合测速方法及其应用 | |
CN105372691A (zh) | 一种模糊度固定的长基线卫星编队gnss相对定位方法 | |
CN101561499B (zh) | 一种单站多普勒测距定位方法 | |
CN103017774A (zh) | 单探测器脉冲星导航方法 | |
Meng et al. | Design and experiment of onboard laser time transfer in Chinese Beidou navigation satellites | |
CN103217161B (zh) | 一种脉冲星导航位置和速度联合估计方法 | |
CN104678371B (zh) | 一种基于时延修正的海面高度测量装置 | |
CN104316048A (zh) | 一种普适性的脉冲星自主导航测量模型构建方法 | |
CN116774264B (zh) | 基于低轨卫星机会信号多普勒的运动目标定位方法 | |
CN103576167A (zh) | 基于hht和支持向量机的周跳探测与修复方法 | |
Tinto et al. | Time-delay interferometric ranging for space-borne gravitational-wave detectors | |
CN103529482A (zh) | 一种高精度确定载体动态加速度的方法 | |
CN103940424A (zh) | 一种x射线脉冲星导航信号整周模糊度检测与估计方法 | |
Bramanto et al. | Long-range single baseline RTK GNSS positioning for land cadastral survey mapping | |
CN106371092B (zh) | 一种基于gps与强震仪观测自适应组合的形变监测方法 | |
CN113031036B (zh) | 基于GNSS 30s采样频率数据的电离层相位闪烁因子构建方法 | |
CN105738931A (zh) | 基于卡尔曼滤波的gps单点定位系统 | |
CN107144274B (zh) | 在轨x射线脉冲星计时模型构建方法 |
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 |