CN110146093B - 双体小行星探测自主协同光学导航方法 - Google Patents
双体小行星探测自主协同光学导航方法 Download PDFInfo
- Publication number
- CN110146093B CN110146093B CN201910529172.0A CN201910529172A CN110146093B CN 110146093 B CN110146093 B CN 110146093B CN 201910529172 A CN201910529172 A CN 201910529172A CN 110146093 B CN110146093 B CN 110146093B
- Authority
- CN
- China
- Prior art keywords
- star
- axis
- detectors
- double
- detector
- 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.)
- Active
Links
Images
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/20—Instruments for performing navigational calculations
-
- 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/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
Abstract
本发明公开的双体小行星探测自主协同光学导航方法,属于自主导航领域。本发明实现方法为:当主星的位置已知,从星的位置不确知,在双小行星系统惯性坐标系下建立探测器的动力学方程;根据探测器所使用的导航系统建立主星和从星的视线观测模型和探测器间的测距测角模型,并设计双体小行星导航的无迹卡尔曼滤波器;以轨道预报给出两颗探测器位置速度和从星位置的先验估计值作为初始状态,结合动力学方程和观测方程,通过无迹卡尔曼滤波器进行估计,实时确定两颗探测器的位置速度和从星的位置。本发明具有实时性好、导航精度高的优点。本发明利用无迹卡尔曼滤波器进行估计,不需要将非线性系统线性化,也不需要计算雅可比矩阵,实现简单、收敛性好。
Description
技术领域
本发明涉及一种双体小行星探测自主协同光学导航方法,属于自主导航领域。
背景技术
小行星探测是当前深空探测领域的研究热点。与单小行星相比,双小行星系统的不规则引力场、自旋以及双星相互作用等因素,导致双星系统动力学更加复杂、探测难度大,对导航系统要求高。自主光学导航是实现双体小行星探测的关键技术之一,为实现双体小行星精准探测,需要研究双体小行星探测器自主精确导航技术,是当前科技人员关注的重点问题之一。
在已发展的双体小行星导航方法中,在先技术[1](参见PeytavíG G,Andert T P,Probst A,et al.3D Imaging for Autonomous Navigation About Small-Bodies[C].AIAA/AAS Astrodynamics Specialist Conference.Long Beach,California,2016),利用三维激光雷达对双小行星系统的主星进行扫描,同时结合扩展卡尔曼滤波进行状态估计,由于双小行星系统引力场建模的不准确和扩展卡尔曼滤波本身的缺陷,导航估计结果发散,无法有效地确定探测器的状态。
在先技术[2](参见Torre F,Vasile M,Serra R,et al.Autonomous Navigationof a Formation of Spacecraft in the Proximity of a Binary Asteroid[C].International Symposium on Space Technology and Science,Ehime,Japan,2017),利用地面站、激光雷达、相机和探测器间的通信获得观测信息,并采用无迹H∞滤波方法进行状态估计,对单颗或多颗探测器在双小行星系统拉格朗日点L4周围的导航和位置保持进行研究。由于导航过程中使用了地面站的观测信息,导航的自主性和实时性不能双体小行星探测的需求。
在先技术[3](参见Vasile M,Torre F,Serra R,et al.Angles-Only Navigationof a Formation in the Proximity of a Binary System[C].2018Space FlightMechanics Meeting,Kissimmee,Florida,2018.),仅利用相机和探测器间的通信获得观测信息,并采用无迹H∞滤波方法进行状态估计,对单颗/两颗探测器在双小行星系统拉格朗日点L4周围的导航和位置保持进行研究。由于该方法中的单颗/两颗探测器导航均是在主星和从星位置已知的情况下进行的,而实际过程中主星和从星星历并不准确,从星相对于主星的位置也因此不准确,导航精度,需要对从星相对于主星的位置进行估计。
发明内容
为解决双体小行星导航过程中,主星和从星星历不准确使从星相对于主星的位置不准确并影响导航精度的问题,本发明公开的双体小行星探测自主协同光学导航方法要解决的技术问题是:利用两颗探测器分别对主星和从星进行视线观测,同时结合探测器间的测距测角,能够确定两颗探测器的位置速度和从星相对于主星的位置,并且具有导航精度高的优点。
本发明目的是通过下述技术方案实现。
本发明公开的双体小行星探测自主协同光学导航方法,当主星的位置已知,从星的位置不确知,在双小行星系统惯性坐标系下建立探测器的动力学方程;根据探测器所使用的导航系统建立主星和从星的视线观测模型和探测器间的测距测角模型,并设计双体小行星导航的无迹卡尔曼滤波器;以轨道预报给出两颗探测器位置速度和从星位置的先验估计值作为初始状态,结合动力学方程和观测方程,通过无迹卡尔曼滤波器进行估计,实时确定两颗探测器的位置速度和从星的位置。
本发明公开的双体小行星探测自主协同光学导航方法,包括步骤如下:
步骤1,在双小行星系统惯性坐标系下建立探测器的动力学方程。
选择双小行星系统的质心作为原点建立双小行星系统惯性坐标系,选择X轴指向春分点,Z轴为沿双小行星互相旋转的角速度方向,Y轴与X轴、Z轴垂直构成右手坐标系。选择从星质心作为原点建立从星固连坐标系,x轴沿从星最小惯量轴方向,z轴沿从星的自旋轴方向,y轴与x轴、z轴垂直构成右手坐标系。
定义rA、rB、rS和rSc代表主星、从星、太阳和探测器分别在双星系统惯性系中的位置,定义探测器相对于主星的位置矢量rASc=rSc-rA,探测器相对于从星的位置矢量rBSc=rSc-rB,探测器相对于太阳的位置矢量rSSc=rSc-rS,则探测器在双小行星系统惯性坐标系下的动力学方程表示为
在动力学方程中,μA为主星的引力常数,μB为从星的引力常数,μS为太阳的引力常数,IRB为由从星固连坐标系到双星系统惯性坐标系的转换矩阵,有
ωB为从星绕Z轴旋转的角速度,t表示时间。
在动力学方程中(1),U20,22为从星的引力场,有
在动力学方程(1)中,aSRP为太阳辐射压力引起的加速度,有
A和mSc分别是探测器的横截面积和质量,CR是反射系数,AU是1天文单位,ρS是1AU的太阳辐射压力。
步骤2,根据光学导航系统建立探测器的观测方程。
在双体小行星探测自主协同光学导航中,探测器通过导航相机和星间测量获取观测信息。选择相机的光心作为原点建立导航相机固连坐标系,ZC轴为相机的光轴,垂直于图像平面,XC轴、YC轴分别与图像平面两边平行,构成右手坐标系。选择探测器质心作为原点建立探测器本体固连坐标系,Zb轴沿探测器的最大惯量主轴,Xb轴、Yb轴指向其它两个惯量主轴,构成右手坐标系。
对于导航相机观测信息,小天体形心的位置在导航相机坐标系中的坐标为(XC,YC,ZC)T,对应的成像点像素为(u,v)T,导航相机的焦距为f。根据小孔成像原理,两点之间的关系如下:
则相机坐标系下小行星形心的单位指向矢量公式为:
则探测器观测主星和从星的视线信息分别为
其中,VCA为主星视线信息,VCB为从星视线信息,主星形心像素为(uA,vA)T,从星形心像素为(uB,vB)T。
对于星间测量,探测器采用激光测距获取两颗探测器间的相对距离,采用视觉导航系统获取两颗探测器间的角度信息。在其中任意一颗探测器的本体系下,两颗探测器间的相对位置为rr=(xr,yr,zr)T,相对距离为dr,视线方位角为俯仰角为ψr,则dr、ψr的表达式分别为:
星间测量由探测器间的相对距离、方位角和俯仰角组成,其表达式如下:
步骤3,设计双体小行星导航的无迹卡尔曼滤波器。
对于双体小行星探测自主协同光学导航,两颗探测器分别使用相机对主星和从星进行观测,同时通过星间测量获取探测间的相对位置dr和角度信息ψr,通过滤波器估计两颗探测器的位置速度和从星的位置。系统的状态方程和观测方程表示为
y=h(x,v)=[VCA1 VCB1 VCA2 VCB2 Zr]T+v
其中rsc1和rsc2分别为两颗探测器的位置,vsc1和vsc2分别为两颗探测器的速度,VCA1和VCA2分别为两颗探测器观测主星的视线信息,VCB1和VCB2分别为两颗探测器观测从星的视线信息。随机变量u、v分别代表过程噪声和观测噪声,u、v均为不相关零均值高斯白噪声序列,均服从分布u~N(0,Q)、v~N(0,R),Q为过程噪声协方差阵,R为观测误差方差阵。
其中χ矩阵由2l+1个向量组成,l为状态向量的维数,控制sigma点分布的参数λ=α2(l+κ)-l,其中κ=3-l,1e-4≤α<1,下标i代表矩阵的第i列。利用状态方程对sigma点进行时间更新,计算加权均值和协方差:
其中β用于减小高阶项误差,可根据状态x分布的先验知识确定,对于高斯分布,β=2最优。利用观测方程和sigma点进行测量更新,计算加权均值和协方差:
计算无迹卡尔曼滤波增益,更新状态向量和协方差矩阵,实现无迹卡尔曼滤波估计:
步骤4,以轨道预报给出两颗探测器位置速度和从星位置的先验估计值作为初始状态,结合动力学方程和观测方程,通过无迹卡尔曼滤波器进行估计,实时确定两颗探测器的位置速度和从星的位置,至此,实现双体小行星探测自主协同光学导航。
有益效果:
1、本发明公开的双体小行星探测自主协同光学导航方法,以轨道预报给出两颗探测器位置速度和从星位置的先验估计值作为初始状态,结合每颗探测器动力学方程和两颗探测器对主星、从星的视线观测和星间测量,通过无迹卡尔曼滤波器进行估计,能够在从星相对于主星位置未知的情况下实时确定两颗探测器的位置速度和从星的位置,即实现双体小行星探测自主协同光学导航,具有实时性好、导航精度高的优点。
2、本发明公开的双体小行星探测自主协同光学导航方法,利用无迹卡尔曼滤波器进行估计,不需要将非线性系统线性化,也不需要计算雅可比矩阵,具有实现简单、收敛性好、且滤波精度高的优点。
附图说明
图1为双体小行星探测自主协同光学导航的示意图;
图2为双体小行星探测自主协同光学导航的流程图;
图3为本发明步骤1中的双小行星系统惯性系下小行星与探测器的位置关系;
图4为本发明步骤2中的导航相机观测模型;
图5为本发明步骤2中的星间测量观测模型;
图6为本发明实例中的双体小行星探测自主协同光学导航的状态误差,其中图6a)为探测器1的位置估计误差曲线,6b)为探测器1的速度估计误差曲线,6c)为探测器2的位置估计误差曲线,6d)为探测器2的速度估计误差曲线,6e)为从星的位置估计误差曲线。
具体实施方式
为了更好地说明本发明的目的与优点,下面结合附图和实例对本发明做进一步说明。
实施例1:
选择双小行星系统65803Didymos作为目标天体,表1列出了双小行星系统65803Didymos的各项参数。令双体小行星绕双小行星系统中心旋转的轨道均与双小行星系统绕太阳旋转的轨道在同一平面内,导航相机分辨率为500×500像素,视场角为90°,导航相机焦距为20mm,相机测量误差为4×10-3rad,星间测距误差为0.1m,星间测角误差为1×10-3rad。探测器质量恒定为500kg,最大横截面面积为20m2,反射系数假定为1。
表1双体小行星65803Didymos的参数
令探测器1的初始位置为双小行星系统拉格朗日点L4,轨道法向量为h1=[0 01]T,则探测器1的初始速度为轨道法向量与初始位置向量积;令探测器2在双小行星系统惯性坐标系下的初始位置为[1000 1000 100]Tm,轨道法向量为h2=[0.9914 -0.01300.1305]T,则探测器2的初始速度为轨道法向量与初始位置向量积。仿真时间与双小行星系统旋转周期相同,约为11.92小时,在从星相对于主星位置未知的情况下,令初始状态估计误差均为实际初始状态的20%。
本实施例公开的双体小行星探测自主协同光学导航方法,具体实现步骤如下:
步骤1,在双小行星系统惯性坐标系下建立探测器的动力学方程。
选择双小行星系统的质心作为原点建立双小行星系统惯性坐标系,选择X轴指向春分点,Z轴为沿双小行星互相旋转的角速度方向,Y轴与X轴、Z轴垂直构成右手坐标系。选择从星质心作为原点建立从星固连坐标系,x轴沿从星最小惯量轴方向,z轴沿从星的自旋轴方向,y轴与x轴、z轴垂直构成右手坐标系。
定义rA、rB、rS和rSc代表主星、从星、太阳和探测器分别在双星系统惯性系中的位置,定义探测器相对于主星的位置矢量rASc=rSc-rA,探测器相对于从星的位置矢量rBSc=rSc-rB,探测器相对于太阳的位置矢量rSSc=rSc-rS,如图3所示,则探测器在双小行星系统惯性坐标系下的动力学方程表示为
在动力学方程中,μA为主星的引力常数,μB为从星的引力常数,μS为太阳的引力常数,IRB为由从星固连坐标系到双星系统惯性坐标系的转换矩阵,有
ωB为从星绕Z轴旋转的角速度,t表示时间。
在动力学方程中(1),U20,22为从星的引力场,有
在动力学方程(1)中,aSRP为太阳辐射压力引起的加速度,有
A和mSc分别是探测器的横截面积和质量,CR是反射系数,AU是1天文单位,ρS是1AU的太阳辐射压力,ρS=4.5605×10-6N/m2。
步骤2,根据光学导航系统建立探测器的观测方程。
在双体小行星探测自主协同光学导航中,探测器通过导航相机和星间测量获取观测信息。选择相机的光心作为原点建立导航相机固连坐标系,ZC轴为相机的光轴,垂直于图像平面,XC轴、YC轴分别与图像平面两边平行,构成右手坐标系。选择探测器质心作为原点建立探测器本体固连坐标系,Zb轴沿探测器的最大惯量主轴,Xb轴、Yb轴指向其它两个惯量主轴,构成右手坐标系。
对于导航相机观测信息,小天体形心的位置在导航相机坐标系中的坐标为(XC,YC,ZC)T,对应的成像点像素为(u,v)T,导航相机的焦距为f。根据小孔成像原理,两点之间的关系如下:
则相机坐标系下小行星形心的单位指向矢量公式为:
则如图4所示,探测器观测主星和从星的视线信息分别为
其中,VCA为主星视线信息,VCB为从星视线信息,主星形心像素为(uA,vA)T,从星形心像素为(uB,vB)T。对于星间测量,探测器采用激光测距获取两颗探测器间的相对距离,采用视觉导航系统获取两颗探测器间的角度信息。如图5所示,在其中任意一颗探测器的本体系下,两颗探测器间的相对位置为rr=(xr,yr,zr)T,相对距离为dr,视线方位角为俯仰角为ψr,则dr、ψr的表达式分别为:
星间测量由探测器间的相对距离、方位角和俯仰角组成,其表达式如下:
步骤3,设计双体小行星导航的无迹卡尔曼滤波器。
对于双体小行星探测自主协同光学导航,两颗探测器分别使用相机对主星和从星进行观测,同时通过星间测量获取探测间的相对位置dr和角度信息ψr,通过滤波器估计两颗探测器的位置速度和从星的位置。系统的状态方程和观测方程表示为
y=h(x,v)=[VCA1 VCB1 VCA2 VCB2 Zr]T+v
其中rsc1和rsc2分别为两颗探测器的位置,vsc1和vsc2分别为两颗探测器的速度,VCA1和VCA2分别为两颗探测器观测主星的视线信息,VCB1和VCB2分别为两颗探测器观测从星的视线信息。随机变量u、v分别代表过程噪声和观测噪声,u、v均为不相关零均值高斯白噪声序列,均服从分布u~N(0,Q)、v~N(0,R),Q为过程噪声协方差阵,R为观测误差方差阵。
其中χ矩阵由2l+1个向量组成,l为状态向量的维数,控制sigma点分布的参数λ=α2(l+κ)-l,其中κ=3-l,α=0.001,下标i代表矩阵的第i列。利用状态方程对sigma点进行时间更新,计算加权均值和协方差:
其中β用于减小高阶项误差,可根据状态x分布的先验知识确定,对于高斯分布,β=2最优。利用观测方程和sigma点进行测量更新,计算加权均值和协方差:
计算无迹卡尔曼滤波增益,更新状态向量和协方差矩阵,实现无迹卡尔曼滤波估计:
步骤4,以轨道预报给出两颗探测器位置速度和从星位置的先验估计值作为初始状态,结合动力学方程和观测方程,通过无迹卡尔曼滤波器进行估计,实时确定两颗探测器的位置速度和从星的位置,至此,实现双体小行星探测自主协同光学导航。
利用本发明方法进行导航滤波估计,其效果如图6所示,由图6可知,两颗探测器的位置确定精度最终均达到10m,速度确定精度均达到0.25m/s,从星的位置确定精度达到2m,在从星相对于主星位置未知的情况下,本发明方法能够准确地确定两颗探测器的位置速度和从星的位置。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.双体小行星探测自主协同光学导航方法,其特征在于:包括步骤如下,
步骤1,在双小行星系统惯性坐标系下建立探测器的动力学方程;
步骤2,根据光学导航系统建立探测器的观测方程;
步骤3,设计双体小行星导航的无迹卡尔曼滤波器;
步骤4,以轨道预报给出两颗探测器位置速度和从星位置的先验估计值作为初始状态,结合动力学方程和观测方程,通过无迹卡尔曼滤波器进行估计,实时确定两颗探测器的位置速度和从星的位置,至此,实现双体小行星探测自主协同光学导航;
其中,步骤1实现方法为,
选择双小行星系统的质心作为原点建立双小行星系统惯性坐标系,选择X轴指向春分点,Z轴为沿双小行星互相旋转的角速度方向,Y轴与X轴、Z轴垂直构成右手坐标系;选择从星质心作为原点建立从星固连坐标系,x轴沿从星最小惯量轴方向,z轴沿从星的自旋轴方向,y轴与x轴、z轴垂直构成右手坐标系;
定义rA、rB、rS和rSc代表主星、从星、太阳和探测器分别在双星系统惯性系中的位置,定义探测器相对于主星的位置矢量rASc=rSc-rA,探测器相对于从星的位置矢量rBSc=rSc-rB,探测器相对于太阳的位置矢量rSSc=rSc-rS,则探测器在双小行星系统惯性坐标系下的动力学方程表示为
在动力学方程中,μA为主星的引力常数,μB为从星的引力常数,μS为太阳的引力常数,IRB为由从星固连坐标系到双星系统惯性坐标系的转换矩阵,有
ωB为从星绕Z轴旋转的角速度,t表示时间;
在动力学方程中(1),U20,22为从星的引力场,有
在动力学方程(1)中,aSRP为太阳辐射压力引起的加速度,有
A和mSc分别是探测器的横截面积和质量,CR是反射系数,AU是1天文单位,ρS是1AU的太阳辐射压力;
步骤2实现方法为,
在双体小行星探测自主协同光学导航中,探测器通过导航相机和星间测量获取观测信息;选择相机的光心作为原点建立导航相机固连坐标系,ZC轴为相机的光轴,垂直于图像平面,XC轴、YC轴分别与图像平面两边平行,构成右手坐标系;选择探测器质心作为原点建立探测器本体固连坐标系,Zb轴沿探测器的最大惯量主轴,Xb轴、Yb轴指向其它两个惯量主轴,构成右手坐标系;
对于导航相机观测信息,小天体形心的位置在导航相机坐标系中的坐标为(XC,YC,ZC)T,对应的成像点像素为(u,v)T,导航相机的焦距为f;根据小孔成像原理,两点之间的关系如下:
则相机坐标系下小行星形心的单位指向矢量公式为:
则探测器观测主星和从星的视线信息分别为
其中,VCA为主星视线信息,VCB为从星视线信息,主星形心像素为(uA,vA)T,从星形心像素为(uB,vB)T;
对于星间测量,探测器采用激光测距获取两颗探测器间的相对距离,采用视觉导航系统获取两颗探测器间的角度信息;在其中任意一颗探测器的本体系下,两颗探测器间的相对位置为rr=(xr,yr,zr)T,相对距离为dr,视线方位角为俯仰角为ψr,则dr、ψr的表达式分别为:
星间测量由探测器间的相对距离、方位角和俯仰角组成,其表达式如下:
步骤3实现方法为,
对于双体小行星探测自主协同光学导航,两颗探测器分别使用相机对主星和从星进行观测,同时通过星间测量获取探测间的相对位置dr和角度信息ψr,通过滤波器估计两颗探测器的位置速度和从星的位置;系统的状态方程和观测方程表示为
y=h(x,v)=[VCA1 VCB1 VCA2 VCB2 Zr]T+v (12)
其中rsc1和rsc2分别为两颗探测器的位置,vsc1和vsc2分别为两颗探测器的速度,VCA1和VCA2分别为两颗探测器观测主星的视线信息,VCB1和VCB2分别为两颗探测器观测从星的视线信息;随机变量u、v分别代表过程噪声和观测噪声,u、v均为不相关零均值高斯白噪声序列,均服从分布u:N(0,Q)、v:N(0,R),Q为过程噪声协方差阵,R为观测误差方差阵;
其中χ矩阵由2l+1个向量组成,l为状态向量的维数,控制sigma点分布的参数λ=α2(l+κ)-l,其中κ=3-l,1e-4≤α<1,下标i代表矩阵的第i列;利用状态方程对sigma点进行时间更新,计算加权均值和协方差:
其中β用于减小高阶项误差,可根据状态x分布的先验知识确定,对于高斯分布,β=2最优;
利用观测方程和sigma点进行测量更新,计算加权均值和协方差:
计算无迹卡尔曼滤波增益,更新状态向量和协方差矩阵,实现无迹卡尔曼滤波估计:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910529172.0A CN110146093B (zh) | 2019-06-19 | 2019-06-19 | 双体小行星探测自主协同光学导航方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910529172.0A CN110146093B (zh) | 2019-06-19 | 2019-06-19 | 双体小行星探测自主协同光学导航方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110146093A CN110146093A (zh) | 2019-08-20 |
CN110146093B true CN110146093B (zh) | 2020-12-15 |
Family
ID=67595761
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910529172.0A Active CN110146093B (zh) | 2019-06-19 | 2019-06-19 | 双体小行星探测自主协同光学导航方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110146093B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110716557B (zh) * | 2019-11-18 | 2021-05-11 | 河北工业大学 | 基于先验动力学知识的机器人参数识别及接触力监测方法 |
CN110672105B (zh) * | 2019-11-22 | 2021-04-20 | 北京理工大学 | 一种小天体接近段双探测器高精度协同光学导航方法 |
CN110889219A (zh) * | 2019-11-22 | 2020-03-17 | 北京理工大学 | 基于器间测距的小天体引力场反演修正方法 |
CN111238484B (zh) * | 2020-02-28 | 2022-04-12 | 上海航天控制技术研究所 | 一种基于球形无迹变换的环火轨道自主导航方法 |
CN111735459B (zh) * | 2020-07-15 | 2022-03-22 | 北京理工大学 | 小天体探测器间协同导航方法 |
CN111680462B (zh) * | 2020-08-11 | 2020-11-10 | 北京控制与电子技术研究所 | 基于空间目标在光学相平面位置变化的制导方法和系统 |
CN115065410B (zh) * | 2022-07-06 | 2023-09-15 | 北京理工大学 | 卫星间高精度激光链路扫描捕获方法 |
CN115877370B (zh) * | 2023-03-08 | 2023-07-07 | 中国西安卫星测控中心 | 一种利用双雷达距离与方位角快速计算航天器轨道的方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB758933A (en) * | 1953-10-01 | 1956-10-10 | Decca Record Co Ltd | Improvements in or relating to navigational computers |
CN105652308A (zh) * | 2014-11-27 | 2016-06-08 | 航天恒星科技有限公司 | 飞行器相对测量方法及系统 |
CN104406583B (zh) * | 2014-12-05 | 2017-02-22 | 中国人民解放军63680部队 | 双星敏感器联合确定载体姿态方法 |
CN104457705B (zh) * | 2014-12-12 | 2016-10-26 | 北京理工大学 | 基于天基自主光学观测的深空目标天体初定轨方法 |
CN106767768A (zh) * | 2016-12-06 | 2017-05-31 | 东南大学 | 一种双星编队的自主导航方法 |
CN107144283A (zh) * | 2017-06-30 | 2017-09-08 | 上海航天控制技术研究所 | 一种用于深空探测器的高可观度光学脉冲星混合导航方法 |
CN108871349B (zh) * | 2018-07-13 | 2021-06-15 | 北京理工大学 | 一种深空探测器光学导航位姿加权确定方法 |
-
2019
- 2019-06-19 CN CN201910529172.0A patent/CN110146093B/zh active Active
Non-Patent Citations (2)
Title |
---|
Control of Spacecraft Trajectories near Collinear Equilibrium Points of Binary Asteroid Systems;Woo,Pamela,Misra,Arun K;《Journal of Guidance Control & Dynamics》;20151231;全文 * |
双小行星系统引力场建模方法研究;步士超等;《深空探测学报》;20170430;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110146093A (zh) | 2019-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110146093B (zh) | 双体小行星探测自主协同光学导航方法 | |
Trebi-Ollennu et al. | Design and analysis of a sun sensor for planetary rover absolute heading detection | |
CN106595674B (zh) | 基于星敏感器和星间链路的heo卫星编队飞行自主导航方法 | |
WO2017113567A1 (zh) | 火星探测器自主导航方法 | |
CN105203101B (zh) | 一种基于目标天体星历修正的深空探测器捕获段天文导航方法 | |
US6622970B2 (en) | Method and apparatus for autonomous solar navigation | |
CN100533065C (zh) | 基于多天体路标的星际巡航自主导航方法 | |
CN103674032A (zh) | 融合脉冲星辐射矢量和计时观测的卫星自主导航系统及方法 | |
CN110412868B (zh) | 一种使用星间光学图像的非合作航天器轨道确定方法 | |
CN109709537A (zh) | 一种基于卫星编队的非合作目标位置速度跟踪方法 | |
CN111102981B (zh) | 一种基于ukf的高精度卫星相对导航方法 | |
CN110146092B (zh) | 基于导航信息评价的双体小行星探测轨迹优化方法 | |
CN111366984B (zh) | 一种基于重力卫星星间激光测距系统确定引力场模型的方法 | |
Liu et al. | X-ray pulsar/starlight Doppler integrated navigation for formation flight with ephemerides errors | |
CN109269510A (zh) | 基于星敏感器和星间链路的heo卫星编队飞行自主导航方法 | |
CN106679653A (zh) | 一种基于星敏感器和星间链路的heo卫星群相对测量方法 | |
Henry et al. | Absolute triangulation algorithms for space exploration | |
CN107883925B (zh) | 一种导航星座星间观测目标卫星图像模拟方法 | |
CN109269508A (zh) | 一种卫星相对小行星视觉自主导航方法 | |
CN106643742B (zh) | 一种卫星自主连续观测小行星的方法 | |
Eisenman et al. | Sun sensing on the Mars exploration rovers | |
Sigel et al. | Star tracker celestial localization system for a lunar rover | |
CN107830856B (zh) | 面向编队飞行的太阳tdoa量测方法及组合导航方法 | |
Pingyuan et al. | Attitude and position determination scheme of lunar rovers basing on the celestial vectors observation | |
Nayak et al. | Design of relative motion and attitude profiles for three-dimensional resident space object imaging with a laser rangefinder |
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 |