捷联惯性导航系统的姿态算法优化设计

   2024-02-18 互联网920
核心提示:摘要:在圆锥运动条件下,具有相同时间间隔的两个角增量的叉乘对圆锥补偿的贡献相等。本文根据这一特点,设计了一种捷联惯导系统

 

摘要:在圆锥运动条件下,具有相同时间间隔的两个角增量的叉乘对圆锥补偿的贡献相等。本文根据这一特点,设计了一种捷联惯导系统姿态算法,在圆锥补偿获得相同精度的情况下,它的计算量较少。同时,利用此特点,推导出利用前一圆锥补偿周期的角增量进行圆锥补偿的算法,提高了补偿精度。本文给出了仿真实验结果。本文为研制激光陀螺捷联惯导系统提供了一种高精度算法。
关键词:捷联惯性导航系统; 姿态算法;圆锥补偿算法

1 引言
  在捷联惯性导航系统中,陀螺仪与加速度计固联于载体上,直接测量载体角运动和线运动,因而要求较高的测量带宽。陀螺仪的输出是积分时间内载体的角增量的量化值。刚体的有限转动不是矢量,其转动次序不可交换。无限小转动是矢量[1]。Bortz的转动向量(orientation vector)微分方程为计算捷联惯性系统的姿态矩阵建立了全面的理论基础[2][3]。转动向量的变化速率是惯性测量的角速率向量(陀螺仪输出)与计算得到的非互易速率向量二者之和。后者是影响捷联惯性导航系统姿态角精度的一个重要因素。因此,在高角速率动态环境中,为防止姿态误差积累,必须对非互易速率向量进行补偿。对非互易速率向量进行补偿的计算一般称为圆锥补偿算法。要提高系统精度,可以有两种选择,一种是对陀螺输出信号进行高速采样,使用简单的圆锥补偿算法;另一种是适当降低对陀螺输出信号进行采样的速率,使用复杂但精确的圆锥补偿算法[4]
  设计补偿精度高、计算量少的圆锥补偿算法一直是人们不断追求得目标。Miller提出了在圆锥运动条件下优化的三子样算法[5],Jang G.Lee等提出了四子样算法[6],这两种算法没有利用前一个圆锥补偿周期的陀螺输出信号,而且计算量较大。Yeon Fuh Jiang等提出了利用前一补偿周期的一类算法[7],包含了圆锥补偿周期内陀螺输出信号的所有可能的叉乘项,并利用前一补偿周期的陀螺输出累加角增量。Musoff提出了圆锥算法的优化指标,分析了圆锥补偿后的算法误差与补偿周期的幂次r的关系[8]。本文利用在经典圆锥运动下,相同时间间隔的角增量向量的叉乘对圆锥误差的贡献相等的特点,推导出一类计算量较少的圆锥补偿算法。同时利用前一圆锥补偿周期的陀螺输出信号,提高了补偿精度。

2 圆锥补偿算法
  载体转动向量的微分方程为

1-1.gif (924 bytes)    (1)
1-2.gif (771 bytes)

  式中,Φ是转动向量,ω是机体角速率向量。方程右边第二项与第三项之和就是非互易误差。对于小Φ,为了简化计算,忽略右边第三项,并用下式代替第二项中的Φ

α=∫ωdt    (2)

  这样,式(1)可简化为

1-3.gif (435 bytes)    (3)
1-4.gif (812 bytes)  (4)

  式中,

1-5.gif (487 bytes)    (5)

  β就是需要进行补偿的非互易向量。随圆锥补偿周期内陀螺仪的采样次数不同,β可以有许多种方式进行计算。
  对于经典圆锥运动,它的转动向量为[9]

1-6.gif (730 bytes)    (6)

  式中,Ω是圆锥运动角频率,a是圆锥运动的幅度。采用上标B来表示机体坐标系,得到

1-7.gif (1369 bytes)    (7)

  假设在t~t+h间隔内,对陀螺输出信号进行N次采样,θi表示第i次采样的陀螺输出的增量角信号。

1-8.gif (348 bytes)    (8)

  其中,

1-9.gif (3217 bytes)    (9)

  则在t~t+h间隔内,转动向量的估计值的计算公式为

1-10.gif (778 bytes)    (10)

根据式(4),圆锥补偿可计算如下

1-11.gif (720 bytes)    (11)

  右边共有N.(N-1)/2项。由式(9),得

1-12.gif (5221 bytes)(12)

  由上式可以看出,在θi×θj中,x分量只与相对时间(i-j)h有关,而与绝对时间无关,y和z分量是绝对时间t的余弦函数。在本文的圆锥运动下,能引起漂移的误差出现在x轴上。因此,具有相同时间间隔的两个角增量向量的叉乘对圆锥误差的贡献是相等的,而可以不考虑它们与绝对时间的关系。利用这个性质,可以将圆锥误差补偿公式简化为

1-13.gif (588 bytes)    (13)

  比较(13)和(11)式可以看出,简化后的算法,计算量大大减小。简化带来的害处是在一般机动条件下,圆锥误差补偿算法的误差会增大。但是即使是在非常大的机动条件下,算法误差也是在可以忽略的水平上,而且这种误差只是在机动期间产生。因此,对圆锥补偿算法的这种简化是合理的。
  如果考虑到可以利用前一圆锥补偿周期的角增量输出,则可以得到以下算式

1-14.gif (907 bytes)    (14)

  ai、bi是加权系数,p是要利用前一圆锥补偿周期的角增量的个数,θN-i+1是前一圆锥补偿周期的陀螺输出信号采样,即角增量。式(14)与文献[7]中的式(9)相比,计算简化了,而且并不是利用前一周期陀螺输出信号采样的累加值。也可以参照文献[7]的式(9),得到利用前一圆锥补偿周期的累加值的圆锥补偿公式:

1-15.gif (713 bytes)    (15)

  式中,θ是前一圆锥补偿周期的累加值。对应采样为2、3、4、5等情况,可得到简化算法式(14)、(15)的各种形式。各种圆锥补偿算法的系数如表1所示。算法的性能可以用经典圆锥运动下在x轴上的算法漂移速率来表征[8]

表1 各种圆锥补偿算法的系数

 算法编号Ga3a2a1b1b2b3b4r
二子样1   -1/3011/15   6
2  1/14013/210323/420   8
3-1/180   32/45   6
三子样4    9/2027/20  6
5   3/28057/140393/280  8
6  -1/4201/40157/4201207/

840

  10
7 1/1848
-31/
4620
61/1540
1607/
4620
13487/
9240
  12
81/3360
   243/560
1539/
1120
  8
四子样9    54/10592/105214/105 8
10   -1/315168/315262/315656/315 10
11
  1/1386
-31/
3465
1 277/
2310
2762/
3465
7321/
3465
 12
12
 -1/
6006
43/
18018
-733/
45045
5717/
10010
69337/
90090
96163/
45045
 14
13
-1/
69300
   8992/
17325
14912/
17325
1696/
825
 10
五子样14    125/25225/24325/2521375/
504
10
15
   5/55441355/
2772
2955/
2772
3455/
2772
15335/
5544
12
16
  -5/
24024
215/
72072
17285/
36036
1090/
1001
87355/
72072
201335/
72072
14
17
 1/20592-19/
24024
295/
48048
16921/
36036
53 321/
48048
28451/
24024
405673/
144144
16
18
1/
1513512
   374375/
756756
1586875/
1513512
241250/
189189
518750/
189189
12

drift=khr

  式中,r是圆锥补偿周期h的幂次,它表征了圆锥补偿算法的精度。对于每一个算法来说,如果圆锥运动幅度a和圆锥运动角频率Ω是确定的,k为一个常数。各个算法的精度表征数r也列于表1中。
  从表1中的系数可以发现,算法3与文献[7]的算法2相同,而算法8与文献[7]的算法3一致,只要后者的系数选择K12=0。算法4与文献[5][7]的优化算法一致,只要后者的系数选择K12=0。由此也可以看出本文这种简化计算的合理性。由表1还表明,利用前一圆锥补偿的角增量输出,能够提高补偿精度。如三采样的算法5与四采样的算法9的r数是一样的,都是8。表1还表明,不管哪个算法,只要增加一个系数,其r数就增加2。因此应该充分利用前一圆锥补偿周期的角增量输出。

3 仿真计算
  仿真程序用Matlab语言编写。其流程图如图1所示。轨迹仿真选择为(6)式的经典圆锥运动。Ω=40 Hz,a=1°,取陀螺仪采样周期为h=0.02、
0.01、0.005 s。计算表1中所有算法及文献[5][6]中的算法的漂移误差。仿真计算结果如表2所示。算法I为不进行圆锥补偿的二采样算法,算法Ⅱ为Jordan预处理算法[5],算法Ⅲ为文献[5]中的三采样算法,算法Ⅳ为文献[5]中的三采样优化算法,算法V为文献[6]中的四采样算法。图2(a)、(b)分别为Ω=40 Hz,a=1°,h=0.01 s时,算法Ⅰ和算法Ⅱ的x轴漂移曲线。

表2 各算法的漂移误差(°/h)(Ω=40 Hz,a=1°)

