发明名称 一种用于捷联惯导软件测试的轨迹数据生成方法
摘要 本发明公开了一种用于捷联惯导软件测试的轨迹数据生成方法;本方法首先分别设计载体的运动路径和姿态演变,得到位置和姿态的变化函数;然后通过上述变化函数及其导数,计算姿态矩阵、速度和加速度等参数;最后,根据上述结果求解惯性敏感单元的输出数据。测试惯导软件时,将惯性敏感单元数据作为输入数据流,并对比惯导软件输出的导航参数与本方法输出的标准导航参数,从而得出惯导软件的计算误差。本方法简单高效,能够灵活产生各种复杂且高精度的仿真轨迹数据,从而为评测捷联惯导软件的正确性、可靠性和算法精度起到重要作用,且在实验室基于光纤陀螺和石英加速度计的捷联惯性导航系统及其导航软件的研究中取得了良好的效果。
申请公布号 CN102607591A 申请公布日期 2012.07.25
申请号 CN201210045683.3 申请日期 2012.02.27
申请人 浙江大学 发明人 楼奇哲;周一览;律新伟;舒晓武
分类号 G01C25/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 周烽
主权项 1.一种用于捷联惯导软件测试的轨迹数据生成方法,其特征在于,它包括位置相关数据生成、姿态相关数据生成和惯性敏感单元数据生成等三部分内容;其中,所述位置相关数据生成的具体步骤如下:(1)设定位置关于时间<img file="877930DEST_PATH_IMAGE002.GIF" wi="9" he="17" />的变化函数并求解其一阶和二阶导数,载体位置即载体相对于地球的高度<img file="2012100456833100001DEST_PATH_IMAGE003.GIF" wi="14" he="20" />、经度<img file="216900DEST_PATH_IMAGE004.GIF" wi="16" he="20" />和纬度<img file="DEST_PATH_IMAGE005.GIF" wi="16" he="18" />;位置变化函数及其导数的一般公式如下:<img file="658377DEST_PATH_IMAGE006.GIF" wi="66" he="76" />;<img file="DEST_PATH_IMAGE007.GIF" wi="69" he="78" />;<img file="632149DEST_PATH_IMAGE008.GIF" wi="72" he="78" />;(2)计算载体相对于地球的速度在导航系中的值<img file="DEST_PATH_IMAGE009.GIF" wi="22" he="26" />采用地球椭球模型时,记地球赤道半径为<img file="694521DEST_PATH_IMAGE010.GIF" wi="20" he="25" />,地球子午面曲率半径为<img file="DEST_PATH_IMAGE011.GIF" wi="24" he="25" />,地球横向曲率半径为<img file="753744DEST_PATH_IMAGE012.GIF" wi="24" he="26" />,地球椭球度为<img file="DEST_PATH_IMAGE013.GIF" wi="13" he="16" />,可得:<img file="61229DEST_PATH_IMAGE014.GIF" wi="169" he="76" />;<img file="DEST_PATH_IMAGE015.GIF" wi="180" he="54" />;采用地球的圆球模型时,只需令<img file="594234DEST_PATH_IMAGE016.GIF" wi="76" he="26" />,且有<img file="DEST_PATH_IMAGE017.GIF" wi="38" he="25" />,下同;(3)计算载体相对于地球的加速度在导航系中的值<img file="328972DEST_PATH_IMAGE018.GIF" wi="24" he="36" />:<img file="DEST_PATH_IMAGE019.GIF" wi="452" he="97" />;<img file="875491DEST_PATH_IMAGE020.GIF" wi="192" he="65" />;(4)计算重力加速度在导航系中的值<img file="DEST_PATH_IMAGE021.GIF" wi="21" he="25" />:<img file="783404DEST_PATH_IMAGE022.GIF" wi="112" he="30" />;重力加速度模型需要考虑纬度和高度的情况,一般公式为:<img file="DEST_PATH_IMAGE023.GIF" wi="338" he="25" />;(5)综合步骤(1)至(4),位置变化函数及其导数、载体相对于地球的速度在导航系中的值<img file="731769DEST_PATH_IMAGE009.GIF" wi="22" he="26" />、载体相对于地球的加速度在导航系中的值<img file="73626DEST_PATH_IMAGE018.GIF" wi="24" he="36" />以及重力加速度在导航系中的值<img file="169758DEST_PATH_IMAGE021.GIF" wi="21" he="25" />共同组成了位置相关数据;所述姿态相关数据生成的具体步骤如下:(1)设定姿态角关于时间<img file="881362DEST_PATH_IMAGE002.GIF" wi="9" he="17" />的变化函数并求解其一阶导数:载体的姿态即载体系相对于导航系的俯仰角<img file="949812DEST_PATH_IMAGE024.GIF" wi="14" he="20" />、横滚角<img file="DEST_PATH_IMAGE025.GIF" wi="14" he="18" />和航向角<img file="26353DEST_PATH_IMAGE026.GIF" wi="17" he="18" />;姿态角变化函数及其导数的一般公式如下:<img file="DEST_PATH_IMAGE027.GIF" wi="69" he="76" />;<img file="771631DEST_PATH_IMAGE028.GIF" wi="73" he="78" />;(2)计算导航系到载体系的变换矩阵<img file="DEST_PATH_IMAGE029.GIF" wi="22" he="26" />:<img file="224609DEST_PATH_IMAGE030.GIF" wi="485" he="76" />;(3)计算地球自转角速率在导航系中的值<img file="DEST_PATH_IMAGE031.GIF" wi="22" he="26" />,其中<img file="147566DEST_PATH_IMAGE032.GIF" wi="22" he="26" />表示地球自转角速率常数:<img file="DEST_PATH_IMAGE033.GIF" wi="108" he="76" />;(4)计算导航系相对于地球的转动速率在导航系中的值<img file="395007DEST_PATH_IMAGE034.GIF" wi="25" he="26" />,其中<img file="DEST_PATH_IMAGE035.GIF" wi="26" he="26" />和<img file="636370DEST_PATH_IMAGE036.GIF" wi="26" he="28" />即为向量<img file="689777DEST_PATH_IMAGE009.GIF" wi="22" he="26" />在导航系x轴和y轴的分量:<img file="DEST_PATH_IMAGE037.GIF" wi="136" he="156" />;(5)求解载体系相对于导航系的转动角速率在载体系中的值<img file="467240DEST_PATH_IMAGE038.GIF" wi="26" he="26" />:<img file="DEST_PATH_IMAGE039.GIF" wi="230" he="76" />;(6)综合步骤(1)至(5),姿态角变化函数及其导数、导航系到载体系的变换矩阵<img file="885583DEST_PATH_IMAGE029.GIF" wi="22" he="26" />、地球自转角速率在导航系中的值<img file="443604DEST_PATH_IMAGE031.GIF" wi="22" he="26" />、导航系相对于地球的转动速率在导航系中的值<img file="972805DEST_PATH_IMAGE034.GIF" wi="25" he="26" />、载体系相对于导航系的转动角速率在载体系中的值<img file="932671DEST_PATH_IMAGE038.GIF" wi="26" he="26" />共同组成了姿态角相关数据;所述惯性敏感单元数据生成的具体步骤如下:(1)陀螺仪输出数据生成:陀螺仪输出角速率<img file="23380DEST_PATH_IMAGE040.GIF" wi="24" he="26" />为载体系相对于惯性系的转动速率在载体系中的投影,计算方法为:<img file="DEST_PATH_IMAGE041.GIF" wi="158" he="26" />;(2)加速度计输出数据生成:加速度计的输出称为比力<img file="740800DEST_PATH_IMAGE042.GIF" wi="22" he="26" />,计算方法为:<img file="DEST_PATH_IMAGE043.GIF" wi="231" he="49" />;(3)综合步骤(1)至(2),陀螺仪输出数据<img file="136009DEST_PATH_IMAGE040.GIF" wi="24" he="26" />和加速度计输出数据<img file="888065DEST_PATH_IMAGE042.GIF" wi="22" he="26" />共同组成了惯性敏感单元数据。
地址 310058 浙江省杭州市西湖区余杭塘路388号