CN102435214B - 光纤光栅谱峰高精度快速寻峰方法 - Google Patents
光纤光栅谱峰高精度快速寻峰方法 Download PDFInfo
- Publication number
- CN102435214B CN102435214B CN201110278661.7A CN201110278661A CN102435214B CN 102435214 B CN102435214 B CN 102435214B CN 201110278661 A CN201110278661 A CN 201110278661A CN 102435214 B CN102435214 B CN 102435214B
- Authority
- CN
- China
- Prior art keywords
- peak
- fiber grating
- threshold
- crest
- point
- 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
- Optical Transform (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种光纤光栅谱峰高精度快速寻峰方法,方法采用了相对阈值和绝对阈值两个参数,寻峰时采用的有效阈值等于反射谱数据的最大值减去相对阈值;当通过相对阈值设定的有效阈值小于绝对阈值时,使有效阈值=绝对阈值;判断波峰是光纤光栅反射峰还是噪声形成的峰时,利用有效阈值分别减、加H得到两个参数,通过比较反射峰两侧两个参数之间数据的关系,判断是否存在上升沿和下降沿,如果上升沿和下降沿都存在,则说明此处为一个波峰,再认定该峰为光纤光栅的反射峰还是噪声形成的波峰,最终进一步得到峰值所在位置。该方法解决了现有方法中寻峰精度低、运算量大的问题。
Description
技术领域
本发明涉及光纤光栅传感领域,尤其涉及一种光纤光栅谱峰高精度快速寻峰方法。
背景技术
随着光纤传感技术的发展,光纤光栅传感器以其抗电磁干扰、耐腐蚀、高绝缘性、测量范围广、便于复用成网等优点,在土木工程、航空航天、石油化工、电力、医疗、船舶工业等领域获得了广泛应用。当外界作用在光纤光栅传感器上的温度或应变发生变化时,光纤光栅的反射谱或透射谱中心波长会发生变化,因此,寻找光纤光栅谱峰峰值的精确位置是光纤光栅传感中的一项关键技术,快速、准确地寻峰方法是保证光纤光栅传感系统实时性和精度的必要条件。目前,常用的寻峰算法有蒙特卡罗算法、直接比较法、二次插值数值微分法、一般多项式拟合法、多项式-高斯公式拟合法等,利用上述方法寻峰,精度低、运算量大,在对传感精度要求比较高时,不能满足要求,另外,某些应用场合需要用微处理器替代PC机实现寻峰,采用上述方法时会严重影响整个传感系统的实时性。因此,研究光纤光栅谱峰高精度快速寻峰方法是测量工作的需要。
发明内容
本发明的目的在于提供一种应用于光纤光栅传感系统中的光纤光栅谱峰高精度快速寻峰方法,实现光纤光栅谱峰峰值位置的精确定位,解决现有寻峰方法中寻峰精度低、运算量大的问题。
本发明通过以下方法与步骤来实现:1)首先根据光纤光栅反射谱多个峰值中的最大值及信噪比,确定寻峰时采用的阈值;2)根据设定的反射谱的宽度和形状特点,判断波峰是光纤光栅反射峰还是噪声形成的峰,并记录光纤光栅反射峰的起始点和结束点的横坐标值;3)找到光纤光栅反射峰两侧比最大值低3dB的左右两个点,然后将两个点对应的横坐标数值求和后除2,得到峰值所在位置;所述确定寻峰时采用的阈值,采用了相对阈值relThreshold和绝对阈值absThreshold两个参数,寻峰时采用的有效阈值threshold等于反射谱数据的最大值减relThreshold;当通过relThreshold设定的有效阈值小于absThreshold时,使有效阈值threshold=absThreshold;所述判断波峰是光纤光栅反射峰还是噪声形成的峰时,利用有效阈值threshold分别减、加迟滞阈值H得到两个参数thresholdLow和thresholdHigh;通过比较反射峰左侧thresholdLow和thresholdHigh之间数据的关系,判断是否存在上升沿,如果反射峰左侧存在上升沿,再采用相同方法判断反射峰右侧thresholdHigh和thresholdLow之间是否存在下降沿;如果上升沿和下降沿都存在,则说明此处为一个波峰,进一步判断当end-start大于波峰宽度width时,则认定该峰为光纤光栅的反射峰,否则,则判定该峰为噪声形成的波峰。
在寻找光纤光栅反射峰两侧比最大值低3dB的左右两个点时,定义左侧的点的坐标为(left3dBPos,left3dB),右侧的点的坐标为(right3dBPos,right3dB),点(left3dBPos,left3dB)两侧的采样点的坐标为(x1-1,f(x1-1))和(x1,f(x1)),点(right3dBPos,right3dB)两侧的采样点的坐标为(x2,f(x2))和(x2+1,f(x2+1)),峰值点的横坐标为peakPos,将光纤光栅反射峰视为对称峰,利用三角形的相似性得到了如下关系式:
可得:
同理可得:
由此推出峰值点的准确位置表达式:
当寻找透射谱中的负峰时,第一步中,确定寻峰时采用的阈值时,采用的有效阈值threshold等于透射谱数据的最小值加上relThreshold,当通过relThreshold设定的有效阈值大于absThreshold时,使有效阈值threshold=absThreshold;第二步中,判断是否为光纤光栅形成的波峰时,应判断波峰左侧thresholdHigh和thresholdLow之间是否存在下降沿,波峰右侧thresholdLow和thresholdHigh之间是否存在上升沿;第三步中,应找到光纤光栅透射峰两侧比最小值高3dB的左右两个点,然后将两个点对应的横坐标数值求和后除2,得到峰值所在位置。
附图说明
图1本发明相对阈值设定图例
图2本发明绝对阈值设定图例
图3本发明峰值区域确定图例
图4本发明峰值位置确定图例
具体实施方式
利用光纤光栅的反射谱和透射谱都可实现温度和应变的传感,反射谱的谱峰为正峰,透射谱的谱峰为负峰,以下结合附图和实施例,以寻找反射谱中的正峰位置为例,对本发明进一步详细说明,本方法主要分三步:
1)首先根据光纤光栅反射谱多个峰值中的最大值及信噪比,确定寻峰时采用的阈值。从反射谱数据的第一数据开始逐个与阈值进行比较,当反射谱数据小于该阈值时,认为当前数据为噪声,当反射谱数据大于该阈值时,初步判断当前数据为光栅反射峰,然后进入下一步进一步判断是旁瓣还是真峰。寻峰时采用的阈值,采用了相对阈值relThreshold和绝对阈值absThreshold两个参数,寻峰时采用的有效阈值threshold等于反射谱数据的最大值减去relThreshold;当通过relThreshold设定的有效阈值小于absThreshold时,使有效阈值threshold=absThreshold。图1中反射谱的最大功率值为-20dBm,当relThreshold设置为7dB时,寻峰时采用的有效阈值threshold为-20dBm-7dB=-27dBm。
absThreshold用于设定寻峰算法采用的最小有效阈值,如图2中所示,由于系统光功率衰减过大等原因,使光纤光栅反射峰的最大值只有-32dBm,如果按照relThreshold=7dB时得到有效阈值threshold=-39dBm进行寻峰,寻到的全部是噪声信号,根据图2中噪声情况,可设置absThreshold=-35dBm,此时,可知absThreshold>-32dBm-relThreshold,则设置有效阈值threshold=absThreshold=-35dBm。
2)根据设定的反射谱的宽度和形状特点,判断波峰是光纤光栅反射峰还是噪声形成的峰,并记录光纤光栅反射峰的起始点和结束点的横坐标值。判断波峰是光纤光栅反射峰还是噪声形成的峰时,利用有效阈值threshold分别减、加迟滞阈值H得到两个参数thresholdLow和thresholdHigh。通过比较反射峰左侧thresholdLow和thresholdHigh之间数据的关系,判断是否存在上升沿,如果反射峰存在上升沿,再采用相同方法判断反射峰右侧thresholdHigh和thresholdLow之间是否存在下降沿。如果上升沿和下降沿都存在,则说明此处为一个波峰,进一步判断当end-start大于波峰宽度width时,则认定该峰为光纤光栅的反射峰,否则,则判定该峰为噪声形成的波峰。图3给出了实例,本实施例中H取3dB,所以,当threshold=-26dBm时,thresholdLow=-29dBm,thresholdHigh=-23dBm。
3)找到光纤光栅反射峰两侧比最大值低3dB的左右两个点,然后将两个点对应的横坐标数值求和后除2,得到峰值所在位置。寻找比光纤光栅反射峰最大值低3dB的左右两个点时,将光纤光栅反射峰视为对称峰,由图4可见,定义左侧的点的坐标为(left3dBPos,left3dB),右侧的点的坐标为(right3dBPos,right3dB),点(left3dBPos,left3dB)两侧的采样点的坐标为(x1-1,f(x1-1))和(x1,f(x1)),点(right3dBPos,right3dB)两侧的采样点的坐标为(x2,f(x2))和(x2+1,f(x2+1)),峰值点的横坐标为peakPos,利用三角形的相似性得到了如下关系式:
由式(1)可得:
同理可得:
找到leff3dB和right3dB点后,可得出峰值点的准确位置:
当利用光纤光栅的透射谱时,应将本实施例中寻找正峰的准确位置改为寻找负峰的准确位置。寻找负峰时,第一步中,应根据多个峰值中的最小值及信噪比确定阈值,有效阈值threshold等于反射谱数据的最小值加上relThreshold,当通过relThreshold设定的有效阈值大于absThreshold时,使有效阈值threshold=absThreshold;第二步中,判断是否为光纤光栅形成的波峰时,应判断波峰左侧thresholdHigh和thresholdLow之间是否存在下降沿,波峰右侧thresholdLow和thresholdHigh之间是否存在上升沿;第三步中,应找到光纤光栅透射峰两侧比最小值高3dB的左右两个点,然后将两个点对应的横坐标数值求和后除2,得到负峰值所在位置。
上述寻峰方法是以反射型光纤光栅为例说明的,相同方法做适当变化也适合于透射型光纤光栅。对于透射型光纤光栅,反射谱表现为负峰,应按上述寻找反射型光纤光栅透射谱负峰的步骤确定峰值的位置;透射谱表现为正峰,应按上述寻找反射型光纤光栅反射谱正峰的步骤确定峰值的位置。
综上所述,本发明实施例通过确定阈值、确定谱峰所在区域和精确查找谱峰位置三步,可精确定位光纤光栅传感系统中光栅谱峰的位置,从而可进一步实现高精度的温度、应变及其他物理量的传感。
Claims (3)
1.一种光纤光栅谱峰高精度快速寻峰方法,其特征在于,寻找反射谱中的正峰时,所述方法包括以下步骤:1)根据光纤光栅反射谱多个峰值中的最大值及信噪比,确定寻峰时采用的阈值;2)根据设定的反射谱的宽度和形状特点,判断波峰是光纤光栅反射峰还是噪声形成的峰,并记录光纤光栅反射峰的起始点和结束点的横坐标值;3)找到光纤光栅反射峰两侧比最大值低3dB的左右两个点,然后将两个点对应的横坐标数值求和后除2,得到峰值所在位置;所述确定寻峰时采用的阈值,采用了相对阈值relThreshold和绝对阈值absThreshold两个参数,寻峰时采用的有效阈值threshold等于反射谱数据的最大值减relThreshold;当通过relThreshold设定的有效阈值小于absThreshold时,使有效阈值threshold=absThreshold;所述判断波峰是光纤光栅反射峰还是噪声形成的峰时,利用有效阈值threshold分别减、加迟滞阈值H得到两个参数thresholdLow和thresholdHigh;通过比较反射峰左侧thresholdLow和thresholdHigh之间数据的关系,判断是否存在上升沿,如果反射峰左侧存在上升沿,再采用相同方法判断反射峰右侧thresholdHigh和thresholdLow之间是否存在下降沿;如果上升沿和下降沿都存在,则说明此处为一个波峰,进一步判断当end-start大于波峰宽度width时,则认定该峰为光纤光栅的反射峰,否则,则判定该峰为噪声形成的波峰。
2.根据权利要求1所述的方法,其特征在于,所述寻找光纤光栅反射峰两侧比最大值低3dB的左右两个点时,定义左侧的点的坐标为(left3dBPos,left3dB),右侧的点的坐标为(right3dBPos,right3dB),点(left3dBPos,left3dB)两侧的采样点的坐标为(x1-1,f(x1-1))和(x1,f(x1)),点(right3dBPos,right3dB)两侧的采样点的坐标为(x2,f(x2))和(x2+1,f(x2+1)),峰值点的横坐标为peakPos将光纤光栅反射峰视为对称峰,利用三角形的相似性得到了如下关系式:
可得:
同理可得:
由此推出峰值点的准确位置表达式:
3.根据权利要求1所述的方法,其特征在于:当寻找透射谱中的负峰时,第一步中,确定寻峰时采用的阈值时,采用的有效阈值threshold等于透射谱数据的最小值加上relThreshold,当通过relThreshold设定的有效阈值大于absThreshold时,使有效阈值threshold=absThreshold;第二步中,判断是否为光纤光栅形成的波峰时,应判断波峰左侧thresholdHigh和thresholdLow之间是否存在下降沿,波峰右侧thresholdLow和thresholdHigh之间是否存在上升沿;第三步中,应找到光纤光栅透射峰两侧比最小值高3dB的左右两个点,然后将两个点对应的横坐标数值求和后除2,得到峰值所在位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110278661.7A CN102435214B (zh) | 2011-09-20 | 2011-09-20 | 光纤光栅谱峰高精度快速寻峰方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110278661.7A CN102435214B (zh) | 2011-09-20 | 2011-09-20 | 光纤光栅谱峰高精度快速寻峰方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102435214A CN102435214A (zh) | 2012-05-02 |
CN102435214B true CN102435214B (zh) | 2014-02-12 |
Family
ID=45983412
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110278661.7A Expired - Fee Related CN102435214B (zh) | 2011-09-20 | 2011-09-20 | 光纤光栅谱峰高精度快速寻峰方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102435214B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103644974B (zh) * | 2013-12-11 | 2016-05-25 | 李智忠 | 一种光纤布拉格光栅反射谱中心波长位移的计测方法 |
CN103913765B (zh) * | 2014-03-24 | 2016-08-31 | 中国船舶重工集团公司第七一九研究所 | 一种核素能谱寻峰方法 |
CN104111082A (zh) * | 2014-07-11 | 2014-10-22 | 中南大学 | 一种高精度fbg传感信号寻峰方法 |
CN104864894B (zh) * | 2015-04-16 | 2017-10-31 | 中国电子科技集团公司第四十一研究所 | 一种botdr布里渊散射谱寻峰方法 |
CN105241592A (zh) * | 2015-09-12 | 2016-01-13 | 承德石油高等专科学校 | 基于多通道dwdm技术的新型波长解调器 |
CN105841723B (zh) * | 2015-10-13 | 2018-04-10 | 北京信息科技大学 | 基于小波奇异性检测的fbg传感解调寻峰方法 |
CN109141827B (zh) * | 2017-06-28 | 2020-07-07 | 福州高意通讯有限公司 | 一种应用于可调谐滤波器光性能检测装置的光谱寻峰方法 |
CN107796421B (zh) * | 2017-10-13 | 2018-08-21 | 北京卫星环境工程研究所 | 快速确定fbg反射光谱异常模式的方法 |
CN108256238B (zh) * | 2018-01-22 | 2021-08-03 | 武汉理工大学 | 一种基于深度学习的光纤光栅波长解调方法及装置 |
CN114047160B (zh) * | 2021-10-28 | 2023-06-13 | 中南大学 | 二次谐波穿线寻峰方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6002472A (en) * | 1997-03-18 | 1999-12-14 | The Furukawa Electric Co., Ltd. | Method of measuring optical fiber drawing tension |
US6052179A (en) * | 1998-04-14 | 2000-04-18 | Physical Optics Corporation | Method and system for determining the wavelength of light transmitted through an optical fiber |
CN101692149A (zh) * | 2009-09-30 | 2010-04-07 | 武汉理工大学 | 波峰寻找方法和波峰寻找系统 |
CN102098101A (zh) * | 2011-01-04 | 2011-06-15 | 武汉理工大学 | 波峰寻找方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001201407A (ja) * | 2000-01-19 | 2001-07-27 | Fuji Electric Co Ltd | 波長測定装置 |
-
2011
- 2011-09-20 CN CN201110278661.7A patent/CN102435214B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6002472A (en) * | 1997-03-18 | 1999-12-14 | The Furukawa Electric Co., Ltd. | Method of measuring optical fiber drawing tension |
US6052179A (en) * | 1998-04-14 | 2000-04-18 | Physical Optics Corporation | Method and system for determining the wavelength of light transmitted through an optical fiber |
CN101692149A (zh) * | 2009-09-30 | 2010-04-07 | 武汉理工大学 | 波峰寻找方法和波峰寻找系统 |
CN102098101A (zh) * | 2011-01-04 | 2011-06-15 | 武汉理工大学 | 波峰寻找方法及系统 |
Non-Patent Citations (7)
Title |
---|
《实验中采集CCD光谱测量原理的初步研究》;旺堆等;《西藏大学学报》;20091031;第24卷(第2期);第122页左栏3.2节第1段、倒数第1段 * |
《提高光纤Bragg光栅波长测量精度的方法》;付建伟等;《光电工程》;20060228;第33卷(第2期);第124页第1节倒数第1段 * |
《高速高精光栅解调器的反馈式相对寻峰算法》;韩屏等;《武汉理工大学学报》;20110228;第33卷(第1期);第10页右栏第2.1节 * |
JP特开2001-201407A 2001.07.27 |
付建伟等.《提高光纤Bragg光栅波长测量精度的方法》.《光电工程》.2006,第33卷(第2期), |
旺堆等.《实验中采集CCD光谱测量原理的初步研究》.《西藏大学学报》.2009,第24卷(第2期), |
韩屏等.《高速高精光栅解调器的反馈式相对寻峰算法》.《武汉理工大学学报》.2011,第33卷(第1期), |
Also Published As
Publication number | Publication date |
---|---|
CN102435214A (zh) | 2012-05-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102435214B (zh) | 光纤光栅谱峰高精度快速寻峰方法 | |
CN103487074B (zh) | 利用三点寻峰算法处理fbg传感信号的方法 | |
CN107430013B (zh) | 使用参考光纤干涉数据配准测量的光纤干涉数据 | |
US20150066390A1 (en) | Error measuring method of gear | |
US10775245B2 (en) | Temperature measurement device, temperature measurement method, and computer-readable non-transitory medium | |
AU2019202787B2 (en) | Temperature measurement device, temperature measurement method, and temperature measurement program | |
CN107036720B (zh) | 一种超短啁啾脉冲时域相位与频域相位测量方法及系统 | |
CN105259548A (zh) | 一种用于fmcw绝对距离测量技术中色散失配校正方法 | |
KR20090025959A (ko) | 거리 측정 방법 및 장치 | |
CN103017887A (zh) | 一种光纤振动传感系统及其检测方法 | |
Ri et al. | Immediate estimation of feedback factor and linewidth enhancement factor from measured self-mixing signals under moderate or strong regime | |
CN105203209A (zh) | 基于光强测量积分时间优化的Stokes矢量测量系统及方法 | |
US10054516B2 (en) | System and method for optical frequency domain reflectometer | |
CN102927925A (zh) | 基于自适应加权融合算法的倾斜光纤光栅应变测量方法 | |
Liu et al. | ESMD-based stability analysis in the progressive collapse of a building model: A case study of a reinforced concrete frame-shear wall model | |
CN108196242A (zh) | 基于边沿检测的激光雷达计时方法及数据处理单元 | |
US10247622B2 (en) | Temperature measurement device and temperature measurement method | |
CN102944696B (zh) | 啁啾光纤光栅传感器及处理方法 | |
CN113343173B (zh) | 一种布里渊频移提取方法 | |
AU2015394728A1 (en) | Temperature measurement device, temperature measurement method, and temperature measurement program | |
CN110471075A (zh) | 雷达测距方法、装置及终端设备 | |
CN106154257A (zh) | 基于FFT与apFFT的精密测量雷达二次测频方法 | |
EP2947446B1 (en) | Method for measuring light physical constants and device for estimating light physical constants | |
CN101710023B (zh) | 一种保偏光纤拍长的测试方法及装置 | |
CN110635841B (zh) | 一种提高混沌光时域反射仪回波信号的方法及装置 |
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 |
Granted publication date: 20140212 Termination date: 20170920 |
|
CF01 | Termination of patent right due to non-payment of annual fee |