CN108896078A - 基于探测器时域响应的光纤布拉格光栅弱信号解调方法 - Google Patents
基于探测器时域响应的光纤布拉格光栅弱信号解调方法 Download PDFInfo
- Publication number
- CN108896078A CN108896078A CN201810302314.5A CN201810302314A CN108896078A CN 108896078 A CN108896078 A CN 108896078A CN 201810302314 A CN201810302314 A CN 201810302314A CN 108896078 A CN108896078 A CN 108896078A
- Authority
- CN
- China
- Prior art keywords
- detector
- sequence
- threshold value
- bragg grating
- fiber bragg
- 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.)
- Granted
Links
- 230000004044 response Effects 0.000 title claims abstract description 43
- 239000000835 fiber Substances 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000003595 spectral effect Effects 0.000 claims abstract description 24
- 238000005316 response function Methods 0.000 claims abstract description 23
- 238000000985 reflectance spectrum Methods 0.000 claims abstract description 8
- 238000001228 spectrum Methods 0.000 claims description 15
- 239000013307 optical fiber Substances 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 230000003287 optical effect Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- BJQHLKABXJIVAM-UHFFFAOYSA-N bis(2-ethylhexyl) phthalate Chemical compound CCCCC(CC)COC(=O)C1=CC=CC=C1C(=O)OCC(CC)CCCC BJQHLKABXJIVAM-UHFFFAOYSA-N 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000003321 amplification Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 206010019133 Hangover Diseases 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000009828 non-uniform distribution Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
- G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
- G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
- G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
- G01D5/353—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
- G01D5/35306—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using an interferometer arrangement
- G01D5/35309—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using an interferometer arrangement using multiple waves interferometer
- G01D5/35316—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using an interferometer arrangement using multiple waves interferometer using a Bragg gratings
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Computing Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Optical Transform (AREA)
Abstract
本发明公开了一种基于光电探测器时域响应的光纤布拉格光栅弱信号解调方法，步骤(1)、确定寻峰区域；步骤(2)、获取采样点与波长值对应序列；步骤(3)、构建探测器响应模型步骤(4)、拟合FBG反射光谱：以高斯函数与探测器冲激响应函数的卷积作为拟合函数拟合方程为两函数、四个参数的卷积结果；步骤(5)、获得对应位置解调波长。与现有技术相比，本发明通过卷积函数对谱型的高精度拟合，提高了解调精度；以光电探测器脉冲响应为基础构建的探测器时域响应模型，模型为包含四个参数的方程，因此十分简便；可普遍用于无畸变及不同畸变程度光纤布拉格光栅反射光谱的解调；可实现温度、应变、压力的高速高精度解调。
Description
技术领域
本发明涉及光纤传感技术领域，特别是涉及一种光纤布拉格光栅反射峰中心波长定位方法。
背景技术
光纤布拉格光栅传感器是一种波长调制型传感器，其工作原理基于对外界物理参量如温度、应变、振动、折射率的敏感响应，光纤布拉格光栅传感器反射的中心波长随着物理参量发生变化产生偏移。
目前，光纤布拉格光栅解调技术大多将光纤布拉格光栅反射峰简化为对称的高斯谱，通过寻峰算法检测波长偏移量实现物理量测量。常见的寻峰算法包括：最大值法、高斯拟合法、质心法、互相关法和快速相位相关法等。但由于光纤布拉格光栅反射光谱一方面受自身光学性质如纤芯折射率分布、非线性啁啾的影响；另一方面受所处应用环境如系统调制频率、探测器性能的影响；因此，存在应变场非均匀分布、传感器间多重反射、光谱遮蔽，易出现反射峰非对称的畸变现象，使系统解调精度及准确性降低。
采用扫描光源的解调方法是当前主流的光纤布拉格光栅传感解调方法。该方法中，光源在一个扫描周期内输出光波长随时间连续变化，即每一时刻光源输出一个特定波长，光源输出的光受光纤布拉格光栅反射作用，光纤布拉格光栅的光谱信息将转化为时域脉冲信号进行输出。因此，探测器对脉冲光信号的响应情况，将直接影响光纤布拉格光栅反射峰谱型。以往方法通常工作于较低扫描速率下，无需考虑探测器脉冲响应导致的畸变，大多将光纤布拉格光栅反射峰简化为对称的高斯谱，对发生畸变的谱型实用性较差。
发明内容
本发明的目的在于解决传统解调方法对光纤布拉格光栅畸变谱峰值位置的解调劣势，提供一种基于光电探测器时域响应的光纤布拉格光栅弱信号解调方法，通过探测器脉冲响应与脉冲函数反卷积，得到不同增益下光电探测器的时域响应函数，以高斯函数和探测器响应函数的卷积作为拟合函数，对不同畸变程度的光纤布拉格光栅反射峰进行高精度拟合，进而提高解调精度。
步骤1、构建基于可调谐F-P滤波器的弱光强解调系统，从而获得不同探测器增益下的归一化光纤布拉格光栅反射光谱，即采样点(1,2,…,N)-幅值序列(A1,A2,…,AN)，设定光谱采样点幅值的对应阈值Φ，截取幅值大于阈值的采样点序列P(Aa+1,Aa+2,…,Aa+n)作为寻峰区域；其中，将阈值Φ的选取将遵循高解调精度、低解调时间的准则，设定为0至1范围内的最佳值；
步骤2、以透射峰值波长作为波长参考，获取采样点(1,2,…,N)-波长值对应序列(λ1,λ2,…,λN)；
步骤3、构建探测器响应模型，仿真得到不同畸变程度光纤布拉格光栅反射峰谱型；表达式如下：
其中，符号*表示卷积运算，fin(t)表示时域输入信号，τ为在0至∞范围内变化的时间值，h(t)表示探测器时域冲激响应系数，通过探测器脉冲响应函数与脉冲函数反卷积得到探测器时域冲激响应系数h(t)，表达式如下:
h(t)＝α·exp(-t/β)
其中，α表示与探测器输出幅值相关的系数，β表示探测器时间响应系数的拟合值；
步骤4、拟合FBG反射光谱：以高斯函数与探测器冲激响应函数的卷积作为拟合函数
其中，高斯函数包含对称轴μ、标准差σ两个参数，探测器响应函数包含幅值参数α、时间系数β两个参数；拟合方程为两函数、四个参数的卷积结果；式中，N表示一个周期内等时间间隔采集的采样点总数，m为在1和N之间变化的变量；由于卷积运算使数据量加倍，即由N个点拓展为2N-1个点，因此n为1至2N-1之间的量；
步骤5、设置参数α、β、μ、σ范围及最佳寻峰阈值，以使得最佳寻峰阈值与步骤(1)中的阈值Φ相等；根据阈值Φ，仿真FBG光谱fFBG(n)中在阈值以上的点构成序列Q(Bb+1,Bb+2,…,Bb+m)，序列Q共包含m个点，其中，b+1为序列Q中超过阈值的起始横坐标，b+m为序列Q中超过阈值的终止横坐标；
通过线性插值，将步骤(1)中序列P转换为序列P’(Ab’+1,Ab’+2,…,Ab’+m)，其中，b’+1为变换后的起始横坐标，b’+m为变换后终止横坐标，且Ab’+1＝Aa+1，Ab’+m＝Aa+n，运用最小二乘法：获得寻峰区域拟合最优参数(μ’,σ’,α’,β’)，其中k从1递增至m；求取拟合结果最大值横坐标b’+k0，k0为1至m之间的值；通过三次样条插值，求取b’+k0在序列(λ1,λ2,…,λN)中的对应波长λ，即为解调波长。
相比于传统的解调方法，本发明具有以下的优点和有益效果：
1、本发明将解调系统探测到的光纤布拉格光栅反射信号看成高斯函数与探测器时域冲激响应函数的卷积结果，通过卷积函数对谱型的高精度拟合，提高了解调精度；
2、本发明以光电探测器脉冲响应为基础构建的探测器时域响应模型，模型为包含四个参数的方程，因此十分简便；
3、本发明可普遍用于无畸变及不同畸变程度光纤布拉格光栅反射光谱的解调，适用范围广泛；
4、本发明可以实现温度、应变、压力的高速高精度解调。
附图说明
图1为现有技术的基于可调谐F-P滤波器的光纤布拉格光栅弱光强解调系统装置示意图；
图2为常见的光纤布拉格光栅畸变谱型曲线示意图；
图3为探测器响应函数及拟合方程曲线图；
图4为光纤布拉格光栅反射光谱曲线及寻峰拟合区域；
图5为图4所示谱型的高斯函数、时间响应函数卷积拟合结果；
图6为光源扫描频率为1600Hz，获取探测器增益为20-70dB/10dB步进时光纤布拉格光栅反射峰谱型；
图7为本发明的基于光电探测器时域响应的光纤布拉格光栅弱信号解调方法整体流程示意图；
图中：1、ASE宽带光源，2、可调谐F-P滤波器，3、标准具，4、光纤布拉格光栅传感器序列，5、增益可调光电探测器，6、采集及处理模块。
具体实施方式
下面将结合附图对本发明的实施方式作进一步的详细描述。
本发明的理论依据如下：
在基于扫描光源的解调系统中，FBG的反射光谱实际是时域上周期性的窄带光信号。光信号由光电探测器接收，然后完成光电转换及放大，再传输至采集卡。本发明将完成光电转换的光电探测器及其增益可调的外围放大电路视为一个探测器系统，则系统的输出表示为时域输入信号fin(t)及探测器时域冲激响应系数h(t)的卷积：
其中，符号*表示卷积运算，h(t)表示探测器时域冲激响应系数，冲激响应为幅值趋于无穷大，脉宽趋于零的信号。由于在实际应用中无法获得极限条件的冲激信号，本发明运用系统脉冲响应的概念，求取冲激响应h(t)，其中τ为在0至∞范围内变化的时间值。方法如下：
通常，系统对单位脉冲信号Pn(t)(幅值为n，脉冲宽度为1/n)的响应称为系统单位脉冲响应hn(t)，根据公式(1)，hn(t)表示为：
其中，ω0为探测器时间响应系数，受探测器增益影响，C为与输入信号幅值相关的常系数。则系统时域冲激响应表示为hn(t)与Pn(t)的反卷积结果：
h(t)＝conv-1(hn(t),pn(t)) (3)
值得注意的是，探测器时间响应系数ω0是描述探测器对入射光响应快慢的物理量，该参数数只受探测器内部耗尽层附近光生载流子扩散时间、耗尽层内光生载流子漂移时间以及外部负载电阻与结电容RC时间等控制，与PN结载流子浓度、结构工艺存在复杂的依赖关系。
为得到时域冲激响应函数h(t)，设定参数n,C,τ0代入式(2)，得探测器响应函数选取衰减的指数函数进行拟合，如图3(1)所示为探测器响应函数及拟合方程曲线图。因此，用指数衰减函数表示探测器响应函数α·exp(-t/β)，其中，参数α为与探测器输出幅值相关的系数，参数β为探测器时间响应系数τ0的拟合值。理想情况下，β是趋近零的值，探测器在无弛豫状态下，可做出瞬时响应，得到输入光强的无失真输出，而实际应用中，β是恒大于0的正实数，导致畸变客观存在；探测器的输出会随β的增大而发生愈加明显的畸变。对于内部半导体结构及外部参数确定的探测器，β是一个仅受系统中放大电路增益影响的时间常数。
不考虑由于FBG传感器内部折射率分布不均匀导致的畸变，则输入信号fin(t)为标准高斯轮廓。在每个光源扫描周期内，输入信号fin(t)表示为：
其中，tλ表示高斯轮廓峰值对应时间，σ表示受光脉冲频率影响，决定轮廓的宽度。因此，探测器系统的时域响应模型表示为：
此模型表示FBG反射光谱是FBG反射光脉冲轮廓与探测器系统冲激响应函数的卷积结果。
基于探测器时域响应模型，本发明进而提出一种同时适用于无畸变及不同畸变程度谱型的光纤布拉格光栅拟合寻峰算法，解决传统高斯拟合法、质心法对非对称谱型的寻峰劣势。该算法的核心是利用探测器响应模型作为拟合函数，再进行寻峰。由于解调系统中，采集卡得到的数据是时域上离散的点，因此公式(5)将转化为：
其中，N表示一个周期内等时间间隔采集的采样点总数，m为在1和N之间变化的变量；由于卷积运算使数据量加倍，即由N个点拓展为2N-1个点，因此n为1至2N-1之间的量。寻峰算法中，只需关注峰值附近区域，因此算法采取最佳阈值，仅对峰值附近，幅值大于阈值的区域进行拟合寻峰计算，减少运算量，提高寻峰速度。绘制光纤布拉格光栅反射光谱曲线及选取的寻峰区域如图4所示。最佳寻峰阈值选取的方法，主要涉及不同阈值下的最大解调误差及解调时间的综合考量。设置不同阈值对几种畸变程度的数据进行解调，随着畸变程度增加，最大解调误差在小范围内波动，解调时间减小。基于确保解调精度、缩短解调时间的原则，设置归一化光谱的解调阈值在0.5至0.7范围内。
用一个决定探测器输出谱型畸变程度的时间系数β构建探测器响应函数的方法，主要涉及通过探测器脉冲响应函数与脉冲函数反卷积，求得探测器时域冲激响应函数。通过验证，探测器响应函数为e指数下降形式h(t)＝a·exp(-t/β)，指数中包含一个决定探测器输出谱型畸变程度的参数。该参数与探测器响应速度直接相关，受探测器增益控制。通过构建探测器时域响应模型可仿真得到不同畸变程度光纤布拉格光栅反射峰谱型。
拟合式(6)为包含对称轴μ、标准差σ、幅值参数α和时间系数β四个拟合参数的方程。其中，对称轴μ与光纤布拉格光栅反射峰峰值横坐标直接相关，标准差σ决定高斯函数展宽大小，α决定拟合函数幅值，β决定畸变程度。本算法采用最小二乘法获取拟合参数，再通过插值得到光纤布拉格光栅反射峰峰值横坐标。
如图1所示，是为实现本发明的基于光电探测器时域响应的光纤布拉格光栅弱信号解调方法而构建的基于光纤可调谐F-P滤波器的光纤布拉格光栅弱光强解调系统的实施例。ASE宽带光源1与可调谐F-P滤波器2组合，形成扫描光源，输出C波段(1525-2565nm)连续变化的扫描光。其中可调谐F-P滤波器频率由信号模块的三角波控制，光源输出扫描光频率与三角波相等。扫描光经由耦合器均分为两路，一路传输至标准具3，作为波长参考，另一路连接环形器，控制光沿逆时针方向单向传输：先经过光纤布拉格光栅传感器序列4，传感器分别选取特定波长进行反射，再经过环形器传输至耦合器1：1分为两路，分别连接增益可调的光电探测器5，完成光电转换及信号放大。数据采集卡实现三路信号并行采集，包括标准具信号、两路探测器输出的传感器反射信号，用于处理及解调。
以具体参数实施说明，设置光源扫描频率为1600Hz，获取探测器增益为20-70dB/10dB步进时光纤布拉格光栅反射峰谱型。如图5所示，曲线(1)～(6)分别为探测器增益为20～70dB时的反射谱型。可见，相同光源扫描频率下，探测器增益是影响探测器输出谱型畸变程度的决定参数。当探测器增益大于50dB时，输出光纤布拉格光栅的谱型失去了明显的对称性。
图1中标准具3的透射峰值波长已知，设置自适应阈值截取标准具阈值以上各峰，通过功率加权法，获得峰值对应采样点序列。再通过标准具Mark点，确定采样点与波长对应关系。在此基础上，截取光纤布拉格光栅归一化反射谱幅值大于峰值0.5至0.7倍范围内的区域，作为寻峰区域。设置初始参数范围，分别用高斯拟合法与基于探测器响应函数的拟合法对寻峰区域进行拟合。图5中曲线(1)～(6)的寻峰区域拟合效果如图6(1)～(6)所示，图中实心圆点为原始采样点序列，虚线对应高斯拟合结果，实线为基于探测器时域响应函数的拟合结果。与传统的高斯拟合法相比，响应函数拟合法平均拟合确定系数R2达到0.988，具有更高普适性，能同时适用于对称性良好的高斯谱型及不同程度的畸变谱型。
如图2所示，常见的光纤布拉格光栅畸变谱型曲线示意图。由于探测器对脉冲光信号的响应情况会直接影响光纤布拉格光栅反射峰谱型，因此，扫描光源频率较高时，探测器的输出光谱将出现对称性下降、谱型展宽和右端拖尾的畸变现象。
Claims (2)
1.一种基于光电探测器时域响应的光纤布拉格光栅弱信号解调方法，其特征在于，该方法包括以下步骤：
步骤(1)、构建基于可调谐F-P滤波器的弱光强解调系统，从而获得不同探测器增益下的归一化光纤布拉格光栅反射光谱，即采样点(1,2,…,N)-幅值序列X(A1,A2,…,AN)，设定光谱采样点幅值的对应阈值Φ，截取幅值大于阈值Φ的采样点序列P(Aa+1,Aa+2,…,Aa+n)作为寻峰区域，序列共包含n个点，其中a+1为序列X中超过阈值的起始采样点坐标，a+n为序列X中超过阈值的终止采样点坐标；
步骤(2)、以透射峰值波长作为波长参考，获取采样点(1,2,…,N)-波长值对应序列(λ1,λ2,…,λN)；
步骤(3)、构建探测器响应模型，仿真得到不同畸变程度光纤布拉格光栅反射峰谱型；表达式如下：
其中，符号*表示卷积运算，fin(t)表示时域输入信号，τ为在0至∞范围内变化的时间值，h(t)表示探测器时域冲激响应系数，通过探测器脉冲响应函数与脉冲函数反卷积得到探测器时域冲激响应系数h(t)，表达式如下:
h(t)＝α·exp(-t/β)
其中，α表示与探测器输出幅值相关的系数，β表示探测器时间响应系数的拟合值；
步骤(4)、拟合FBG反射光谱：以高斯函数与探测器冲激响应函数的卷积作为拟合函数，表达式如下：
其中，高斯函数包含对称轴μ、标准差σ两个参数，探测器响应函数包含幅值参数α、时间系数β两个参数；拟合方程为两函数、四个参数的卷积结果；式中，N表示一个周期内等时间间隔采集的采样点总数，m为在1和N之间变化的变量，n为1至2N-1之间的量；
步骤(5)、设置参数α、β、μ、σ范围及最佳寻峰阈值，以使得最佳寻峰阈值与步骤(1)中的阈值Φ相等；根据阈值Φ，仿真FBG光谱fFBG(n)中在阈值以上的点构成序列Q(Bb+1,Bb+2,…,Bb+m)，序列Q共包含m个点，其中，b+1为序列Q中超过阈值的起始横坐标，b+m为序列Q中超过阈值的终止横坐标；
通过线性插值，将步骤(1)中序列P转换为序列P’(Ab’+1,Ab’+2,…,Ab’+m)，其中，b’+1为变换后的起始横坐标，b’+m为变换后终止横坐标，且Ab’+1＝Aa+1，Ab’+m＝Aa+n，运用最小二乘法：获得寻峰区域拟合最优参数(μ’,σ’,α’,β’)，其中k从1递增至m；求取拟合结果最大值横坐标b’+k0，k0为1至m之间的值；通过三次样条插值，求取b’+k0在序列(λ1,λ2,…,λN)中的对应波长λ，即为解调波长。
2.如权利要求1所述的一种基于光电探测器时域响应的光纤布拉格光栅弱信号解调方法，其特征在于，所述阈值Φ的选取遵循高解调精度、低解调时间的准则，设定为0至1范围内的最佳值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810302314.5A CN108896078B (zh) | 2018-04-04 | 2018-04-04 | 基于探测器时域响应的光纤布拉格光栅弱信号解调方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810302314.5A CN108896078B (zh) | 2018-04-04 | 2018-04-04 | 基于探测器时域响应的光纤布拉格光栅弱信号解调方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108896078A true CN108896078A (zh) | 2018-11-27 |
| CN108896078B CN108896078B (zh) | 2020-04-03 |
Family
ID=64342315
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810302314.5A Active CN108896078B (zh) | 2018-04-04 | 2018-04-04 | 基于探测器时域响应的光纤布拉格光栅弱信号解调方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108896078B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110260898A (zh) * | 2019-06-24 | 2019-09-20 | 武汉理工光科股份有限公司 | 抗干扰的光栅波长解调方法及系统 |
| CN111854621A (zh) * | 2020-06-05 | 2020-10-30 | 北京航空航天大学 | 机载分布式pos用光纤光栅传感器数据拟合方法和装置 |
| CN114518163A (zh) * | 2022-02-21 | 2022-05-20 | 无边界(苏州)新材料科技有限公司 | 一种基于Gaussian-LM算法进行全声态光纤监测的方法 |
| CN114861723A (zh) * | 2022-05-07 | 2022-08-05 | 重庆邮电大学 | 一种应用于光纤布拉格光栅解调的系统及方法 |
| CN116907556A (zh) * | 2023-09-11 | 2023-10-20 | 武汉理工大学 | 分布式光纤传感多特征混合解调系统及方法 |
| CN117490740A (zh) * | 2023-12-29 | 2024-02-02 | 江西飞尚科技有限公司 | 一种光纤光栅调解方法及系统 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110255078A1 (en) * | 2007-10-23 | 2011-10-20 | Us Sensor Systems, Inc. | Interrogator for a plurality of sensor fiber optic gratings |
| US20140152995A1 (en) * | 2012-11-27 | 2014-06-05 | Sentek Instrument LLC | Serial weak fbg interrogator |
| WO2016086310A1 (en) * | 2014-12-04 | 2016-06-09 | Hifi Engineering Inc. | Optical interrogator for performing interferometry using fiber bragg gratings |
| CN105783953A (zh) * | 2016-03-25 | 2016-07-20 | 武汉理工大学 | 应用于光纤布拉格光栅波长解调的快速高斯拟合方法 |
| CN105973282A (zh) * | 2016-05-20 | 2016-09-28 | 武汉理工大学 | 光纤f-p传感器腔长小波相位提取解调方法 |
| US20170219390A1 (en) * | 2014-02-28 | 2017-08-03 | Hitachi, Ltd. | Optical fiber sensor device |
| CN107490397A (zh) * | 2016-09-14 | 2017-12-19 | 北京卫星环境工程研究所 | 高精度自适应滤波fbg光谱快速寻峰方法 |
| CN107560645A (zh) * | 2017-08-29 | 2018-01-09 | 北京航空航天大学 | 一种光纤布拉格光栅传感器波长解调寻峰方法 |
-
2018
- 2018-04-04 CN CN201810302314.5A patent/CN108896078B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110255078A1 (en) * | 2007-10-23 | 2011-10-20 | Us Sensor Systems, Inc. | Interrogator for a plurality of sensor fiber optic gratings |
| US20140152995A1 (en) * | 2012-11-27 | 2014-06-05 | Sentek Instrument LLC | Serial weak fbg interrogator |
| US20170219390A1 (en) * | 2014-02-28 | 2017-08-03 | Hitachi, Ltd. | Optical fiber sensor device |
| WO2016086310A1 (en) * | 2014-12-04 | 2016-06-09 | Hifi Engineering Inc. | Optical interrogator for performing interferometry using fiber bragg gratings |
| CN105783953A (zh) * | 2016-03-25 | 2016-07-20 | 武汉理工大学 | 应用于光纤布拉格光栅波长解调的快速高斯拟合方法 |
| CN105973282A (zh) * | 2016-05-20 | 2016-09-28 | 武汉理工大学 | 光纤f-p传感器腔长小波相位提取解调方法 |
| CN107490397A (zh) * | 2016-09-14 | 2017-12-19 | 北京卫星环境工程研究所 | 高精度自适应滤波fbg光谱快速寻峰方法 |
| CN107560645A (zh) * | 2017-08-29 | 2018-01-09 | 北京航空航天大学 | 一种光纤布拉格光栅传感器波长解调寻峰方法 |
Non-Patent Citations (4)
| Title |
|---|
| 张红霞等: "光纤布拉格光栅传感器的一种波长解调方法", 《天津大学学报》 * |
| 朱梅等: "光纤布拉格光栅中心波长检测中的寻峰算法", 《光通信研究》 * |
| 王梓蒴: "光纤光栅传感信号解调方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
| 范伟凯: "光纤光栅传感器波长检测与解调系统的研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110260898A (zh) * | 2019-06-24 | 2019-09-20 | 武汉理工光科股份有限公司 | 抗干扰的光栅波长解调方法及系统 |
| CN110260898B (zh) * | 2019-06-24 | 2021-07-06 | 武汉理工光科股份有限公司 | 抗干扰的光栅波长解调方法及系统 |
| CN111854621A (zh) * | 2020-06-05 | 2020-10-30 | 北京航空航天大学 | 机载分布式pos用光纤光栅传感器数据拟合方法和装置 |
| CN111854621B (zh) * | 2020-06-05 | 2021-10-15 | 北京航空航天大学 | 机载分布式pos用光纤光栅传感器数据拟合方法和装置 |
| CN114518163A (zh) * | 2022-02-21 | 2022-05-20 | 无边界(苏州)新材料科技有限公司 | 一种基于Gaussian-LM算法进行全声态光纤监测的方法 |
| CN114518163B (zh) * | 2022-02-21 | 2024-03-19 | 无边界(苏州)新材料科技有限公司 | 一种基于Gaussian-LM算法进行全声态光纤监测的方法 |
| CN114861723A (zh) * | 2022-05-07 | 2022-08-05 | 重庆邮电大学 | 一种应用于光纤布拉格光栅解调的系统及方法 |
| CN116907556A (zh) * | 2023-09-11 | 2023-10-20 | 武汉理工大学 | 分布式光纤传感多特征混合解调系统及方法 |
| CN116907556B (zh) * | 2023-09-11 | 2024-04-16 | 武汉理工大学 | 分布式光纤传感多特征混合解调系统及方法 |
| CN117490740A (zh) * | 2023-12-29 | 2024-02-02 | 江西飞尚科技有限公司 | 一种光纤光栅调解方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108896078B (zh) | 2020-04-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108896078A (zh) | 基于探测器时域响应的光纤布拉格光栅弱信号解调方法 | |
| CN106940444B (zh) | 基于微波差分增益的相干多普勒测风激光雷达 | |
| US7787779B2 (en) | Photonic time-domain electromagnetic signal generator and system using the same | |
| CN107219002B (zh) | 一种超高分辨率光谱测量方法及系统 | |
| CN103176173B (zh) | 基于光纤采样技术的lfmcw激光雷达调频的非线性校正方法 | |
| CN106907997B (zh) | 一种基于光纤位移传感系统的位移测量信号分析方法 | |
| GB2517100A (en) | Method and apparatus for optical sensing | |
| CN103592261A (zh) | 全光纤温度补偿型气体传感器及其补偿方法 | |
| CN109883458A (zh) | 一种采用新型光学微波鉴频器和新型扰偏器的布里渊传感系统 | |
| CN105577280B (zh) | 一种光载微波信号动态宽频实时数字解调系统 | |
| CN102589588A (zh) | 利用光纤光栅解调Fabry–Pérot腔腔长的方法 | |
| CN103414513B (zh) | 一种具有高动态范围的脉冲光动态消光比测量装置及方法 | |
| Laghezza et al. | Field evaluation of a photonics‐based radar system in a maritime environment compared to a reference commercial sensor | |
| CN101241029A (zh) | 光纤Bragg光栅传感器解调器 | |
| CN101419317B (zh) | 一种基于光纤Bragg光栅的双边缘滤波器 | |
| CN113391136A (zh) | 一种基于固定低频检测的微波光子频率测量装置及方法 | |
| CN116839758A (zh) | 高信噪比高精度的光纤传感解调系统及其实现方法 | |
| CN110082068A (zh) | 一种具有波长修正功能的光纤光栅波长解调系统及方法 | |
| CN103398736A (zh) | 光电探测器频率响应的测量系统 | |
| CN105606345A (zh) | 基于波长编码技术光电探测器频响的测试装置及测试方法 | |
| CN102607702A (zh) | 宽带参考光源光频域游标法光谱仪 | |
| CN104848879A (zh) | 基于线性工作匹配光栅法的光纤布拉格光栅传感器信号解调方法 | |
| CN108759879B (zh) | 一种基于光栅传感器的波长解析装置 | |
| CN203323891U (zh) | 一种基于awg和光开关的光波长计 | |
| CN113567955B (zh) | 基于单腔双工作波长fpi的水体探测激光雷达 |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |