CN107063306A - 一种基于改进的eemd和排列熵的光纤陀螺振动补偿算法 - Google Patents
一种基于改进的eemd和排列熵的光纤陀螺振动补偿算法 Download PDFInfo
- Publication number
- CN107063306A CN107063306A CN201710243138.8A CN201710243138A CN107063306A CN 107063306 A CN107063306 A CN 107063306A CN 201710243138 A CN201710243138 A CN 201710243138A CN 107063306 A CN107063306 A CN 107063306A
- Authority
- CN
- China
- Prior art keywords
- imf
- signal
- vibration
- optical fibre
- entropy
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
- G01C25/005—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Gyroscopes (AREA)
Abstract
本发明公开了一种基于改进的集合经验模态分解(EEMD)和排列熵的光纤陀螺振动补偿算法,属于数字信号处理领域。为了消除传统经验模态分解的模态混叠缺陷,有效分离噪声和陀螺振动信号,进而对陀螺输出进行补偿,达到消除噪声和振动的目的。算法首先通过对人工白噪声进行EMD分解,多次利用其第一个模态分量对陀螺信号进行分解;然后利用排列熵对各个IMF求取熵值,基于熵值找出噪声相关IMF和振动相关IMF;去掉噪声相关IMF和振动相关IMF后,剩下的IMF当做光纤陀螺的最后输出信号,从而达到去噪和消除振动影响的目的。
Description
技术领域
本发明涉及一种基于改进的集合经验模态分解(EEMD)和排列熵的光纤陀螺振动补偿算法,属于数字信号处理领域。
背景技术
光纤陀螺是一种高精度惯性传感器,但容易受到环境的干扰,对外界振动情况较为敏感,也是影响陀螺精度及其应用的主要因素之一。如何检测并消除振动的影响成为工程应用中的一个重点和难点。
角陀螺的振动光纤陀螺的一种主要振动形式,它是一种旋转振动,会直接导致陀螺的输出误差。载体的振动上分为周期振动和随机振动,前者由周期的扰动和不平衡力造成,后者则是由载体的航速、航向、操纵和航行环境造成的。相比而言,周期垂向振动占主导作用。要消除振动的干扰,首先要准确地分理出振动信号。因此在理论和工程中,如何检测并消除陀螺的振动信息是提高陀螺精度的关键,也是当前的研究热点和难点。除了传感器层面的影响之外,导航初始对准和导航解算都会受到陀螺振动的影响,因此,陀螺的振动问题不止是器件问题,更是***问题。消除振动问题,主要有物理方法和软件算法,物理方法虽效果较好,却带来了成本、体积、散热等问题;相比而言,软件具有更高的灵活性及经济性。经验模态分解(EMD)是一种基于时间尺度的自适应分解方法,近些年得到了广泛的应用。但是EMD存在模态混叠和端点效应等问题。集合经验模态分解(Ensemble EmpiricalMode Decomposition,EEMD)以改进的噪声辅助为手段,较好地解决了EMD的问题。EEMD也有一些缺点,比如噪声幅值的选择和迭代的次数都是经验化的。当混合信号的频率较为接近时(尤其是频率的倍数小于2),EEMD依然会有信号混叠的问题。为了提高EEMD消除模态混叠的性能,噪声对分解效果的影响需要深入研究。排列熵算法是一种分析***非线性的新方法,能够以较高的分辨率检测信号的非线性度。通过排列熵可以分辨出噪声和信号,通过这种方法可以达到去噪的作用。
基于上述陀螺振动的问题,更多地关注振动的软件算法,可以在不增加成本、体积等问题的前提下,提高解决问题的精度和实时性。
发明内容
发明目的:为了解决上述问题,本发明目的在于提供一种基于改进的集合经验模态分解和排列熵的光纤陀螺振动补偿算法,根据陀螺振动的特征,结合创新的经验模态分解算法来提取和补偿陀螺的振动误差,能够消除传统经验模态分解的模态混叠缺陷,有效分离噪声和陀螺振动信号,进而对陀螺输出进行补偿,达到消除噪声和振动的目的。
技术方案:为实现上述发明目的,本发明采用如下技术方案:
一种基于改进的集合经验模态分解和排列熵的光纤陀螺振动补偿算法,包括如下步骤:
(1)构造l个人工白噪声并用经验模态分解(EMD)进行分解,得到第i个人造白噪声w(i)的第一个模态分量(IMF)E1(w(i)),i=1,...,l,l≥2;
(2)使用步骤(1)得到的E1(w(i))对光纤陀螺信号用集合经验模态分解(EEMD)进行分解,得到光纤陀螺信号的前m个IMF,m≥2;
(3)对于光纤陀螺剩余的信号用EMD分解,进而得到光纤陀螺信号的完整分解;
(4)计算光纤陀螺信号各个IMF的排列熵的值;
(5)将排列熵值大于设定阈值的IMF作为噪声相关IMF,排列熵值曲线上波谷对应的IMF作为振动相关IMF,去除噪声相关IMF和振动相关IMF后,剩下的作为光纤陀螺最终输出信号。
进一步地,所述步骤(2)中包括:
(2.1)设光纤陀螺原始信号为x,构造分解信号为x(i)=x+β(E1(w(i))),其中β>0,i=1,...,l,分别使用EMD对x(i)进行分解,再对得到的l个第一IMF求均值,得到光纤陀螺原始信号的第一个IMF,即分解得到IMF1,剩余的信号作为残留信号r1=x-IMF1;
(2.2)以残留信号替代步骤(2.1)中的陀螺原始信号,进行构造分解信号、分解及求均值操作,重复步骤(2.1)m-1次,得到第m个模态分量剩余信号
有益效果:本发明公开的一种基于改进的集合经验模态分解和排列熵的光纤陀螺振动补偿算法,通过对人工白噪声进行EMD分解,多次利用其第一个模态分量(IMF1)对陀螺信号进行分解,可以有效地消除模态混叠,从而分离出振动信号。然后利用排列熵对各个IMF求取熵值,基于排列熵值判断出噪声相关IMF和振动相关IMF,去掉噪声相关IMF和振动相关IMF后,剩下的IMF当做光纤陀螺的最后输出信号,从而达到去噪和消除振动影响的目的。本发明算法,根据陀螺振动的特征,结合创新的经验模态分解算法来提取和补偿陀螺的振动误差,能够消除传统经验模态分解的模态混叠缺陷,有效分离噪声和陀螺振动信号。
附图说明
图1为本发明的算法的流程图。
具体实施方式
下面结合附图和实施例对本发明做进一步详细说明。应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
如图1所示,本发明实施例公开的一种基于改进的集合经验模态分解和排列熵的光纤陀螺振动补偿算法。首先,对人工白噪声进行EMD分解,只取其第一个IMF作为分解噪声带入EEMD算法得到前m个IMF(其中m和信号的噪声水平相关,一般取m大于2)。然后,剩下的残留信号使用EMD进行分解,进而得到光纤陀螺振动信号的完整分解结果。第三步,对各个IMF求排列熵,然后对IMF进行分类,具体为a)大于经验阈值的为噪声。b)波谷的为振动信号。C)其余的为陀螺输出信号。滤除上述的噪声IMF和振动IMF后,得到的最后信号即为光纤陀螺的输出,进而达到滤波和消除振动影响的目的。详细的算法流程如下:
第一步,利用EMD分解人工白噪声。
假设w(i)是第i个人造白噪声,E1(w(i))表示用EMD方法对w(i)分解得到的第1个IMF,其中i=1,...,l,l≥2,一般可取100的整数倍,如100或200。
第二步,使用第一步得到的E1(w(i))分解得到前m个IMF。
假设光纤陀螺原始输出为x,构造分解信号为
x(i)=x+β(E1(w(i))) (2a)
其中β>0,为噪声幅值。分别使用EMD对x(i)进行分解,得到各个IMFk (i),其中k=1,...,K。再对IMF1 (i)求均值,即可得到光纤陀螺原始信号的第一个IMF,即
分解得到第一个IMF后,剩余的信号作为残留信号r1,即有
r1=x-IMF1 (2c)
以残留信号替代式(2a)中的陀螺原始输出信号,构造分解信号,进行多次分解,求均值操作,重复共m-1次,其中m≥2,共得到光纤陀螺原始信号的前m个IMF,以及残留信号
与第一个IMF类似,第m个IMF由下式求得
第三步,对于剩余信号rm运用EMD分解得到剩下的IMF。
第四步,通过计算各个IMF分离出噪声、振动和陀螺输出信号。
计算各个IMF的排列熵的值,大于设定阈值(阈值为经验值,一般可在[0.9,1]范围内取值)的IMF作为噪声,记为IMFnoise;在排列熵值的曲线上寻找波谷,求得其所在列的索引值为col,其对应的IMF记为IMFcol,则光纤陀螺在消除噪声和振动后的输出xout为
xout=x-∑IMFnoise-∑IMFcol (4a)。
Claims (3)
1.一种基于改进的集合经验模态分解和排列熵的光纤陀螺振动补偿算法,其特征在于,包括如下步骤:
(1)构造l个人工白噪声并用经验模态分解(EMD)进行分解,得到第i个人造白噪声w(i)的第一个模态分量(IMF)E1(w(i)),i=1,...,l,l≥2;
(2)使用步骤(1)得到的E1(w(i))对光纤陀螺信号用集合经验模态分解(EEMD)进行分解,得到光纤陀螺信号的前m个IMF,m≥2;
(3)对于光纤陀螺剩余的信号用EMD分解,进而得到光纤陀螺信号的完整分解;
(4)计算光纤陀螺信号各个IMF的排列熵的值;
(5)将排列熵值大于设定阈值的IMF作为噪声相关IMF,排列熵值曲线上波谷对应的IMF作为振动相关IMF,去除噪声相关IMF和振动相关IMF后,剩下的作为光纤陀螺最终输出信号。
2.根据权利要求1所述的一种基于改进的集合经验模态分解和排列熵的光纤陀螺振动补偿算法,其特征在于,所述步骤(2)包括:
(2.1)设光纤陀螺原始信号为x,构造分解信号为x(i)=x+β(E1(w(i))),其中β>0,i=1,...,l,分别使用EMD对x(i)进行分解,再对得到的l个第一IMF求均值,得到光纤陀螺原始信号的第一个IMF,即分解得到IMF1,剩余的信号作为残留信号r1=x-IMF1;
(2.2)以残留信号替代步骤(2.1)中的陀螺原始信号,进行构造分解信号、分解及求均值操作,重复步骤(2.1)m-1次,得到第m个模态分量剩余信号
3.根据权利要求1所述的一种基于改进的集合经验模态分解和排列熵的光纤陀螺振动补偿算法,其特征在于,所述步骤(5)中将排列熵值大于0.9的IMF作为噪声相关IMF。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710243138.8A CN107063306A (zh) | 2017-04-14 | 2017-04-14 | 一种基于改进的eemd和排列熵的光纤陀螺振动补偿算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710243138.8A CN107063306A (zh) | 2017-04-14 | 2017-04-14 | 一种基于改进的eemd和排列熵的光纤陀螺振动补偿算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107063306A true CN107063306A (zh) | 2017-08-18 |
Family
ID=59599504
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710243138.8A Pending CN107063306A (zh) | 2017-04-14 | 2017-04-14 | 一种基于改进的eemd和排列熵的光纤陀螺振动补偿算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107063306A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107885940A (zh) * | 2017-11-10 | 2018-04-06 | 吉林大学 | 一种用于分布式光纤振动传感***的信号特征提取方法 |
CN108506736A (zh) * | 2018-03-19 | 2018-09-07 | 浙江理工大学 | 一种入侵信号诊断***和方法 |
CN108759871A (zh) * | 2018-07-03 | 2018-11-06 | 哈尔滨工业大学 | 一种基于改进emd预处理算法的捷联惯性导航***粗对准方法 |
CN109443393A (zh) * | 2018-12-11 | 2019-03-08 | 中国人民解放***箭军工程大学 | 一种基于盲分离算法的捷联惯导信号提取方法及*** |
CN113702666A (zh) * | 2021-08-03 | 2021-11-26 | 哈尔滨工程大学 | 一种光纤陀螺惯性测量单元信号联合降噪方法 |
CN114088077A (zh) * | 2021-12-10 | 2022-02-25 | 哈尔滨工业大学 | 一种改进的半球谐振陀螺信号去噪方法 |
CN114896559A (zh) * | 2022-04-29 | 2022-08-12 | 盐城工学院 | 一种基于排列熵的自适应噪声完备集合经验模式分解方法 |
CN115031759A (zh) * | 2022-02-25 | 2022-09-09 | 长光卫星技术股份有限公司 | 基于等效噪声带宽法的在轨光纤陀螺噪声获取方法、可读存储介质及电子设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101566483A (zh) * | 2009-05-22 | 2009-10-28 | 哈尔滨工程大学 | 光纤陀螺捷联惯性测量***振动误差补偿方法 |
US20100092028A1 (en) * | 2008-10-10 | 2010-04-15 | National Central University | Data Decomposition Method and Computer System Therefrom |
CN103091096A (zh) * | 2013-01-23 | 2013-05-08 | 北京信息科技大学 | 基于eemd和小波包变换的早期故障敏感特征提取方法 |
CN104019831A (zh) * | 2014-06-20 | 2014-09-03 | 哈尔滨工业大学 | 基于emd和熵权的陀螺仪故障诊断方法 |
CN104931197A (zh) * | 2015-06-29 | 2015-09-23 | 吉林大学 | 基于eemd的自动平衡机振动信号处理方法 |
CN105371836A (zh) * | 2015-12-18 | 2016-03-02 | 哈尔滨工业大学 | 基于eemd和fir的混合型光纤陀螺信号滤波方法 |
-
2017
- 2017-04-14 CN CN201710243138.8A patent/CN107063306A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100092028A1 (en) * | 2008-10-10 | 2010-04-15 | National Central University | Data Decomposition Method and Computer System Therefrom |
CN101566483A (zh) * | 2009-05-22 | 2009-10-28 | 哈尔滨工程大学 | 光纤陀螺捷联惯性测量***振动误差补偿方法 |
CN103091096A (zh) * | 2013-01-23 | 2013-05-08 | 北京信息科技大学 | 基于eemd和小波包变换的早期故障敏感特征提取方法 |
CN104019831A (zh) * | 2014-06-20 | 2014-09-03 | 哈尔滨工业大学 | 基于emd和熵权的陀螺仪故障诊断方法 |
CN104931197A (zh) * | 2015-06-29 | 2015-09-23 | 吉林大学 | 基于eemd的自动平衡机振动信号处理方法 |
CN105371836A (zh) * | 2015-12-18 | 2016-03-02 | 哈尔滨工业大学 | 基于eemd和fir的混合型光纤陀螺信号滤波方法 |
Non-Patent Citations (5)
Title |
---|
XIAOYUAN ZHANG,ET AL.: "A novel bearing fault diagnosis model integrated permutation entropy, ensemble empirical mode decomposition and optimized SVM", 《MEASUREMENT》 * |
XIYUAN CHEN ET AL.: "Efficient modeling of fiber optic gyroscope drift using improved EEMD and extreme learning machine", 《SIGNAL PROCESSING》 * |
刘涛 等: "基于改进的集合经验模态方法振动信号分解", 《上海交通大学学报》 * |
李永军 等: "一种振动信号的EMD去噪方法及其仿真分析", 《自动化仪表》 * |
杨宏: "经验模态分解及其在水声信号处理中的应用", 《中国博士学位论文全文数据库信息科技辑》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107885940A (zh) * | 2017-11-10 | 2018-04-06 | 吉林大学 | 一种用于分布式光纤振动传感***的信号特征提取方法 |
CN108506736A (zh) * | 2018-03-19 | 2018-09-07 | 浙江理工大学 | 一种入侵信号诊断***和方法 |
CN108506736B (zh) * | 2018-03-19 | 2019-11-29 | 浙江理工大学 | 一种入侵信号诊断***和方法 |
CN108759871A (zh) * | 2018-07-03 | 2018-11-06 | 哈尔滨工业大学 | 一种基于改进emd预处理算法的捷联惯性导航***粗对准方法 |
CN109443393A (zh) * | 2018-12-11 | 2019-03-08 | 中国人民解放***箭军工程大学 | 一种基于盲分离算法的捷联惯导信号提取方法及*** |
CN113702666A (zh) * | 2021-08-03 | 2021-11-26 | 哈尔滨工程大学 | 一种光纤陀螺惯性测量单元信号联合降噪方法 |
CN114088077A (zh) * | 2021-12-10 | 2022-02-25 | 哈尔滨工业大学 | 一种改进的半球谐振陀螺信号去噪方法 |
CN115031759A (zh) * | 2022-02-25 | 2022-09-09 | 长光卫星技术股份有限公司 | 基于等效噪声带宽法的在轨光纤陀螺噪声获取方法、可读存储介质及电子设备 |
CN114896559A (zh) * | 2022-04-29 | 2022-08-12 | 盐城工学院 | 一种基于排列熵的自适应噪声完备集合经验模式分解方法 |
CN114896559B (zh) * | 2022-04-29 | 2023-03-07 | 盐城工学院 | 一种基于排列熵的自适应噪声完备集合经验模式分解方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107063306A (zh) | 一种基于改进的eemd和排列熵的光纤陀螺振动补偿算法 | |
CN101087365B (zh) | 一种滤除图像混合噪声的方法 | |
CN109883392B (zh) | 一种基于相位补偿的捷联惯导升沉测量方法 | |
CN101840529B (zh) | 基于局部可变集成神经网络的光纤陀螺随机漂移建模方法 | |
CN108759871B (zh) | 一种基于改进emd预处理算法的捷联惯性导航***粗对准方法 | |
CN106409351A (zh) | 一种核电厂松动部件振动信号的处理方法及装置 | |
CN110940524A (zh) | 一种基于稀疏理论的轴承故障诊断方法 | |
CN117490834A (zh) | 振动频谱峰值捕获方法 | |
CN109425897B (zh) | 消除地震资料野值干扰的方法及*** | |
CN115743101A (zh) | 车辆轨迹预测方法、轨迹预测模型训练方法和装置 | |
CN113702666A (zh) | 一种光纤陀螺惯性测量单元信号联合降噪方法 | |
CN116380462A (zh) | 一种多尺度均值变分模态分解的轴承振动信号降噪方法 | |
CN112183407B (zh) | 一种基于时频域谱减法的隧道地震波数据去噪方法及*** | |
CN104502992B (zh) | 一种基于空时过采样扫描的弱小点目标精确定位方法及*** | |
CN111596358B (zh) | 一种多次波的压制方法、设备以及*** | |
CN106840202A (zh) | 一种陀螺振动信号提取与补偿方法 | |
CN107301431A (zh) | 一种基于权重的非极大值抑制的方法及装置 | |
CN102073422B (zh) | 电容式触控板的控制电路及方法 | |
CN105301655B (zh) | 一种共成像点道集线性噪声去除方法及装置 | |
CN114444548A (zh) | 一种基于鲁棒性主成分分析的小波声音降噪方法及*** | |
CN114088077B (zh) | 一种改进的半球谐振陀螺信号去噪方法 | |
CN116339545B (zh) | 一种对触摸屏干扰信号处理方法与*** | |
CN103528608A (zh) | 一种基于ccd分区定位的多模光纤散斑传感*** | |
CN116626650B (zh) | 一种单光子探测的运动目标全波形恢复方法与*** | |
CN203758528U (zh) | 一种基于ccd分区定位的多模光纤散斑传感*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170818 |