CN102937448B - 基于斜率的去除光纤陀螺脉冲噪声的方法 - Google Patents
基于斜率的去除光纤陀螺脉冲噪声的方法 Download PDFInfo
- Publication number
- CN102937448B CN102937448B CN201210417476.6A CN201210417476A CN102937448B CN 102937448 B CN102937448 B CN 102937448B CN 201210417476 A CN201210417476 A CN 201210417476A CN 102937448 B CN102937448 B CN 102937448B
- Authority
- CN
- China
- Prior art keywords
- data
- group
- slope
- groups
- fiber
- 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.)
- Active
Links
Landscapes
- Gyroscopes (AREA)
Abstract
本发明公开了一种基于斜率的去除光纤陀螺脉冲噪声的方法,属于光纤陀螺信号处理技术领域。所述方法首先对采样数据均分为三组,在分组的过程中,组与组之间有间隙;然后将各个组内数据的幅值相加,计算累加值;最后根据累加值之间的差值计算两两之间的斜率,比较斜率绝对值的大小,保留组成斜率最小的两组数据即为需要的采集数据。本发明在保证精度的前提下,由于只是差值计算,计算较为快速,减少了工作量。在硬件主要有FPGA实现的情况下,斜率的计算对硬件资源的占用较少,便于硬件更加深度的应用;数据处理实时性更好,改进空间大,差值计算斜率后,对舍去的一组数据进行分析,便于有针对性的找出原因,有利于提高光纤陀螺的精度。
Description
技术领域
本发明主要涉及光纤陀螺信号处理技术领域,特别适用于光纤陀螺信号处理去除瞬时尖峰脉冲噪声的影响,具体地说,是指一种基于斜率的去除光纤陀螺脉冲噪声的方法。
背景技术
光纤陀螺是一种用来检测载体运动角速率的传感器,利用光学上的萨格奈克效应与光纤传输原理,克服了机电陀螺与激光陀螺的一系列缺点,具有全固态、可靠性高、动态范围大、频带宽、体积小、振动性能好等优点,被广泛应用于航天、航空、航海、石油勘探等领域。
在光纤陀螺信号处理部分采用探测器将光信号转换成电信号处理,通过电信号计算出相位差,进而得到载体运动角速率。因此光纤陀螺的信号处理及其噪声的抑制显得尤为重要。
在光电转换过程中,探测器输出的信号中有用信号是微弱信号,噪声比较大,有用信号常常淹没在噪声中,现有技术中一般采用自相关平均法处理信号,即基于信号和噪声的统计特性进行检测。
发明内容
本发明的目的是提供一种利用斜率简单快速剔除光纤陀螺信号处理尖峰噪声的方法,并保证光纤陀螺的测量精度。
在光电转换的过程中,由于环境干扰和电路串扰的影响,探测器输出的信号会有尖峰脉冲噪声,因此在计算自相关平均时需要剔除尖峰脉冲的影响,本发明提供基于斜率的去除光纤陀螺脉冲噪声的方法,具体的处理方法步骤为:
第一步,对探测器的输出信号采样后,把采样得到的数字信号数据均分为三组。在均分的过程中,为了避免尖峰脉冲同时落在两组内的情况,各组之间不连续,即组与组之间有间隙。
所述的间隙大于尖峰脉冲的宽度。
第二步,将各个组内数据的幅值相加,计算累加值。
第三步,在第二步的基础上,根据计算得到的三个累加值,计算这各组两两之间的累加值的差值,进而计算斜率,比较斜率绝对值的大小,保留组成斜率绝对值最小的两组数据对应的组内数据即为需要的采集数据。
在探测器采集信号的过程中,采集到的数据中包含了有用信号和噪声,有用信号特别微弱,噪声主要是随机白噪声,从统计学来看白噪声是随机不相关的,对其幅值时间平均,组与组之间幅值相差特别小,其幅值是均匀的,如果有尖峰脉冲噪声的话,幅值平均就会明显偏大。因此上述方法中采用的斜率方式可以有效去除光纤陀螺脉冲噪声。
本发明用于光纤陀螺信号处理去除尖峰脉冲噪声的方法,与现有的方法相比,优点在于:
(1)由于尖峰脉冲幅值很大,把采集到的所有数据平均处理会由于尖峰脉冲的影响而使均值偏大,基于斜率的处理去除了尖峰脉冲的影响,结果更准确。
(2)数据分组不连续,排除了尖峰脉冲同时在两组内的情况,分组简单可行。
(3)在保证精度的前提下,由于只是差值计算,计算较为快速,减少了工作量。
(4)在硬件主要有FPGA实现的情况下,斜率的计算对硬件资源的占用较少,便于硬件更加深度的应用。
(5)有效数据的选择是基于“三选二”的原则,可靠性好,数据处理实时性好。
(6)数据处理时,组与组之间的间隙可以根据理论计算或者实验结果调整,数据的弹性更大,处理更为灵活。
(7)根据差值计算斜率后,对舍去的各组数据进行分析,便于有针对性的找出原因,有利于提高光纤陀螺的精度,改进空间大。
附图说明
图1是正常情况下探测器输出的信号;
图2是有尖峰脉冲噪声时的信号;
图3A是本发明信号处理分组方法示意图;
图3B是本发明中进行各组间斜率计算示意图;
图4是本发明提供方法的流程图;
图5是数据分组分界点落在尖峰脉冲内时的示意图;
图6是组间分界点落在尖峰脉冲内时的斜率计算示意图。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明提供一种基于斜率的去除光纤陀螺脉冲噪声的方法,如图4所示流程,包括如下步骤:
第一步,对探测器的输出信号进行采样,把采样得到的数字信号数据均分为三组A1、A2、A3,如图3A。在分组的过程中,为了避免尖峰脉冲同时落在两组内,三个分组的数据之间不连续,即组与组之间有间隙。所述的间隙大小要满足大于尖峰脉冲的宽度。
如图1所示,正常情况下没有尖峰脉冲噪声时的信号,如图2所示,探测器输出有尖峰脉冲噪声时的信号,本发明提供的方法就是要去除自相关平均时的尖峰噪声。
第二步,将各个组内数字信号的幅值相加,计算三个幅值累加值SA1、SA2、SA3。
第三步,在第二步的基础上,根据计算得到的三个幅值累加值SA1、SA2、SA3,计算这三组两两之间的幅值累加值的差值,根据所述的差值计算斜率,得到三个斜率值K13、K12、K23:
K13=(SA1-SA3)/t;
K12=(SA1-SA2)/t;
K23=(SA2-SA3)/t;
其中t代表时间。
比较斜率K13、K12、K23的绝对值|K13|、|K12|、|K23|的大小,如图3B所示,保留组成斜率的绝对值最小的两组幅值累加值对应的组内数据即为需要的采集数据,即数据处理部分的:“三选二”。本实施例中如图3B所示,斜率|K12|最小,所以A1组和A2组的数据即为所需要的数据,也就是这两组数据中不包含尖峰噪声。
假如在分组的过程中,两组的分界点同时落在尖峰脉冲上,即图5所示,A1组和A2组的分界点同时落在尖峰脉冲上,此时斜率计算如图6所示,则组成斜率绝对值最小的两组为A2组和A3组,但数据明显不符合要求,所以在分组的过程中,组与组之间不连续,以防止分界点同时落在尖峰脉冲上。
如果没有尖峰脉冲的情况下,数据内噪声主要是白噪声,统计平均发现白噪声的幅值是均匀的,如果组内有尖峰脉冲,则组内数据的累加值会明显偏大,那么包含尖峰脉冲的一组数据与其他组数据做斜率计算就会明显大于不包含尖峰脉冲噪声的任意两组之间的斜率,这样就可以去除含有尖峰脉冲的一组数据。
Claims (1)
1.基于斜率的去除光纤陀螺脉冲噪声的方法,其特征在于,所述的方法包括如下步骤:
第一步,探测器输出信号进行采样,把采样得到的数据均分为三组;数据均分的过程中,均分的各组数据之间不连续,即组与组之间有间隙;所述的间隙大于尖峰脉冲的宽度;
第二步,将各个组内数据的幅值相加,计算幅值累加值;
第三步,根据计算得到的三个幅值累加值,计算这三组两两之间的幅值累加值差值,进而计算斜率,比较斜率绝对值的大小,保留组成斜率绝对值最小的两组数据即为需要的采集数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210417476.6A CN102937448B (zh) | 2012-10-26 | 2012-10-26 | 基于斜率的去除光纤陀螺脉冲噪声的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210417476.6A CN102937448B (zh) | 2012-10-26 | 2012-10-26 | 基于斜率的去除光纤陀螺脉冲噪声的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102937448A CN102937448A (zh) | 2013-02-20 |
CN102937448B true CN102937448B (zh) | 2015-04-15 |
Family
ID=47696359
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210417476.6A Active CN102937448B (zh) | 2012-10-26 | 2012-10-26 | 基于斜率的去除光纤陀螺脉冲噪声的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102937448B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106885565B (zh) * | 2017-02-14 | 2020-02-07 | 北京航空航天大学 | 一种基于中值滤波的干涉式光纤陀螺输出信号附加零偏的抑制方法 |
CN113031524B (zh) * | 2021-02-07 | 2022-06-17 | 南京航空航天大学 | 一种基于三次样条的压装力包络线生成方法 |
CN113739783B (zh) * | 2021-11-08 | 2022-03-01 | 华中光电技术研究所(中国船舶重工集团公司第七一七研究所) | 一种光纤陀螺振铃效应误差抑制系统及抑制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5684589A (en) * | 1995-08-28 | 1997-11-04 | Litton Systems, Inc. | Loop controller for fiber optic gyro with distributed data processing |
US5854678A (en) * | 1996-06-28 | 1998-12-29 | Honeywell Inc. | Three-axis fiber optic gyroscope having a single source and multi-coupler configuration |
CN102141403A (zh) * | 2010-12-17 | 2011-08-03 | 北京航空航天大学 | 基于小波阈值去噪、中值滤波和均值滤波的实时混合去噪方法 |
-
2012
- 2012-10-26 CN CN201210417476.6A patent/CN102937448B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5684589A (en) * | 1995-08-28 | 1997-11-04 | Litton Systems, Inc. | Loop controller for fiber optic gyro with distributed data processing |
US5854678A (en) * | 1996-06-28 | 1998-12-29 | Honeywell Inc. | Three-axis fiber optic gyroscope having a single source and multi-coupler configuration |
CN102141403A (zh) * | 2010-12-17 | 2011-08-03 | 北京航空航天大学 | 基于小波阈值去噪、中值滤波和均值滤波的实时混合去噪方法 |
Non-Patent Citations (2)
Title |
---|
光纤陀螺漂移误差动态Allan方差分析;李颖等;《光电子激光》;20080228;第19卷(第2期);第183-186页 * |
基于动态Allan 方差的机抖激光陀螺随机误差研究;魏国等;《中国激光》;20101231;第37卷(第12期);第2975-2979页 * |
Also Published As
Publication number | Publication date |
---|---|
CN102937448A (zh) | 2013-02-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102937448B (zh) | 基于斜率的去除光纤陀螺脉冲噪声的方法 | |
CN104008622B (zh) | 基于短时能量和过零率的光纤周界安防系统端点检测方法 | |
CN102928873B (zh) | 基于四维能量聚焦的地面微地震定位方法 | |
CN104467969B (zh) | 分数阶傅里叶变换测量光纤链路色散的方法 | |
CN102353393A (zh) | 基于π/2相位调制的干涉型光传感器的正交解调装置 | |
CN108955734A (zh) | 一种光纤f-p温度/压力复合传感器的腔长解调方法 | |
CN104132693A (zh) | 相位otdr系统中振动信号位置和频率的同时提取方法 | |
CN102650526A (zh) | 一种基于相位比较调频连续波光纤陀螺的开环检测电路 | |
CN106885565B (zh) | 一种基于中值滤波的干涉式光纤陀螺输出信号附加零偏的抑制方法 | |
CN102183249B (zh) | 一种光纤陀螺仪的萨格奈克相移跟踪方法 | |
CN103063237B (zh) | 一种使编码器抗干扰的方法及装置 | |
CN104122412B (zh) | 一种基于北斗二代速度信息的加速度计标定方法 | |
RU2405160C1 (ru) | Устройство для измерения ускорений | |
CN107228684B (zh) | 基于自相关滤波技术的psd输出有用信号提取方法 | |
CN113219539B (zh) | 光纤声波传感地震数据中井筒波干扰的去除方法及装置 | |
CN109059908B (zh) | 一种运动平台光电跟踪系统惯性传感器信号中扰动信息提取方法 | |
CN108444464B (zh) | 一种光纤陀螺仪中的光源相对强度噪声抑制方法 | |
CN103197161A (zh) | 一种微弱直流信号中的脉冲干扰检测方法 | |
WO2014063989A1 (en) | Motion sensor assembly for determining the angular velocity of a moving contrast in its field of view with a high accuracy | |
CN115420272A (zh) | 一种光纤陀螺光源相对强度噪声自适应抑制实现方法 | |
CN110530496A (zh) | 一种基于数字滤波和二元脉冲调制的光纤传感解调方法 | |
CN110285881A (zh) | 一种基于全相位滤波的密集谱频率估计法 | |
CN102082606A (zh) | 一种光纤电流互感器相差检测及处理装置 | |
Yu et al. | OSNR monitoring by using single sampling channel generated 2-D phase portrait | |
US20080175338A1 (en) | Apparatus and method for detecting null symbols |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |