[发明专利]一种适用于筒装斜发射导弹的动基座传递对准方法有效
申请号: | 201911380922.9 | 申请日: | 2019-12-27 |
公开(公告)号: | CN111089588B | 公开(公告)日: | 2022-02-01 |
发明(设计)人: | 罗伟;曾奎;张力;杨帆;韩根森;张如飞;杨登峰 | 申请(专利权)人: | 湖北航天技术研究院总体设计所 |
主分类号: | G01C21/16 | 分类号: | G01C21/16;G01C21/20;G01C25/00 |
代理公司: | 武汉智权专利代理事务所(特殊普通合伙) 42225 | 代理人: | 唐勇 |
地址: | 430040 湖*** | 国省代码: | 湖北;42 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 适用于 筒装斜 发射 导弹 基座 传递 对准 方法 | ||
1.一种适用于筒装斜发射导弹的动基座传递对准方法,其特征在于,包括以下步骤:
将主惯导的初始位置、速度和姿态信息作为子惯导的初值,并进行惯导捷联解算,所述进行惯导捷联解算,具体包括以下步骤:
进行姿态更新计算:
qk=[q0(tk) q1(tk) q2(tk) q3(tk)]
其中:q0(tk),q1(tk),q2(tk),q3(tk)为递推更新前的姿态四元数,q0(tk+1),q1(tk+1),q2(tk+1),q3(tk+1)为递推更新后的姿态四元数,为三轴角增量,ΔT为子惯组采样周期,为三轴附加角速度,Fv→q(v)表示将选择矢量转换为四元数,Fq×v表示四元数与向量相乘;
进行速度更新计算:
Δv=Fq×v(qk,Δv)
其中:Δv为加速度计测量值,为递推下一步的北向、东向与天向速度,单位均为:米/秒,为当前位置的维度,单位:弧度,为子午圈半径和卯酉圈半径,单位:米,ωie为地球自转角速率,Δvsfm为速度旋转效应补偿项;
进行位置更新计算:
其中:Lk,λk,hk为当前步的地理纬度、经度与高度值,单位分别为:弧度、弧度与米,Lk+1,λk+1,hk+1为递推下一步的地理纬度、经度与高度值,单位分别为:弧度、弧度与米,为当前步的北向、东向与天向速度,单位均为:米/秒;
建立包括子惯导的三轴加表零偏和天向速度偏差在内的15维状态向量Xk,并根据导弹的起竖角度,调整测量噪声方差矩阵Rk,进行子惯导的对准解算,得到子惯导实时的姿态失准角和速度;
分别对子惯导实时的姿态和速度进行校正计算,并将校正结果作为下一周期导航计算的初值。
2.如权利要求1所述的一种适用于筒装斜发射导弹的动基座传递对准方法,其特征在于:所述根据导弹的起竖角度,调整测量噪声方差矩阵Rk,进行子惯导的对准解算,包括以下步骤:依次进行导航系速度增量计算、状态转移矩阵更新计算、滤波时间更新计算、观测向量更新计算、观测矩阵更新计算以及滤波测量更新计算。
3.如权利要求2所述的一种适用于筒装斜发射导弹的动基座传递对准方法,其特征在于:所述状态转移矩阵更新计算中的初始矩阵为15×15维的单位矩阵,其实时计算的相关项如下:
其中:为附加角速度,Ak(1:3,1:3)为矩阵Ak的1至3行与1至3列构成的3×3矩阵,Ak(4:6,1:3)为矩阵Ak的4至6行与1至3列构成的3×3矩阵,(·×)为列向量的反对称矩阵,Δvsfm为速度旋转效应补偿项。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于湖北航天技术研究院总体设计所,未经湖北航天技术研究院总体设计所许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201911380922.9/1.html,转载请声明来源钻瓜专利网。