CN108645400B - 用于空间非合作目标相对导航的惯性参数辨识方法及系统 - Google Patents
用于空间非合作目标相对导航的惯性参数辨识方法及系统 Download PDFInfo
- Publication number
- CN108645400B CN108645400B CN201810268150.9A CN201810268150A CN108645400B CN 108645400 B CN108645400 B CN 108645400B CN 201810268150 A CN201810268150 A CN 201810268150A CN 108645400 B CN108645400 B CN 108645400B
- Authority
- CN
- China
- Prior art keywords
- relative
- calculating
- maximum value
- inertial parameter
- minimum value
- 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/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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Image Analysis (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种用于空间非合作目标相对导航的惯性参数辨识方法及系统,其中,该方法包括如下步骤:步骤一:根据相机获得的不同时刻的相邻两张图像计算角速度ω,并计算动力学参数α;步骤二:根据步骤一中的角速度ω计算|ω|2的周期T、最大值β和最小值γ;步骤三:根据步骤二中得到的周期T、最大值β和最小值γ计算模数k;步骤四:根据步骤一中的动力学参数α、步骤二中的最大值β和最小值γ、步骤三中的模数k得到三个相对惯性参数Jx,Jy,Jz。本发明只利用视觉信息,在无外力和外力矩情况下估计非合作目标的惯性参数,进一步能够处理质量分布不均匀、惯量和几何参数完全未知的非合作目标。
Description
技术领域
本发明属于航天器自主导航与控制技术领域,尤其涉及一种用于空间非合作目标相对导航的惯性参数辨识方法及系统。
背景技术
空间非合作目标抓捕动力学与导航辨识是未来太空任务的重要方向,其中最关键技术之一就是识别非合作目标的运动状态和惯性参数。合作目标通常带有辅助标识和应答装置,可以通过指令进行动作配合。这类合作目标的状态感知、测量和抓捕操作相关技术已经非常成熟。但是对于失效卫星、敌方卫星、空间碎片和小行星等为代表的空间非合作目标,其共同特征就是无特定标识、无应答装置、无法执行动作配合,而且往往处于失效滚动状态。因此,确定这类非合作目标的运动状态和惯性参数是对其实施导航与抓捕操作的前提。
对于在轨目标进行识别的主要信息源就是视觉信息。传统的通过视觉信息估计空间运动目标位置和速度的方法都是基于惯性参数和几何参数已知的情况,也就是合作目标,因此不能用来处理惯性参数和几何参数未知的非合作目标,例如空间碎片,以及由于毁坏或变形导致参数改变的目标。目前,通过图像获得合作目标的运动状态和结构外形研究较多。针对非合作目标的惯性参数辨识研究很少,少量的研究也仅限于假设角动量保持常量的对称刚体模型。而传统的基于力的惯性参数辨识方法,大都以推进器产生的反作用推力为系统的输入,通过测量加速度信息以求得飞行器的惯性参数,需要已知外部激励,不适用于非合作目标的惯性参数辨识。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供了一种用于空间非合作目标相对导航的惯性参数辨识方法及系统,只利用视觉信息,在无外力和外力矩情况下估计非合作目标的惯性参数,进一步能够处理质量分布不均匀、惯量和几何参数完全未知的非合作目标。
本发明目的通过以下技术方案予以实现:根据本发明的一个方面,提供了一种用于空间非合作目标相对导航的惯性参数辨识方法,所述方法包括如下步骤:
步骤一:建立目标刚体的笛卡尔坐标系Oxyz,原点O过目标刚体质心,根据相机获得的不同时刻的相邻两张目标刚体的图像计算角速度ω,并计算动力学参数α;
步骤二:根据步骤一中的角速度ω计算|ω|2的周期T、最大值β和最小值γ;
步骤三:根据步骤二中得到的周期T、最大值β和最小值γ计算模数k;
步骤四:根据步骤一中的动力学参数α、步骤二中的最大值β和最小值γ、步骤三中的模数k得到三个相对惯性参数Jx,Jy,Jz。
上述用于空间非合作目标相对导航的惯性参数辨识方法中,在步骤一中,角速度ω为:
其中,v为单位矢量,θ为角度,Δt=t2-t1是从t1到t2时刻的时间间隔。
上述用于空间非合作目标相对导航的惯性参数辨识方法中,在步骤一中,动力学参数α为:
上述用于空间非合作目标相对导航的惯性参数辨识方法中,在步骤二中,周期T通过|ω|2的FFT变换得到。
上述用于空间非合作目标相对导航的惯性参数辨识方法中,在步骤三中,模数k通过下式求解:
上述用于空间非合作目标相对导航的惯性参数辨识方法中,在步骤四中,三个相对惯性参数Jx,Jy,Jz的计算步骤如下:
设f(ξ)=C2ξ2+C1ξ+C0=0的解ξ1<ξ2,其中C2=k2γ(αβ-1)-(β-γ),C1=-2{k2(αβ-1)-α(β-γ)},C0=α{k2(αβ-1)-α(β-γ)},C2、C1、C0均为方程系数;
当C2>0,α>Jy时,Jx=ξ1为负数,不符合实际情况,舍去;
根据本发明的另一方面,还提供了一种用于空间非合作目标相对导航的惯性参数辨识系统,包括:第一模块,用于根据相机获得的不同时刻的相邻两张图像计算角速度ω,并计算动力学参数α;第二模块,用于根据角速度ω计算|ω|2的周期T、最大值β和最小值γ;第三模块,用于根据周期T、最大值β和最小值γ计算模数k;第四模块,用于根据动力学参数α、最大值β、最小值γ、模数k得到三个相对惯性参数Jx,Jy,Jz。
其中,v为单位矢量,θ为角度,Δt=t2-t1是从t1到t2时刻的时间间隔;
上述用于空间非合作目标相对导航的惯性参数辨识系统中,周期T通过|ω|2的FFT变换得到;最大值β为:β=A;最小值γ为γ=A-B;其中, Jx为绕x轴的相对惯性参数,Jy为绕y轴的相对惯性参数,Jz为绕z轴的相对惯性参数。
上述用于空间非合作目标相对导航的惯性参数辨识系统中,模数k通过下式求解:
上述用于空间非合作目标相对导航的惯性参数辨识系统中,三个相对惯性参数Jx,Jy,Jz的计算步骤如下:
设f(ξ)=C2ξ2+C1ξ+C0=0的解ξ1<ξ2,其中C2=k2γ(αβ-1)-(β-γ),C1=-2{k2(αβ-1)-α(β-γ)},C0=α{k2(αβ-1)-α(β-γ)},C2、C1、C0均为方程系数;
当C2>0,α>Jy时,Jx=ξ1为负数,不符合实际情况,舍去;
本发明与现有技术相比具有如下有益效果:
(1)本发明提出了一种基于欧拉公式从一系列角速度描述的转动来计算动力学非合作目标相对惯性参数的方法,对非合作目标的在轨动力学参数辨识是一种有益的补充;
(2)本发明计算简单,概念明确,只利用视觉信息,可有效的在无外力和外力矩情况下估计非合作目标的相对惯性参数。该方法能够处理质量分布不均匀、惯量和几何参数完全未知的非合作目标。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1是本发明实施例提供的用于空间非合作目标近距离相对导航的惯性参数辨识方法的流程图;
图2是本发明实施例提供的刚体运动示意图;
图3是本发明实施例提供的角速度示意图。
具体实施方式
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例,然而应当理解,可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本公开,并且能够将本公开的范围完整的传达给本领域的技术人员。需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
图1是本发明实施例提供的用于空间非合作目标近距离相对导航的惯性参数辨识方法的流程图。如图1所示,用于空间非合作目标近距离相对导航的惯性参数辨识方法包括以下步骤:
步骤一:建立目标刚体的笛卡尔坐标系Oxyz,原点O过目标刚体质心,根据相机获得的不同时刻的相邻两张目标刚体的图像计算角速度ω,并计算动力学参数α;
步骤二:根据步骤一中的角速度ω计算|ω|2的周期T、最大值β和最小值γ;
步骤三:根据步骤二中得到的T、β和γ计算模数k;
步骤四:根据步骤一、步骤二和步骤三中得到的动力学参数α、β、γ和k计算相对惯性参数Jx,Jy,Jz。
在步骤一中,具体的,建立目标刚体的笛卡尔坐标系Oxyz,原点O过目标刚体质心,假设相机固定在惯性空间里,通过计算一系列角速度ω来估计未知刚体的运动。在t1和t2时刻的两幅图像中,设x1和x2分别是刚体上某一点在t1和t2时刻的空间位置矢量,则刚体从t1到t2时刻的运动可表示为
x2=Rx1+T (1)
其中,R∈R3×3是转动矩阵,T∈R3是平动向量。
转动矩阵R可表示为单位矢量v和角度θ的乘积。如果时间间隔Δt=t2-t1足够小,可以近似计算出角速度
在步骤一中,具体的,参见图2,角速度矢量总是位于固定面内,有下式成立
通过最小二乘法可得
在步骤二中,具体的,参见图3,周期T可通过|ω|2的FFT变换得到。
参见图2,刚体运动欧拉运动方程可表示为
其中,Ix为绕x轴的惯性参数,Iy为绕y轴的惯性参数,Iz为绕z轴的惯性参数,ωx为角速度ω的x分量,ωy为角速度ω的y分量,ωz为角速度ω的z分量。
求解式(6)和(7),可得
当α>Jy时
当α<Jy时
其中,sn(k,λt-δ),cn(k,λt-δ),dn(k,λt-δ)是雅克比椭圆函数。
由式(8)-(13)可得
|ω|2=A-Bsn2(k,λt-δ) (14)
通过最小二乘法,求得最大值β和γ为
β=A,γ=A-B (15)
在步骤三中,具体的,模数k可通过下式求解
f(ξ)=C2ξ2+C1ξ+C0=0 (18)
其中C2=k2γ(αβ-1)-(β-γ),C1=-2{k2(αβ-1)-α(β-γ)},C0=α{k2(αβ-1)-α(β-γ)}为方程系数。
设(18)式的解ξ1<ξ2。
当C2>0,α<Jy时,
当C2>0,α>Jy时,Jx=ξ1为负数,不符合实际情况,舍去。
当C2<0,α<Jy时,
当C2<0,α>Jy时,
到此就获得了非合作目标相对惯性参数Jx,Jy,Jz。
本实施例还提供了一种用于空间非合作目标相对导航的惯性参数辨识系统,包括:第一模块,用于根据相机获得的不同时刻的相邻两张图像计算角速度ω,并计算动力学参数α;第二模块,用于根据角速度ω计算|ω|2的周期T、最大值β和最小值γ;第三模块,用于根据周期T、最大值β和最小值γ计算模数k;第四模块,用于根据动力学参数α、最大值β、最小值γ、模数k得到三个相对惯性参数Jx,Jy,Jz。
上述实施例中,三个相对惯性参数Jx,Jy,Jz的计算步骤如下:设f(ξ)=C2ξ2+C1ξ+C0=0的解ξ1<ξ2,其中C2=k2γ(αβ-1)-(β-γ),C1=-2{k2(αβ-1)-α(β-γ)},C0=α{k2(αβ-1)-α(β-γ)},C2、C1、C0均为方程系数;当C2>0,α<Jy时,当C2>0,α>Jy时,Jx=ξ1为负数,不符合实际情况,舍去;当C2<0,α<Jy时,当C2<0,α>Jy时,
本实施例提出了一种基于欧拉公式从一系列角速度描述的转动来计算动力学非合作目标相对惯性参数的方法,对非合作目标的在轨动力学参数辨识是一种有益的补充;并且本发明计算简单,概念明确,只利用视觉信息,可有效的在无外力和外力矩情况下估计非合作目标的相对惯性参数。该方法能够处理质量分布不均匀、惯量和几何参数完全未知的非合作目标。
以上所述的实施例只是本发明较优选的具体实施方式,本领域的技术人员在本发明技术方案范围内进行的通常变化和替换都应包含在本发明的保护范围内。
Claims (4)
1.一种用于空间非合作目标相对导航的惯性参数辨识方法,其特征在于,所述方法包括如下步骤:
步骤一:建立目标刚体的笛卡尔坐标系Oxyz,原点O过目标刚体质心,根据相机获得的不同时刻的相邻两张目标刚体的图像计算角速度ω,并计算动力学参数α;
步骤二:根据步骤一中的角速度ω计算|ω|2的周期T、最大值β和最小值γ;
步骤三:根据步骤二中得到的周期T、最大值β和最小值γ计算模数k;
步骤四:根据步骤一中的动力学参数α、步骤二中的最大值β和最小值γ、步骤三中的模数k得到三个相对惯性参数Jx,Jy,Jz;
在步骤一中,动力学参数α为:
在步骤三中,模数k通过下式求解:
在步骤四中,三个相对惯性参数Jx,Jy,Jz的计算步骤如下:
设f(ξ)=C2ξ2+C1ξ+C0=0的解ξ1<ξ2,其中C2=k2γ(αβ-1)-(β-γ),C1=-2{k2(αβ-1)-α(β-γ)},C0=α{k2(αβ-1)-α(β-γ)},C2、C1、C0均为方程系数;
当C2>0,α>Jy时,Jx=ξ1为负数,不符合实际情况,舍去;
3.根据权利要求1所述的用于空间非合作目标相对导航的惯性参数辨识方法,其特征在于:在步骤二中,周期T通过|ω|2的FFT变换得到。
4.一种用于空间非合作目标相对导航的惯性参数辨识系统,其特征在于包括:
第一模块,用于根据相机获得的不同时刻的相邻两张图像计算角速度ω,并计算动力学参数α;
第二模块,用于根据角速度ω计算|ω|2的周期T、最大值β和最小值γ;
第三模块,用于根据周期T、最大值β和最小值γ计算模数k;
第四模块,用于根据动力学参数α、最大值β、最小值γ、模数k得到三个相对惯性参数Jx,Jy,Jz;
其中,v为单位矢量,θ为角度,Δt=t2-t1是从t1到t2时刻的时间间隔;
模数k通过下式求解:
三个相对惯性参数Jx,Jy,Jz的计算步骤如下:
设f(ξ)=C2ξ2+C1ξ+C0=0的解ξ1<ξ2,其中C2=k2γ(αβ-1)-(β-γ),C1=-2{k2(αβ-1)-α(β-γ)},C0=α{k2(αβ-1)-α(β-γ)},C2、C1、C0均为方程系数;
当C2>0,α>Jy时,Jx=ξ1为负数,不符合实际情况,舍去;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810268150.9A CN108645400B (zh) | 2018-03-29 | 2018-03-29 | 用于空间非合作目标相对导航的惯性参数辨识方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810268150.9A CN108645400B (zh) | 2018-03-29 | 2018-03-29 | 用于空间非合作目标相对导航的惯性参数辨识方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108645400A CN108645400A (zh) | 2018-10-12 |
CN108645400B true CN108645400B (zh) | 2020-08-14 |
Family
ID=63745181
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810268150.9A Active CN108645400B (zh) | 2018-03-29 | 2018-03-29 | 用于空间非合作目标相对导航的惯性参数辨识方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108645400B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110081906B (zh) * | 2019-03-28 | 2022-11-22 | 西北工业大学 | 基于吸附过程的非合作目标惯性特征参数的两步辨识方法 |
CN110361001B (zh) * | 2019-06-20 | 2023-04-07 | 河北工业大学 | 一种用于空间碎片运动测量系统及标定方法 |
CN112497240B (zh) * | 2020-11-24 | 2022-12-06 | 西北工业大学 | 非刚性抓捕的目标状态观测方法、计算机设备及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102636150A (zh) * | 2012-04-28 | 2012-08-15 | 北京航空航天大学 | 一种快速确定空间轴对称刚体目标姿态角的方法 |
CN103970964A (zh) * | 2014-05-23 | 2014-08-06 | 哈尔滨工业大学 | 一种挠性卫星模态参数在轨辨识方法 |
WO2014199073A1 (fr) * | 2013-06-12 | 2014-12-18 | Star Nav | Sextant électronique a système inertiel trois axes et procédé de détermination de la position |
-
2018
- 2018-03-29 CN CN201810268150.9A patent/CN108645400B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102636150A (zh) * | 2012-04-28 | 2012-08-15 | 北京航空航天大学 | 一种快速确定空间轴对称刚体目标姿态角的方法 |
WO2014199073A1 (fr) * | 2013-06-12 | 2014-12-18 | Star Nav | Sextant électronique a système inertiel trois axes et procédé de détermination de la position |
CN103970964A (zh) * | 2014-05-23 | 2014-08-06 | 哈尔滨工业大学 | 一种挠性卫星模态参数在轨辨识方法 |
Non-Patent Citations (2)
Title |
---|
"Estimation of state, shape, and inertial parameters of space objects from sequences of range images";Matthew D. Lichter;《Intelligent Robots and Computer Vision XXI: Algorithms, Techniques, and Active Vision》;20031231;正文第194-204页 * |
空间非合作目标惯性参数的Adaline网络辨识方法;孙俊等;《航空学报》;20161231(第09期);正文第2799-2806页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108645400A (zh) | 2018-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108645400B (zh) | 用于空间非合作目标相对导航的惯性参数辨识方法及系统 | |
Xu et al. | Autonomous rendezvous and robotic capturing of non-cooperative target in space | |
CN104457761B (zh) | 基于多目视觉的相对位置和姿态的特征接力方法 | |
CN108376411B (zh) | 一种基于双目视觉的非合作目标相对状态解算方法 | |
CN106595674B (zh) | 基于星敏感器和星间链路的heo卫星编队飞行自主导航方法 | |
WO2017113567A1 (zh) | 火星探测器自主导航方法 | |
Fourie et al. | Flight results of vision-based navigation for autonomous spacecraft inspection of unknown objects | |
CN110567461B (zh) | 一种考虑无陀螺仪的非合作航天器姿态和参数估计方法 | |
CN109612438B (zh) | 一种虚拟共面条件约束下的空间目标初轨确定方法 | |
CN110146224B (zh) | 一种辨识组合体航天器质量、质心位置和惯性张量的方法 | |
CN109178345B (zh) | 一种用于空中动目标跟踪的云台指向及星体姿态协同控制方法 | |
CN103217544B (zh) | 用星敏感器星点位置变化估算星体角速度的方法及系统 | |
Mumm et al. | Miniature control moment gyroscope development | |
Ivanov et al. | Relative pose and inertia determination of unknown satellite using monocular vision | |
Zhao et al. | Image-based control for rendezvous and synchronization with a tumbling space debris | |
CN110567462B (zh) | 一种近似自旋非合作航天器三轴转动惯量比的辨识方法 | |
Benninghoff et al. | Autonomous navigation for on-orbit servicing | |
Wang et al. | Ego-motion estimation of a quadrotor based on nonlinear observer | |
Hosen et al. | Vision-aided nonlinear observer for fixed-wing unmanned aerial vehicle navigation | |
CN110186465A (zh) | 一种基于单目视觉的空间非合作目标相对状态估计方法 | |
Du et al. | Pose estimation of large non-cooperative spacecraft based on extended PnP model | |
Meng et al. | Estimate of all the inertial parameters of a free-floating object in orbit | |
Sawada et al. | Shape estimation of IKAROS's solar power sail by images of monitor cameras | |
Tippetts et al. | Real-time vision sensor for an autonomous hovering micro unmanned aerial vehicle | |
Corke et al. | Sensing and control on the sphere |
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 |