二采样算法1算法2算法3算法Ⅰ算法Ⅱ  
h=0.02 s2.00.273.44.7e+215.4  
h=0.01 s3.4e-21.3e-36.0e-21.3e+21.03  
h=0.005 s5.6e-4-1.5e-61.0e-333.16.6e-2  
三采样算法4算法5算法6算法7算法8算法Ⅲ算法Ⅳ
h=0.02 s2.00.296.3e-23.2e-20.56-7.82.0
h=0.01 s3.5e-22.0e-38.3e-48.1e-43.3e-3-0.680.035
h=0.005 s5.7e-41.1e-54.4e-64.4e-61.7e-5-4.9e-25.7e-4
四采样算法9算法10算法11算法12算法13算法V 
h=0.02 s0.370.140.110.100.180.37 
h=0.01 s4.6e-33.5e-33.4e-33.4e-33.5e-34.1e-3 
h=0.005 s1.0e-41.0e-41.1e-49.8e-59.7e-55.4e-5 
五采样算法14算法15算法16算法17算法18  
h=0.02 s0.260.230.220.220.23  
h=0.01 s1.2e-21.2e-21.2e-21.2e-21.2e-2  
h=0.005 s3.8e-44.0e-44.0e-44.0e-43.8e-4  
  由表2看出,在所给仿真条件下,三采样算法7的漂移速率最小,精度最高。当圆锥补偿周期的采样数小于4时,遵循表1的r数越大,算法精度越高的原则。当圆锥补偿周期的采样数大于、等于4时,不再遵循表1的r数越大,算法精度越高的原则。

101.gif (3064 bytes)

图1 仿真程序流程图

102.gif (2386 bytes)

(a)

103.gif (5106 bytes)

(b)
图2 x轴姿态角漂移曲线(Ω=40 Hz,a=1°,h=0.01 s):(a) 算法Ⅰ,(b) 算法Ⅱ

4 结论
  在所选定的载体圆锥运动条件下的仿真结果表明:
  ① 这类算法比相同精度的其它算法计算量少。如算法4和算法Ⅳ计算精度相同,计算量却少三分之一。算法9与算法V计算精度基本相同,但计算量要少50%。
  ② 所有算法的精度都随着陀螺输出信号的采样频率提高而提高,说明在激光陀螺捷联惯导系统中,要努力提高激光陀螺仪的采样频率。
  ③ 本文提出的三采样算法7在所有仿真算法中,漂移为最小,表明在陀螺输出信号采样频率一定时,一个圆锥补偿周期内的采样数不是越多越好。
  ④ 圆锥补偿周期的采样数大于、等于4时,利用前一补偿周期的采样,对增加精度效果不明显。圆锥补偿周期的采样数小于4时,分别利用前一圆锥补偿周期陀螺输出信号的各次采样,而不是利用前一补偿周期各次采样输出的累加量,补偿精度将更高。
  ⑤ 圆锥补偿周期的采样数小于4时,在圆锥补偿算法中,利用前一圆锥补偿周期陀螺输出信号,相当于提高了陀螺输出信号的采样频率。

基金项目:国防科工委“九五”重点预研项目9.2.7

作者简介:刘巧光,男,清华大学博士生,从事惯性技术研究。

作者单位:清华大学,北京 100084

参考文献:
[1] 陈哲.捷联惯性导航系统原理.北京:宇航出版社,1986.
[2] Bortz J E. A new mathematical formulation for strapdown inertial navigation [J].IEEE Transacfions on AerosPACe and ElectroNIc Systems,1971,7(1):61~66.
[3] Ignagni M B. On the orientation vector differential equation in strapdown inertial systems [J].IEEE Transactions on Aerospace and Electronic Systerms,1994,30(4):1076~1081.
[4] Savage P G.Strapdown inertial navigation integration algorithm design,Part l:Attitude algorithm [J].Journal of Guidance,Control,and Dynamics,1998,21(1):19~28.
[5] Miller R B. A new strapdown attitude algorithm [J]. Journal of Guidance,Control,and Dynamics,1983,6(4):287~291.
[6] Lee J G, Yoon Y J,Mark J G.Extension of strapdown attitude algorithm for highfrequency base motion [J].Journal of Guidance,Control,and Dynamics,1990,13(4):738~743.
[7] Jiang Y F,Lin Y P.Improved strapdown coning algorithm [J].IEEE Transactions on Aerospace and Electronic Systems,1992,28(2):484~490.
[8] Musoff H,Murphy J H.Study of strapdown navigation attitude algorithms [J].Journal of Guidance,Control,and Dynamics,1995,18(2):287~290.
[9] Jiang Y F,Lin Y P.On the rotation vector differential equation [J].IEEE Transactions on Aerospace and Electronic Systems,1991,27(1):181~183.


 
举报收藏 0打赏 0评论 0
 
更多>同类资讯
推荐图文
推荐资讯
点击排行
网站首页  |  关于我们  |  联系方式  |  使用协议  |  版权隐私  |  隐私政策  |  网站地图  |  排名推广  |  广告服务  |  积分换礼  |  RSS订阅