基于小波方差的光纤陀螺随机误差特性提取方法技术

本发明专利技术公开了一种基于小波方差的光纤陀螺随机误差特性提取方法，该方法通过对光纤陀螺输出的角速度信息ｘｔ－０首先进行野值剔除处理，从而得到不含有野值的角速度信息ｘｔ，然后采用离散小波变换对ｘｔ进行处理得到不同分解尺度下的离散小波变换系数ｄｊ，τ，然后对ｄｊ，τ进行方差处理，从而获得各尺度因子τ下的从低频带至高频带的小波方差σｘ２（τ），最后采用随机误差源方差平方和关系进行小波方差σｘ２（τ）与尺度因子τ的双对数曲线拟合，从而提取出光纤陀螺输出的角速度信息中的五项随机误差系数Ｎ，Ｂ，Ｋ，Ｒ，Ｑ。本发明专利技术能够实现更加精确的光纤陀螺随机误差特性分析，从而为提高光纤陀螺的性能设计和补偿随机误差提供指导。

【技术实现步骤摘要】

本专利技术涉及一种对光纤陀螺输出的角速度信息xt-0进行性能参数的提取方法，更特别地说，是指一种采用小波方差法来对角速度信息xt-0进行处理后获得光纤陀螺的五项随机误差系数(N，B，K，R，Q)，根据这五项随机误差系数能够评价光纤陀螺在工作环境下的性能。
技术介绍
光纤陀螺(FOG)是新一代很有发展前途的惯性元件，随机误差是影响光纤陀螺精度的一个关键因素。光纤陀螺随机误差主要包括量化噪声、角度随机游走、偏置不稳定性、速率随机游走和速率斜坡等五个噪声项。 传统的分析光纤陀螺随机误差频率稳定性的方法是Allan方差法，它是基于数据分析的时域分析技术，也是IEEE标准化组织推荐的光学陀螺随机性能分析方法，这个方法的主要特点是它能非常容易地对光纤陀螺随机误差N，B，K，R，Q(N表示光纤陀螺输出的角速度信息xt中的角度随机游走系数、B表示光纤陀螺输出的角速度信息xt中的偏置不稳定性系数、K表示光纤陀螺输出的角速度信息xt中的速率随机游走系数、R表示光纤陀螺输出的角速度信息xt中的速率斜坡系数、Q表示光纤陀螺输出的角速度信息xt中的量化噪声系数)进行细致的表征和辨识，是测量和评价角速度信息xt中各类误差和噪声特性的一种重要手段，利用这种方法能够确定产生数据噪声的基本随机过程的特性，有助于识别数据中给定噪声项的来源。然而，Allan方差一般用来分析误差信号的稳定性，它对理想的时变、平稳信号是有效的，由于Allan方差的转移函数有基本的泄露现象，类似于无窗的离散傅立叶变化，因此Allan方差不能精确的表征非平稳信号在不同频率范围的方差波动情况。 从信号分析的角度来看，小波分析是时间频率二维平面分析的一种，是短时傅里叶变换(STFT)和Wingner-Ville分布(WD)的一个发展。小波分析实际上是一组带通滤波器，它使不同尺度、不同频率的信号通过不同的频带通道分离出来。小波分析在高频部分具有较高的时间分辨率和较低的频率分辨率，在低频部分具有较低的时间分辨率和较高的频率分辨率，因此小波分析被誉为“数学显微镜”。 实际的光纤陀螺输出信号，尤其是特定环境条件(如空间辐照)下的误差信号，即使在很短的时间内信号也表现出不稳定性。光纤陀螺的各项噪声分布于不同的频段上，如量化噪声对应高频段，而随机游走一般位于小于10Hz的频段，不稳定性噪声在更低的频段，各噪声源的频率特性各不相同，因此本专利技术用小波分析的方法把信号在不同的尺度分解，提取不同频率分量，并求取方差来表征噪声在不同频率分量的变化情况。
技术实现思路
本专利技术的目的是提出一种，该方法通过对光纤陀螺输出的角速度信息xt-0首先进行野值剔除处理，从而得到不含有野值的角速度信息xt，然后采用离散小波变换(DWT)对xt进行处理得到不同分解尺度下的离散小波变换系数dj，τ，然后对各尺度下的离散小波变换系数dj，τ进行方差处理，从而获得各尺度因子τ下的从低频带至高频带的小波方差σx2(τ)，最后采用随机误差源方差平方和关系 进行小波方差σx2(τ)与尺度因子τ的双对数曲线拟合，从而提取出光纤陀螺输出的角速度信息中的五项随机误差系数N，B，K，R，Q。本专利技术能够实现更加精确的光纤陀螺随机误差特性分析，从而为提高光纤陀螺的性能设计和补偿随机误差提供指导。 本专利技术的具有以下优点 (1)采用小波方差取代传统的Allan方差，实现了光纤陀螺随机误差项系数的提取和随机误差特性的分析。 (2)本专利技术采用的小波方差克服了Allan方差法存在的能量泄露的缺点，从而更加精确地进行噪声项的提取，并准确地反映光纤陀螺所处工作环境下的各项噪声的变化情况，为提高光纤陀螺的性能设计和补偿随机误差提供指导。 (3)在同等数据量的实施验证条件下，小波方差的计算效率比Allan方差的计算效率高，其算法运行速度明显增加、稳定性好。 附图说明 图1是本专利技术的流程图。 图2是整个辐照实验过程中的光纤陀螺输出信号图。 图3A是在辐照剂量为13.5Krad本专利技术所采用的算法分析得到的小波方差与相关时间的双对数曲线图。 图3B是在辐照剂量为31.5Krad本专利技术所采用的算法分析得到的小波方差与相关时间的双对数曲线图。 图3C是在辐照剂量为49.95Krad本专利技术所采用的算法分析得到的小波方差与相关时间的双对数曲线图。 图4A采用Allan方差法提取辐照环境下光纤陀螺量化噪声系数Q。 图4B采用Allan方差法提取辐照环境下光纤陀螺角度随机游走系数N。 图4C采用Allan方差法提取辐照环境下光纤陀螺偏置不稳定性系数B。 图4D采用Allan方差法提取辐照环境下光纤陀螺速率随机游走系数K。 图4E采用Allan方差法提取辐照环境下光纤陀螺速率斜坡系数R。 图5A采用本专利技术的基于小波方差法提取辐照环境下光纤陀螺量化噪声系数Q。 图5B采用本专利技术的基于小波方差法提取辐照环境下光纤陀螺角度随机游走系数N。 图5C采用本专利技术的基于小波方差法提取辐照环境下光纤陀螺偏置不稳定性系数B。 图5D采用本专利技术的基于小波方差法提取辐照环境下光纤陀螺速率随机游走系数K。 图5E采用本专利技术的基于小波方差法提取辐照环境下光纤陀螺速率斜坡系数R。 具体实施例方式 下面将结合附图和实施例对本专利技术做进一步的详细说明。 参见图1所示，本专利技术是一种，该提取方法具有下列步骤 步骤1剔除野值处理 采用莱特法则对光纤陀螺输出的角速度信息xt-0进行剔除野值，获得不含野值的有效的角速度信息xt，简称为野值剔除信息xt。 在本专利技术中，是以每组光纤陀螺输出的角速度信息xt-0的剩余误差的绝对值超过3倍的标准差作为条件，如果条件成立则该数据为野值加以剔除，从而得到完整的有效的光纤陀螺输出角速度信息xt。 在本专利技术中，莱特法则的数学表达形式可以为 以x1，x2，……，xn为一组在采样时间t内光纤陀螺输出的角速度信息，n表示采样个数。 计算在采样时间t内角速度信息的均值 计算在采样时间t内角速度信息的标准差 计算在采样时间t内角速度信息的野值判别|xi-x|≥3S。 步骤2离散小波变换处理 将野值剔除信息xt进行离散小波变换DWT获得各尺度因子τ下的从低频带至高频带的离散小波变换系数dj，τ(j表示时间偏移量，τ表示尺度因子)，dj，τ表示一组从低频带至高频带的离散小波变换系数，简称为单组小波系数。 在本专利技术中，尺度因子τ是指在采样时间t内采样个数n的偶数，数学表达为 离散小波变换(DWT)是将野值剔除信息xt进行基本小波ψ(t)的尺度伸缩和平移作处理，基本小波ψ(t)的尺度伸缩和平移为 ψj，τ(t)表示基本小波函数的尺度伸缩和平移量，j表示时间偏移量，τ表示尺度因子，t表示采样时间。 在本专利技术中，离散小波变换(DWT)系数 dj，τ表示离散小波变换系数，j表示时间偏移量，τ表示尺度因子，t表示采样时间，xt表示野值剔除信息，ψj，τ(t)表示基本小波函数的尺度伸缩和平移量。 步骤3求取方差处理 对单组小波系数dj，τ进行方差处理，获得各尺度因子τ下的从低频带至高频带的小波方差 σx2(τ)表示小波方差，τ表示尺度因子，dj，τ表示单组小波系数。该方差本文档来自技高网...

