基于改进GM(1,1)模型的激光陀螺仪随机误差预测 下载: 1200次
1 引言
激光陀螺仪作为高精度的惯性器件在激光捷联惯导领域应用广泛[1]。很多学者已对激光陀螺仪误差模型进行深入研究,通过对误差参数进行辨识、标定及补偿,减小激光陀螺仪的误差,提高激光捷联惯导的导航精度。激光陀螺仪的输出信号包含很多确定误差项和不确定的随机误差项[2],随机误差包括量化噪声、角度随机游走、零偏不稳定性、角速率随机游走、速率斜坡等,其中量化噪声、角度随机游走、零偏不稳定性是激光陀螺仪主要的性能指标[3]。对上述的误差进行分析的方法有常规的样本均值、自相关函数、功率率函数、Allan方差辨识等方法。样本均值法和方差法都无法知道误差的来源,自相关函数和功率谱函数法虽然是从时域和频域进行统计特性分析的,但是很难将随机误差从数据中分离出来[4]。静态下的Allan方差可以有效地辨识随机误差,但是在实际应用中的动态环境下,随机误差尤其是量测噪声,会随着工况、环境等因素的不同发生变化,因此研究动态Allan方差分析方法显得尤为重要。
Allan方差分析法是20世纪60年代由国家标准局提出来的一种时域分析方法[5],它可以用来分析陀螺仪等惯性器件或者其他精密仪器的随机误差,是一种公认的陀螺仪参数分析方法,Allan方差分析法可以对各种误差源和噪声统计特性进行辨识[6]。动态Allan方差可以看作是Allan方差的扩展,最初由Galleani等[7-8]提出,用于测量原子钟的时变特性,同时由美国海军用于时钟监控[8]。
灰色系统理论由我国学者邓聚龙1982年建立[9],其主要思想是信息可以分为白色、黑色、灰色,分别对应数据全部已知、数据未知、数据已知一部分的状态[10],当数据量比较小时将其视为小样本,并将小样本考虑为大样本的已知部分,最后将整个数据样本作为一个灰色系统进行预测,根据样本发展趋势拓展样本数量,得到整个样本未知的统计特性。灰色预测理论是一种基于灰色理论的有效数据预测的方法,其中的灰色GM(1,1)预测模型可以应用于实际生产和科学研究的诸多领域[11]。灰色GM(1,1)预测模型可以在已知一部分先验随机误差参数的情况下,对后续的随机误差进行灰色预测和残差检验,通过对随机误差的动态预测来补偿激光陀螺仪随机误差[12],最终达到提高惯性器件的数据精度的目的。
小波变换是信号处理和分析领域的重要方法[13],其中通过小波函数进行信号降噪和信号检测是研究的热点[14-15]。小波变换可以将信号分解为近似系数和细节系数,基于小波变换对信号进行小波分解可以得到信号的高频部分和低频部分,近似系数经过多层小波分解隔离掉细节抖动后,可以从近似系数曲线中检测得到信号是否存在间断点和坏值。小波分解后的细节系数主要是表征噪声干扰抖动,抑制信号分离的细节系数可以很好地降低信号受到的噪声干扰。
本文基于改进的灰色GM(1,1)预测模型,对动态工况下的激光陀螺仪随机误差进行预测,通过预测激光陀螺仪的随机误差,为激光陀螺仪的误差补偿提供依据。首先,提出一种新的动态工况下的Allan方差分析方法,建立灰色GM(1,1)预测模型,通过先验序列预测后续同种工况下的随机误差;接着,进行实测数据的实验验证。由实测数据直接建立灰色GM(1,1)预测模型的预测精度较低,只能达到精度3级,预测效果不佳[16-17]。本研究进一步采用小波滤波结合残差模型对传统GM(1,1)预测模型进行改进,采用自适应小波对实测随机误差系数进行平滑处理,使用残差修正GM(1,1)模型提高模型的精度[18-19]。
实验分析结果表明,动态Allan方差分析可以很好地辨识出激光陀螺仪在等时间间隔内的随机误差参数,从时域动态Allan方差分析曲线上可以看出,随机误差参数存在着波动上升的趋势。从Allan方差需求数据量大、得出的随机误差系数量小的特点来看,动态Allan方差分析符合灰色理论应用条件,在灰色GM(1,1)预测模型的基础上,基于改进GM(1,1)预测模型,预测出后续的随机误差参数,其预测结果精度比传统GM(1,1)预测模型精度提高了2个精度等级。
2 动态Allan方差分析
2.1 动态Allan方差分析
提出一种新的动态工况下激光陀螺仪信号的Allan方差分析方法,新动态Allan方差分析方法是利用矩形时间框对整个先验数据进行采样,由于静态Allan方差分析需要大量数据才能辨识随机误差参数[20],假设每段矩形框选择的数据时间段内随机误差参数缓慢变化,通过静态Allan方差分析得出所有矩形框的随机误差,并将随机误差按时间点列成灰色预测序列后进行灰色GM(1,1)预测模型预测,选择时间节点t=t1作为时间矩形框中心点,矩形框时间跨度为T,矩形框点所包含的时间范围为t1∈[t1-T/2,t1+T/2],并且假设在T时间内随机误差变化缓慢即静态Allan方差分析辨识误差系数作为该段时间内的误差系数值,计算得到的随机误差在t1时刻值为x(t1),假设先验数据框数目为n,动态Allan方差辨识出随机误差序列X=
3 GM(1,1)预测模型
3.1 建立GM(1,1)预测模型
灰色量的处理过程是利用灰色理论建立预测模型,采用数据生成的方法规律性地预测动态随机误差,常用的单序列一阶线性灰色模型记作GM(1,1)。
GM(1,1)建模的原始数据为
式中:x(0)(ti)为ti时刻动态Allan方差辨识的随机误差参数。
利用一次累加生成1-AGO单增序列,即
式中:x(1)(t1)=x(0)(t1),x(1)(tk)=
对于非负的数据序列,累加生成操作可以弱化随机性,增加规律性,使生成序列呈指数增长规律;对于产生的阶跃数据序列,采取小波检测间断点的方式识别出间断点,并用自适应小波对波动干扰数据进行光滑;最后利用处理后的数据建立GM(1,1)预测模型,步骤如下。
1) 建立GM(1,1)预测模型
设Z(1)为X(1)的邻近均值生成序列,即
式中:z(1)(tk)=[x(1)(tk)+x(1)(tk-1)]/2。则GM(1,1)的灰色微分方程为
式中:x(0)(tk)为灰导数;a为发展系数;z(1)(tk)为背景值;b为灰作用量,k=2,3,…,n。写成矩阵形式为
式中:Y为激光陀螺仪随机误差数据向量,Y=
利用最小二乘法求解,可得
将P=
因此,灰色微分方程的时间响应序列为
式中:tk+1时刻前数据为先验建模序列;tk+1时刻以后的数据为预测序列。
2)预测tk+1时刻的原始数据
灰色模型实际上是生成数列的模型,模型预测的数据需要经过逆生成还原后才可使用。还原值为
3)模型精度检验
检验GM(1,1)模型的精度,采用后验差进行检验。
预测误差定义为
为求解,首先设置一个估计误差函数,为
原始数列的均值和标准差分别定义为
预测误差的标准差定义为
方差比定义为
3.2 灰色预测模型的改进
原始数据的光滑特性是影响模型精度的主要因素,本研究结合激光陀螺仪的随机误差参数特点对GM(1,1)模型进行改进,改进方法如下。
1) 预处理原始数据序列,使其满足序列级比判断,引入小波滤波函数检测原始数据中的间断点与坏值。
2) 对原始数据进行光滑处理,利用自适应小波函数对数据的波动干扰进行平滑处理,增加数据的光滑性,对自适应小波滤波平滑后的预测序列进行GM(1,1)预测模型建模,建模结果的平均相对误差明显减小,提高了预测精度。
3) GM(1,1)预测模型的预测初始x(0)值与后续预测值有间断点,造成了曲线拟合效果差、预测精度差的问题。采用初始序列动态预测的方法[21],放弃初始先验序列x(0),加入x(k+1)构成新序列(k个数据为建模序列,得到第k+1个数据,也就是模型中的第n+1个递推值),也就是
4) 建立残差修正模型,修正原点附近的数据。
已知原有的预测模型(8)式,对应数据为
令k=i,i+1,…,n,得到残差序列q(0)(tk)=
对q(0)建立GM(1,1)模型,时间响应为
将上述残差的GM(1,1)模型加入原模型可得
式中:
改进GM(1,1)预测模型算法运行具体步骤如
4 实验分析
实验采用实验车、带有三轴激光陀螺仪的激光捷联惯导系统(LSINS)装置进行动态跑车测试实验,实验装置如
传统灰色GM(1,1)预测模型直接利用原始序列进行递推预测,结果如
表 2. GM(1,1)模型精度分级
Table 2. Accuracy classification of GM (1,1) model
|
表 1. 传统GM(1,1)模型随机误差预测方差比
Table 1. Variance ratios of random error prediction by traditional GM (1,1) model
|
表 3. 改进GM(1,1)模型随机误差预测方差比
Table 3. Variance ratios of random error prediction by modified GM (1,1) model
|
图 3. 传统灰色预测算法量测噪声系数预测图。 (a) X轴量测噪声系数;(b) Y轴量测噪声系数;(c) Z轴量测噪声系数
Fig. 3. Noise coefficient prediction graphs measured by traditional grey prediction algorithm. (a) X-axis measuremental noise coefficient; (b) Y-axis measuremental noise coefficient; (c) Z-axis measuremental noise coefficient
原始序列经过小波滤波平滑处理,处理后的曲线如
将平滑处理后的序列用于改进GM(1,1)预测模型建模,以量测噪声系数为例,利用建模数据和检验数据作图,如
图 5. 改进算法量测噪声系数预测图。 (a) X轴量测噪声系数;(b) Y轴量测噪声系数; (c) Z轴量测噪声系数
Fig. 5. Noise coefficient prediction graphs measured improved coefficient. (a) X-axis measuremental noise coefficient; (b) Y-axis measuremental noise coefficient; (c) Z-axis measuremental noise coefficient
图 6. 激光陀螺仪随机误差预测曲线图。 (a)量测噪声系数;(b)角度随机游走系数; (c)零偏不稳定性系数
Fig. 6. Random error prediction curves of laser gyroscope. (a) Measuremental noise coefficient; (b) angle random walk coefficient; (c) coefficient of zero deviation instability
通过实测数据计算改进GM(1,1)模型随机误差预测方差比如
5 结论
基于改进GM(1,1)预测模型预测激光陀螺仪的随机误差参数,采用矩形时间框等间隔采集动态Allan方差辨识随机误差参数。基于小波变换检测GM(1,1)预测模型序列间断点,使用小波滤波平滑GM(1,1)预测模型序列,对模型进行残差修正,最后利用多组实测数据进行模型的建模和验证。实验结果表明,改进的灰色GM(1,1)预测模型可以预测激光陀螺仪的随机误差参数,预测精度方差比相比于传统灰色GM(1,1)预测模型提高了2个精度等级,达到了1级的精度等级,对于激光陀螺仪的随机误差的补偿提供了重要依据。
[1] 战德军, 秦石乔, 王省书, 等. 速率偏频激光陀螺标度因数正反转不对称性的研究[J]. 中国激光, 2009, 36(8): 1995-1998.
[2] 杜少林, 陈书钊, 陈鹏光, 等. 基于Allan方差的MEMS陀螺仪噪声分析[J]. 仪表技术与传感器, 2018( 5): 20- 22, 27.
Du SL, Chen SZ, Chen PG, et al. Analysis of MEMS gyroscope noise based on Allan variance[J]. Instrument Technique and Sensor, 2018( 5): 20- 22, 27.
[3] 马群, 王庆, 阳媛, 等. 基于Allan方差的MEMS陀螺仪随机误差辨识与抑制[J]. 传感器与微系统, 2019, 38(6): 62-65.
Ma Q, Wang Q, Yang Y, et al. Random error identification and suppression of MEMS gyroscope based on Allan variance[J]. Transducer and Microsystem Technologies, 2019, 38(6): 62-65.
[4] 魏国, 龙兴武. 基于动态Allan方差的机抖激光陀螺随机误差研究[J]. 中国激光, 2010, 37(12): 2975-2979.
[5] El-Sheimy N, Hou H Y, Niu X J. Analysis and modeling of inertial sensors using Allan variance[J]. IEEE Transactions on Instrumentation and Measurement, 2008, 57(1): 140-149.
[6] 张谦, 王玮, 王蕾, 等. 基于动态Allan方差的光纤陀螺随机误差分析及算法改进[J]. 光学学报, 2015, 35(4): 0406003.
[7] GalleaniL, TavellaP. Tracking nonstationarities in clock noises using the dynamic Allan variance[C]∥Proceedings of the 2005 IEEE International Frequency Control Symposium and Exposition, August 29-31, 2005, Vancouver, BC, Canada. New York: IEEE, 2005: 392- 396.
[8] GalleaniL, TavellaP. The characterization of clock behavior with the dynamic Allan variance[C]∥IEEE International Frequency Control Symposium and PDA Exhibition Jointly with the 17th European Frequency and Time Forum, May 4-8, 2003, Tampa, FL, USA. New York: IEEE, 2003: 8029712.
[9] 罗党, 韦保磊. 一类离散灰色预测模型的统一处理方法及应用[J]. 系统工程理论与实践, 2019, 39(2): 451-462.
Luo D, Wei B L. A unified treatment approach for aclass of discrete grey forecasting models and its application[J]. Systems Engineering: Theory & Practice, 2019, 39(2): 451-462.
[10] 刘思峰. 灰色系统理论的产生与发展[J]. 南京航空航天大学学报, 2004, 36(2): 267-272.
Liu S F. Emergence and development of grey system theory and its forward trends[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2004, 36(2): 267-272.
[11] 杨书霞, 胡艳. 基于动态灰色预测的交通事故损失预测[J]. 武汉理工大学学报(信息与管理工程版), 2017, 39(6): 674-678.
Yang S X, Hu Y. Prediction of traffic accident loss based on dynamic grey prediction[J]. Journal of Wuhan University of Technology (Information & Management Engineering), 2017, 39(6): 674-678.
[12] 唐天国, 万星, 刘浩吾. 高边坡安全监测的改进GM模型预测研究[J]. 岩石力学与工程学报, 2005, 24(2): 307-312.
Tang T G, Wan X, Liu H W. Improved GM model for safety monitoring of high rock slopes[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(2): 307-312.
[13] 李家垒, 许化龙, 何婧. 光纤陀螺信号的小波包去噪及改进[J]. 光学学报, 2010, 30(8): 2224-2228.
[14] 孟鑫, 刘磊, 江升, 等. 基于小波变换的干涉光谱信号检测与校正方法[J]. 光学学报, 2019, 39(9): 0930007.
[15] 王行芳, 金施群, 侯少阳. 基于小波阈值的激光探测声音信号去噪研究[J]. 大气与环境光学学报, 2018, 13(5): 388-394.
[16] 方文辉, 卢伟, 洪德林, 等. 基于多尺度小波变换和灰色神经网络的稻种发芽率红外热预测模型的研究[J]. 光学学报, 2015, 35(10): 1012001.
[17] 程宇琼, 卢伟, 罗慧, 等. 基于连续偏振光谱技术与嵌入型灰色神经网络的稻种发芽率检测方法研究[J]. 光学学报, 2015, 35(12): 1230001.
[18] 鞠荟荟, 刘志刚, 姜江军, 等. 基于低通滤波残差图的高光谱条带噪声去除[J]. 光学学报, 2018, 38(12): 1228002.
[19] 周朴, 王小林, 马阎星, 等. 基于随机并行梯度下降算法自适应光学系统的校正残差分析[J]. 光学学报, 2010, 30(3): 613-617.
[20] 张娜, 李绪友. 动态Allan方差的理论改进及其应用研究[J]. 光学学报, 2011, 31(11): 1106003.
[21] 王立冬, 车琳, 鲁军, 等. 基于改进型灰色预测模型的SINS/GPS组合导航系统[J]. 中国惯性技术学报, 2015, 23(2): 248-251, 257.
Wang L D, Che L, Lu J, et al. SINS/GPS integrated navigation system based on improved grey forecasting model[J]. Journal of Chinese Inertial Technology, 2015, 23(2): 248-251, 257.
Article Outline
李想, 汪立新, 沈强. 基于改进GM(1,1)模型的激光陀螺仪随机误差预测[J]. 光学学报, 2020, 40(12): 1204001. Xiang Li, Lixin Wang, Qiang Sheng. Prediction of the Random Error of a Laser Gyroscope Using the Modified GM (1, 1) Model[J]. Acta Optica Sinica, 2020, 40(12): 1204001.