CN101216319A - 基于联邦ukf算法的低轨卫星多传感器容错自主导航方法 - Google Patents
基于联邦ukf算法的低轨卫星多传感器容错自主导航方法 Download PDFInfo
- Publication number
- CN101216319A CN101216319A CNA2008100191018A CN200810019101A CN101216319A CN 101216319 A CN101216319 A CN 101216319A CN A2008100191018 A CNA2008100191018 A CN A2008100191018A CN 200810019101 A CN200810019101 A CN 200810019101A CN 101216319 A CN101216319 A CN 101216319A
- Authority
- CN
- China
- Prior art keywords
- earth
- satellite
- equation
- measurement
- chi
- 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
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 22
- 238000005259 measurement Methods 0.000 claims abstract description 84
- 239000011159 matrix material Substances 0.000 claims abstract description 32
- 230000004927 fusion Effects 0.000 claims abstract description 18
- 238000005070 sampling Methods 0.000 claims abstract description 9
- 238000001914 filtration Methods 0.000 claims description 40
- 238000010998 test method Methods 0.000 claims description 8
- 230000005358 geomagnetic field Effects 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 230000009897 systematic effect Effects 0.000 claims description 6
- 238000002955 isolation Methods 0.000 claims description 5
- 239000000654 additive Substances 0.000 claims description 3
- 230000000996 additive effect Effects 0.000 claims description 3
- 230000003094 perturbing effect Effects 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims description 3
- 230000026676 system process Effects 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 3
- 230000007257 malfunction Effects 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000010365 information processing Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000033772 system development Effects 0.000 description 1
- 230000035899 viability Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Landscapes
- Navigation (AREA)
Abstract
一种基于联邦UKF算法的低轨卫星多传感器容错自主导航方法,属卫星自主导航方法。该方法包括:建立直角坐标系的地球卫星轨道动力学方程;建立以星敏感器和红外地球敏感器的输出值为测量量的子***量测方程;建立以磁强计和雷达高度计的输出值为量测量的子***量测方程;建立以紫外敏感器的输出值为量测量的子***量测方程;选择Sigma采样点;建立离散型UKF算法的预测方程和更新方程;各子***分别独立进行Sigma采样点计算及预测更新和量测更新;根据预测滤波残差判断各子滤波器输出是否有效,如有故障将其隔离,否则将滤波结果输入主滤波器进行信息融合;建立基于UKF算法的无复位联邦UKF滤波方程;按照上述步骤,输出为地球卫星状态估计值及方差矩阵P。
Description
技术领域
本发明属于卫星多传感器信息融合自主导航技术领域,涉及一种基于联邦UKF算法的卫星多传感器容错自主导航方法,可用于近地轨道地球卫星的自主导航定位。
背景技术
卫星自主导航是指卫星在不依赖地面站的条件下,仅依靠卫星上的测量设备实时确定卫星的位置和速度,在轨完成的飞行任务要求的功能或操作。航天器通常采用的依靠地面站无线电测控进行定位导航的方法,由于受地理条件的限制,难以实现对中低轨道卫星整个轨道的定位导航;同时,随着空间运行的各种航天器的数目的增加,信息传输量剧增,完全依靠地面站测控,将会引起测控***的信息阻塞和地面站负荷过重。因此,采用自主导航技术,提高卫星的自主运行、自主管理和在轨生存能力,已经成为现代卫星导航***发展的必然趋势。
在卫星自主导航***中,为了提高可靠性和精度,通常采用多种传感器同时进行定轨。目前,星上常用的自主导航传感器有以下几类:利用地球物理特性的敏感器,如红外地平仪、磁强计、雷达高度计等;利用天***置的敏感器,如太阳敏感器、星敏感器、紫外敏感器等;利用惯性的敏感器,如陀螺仪、角加速度计等。这些传感器的测量方程都具有较强的非线性,如何将其信息有效的组织并充分利用,是卫星自主导航的关键问题。传统的信息处理方法是采用集中式扩展卡尔曼滤波器,但是这种方法存在两个缺陷:一是采用集中式处理方式,导致计算负担过重,容错性差、通信负担重;二是这种方法的假设前提十分苛刻,当***非线性特性较强或对噪声估计不准确时,将会造成滤波器发散,导致导航精度下降。
Unscented卡尔曼滤波(UKF)是一种新的信息处理方法,它采用卡尔曼线性滤波框架,使用UT变换来处理一步预测方程中的均值和协方差的非线性传递,不需要求导计算Jacobian矩阵,算法简单,容易实现,利用UKF算法可以获得比传统的扩展卡尔曼滤波(EKF)更高的估计精度。近年来,基于UKF算法的轨道及姿态确定方法研究得到迅速发展。如何采用信息分配原理,将Unscented卡尔曼滤波算法和信息融合技术相结合,实现基于UKF的多传感器非线性滤波算法,使***能够充分运用各导航***的信息进行信息互补和信息融合,并具有故障检测、隔离和重构的功能,对卫星自主导航技术发展有着非常重要的意义。
发明内容
发明目的:本发明的目的是提供一种具有容错功能的近地轨道卫星多传感器自主导航定位方法,该方法能够充分利用卫星上的测量设备提供的导航信息实现自主定轨,并解决了如何用联邦滤波技术对Unscented卡尔曼滤波结果进行信息融合的问题。
技术方案:为达到上述发明目的,本发明设计了一种基于多传感器信息融合的联邦UKF滤波器,该滤波器采用地球卫星的轨道动力学方程作为状态方程,分别以星敏感器+红外地球敏感器、磁强计+雷达高度计、紫外敏感器作为测量器件建立三个子滤波器;用无复位联邦滤波技术对各子滤波器的Unscented卡尔曼滤波结果进行信息融合;用残差检验技术对子滤波器进行故障检测、隔离和重构。本发明可在测量***非线性较强、对***噪声和量测噪声估计不准确的情况下获得较高的定轨精度。各个子滤波器独立工作,相互之间没有干扰,既保证了输出导航参数的连续性和稳定性,也便于工程实现,并可通过主滤波器迅速实现***重构,使因部分导航传感器故障导致子***失效后,组合导航***仍能继续工作,具容错性。大大提高了地球卫星导航***的适用范围。本发明具体步骤如下:
(1)***初始化;
读取轨道数据初值,状态初始条件为:
状态矢量X0=[x0,y0,z0,vx0,vy0,vz0]T,x0,y0,z0,vx0,vy0,vz0,分别为卫星在X,Y,Z三个方向的初始位置和速度;并设***过程噪声w(k)和量测噪声v(k)为加性噪声,Qk和Rk分别为***和量测噪声协方差阵。
(2)建立基于直角坐标系的地球卫星轨道运动学方程;
选取历元(J2000.0)地球赤道惯性坐标系,建立基于直角坐标系的地球卫星轨道运动学方程,即状态方程,通过四阶龙格-库塔方法解微分方程计算出地球卫星的位置(x,y,z)和速度(vx,vy,vz)信息。简化的***状态方程为:
式中,Re为地球参考赤道半径,状态矢量X=[x,y,z,vx,vy,vz]T,x,y,z,vx,vy,vz分别为卫星在X,Y,Z三个方向的位置和速度, 为地心距,μ是地球引力常数,J2为地球引力二阶带谐项系数,ΔFx,ΔFy,ΔFz为地球非球形摄动的高阶摄动项和日、月摄动以及太阳光压摄动和大气摄动等摄动力的影响,在简化模型中这些摄动力的影响用***过程噪声w(k)表示。
对于***过程噪声为零均值白噪声有:
E{w(k)}=0,E{w(k)wT(j)}=Qkδkj,k,j=1,2,3…
式中,δkj为克朗涅克笛耳塔(Kronecker Delta)函数。
(3)建立以星敏感器和红外地球敏感器的输出值为量测量的子***量测方程;
以星敏感器和红外地球敏感器的子***的量测量为星光角距,星光角距是指从卫星上观测到的导航恒星星光的矢量方向与地心矢量方向之间的夹角,星光角距a的量测方程为:
设量测噪声为零均值白噪声:
E{va(k)}=0, k,j=1,2,3…
(4)建立以磁强计和雷达高度计的输出值为量测量的子***量测方程;
磁强计和雷达高度计构成的子***的量测量为地磁场强度矢量和卫星轨道高度,由于通过磁强计对磁场强度的测量即可单独完成定轨功能,雷达高度计起辅助定位、提高***精度的作用,因此本子***对雷达高度计设计了故障检测阈值,当雷达高度计测量值和计算值的差值超出阈值时,将雷达高度计隔离,仅使用磁强计进行定轨,保证了在雷达高度计失效时子***2仍能正常工作。
三轴磁强计的测量值为卫星所在位置的地磁场矢量在卫星本体系中的三个分量,为简化量测模型,取地磁场强度矢量的模作为磁测自主导航观测量,比较此值与国际地磁场模型(IGRF)之间的差值来提供导航信息。量测方程为:
设量测噪声为零均值白噪声:E{vB(k)}=0, k,j=1,2,3…
星载雷达高度计的测量值为卫星到星下点实际海平面的距离信息,由于海平面很接近于大地水准面,所以可认为雷达高度计测得的卫星高度是卫星至星下点处大地水准面的距离。取测量模型为:
其中,Re为地球参考赤道半径,为地心纬度,f为地球参考椭球的扁率,vH(t)为量测噪声。
设量测噪声为零均值白噪声:
E{vH(k)}=0, k,j=1,2,3…
从式5和式6可以看出,在地磁场强度和轨道高度的测量模型中,磁场强度和轨道高度是根据卫星在地球固联坐标系下的位置信息计算的,而卫星的位置矢量一般是在地心惯性坐标系中给出的,因此,需建立地心惯性坐标系和地球固联坐标系之间的转换矩阵。从地心赤道惯性坐标系Si到地球固连坐标系Se的坐标变换矩阵为:
式中:aG为Greenwich赤经。
(5)建立以紫外敏感器的输出值为量测量的子***量测方程;
紫外敏感器的量测值为卫星的位置矢量。紫外敏感器工作在紫外波段,它可以同时观测多个天体目标,提供三轴的航天器姿态信息,定姿精度可达0.05°。另外,通过对地球圆盘的图像处理,可提取星体坐标系中的地心方向和地心距离信息。结合定姿过程得到的姿态信息,可计算出卫星在惯性系中的位置矢量。子***的量测方程为:
设量测噪声为零均值白噪声:
E{vu(k)}=0, k,j=1,2,3…
(6)依照最小偏度采样策略选择Sigma采样点;
①选择0≤W0<1,Sigma权值为:
式中,W为权值,n为状态方程维数,Wi m为均值加权值,Wi c为协方差加权值。
②对应状态为1维情况,迭代初始向量:
③对于输入维数j=2,…n时,迭代公式为:
(7)根据最小方差估计原则建立离散型UKF算法的预测方程和更新方程,各子***分别独立进行Sigma采样点计算以及预测更新和量测更新;
①预测方程:
对所生成的Sigma点加入状态矢量X的均值和协方差信息:
χi(k+1|k)=f[χi(k|k)] (13)
Zi(k+1|k)=h[χi(k+1|k)] (16)
②更新方程:
P(k+1|k+1)=P(k+1|k)-W(k+1)Pyy(k+1|k)WT(k+1) (22)
(8)利用残差检验法,根据预测滤波残差判断各子滤波器输出值是否有效,如出现故障则将其隔离,否则将滤波结果输入主滤波器进行信息融合;
残差检验法的故障函数为:
其中, 即χm 2(k)为具有m个自由度的χ2变量,m为量测矢量Z的维数;vk为残差,
对于正常滤波,取误警率 其中km(·)为χm 2的概率密度函数,a为显著性水平,TD为由误警率确定的门限值。
(9)建立基于UKF算法的无复位联邦UKF滤波方程;
联邦滤波器的信息融合方式采用无复位模式,仅对状态量进行融合和反馈,各子***的协方差阵独立进行时间更新和预测更新,不参加主滤波器的融合:
(10)按照上述步骤(1)-(9),输出为地球卫星状态矢量估计值及其方差矩阵Pg,其中状态估计值包括地球卫星位置和速度矢量[x,y,z,vx,vy,vz]T,状态估计方差矩阵Ps包括地球卫星位置和速度估计方差[Px g,Py g,Pz g,Pvx g,Pvy g,Pvz g]T。
3、有益效果:本发明的方法有如下优点:(1)采用星敏感器+红外地球敏感器、磁强计+雷达高度计、紫外敏感器的组合,能够充分利用星上的测量设备提供的信息进行导航定位,具有完全的自主性;(2)子***2采用基于UKF技术的磁强计和雷达高度计组合定轨,并设计故障检测阈值,保证在雷达高度计失效后***仍然能继续工作,***具较高的精度和可靠性;(3)采用UKF算法对子滤波器进行状态估计,在测量***非线性较强、或对***噪声和量测噪声估计不准确的情况下均能够获得较高的定轨精度;(4)UKF算法采用最小偏度采样策略,减少了采样点数,在保证精度的前提下提高了计算效率;(5)用无复位联邦滤波技术对各子滤波器的UKF滤波结果进行信息融合,对每个子***单独计算Sigma采样点,各子滤波器独立工作,相互之间没有干扰,既保证了输出导航参数的连续性和稳定性,也便于工程实现,(6)用残差检验技术对子滤波器进行故障检测,并可在检测到故障时,通过主滤波器将其隔离,迅速实现***重构,使因部分导航传感器故障导致子***失效后,组合导航***仍能继续工作,具容错性。
附图说明:
图1卫星多传感器容错联邦UKF***流程图;
图2卫星多传感器容错联邦UKF***结构原理图;
图3标称轨迹示意图;
图4仿真轨迹示意图;
图5标称轨迹位置变化曲线;
图6标称轨迹速度变化曲线;
图7仿真轨迹位置变化曲线;
图8仿真轨迹速度变化曲线;
图9理想情况下联邦UKF滤波误差曲线;
图10对***噪声估计不足情况下的联邦UKF滤波误差曲线;
图11对量测噪声估计不足情况下的联邦UKF滤波误差曲线;
图12子***1故障情况下的联邦UKF误差曲线;
图13子***2故障情况下的联邦UKF误差曲线;
图14子***3故障情况下的联邦UKF误差曲线;
具体实施方式
本实例的具体组成图如图1所示,实施的流程如图2所示,轨道数据由STK(卫星仿真工具包)软件产生,具体仿真条件如下;
①仿真时间2008年6月1日12:00-2008年6月1日24:00,采样周期T为10s。标称轨道数据为:半长轴:a=6823.082km;偏心率:e=1.83×10-16;轨道倾角:i=97.2°;升交点赤经:Ω=157.5°;近升角距:ω=0°。
②为了模拟卫星在轨运行的真实运动情况,在计算卫星轨道时,使用力学模型JGM-3,并考虑如下摄动因素:地球非球形摄动;太阳引力、月球引力;太阳光压;大气阻力。
③测量仪器精度:导航星采用分布在天球上的80颗最亮的恒星(星等≤2.3),星敏感器的视场是25°×25°,精度为3″;红外地平仪的精度为0.02°;磁强计测量噪声误差取为50nT,地磁场模型误差取为150nT;雷达高度计测量误差取为200m。紫外敏感器测距误差500m,地心方向测量误差0.01°。
④轨道初始位置误差为5km,初始速度误差取为10m/s。
本实例以3个子***的量测输出值作为量测量,结合地球卫星的状态方程,分别利用基于最小偏度采样的UKF算法对各子***单独进行时间更新和量测更新,并利用残差检验法对3个子***的滤波结果进行分析,判断子***输出是否有效,如果有效则将该结果输出给主滤波器,利用无复位联邦滤波器在主滤波器中对有效的子滤波器的状态估计值进行信息融合,得到地球卫星的位置、速度信息。
具体的实施步骤如下:
(1)***初始化;
读取轨道数据初值,状态初始条件为:
状态矢量X0=[x0,y0,z0,vx0,vy0,vz0]T,x0,y0,z0,vx0,vy0,vz0,分别为卫星在X,Y,Z三个方向的初始位置和速度;并设***过程噪声w(k)和量测噪声v(k)为加性噪声,Qk和Rk分别为***和量测噪声协方差阵。
(2)建立基于直角坐标系的地球卫星轨道运动学方程;
选取历元(J2000.0)地球赤道惯性坐标系,建立基于直角坐标系的地球卫星轨道运动学方程,即状态方程,通过四阶龙格-库塔方法解微分方程计算出地球卫星的位置(x,y,z)和速度(vx,vy,vz)信息。简化的***状态方程为:
式中,Re为地球参考赤道半径,状态矢量X=[x,y,z,vx,vy,vz]T,x,y,z,vx,vy,vz分别为卫星在X,Y,Z三个方向的位置和速度, 为地心距,μ是地球引力常数,J2为地球引力二阶带谐项系数,ΔFx,ΔFy,ΔFz为地球非球形摄动的高阶摄动项和日、月摄动以及太阳光压摄动和大气摄动等摄动力的影响,在简化模型中这些摄动力的影响用***过程噪声w(k)表示。
对于***过程噪声为零均值白噪声有:
E{w(k)}=0,E{w(k)wT(j)}=Qkδkj,k,j=1,2,3…式中,δkj为克朗涅克笛耳塔(Kronecker Delta)函数。
(3)建立以星敏感器和红外地球敏感器的输出值为量测量的子***量测方程;
以星敏感器和红外地球敏感器的子***的量测量为星光角距,星光角距是指从卫星上观测到的导航恒星星光的矢量方向与地心矢量方向之间的夹角,星光角距a的量测方程为:
设量测噪声为零均值白噪声:
E{va(k)}=0, k,j=1,2,3…
(4)建立以磁强计和雷达高度计的输出值为量测量的子***量测方程;
磁强计和雷达高度计构成的子***的量测量为地磁场强度矢量和卫星轨道高度,由于通过磁强计对磁场强度的测量即可单独完成定轨功能,雷达高度计起辅助定位、提高***精度的作用,因此本子***对雷达高度计设计了故障检测阈值,当雷达高度计测量值和计算值的差值超出阈值时,将雷达高度计隔离,仅使用磁强计进行定轨,保证了在雷达高度计失效时子***2仍能正常工作。三轴磁强计的测量值为卫星所在位置的地磁场矢量在卫星本体系中的三个分量,为简化量测模型,取地磁场强度矢量的模作为磁测自主导航观测量,比较此值与国际地磁场模型(IGRF)之间的差值来提供导航信息。量测方程为:
设量测噪声为零均值白噪声:
E{vB(k)}=0, k,j=1,2,3…
星载雷达高度计的测量值为卫星到星下点实际海平面的距离信息,由于海平面很接近于大地水准面,所以可认为雷达高度计测得的卫星高度是卫星至星下点处大地水准面的距离。取测量模型为:
其中,Re为地球参考赤道半径,为地心纬度,f为地球参考椭球的扁率,vH(t)为量测噪声。
设量测噪声为零均值白噪声:
E{vH(k)}=0, k,j=1,2,3…
从式5和式6可以看出,在地磁场强度和轨道高度的测量模型中,磁场强度和轨道高度是根据卫星在地球固联坐标系下的位置信息计算的,而卫星的位置矢量一般是在地心惯性坐标系中给出的,因此,需建立地心惯性坐标系和地球固联坐标系之间的转换矩阵。从地心赤道惯性坐标系Si到地球固连坐标系Se的坐标变换矩阵为:
式中:aG为Greenwich赤经。
(5)建立以紫外敏感器的输出值为量测量的子***量测方程;
紫外敏感器的量测值为卫星的位置矢量。紫外敏感器工作在紫外波段,它可以同时观测多个天体目标,提供三轴的航天器姿态信息,定姿精度可达0.05°。另外,通过对地球圆盘的图像处理,可提取星体坐标系中的地心方向和地心距离信息。结合定姿过程得到的姿态信息,可计算出卫星在惯性系中的位置矢量。子***的量测方程为:
设量测噪声为零均值白噪声:
E{vu(k)}=0, k,j=1,2,3…
(6)依照最小偏度采样策略选择Sigma采样点;
①选择0≤W0<1,Sigma权值为:
式中,W为权值,n为状态方程维数,Wi m为均值加权值,Wi c为协方差加权值。
②对应状态为1维情况,迭代初始向量:
③对于输入维数j=2,…n时,迭代公式为:
(7)根据最小方差估计原则建立离散型UKF算法的预测方程和更新方程,各子***分别独立进行Sigma采样点计算以及预测更新和量测更新;
①预测方程:
对所生成的Sigma点加入状态矢量X的均值和协方差信息:
χi(k+1|k)=f[χi(k/k)] (13)
Zi(k+1|k)=h[χi(k+1)] (16)
②更新方程:
P(k+1|k+1)=P(k+1|k)-W(k+1)Pyy(k+1|k)WT(k+1) (22)
(8)利用残差检验法,根据预测滤波残差判断各子滤波器输出值是否有效,如出现故障则将其隔离,否则将滤波结果输入主滤波器进行信息融合;
残差检验法的故障函数为:
其中, 即χm 2(k)为具有m个自由度的χ2变量,m为量测矢量Z的维数;vk为残差,
对于正常滤波,取误警率 其中km(·)为χm 2的概率密度函数,a为显著性水平,TD为由误警率确定的门限值。
(9)建立基于UKF算法的无复位联邦UKF滤波方程;
联邦滤波器的信息融合方式采用无复位模式,仅对状态量进行融合和反馈,各子***的协方差阵独立进行时间更新和预测更新,不参加主滤波器的融合:
(10)按照上述步骤(1)-(9),输出为地球卫星状态矢量估计值及其方差矩阵Pg,其中状态估计值包括地球卫星位置和速度矢量[x,y,z,vx,vy,vz]T,状态估计方差矩阵Pg包括地球卫星位置和速度估计方差[Px g,Py g,Pz g,Pvx g,Pvy g,Pvz g]T。
对本发明的有益效果说明如下:
(1)本发明的***组成图如图1所示,程序流程图如图2所示;
(2)设计轨道为太阳同步轨道,轨道高度500km,降交点地方时上午10:30,轨道倾角97.2°,仿真轨迹如图3-8所示;
(3)为充分说明本发明方法对导航结果的修正作用,分别给出理想状态下的联邦UKF、对***噪声估计不足情况下的联邦UKF、对量测噪声估计不足情况下的联邦UKF滤波的滤波结果,分别如下所示:
①理想情况下联邦UKF滤波误差曲线如图9所示;
②对***噪声估计不足情况下的联邦UKF滤波误差曲线如图10所示;
③对量测噪声估计不足情况下的联邦UKF滤波误差曲线如图11所示;
分别对比3种情况下的联邦UKF滤波结果,可以看出:
①在理想情况下,组合导航***能够迅速收敛,且没有发散现象;在***收敛后,联邦自适应UKF的位置误差不超过400m,速度误差不超过0.8m/s,能够很好的完成定轨功能。
②当对***噪声估计不足时,联邦UKF滤波结果所受影响不大,仍能很好的完成定轨功能。
③当对量测噪声估计不足时,联邦UKF滤波结果所受影响不大,仍能很好的完成定轨功能。
④三种情况下的滤波精度如表1所示。
表1不同情况下的联邦UKF在定轨精度比较
滤波结果(1σ) | 理想状态 | ***噪声估计不足 | 量测噪声估计不足 |
位置精度m速度精度m/s | 238.80.65 | 257.70.73 | 257.20.73 |
从以上分析可以看出,本发明的方法本文采用的联邦UKF万法对低轨卫星位置和速度的误差的修正效果良好,该方法适合于***非线性较强、***状态噪声阵或量测噪声阵未知的情况。
本发明的采用容错联邦UKF算法来提高非线性***的稳定性和滤波精度,因此这里对这种手段所起的作用进行说明。
表2子***故障情况下定轨精度比较
滤波结果(1σ) 无故障 子***1故障 子***2故障 子***3故障
位置精度m 238.8 263.9 264.5 273.3
速度精度m/s 0.65 0.69 0.85 0.81
子***故障情况下联邦UKF滤波结果如图12-14和表2所示,分别在不同时间对3个子***加入突变故障,可以看出,当单个子***故障时,主***可通过残差检验法迅速检验出故障,并进行***重构,子***的故障不会对主***造成污染,滤波结果变化较小。
本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (1)
1.一种基于联邦UKF算法的低轨卫星多传感器容错自主导航方法,其特征在于,包括下列步骤:
(一)***初始化;
读取轨道数据初值,状态初始条件为:
状态矢量X0=[x0,y0,z0,vx0,vy0,vz0]T,x0,y0,z0,vx0,vy0,vz0,分别为卫星在X,Y,Z三个方向的初始位置和速度;并设***过程噪声w(k)和量测噪声v(k)为加性噪声,Qk和Rk分别为***和量测噪声协方差阵;
(二)建立基于直角坐标系的地球卫星轨道运动学方程;
选取历元(J2000.0)地球赤道惯性坐标系,建立基于直角坐标系的地球卫星轨道运动学方程,即状态方程,通过四阶龙格-库塔方法解微分方程计算出地球卫星的位置(x,y,z)和速度(vx,vy,vz)信息,简化的***状态方程为:
式中,Re为地球参考赤道半径,状态矢量X=[x,y,z,vx,vy,vz]T,x,y,z,vx,vy,vz分别为卫星在X,Y,Z三个方向的位置和速度, 为地心距,μ是地球引力常数,J2为地球引力二阶带谐项系数,ΔFx,ΔFy,ΔFz为地球非球形摄动的高阶摄动项和日、月摄动以及太阳光压摄动和大气摄动等摄动力的影响,在简化模型中这些摄动力的影响用***过程噪声w(k)表示;
对于***过程噪声为零均值白噪声有:
E{w(k)}=0,E{w(k)wT(j)}=Qkδkj,k,j=1,2,3…;
式中,δkj为克朗涅克笛耳塔(Kronecker Delta)函数;
(三) 建立以星敏感器和红外地球敏感器的输出值为量测量的子***量测方程;
以星敏感器和红外地球敏感器的子***的量测量为星光角距,星光角距是指从卫星上观测到的导航恒星星光的矢量方向与地心矢量方向之间的夹角,星光角距a的量测方程为:
设量测噪声为零均值白噪声:
E{va(k)}=0, k,j=1,2,3…;
(四)建立以磁强计和雷达高度计的输出值为量测量的子***量测方程;三轴磁强计的测量值为卫星所在位置的地磁场矢量在卫星本体系中的三个分量,为简化量测模型,取地磁场强度矢量的模作为磁测自主导航观测量,比较此值与国际地磁场模型(IGRF)之间的差值来提供导航信息,量测方程为:
设量测噪声为零均值白噪声:
E{vB(k)}=0, k,j=1,2,3…;
星载雷达高度计的测量值为卫星到星下点实际海平面的距离信息,取测量模型为:
其中,Re为地球参考赤道半径,为地心纬度,f为地球参考椭球的扁率;vH(t)为量测噪声;
设量测噪声为零均值白噪声:
E{vH(k)}=0, k,j=1,2,3…;
从式5和式6看出,需建立地心惯性坐标系和地球固联坐标系之间的转换矩阵,从地心赤道惯性坐标系Si到地球固连坐标系Se的坐标变换矩阵为:
式中:aG为Greenwich赤经;
(五)建立以紫外敏感器的输出值为量测量的子***量测方程;
紫外敏感器的量测值为卫星的位置矢量,紫外敏感器工作在紫外波段,同时观测多个天体目标,提供三轴的航天器姿态信息,定姿精度可达0.05°;另外,通过对地球圆盘的图像处理,提取星体坐标系中的地心方向和地心距离信息;结合定姿过程得到的姿态信息,计算出卫星在惯性系中的位置矢量,子***的量测方程为:
设量测噪声为零均值白噪声:
E{vu(k)}=0, k,j=1,2,3…;
(六) 依照最小偏度采样策略选择Sigma采样点;
①选择0≤W0<1,Sigma权值为:
式中,W为权值,n为状态方程维数,Wi m为均值加权值,Wi c为协方差加权值;
②对应状态为1维情况,迭代初始向量:
③对于输入维数j=2,…n时,迭代公式为:
(七)根据最小方差估计原则建立离散型UKF算法的预测方程和更新方程,各子***分别独立进行Sigma采样点计算以及预测更新和量测更新;
①预测方程:
对所生成的Sigma点加入状态矢量X的均值和协方差信息:
χi(k+1|k)=f[χi(k|k)] (13)
Zi(k+1|k)=h[χi(k+1|k)] (16)
②更新方程:
P(k+1|k+1)=P(k+1|k)-W(k+1)Pyy(k+1|k)WT(k+1) (22)
(八)利用残差检验法,根据预测滤波残差判断各子滤波器输出值是否有效,如出现故障则将其隔离,否则将滤波结果输入主滤波器进行信息融合;残差检验法的故障函数为:
其中, 即χm 2(k)为具有m个自由度的χ2变量,m为量测矢量Z的维数;vk为残差,
对于正常滤波,取误警率 其中km(·)为χm 2的概率密度函数,a为显著性水平,TD为由误警率确定的门限值;
(九)建立基于UKF算法的无复位联邦UKF滤波方程;
联邦滤波器的信息融合方式采用无复位模式,仅对状态量进行融合和反馈,各子***的协方差阵独立进行时间更新和预测更新,不参加主滤波器的融合:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100191018A CN101216319B (zh) | 2008-01-11 | 2008-01-11 | 基于联邦ukf算法的低轨卫星多传感器容错自主导航方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100191018A CN101216319B (zh) | 2008-01-11 | 2008-01-11 | 基于联邦ukf算法的低轨卫星多传感器容错自主导航方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101216319A true CN101216319A (zh) | 2008-07-09 |
CN101216319B CN101216319B (zh) | 2012-01-11 |
Family
ID=39622799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008100191018A Expired - Fee Related CN101216319B (zh) | 2008-01-11 | 2008-01-11 | 基于联邦ukf算法的低轨卫星多传感器容错自主导航方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101216319B (zh) |
Cited By (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101520328B (zh) * | 2009-04-01 | 2010-08-25 | 西北工业大学 | 地磁场线图匹配自主导航方法 |
CN101608921B (zh) * | 2009-07-21 | 2011-08-24 | 华中科技大学 | 一种脉冲星/cns组合导航方法 |
CN102175260A (zh) * | 2010-12-31 | 2011-09-07 | 北京控制工程研究所 | 一种自主导航***误差校正方法 |
CN102288200A (zh) * | 2011-07-06 | 2011-12-21 | 清华大学 | 用于星敏感器的精度测量*** |
CN102597701A (zh) * | 2009-10-15 | 2012-07-18 | 纳夫科姆技术公司 | 用于补偿错误测量的***和方法 |
CN101556155B (zh) * | 2009-05-20 | 2012-07-18 | 上海微小卫星工程中心 | 小卫星姿态确定***及方法 |
CN102608631A (zh) * | 2011-10-28 | 2012-07-25 | 北京航空航天大学 | 基于模糊逻辑的自适应强跟踪ukf定位滤波算法 |
CN102654582A (zh) * | 2012-04-16 | 2012-09-05 | 东莞市泰斗微电子科技有限公司 | 一种联合导航***及方法 |
CN102735265A (zh) * | 2012-06-18 | 2012-10-17 | 北京控制工程研究所 | 一种基于陀螺漂移估计值的星敏感器周期性故障检测方法 |
CN102735238A (zh) * | 2012-06-18 | 2012-10-17 | 北京控制工程研究所 | 一种基于相关性模型的红外地球敏感器测点优化配置方法 |
CN102914308A (zh) * | 2012-10-24 | 2013-02-06 | 南京航空航天大学 | 一种基于新息正交性的抗野值联邦滤波方法 |
CN102927987A (zh) * | 2012-10-24 | 2013-02-13 | 北京控制工程研究所 | 一种地球敏感器探头在轨故障情况下的处理方法 |
CN103309348A (zh) * | 2013-06-28 | 2013-09-18 | 哈尔滨工业大学 | 一种利用二阶卡尔曼滤波算法估计卫星姿态控制***执行机构加性故障大小的方法 |
CN103411627A (zh) * | 2013-08-07 | 2013-11-27 | 北京航空航天大学 | 火星动力下降段非线性三步滤波方法 |
CN103759733A (zh) * | 2013-10-29 | 2014-04-30 | 南昌大学 | 基于联邦滤波的ins/vkm/vdm车载导航*** |
CN103776449A (zh) * | 2014-02-26 | 2014-05-07 | 北京空间飞行器总体设计部 | 一种提高鲁棒性的动基座初始对准方法 |
CN103810382A (zh) * | 2014-01-27 | 2014-05-21 | 中国电子科技集团公司第十研究所 | 机载分布式多传感器两级数据融合策略选择方法 |
CN104020671A (zh) * | 2014-05-30 | 2014-09-03 | 哈尔滨工程大学 | 一种量测干扰下用于飞行器姿态估计的鲁棒递推滤波方法 |
CN105139406A (zh) * | 2015-09-08 | 2015-12-09 | 哈尔滨工业大学 | 一种基于序列图像的跟踪精度反演方法 |
CN106597498A (zh) * | 2017-01-18 | 2017-04-26 | 哈尔滨工业大学 | 多传感器融合***空时偏差联合校准方法 |
CN106646543A (zh) * | 2016-12-22 | 2017-05-10 | 成都正扬博创电子技术有限公司 | 基于主从式aukf算法的高动态卫星导航信号载波追踪方法 |
CN106643666A (zh) * | 2016-09-13 | 2017-05-10 | 首都师范大学 | 基于同平台多个传感器的运动目标探测方法和*** |
CN106767998A (zh) * | 2017-01-05 | 2017-05-31 | 清华大学 | 空间信息网络中多卫星协同的探测方法和装置 |
CN107421541A (zh) * | 2017-05-25 | 2017-12-01 | 西北工业大学 | 一种容错非接触式失效卫星的形态参数测算方法 |
CN107728123A (zh) * | 2017-09-08 | 2018-02-23 | 中国电子科技集团公司信息科学研究院 | 雷达极‑直坐标转换观测精度分析方法、装置和*** |
CN108469622A (zh) * | 2018-06-22 | 2018-08-31 | 北京电子工程总体研究所 | 基于自适应阈值的卫星导航数据完好性检测***及方法 |
CN108981708A (zh) * | 2018-08-02 | 2018-12-11 | 南京航空航天大学 | 四旋翼扭矩模型/航向陀螺/磁传感器容错组合导航方法 |
CN109211270A (zh) * | 2018-08-17 | 2019-01-15 | 中国航空工业集团公司西安飞行自动控制研究所 | 一种惯性天文卫星组合导航装置的故障检测*** |
CN109471143A (zh) * | 2018-12-11 | 2019-03-15 | 北京交通大学 | 自适应容错的列车组合定位方法 |
CN110017837A (zh) * | 2019-04-26 | 2019-07-16 | 沈阳航空航天大学 | 一种姿态抗磁干扰的组合导航方法 |
CN110196068A (zh) * | 2019-05-27 | 2019-09-03 | 哈尔滨工程大学 | 一种极区集中滤波组合导航***残差向量故障检测与隔离方法 |
CN110285815A (zh) * | 2019-05-28 | 2019-09-27 | 山东航天电子技术研究所 | 一种可在轨全程应用的微纳卫星多源信息姿态确定方法 |
CN110502023A (zh) * | 2019-07-18 | 2019-11-26 | 南京航空航天大学 | 一种基于分布式智能敏感器的航天器姿态确定实现方法 |
CN110779532A (zh) * | 2019-11-18 | 2020-02-11 | 河南工业大学 | 一种应用于近地轨道卫星的地磁导航***及方法 |
CN110865650A (zh) * | 2019-11-19 | 2020-03-06 | 武汉工程大学 | 基于主动视觉的无人机位姿自适应估计方法 |
CN111339629A (zh) * | 2019-11-22 | 2020-06-26 | 北京理工大学 | 一种用于天基观测的空间目标机动轨道确定方法 |
CN111811512A (zh) * | 2020-06-02 | 2020-10-23 | 北京航空航天大学 | 基于联邦平滑的mpos离线组合估计方法和装置 |
CN112212860A (zh) * | 2020-08-28 | 2021-01-12 | 山东航天电子技术研究所 | 具有故障容错的分布式滤波微纳卫星姿态确定方法 |
CN117315921A (zh) * | 2023-10-25 | 2023-12-29 | 招商局检测车辆技术研究院有限公司 | 一种基于融合感知的斑马线预警方法及*** |
-
2008
- 2008-01-11 CN CN2008100191018A patent/CN101216319B/zh not_active Expired - Fee Related
Cited By (64)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101520328B (zh) * | 2009-04-01 | 2010-08-25 | 西北工业大学 | 地磁场线图匹配自主导航方法 |
CN101556155B (zh) * | 2009-05-20 | 2012-07-18 | 上海微小卫星工程中心 | 小卫星姿态确定***及方法 |
CN101608921B (zh) * | 2009-07-21 | 2011-08-24 | 华中科技大学 | 一种脉冲星/cns组合导航方法 |
CN102597701A (zh) * | 2009-10-15 | 2012-07-18 | 纳夫科姆技术公司 | 用于补偿错误测量的***和方法 |
CN102597701B (zh) * | 2009-10-15 | 2015-06-03 | 纳夫科姆技术公司 | 用于补偿错误测量的***和方法 |
CN102175260A (zh) * | 2010-12-31 | 2011-09-07 | 北京控制工程研究所 | 一种自主导航***误差校正方法 |
CN102288200A (zh) * | 2011-07-06 | 2011-12-21 | 清华大学 | 用于星敏感器的精度测量*** |
CN102288200B (zh) * | 2011-07-06 | 2012-07-18 | 清华大学 | 用于星敏感器的精度测量*** |
CN102608631A (zh) * | 2011-10-28 | 2012-07-25 | 北京航空航天大学 | 基于模糊逻辑的自适应强跟踪ukf定位滤波算法 |
CN102654582A (zh) * | 2012-04-16 | 2012-09-05 | 东莞市泰斗微电子科技有限公司 | 一种联合导航***及方法 |
CN102735238B (zh) * | 2012-06-18 | 2014-08-27 | 北京控制工程研究所 | 一种基于相关性模型的红外地球敏感器测点优化配置方法 |
CN102735238A (zh) * | 2012-06-18 | 2012-10-17 | 北京控制工程研究所 | 一种基于相关性模型的红外地球敏感器测点优化配置方法 |
CN102735265A (zh) * | 2012-06-18 | 2012-10-17 | 北京控制工程研究所 | 一种基于陀螺漂移估计值的星敏感器周期性故障检测方法 |
CN102735265B (zh) * | 2012-06-18 | 2014-12-17 | 北京控制工程研究所 | 一种基于陀螺漂移估计值的星敏感器周期性故障检测方法 |
CN102914308B (zh) * | 2012-10-24 | 2014-12-10 | 南京航空航天大学 | 一种基于新息正交性的抗野值联邦滤波方法 |
CN102927987A (zh) * | 2012-10-24 | 2013-02-13 | 北京控制工程研究所 | 一种地球敏感器探头在轨故障情况下的处理方法 |
CN102914308A (zh) * | 2012-10-24 | 2013-02-06 | 南京航空航天大学 | 一种基于新息正交性的抗野值联邦滤波方法 |
CN102927987B (zh) * | 2012-10-24 | 2015-07-08 | 北京控制工程研究所 | 一种地球敏感器探头在轨故障情况下的处理方法 |
CN103309348A (zh) * | 2013-06-28 | 2013-09-18 | 哈尔滨工业大学 | 一种利用二阶卡尔曼滤波算法估计卫星姿态控制***执行机构加性故障大小的方法 |
CN103411627A (zh) * | 2013-08-07 | 2013-11-27 | 北京航空航天大学 | 火星动力下降段非线性三步滤波方法 |
CN103411627B (zh) * | 2013-08-07 | 2016-12-07 | 北京航空航天大学 | 火星动力下降段非线性三步滤波方法 |
CN103759733B (zh) * | 2013-10-29 | 2016-10-12 | 南昌大学 | 基于联邦滤波的ins/vkm/vdm车载导航*** |
CN103759733A (zh) * | 2013-10-29 | 2014-04-30 | 南昌大学 | 基于联邦滤波的ins/vkm/vdm车载导航*** |
CN103810382B (zh) * | 2014-01-27 | 2017-01-25 | 中国电子科技集团公司第十研究所 | 机载分布式多传感器两级数据融合策略选择方法 |
CN103810382A (zh) * | 2014-01-27 | 2014-05-21 | 中国电子科技集团公司第十研究所 | 机载分布式多传感器两级数据融合策略选择方法 |
CN103776449A (zh) * | 2014-02-26 | 2014-05-07 | 北京空间飞行器总体设计部 | 一种提高鲁棒性的动基座初始对准方法 |
CN103776449B (zh) * | 2014-02-26 | 2016-05-04 | 北京空间飞行器总体设计部 | 一种提高鲁棒性的动基座初始对准方法 |
CN104020671B (zh) * | 2014-05-30 | 2017-01-11 | 哈尔滨工程大学 | 一种量测干扰下用于飞行器姿态估计的鲁棒递推滤波方法 |
CN104020671A (zh) * | 2014-05-30 | 2014-09-03 | 哈尔滨工程大学 | 一种量测干扰下用于飞行器姿态估计的鲁棒递推滤波方法 |
CN105139406A (zh) * | 2015-09-08 | 2015-12-09 | 哈尔滨工业大学 | 一种基于序列图像的跟踪精度反演方法 |
CN105139406B (zh) * | 2015-09-08 | 2018-02-23 | 哈尔滨工业大学 | 一种基于序列图像的跟踪精度反演方法 |
US10970857B2 (en) | 2016-09-13 | 2021-04-06 | Capital Normal University | Moving object detection method and system based on multiple sensors on same platform |
CN106643666A (zh) * | 2016-09-13 | 2017-05-10 | 首都师范大学 | 基于同平台多个传感器的运动目标探测方法和*** |
CN106643666B (zh) * | 2016-09-13 | 2019-12-06 | 首都师范大学 | 基于同平台多个传感器的运动目标探测方法和*** |
WO2018049740A1 (zh) * | 2016-09-13 | 2018-03-22 | 首都师范大学 | 基于同平台多个传感器的运动目标探测方法和*** |
CN106646543A (zh) * | 2016-12-22 | 2017-05-10 | 成都正扬博创电子技术有限公司 | 基于主从式aukf算法的高动态卫星导航信号载波追踪方法 |
CN106767998B (zh) * | 2017-01-05 | 2019-06-04 | 清华大学 | 空间信息网络中多卫星协同的探测方法和装置 |
CN106767998A (zh) * | 2017-01-05 | 2017-05-31 | 清华大学 | 空间信息网络中多卫星协同的探测方法和装置 |
CN106597498A (zh) * | 2017-01-18 | 2017-04-26 | 哈尔滨工业大学 | 多传感器融合***空时偏差联合校准方法 |
CN107421541B (zh) * | 2017-05-25 | 2020-05-12 | 西北工业大学 | 一种容错非接触式失效卫星的形态参数测算方法 |
CN107421541A (zh) * | 2017-05-25 | 2017-12-01 | 西北工业大学 | 一种容错非接触式失效卫星的形态参数测算方法 |
CN107728123A (zh) * | 2017-09-08 | 2018-02-23 | 中国电子科技集团公司信息科学研究院 | 雷达极‑直坐标转换观测精度分析方法、装置和*** |
CN108469622A (zh) * | 2018-06-22 | 2018-08-31 | 北京电子工程总体研究所 | 基于自适应阈值的卫星导航数据完好性检测***及方法 |
CN108981708A (zh) * | 2018-08-02 | 2018-12-11 | 南京航空航天大学 | 四旋翼扭矩模型/航向陀螺/磁传感器容错组合导航方法 |
CN109211270A (zh) * | 2018-08-17 | 2019-01-15 | 中国航空工业集团公司西安飞行自动控制研究所 | 一种惯性天文卫星组合导航装置的故障检测*** |
CN109211270B (zh) * | 2018-08-17 | 2022-03-15 | 中国航空工业集团公司西安飞行自动控制研究所 | 一种惯性天文卫星组合导航装置的故障检测*** |
CN109471143A (zh) * | 2018-12-11 | 2019-03-15 | 北京交通大学 | 自适应容错的列车组合定位方法 |
CN109471143B (zh) * | 2018-12-11 | 2023-05-05 | 北京交通大学 | 自适应容错的列车组合定位方法 |
CN110017837A (zh) * | 2019-04-26 | 2019-07-16 | 沈阳航空航天大学 | 一种姿态抗磁干扰的组合导航方法 |
CN110017837B (zh) * | 2019-04-26 | 2022-11-25 | 沈阳航空航天大学 | 一种姿态抗磁干扰的组合导航方法 |
CN110196068A (zh) * | 2019-05-27 | 2019-09-03 | 哈尔滨工程大学 | 一种极区集中滤波组合导航***残差向量故障检测与隔离方法 |
CN110285815B (zh) * | 2019-05-28 | 2023-10-20 | 山东航天电子技术研究所 | 一种可在轨全程应用的微纳卫星多源信息姿态确定方法 |
CN110285815A (zh) * | 2019-05-28 | 2019-09-27 | 山东航天电子技术研究所 | 一种可在轨全程应用的微纳卫星多源信息姿态确定方法 |
CN110502023A (zh) * | 2019-07-18 | 2019-11-26 | 南京航空航天大学 | 一种基于分布式智能敏感器的航天器姿态确定实现方法 |
CN110779532B (zh) * | 2019-11-18 | 2023-03-31 | 河南工业大学 | 一种应用于近地轨道卫星的地磁导航***及方法 |
CN110779532A (zh) * | 2019-11-18 | 2020-02-11 | 河南工业大学 | 一种应用于近地轨道卫星的地磁导航***及方法 |
CN110865650A (zh) * | 2019-11-19 | 2020-03-06 | 武汉工程大学 | 基于主动视觉的无人机位姿自适应估计方法 |
CN110865650B (zh) * | 2019-11-19 | 2022-12-20 | 武汉工程大学 | 基于主动视觉的无人机位姿自适应估计方法 |
CN111339629A (zh) * | 2019-11-22 | 2020-06-26 | 北京理工大学 | 一种用于天基观测的空间目标机动轨道确定方法 |
CN111811512A (zh) * | 2020-06-02 | 2020-10-23 | 北京航空航天大学 | 基于联邦平滑的mpos离线组合估计方法和装置 |
CN111811512B (zh) * | 2020-06-02 | 2023-08-01 | 北京航空航天大学 | 基于联邦平滑的mpos离线组合估计方法和装置 |
CN112212860B (zh) * | 2020-08-28 | 2023-03-03 | 山东航天电子技术研究所 | 具有故障容错的分布式滤波微纳卫星姿态确定方法 |
CN112212860A (zh) * | 2020-08-28 | 2021-01-12 | 山东航天电子技术研究所 | 具有故障容错的分布式滤波微纳卫星姿态确定方法 |
CN117315921A (zh) * | 2023-10-25 | 2023-12-29 | 招商局检测车辆技术研究院有限公司 | 一种基于融合感知的斑马线预警方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN101216319B (zh) | 2012-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101216319A (zh) | 基于联邦ukf算法的低轨卫星多传感器容错自主导航方法 | |
CN109556632B (zh) | 一种基于卡尔曼滤波的ins/gnss/偏振/地磁组合导航对准方法 | |
CN107065025B (zh) | 一种基于重力场梯度不变量的轨道要素估计方法 | |
CN100575877C (zh) | 基于多信息融合的航天器组合导航方法 | |
CN100501331C (zh) | 基于x射线脉冲星的导航卫星自主导航***与方法 | |
CN106338753A (zh) | 一种基于地面站/星间链路/gnss联合测量的地球同步轨道星座定轨方法 | |
CN1987355A (zh) | 一种基于自适应扩展卡尔曼滤波的地球卫星自主天文导航方法 | |
CN101246011B (zh) | 一种基于凸优化算法的多目标多传感器信息融合方法 | |
CN104406605A (zh) | 机载多导航源综合导航仿真*** | |
CN103018783B (zh) | 重力卫星编队轨道稳定性优化设计和精密反演地球重力场方法 | |
CN104848862B (zh) | 一种环火探测器精密同步定位守时方法及*** | |
CN103868514A (zh) | 一种在轨飞行器自主导航*** | |
CN112161632B (zh) | 一种基于相对位置矢量测量的卫星编队初始定位方法 | |
CN105737858A (zh) | 一种机载惯导***姿态参数校准方法与装置 | |
CN106643744A (zh) | 一种基于四程中继跟踪模式的远月面着陆器精密定位方法 | |
CN1987356A (zh) | 一种航天器天文/多普勒组合导航方法 | |
CN103076639B (zh) | 基于残余星间速度原理反演地球重力场的方法 | |
Xinlong et al. | An autonomous navigation scheme based on geomagnetic and starlight for small satellites | |
Roh et al. | Orbit determination using the geomagnetic field measurement via the unscented Kalman filter | |
Jiancheng et al. | Installation direction analysis of star sensors by hybrid condition number | |
CN101608921B (zh) | 一种脉冲星/cns组合导航方法 | |
Gao et al. | Tightly coupled INS/CNS/spectral redshift integrated navigation system with the aid of redshift error measurement | |
Matsumoto et al. | Impact of far-side satellite tracking on gravity estimation in the SELENE project | |
Vigrahala et al. | Attitude, Position and Velocity determination using Low-cost Inertial Measurement Unit for Global Navigation Satellite System Outages | |
CN112526624A (zh) | 重力卫星东西方向差分观测数据构建及反演方法和*** |
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: 20120111 Termination date: 20150111 |
|
EXPY | Termination of patent right or utility model |