CN101561281A - 一种捷联磁惯性组合系统的工作方法 - Google Patents

一种捷联磁惯性组合系统的工作方法 Download PDF

Info

Publication number
CN101561281A
CN101561281A CNA2009100847820A CN200910084782A CN101561281A CN 101561281 A CN101561281 A CN 101561281A CN A2009100847820 A CNA2009100847820 A CN A2009100847820A CN 200910084782 A CN200910084782 A CN 200910084782A CN 101561281 A CN101561281 A CN 101561281A
Authority
CN
China
Prior art keywords
module
signal
resolves
angle
omega
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
Application number
CNA2009100847820A
Other languages
English (en)
Other versions
CN101561281B (zh
Inventor
苏中
李兴城
施国兴
杨树兴
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
BEIJING XINGJIAN CHANGKONG OBSERVATION AND CONTROL TECHNOLOGY Co Ltd
Beijing Institute of Technology BIT
Original Assignee
BEIJING XINGJIAN CHANGKONG OBSERVATION AND CONTROL TECHNOLOGY Co Ltd
Beijing Institute of Technology BIT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by BEIJING XINGJIAN CHANGKONG OBSERVATION AND CONTROL TECHNOLOGY Co Ltd, Beijing Institute of Technology BIT filed Critical BEIJING XINGJIAN CHANGKONG OBSERVATION AND CONTROL TECHNOLOGY Co Ltd
Priority to CN2009100847820A priority Critical patent/CN101561281B/zh
Publication of CN101561281A publication Critical patent/CN101561281A/zh
Application granted granted Critical
Publication of CN101561281B publication Critical patent/CN101561281B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Gyroscopes (AREA)
  • Navigation (AREA)

Abstract

本发明涉及一种捷联磁惯性组合系统的工作方法。该方法中,信号调理模块对角速率信号和正弦信号进行滤波和放大;采集模块对信号调理模块调理后的角速率信号和正弦信号进行模/数变换,并将正弦信号发送至第一解算模块,将角速率信号输出至第二解算模块;第一解算模块根据正弦信号解算滚转角γ以及滚转角速率ωx;第二解算模块根据角速率信号和滚转角速率ωx解算方位角ψ和俯仰角θ;输出模块,用于输出滚转角γ、方位角ψ和俯仰角θ。本发明对于横滚角速率的测量不再使用原有的角速率陀螺,而是直接使用地磁滚转角测量单元来解算旋转弹的横滚角及横滚角速率,由于可以直接获取横滚角位置,从而避免了由积分获取位置量所带来的积分积累误差。

Description

一种捷联磁惯性组合系统的工作方法
技术领域
本发明属于导航技术测量领域,涉及一种捷联磁惯性组合系统的工作方法。
背景技术
捷联惯导系统从上世纪60年代初得到成熟的探索,尤其在1969年,捷联惯导系统作为“阿波罗”-13号登月飞船的应急备份装置,在其服务舱发生爆炸时将飞船成功的引导到返回地球的轨道上时起到了决定性的作用,成为捷联式惯导系统发展中的一个里程碑。从本世纪七十年代以来,在惯性传感器及惯性导航系统技术方面又取得了长足的进步,成功的研制了各种类型的惯性传感器(如激光陀螺、光纤陀螺、动态调谐陀螺等),并得到了广泛的应用。在系统技术方面,捷联式惯导系统也得到了极大的发展,除了精度要求比较高的战略应用(战略导弹、战略核潜艇)之外,捷联惯导系统几乎取代了所有的平台式惯导系统。据资料报道,美国军用惯导系统1984年前全部为平台惯导,到1989年一半改为捷联惯导,1994年捷联式已经占有90%。
目前用于旋转弹的捷联惯性组合,由于缺少敏感高速横滚角速率的陀螺仪,都在采取各种方法避免使捷联惯性测量单元承受弹体的旋转环境。如:美国为首的5国正在联合研制的GMLRS227mm远程制导火箭弹,为了避免惯性测量单元承受旋转环境,采取了增加一个滑动轴承,隔离发动机的旋转,通过空气舵使制导舱保持倾斜稳定。这种方式的缺点是系统复杂、增加了无效重量,牺牲了一部分舵资源。美国正在研究能够测量范围达3600°/s、标度因数误差小于5PPM的角速率陀螺仪,但是成本问题也难以在短期内得到解决。
发明内容
为了解决上述的技术问题,本发明提出了一种捷联惯性组合系统的工作方法,其目的在于解决对转速在10r/s~20r/s的旋转弹下的姿态测量,能实时解算旋转弹的姿态,具有很强的环境适用性,以有效解决现有旋转弹的制导改造需要。
本发明提供了一种捷联磁惯性组合系统,两只单轴角速率微机械陀螺、一个双轴地磁滚转角测量单元、第一解算模块、第二解算模块、信号调理模块、采集模块、输出模块和供电模块以捷联方式直接固连在弹体上;微机械陀螺与双轴地磁滚转角测量单元相互正交;
供电模块,用于为两只单轴角速率微机械陀螺、一个双轴地磁滚转角测量单元、第一解算模块、第二解算模块、信号调理模块、采集模块和输出模块供电;
信号调理模块,用于接收单轴角速率微机械陀螺输出的角速率信号并进行滤波和放大;用于接收双轴地磁滚转角测量单元输出的正弦信号并进行滤波和放大;
采集模块,用于对信号调理模块调理后的单轴角速率微机械陀螺输出的角速率信号和双轴地磁滚转角测量单元输出的正弦信号进行模/数变换,并将模/数变换后的正弦信号发送至第一解算模块,将模/数变换后的角速率信号输出至第二解算模块;
第一解算模块,用于根据模/数变换后的正弦信号解算滚转角γ以及滚转角速率ωx
第二解算模块,用于根据模/数变换后的角速率信号和第一解算模块发送的滚转角速率ωx解算方位角ψ和俯仰角θ;
输出模块,用于输出滚转角γ、方位角ψ和俯仰角θ。
第一解算模块根据以下公式解算滚转角:
当Hy>0时,γ=90-[arctan(Hz/Hy)]×180/π;
当Hy<0时,γ=270-[arctan(Hz/Hy)]×180/π;
当Hy=0且Z<0时,γ=180;
当Hy=0且Z>0时,γ=0;
其中,Hy和Hz分别为地磁敏感双轴信号。
第二解算模块根据以下方法解算方位角ψ和俯仰角θ:
根据公式 q 0 q 1 q 2 q 3 = 0 - ω x - ω y ω z ω x 0 ω z - ω y ω y - ω z 0 ω x ω z ω y - ω x 0 得到四元数:q0,q1,q2,q3
将游动方位坐标系与弹体坐标系的方向余弦矩阵式为等效为四元数姿态
T = C b n = C n b ′
矩阵: = q 0 2 + q 1 2 - q 2 2 - q 3 2 2 ( q 1 q 2 - q 0 q 3 ) 2 ( q 1 q 3 + q 0 q 2 ) 2 ( q 1 q 2 + q 0 q 3 ) q 0 2 - q 1 2 + q 2 2 - q 3 2 2 ( q 2 q 3 - q 0 q 1 ) 2 ( q 1 q 3 - q 0 q 2 ) 2 ( q 2 q 3 + q 0 q 1 ) q 0 2 - q 1 2 - q 2 2 + q 3 2 ;
= T 11 T 12 T 13 T 21 T 22 T 23 T 31 T 32 T 33
θ=arcsin(T32)
根据公式ψab=arctg(-T12/T22)得到方位角ψ和俯仰角θ;
ψ=ψab
其中,α为初始方位,ωy和ωz分别为两只单轴角速率微机械陀螺测得的角速率信号进行模/数变换后的角速率信号。
双轴地磁滚转角测量单元由双轴磁阻芯片实现。
信号调理模块由巴特沃斯滤波器和仪表放大器组成。
捷联磁惯性组合系统还包括相互正交的三只微机械加速度计;
信号调理模块,还用于接收微机械加速度计输出的加速度信号并进行滤波和放大;
采集模块,还用于对信号调理模块调理后的微机械加速度计输出的加速度信号进行模/数变换,并将模/数变换后的加速度信号输出至第二解算模块;
供电模块,还用于为微机械加速度计供电。
捷联磁惯性组合系统具有外壳、第一层框架板、第二层框架板和第三层框架板;
双轴地磁滚转角测量单元单独灌封成整体安装在第三层框架板上;
单轴角速率微机械陀螺和微机械加速度计灌封成整体,调整好交叉耦合后,安装在第二框架板上;供电模块安装在第二框架板上;
第一解算模块、第二解算模块、信号调理模块、采集模块和输出模块安装在第一层框架板上;
外壳上设置有安装座,通过安装孔与弹体固定连接。
外壳、第一层框架板、第二层框架板和第三层框架板均由无磁铝构成。
双轴地磁滚转角测量单元与第三层框架板之间垫有阻尼片;
单轴角速率微机械陀螺和微机械加速度计与第二框架板之间垫有阻尼片;供电模块与第二框架板之间垫有阻尼片;
第一解算模块、第二解算模块、信号调理模块、采集模块和输出模块与第一层框架板之间垫有阻尼片。
双轴地磁滚转角测量单元输出的正弦信号直接输入采集单元。
本发明提供了一种捷联磁惯性组合系统的工作方法,包括:
步骤S1,信号调理模块接收单轴角速率微机械陀螺输出的角速率信号并进行滤波和放大,以及接收双轴地磁滚转角测量单元输出的正弦信号并进行滤波和放大;
步骤S2,采集模块对信号调理模块调理后的单轴角速率微机械陀螺输出的角速率信号和双轴地磁滚转角测量单元输出的正弦信号进行模/数变换,并将双轴地磁滚转角测量单元输出的正弦信号发送至第一解算模块,模/数变换将单轴角速率微机械陀螺输出的角速率信号输出至第二解算模块;
步骤S3,第一解算模块根据模/数变换后的正弦信号解算滚转角γ以及滚转角速率ωx
步骤S4,第二解算模块根据模/数变换后的角速率信号和第一解算模块发送的滚转角速率ωx解算方位角ψ和俯仰角θ;
步骤S5,输出模块输出滚转角γ、方位角ψ和俯仰角θ;
捷联磁惯性组合系统中的两只单轴角速率微机械陀螺、一个双轴地磁滚转角测量单元、第一解算模块、第二解算模块、信号调理模块、采集模块、输出模块和供电模块以捷联方式直接固连在弹体上;微机械陀螺与双轴地磁滚转角测量单元相互正交;
供电模块为两只单轴角速率微机械陀螺、一个双轴地磁滚转角测量单元、第一解算模块、第二解算模块、信号调理模块、采集模块和输出模块供电。
步骤S3中,第一解算模块根据以下公式解算滚转角:
当Hy>0时,γ=90-[arctan(Hz/Hy)]×180/π;
当Hy<0时,γ=270-[arctan(Hz/Hy)]×180/π;
当Hy=0且Z<0时,γ=180;
当Hy=0且Z>0时,γ=0;
其中,Hy和Hz分别为地磁敏感双轴信号。
步骤S4中,第二解算模块根据以下方法解算方位角ψ和俯仰角θ:
根据公式 q 0 q 1 q 2 q 3 = 0 - ω x - ω y ω z ω x 0 ω z - ω y ω y - ω z 0 ω x ω z ω y - ω x 0 得到四元数:q0,q1,q2,q3
将游动方位坐标系与弹体坐标系的方向余弦矩阵式为等效为四元数姿态
T = C b n = C n b ′
矩阵: = q 0 2 + q 1 2 - q 2 2 - q 3 2 2 ( q 1 q 2 - q 0 q 3 ) 2 ( q 1 q 3 + q 0 q 2 ) 2 ( q 1 q 2 + q 0 q 3 ) q 0 2 - q 1 2 + q 2 2 - q 3 2 2 ( q 2 q 3 - q 0 q 1 ) 2 ( q 1 q 3 - q 0 q 2 ) 2 ( q 2 q 3 + q 0 q 1 ) q 0 2 - q 1 2 - q 2 2 + q 3 2 ;
= T 11 T 12 T 13 T 21 T 22 T 23 T 31 T 32 T 33
θ=arcsin(T32)
根据公式ψab=arctg(-T12/T22)得到方位角ψ和俯仰角θ;
ψ=ψab
其中,α为初始方位,ωy和ωz分别为两只单轴角速率微机械陀螺测得的角速率信号进行模/数变换后的角速率信号。
双轴地磁滚转角测量单元由双轴磁阻芯片实现。
信号调理模块由巴特沃斯滤波器和仪表放大器组成。
捷联磁惯性组合系统还包括相互正交的三只微机械加速度计;
步骤S1中,信号调理模块还接收微机械加速度计输出的加速度信号并进行滤波和放大;
步骤S2中,采集模块还对信号调理模块调理后的微机械加速度计输出的加速度信号进行模/数变换,并将模/数变换后的加速度信号输出至第二解算模块;
供电模块,还用于为微机械加速度计供电。
捷联磁惯性组合系统具有外壳、第一层框架板、第二层框架板和第三层框架板;
双轴地磁滚转角测量单元单独灌封成整体安装在第三层框架板上;
单轴角速率微机械陀螺和微机械加速度计灌封成整体,调整好交叉耦合后,安装在第二框架板上;供电模块安装在第二框架板上;
第一解算模块、第二解算模块、信号调理模块、采集模块和输出模块安装在第一层框架板上;
外壳上设置有安装座,通过安装孔与弹体固定连接。
外壳、第一层框架板、第二层框架板和第三层框架板均由无磁铝构成。
双轴地磁滚转角测量单元与第三层框架板之间垫有阻尼片;
单轴角速率微机械陀螺和微机械加速度计与第二框架板之间垫有阻尼片;供电模块与第二框架板之间垫有阻尼片;
第一解算模块、第二解算模块、信号调理模块、采集模块和输出模块与第一层框架板之间垫有阻尼片。
双轴地磁滚转角测量单元输出的正弦信号直接输入采集单元。
采用本发明的装置,可以实现如下有益效果:1)体积小,抗高过载冲击;2)对于横滚角速率的测量不再使用原有的角速率陀螺,而是直接使用地磁滚转角测量单元来解算旋转弹的横滚角及横滚角速率,由于可以直接获取横滚角位置,从而避免了由积分获取位置量所带来的积分积累误差。解算装置通过所测的横滚角及横滚角速率来完成对其它姿态信号的解旋测量。
附图说明
图1是磁惯性组合系统各模块结构图;
图2是磁惯性组合系统的内部结构示意图;
图3是磁惯性组合的产品图;
图4是利用地磁敏感双轴信号来解算滚转角的坐标示意图;
图5是坐标系的相互转换关系示意图;
图6是磁惯性组合系统的工作流程示意图。
具体实施方式
本发明的磁惯性组合系统是按照如下的技术方案实施的:主要由两只单轴角速率微机械陀螺、三只微机械加速度计、一个双轴地磁滚转角测量单元、弹载计算机及电源电路以捷联方式直接固连在结构本体上。所有元器件,MEMS传感器都采用固体芯片,分层通过阻尼片装在板级上,利用阻尼胶灌封在壳体内,成为整体。工作流程:装订旋转弹的射向、所处当地磁倾角、磁偏角,随着弹体的转动,地磁滚转角测量单元获取正弦变化的磁信号,然后根据这些装订参数,解算出主轴的滚转角和滚转速率。俯仰轴和偏航轴的速率可由安装在对应轴上的两只单轴角速率微机械陀螺测出,然后通过解旋计算所得。类似的三轴加速度可由正交安装的三只微机械加速度计所测的。
下面结合附图,对本发明做进一步的详细描述。
图1给出了本发明的系统模块结构图。它包括双轴地磁滚转角测量单元1,MEMS(微机电系统)陀螺2,MEMS加速度计3,信号调理模块4,采集模块5,单片机解算模块6,DSP姿态解算模块7和输出模块8。双轴地磁滚转角测量单元1由双轴磁阻传感器及其配套的信号处理电路组成。MEMS陀螺2由单轴微机械陀螺组成。MEMS加速度计3由三轴正交的微机械加速度计组成。信号调理模块4为差分放大电路和四阶巴特沃斯低通滤波器构成,滤波器截止频率100H。采集模块5由高速ADC板构成,它是12通道,采样精度16bit,单路采样速率50KHz的ADC板。单片计算机解算模块6中的单片计算机是采用Cygnal公司的C8051F020,最高运行频率可达到25MHz,同时具有比较大的FLASH和RAM。DSP姿态解算模块7采用TI公司的DSP芯片TMS320F2801及电源芯片TPS70151。输出模块8主要由串口芯片Max3488组成。电源模块9采用DC-DC模块和78系列稳压芯片为系统提供3.3V,5V,15V等电压。其中滚转角、滚转速率解算程序装载在单片计算机解算模块中,姿态解算算法装载在DSP姿态解算模块7中。
图1中的MEMS加速度计并不是本发明必须的模块,MEMS加速度计测得的信号经过调理和采集之后输入DSP姿态解算模块并处理后,由输出模块输出,它主要是用来为后续的导航提供依据。
工作时,由弹上电源通过供电模块9供电,利用双轴地磁滚转角测量模块1对地磁场进行敏感,利用MEMS陀螺2和MEMS加速度计3对转动速率和三轴及加速度进行敏感,随着弹体的转动,输出正弦的敏感信号,这些敏感信号经过信号调理模块4,经过采集模块5进行AD采样,分别输出给单片计算解算模块7和DSP解算模块8进行三个姿态角的解算,最后通过输出模块8输出给后续的控制系统,进行惯性制导。
图1中滚转角测量单元用于测出随弹体滚转变化的双轴正弦信号,它是采用双轴磁阻芯片HMC1022芯片实现的。
信号调理模块4主要用于信号滤波和放大,它可以采用二到四阶巴特沃斯滤波器进行滤波,能很好的滤除100Hz以上的干扰信号。同时利用仪表放大器对滤波后的信号进行放大,以供给采集模块5。
采集模块5是将模拟信号转换为数字信号,把地磁信号提供给单片机解算模块,陀螺信号和加速度信号提供给DSP姿态解算模块。
单片机解算模块利用地磁敏感双轴信号来解算滚转角,具体方法如下:
如图4所示:设地磁敏感双轴信号为Hy,Hz,当传感器位于不同滚转角指向时,Hy,Hz会随之发生变化,通过如下公式,即可解算出0°~360°的滚转角γ。根据各个时刻解算的滚转角γ,能实时给出滚转角速率。
由于需要解算的滚转角γ是在0°~360°内,而反正切函数解算的角度范围只有-90~90,故需要分段处理,如下式:
当Hy>0时:
γ=90-[arctan(z/y)]×180/π             (1)
当Hy<0时:
γ=270-[arctan(z/y)]×180/π            (2)
当Hy=0,Z<0时
γ=180                                  (3)
当Hy=0,Z>0时
γ=0                                    (4)
设滚转角速率为ωx,可由如下公式得到:
ωx=dγ/dt                              (5)
姿态解算模块采用游动方位坐标系的捷联式惯性导航数学平台,采用四元数法与等效转动矢量的方法进行姿态矩阵的解算以从中提取姿态角信息。
载体的姿态可用方位角ψ、俯仰角θ和滚转角γ表示。
如图4,游动方位坐标系与弹体坐标系的方向余弦矩阵式为:
C n b = cos γ cos ψ ab + sin γ sin θ sin ψ ab - cos γ sin ψ ab + sin γ sin θ sin ψ ab - sin γ cos θ cos θ sin ψ ab cos θ cos ψ ab sin θ sin γ cos ψ ab - cos γ sin θ sin ψ ab - sin γ sin ψ ab - cos γ sin θ cos ψ ab cos γ cos θ - - - ( 5 )
由四元法可知:
q 0 q 1 q 2 q 3 = 0 - ω x - ω y ω z ω x 0 ω z - ω y ω y - ω z 0 ω x ω z ω y - ω x 0 - - - ( 6 )
由(6)可得到四元数:q0,q1,q2,q3。Wy就是两只单轴MEMS陀螺测得的角速率。
则式(5)的方向余弦矩阵可等效为四元数姿态矩阵:
T = C b n = C n b ′
= q 0 2 + q 1 2 - q 2 2 - q 3 2 2 ( q 1 q 2 - q 0 q 3 ) 2 ( q 1 q 3 + q 0 q 2 ) 2 ( q 1 q 2 + q 0 q 3 ) q 0 2 - q 1 2 + q 2 2 - q 3 2 2 ( q 2 q 3 - q 0 q 1 ) 2 ( q 1 q 3 - q 0 q 2 ) 2 ( q 2 q 3 + q 0 q 1 ) q 0 2 - q 1 2 - q 2 2 + q 3 2 - - - ( 8 )
= T 11 T 12 T 13 T 21 T 22 T 23 T 31 T 32 T 33
对式(8)实时求解q0,q1,q2,q3的值,即可实时求出俯仰角和方位角为:
θ=arcsin(T32)
ψab=arctg(-T12/T22)                                 (9)
ψ=ψab
其中α为初始方位,其中滚转角通过地磁滚转角测量单元给出。
双轴地磁滚转角测量模块1输出的双轴正弦信号也可以不经过信号调理模块4,而直接输入到采集模块5;其余各模块之间的连接关系如图1所示,同样也可以实现本发明。
图2给出了本发明磁惯性组合的内部安装示意图;本发明磁惯性组合所有元器件都采用固体芯片,分装在3层框架板10~12,每块板级和固体芯片之间都垫有阻尼片。其中磁滚转角测量单元1单独灌封成整体,安装在框架板12上,MEMS陀螺2和MEMS加速度计3灌封成整体,调整好交叉耦合后,安装在框架板11上,框架板11的富余空间上安装供电模块9,计算机模块包括单片计算机解算模块6,DSP姿态解算模块7,安装在框架板10上。其中采集模块,信号调理模块和输出模块也安装在框架板10上。壳体13有安装座,通过安装孔14与弹体固定连接,壳体上有坐标线及产品编号(图中未画出),整个壳体7和框架板10~12都是采用无磁铝构成。
图3给出了本发明磁惯性组合产品外观图。
整个系统设计采用模块化处理,各个模块功能单一,方便调试。
图6是本发明提供的磁惯性组合系统的工作流程示意图,包括:
步骤601,信号调理模块接收单轴角速率微机械陀螺输出的角速率信号并进行滤波和放大,以及接收双轴地磁滚转角测量单元输出的正弦信号并进行滤波和放大;
步骤602,采集模块对信号调理模块调理后的单轴角速率微机械陀螺输出的角速率信号和双轴地磁滚转角测量单元输出的正弦信号进行模/数变换,并将双轴地磁滚转角测量单元输出的正弦信号发送至第一解算模块,模/数变换将单轴角速率微机械陀螺输出的角速率信号输出至第二解算模块;
步骤603,第一解算模块根据模/数变换后的正弦信号解算滚转角γ以及滚转角速率ωx
步骤604,第二解算模块根据模/数变换后的角速率信号和第一解算模块发送的滚转角速率ωx解算方位角ψ和俯仰角θ;
步骤605,输出模块输出滚转角γ、方位角ψ和俯仰角θ。
本领域的技术人员在不脱离权利要求书确定的本发明的精神和范围的条件下,还可以对以上内容进行各种各样的修改。因此本发明的范围并不仅限于以上的说明,而是由权利要求书的范围来确定的。

