CN111220182A - 一种火箭传递对准方法及*** - Google Patents
一种火箭传递对准方法及*** Download PDFInfo
- Publication number
- CN111220182A CN111220182A CN202010210484.8A CN202010210484A CN111220182A CN 111220182 A CN111220182 A CN 111220182A CN 202010210484 A CN202010210484 A CN 202010210484A CN 111220182 A CN111220182 A CN 111220182A
- Authority
- CN
- China
- Prior art keywords
- inertial navigation
- matrix
- sub
- attitude
- error
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
- G01C25/005—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Navigation (AREA)
Abstract
本申请提供一种火箭传递对准方法及***,该方法包括如下步骤:对子惯导进行惯性解算,更新子惯导的姿态矩阵和速度;根据更新后的子惯导的姿态矩阵计算预先构建的传递对准滤波模型中的状态一步转移矩阵和量测矩阵;根据计算的状态一步转移矩阵和量测矩阵,对预先构建的传递对准滤波模型中的各个状态变量进行卡尔曼滤波估计;对各个状态变量进行收敛性判断,若各个状态变量的均方误差小于预设阈值,则收敛完成,对准结束,否则,继续返回执行惯性解算步骤。本申请在主惯导的对准精度水平前提下,大幅降低设备成本,实现了动态环境下火箭的初始对准。
Description
技术领域
本申请涉及火箭技术领域,尤其涉及一种火箭传递对准方法及***。
背景技术
目前,传递对准是把基准坐标系传递给惯性导航***以进行初始对准的方法。火箭对准方法主要是光瞄,也有自主对准,其中,光瞄的方法是采用光学瞄准确定初始方位角,用捷联惯组加速度计输出计算水平姿态角,如《箭载捷联惯导***水平自对准的两种实用方法》(中国惯性技术学报,2009 17(4))。自主对准常采用以惯性系为参考坐标系的解析粗对准和卡尔曼滤波精对准的方法,如《一种可用于运载火箭的SINS/GNSS自主导航方案》(中国惯性技术学报,2017 25(1))。
然而,光瞄的方法存在如下缺陷:光瞄是采用远距离斜瞄或近距离平瞄等外部光学瞄准的方法,需要更专业的测试开发人员,配置与测试开发流程较为复杂,对火箭结构也有诸多要求,不具备全天候、无人化发射的能力,给火箭瞄准工作带来了许多限制,且光瞄通常用于静止或小幅晃动的情况,一旦振动幅度过大则会引起失效,无法满足诸如海上发射等环境下的对准要求。
自主对准的方法不需要借助外部设备进行对准,但其对惯导***的本身的精度提出了很高的要求,例如惯导器件中的陀螺组件精度通常在0.01度/小时以上,导致整体导航***的成本很高,传递对准常用于机载或弹上,需要机动条件进行误差分离,现有火箭发射中不存在类似的对准方案。
发明内容
本申请的目的在于提供一种火箭传递对准方法及***,该对准方法目的是在实现主惯导的对准精度水平前提下,大幅降低设备成本,实现在动态环境下火箭的初始对准。
为达到上述目的,本申请提供一种火箭传递对准方法,该方法包括如下步骤:
对子惯导进行惯性解算,更新子惯导的姿态矩阵和速度;
根据更新后的子惯导的姿态矩阵计算预先构建的传递对准滤波模型中的状态一步转移矩阵和量测矩阵;
根据计算的状态一步转移矩阵和量测矩阵,对预先构建的传递对准滤波模型中的各个状态变量进行卡尔曼滤波估计;
根据卡尔曼滤波估计的结果对各个状态变量进行收敛性判断,若各个状态变量的均方误差小于预设阈值,则收敛完成,对准结束,否则,继续更新子惯导的姿态矩阵和速度。
如上的,其中,构建传递对准滤波模型为:
其中,Xk是tk时刻15×1维的状态向量;Xk-1是tk-1时刻15×1维的状态向量;Φk/k- 115×15维的状态一步转移矩阵;Wk-1为tk-1时刻6×1维的***噪声向量;Zk是6×1维的量测向量;Hk表示tk时刻6×15维的量测矩阵;Vk表示tk时刻6×1维的量测噪声向量。
如上的,其中,将传递对准模型中的状态向量Xk定义为:
如上的,其中,将传递对准滤波模型中的观测向量Zk定义为:
Zk=[φ δv]T;
其中,φ表示子惯导与主惯导的姿态误差,δv表示子惯导与主惯导的速度误差;T表示矩阵的转置。
如上的,其中,子惯导与主惯导的姿态误差φ的计算方法如下:
φ=[atan2(C21,C11) -asin(C31) atan2(C32,C33)];
其中,atan()表示反正切函数;asin()表示反正弦函数。
如上的,其中,子惯导与主惯导的速度误差δv计算方法如下:
如上的,其中,计算状态一步转移矩阵的公式如下:
其中,Φk/k-1表示状态一步转移矩阵;符号×代表叉乘法;表示地球自转角速度;表示的反对称阵;fn表示导航系下加速度值;fn×表示fn的反对称阵;表示子惯导的姿态矩阵;T表示一个滤波周期;I3×3为3×3的单位阵;03×3表示3行3列0矩阵。
如上的,其中,计算量测矩阵Hk的公式如下:
如上的,其中,***噪声向量Wk-1表示为:
Wk-1=[wgN wgU wgE waN waU waE]T;
其中,wgi(i=N,U,E)为等效陀螺噪声;wai(i=N,U,E)为等效加速度计噪声;T表示矩阵的转置。
一种火箭传递对准***,包括:安装在发射架上的主惯导,安装在火箭上的子惯导;
更新模块,用于对子惯导进行惯性解算,更新子惯导的姿态矩阵和速度;
计算处理模块,用于根据更新后的子惯导的姿态矩阵计算预先构建的传递对准滤波模型中的状态一步转移矩阵和量测矩阵;
卡尔曼滤波估计模块,用于根据状态一步转移矩阵和量测矩阵,对预先构建的传递对准滤波模型中的各个状态变量进行卡尔曼滤波估计;
判断模块,对各个状态变量进行收敛性判断,若各个状态变量的均方误差小于预设阈值,则收敛完成,对准结束,否则,继续更新子惯导的姿态矩阵和速度。
本申请实现的有益效果如下:
(1)本发明用于舰船摇摆晃动条件下火箭初始姿态确定。箭上安装有子惯导用于敏感角速度和加速度信息,将另一高精度主惯导布置在火箭发射架上且靠近子惯导处,利用高精度的主惯导信息进行卡尔曼滤波估计,完成箭上子惯导的对准,提高了子惯导对准的精度和准确度。
(2)本发明对准方法实现主惯导的对准精度水平前提下,对器件精度的选择主要考虑组合导航的精度要求即可,极大降低了器件的成本,进而大幅降低设备成本,实现在动态环境下火箭的初始对准。
(3)本发明在大晃动条件下具有良好的效果,其考虑了杆臂误差、安装误差等误差,以及充分的机动条件下,子惯导能够达到接近主惯导的初始对准效果,本发明对子惯导的速度误差进行杆臂误差补偿,尽量使得主、子惯导运动保持一致。本发明与光瞄的对准相比,具有结构简单,仅需要在发射架(或者其他与子***由刚性连接的活动结构件)上安装主惯导,对准流程快捷,无需人工参与和可靠性高的优点。
(4)本发明以主惯导输出的信息作为信息源,进行惯性信息匹配,从而估计出子惯导的姿态误差,实现初始对准,并估计出部分器件误差数据。
(5)本发明结合要求的对准时间对滤波状态进行判断,对各个状态变量的收敛性判断进行状态固化,并在对准结束后完成误差校正。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域技术人员来讲,还可以根据这些附图获得其他的附图。
图1为本发明实施例的一种火箭传递对准方法的流程图。
图2为本发明实施例的更新子惯导的姿态矩阵和速度的流程图。
图3为本发明实施例的计算子惯导与主惯导姿态误差的流程图。
图4为本发明实施例的计算子惯导与主惯导速度误差的流程图。
图5为本发明实施例的一种火箭传递对准***的结构示意图。
图6为本发明实施例的主惯导和子惯导的安装示意图。
附图标记:1-发射架;2-火箭;10-主惯导;20-子惯导;30-初始化模块;40-更新模块;50-计算处理模块;60-卡尔曼滤波估计模块;70-判断模块;100-火箭传递对准***。
具体实施方式
下面结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
实施例一
如图1所示,本发明提供了一种火箭传递对准方法,该方法包括如下步骤:
步骤S1,利用主惯导信息对子惯导的姿态和速度进行初始化,获得子惯导的初始状态下的姿态矩阵和速度。
根据本发明的一个具体实施例,主惯导与子惯导的坐标系相同,则采用主惯导的姿态和速度的输出值分别作为子惯导的姿态和速度的初始值。
根据本发明的另一个具体实施例,主惯导与子惯导的坐标系不同,利用现有的导航中坐标系转换方法将主惯导的姿态和速度按照子惯导的坐标系定义进行转换怎么进行转换的?,将转换后的主惯导的姿态阵和速度作为初始值赋予子惯导,作为子惯导的姿态和速度的初始值。
步骤S2,对子惯导进行惯性解算,更新子惯导的姿态矩阵和速度。
如图2所示,步骤S2包括如下子步骤:
步骤S210包括:
步骤S212,根据子惯导的姿态变化矩阵更新子惯导的姿态矩阵,公式如下:
更新子惯导的速度的公式如下:
步骤S3,根据更新后的子惯导的姿态矩阵计算预先构建的传递对准滤波模型中的状态一步转移矩阵和量测矩阵。
根据本发明的一个具体实施例,构建传递对准滤波模型为:
其中,Xk是tk时刻15×1维的状态向量;Xk-1是tk-1时刻15×1维的状态向量;Φk/k- 115×15维的状态一步转移矩阵;Wk-1为tk-1时刻6×1维的***噪声向量;Zk是6×1维的量测向量;Hk表示tk时刻6×15维的量测矩阵;Vk表示tk时刻6×1维的量测噪声向量。
根据本发明的一个具体实施例,将传递对准滤波模型中的状态向量Xk定义为:
δv表示子惯导与主惯导的速度误差,速度误差包括北速度误差δvn、天速度误差δvu和东速度误差δve。
ε表示三轴陀螺常漂,三轴陀螺常漂包括X轴陀螺常漂εx、Y轴陀螺常漂εy和Z轴陀螺常漂εz。
λ表示子惯导相对主惯导安装误差,子惯导相对主惯导安装误差包括X轴安装误差λx、Y轴安装误差λy和Z轴安装误差λz。
其中,T表示矩阵的转置。
根据本发明的一个具体实施例,将传递对准滤波模型中的观测向量Zk定义为:
Zk=[φ δv]T; (6)
其中,φ表示子惯导与主惯导的姿态误差,δv表示子惯导与主惯导的速度误差;T表示矩阵的转置。
如图3所示,根据本发明的一个具体实施例,子惯导与主惯导的姿态误差φ的计算方法如下:
其中,表示姿态误差的转换矩阵;表示子惯导的姿态矩阵;表示主惯导的姿态矩阵;表示的转置;C11-C33为转换矩阵的各个元素;C11表示转换矩阵第1行第1列的元素;C12表示转换矩阵第1行第2列的元素;C11、C12、C13、C21、C22、C23、C31、C32和C33无实际意义。
具体的,子惯导与主惯导的姿态误差φ的计算公式为:
φ=[atan2(C21,C11) -asin(C31) atan2(C32,C33)]; (8)
其中,atan()表示反正切函数;asin()表示反正弦函数;C21表示转换矩阵第2行第1列的元素;C11表示转换矩阵第1行第1列的元素;C31表示转换矩阵第3行第1列的元素;C32表示转换矩阵第3行第2列的元素;C33表示转换矩阵第3行第3列的元素。
根据子惯导与主惯导的姿态误差可知:
如图4所示,根据本发明的一个具体实施例,子惯导与主惯导的速度误差δv计算方法如下:
步骤P1,计算杆臂速度误差。
如图6所示,杆臂δlb=lm-ls;lm为主惯导在载体系下(动坐标系,按X-Y-Z的轴向定义)坐标;ls为子惯导在载体系下(动坐标系,按X-Y-Z的轴向定义)坐标。
步骤P2,计算子惯导与主惯导的速度误差。
子惯导与主惯导的速度误差的计算公式如下:
步骤S3包括:
步骤S310,每个滤波周期内对***的状态一步转移矩阵进行计算。Φk/k-1为15×15维的状态一步转移矩阵,其计算方法如下:
其中,符号×代表叉乘法,表示地球自转角速度;表示的反对称阵;fn表示导航系下加速度值;fn×表示fn的反对称阵;表示子惯导的姿态矩阵;T表示一个滤波周期;I3×3为3×3的单位阵;;03×3表示3行3列0矩阵。
导航系下加速度值fn的计算公式如下:
步骤S320,每个滤波周期内对量测矩阵Hk进行计算。
根据本发明的具体实施例,6×1维的***噪声向量Wk-1如下:
Wk-1=[wgN wgU wgE waN waU waE]T; (14)
其中,wgi(i=N,U,E)为等效陀螺噪声;wai(i=N,U,E)为等效加速度计噪声;T表示矩阵的转置;N、U、E分别代表东北天坐标系中的北向、天向和东向;wgN表示北向等效陀螺噪声;wgU表示天向等效陀螺噪声;wgE表示东向等效陀螺噪声;waN表示北向加速度计噪声;waU表示天向等效加速度计噪声;waE表示东向等效加速度计噪声。
6×1维的***噪声向量Wk-1和6×1维的量测噪声向量Vk都是零均值的高斯白噪声向量序列(服从正态分布),且它们之间互不相关,即满足:
步骤S4,根据状态一步转移矩阵和量测矩阵,对预先构建的传递对准滤波模型中的各个状态变量进行卡尔曼滤波估计。
根据本发明的一个具体实施例,按照如下五个滤波方程进行卡尔曼滤波估计:
状态一步预测:
可知,这一步的状态向量由上一步的状态向量估计值和这一步的状态一步转移矩阵相乘获得的。
状态一步预测均方误差:
滤波增益方程为:
状态估计:
状态估计均方误差:
Pk=(I-KkHk)Pk/k-1; (19)
其中,I表示15×15维的单位阵;Hk表示滤波增益;Hk表示6×15维的量测矩阵;Pk表示tk时刻的状态变量的估计均方误差;Pk/k-1表示一个滤波周期内状态一步预测均方误差。
子惯导根据主惯导给出的姿态和速度进行传递对准,工作在闭环修正的导航状态,修正量来自卡尔曼滤波器。
步骤S5,对各个状态变量进行收敛性判断,若各个状态变量的均方误差小于预设阈值,则收敛完成,对准结束,否则,继续返回执行步骤S2。
采用范围判断和收敛趋势判断相结合的方法进行。
对状态变量的估计均方误差分别进行平稳性判断,当十五个状态变量的估计均方误差减小到设定阈值内时,则认为该状态变量收敛完成,并记录估计值。
步骤S6,对准结束后,误差进行一次修正,主惯导分离,火箭准备发射。
误差修正包括:姿态误差补偿,安装误差角步长和器件误差补偿。
误差修正完成后,主惯导固定在发射架1上,并与箭体分离。
实施例二
如图5和图6所示,一种火箭传递对准***100,包括安装在发射架1上的主惯导10,安装在火箭2上的子惯导20,主惯导10和子惯导20之间可视为刚性连接;
还包括初始化模块30,利用主惯导信息对子惯导的姿态矩阵和速度进行初始化;
更新模块40,用于对子惯导进行惯性解算,更新子惯导的姿态矩阵和速度;
计算处理模块50,用于根据更新后的子惯导的姿态矩阵计算预先构建的传递对准滤波模型中的状态一步转移矩阵和量测矩阵;
卡尔曼滤波估计模块60,用于根据状态一步转移矩阵和量测矩阵,对预先构建的传递对准滤波模型中的各个状态变量进行卡尔曼滤波估计;
判断模块70,对各个状态变量进行收敛性判断,若各个状态变量的均方误差小于预设阈值,则收敛完成,对准结束,否则,继续对子惯导进行对准。
本申请实现的有益效果如下:
(1)本发明用于舰船摇摆晃动条件下火箭初始姿态确定。箭上安装有子惯导用于敏感角速度和加速度信息,将另一高精度主惯导布置在火箭发射架上且靠近子惯导处,利用高精度的主惯导信息进行卡尔曼滤波估计,完成箭上子惯导的对准,提高了子惯导对准的精度和准确度。
(2)本发明对准方法实现主惯导的对准精度水平前提下,对器件精度的选择主要考虑组合导航的精度要求即可,极大降低了器件的成本,进而大幅降低设备成本,实现在动态环境下火箭的初始对准。
(3)本发明在大晃动条件下具有良好的效果,其考虑了杆臂误差、安装误差等误差,以及充分的机动条件下,子惯导能够达到接近主惯导的初始对准效果,本发明对子惯导的速度误差进行杆臂误差补偿,尽量使得主、子惯导运动保持一致。本发明与光瞄的对准相比,具有结构简单,仅需要在发射架(或者其他与子***由刚性连接的活动结构件)上安装主惯导,对准流程快捷,无需人工参与和可靠性高的优点。
(4)本发明以主惯导输出的信息作为信息源,进行惯性信息匹配,从而估计出子惯导的姿态误差,实现初始对准,并估计出部分器件误差数据。
(5)本发明结合要求的对准时间对滤波状态进行判断,对各个状态变量的收敛性判断进行状态固化,并在对准结束后完成误差校正。
Claims (10)
1.一种火箭传递对准方法,其特征在于,该方法包括如下步骤:
对子惯导进行惯性解算,更新子惯导的姿态矩阵和速度;
根据更新后的子惯导的姿态矩阵计算预先构建的传递对准滤波模型中的状态一步转移矩阵和量测矩阵;
根据计算的状态一步转移矩阵和量测矩阵,对预先构建的传递对准滤波模型中的各个状态变量进行卡尔曼滤波估计;
根据卡尔曼滤波估计的结果对各个状态变量进行收敛性判断,若各个状态变量的均方误差小于预设阈值,则收敛完成,对准结束,否则,继续更新子惯导的姿态矩阵和速度。
4.根据权利要求2所述的火箭传递对准方法,其特征在于,将传递对准滤波模型中的观测向量Zk定义为:
Zk=[φ δv]T;
其中,φ表示子惯导与主惯导的姿态误差,δv表示子惯导与主惯导的速度误差;T表示矩阵的转置。
9.根据权利要求2所述的火箭传递对准方法,其特征在于,***噪声向量Wk-1表示为:
Wk-1=[wgN wgU wgE waN waU waE]T;
其中,wgi(i=N,U,E)为等效陀螺噪声;wai(i=N,U,E)为等效加速度计噪声;T表示矩阵的转置。
10.一种火箭传递对准***,其特征在于,包括:安装在发射架上的主惯导,安装在火箭上的子惯导;
更新模块,用于对子惯导进行惯性解算,更新子惯导的姿态矩阵和速度;
计算处理模块,用于根据更新后的子惯导的姿态矩阵计算预先构建的传递对准滤波模型中的状态一步转移矩阵和量测矩阵;
卡尔曼滤波估计模块,用于根据状态一步转移矩阵和量测矩阵,对预先构建的传递对准滤波模型中的各个状态变量进行卡尔曼滤波估计;
判断模块,对各个状态变量进行收敛性判断,若各个状态变量的均方误差小于预设阈值,则收敛完成,对准结束,否则,继续更新子惯导的姿态矩阵和速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010210484.8A CN111220182B (zh) | 2020-03-23 | 2020-03-23 | 一种火箭传递对准方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010210484.8A CN111220182B (zh) | 2020-03-23 | 2020-03-23 | 一种火箭传递对准方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111220182A true CN111220182A (zh) | 2020-06-02 |
CN111220182B CN111220182B (zh) | 2022-02-11 |
Family
ID=70826471
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010210484.8A Active CN111220182B (zh) | 2020-03-23 | 2020-03-23 | 一种火箭传递对准方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111220182B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113175942A (zh) * | 2021-04-02 | 2021-07-27 | 上海机电工程研究所 | 动态环境下的杆臂速度补偿方法、***及介质 |
WO2022042241A1 (zh) * | 2020-08-31 | 2022-03-03 | 蓝箭航天空间科技股份有限公司 | 回收式运载火箭分布式导航***初始对准方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101963513A (zh) * | 2010-09-03 | 2011-02-02 | 哈尔滨工程大学 | 消除水下运载体捷联惯导***杆臂效应误差的对准方法 |
CN103047999A (zh) * | 2012-12-18 | 2013-04-17 | 东南大学 | 一种舰载主/子惯导传递对准过程中的陀螺误差快速估计方法 |
US20130231860A1 (en) * | 2012-03-01 | 2013-09-05 | Honeywell International Inc. | Systems and methods to incorporate master navigation system resets during transfer alignment |
CN105783943A (zh) * | 2016-04-26 | 2016-07-20 | 哈尔滨工程大学 | 一种基于无迹卡尔曼滤波的极区环境下舰船大方位失准角传递对准方法 |
CN109724626A (zh) * | 2019-03-08 | 2019-05-07 | 哈尔滨工程大学 | 一种针对极区传递对准动态挠曲杆臂效应的模型补偿方法 |
CN110388941A (zh) * | 2019-08-28 | 2019-10-29 | 北京机械设备研究所 | 一种基于自适应滤波的车辆姿态对准方法 |
CN110646012A (zh) * | 2018-06-27 | 2020-01-03 | 北京自动化控制设备研究所 | 一种惯导***单位置初始对准最优方法 |
-
2020
- 2020-03-23 CN CN202010210484.8A patent/CN111220182B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101963513A (zh) * | 2010-09-03 | 2011-02-02 | 哈尔滨工程大学 | 消除水下运载体捷联惯导***杆臂效应误差的对准方法 |
US20130231860A1 (en) * | 2012-03-01 | 2013-09-05 | Honeywell International Inc. | Systems and methods to incorporate master navigation system resets during transfer alignment |
CN103047999A (zh) * | 2012-12-18 | 2013-04-17 | 东南大学 | 一种舰载主/子惯导传递对准过程中的陀螺误差快速估计方法 |
CN105783943A (zh) * | 2016-04-26 | 2016-07-20 | 哈尔滨工程大学 | 一种基于无迹卡尔曼滤波的极区环境下舰船大方位失准角传递对准方法 |
CN110646012A (zh) * | 2018-06-27 | 2020-01-03 | 北京自动化控制设备研究所 | 一种惯导***单位置初始对准最优方法 |
CN109724626A (zh) * | 2019-03-08 | 2019-05-07 | 哈尔滨工程大学 | 一种针对极区传递对准动态挠曲杆臂效应的模型补偿方法 |
CN110388941A (zh) * | 2019-08-28 | 2019-10-29 | 北京机械设备研究所 | 一种基于自适应滤波的车辆姿态对准方法 |
Non-Patent Citations (2)
Title |
---|
YEON FUH JIANG 等: ""Error estimation of INS ground alignment through observability analysis"", 《 IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》 * |
王东升 等: ""舰载平台式惯导***的传递对准"", 《中国惯性技术学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022042241A1 (zh) * | 2020-08-31 | 2022-03-03 | 蓝箭航天空间科技股份有限公司 | 回收式运载火箭分布式导航***初始对准方法 |
CN113175942A (zh) * | 2021-04-02 | 2021-07-27 | 上海机电工程研究所 | 动态环境下的杆臂速度补偿方法、***及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN111220182B (zh) | 2022-02-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107525503B (zh) | 基于双天线gps和mimu组合的自适应级联卡尔曼滤波方法 | |
CN109596018B (zh) | 基于磁测滚转角速率信息的旋转弹飞行姿态高精度估计方法 | |
CN101413800B (zh) | 导航/稳瞄一体化***的导航、稳瞄方法 | |
CN100516775C (zh) | 一种捷联惯性导航***初始姿态确定方法 | |
CN110207697A (zh) | 基于角加速度计/陀螺/加速度计的惯性导航解算方法 | |
CN112611394B (zh) | 一种在发射坐标系下的飞行器姿态对准方法及*** | |
CN111121766B (zh) | 一种基于星光矢量的天文与惯性组合导航方法 | |
CN106979781B (zh) | 基于分布式惯性网络的高精度传递对准方法 | |
WO2004063669A2 (en) | Attitude change kalman filter measurement apparatus and method | |
CN112504275B (zh) | 一种基于级联卡尔曼滤波算法的水面舰船水平姿态测量方法 | |
CN111323050A (zh) | 一种捷联惯导和多普勒组合***标定方法 | |
CN110926468A (zh) | 基于传递对准的动中通天线多平台航姿确定方法 | |
CN108458709B (zh) | 基于视觉辅助测量的机载分布式pos数据融合方法和装置 | |
CN111220182B (zh) | 一种火箭传递对准方法及*** | |
CN113503894B (zh) | 基于陀螺基准坐标系的惯导***误差标定方法 | |
CN109708663B (zh) | 基于空天飞机sins辅助的星敏感器在线标定方法 | |
CN111207745A (zh) | 一种适用于大机动无人机垂直陀螺仪的惯性测量方法 | |
CN111189442A (zh) | 基于cepf的无人机多源导航信息状态预测方法 | |
CN113551668A (zh) | 一种航天器惯性/恒星星光矢量/星光折射组合导航方法 | |
CN112880669A (zh) | 一种航天器星光折射和单轴旋转调制惯性组合导航方法 | |
CN111141285B (zh) | 一种航空重力测量装置 | |
CN111207734B (zh) | 一种基于ekf的无人机组合导航方法 | |
Hong et al. | Application of EKF for missile attitude estimation based on “SINS/CNS” integrated guidance system | |
CN114061575B (zh) | 大失准角条件下的导弹姿态角精对准方法及*** | |
CN107702718B (zh) | 一种基于瞬间可观测度模型的机载pos机动优化方法与装置 |
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 |