CN102098101B - 波峰寻找方法及系统 - Google Patents
波峰寻找方法及系统 Download PDFInfo
- Publication number
- CN102098101B CN102098101B CN2011100005764A CN201110000576A CN102098101B CN 102098101 B CN102098101 B CN 102098101B CN 2011100005764 A CN2011100005764 A CN 2011100005764A CN 201110000576 A CN201110000576 A CN 201110000576A CN 102098101 B CN102098101 B CN 102098101B
- Authority
- CN
- China
- Prior art keywords
- data
- crest
- current
- state
- attitude
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Mobile Radio Communication Systems (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种波峰寻找方法,包括如下步骤:(1)对光栅反射波形信号进行高速A/D采样,获取一组离散波形数据;(2)根据离散波形数据的当前数据和前一数据实时切换系统当前状态,直至系统当前状态为完成态或错误态,状态转换过程中获取离散波形数据中的有效数据;(3)当系统当前状态为错误态时,放弃当前采样的离散波形数据;(4)当系统当前状态为完成态时,对有效数据进行数学期望和分位数计算,获得波峰位置值。本方法能完成对代表光栅反射波长信号的位移量的离散数字信号的分析和处理,以获得具有高重复性、高稳定度的光栅解调结果。另外本方法还包括对波峰位置值进行有效审核的步骤。本发明同时公开了一种波峰寻找系统。
Description
技术领域
本发明涉及光纤光栅数据传输系统,尤其涉及光纤光栅数据传输系统中的高重复性、稳定性寻找波峰的方法和系统。
背景技术
近年来,光纤光栅技术在各行各业中的应用越来越广泛,因此对光栅解调仪的需求也随之扩大。
光栅传感器用于感知物体的表面由应力、温度、振动等导致的形变等物理参数,根据各种物理参数产生光信号(光栅反射波长信号),其波长信号位移量的变化程度和速度由物理参数决定。
光栅解调仪将光栅传感器产生的光信号转换为电子设备能够识别的电信号,并将模拟电信号经过A/D采样后转换为数字电信号,所述数字电信号代表了光栅反射波长信号的位移量。
波峰寻找方法分析和处理该数字电信号,获得具有高重复性、高稳定度的光栅解调结果。
已有的各种波峰寻找方法大多采用多项式-高斯公式拟合法指令或者高斯公式拟合法指令,这种方法一方面需要有足够的内存来存储采集的一个周期的数据,并且所述指令需要消耗很多的指令周期,不仅影响了波峰寻找的速度,而且限制了光纤光栅数据解调系统的速度;另一方面多以寻找波峰的峰值点的位置来代表波峰的位置,在实际使用过程中,由于各种干扰因素的存在会对波峰的局部造成影响,进而会造成峰值点的波动和起伏,使得对同一波形信号的多次寻峰结果会存在差异,降低了光栅解调仪的重复性和稳定性。
因此,有必要提供一种改进的波峰寻找方法和系统来克服现有技术的缺陷。
发明内容
本发明的目的是提供一种波峰寻找方法及系统,完成对代表光栅反射波长信号的位移量的离散数字信号的分析和处理,以获得具有高重复性、高稳定度的光栅解调结果。
为了实现上述目的,本发明提供了一种波峰寻找方法,包括如下步骤:
(1)对光栅反射波形信号进行高速A/D采样,获取一组离散波形数据;
(2)根据所述离散波形数据的当前数据和前一数据实时切换系统当前状态,直至系统当前状态为完成态或错误态,状态转换过程中获取离散波形数据中的有效数据;
(3)当系统当前状态为错误态时,放弃当前采样的离散波形数据;
(4)当系统当前状态为完成态时,对所述有效数据进行数学期望和分位数计算,获得波峰位置值。
在本发明的一个实施例中,所述步骤(2)具体为:
当当前采样值不大于阀值时,前一状态是初始态,将系统当前状态设置为是初始态;
当当前数据和前一数据均大于采样阀值,且当前数据小于前一数据、前一数据不等于采样阀值时,将系统当前状态设置为下降态;
当当前数据和前一数据均大于阀值,前一状态是下降态,且当前数据大于前一数据或当前数据等于前一数据时,将系统当前状态设置为错误态;
当当前数据和前一数据均大于阀值,当前数据大于前一数据,且前一状态不是下降态时,将系统当前状态设置为上升态;
当当前数据和前一数据均大于阀值,当前数据等于前一数据,且前一状态不是下降态时,将系统当前状态设置为平稳态;
当当前数据小于阀值,前一数据大于阀值,且前一状态不是错误态时,将系统当前状态设置为完成态,
其中,当系统当前状态为下降态、上升态或平稳态时,转入对当前数据的下一数据的判断,直至当前状态为完成态或错误态。
在本发明的另一实施例中,所述步骤(4)具体为:
对所述有效数据进行数学期望计算,获得数学期望值;
选择分位数,对所述有效数据进行分位数计算,获得分位数值;
将所述数学期望值和所述分位数值相加平均后得到波峰质心值;
取在所述离散波形数据中位于所述波峰质心值左右两侧且与波峰质心值最接近的两个采样点的横坐标值;
取所述两个采样点的横坐标值的中间位置为波峰位置。
在本发明的又一实施例中,所述步骤(2)还包括:
在切换系统当前状态过程中,当当前数据大于阀值且前一数据小于阀值时,将波形开始位置更新为当前数据的横坐标;
在切换系统当前状态过程中,当当前数据小于阀值且前一数据大于阀值时,将波形结束位置更新为当前数据的横坐标。
在本发明的再一实施例中,所述方法还包括:
当波峰位置值不在所述起始位置值与所述结束位置值之间时,根据当前数据的下一数据判断系统当前状态,直至确定的波峰位置值在起始位置值与结束位置值之间。
本发明还提供了一种波峰寻找系统,包括:
波形数据采集模块,用于对光栅反射波形信号进行高速A/D采样,获取一组离散波形数据;
状态转换模块,用于根据所述波形数据采集模块获取的离散波形数据的当前数据和前一数据实时切换系统当前状态,直至系统当前状态为完成态或错误态,状态转换过程中获取离散波形数据中的有效数据;
错误态处理模块,用于当所述状态转换模块转换后的系统当前状态为错误态时,放弃当前采样的离散波形数据;
完成态处理模块,对所述有效数据进行数学期望和分位数计算,获得波峰位置值。
在本发明的一个实施例中,所述状态转换模块包括:
初始态确定子模块,用于当所述波形数据采集模块获取的离散波形数据的当前采样值不大于阀值时,前一状态是初始态,将系统当前状态设置为是初始态;
下降态确定子模块,用于当所述波形数据采集模块获取的离散波形数据的当前数据和前一数据均大于采样阀值,且当前数据小于前一数据、前一数据不等于采样阀值时,将系统当前状态设置为下降态;
错误态确定子模块,用于当所述波形数据采集模块获取的离散波形数据的当前数据和前一数据均大于阀值,前一状态是下降态,且当前数据大于前一数据或当前数据等于前一数据时,将当前状态设置为错误态;
上升态确定子模块,用于当所述波形数据采集模块获取的离散波形数据的当前数据和前一数据均大于阀值,当前数据大于前一数据,且前一状态不是下降态时,将当前状态设置为上升态;
平稳态确定子模块,用于当所述波形数据采集模块获取的离散波形数据的当前数据和前一数据均大于阀值,当前数据等于前一数据,且前一状态不是下降态时,将当前状态设置为平稳态;
完成态确定子模块,用于当所述波形数据采集模块获取的离散波形数据的当前数据小于阀值,前一数据大于阀值,且当前状态不是错误态时,将当前状态设置为完成态,
其中,当所述下降态确定子模块确定的系统当前状态为下降态、或所述上升态确定子模块确定的系统状态为上升态,或所述平稳态确定子模块确定的系统状态为平稳态时,转入对所述波形数据采集模块获取的离散波形数据的当前数据的下一数据的判断,直至系统当前状态为完成态或错误态。
在本发明的另一实施例中,所述完成态处理模块包括:
数学期望值计算子模块,用于对所述状态转换模块确定的有效数据进行数学期望计算,获得数学期望值;
分位数值计算子模块,用于选择分位数,并对所述状态转换模块确定的波形数据进行分位数计算,获得分位数值;
波峰质心值获取子模块,用于将所述数学期望值计算子模块确定的数学期望值和所述分位数值计算子模块确定的分位数值相加平均后得到波峰质心值;
波峰位置确定子模块,用于取在离散波形数据中位于所述波峰质心值获取子模块确定的波峰质心值左右两侧且与波峰质心值最接近的两个采样点的横坐标值;取所述两个采样点的横坐标值的中间位置为波峰位置。
在本发明的又一实施例中,所述系统还包括:
开始位置值更新模块,用于在所述状态转换模块切换系统当前状态过程中,当当前数据大于阀值且前一数据小于阀值时,将波形开始位置更新为当前数据的横坐标;
结束位置值更新模块,用于在所述状态转换模块切换系统当前状态过程中,当当前数据小于阀值且前一数据大于阀值时,将波形结束位置更新为当前数据的横坐标。
在本发明的再一实施例中,所述系统还包括:
波峰位置值审核模块,用于判断所述完成态处理模块确定的波峰位置值是否在所述开始位置值更新模块确定的起始位置值与所述结束位置值更新模块确定的结束位置值之间,
其中,当波峰位置值审核模块判断出波峰位置值不在起始位置值与结束位置值之间时,所述状态转换模块根据当前数据的下一数据判断系统当前状态,直至波峰位置值审核模块判断所述完成态处理模块确定的波峰位置值在起始位置值与结束位置值之间。
与现有技术相比,本发明波峰寻找方法及系统具有如下优点:
1)状态转换模块将采样数据导入状态转换表,通过系统状态来实时跟踪并确定有效波形数据,通过系统状态(初始态、上升态、平稳态、下降态、完成态、错误态)的转换,一方面使系统始终与当前状态同步,因此无需大量存储空间和众多指令周期,能够从采样的离散数据中快速确定有效波形数据,降低了处理时间且使判断结果具有很高的可靠性。
2)通过数学期望值计算子模块和分位数值计算子模块来计算数学期望和分位数对波形数据进行整体处理,采用波峰质心值获取子模块和波峰位置确定子模块来获得波形位置,降低了干扰因素对波形的局部影响,具有很强的抗干扰能力和稳定性,能够保证寻峰结果的稳定性和重复性,适合在高速高精光栅解调仪上应用。
通过以下的描述并结合附图,本发明将变得更加清晰,这些附图用于解释本发明的实施例。
附图说明
图1为本发明波峰寻找方法的流程图。
图1a为图1所述波峰寻找方法中系统当前状态转换的流程图。
图1b为图1所述波峰寻找方法中确定波峰位置的流程图。
图1c为图1a所示当前状态转换的示意图。
图2为本发明波峰寻找系统的结构框图。
图2a为图2所述波峰寻找系统的状态转换模块的结构框图。
图2b为图2所述波峰寻找系统的完成态处理模块的结构框图。
具体实施方式
现在参考附图描述本发明的实施例,附图中类似的元件标号代表类似的元件。
图1为本发明波峰寻找方法的流程图。如图1所示,所述峰寻找方法包括如下步骤:
步骤S10:对光栅反射波形信号进行高速A/D采样,获取一组代表波形信号的离散波形数据;
步骤S20:判断所述离散波形数据中80%以上的数据是否高于采样阀值(所述采样阀值可以根据实际环境动态地设置,阀值的设置方法是:光源发出的光入射光栅传感器,光栅传感器会反射一个中心波长为布拉格波长的波峰信号,如果光栅传感器所测的物理量发生了变化,则该波峰信号会相对于其中心波长发生位移,寻峰算法就是要获取波峰的位移量。波峰的光信号经过光电转换后成为连续的模拟电信号,对该模拟电信号进行前端采样和放大后,根据所得到的电信号的信噪比的大小,由软件动态计算得到采样阀值。大于采样阀值的数据认为是有效数据,小于采样阀值的数据认为是无效数据。所述判断步骤是为了过滤掉系统由于电源,布线等产生的系统干扰),如果否(即当前数据可以认为是来自外界的干扰,丢弃当前采样数据,重新进行采样),转步骤S10,如果是,继续下一步;
步骤S30:将所述离散波形数据导入状态转换表,根据所述离散波形数据的当前数据和前一数据实时切换系统当前状态,直至系统当前状态为完成态或错误态,状态转换过程中获取起始位置值、结束位置值、有效数据、无效数据;
步骤S40:判断系统的当前状态为错误态还是完成态,当系统的当前状态为错误态时(说明先前系统采集到了大于采样阀值的数据,由于种种原因波峰形状出现了畸变,此时波形确定失败),转步骤S60,当系统的当前状态为完成态时(波形确定成功),转步骤S60;
步骤S50:放弃当前采样的离散波形数据,转步骤S10;
步骤S60:根据所述有效数据确定波峰位置值;
步骤S70:判断所述波峰位置值是否在所述起始位置值与所述结束位置值之间,如果是,继续下一步,如果否,转步骤S10;
步骤S80:系统记录波峰位置值,转入步骤S10。
如图1b,所述步骤S60具体为:
步骤S801,对所述有效数据进行数学期望计算,获得所述有效数据的数学期望值,数学期望计算公式为:
其中,Xi代表第i个有效数据,Pi是有效数据个数的倒数;
步骤S802,选择分位数,对所述有效数据进行分位数计算,获得分位数值(根据实际需要选择分位数,通常为75%分位数,则进行分位数计算式从有效数据中取一个值,使所述有效数据中75%的数值小于等于该值);
步骤S803,将数学期望值和分位数值相加平均后得到波峰质心值Yavr;
步骤S804,取在采集的离散波形数据中位于波峰质心值左右两侧且与“波峰质心值”Yavr最接近的两个采样点(离散采样值)的横坐标值,即位置值Xn和Xm;
步骤S805,取Xn和Xm的中间位置为“波峰位置”Xavr,并存入“波峰位置”寄存器中。
如图1a,所述步骤S30具体为:
步骤S301,将系统当前状态初始化为初始态;
步骤S302,判断对离散波形数据的所有数据是否都处理完毕,如果是,继续下一步,如果否,转步骤S304;
步骤S303,置“数据无效”标志位,转步骤S10(重新进行数据采样),如果否,继续步骤S400(进行下面的处理);
步骤S304,从离散数据中的当前数据转向下一个数据;
步骤S305,判断当前数据是否大于采样阀值,如果是,继续下一步,如果否,转步骤S317;
步骤S306,判断前一数据是否大于采样阀值,如果是,继续下一步,如果否,转步骤S316;
步骤S307,判断当前数据是否大于前一数据,如果是,继续下一步,如果否,转步骤S614;
步骤S308,判断系统的前一状态是否处于下降态,如果是,继续下一步,如果否,转步骤S310;
步骤S309,当前数据确定为无效数据,将系统当前状态设置为错误态,转步骤S450;
步骤S310,将系统当前状态设置为上升态(此时“数据有效”标志已置位),转步骤S302;
步骤S311,判断当前数据是否等于前一数据,如果是,继续下一步,如果否,转步骤S314;
步骤S312,判断系统的前一状态是否处于下降态,如果是,转步骤S309,如果否,继续下一步;
步骤S313,将系统当前状态设置为平稳态(此时“数据有效”标志已置位),转步骤S302;
步骤S314,判断前一数据是否等于采样阀值,如果是,转步骤S623;如果否,继续下一步;
步骤S315,将系统当前状态设置为下降态(此时“数据有效”标志已置位),转步骤S302;
步骤S316,置位“数据有效”标志,将当前数据的横坐标记录到“开始位置”寄存器中,作为有效波形数据的开始位置,转步骤S304;
步骤S317,判断当前数据的前一数据是否小于采样阀值,如果是,转步骤S302,如果否,继续下一步;
步骤S318,判断系统的前一状态是否为错误态,如果是,转步骤S302,如果否,继续下一步;
步骤S319,当前数据确定为有效数据,将系统当前状态设置为完成态,将当前数据的横坐标记录到“结束位置”寄存器中,作为有效波形数据的结束位置,转步骤S40.
由上可以看出,本实施例通过系统状态来实时跟踪所采集的离散数据,通过系统状态(初始态、上升态、平稳态、下降态、完成态、错误态)的转换,使系统始终与当前状态同步。通过对系统状态变化的判断,能够了解波形存在的畸变情况,由此来决定是否丢弃该次的采样离散数据。
本方法实时记录有效波形数据的开始位置、结束位置。当当前数据大于阀值且前一数据小于阀值时,更新“开始位置”寄存器;当当前数据小于阀值且前一数据大于阀值时,更新“结束位置”寄存器。
图1c为图1a所示系统状态转换示意图。如图1c所示,系统状态可以是:当系统进入数据分析阶段开始处理采集到的离散数据,系统进入初始态。当当前数据大于前一数据且系统处于初始态时,系统从初始态转换到上升态,如图中箭头20所示。当当前数据小于前一数据且系统处于初始态时,系统从初始态转换到下降态,如图中箭头21所示。当当前数据小于前一数据且系统处于上升态时,系统从上身态转换到下降态,并进入数据审核,如图中箭头22所示。当当前数据大于前一数据且系统处于下降态时,系统从下降态转换到上升态,如图中箭头23所示。当当前数据等于前一数据且系统处于上升态时,系统从上升态转换到平稳态,如图中箭头24所示。当当前数据大于前一数据且系统处于平稳态时,系统从平稳态转换到上升态,如图中箭头25所示。当当前数据等于前一数据且系统处于下降态时,系统从下降态转换到平稳态,如图中箭头26所示。当当前数据小于前一数据且系统处于平稳态时,系统从平稳态转换到下降态,如图中箭头27所示。
当系统由上升态进入数据审核,如果当前数据小于阀值,则重新进入上升态;如果当前数据大于阀值,系统前一状态处于上升态或平稳态而当前状态处于上升态,则系统重新进入上升态;如果当前数据大于阀值,系统前一状态处于下降态而当前状态处于上升态,则系统进入错误态,如图中双向箭头32所示。
当系统由平稳态进入数据审核,如果当前数据小于阀值,则重新进入平稳态;如果当前数据大于阀值,系统前一状态为平稳态,当前状态为平稳态,则系统重新进入平稳态;如果当前数据大于阀值,系统前一状态为下降态,当前状态为平稳态,而下一状态为上升态,则系统进入错误态,如图中上向箭头31所示。
当系统由下降态进入数据审核,如果当前数据小于阀值,且前一数据也小于阀值,则重新进入下降态;如果当前数据大于阀值,且前一数据也大于阀值,重新进入下降态;如果当前数据小于阀值,且前一数据大于阀值,则系统进入完成态,如图中双向箭头30所示。
当系统处于错误态,则丢弃当前处理的离散采样数据,开始新一轮采样,系统进入初始态,如图中箭头29所示。
当系统处于完成态,则获取有效波形数据的开始位置和结束位置后开始新一轮采样,系统进入初始态,如图中箭头28所示。
如图2,本发明波峰寻找方法系统,包括:
波形数据采集模块100,用于对光栅反射波形信号进行高速A/D采样,获取一组代表波形信号的离散波形数据;
状态转换模块200,用于根据所述波形数据采集模块100获取的离散波形数据的当前数据和前一数据实时切换系统当前状态,直至系统当前状态为完成态或错误态,状态转换过程中获取有效波形数据的起始位置值和结束位置值、以及离散波形数据中的有效数据、无效数据;
错误态处理模块300,用于当所述状态转换模块200转换后的系统当前状态为错误态时,放弃采样数据,所述波形数据采集模块100重新对光栅反射波形信号进行高速A/D采样,获取另一组代表波形信号的离散波形数据;
完成态处理模块400,用于当所述状态转换模块200状态转换过程中获取的有效数据确定波峰位置。
开始位置值更新模块600,用于在所述状态转换模块200切换系统当前状态过程中,当当前数据大于阀值且前一数据小于阀值时,将波形开始位置更新为当前数据的横坐标;
结束位置值更新模块700,用于在所述状态转换模块200切换系统当前状态过程中,当当前数据小于阀值且前一数据大于阀值时,将波形结束位置更新为当前数据的横坐标。
波峰位置值审核模块800,用于根据所述开始位置值更新模块600确定的起始位置值、所述结束位置值更新模块700确定的结束位置值判断所述完成态处理模块400确定的波峰位置值是否有效,当所述波峰位置值无效时,所述状态转换模块200根据当前数据的下一数据判断系统当前状态,直至所述完成态处理模块400确定的波峰位置值有效。
其中,如图2a,所述状态转换模块200包括:
初始态确定子模块,用于当所述波形数据采集模块获取的离散波形数据的当前数据和前一数据均小于阀值时,当前状态是前一状态是初始态,将当前状态设置为初始态;
下降态确定子模块210,用于当所述波形数据采集模块获取的离散波形数据的当前数据和前一数据均大于阀值,且当前数据小于前一数据时,将当前状态设置为下降态;
错误态确定子模块220,用于当所述波形数据采集模块获取的离散波形数据的当前数据和前一数据均大于阀值,前一状态是下降态,且当前数据大于前一数据或当前数据等于前一数据时,将当前状态设置为错误态;
上升态确定子模块230,用于当所述波形数据采集模块获取的离散波形数据的当前数据和前一数据均大于阀值,当前数据大于前一数据,且当前状态不是下降态时,将当前状态设置为上升态;
平稳态确定子模块240,用于当所述波形数据采集模块获取的离散波形数据的当前数据和前一数据均大于阀值,当前数据等于前一数据,且当前状态不是下降态时,将当前状态设置为平稳态;
完成态确定子模块250,用于当所述波形数据采集模块获取的离散波形数据的当前数据小于阀值,前一数据大于阀值,且当前状态不是错误态时,将当前状态设置为完成态,
其中,当所述下降态确定子模块210确定的系统当前状态为下降态、或所述上升态确定子模块230确定的系统状态为上升态,或所述平稳态确定子模块240确定的系统状态为平稳态时,转入对所述波形数据采集模块100获取的离散波形数据的当前数据的下一数据的判断,直至系统当前状态为完成态或错误态。
其中,如图2b,所述完成态处理模块400包括:
数学期望值计算子模块410,用于对所述状态转换模块200切换系统状态过程中确定的有效数据进行数学期望计算,获得数学期望值;
分位数值计算子模块420,用于选择分位数,并对所述状态转换模块200切换系统状态过程中确定的有效数据进行分位数计算,获得分位数值;
波峰质心值获取子模块430,用于将所述数学期望值计算子模块410确定的数学期望值和所述分位数值计算子模块420确定的分位数值相加平均后得到波峰质心值;
波峰位置确定子模块440,用于取在离散波形数据中位于所述波峰质心值获取子模块430确定的波峰质心值左右两侧且与波峰质心值最接近的两个采样点的横坐标值;取所述两个采样点的横坐标值的中间位置为波峰位置。
以上结合最佳实施例对本发明进行了描述,但本发明并不局限于以上揭示的实施例,而应当涵盖各种根据本发明的本质进行的修改、等效组合。
Claims (6)
1.一种波峰寻找方法,包括如下步骤:
(1)对光栅反射波形信号进行高速A/D采样,获取一组离散波形数据;
(2)根据所述离散波形数据的当前数据和前一数据实时切换系统当前状态,直至系统当前状态为完成态或错误态,状态转换过程中获取离散波形数据中的有效数据;
(3)当系统当前状态为错误态时,放弃当前采样的离散波形数据;
(4)当系统当前状态为完成态时,对所述有效数据进行数学期望和分位数计算,获得波峰位置值,
其特征在于,所述步骤(4)具体为:
对所述有效数据进行数学期望计算,获得数学期望值;
选择分位数,对所述有效数据进行分位数计算,获得分位数值;
将所述数学期望值和所述分位数值相加平均后得到波峰质心值;
取在所述离散波形数据中位于所述波峰质心值左右两侧且与波峰质心值最接近的两个采样点的横坐标值;
取所述两个采样点的横坐标值的中间位置为波峰位置。
2.如权利要求1所述的波峰寻找方法,其特征在于,所述步骤(2)还包括:
在切换系统当前状态过程中,当当前数据大于阀值且前一数据小于阀值时,将波形开始位置值更新为当前数据的横坐标;
在切换系统当前状态过程中,当当前数据小于阀值且前一数据大于阀值时,将波形结束位置值更新为当前数据的横坐标。
3.如权利要求2所述的波峰寻找方法,其特征在于,还包括:
当波峰位置值不在所述开始位置值与所述结束位置值之间时,根据当前数据的下一数据判断系统当前状态,直至确定的波峰位置值在起始位置值与结束位置值之间。
4.一种波峰寻找系统,包括:
波形数据采集模块,用于对光栅反射波形信号进行高速A/D采样,获取一组离散波形数据;
状态转换模块,用于根据所述波形数据采集模块获取的离散波形数据的当前数据和前一数据实时切换系统当前状态,直至系统当前状态为完成态或错误态,状态转换过程中获取离散波形数据中的有效数据;
错误态处理模块,用于当所述状态转换模块转换后的系统当前状态为错误态时,放弃当前采样的离散波形数据;
完成态处理模块,对所述有效数据进行数学期望和分位数计算,获得波峰位置值,其特征在于,
所述完成态处理模块包括:
数学期望值计算子模块,用于对所述状态转换模块确定的有效数据进行数学期望计算,获得数学期望值;
分位数值计算子模块,用于选择分位数,并对所述状态转换模块确定的波形数据进行分位数计算,获得分位数值;
波峰质心值获取子模块,用于将所述数学期望值计算子模块确定的数学期望值和所述分位数值计算子模块确定的分位数值相加平均后得到波峰质心值;
波峰位置确定子模块,用于取在离散波形数据中位于所述波峰质心值获取子模块确定的波峰质心值左右两侧且与波峰质心值最接近的两个采样点的横坐标值;取所述两个采样点的横坐标值的中间位置为波峰位置。
5.如权利要求4所述的波峰寻找系统,其特征在于,还包括:
开始位置值更新模块,用于在所述状态转换模块切换系统当前状态过程中,当当前数据大于阀值且前一数据小于阀值时,将波形开始位置值更新为当前数据的横坐标;
结束位置值更新模块,用于在所述状态转换模块切换系统当前状态过程中,当当前数据小于阀值且前一数据大于阀值时,将波形结束位置值更新为当前数据的横坐标。
6.如权利要求5所述的波峰寻找系统,其特征在于,还包括:
波峰位置值审核模块,用于判断所述完成态处理模块确定的波峰位置值是否在所述开始位置值更新模块确定的开始位置值与所述结束位置值更新模块确定的结束位置值之间,
其中,当波峰位置值审核模块判断出波峰位置值不在开始位置值与结束位置值之间时,所述状态转换模块根据当前数据的下一数据判断系统当前状态,直至波峰位置值审核模块判断所述完成态处理模块确定的波峰位置值在开始位置值与结束位置值之间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100005764A CN102098101B (zh) | 2011-01-04 | 2011-01-04 | 波峰寻找方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100005764A CN102098101B (zh) | 2011-01-04 | 2011-01-04 | 波峰寻找方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102098101A CN102098101A (zh) | 2011-06-15 |
CN102098101B true CN102098101B (zh) | 2013-05-08 |
Family
ID=44130986
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011100005764A Expired - Fee Related CN102098101B (zh) | 2011-01-04 | 2011-01-04 | 波峰寻找方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102098101B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102435214B (zh) * | 2011-09-20 | 2014-02-12 | 华北电力大学(保定) | 光纤光栅谱峰高精度快速寻峰方法 |
CN102706372B (zh) * | 2012-03-28 | 2014-07-16 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种光纤波长解调光谱峰值定位方法 |
CN102708394B (zh) * | 2012-04-17 | 2016-02-17 | 重庆大学 | 基于saw的无源温度标签及其阅读器 |
CN103822646B (zh) * | 2014-03-03 | 2016-03-30 | 武汉理工大学 | 光纤光栅解调系统中长距离光延迟效应的解调方法 |
CN103913765B (zh) * | 2014-03-24 | 2016-08-31 | 中国船舶重工集团公司第七一九研究所 | 一种核素能谱寻峰方法 |
CN104864894B (zh) * | 2015-04-16 | 2017-10-31 | 中国电子科技集团公司第四十一研究所 | 一种botdr布里渊散射谱寻峰方法 |
CN107560645A (zh) * | 2017-08-29 | 2018-01-09 | 北京航空航天大学 | 一种光纤布拉格光栅传感器波长解调寻峰方法 |
CN109087476A (zh) * | 2018-02-06 | 2018-12-25 | 上海光栅信息技术有限公司 | 一种光纤光栅线型感温火灾探测器解调程序运行方法 |
CN110260898B (zh) * | 2019-06-24 | 2021-07-06 | 武汉理工光科股份有限公司 | 抗干扰的光栅波长解调方法及系统 |
CN113283413B (zh) * | 2021-07-26 | 2021-09-24 | 枫树谷(成都)科技有限责任公司 | 脉冲波形模板库的创建方法、系统、存储介质及设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1614359A (zh) * | 2004-12-07 | 2005-05-11 | 天津大学 | 实现多通道光纤光栅传感装置高灵敏度测量的方法 |
CN101241029A (zh) * | 2007-02-09 | 2008-08-13 | 中国计量学院 | 光纤Bragg光栅传感器解调器 |
CN101692149A (zh) * | 2009-09-30 | 2010-04-07 | 武汉理工大学 | 波峰寻找方法和波峰寻找系统 |
-
2011
- 2011-01-04 CN CN2011100005764A patent/CN102098101B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1614359A (zh) * | 2004-12-07 | 2005-05-11 | 天津大学 | 实现多通道光纤光栅传感装置高灵敏度测量的方法 |
CN101241029A (zh) * | 2007-02-09 | 2008-08-13 | 中国计量学院 | 光纤Bragg光栅传感器解调器 |
CN101692149A (zh) * | 2009-09-30 | 2010-04-07 | 武汉理工大学 | 波峰寻找方法和波峰寻找系统 |
Also Published As
Publication number | Publication date |
---|---|
CN102098101A (zh) | 2011-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102098101B (zh) | 波峰寻找方法及系统 | |
CN101604287A (zh) | 一种基于硬件计数器获取性能数据实现动态优化服务器性能的方法 | |
CN101701869B (zh) | 一种适于控制锥齿轮啮合质量的便携式传动噪声测试仪 | |
WO2023123970A1 (zh) | 光纤传感数据的全局相位正交解调方法及装置 | |
CN108982941A (zh) | 一种周期频率的实时测量方法及装置、示波器、存储介质 | |
CN102761360A (zh) | 一种光信号处理方法及装置 | |
CN106526380B (zh) | 电磁发射要素检测分析系统 | |
US20170363671A1 (en) | Noise spectrum analysis for electronic device | |
CN116706884A (zh) | 光伏发电量的预测方法、装置、终端及存储介质 | |
CN109446572B (zh) | 一种量子测控数据的处理方法及装置 | |
CN106197803A (zh) | 跌落数据的采集方法及终端设备 | |
CN111858703B (zh) | 一种benchmarkSQL性能测试数据统计方法及装置 | |
CN113223499B (zh) | 一种音频负样本的生成方法及装置 | |
FR2758639A1 (fr) | Procede pour transmettre les donnees de pointage d'une commande a distance et procede pour traiter les donnees recues | |
CN110991641B (zh) | 一种油藏类型分析方法、装置及电子设备 | |
CN115222018A (zh) | 非侵入式负荷分解方法及装置 | |
CN111221821A (zh) | 一种ai模型迭代更新方法、电子设备及存储介质 | |
CN201017228Y (zh) | 数控机床加工动力学特性测试分析系统 | |
CN102236261B (zh) | 一种基于正交化模型的离轴信号处理方法、装置及系统 | |
Kneller | Measurement, control and other processes: to the problem of knowledge systematization | |
CN117349661B (zh) | 柱塞泵振动信号特征提取方法、装置、设备和存储介质 | |
CN116297883B (zh) | 一种基于敲击声的结构识别方法、装置、系统及终端设备 | |
CN113407435B (zh) | 一种测试用例的生成方法、装置及控制设备 | |
CN101806820B (zh) | 软件自动调零方法 | |
CN111198789B (zh) | 一种fft硬件实现模块的验证方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130508 Termination date: 20210104 |