专利名称：基于角速度的欧拉角埃米特近似输出方法
技术领域：
本发明涉及一种飞行器机动飞行姿态确定方法，特别是涉及一种基于角速度的欧拉角埃米特近似输出方法。
背景技术：
惯性设备在运动体导航和控制中具有重要作用；刚体运动的加速度、角速度和姿态等通常都依赖于惯性设备输出，因此提高惯性设备的输出精度具有明确的实际意义；在惯性设备中，加速度采用加速度计、角速度采用角速率陀螺直接测量方式，刚体的姿态精度要求很高时如飞行试验等采用姿态陀螺测量，但在很多应用领域都有角速度等测量直接解算输出；主要原因是由于动态姿态传感器价格昂贵、体积大，导致很多飞行器采用角速率陀螺等解算三个欧拉角，使得姿态时间更新输出成为导航等核心内容，也使其成为影响惯导系统精度的主要因素之一，因此设计和采用合理的姿态时间更新输出方法就成为研究的热点课题；从公开发表的文献中对姿态输出主要基于角速度采用欧拉方程直接近似法或采用近似龙格库塔方法解算(孙丽、秦永元，捷联惯导系统姿态算法比较，中国惯性技术学报， 2006, Vol. 14(3) :6-10 ；Pu Li, Wang TianMiao,Liang JianHong, Wang Song, An Attitude Estimate Approach using MEMS Sensors for Small UAVs,2006, IEEE International Conference on Industrial Informatics, 1113-1117)；由于欧拉方程中三个欧拉角互相耦合，属于非线性微分方程，在不同初始条件和不同飞行状态下的误差范围不同，难以保证实际工程要求的精度。发明内容
为了克服现有的飞行器机动飞行时欧拉角输出精度差的问题，本发明提供一种基于角速度的欧拉角埃米特近似输出方法。该方法通过引入多个参数并将滚转、俯仰、偏航角速度按照埃米特正交多项式展开，通过按照依次求解俯仰角、滚转角、偏航角，直接对欧拉角的表达式进行高阶逼近积分，使得欧拉角的求解按照超线性逼近，从而可以保证确定欧拉角的时间更新迭代计算精度和惯性单元的输出精度。
本发明解决其技术问题所采用的技术方案是一种基于角速度的欧拉角埃米特近似输出方法，其特点是包括以下步骤
1、(a)根据欧拉方程
权利要求
1. 一种基于角速度的欧拉角埃米特近似输出方法，其特征在于包括以下步骤 步骤1、(a)根据欧拉方程\β^( βα\φ+ηχ^,φ)Ι cos; < 表 φ^φ-ν ηφ￠^ /^+(gsin 炉+AXIS 炉)tan;式中φ,θ,ψ分别指滚转、俯仰、偏航角；p，q，r分别为滚转、俯仰、偏航角速度；全文参数定义相同；这三个欧拉角的计算按照依次求解俯仰角、滚转角、偏航角的步骤进行；滚转、俯仰、偏航角速度P，q，r的展开式分别为p(t) = [p0 P1 L Pn^1 pn] [l0(t) I1WL ξ η_! (t) ln(t)]T q(t) = [q0 Q1 L qn_! qn] [ξ0( ) I1^L ξ (t) ln(t)]T r(t) = [r0 Γι L Iv1 rn] [ξ0( ) I1WL ξ (t) ln(t)]T 其中^1(I) = 2t ξ2{ ) = Α 2-1 4(/) = 8/3-12/= 16/4-48/2+ 12i = 2，3，L，n-lξ5( ) = 32/5 -160/3 +120/ ξ6 (t) = 64/6 - 480/4 + 120t-120 M^,,(/) = 2^,.(/)-2^(/)为埃米特正交多项式的递推形式； (b)俯仰角的时间更新求解式为式中 1 = 1+ [p0 P1 L P^1 Ρ ]|^ + 1>Τ[ξ( )ξΤ( )]^ Η T [p0 P1 L ρη_, ρ ]τ+ [9。 L L qn_, 9 ]|^ + 1>Τ[ξ( )ξΤ( )]^ Η τ Iq0 q, L qn_, qn ]Τ+ |>。、 L r _t rn + [l(t)lT (t)]dtH T [r0 rt L r _t r ]T-[P。P1 L pn_, p ]H [“⑴ ξ,α) L ξη( ) “ + 1 ⑴]T.[“⑴L ξη( ) “ + 1( )] Ητ[ρ。P1 L ρη_, ρ ]Τ-[9ο 9ι L qn_, g ]H [“⑴ ξ,α) L ξη( )+ 1(0]Τ [^+"7 '[“⑴L ξη( ) +1( )] H τ Iq0 q, L qn_, qn f- |>。L r _t r ]H [f0(i) ξ,α) L ξη( )+ 1( )Γ [" + 1)Τ ■[“⑴ ξλ( ) L ξ Jt)+ 1(0] H τ [r0 rt L r _t r ]τ+ 0.25{[p0 P1 L Prl] ρ ]Η[“( )L “⑴+ 1( )Γ [* + 1)Τ}2+ 0.25{[g。qi L I1 g ]H [“⑴。⑴ L “⑴ “ + 1 ⑴]Τ}2-0.2 5{[r0 rt L r _t r ] H [f 0 ( )。⑴ L “⑴+ 1(0]Τ [7"'τ (k+l)T 2 =L q“ qn]H^0(t) m L ξ州 ξη+ι( )]虹
全文摘要
本发明公开了一种基于角速度的欧拉角埃米特近似输出方法，用于解决现有的飞行器机动飞行时欧拉角输出精度差的技术问题。技术方案是通过引入多个参数并将滚转、俯仰、偏航角速度按照埃米特正交多项式展开，按照依次求解俯仰角、滚转角、偏航角，直接对欧拉角的表达式进行高阶逼近积分，使得欧拉角的求解按照超线性逼近，保证了确定欧拉角的时间更新迭代计算精度，从而提高了惯性设备输出飞行姿态的准确性。
文档编号G06F17/11GK102495827SQ20111038797
公开日2012年6月13日 申请日期2011年11月30日 优先权日2011年11月30日
发明者史忠科 申请人:西北工业大学