阅读说明：本技术 一种光纤布拉格光栅传感系统的信号解调方法 (Signal demodulation method of optical fiber Bragg grating sensing system ) 是由 尚秋峰 秦文婕 于 2019-10-10 设计创作，主要内容包括：一种光纤布拉格光栅传感系统的信号解调方法,所述方法采用兼容平滑与边缘保持的降噪方法,去除FBG传感信号中的噪声；然后对降噪后的数据进行迭代拟合区域分割,确定所有单峰区域；最后以修正后谱曲线的斜率比归1度为每次迭代的目标函数,利用加权的一阶导数对单峰区域的非对称谱进行锐化修正,对修正后的对称谱进行谱峰解调,得到FBG信号的中心波长位置。本发明采用兼容平滑与边缘保持的降噪方法,在去除噪声干扰的同时还能保护信号的特征信息不丢失；迭代拟合区域分割克服了传统方法复杂的递回,分割区域不对称等缺点,提高了分割效率和分割准确性；通过锐化修正非对称谱提高谱峰解调精度,从而提高了光纤布拉格光栅传感系统的检测精度。(A signal demodulation method of a fiber Bragg grating sensing system adopts a noise reduction method compatible with smoothness and edge preservation to remove noise in an FBG sensing signal; then, carrying out iterative fitting region segmentation on the denoised data to determine all single-peak regions; and finally, taking the slope ratio of the corrected spectral curve to 1 degree as a target function of each iteration, carrying out sharpening correction on the asymmetric spectrum of the single-peak region by using the weighted first derivative, and carrying out spectral peak demodulation on the corrected symmetric spectrum to obtain the central wavelength position of the FBG signal. The invention adopts a noise reduction method compatible with smoothness and edge maintenance, which can remove noise interference and simultaneously protect the characteristic information of the signal from losing; the iterative fitting region segmentation overcomes the defects of complex recursion, asymmetric segmentation region and the like of the traditional method, and improves the segmentation efficiency and the segmentation accuracy; the spectrum peak demodulation precision is improved by sharpening and correcting the asymmetric spectrum, so that the detection precision of the fiber Bragg grating sensing system is improved.)

1. a signal demodulation method of a fiber Bragg grating sensing system is characterized in that a noise reduction method compatible with smoothness and edge preservation is adopted to remove high-frequency noise and low-frequency noise in an FBG sensing signal; then, carrying out iterative fitting region segmentation on the denoised data to determine all single-peak regions; and finally, taking the slope ratio of the corrected spectral curve to 1 degree as a target function of each iteration, carrying out sharpening correction on the asymmetric spectrum of the single-peak region by using the weighted first derivative, and carrying out spectral peak demodulation on the corrected symmetric spectrum to obtain the central wavelength position of the FBG signal.

2. The method of claim 1, wherein the signal demodulation method comprises the following steps:

a. noise reduction processing compatible with smoothing and edge preservation

In a plane coordinate system, a collected signal is a set of a column of discrete data, the width of a noise reduction sliding window is set as n, the length of a sliding array of noise reduction processing is set as w ═ 2n +1, the sliding window contains data points x ═ n, -n +1,, 0,1, … n-1, n, and each data point in the width is fitted by using a high-order polynomial to obtain: wherein q is the polynomial degree, ai is the coefficient of the ith polynomial of the polynomial, w equations can be written out by moving the sliding window once, and if the residual error of least square fitting is set as e, then:

that is, f is taken as an original signal, X is a data point matrix, a is a coefficient matrix, the least square solution of a is (XT · X) -1 · XT · f, E is a residual matrix of least square fitting, and the signal after noise reduction is f '═ X · a' ═ X · (XT · X) -1 · XT · f;

b. Iterative fitting region segmentation on denoised data

Setting a sensing channel consisting of R FBGs connected in series, firstly carrying out maximum value search on a signal f ' subjected to noise reduction to obtain R maximum values, starting from a first maximum value point, respectively selecting fitting data point numbers Nj on the left and the right of the signal f ', wherein Nj is used as a control factor of a fitting area and changes every iteration, wherein j (j is more than or equal to 0) is iteration times, N0 is an initial value when j is 0, enabling Ns to be the total number of data points in the signal f ', carrying out Gaussian fitting once when 2Nj +1 changes in a range of [ N0, Nmax ], finishing fitting every time to calculate an error epsilon, finding out the value of Nj when epsilon is minimum, and determining a corresponding unimodal area;

