CN102359790B - 一种刚体空间运动状态的傅里埃输出方法 - Google Patents
一种刚体空间运动状态的傅里埃输出方法 Download PDFInfo
- Publication number
- CN102359790B CN102359790B CN 201110280706 CN201110280706A CN102359790B CN 102359790 B CN102359790 B CN 102359790B CN 201110280706 CN201110280706 CN 201110280706 CN 201110280706 A CN201110280706 A CN 201110280706A CN 102359790 B CN102359790 B CN 102359790B
- Authority
- CN
- China
- Prior art keywords
- sin
- centerdot
- cos
- nωt
- state
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
本发明公开了一种刚体空间运动状态的傅里埃输出方法,该方法通过定义三元数,使得机体轴系三个速度分量和三元数构成线性微分方程组,并采用傅里埃级数对滚转、俯仰、偏航角速度p,q,r进行近似逼近描述,可以按照任意阶保持器的方式求解系统的状态转移矩阵,进而得到刚体运动离散状态方程的表达式,避免了姿态方程奇异问题,从而得到刚体主要运动状态;本发明通过引入三元数使得状态转移矩阵为分块上三角形式,可以降阶求解状态转移矩阵,大大简化了计算复杂度,便于工程使用。
Description
技术领域
本发明涉及空间运动刚体模型,特别涉及飞行器大机动飞行状态输出问题。
背景技术
机体轴系刚体运动微分方程是描述飞行器、鱼雷、航天器等空间运动的基本方程。通常,在数据处理等应用中,体轴系的状态变量主要包含3个速度分量、三个欧拉角、以及地面坐标系的XE,YE,ZE等,由于ZE定义为垂直地面指向地球中心,因此ZE实际为负的飞行高度;XE,YE通常主要依赖GPS、GNSS、北斗等直接给出;欧拉角表示刚体空间运动姿态,而刻画刚体姿态的微分方程又是其中的核心,通常以三个欧拉角即俯仰、滚转和偏航角来描述。当刚体的俯仰角为±90°时,滚转角和偏航角无法定值,同时临近该奇点的区域求解误差过大,导致工程上不可容忍的误差而不能使用;为了避免这一问题,人们首先采用限制俯仰角取值范围的方法,这使得方程式退化,不能全姿态工作,因而难以广泛用于工程实践。随着对飞行器极限飞行的研究,人们又相继采用了方向余弦法、等效转动矢量法、四元数法等推算刚体运动姿态。
方向余弦法避免了欧拉角描述方法的“奇异”现象,用方向余弦法计算姿态矩阵没有方程退化问题,可以全姿态工作,但需要求解9个微分方程,计算量较大,实时性较差,无法满足工程实践要求。等效转动矢量法如单子样递推、双子样转动矢量、三子样转动矢量和四子样旋转矢量法以及在此基础上的各种修正算法和递推算法等。文献中研究旋转矢量时,都是基于速率陀螺输出为角增量的算法。然而在实际工程中,一些陀螺的输出是角速率信号,如光纤陀螺、动力调谐陀螺等。当速率陀螺输出为角速率信号时,旋转矢量法的算法误差明显增大。四元数法是定义4个欧拉角的函数来计算航姿,能够有效弥补欧拉角描述方法的奇异性,只要解4个一阶微分方程式组即可,比方向余弦姿态矩阵微分方程式计算量有明显的减少,能满足工程实践中对实时性的要求。其常用的计算方法有毕卡逼近法、二阶、四阶龙格-库塔法和三阶泰勒展开法等。毕卡逼近法实质是单子样算法,对有限转动引起的不可交换误差没有补偿,在高动态情况下姿态解算中的算法漂移会十分严重。采用四阶龙格-库塔法求解四元数微分方程时,随着积分误差的不断积累,会出现三角函数取值超出±1的现象,从而导致计算发散;泰勒展开法也因计算精度的不足而受到制约。当刚体大机动时,角速率较大导致上述方法的误差更大;不仅如此,姿态估计的误差常常会导致速度4个分量、高度输出的误差急剧增大。
发明内容
为了克服现有刚体运动模型输出误差大的问题,本发明提供一种刚体空间运动状态的傅里埃输出方法,该方法通过定义三元数,使得机体轴系三个速度分量和三元数构成线性微分方程组,并采用傅里埃级数对滚转、俯仰、偏航角速度p,q,r进行近似逼近描述,可以按照任意阶保持器的方式求解系统的状态转移矩阵,进而得到刚体运动离散状态方程的表达式,避免了姿态方程奇异问题,从而得到刚体主要运动状态。
本发明解决其技术问题采用的技术方案是,一种刚体空间运动状态的傅里埃输出方法,其特征包括以下步骤:
1、机体轴系三个速度分量输出为:
其中:u,v,w分别为沿刚体机体轴系x,y,z轴的速度分量,nx,ny,nz分别为沿x,y,z轴的过载,g为重力加速度,s1、s2、s3为定义的三元数,且
p,q,r分别为滚转、俯仰、偏航角速度,T为采样周期;全文参数定义相同;
ξA1(t)=[t sin(ωt)…sin[(n-1)ωt]sin(nωt)]T
ξB1(t)=[cos(ωt)cos(2ωt)…cos[(n-1)ωt]cos(nωt)]T
ω为角频率,滚转、俯仰、偏航角速度p,q,r的n阶展开式分别为
p(t)=[pa0 pa1…pa(n-1)pan][1cos(ωt)…cos[(n-1)ωt]cos(nωt)]T
+[pb1 pb2…pb(n-1) pbn][sin(ωt)sin(2ωt)…sin[(n-1)ωt]sin(nωt)]T
q(t)=[qa0 qa1…qa(n-1) qan][1cos(ωt)…cos[(n-1)ωt]cos(nωt)]T
+[qb1 qb2…qb(n-1) qbn][sin(ωt)sin(2ωt)…sin[(n-1)ωt]sin(nωt)]T
r(t)=[ra0 ra1…ra(n-1) ran][cos(ωt)cos(2ωt)…cos[(n-1)ωt]cos(nωt)]T
+[rb1 rb2…rb(n-1) rbn][sin(ωt)sin(2ωt)…sin[(n-1)ωt]sin(nωt)]T
2、高度输出为:
其中:h为高度;
3、姿态角的输出为:
其中:θ,ψ分别表示滚转、俯仰、偏航角,
本发明的有益效果是:通过引入三元数使得状态转移矩阵为分块上三角形式,可以降阶求解状态转移矩阵,大大简化了计算复杂度,便于工程使用。
下面结合实施例对本发明作详细说明。
具体实施方式
1、机体轴系三个速度分量输出为:
其中:u,v,w分别为沿刚体机体轴系x,y,z轴的速度分量,nx,ny,nz分别为沿x,y,z轴的过载,g为重力加速度,s1、s2、s3为定义的三元数,且
p,q,r分别为滚转、俯仰、偏航角速度,T为采样周期;全文参数定义相同;
ξA1(t)=[t sin(ωt)…sin[(n-1)ωt]sin(nωt)]T
ξB1(t)=[cos(ωt)cos(2ωt)…cos[(n-1)ωt]cos(nωt)]T
ω为角频率,滚转、俯仰、偏航角速度p,q,r的n阶展开式分别为
p(t)=[pa0 pa1 … pa(n-1) pan][1cos(ωt)…cos[(n-1)ωt]cos(nωt)]T
+[pb1 pb2 … pb(n-1) pbn][sin(ωt)sin(2ωt)…sin[(n-1)ωt]sin(nωt)]T
q(t)=[qa0 qa1 … qa(n-1) qan][1cos(ωt)…cos[(n-1)ωt]cos(nωt)]T
+[qb1 qb2 … qb(n-1) qbn][sin(ωt)sin(2ωt)…sin[(n-1)ωt]sin(nωt)]T
r(t)=[ra0 ra1 … ra(n-1) ran][cos(ωt)cos(2ωt)…cos[(n-1)ωt]cos(nωt)]T
+[rb1 rb2 … rb(n-1) rbn][sin(ωt)sin(2ωt)…sin[(n-1)ωt]sin(nωt)]T
2、高度输出为:
其中:h为高度;
3、姿态角的输出为:。
其中:θ,ψ分别表示滚转、俯仰、偏航角,
Claims (1)
1.一种刚体空间运动状态的傅里埃输出方法,其特征包括以下步骤:
机体轴系三个速度分量输出为:
其中:u,v,w分别为沿刚体机体轴系x,y,z轴的速度分量,nx,ny,nz分别为沿x,y,z轴的过载,g为重力加速度,s1、s2、s3为定义的三元数,且
p,q,r分别为滚转、俯仰、偏航角速度,T为采样周期;全文参数定义相同;
ξA1(t)=[t sin(ωt)…sin[(n-1)ωt]sin(nωt)]T
ξB1(t)=[cos(ωt)cos(2ωt)…cos[(n-1)ωt]cos(nωt)]T
ω为角频率,滚转、俯仰、偏航角速度p,q,r的n阶展开式分别为
p(t)=[pa0 pa1…pa(n-1) pan][1 cos(ωt)…cos[(n-1)ωt]cos(nωt)]T
+[pb1 pb2…pb(n-1) pbn][sin(ωt)sin(2ωt)…sin[(n-1)ωt]sin(nωt)]T
q(t)=[qa0 qa1…qa(n-1) qan][1 cos(ωt)…cos[(n-1)ωt]cos(nωt)]T
+[qb1 qb2…qb(n-1) qbn][sin(ωt)sin(2ωt)…sin[(n-1)ωt]sin(nωt)]T
r(t)=[ra0 ra1…ra(n-1) ran][cos(ωt)cos(2ωt)…cos[(n-1)ωt]cos(nωt)]T
+[rb1 rb2…rb(n-1) rbn][sin(ωt)sin(2ωt)…sin[(n-1)ωt]sin(nωt)]T
高度输出为:
其中:h为高度;
姿态角的输出为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110280706 CN102359790B (zh) | 2011-09-20 | 2011-09-20 | 一种刚体空间运动状态的傅里埃输出方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110280706 CN102359790B (zh) | 2011-09-20 | 2011-09-20 | 一种刚体空间运动状态的傅里埃输出方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102359790A CN102359790A (zh) | 2012-02-22 |
CN102359790B true CN102359790B (zh) | 2013-06-05 |
Family
ID=45585154
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110280706 Active CN102359790B (zh) | 2011-09-20 | 2011-09-20 | 一种刚体空间运动状态的傅里埃输出方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102359790B (zh) |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04193699A (ja) * | 1990-11-27 | 1992-07-13 | Mitsubishi Electric Corp | 三次元位置姿勢計測用マーク |
US6377906B1 (en) * | 2000-02-03 | 2002-04-23 | Independence Technology, L.L.C. | Attitude estimation in tiltable body using modified quaternion data representation |
ES2238936B1 (es) * | 2004-02-27 | 2006-11-16 | INSTITUTO NACIONAL DE TECNICA AEROESPACIAL "ESTEBAN TERRADAS" | Sistema y metodo de fusion de sensores para estimar posicion, velocidad y orientacion de un vehiculo, especialmente una aeronave. |
RU2282826C1 (ru) * | 2005-03-28 | 2006-08-27 | Федор Антонович Шаповалов | Способ автономного измерения вектора угловой скорости |
JP4743903B2 (ja) * | 2008-01-25 | 2011-08-10 | Necエンジニアリング株式会社 | スタートラッカ用ノイズ除去/グルーピング回路 |
DE102009059694B3 (de) * | 2009-12-18 | 2011-03-17 | Astrium Gmbh | Verfahren zur Lagebestimmung eines Raumfahrzeugs, insbesondere eines Satelliten, Computerprogramm zur Durchführung des Verfahrens, Datenträger und Raumfahrzeug |
-
2011
- 2011-09-20 CN CN 201110280706 patent/CN102359790B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN102359790A (zh) | 2012-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102436437B (zh) | 基于角速度的飞行器极限飞行时四元数傅里埃近似输出方法 | |
CN102589553B (zh) | 建立飞行器运动模型的切换方法 | |
CN113108781A (zh) | 一种应用于无人船行进间的改进粗对准算法 | |
CN102445202B (zh) | 一种刚体空间运动状态的拉盖尔输出方法 | |
CN102323990B (zh) | 一种刚体空间运动气动模型的建模方法 | |
CN102508819B (zh) | 基于角速度的飞行器极限飞行时四元数勒让德近似输出方法 | |
CN102359790B (zh) | 一种刚体空间运动状态的傅里埃输出方法 | |
CN102495830B (zh) | 基于角速度的飞行器极限飞行时四元数Hartley近似输出方法 | |
CN102495829B (zh) | 基于角速度的飞行器极限飞行时四元数沃尔什近似输出方法 | |
CN102384747A (zh) | 一种刚体空间运动状态的Hartley输出方法 | |
CN102346729B (zh) | 一种刚体空间运动状态的勒让德输出方法 | |
CN102495831B (zh) | 基于角速度的飞行器极限飞行时四元数埃米特近似输出方法 | |
CN102359789B (zh) | 一种刚体空间运动状态的任意阶输出方法 | |
CN102508818B (zh) | 一种刚体空间运动状态的任意步长正交级数输出模型建模方法 | |
CN102323992B (zh) | 一种刚体空间运动状态多项式类输出模型的建模方法 | |
CN102445203B (zh) | 一种刚体空间运动状态的埃米特输出方法 | |
CN102384746B (zh) | 一种刚体空间运动状态的切比雪夫输出的建模方法 | |
CN102506866B (zh) | 基于角速度的飞行器极限飞行时四元数切比雪夫近似输出方法 | |
CN102323991B (zh) | 一种刚体空间运动状态沃尔什输出模型的建模方法 | |
CN102495825B (zh) | 基于角速度的飞行器极限飞行时四元数超线性输出方法 | |
CN102506864B (zh) | 飞行器极限飞行时四元数任意步长正交级数近似输出方法 | |
CN102375803B (zh) | 一种刚体空间运动的气流轴系模型的建立方法 | |
CN102508821B (zh) | 一种刚体空间运动的状态输出模型建模方法 | |
CN102506865B (zh) | 基于角速度的飞行器极限飞行时四元数多项式类近似输出方法 | |
CN102494688B (zh) | 基于角速度的飞行器极限飞行时四元数拉盖尔近似输出方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |