发明名称 基于角速度的欧拉角任意步长正交级数指数型近似输出方法
摘要 本发明公开了一种基于角速度的欧拉角任意步长正交级数指数型近似输出方法,用于解决现有的飞行器机动飞行时欧拉角输出精度差的技术问题。技术方案是通过引入多个参数并将滚转、俯仰、偏航角速度采用改进的类似切比雪夫正交多项式的递推形式展开逼近,按照依次求解俯仰角、滚转角、偏航角,直接对欧拉角的表达式进行高阶逼近积分,使得欧拉角的求解按照超线性逼近,保证了确定欧拉角的时间更新迭代计算精度,从而提高了惯性设备输出飞行姿态的准确性。
申请公布号 CN102435192B 申请公布日期 2013.10.09
申请号 CN201110380055.6 申请日期 2011.11.25
申请人 西北工业大学 发明人 史忠科
分类号 G01C21/16(2006.01)I 主分类号 G01C21/16(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 1.一种基于角速度的欧拉角任意步长正交级数指数型近似输出方法,其特征在于包括以下步骤: 步骤1、(a)根据欧拉方程: <img file="FDA00003351809300011.GIF" wi="605" he="215" />式中:<img file="FDA00003351809300012.GIF" wi="41" he="42" /><img file="1.GIF" wi="16" he="20" />,Ψ分别指滚转角、俯仰角、偏航角;p,q,r分别为滚转角速度、俯仰角速度、偏航角速度;这三个欧拉角的计算按照依次求解俯仰角、滚转角、偏航角的步骤进行;滚转角速度、俯仰角速度、偏航角速度p,q,r的展开式分别为p(t)=pξ,q(t)=qξ,r(t)=rξ 其中 p=[p<sub>0</sub>  p<sub>1</sub>  …  p<sub>n-1</sub>  p<sub>n</sub>]q=[q<sub>0</sub>  q<sub>1</sub>  …  q<sub>n-1</sub>  q<sub>n</sub>] r=[r<sub>0</sub>  r<sub>1</sub>  …  r<sub>n-1</sub>  r<sub>n</sub>]ξ=[ξ<sub>0</sub>(t)  ξ<sub>1</sub>(t)  …  ξ<sub>n-1</sub>(t)  ξ<sub>n</sub>(t)]<sup>T</sup><img file="FDA00003351809300013.GIF" wi="1495" he="365" />为改进的类似切比雪夫正交多项式的递推形式,a为任意实数,T为采样周期; (b)俯仰角的时间更新求解式为: <img file="FDA00003351809300014.GIF" wi="1025" he="150" /><img file="FDA00003351809300015.GIF" wi="1177" he="150" /><img file="FDA00003351809300016.GIF" wi="1163" he="149" />式中: a<sub>1</sub>=(qζ)<sup>2</sup>+(rζ)<sup>2</sup>-(pζ)<sup>2</sup>a<sub>2</sub>=pΩr<sup>T</sup>a<sub>3</sub>=pΩq<sup>T</sup><img file="FDA00003351809300017.GIF" wi="576" he="80" /><img file="FDA00003351809300021.GIF" wi="710" he="90" /><img file="FDA00003351809300022.GIF" wi="985" he="90" /><img file="FDA00003351809300023.GIF" wi="742" he="111" /><img file="FDA00003351809300024.GIF" wi="803" he="109" /><img file="FDA00003351809300025.GIF" wi="937" he="90" /><img file="FDA00003351809300026.GIF" wi="585" he="90" /><img file="FDA00003351809300027.GIF" wi="1693" he="111" /><img file="FDA00003351809300028.GIF" wi="1651" he="108" /><img file="FDA00003351809300029.GIF" wi="778" he="110" /><img file="FDA000033518093000210.GIF" wi="1343" he="109" />步骤2、在已知俯仰角的情况下,滚转角的时间更新求解式为: <img file="FDA000033518093000211.GIF" wi="1515" he="150" /><img file="FDA000033518093000212.GIF" wi="876" he="149" /><img file="FDA000033518093000213.GIF" wi="1099" he="149" />其中 a<sub>4</sub>=(pζ)<sup>2</sup>+(rζ)<sup>2</sup>-(qζ)<sup>2</sup>a<sub>5</sub>=qΩp<sup>T</sup>a<sub>6</sub>=qΩr<sup>T</sup>步骤3、在俯仰角、滚转角已知情况下,偏航角的求解式为: <img file="FDA000033518093000214.GIF" wi="663" he="87" />式中: <img file="FDA000033518093000215.GIF" wi="830" he="121" /><img file="FDA000033518093000216.GIF" wi="393" he="121" />
地址 710072 陕西省西安市友谊西路127号