发明名称 一种用于天线跟踪稳定平台的捷联惯导初始对准方法
摘要 本发明属于惯性导航技术领域,为进一步提高捷联惯导系统初始对准的精度与速度,本发明提出了一种应用于天线跟踪稳定平台的捷联惯导初始对准方法,并建立以速度误差为状态变量的卡尔曼滤波状态方程及速度误差为观测量的观测方程,通过卡尔曼滤波技术估计出载体失准角并反馈到系统中完成系统的初始对准。本发明主要优点在于在不改变系统现有硬件结构条件下,充分利用了天线跟踪稳定平台的结构特点进行转动基座初始对准,有利于提高捷联惯导系统的可观测性,从而有效地提高了捷联惯导系统初始对准的对准精度。
申请公布号 CN103557876B 申请公布日期 2016.01.20
申请号 CN201310569032.9 申请日期 2013.11.15
申请人 山东理工大学 发明人 张华强;许敬;李东兴
分类号 G01C25/00(2006.01)I;G01C21/16(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 代理人
主权项 一种用于天线跟踪稳定平台的捷联惯导初始对准方法,其特征在于:该方法的实现包括以下步骤:步骤一、将捷联惯导安装到天线座旋转底盘上,捷联惯导到天线座底盘旋转中心的距离为r;步骤二、利用全球定位系统GPS确定载体所在位置的经度λ、纬度<i>L</i>;步骤三、根据所确定的经度、纬度信息以及惯导系统惯性器件的输出进行粗对准,得到惯性测量单元姿态矩阵的近似估计值;步骤四、启动天线座底盘按照逆时针方向旋转,旋转角速率记为ω,并控制天线座底盘旋转整数圈,旋转圈数记为n;步骤五、按照步骤三所确定的旋转方案旋转的同时,利用粗对准得到的姿态矩阵、惯性测量单元中陀螺实时输出的角速度信息和加速度计实时输出的比力信息,以更新惯性测量单元的姿态矩阵;步骤六、通过卡尔曼滤波完成捷联惯导转动基座对准,具体步骤如下:<img file="dest_path_image001.GIF" wi="16" he="24" />转动基座捷联惯导误差建模首先,建立旋转捷联惯导系统的平台误差角方程,旋转捷联惯导从惯性测量单元中获取角速度和比力信息,经过<img file="dest_path_image002.GIF" wi="24" he="26" />可实现惯性器件敏感系s到载体坐标系b的转换,转换后的捷联解算流程与常规的捷联惯导系统是一致的,因此,其平台误差角方程形式仍为:<img file="dest_path_image003.GIF" wi="269" he="26" />式中,<img file="dest_path_image004.GIF" wi="24" he="21" />为失准角矢量,<img file="dest_path_image005.GIF" wi="24" he="26" />为地球自转角速率在导航坐标系下的投影,<img file="dest_path_image006.GIF" wi="26" he="26" />为导航坐标系相对于地球坐标系在导航坐标系下的投影,<img file="dest_path_image007.GIF" wi="33" he="22" />为陀螺等效到导航坐标系下的误差;其次,建立旋转捷联惯导系统的速度误差方程,由比力方程可得速度误差方程为:<img file="dest_path_image008.GIF" wi="408" he="26" />式中,<img file="dest_path_image009.GIF" wi="24" he="22" />为载体在导航坐标系下的速度,<img file="dest_path_image010.GIF" wi="24" he="25" />为加速度计输出的比力;最后,建立旋转捷联惯导系统的位置误差方程,惯导系统位置误差方程为:<img file="dest_path_image011.GIF" wi="272" he="150" />式中,<img file="dest_path_image012.GIF" wi="28" he="25" />和<img file="dest_path_image013.GIF" wi="29" he="26" />分别为东向和北向速度误差,<img file="dest_path_image014.GIF" wi="18" he="25" />和<img file="dest_path_image015.GIF" wi="19" he="26" />分别为东向和北向速度,<img file="dest_path_image016.GIF" wi="25" he="25" />为地球子午圈主曲率半径,<img file="dest_path_image017.GIF" wi="24" he="25" />为地球卯酉圈主曲率半径,<img file="dest_path_image018.GIF" wi="14" he="19" />为高度,<img file="dest_path_image019.GIF" wi="16" he="18" />为载体所在地球纬度;对于旋转基座捷联惯导系统初始对准来讲,可以不考虑<img file="dest_path_image020.GIF" wi="14" he="14" />轴方向的加速度和速度,即可以认为<img file="dest_path_image021.GIF" wi="44" he="26" />,<img file="dest_path_image022.GIF" wi="44" he="25" />,同时捷联惯导系统本身不存在相对于地球的线速度,因此可以略去以上惯导系统误差方程中的有关项,可得旋转基座捷联惯导初始对准系统的误差方程如下:捷联惯导平台误差角简化方程:<img file="dest_path_image023.GIF" wi="326" he="145" />捷联惯导速度误差简化方程:<img file="dest_path_image024.GIF" wi="265" he="54" />捷联惯导位置误差简化方程:<img file="dest_path_image025.GIF" wi="272" he="150" />式中,<img file="dest_path_image026.GIF" wi="18" he="25" />、<img file="dest_path_image027.GIF" wi="18" he="26" />和<img file="dest_path_image028.GIF" wi="18" he="25" />分别为东向、北向和天向平台失准角,<img file="dest_path_image029.GIF" wi="22" he="25" />为地球自转角速率,<img file="dest_path_image030.GIF" wi="32" he="25" />、<img file="dest_path_image031.GIF" wi="33" he="26" />、<img file="dest_path_image032.GIF" wi="32" he="25" />等效为导航坐标系东北天方向陀螺的误差;惯性器件随机误差在导航坐标系以及惯性器件敏感坐标系之间的关系如下:<img file="dest_path_image033.GIF" wi="112" he="26" />、<img file="dest_path_image034.GIF" wi="106" he="26" />,其中,<img file="253521dest_path_image007.GIF" wi="33" he="22" />和<img file="dest_path_image035.GIF" wi="31" he="25" />分别为陀螺和加速度计在导航坐标系下的随机误差,<img file="dest_path_image036.GIF" wi="24" he="26" />为载体坐标系到导航坐标系的转换矩阵,<img file="dest_path_image037.GIF" wi="24" he="26" />为惯性器件敏感系到载体坐标系的转换矩阵,<img file="dest_path_image038.GIF" wi="32" he="22" />和<img file="dest_path_image039.GIF" wi="29" he="25" />分别为陀螺和加速度计在惯性器件敏感坐标系下的随机误差;<img file="dest_path_image040.GIF" wi="16" he="24" />卡尔曼滤波器状态方程与观测方程的建立捷联惯导在初始对准时,移动载体相对于地理坐标系是静止的,由于地理位置信息精确已知,因此可以忽略误差模型中的位置误差;同时考虑到惯导系统的垂直通道不稳定,因此初始对准状态方程中只采用速度误差和加速度误差的水平分量作为状态量,此时系统状态变量的维数可降为10维,惯导系统在导航坐标系下的初始对准误差模型为:<img file="dest_path_image041.GIF" wi="92" he="22" />式中,<img file="dest_path_image042.GIF" wi="19" he="18" />为系统的状态向量,<img file="dest_path_image043.GIF" wi="18" he="18" />为状态转移矩阵,<img file="dest_path_image044.GIF" wi="22" he="20" />为系统噪声向量;将这些参数写成矩阵形式如下:<img file="dest_path_image045.GIF" wi="364" he="26" /><img file="dest_path_image046.GIF" wi="430" he="150" />式中,<img file="dest_path_image047.GIF" wi="29" he="25" />和<img file="dest_path_image048.GIF" wi="31" he="25" />分别为东向和北向速度误差,<img file="dest_path_image049.GIF" wi="19" he="25" />、<img file="dest_path_image050.GIF" wi="21" he="25" />和<img file="dest_path_image051.GIF" wi="21" he="25" />分别为东向、北向和天向平台失准角,<img file="dest_path_image052.GIF" wi="22" he="25" />和<img file="dest_path_image053.GIF" wi="24" he="26" />分别为<img file="dest_path_image054.GIF" wi="14" he="16" />轴和<img file="dest_path_image055.GIF" wi="16" he="18" />轴加速度计偏置,<img file="dest_path_image056.GIF" wi="18" he="25" />、<img file="dest_path_image057.GIF" wi="19" he="26" />和<img file="dest_path_image058.GIF" wi="18" he="25" />分别为<img file="613456dest_path_image054.GIF" wi="14" he="16" />轴、<img file="862035dest_path_image055.GIF" wi="16" he="18" />轴和<img file="dest_path_image059.GIF" wi="14" he="14" />轴陀螺常值漂移,<img file="dest_path_image060.GIF" wi="93" he="25" />,<img file="dest_path_image061.GIF" wi="96" he="25" />,<img file="dest_path_image062.GIF" wi="22" he="26" />为姿态矩阵<img file="dest_path_image063.GIF" wi="41" he="26" />元素,即<img file="dest_path_image064.GIF" wi="122" he="28" />;取和作为卡尔曼滤波器的观测量时,系统观测方程为:<img file="dest_path_image065.GIF" wi="86" he="20" />式中,<img file="dest_path_image066.GIF" wi="18" he="18" />为系统观测向量,<img file="dest_path_image067.GIF" wi="109" he="28" />为观测矩阵,<img file="dest_path_image068.GIF" wi="17" he="20" />为观测噪声,服从<img file="dest_path_image069.GIF" wi="56" he="22" />分布,<img file="dest_path_image070.GIF" wi="18" he="18" />为观测噪声协方差阵;将所求得的状态方程和观测方程进行卡尔曼滤波,估计出姿态角的误差和惯性器件的零偏,即可完成初始对准。
地址 255086 山东省淄博市高新技术产业开发区高创园D座1012室
您可能感兴趣的专利