CN107702709B - 一种非合作目标运动与惯性参数的时频域混合辨识方法 - Google Patents
一种非合作目标运动与惯性参数的时频域混合辨识方法 Download PDFInfo
- Publication number
- CN107702709B CN107702709B CN201710772043.5A CN201710772043A CN107702709B CN 107702709 B CN107702709 B CN 107702709B CN 201710772043 A CN201710772043 A CN 201710772043A CN 107702709 B CN107702709 B CN 107702709B
- Authority
- CN
- China
- Prior art keywords
- target
- formula
- equation
- time
- quaternion
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M1/00—Testing static or dynamic balance of machines or structures
- G01M1/12—Static balancing; Determining position of centre of gravity
- G01M1/122—Determining position of centre of gravity
- G01M1/125—Determining position of centre of gravity of aircraft
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Feedback Control In General (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Navigation (AREA)
Abstract
本发明涉及一种非合作目标运动与惯性参数的时频域混合辨识方法,首先通过视觉相机获得目标的位姿信息。本文提出一种在视觉测量目标位姿的基础上,解耦运动过程中的平动与转动,时域与频域同步辨识目标的状态与动力学参数,频域上采用DFT辨识与平动相关参量,时域上线性化空间漂浮目标动力学方程,采用EFIR辨识与转动相关参量。
Description
技术领域
本发明属于航天器参数辨识领域,涉及一种非合作目标运动与惯性参数的时频域混合辨识方法。
背景技术
空间操作作为未来航天技术发展的方向,越来越受到航天大国的重视。由平台基座、操作机械臂组成的空间机械臂机器人可广泛应用于在轨加注、故障卫星维护、轨道垃圾清理、辅助变轨等在轨服务任务。空间机械臂机器人为空间机器人的一种,其具有较大的工作空间和灵巧的操作性能。已经成为当前空间操作的主要手段。
空间操作的前提是对目标的抓捕,由于被抓捕对象为非合作目标,其运动状态与惯性参数未知,这对于抓捕及抓捕后复合体的稳定控制造成了很大的影响,已有的研究成果及实验都是针对参数已知的合作目标,或者采用激励的方法进行辨识,需要与目标物理连接,这在目标完全未知的条件下是有风险的。而已有的基于视觉的参数辨识都是在时域进行处理,需要提前已知或假设量测噪声与过程噪声,这在空间复杂背景条件下又是极其苛刻的。并且由于待估计参数较多,需要更长的处理时间,这对辨识精度与燃料消耗都是极其不利的。因此有必要发明一种更准确更快速的非合作目标参数辨识方法,以获得目标的运动速度、角速度以及惯性参量等参数。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,本发明提出一种非合作目标运动与惯性参数的时频域混合辨识方法,主要针对具备视觉观测系统的空间智能机器人对非合作目标的目标位置、平动速度、转动角速度以及惯性参量等运动参量及惯性参量的辨识方法.
技术方案
一种非合作目标运动与惯性参数的时频域混合辨识方法,其特征在于步骤如下:
步骤1:由视觉相机获得的非合作目标相对位姿信息;
步骤2、根据视觉相机获得的目标位姿信息以及航天器运动方程,采用频域最优滤波器对非合作目标的质心位置及平动速度进行估计:
2-1)设计频域滤波器,对非合作目标位置序列进行DFT变换,相当于把这个序列通过一个FIR数字滤波器,位置序列rs(n)的离散付里叶变换式
式中,N为变换点数目,rs(n)为n时刻目标的位置,由视觉相机获得,如图1所示。这相当于使序列rs(n)中的每个样本有相位滞后最后在n=N-1时刻把经过相移的各个样本叠加起来,就得出了序列rs(n)中频率为kΩ的频谱分量的付里叶系数X(k);
2-2)由于目标存在平动与转动,将方程进行平动与转动分解,
rs=r+R(q)ρ=r0+vt+R(q)ρ
其中:ρ为卫星质心到目标点位置矢量,ω为目标旋转角速度,ρ,ω均为在相机所在坐标系下的表达,r为目标质心位置,r0为初始时刻目标质心位置;相机所在坐标系相对于惯性坐标系的旋转通过四元数q表示;
角速度与四元数的运动学关系表示为:
质心位置r0与平动线速度v均为目标位置序列及其差分的直流分量;
步骤3、辨识质心位置r0与平动线速度v:根据视觉相机获得的目标位姿信息以及航天器运动方程,在时域上采用扩展有限冲击相应滤波器对非合作目标的运动角速度以及惯性参量进行估计:
3-1)设计时域滤波器,
EFIR滤波器首先假定窗口长度N,通过计算量测偏差的均方值在线估计得出;当获得第n时刻量测值时,从时刻m=n-N+1到时刻n之间的N个量测值为有效量测值,目标状态的估计值表示为
式中Gl为广义噪声功率增益,通过下式迭代计算得出
所述求解H和Φ的步骤如下:
根据航天器动力学方程,设计需要估计的状态变量x、推导算法所需要的观测矩阵H以及状态转移矩阵Φ表达式:
建立通过欧拉方程描述的目标姿态动力学方程:
其中,I为转动惯量,τ为扰动力矩;
引入惯量参数表示:
式中,Ixx,Iyy,Izz为目标主惯量,满足以下条件:
Ixx+Iyy>Izz
Iyy+Izz>Ixx
Izz+Ixx>Iyy
进而求得惯量参数的约束条件
px>-1,py>-1,pz>-1
式中
Ic=diag(Ixx,Iyy,Izz),tr(Ic)表示Ic的迹,ετ用来描述太阳能帆板支架、重力梯度产生的阻尼效果,选择采用过程噪声来描述;
视觉测量系统获得带有噪声的目标点相对位姿并解算到惯性坐标系下,目标位置由矢量rs表示,目标所在目标坐标系{C}相对于惯性坐标系{A}的姿态由四元数η表示,测量向量Z表示为:
式中:为测量噪声,对于空间自旋目标,测量得到的目标姿态η由预先未知的两个姿态变换得到,设目标坐标系{C}相对于相机坐标系{B}的姿态由四元数μ表示;相机坐标系{B}相对于惯性坐标系{A}的姿态由四元数q表示,则有下式成立:
x=[qv T ωT pT ρT μv T]T
式中:qv,μv分别为四元数的矢量部分;则观测方程中转动部分表示为:
式中
描述微小旋转的四元数方程:
式中:
至此,
X(n+1)=Φ(n)X(n)+V(n)
其中,V(n)表示离散的过程噪声,Φ(n)表示离散状态转移矩阵;
有益效果
本发明提出的一种非合作目标运动与惯性参数的时频域混合辨识方法,包括以下步骤:1)由视觉相机获得的非合作目标相对位姿信息;2)根据视觉相机获得的目标位姿信息以及航天器运动方程,采用频域最优滤波器对非合作目标的质心位置及平动速度进行估计;3)根据视觉相机获得的目标位姿信息以及航天器运动方程,在时域上采用扩展有限冲击相应滤波器对非合作目标的运动角速度以及惯性参量进行估计。
与现有技术相比,本发明具有如下有益效果:
本发明的基本思想是将航天器运动分解为平动与转动,采用数据的时域与频域同步处理,针对每一时刻检测获得的数据,采用扩展有限冲击响应滤波器对参数进行估计,获得目标的转动参数。对于成块数据,采用离散傅立叶变换加以处理,从而获得目标的平动参数。时频域的同步处理最大限度的使用了采样数据所提供的信息,降低了滤波器的估计维度,提高了处理速度,缩短了估计时间。加窗滤波器的引入提高了估计算法的鲁棒性,增加了对噪声的抑制能力。本发明具体拥有以下优点:
1.估计速度快
本发明将航天器运动分解为平动与转动,采用数据的时域与频域同步处理,针对每一时刻检测获得的数据,采用扩展有限冲击响应滤波器对参数进行估计,获得目标的转动参数。对于成块数据,采用离散傅立叶变换加以处理,从而获得目标的平动参数。时频域的同步处理最大限度的使用了采样数据所提供的信息,降低了滤波器的估计维度,提高了处理速度,缩短了估计时间。
2.鲁棒性强
本发明针对对传统的卡尔曼滤波器中对测量噪声与过程噪声估计不准确所造成的滤波发散问题,采用滤波器加窗处理,增强了系统对噪声的抑制能力,提高了系统的鲁棒性。
3.实时性高
本发明采用时频域同时处理的方法,减少了时域滤波器的状态维度,减小了计算量,从而缩短了运算时间,提高了系统的实时性。
附图说明
图1为目标位置测量示意图
图2为质心位置与线速度频率域估计
图3为目标位置角速度及惯性参数时域估计
具体实施方式
现结合实施例、附图对本发明作进一步描述:
为实现上述目的,本发明所采用的技术方案包括以下步骤:
1)由视觉相机获得的非合作目标相对位姿信息
2)根据视觉相机获得的目标位姿信息以及航天器运动方程,采用频域最优滤波器对非合作目标的质心位置及平动速度进行估计。
3)根据视觉相机获得的目标位姿信息以及航天器运动方程,在时域上采用扩展有限冲击相应滤波器对非合作目标的运动角速度以及惯性参量进行估计。
所述步骤1)中,计算目标惯性系坐标的具体步骤如下:
1-1):设相机坐标系中目标位置为P(xc,yc,zc),目标在相机坐标系中成像位置(u,v),则有以下公式成立:
式中(u0,v0)为图像坐标系的中心,(ax,ay)为有效焦距长度。
1-2):待求惯性系下坐标P(xw,yw,zw)便可以表示为:
其中旋转矩阵R、平移矩阵T分别为:
T=[tx,ty,tz]T (3)
1-3):将公式(3),(1)代入(2),可得到
其中B=[xw,yw,zw,1]T,待求惯性系下的坐标P(xw,yw,zw)可通过对相机坐标系中成像位置(u,v)的实时检测及公式(4)解算获得,坐标点的集合便构成目标运动轨迹。
所述步骤2)中,根据视觉相机获得的目标位姿信息以及航天器运动方程,采用频域最优滤波器对非合作目标的质心位置及平动速度进行估计的具体方法如下:
2-1)设计频域滤波器,对目标位置序列进行DFT变换,相当于把这个序列通过一个FIR数字滤波器,位置序列rs(n)的离散付里叶变换式
式中,N为变换点数目,rs(n)为n时刻目标的位置,由视觉相机获得,如图1所示。这相当于使序列rs(n)中的每个样本有相位滞后最后在n=N-1时刻把经过相移的各个样本叠加起来,就得出了序列rs(n)中频率为kΩ的频谱分量的付里叶系数X(k)。
2-2)推导目标位置方程,由于目标存在平动与转动,将方程进行平动与转动分解。
rs=r+R(q)ρ=r0+vt+R(q)ρ (6)
如图1所示,ρ为卫星质心到目标点位置矢量,ω为目标旋转角速度,ρ,ω均为在坐标系{B}下的表达,r为目标质心位置,r0为初始时刻目标质心位置。坐标系{B}相对于{A}的旋转通过四元数q表示,角速度与四元数的运动学关系可表示为:
根据公式(6)、(7),质心位置r0与平动线速度v均为目标位置序列及其差分的直流分量,因此可通过DFT滤波器准确估计出,并且由于估计值为直流分量,因而对高频观测噪声以及信号的非等间隔采样具有很强抑制性,算法鲁棒性好。估计结果入图2所示。
所述步骤3)中,根据视觉相机获得的目标位姿信息以及航天器运动方程,在时域上采用扩展有限冲击相应滤波器对非合作目标的运动角速度以及惯性参量进行估计的具体方法如下:
3-1)设计时域滤波器,由于空间环境复杂的特点,传统卡尔曼滤波器所要求的提前预知量测噪声与过程噪声无法满足,在噪声估计不准确的情况下,卡尔曼滤波器会产生误差累计进而发散,故本发明采用加窗处理的扩展有限冲击响应(EFIR)滤波器,本发明将其首次应用于航天器惯性参数及角速度估计。不同于采用递归算法的卡尔曼滤波器,EFIR滤波器采用迭代算法,因为在算法中做了加窗处理,不会产生因持续的偏差累计而导致滤波器发散的现象,鲁棒性好。
EFIR滤波器首先假定窗口长度N,通过计算量测偏差的均方值在线估计得出。当获得第n时刻量测值时,从时刻m=n-N+1到时刻n之间的N个量测值为有效量测值,这使得滤波器不需要预先知道噪声特性参数。目标状态的估计值表示为
式中Gl为广义噪声功率增益,通过下式迭代计算得出
可见,Gl的计算仅需已知矩阵Hl与Φl,Hl为l时刻观测矩阵,Φl为l时刻状态转移矩阵,在滤波器增益计算中,不需要预测噪声统计特性,该特点使算法适合噪声特性未知的复杂条件下的空间任务。算法流程如下所示。
3-2)根据航天器动力学方程,设计需要估计的状态变量x、推导算法所需要的观测矩阵H以及状态转移矩阵Φ表达式
为了估计目标的运动状态与参数,建立通过欧拉方程描述的目标姿态动力学方程:
其中,I为转动惯量,τ为扰动力矩,考虑到绕坐标轴三个方向转动的耦合性,无法通过目标运动轨迹分别独立估计出三个方向转动惯量,为解决这个问题,引入惯量参数表示:
式中,Ixx,Iyy,Izz为目标主惯量,应满足以下条件:
Ixx+Iyy>Izz
Iyy+Izz>Ixx
Izz+Ixx>Iyy
进而可求得惯量参数的约束条件
px>-1,py>-1,pz>-1
用惯性参数改写动力学方程为:
式中
Ic=diag(Ixx,Iyy,Izz),tr(Ic)表示Ic的迹,ετ用来描述太阳能帆板支架、重力梯度等产生的阻尼效果,可以选择采用过程噪声来描述。
假定视觉测量系统可以获得带有噪声的目标点相对位姿并可解算到惯性坐标系下,如图1所示,目标位置由矢量rs表示,目标所在坐标系{C}相对于惯性坐标系{A}的姿态由四元数η表示,测量向量Z可表示为:
式中:为测量噪声,对于空间自旋目标,测量得到的目标姿态η由预先未知的两个姿态变换得到,设坐标系{C}相对于坐标系{B}的姿态由四元数μ表示。坐标系{B}相对于{A}的姿态由四元数q表示,则有下式成立:
x=[qv T ωT pT ρT μv T]T (19)
式中:qv,μv分别为四元数的矢量部分。则观测方程中转动部分可表示为:
式中
该方程为四元数q,μ的非线性方程,为满足线性滤波器的要求,需要对观测方程线性化。考虑到在采样时间较短的条件下,目标有很小的位姿变化,从而有||δqv||<<1,||δμv||<<1,δq0≈1,观测矩阵中位置部分的线性化可由公式(10)得一阶近似公式:
描述微小旋转的四元数方程:
将公式(22)代入公式(18)并忽略高阶小项可得:
式中:
根据公式(8)、(16)所描述的目标动力学方程,将连续时间状态空间模型描述为:
对于方程(16)的线性化可直接由全微分公式给出:
将方程(26)写成离散形式,有
X(n+1)=Φ(n)X(n)+V(n) (28)
其中,V(n)表示离散的过程噪声,Φ(n)表示离散状态转移矩阵。对于时常系统有
本方法用于空间在轨服务任务中,具有不需要提前已知量测噪声与过程噪声,时域与频域同步处理的特点。本发明具有鲁棒性强,辨识速度快,辨识精度高等特点。
Claims (2)
1.一种非合作目标运动与惯性参数的时频域混合辨识方法,其特征在于步骤如下:
步骤1:由视觉相机获得目标位姿信息;
步骤2、根据视觉相机获得的目标位姿信息以及航天器运动方程,采用频域最优滤波器对非合作目标的质心位置及平动速度进行估计:
2-1)设计频域滤波器,对非合作目标位置序列进行DFT变换,相当于把这个序列通过一个FIR数字滤波器,位置序列rs(n)的离散付里叶变换式
式中,N为变换点数目,rs(n)为n时刻目标的位置,由视觉相机获得,这相当于使序列rs(n)中的每个样本有相位滞后最后在n=N-1时刻把经过相移的各个样本叠加起来,就得出了序列rs(n)中频率为kΩ的频谱分量的付里叶系数X(k);
2-2)由于目标存在平动与转动,将方程进行平动与转动分解,
rs=r+R(q)ρ=r0+vt+R(q)ρ
其中:ρ为卫星质心到目标点位置矢量,ω为目标旋转角速度,ρ,ω均为在相机所在坐标系下的表达,r为目标质心位置,r0为初始时刻目标质心位置;相机所在坐标系相对于惯性坐标系的旋转通过四元数q表示;
角速度与四元数的运动学关系表示为:
式中,E为单位矢量,qv为四元数矢量部分,q0为四元数标量部分,四元数表示为q=[qv Tq0]T;
R(q)=(2q0 2-1)E3+2q0[qv×]+2qvqv T
质心位置r0与平动线速度v均为目标位置序列及其差分的直流分量;
步骤3、辨识质心位置r0与平动线速度v:根据视觉相机获得的目标位姿信息以及航天器运动方程,在时域上采用扩展有限冲击相应滤波器对非合作目标的运动角速度以及惯性参量进行估计:
3-1)设计时域滤波器,
EFIR滤波器首先假定窗口长度N,通过计算量测偏差的均方值在线估计得出;当获得第n时刻量测值时,从时刻m=n-N+1到时刻n之间的N个量测值为有效量测值,目标状态的估计值表示为
式中Gl为广义噪声功率增益,通过下式迭代计算得出
对与转动有关的状态矢量x=[qv T ωT pT ρT μv T]T的辨识流程如下所示:
2.根据权利要求1所述非合作目标运动与惯性参数的时频域混合辨识方法,其特征在于:求解H和Φ的过程如下:
根据航天器动力学方程,设计需要估计的状态变量x、推导算法所需要的观测矩阵H以及状态转移矩阵Φ表达式:
建立通过欧拉方程描述的目标姿态动力学方程:
其中,I为转动惯量,τ为扰动力矩;
引入惯量参数表示:
式中,Ixx,Iyy,Izz为目标主惯量,满足以下条件:
Ixx+Iyy>Izz
Iyy+Izz>Ixx
Izz+Ixx>Iyy
进而求得惯量参数的约束条件
px>-1,py>-1,pz>-1
式中
Ic=diag(Ixx,Iyy,Izz),tr(Ic)表示Ic的迹,ετ用来描述太阳能帆板支架、重力梯度产生的阻尼效果,选择采用过程噪声来描述;
视觉测量系统获得带有噪声的目标点相对位姿并解算到惯性坐标系下,目标位置由矢量rs表示,目标所在目标坐标系{C}相对于惯性坐标系{A}的姿态由四元数η表示,测量向量Z表示为:
式中:为测量噪声,对于空间自旋目标,测量得到的目标姿态η由预先未知的两个姿态变换得到,设目标坐标系{C}相对于相机坐标系{B}的姿态由四元数μ表示;相机坐标系{B}相对于惯性坐标系{A}的姿态由四元数q表示,则有下式成立:
x=[qv T ωT pT ρT μv T]T
式中:qv,μv分别为四元数的矢量部分;则观测方程中转动部分表示为:
式中
该方程为四元数q,μ的非线性方程,为满足线性滤波器的要求,需要对观测方程线性化;有||δqv||<<1,||δμv||<<1,δq0≈1,观测矩阵中位置部分的线性化由公式R(q)=(2q0 2-1)E3+2q0[qv×]+2qvqv T得一阶近似公式:
描述微小旋转的四元数方程:
式中:
至此,
X(n+1)=Φ(n)X(n)+V(n)
其中,V(n)表示离散的过程噪声,Φ(n)表示离散状态转移矩阵;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710772043.5A CN107702709B (zh) | 2017-08-31 | 2017-08-31 | 一种非合作目标运动与惯性参数的时频域混合辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710772043.5A CN107702709B (zh) | 2017-08-31 | 2017-08-31 | 一种非合作目标运动与惯性参数的时频域混合辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107702709A CN107702709A (zh) | 2018-02-16 |
CN107702709B true CN107702709B (zh) | 2020-09-25 |
Family
ID=61171435
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710772043.5A Active CN107702709B (zh) | 2017-08-31 | 2017-08-31 | 一种非合作目标运动与惯性参数的时频域混合辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107702709B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109062051B (zh) * | 2018-08-28 | 2021-07-23 | 苏州艾利特机器人有限公司 | 一种提高机器人动力学参数辨识精度的方法 |
CN110470297A (zh) * | 2019-03-11 | 2019-11-19 | 北京空间飞行器总体设计部 | 一种空间非合作目标的姿态运动与惯性参数估计方法 |
CN110081906B (zh) * | 2019-03-28 | 2022-11-22 | 西北工业大学 | 基于吸附过程的非合作目标惯性特征参数的两步辨识方法 |
CN110203424B (zh) * | 2019-05-05 | 2021-04-20 | 中国人民解放军63921部队 | 利用测速数据估计航天器自旋运动的方法和设备 |
CN110132272A (zh) * | 2019-06-20 | 2019-08-16 | 河北工业大学 | 一种用于空间碎片运动参数的测量方法及系统 |
CN110186465B (zh) * | 2019-07-03 | 2022-08-05 | 西北工业大学 | 一种基于单目视觉的空间非合作目标相对状态估计方法 |
CN110906922A (zh) * | 2019-11-08 | 2020-03-24 | 沈阳无距科技有限公司 | 无人机位姿信息的确定方法及装置、存储介质、终端 |
CN111307176B (zh) * | 2020-03-02 | 2023-06-16 | 北京航空航天大学青岛研究院 | 一种vr头戴显示设备中视觉惯性里程计的在线标定方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104330799A (zh) * | 2014-11-17 | 2015-02-04 | 西安电子科技大学 | 一种基于粒子群滤波优化的isar成像方法 |
CN104656666A (zh) * | 2015-03-11 | 2015-05-27 | 哈尔滨工业大学 | 针对空间非合作目标的相对轨道设计及高精度姿态指向控制方法 |
CN105486312A (zh) * | 2016-01-30 | 2016-04-13 | 武汉大学 | 一种星敏感器与高频角位移传感器组合定姿方法及系统 |
CN106548475A (zh) * | 2016-11-18 | 2017-03-29 | 西北工业大学 | 一种适用于空间非合作自旋目标运动轨迹的预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9476715B2 (en) * | 2015-02-26 | 2016-10-25 | Space Systems/Loral, Llc | Navigational route selection to mitigate probability mobile terminal loses communication capability |
-
2017
- 2017-08-31 CN CN201710772043.5A patent/CN107702709B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104330799A (zh) * | 2014-11-17 | 2015-02-04 | 西安电子科技大学 | 一种基于粒子群滤波优化的isar成像方法 |
CN104656666A (zh) * | 2015-03-11 | 2015-05-27 | 哈尔滨工业大学 | 针对空间非合作目标的相对轨道设计及高精度姿态指向控制方法 |
CN105486312A (zh) * | 2016-01-30 | 2016-04-13 | 武汉大学 | 一种星敏感器与高频角位移传感器组合定姿方法及系统 |
CN106548475A (zh) * | 2016-11-18 | 2017-03-29 | 西北工业大学 | 一种适用于空间非合作自旋目标运动轨迹的预测方法 |
Non-Patent Citations (2)
Title |
---|
利用频域连续性检测非合作突发通信信号;胡亚 等;《北京理工大学学报》;20121030;第32卷(第10期);第1071-1076页 * |
基于一的非合作水声脉冲信号检测方法;王晓燕 等;《信号处理》;20110831;第27卷(第8期);第1271-1278页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107702709A (zh) | 2018-02-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107702709B (zh) | 一种非合作目标运动与惯性参数的时频域混合辨识方法 | |
CN108731670B (zh) | 基于量测模型优化的惯性/视觉里程计组合导航定位方法 | |
CN112649016B (zh) | 一种基于点线初始化的视觉惯性里程计方法 | |
CN108376411B (zh) | 一种基于双目视觉的非合作目标相对状态解算方法 | |
CN108645416B (zh) | 基于视觉测量系统的非合作目标相对导航仿真验证方法 | |
Piepmeier et al. | Tracking a moving target with model independent visual servoing: A predictive estimation approach | |
Wu et al. | A self-adaptive unscented Kalman filtering for underwater gravity aided navigation | |
CN108917772A (zh) | 基于序列图像的非合作目标相对导航运动估计方法 | |
CN110941183A (zh) | 一种基于神经网络的工业机器人动力学辨识方法 | |
CN110186465B (zh) | 一种基于单目视觉的空间非合作目标相对状态估计方法 | |
CN107167145B (zh) | 一种自适应非接触式失效卫星的形态参数测算方法 | |
CN110146092B (zh) | 基于导航信息评价的双体小行星探测轨迹优化方法 | |
CN113129377B (zh) | 一种三维激光雷达快速鲁棒slam方法和装置 | |
CN108645400B (zh) | 用于空间非合作目标相对导航的惯性参数辨识方法及系统 | |
CN110861081A (zh) | 一种欠约束索并联机器人末端执行器的自主定位方法 | |
McCann et al. | Rigid body pose estimation on TSE (3) for spacecraft with unknown moments of inertia | |
CN112407344B (zh) | 空间非合作目标的位姿预测方法和装置 | |
Wu et al. | Correspondence matching and time delay estimation for hand-eye calibration | |
Feng et al. | Pose and motion estimation of unknown tumbling spacecraft using stereoscopic vision | |
Meng et al. | Estimate of all the inertial parameters of a free-floating object in orbit | |
Aghili et al. | Adaptive motion estimation of a tumbling satellite using laser-vision data with unknown noise characteristics | |
CN113252029B (zh) | 一种基于光学陀螺量测信息的天文导航姿态传递方法 | |
Wu et al. | Adaptive robust Kalman filter for vision-based pose estimation of industrial robots | |
CN114764830A (zh) | 一种基于四元数ekf和未标定手眼系统的物体位姿估算方法 | |
Han et al. | Trajectory prediction of space robot for capturing non-cooperative target |
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 |