c. sharpening correction for asymmetric spectrum of single-peak region

Let S0 represent a single peak, S0 is corrected using the first derivative S0' of S0: s1 ═ S0+ C1 · S0', where S1 is the corrected signal, C1 is the iteratively obtained weighting factor, and the absolute value Q of the slope ratio on both sides of the peak is calculated for each correction:

wherein ρ l is the slope of the left side of the peak at S1, Sl 'and S' l are the first and second derivatives, respectively, of the spectral data left side of the peak at S1, ρ r is the slope of the right side of the peak at S1, and Sr 'and S' r are the first and second derivatives, respectively, of the spectral data right side of the peak at S1;

Comparing Q with 1, if Q is greater than 1, increasing C1, and if Q is less than 1, decreasing C1 until Q is 1, and obtaining a corrected symmetric spectrum S1;

and performing spectral peak demodulation on the S1 to obtain the central wavelength positions of all the FBGs in one sensing channel.

3. The method of claim 2, wherein the sliding window width w is greater than the polynomial order q.

Technical Field

The invention relates to a signal demodulation method of a fiber Bragg grating sensing system, which is not only suitable for symmetrical Gaussian spectrum type sensing signals with high signal-to-noise ratio, but also suitable for asymmetrical Gaussian spectrum type sensing signals with serious noise influence, and belongs to the technical field of measurement.

Background

At present, a sensing system based on Fiber Bragg Gratings (FBGs) is widely applied to measurement and monitoring of physical parameters such as strain, temperature and acceleration. The accurate central wavelength position is obtained through an FBG demodulation algorithm, and the FBG demodulation algorithm plays an important role in the measurement accuracy of the sensing system. The traditional demodulation algorithm extracts the central wavelength characteristic, has a good effect on a symmetrical Gaussian spectrum type sensing signal with high signal-to-noise ratio, but is not suitable for an asymmetrical Gaussian spectrum type sensing signal seriously influenced by noise.

For the FBG sensing system with a large number of nodes, the signal-to-noise ratio of the sensing signal is low, and the traditional demodulation method has the problems of FBG signal characteristic point fuzzification, spectrum type detail information loss, single-peak area asymmetry, long operation time, poor measurement accuracy of asymmetric Gaussian spectrum type signals and the like.

Disclosure of Invention

The invention aims to provide a signal demodulation method of a fiber Bragg grating sensing system aiming at the defects of the prior art so as to improve the detection precision of the fiber Bragg grating sensing system.

The problems of the invention are solved by the following technical scheme:

a signal demodulation method of a fiber Bragg grating sensing system comprises three links of noise reduction, region segmentation and spectral peak demodulation; the method comprises the steps that firstly, a noise reduction method compatible with smoothing and edge keeping is adopted to remove high-frequency noise and low-frequency noise in FBG sensing signals; then, carrying out iterative fitting region segmentation on the denoised data to determine all single-peak regions; and finally, taking the slope ratio of the corrected spectral curve to 1 degree as a target function of each iteration, carrying out sharpening correction on the asymmetric spectrum of the single-peak region by using the weighted first derivative, and carrying out spectral peak demodulation on the corrected symmetric spectrum to obtain the central wavelength position of the FBG signal.

The signal demodulation method of the fiber Bragg grating sensing system comprises the following steps:

a. Noise reduction processing compatible with smoothing and edge preservation

in a plane coordinate system, a collected signal is a set of a column of discrete data, the width of a noise reduction sliding window is set as n, the length of a sliding array of noise reduction processing is set as w ═ 2n +1, the sliding window contains data points x ═ n, -n +1,, 0,1, … n-1, n, and each data point in the width is fitted by using a high-order polynomial to obtain: wherein q is the polynomial degree, ai is the coefficient of the ith polynomial of the polynomial, w equations can be written out by moving the sliding window once, and if the residual error of least square fitting is set as e, then:

