CN106839972A - A kind of interference signal processing method of full light fiber white light interference instrument - Google Patents
A kind of interference signal processing method of full light fiber white light interference instrument Download PDFInfo
- Publication number
- CN106839972A CN106839972A CN201710020296.7A CN201710020296A CN106839972A CN 106839972 A CN106839972 A CN 106839972A CN 201710020296 A CN201710020296 A CN 201710020296A CN 106839972 A CN106839972 A CN 106839972A
- Authority
- CN
- China
- Prior art keywords
- data
- interference signal
- direct current
- interference
- group
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02055—Reduction or prevention of errors; Testing; Calibration
- G01B9/02062—Active error reduction, i.e. varying with time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/0209—Low-coherence interferometers
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Instruments For Measurement Of Length By Optical Means (AREA)
Abstract
A kind of interference signal processing method of full light fiber white light interference instrument, is related to optical interference signal data processing method, belongs to optical measurement and technical field of optical fiber sensing.The present invention provide method be:Assuming that containing N number of data altogether by the interference signal that full light fiber white light interference instrument is collected, N is positive integer;IiRepresent i-th numerical value of data, i ∈ (0, N 1);Data processing need to be as shown in Figure 1 by 4 steps, algorithm flow;The original interference signal that this method can mutually fluctuate larger to direct current is processed, and is that the calculating of follow-up associated arguments lays the foundation so as to obtain stable interference signal.
Description
Technical field
The present invention relates to optical interference signal data processing method, belong to optical measurement and technical field of optical fiber sensing, especially
It is related to a kind of white light interference signal processing method based on all -fiber Michelson's interferometer.
Background technology
Measuring system based on interference has turned into one of system the most accurate during physical quantity is detected.Conventional interference system
The interference system being made up of the mixing of Michelson's interferometer, Mach Zehnder interferometer, sagnac interferometer and each interferometer.
Michelson's interferometer is because with simple structure, fringe contrast is good, signal to noise ratio is high, the calculating of striped counting and tested displacement
The advantages of relation is simple and widely studied and applied.Michelson's interferometer can be divided into spatial light according to its structure type
Two kinds of Michelson's interferometer and all -fiber formula Michelson's interferometer.Spatial light Michelson's interferometer optical path adjusting is more multiple
Miscellaneous, the stability of a system and mobility are poor.Although all -fiber formula Michelson's interferometer is in the stability of a system, mobility and light
Road regulation complexity on be better than spatial light Michelson's interferometer, but haves the shortcomings that to go wrong do not allow easy-to-search and adjust.
White light spectral range is wide, continuous, coherence length is short, only optical path difference very little when, can just interfere.Work as optical path difference
When being zero, each spectral line dual-beam is completely superposed in white-light spectrum, and various wavelength lights are overlapped, the maximum center zero of degree of being contrasted
Level striped, i.e., optimal interference position is realized to parameter high-acruracy survey by interference.White light source can be stepped as all -fiber
The light source of Ke Erxun interferometers.
The quality of interference signal directly decides the precision of measurement.Because optical fiber type Michelson's interferometer and spatial light are stepped
Ke Erxun interferometers are all to compensate optical path difference by mechanical scanning arm, and the mechanical oscillation for scanning guide rail will reduce the precision of measurement.
Usual mechanical oscillation are prevalent in all kinds of instrument and meters and equipment, and the performance that vibration directly influences all kinds of instrument and meters refers to
Mark, the influence for reducing vibration signal turns into one of key of raising instrument and meter and equipment performance.Such as Publication No.
The signal of patent of invention " the abnormal interference elimination method of mechanical oscillation signal detection " successively to collecting of 103424183B does 1 rank
After the computings such as difference, high-order smoothing processing, reconstruction signal, exclusion impulse disturbances, it is possible to decrease the interference of vibration signal, but the method
Calculating process is more complicated.A kind of and for example patent of invention " vibration for eliminating velocity hydrophone of Publication No. 102494758A
The device of noise " is counteracted to underwater platform by using vibration monitoring sensor, signal conditioning circuit and differential amplifier circuit
Interference of the vibration noise to velocity hydrophone, which thereby enhance the accuracy of measurement of velocity hydrophone, but the method is applicable model
It is with limit.
The content of the invention
The main object of the present invention is directed to because all -fiber formula white light Michelson's interferometer scans the non-linear of guide rail
Or interference signal direct current is mutually unstable caused by vibration, so as to influence to interfere the certainty of measurement of primary maximum and secondary maximum, there is provided one
Full light fiber white light interference instrument experimental provision is planted, a kind of data processing method is proposed.This method can mutually fluctuate larger original to direct current
Beginning interference signal is processed, and is that the calculating of follow-up associated arguments lays the foundation so as to obtain stable interference signal.
The present invention proposes a kind of interference signal processing method for full light fiber white light interference instrument.Assuming that white by all -fiber
The interference signal that optical interferometer is collected contains N number of data altogether (N is positive integer).IiRepresent i-th numerical value of data, i ∈ (0,
N-1).Data processing need to be as shown in Figure 1 by following 4 steps, algorithm flow.
Step 1:The data that will be collected are grouped with T adjacent data as one group, it is assumed that total G groups (0,1,2 ...
G-1), adding window is carried out to data equivalent in time domain.The initial data of g (G-2 of g=0,1,2 ...) group data is to collect
N (since 0 count) individual initial data the gT data, last data is (g+1) T- for the N number of data for collecting
1 data, the initial data of G-1 group data is (G-1) T data, and end data are the N-1 data.
Step 2:Calculate g (G-1 of g=0,1,2 ...) group data direct current phase numerical value beDue to every group of collection T number
According to the distance of scanning guide rail movement is very short, and negligible guide rail is non-linear to interference signal direct current phase in the segment distance
Influence, g groups data interference direct current is mutually represented by:
It is the average of g group data, represents this group of size of data direct current phase.
Step 3:Recalculate interference signal direct current phase.Interference signal direct current is mutually represented by:
That is M is the average of institute's gathered data, represents interference signal direct current phase.
Step 4:Reconstitution experiments data.
Represent the experimental data after reconstruct.It is equal to raw experimental data IiSubtract the direct current phase of group where dataAgain
Plus overall interference signal direct current phase M.
Raw experimental data by after above-mentioned 4 step process, its stationarity must to greatly raising, effect is significant,
For the calculating of the associated arguments such as subsequent coupling intensity is laid a good foundation.
Experimental data is processed based on the data processing algorithm flow chart shown in Fig. 1.Experiment acquires 5013000 altogether
Individual data, as shown in Figure 3.Multiplying window is carried out to initial data by window width of T=125 in 1 the step of technical scheme, is used
The average of the data in different windowRepresent the DC component of signal in the window.The direct current phase of interference signal is calculated in step 3
Equal to M=0.0748.Step 4 is reconstructed to experimental dataExperimental data after reconstruct is as shown in Figure 4.Figure
It is that partial enlargement is carried out to Fig. 3 experimental datas shown in 5, is 310 meters of optical fiber incidence end Coupling point raw experimental datas, by above-mentioned
Method treatment after experimental data it is as shown in Figure 6.From Fig. 4 and Fig. 6, the data processing algorithm greatly improves data
Stationarity, effect is significant.
A kind of interference signal processing method for full light fiber white light interference instrument proposed by the present invention, by time-domain signal
Multiplying window is carried out, the DC component of different window internal interference signal is calculated.Raw experimental data is subtracted in time window where it
Interference DC component, along with overall interference direct current phase.Interference direct current caused by the non-linear or vibration of scanning guide rail can be eliminated
Influence of the signal fluctuation to signal, greatly improves the stationarity of data, and effect is obvious.
Brief description of the drawings
Fig. 1 data processing algorithm flow charts
Fig. 2 full light fiber white light interference instrument experimental provision structure charts
Fig. 3 optical fiber scannings gained raw experimental data
Fig. 4 optical fiber scannings data after treatment
The experimental data of Fig. 5 optical fiber incidence end Coupling point before processings
Experimental data after the treatment of Fig. 6 optical fiber incidence ends Coupling point
Specific embodiment
Embodiment 1:
Fig. 2 is full light fiber white light interference instrument experimental provision.The device is by light source module, the polarizer, polarization maintaining optical fibre, 45 degree of inspections
Inclined device module, Michelson's interferometer module, photodetector, data collecting card, 8 parts of computer are constituted.
The light that SLD white light sources send is changed into line polarisation after the polarizer, and the line polarisation is aligned with polarization maintaining optical fibre fast axle, enters
It is mapped in polarization maintaining optical fibre;When there is external force to act on polarization maintaining optical fibre, line polarisation is transmitted to external force and makees use, has polarization coupled
Phenomenon occurs, and excites the portion of energy crosstalk of mould to slow axis along fast axis propagation, forms coupled mode;Due to fiber mode birefringence
ΔnbPresence, excite mould and coupled mode to generate certain optical path difference at fiber exit end;By after analyzer, inciding wheat
Ke Erxun interferometers.Photodetector converts light signals into electric signal, then carries out signal acquisition by data collecting card, is input to
Computer carries out signal transacting.In fig. 2, light source module is made up of SLD light sources (centre wavelength is 1550nm) and isolator.Rise
Inclined device is the polarizer of fiber type.Polarization maintaining optical fibre to be measured works in 1550nm, and length is 310m.45 degree of analyzing modules are by polarization-maintaining
Wire jumper (centre is blocked, fast axle and slow axis welding again in angle of 45 degrees) and analyzer are concatenated and formed.Coupler is 2 × 2 types
Fiber coupler.Variable optical delay line is produced by AM General photoelectricity, its model MDL -002, works in 1550nm, is postponed
Scope is 0~560ps.Two speculums are faraday's rotating mirror of fibre-optical splice.Photodetector is Thorlabs's
PDA10CS-EC types.Data collecting card uses NI USB6251.Software is write using Labview.
The interference signal processing method for full light fiber white light interference instrument of the present embodiment is:Assuming that passing through full light fiber white light
The interference signal that interferometer is collected contains N number of data altogether (N is positive integer).IiRepresent i-th numerical value of data, i ∈ (0, N-
1).Data processing need to be as shown in Figure 1 by following 4 steps, algorithm flow.
Step 1:The data that will be collected are grouped with T adjacent data as one group, it is assumed that total G groups (0,1,2 ...
G-1), adding window is carried out to data equivalent in time domain.The initial data of g (G-2 of g=0,1,2 ...) group data is to collect
N (since 0 count) individual initial data the gT data, last data is (g+1) T- for the N number of data for collecting
1 data, the initial data of G-1 group data is (G-1) T data, and end data are the N-1 data.
Step 2:Calculate g (G-1 of g=0,1,2 ...) group data direct current phase numerical value beDue to every group of collection T number
According to the distance of scanning guide rail movement is very short, and negligible guide rail is non-linear to interference signal direct current phase in the segment distance
Influence, g groups data interference direct current is mutually represented by:
It is the average of g group data, represents this group of size of data direct current phase.
Step 3:Recalculate interference signal direct current phase.Interference signal direct current is mutually represented by:
That is M is the average of institute's gathered data, represents interference signal direct current phase.
Step 4:Reconstitution experiments data.
Represent the experimental data after reconstruct.It is equal to raw experimental data IiSubtract the direct current phase of group where dataAgain
Plus overall interference signal direct current phase M.
Experimental data is processed based on the data processing algorithm flow chart shown in Fig. 1.Experiment acquires 5013000 altogether
Individual data, as shown in Figure 3.Multiplying window is carried out to initial data by window width of T=125 in 1 the step of technical scheme, is used
The average of the data in different windowRepresent the DC component of signal in the window.The direct current phase of interference signal is calculated in step 3
Equal to M=0.0748.Step 4 is reconstructed to experimental dataExperimental data after reconstruct is as shown in Figure 4.Fig. 5
Shown is to carry out partial enlargement to Fig. 3 experimental datas, is 310 meters of optical fiber incidence end Coupling point raw experimental datas, by above-mentioned
Method treatment after experimental data it is as shown in Figure 6.From Fig. 4 and Fig. 6, the data processing algorithm greatly improves data
Stationarity, effect is significant.
Claims (1)
1. a kind of interference signal processing method of full light fiber white light interference instrument, it is characterised in that the method is:Assuming that passing through full light
The interference signal that fiber white light interference instrument is collected contains N number of data altogether, and N is positive integer;IiRepresent i-th numerical value of data, i ∈
(0, N-1);Data processing need to be as shown in Figure 1 by following 4 steps, algorithm flow;
Step 1:The data that will be collected are grouped with T adjacent data as one group, it is assumed that total G groups (0,1,2 ... G-
1) adding window, is carried out to data equivalent in time domain;The initial data of g (G-2 of g=0,1,2 ...) group data is the N for collecting
The gT data of (being counted since 0) individual initial data, last data are (g+1) T-1 of the N number of data for collecting
Data, the initial data of G-1 group data is (G-1) T data, and end data are the N-1 data;
Step 2:Calculate g (G-1 of g=0,1,2 ...) group data direct current phase numerical value beDue to T data of every group of collection, sweep
The distance for retouching guide rail movement is very short, the non-linear influence to interference signal direct current phase of negligible guide rail in the segment distance,
G groups data interference direct current is mutually represented by:
It is the average of g group data, represents this group of size of data direct current phase;
Step 3:Recalculate interference signal direct current phase;Interference signal direct current is mutually represented by:
That is M is the average of institute's gathered data, represents interference signal direct current phase;
Step 4:Reconstitution experiments data;
Represent the experimental data after reconstruct.It is equal to raw experimental data IiSubtract the direct current phase of group where dataAdd
Overall interference signal direct current phase M;
Raw experimental data is by after above-mentioned 4 step process, its stationarity will must be improved to great, effect is significant, after being
The calculating of the associated arguments such as continuous stiffness of coupling is laid a good foundation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710020296.7A CN106839972B (en) | 2017-01-11 | 2017-01-11 | A kind of interference signal processing method of full light fiber white light interference instrument |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710020296.7A CN106839972B (en) | 2017-01-11 | 2017-01-11 | A kind of interference signal processing method of full light fiber white light interference instrument |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106839972A true CN106839972A (en) | 2017-06-13 |
CN106839972B CN106839972B (en) | 2019-05-10 |
Family
ID=59123194
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710020296.7A Active CN106839972B (en) | 2017-01-11 | 2017-01-11 | A kind of interference signal processing method of full light fiber white light interference instrument |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106839972B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113739708A (en) * | 2021-08-20 | 2021-12-03 | 哈尔滨工程大学 | Double-path white light interference differential measurement device and method |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4889986A (en) * | 1988-08-18 | 1989-12-26 | The United States Of America As Represented By The Secretary Of The Navy | Serial interferometric fiber-optic sensor array |
JP3630581B2 (en) * | 1999-03-04 | 2005-03-16 | アルプス電気株式会社 | Spread modulation signal receiver |
CN102045120A (en) * | 2010-10-29 | 2011-05-04 | 成都九洲电子信息系统有限责任公司 | Vibration signal identification method for optical fiber perimeter system |
CN102494758A (en) * | 2011-12-29 | 2012-06-13 | 中国船舶重工集团公司第七一○研究所 | Apparatus for eliminating interference of vibration noise on velocity hydrophone |
CN102636196A (en) * | 2012-04-09 | 2012-08-15 | 天津大学 | Distributed disturbance sensor on basis of Rayleigh scattering spectrum related coefficient and demodulating method thereof |
CN102832620A (en) * | 2012-08-31 | 2012-12-19 | 天津理工大学 | System and method for harmonic detection and control of APF (active power filter) based on windowed all-phase FFT (fast Fourier transform) |
CN103424183A (en) * | 2013-08-06 | 2013-12-04 | 常州大学 | Method for eliminating abnormal interference on detection for mechanical vibration signals |
CN103759804A (en) * | 2014-01-23 | 2014-04-30 | 安徽大学 | Method and device for optical fiber white light differential interference non-contact vibration measurement |
CN104126103A (en) * | 2011-12-05 | 2014-10-29 | 直观外科手术操作公司 | Method and apparatus for motion compensation in interferometric sensing systems |
CN104766281A (en) * | 2015-03-30 | 2015-07-08 | 中国资源卫星应用中心 | Correction method for uneven direct current component in interference image |
-
2017
- 2017-01-11 CN CN201710020296.7A patent/CN106839972B/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4889986A (en) * | 1988-08-18 | 1989-12-26 | The United States Of America As Represented By The Secretary Of The Navy | Serial interferometric fiber-optic sensor array |
JP3630581B2 (en) * | 1999-03-04 | 2005-03-16 | アルプス電気株式会社 | Spread modulation signal receiver |
CN102045120A (en) * | 2010-10-29 | 2011-05-04 | 成都九洲电子信息系统有限责任公司 | Vibration signal identification method for optical fiber perimeter system |
CN104126103A (en) * | 2011-12-05 | 2014-10-29 | 直观外科手术操作公司 | Method and apparatus for motion compensation in interferometric sensing systems |
CN102494758A (en) * | 2011-12-29 | 2012-06-13 | 中国船舶重工集团公司第七一○研究所 | Apparatus for eliminating interference of vibration noise on velocity hydrophone |
CN102636196A (en) * | 2012-04-09 | 2012-08-15 | 天津大学 | Distributed disturbance sensor on basis of Rayleigh scattering spectrum related coefficient and demodulating method thereof |
CN102832620A (en) * | 2012-08-31 | 2012-12-19 | 天津理工大学 | System and method for harmonic detection and control of APF (active power filter) based on windowed all-phase FFT (fast Fourier transform) |
CN103424183A (en) * | 2013-08-06 | 2013-12-04 | 常州大学 | Method for eliminating abnormal interference on detection for mechanical vibration signals |
CN103759804A (en) * | 2014-01-23 | 2014-04-30 | 安徽大学 | Method and device for optical fiber white light differential interference non-contact vibration measurement |
CN104766281A (en) * | 2015-03-30 | 2015-07-08 | 中国资源卫星应用中心 | Correction method for uneven direct current component in interference image |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113739708A (en) * | 2021-08-20 | 2021-12-03 | 哈尔滨工程大学 | Double-path white light interference differential measurement device and method |
Also Published As
Publication number | Publication date |
---|---|
CN106839972B (en) | 2019-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102798411B (en) | System and method for distributed optical fibre sensing measurement based on Brillouin scattering | |
CN102332956B (en) | Dispersion compensation method for broadband light source | |
CN104279959B (en) | A kind of new method of the fine length of use vector network analyzer precise measuring | |
CN103017934B (en) | Self-correction method capable of eliminating wavelength dispersion of distributed Raman temperature measurement system | |
CN102288388B (en) | Device and method for improving polarization-maintaining optical fiber polarization coupling measurement precision and symmetry | |
CN103364775A (en) | Optical frequency comb calibration-based dual-color laser scanning absolute distance measuring device and method | |
CN105571830B (en) | The method for measuring super-narrow line width laser device laser line width | |
CN102914421B (en) | A kind of method for measuring polarization interference in optical birefringence medium and equipment thereof | |
CN107894327B (en) | Measuring device for optical fiber dispersion coefficient | |
CN104006948B (en) | Based on the method that multimodal division cycle demodulates polarization maintaining optical fibre polarization coupled point position | |
CN113804404B (en) | Light source sweep frequency nonlinear correction method for optical frequency domain polarization crosstalk measurement | |
CN101592551A (en) | A kind of beat length of polarization maintaining optical fiber method of testing and proving installation based on the Sagnac interferometer | |
CN104776871B (en) | Optical fiber Brillouin distributed measurement light path, apparatus and method | |
CN103968864B (en) | For the maximal phase seemingly Match Analysis of the frequency displacement of Measurement accuracy Brillouin spectrum | |
CN101660998B (en) | Method for measuring group delay by using wavelet transformation | |
CN106989904A (en) | A kind of measuring method of polarization maintaining optical fibre extinction ratio | |
CN105181152A (en) | Calculation method for frequency shift of distributed Brillouin scattered spectrum | |
CN102419247B (en) | High-precision detection device and method of reflection type optical fiber turbulence | |
CN102706538B (en) | Polarization maintaining optical fiber extinction-ratio self-adaptive demodulation method | |
CN103411689B (en) | Based on optical maser wavelength direct measuring method and the device of the orthogonal line polarized light of single-frequency | |
CN113218320B (en) | OFDR (offset-field-of-view) large strain measurement method based on distance domain compensation | |
CN104006950B (en) | A kind of polarization maintaining optical fibre birefringence dispersion measuring method | |
CN105953919A (en) | All-fiber Fourier expansion spectrum analyzer | |
CN103822646A (en) | Method for demodulating long distance light delayed effect in fiber bragg grating demodulation system | |
CN106839972B (en) | A kind of interference signal processing method of full light fiber white light interference instrument |
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 |