CN111947653A - 一种月面巡视探测器双模式惯性/视觉/天文导航方法 - Google Patents
一种月面巡视探测器双模式惯性/视觉/天文导航方法 Download PDFInfo
- Publication number
- CN111947653A CN111947653A CN202010811874.0A CN202010811874A CN111947653A CN 111947653 A CN111947653 A CN 111947653A CN 202010811874 A CN202010811874 A CN 202010811874A CN 111947653 A CN111947653 A CN 111947653A
- Authority
- CN
- China
- Prior art keywords
- error
- inertial
- measurement
- vector
- lunar
- 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.)
- Pending
Links
- 230000000007 visual effect Effects 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000007689 inspection Methods 0.000 title claims abstract description 17
- 239000013598 vector Substances 0.000 claims abstract description 85
- 241001061260 Emmelichthys struhsakeri Species 0.000 claims abstract description 46
- 230000003068 static effect Effects 0.000 claims abstract description 18
- 239000000523 sample Substances 0.000 claims abstract description 12
- 238000005259 measurement Methods 0.000 claims description 104
- 239000011159 matrix material Substances 0.000 claims description 50
- 230000009466 transformation Effects 0.000 claims description 26
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims description 23
- 238000001914 filtration Methods 0.000 claims description 12
- 230000001133 acceleration Effects 0.000 claims description 10
- 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 claims description 4
- 238000009434 installation Methods 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 description 3
- 241000283973 Oryctolagus cuniculus Species 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
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/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
- G01C21/16—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
- G01C21/165—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
-
- 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
- G01C21/025—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means with the use of startrackers
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)
- Navigation (AREA)
Abstract
本发明涉及一种月面巡视探测器双模式惯性/视觉/天文导航方法,当月面巡视探测器运动时,采用惯性/视觉/天文组合导航模式,利用天文导航提供的星光矢量及惯性导航***和视觉导航***提供的相对位置和姿态误差,对月面巡视探测器的姿态误差和惯性器件误差进行修正。当月面巡视探测器静止时,导航***工作在惯性/天文组合导航模式,可由天文导航***提供的星光矢量、星光高度和惯性导航***提供的水平速度误差修正导航***的位置和姿态误差。本发明可为月面巡视探测器提供高精度的位置、速度及姿态信息,对深空探测器自主导航具有重要的实际意义。
Description
技术领域
本发明属于航天器自主导航领域,涉及一种适用于月面巡视探测器“走走停停”运动模式的双模式自主惯性/视觉/天文导航方法。
背景技术
深空探测技术作为一个国家综合国力、军事应用价值和科学技术发展水平的重要特征与标志,已引起世界各国的极大关注。随着我国深空探测技术的不断发展,我国的深空探测能力不断提升,我国已经具备探测月球的能力。
月面巡视探测器如“玉兔二号”位于月球背面,无法靠地面站直接进行测控,虽然可通过中继星中转,但由于中继星的能力有限,会导致时延非常长,因此就要求月面巡视探测器必须具备长时间、长距离的高精度、高可靠、自主导航能力,但惯性导航误差随时间发散,视觉导航误差随距离发散,都无法达到这一要求,而天文导航的误差不受时间和距离的影响,可有效修正惯性和视觉导航的误差。因此,惯性导航、视觉导航和天文导航具有很好的优势互补特性,对于长时间、长距离的探测器自主导航任务,利用天文信息辅助惯性导航和视觉导航实现月面巡视探测器的自主导航是一种可行的方法。由于月面巡视探测器在运动时受水平加速度的影响,无法获取高精度的水平信息,此时天文导航不能提供位置信息,只能提供姿态信息。而对于处于静止状态的月面巡视探测器,由于没有水平运动加速度的干扰,利用加速度计等重力敏感器件等可获得高精度的水平基准,结合星敏感器的测量信息,就可获得天体高度,并进一步获得月面巡视探测器的位置。因此,在静止条件下,天文导航可以辅助惯性导航实现位置和姿态的获取。
发明内容
本发明要解决的技术问题是:克服单独使用某种导航方法存在的缺点,为月面巡视探测器提供一种将惯性导航、视觉导航和天文导航结合的适用于月面巡视探测器运动和静止两种工作模式的双模式自主导航方法,提供高精度的位置、速度和姿态信息。
本发明解决其技术问题所采用的技术方案为:一种月面巡视探测器双模式惯性/视觉/天文导航方法,包括如下步骤:
步骤1:获取星敏感器测量的星光矢量、天体高度,惯性导航***和视觉导航***测量的相对位置和姿态误差,以及惯性导航***测量的水平速度误差;
步骤2:当月面巡视探测器运动时,首先根据月球固连坐标系下捷联惯导的力学编排,以惯性导航误差方程作为状态模型;再分别以惯性导航***和视觉导航***获得的相对位置和姿态误差及星敏感器测量的星光矢量作为量测量,之后根据这些量测量分别建立惯性导航***和视觉导航***获得的相对位置和姿态误差的量测模型及星光矢量的量测模型;并使用无迹卡尔曼滤波UKF滤波估计月面巡视探测器的位置、速度和姿态,对月面巡视探测器的姿态误差和惯性器件误差进行修正;
步骤3:当月面巡视探测器静止时,速度为零,首先建立简化的惯导误差方程作为状态模型;再分别以惯性导航***获得的水平速度误差及星敏感器测量的星光矢量和星光高度作为量测量,之后根据这些量测量分别建立水平速度误差的量测模型及星光矢量、星光高度的量测模型;并使用UKF滤波估计月面巡视探测器的位置、速度和姿态,对月面巡视探测器的位置和姿态误差以及惯性器件误差进行修正。
进一步的,所述步骤2中,当月面巡视探测器运动时的具体导航步骤如下:
步骤2.1建立月面巡视探测器运动时惯性/视觉/天文组合导航***状态模型
以月球固连坐标系下捷联惯导的误差方程作为***状态模型:
式中,φ=[φE φN φU]T是姿态误差角,φE φN φU是地理坐标系下每个轴向的姿态误差角,fn是n系下的比力矢量,是速度误差,是地理坐标系下各轴向的速度误差,δrn=[δL δλ]T是位置误差,δL是经度误差δλ是纬度误差,是的误差,是巡视器在n系下的月球角速率,是的误差,是n系相对m系的旋转角速率在n系中的表示,是的误差,为n系相对惯性空间的转动角速度矢量在n系下的投影,ε=(εx εy εz)T是陀螺仪的常值漂移,εx εy εz是陀螺仪在本体坐标系下三个轴的常值漂移,是加速度计的常值偏置,是加计在本体坐标系下三个轴的常值偏差,变量上方的圆点表示导数;
步骤2.2建立惯性导航***和视觉导航***获得的相对位置和姿态误差及星光矢量的量测模型
首先,建立惯性导航***和视觉导航***获得的相对位置和姿态误差的量测模型的过程如下:
假设在tk时刻巡视器的真实姿态矩阵,Rerr是的估计误差,定义为rk是巡视器在tk时刻月球固连坐标系下的真实位置;δr(k)=rINS(k)-r(k)是rINS(k)的估计误差;姿态矩阵和位置rINS(k)表示为:
假设校正后tk-1时刻的估计误差很小,满足:
相对旋转矩阵RINS和位置矢量TINS可以通过把方程式(3)代入方程式(4)获得:
相对旋转矩阵RVNS位置矢量TVNS通过VNS获得可以表示为:.
式中,ΔRVNS和ΔTVNS表示RVNS和TVNS的误差
从等式(5)和等式(6)得到:
当姿态误差角φ小于预定常数,将φ表示为Rerr的函数:
φ=g(Rerr) (8)
式中,Δφ是φ的误差;
INS和VNS之间的位置矢量误差表示为:
λ(k)=λINS(k)+Δλ(k)
L(k)=LINS(k)+ΔL(k) (12)
式中,λINS(k)和LINS(k)是INS在tk时刻获得的经度和纬度;
从等式(9)和等式(10)中获得的相对运动误差的量测模型表示为:
式中V1=[Δφ,ΔTVNS]是量测噪声;
建立基于星光矢量的量测模型,Si和Ss之间的关系如下:
式中,ss是星敏感器坐标系下的恒星星光方向矢量,si是惯性坐标系下所观测恒星星光方向矢量,是从n系到m系的变换矩阵,是从b系到n系的变换矩阵,也是姿态矩阵,Rerr是的估计误差,是从m系到i系的变换矩阵,是星敏感器的安装矩阵,星敏感器坐标系相对于惯性坐标系的转换矩阵;
根据公式(14),星光矢量的量测模型简写为:
Z2=[ss]=h2(X)+V2 (15)
式中V2是星敏感器的量测噪声;
步骤2.3量测模型(13)、(15)都是非线性的,使用非线性UKF滤波方法,估计月面巡视探测器的位置、速度和姿态,并修正月面巡视探测器的姿态误差和惯性器件误差。
进一步的,所述步骤3中,当月面巡视探测器静止时的具体导航步骤如下:
步骤3.1月面巡视探测器静止时惯性/天文导航***状态模型
当月面巡视探测器静止时,其真实速度为零,简化的INS误差模型从公式(1)得到,如下所示:
步骤3.2量测模型
静止时的惯性/天文组合导航模式下的量测量包括惯导***获得的水平速度误差、星光矢量和星光高度;
首先,惯性导航***获得的水平速度误差的量测模型如下:
月面巡视探测器静止时惯导解算得到的水平速度即为速度误差δv,相应的量测方程为:
天文导航是以已知准确空间位置的自然天体为基准,通过天体测量仪器测得的天体高度和星光方位矢量信息进行导航;
其次,基于星光方位矢量的量测模型如下:
静止时的惯性/天文组合导航模式下星光矢量的量测方程如下所示:
方程(18)是X1的非线性方程;使用Z1=[ss],量测方程表示为:
Z1=h(X1)+V1 (19)
式中V1是量测噪声;
最后,星光高度的量测模型如下:
星光高度是天体投影点与月面巡视探测器位置(L,λ)的函数,它们之间的关系为:
H=arcsin(sinLsinΔ+cosLcosΔcos(λ+tSHA-RA)) (20)
式中,Δ、RA分别为天体的赤纬和赤经;
天体高度利用加速度计提供的水平基准和星敏感器提供的星光方位矢量获取;当加速度计不存在测量误差时:
式中,fbr为不存在测量误差时的加速度计输出,sb为本体坐标系下的星光方向矢量;当加速度存在测量误差时:
将式(22)代入式(21)得到:
式中,gn=[0 0 -gm]T,gm=1.618m/s2是月球重力加速度,通过公式(17)、(18)、(23),以Z2=[δv ss sinHm]作为量测量,得量测方程为:
Z2=h(X2)+V2 (24)
式中,V2为量测噪声;
步骤3.3由于量测模型(19)、(24)都是非线性的,因此使用非线性UKF滤波方法,估计月面巡视探测器的位置、速度和姿态,并修正月面巡视探测器的姿态误差和惯性器件误差。
本发明与现有技术相比的优点在于:
(1)利用多种导航方法并基于信息融合技术进行组合导航,实现它们之间的取长补短和优势互补,最大限度地提高导航的精度和可靠性,实现对航天器的高精度自主导航。
(2)月面巡视探测器运动时,引入天体方向矢量量测信息,大幅度提高了姿态精度,尤其是航向角精度,最大航向误差减小了97.3%,终点航向误差减小了99.3%。通过航向角精度的提高,又间接提高了位置精度,最大位置误差减小了41.7%,终点误差减小了82.05%。
附图说明
图1为本发明月面巡视探测器双模式惯性/视觉/天文导航方法***流程图;
图2为本发明中月面巡视探测器运动时以相对运动参数差和星光矢量为量测量的惯性/视觉/天文导航方法***流程图;
图3为本发明中月面巡视探测器静止时以水平速度误差、星光矢量和星光高度为量测量的惯性/天文导航方法***流程图;
图4为本发明中导航坐标系示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例仅为本发明的一部分实施例,而不是全部的实施例,基于本发明中的实施例,本领域的普通技术人员在不付出创造性劳动的前提下所获得的所有其他实施例,都属于本发明的保护范围。
图1给出了月面巡视探测器双模式惯性/视觉/天文导航方法***流程图。接下来详细说明本发明的具体实施过程:
月面巡视探测器运动时的惯性/视觉/天文导航模式:
①建立月面巡视探测器运动时惯性/视觉/天文导航***状态模型
以月球固连坐标系下捷联惯导的误差方程作为***状态模型:
式中,φ=[φE φN φU]T是姿态误差角,φE φN φU是地理坐标系下每个轴向的姿态误差角,fn是n系下的比力矢量,是速度误差,是地理坐标系下各轴向的速度误差,δrn=[δL δλ]T是位置误差,δL是经度误差δλ是纬度误差,是的误差,是巡视器在n系下的月球角速率,是的误差,是n系相对m系的旋转角速率在n系中的表示,是的误差,为n系相对惯性空间的转动角速度矢量在n系下的投影,ε=(εx εy εz)T是陀螺仪的常值漂移,εx εy εz是陀螺仪在本体坐标系下三个轴的常值漂移,是加速度计的常值偏置,是加计在本体坐标系下三个轴的常值偏差,变量上方的圆点表示导数。
②建立惯性导航***和视觉导航***获得的相对位置和姿态误差及星光矢量的量测模型
首先介绍惯性导航***和视觉导航***获得的相对位置和姿态误差的量测模型。假设在tk时刻巡视器的真实姿态矩阵,Rerr是的估计误差,定义为r(k)是巡视器在tk时刻月球固连坐标系下的真实位置。δr(k)=rINS(k)-r(k)是rINS(k)的估计误差。姿态矩阵和位置rINS(k)可以表示为:
假设校正后tk-1时刻的估计误差很小,满足:
相对旋转矩阵RINS和位置矢量TINS可以通过把方程式(3)代入方程式(4)获得:
相对旋转矩阵RVNS位置矢量TVNS通过VNS获得可以表示为:.
式中,ΔRVNS和ΔTVNS表示RVNS和TVNS的误差
从等式(5)和等式(6)可以得到:
当姿态误差角φ非常小时,可以将φ表示为Rerr的函数:
φ=g(Rerr) (8)
式中,Δφ是φ的误差。
INS和VNS之间的位置矢量误差可以表示为:
λ(k)=λINS(k)+Δλ(k)
L(k)=LINS(k)+ΔL(k) (12)
式中,λINS(k)和LINS(k)是INS在tk时刻获得的经度和纬度。
从等式(9)和等式(10)中获得的相对运动误差的量测模型可以简单表示为:
式中V1=[Δφ,ΔTVNS]是量测噪声。
接下来介绍基于星光矢量的量测模型。Si和Ss之间的关系如下:
式中,ss是星敏感器坐标系下的恒星星光方向矢量,si是惯性坐标系下所观测恒星星光方向矢量,是从n系到m系的变换矩阵,是从b系到n系的变换矩阵,也是姿态矩阵,Rerr是的估计误差,是从m系到i系的变换矩阵,是星敏感器的安装矩阵,星敏感器坐标系相对于惯性坐标系的转换矩阵。
根据公式(14),星光矢量的量测模型可以简写为:
Z2=[ss]=h2(X)+V2 (15)
式中V2是星敏感器的量测噪声。
③由于量测模型(13)、(15)都是非线性的,因此使用了非线性UKF滤波方法。估计月面巡视探测器的位置、速度和姿态,并修正月面巡视探测器的姿态误差和惯性器件误差。
月面巡视探测器静止时的惯性/天文导航模式:
①月面巡视探测器静止时惯性/天文导航***状态模型
当月面巡视探测器静止时,其真实速度为零,简化的INS误差模型可以从公式(1)得到,如下所示:
②量测模型
静止时的惯性/天文组合导航模式下的量测量包括惯导***获得的的水平速度误差、星光矢量和星光高度。
首先介绍惯性导航***获得的水平速度误差的量测模型。
月面巡视探测器静止时惯导解算得到的水平速度即为速度误差δv,相应的量测方程为:
天文导航是以已知准确空间位置的自然天体为基准,通过天体测量仪器测得的天体高度和星光方位矢量等信息进行导航。
接下来介绍基于星光方位矢量的量测模型。
静止时的惯性/天文导航模式下星光矢量的量测方程如下所示:
方程(18)是X1的非线性方程。使用Z1=[ss],量测方程可以表示为:
Z1=h(X1)+V1 (19)
式中V1是量测噪声。
最后介绍星光高度的量测模型
星光高度是天体投影点与月面巡视探测器位置(L,λ)的函数,它们之间的关系为:
H=arcsin(sinLsinΔ+cosLcosΔcos(λ+tSHA-RA)) (20)
式中,Δ、RA分别为天体的赤纬和赤经。
天体高度利用加速度计提供的水平基准和星敏感器提供的星光方位矢量获取。当加速度计不存在测量误差时:
式中,fbr为不存在测量误差时的加速度计输出,sb为本体坐标系下的星光方向矢量。当加速度存在测量误差时:
将式(22)代入式(21)可以得到:
式中,gn=[0 0 -gm]T,gm=1.618m/s2是月球重力加速度,通过公式(17)、(18)、(23),以Z2=[δv ss sinHm]作为量测量,可得量测方程为:
Z2=h(X2)+V2 (24)
式中,V2为量测噪声。
图2给出了月面巡视探测器运动时以相对运动参数差和星光矢量为量测量的惯性/视觉/天文导航方法***流程图。接下来详细说明本发明的具体实施过程:
运动时的惯性/视觉/天文导航模式:
①获取星敏感器测量的星光矢量,惯性导航***和视觉导航***测量的相对位置和姿态误差;
②建立月面巡视探测器运动时惯性/视觉/天文导航***状态模型。
以月球固连坐标系下捷联惯导的误差方程作为***状态模型:
式中,φ=[φE φN φU]T是姿态误差角,φE φN φU是地理坐标系下每个轴向的姿态误差角,fn是n系下的比力矢量,是速度误差,是地理坐标系下各轴向的速度误差,δrn=[δL δλ]T是位置误差,δL是经度误差δλ是纬度误差,是的误差,是巡视器在n系下的月球角速率,是的误差,是n系相对m系的旋转角速率在n系中的表示,是的误差,为n系相对惯性空间的转动角速度矢量在n系下的投影,ε=(εx εy εz)T是陀螺仪的常值漂移,εx εy εz是陀螺仪在本体坐标系下三个轴的常值漂移,是加速度计的常值偏置,是加计在本体坐标系下三个轴的常值偏差,变量上方的圆点表示导数。
③建立惯性导航***和视觉导航***获得的相对位置和姿态误差及星光矢量的量测模型
首先介绍惯性导航***和视觉导航***获得的相对位置和姿态误差的量测模型。假设在tk时刻巡视器的真实姿态矩阵,Rerr是的估计误差,定义为rk是巡视器在tk时刻月球固连坐标系下的真实位置。δr(k)=rINS(k)-r(k)是rINS(k)的估计误差。姿态矩阵和位置rINS(k)可以表示为:
假设校正后tk-1时刻的估计误差很小,满足:
相对旋转矩阵RINS和位置矢量TINS可以通过把方程式(3)代入方程式(4)获得:
相对旋转矩阵RVNS位置矢量TVNS通过VNS获得可以表示为:.
式中,ΔRVNS和ΔTVNS表示RVNS和TVNS的误差
从等式(5)和等式(6)可以得到:
当姿态误差角φ非常小时,可以将φ表示为Rerr的函数:
φ=g(Rerr) (8)
式中,Δφ是φ的误差。
INS和VNS之间的位置矢量误差可以表示为:
λ(k)=λINS(k)+Δλ(k)
L(k)=LINS(k)+ΔL(k) (12)
式中,λINS(k)和LINS(k)是INS在tk时刻获得的经度和纬度。
从等式(9)和等式(10)中获得的相对运动误差的量测模型可以简单表示为:
式中V1=[Δφ,ΔTVNS]是量测噪声。
接下来介绍基于星光矢量的量测模型。Si和Ss之间的关系如下:
式中,ss是星敏感器坐标系下的恒星星光方向矢量,si是惯性坐标系下所观测恒星星光方向矢量,是从n系到m系的变换矩阵,是从b系到n系的变换矩阵,也是姿态矩阵,Rerr是的估计误差,是从m系到i系的变换矩阵,是星敏感器的安装矩阵,星敏感器坐标系相对于惯性坐标系的转换矩阵。
根据公式(14),星光矢量的量测模型可以简写为:
Z2=[ss]=h2(X)+V2 (15)
式中V2是星敏感器的量测噪声。
④由于量测模型(13)、(15)都是非线性的,因此使用了非线性UKF滤波方法。估计月面巡视探测器的位置、速度和姿态,并修正月面巡视探测器的姿态误差和惯性器件误差。
图3给出了月面巡视探测器静止时以水平速度误差、星光矢量和星光高度为量测量的惯性/天文组合导航方法***流程图。接下来详细说明本发明的具体实施过程:
静止时的惯性/天文导航模式:
①获取星敏感器测量的星光矢量、天体高度,惯性导航***测量的水平速度误差;
②月面巡视探测器静止时惯性/天文导航***状态模型
当月面巡视探测器静止时,其真实速度为零,简化的INS误差模型可以从公式(1)得到,如下所示:
③量测模型
静止时的惯性/天文导航模式下的量测量包括惯导***获得的的水平速度误差、星光矢量和星光高度。
首先介绍惯性导航***获得的水平速度误差的量测模型。
月面巡视探测器静止时惯导解算得到的水平速度即为速度误差δv,相应的量测方程为:
天文导航是以已知准确空间位置的自然天体为基准,通过天体测量仪器测得的天体高度和星光方位矢量等信息进行导航。
接下来介绍基于星光方位矢量的量测模型。
静止时的惯性/天文导航模式下星光矢量的量测方程如下所示:
方程(18)是X1的非线性方程。使用Z1=[ss],量测方程可以表示为:
Z1=h(X1)+V1 (19)
式中V1是量测噪声。
最后介绍星光高度的量测模型。
星光高度是天体投影点与月面巡视探测器位置(L,λ)的函数,它们之间的关系为:
H=arcsin(sinLsinΔ+cosLcosΔcos(λ+tSHA-RA)) (20)
式中,Δ、RA分别为天体的赤纬和赤经,tSHA是春分点的格林时角,可以通过观测时间查找天文历表获得。
天体高度利用加速度计提供的水平基准和星敏感器提供的星光方位矢量获取。当加速度计不存在测量误差时:
式中,fbr为不存在测量误差时的加速度计输出,sb为本体坐标系下的星光方向矢量。当加速度存在测量误差时:
式中,fb是加速度计的真实输出。
将式(22)代入式(21)可以得到:
式中,gn=[0 0 -gm]T,gm=1.618m/s2是月球重力加速度,通过公式(17)、(18)、(23),以Z2=[δv ss sinHm]作为量测量,可得量测方程为:
Z2=h(X2)+V2 (24)
式中,V2为量测噪声。
图4给出了导航坐标系示意图,介绍了惯性导航和视觉导航中涉及的常用坐标系。本文涉及四个坐标系:惯性坐标系(i系),月球固连坐标系(m系),地理坐标系(n系),巡视器本体坐标系(b系)。
①惯性坐标系(i系)
惯性坐标系的原点位于月球质心Om,zi轴垂直于赤道平面,且与月球自转角速率矢量的方向一致,xi轴在赤道平面内并指向春分点方向,yi轴与xi轴和zi轴均垂直,并构成右手直角坐标系。
②月球固连坐标系(m系)
月球固连坐标系的坐标原点位于月球质心Om,轴垂直于赤道平面,且与月球自转角速率矢量的方向一致,xm轴在赤道平面内并指向零度经线方向,ym轴垂直于xm轴和轴所在的平面,并构成右手直角坐标系。从月球惯性坐标系到月球固连坐标系的坐标转换矩阵为:
式中,tSHA是春分点格林时角,可通过观测时间查找天文历表获得。
③地理坐标系(n系)
地理坐标系为当地东北天坐标系,其原点位于巡视器质心O,xn轴沿月球参考椭球的卯酉圈方向并指向东向,yn轴沿月球参考椭球的子午圈方向并指向北向,zn轴沿月球参考椭球的外法线方向指向天顶。从地理坐标系到月球固连坐标系的坐标转换矩阵为:
式中,Rx(90°-L)和Rz(90°+L)分别表示绕x轴和z轴的旋转矩阵,L,λ分别为巡视器所在位置的纬度和经度。
④巡视器本体坐标系(b系)
巡视器本体坐标系与巡视器刚性固连,坐标原点位于巡视器质心O,xboyb平面与巡视器甲板平面平行,xb轴指向巡视器右方,yb轴沿巡视器对称轴方向指向车体前方,zb轴垂直于xboyb平面并指向上方,构成右手直角坐标系。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,且应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (3)
1.一种月面巡视探测器双模式惯性/视觉/天文导航方法,其特征在于,包括如下步骤:
步骤1:获取星敏感器测量的星光矢量、天体高度,惯性导航***和视觉导航***测量的相对位置和姿态误差,以及惯性导航***测量的水平速度误差;
步骤2:当月面巡视探测器运动时,首先根据月球固连坐标系下捷联惯导的力学编排,以惯性导航误差方程作为状态模型;再分别以惯性导航***和视觉导航***获得的相对位置和姿态误差及星敏感器测量的星光矢量作为量测量,之后根据这些量测量分别建立惯性导航***和视觉导航***获得的相对位置和姿态误差的量测模型及星光矢量的量测模型;并使用无迹卡尔曼滤波UKF滤波估计月面巡视探测器的位置、速度和姿态,对月面巡视探测器的姿态误差和惯性器件误差进行修正;
步骤3:当月面巡视探测器静止时,速度为零,首先建立简化的惯导误差方程作为状态模型;再分别以惯性导航***获得的水平速度误差及星敏感器测量的星光矢量和星光高度作为量测量,之后根据这些量测量分别建立水平速度误差的量测模型及星光矢量、星光高度的量测模型;并使用UKF滤波估计月面巡视探测器的位置、速度和姿态,对月面巡视探测器的位置和姿态误差以及惯性器件误差进行修正。
2.根据权利要求1所述的一种月面巡视探测器双模式惯性/视觉/天文导航方法,其特征在于:所述步骤2中,当月面巡视探测器运动时的具体导航步骤如下:
步骤2.1建立月面巡视探测器运动时惯性/视觉/天文组合导航***状态模型
以月球固连坐标系下捷联惯导的误差方程作为***状态模型:
式中,φ=[φE φN φU]T是姿态误差角,φE φN φU是地理坐标系下每个轴向的姿态误差角,fn是n系下的比力矢量,是速度误差,是地理坐标系下各轴向的速度误差,δrn=[δL δλ]T是位置误差,L为经度,λ是纬度,δL是经度误差,δλ是纬度误差,是的误差,是巡视器在n系下的月球角速率,是的误差,是n系相对m系的旋转角速率在n系中的表示,是的误差,为n系相对惯性空间的转动角速度矢量在n系下的投影,ε=(εx εy εz)T是陀螺仪的常值漂移,εxεyεz是陀螺仪在本体坐标系下三个轴的常值漂移,▽=(▽x ▽y ▽z)T是加速度计的常值偏置,▽x ▽y ▽z是加计在本体坐标系下三个轴的常值偏差,变量上方的圆点表示导数;
状态量:X=[φT δvT δrT εT ▽T]T;
步骤2.2建立惯性导航***和视觉导航***获得的相对位置和姿态误差及星光矢量的量测模型
首先,建立惯性导航***和视觉导航***获得的相对位置和姿态误差的量测模型的过程如下:
假设在tk时刻巡视器的真实姿态矩阵,Rerr是的估计误差,定义为rk是巡视器在tk时刻月球固连坐标系下的真实位置;δr(k)=rINS(k)-r(k)是rINS(k)的估计误差;姿态矩阵和位置rINS(k)表示为:
假设校正后tk-1时刻的估计误差很小,满足:
相对旋转矩阵RINS和位置矢量TINS可以通过把方程式(3)代入方程式(4)获得:
相对旋转矩阵RVNS位置矢量TVNS通过视觉导航***VNS获得表示为:.
式中,ΔRVNS和ΔTVNS表示RVNS和TVNS的误差
从等式(5)和等式(6)得到:
当姿态误差角φ小于预定常数,将φ表示为Rerr的函数:
φ=g(Rerr) (8)
式中,Δφ是φ的误差;
INS表示惯性导航***,VNS表示视觉导航***,INS和VNS之间的位置矢量误差表示为:
λ(k)=λINS(k)+Δλ(k)
L(k)=LINS(k)+ΔL(k) (12)
式中,λINS(k)和LINS(k)是惯性导航***INS在tk时刻获得的经度和纬度;
从等式(9)和等式(10)中获得的相对运动误差的量测模型表示为:
式中V1=[Δφ,ΔTVNS]是量测噪声;
建立基于星光矢量的量测模型,Si和Ss之间的关系如下:
式中,ss是星敏感器坐标系下的恒星星光方向矢量,si是惯性坐标系下所观测恒星星光方向矢量,是从n系到m系的变换矩阵,是从b系到n系的变换矩阵,也是姿态矩阵,Rerr是的估计误差,是从m系到i系的变换矩阵,是星敏感器的安装矩阵,星敏感器坐标系相对于惯性坐标系的转换矩阵;
根据公式(14),星光矢量的量测模型简写为:
Z2=[ss]=h2(X)+V2 (15)
式中V2是星敏感器的量测噪声;
步骤2.3量测模型(13)、(15)都是非线性的,使用非线性UKF滤波方法,估计月面巡视探测器的位置、速度和姿态,并修正月面巡视探测器的姿态误差和惯性器件误差。
3.根据权利要求1所述的一种月面巡视探测器双模式惯性/视觉/天文导航方法,其特征在于:
所述步骤3中,当月面巡视探测器静止时的具体导航步骤如下:
步骤3.1月面巡视探测器静止时惯性/天文导航***状态模型
当月面巡视探测器静止时,其真实速度为零,简化的INS误差模型从公式(1)得到,如下所示:
状态量:X1=[φT δvT δrT εT ▽T]T
步骤3.2量测模型
静止时的惯性/天文组合导航模式下的量测量包括惯导***获得的水平速度误差、星光矢量和星光高度;
首先,惯性导航***获得的水平速度误差的量测模型如下:
月面巡视探测器静止时惯导解算得到的水平速度即为速度误差δv,相应的量测方程为:
天文导航是以已知准确空间位置的自然天体为基准,通过天体测量仪器测得的天体高度和星光方位矢量信息进行导航;
其次,基于星光方位矢量的量测模型如下:
静止时的惯性/天文组合导航模式下星光矢量的量测方程如下所示:
方程(18)是X1的非线性方程;使用Z1=[ss],量测方程表示为:
Z1=h(X1)+V1 (19)
式中V1是量测噪声;
最后,星光高度的量测模型如下:
星光高度是天体投影点与月面巡视探测器位置(L,λ)的函数,它们之间的关系为:
H=arcsin(sinLsinΔ+cosLcosΔcos(λ+tSHA-RA)) (20)
式中,Δ、RA分别为天体的赤纬和赤经;
天体高度利用加速度计提供的水平基准和星敏感器提供的星光方位矢量获取;当加速度计不存在测量误差时:
式中,fbr为不存在测量误差时的加速度计输出,sb为本体坐标系下的星光方向矢量;当加速度存在测量误差时:
将式(22)代入式(21)得到:
式中,gn=[0 0 -gm]T,gm=1.618m/s2是月球重力加速度,通过公式(17)、(18)、(23),以Z2=[δv ss sinHm]作为量测量,得量测方程为:
Z2=h(X2)+V2 (24)
式中,V2为量测噪声;
步骤3.3由于量测模型(19)、(24)都是非线性的,因此使用非线性UKF滤波方法,估计月面巡视探测器的位置、速度和姿态,并修正月面巡视探测器的姿态误差和惯性器件误差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010811874.0A CN111947653A (zh) | 2020-08-13 | 2020-08-13 | 一种月面巡视探测器双模式惯性/视觉/天文导航方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010811874.0A CN111947653A (zh) | 2020-08-13 | 2020-08-13 | 一种月面巡视探测器双模式惯性/视觉/天文导航方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111947653A true CN111947653A (zh) | 2020-11-17 |
Family
ID=73332505
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010811874.0A Pending CN111947653A (zh) | 2020-08-13 | 2020-08-13 | 一种月面巡视探测器双模式惯性/视觉/天文导航方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111947653A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113030517A (zh) * | 2021-02-18 | 2021-06-25 | 北京控制工程研究所 | 一种火星着陆过程利用测速敏感器的姿态修正方法 |
CN113551667A (zh) * | 2021-07-21 | 2021-10-26 | 北京航空航天大学 | 航天器惯性/恒星星光矢量/太阳多普勒速度组合导航方法 |
CN113551668A (zh) * | 2021-07-21 | 2021-10-26 | 北京航空航天大学 | 一种航天器惯性/恒星星光矢量/星光折射组合导航方法 |
CN113849003A (zh) * | 2021-10-13 | 2021-12-28 | 西安尹纳数智能科技有限公司 | 一种动中通天线运动隔离的控制方法 |
CN114111723A (zh) * | 2021-11-15 | 2022-03-01 | 北京航天计量测试技术研究所 | 一种基于动态小视场的天体捕获方法 |
CN114509071A (zh) * | 2022-04-20 | 2022-05-17 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种风洞试验模型姿态测量方法 |
CN114526726A (zh) * | 2022-02-09 | 2022-05-24 | 中国人民解放***箭军研究院科技创新研究中心 | 基于可观性分析的星光折射导航观星方案优化设计方法 |
CN114608564A (zh) * | 2022-05-11 | 2022-06-10 | 北京航空航天大学 | 一种基于夜间月光偏振-星光信息融合的组合定位方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060085130A1 (en) * | 2004-10-18 | 2006-04-20 | Trex Enterprises Corporation | Daytime stellar imager for attitude determination |
CN103033189A (zh) * | 2012-12-26 | 2013-04-10 | 北京航空航天大学 | 一种深空探测巡视器惯性/视觉组合导航方法 |
CN106932804A (zh) * | 2017-03-17 | 2017-07-07 | 南京航空航天大学 | 天文辅助的惯性/北斗紧组合导航***及其导航方法 |
CN107045136A (zh) * | 2017-03-17 | 2017-08-15 | 南京航空航天大学 | 可配置的惯性/天文/北斗多组合导航***及其导航方法 |
CN107576327A (zh) * | 2017-08-07 | 2018-01-12 | 西南技术物理研究所 | 基于可观测度分析的可变结构综合导航***设计方法 |
CN108871326A (zh) * | 2018-07-09 | 2018-11-23 | 北京航空航天大学 | 一种单轴旋转调制惯性-天文深组合导航方法 |
CN109631876A (zh) * | 2019-01-18 | 2019-04-16 | 辽宁工程技术大学 | 一种基于单相机导航影像的巡视探测器定位方法 |
CN111121766A (zh) * | 2019-12-17 | 2020-05-08 | 南京理工大学 | 一种基于星光矢量的天文与惯性组合导航方法 |
CN111238467A (zh) * | 2020-02-07 | 2020-06-05 | 西北工业大学 | 一种仿生偏振光辅助的无人作战飞行器自主导航方法 |
-
2020
- 2020-08-13 CN CN202010811874.0A patent/CN111947653A/zh active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060085130A1 (en) * | 2004-10-18 | 2006-04-20 | Trex Enterprises Corporation | Daytime stellar imager for attitude determination |
CN103033189A (zh) * | 2012-12-26 | 2013-04-10 | 北京航空航天大学 | 一种深空探测巡视器惯性/视觉组合导航方法 |
CN106932804A (zh) * | 2017-03-17 | 2017-07-07 | 南京航空航天大学 | 天文辅助的惯性/北斗紧组合导航***及其导航方法 |
CN107045136A (zh) * | 2017-03-17 | 2017-08-15 | 南京航空航天大学 | 可配置的惯性/天文/北斗多组合导航***及其导航方法 |
CN107576327A (zh) * | 2017-08-07 | 2018-01-12 | 西南技术物理研究所 | 基于可观测度分析的可变结构综合导航***设计方法 |
CN108871326A (zh) * | 2018-07-09 | 2018-11-23 | 北京航空航天大学 | 一种单轴旋转调制惯性-天文深组合导航方法 |
CN109631876A (zh) * | 2019-01-18 | 2019-04-16 | 辽宁工程技术大学 | 一种基于单相机导航影像的巡视探测器定位方法 |
CN111121766A (zh) * | 2019-12-17 | 2020-05-08 | 南京理工大学 | 一种基于星光矢量的天文与惯性组合导航方法 |
CN111238467A (zh) * | 2020-02-07 | 2020-06-05 | 西北工业大学 | 一种仿生偏振光辅助的无人作战飞行器自主导航方法 |
Non-Patent Citations (2)
Title |
---|
XIAOLIN NING,AT AL: "A Two-Mode INS/CNS Navigation Method for Lunar Rovers", 《IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT》 * |
XIAOLIN NING,AT AL: "INS/VNS/CNS integrated navigation method for planetary rovers", 《AEROSPACE SCIENCE AND TECHNOLOGY》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113030517B (zh) * | 2021-02-18 | 2022-10-28 | 北京控制工程研究所 | 一种火星着陆过程利用测速敏感器的姿态修正方法 |
CN113030517A (zh) * | 2021-02-18 | 2021-06-25 | 北京控制工程研究所 | 一种火星着陆过程利用测速敏感器的姿态修正方法 |
CN113551667A (zh) * | 2021-07-21 | 2021-10-26 | 北京航空航天大学 | 航天器惯性/恒星星光矢量/太阳多普勒速度组合导航方法 |
CN113551668A (zh) * | 2021-07-21 | 2021-10-26 | 北京航空航天大学 | 一种航天器惯性/恒星星光矢量/星光折射组合导航方法 |
CN113551668B (zh) * | 2021-07-21 | 2024-05-28 | 北京航空航天大学 | 一种航天器惯性/恒星星光矢量/星光折射组合导航方法 |
CN113551667B (zh) * | 2021-07-21 | 2024-05-24 | 北京航空航天大学 | 航天器惯性/恒星星光矢量/太阳多普勒速度组合导航方法 |
CN113849003A (zh) * | 2021-10-13 | 2021-12-28 | 西安尹纳数智能科技有限公司 | 一种动中通天线运动隔离的控制方法 |
CN113849003B (zh) * | 2021-10-13 | 2024-04-26 | 复远芯(上海)科技有限公司 | 一种动中通天线运动隔离的控制方法 |
CN114111723A (zh) * | 2021-11-15 | 2022-03-01 | 北京航天计量测试技术研究所 | 一种基于动态小视场的天体捕获方法 |
CN114526726A (zh) * | 2022-02-09 | 2022-05-24 | 中国人民解放***箭军研究院科技创新研究中心 | 基于可观性分析的星光折射导航观星方案优化设计方法 |
CN114509071A (zh) * | 2022-04-20 | 2022-05-17 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种风洞试验模型姿态测量方法 |
CN114608564A (zh) * | 2022-05-11 | 2022-06-10 | 北京航空航天大学 | 一种基于夜间月光偏振-星光信息融合的组合定位方法 |
CN114608564B (zh) * | 2022-05-11 | 2022-07-29 | 北京航空航天大学 | 一种基于夜间月光偏振-星光信息融合的组合定位方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111947653A (zh) | 一种月面巡视探测器双模式惯性/视觉/天文导航方法 | |
CN111947652B (zh) | 一种适用于月球着陆器的惯性/视觉/天文/激光测距组合导航方法 | |
CN108051866B (zh) | 基于捷联惯性/gps组合辅助水平角运动隔离的重力测量方法 | |
CN101788296B (zh) | 一种sins/cns深组合导航***及其实现方法 | |
CN103994763B (zh) | 一种火星车的sins/cns深组合导航***及其实现方法 | |
CN101881619B (zh) | 基于姿态测量的船用捷联惯导与天文定位方法 | |
CN106289246B (zh) | 一种基于位置和姿态测量***的柔性杆臂测量方法 | |
EP2557394B1 (en) | System for processing pulse signals within an inertial navigation system | |
CN102879011B (zh) | 一种基于星敏感器辅助的月面惯导对准方法 | |
CN103900576B (zh) | 一种深空探测自主导航的信息融合方法 | |
CN110221332A (zh) | 一种车载gnss/ins组合导航的动态杆臂误差估计和补偿方法 | |
CN103900611B (zh) | 一种惯导天文高精度复合两位置对准及误差标定方法 | |
CN104880192B (zh) | 一种基于偏振罗盘的载体航向角计算方法 | |
CN112595350B (zh) | 一种惯导***自动标定方法及终端 | |
CN110296719B (zh) | 一种在轨标定方法 | |
CN113551668B (zh) | 一种航天器惯性/恒星星光矢量/星光折射组合导航方法 | |
CN107677292B (zh) | 基于重力场模型的垂线偏差补偿方法 | |
CN115327588A (zh) | 一种基于网络rtk的无人自动化作业特种车高精度定位方法 | |
CN112880669A (zh) | 一种航天器星光折射和单轴旋转调制惯性组合导航方法 | |
CN116105730A (zh) | 基于合作目标卫星甚短弧观测的仅测角光学组合导航方法 | |
CN109084755B (zh) | 一种基于重力视速度与参数辨识的加速度计零偏估计方法 | |
Nebot et al. | Initial calibration and alignment of an inertial navigation | |
CN112525204A (zh) | 一种航天器惯性和太阳多普勒速度组合导航方法 | |
CN113465570B (zh) | 一种基于高精度imu的气浮台初始对准方法及*** | |
Ploen et al. | Bias compensated inertial navigation for Venus balloon missions |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20201117 |
|
RJ01 | Rejection of invention patent application after publication |