Claims (10)

1.一种捷联磁惯性组合系统的工作方法,其特征在于,包括:
步骤S1,信号调理模块接收单轴角速率微机械陀螺输出的角速率信号并进行滤波和放大,以及接收双轴地磁滚转角测量单元输出的正弦信号并进行滤波和放大;
步骤S2,采集模块对信号调理模块调理后的单轴角速率微机械陀螺输出的角速率信号和双轴地磁滚转角测量单元输出的正弦信号进行模/数变换,并将双轴地磁滚转角测量单元输出的正弦信号发送至第一解算模块,模/数变换将单轴角速率微机械陀螺输出的角速率信号输出至第二解算模块;
步骤S3,第一解算模块根据模/数变换后的正弦信号解算滚转角γ以及滚转角速率ωx
步骤S4,第二解算模块根据模/数变换后的角速率信号和第一解算模块发送的滚转角速率ωx解算方位角ψ和俯仰角θ;
步骤S5,输出模块输出滚转角γ、方位角ψ和俯仰角θ;
捷联磁惯性组合系统中的两只单轴角速率微机械陀螺、一个双轴地磁滚转角测量单元、第一解算模块、第二解算模块、信号调理模块、采集模块、输出模块和供电模块以捷联方式直接固连在弹体上;微机械陀螺与双轴地磁滚转角测量单元相互正交;
供电模块为两只单轴角速率微机械陀螺、一个双轴地磁滚转角测量单元、第一解算模块、第二解算模块、信号调理模块、采集模块和输出模块供电。
2.如权利要求1所述的捷联磁惯性组合系统的工作方法,其特征在于,步骤S3中,第一解算模块根据以下公式解算滚转角:
当Hy>0时,γ=90-[arctan(Hz/Hy)]×180/π;
当Hy<0时,γ=270-[arctan(Hz/Hy)]×180/π;
当Hy=0且Z<0时,γ=180;
当Hy=0且Z>0时,γ=0;
其中,Hy和Hz分别为地磁敏感双轴信号。
3.如权利要求2所述的捷联磁惯性组合系统的工作方法,其特征在于,步骤S4中,第二解算模块根据以下方法解算方位角ψ和俯仰角θ:
根据公式 q 0 q 1 q 2 q 3 = 0 - ω x - ω y ω z ω x 0 ω z - ω y ω y - ω z 0 ω x ω z ω y - ω x 0 得到四元数:q0,q1,q2,q3
将游动方位坐标系与弹体坐标系的方向余弦矩阵式为等效为四元数姿态矩阵:
T = C b n = C n b ′
= q 0 2 + q 1 2 - q 2 2 - q 3 2 2 ( q 1 q 2 - q 0 q 3 ) 2 ( q 1 q 3 + q 0 q 2 ) 2 ( q 1 q 2 + q 0 q 3 ) q 0 2 - q 1 2 + q 2 2 - q 3 2 2 ( q 2 q 3 - q 0 q 1 ) 2 ( q 1 q 3 - q 0 q 2 ) 2 ( q 2 q 3 + q 0 q 1 ) q 0 2 - q 1 2 - q 2 2 + q 3 2 ;
= T 11 T 12 T 13 T 21 T 22 T 23 T 31 T 32 T 33
θ=arcsin(T32)
根据公式ψab=arctg(-T12/T22)得到方位角ψ和俯仰角θ:
ψ=ψab
其中,α为初始方位,ωy和ωz分别为两只单轴角速率微机械陀螺测得的角速率信号进行模/数变换后的角速率信号。
4.如权利要求1所述的捷联磁惯性组合系统的工作方法,其特征在于,双轴地磁滚转角测量单元由双轴磁阻芯片实现。
5.如权利要求1所述的捷联磁惯性组合系统的工作方法,其特征在于,信号调理模块由巴特沃斯滤波器和仪表放大器组成。
6.如权利要求1所述的捷联磁惯性组合系统的工作方法,其特征在于,捷联磁惯性组合系统还包括相互正交的三只微机械加速度计;
步骤S1中,信号调理模块还接收微机械加速度计输出的加速度信号并进行滤波和放大;
步骤S2中,采集模块还对信号调理模块调理后的微机械加速度计输出的加速度信号进行模/数变换,并将模/数变换后的加速度信号输出至第二解算模块;
供电模块,还用于为微机械加速度计供电。
7.如权利要求6所述的捷联磁惯性组合系统的工作方法,其特征在于,捷联磁惯性组合系统具有外壳、第一层框架板、第二层框架板和第三层框架板;
双轴地磁滚转角测量单元单独灌封成整体安装在第三层框架板上;
单轴角速率微机械陀螺和微机械加速度计灌封成整体,调整好交叉耦合后,安装在第二框架板上;供电模块安装在第二框架板上;
第一解算模块、第二解算模块、信号调理模块、采集模块和输出模块安装在第一层框架板上;
外壳上设置有安装座,通过安装孔与弹体固定连接。
8.如权利要求7所述的捷联磁惯性组合系统的工作方法,其特征在于,外壳、第一层框架板、第二层框架板和第三层框架板均由无磁铝构成。
9.如权利要求7所述的捷联磁惯性组合系统的工作方法,其特征在于,
双轴地磁滚转角测量单元与第三层框架板之间垫有阻尼片;
单轴角速率微机械陀螺和微机械加速度计与第二框架板之间垫有阻尼片;供电模块与第二框架板之间垫有阻尼片;
第一解算模块、第二解算模块、信号调理模块、采集模块和输出模块与第一层框架板之间垫有阻尼片。
10.如权利要求7所述的捷联磁惯性组合系统的工作方法,其特征在于,双轴地磁滚转角测量单元输出的正弦信号直接输入采集单元。
CN2009100847820A 2009-05-19 2009-05-19 一种捷联磁惯性组合系统的工作方法 Expired - Fee Related CN101561281B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100847820A CN101561281B (zh) 2009-05-19 2009-05-19 一种捷联磁惯性组合系统的工作方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100847820A CN101561281B (zh) 2009-05-19 2009-05-19 一种捷联磁惯性组合系统的工作方法

