发明名称 提高捷联惯导系统姿态测量输出频率求姿态角的方法
摘要 本发明提出了一种通过提高捷联惯性导航系统姿态测量输出频率,利用欧拉角法求得姿态角的新方法。本发明提高了捷联惯性导航系统姿态测量输出频率,利用欧拉角法以求得姿态角。
申请公布号 CN103411608A 申请公布日期 2013.11.27
申请号 CN201310156778.7 申请日期 2013.04.26
申请人 哈尔滨工程大学 发明人 奔粤阳;杨娇龙;鲍桂清;高伟;李倩;于飞
分类号 G01C21/16(2006.01)I 主分类号 G01C21/16(2006.01)I
代理机构 代理人
主权项 1.一种提高捷联惯性导航系统姿态测量输出频率的方法,其特征是:一种提高捷联惯性导航系统姿态测量输出频率求姿态角的方法,其特征是: 步骤1、通过外部设备确定载体的初始位置参数与初始速度值; 步骤2、光纤陀螺捷联惯性导航系统进行初始对准,确定载体坐标系相对导航坐标系在k时刻的初始姿态<img file="DEST_PATH_FSB0000115082430000011.GIF" wi="108" he="52" />θ(k),γ(k);步骤3、利用步骤2所确定的初始姿态<img file="DEST_PATH_FSB0000115082430000012.GIF" wi="108" he="61" />θ(k),γ(k),建立载体坐标系相对于导航坐标系在k时刻的捷联矩阵T(k);其中,<img file="DEST_PATH_FSB0000115082430000013.GIF" wi="1716" he="288" />步骤4、根据光纤陀螺输出的载体坐标系相对于惯性坐标系在k时刻的角速度<img file="DEST_PATH_FSB0000115082430000014.GIF" wi="334" he="62" />利用欧拉角微分方程求解;其中,欧拉角微分方程为<img file="DEST_PATH_FSB0000115082430000015.GIF" wi="1680" he="309" />而<img file="DEST_PATH_FSB0000115082430000016.GIF" wi="101" he="215" />为机体坐标系相对于平台坐标系在k时刻的角速度向量<img file="DEST_PATH_FSB0000115082430000017.GIF" wi="81" he="67" /><img file="DEST_PATH_FSB0000115082430000018.GIF" wi="574" he="96" />步骤5、利用一阶解欧拉算法求得在<img file="DEST_PATH_FSB0000115082430000019.GIF" wi="137" he="105" />时刻的姿态<img file="DEST_PATH_FSB00001150824300000110.GIF" wi="541" he="127" />其中,一阶欧拉算法为:y<sub>i</sub>+1=y<sub>i</sub>+τf(ω<sub>i</sub>,y<sub>i</sub>,t<sub>i</sub>)。根据此公式,可得出<img file="DEST_PATH_FSB0000115082430000021.GIF" wi="725" he="96" /><img file="DEST_PATH_FSB0000115082430000022.GIF" wi="683" he="86" /><img file="DEST_PATH_FSB0000115082430000023.GIF" wi="730" he="93" />步骤6、根据步骤5所求得的<img file="DEST_PATH_FSB0000115082430000024.GIF" wi="585" he="100" />重复步骤2,再次建立载体坐标系相对于导航坐标系在<img file="DEST_PATH_FSB0000115082430000025.GIF" wi="111" he="98" />时刻的捷联矩阵<img file="DEST_PATH_FSB0000115082430000026.GIF" wi="197" he="97" />其中,<img file="DEST_PATH_FSB0000115082430000027.GIF" wi="1992" he="349" />步骤7、根据光纤陀螺输出载体坐标系相对于惯性坐标系在<img file="DEST_PATH_FSB0000115082430000028.GIF" wi="111" he="100" />时刻的角速度<img file="DEST_PATH_FSB0000115082430000029.GIF" wi="570" he="99" />重复步骤3,利用欧拉角微分方程求解。其中,欧拉角微分方程为: <img file="DEST_PATH_FSB00001150824300000210.GIF" wi="1770" he="422" />步骤8、利用二阶龙格-库塔法,将<img file="DEST_PATH_FSB00001150824300000211.GIF" wi="597" he="105" />的数值,更新为<img file="DEST_PATH_FSB00001150824300000212.GIF" wi="604" he="109" />其中,二阶龙格-库塔法为:在步骤中,所用的一阶欧拉法求得点(x<sub>i</sub>,y<sub>i</sub>)的斜率K<sub>1</sub>=f(ω<sub>i</sub>,y<sub>i</sub>,t<sub>i</sub>)的基础上,可得点t<sub>i+1</sub>处的预报值,即<img file="DEST_PATH_FSB00001150824300000213.GIF" wi="479" he="57" />及t<sub>i</sub>点得斜率值K<sub>1</sub>=f(ω<sub>i</sub>,y<sub>i</sub>,t<sub>i</sub>),再由点t<sub>i+1</sub>的预报值<img file="DEST_PATH_FSB00001150824300000214.GIF" wi="83" he="51" />求t<sub>i+1</sub>点的斜率值K<sub>2</sub>=f(ω<sub>i+1</sub>,y<sub>i+1</sub>,t<sub>i+1</sub>),然后求平均斜率<img file="DEST_PATH_FSB00001150824300000215.GIF" wi="351" he="110" />从而,<img file="DEST_PATH_FSB00001150824300000216.GIF" wi="475" he="99" />则由二阶龙格-库塔法可得到 <img file="DEST_PATH_FSB0000115082430000031.GIF" wi="1095" he="106" /><img file="DEST_PATH_FSB0000115082430000032.GIF" wi="1081" he="105" /><img file="DEST_PATH_FSB0000115082430000033.GIF" wi="1075" he="106" />步骤9、再利用步骤8所更新后,所得到的<img file="DEST_PATH_FSB0000115082430000034.GIF" wi="575" he="101" />数值,建立载体坐标系相对于导航坐标系在<img file="DEST_PATH_FSB0000115082430000035.GIF" wi="110" he="97" />时刻的捷联矩阵<img file="DEST_PATH_FSB0000115082430000036.GIF" wi="211" he="97" />重复步骤3至步骤8,以不断循环更新下一时刻的姿态和捷联矩阵,从而通过提高捷联惯性导航系统姿态测量输出频率能够求得姿态角。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号