that is, f is taken as an original signal, X is a data point matrix, a is a coefficient matrix, the least square solution of a is (XT · X) -1 · XT · f, E is a residual matrix of least square fitting, and the signal after noise reduction is f '═ X · a' ═ X · (XT · X) -1 · XT · f;

b. Iterative fitting region segmentation on denoised data

setting a sensing channel consisting of R FBGs connected in series, firstly carrying out maximum value search on a signal f ' subjected to noise reduction to obtain R maximum values, starting from a first maximum value point, respectively selecting fitting data point numbers Nj on the left and the right of the signal f ', wherein Nj is used as a control factor of a fitting area and changes every iteration, wherein j (j is more than or equal to 0) is iteration times, N0 is an initial value when j is 0, enabling Ns to be the total number of data points in the signal f ', carrying out Gaussian fitting once when 2Nj +1 changes in a range of [ N0, Nmax ], finishing fitting every time to calculate an error epsilon, finding out the value of Nj when epsilon is minimum, and determining a corresponding unimodal area;

c. Sharpening correction for asymmetric spectrum of single-peak region

Let S0 represent a single peak, S0 is corrected using the first derivative S0' of S0: s1 ═ S0+ C1 · S0', where S1 is the corrected signal, C1 is the iteratively obtained weighting factor, and the absolute value Q of the slope ratio on both sides of the peak is calculated for each correction:

wherein ρ l is the slope of the left side of the peak at S1, Sl 'and Sl "are the first and second derivatives, respectively, of the spectral data left side of the peak at S1, ρ r is the slope of the right side of the peak at S1, and Sr' and Sr" are the first and second derivatives, respectively, of the spectral data right side of the peak at S1;

Comparing Q with 1, if Q is greater than 1, increasing C1, and if Q is less than 1, decreasing C1 until Q is 1, and obtaining a corrected symmetric spectrum S1; and performing spectral peak demodulation on the S1 to obtain the central wavelength positions of all the FBGs in one sensing channel.

According to the signal demodulation method of the fiber Bragg grating sensing system, the width w of the sliding window is larger than the polynomial order q.

The invention adopts a noise reduction method compatible with smoothness and edge maintenance, which can remove noise interference and simultaneously protect the characteristic information of the signal from losing; the iterative fitting region segmentation overcomes the defects of complex recursion, asymmetric segmentation region and the like of the traditional method, and improves the segmentation efficiency and the segmentation accuracy; the spectrum peak demodulation precision is improved by sharpening and correcting the asymmetric spectrum, so that the detection precision of the fiber Bragg grating sensing system is greatly improved.

Drawings

the present invention will be described in further detail with reference to the accompanying drawings.

FIG. 1 is a signal flow diagram;

FIG. 2 is a flow chart of iterative fitting region segmentation;

FIG. 3 is a schematic diagram of a sharpening modified asymmetric spectrum.

the symbols used herein are: n is the width of the noise reduction sliding window, w is the sliding array length of the noise reduction process, X is the data points contained in the sliding window, X is the data point matrix, fw is the polynomial fitted to the data points, Q is the polynomial degree, ai is the coefficient of the ith term of the polynomial, X is the unknown number of the polynomial, f is the collected original signal, f 'is the noise-reduced signal, Ns is the total number of data points in the signal f', E is the residual error of the least square fit, a is the coefficient matrix, a 'is the least square solution of a, E is the residual matrix of the least square fit, Nj is the data point number, N0 is the initial value of Nj when j is 0, epsilon is the error, S0 represents a single peak, S0' represents the first derivative of S0, S1 represents the modified symmetric spectrum, C1 is the weighting factor obtained by iteration, Q is the absolute value of the slope ratio of both sides of the S1 peak, ρ l is the slope ratio of the left side of the S1 peak, sl 'and Sl "are the first and second derivatives, respectively, of the spectral data to the left of the S1 peak, ρ r is the slope to the right of the S1 peak, Sr' and Sr" are the first and second derivatives, respectively, of the spectral data to the right of the S1 peak, and β is the iteration step.

Detailed Description

