第16卷第lo期2008年10月光学精密工程OpticsandPrecisionEngineeringV01.16NO.10Oct.2008文章编号1004—924X(2008)10—1956—07捷联惯导系统姿态解算模块的实现杜海龙1’2,张荣辉1’2,刘平3,郑喜凤1,贾宏光1,马海涛4’1(1.中国科学院长春光学精密机械与物理研究所,吉林长春130033;2.中国科学院研究生院,北京100039;3.大庆职业学院,黑龙江大庆163255;4.吉林大学通信工程学院,吉林长春130012)摘要:提出了一种实现捷联惯性导航系统姿态解算模块的方法。基于DSP的硬件平台和四元数的数学平台,设计了捷联惯性导航系统的姿态解算模块。介绍了捷联惯性导航系统的工作原理和姿态解算的基本算法,并给出了四元数法的四阶龙格一库塔数值解法。设计了姿态解算模块的硬件电路和软件实现程序。实验测试结果表明,在增量角<5。的情况下,用四阶龙格一库塔法进行姿态解算,误差<0.0053%;应用TMS320C6713B进行硬件电路设计,每次解算时间<36肚s,能够满足捷联惯性导航系统对精度和速度的要求。关键词:捷联系统;惯性导航;姿态解算;四元数中图分类号:V448.22文献标识码:ARealizationofattitudealgorithmmoduleinstrapdowninertialguidancesystemDUHai—lon91'2,ZHANGRong—huil’2,IAUPing”,ZHENGXi—fen91,JIAHong—guan91,MAHal—ta04(1.ChangchunInstituteol厂Optics,FineChineseAcademyMechanicsandPhysi(‘S,ofSciences,Changchun130033,China;2.GraduateUniversityoJ’ChineseAcademyoJ’Sciences,Beijing100039,China;3.DaqingVocationalCollege,Daqing163255,China;4.SchoolofCornmunic’ationEngineering,JilinUniversity,Changc’hun130012,Chi*la)Abstract:AmethodBasedontorealizeattitudealgorithmmoduleinstrapdowninertialguidanceispresented.thedigitalsignalprocessorandquaternions,an.attitudealgorithmmoduleinstrapdowniner—tialguidanceisdesigned.Theworkingprincipleoftrapdowninertialguidance,theattitudealgorithmandthenumericalmethodThetestintroduced.Thecircuitandsoftwareofattitudealgorithmerrordesigned.resultsshowthattherelativecaniS1essthan0.0053%with4orderRunge~Kuttawhentheincrementislessthan5。;itmeetthedemandsoftrapdowninertialguidance,whenthealgorithmtimeislessthan36ptswithTMS320C6713B.Keywords:strapdownsystem;inertialguidance;attitudealgorithm;quaternion收稿日期:2008—08—07;修订日期:2008—09—05.基金项目:吉林省杰出青年基金资助项目(No.20060115);中国科学院“三期创新”平台资助项目万方数据 第10期杜海龙,等:捷联惯导系统姿态解算模块的实现引言捷联惯性导航技术是一门综合性技术,用于对运动体的姿态、速度和位置参数的确定,该项技术广泛应用于航天、航空、航海和大地测量等领域。由于捷联惯性导航不受外界电磁和气象环境干扰,不依赖于外界条件,并且具有良好的隐蔽性,真正地实现了自主导航,成为自主式测量的最佳手段,使得惯性导航技术在军事上具有特殊的应用价值。目前,在发达国家,捷联惯性导航系统广泛用于机器人、无人飞机和精确制导炸弹等飞行器的导航控制中。捷联惯导系统是通过提取陀螺和加速度计的测量值实时地计算姿态矩阵,进而从姿态矩阵的元素中获得运动载体的姿态和航向信息来进行导航计算的。姿态解算算法是捷联惯导系统算法的核心,也是影响捷联惯导系统精度的主要因素之一,因此设计和采用合理的姿态解算算法就成为需要研究的课题。目前,国内的飞行器采用的计算机大多是冯·诺伊曼结构的通用型微处理器,这些通用型微处理器的乘法用软件实现,常常需要若干个机器周期才能完成,而且由于同时运行多个进程,使数据处理速度较慢。并且采用这些微处理器的器件体积大,功耗大,电路结构复杂,可靠性低,成本高,无法更好地满足当代捷联惯性导航系统对小型化、高精度、实时性的要求。对于嵌入式数字信号处理系统,基于数字信号处理器(DSP)的实时信号处理系统具有速度快、可靠性高、成本低、开发周期短等特点。DSP是专为高速数字信号处理而设计的微处理器,其改进的哈佛结构、先进的多总线和多级流水线机制、专用的硬件乘法器、高效的指令集,使其易于处理嵌入式各种信号的滤波、信息融合算法,而不会牺牲系统的实时性。此外,DSP芯片中集成了RAM存储器和各种外设接口,便于系统的设计和调试。综上所述,DSP系统具有速度快、精度高、体积小、成本低、可靠性高以及抗干扰能力强等优点,可以满足捷联式惯性导航系统对硬件平台的要求H“。万 方数据2捷联惯性导航系统的工作原理捷联惯性导航系统的工作是把捷联惯性测量装置测量的运载体信息处理成导航所需要的各种信息,形成控制指令,控制运载体按预定轨道运行。惯性导航的基本原理是以牛顿力学定律为基础,在运载体内用加速度表和陀螺测量运载体加速度和角速度,通过积分运算获得运载体姿态、速度和位置的信息,利用这些信息进行导航。捷联惯性导航系统的关键是高精度的惯性测量仪表和高性能的计算机以及软件设计。捷联惯性导航系统的特点是把惯性测量元件加速度表和陀螺直接固连在运载体上,这些元件测量出沿载体坐标系三轴的运载体的角速度和线加速度,计算机实时计算出姿态矩阵,通过姿态矩阵把加速度计测量的载体坐标系的轴向加速度信息变换到导航坐标系。捷联惯性导航原理框图如图1所示。加速度计H姿态矩阵H导航计算卜.一载体位置速度陀螺仪H姿态矩阵计算姿态角计算卜载体姿态角图1捷联惯性导航原理框图Fig.1Schematicdiagramofstrapdowninertialguidance3捷联惯性导航姿态解算的基本算法捷联惯性导航的姿态解算主要是完成坐标系的变换。即将运载体相对于载体坐标系下测得的加速度和角速率转换到导航坐标系下,同时计算出运载体的姿态、速度和位置。常用的算法有欧拉角法、方向余弦法和四元数法。由于欧拉角法求解时,方程中存在奇点,所以不能用于全姿态飞行器上。方向余弦法虽然可以全姿态解算,但是由于计算量大,所有也不能用于飞行器姿态实时解算。现在常用的方法是四元数法,由于其计算无奇点,计算量小,所以常用于飞行器实时解算。1958光学精密工程第16卷3.1四元数法四元数(四维数)的概念是1843年由哈密顿首先提出的,他是代数学中的内容之一。近些年来,随着控制理论、惯性技术、计算技术,特别是捷联惯性导航技术的发展,为了更简便地描述刚体的角运动,设计控制系统,人们采用了四元数这个数学工具,用它来描述刚体角运动的3个欧拉角参数在设计控制系统时的不足∞]。四元数是指由1个实数单位A和3个虚数单位P,、P。、P。组成的数,表示为:g—A+Pzi+PzJ『+P。足。其中i,J『,七为3个虚数单位,也可以看成是三维空间的单位矢量。3.2基于四元数法的姿态解算算法‘4。83(1)在捷联惯性导航中,初始四元数确定计算如下:cOs譬cOs鲁cOs警+sin譬sin鲁sin警匿cOs警cOs鲁sin警一sin譬sin鲁cOs等cos譬s·’nioocos警+sin警cos譬sin等sin譬COs鲁cOs誓一cOs譬sin鲁sin警(1)式中,A、P.、P2、P。为四元数;A。、P1(o)、P:(0)、P。(o)为初始四元数;瓯、幽、托为初始姿态角。(2)四元数微分方程的四阶龙格一库塔数值解法如下:KIo—i1[一(u。(f)P1(f)一哟(f)P2(f)一叫:(f)P3(£)]Kll=百1[∽(£)A(f)+叫=(f)P2(£)一叫。(f)P{(f)]K12=i1[吣(f)A(f)~叫:(f)Pl(f)+∽(f)P{(f)]K13=百1[∽(f)A(f)+吣(f)Pl(f)一∽(f)P2(f)]r9、A。一A(£)+鉴#A。一P.(£)+鉴#(3)A。一P。(})+丝粤A。一Pm)+掣万 方数据‰÷[一叫√f+号)A,一山,(f+丢)A:(f)一叫:(件詈)A。(£)]%÷[∽(f+手)A,,+Ⅲ:(计手)A。(f)一叫,(f+等)A。(f)]%÷[甜,(,+手)一。一叫:(,T等)月,(,)+甜。(,+等)4。(,)]阮专[叫m+吾地,+吣(件等)Am)一叭(件等)A:(f)](4)Bo刊卅等B1_P1(卅等B2_P2(卅等(5)Ba=P=;㈤+等K。。=可1[一叫√f+等)B。一叫,(升号)B。(f)一叫:(件号)B。(f)]K。。=÷[∽(,+吾)B。,+甜:(}+号)B:(r)一甜,(f+吾)B。(f)]托:=丁1[叫沁+号)Bf,一甜m+号)B】(f)+∽(f+号)疡(f)]尥。=T1[埘m+号)BfJ+哟(f+等)B1(f)一∽(件丢)B:(f)](6)fCt,一A(£)+K。c,丁JCl—P“幻+K引T(7)lC2一P2(t)+K2Tlc。一P。(£)+K。。T虬。=÷[一∽(,+了1)c:一蛳(z+丁)G(,)一甜=(外丁)G(,)]‰=÷‰(f+眦。+“件丁)c加)一%(汁丁)C3“)]K舵:{[甜,(件丁)c。,一([J:(f+丁)cl(f)+∽(f+丁)c。(,)]’K。。=虿1㈧f+丁)c。+啪(f+丁)啪)一。(升丁)c2(f)](8)A(£+l、)一A(f)+÷(Kl。+2K2。,+2K={()+K¨)P。(f+T)一P。(f)+÷(Kll+2K21+2K3l+K。,)P2(£+T)一P2(£)+÷(Klb2+2K22+2K32+K;2)P。(f+T)一P。(£)+÷(K,。+2Kb23+2K。。+K。。)第10期杜海龙,等:捷联惯导系统姿态解算模块的实现FTllT1式中~,(f)为38轴行时刻角速度;叫。(t)为y轴”时刻角速度;叫:(t)为2轴”时刻角速度。叫。(t+TT13]7’33C鼍一l1、z1T22,ll,(11)÷)为T轴”时刻与”+1时刻角速度的中值;叫。『瓦。臼一T{。T;。J(12)T式中,b表示弹体坐标系,E表示导航坐标系。(4)飞行姿态角的计算如下:arcsin(Tl3(卵)),(£+i1)为y轴,2时刻与”+1时刻角速度的中T值;叫:(£+÷)为之轴”时刻与”+1时刻角速度式中,0为俯仰角。的中值。∽(f+T)为z轴”4-l时刻角速度;叫,(£+T)为y轴,z+1时刻角速度;叫:(£+T)为2轴钾+l时刻角速度。(3)姿态矩阵的计算如下:阱+薪一霞一篇式中,函为偏航角。2晒鹿柏众)芹+程一荔一篇2铂P。--AA)]2(忽A斗AA)l,芹+痨一西一鹾J(10)卿rctan(黑),一rctan(搿),T2一”T1,(13)(14)式中,),为滚转角。(5)飞行速度、位置的计算如下:(15)Q—l2(向忽一A)2(AA—hA)2(AA—向)式中,11.为采样周期;T:导航计算周期。简易表示为:陬2+Pi—Pi—Pi2(Pl户2--2p。)A2(P1P2+2p;)△∞2+P;一户;一p;2(P2P:{一Apl)L2(P·p。+Ap2)∞趣印旷时㈠趣m一△∞(16),,疰△∞式中,叫。为惯性器件的数字输出信号。速度计算:4系统设计与分析捷联惯性导航系统由数据采集模块、姿态解∑幽彬(j)∑△叫业(J)∑△毗(j)汁靴,(17)算模块、通信模块、电源模块和数据存储模块等组成‘“,其结构框图如图2所示。鹾k㈦憧置司l叫I川“k一㈦∥十T=}j+嘲∑∑△∽F(j)一寺∑△∽F(j)∑∑△∞妒(j)一寺∑△叫。。(j)∑∑△叫出(j)一寺∑△∞:,(j)Fig.2图2捷联惯性导航系统结构框图Structurechartofstrapdowninertialguidance(18)万方数据 光学精密工程第16卷4.1姿态解算模块的功能4.1.1模块的启动和自检测模块启动后,各部分的工作是否正常,要通过自检测程序加以检测,其中包括电源、惯性仪表和串行口等。4.1.2系统初始化系统初始化包括3项任务:(1)给定运载体的初始位置和初始姿态值;(2)在捷联惯性导航系统中,初始对准,确定姿态矩阵的初始值;(3)惯性仪表校正,输入运载体上惯性仪表标定系数,并将陀螺、加速度计的标定系数及补偿系数读入计算机,这是提高系统精度的重要保证。4.1.3惯性仪表的误差补偿为提高系统精度,首先对读人的数据进行中值滤波,然后对惯性仪表的误差进行补偿,计算机内通过软件实现误差补偿。4.1.4姿态矩阵计算姿态矩阵计算是捷联惯性导航算法中最重要的一部分,也是捷联系统所特有的。4.1.5导航计算导航计算就是把惯性测量装置测得的载体坐标系下的量值,变换到导航坐标系中,然后计算出运载体的姿态、速度和位置等导航所需要的信息。4.1.6导航和控制计算导航计算指的是按照导航规律算出控制运载体运动的控制指令;控制计算就是根据稳定控制系统设计所需要的控制算法,形成控制信号。4.2姿态解算模块的硬件电路设计姿态解算模块是捷联惯性导航系统的核心模块,其功能是将测得的载体坐标系下的加速度和图3姿态解算模块结构框图Fig.3Structurechartofattitudealgorithmmodule万 方数据角速度信息,变换到导航坐标系下,并计算出运载体的姿态、速度和位置信息,用于导航,其结构框图如图3所示。其中,DSP完成姿态解算的运算部分,FI.ASH完成程序的存储,EEPROM内部存放运载体的初始姿态和位置,状态显示部分完成系统自检及初始化的显示。4.3捷联惯性导航姿态解算模块时间分配根据姿态解算模块的功能,将其工作过程分为数据读入、数据解码、姿态解算、数据编码和数据输出,时间构成如表l所示。表1姿态解算模块时间构成(单位:ms)Timecompositionofattitudealgorithm(Unit:ms)数据读入时间:(19Byte×9bit/Byte)/38400bits/s≈4.5ms,数据输出时间:(36Byte×9bit/Byte)/1.5Mbit/s≈0.3ms.4.4捷联惯性导航姿态解算软件设计根据姿态解算模块的功能跏,设计了姿态解算的软件,流程图如图4所示。图4姿态解算流程图Fig.4FlowchartofattitudealgorithmTab.1第10期杜海龙,等:捷联惯导系统姿态解算模块的实现令整个解算运行120s,并计算解算的次数,5实验过程、结果及误差分析然后用所得到的时间除以循环次数,求得每次姿态解算所用时间。5.1求姿态解算的最大相对误差及解算时间5.2姿态解算的结果及误差分析设置时间问隔为0.01S,角速度分别为(1)增量角为5。,T一0.01S,叫一500。/s,如500。/s和50。/s,并计算相对误差,然后输出整个表2所示。计算过程中的最大误差,从而得到解算过程中的(2)增量角为0.5。,T一0.01s,叫一50。/s,如最大误差。表3所示。表2△=5。姿态解算结果Tab.2Resultsofattitudealgorithmas△一5。表3△=0.5。姿态解算结果Tab.3ResultsofattitudealgorithmA一0.5。其中:T为采样时间间隔,∞为旋转角速率。解算的实现方法。实验结果表明,四阶龙格一库塔通过数据对比,可以看出,提高计算阶数和减法进行姿态解算,误差d0.0053%,能够满足捷小解算步长,都可以减小解算误差。提高解算阶联惯性导航系统的精度要求;应用数,增加了相应的解算时间;减小步长,需要提高TMS320C6713B进行硬件电路设计,每次解算时数据的采样频率。间<36“s,能够满足捷联惯性导航的实时性要求。提高计算阶数和提高采样频率,可以减小误6结论差,但是增加了DSP的计算量,必须选择更高主频的DSP,增加了成本。本文提出了一种捷联惯性导航系统的姿态参考文献:[1]陆元九.惯性器件[M].北京:中国宇航出版社,1990.I。uYJ.InertialDevice[M].Beijing:ChinaAerospacePress,1990.(inChinese)[23SPRS249B.TMS320C6713BFloating—pointDigitalSignalProcessor[G].TexasInstrument,2006.[3]邓正隆.惯性技术[M].哈尔滨:哈尔滨工业大学出版社,2006.DENGZHI..Inertial孔(hnology[M].Harbin:HarbinInstituteofTechnologyPress,2006.(inChinese)E43于本水.防空导弹总体设计[M].北京:中国宇航出版社,1982.YUBSH.OverallDesignoJ、Anti—aircraJ、tMissile[M].Beijing:ChinaAerospacePress,1982.(inChinese)[5]王玫。张炎华.捷联惯性导航的工程实现[J].哈尔滨商业大学学报,2005,21(4):470—474.WANGM,ZHANGYH.Studyengineeringimplementationofstrapdowninertialnavigation[J].Journal0,万 方数据1962光学精密工程第16卷HarbinUniversitY0rCOl71D!PF(P(NaturalS(、iem‘P^Edition),2005,21(4):470—474.(inChinese)[6]赵峰,袁东风,张海霞,等.多DSP图像压缩实时并行处理系统[J].光学精密工程,2007,15(9):1451-1455.ZHA()F,YUANDF.ZHANGHX,ela1..Multi—DSPreal—timeparallelprocessingsystemforimagecompression[J].Opt.Precision[7]Eng.,2007,15(9):14511455.(inChinese)崔中兴.惯性导航系统[M].北京:国防工业出版礼,1982.CUIZHx.InertialNavigationSystem[M].Bcijing:NationalDefenseIndustryPress,l982.(inChinese)4(3):490491.[8]翟林培,刘明,修吉宏.考虑飞机姿态角时倾斜航空相机像移速度计算[J].光学精密工程,2006,1ZAII.P,IAUM,XIUJH.CalculationofimagemotionvelocityconsideringairplaneEng.,2006,14(3):490—494.(inChinese)gestureangleinobliqueaerialcameraEJ].Opt.Pr“ision[9]张鹏飞,汤建勋.龙兴武.基于DSI,的捷联惯导计算机系统的设计与开发[J].计算机工程,2007,33(2):244ZHANGPF,TANGjX,I,()NGforX248.DSPW.DesignanddevelopmentofthenavigationcomputersystembasedChinese)SINS[J].ComputerEnl4ineering,2007,33(2):244—248.(in[io]刘延斌,许晖,徐长吉,等.四元数矩阵在航天相机像面位置求解中的应用[J].光学精密工程,2004,12(2):136140.I。IUYB,XUH,XUCHJ,ela1..ApplicationofquaternionmatrixinresolutionofaerospacecameraimEng.,2004,12(2):136140.(inChinese)ageplaneposition[J].Opt.PrP(ision作者简介:杜海龙(1983一),男,博=[_=研究生,主要从事数字信号处理方面的研究。E-mail:dhlmao@126.com导师简介:贾宏光(1971一),男,研究员,博上生导师,主要从事数字信号处理及精密机械方面的研究。E—mail:jiahg@Comp.ac.cn●下期预告基于空间编码投影法测量三维大物体面形刘晨,金施群,费业泰,唐黎莺,姜晨(合肥工业大学仪器科学与光电工程学院,安徽合肥230009)三维图像拼接技术是实现大型物体形貌测量的关键技术。在图像控制点约束的基础上,提出了一种实现三维图像拼接的新方法,即将一伪随机空间编码投射到一大型物体表面,然后用CCI)在不同视角拍摄不同视场中的物体,利用重叠区域相同的某一窗口特征点进行配准,采用四元组法求取坐标转换矩阵,从而实现大范围自由曲面三维形貌的测量。实验结果表明:基于空问编码的图像拼接新方法,可以实现大尺度三维面形的测量,拼接数据相对误差<o.8%。万方数据 捷联惯导系统姿态解算模块的实现
作者:作者单位:
杜海龙, 张荣辉, 刘平, 郑喜凤, 贾宏光, 马海涛, DU Hai-long, ZHANG Rong-hui, LIU Ping, ZHENG Xi-feng, JIA Hong-guang, MA Hai-tao
杜海龙,张荣辉,DU Hai-long,ZHANG Rong-hui(中国科学院,长春光学精密机械与物理研究所,吉林,长春,130033;中国科学院,研究生院,北京,100039), 刘平,LIU Ping(大庆职业学院,黑龙江,大庆,163255), 郑喜凤,贾宏光,ZHENG Xi-feng,JIA Hong-guang(中国科学院,长春光学精密机械与物理研究所,吉林,长春,130033), 马海涛,MA Hai-tao(吉林大学,通信工程学院,吉林,长春,130012;中国科学院,长春光学精密机械与物理研究所,吉林,长春,130033)
光学精密工程
OPTICS AND PRECISION ENGINEERING2008,16(10)1次
刊名:英文刊名:年,卷(期):被引用次数:
1.陆元九 惯性器件 1990
2.SPRS24 9B.TMS320C6713B Floating-point Digital Signal Processor 20063.邓正隆 惯性技术 20064.于本水 防空导弹总体设计 1982
5.王玫;张炎华 捷联惯性导航的工程实现[期刊论文]-哈尔滨商业大学学报 2005(04)6.赵峰;袁东风;张海霞 多DSP图像压缩实时并行处理系统[期刊论文]-光学精密工程 2007(09)7.崔中兴 惯性导航系统 1982
8.翟林培;刘明;修吉宏 考虑飞机姿态角时倾斜航空相机像移速度计算[期刊论文]-光学精密工程 2006(03)9.张鹏飞;汤建勋;龙兴武 基于DSP的捷联惯导计算机系统的设计与肝发[期刊论文]-计算机工程 2007(02)10.刘延斌;许晖;徐长吉 四元数矩阵在航天相机像面位置求解巾的应用[期刊论文]-光学精密工程 2004(02)
1. 张荣辉.贾宏光.陈涛.张跃.ZHANG Rong-hui.JIA Hong-guang.CHEN Tao.ZHANG Yue 基于四元数法的捷联式惯性导航系统的姿态解算[期刊论文]-光学精密工程2008,16(10)
2. 项凤涛.王正志.吴第旻.岳达 捷联系统四元数姿态解算的精细积分法[期刊论文]-四川兵工学报2010,31(5)3. 邢质皙.王爱民.Xing Zhixi.Wang Aimin 捷联惯导系统姿态算法比较[期刊论文]-舰船电子工程2009,29(7)4. 张春慧.吴简彤.何昆鹏.郭新伟 四阶龙格-库塔法在捷联惯导系统姿态解算中的应用[期刊论文]-声学与电子工程2005(1)
5. 张春慧.吴简彤.何昆鹏.周雪梅.ZHANG Chun-hui.WU Jian-tong.HE Kun-peng.ZHOU Xue-mei 四阶龙格-库塔法在捷联惯导系统姿态解算中的应用[期刊论文]-应用科技2005,32(6)
1.小型无人机CAN总线信号记录仪[期刊论文]-数据采集与处理 2009(z1)
本文链接:http://d.g.wanfangdata.com.cn/Periodical_gxjmgc200810028.aspx