CN104034334B - 一种小视场星敏感器的单星及双星定姿方法 - Google Patents
一种小视场星敏感器的单星及双星定姿方法 Download PDFInfo
- Publication number
- CN104034334B CN104034334B CN201410246926.9A CN201410246926A CN104034334B CN 104034334 B CN104034334 B CN 104034334B CN 201410246926 A CN201410246926 A CN 201410246926A CN 104034334 B CN104034334 B CN 104034334B
- Authority
- CN
- China
- Prior art keywords
- star
- attitude
- carrier
- coordinate system
- sensor
- 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 40
- 239000013598 vector Substances 0.000 claims abstract description 25
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 23
- 230000000007 visual effect Effects 0.000 claims abstract description 18
- 238000001914 filtration Methods 0.000 claims abstract description 8
- 230000009191 jumping Effects 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 20
- 241000612151 Dodecatheon Species 0.000 claims description 19
- 235000008166 Dodecatheon meadia Nutrition 0.000 claims description 19
- 238000005259 measurement Methods 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 8
- 238000009434 installation Methods 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000012795 verification Methods 0.000 claims description 2
- 210000004087 cornea Anatomy 0.000 abstract 2
- 238000004364 calculation method Methods 0.000 description 5
- 238000000605 extraction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
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/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
本发明提供的是一种小视场星敏感器的单星及双星定姿方法。一:判断视场内恒星的数目,若视场内恒星数目不小于3,利用三角形算法进行识别,若识别成功,则计算并输出姿态,否则跳到步骤四;二:若视场内恒星数目为1,利用上一时刻的姿态信息和陀螺输出信息进行星点识别,若识别成功,则结合卡尔曼滤波实时估计当前姿态,否则跳到步骤四;三:若视场内恒星数目为2,利用上一时刻的姿态信息、陀螺输出信息和星点的亮度和位置信息进行星图识别,若识别成功,用双矢量算法计算姿态,否则跳到步骤四;四:星图识别失败后利用陀螺信息估计当前时刻姿态。本发明在视场内只有一、两颗星的情况下仍然能够进行星图识别并计算载体姿态。适用于小视场的星敏感器。
Description
技术领域
本发明涉及的是一种星图识别和姿态确定方法,具体地说是一种小视场星敏感器的单星及双星定姿方法。
背景技术
精确的姿态信息是航天器完成航天任务的基础。能够提供姿态信息的测量器件有很多,但是精度最高应用最广泛的是星敏感器。星敏感器的姿态确定过程分为:星空成像、星点提取、星图识别和姿态计算四个步骤。
星空成像和星点提取属于星图的预处理阶段,通过阈值分割算法可以将星点信息与星图的背景分离开来,利用连通域法或者多阈值聚类的方法将星点目标区分开来,最后使用质心法或者高斯曲面法等质心提取方法将星点的位置和亮度等信息提取出来。星点的位置和亮度等信息是星图识别阶段最重要的信息,星图识别的结果直接影响到姿态解算的精度。随着星敏感器的发展,出现了许多的星图识别算法,例如三角形算法、多边形算法、栅格算法等。然而目前的星图识别算法都要求视场内的星点达到一定的数目,即使是要求最低的三角形算法也至少需要三颗星才能够进行星图识别,对于视场内只有两颗星或者一颗星的情况无能为力,这就限制了这些算法在小视场星敏感器中的应用。但是在小型卫星上,为了降低成本和体积,小视场星敏感器的应用还是很多的,如果使用传统星图识别和姿态计算方法,当视场内只有一颗或者两颗星时不能够进行星图识别,此时就不能实时输出姿态信息。
发明内容
本发明的目的在于提供一种视场内只有一颗或者两颗星的情况下也能够实时估计载体的姿态的小视场星敏感器的单星及双星定姿方法。
本发明的目的是这样实现的:
步骤一:判断视场内恒星的数目,若视场内恒星数目不小于3,利用三角形算法进行识别,若识别成功,则计算并输出姿态,否则跳到步骤四;
步骤二:若视场内恒星数目为1,利用上一时刻的姿态信息和陀螺输出信息进行星点识别,若识别成功,则结合卡尔曼滤波实时估计当前姿态,否则跳到步骤四;
步骤三:若视场内恒星数目为2,利用上一时刻的姿态信息、陀螺输出信息和星点的亮度和位置信息进行星图识别,若识别成功,用双矢量算法计算姿态,否则跳到步骤四;
步骤四:星图识别失败后利用陀螺信息估计当前时刻姿态。
本发明基于卡尔曼滤波对视场内单星及双星定姿技术进行研究。当视场内只有一颗星时,利用上一时刻的姿态信息和陀螺输出信息对其进行识别,识别成功后根据四元数姿态动力学建立状态方程,然后以单星位置信息作为观测量并使用卡尔曼滤波实时估计当前的姿态信息。当视场内有两颗星时,利用上一时刻的姿态信息、陀螺输出信息以及两颗星的角距及亮度等信息进行星图识别,识别成功以后直接使用双矢量定姿的方法计算载体的姿态。
本发明的优点主要体现在:
1、本发明在视场内只有一颗星和两颗星的情况下仍然能够进行星图识别并计算载体姿态,打破了传统算法的限制,提高了姿态计算的精度。
2、本发明适用于小视场的星敏感器,对于小型、低成本卫星的姿态确定有着重要的意义。
附图说明
图1是单星识别示意图;
图2是双星识别示意图;
图3星敏感器视场为6°×6°时航向角的对比曲线;
图4星敏感器视场为4°×4°时航向角的对比曲线;
图5星敏感器视场为3°×3°时航向角的对比曲线;
图6是本发明的流程图。
具体实施方式
下面将结合附图和实施实例对本发明作进一步的详细说明。结合图6本发明具体包括如下步骤:
步骤一:判断视场内恒星的数目,若视场内恒星数目不小于3,利用三角形算法进行识别,若识别成功,则计算并输出姿态,否则跳到步骤四;
步骤二:若视场内恒星数目为1,利用上一时刻的姿态信息和陀螺输出信息进行星点识别,若识别成功,则结合卡尔曼滤波实时估计当前姿态,否则跳到步骤四;
具体是实现步骤为:
步骤A:根据载体上一时刻的姿态信息和陀螺的输出信息计算载体当前时刻的姿态;
步骤B:利用星敏感器与载体之间的安装矩阵计算星敏感器视轴的方向矢量P;
使用下标b代表载体坐标系,上标i代表惯性坐标系,当前时刻由陀螺计算得到的载体系到惯性坐标系的转换矩阵为用下标s代表星敏感器坐标系,则星敏感器坐标系到惯性坐标系的转换矩阵为星敏感器的视轴方向即为星敏感器坐标系的z轴方向,方向矢量P即为的第三列元素;
步骤C:利用P从星表中选择落入视场中的星,并将其投影到拍摄星图中;
一般情况下,生成的投影星图中只有两颗星,拍摄星和投影星。此时如果盲目的将这颗投影星作为拍摄星的对应星,则之后可能输出错误的姿态,这是姿态确定中所不能忍受的。若陀螺和星敏感器都不存在噪声,则拍摄星和其对应的投影星将会重合。陀螺和星敏感器的噪声在星图中体现为星点的位置噪声,从而使得拍摄星和其投影星不重合。假设陀螺按照载体三轴安装,其漂移为1°/h,星敏感器更新周期为1秒,则由陀螺漂移引入的噪声可近似为星敏感器噪声为3",则总的噪声约为而在星图预处理时,星点之间的最小间距一般要大于5个像素,当视场大小为4°×4°,像素为512×512时,此间距大约对应0.04°>>0.0013°,即陀螺漂移和星敏感器噪声产生的星点位置噪声要远小于星图中星点之间的间距,实际上现在的陀螺和星敏感器的精度要高于上述的假设值,那么陀螺漂移和星敏感器噪声产生的星点位置噪声与星图中星点之间的间距差距就更加明显。因此,可以将星点最小间距值d作为验证的基础条件。若只有拍摄星和投影星时,计算两颗星的欧拉距离,若其小于d,则认为投影星为该拍摄星的对应星。由于噪声的原因,偶尔会有星图中出现多颗投影星的情况,如图1所示,此时在投影星图中计算拍摄星和每颗投影星之间的欧拉距离,距离最小的星若满足上述基础条件,则它就是拍摄星的对应星;若星图中没有投影星则此次星图识别失败;
步骤D:利用识别成功星点的位置信息实时估计当前姿态;
具体实现方法为:
根据四元数姿态动力学,建立如下的状态方程
其中,q=[q0 q1 q2 q3]为载体的姿态四元数,M(ω)具体形式如下
其中,ωx、ωyωz为陀螺的测量值;假设拍摄星点在星敏感器坐标系中的方向矢量为ps,其对应的投影星在地心赤道坐标系的方向矢量为pi,建立如下的量测方程
其中,为星敏感器的安装矩阵,C(q)为姿态四元数表示的姿态矩阵,其具体表示如下
根据上述建立的状态方程和量测方程,利用扩展卡尔曼滤波可以实时估计载体当前时刻的姿态。
步骤三:若视场内恒星数目为2,利用上一时刻的姿态信息、陀螺输出信息和星点的亮度和位置等信息进行星图识别,若识别成功,用双矢量算法计算姿态,否则跳到步骤四;
具体实现方法为:
步骤A:根据载体上一时刻的姿态信息和陀螺的输出信息计算载体当前时刻的姿态;
步骤B:利用星敏感器与载体之间的安装矩阵计算星敏感器视轴的方向矢量P;
使用下标b代表载体坐标系,上标i代表惯性坐标系,当前时刻由陀螺计算得到的载体系到惯性坐标系的转换矩阵为用下标s代表星敏感器坐标系,则星敏感器坐标系到惯性坐标系的转换矩阵为星敏感器的视轴方向即为星敏感器坐标系的z轴方向,方向矢量P即为的第三列元素;
步骤C:利用P从星表中选择落入视场中的星,并将其投影到拍摄星图中;
若投影星只有一颗,则计算这颗投影星与两颗拍摄星之间的欧拉距离,欧拉距离小且满足步骤二基础条件的就是其对应的拍摄星,此时利用步骤二的方法估计载体的姿态;若没有投影星,则星图识别失败;若有两颗或多颗投影星,如图2所示,根据单星识别的原理,计算每颗拍摄星与每颗投影星之间的欧拉距离,从中选取欧拉距离最小且满足基础条件的星作为拍摄星的对应星。
步骤D:星图识别成功以后利用双矢量算法计算载体姿态;
若步骤C中将两颗星都识别成功,则可以得到两颗星在惯性坐标系的方向矢量U、V,利用这两个方向矢量可建立三个单位矢量Xn、Yn、Zn,并构成一个新的正交坐标基,即
令Fn=[Xn Yn Zn]T,利用同样的方法可以建立星敏感器坐标系下的坐标基Fs=[Xs Ys Zs]T,则有由于式中Fs、Fn为3×3方阵,因此存在逆矩阵,则载体的姿态矩阵为
利用公式(4)就可以计算出载体的姿态;
步骤四:星图识别失败后利用陀螺信息估计当前时刻姿态;
由于各种噪声的存在,星敏感器不能够完美的呈现星空。因此,为了模拟真实情况下的星图识别过程,用两种形式的高斯白噪声对所有的噪声进行模拟。第一种是星等噪声MN(Magnitude Noise),由于恒星亮度的变化引起。第二种是星点的位置噪声PN(PositionNoise),由星敏感器噪声和陀螺测量噪声引起,测量单位为像素。另外,由于恒星的自行、测量误差等因素的存在,星表中也存在噪声。对于大视场星敏感器而言,与其它噪声相比星表噪声非常小可以忽略不计。然而,对于本文使用的小视场星敏感器,星表噪声是需要考虑的一个因素,因此在本文中星表噪声也用MN和PN来表征。
仿真中,使用的星敏感器的精度为3",星等敏感极限为6.5Mv,以卫星为载体对比传统的三角形算法和本发明方法,图3-图5为以航向角为例视场大小为6°×6°,4°×4°,3°×3°的对比曲线。
如图3所示,当视场大小为6°×6°时两种方法的姿态测量精度相当。这是因为当视场大小为6°×6°时,星敏感器能够捕获的星点数目一般多大于3颗,此时本发明方法在大部分情况下应用的依然是传统的三角形算法。根据图4和图5的对比结果,随着星敏感器视场的减小,新方法的优越性渐渐体现出来,当视场大小为3°×3°时,新方法要明显由于传统算法。因为此时星敏感器捕获到的星点数目出现小于三颗的情况大大增加,传统算法对这种情况无能为力,只能依靠陀螺来对姿态进行解算,由于陀螺漂移的存在解算精度降低。同时从图中还可以看出,即使用本发明的方法,当星敏感器视场变小时姿态测量的精度依然要低于大视场情况下的测量精度,这时因为当只有一颗拍摄星时***姿态的状态量不完全可观测,此时只能利用得到的恒星位置信息来补偿一部分陀螺误差,却不能够完全补偿。而且,当视场变小时星图识别的失败率也会相应增加,此时的姿态测量完全依靠陀螺,从而导致精度的较低。但是从总体来看,本文的方法能够借助识别成功的一颗或者两颗恒星来对姿态进行更新,相对于传统方法很大程度上提高了姿态计算的精度。因此,本发明的方法更加适用于小视场的星敏感器,对于低成本、小型卫星的姿态确定有着较大的意义。
Claims (5)
1.一种小视场星敏感器的单星及双星定姿方法,其特征是:
步骤一:判断视场内恒星的数目,若视场内恒星数目不小于3,利用三角形算法进行识别,若识别成功,则计算并输出姿态,否则跳到步骤四;
步骤二:若视场内恒星数目为1,利用上一时刻的姿态信息和陀螺输出信息进行星点识别,若识别成功,则结合卡尔曼滤波实时估计当前姿态,否则跳到步骤四;所述利用上一时刻的姿态信息和陀螺输出信息进行星点识别具体包括:
步骤A:根据载体上一时刻的姿态信息和陀螺的输出信息计算载体当前时刻的姿态;
步骤B:利用星敏感器与载体之间的安装矩阵计算星敏感器视轴的方向矢量P;具体包括:使用下标b代表载体坐标系,上标i代表惯性坐标系,当前时刻由陀螺计算得到的载体系到惯性坐标系的转换矩阵为用下标s代表星敏感器坐标系,则星敏感器坐标系到惯性坐标系的转换矩阵为星敏感器的视轴方向即为星敏感器坐标系的z轴方向,方向矢量P即为的第三列元素;
步骤C:利用P从星表中选择落入视场中的星,并将其投影到拍摄星图中;
步骤D:利用识别成功星点的位置信息实时估计当前姿态;
步骤三:若视场内恒星数目为2,利用上一时刻的姿态信息、陀螺输出信息和星点的亮度和位置信息进行星图识别,若识别成功,用双矢量算法计算姿态,否则跳到步骤四;所述利用上一时刻的姿态信息、陀螺输出信息和星点的亮度和位置信息进行星图识别具体包括:
步骤E:根据载体上一时刻的姿态信息和陀螺的输出信息计算载体当前时刻的姿态;
步骤F:利用星敏感器与载体之间的安装矩阵计算星敏感器视轴的方向矢量P;
步骤G:利用P从星表中选择落入视场中的星,并将其投影到拍摄星图中;具体包括:若投影星只有一颗,则计算这颗投影星与两颗拍摄星之间的欧拉距离,欧拉距离最小且视场内恒星数目为1的就是其对应的拍摄星,此时利用步骤二的方法估计载体的姿态;若没有投影星,则星图识别失败;若有两颗或多颗投影星,根据单星识别的原理,计算每颗拍摄星与每颗投影星之间的欧拉距离,从中选取欧拉距离最小且满足基础条件的星作为拍摄星的对应星;
步骤H:星图识别成功以后利用双矢量算法计算载体姿态;
步骤四:星图识别失败后利用陀螺信息估计当前时刻姿态。
2.根据权利要求1所述的小视场星敏感器的单星及双星定姿方法,其特征是步骤C具体包括:生成的投影星图中只有两颗星,拍摄星和投影星;若陀螺和星敏感器都不存在噪声,则拍摄星和其对应的投影星重合,陀螺和星敏感器的噪声在星图中体现为星点的位置噪声,从而使得拍摄星和其投影星不重合;设陀螺按照载体三轴安装,其漂移为1°/h,星敏感器更新周期为1秒,则由陀螺漂移引入的噪声为星敏感器噪声为3",则总的噪声为而在星图预处理时,星点之间的最小间距大于5个像素,当视场大小为4°×4°,像素为512×512时,此间距大约对应0.04°>>0.0013°,即陀螺漂移和星敏感器噪声产生的星点位置噪声要小于星图中星点之间的间距,将星点最小间距值d作为验证的基础条件;若只有拍摄星和投影星时,计算两颗星的欧拉距离,若其小于d,则认为投影星为该拍摄星的对应星;由于噪声的原因,偶尔会出现星图中出现多颗投影星的情况,此时在投影星图中计算拍摄星和每颗投影星之间的欧拉距离,距离最小的星若满足上述基础条件,则它就是拍摄星的对应星;若星图中没有投影星则此次星图识别失败。
3.根据权利要求2所述的小视场星敏感器的单星及双星定姿方法,其特征是步骤D具体包括:根据四元数姿态动力学,建立如下的状态方程
其中,q=[q0 q1 q2 q3]为载体的姿态四元数,M(ω)具体形式如下
其中,ωx、ωy ωz为陀螺的测量值;假设拍摄星点在星敏感器坐标系中的方向矢量为ps,其对应的投影星在地心赤道坐标系的方向矢量为pi,建立如下的量测方程
其中,为星敏感器的安装矩阵,C(q)为姿态四元数表示的姿态矩阵,其具体表示如下
根据上述建立的状态方程和量测方程,利用扩展卡尔曼滤波实时估计载体当前时刻的姿态。
4.根据权利要求3所述的小视场星敏感器的单星及双星定姿方法,其特征是步骤F具体包括:使用下标b代表载体坐标系,上标i代表惯性坐标系,当前时刻由陀螺计算得到的载体系到惯性坐标系的转换矩阵为用下标s代表星敏感器坐标系,则星敏感器坐标系到惯性坐标系的转换矩阵为星敏感器的视轴方向即为星敏感器坐标系的z轴方向,方向矢量P即为的第三列元素。
5.根据权利要求4所述的小视场星敏感器的单星及双星定姿方法,其特征是步骤H具体包括:若步骤F中将两颗星都识别成功,则得到两颗星在惯性坐标系的方向矢量U、V,利用这两个方向矢量建立三个单位矢量Xn、Yn、Zn,并构成一个新的正交坐标基,即
令Fn=[Xn Yn Zn]T,利用同样的方法建立星敏感器坐标系下的坐标基Fs=[Xs Ys Zs]T,则有由于式中Fs、Fn为3×3方阵,因此存在逆矩阵,则载体的姿态矩阵为
利用公式就计算出载体的姿态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410246926.9A CN104034334B (zh) | 2014-06-05 | 2014-06-05 | 一种小视场星敏感器的单星及双星定姿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410246926.9A CN104034334B (zh) | 2014-06-05 | 2014-06-05 | 一种小视场星敏感器的单星及双星定姿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104034334A CN104034334A (zh) | 2014-09-10 |
CN104034334B true CN104034334B (zh) | 2016-09-14 |
Family
ID=51465199
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410246926.9A Expired - Fee Related CN104034334B (zh) | 2014-06-05 | 2014-06-05 | 一种小视场星敏感器的单星及双星定姿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104034334B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104567865B (zh) * | 2014-12-29 | 2017-05-10 | 北京控制工程研究所 | 一种星敏感器空间粒子干扰条件下的姿态捕获方法 |
CN105241451B (zh) * | 2015-10-30 | 2018-01-05 | 中国电子科技集团公司第二十研究所 | 一种小视场星***天文观测信息异常检测方法 |
CN107246883B (zh) * | 2017-08-07 | 2020-11-03 | 上海航天控制技术研究所 | 一种高精度星敏感器安装矩阵在轨实时校准方法 |
CN110487260B (zh) * | 2019-07-18 | 2021-03-05 | 南京航空航天大学 | 一种航天器姿态确定智能星敏感器及其姿态确定方法 |
CN111854764B (zh) * | 2020-07-20 | 2022-05-27 | 中国科学院微小卫星创新研究院 | 基于星间测量信息的航天器姿态确定方法及*** |
CN113670299B (zh) * | 2021-08-16 | 2022-10-25 | 北京航空航天大学 | 一种用于星敏感器的串行并行组合式视场选通成像方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101196398A (zh) * | 2007-05-25 | 2008-06-11 | 北京航空航天大学 | 一种基于Euler-q算法和DD2滤波的航天器姿态确定方法 |
CN101830290A (zh) * | 2010-02-12 | 2010-09-15 | 哈尔滨工业大学 | 小天体撞击探测自主导航与制导控制规划调度方法 |
CN101846510A (zh) * | 2010-05-28 | 2010-09-29 | 北京航空航天大学 | 一种基于星敏感器和陀螺的高精度卫星姿态确定方法 |
CN102155945A (zh) * | 2011-03-08 | 2011-08-17 | 哈尔滨工业大学 | 一种提高ccd星敏感器动态性能的方法 |
CN103076015A (zh) * | 2013-01-04 | 2013-05-01 | 北京航空航天大学 | 一种基于全面最优校正的sins/cns组合导航***及其导航方法 |
CN103616028A (zh) * | 2013-11-29 | 2014-03-05 | 哈尔滨工程大学 | 一种基于单星敏感器的星光折射卫星自主导航方法 |
CN103630109A (zh) * | 2013-12-08 | 2014-03-12 | 北京航空航天大学 | 一种基于星光折射确定地心矢量的方法 |
CN103674021A (zh) * | 2013-11-25 | 2014-03-26 | 哈尔滨工业大学 | 基于捷联惯导与星敏感器的组合导航***及方法 |
-
2014
- 2014-06-05 CN CN201410246926.9A patent/CN104034334B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101196398A (zh) * | 2007-05-25 | 2008-06-11 | 北京航空航天大学 | 一种基于Euler-q算法和DD2滤波的航天器姿态确定方法 |
CN101830290A (zh) * | 2010-02-12 | 2010-09-15 | 哈尔滨工业大学 | 小天体撞击探测自主导航与制导控制规划调度方法 |
CN101846510A (zh) * | 2010-05-28 | 2010-09-29 | 北京航空航天大学 | 一种基于星敏感器和陀螺的高精度卫星姿态确定方法 |
CN102155945A (zh) * | 2011-03-08 | 2011-08-17 | 哈尔滨工业大学 | 一种提高ccd星敏感器动态性能的方法 |
CN103076015A (zh) * | 2013-01-04 | 2013-05-01 | 北京航空航天大学 | 一种基于全面最优校正的sins/cns组合导航***及其导航方法 |
CN103674021A (zh) * | 2013-11-25 | 2014-03-26 | 哈尔滨工业大学 | 基于捷联惯导与星敏感器的组合导航***及方法 |
CN103616028A (zh) * | 2013-11-29 | 2014-03-05 | 哈尔滨工程大学 | 一种基于单星敏感器的星光折射卫星自主导航方法 |
CN103630109A (zh) * | 2013-12-08 | 2014-03-12 | 北京航空航天大学 | 一种基于星光折射确定地心矢量的方法 |
Non-Patent Citations (6)
Title |
---|
基于CCD星图的空间目标精确光学定位方法研究;谌钊;《中国优秀硕士学位论文全文数据库基础科学辑》;20111215(第S2期);正文第22-56页 * |
基于双星敏感器的船体姿态测量***设计;郭敬明等;《光电子技术》;20140331;第34卷(第1期);5-10,20 * |
基于星敏感器双矢量观测信息的卫星姿态确定算法研究;***等;《科学技术与工程》;20111130;第11卷(第33期);8252-8255 * |
星图匹配制导中的关键技术;刘朝山等;《宇航学报》;20060130(第01期);31-35 * |
星图识别算法及空间飞行器姿态确定研究;史话;《中国优秀博硕士论文全文数据库(硕士)工程科技II辑》;20070415(第4期);正文第30-48页 * |
组合星图技术中的双星敏感器安装角度分析;崔祥祥等;《光电工程》;20131231;第40卷(第12期);49-53 * |
Also Published As
Publication number | Publication date |
---|---|
CN104034334A (zh) | 2014-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104034334B (zh) | 一种小视场星敏感器的单星及双星定姿方法 | |
US20190323845A1 (en) | Method and System for Accurate Long Term Simultaneous Localization and Mapping with Absolute Orientation Sensing | |
US20150201180A1 (en) | 3-d motion estimation and online temporal calibration for camera-imu systems | |
Yu et al. | Stereo‐vision‐based relative pose estimation for the rendezvous and docking of noncooperative satellites | |
CN105509750A (zh) | 一种天文测速与地面无线电组合的火星捕获段导航方法 | |
Van Pham et al. | Vision‐based absolute navigation for descent and landing | |
Ivanov et al. | Relative pose and inertia determination of unknown satellite using monocular vision | |
Li et al. | Autonomous navigation and environment modeling for MAVs in 3-D enclosed industrial environments | |
CN104931071A (zh) | 一种基于选权迭代的星敏感器在轨几何定标方法及*** | |
CN109211231A (zh) | 一种基于牛顿迭代法的炮弹姿态估计方法 | |
CN116184430A (zh) | 一种激光雷达、可见光相机、惯性测量单元融合的位姿估计算法 | |
EP3166849B1 (en) | System for locating the barycenter of at least one object orbiting in space and related process of physical and mechanical characterization of the identified object | |
CN103940432B (zh) | 一种星敏感器的姿态确定方法 | |
Toloei et al. | A new composite algorithm for identifying the stars in the star tracker | |
CN107727102A (zh) | 天文测速与地面无线电组合的火星捕获段导航方法 | |
Arani et al. | A geometric star identification algorithm based on triple triangle pattern | |
Olivares-Mendez et al. | Autonomous landing of an unmanned aerial vehicle using image-based fuzzy control | |
CN105424047B (zh) | 基于路标信息的航天器姿态指向误差辨识方法 | |
Yan et al. | Horizontal velocity estimation via downward looking descent images for lunar landing | |
CN110514199B (zh) | 一种slam***的回环检测方法及装置 | |
Luckett | Comparison of three machine vision pose estimation systems based on corner, line, and ellipse extraction for satellite grasping | |
CN105445763A (zh) | 基于跟瞄信息的目标重构方法 | |
Song et al. | ETA-INIT: Enhancing the Translation Accuracy for Stereo Visual-Inertial SLAM Initialization | |
Qiuying et al. | Research on pedestrian location based on dual MIMU/magnetometer/ultrasonic module | |
Arai et al. | Fast vision-based localization for a mars airplane |
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: 20160914 |
|
CF01 | Termination of patent right due to non-payment of annual fee |