【技术保护点】
一种基于小波方差的光纤陀螺随机误差特性提取方法，其特征在于提取方法具有下列步骤：步骤１：剔除野值处理采用莱特法则对光纤陀螺输出的角速度信息ｘ↓［ｔ－０］进行剔除野值，从而获得野值剔除信息ｘ↓［ｔ］；步骤２：离散小波变换处理将野值剔除信息ｘ↓［ｔ］进行离散小波变换ＤＷＴ获得各尺度因子τ下的单组小波系数ｄ↓［ｊ，τ］＝＊ｘ↓［ｉ］ψ↓［ｊ，τ］（ｔ），ｊ表示时间偏移量，τ表示尺度因子，ｔ表示采样时间，ｘ↓［ｔ］表示野值剔除信息，ψ↓［ｊ，τ］（ｔ）表示基本小波函数的尺度伸缩和尺度平移量；所述的尺度因子τ是指在采样时间ｔ内采样个数ｎ的偶数，数学表达为τ＝１，２，４，．．．，ｎ／２；步骤３：求取方差处理对单组小波系数ｄ↓［ｊ，τ］进行方差处理，获得各尺度因子τ下的从低频带至高频带的小波方差σ↓［ｘ］↑［２］（τ）≡ｖａｒ｛ｄ↓［ｊ，τ］｝／τ，σ↓［ｘ］↑［２］（τ）表示小波方差，τ表示尺度因子，ｄ↓［ｊ，τ］表示单组小波系数；该方差σ↓［ｘ］↑［２］（τ）能够表征不同频率成分的信号波动性；根据离散小波变换ＤＷＴ系数ｄ↓［ｊ，τ］与尺度因子τ的关系，获得小波方差的估计值＊↓［ｘ］↑［２］（τ）＝１／τ１／ｎ／２τ＊ｄ↓［ｊ，τ］↑［２］，＊↓［ｘ］↑［２］（τ）表示在尺度因子τ下的小波方差估计值，τ表示尺度因子，ｎ表示光纤陀螺输出信号在采样时间ｔ内的采样个数，ｊ表示时间偏移量，ｄ↓［ｊ，τ］↑［２］表示单组小波系数的平方；步骤４：五项随机误差提取采用随机误差源方差平方和关系σ↓［ｘ］↑［２］（τ）＝Ｒ↑［２］τ↑［２］／２＋Ｋ↑［２］τ／３＋Ｂ↑［２］［２／π］ｌｎ２＋Ｎ↑［２］／τ＋３Ｑ↑［２］／τ↑［２］进行小波方差σ↓［ｘ］↑［２］（τ）与尺度因子τ的双对数曲线拟合，小波方差σ↓［ｘ］↑［２］（τ）为纵坐标的数值，尺度因子τ为横坐标的数值；从而提取出光纤陀螺输出的角速度信息中的五项随机误差系数Ｎ，Ｂ，Ｋ，Ｒ，Ｑ；Ｒ↑［２］τ↑［２］／２表示速率斜坡系数Ｒ噪声源的方差；Ｋ↑［２］τ／３表示速率随机游走系数Ｋ噪声源的方差；Ｂ↑［２］［２／π］ｌｎ２表示偏置不稳定性系数Ｂ噪声源的方差；Ｎ↑［２］／τ表示角度随机游走系数Ｎ噪声源的方差；３Ｑ↑［２］／τ↑［２］表示量化噪声系数Ｑ噪声源的方差。...

