技术领域
[0001] 本
发明涉及一种基于交叠M带
离散小波变换(Overlap M-band Discrete Wavelet Transform，OMDWT)的光纤陀螺信号去噪方法，属于惯性导航技术领域。
背景技术
[0002] 光纤陀螺是以Sagnac效应为
基础而发展起来的新型全固态
陀螺仪，受到世界各国的高度重视，已成为新一代中低
精度捷联惯性导航、制导系统中的理想惯性器件。捷联式惯导系统中的一个关键技术是如何完全自主地实现快速初始对准，即快速自主地确定载体的初始方位
角和
水平
姿态角。影响方位对准(寻北)精度的主要因素是陀螺漂移，有效地消除陀螺漂移是保证方位对准精度的关键。
[0003] 光纤陀螺漂移可分为系统性漂移和随机漂移两种类型。由于随机漂移是弱非线性的、慢时变的，同时受外部环境等不确定性因素的影响，无法建立其精确系统模型，在惯导系统中不能用简单的方法加以补偿，因此随机漂移成为衡量陀螺仪精度的重要指标。为了提高初始对准精度，必须采用有效的
信号处理方法对陀螺信号处理，来抑制这些随机漂移。目前对陀螺信号的处理有两种思想：(1)建立陀螺随机漂移模型，使用Kalman滤波、Wiener滤波等方法进行补偿；(2)直接对陀螺
输出信号进行消噪处理。常用方法有数字低通滤波、自适应滤波和小波变换
阈值滤波等。鉴于随机漂移的特性，事先不可能得到精确的统计特性，需要采用直接滤波方法对陀螺信号进行去噪处理。
[0004] 小波具有较好的时频分辨特性，因此被应用到光纤陀螺信号的去噪中。现有的基于小波理论的光纤陀螺信号分析研究中，都是基于2带小波的，去噪结果具有一定的局限性，主要体现在两方面：一方面，随着分解的进行，小波系数的长度以2倍关系递减，分解速度慢；另一方面，光纤陀螺信号循环移动一定
节拍后，其离散小波变换(Discrete Wavelet Transform，DWT)小波系数和尺度系数不能循环移动同样节拍，去噪效果差。
发明内容
[0005] 本发明的技术解决问题是：本发明的技术是针对已有的光纤陀螺信号分析方法中的不足，将交叠M带离散小波变换引入陀螺信号去噪。一方面，克服了基于2带小波的光纤陀螺信号处理分解速度慢的不足，M带小波理论为小波基的选取提供了更广泛的选择余地；在分解子带数相同条件下，M带离散小波变换(M-band Discrete Wavelet Transform，MDWT)中信号是按多通道进行分解的；另一方面，克服了基于2带小波进行光纤陀螺信号去噪时，光纤陀螺信号循环移动一定节拍后，DWT小波和尺度系数不能循环移动同样节拍的不足。
[0006] 本发明的技术解决方案是：本发明提出了一种基于交叠M带离散小波变换的光纤陀螺信号去噪方法，给出了计算有效的光纤陀螺信号分解、
门限阈值处理和重构
算法。具体步骤如下：
[0007] (1)建立光纤陀螺输出信号模型
[0008]
[0009] 其中，Xt为光纤陀螺输出信号；t＝0，...，N-1，N表示光纤陀螺输出信号的长度；ωie为地球自转角速率，ωie＝7.27×10-5rad/s；为测试点的地理纬度；εt为陀螺漂移，陀螺漂移由常值分量、周期分量和白噪声组成，即：
[0010] εt＝εd+Ωdsin(2πfd+θ0)+Wt
[0011] 其中，εd为常值分量，短时间内近似为一个常数；Ωd为周期分量的幅值；fd为周期分量的
频率；θ0为初始
相位；Wt为零均值高斯白噪声。式中的Ωdsin(2πfd+θ0)和Wt两项构成了光纤陀螺输出的随机漂移项。
[0012] (2)对光纤陀螺输出信号作OMDWT小波分解
[0013] 对步骤(1)得到的长度N的光纤陀螺输出信号Xt作OMDWT变换，将OMDWT变换所得到的第j尺度的OMDWT小波系数和尺度系数分别用N维矢量 和 表示，则 对应的元素 和 对应的元素 分别为：
[0014]
[0015] 式中，t＝0，...，N-1，N表示光纤陀螺信号的长度；l＝1，...，L，L＝2M表示第J0尺度上MDWT
滤波器的长度，M表示小波带数，即M带小波；i＝1，...，M-1；j＝J0，J0+1，...，J，J0表示小波变换的初始尺度，J表示小波变换的终止尺度； 和 分别为第j尺度上OMDWT的第i个小波滤波器系数和尺度滤波器系数；lmodN表示l对N取余。
[0016]
[0017] 式中， 和{gj，l}分别表示第j尺度上MDWT的第i个小波滤波器系数和尺度滤j/2 j波器系数。M 表示M的j/2次方；第j尺度上滤波器 的长度是Lj≡(M-1)(L-1)+1。
[0018] 对光纤陀螺信号Xt分别用周期滤波器进行滤波，可利用光纤陀螺信号在第j-1尺度上的尺度系数矢量 计算出光纤陀螺信号在第j尺度上的小波系数矢量 和尺度系数矢量 的元素。 和 对应的元素分别用 和 表示，则：
[0019]
[0020]
[0021] 式中， 和 分别表示OMDWT的第i个小波滤波器系数和尺度滤波器系数；j-1
表示t-M lmodN时刻第j-1尺度的尺度系数，j＝J0，J0+1，...，J。
[0022] 在上述变换中，OMDWT小波滤波器系数 和尺度滤波器系数 的求解如下：
[0023]
[0024] 式中，l＝1，...，L； 和{gl}分别为MDWT的第i个小波滤波器和尺度滤波器系数。
[0025] 步骤(2)对光纤陀螺信号作OMDWT小波分解时，按不同频带分解从而速度较快，且对光纤陀螺信号有更细的频带划分；光纤陀螺信号循环移动后，其OMDWT小波系数和尺度系数会循环移动相应节拍，有利于光纤陀螺的信号去噪。
[0026] (3)对OMDWT小波系数作门限阈值处理
[0027] 对步骤(2)中光纤陀螺信号分解后得到的小波系数 作门限阈值处理。利用门限阈值处理，将绝对值小于门限阈值的小波系数置零，但对绝对值大于门限阈值的小波系数作了收缩。
[0028]
[0029] 式中，Ts是门限阈值，10-5＜Ts＜10-3； 表示在尺度j上经过门限阈值处理后的第i个小波系数。
[0030] (4)对门限阈值处理后的光纤陀螺输出信号进行OMDWT小波重构
[0031]
[0032] 式中， 表示t时刻第j-1尺度的尺度系数； 和 分别表示第j-1
t+M lmodN时刻尺度j的尺度系数和门限阈值处理后的第i个小波系数。
[0033] 利用上式经过计算，最后在第J0尺度上得到重构的去噪后的光纤陀螺信号[0034] 本发明与
现有技术相比的优点在于：本发明在进行去噪时，按不同频带分解从而速度较快，且对光纤陀螺信号有更细的频带划分；(2)光纤陀螺信号循环移动后，其OMDWT小波系数和尺度系数会循环移动相应节拍，有利于光纤陀螺的信号去噪。
附图说明
具体实施方式
[0036] 如图1所示，本发明的实现方法如下：
[0037] 1、建立光纤陀螺输出信号模型
[0038]
[0039] 其中，Xt为光纤陀螺输出信号；t＝0，...，N-1，N表示光纤陀螺输出信号的长度；-5
ωie为地球自转角速率，ωie＝7.27×10 rad/s；为测试点的地理纬度；εt为陀螺漂移，陀螺漂移由常值分量、周期分量和白噪声组成，即
[0040] εt＝εd+Ωdsin(2πfd+θ0)+Wt (2)
[0041] 式中，εd为常值分量，短时间内近似为一个常数；Ωd为周期分量的幅值；fd为周期分量的频率；θ0为初始相位；Wt为零均值高斯白噪声。式中的Ωdsin(2πfd+θ0)和Wt两项构成了光纤陀螺输出的随机漂移项。
[0042] 2、对光纤陀螺输出信号作OMDWT小波分解
[0043] a.对光纤陀螺输出信号作OMDWT一步分解
[0044] 令
[0045]
[0046] 定义
[0047]
[0048]
[0049] OMDWT算 法 的 一 步 分 解 过 程 为：用 小 波 系 数 分 别 对Xt，进行MDWT分解，得到 和V1，并按式(5)排列。也就是说，进行M次MDWT分解后，把M组系数排列起来，即
[0050]
[0051] 其中，中的上标i表示此矩阵由M带小波变换的第i个小波系数组成。
[0052]
[0053] 例如，当M＝3，L＝6时
[0054]
[0055] 把式(6)中的小波滤波器系数 替换为gl，得到
[0056] 和V1的元素可表示为：
[0057]
[0058]
[0059] 可以看出，用小波滤波器 和尺度滤波器{gl}与时间序列Xt循环卷积，得到和V1。 和V1分别由长度为N的序列 和 的第M-1，2M-1，...，N-1个元素组成。
[0060] 把式(5)中的Xt替换为循环移动矢量 和 则有
[0061]
[0062] 定义
[0063]
[0064]
[0065] .
[0066] .
[0067] .
[0068]
[0069] 则
[0070]
[0071] 例如，当M＝3，L＝6时，有
[0072]
[0073]
[0074] 把 替换为gl，可得到 和 的元素可表示为：
[0075]
[0076] 的元素可表示为：
[0077]
[0078] 同样， 的元素可表示为：
[0079]
[0080] 的元素可表示为：
[0081]
[0082] 定义
[0083]
[0084]1/M
[0085] 将 和 的元素除以系数M ，并将它们按式(21)排列，可得到 用类似方法由V1， 和 得到。这里， 和 的元素
实际上是用OMDWT的滤波器 和 对Xt滤波后，得到的输出 和 通过把矩阵
和 中的每个 替换为 并将这些矩阵中的行按式(23)排列，得到
N×N矩阵
[0086]
[0087] 例如，当M＝3，L＝6时
[0088]
[0089] 故
[0090]
[0091] 把式(23)中的 替换为 可得到 则 最后，把OMDWT算法的一步分解表示为：
[0092] 其中
[0093] 通过这个过程，对光纤陀螺信号作OMDWT小波分解时，按不同频带分解从而速度较快，且对光纤陀螺信号有更细的频带划分；光纤陀螺信号循环移动后，其OMDWT小波系数和尺度系数会循环移动相应节拍。
[0094] b.求光纤陀螺信号Xt在第j尺度上的OMDWT变换
[0095] 对步骤1得到的长度N的光纤陀螺信号Xt作OMDWT变换，将第j尺度的OMDWT小波系数和尺度系数分别用N维矢量 和 表示，则 对应的元素 和 对应的元素 分别为：
[0096]
[0097] 式中，t＝0，...，N-1，N表示光纤陀螺信号的长度；l＝1，...，L，j＝J0，J0+1，...，J，J0表示小波变换的初始尺度，J表示小波变换的终止尺度； 和 分别为第j尺度上OMDWT的第i个小波滤波器系数和尺度滤波器系数；lmodN表示l对N取余。
[0098]
[0099] 和{gj，l}分别表示第j尺度上MDWT的第i个小波滤波器系数和尺度滤波器j系数。可以推得，滤波器 的长度是Lj≡(M-1)(L-1)+1。
[0100] 利用光纤陀螺信号在第j-1尺度上的尺度系数矢量 计算出光纤陀螺信号在第j尺度上的小波系数矢量 和尺度系数矢量 的元素。 和 对应的元素分别用和 表示，则：
[0101]
[0102]
[0103] 式中， 和 分别表示OMDWT的第i个小波滤波器系数和尺度滤波器系数；j-1
表示t-M lmodN时刻第j-1尺度的尺度系数，j＝J0，J0+1，...，J。
[0104] c.求OMDWT小波滤波器系数 和尺度滤波器系数
[0105]
[0106] 式中，l＝1，...，L；i＝1，...，M-1； 和{gl}分别为MDWT的第i个小波滤波器和尺度滤波器系数。
[0107] 3、对OMDWT小波系数作门限阈值处理
[0108] 对步骤2中光纤陀螺信号分解后得到的小波系数 作门限阈值处理。利用门限阈值处理，将绝对值小于门限阈值的小波系数置零，但对绝对值大于门限阈值的小波系数作了收缩。
[0109]
[0110] 式中，Ts是门限阈值，10-5＜Ts＜10-3； 表示在尺度j上经过门限阈值处理后的第i个小波系数。
[0111] 4、对门限阈值处理之后的光纤陀螺输出信号进行OMDWT小波重构光纤陀螺输出信号的OMDWT重构算法用下式描述：
[0112]
[0113] 式中，t＝0，1，...，N-1； 和 分别表示第t+Mj-1lmodN时刻尺度j的尺度系数和门限阈值处理后的第i个小波系数。
[0114] 利用上式经过计算，最后在第J0尺度上得到重构的去噪后的光纤陀螺信号[0115] 本发明基于OMDWT对光纤陀螺信号进行去噪，与基于2带小波的光纤陀螺信号去噪方法相比，计算速度快，噪声压缩比大。
[0116] 本发明
说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