the invention provides a novel demodulation method of an FBG sensing signal with a low signal-to-noise ratio, which can effectively improve the detection precision of an FBG sensing system.

the method comprises three parts of a noise reduction method compatible with smoothing and edge preservation, iterative fitting region segmentation and sharpening correction asymmetric spectrum.

The method is compatible with a smoothing and edge-keeping noise reduction method, removes high-frequency noise, emphasizes on the elimination of low-frequency noise, protects FBG signal feature point data from being eliminated, and can ensure that the shape and the width of a signal are unchanged. The core idea is a weighted average algorithm of a moving window, but the weighting coefficient is not a simple constant, but is obtained by performing high-order polynomial fitting on data points in a window sliding for a certain length. By adjusting the sliding window and the order of the polynomial, the noise reduction result compatible with the smoothness and the edge preservation can be realized.

And the iterative fitting region segmentation is to iteratively obtain the optimal solution of the region segmentation by comparing the fitting result with the standard deviation of the original data. Since the unimodal spectrum of the FBG sensing signal is not always symmetric, the accuracy of the spectrum demodulation is limited. And the iterative fitting region segmentation firstly carries out maximum value search on the data after noise reduction, iteratively selects the number of data near the maximum value after the position of the maximum value is detected, respectively fits the data points near the maximum value, and when the fitting error is minimum, the data are the best single-peak regions, and the purpose of region segmentation is achieved after all the single-peak regions are determined to be finished.

Sharpening the modified asymmetric spectrum is mainly characterized in that: for the asymmetric spectrum, sharpening modification takes the slope ratio of a modified spectral curve to 1 degree as a target function of each iteration, the original spectral curve is continuously modified by calculating the derivative of the original spectral curve, spectral characteristic parameters (spectral height, spectral width and spectral area) are adaptively re-calibrated, the processed signal is symmetrical left and right, the spectral width is narrowed, and the detection precision of the processed FBG sensing signal with low signal-to-noise ratio is improved. Sharpening modifies the asymmetric spectrum, the original spectrum is modified by solving a derivative of the unimodal original curve as a modification factor, and the existence of the derivative modification factor enables the modified spectrum width to be narrowed, so that the sharpening effect is achieved. The correction weight of the derivative correction factor is adaptively and iteratively adjusted in the invention until the final spectral slope ratio is detected to be 1, and the correction is finished.

The invention comprises the following 3 steps:

Step 1: and noise reduction processing compatible with smoothing and edge maintaining is carried out, and the signal to noise ratio of the FBG sensing signal is improved. In a planar coordinate system, the acquired signal is a set of a series of discrete data. Assuming that the width of the noise reduction sliding window is n, the length of the sliding array of the noise reduction processing is w ═ 2n +1, and the data points included in the sliding window are x ═ (-n, -n +1,, 0,1, … n-1, n), fitting each data point in the width by using a high-order polynomial to obtain: and q is polynomial degree, the smoothing intensity is controlled, the smaller q, the smoother the fitted curve is, the poorer the edge retentivity is, and conversely, the larger q, the better the marginality is, and the poorer the smoothness is. The sliding window moves once to write w equations, and if the residual error of least square fitting is set as e, then:

That is, f is X · a + E, the least-squares solution of a is a ' ═ (XT · X) -1 · XT · f, f is regarded as the collected original signal, and is denoised to obtain f ' ═ X · a ' ═ X · (XT · X) -1 · XT · f. The purpose of compatible smoothing and edge keeping can be achieved by adjusting the sizes of w and q, when the value of w is larger, the processed signal is smoother, the edge performance is poorer, but the larger the value of q is, the better the edge performance is, and the smoother is. Meanwhile, if the equation set has a solution, w is more than or equal to q, w is generally selected to be more than q, in the selected embodiment of the invention, the sliding window w is 5, and the polynomial order q is 3.

fig. 2 shows a flow of iterative fitting region segmentation. Aiming at the problem of single-peak region asymmetry after Hilbert transform segmentation, a comprehensive iteration and Gaussian fitting solving technology is provided. The technology dynamically selects the fitting range of each single peak according to the principle of minimum error so as to optimize the fitting of the single peak area and avoid the area coincidence between two peaks.

