CN103389099B - 基于x射线脉冲星的航天器姿态和位置测量系统及方法 - Google Patents
基于x射线脉冲星的航天器姿态和位置测量系统及方法 Download PDFInfo
- Publication number
- CN103389099B CN103389099B CN201310320673.0A CN201310320673A CN103389099B CN 103389099 B CN103389099 B CN 103389099B CN 201310320673 A CN201310320673 A CN 201310320673A CN 103389099 B CN103389099 B CN 103389099B
- Authority
- CN
- China
- Prior art keywords
- pulsar
- satellite
- borne
- ray detector
- attitude
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 230000005855 radiation Effects 0.000 claims abstract description 57
- 238000005259 measurement Methods 0.000 claims description 24
- 230000004907 flux Effects 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000001514 detection method Methods 0.000 claims description 6
- 230000002708 enhancing effect Effects 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 230000014509 gene expression Effects 0.000 claims description 4
- 230000003247 decreasing effect Effects 0.000 claims description 2
- 230000004927 fusion Effects 0.000 abstract 1
- 230000032258 transport Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 7
- 238000011160 research Methods 0.000 description 7
- 230000009466 transformation Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- MKZGVLPHKXXSSG-UHFFFAOYSA-N ethyl n-[4-[benzyl(2-phenylethyl)amino]-2-[4-(trifluoromethyl)phenyl]-1h-imidazo[4,5-c]pyridin-6-yl]carbamate Chemical compound N=1C(NC(=O)OCC)=CC=2NC(C=3C=CC(=CC=3)C(F)(F)F)=NC=2C=1N(CC=1C=CC=CC=1)CCC1=CC=CC=C1 MKZGVLPHKXXSSG-UHFFFAOYSA-N 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- RKWPMPQERYDCTB-UHFFFAOYSA-N ethyl n-[4-[benzyl(2-phenylethyl)amino]-2-(4-nitrophenyl)-1h-imidazo[4,5-c]pyridin-6-yl]carbamate Chemical compound N=1C(NC(=O)OCC)=CC=2NC(C=3C=CC(=CC=3)[N+]([O-])=O)=NC=2C=1N(CC=1C=CC=CC=1)CCC1=CC=CC=C1 RKWPMPQERYDCTB-UHFFFAOYSA-N 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
Landscapes
- Navigation (AREA)
Abstract
本发明公开了一种基于X射线脉冲星的航天器姿态和位置测量系统及方法,所述系统利用星载X射线探测器配合准直器获取脉冲到达时间和精确获取脉冲星辐射方向矢量,并利用红外或紫外地平敏感器获取地心方向矢量;利用脉冲星脉冲到达时间、辐射方向矢量和地心方向矢量,导航计算机分别完成地心坐标系下脉冲星计时观测导航算法、基于角距测量的脉冲星定位算法和基于角距测量的脉冲星定姿算法,并通过几种算法的融合实时获取航天器位置、速度、姿态以及在轨运动方向,完成轨道卫星自主导航,输出高精度的位置、姿态和速度信息。本发明适用于地球或其它行星的轨道卫星或星座的自主运控。
Description
技术领域
本发明涉及导航技术领域,具体涉及基于X射线脉冲星的航天器姿态和位置测量系统及其方法,所述方法用于为近地轨道航天器进行高精度自主导航以及姿态测量提供服务。
背景技术
X射线脉冲星导航(XPNAV)是通过测量脉冲星辐射的X射线光子到达时间和脉冲星影像角位置来更新航天器位置、速度、时间和姿态等导航参数,根据光子到达时间计算得到导航过程中的基本信息——脉冲到达时间(TOA),利用这一信息可以计算航天器的位置、速度,并对星上的时钟进行校正。根据光子到达星载X射线探测器的位置,则可以确定脉冲星的方位信息,从而完成对航天器姿态的测量。
1974年,美国喷气推进实验室的Dr.Downs首次提出基于射电脉冲星的行星际飞行航天器自主轨道确定方法,这是一种建立在脉冲星无线电信号基础上的在轨航天器导航方法。随着天文观测能力的增强以及对X射线脉冲星特性研究的深入,1981年美国通信系统研究所的T.J.Chester和S.A.Butman提出了使用X射线波段的脉冲星为航天器导航的构想。1982年,第一颗毫秒级脉冲星被确认,人们发现利用某些毫秒脉冲星的频率稳定性,可以将其用作天然高精度的频率标准源。1999年,美国ARGOS卫星发射升空,初步验证了脉冲星自主导航的可能性,为其发展与应用奠定了重要的基础。2004年,美国航空航天局和海军天文台等多家单位着手拟定和启动脉冲星导航的研究计划,同时X射线脉冲星导航已纳入国防部长期发展战略规划纲要,并逐年增加项目研究经费,持续开展脉冲星导航的理论研究方法、关键技术攻关和原理样机研制等方面的研究工作。近年来,我国在X射线脉冲星自主导航方面的研究工作发展很快,从最初跟踪国外的科研项目,到现在能够自主创新、提出越来越多的观点和方法。2006年1月,中国科学院高能物理研究所和清华大学天体物理研究中心共同设计研制了硬X射线调制望远镜HXMT,目标在于实现宽波段X射线(1~250keV)巡天,探索利用X射线脉冲星实现航天器自主导航的技术和原理。
但是,脉冲星信号极其微弱,而从目前的脉冲星导航原理来看,在太阳系质心坐标系中,X射线脉冲星导航的基本观测量是脉冲星计时模型的预报时间与卫星观测时间之差,换算为距离之后进行定位,这种方法受诸多因素的影响,定位精度不高。并且传统的脉冲星定姿方案中,需要使用两个星载X射线探测器扫描同一颗星或者使用一个星载X射线探测器观测两个脉冲星。基于这一原理,在实际中搭载两个星载X射线探测器的方法成本太高,必须考虑单星载X射线探测器定姿。
发明内容
针对上述问题,本发明的主要目的在于克服上述已有技术的不足,提出一种基于X射线脉冲星的航天器姿态和位置测量系统及其方法,从而实现单星载X射线探测器观测单颗脉冲星完成航天器自主姿态和位置的联合测量,以克服传统利用脉冲星对三轴稳定航天器定姿中单星载X射线探测器必须分时观测两颗脉冲星而引入的附加噪声和频繁调整星载X射线探测器指向等问题。本发明还能在姿态测量同时,完成航天器指向脉冲星辐射矢量与地平观测矢量夹角的测量,从而进一步进行脉冲星定位,此定位方法弥补了计时观测在航天器定位中精度差的缺点。因此,本发明所提出的联合测量方法能充分利用脉冲星的导航能力,降低成本,实现航天器自主高精度姿态测量和轨道确定。
为了实现上述目的,本发明采用的技术方案如下:
基于X射线脉冲星的航天器姿态和位置测量系统,其特征在于,所述系统包括:星载X射线探测器、准直器、星载原子钟、X射线脉冲星特征参数数据库、X射线脉冲星辨识算法库、太阳系行星参数数据库、星载计算机、地平敏感仪,导航算法库、星载万向支架等;其中
所述星载计算机根据当前所记录的航天器姿态和所述X射线脉冲星特征参数数据库,控制所述星载万向支架支撑星载X射线探测器指向目标脉冲星;
所述星载X射线探测器配合所述准直器,测量到达所述星载X射线探测器的X射线脉冲星信号的光子流量强度,从而计算得出所述星载X射线探测器指向与脉冲星辐射矢量的偏差,将该偏差反馈到所述星载计算机用于进行航天器姿态的控制增强,不断调整所述星载X射线探测器指向控制,以获得精确地脉冲星辐射方向矢量;
所述地平敏感仪敏感地平得到地平线张角,从而计算地心方向矢量,将该地心方向矢量和所述脉冲星辐射方向矢量共同送至所述星载计算机中的脉冲星定姿算法,同时将所述地平线张角和所述脉冲星辐射方向矢量共同送至脉冲星定位算法;
所述星载计算机调用所述导航算法库中的联合定位算法对定姿算法和定位算法进行结合,输出高精度的自主导航参数信息,包括姿态和位置。
需要说明的是,所述星载X射线探测器与所述准直器配合用于接收脉冲星辐射的X射线,所述准直器安装于所述星载X射线探测器的前部;所述准直器由若干空心圆柱体互相紧密排列组合而成。
一种使用基于X射线脉冲星的航天器姿态和位置联合测量系统测量的方法,所述方法包括以下步骤:
(1)探测目标脉冲星:所述星载计算机根据当前所记录的航天器姿态和所述X射线脉冲星特征参数数据库,控制所述星载万向支架支撑星载X射线探测器指向目标脉冲星;
(2)调整所述星载X射线探测器方向:所述星载X射线探测器配合所述准直器,测量X射线脉冲星信号的光子流量强度,并计算得出所述星载X射线探测器指向与脉冲星辐射矢量的偏差;
(3)校正目标脉冲星位置:将步骤(2)获得的偏差反馈到所述星载计算机用于进行航天器姿态的控制增强,从而不断调整所述星载X射线探测器指向控制,以获得精确地脉冲星辐射方向矢量;
(4)确定航天器姿态与位置:所述地平敏感仪敏感地平得到地平线张角,从而计算地心方向矢量,将该地心方向矢量和所述脉冲星辐射方向矢量共同送至所述星载计算机中的脉冲星定姿算法,同时将所述地平线张角和所述脉冲星辐射方向矢量共同送至脉冲星定位算法;所述星载计算机调用所述导航算法库中的联合定位算法对定姿算法和定位算法进行结合,输出高精度的自主导航参数信息,包括姿态和位置。
需要说明的是,所述步骤(2)中的偏差获取方法如下:首先定义两个坐标系,分别是B-系和C-系,B-系的坐标原点为航天器的质心,其X轴指向航天器的前进方向,Z轴垂直于偏航平面,Y轴根据右手定则确定,即(XB,YB,ZB);C-系(XC,YC,ZC)中,YC平行于支架的水平轴,ZC指向星载X射线探测器的中心轴,XC垂直于YC和ZC的平面;明显地,所述星载X射线探测器指向控制矢量ZC与脉冲星辐射矢量的夹角θ是所要求的偏差,其中,C-系与B-系的坐标原点相同。
需要说明的是,所述星载X射线探测器指向控制的实现方法如下:
根据三角关系,圆筒形准直器的有效探测面积为:
使得:
其中,θ为星载X射线探测器指向与脉冲星辐射矢量的夹角,H为准直器的高,R为准直器横截圆面的半径。
需要说明的是,所述脉冲星辐射方向矢量的搜索方法如下:所述星载计算机控制所述星载万向支架的方向,同时控制所述星载X射线探测器指向,所述星载X射线探测器记录敏感到的脉冲星辐射光子流量,通过改变其的指向ZC来减小θ,当光子流量强度最大时,即可认为所述星载X射线探测器指向与脉冲星辐射方向矢量平行;所述星载万向支架的水平轴XC只能在(XC,YC)平面内以ZB为轴进行旋转,α用来表示旋转角度;由于ZC与XC是垂直的,那么ZC以XC为轴进行旋转时,两个轴在线性独立方向的旋转角度都是θ,同时要求脉冲星光子流量强度作为目标函数是线性递减的;其中包括:
(1)根据姿态信息的先验知识,所述星载X射线探测器根据参考星历对目标脉冲星在B-系中的位置进行粗估计,所述星载X射线探测器旋转中心轴使其指向脉冲星以获取脉冲星信号。同时所述星载计算机使用FFT或是贝叶斯定理通过捕获到的光子序列Ph(θ)来检测目标信号,检测成功后进行步骤(2);
(2)在B-系中,XC被表示为矢量(cosα,sinα,0)Τ,如果XC以ZB为轴旋转了Δα,则新的矢量可以表示为(cos(α+Δα),sin(α+Δα),0)Τ;设β为ZC以XC为轴旋转的角度,那么
设为ZC的单位向量,则ZC绕XC旋转角度β可以表示为 具体是:
(2a)设置初始点x(0)=(α0,β0)以及两个线性独立方向α1,β1,允许误差ε≥0且k=1;
(2b)设置x(k,0)=x(k-1)。从x(k,0)沿着αk方向使用线性搜索方法:
结果为x(k,1),同样通过沿着βk方向搜索可以得到xk点;
(2c)如果||ph(θ)-phmax||≤ε,则进行步骤(3),其中phmax是所述星载X射线探测器记录的平均最大流量强度;否则,设置x0=xk,返回步骤(2a)。
(3)所述星载X射线探测器持续观测目标脉冲星,星载计算机计算得到平均流量强度;如果||ph(θ)-phmax||>ε,进行步骤(2)与(2a),并再一次调整指向。
需要说明的是,所述航天器确定姿态实现方法如下:将获得的地心方向矢量和脉冲星辐射方向矢量,在惯性坐标系中,利用欧拉角建立姿态变换矩阵,即:
需要说明的是,所述航天器确定位置实现方法如下:所述星载计算机接收所述星载X射线探测器输出的脉冲星辐射方向矢量和所述地平敏感仪输出的地平线张角和地心方向矢量,建立两个关系式A和B;其中关系式A为地平线张角与卫星位置矢量的函数,即关系式B为地心方向矢量与脉冲星方向矢量之间的夹角以及卫星位置矢量的函数,即同时,利用脉冲星计时观测通过累计轮廓与标准相位轮廓比相而得到相差,从而计算距离信息;将该距离信息与前述的关系式A和B一起送入导航滤波算法库,完成定位。
需要说明的是,所述导航滤波算法实现方法为:设地平敏感仪观测到的地平线张角为α,卫星地心方向矢量与脉冲星辐射方向矢量的夹角为β,建立观测方程:
其中ns包括泰勒展开的高阶项和噪声,nh包括泰勒展开的高阶项和地心矢量观测噪声,E[ns(t)]=0,E[nh(t)]=0,同时脉冲星计时观测建立的观测方程为:
综合(1)、(2)和(3)式建立组合观测方程
其中,H为观测矩阵,η(t)为零均值高斯观测噪声且η(t)=[ns(t),nh(t),ωN(t)]Τ, 其中ωN(t)为第N颗脉冲星的TOA量测噪声。
一种利用基于X射线脉冲星的航天器姿态和位置联合测量的方法,适用于地球外的其它近天体轨道航天器。
本发明有益效果在于:实现单星载X射线探测器观测单颗脉冲星完成航天器自主姿态和位置的联合测量,以克服传统利用脉冲星对三轴稳定航天器定姿中单星载X射线探测器必须分时观测两颗脉冲星而引入的附加噪声和频繁调整星载X射线探测器指向等问题。本发明还能在姿态测量同时,完成航天器指向脉冲星辐射矢量与地平观测矢量夹角的测量,从而进一步进行脉冲星定位,此定位方法弥补了计时观测在航天器定位中精度差的缺点。
附图说明
图1为本发明的利用X射线脉冲星的航天器姿态和位置的联合测量的几何方法系统结构框图;
图2为本发明在地球质心参考坐标系中基于计时观测和仰角观测的导航算法原理图;
图3为本发明的脉冲星辐射矢量搜索方法所建立坐标系的示意图;
图4为本发明的准直器结构示意图和有效使用面积与脉冲星辐射方向矢量和准直器指向夹角之间的关系;
图5为星载X射线探测器有效使用面积计算方法示意图。
具体实施方式
下面将结合附图对本发明作进一步的描述。
如图1所示,本发明为基于X射线脉冲星的航天器姿态和位置联合测量系统,包括星载X射线探测器4、准直器3、星载原子钟、X射线脉冲星特征参数数据库6、X射线脉冲星辨识算法库7、太阳系行星参数数据库、红外(或紫外)地平敏感仪,导航算法库、卫星姿态和位置控制平台等;脉冲星辐射出的X射线光子流1通过准直器3的选通处理,当准直器3的指向与脉冲星辐射方向一致或在选通角范围之内时,X射线光子流进入星载X射线探测器4,否则脉冲星辐射的X射线光子流无法进入星载X射线探测器4进行处理。星载X射线探测器4将接收到的X射线光子流分别送入脉冲星信号辨识模块8和记录光子到达模块9。脉冲星信号辨识模块8利用X射线脉冲星特征参数库6和X射线脉冲星辨识算法库7完成对所接收到的脉冲星信号的辨识,由于X射线脉冲星信号具有稳定的脉冲辐射周期和唯一可辨识的特征轮廓,因此具有非常好的可辨识性和可检验性。统计脉冲星信号辨识模块8输出的光子信号的流量强度,当光子流量强度最大时认为准直器3与X射线脉冲星辐射方向矢量平行,即准直器3此时所指方向为脉冲星在天球坐标系中的辐射矢量;同时,地平敏感仪敏感地平,得到地平线张角11和航天器相对于地心方向的位置矢量10,将该地心方向矢量10与脉冲星辐射方向矢量12共同送入姿态矩阵确定模块13,在惯性坐标系中,利用欧拉角建立姿态变换矩阵,用于进行航天器姿态确定;星载X射线探测器4指向控制模块14不断改变星载X射线探测器4在天球坐标系中的指向,使单位时间内到达星载X射线探测器4的X射线信号光子流量最大;本发明在完成定姿要求的同时,可以完成航天器的定位功能。基于计时观测的导航算法库15综合前述的地平线张角11、地心方向矢量、脉冲星辐射方向矢量和光子到达时间,建立观测方程,结合航天器轨道动力学完成近地轨道卫星的高精度导航;航天器姿态信息和高精度的导航信息送至卫星姿态、位置控制平台。
如图3所示,所述的星载X射线探测器4指向与脉冲星辐射方向矢量的偏差的获取方法为:首先定义两个坐标系,B-系和C-系。B-系的坐标原点为航天器的质心,其X轴指向航天器的前进方向,Z轴垂直于偏航平面,Y轴根据右手定则确定,即(XB,YB,ZB)。C-系(XC,YC,ZC)中,YC平行于支架的水平轴,ZC指向星载X射线探测器4的中心轴,XC垂直于YC和ZC的平面。显然,星载X射线探测器4指向控制矢量Zx与脉冲星辐射矢量的夹角θ就是所要求的偏差。C-系与B-系的坐标原点相同。
如图5所示,所述的星载X射线探测器4指向控制的实现方法为:根据三角关系,圆筒形准直器3的有效探测面积为
使得:
其中,θ为星载X射线探测器4指向与脉冲星辐射矢量的夹角,H为准直器的高,R为准直器横截圆面的半径。
如图3所示,所述的脉冲星辐射方向矢量的搜索方法为:计算机控制星载万向支架的方向,同时控制星载X射线探测器4指向,星载X射线探测器4记录敏感到的脉冲星辐射光子流量,通过改变星载X射线探测器4的指向ZC来减小θ,当光子流量强度最大时,即可认为星载X射线探测器4指向与脉冲星辐射方向矢量平行。支架的水平轴XC只能在(XC,YC)平面内以ZB为轴进行旋转,α用来表示旋转角度。由于ZC与XC是垂直的,那么ZC以XC为轴进行旋转时,两个轴在线性独立方向的旋转角度都是θ,同时要求脉冲星光子流量强度作为目标函数是线性递减的。具体的搜索方法为:
1、根据姿态信息的先验知识,星载X射线探测器4根据参考星历对目标脉冲星在B-系中的位置进行粗估计,星载X射线探测器旋转中心轴使其指向脉冲星以获取脉冲星信号。同时星载计算机使用FFT或是贝叶斯定理通过捕获到的光子序列Ph(θ)来检测目标信号。检测成功那么进行下一步骤;
2、在B-系中,XC被表示为矢量(cosα,sinα,0)Τ。如果XC以ZB为轴旋转了Δα,则新的矢量可以表示为(cos(α+Δα),sin(α+Δα),0)Τ。设β为ZC以XC为轴旋转的角度,那么:
设为ZC的单位向量,则ZC绕XC旋转角度β可以表示为具体的搜索方法表示如下:
(2a)设置初始点x(0)=(α0,β0)以及两个线性独立方向α1,β1,允许误差ε≥0且k=1;
(2b)设置x(k,0)=x(k-1)。从x(k,0)沿着αk方向使用线性搜索方法
结果为x(k,1)。同样通过沿着βk方向搜索可以得到xk点;
(2c)如果||ph(θ)-phmax||≤ε,则进行步骤3,其中phmax是星载X射线探测器记录的平均最大流量强度。否则,设置x0=xk,返回步骤(2a);
3、星载X射线探测器持续观测目标脉冲星,计算机计算得到平均流量强度。如果||ph(θ)-phmax||>ε,进行步骤(2)与(2a),并再一次调整指向。
如图4所示,以准直器3的通孔截面为圆形时为例,给出了一个准直器3单元的结构示意图,实际使用的准直器3是由多个这样的单元规则排列而成。其中,和分别为准直器3轴心指向方向矢量和脉冲星辐射方向矢量,θ为脉冲星辐射的X射线与准直器3轴心方向的夹角,H为准直器3的高度,R为准直器3截面圆的半径。只有当准直器3的轴心方向与脉冲星辐射方向一致时,全部的脉冲星信号能够无遮挡的被星载X射线探测器4接收,此时星载X射线探测器4接收到的光子流强度最大。当准直器轴心方向与脉冲星辐射方向有一个夹角为θ,则部分脉冲星信号能无遮挡的被星载X射线探测器4接收,此时星载X射线探测器4接收到的光子流量强度就是关于θ的函数。因此,当θ=0时,对应的星载X射线探测器4的指向就是脉冲星辐射方向矢量。
如图1所示,姿态矩阵确定13是利用脉冲星辐射方向矢量和地心方向矢量,在惯性坐标系内,以欧拉角表示的姿态变换矩阵,即
如图2所示,为本发明的基于地球质心参考坐标系中的脉冲星、航天器、太阳质心、地球和地平线张角之间的几何关系以及基于计时观测和仰角观测的导航算法原理图。地平线张角可以通过红外或紫外地平敏感仪测量得到,并且所观测的地心方向矢量和脉冲星矢量处于同一个平面内,如图2所示,为航天器相对于地心方向矢量,它与脉冲星辐射方向矢量的夹角为β,地平敏感仪敏感到的地平线张角为α,敏感仪的地平观测旋转α2角度可以得到地心方向矢量,地球半径为Re,则可以得到式(5)和式(6):
其中是星载X射线探测器对脉冲星的辐射矢量观测得到的。
如图2所示,图1中的导航算法库15的实现方法是建立在计时观测和仰角观测基础之上的。基于仰角观测的导航方法,其观测角为:
对其进行一阶泰勒展开,
其中ns包括泰勒展开的高阶项和噪声,,E[ns(t)]=0,同样对(6)式做一阶泰勒展开,
其中,nh包括泰勒展开的高阶项和地心矢量观测噪声,E[nh(t)]=0,
如图2所示,基于计时观测的脉冲星导航,在原子钟足够稳定的前提下,至少需要观测3颗脉冲星信号。为了实现单颗脉冲星导航,本发明使用脉冲星计时观测和仰角观测结合的方式,建立观测方程。计时观测的TOA观测方程为:
其中ωN(t)为第N颗脉冲星的TOA量测噪声。综合(8)式、(9)式、(10)式可建立观测方程为:
其中,H为观测矩阵,η(t)为零均值高斯观测噪声且:
η(t)=[ns(t),nh(t),ωN(t)]Τ,
本发明以地球卫星为例对所发明的系统和方法进行了阐述,而从理论上讲,本发明的利用X射线脉冲星的航天器姿态和位置的联合测量的几何方法完全适用于出地球外的其它近天体轨道航天器,仅需要将时间和坐标系统转换到相应的目标惯性系下即可。本专业领域人员在不背离本发明权利要求范围和主旨的前提下可以实现多种显而易见的改进,因此上述内容只是本发明借以阐述的实例,本发明的权利要求范围并不限于以上论述。
本说明未做详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (9)
1.基于X射线脉冲星的航天器姿态和位置测量系统,其特征在于,所述系统包括:星载X射线探测器、准直器、星载原子钟、X射线脉冲星特征参数数据库、X射线脉冲星辨识算法库、太阳系行星参数数据库、星载计算机、地平敏感仪,导航算法库、星载万向支架;其中
所述星载计算机根据当前所记录的航天器姿态和所述X射线脉冲星特征参数数据库,控制所述星载万向支架支撑星载X射线探测器指向目标脉冲星;
所述星载X射线探测器配合所述准直器,测量到达所述星载X射线探测器的X射线脉冲星信号的光子流量强度,从而计算得出所述星载X射线探测器指向与脉冲星辐射矢量的偏差,将该偏差反馈到所述星载计算机用于进行航天器姿态的控制增强,不断调整所述星载X射线探测器指向控制,以获得精确地脉冲星辐射方向矢量;
所述地平敏感仪敏感地平得到地平线张角,从而计算地心方向矢量,将该地心方向矢量和所述脉冲星辐射方向矢量共同送至所述星载计算机中的脉冲星定姿算法,同时将所述地平线张角和所述脉冲星辐射方向矢量共同送至脉冲星定位算法;
所述星载计算机调用所述导航算法库中的联合定位算法对定姿算法和定位算法进行结合,输出高精度的自主导航参数信息,包括姿态和位置。
2.根据权利要求1所述的航天器姿态和位置测量系统,其特征在于,所述星载X射线探测器与所述准直器配合用于接收脉冲星辐射的X射线,所述准直器安装于所述星载X射线探测器的前部;所述准直器由若干空心圆柱体互相紧密排列组合而成。
3.一种利用权利要求1所述的系统对航天器姿态和位置进行测量的方法,其特征在于,所述方法包括以下步骤:
(1)探测目标脉冲星:所述星载计算机根据当前所记录的航天器姿态和所述X射线脉冲星特征参数数据库,控制所述星载万向支架支撑星载X射线探测器指向目标脉冲星;
(2)调整所述星载X射线探测器方向:所述星载X射线探测器配合所述准直器,测量X射线脉冲星信号的光子流量强度,并计算得出所述星载X射线探测器指向与脉冲星辐射矢量的偏差;
(3)校正目标脉冲星位置:将步骤(2)获得的偏差反馈到所述星载计算机用于进行航天器姿态的控制增强,从而不断调整所述星载X射线探测器指向控制,以获得精确地脉冲星辐射方向矢量;
(4)确定航天器姿态与位置:所述地平敏感仪敏感地平得到地平线张角,从而计算地心方向矢量,将该地心方向矢量和所述脉冲星辐射方向矢量共同送至所述星载计算机中的脉冲星定姿算法,同时将所述地平线张角和所述脉冲星辐射方向矢量共同送至脉冲星定位算法;所述星载计算机调用所述导航算法库中的联合定位算法对定姿算法和定位算法进行结合,输出高精度的自主导航参数信息,包括姿态和位置。
4.根据权利要求3所述的对航天器姿态和位置进行测量的方法,其特征在于,所述步骤(2)中的偏差获取方法如下:首先定义两个坐标系,分别是B-系和C-系,B-系的坐标原点为航天器的质心,其X轴指向航天器的前进方向,Z轴垂直于偏航平面,Y轴根据右手定则确定,即(XB,YB,ZB);C-系(XC,YC,ZC)中,YC平行于支架的水平轴,ZC指向探测器的中心轴,XC垂直于YC和ZC的平面;明显地,所述星载X射线探测器指向控制矢量ZC与脉冲星辐射矢量的夹角θ是所要求的偏差,其中,C-系与B-系的坐标原点相同。
5.根据权利要求3或4所述的对航天器姿态和位置进行测量的方法,其特征在于,所述星载X射线探测器指向控制的实现方法如下:
根据三角关系,圆筒形准直器的有效探测面积为:
使得:
其中,θ为探测器指向与脉冲星辐射矢量的夹角,H为准直器的高,R为准直器横截圆面的半径。
6.根据权利要求3或4所述的对航天器姿态和位置进行测量的方法,其特征在于,所述脉冲星辐射方向矢量的搜索方法如下:所述星载计算机控制所述星载万向支架的方向,同时控制所述星载X射线探测器指向,所述星载X射线探测器记录敏感到的脉冲星辐射光子流量,通过改变其的指向ZC来减小θ,当光子流量强度最大时,即可认为所述星载X射线探测器指向与脉冲星辐射方向矢量平行;所述星载万向支架的水平轴XC只能在(XC,YC)平面内以ZB为轴进行旋转,α用来表示旋转角度;由于ZC与XC是垂直的,那么ZC以XC为轴进行旋转时,两个轴在线性独立方向的旋转角度都是α,同时要求脉冲星光子流量强度作为目标函数是线性递减的;其中包括:
(1)根据姿态信息的先验知识,所述星载X射线探测器根据参考星历对目标脉冲星在B-系中的位置进行粗估计,所述星载X射线探测器旋转中心轴使其指向脉冲星以获取脉冲星信号;同时所述星载计算机使用FFT或是贝叶斯定理通过捕获到的光子序列Ph(θ)来检测目标信号,检测成功后进行步骤(2);
(2)在B-系中,XC被表示为矢量(cosα,sinα,0)Τ,如果XC以ZB为轴旋转了Δα,则新的矢量表示为(cos(α+Δα),sin(α+Δα),0)Τ;设β为ZC以XC为轴旋转的角度,那么
设为ZC的单位向量,则ZC绕XC旋转角度β表示为具体是:
(2a)设置初始点x(0)=(α0,β0)以及两个线性独立方向α1,β1,允许误差ε≥0且k=1;
(2b)设置x(k,0)=x(k-1),从x(k,0)沿着αk方向使用线性搜索方法:
结果为x(k,1),同样通过沿着βk方向搜索可以得到xk点;
(2c)如果||ph(θ)-phmax||≤ε,则进行步骤(3),其中phmax是所述星载X射线探测器记录的平均最大流量强度;否则,设置x0=xk,返回步骤(2a);
(3)所述星载X射线探测器持续观测目标脉冲星,星载计算机计算得到平均流量强度;如果||ph(θ)-phmax||>ε,进行步骤(2)与(2a),并再一次调整指向。
7.根据权利要求3所述的对航天器姿态和位置进行测量的方法,其特征在于,所述航天器确定位置实现方法如下:所述星载计算机接收所述星载X射线探测器输出的脉冲星辐射方向矢量和所述地平敏感仪输出的地平线张角和地心方向矢量,建立两个关系式A和B;其中关系式A为地平线张角与卫星位置矢量的函数,即关系式B为地心方向矢量与脉冲星方向矢量之间的夹角以及卫星位置矢量的函数,即同时,利用脉冲星计时观测通过累计轮廓与标准相位轮廓比相而得到相差,从而计算距离信息;将该距离信息与前述的关系式A和B一起送入导航滤波算法库,完成定位。
8.根据权利要求3所述的对航天器姿态和位置进行测量的方法,其特征在于,所述导航滤波算法实现方法为:设地平敏感仪观测到的地平线张角为α,卫星地心方向矢量与脉冲星辐射方向矢量的夹角为β,建立观测方程:
α-α=FsδX+ns(t)=Δα (1)
β-β=FhδX+nh(t)=Δβ (2)
其中ns包括泰勒展开的高阶项和噪声,nh包括泰勒展开的高阶项和地心矢量观测噪声,Ens(t)]=0,E[nh(t)]=0,同时脉冲星计时观测建立的观测方程为:
Δτ=ΔtTOA-ΔtTOA=GNδX+ωN(t) (3)
综合(1)、(2)和(3)式建立组合观测方程
其中,H为观测矩阵,η(t)为零均值高斯观测噪声且η(t)=[ns(t),nh(t),ωN(t)]Τ,其中ωN(t)为第N颗脉冲星的TOA量测噪声。
9.一种利用权利要求3所述对航天器姿态和位置进行测量的方法,适用于地球外的其它近天体轨道航天器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310320673.0A CN103389099B (zh) | 2013-07-26 | 2013-07-26 | 基于x射线脉冲星的航天器姿态和位置测量系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310320673.0A CN103389099B (zh) | 2013-07-26 | 2013-07-26 | 基于x射线脉冲星的航天器姿态和位置测量系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103389099A CN103389099A (zh) | 2013-11-13 |
CN103389099B true CN103389099B (zh) | 2016-12-28 |
Family
ID=49533444
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310320673.0A Expired - Fee Related CN103389099B (zh) | 2013-07-26 | 2013-07-26 | 基于x射线脉冲星的航天器姿态和位置测量系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103389099B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104567879B (zh) * | 2015-01-27 | 2018-08-21 | 北京控制工程研究所 | 一种组合视场导航敏感器地心方向提取方法 |
CN104864875B (zh) * | 2015-04-03 | 2018-01-05 | 北京控制工程研究所 | 一种基于非线性h∞滤波的航天器自主定位方法 |
CN105224731B (zh) * | 2015-09-17 | 2018-04-27 | 南京大学 | 地球静止卫星紫外成像敏感器的辐射模拟仿真方法 |
CN105424048B (zh) * | 2015-11-03 | 2018-12-21 | 北京空间飞行器总体设计部 | 一种基于周期变星的航天器自主导航方法 |
CN105571597B (zh) * | 2015-12-16 | 2019-02-15 | 中国空间技术研究院 | 一种超深空x射线主动导航系统 |
CN108020242B (zh) * | 2017-12-06 | 2019-12-13 | 中国人民解放军国防科技大学 | 脉冲星探测器指向误差在轨标定方法、处理器及存储介质 |
CN109631877B (zh) * | 2018-12-05 | 2020-06-09 | 北京控制工程研究所 | 一种指向驻留的螺旋分段式脉冲星搜索方法 |
CN109683208B (zh) * | 2019-01-25 | 2021-02-09 | 北京空间飞行器总体设计部 | 一种适应空间x射线源定位精度分析方法 |
CN110186464B (zh) * | 2019-05-30 | 2023-03-28 | 西安电子科技大学 | 一种基于贝叶斯估计的x射线脉冲星导航toa估计方法 |
CN111142358B (zh) * | 2020-01-20 | 2022-01-07 | 山东科技大学 | 一种基于自然x射线源观测的深空自主守时方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102168981A (zh) * | 2011-01-13 | 2011-08-31 | 北京航空航天大学 | 一种深空探测器火星捕获段自主天文导航方法 |
CN102243311A (zh) * | 2011-04-15 | 2011-11-16 | 中国人民解放军国防科学技术大学 | 一种x射线脉冲星导航使用的选星方法 |
CN102829787A (zh) * | 2012-08-30 | 2012-12-19 | 西安电子科技大学 | X射线脉冲星成像仪及其编码成像方法 |
CN103048000A (zh) * | 2012-12-29 | 2013-04-17 | 中国空间技术研究院 | X射线脉冲星导航地面试验系统 |
CN103105811A (zh) * | 2012-12-29 | 2013-05-15 | 中国空间技术研究院 | X射线脉冲星导航地面试验信号控制系统 |
CN103217162A (zh) * | 2013-03-21 | 2013-07-24 | 西安电子科技大学 | 采用稀疏表示的脉冲星累积脉冲轮廓时间延迟测量方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101285359B1 (ko) * | 2006-11-14 | 2013-07-11 | 삼성전자주식회사 | 동작 추적 장치 및 방법 |
-
2013
- 2013-07-26 CN CN201310320673.0A patent/CN103389099B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102168981A (zh) * | 2011-01-13 | 2011-08-31 | 北京航空航天大学 | 一种深空探测器火星捕获段自主天文导航方法 |
CN102243311A (zh) * | 2011-04-15 | 2011-11-16 | 中国人民解放军国防科学技术大学 | 一种x射线脉冲星导航使用的选星方法 |
CN102829787A (zh) * | 2012-08-30 | 2012-12-19 | 西安电子科技大学 | X射线脉冲星成像仪及其编码成像方法 |
CN103048000A (zh) * | 2012-12-29 | 2013-04-17 | 中国空间技术研究院 | X射线脉冲星导航地面试验系统 |
CN103105811A (zh) * | 2012-12-29 | 2013-05-15 | 中国空间技术研究院 | X射线脉冲星导航地面试验信号控制系统 |
CN103217162A (zh) * | 2013-03-21 | 2013-07-24 | 西安电子科技大学 | 采用稀疏表示的脉冲星累积脉冲轮廓时间延迟测量方法 |
Non-Patent Citations (2)
Title |
---|
《X射线脉冲星信号检测及位相测量技术研究》;张华;《万方学位论文》;20120731;全文 * |
《基于光量子探测的XPNAV半物理仿真及建模》;张华等;《光电子 激光》;20110630;第22卷(第6期);905-910 * |
Also Published As
Publication number | Publication date |
---|---|
CN103389099A (zh) | 2013-11-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103389099B (zh) | 基于x射线脉冲星的航天器姿态和位置测量系统及方法 | |
CN102168981B (zh) | 一种深空探测器火星捕获段自主天文导航方法 | |
CN104332707B (zh) | 一种用于低轨星载天线跟踪地面站的方法 | |
CN102175241B (zh) | 一种火星探测器巡航段自主天文导航方法 | |
CN103674032B (zh) | 融合脉冲星辐射矢量和计时观测的卫星自主导航系统及方法 | |
CN104567819B (zh) | 一种星载相机全视场偏流角确定与补偿方法 | |
CN103983254A (zh) | 一种新型敏捷卫星机动中成像方法 | |
CN103744052B (zh) | 一种应用于空中目标定位的双星测时差测向方法及装置 | |
CN107380485B (zh) | 一种微小卫星大面阵广域多模凝视成像控制方法 | |
CN102116628A (zh) | 一种着陆或附着深空天体探测器的高精度导航方法 | |
CN106595674A (zh) | 基于星敏感器和星间链路的heo卫星编队飞行自主导航方法 | |
CN104764449B (zh) | 一种基于星历修正的捕获段深空探测器自主天文导航方法 | |
CN105184002A (zh) | 一种数传天线指向角度的仿真分析方法 | |
CN101858747A (zh) | 一种有效利用地球辐照能的卫星帆板对日定向目标姿态的解析确定方法 | |
CN102514734B (zh) | 基于日地系统Halo轨道探测器构型与姿态指向的姿态递推方法 | |
CN102944238B (zh) | 一种行星探测器接近目标过程中相对位置确定方法 | |
Wang et al. | Absolute navigation for Mars final approach using relative measurements of X-ray pulsars and Mars orbiter | |
CN111427003A (zh) | 地面测站天线对卫星的指向导引系统 | |
CN111427001A (zh) | 适用于地面测站天线对卫星指向的目标定位方法 | |
CN106289156A (zh) | 一种卫星以任意姿态成像时获取摄影点太阳高度角的方法 | |
Yang et al. | A real-time calibration method for the systematic errors of a star sensor and gyroscope units based on the payload multiplexed | |
CN103950555A (zh) | 超近距离的高精度相对位置保持控制方法 | |
Cui et al. | Real-time navigation for Mars final approach using X-ray pulsars | |
CN102519454B (zh) | 一种日地月导航的月心方向修正方法 | |
Xu et al. | A novel X-ray pulsar integrated navigation method for ballistic aircraft |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20161228 |