CN102621565B - 一种机载分布式pos的传递对准方法 - Google Patents
一种机载分布式pos的传递对准方法 Download PDFInfo
- Publication number
- CN102621565B CN102621565B CN2012101133957A CN201210113395A CN102621565B CN 102621565 B CN102621565 B CN 102621565B CN 2012101133957 A CN2012101133957 A CN 2012101133957A CN 201210113395 A CN201210113395 A CN 201210113395A CN 102621565 B CN102621565 B CN 102621565B
- Authority
- CN
- China
- Prior art keywords
- wing
- beta
- deflection deformation
- measurement unit
- inertial measurement
- 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
Landscapes
- Navigation (AREA)
Abstract
一种机载分布式POS的传递对准方法,利用力学的方法对载机机翼的挠曲运动进行建模,并将由机翼挠曲运动产生的挠曲变形角和挠曲变形角速度增广为卡尔曼滤波的状态变量。在此基础上,采用“速度+姿态”的匹配方式设计卡尔曼滤波器,利用主惯性测量单元与GPS进行信息融合后得到的速度、姿态基准信息对子惯性测量单元进行传递对准,最终获得每个子惯性测量单元安装点的经校正后的速度和姿态信息。本发明具有自主性强、精度高的特点,可用于提高载机存在挠曲变形时分布式POS的传递对准精度。
Description
技术领域
本发明涉及一种机载分布式POS的传递对准方法,可用于提高载机存在挠曲变形时分布式POS的传递对准精度。
背景技术
航空遥感是以飞机为观测平台,利用各种成像载荷获取地球表面或表层的大范围、高精度图像的战略技术,对国家经济建设和国家安全具有重大意义。随着我国经济的发展,航空遥感不断向着高分辨率、高精度、多遥感载荷集成及阵列载荷方向发展。
为实现遥感载荷的高精度成像,飞机需做匀速直线运动,但飞机在实际飞行中受气流扰动、飞控***误差等因素影响,必然偏离理想的匀速直线运动,导致遥感成像分辨率和精度严重下降。因此,高精度航空遥感必须进行运动补偿。位置姿态测量***(Position and Orientation System,POS)是一种特殊的惯性/卫星组合测量***,可精确测量遥感载荷中心的位置、速度和姿态等运动参数,是高分辨率航空遥感***实现运动补偿的关键设备。
但是对于装备了多个观测载荷的高性能航空遥感***,由于多个观测载荷安装在飞机的不同位置,飞机弹性变形导致载荷间的空间相对关系发生变化。此时,采用传统的单一POS无法实现多点的高精度位置姿态测量。因此,必须建立起高精度分布式时空基准***(分布式POS)为高性能航空遥感***中所有载荷提供高精度的时间、空间信息。分布式POS由主导航***和子惯性测量单元组成。其中主导航***即为采用高精度惯性测量单元的POS,可精确测量载体的运动参数;子惯性测量单元采用低精度的惯性测量单元。主惯性测量单元(导航解算后称为主惯导)通常安装在机腹中央,子惯性测量单元(导航解算后称为子惯导)安装在机翼两侧的不同载荷附近,用于测量载荷中心的运动参数。
但是,子惯导精度不高,测量误差随时间累积,需要主惯导不断对其进行传递对准,以达到高精度测量。在传递对准过程中,主惯导与GPS信息进行融合后所提供的信息可视为高精度的参考信息。利用精度较高的主惯导来校准子惯导的传递对准技术就是设法估计出失准角并消除其影响,提高对准精度。因此传递对准的核心是确定主惯性测量单元与子惯性测量单元之间的姿态失准角。此时,在飞机飞行过程中,主惯性测量单元与子惯性测量单元之间的姿态失准角包括两部分:固定安装误差角和机翼弹性运动引起的挠曲变形角。其中,固定安装误差角是常值,并且不受外界干扰的影响,容易确定。而挠曲变形角是由载机机动、载机内部振源或阵风等使载机机翼产生挠曲变形而产生的动态变形角,不易确定,对传递对准的精度影响较大。现有的对挠曲变形角进行补偿的方法有两种:一种是用经验建模的方法,将机翼的挠曲变形考虑为二阶或三阶的马尔科夫过程,并将挠曲变形角增广为卡尔曼滤波的状态变量,通过卡尔曼滤波估计出该挠曲变形角并进行补偿。用此方法时,模型的有关参数全凭经验设定。并且,对于安装在机翼不同位置的子惯性测量单元,经验模型和模型参数也完全相同。而对于安装在机翼不同位置的子惯性测量单元来说,挠曲运动变形是有明显差别的。靠近机腹处的挠曲变形程度较低,而靠近机翼尖端处的挠曲变形程度较高。这种不加区分的对整个机翼建立统一的挠曲变形运动模型是不够精确的,必然影响了传递对准的精度;另一种补偿方法是采用加大卡尔曼滤波器过程噪声矩阵来减小机翼变形和振动的影响。此方法是依据真实模型的协方差分析结果来确定注入白噪声的强度,即通过加大过程噪声来补偿建模的挠曲变形。该方法除了可以补偿挠曲运动外,还可以增加滤波器的鲁棒性,但是以降低传递对准的精度为代价的。对于机载分布式POS来说,其用途是运动补偿,因此对传递对准精度要求较高,现有的对机翼挠曲变形进行补偿的两种方法都很难保证对准精度。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提出一种机载分布式POS的传递对准方法,该方法可提高载机机翼存在挠曲变形时分布式POS的传递对准精度。
本发明的技术解决方案为:一种机载分布式POS的传递对准方法,具体步骤如下:
机载分布式POS的传递对准是利用卡尔曼滤波技术估计出主、子惯性测量单元间的失准角,包括:固定安装误差角λ和挠曲变形角μ。利用主惯性测量单元与GPS进行信息融合后得到的基准信息,对主、子惯性测量单元的失准角进行校正,最终获得每个子惯性测量单元安装点的高精度速度、姿态信息,完成传递对准。
(2)采用卡尔曼滤波技术进行主惯性测量单元与GPS的信息融合,获得主惯性测量单元安装点的速度和姿态基准信息;
(3)采用“速度+姿态”匹配方式建立卡尔曼滤波模型,将机翼由于挠曲变形运动产生的挠曲变形角μ和挠曲变形角速度,增广为卡尔曼滤波的状态变量。通过卡尔曼滤波估计出主、子惯性测量单元间速度误差,固定安装误差角λ和挠曲变形角μ;
(4)利用步骤(2)得到的基准信息和步骤(3)估计出的主、子惯性测量单元间的姿态失准角,对每个子惯性测量单元的速度和姿态进行校正,最终获得每个子惯性测量单元安装点的高精度的速度和姿态信息。
上述所采用的对机翼挠曲变形进行建模的方法为一种基于力学的建模方法,该方法的具体步骤为:
(1)利用ANSYS软件对机翼结构进行模态分析,得到机翼的第一阶弯曲模态频率ωm和扭转模态频率ωn;
(2)建立机翼弯曲和扭转运动与时间相关项的数学方程为:
式中,m(t)是弯曲运动与时间相关项;是对m(t)时间t的一阶导数;是m(t)对时间t的二阶导数;n(t)是扭转运动与时间相关项;是n(t)对时间t的一阶导数;是n(t)对时间t的二阶导数;ξm、ξn分别是弯曲和扭转的模态阻尼系数;ωm、ωn分别是弯曲和扭转的模态频率;fm(t)、fn(t)是白噪声;
(3)由于用于对地观测成像的飞机一般为平直翼飞机,且处于亚音速飞行状态,建立机翼弯曲和扭转运动与位置相关项的数学方程为:
式中:p(x)是机翼弯曲运动与空间位置相关项,q(x)是机翼扭转运动与空间位置相关项,自变量x为机翼上的坐标点,L为机翼长度,C1,C2为比例常数,是机翼结构特性固有参数,ω是机翼自由弯曲振动时的振动圆频率,是与机翼刚度和线质量相关的固有参数,EI是机翼的抗弯刚度,m代表机翼的线质量;
运用麦夸特曲线拟合的方法拟合得到p(x),q(x),即确定上式中的系数C1,C2,β;
(4)结合步骤(2)得到的m(t),n(t)以及步骤(3)得到的p(x),q(x),可建立载机机翼弹性运动方程为:
式中,γc(x,t)是机翼沿z方向的弯曲位移,γs(x,t)是机翼绕ox轴的扭转角位移,γc(x,t)和γs(x,t)均是随时间变化的;
(5)建立机翼挠曲变形角μ(x,t)=[μx(x,t) μy(x,t) μz(x,t)]T与扭转和弯曲位移关系为:
式中,μx(x,t)是机翼沿x方向的挠曲变形角,μy(x,t)是机翼沿y方向的挠曲变形角,μz(x,t)是机翼沿z方向的挠曲变形角。μx(x,t),μy(x,t),μz(x,t)均是随时间变化的。μ(x,t)的具体表达式为:
机翼挠曲变形角速度 为:
本发明的原理是:如图4所示,通过对机翼挠曲变形运动进行动力学建模,将由机翼挠曲运动引起的挠曲变形角和挠曲变形角速度增广为卡尔曼滤波的状态变量。利用卡尔曼滤波估计出主、子惯性测量单元间的固定安装误差角λ和挠曲变形角μ。利用主惯性测量单元与GPS进行信息融合后得到的基准信息,对子惯性测量单元速度和姿态进行校正,最终获得每个子惯性测量单元安装点的高精度速度、姿态信息,完成传递对准。
如图1所示,在对机翼挠曲变形运动进行建模时,充分考虑了机翼挠曲运动受时间因素和空间因素两方面影响。时间因素是指在不同时刻,机翼的挠曲变形程度不同。因此在建模时建立了机翼挠曲运动(包括弯曲和扭转)与时间相关项的数学方程;空间因素是指机翼不同位置的挠曲变形程度不同。比如,靠近机腹处的挠曲变形程度较低,而靠近机翼尖端处的挠曲变形程度较高。因此在建模时建立了机翼挠曲运动(包括弯曲和扭转)与位置相关项的数学方程。最后,将两方面因素综合到一起,建立了机翼挠曲变形运动整体的数学模型。数学模型中包括了机翼挠曲运动引起的挠曲变形角和挠曲变形角速度的模型,基于该模型可建立相应的卡尔曼滤波方程。通过卡尔曼滤波可估计出挠曲变形角。在此基础上,利用基准信息完成姿态校正,从而完成传递对准。
本发明与现有技术相比的优点在于:充分考虑了机翼挠曲变形对传递对准精度的影响,并通过将挠曲变形角及变形角速度增广为卡尔曼滤波状态变量的方法对其进行补偿。在建立卡尔曼滤波模型时,利用力学的方法对载机机翼挠曲变形进行建模,是在充分考虑载机机翼的结构特性的基础上,从力学的角度精确的建立机翼挠曲变形运动的数学模型,模型中考虑了机翼不同位置处挠曲变形的不同情况。在此基础上,利用卡尔曼滤波可精确的估计出安装误差角和挠曲变形角,从而完成子惯性测量单元的姿态校正,实现高精度的传递对准。
附图说明
图1为本发明的流程图;
图2为本发明的机翼挠曲变形运动建模的流程图;
图3卡尔曼滤波算法基本编排图;
图4为本发明的“速度+姿态”匹配传递对准原理图;
具体实施方式
由于分布式POS的主惯性测量单元通常安装在机腹中央,子惯性测量单元安装在机翼两侧的不同载荷附近,在大气湍流的影响下,机翼较机腹更容易发生挠曲变形,因此本发明对机翼的挠曲变形进行建模。
如图1所示,本发明的具体方法实施如下:
1、利用力学的方法建立载机机翼挠曲变形的运动模型,得到挠曲变形角μ和挠曲变形角速度的运动方程;
(1)利用ANSYS软件对机翼结构进行模态分析,得到机翼的第一阶弯曲模态频率ωm和扭转模态频率ωn;
(2)机翼弯曲和扭转运动与时间相关项的数学方程为:
式中,m(t)是弯曲运动与时间相关项;是对m(t)时间t的一阶导数;是m(t)对时间t的二阶导数;n(t)是扭转运动与时间相关项;是n(t)对时间t的一阶导数;是n(t)对时间t的二阶导数;ξm、ξn分别是弯曲和扭转的模态阻尼系数;ωm、ωn分别是弯曲和扭转的模态频率;fm(t)、fn(t)是白噪声;
(3)由于用于对地观测成像的飞机一般为平直翼飞机,且处于亚音速飞行状态,建立机翼弯曲和扭转运动与位置相关项的数学方程为:
式中:p(x)是机翼弯曲运动与空间位置相关项,q(x)是机翼扭转运动与空间位置相关项,自变量x为机翼上的坐标点,L为机翼长度,C1,C2为比例常数,是机翼结构特性固有参数,ω是机翼自由弯曲振动时的振动圆频率,是与机翼刚度和线质量相关的固有参数,EI是机翼的抗弯刚度,m代表机翼的线质量。
通过ANSYS对机翼结构的分析,可以得到x-p(x),x-q(x)曲线图。运用麦夸特曲线拟合的方法拟合得到p(x),q(x)即确定(2)式中的系数C1,C2,β;
(4)结合步骤(2)得到的m(t)和n(t)以及步骤(3)得到的p(x)和q(x),可建立载机机翼弹性运动方程为:
式中,γc(x,t)是机翼沿z方向的弯曲位移,γs(x,t)是机翼绕ox轴的扭转角位移,γc(x,t)和γs(x,t)均是随时间变化的;
(5)建立机翼挠曲变形角μ(x,t)=[μx(x,t) μy(x,t) μz(x,t)]T与扭转和弯曲位移关系为:
将式(2),式(3)代入式(4),得到:
式中,μx(x,t)是机翼沿x方向的挠曲变形角,μy(x,y)是机翼沿y方向的挠曲变形角,μz(x,t)是机翼沿z方向的挠曲变形角。μx(x,t),μy(x,t),μz(x,t)均是随时间变化的。p′(x)是机翼弯曲运动与空间位置相关项p(x)的一阶导数。
2、采用卡尔曼滤波进行主惯性测量单元与GPS的信息融合,获得主惯性测量单元安装点的速度和姿态基准信息,具体公式如下:
(1)建立***状态方程:
式中,X为***状态矢量,W为***噪声矢量,F为***转移矩阵,G为噪声转换矩阵:
其中,φx、φy和φz为平台失准角,δvx、δvy和δvz为速度误差,δL、δλ和δh分别为纬度误差、经度误差和高度误差,εx、εy和εz为陀螺仪随机常值漂移误差,和为加速度计随机常值偏置误差,FINS为惯性***矩阵,为***姿态转换矩阵;
(2)建立***的量测方程:
Z=HX+η (8)
其中,Z为观测矢量,H为观测矩阵,η为量测噪声,I为单位矩阵:
Z=[δL δλ δh δVE δVN δVU]T
(3)卡尔曼滤波基本算法编排,该算法的流程图如图3所示。
状态一步预测方程:
状态估值计算方程;
滤波增量方程:
其中,Pk/k-1、Rk分别为k时刻***协方差阵的一步预测、k时刻***量测噪声矩阵;
一步预测均方误差方程:
其中,Pk-1、Qk-1、Гk-1分别为k-1时刻***协方差阵、k-1时刻***噪声矩阵、k-1时刻***噪声驱动矩阵;
估计均方误差方程:
其中,Pk为k时刻***状态协方差阵,I为单位阵;
3、采用“速度+姿态”的匹配方式建立卡尔曼滤波模型,在考虑机翼挠曲变形的情况下,结合步骤1的建模,将机翼因挠曲变形运动产生的挠曲变形角μ和挠曲变形角速度增广为卡尔曼滤波的状态变量,即在确定出机翼挠曲变形角和挠曲变形角速度的运动模型后,将x,y方向的挠曲变形角μx,μy和将x,y方向的挠曲变形角速度增广为卡尔曼滤波的状态变量(z方向的挠曲变形角μz和挠曲变形角速度均为零);结合步骤2得到的主惯性测量单元安装点的速度和姿态的基准信息,通过卡尔曼滤波估计出主、子惯性测量单元间的姿态失准角;其中卡尔曼滤波器的模型包括状态方程和量测方程,分别如式(14)和式(15)所示。
(1)建立***状态方程:
式中:X为***状态矢量,W为***噪声矢量,F为***转移矩阵,G为噪声转换矩阵:
其中,φE是子惯性测量单元东向平台失准角,φN是子惯性测量单元东向平台失准角,φU是子惯性测量单元天向平台失准角;δVE是主、子惯性测量单元间的东向速度误差,δVN是主、子惯性测量单元间的北向速度误差,δVU是主、子惯性测量单元间的天向速度误差;εx是子惯性测量单元载体坐标系x轴陀螺常值漂移,εy是子惯性测量单元载体坐标系y轴陀螺常值漂移,εz是子惯性测量单元载体坐标系z轴陀螺常值漂移;是子惯性测量单元载体坐标系x轴加计常值偏置,是子惯性测量单元载体坐标系y轴加计常值偏置,是子惯性测量单元载体坐标系z轴加计常值偏置;λx是子惯性测量单元载体坐标系x轴固定安装误差角,λy是子惯性测量单元载体坐标系y轴固定安装误差角,λz是子惯性测量单元载体坐标系z轴固定安装误差角;μx是子惯性测量单元载体坐标系x轴挠曲变形角,μy是子惯性测量单元载体坐标系y轴挠曲变形角;是子惯性测量单元载体坐标系x轴挠曲变形角速度,是子惯性测量单元载体坐标系y轴挠曲变形角速度。
wεx是子惯性测量单元载体坐标系x轴陀螺随机误差,wεy是子惯性测量单元载体坐标系y轴陀螺随机误差,wεz是子惯性测量单元载体坐标系z轴陀螺随机误差;是子惯性测量单元载体坐标系x轴加计随机误差,是子惯性测量单元载体坐标系y轴加计随机误差,是子惯性测量单元载体坐标系z轴加计随机误差;fm(t)是机翼弯曲运动白噪声,fn(t)是机翼扭转运动白噪声。
其中,是主惯性测量单元的导航坐标系相对惯性坐标系的旋转角速度在导航坐标系上的投影分量组成的反对称阵;[fn×]是子惯性测量单元加计测得的比例在导航坐标系上的投影分量组成的反对称阵;是子惯性测量单元载体坐标系到导航坐标系的转换矩阵;ξm、ξn分别是弯曲和扭转的模态阻尼系数;ωm、ωn分别是弯曲和扭转的模态频率。
其中,RM是地球子午圈主曲率半径,RN是地球卯酉圈主曲率半径,H是高度。
其中,VE,VN,VU分别是惯性测量单元与GPS信息进行融合后东、北、天方向的基准速度,L是基准纬度,ωie是地球自转角速度。
(2)建立***量测方程:
当采用“速度+姿态”匹配方法,即量测量选为主、子惯性测量单元的速度误差和姿态误差进行分布式POS的传递对准时,量测值为主、子惯性测量单元的姿态角误差和速度误差。***量测方程为:
Z=HX+η (15)
式中:Z为量测变量,H为量测矩阵,η是量测噪声,假定其为零均值的高斯白噪声;
Z=[δψ δθ δγ δVE δVN δVU]T
其中,δψ、δθ、δγ分别是是主、子惯性测量单元的航向角误差、俯仰角误差和横滚角误差,即三个方向的姿态误差。δVE、δVN、δVU分别是是主、子惯性测量单元的东向、北向、天向速度误差,即三个方向的速度误差。
其中,p′(x)是机翼弯曲运动与空间位置相关项p(x)的一阶导数。q(x)是机翼扭转运动与空间位置相关项。
(3)卡尔曼滤波基本算法编排,根据上述***状态方程和量测方程,可建立离散卡尔曼滤波方程如下:
状态一步预测方程:
状态估值计算方程:
滤波增量方程:
其中,Pk/k-1、Rk分别为k时刻***协方差阵的一步预测、k时刻***量测噪声矩阵;
一步预测均方误差方程:
其中,Pk-1、Qk-1、Гk-1分别为k-1时刻***协方差阵、k-1时刻***噪声矩阵、k-1时刻***噪声驱动矩阵;
估计均方误差方程:
其中,Pk为k时刻***状态协方差阵,I是单位阵;
4、利用步骤2得到的基准信息和步骤3估计出的主、子惯性测量单元间的姿态失准角,对每个子惯性测量单元计算出速度、姿态进行校正,最终获得每个子惯性测量单元安装点的经校正后的速度和姿态信息。
(1)速度校正:
速度校正是直接利用子惯性测量单元量测得的速度减去在步骤3中利用卡尔曼滤波估计出的主、子惯性测量单元间的速度误差,即可得到子惯性测量单元被校正后的速度信息。
Vs_new=Vs_old-δV (21)
其中,Vs_new是子惯性测量单元被更新后的速度,Vs_old是子惯性测量单元测得到的速度,δV是用卡尔曼滤波估计出的主、子惯性测量单元间的速度误差。
(2)姿态校正:
利用估计出的固定安装误差角λ和挠曲变形角μ,计算主、子惯性测量单元的相对姿态矩阵为:
式中:(λ+μ)×为子***相对主***的失准角组成的反对称阵。
经展开计算后得到:
又因为
因此,由(23)式和(24)式,可以确定出航向角俯仰角θs和横滚角γs的主值,即
θs主=arcsin(T32) (26)
θs=θs主 (28)
通过对子惯性测量单元的速度、姿态进行校正,能够得到子惯性测量单元安装点的高精度速度、姿态信息,完成传递对准。
本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (1)
1.一种机载分布式POS的传递对准方法,具体步骤为:
(1.1)利用有限元软件对机翼结构进行模态分析,得到机翼的第一阶弯曲模态频率ωm和扭转模态频率ωn;
(1.2)建立机翼弯曲和扭转运动与时间相关项的数学方程为:
式中,m(t)是弯曲运动与时间相关项;是对m(t)时间t的一阶导数;是m(t)对时间t的二阶导数;n(t)是扭转运动与时间相关项;是n(t)对时间t的一阶导数;是n(t)对时间t的二阶导数;ξm、ξn分别是弯曲和扭转的模态阻尼系数;ωm、ωn分别是弯曲和扭转的模态频率;fm(t)、fn(t)是白噪声;
(1.3)机翼弯曲和扭转运动与位置坐标相关项的数学方程为:
式中:p(x)是机翼弯曲运动与空间位置相关项,q(x)是机翼扭转运动与空间位置相关项,自变量x为机翼上的坐标点,L为机翼长度,C1,C2为比例常数,是机翼结构特性固有参数,ω是机翼自由弯曲振动时的振动圆频率,是与机翼刚度和线质量相关的固有参数,EI是机翼的抗弯刚度,m代表机翼的线质量;
运用麦夸特曲线拟合的方法拟合得到p(x),q(x)即确定(2)式中的系数C1,C2,β;
(1.4)结合步骤(1.2)得到的m(t),n(t)以及步骤(1.3)得到的p(x),q(x),建立载机机翼完整弹性运动方程为:
式中,γc(x,t)是机翼沿z方向的弯曲位移,γs(x,t)是机翼绕ox轴的扭转角位移,γc(x,t)和γs(x,t)均是随时间变化的;
(1.5)建立机翼挠曲变形角μ(x,t)=[μx(x,t) μy(x,t) μz(x,t)]T与扭转和弯曲位移关系为:
式中,μx(x,t)是机翼沿x方向的挠曲变形角,μy(x,t)是机翼沿y方向的挠曲变形角,μz(x,t)是机翼沿z方向的挠曲变形角,μx(x,t),μy(x,t),μz(x,t)均是随时间变化的,μ(x,t)的具体表达式为:
机翼挠曲变形角速度 为:
(2)采用卡尔曼滤波技术进行主惯性测量单元与GPS的信息融合,获得主惯性测量单元安装点的速度和姿态基准信息;
(4)利用步骤(2)得到的基准信息和步骤(3)估计出的主、子惯性测量单元间的姿态失准角,对每个子惯性测量单元计算出速度、姿态进行校正,最终获得每个子惯性测量单元安装点的经校正后的速度和姿态信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012101133957A CN102621565B (zh) | 2012-04-17 | 2012-04-17 | 一种机载分布式pos的传递对准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012101133957A CN102621565B (zh) | 2012-04-17 | 2012-04-17 | 一种机载分布式pos的传递对准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102621565A CN102621565A (zh) | 2012-08-01 |
CN102621565B true CN102621565B (zh) | 2013-12-04 |
Family
ID=46561591
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012101133957A Active CN102621565B (zh) | 2012-04-17 | 2012-04-17 | 一种机载分布式pos的传递对准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102621565B (zh) |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103196448B (zh) * | 2013-03-22 | 2016-06-08 | 南京理工大学 | 一种机载分布式惯性测姿***及其传递对准方法 |
CN104950688B (zh) * | 2014-03-27 | 2017-11-14 | 南京理工大学 | 分布式传递对准半实物仿真***及其方法 |
CN104949687A (zh) * | 2014-03-31 | 2015-09-30 | 北京自动化控制设备研究所 | 一种长时间导航***全参数精度评估方法 |
CN103995918A (zh) * | 2014-04-17 | 2014-08-20 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种机翼变形和振动对飞机传递对准影响的分析方法 |
CN103913181B (zh) * | 2014-04-24 | 2017-03-29 | 北京航空航天大学 | 一种基于参数辨识的机载分布式pos传递对准方法 |
CN103983280B (zh) * | 2014-05-27 | 2017-02-15 | 上海新跃仪表厂 | 一种用速度差信息的天基惯性基准传递方法 |
CN104567930A (zh) * | 2014-12-30 | 2015-04-29 | 南京理工大学 | 一种能够估计和补偿机翼挠曲变形的传递对准方法 |
CN104655152B (zh) * | 2015-02-11 | 2017-06-27 | 北京航空航天大学 | 一种基于联邦滤波的机载分布式pos实时传递对准方法 |
CN106289246B (zh) * | 2016-07-25 | 2018-06-12 | 北京航空航天大学 | 一种基于位置和姿态测量***的柔性杆臂测量方法 |
CN107796388B (zh) * | 2016-09-06 | 2021-03-16 | 北京自动化控制设备研究所 | 一种基于惯性技术的相对姿态测量方法 |
CN107883952A (zh) * | 2016-09-30 | 2018-04-06 | 北京自动化控制设备研究所 | 一种姿态匹配过程中抗挠曲误差的自适应滤波方法 |
CN106595649B (zh) * | 2016-11-22 | 2019-10-22 | 北京航天自动控制研究所 | 一种飞行中惯性初始基准偏差补偿方法 |
CN106814383B (zh) * | 2017-01-16 | 2018-01-12 | 立得空间信息技术股份有限公司 | 一种适用于多种模式下的高精度pos快速对准方法 |
CN107747944B (zh) * | 2017-09-18 | 2020-03-17 | 北京航空航天大学 | 基于融合权重矩阵的机载分布式pos传递对准方法和装置 |
CN107765244B (zh) * | 2017-09-18 | 2020-10-27 | 北京航空航天大学 | 基于机载双天线InSAR基线测量方法和装置 |
CN108398130B (zh) * | 2018-02-22 | 2020-08-14 | 北京航空航天大学 | 挠曲形变测量网络的分布式pos传递对准建模方法和装置 |
CN108387227B (zh) * | 2018-02-22 | 2020-03-24 | 北京航空航天大学 | 机载分布式pos的多节点信息融合方法及*** |
CN109141476B (zh) * | 2018-09-27 | 2019-11-08 | 东南大学 | 一种动态变形下传递对准过程中角速度解耦合方法 |
CN110371318B (zh) * | 2019-05-17 | 2020-12-11 | 东南大学 | 一种动态变形下基于双重滤波器的传递对准方法 |
CN110514222B (zh) * | 2019-08-16 | 2023-05-23 | 中国航空工业集团公司西安飞行自动控制研究所 | 一种基于协方差分析的弹类惯导导航误差分析方法及装置 |
CN112100735B (zh) * | 2020-08-03 | 2022-11-11 | 东南大学 | 一种基于机翼形变的机载imu高精度参考基准获取方法 |
CN111811543B (zh) * | 2020-08-31 | 2020-12-11 | 蓝箭航天空间科技股份有限公司 | 回收式航天运载器分布式导航***初始对准方法 |
CN112326162B (zh) * | 2020-09-17 | 2021-07-06 | 北京航空航天大学 | 一种机载分布式pos用机翼弹性变形测量方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU642638B2 (en) * | 1989-12-11 | 1993-10-28 | Caterpillar Inc. | Integrated vehicle positioning and navigation system, apparatus and method |
US6580389B2 (en) * | 2000-08-11 | 2003-06-17 | The Regents Of The University Of California | Attitude determination using a global positioning system |
CN100386647C (zh) * | 2005-01-25 | 2008-05-07 | 中国科学院电子学研究所 | 基于pos的机载干涉sar***及方法 |
CN102322873B (zh) * | 2011-08-23 | 2014-02-26 | 北京航空航天大学 | 一种分布式pos地面演示验证*** |
-
2012
- 2012-04-17 CN CN2012101133957A patent/CN102621565B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN102621565A (zh) | 2012-08-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102621565B (zh) | 一种机载分布式pos的传递对准方法 | |
CN103913181B (zh) | 一种基于参数辨识的机载分布式pos传递对准方法 | |
CN103245359B (zh) | 一种惯性导航***中惯性传感器固定误差实时标定方法 | |
CN108413887B (zh) | 光纤光栅辅助分布式pos的机翼形变测量方法、装置和平台 | |
CN102608596B (zh) | 一种用于机载惯性/多普勒雷达组合导航***的信息融合方法 | |
CN109974697A (zh) | 一种基于惯性***的高精度测绘方法 | |
CN100585602C (zh) | 惯性测量***误差模型验证试验方法 | |
CN106289246A (zh) | 一种基于位置和姿态测量***的柔性杆臂测量方法 | |
CN103674034B (zh) | 多波束测速测距修正的鲁棒导航方法 | |
CN102680004B (zh) | 一种挠性陀螺位置姿态测量***pos的标度因数误差标定与补偿方法 | |
CN101660914B (zh) | 耦合惯性位置误差的机载星光和惯性组合的自主导航方法 | |
CN101706284B (zh) | 提高船用光纤陀螺捷联惯导***定位精度的方法 | |
CN108398130B (zh) | 挠曲形变测量网络的分布式pos传递对准建模方法和装置 | |
CN103278163A (zh) | 一种基于非线性模型的sins/dvl组合导航方法 | |
CN103674030A (zh) | 基于天文姿态基准保持的垂线偏差动态测量装置和方法 | |
CN104019828A (zh) | 高动态环境下惯性导航***杆臂效应误差在线标定方法 | |
CN103076025B (zh) | 一种基于双解算程序的光纤陀螺常值误差标定方法 | |
CN103852760B (zh) | 一种基于刚性和柔性基线组合的多基线测量方法 | |
CN101949703A (zh) | 一种捷联惯性/卫星组合导航滤波方法 | |
CN104697526A (zh) | 用于农业机械的捷联惯导***以及控制方法 | |
CN103076026B (zh) | 一种捷联惯导***中确定多普勒计程仪测速误差的方法 | |
CN103256942A (zh) | 传递对准中考虑杆臂补偿的变形角测量方法 | |
CN107764261B (zh) | 一种分布式pos传递对准用模拟数据生成方法和*** | |
CN102519470A (zh) | 多级嵌入式组合导航***及导航方法 | |
CN104764467A (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 |