Step 2: the noise-reduced signal f' is first subjected to a maximum search. A sensing channel is formed by R FBGs connected in series, so that the signal f' contains R maxima, one maximum corresponding to a single peak. Starting from the first maximum value point, selecting the number Nj of fitting data points on the left and right of the first maximum value point respectively, wherein the Nj is used as a control factor of a fitting area and changes every iteration, and j (j is more than or equal to 0) is the iteration frequency. N0 is the starting value when j is 0. The starting value N0 and the iteration step β of Nj are selected according to the sampling rate of the original data f and engineering experience, and in the embodiment of the present invention, N0 is 10, and β is 10. Let Ns be the total number of data points in signal f'. And 2Nj +1 is subjected to Gaussian fitting once every time when the variation is within the range of [ N0, Nmax ], the calculation error epsilon is fitted every time, the value of Nj when epsilon is minimum is found out, the region corresponding to a single peak can be determined, and for R maximum value points of the signal f', the region segmentation is completed by determining R single peak regions in total.

FIG. 3 shows a flow of sharpening the modified asymmetric spectrum. Due to the spectral asymmetry of the FBG sensing signal, the accuracy of the demodulation is greatly limited. The sharpening correction asymmetric spectrum technology corrects the spectrum type into a symmetric spectrum so as to accurately demodulate a target spectrum peak. Only the derivative factor correction is needed to be carried out on the original spectrum function, and a sharper spectrum peak can be obtained by utilizing the edge enhancement technology. Sharpening correction of the asymmetric spectrum is an iterative process, the slope ratio of a corrected spectrum curve is reduced to 1 degree as a target function of each iteration, and the symmetric spectrum is obtained through continuous correction.

And step 3: the method is a unimodal spectrum processing after the signal f 'obtained by noise reduction is subjected to region segmentation, wherein S0 is one of the unimodal (the unimodal is not necessarily a symmetric spectrum), S1 is a modified symmetric spectrum, if a spectrum curve is asymmetric, absolute values of slopes on two sides of the spectrum curve are not equal, and the absolute value of the slope on one side is larger than that on the other side, then the asymmetric spectrum can be modified by weighted addition of the first derivative term S0' of S0. Since the first derivative of the spectral curve is antisymmetric (positive on the one hand and negative on the other hand) as shown in fig. 3, the slope is lower on the right side of the original spectral curve S0 (solid line) than on the left side, and the first derivative S0' (dashed line) has a positive lobe on the left side and a wider but smaller negative lobe on the right side. When the original spectral curve is added with the weighted summation of the first derivative, the positive lobe of the derivative reinforces the front edge, and the negative lobe restrains the back edge, the processed S1 is the symmetric spectral curve, and the introduction of the derivative term narrows the width of the spectral curve to achieve the sharpening effect. The process using the first derivative is: s1 is S0+ C1 · S0', where C1 is an iteratively obtained weighting factor, and the criterion for whether the correction is completed is to see whether the absolute value Q of the slope ratio on both sides of the peak is 1, and if Q is greater than 1, increase C1, and vice versa. And finally, demodulating the spectrum peak of the corrected symmetrical spectrum to demodulate the central wavelength positions of all the FBGs of one sensing channel.

wherein the content of the first and second substances,

ρ l is the slope of the left side of the peak at S1, Sl 'and Sl "are the first and second derivatives, respectively, of the spectral data left side of the peak at S1, and similarly ρ r is the slope of the right side of the peak at S1, and Sr' and Sr" are the first and second derivatives, respectively, of the spectral data right side of the peak at S1.

The invention has the following advantages:

1. The method is compatible with the noise reduction processing of smoothing and edge keeping, removes noise interference and simultaneously can protect the characteristic information of the signal from being lost.

2. the iterative fitting region segmentation overcomes the defects of complex recursion, asymmetric segmentation region and the like of the traditional method, and is accurate and efficient in segmentation.

3. The asymmetric spectrum can be effectively corrected by sharpening, so that a symmetric spectrum can be obtained, and the spectrum peak demodulation precision is improved.