Publications (2)

Publication Number Publication Date
CN101561281A true CN101561281A (zh) 2009-10-21
CN101561281B CN101561281B (zh) 2011-04-13

Family

ID=41220180

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100847820A Expired - Fee Related CN101561281B (zh) 2009-05-19 2009-05-19 一种捷联磁惯性组合系统的工作方法

Country Status (1)

Country Link
CN (1) CN101561281B (zh)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101900559A (zh) * 2009-11-06 2010-12-01 北京自动化控制设备研究所 一种捷联惯性导航系统双轴旋转调制方法
CN102322858A (zh) * 2011-08-22 2012-01-18 南京航空航天大学 用于地磁/捷联惯导组合导航系统的地磁匹配导航方法
CN103712622A (zh) * 2013-12-31 2014-04-09 清华大学 基于惯性测量单元旋转的陀螺漂移估计补偿方法及装置
CN103925917A (zh) * 2014-05-05 2014-07-16 上海新跃仪表厂 一种运载火箭姿态角速率信号的测量系统及方法
CN104764451A (zh) * 2015-04-23 2015-07-08 北京理工大学 一种基于惯性和地磁传感器的目标姿态跟踪方法
CN105737793A (zh) * 2016-05-11 2016-07-06 西安中星测控有限公司 滚转角测量单元及测量方法
CN106403935A (zh) * 2016-08-24 2017-02-15 易文俊 一种弹载磁组合姿态测量装置的应用
CN107063254A (zh) * 2016-12-27 2017-08-18 南京理工大学 一种陀螺地磁组合的姿态解算方法
CN107883944A (zh) * 2016-09-29 2018-04-06 北京航空航天大学 一种捷联星敏感器实现间接敏感地平的导弹姿态机动方法
CN108196570A (zh) * 2017-12-26 2018-06-22 深圳市道通智能航空技术有限公司 一种无人机航向修正方法、装置和无人机
CN110017808A (zh) * 2019-03-26 2019-07-16 北京理工大学 利用地磁信息和加速计解算飞行器姿态的方法
CN110765593A (zh) * 2019-10-09 2020-02-07 上海机电工程研究所 适用于便携式导弹初始装订信息的评估方法及系统
CN111366188A (zh) * 2018-12-26 2020-07-03 北京信息科技大学 一种面向动态环境力测量领域的数据采集存储装置
CN111443735A (zh) * 2020-03-25 2020-07-24 浙江大华技术股份有限公司 一种车载云台摄像机姿态保持的方法、装置及系统
CN111504256A (zh) * 2020-04-29 2020-08-07 中国北方工业有限公司 一种基于最小二乘法的滚转角实时估计方法
CN113587746A (zh) * 2021-10-08 2021-11-02 北京信息科技大学 基于地磁信息测量弹体的大跨度转速的方法、装置及系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7219033B2 (en) * 2005-02-15 2007-05-15 Magneto Inertial Sensing Technology, Inc. Single/multiple axes six degrees of freedom (6 DOF) inertial motion capture system with initial orientation determination capability
CN100489459C (zh) * 2006-07-17 2009-05-20 北京航空航天大学 适用于全光纤数字测斜仪的捷联惯性组合测量控制装置

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101900559A (zh) * 2009-11-06 2010-12-01 北京自动化控制设备研究所 一种捷联惯性导航系统双轴旋转调制方法
CN102322858A (zh) * 2011-08-22 2012-01-18 南京航空航天大学 用于地磁/捷联惯导组合导航系统的地磁匹配导航方法
CN103712622A (zh) * 2013-12-31 2014-04-09 清华大学 基于惯性测量单元旋转的陀螺漂移估计补偿方法及装置
CN103712622B (zh) * 2013-12-31 2016-07-20 清华大学 基于惯性测量单元旋转的陀螺漂移估计补偿方法及装置
CN103925917A (zh) * 2014-05-05 2014-07-16 上海新跃仪表厂 一种运载火箭姿态角速率信号的测量系统及方法
CN104764451A (zh) * 2015-04-23 2015-07-08 北京理工大学 一种基于惯性和地磁传感器的目标姿态跟踪方法
CN105737793B (zh) * 2016-05-11 2018-06-15 西安中星测控有限公司 滚转角测量单元及测量方法
CN105737793A (zh) * 2016-05-11 2016-07-06 西安中星测控有限公司 滚转角测量单元及测量方法
CN106403935B (zh) * 2016-08-24 2019-03-29 易文俊 一种弹载磁组合姿态测量装置的应用
CN106403935A (zh) * 2016-08-24 2017-02-15 易文俊 一种弹载磁组合姿态测量装置的应用
CN107883944B (zh) * 2016-09-29 2021-03-09 北京航空航天大学 一种捷联星敏感器实现间接敏感地平的导弹姿态机动方法
CN107883944A (zh) * 2016-09-29 2018-04-06 北京航空航天大学 一种捷联星敏感器实现间接敏感地平的导弹姿态机动方法
CN107063254A (zh) * 2016-12-27 2017-08-18 南京理工大学 一种陀螺地磁组合的姿态解算方法
CN107063254B (zh) * 2016-12-27 2023-08-25 南京理工大学 一种陀螺地磁组合的姿态解算方法
CN108196570A (zh) * 2017-12-26 2018-06-22 深圳市道通智能航空技术有限公司 一种无人机航向修正方法、装置和无人机
CN111366188A (zh) * 2018-12-26 2020-07-03 北京信息科技大学 一种面向动态环境力测量领域的数据采集存储装置
CN110017808A (zh) * 2019-03-26 2019-07-16 北京理工大学 利用地磁信息和加速计解算飞行器姿态的方法
CN110017808B (zh) * 2019-03-26 2021-08-27 北京理工大学 利用地磁信息和加速计解算飞行器姿态的方法
CN110765593B (zh) * 2019-10-09 2023-08-29 上海机电工程研究所 适用于便携式导弹弹上设备的评估方法及系统
CN110765593A (zh) * 2019-10-09 2020-02-07 上海机电工程研究所 适用于便携式导弹初始装订信息的评估方法及系统
CN111443735A (zh) * 2020-03-25 2020-07-24 浙江大华技术股份有限公司 一种车载云台摄像机姿态保持的方法、装置及系统
CN111443735B (zh) * 2020-03-25 2023-10-24 浙江大华技术股份有限公司 一种车载云台摄像机姿态保持的方法、装置及系统
CN111504256A (zh) * 2020-04-29 2020-08-07 中国北方工业有限公司 一种基于最小二乘法的滚转角实时估计方法
CN113587746A (zh) * 2021-10-08 2021-11-02 北京信息科技大学 基于地磁信息测量弹体的大跨度转速的方法、装置及系统