【技术特征摘要】
1.一种基于小波方差的光纤陀螺随机误差特性提取方法，其特征在于提取方法具有下列步骤步骤1剔除野值处理采用莱特法则对光纤陀螺输出的角速度信息xt-0进行剔除野值，从而获得野值剔除信息xt；步骤2离散小波变换处理将野值剔除信息xt进行离散小波变换DWT获得各尺度因子τ下的单组小波系数j表示时间偏移量，τ表示尺度因子，t表示采样时间，xt表示野值剔除信息，ψj，τ(t)表示基本小波函数的尺度伸缩和尺度平移量；所述的尺度因子τ是指在采样时间t内采样个数n的偶数，数学表达为步骤3求取方差处理对单组小波系数dj，τ进行方差处理，获得各尺度因子τ下的从低频带至高频带的小波方差σx2(τ)表示小波方差，τ表示尺度因子，dj，τ表示单组小波系数；该方差σx2(τ)能够表征不同频率成分的信号波动性；根据离散小波变换DWT系数dj，τ与尺度因子τ的关系，获得小波方差的估计值表示在尺度因子τ下的小波方差估计值，τ表示尺度因子，n表示光纤陀螺输出信号在采样时间t内的采样个数，j表示时间偏移量，dj，τ2表示单组小波系数的平方；步骤4五项随机误差提取采用随机误差源方差平方和关系进行小波方差σx2(τ)与尺度因子τ的双对数曲...

【专利技术属性】
技术研发人员：宋凝芳，陈婧，金靖，李敏，张晨，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：11[中国|北京]