CN106643835A - Fiber fabry-perot cavity demodulation method, device thereof and fiber fabry-perot interferometer - Google Patents
Fiber fabry-perot cavity demodulation method, device thereof and fiber fabry-perot interferometer Download PDFInfo
- Publication number
- CN106643835A CN106643835A CN201510740349.3A CN201510740349A CN106643835A CN 106643835 A CN106643835 A CN 106643835A CN 201510740349 A CN201510740349 A CN 201510740349A CN 106643835 A CN106643835 A CN 106643835A
- Authority
- CN
- China
- Prior art keywords
- chamber
- value
- fabry
- optical fiber
- length
- 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
- 238000000034 method Methods 0.000 title claims abstract description 64
- 239000000835 fiber Substances 0.000 title claims abstract description 22
- 230000003595 spectral effect Effects 0.000 claims abstract description 139
- 239000013307 optical fiber Substances 0.000 claims description 139
- 238000010276 construction Methods 0.000 claims description 32
- 230000002452 interceptive effect Effects 0.000 claims 2
- 238000001228 spectrum Methods 0.000 description 21
- 238000010586 diagram Methods 0.000 description 7
- 230000009466 transformation Effects 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 230000003287 optical effect Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- PNEYBMLMFCGWSK-UHFFFAOYSA-N Alumina Chemical compound [O-2].[O-2].[O-2].[Al+3].[Al+3] PNEYBMLMFCGWSK-UHFFFAOYSA-N 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000009530 blood pressure measurement Methods 0.000 description 2
- 238000011088 calibration curve Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000001831 conversion spectrum Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005538 encapsulation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 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/35312—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 Fabry Perot
Abstract
The present invention discloses a fiber fabry-perot cavity demodulation method, a device thereof and a fiber fabry-perot interferometer. The fiber fabry-perot cavity demodulation method comprises the steps of acquiring a spectral signal of an interference fiber fabry-perot cavity; subjecting the spectral signal to first demodulation through adopting non-equally dispersed fourier transform so as to obtain a first cavity length value of the fiber fabry-perot cavity, wherein the accuracy of the first cavity length value is greater than a first threshold value and the cavity length demodulation of the fiber fabry-perot cavity is free of the half-wave hopping phenomenon when the accuracy of the first cavity length value is greater than the first threshold value; and according to the first cavity length value, subjecting the spectral signal to second demodulation so as to obtain the second cavity length value of the fiber fabry-perot cavity, wherein the accuracy of the second cavity length value of the fiber fabry-perot cavity is greater than the accuracy of the first cavity length value. According to the technical scheme of the invention, the problem in the prior art that the cavity length demodulation accuracy of the fiber fabry-perot cavity is not high can be solved. Therefore, the demodulation accuracy of the cavity length of the fiber fabry-perot cavity is improved.
Description
Technical field
The present invention relates to Fabry-perot optical fiber interferometer field, in particular to a kind of Fabry-perot optical fiber chamber demodulation method and dress
Put and Fabry-perot optical fiber interferometer.
Background technology
Fabry-perot optical fiber interferometer (i.e. Fiber-optical Fabry-Perot Interferometric Cavity, referred to as optical fiber FP interferometers), as one kind
The advantages of fiber optic sensing device, good stability simple due to making, high precision, is widely used to various components or ring
Detect in border.According to practical application needs, with reference to different encapsulation, optical fiber FP interferometers can become temperature sensor, press
Strong sensor, stress/strain sensor and vibrating sensor etc..As respective environment parameter (such as temperature, pressure, is answered
Become etc.) change can cause the change of FP chambers chamber length, therefore it is long to demodulate chamber by interference spectrum, just can obtain and treat
Detection limit information.
The major way of FP chambers demodulation has intensity demodulation and phase demodulating, and intensity demodulation is the technology for using earliest, simply
Easy but error is larger, what is be currently mainly used is the higher phase demodulating of precision.For wide range optical interface fiber FP
Chamber, conventional spectral analysis method have Fringe-counting method, Wavelength-tracking method, fourier transform method and approximating method and this
The improvement and combination of a little methods.Fringe-counting method obtains the absolute of chamber length according to the relation of fringe period and wavelength
Value, but precision is very low；Wavelength-tracking method precision is of a relatively high, but the only dynamic range of half-wavelength, and has half
Ripple saltus step, wherein, half-wave saltus step refer to have two as algorithm factor causes the chamber long value for demodulating and actual value/
The deviation of one wavelength；Periodic signal is transformed to its frequency domain by fourier transform method, that is, correspond to chamber long message, but its mistake
Difference is not also little.Approximating method is most connect by building accurate demodulation model, scanning construction signal and actual signal
Chamber length output during nearly actual signal, approximating method mainly have method of correlation and method of least square etc., and this kind of demodulation method can
Sub- nano-precision is reached, but the big speed of operand is slower, and half-wave hopping phenomenon also occurs, affect cavity length demodulating
Precision.
For the not high problem of Fabry-perot optical fiber chamber cavity length demodulating precision in correlation technique, not yet propose at present effectively to solve
Scheme.
The content of the invention
Present invention is primarily targeted at a kind of Fabry-perot optical fiber chamber demodulation method and device and Fabry-perot optical fiber interferometer are provided,
To solve the problems, such as that Fabry-perot optical fiber chamber cavity length demodulating precision is not high in correlation technique.
To achieve these goals, according to an aspect of the invention, there is provided a kind of Fabry-perot optical fiber chamber demodulation method.
Fabry-perot optical fiber chamber of the invention demodulation method includes：Obtain the spectral signal of interference optical fiber Fa-Po cavity；Using
The discrete Fourier transform pairs spectral signal of Unequal distance carries out the first demodulation, obtains the first chamber long value in Fabry-perot optical fiber chamber, its
In, the precision of the first chamber long value is more than first threshold, the Fabry-perot optical fiber when the precision of the first chamber long value is more than first threshold
The cavity length demodulating in chamber is not in half-wave saltus step；And the second demodulation is carried out to spectral signal according to the first chamber long value, obtain
To the second chamber long value in Fabry-perot optical fiber chamber, wherein, the precision of the second chamber long value is more than the precision of the first chamber long value.
Further, the first demodulation is carried out to spectral signal, the first chamber long value for obtaining Fabry-perot optical fiber chamber includes：Using
The discrete Fourier transform pairs spectral signal of Unequal distance carries out the first demodulation, obtains the first chamber long value in Fabry-perot optical fiber chamber.
Further, the first demodulation is being carried out using the discrete Fourier transform pairs spectral signal of Unequal distance, is obtaining Fiber Optic Sensor
Before the first chamber long value in amber chamber, method also includes：Spectral signal is multiplied with window function, the first spectral signal is obtained,
The first demodulation is carried out using the discrete Fourier transform pairs spectral signal of Unequal distance includes：Become using Unequal distance discrete fourier
Changing carries out the first demodulation to the first spectral signal, obtains the first chamber long value in Fabry-perot optical fiber chamber.
Further, the second demodulation is carried out to spectral signal according to the first chamber long value, obtains second chamber in Fabry-perot optical fiber chamber
Long value includes：Obtain the construction signal of spectral signal；Spectral signal is calculated according to the first chamber long value and signal is constructed most
Little mean square deviation estimated value；And the second chamber long value in Fabry-perot optical fiber chamber is determined according to nonlinear IEM model value.
Further, the nonlinear IEM model value for calculating spectral signal and construction signal according to the first chamber long value includes：
Obtain the first step-length；Sweep limitss are determined according to the first chamber long value and the first step-length；The second step-length is obtained, wherein, the
Two step-lengths are less than the first step-length；Value is carried out every second step length from sweep limitss, multiple values are obtained；Difference root
According to multiple exploitation spectral signals and the mean square deviation estimated value of construction signal, multiple mean square deviation estimated values are obtained；Obtain
In multiple mean square deviation estimated values, minimum mean square deviation estimated value is used as nonlinear IEM model value.
Further, method also includes：Judge the second step-length whether more than Second Threshold；Judging that second step grows up
When Second Threshold, continue to obtain the 3rd step-length, wherein, the 3rd step-length is less than the second step-length；According to multiple mean square deviations
In estimated value, minimum mean square deviation estimated value determines the 3rd chamber long value in Fabry-perot optical fiber chamber；And according to the 3rd chamber long value and
Second step-length determines sweep limitss, wherein, value is carried out every the 3rd step-length from sweep limitss.
To achieve these goals, according to a further aspect in the invention, there is provided a kind of Fabry-perot optical fiber chamber demodulating equipment.
The Fabry-perot optical fiber chamber demodulating equipment includes：First acquisition unit, for obtaining the spectral signal of interference optical fiber Fa-Po cavity；
First demodulating unit, for carrying out the first demodulation using the discrete Fourier transform pairs spectral signal of Unequal distance, obtains optical fiber
First chamber long value of Fa-Po cavity, wherein, the precision of the first chamber long value is more than first threshold, in the precision of the first chamber long value
Be not in half-wave saltus step more than the cavity length demodulating in Fabry-perot optical fiber chamber during first threshold；And second demodulating unit, it is used for
Second demodulation is carried out to spectral signal according to the first chamber long value, the second chamber long value in Fabry-perot optical fiber chamber is obtained, wherein, the
Precision of the precision of two chamber long values more than the first chamber long value.
Further, the first demodulating unit is for carrying out the first solution using the discrete Fourier transform pairs spectral signal of Unequal distance
Adjust, obtain the first chamber long value in Fabry-perot optical fiber chamber.
Further, device also includes：Arithmetic element, for spectral signal is multiplied with window function, obtains the first light
Spectrum signal, the first demodulating unit for carrying out the first demodulation using discrete the first spectral signal of Fourier transform pairs of Unequal distance,
Obtain the first chamber long value in Fabry-perot optical fiber chamber.
Further, the second demodulating unit includes：Acquisition module, for obtaining the construction signal of spectral signal；Calculate
Module, for the nonlinear IEM model value of spectral signal and construction signal is calculated according to the first chamber long value；And determine
Module, for the second chamber long value in Fabry-perot optical fiber chamber is determined according to nonlinear IEM model value.
Further, computing module includes：First acquisition submodule, for obtaining the first step-length；First determines submodule
Block, for determining sweep limitss according to the first chamber long value and the first step-length；Second acquisition submodule, for obtaining second
Step-length, wherein, the second step-length is less than the first step-length；Value submodule, for from sweep limitss every the second step-length
Value is carried out, multiple values are obtained；Calculating sub module, for respectively according to multiple exploitation spectral signals and construction
The mean square deviation estimated value of signal, obtains multiple mean square deviation estimated values；Obtain minimum mean square in multiple mean square deviation estimated values
Difference estimated value is used as nonlinear IEM model value.
Further, device also includes：Judging unit, for judging the second step-length whether more than Second Threshold；Second
Acquiring unit, for when judging that second step length is more than Second Threshold, obtaining the 3rd step-length, wherein, the 3rd step-length
Less than the second step-length；First determining unit, it is true for the mean square deviation estimated value according to minimum in multiple mean square deviation estimated values
Determine the 3rd chamber long value in Fabry-perot optical fiber chamber；And second determining unit, for true according to the 3rd chamber long value and the second step-length
Determine sweep limitss, wherein, value is carried out every the 3rd step-length from sweep limitss.
According to a further aspect of the invention, there is provided a kind of Fabry-perot optical fiber interferometer, Fabry-perot optical fiber interferometer include
State any one Fabry-perot optical fiber chamber demodulating equipment.
By the present invention, using the spectral signal for obtaining interference optical fiber Fa-Po cavity；Using Unequal distance discrete Fourier transform
First demodulation is carried out to spectral signal, the first chamber long value in Fabry-perot optical fiber chamber is obtained, wherein, the precision of the first chamber long value
More than first threshold, when the precision of the first chamber long value is more than first threshold, the cavity length demodulating in Fabry-perot optical fiber chamber is not in
Half-wave saltus step；And the second demodulation is carried out to spectral signal according to the first chamber long value, obtain second chamber in Fabry-perot optical fiber chamber
Long value, wherein, the precision of the second chamber long value solves Fabry-perot optical fiber in correlation technique more than the precision of the first chamber long value
The not high problem of chamber cavity length demodulating precision, and then reached the effect for improving Fabry-perot optical fiber chamber cavity length demodulating precision.
Description of the drawings
The accompanying drawing for constituting the part of the application is used for providing a further understanding of the present invention, the schematic reality of the present invention
Apply example and its illustrate, for explaining the present invention, not constituting inappropriate limitation of the present invention.In the accompanying drawings：
Fig. 1 is the flow chart of Fabry-perot optical fiber chamber demodulation method according to embodiments of the present invention；
Fig. 2 is the schematic diagram of according to embodiments of the present invention three kind Fourier transform spectrum；
Fig. 3 is the test system schematic diagram of Fabry-perot optical fiber cavity pressure sensor according to embodiments of the present invention；
Fig. 4 is the calibration curve figure of Fabry-perot optical fiber cavity pressure sensor according to embodiments of the present invention；And
Fig. 5 is the schematic diagram of Fabry-perot optical fiber chamber demodulating equipment according to embodiments of the present invention.
Specific embodiment
It should be noted that in the case where not conflicting, the feature in embodiment and embodiment in the application can phase
Mutually combine.Below with reference to the accompanying drawings and in conjunction with the embodiments describing the present invention in detail.
In order that those skilled in the art more fully understand application scheme, below in conjunction with the embodiment of the present application
Accompanying drawing, is clearly and completely described to the technical scheme in the embodiment of the present application, it is clear that described embodiment
The only embodiment of the application part, rather than the embodiment of whole.Based on the embodiment in the application, ability
The every other embodiment obtained under the premise of creative work is not made by domain those of ordinary skill, should all belong to
The scope of the application protection.
It should be noted that the description and claims of this application and the term " first " in above-mentioned accompanying drawing, "
Two " it is etc. for distinguishing similar object, without for describing specific order or precedence.It should be appreciated that this
The data that sample is used can be exchanged in the appropriate case, so as to embodiments herein described herein.Additionally, term
" comprising " and " having " and their any deformation, it is intended that covering is non-exclusive to be included, for example, comprising
The process of series of steps or unit, method, system, product or equipment are not necessarily limited to those steps clearly listed
Rapid or unit, but may include clearly not listing or intrinsic for these processes, method, product or equipment
Other steps or unit.
A kind of Fabry-perot optical fiber chamber demodulation method is embodiments provided, Fig. 1 is optical fiber according to embodiments of the present invention
The flow chart of Fa-Po cavity demodulation method.
As shown in figure 1, the method includes steps S102 to step S106：
Step S102：Obtain the spectral signal of interference optical fiber Fa-Po cavity.
After using wide range optical interface fiber Fa-Po cavity, spectrogrph can be passed through or (FBG) demodulator obtains the wide range interference of light
The spectral signal that Fabry-perot optical fiber chamber obtains, specifically, spectral signal can wide range optical interface fiber Fa-Po cavity obtain
Interference fringe.There is one-to-one relation in the chamber length in spectral signal and Fabry-perot optical fiber chamber, by demodulating the spectral signal
The chamber that Fabry-perot optical fiber chamber can be obtained is long.
Step S104：First demodulation is carried out using the discrete Fourier transform pairs spectral signal of Unequal distance, Fabry-perot optical fiber is obtained
The first chamber long value in chamber, wherein, the precision of the first chamber long value is more than first threshold, is more than in the precision of the first chamber long value
During first threshold, the cavity length demodulating in Fabry-perot optical fiber chamber is not in half-wave saltus step.
Pre-demodulating (the i.e. first demodulation) is carried out to spectral signal, obtains estimating chamber long value (i.e. the first chamber long value).It is pre- to solve
Tune be in order to prepare to essence demodulation is subsequently carried out to spectral signal, therefore, the required precision to estimating chamber long value is not high,
The precision for estimating chamber long value need to be only met more than first threshold, and first threshold can be to suppress to reach needed for half-wave saltus step
Accuracy value, specifically, first threshold could be arranged to quarter-wave.
Specifically, the embodiment of the present invention carries out the first demodulation using the discrete Fourier transform pairs spectral signal of Unequal distance, obtains
To the first chamber long value in Fabry-perot optical fiber chamber.
In field of signal processing, Fourier transformation is usually used in for time-domain signal being transformed to frequency-region signal being analyzed.In profit
When spectral signal analysis is carried out with Fourier transformation, as the spectral signal obtained from spectrogrph or (FBG) demodulator is typically closed
In wavelength equidistant (i.e. wave number or frequency unequal-interval), and directly using Fourier transformation generally need wave number or
The equidistant spectrum of frequency, the error of the chamber long value in the Fabry-perot optical fiber chamber for otherwise measuring are larger, typically will in prior art
Original spectrum signal differential to frequency it is equidistant after carry out Fourier transformation again, but so can be because changing original spectrum letter
Number and introduce error, additionally, when single frames spectral signature data point (i.e. the length of a frame spectrum number) is insufficient, adopting
The characteristic peak points obtained with fast Fourier transform are less, and peak information error is big.The embodiment of the present invention is according in Fu
The spectral information of original acquisition is carried out cross-correlation with corresponding sine and cosine sequence and obtains eigenvalue by leaf transformation original definition,
Fourier transform spectrum is obtained, thus obtained characteristic peak information point is more, peak error reduces.Preferably, the present invention
Embodiment carries out the first demodulation using the discrete Fourier transform pairs spectral signal of Unequal distance, wherein, Unequal distance direct computation of DFT
Leaf transformation refers to that the direct spectral signal to wave number or frequency unequal-interval carries out discrete Fourier transform, non-equally
Scattered Fourier transform formula is as follows：
Wherein, d represents the chamber long value in Fabry-perot optical fiber chamber, and N represents the spectral sequence columns of original spectral signal, λ tables
Oscillography is long, and x (λ) represents original spectral signal, and nDFT (d) is represented and obtained through Unequal distance discrete Fourier transform
Spectral signal, the referred to as discrete Fourier transform spectrum of Unequal distance.
Preferably, in order to the Unequal distance discrete fourier for effectively suppressing the above-mentioned Unequal distance discrete Fourier transforms of Jing to obtain becomes
The secondary lobe (other peaks i.e. by the discrete Fourier transform spectrum main peak of Unequal distance) exchanged cards with personal and family details when becoming sworn brothers, can become in above-mentioned discrete fourier
To original spectral signal adding window before changing.Preferably, the embodiment of the present invention is using the discrete Fourier transform pairs of Unequal distance
Spectral signal carries out the first demodulation, and before obtaining the first chamber long value in Fabry-perot optical fiber chamber, the method also includes：By spectrum
Signal is multiplied with window function, obtains the first spectral signal, is carried out using the discrete Fourier transform pairs spectral signal of Unequal distance
First demodulation includes：First demodulation is carried out using discrete the first spectral signal of Fourier transform pairs of Unequal distance, optical fiber is obtained
First chamber long value of Fa-Po cavity.
Spectral signal is multiplied with window function Ji Wei original spectral signal adding window, it is preferable that the window function is Bu Laike
Graceful window function.Unequal distance discrete Fourier transform is carried out to the spectral signal (i.e. the first spectral signal) after adding window can have
Effect obtains more accurately peak information, reduces pre-demodulating error.Fig. 2 is according to embodiments of the present invention three kind Fourier
The schematic diagram of conversion spectrum, as shown in Fig. 2 abscissa represents that chamber is long, unit is micron (um), and vertical coordinate is represented returns
One changes intensity, and (i.e. relative intensity, each intensity level divided by maximum intensity value, maximum is that 1), FFT represents quick
Fourier transform spectrum, nDFT represent Unequal distance discrete fourier change spectrum, and win-nDFT represents that adding window Unequal distance is discrete
Fourier's change spectrum (carries out the spectral signal that the change of Unequal distance discrete fourier is obtained) to the spectral signal after adding window.
Figure it is seen that adding window Unequal distance discrete fourier change spectrum signature peak information is significantly more than fast Fourier transform
Spectrum, and secondary lobe information is effectively inhibited, it is more accurate so as to obtain peak information, reduce pre-demodulating error.
Step S106：Second demodulation is carried out to spectral signal according to the first chamber long value, second chamber in Fabry-perot optical fiber chamber is obtained
Long value, wherein, the precision of the second chamber long value is more than the precision of the first chamber long value.
The embodiment of the present invention carries out pre-demodulating by above-mentioned steps and obtains estimating chamber long value that (i.e. the first chamber is long to spectral signal
Value), it is based further on this and estimates chamber long value carrying out spectral signal smart demodulation (i.e. second demodulation), obtains a precision
The chamber long value (i.e. the second chamber long value) in higher Fabry-perot optical fiber chamber, specifically, the embodiment of the present invention can be based on minimum
Square law carries out smart demodulation, or carries out smart demodulation to spectral signal based on cross-correlation method to spectral signal.This
Bright embodiment obtains the first chamber long value by the first demodulation is carried out to spectral signal, and ensures the precision of the first chamber long value
More than first threshold, carry out the second demodulation further to improve Fabry-perot optical fiber chamber to spectral signal based on the first chamber long value
The demodulation accuracy of chamber long value.The embodiment of the present invention is demodulated to spectral signal by demodulating with reference to the first demodulation second,
The precision of Fabry-perot optical fiber chamber cavity length demodulating is improve, and can effectively suppress half-wave hopping phenomenon, solve correlation technique
The not high problem of middle Fabry-perot optical fiber chamber cavity length demodulating precision.
Preferably, the second demodulation is carried out to spectral signal according to the first chamber long value, the second chamber for obtaining Fabry-perot optical fiber chamber is long
Value includes：Obtain the construction signal of spectral signal；The minimum of spectral signal and construction signal is calculated according to the first chamber long value
Mean square deviation estimated value；And the second chamber long value in Fabry-perot optical fiber chamber is determined according to nonlinear IEM model value.
The embodiment of the present invention can set up the construction signal of above-mentioned spectral signal according to spectral model, can be by calculating light
The mean square deviation estimated value of spectrum signal and construction signal determines the second chamber long value in Fabry-perot optical fiber chamber, specifically, spectral signal
It is long with the nonlinear IEM model value corresponding Fabry-perot optical fiber chamber chamber of construction signal, as the second chamber long value.Spectral signal
It is as follows with the computing formula of the mean square deviation estimated value of construction signal：
Wherein, d, d0Represent that the Fabry-perot optical fiber chamber chamber of the Fabry-perot optical fiber chamber chamber length and reality of construction is long respectively, x (n, d)
Represent construction signal, x (n, d0) spectral signal is represented, N represents the spectral series columns of spectral signal, and mse (d) is represented
Mean square deviation estimated value.
Alternatively, can be in the hope of spectral signal and the lowest mean square of construction signal by method of least square or cross-correlation method
Difference estimated value, i.e. second chamber long value long so as to obtain nonlinear IEM model value corresponding Fabry-perot optical fiber chamber chamber.But,
Amount of calculation when calculating spectral signal and construct the nonlinear IEM model value of signal using method of least square or cross-correlation method
Big speed is slow, less efficient.
In order to improve the efficiency of the nonlinear IEM model value for asking for spectral signal and construction signal, the embodiment of the present invention is adopted
With a kind of Algorithm for Solving spectral signal for gradually refining and the nonlinear IEM model value of construction signal, i.e., by gradually contracting
The mode of little scanning step solves nonlinear IEM model value.Preferably, the embodiment of the present invention is according to the first chamber long value meter
The nonlinear IEM model value for calculating spectral signal and construction signal includes：Obtain the first step-length；According to the first chamber long value and
First step-length determines sweep limitss；The second step-length is obtained, wherein, the second step-length is less than the first step-length；From sweep limitss
It is interior to carry out value every second step length, obtain multiple values；Respectively according to multiple exploitation spectral signals and construction letter
Number mean square deviation estimated value, obtain multiple mean square deviation estimated values；Obtain the mean square deviation of minimum in multiple mean square deviation estimated values
Estimated value is used as nonlinear IEM model value.
First step-length can be configured according to the error of the first chamber long value, for example, when the first chamber that the first demodulation is obtained
When the error of long value is 100nm, a length of 100nm of the first step can be set, it is also possible to arrange and must be slightly larger than 100nm,
For example, a length of 150nm of the first step is set.In practical situation, the error of the first chamber long value can be calculated by analyzing demodulation
Method is obtained, and the error of the chamber long value in the Fabry-perot optical fiber chamber obtained with same demodulating algorithm in same device is steady
Fixed, for example, the chamber for obtaining is demodulated in same Fabry-perot optical fiber cavity pressure sensor using discrete Fourier transform long
The error of value is essentially identical.After the first step-length is got, sweep limitss are determined according to the first step-length and the first chamber long value,
For example, a length of s1 of the first step, the first chamber long value are d1, then sweep limitss could be arranged to [d1-s1, d1+s1].The
Two step-lengths are used for the span for determining mean square deviation estimation scan, and the second step-length is less than the first step-length, specifically, can arrange
0.1 times of a length of first step-length of second step, it is also possible to arrange second step it is a length of other less than the first step-length value, for example,
The a length of 100nm of the first step, the second step-length could be arranged to 10nm, 20nm or 5nm etc..It is determined that second step
After long and sweep limitss, value is carried out every second step length from sweep limitss, obtain multiple values, for example, scanning
Scope is [d1-s1, d1+s1], and a length of s2 of second step is from the beginning of d1-s1, every from sweep limitss [d1-s1, d1+s1]
Value is carried out every the second step-length s2, multiple values are obtained, using each value in above-mentioned multiple values as above-mentioned
The chamber long value d of construction signal substitutes into above formula formula and calculates mean square deviation estimated value, can obtain multiple mean square deviation estimated values,
The above-mentioned multiple mean square deviation estimated values of comparison, using wherein minimum mean square deviation estimated value as nonlinear IEM model value.
Preferably, in order to further improve demodulation accuracy, the method also includes：Judge the second step-length whether more than second
Threshold value；When judging that second step length is more than Second Threshold, continue to obtain the 3rd step-length, wherein, the 3rd step-length is less than
Second step-length；Determine that the 3rd chamber in Fabry-perot optical fiber chamber is long according to mean square deviation estimated value minimum in multiple mean square deviation estimated values
Value；And sweep limitss are determined according to the 3rd chamber long value and the second step-length, wherein, from sweep limitss every the 3rd step
Length carries out value.
The Second Threshold of the embodiment of the present invention is usually arranged as the demodulation accuracy for finally needing to reach, and for example, arranges second
Threshold value is 0.1nm.If second step length is more than 0.1nm, illustrate that the chamber long value in currently available Fabry-perot optical fiber chamber does not reach
To default demodulation accuracy, now, continuation is demodulated using said method, redefines sweep limitss and scanning step
It is long, it is preferable that corresponding chamber long value (i.e. the 3rd chamber long value) is determined according to minimum mean square deviation estimated value obtained above,
And sweep limitss are redefined according to the second step-length and the 3rd chamber long value, and continue to obtain the 3rd step-length, wherein, the 3rd
Step-length is less than the second step-length, for example, can arrange that the 3rd step-length is the second step-length 0.1 times.Redefine from this
In sweep limitss, each the 3rd step-length value substitutes into above formula formula and calculates mean square deviation estimated value.Constantly proceed as described above,
Until scanning step reaches default demodulation accuracy, then the mean square deviation estimated value of the minimum for obtaining last time is used as most
Whole nonlinear IEM model value, and determine that the second chamber is long according to the nonlinear IEM model value.The embodiment of the present invention is led to
Cross, compared to existing
Method of least square solve nonlinear IEM model value and demodulate Fabry-perot optical fiber chamber chamber long value, ensureing higher demodulation essence
Under the conditions of degree, demodulation efficiency is greatly improved.
Application with example to the Fabry-perot optical fiber chamber demodulation method of the embodiment of the present invention below is illustrated：
Fabry-perot optical fiber chamber cavity length demodulating is the important component part of whole optical fiber pressure measuring system, directly affects and demodulates
The resolution of instrument, stability and certainty of measurement.Fig. 3 is Fabry-perot optical fiber cavity pressure sensing according to embodiments of the present invention
The test system schematic diagram of device.As shown in figure 3, by 1 interference light circulator 2 of light source, the collection spectrum of spectrogrph 3
Signal, specifically, spectra collection use Bayspec demodulation modules, its wave-length coverage be 1510～1590nm, single frames
Figures are 512, and spectral signal sends computer 4, the Fabry-perot optical fiber chamber according to embodiments of the present invention of computer 4 to
Demodulation method carries out cavity length demodulating according to spectral signal.In experimentation, by the spy of Fabry-perot optical fiber cavity pressure sensor
Head is placed in barochamber 5, barochamber 5 is pressurizeed using compression pump 6, and digital precision pressure meter 7 carries out pressure measurement.Figure
4 is the calibration curve figure of Fabry-perot optical fiber cavity pressure sensor according to embodiments of the present invention, in the diagram, abscissa table
Show actual pressure value, unit is MPa (MPa), vertical coordinate represents that demodulation chamber is long, and unit is micron (um), in 70MPa
In pressure limit, the curve linear is fitted degree of association up to 0.999999, and pressure measurement accuracy is 0.03%F.S.Thus may be used
To find out, the Fabry-perot optical fiber chamber demodulation method of the embodiment of the present invention has higher demodulation accuracy.
As can be seen from the above description, present invention achieves following technique effect：
Spectral signal of the embodiment of the present invention by acquisition interference optical fiber Fa-Po cavity；First demodulation is carried out to spectral signal,
The first chamber long value in Fabry-perot optical fiber chamber is obtained, wherein, the precision of the first chamber long value is more than first threshold；And according to
One chamber long value carries out the second demodulation to spectral signal, obtains the second chamber long value in Fabry-perot optical fiber chamber, wherein, the second chamber is long
The precision of value is demodulated to spectral signal by demodulating with reference to the first demodulation second more than the precision of the first chamber long value,
The precision of Fabry-perot optical fiber chamber cavity length demodulating is improve, and can effectively suppress half-wave hopping phenomenon, solve correlation technique
The not high problem of middle Fabry-perot optical fiber chamber cavity length demodulating precision.Additionally, being carried out using the discrete Fourier transform of adding window Unequal distance
The algorithm that pre-demodulating combination are gradually refined carries out smart demodulation, is capable of achieving the long Larger Dynamic in wide range optical interface fiber Fa-Po cavity chamber
Scope, high accuracy, the fast demodulation of absolute chamber length, it is possible to be prevented effectively from demodulation result and half-wave hopping phenomenon occur.
It should be noted that can be in such as one group of computer executable instructions the step of the flow process of accompanying drawing is illustrated
Perform in computer system, and, although show logical order in flow charts, but in some cases, can
With with different from the shown or described step of order execution herein.
A kind of another aspect according to embodiments of the present invention, there is provided Fabry-perot optical fiber chamber demodulating equipment, the Fabry-perot optical fiber chamber
Demodulating equipment can be used for the Fabry-perot optical fiber chamber demodulation method for performing the embodiment of the present invention, the Fiber Optic Sensor of the embodiment of the present invention
Amber chamber demodulation method can also be performed by the Fabry-perot optical fiber chamber demodulating equipment of the embodiment of the present invention.
Fig. 5 is the schematic diagram of Fabry-perot optical fiber chamber demodulating equipment according to embodiments of the present invention, as shown in figure 5, the device
Including：First acquisition unit 10, the first demodulating unit 20 and the second demodulating unit 30.
First acquisition unit 10, for obtaining the spectral signal of interference optical fiber Fa-Po cavity.
After using wide range optical interface fiber Fa-Po cavity, spectrogrph can be passed through or (FBG) demodulator obtains the wide range interference of light
The spectral signal that Fabry-perot optical fiber chamber obtains, specifically, spectral signal can wide range optical interface fiber Fa-Po cavity obtain
Interference fringe.There is one-to-one relation in the chamber length in spectral signal and Fabry-perot optical fiber chamber, by demodulating the spectral signal
The chamber that Fabry-perot optical fiber chamber can be obtained is long.
First demodulating unit 20, for carrying out the first demodulation using the discrete Fourier transform pairs spectral signal of Unequal distance, obtains
To the first chamber long value in Fabry-perot optical fiber chamber, wherein, the precision of the first chamber long value is more than first threshold, in the first chamber long value
Precision be more than first threshold when Fabry-perot optical fiber chamber cavity length demodulating be not in half-wave saltus step.
Pre-demodulating (the i.e. first demodulation) is carried out to spectral signal, obtains estimating chamber long value (i.e. the first chamber long value).It is pre- to solve
Tune be in order to prepare to essence demodulation is subsequently carried out to spectral signal, therefore, the required precision to estimating chamber long value is not high,
The precision for estimating chamber long value need to be only met more than first threshold, and first threshold can be to suppress the required of half-wave saltus step to reach
Accuracy value, specifically, first threshold could be arranged to quarter-wave.
The embodiment of the present invention according to Fourier transformation original definition, by the spectral information of original acquisition and corresponding sine and cosine sequence
Row carry out cross-correlation and obtain eigenvalue, that is, obtain Fourier transform spectrum, and thus obtained characteristic peak information point is more, peak value
Error reduces.Above-mentioned Unequal distance discrete Fourier transform is directly entered to the spectral signal of wave number or frequency unequal-interval
Row discrete Fourier transform.
Preferably, the device also includes：Arithmetic element, for spectral signal is multiplied with window function, obtains the first light
Spectrum signal, the first demodulating unit for carrying out the first demodulation using discrete the first spectral signal of Fourier transform pairs of Unequal distance,
Obtain the first chamber long value in Fabry-perot optical fiber chamber.
Spectral signal is multiplied with window function Ji Wei original spectral signal adding window, it is preferable that the window function is Bu Laike
Graceful window function.Unequal distance discrete Fourier transform is carried out to the spectral signal (i.e. the first spectral signal) after adding window can have
Effect obtains more accurately peak information, reduces pre-demodulating error.
Second demodulating unit 30, for carrying out second demodulation to spectral signal according to the first chamber long value, obtains Fabry-perot optical fiber
The second chamber long value in chamber, wherein, the precision of the second chamber long value is more than the precision of the first chamber long value.
Specifically, the embodiment of the present invention can carry out second demodulation to spectral signal based on method of least square, it is also possible to
It is that the second demodulation is carried out to spectral signal based on cross-correlation method.
The embodiment of the present invention obtains the spectral signal of interference optical fiber Fa-Po cavity by first acquisition unit 10；First demodulation is single
First 20 pairs of spectral signals carry out the first demodulation, obtain the first chamber long value in Fabry-perot optical fiber chamber, wherein, the first chamber long value
Precision is more than first threshold；And second demodulating unit 30 second demodulation is carried out to spectral signal according to the first chamber long value,
The second chamber long value in Fabry-perot optical fiber chamber is obtained, wherein, the precision of the second chamber long value is more than the precision of the first chamber long value.It is logical
Cross and spectral signal is demodulated with reference to the first demodulating unit 20 and the second demodulating unit 30, improve Fabry-perot optical fiber chamber
The precision of cavity length demodulating, and can effectively suppress half-wave hopping phenomenon, in solving correlation technique, Fabry-perot optical fiber chamber chamber is long
The not high problem of demodulation accuracy.
Preferably, the second demodulating unit 30 includes：Acquisition module, for obtaining the construction signal of spectral signal；Calculate
Module, for the nonlinear IEM model value of spectral signal and construction signal is calculated according to the first chamber long value；And determine
Module, for the second chamber long value in Fabry-perot optical fiber chamber is determined according to nonlinear IEM model value.
Preferably, it is in order to improve the efficiency of the nonlinear IEM model value for asking for spectral signal and construction signal, of the invention
The computing module of embodiment includes：First acquisition submodule, for obtaining the first step-length；First determination sub-module, uses
In determining sweep limitss according to the first chamber long value and the first step-length；Second acquisition submodule, for obtaining the second step-length,
Wherein, the second step-length is less than the first step-length；Value submodule, for being taken every second step length from sweep limitss
Value, obtains multiple values；Calculating sub module, for respectively according to multiple exploitation spectral signals and construction signal
Mean square deviation estimated value, obtains multiple mean square deviation estimated values；In obtaining multiple mean square deviation estimated values, minimum mean square deviation is estimated
Value is used as nonlinear IEM model value.
Preferably, the device also includes：Judging unit, for judging the second step-length whether more than Second Threshold；Second
Acquiring unit, for when judging that second step length is more than Second Threshold, obtaining the 3rd step-length, wherein, the 3rd step-length
Less than the second step-length；First determining unit, it is true for the mean square deviation estimated value according to minimum in multiple mean square deviation estimated values
Determine the 3rd chamber long value in Fabry-perot optical fiber chamber；And second determining unit, for true according to the 3rd chamber long value and the second step-length
Determine sweep limitss, wherein, value is carried out every the 3rd step-length from sweep limitss.
The embodiment of the present invention additionally provides a kind of Fabry-perot optical fiber interferometer, and the Fabry-perot optical fiber interferometer includes above-mentioned any one
Individual Fabry-perot optical fiber chamber demodulating equipment.
Obviously, those skilled in the art should be understood that each module or each step of the above-mentioned present invention can be with general
Computing device realizing, they can be concentrated on single computing device, or are distributed in multiple computing device institutes
On the network of composition, alternatively, they can be realized with the executable program code of computing device, it is thus possible to
It is stored in storage device by computing device performing, or they is fabricated to each integrated circuit die respectively
Block, or the multiple modules or step in them are fabricated to single integrated circuit module to realize.So, the present invention
It is not restricted to any specific hardware and software to combine.
The preferred embodiments of the present invention are the foregoing is only, the present invention is not limited to, for the skill of this area
For art personnel, the present invention can have various modifications and variations.It is all within the spirit and principles in the present invention, made
Any modification, equivalent substitution and improvements etc., should be included within the scope of the present invention.
Claims (11)
1. a kind of Fabry-perot optical fiber chamber demodulation method, it is characterised in that include：
Obtain the spectral signal for interfering the Fabry-perot optical fiber chamber；
First demodulation is carried out using spectral signal described in the discrete Fourier transform pairs of Unequal distance, the Fiber Optic Sensor is obtained
The first chamber long value in amber chamber, wherein, the precision of first chamber long value is more than first threshold, in first chamber
When the precision of long value is more than the first threshold, the cavity length demodulating in the Fabry-perot optical fiber chamber is not in half-wave saltus step；
And
Second demodulation is carried out to the spectral signal according to first chamber long value, the Fabry-perot optical fiber chamber is obtained
Second chamber long value, wherein, the precision of second chamber long value is more than the precision of first chamber long value.
2. Fabry-perot optical fiber chamber according to claim 1 demodulation method, it is characterised in that
First demodulation is being carried out using spectral signal described in the discrete Fourier transform pairs of Unequal distance, the optical fiber is being obtained
Before first chamber long value of Fa-Po cavity, methods described also includes：The spectral signal is multiplied with window function, is obtained
To the first spectral signal,
The first demodulation is carried out using spectral signal described in the discrete Fourier transform pairs of Unequal distance includes：Using Unequal distance
Discrete Fourier transform carries out the first demodulation to first spectral signal, obtains the first of the Fabry-perot optical fiber chamber
Chamber long value.
3. Fabry-perot optical fiber chamber according to claim 1 demodulation method, it is characterised in that according to first chamber long value
Second demodulation is carried out to the spectral signal, the second chamber long value for obtaining the Fabry-perot optical fiber chamber includes：
Obtain the construction signal of the spectral signal；
The nonlinear IEM model of the spectral signal and the construction signal is calculated according to first chamber long value
Value；And
The second chamber long value in the Fabry-perot optical fiber chamber is determined according to the nonlinear IEM model value.
4. Fabry-perot optical fiber chamber according to claim 3 demodulation method, it is characterised in that according to first chamber long value
The nonlinear IEM model value for calculating the spectral signal and the construction signal includes：
Obtain the first step-length；
Sweep limitss are determined according to first chamber long value and first step-length；
The second step-length is obtained, wherein, second step-length is less than first step-length；
Value is carried out every the second step length from the sweep limitss, multiple values are obtained；
Respectively according to the plurality of exploitation spectral signal and the construction signal mean square deviation estimated value,
Obtain multiple mean square deviation estimated values；
The mean square deviation estimated value of minimum in the plurality of mean square deviation estimated value is obtained as nonlinear IEM model value.
5. Fabry-perot optical fiber chamber according to claim 4 demodulation method, it is characterised in that methods described also includes：
Judge second step-length whether more than Second Threshold；
When judging that the second step length is more than the Second Threshold, continue to obtain the 3rd step-length, wherein, institute
The 3rd step-length is stated less than second step-length；
The of the Fabry-perot optical fiber chamber is determined according to mean square deviation estimated value minimum in the plurality of mean square deviation estimated value
Three chamber long values；And
The sweep limitss are determined according to the 3rd chamber long value and second step-length, wherein, from the scanning
In the range of carry out value every the 3rd step-length.
6. a kind of Fabry-perot optical fiber chamber demodulating equipment, it is characterised in that include：
First acquisition unit, for obtaining the spectral signal for interfering the Fabry-perot optical fiber chamber；
First demodulating unit, for carrying out the first solution using spectral signal described in the discrete Fourier transform pairs of Unequal distance
Adjust, obtain the first chamber long value in the Fabry-perot optical fiber chamber, wherein, the precision of first chamber long value is more than first
Threshold value, the cavity length demodulating in the Fabry-perot optical fiber chamber when precision of long value is more than the first threshold in first chamber
It is not in half-wave saltus step；And
Second demodulating unit, for carrying out second demodulation to the spectral signal according to first chamber long value, obtains
To the second chamber long value in the Fabry-perot optical fiber chamber, wherein, the precision of second chamber long value is more than first chamber
The precision of long value.
7. Fabry-perot optical fiber chamber according to claim 6 demodulating equipment, it is characterised in that
Described device also includes：Arithmetic element, for the spectral signal is multiplied with window function, obtains first
Spectral signal,
First demodulating unit is for being carried out using the first spectral signal described in the discrete Fourier transform pairs of Unequal distance
First demodulation, obtains the first chamber long value in the Fabry-perot optical fiber chamber.
8. Fabry-perot optical fiber chamber according to claim 6 demodulating equipment, it is characterised in that the second demodulating unit bag
Include：
Acquisition module, for obtaining the construction signal of the spectral signal；
Computing module, for calculating the spectral signal and the construction signal most according to first chamber long value
Little mean square deviation estimated value；And
Determining module, the second chamber for the Fabry-perot optical fiber chamber is determined according to the nonlinear IEM model value are long
Value.
9. Fabry-perot optical fiber chamber according to claim 8 demodulating equipment, it is characterised in that the computing module includes：
First acquisition submodule, for obtaining the first step-length；
First determination sub-module, for determining sweep limitss according to first chamber long value and first step-length；
Second acquisition submodule, for obtaining the second step-length, wherein, second step-length is less than the first step
It is long；
Value submodule, for carrying out value every the second step length from the sweep limitss, obtains multiple
Value；
Calculating sub module, for spectral signal and the construction signal according to the plurality of exploitation respectively
Mean square deviation estimated value, obtain multiple mean square deviation estimated values；
The mean square deviation estimated value of minimum in the plurality of mean square deviation estimated value is obtained as nonlinear IEM model value.
10. Fabry-perot optical fiber chamber according to claim 9 demodulating equipment, it is characterised in that described device also includes：
Judging unit, for judging second step-length whether more than Second Threshold；
Second acquisition unit, for when judging that the second step length is more than the Second Threshold, obtaining the 3rd
Step-length, wherein, the 3rd step-length is less than second step-length；
First determining unit, for being determined according to mean square deviation estimated value minimum in the plurality of mean square deviation estimated value
The 3rd chamber long value in the Fabry-perot optical fiber chamber；And
Second determining unit, for determining the sweep limitss according to the 3rd chamber long value and second step-length,
Wherein, value is carried out every the 3rd step-length from the sweep limitss.
11. a kind of Fabry-perot optical fiber interferometers, it is characterised in that including the Fiber Optic Sensor any one of claim 6 to 10
Amber chamber demodulating equipment.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510740349.3A CN106643835B (en) | 2015-11-03 | 2015-11-03 | Fabry-perot optical fiber chamber demodulation method and device and Fabry-perot optical fiber interferometer |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510740349.3A CN106643835B (en) | 2015-11-03 | 2015-11-03 | Fabry-perot optical fiber chamber demodulation method and device and Fabry-perot optical fiber interferometer |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106643835A true CN106643835A (en) | 2017-05-10 |
| CN106643835B CN106643835B (en) | 2019-07-05 |
Family
ID=58850935
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510740349.3A Active CN106643835B (en) | 2015-11-03 | 2015-11-03 | Fabry-perot optical fiber chamber demodulation method and device and Fabry-perot optical fiber interferometer |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106643835B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108534811A (en) * | 2018-04-13 | 2018-09-14 | 西安工业大学 | A kind of cavity length demodulating algorithm of short cavity optical fiber Fabry-Perot sensor |
| CN108955734A (en) * | 2018-06-08 | 2018-12-07 | 武汉理工大学 | A kind of cavity length demodulating method of fiber F-P temperature/pressure compound sensor |
| CN110595516A (en) * | 2019-09-18 | 2019-12-20 | 华中科技大学 | FPI cavity length demodulation method and system |
| CN111707303A (en) * | 2020-07-17 | 2020-09-25 | 武汉理工大学 | Variable-step-length rapid high-precision signal demodulation method for optical fiber F-P sensor |
| CN112050942A (en) * | 2020-07-30 | 2020-12-08 | 东北大学 | Optical fiber interference spectrum cavity length correction method based on phase compensation |
| CN112081587A (en) * | 2020-10-23 | 2020-12-15 | 西南石油大学 | Horizontal well output profile calculation method based on optical fiber noise data |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070086018A1 (en) * | 2005-10-14 | 2007-04-19 | Board Of Regents, The University Of Texas System | Broadband cavity spectrometer apparatus and method for determining the path length of an optical structure |
| CN101158602A (en) * | 2007-11-21 | 2008-04-09 | 北京理工大学 | Optical fibre fourier transformation white light interferometric method |
| CN101958729A (en) * | 2010-09-30 | 2011-01-26 | 哈尔滨工程大学 | Quick chirp modulation interference detection and suppression method for direct-sequence spread-spectrum communication system |
| CN102607612A (en) * | 2012-03-22 | 2012-07-25 | 天津大学 | Low-coherence interference demodulation method based on monochrome frequency absolute phase |
-
2015
- 2015-11-03 CN CN201510740349.3A patent/CN106643835B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070086018A1 (en) * | 2005-10-14 | 2007-04-19 | Board Of Regents, The University Of Texas System | Broadband cavity spectrometer apparatus and method for determining the path length of an optical structure |
| CN101158602A (en) * | 2007-11-21 | 2008-04-09 | 北京理工大学 | Optical fibre fourier transformation white light interferometric method |
| CN101958729A (en) * | 2010-09-30 | 2011-01-26 | 哈尔滨工程大学 | Quick chirp modulation interference detection and suppression method for direct-sequence spread-spectrum communication system |
| CN102607612A (en) * | 2012-03-22 | 2012-07-25 | 天津大学 | Low-coherence interference demodulation method based on monochrome frequency absolute phase |
Non-Patent Citations (1)
| Title |
|---|
| 周新磊: "高精度光纤传感系统及其在油田测井中的应用", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108534811A (en) * | 2018-04-13 | 2018-09-14 | 西安工业大学 | A kind of cavity length demodulating algorithm of short cavity optical fiber Fabry-Perot sensor |
| CN108534811B (en) * | 2018-04-13 | 2020-11-27 | 西安工业大学 | Cavity length demodulation algorithm of short-cavity optical fiber Fabry-Perot sensor |
| CN108955734A (en) * | 2018-06-08 | 2018-12-07 | 武汉理工大学 | A kind of cavity length demodulating method of fiber F-P temperature/pressure compound sensor |
| CN110595516A (en) * | 2019-09-18 | 2019-12-20 | 华中科技大学 | FPI cavity length demodulation method and system |
| CN110595516B (en) * | 2019-09-18 | 2020-12-18 | 华中科技大学 | FPI cavity length demodulation method and system |
| CN111707303A (en) * | 2020-07-17 | 2020-09-25 | 武汉理工大学 | Variable-step-length rapid high-precision signal demodulation method for optical fiber F-P sensor |
| CN112050942A (en) * | 2020-07-30 | 2020-12-08 | 东北大学 | Optical fiber interference spectrum cavity length correction method based on phase compensation |
| CN112050942B (en) * | 2020-07-30 | 2021-06-15 | 东北大学 | Optical fiber interference spectrum cavity length correction method based on phase compensation |
| CN112081587A (en) * | 2020-10-23 | 2020-12-15 | 西南石油大学 | Horizontal well output profile calculation method based on optical fiber noise data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106643835B (en) | 2019-07-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106643835A (en) | Fiber fabry-perot cavity demodulation method, device thereof and fiber fabry-perot interferometer | |
| CN104634460B (en) | Multi-peak self-adaption accurate peak searching method for distributed FBG (Fiber Bragg Grating) sensing network | |
| JP6817209B2 (en) | Alignment of measured interference data of optical fiber with reference interference data of optical fiber | |
| CN108955734B (en) | Cavity length demodulation method of optical fiber F-P temperature/pressure composite sensor | |
| CN108534811B (en) | Cavity length demodulation algorithm of short-cavity optical fiber Fabry-Perot sensor | |
| CN105973282A (en) | Fiber F-P sensor cavity length wavelet phase extraction demodulation method | |
| CN109916434B (en) | Optical fiber Fabry-Perot sensor cavity length multiple frequency related demodulation method | |
| US11385130B2 (en) | Multi-spectral feature sensing techniques and sensors for determining a phase shift by comparing a measured narrowband response to a known wideband response | |
| CN105841723A (en) | FBG sensor demodulation peak-searching method based on wavelet singularity detection | |
| CN108896274B (en) | Distributed optical fiber strain demodulation method based on subset window length optimization algorithm | |
| US10445403B2 (en) | Method and system for tracking the centre of a peak from a plurality of sample points in an optical system | |
| CN112050942B (en) | Optical fiber interference spectrum cavity length correction method based on phase compensation | |
| CN104061875A (en) | High-precision fiber bragg grating demodulation method based on Hilbert transformation and bispectrum estimation | |
| Zhang et al. | Segmented noise reduction based on Brillouin-spectrum-partition in Brillouin optical time domain sensors | |
| CN108645431B (en) | Fitting peak searching method for cavity length correlation demodulation of optical fiber Fabry-Perot sensor | |
| CN115597641A (en) | Dynamic cavity length demodulation method and system of optical fiber Fabry-Perot sensor | |
| CN107250788B (en) | Noise level estimation method and determination data processing unit | |
| Gerner et al. | Improving self‐noise estimates of broadband seismometers by 3D trace rotation | |
| Harasim et al. | Accuracy improvement of Fiber Bragg Grating peak wavelength demodulation using wavelet transform and various center wavelength detection algorithms | |
| JP2007322437A (en) | Device and method for performing measurement processing of fiber bragg grating element reflection light wavelength | |
| CN112986180B (en) | Spectrum type gas sensing data processing method and system | |
| CN114518182B (en) | Method and system for simultaneously extracting temperature and strain information in brillouin scattering spectrum image | |
| CN117168515A (en) | High-spatial-resolution OFDR data processing method based on blind source separation | |
| Vila et al. | Improvements of Prakash's method for analyzing ITC peaks | |
| Rao et al. | A novel algorithm for 2-dimensional frequency estimation |
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 |