摘要:在圆锥运动条件下,具有相同时间间隔的两个角增量的叉乘对圆锥补偿的贡献相等。本文根据这一特点,设计了一种捷联惯导系统姿态算法,在圆锥补偿获得相同精度的情况下,它的计算量较少。同时,利用此特点,推导出利用前一圆锥补偿周期的角增量进行圆锥补偿的算法,提高了补偿精度。本文给出了仿真实验结果。本文为研制激光陀螺捷联惯导系统提供了一种高精度算法。
关键词:捷联惯性导航系统; 姿态算法;圆锥补偿算法
1 引言
在捷联惯性导航系统中,陀螺仪与加速度计固联于载体上,直接测量载体角运动和线运动,因而要求较高的测量带宽。陀螺仪的输出是积分时间内载体的角增量的量化值。刚体的有限转动不是矢量,其转动次序不可交换。无限小转动是矢量[1]。Bortz的转动向量(orientation vector)微分方程为计算捷联惯性系统的姿态矩阵建立了全面的理论基础[2][3]。转动向量的变化速率是惯性测量的角速率向量(陀螺仪输出)与计算得到的非互易速率向量二者之和。后者是影响捷联惯性导航系统姿态角精度的一个重要因素。因此,在高角速率动态环境中,为防止姿态误差积累,必须对非互易速率向量进行补偿。对非互易速率向量进行补偿的计算一般称为圆锥补偿算法。要提高系统精度,可以有两种选择,一种是对陀螺输出信号进行高速采样,使用简单的圆锥补偿算法;另一种是适当降低对陀螺输出信号进行采样的速率,使用复杂但精确的圆锥补偿算法[4]。
设计补偿精度高、计算量少的圆锥补偿算法一直是人们不断追求得目标。Miller提出了在圆锥运动条件下优化的三子样算法[5],Jang G.Lee等提出了四子样算法[6],这两种算法没有利用前一个圆锥补偿周期的陀螺输出信号,而且计算量较大。Yeon Fuh Jiang等提出了利用前一补偿周期的一类算法[7],包含了圆锥补偿周期内陀螺输出信号的所有可能的叉乘项,并利用前一补偿周期的陀螺输出累加角增量。Musoff提出了圆锥算法的优化指标,分析了圆锥补偿后的算法误差与补偿周期的幂次r的关系[8]。本文利用在经典圆锥运动下,相同时间间隔的角增量向量的叉乘对圆锥误差的贡献相等的特点,推导出一类计算量较少的圆锥补偿算法。同时利用前一圆锥补偿周期的陀螺输出信号,提高了补偿精度。
2 圆锥补偿算法
载体转动向量的微分方程为
式中,Φ是转动向量,ω是机体角速率向量。方程右边第二项与第三项之和就是非互易误差。对于小Φ,为了简化计算,忽略右边第三项,并用下式代替第二项中的Φ
α=∫ωdt (2)
这样,式(1)可简化为
(5)
β就是需要进行补偿的非互易向量。随圆锥补偿周期内陀螺仪的采样次数不同,β可以有许多种方式进行计算。
对于经典圆锥运动,它的转动向量为[9]
(6)
式中,Ω是圆锥运动角频率,a是圆锥运动的幅度。采用上标B来表示机体坐标系,得到
(7)
假设在t~t+h间隔内,对陀螺输出信号进行N次采样,θi表示第i次采样的陀螺输出的增量角信号。
(8)
其中,
(9)
则在t~t+h间隔内,转动向量的估计值的计算公式为
(10)
根据式(4),圆锥补偿可计算如下
(11)
右边共有N.(N-1)/2项。由式(9),得
(12)
由上式可以看出,在θi×θj中,x分量只与相对时间(i-j)h有关,而与绝对时间无关,y和z分量是绝对时间t的余弦函数。在本文的圆锥运动下,能引起漂移的误差出现在x轴上。因此,具有相同时间间隔的两个角增量向量的叉乘对圆锥误差的贡献是相等的,而可以不考虑它们与绝对时间的关系。利用这个性质,可以将圆锥误差补偿公式简化为
(13)
比较(13)和(11)式可以看出,简化后的算法,计算量大大减小。简化带来的害处是在一般机动条件下,圆锥误差补偿算法的误差会增大。但是即使是在非常大的机动条件下,算法误差也是在可以忽略的水平上,而且这种误差只是在机动期间产生。因此,对圆锥补偿算法的这种简化是合理的。
如果考虑到可以利用前一圆锥补偿周期的角增量输出,则可以得到以下算式
(14)
ai、bi是加权系数,p是要利用前一圆锥补偿周期的角增量的个数,θ′N-i+1是前一圆锥补偿周期的陀螺输出信号采样,即角增量。式(14)与文献[7]中的式(9)相比,计算简化了,而且并不是利用前一周期陀螺输出信号采样的累加值。也可以参照文献[7]的式(9),得到利用前一圆锥补偿周期的累加值的圆锥补偿公式:
(15)
式中,θ′是前一圆锥补偿周期的累加值。对应采样为2、3、4、5等情况,可得到简化算法式(14)、(15)的各种形式。各种圆锥补偿算法的系数如表1所示。算法的性能可以用经典圆锥运动下在x轴上的算法漂移速率来表征[8]
表1 各种圆锥补偿算法的系数
| 算法编号 | G | a3 | a2 | a1 | b1 | b2 | b3 | b4 | r | |
| 二子样 | 1 | -1/30 | 11/15 | 6 | ||||||
| 2 | 1/140 | 13/210 | 323/420 | 8 | ||||||
| 3 | -1/180 | 32/45 | 6 | |||||||
| 三子样 | 4 | 9/20 | 27/20 | 6 | ||||||
| 5 | 3/280 | 57/140 | 393/280 | 8 | ||||||
| 6 | -1/420 | 1/40 | 157/420 | 1207/ 840 | 10 | |||||
| 7 | 1/1848 | -31/ 4620 | 61/1540 | 1607/ 4620 | 13487/ 9240 | 12 | ||||
| 8 | 1/3360 | 243/560 | 1539/ 1120 | 8 | ||||||
| 四子样 | 9 | 54/105 | 92/105 | 214/105 | 8 | |||||
| 10 | -1/315 | 168/315 | 262/315 | 656/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/252 | 25/24 | 325/252 | 1375/ 504 | 10 | ||||
| 15 | 5/5544 | 1355/ 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中。 3 仿真计算 表2 各算法的漂移误差(°/h)(Ω=40 Hz,a=1°) |
| 二采样 | 算法1 | 算法2 | 算法3 | 算法Ⅰ | 算法Ⅱ | ||
| h=0.02 s | 2.0 | 0.27 | 3.4 | 4.7e+2 | 15.4 | ||
| h=0.01 s | 3.4e-2 | 1.3e-3 | 6.0e-2 | 1.3e+2 | 1.03 | ||
| h=0.005 s | 5.6e-4 | -1.5e-6 | 1.0e-3 | 33.1 | 6.6e-2 | ||
| 三采样 | 算法4 | 算法5 | 算法6 | 算法7 | 算法8 | 算法Ⅲ | 算法Ⅳ |
| h=0.02 s | 2.0 | 0.29 | 6.3e-2 | 3.2e-2 | 0.56 | -7.8 | 2.0 |
| h=0.01 s | 3.5e-2 | 2.0e-3 | 8.3e-4 | 8.1e-4 | 3.3e-3 | -0.68 | 0.035 |
| h=0.005 s | 5.7e-4 | 1.1e-5 | 4.4e-6 | 4.4e-6 | 1.7e-5 | -4.9e-2 | 5.7e-4 |
| 四采样 | 算法9 | 算法10 | 算法11 | 算法12 | 算法13 | 算法V | |
| h=0.02 s | 0.37 | 0.14 | 0.11 | 0.10 | 0.18 | 0.37 | |
| h=0.01 s | 4.6e-3 | 3.5e-3 | 3.4e-3 | 3.4e-3 | 3.5e-3 | 4.1e-3 | |
| h=0.005 s | 1.0e-4 | 1.0e-4 | 1.1e-4 | 9.8e-5 | 9.7e-5 | 5.4e-5 | |
| 五采样 | 算法14 | 算法15 | 算法16 | 算法17 | 算法18 | ||
| h=0.02 s | 0.26 | 0.23 | 0.22 | 0.22 | 0.23 | ||
| h=0.01 s | 1.2e-2 | 1.2e-2 | 1.2e-2 | 1.2e-2 | 1.2e-2 | ||
| h=0.005 s | 3.8e-4 | 4.0e-4 | 4.0e-4 | 4.0e-4 | 3.8e-4 |
| 由表2看出,在所给仿真条件下,三采样算法7的漂移速率最小,精度最高。当圆锥补偿周期的采样数小于4时,遵循表1的r数越大,算法精度越高的原则。当圆锥补偿周期的采样数大于、等于4时,不再遵循表1的r数越大,算法精度越高的原则。 图1 仿真程序流程图 (a) (b) 4 结论 基金项目:国防科工委“九五”重点预研项目9.2.7 作者简介:刘巧光,男,清华大学博士生,从事惯性技术研究。 作者单位:清华大学,北京 100084 参考文献: |