Also Published As

Publication number Publication date
CN101561281B (zh) 2011-04-13

Similar Documents

Publication Publication Date Title
CN101561281B (zh) 一种捷联磁惯性组合系统的工作方法
CN101561280B (zh) 一种捷联磁惯性组合系统
CN107655476A (zh) 基于多信息融合补偿的行人高精度足部导航算法
CN102980577B (zh) 一种微型捷联航姿系统及其工作方法
CN103090870B (zh) 一种基于mems传感器的航天器姿态测量方法
CN102445200A (zh) 微小型个人组合导航系统及其导航定位方法
CN103900571B (zh) 一种基于惯性坐标系旋转型捷联惯导系统的载体姿态测量方法
CN202904027U (zh) 一种自主室内定位系统
CN105890593A (zh) 一种mems惯性导航系统及基于该系统的轨迹重构方法
CN104864874B (zh) 一种低成本单陀螺航位推算导航方法及系统
CN1932444B (zh) 适用于高速旋转体的姿态测量方法
CN107390247A (zh) 一种导航方法、系统及导航终端
CN201955092U (zh) 一种基于地磁辅助的平台式惯性导航装置
CN202974288U (zh) 一种微型捷联航姿系统
CN104697526A (zh) 用于农业机械的捷联惯导系统以及控制方法
CN105973271A (zh) 一种混合式惯导系统自标定方法
CN106052682A (zh) 一种混合式惯性导航系统及导航方法
CN109282804B (zh) 一种单轴光纤陀螺寻北算法
CN103900611A (zh) 一种惯导天文高精度复合两位置对准及误差标定方法
CN103017764A (zh) 高速列车自主导航及姿态测量装置
CN102564452A (zh) 一种基于惯性导航系统的在线自主标定方法
CN105674815A (zh) 一种应用于高速旋转弹体的姿态测量装置
CN104697521A (zh) 一种采用陀螺冗余斜交配置方式测量高速旋转体姿态和角速度的方法
CN102072729A (zh) 用于测量旋转飞行器姿态的测量装置及其测量方法
CN102174885B (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
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110413

Termination date: 20130519