CN112946611B - Sweep frequency nonlinear correction distance measurement method based on similar triangular interpolation sampling - Google Patents

Sweep frequency nonlinear correction distance measurement method based on similar triangular interpolation sampling Download PDF

Info

Publication number
CN112946611B
CN112946611B CN202110154331.0A CN202110154331A CN112946611B CN 112946611 B CN112946611 B CN 112946611B CN 202110154331 A CN202110154331 A CN 202110154331A CN 112946611 B CN112946611 B CN 112946611B
Authority
CN
China
Prior art keywords
interferometer
measuring
signal
sweep
frequency
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.)
Active
Application number
CN202110154331.0A
Other languages
Chinese (zh)
Other versions
CN112946611A (en
Inventor
甘雨
刘国栋
段长昊
陈凤东
刘炳国
卢丙辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN202110154331.0A priority Critical patent/CN112946611B/en
Publication of CN112946611A publication Critical patent/CN112946611A/en
Application granted granted Critical
Publication of CN112946611B publication Critical patent/CN112946611B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/497Means for monitoring or calibrating
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/02Systems using the reflection of electromagnetic waves other than radio waves
    • G01S17/06Systems determining position data of a target
    • G01S17/08Systems determining position data of a target for measuring distance only
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/483Details of pulse systems
    • G01S7/486Receivers
    • G01S7/4861Circuits for detection, sampling, integration or read-out

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Electromagnetism (AREA)
  • Instruments For Measurement Of Length By Optical Means (AREA)

Abstract

A frequency sweep nonlinear correction ranging method based on similar triangular interpolation sampling relates to the technical fields of Frequency Sweep Interferometry (FSI) measurement, FMCW laser radar and the like, aims at solving the problems that an accurate resampling sequence cannot be provided by the existing method, and the frequency sweep nonlinearity is not thoroughly eliminated to cause the reduction of the frequency sweep interferometry precision, is used for eliminating the frequency spectrum broadening effect generated by the frequency sweep nonlinearity of a frequency sweep laser, and can effectively improve the excellent property of measuring the full width at half maximum of a frequency spectrum. Especially for measuring extreme distances, the error is smaller than taking a near zero point as a sampling point. Meanwhile, similar interpolation algorithms are faster in signal processing speed due to the use of matrix multiplication. The invention ensures the high-precision and high-real-time processing of the measurement result when the absolute distance measurement system carries out remote measurement.

Description

Sweep frequency nonlinear correction distance measurement method based on similar triangular interpolation sampling
Technical Field
The invention relates to the technical fields of Frequency Sweep Interferometry (FSI) measurement, FMCW laser radar and the like, in particular to a frequency sweep nonlinear correction ranging system and method based on similar triangular interpolation sampling.
Background
The sweep frequency interference measurement has the advantages of low emission power, no range finding ambiguity, no guide rail matching, applicability to cooperative target measurement, capability of realizing high-resolution measurement and the like, so that the sweep frequency interference measurement is widely applied to the field of high-precision absolute distance measurement, such as the fields of frequency modulation continuous laser radar, optical frequency meters, optical coherence tomography and the like. The basic principle is that a fixed frequency difference is formed between an emitted measuring optical signal and a swept frequency measuring optical signal reflected by a measured target to determine the absolute distance of the target, and the absolute distance parameter of the target is reflected by means of signal frequency domain characteristics which are easy to process. High distance resolution and measurement accuracy can be obtained by high-linearity broadband frequency sweep measurement light, but a linear tuned laser in the current market cannot meet the condition of high-linearity linear frequency sweep, nonlinear frequency sweep brings difficulty to signal frequency domain characteristic analysis, and a non-single frequency beat signal not only can reduce the distance resolution, but also greatly reduces the measurement accuracy. In the existing method, an auxiliary interferometer is adopted to sample sweep frequency optical signals, and a sampling sequence is formed by extracting the zero crossing point of the auxiliary interferometer signals to collect measurement interferometer signals. However, the zero crossing point of the auxiliary interferometer signal is not necessarily positioned in the auxiliary interferometer signal sampling, so that the resampling sequence cannot be accurately determined, the sweep nonlinearity cannot be completely eliminated, and the precision of sweep interference measurement is reduced.
Disclosure of Invention
The purpose of the invention is: aiming at the problems that the existing method can not provide an accurate resampling sequence and sweep frequency nonlinearity elimination is not thorough, so that sweep frequency interference measurement precision is reduced, a sweep frequency nonlinearity correction distance measurement method based on similar triangular interpolation sampling is provided.
The technical scheme adopted by the invention to solve the technical problems is as follows:
a sweep frequency nonlinear correction ranging method based on similar triangular interpolation sampling is realized by utilizing a sweep frequency nonlinear correction device;
the frequency sweep nonlinear correction device comprises: the system comprises a measuring interferometer, an auxiliary interferometer and an acquisition card;
the measuring interferometer comprises a PBS, a focusing optical system, a 1/4 wave plate and a balanced detector;
the auxiliary interferometer comprises two 3dB optical fiber couplers;
the sweep laser beam splitting ratio is 95: after the optical fiber coupler of 5 splits light, 95% sweep light and 5% sweep light are obtained, and the splitting ratio of the 95% sweep light is 99: the optical fiber coupler of 1 is divided into measuring light of a measuring interferometer and reference light of the measuring interferometer, the measuring light of the measuring interferometer sequentially passes through a PBS (polarization beam splitter), a focusing optical system and a 1/4 wave plate and then is emitted, and after being reflected by a measuring target, an original path returns and is mixed with the reference light of the measuring interferometer on a detection plane of a balance detector to obtain a beat frequency signal of the measuring interferometer, namely an optical signal of the measuring interferometer;
the 5% sweep light is divided into 50% by a 3dB fiber coupler: measuring light of an auxiliary interferometer and reference light of the auxiliary interferometer of 50, combining the measuring light of the auxiliary interferometer and the reference light of the auxiliary interferometer by a 3dB optical fiber coupler, and mixing on a balance detector to obtain a beat frequency signal of the auxiliary interferometer, namely an optical signal of the auxiliary interferometer;
the acquisition card is used for acquiring and storing a beat frequency signal of the measuring interferometer and a beat frequency signal of the auxiliary interferometer;
the sweep frequency nonlinear correction ranging method based on similar triangular interpolation sampling comprises the following specific steps:
the method comprises the following steps: obtaining a measurement interferometer optical signal ImcAnd auxiliary interferometer optical signal IdcAnd obtaining a measurement interferometer signal sequence I according to the measurement interferometer optical signal and the auxiliary interferometer optical signalmAnd an auxiliary interferometer signal sequence Id
Step two: extracting a signal sequence I of the measuring interferometermAnd an auxiliary interferometer signal sequence IdIn odd order1m、I1dAnd even-order data sequences I2m、I2d
Step three: will I1dAnd I2dCorresponding multiplication is carried out, the multiplication results are screened to obtain a plurality of groups of adjacent data points with zero points between two points, then data point labels are recorded, and a set I is formedo
Step four: will be set IoPerforming similar fitting on each group of adjacent data points to obtain a zero horizontal coordinate value set A between each group of adjacent data pointsoAfter that, the set A isoAfter arrangement, determining a resampling sequence Bo
Step five: using set IoData point index pair I of two adjacent points in which zero exists1mAnd I2mReading data to obtain a resampling sequence BoMeasuring interferometer data points I corresponding to two sides1mo、I2moThen use I1mo、I2moAnd a resampling sequence BoTo I1mo、I2moPerforming similar interpolation to obtain a corrected signal of the measuring interferometer;
step six: and carrying out spectrum analysis on the corrected measuring interferometer signal, searching for the maximum value of the spectrum peak value, and carrying out corresponding linear transformation to obtain a corrected distance measuring result.
Further, the swept laser is obtained by an external cavity laser.
Further, the frequency of the beat frequency signal of the measuring interferometer is determined according to the arm length difference of the measuring interferometer and the measuring distance.
Further, the measurement interferometer optical signal is represented as:
Figure BDA0002933990470000031
wherein, AmFor measuring interferometer optical signal amplitude, 2RmTo measure the interferometer arm length difference, Δ f (t) is the sweep frequency variation, and c is the speed of light.
Further, the auxiliary interferometer light signal is represented as:
Figure BDA0002933990470000032
wherein A is0To assist the interferometer in amplitude of the optical signal, R0To assist interferometer arm length difference, Δ f (t) is sweep frequency variation, and c is speed of light.
Further, the sweep variation Δ f (t) is expressed as:
Figure BDA0002933990470000033
wherein k is the resample sequence sample sequence index.
Further, the interferometer signal measured in the fifth step is expressed as:
Figure BDA0002933990470000034
the beneficial effects of the invention are:
the method is used for eliminating the spectrum broadening effect generated by the non-linearity of the sweep frequency laser, and can effectively improve the excellent property of measuring the full width at half maximum of the spectrum. Especially for measuring extreme distances, the error is smaller than taking a near zero point as a sampling point. Meanwhile, similar interpolation algorithms are faster in signal processing speed due to the use of matrix multiplication. The invention ensures the high-precision and high-real-time processing of the measuring result when the absolute distance measuring system carries out remote measurement.
Drawings
FIG. 1 is a diagram of the optical path of the hardware system of the present application;
FIG. 2 is a schematic diagram of a linear fit similarity difference algorithm;
FIG. 3 is a flow chart of swept nonlinear correction ranging with similar triangular interpolation sampling;
FIG. 4 is a schematic diagram of a 0.1s sweep time simulation;
FIG. 5 is a schematic diagram of a 0.105s sweep time simulation;
FIG. 6 is a schematic diagram of a sweep time 0.110s simulation;
FIG. 7 is a schematic diagram of a 0.115s sweep time simulation;
FIG. 8 is a schematic diagram of a sweep time 0.12s simulation;
FIG. 9 is a schematic diagram of a sweep time 0.125s simulation;
FIG. 10 is a schematic diagram of a sweep time 0.13s simulation;
FIG. 11 is a schematic diagram of a sweep time of 0.135s simulation;
FIG. 12 is a schematic diagram of a 0.14s sweep time simulation;
FIG. 13 is a schematic diagram of a sweep time 0.145s simulation;
FIG. 14 is a schematic diagram of a 0.15s sweep time simulation;
FIG. 15 is a schematic diagram of a sweep time 0.155s simulation;
FIG. 16 is a schematic diagram of a 0.16s sweep time simulation;
FIG. 17 is a schematic diagram of a sweep time 0.165s simulation;
FIG. 18 is a schematic diagram of a 0.17s sweep time simulation;
FIG. 19 is a schematic diagram of a sweep time 0.175s simulation;
FIG. 20 is a schematic diagram of a sweep time 0.18s simulation;
FIG. 21 is a schematic diagram of a sweep time 0.185s simulation;
FIG. 22 is a schematic diagram of a sweep time 0.19s simulation;
fig. 23 is a simulation diagram of sweep time 0.195s.
Detailed Description
It should be noted that, in the case of conflict, the various embodiments disclosed in the present application may be combined with each other.
The first specific implementation way is as follows: specifically describing the present embodiment with reference to fig. 1, fig. 2 and fig. 3, the method for distance measurement based on sweep frequency nonlinear correction of similar triangular interpolation sampling according to the present embodiment is implemented by using a sweep frequency nonlinear correction device;
the frequency sweep nonlinear correction device comprises: the system comprises a measuring interferometer, an auxiliary interferometer and an acquisition card;
the measuring interferometer comprises a PBS, a focusing optical system, a 1/4 wave plate and a balanced detector;
the auxiliary interferometer comprises two 3dB optical fiber couplers;
the sweep laser beam splitting ratio is 95: after the optical fiber coupler of 5 splits light, 95% sweep light and 5% sweep light are obtained, and the splitting ratio of the 95% sweep light is 99: the optical fiber coupler of 1 is divided into measuring light of a measuring interferometer and reference light of the measuring interferometer, the measuring light of the measuring interferometer sequentially passes through a PBS (polarization beam splitter), a focusing optical system and a 1/4 wave plate and then is emitted, and after being reflected by a measuring target, an original path returns and is mixed with the reference light of the measuring interferometer on a detection plane of a balance detector to obtain a beat frequency signal of the measuring interferometer, namely an optical signal of the measuring interferometer;
the 5% sweep light is divided into 50% by a 3dB fiber coupler: measuring light of an auxiliary interferometer and reference light of the auxiliary interferometer of 50, combining the measuring light of the auxiliary interferometer and the reference light of the auxiliary interferometer by a 3dB optical fiber coupler, and mixing on a balance detector to obtain a beat frequency signal of the auxiliary interferometer, namely an optical signal of the auxiliary interferometer;
the acquisition card is used for acquiring and storing the beat frequency signal of the measurement interferometer and the beat frequency signal of the auxiliary interferometer;
the sweep frequency nonlinear correction ranging method based on similar triangular interpolation sampling comprises the following specific steps:
the method comprises the following steps: obtaining measurement interferometer optical signals ImcAnd auxiliary interferometer optical signal IdcAnd obtaining a measurement interferometer signal sequence I according to the measurement interferometer optical signal and the auxiliary interferometer optical signalmAnd an auxiliary interferometer signal sequence Id
Step two: extracting a signal sequence I of the measuring interferometermAnd auxiliary interferometer signal sequence IdIn odd order1m、I1dAnd even-order data sequences I2m、I2d
Step three: will I1dAnd I2dCorresponding multiplication is carried out, the multiplication results are screened to obtain a plurality of groups of adjacent data points with zero point between two points, then data point labels are recorded, and a set I is formedo
Step four: set IoPerforming similar fitting on each group of adjacent data points to obtain a zero horizontal coordinate value set A between each group of adjacent data pointsoAfter that, the set A isoDetermining a resampling sequence B after the arrangemento
Step five: using set IoData point index pair I of two adjacent points in which zero exists1mAnd I2mReading data to obtain a resample sequence BoMeasuring interferometer data points I corresponding to two sides1mo、I2moThen use I1mo、I2moAnd a resampling sequence BoTo I1mo、I2moPerforming similar interpolation to obtain a corrected signal of the measuring interferometer;
step six: and carrying out spectrum analysis on the corrected measuring interferometer signal, searching for the maximum value of the spectrum peak value, and carrying out corresponding linear transformation to obtain a corrected distance measuring result.
According to the method and the device, the sampling points of the auxiliary interferometer are found, the sampling points are fine enough, linear fitting can be conducted on the sampling points on two sides of the zero crossing point, and the time domain coordinate of the zero crossing point of the auxiliary interferometer is found accurately through corresponding operation of similar triangles. Compared with the traditional method, the method provided by the patent improves the measurement precision of the signal of the measuring interferometer. In simulation, the fact that when the optical path difference of the auxiliary interferometer is 2 times of the optical path difference of the measuring interferometer, the spectrum peak of the measuring signal changes after the zero-crossing point resampling is not accurately positioned at the moment, and finally micrometer-level deviation is generated in absolute distance measurement.
The method is used for eliminating the spectrum broadening effect generated by the non-linearity of the sweep frequency laser, and can effectively improve the excellent property of measuring the full width at half maximum of the spectrum. Especially for measuring extreme distances, the error is smaller than taking a near zero point as a sampling point. Meanwhile, similar interpolation algorithms are faster in signal processing speed due to the use of matrix multiplication. The high-precision and high-real-time processing of the measurement result when the absolute distance measurement system carries out remote measurement is ensured.
The hardware part of the invention is mainly divided into a sweep frequency interference distance measuring light path, an auxiliary interferometer 1 and an acquisition card. The frequency sweep interference distance measuring optical path mainly comprises a Mach-Zehnder interferometer formed by a plurality of fiber-coupled PBS, WDM, a focusing optical system, a 1/4 wave plate and a balance detector, the frequency sweep nonlinearity is measured by an auxiliary interferometer, and the signal of the measuring interferometer is resampled by adopting a frequency sampling method of software to eliminate the nonlinear phase introduced by frequency sweep laser. Since the sampling points hardly coincide with the theoretical zero crossings, the sequence of frequency samples generated on the basis of the auxiliary interferometer cannot be determined accurately when the auxiliary interferometer signal frequency drops, and at the same time re-sampling on the basis of this sequence will not completely eliminate the non-linearity of the measuring interferometer. Therefore, the patent proposes to judge the real zero crossing point of the signal based on the triangular similarity, and further improve the nonlinear elimination effect of the frequency sampling method.
Two paths of optical signals are generated in the optical path, wherein one path of optical signal passes through a measurement target, and a return optical beat frequency signal of the optical signal is called as a measurement interferometer optical signal; and the other optical signal passes through the optical path of the auxiliary interferometer, the optical beat frequency signal received by the detector is called as the optical signal of the auxiliary interferometer, and the two optical signals are expressed mathematically as follows:
Figure BDA0002933990470000061
Figure BDA0002933990470000062
in the formula Im(t)、I0(t) measuring the interferometer light signal on the auxiliary interferometer light signal, Am、A0For measuring the amplitude of the optical signal of the interferometer at the amplitude of the optical signal of the auxiliary interferometer, 2Rm、R0In order to measure the interferometer arm length difference and the auxiliary interferometer arm length difference, Δ f (t) is the sweep frequency variation.
Figure BDA0002933990470000063
Figure BDA0002933990470000064
To eliminate the sweep nonlinearity caused by Δ f (t), the signal processing procedure is called frequency sampling method by determining the zero crossings of equation 2 and then resampling equation 1 by the sampling sequence generated by these zero crossings. The sampling mode of the frequency sampling method is shown in formula 3, in the sampling process, the sweep frequency nonlinearity delta f (t) in the formula can be represented by a linear sequence, and a measurement interferometer signal after the frequency sampling method is shown in formula 4, so that a single-frequency signal related to the measurement distance is obtained theoretically. Fig. 3 shows a flow chart of the method adopted by the patent.
The determination of zero crossings by linear fitting will prove to be more effective than the determination of sampling points near zero in the following simulations, which can be illustrated by the following simulations.
The simulation conditions are set such that the measurement interferometer arm length difference setting range is 1m to 54m, the auxiliary interferometer arm length difference is 109m, the sampling frequency of the signal in the simulation is set to 10MHz, and the sweep time is set to 0.1s, 0.105s, 0.11s, 0.115s, 0.120s, 0.125s, 0.130s, 0.135s, 0.140s, 0.145s, 0.150s, 0.155s, 0.160s, 0.165s, 0.170s, 0.175s, 0.180s, 0.185s, 0.190s, 0.195s. The simulation shows that the measurement error generated when the measurement distance approaches the arm length difference of the auxiliary interferometer can be eliminated to a great extent, and the algorithm value deviation caused by data processing is eliminated. In the simulation diagram, linear refers to simulation data points by adopting a Linear fitting method of similar triangles, and Near refers to simulation data points by adopting a sampling point method Near a zero point.
Simulation description:
according to the description of fig. 4 to fig. 23, simulation is performed, the distance from the measurement target to the balanced detector in the simulation is a variable from 1m to 54m, and it can be seen from each simulation that when the measurement interferometer arm length difference gradually approaches to half of the auxiliary interferometer arm length difference, the sampling method using the proximity point as the resampling sequence cannot completely eliminate the influence of the sweep frequency nonlinearity on the measurement interferometer frequency spectrum, and the nonlinear correction method using the similarity fitting to determine the resampling sequence has smaller measurement error. The method adopting the near point as the resampling sequence is sensitive to the parameter of the sequence length, and the distance demodulation error and the sequence length are in a nonlinear relation in simulation. As demonstrated in the above simulations, the sampling method employed in this patent can ensure the effect on the cancellation of non-linear errors in fast distance measurements, attenuating the errors introduced by clock jitter. In the actual measurement process, the method adopted by the patent effectively avoids errors introduced by an algorithm.
It should be noted that the detailed description is only for explaining and explaining the technical solution of the present invention, and the scope of protection of the claims is not limited thereby. It is intended that all such modifications and variations be included within the scope of the invention as defined in the following claims and the description.

Claims (7)

1. The sweep frequency nonlinear correction distance measurement method based on similar triangular interpolation sampling is characterized in that the method is realized by utilizing a sweep frequency nonlinear correction device;
the frequency sweep nonlinear correction device comprises: the system comprises a measuring interferometer, an auxiliary interferometer and an acquisition card;
the measuring interferometer comprises a PBS, a focusing optical system, a 1/4 wave plate and a balanced detector;
the auxiliary interferometer comprises two 3dB optical fiber couplers;
the sweep laser beam splitting ratio is 95: after the optical fiber coupler of 5 splits light, 95% sweep light and 5% sweep light are obtained, and the splitting ratio of the 95% sweep light is 99: the optical fiber coupler of 1 is divided into measuring light of a measuring interferometer and reference light of the measuring interferometer, the measuring light of the measuring interferometer sequentially passes through a PBS (polarization beam splitter), a focusing optical system and a 1/4 wave plate and then is emitted, and after being reflected by a measuring target, an original path returns and is mixed with the reference light of the measuring interferometer on a detection plane of a balance detector to obtain a beat frequency signal of the measuring interferometer, namely an optical signal of the measuring interferometer;
the 5% sweep light is divided into 50% by a 3dB fiber coupler: 50, combining the measuring light of the auxiliary interferometer and the reference light of the auxiliary interferometer by a 3dB optical fiber coupler, and mixing on a balance detector to obtain an auxiliary interferometer beat frequency signal, namely an auxiliary interferometer optical signal;
the acquisition card is used for acquiring and storing the beat frequency signal of the measurement interferometer and the beat frequency signal of the auxiliary interferometer;
the sweep frequency nonlinear correction ranging method based on similar triangular interpolation sampling comprises the following specific steps:
the method comprises the following steps: obtaining a measurement interferometer optical signal ImcAnd auxiliary interferometer optical signal IdcAnd obtaining a measurement interferometer signal sequence I according to the measurement interferometer optical signal and the auxiliary interferometer optical signalmAnd an auxiliary interferometer signal sequence Id
Step two: extracting a signal sequence I of the measuring interferometermAnd an auxiliary interferometer signal sequence IdIn odd order1m、I1dAnd even-order data sequences I2m、I2d
Step three: will I1dAnd I2dCorresponding multiplication is carried out, the multiplication results are screened to obtain a plurality of groups of adjacent data points with zero point between two points, then data point labels are recorded, and a set I is formedo
Step four: will be set IoPerforming similar fitting on each group of adjacent data points to obtain a zero horizontal coordinate value set A between each group of adjacent data pointsoAfter that, the set A isoDetermining a resampling sequence B after the arrangemento
Step five: using set IoData point index pair I of two adjacent points in which zero exists1mAnd I2mReading data to obtain a resample sequence BoMeasuring interferometer data points I corresponding to two sides1mo、I2moThen use I1mo、I2moAnd a resampling sequence BoTo I1mo、I2moPerforming similar interpolation to obtain a corrected signal of the measuring interferometer;
step six: and carrying out spectrum analysis on the corrected measuring interferometer signal, searching for the maximum value of the spectrum peak value, and carrying out corresponding linear transformation to obtain a corrected distance measuring result.
2. A swept frequency nonlinear correction ranging method based on similar triangular interpolation sampling according to claim 1, wherein the swept frequency laser is obtained by an external cavity laser.
3. A swept frequency nonlinear correction ranging method based on similar triangular interpolation sampling according to claim 2, wherein the frequency of the beat signal of the measuring interferometer is determined according to the arm length difference and the measuring distance of the measuring interferometer.
4. A swept frequency nonlinear correction ranging method based on similar triangular interpolation sampling according to claim 1, characterized in that the measurement interferometer optical signal is represented as:
Figure FDA0002933990460000021
wherein A ismFor measuring interferometer optical signal amplitude, 2RmTo measure the interferometer arm length difference, Δ f (t) is the sweep frequency variation, and c is the speed of light.
5. A swept frequency nonlinear correction ranging method based on similar triangular interpolation sampling according to claim 1, characterized in that the auxiliary interferometer optical signal is represented as:
Figure FDA0002933990460000022
wherein A is0To assist the interferometer in amplitude of the optical signal, R0To assist the interferometer arm length difference, Δ f (t) is the sweep frequency variation, and c is the speed of light.
6. A swept frequency nonlinear correction ranging method based on similar triangular interpolation sampling according to claim 4 or 5, characterized in that the sweep frequency variation Δ f (t) is expressed as:
Figure FDA0002933990460000023
wherein k is the resample sequence sample sequence index.
7. A swept frequency nonlinear correction ranging method based on similar triangular interpolation sampling according to claim 6, wherein the measured interferometer signal in the fifth step is represented as:
Figure FDA0002933990460000031
CN202110154331.0A 2021-02-04 2021-02-04 Sweep frequency nonlinear correction distance measurement method based on similar triangular interpolation sampling Active CN112946611B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110154331.0A CN112946611B (en) 2021-02-04 2021-02-04 Sweep frequency nonlinear correction distance measurement method based on similar triangular interpolation sampling

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110154331.0A CN112946611B (en) 2021-02-04 2021-02-04 Sweep frequency nonlinear correction distance measurement method based on similar triangular interpolation sampling

Publications (2)

Publication Number Publication Date
CN112946611A CN112946611A (en) 2021-06-11
CN112946611B true CN112946611B (en) 2022-11-01

Family

ID=76243832

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110154331.0A Active CN112946611B (en) 2021-02-04 2021-02-04 Sweep frequency nonlinear correction distance measurement method based on similar triangular interpolation sampling

Country Status (1)

Country Link
CN (1) CN112946611B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113776445A (en) * 2021-07-20 2021-12-10 重庆大学 Single-frequency-interference rotor and stator axial clearance high-speed dynamic measurement method
CN115327514B (en) * 2022-08-10 2023-05-05 哈尔滨工业大学 Sweep frequency interference dynamic measurement system and measurement method based on phase transfer
CN115327565B (en) * 2022-08-10 2023-05-05 哈尔滨工业大学 Phase transfer-based nonlinear correction system and correction method for sweep frequency interferometry beat signals
CN115327515B (en) * 2022-08-10 2023-04-18 哈尔滨工业大学 Double-sweep frequency interference dynamic measurement system and measurement method based on phase transmission
CN116659394B (en) * 2023-06-07 2024-03-01 哈尔滨工业大学 Time-delay resampling-based nonlinear synchronous error correction method for frequency sweep interferometry
CN116659395B (en) * 2023-06-07 2024-02-23 哈尔滨工业大学 Phase compensation-based nonlinear synchronous error correction method for frequency sweep interferometry

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105259548A (en) * 2015-10-29 2016-01-20 哈尔滨工业大学 Dispersion mismatch correction method used in FMCW absolute distance measurement technology
CN105954735A (en) * 2016-07-15 2016-09-21 哈尔滨工业大学 Improved high-speed dispersion mismatch correction method for FMCW absolute distance measuring technique
CN106226775A (en) * 2016-07-15 2016-12-14 哈尔滨工业大学 A kind of absolute distance dynamic measurement system based on swept frequency interferometer and measuring method thereof
CN109029271A (en) * 2018-09-05 2018-12-18 天津大学 DFB array swept light source optical fiber frequency domain interfeerometry ranging system and method
CN109188453A (en) * 2018-09-11 2019-01-11 哈尔滨工业大学 Dynamic frequency scanning interfeerometry ranging system and distance measuring method based on locking phase gamma correction
CN109188454A (en) * 2018-09-11 2019-01-11 哈尔滨工业大学 Dynamic frequency scanning interfeerometry ranging system and method based on digital servo-control gamma correction
CN110716208A (en) * 2019-10-18 2020-01-21 中国科学院光电技术研究所 Frequency sweep interference dynamic absolute distance measuring system based on acousto-optic modulator
CN111397644A (en) * 2020-03-26 2020-07-10 南京大学 Laser nonlinear tuning effect compensation system and compensation method for optical frequency domain reflectometer
CN111694008A (en) * 2020-06-08 2020-09-22 中国科学院光电技术研究所 Method for eliminating laser mode hopping influence in frequency sweep coherent ranging
CN111948664A (en) * 2020-08-18 2020-11-17 中国科学院光电技术研究所 Dispersion compensation method of frequency modulation continuous wave laser radar based on dispersion coefficient modulation
CN112114326A (en) * 2020-09-21 2020-12-22 哈尔滨工业大学 Frequency sweep signal splicing method and device for FMCW distance measurement

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102322880B (en) * 2011-08-18 2013-06-05 天津大学 Polarization sensitive distributive optical frequency domain reflection disturbance sensor and demodulation method
CN106248121B (en) * 2016-08-11 2018-03-06 天津大学 The fiber grating sensing demodulation device and demodulation method of suppression are fluctuated under environment alternating temperature

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105259548A (en) * 2015-10-29 2016-01-20 哈尔滨工业大学 Dispersion mismatch correction method used in FMCW absolute distance measurement technology
CN105954735A (en) * 2016-07-15 2016-09-21 哈尔滨工业大学 Improved high-speed dispersion mismatch correction method for FMCW absolute distance measuring technique
CN106226775A (en) * 2016-07-15 2016-12-14 哈尔滨工业大学 A kind of absolute distance dynamic measurement system based on swept frequency interferometer and measuring method thereof
CN109029271A (en) * 2018-09-05 2018-12-18 天津大学 DFB array swept light source optical fiber frequency domain interfeerometry ranging system and method
CN109188453A (en) * 2018-09-11 2019-01-11 哈尔滨工业大学 Dynamic frequency scanning interfeerometry ranging system and distance measuring method based on locking phase gamma correction
CN109188454A (en) * 2018-09-11 2019-01-11 哈尔滨工业大学 Dynamic frequency scanning interfeerometry ranging system and method based on digital servo-control gamma correction
CN110716208A (en) * 2019-10-18 2020-01-21 中国科学院光电技术研究所 Frequency sweep interference dynamic absolute distance measuring system based on acousto-optic modulator
CN111397644A (en) * 2020-03-26 2020-07-10 南京大学 Laser nonlinear tuning effect compensation system and compensation method for optical frequency domain reflectometer
CN111694008A (en) * 2020-06-08 2020-09-22 中国科学院光电技术研究所 Method for eliminating laser mode hopping influence in frequency sweep coherent ranging
CN111948664A (en) * 2020-08-18 2020-11-17 中国科学院光电技术研究所 Dispersion compensation method of frequency modulation continuous wave laser radar based on dispersion coefficient modulation
CN112114326A (en) * 2020-09-21 2020-12-22 哈尔滨工业大学 Frequency sweep signal splicing method and device for FMCW distance measurement

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
OFDR光源非线性调谐引起的距离偏移校正研究;唐超 等;《光传感》;20201231(第4期);30-32 *
基于相位补偿的调频连续波大长度测距中的色散校正;史春钊 等;《红外与毫米波学报》;20181031;第37卷(第5期);642-648 *

Also Published As

Publication number Publication date
CN112946611A (en) 2021-06-11

Similar Documents

Publication Publication Date Title
CN112946611B (en) Sweep frequency nonlinear correction distance measurement method based on similar triangular interpolation sampling
CN110068828B (en) Device and dispersion compensation method based on laser frequency modulation continuous wave long-distance measurement
CN102332956B (en) Dispersion compensation method for broadband light source
CN112923960B (en) Optical fiber parameter measuring device for correcting nonlinear tuning effect
CN110487313A (en) Light source frequency sweep Nonlinear Self-tuning method in optical frequency domain reflection technology
CN112051583B (en) Beat frequency signal nonlinear correction method in FMCW distance measurement system
CN108663684A (en) A kind of phase difference ranging method based on equal optical frequency intervals resampling
CN105259548A (en) Dispersion mismatch correction method used in FMCW absolute distance measurement technology
WO2021082377A1 (en) Phase derivation-based optical component delay measurement method, and device
CN102420650A (en) Device and method for inhibiting nonlinear scanning of laser of optical frequency domain reflectometer
CN111694008A (en) Method for eliminating laser mode hopping influence in frequency sweep coherent ranging
CN115327514B (en) Sweep frequency interference dynamic measurement system and measurement method based on phase transfer
CN113029034A (en) Nonlinear correction and range expansion device and method for laser frequency sweep interferometry
CN115327565B (en) Phase transfer-based nonlinear correction system and correction method for sweep frequency interferometry beat signals
CN111948664A (en) Dispersion compensation method of frequency modulation continuous wave laser radar based on dispersion coefficient modulation
CN101660998B (en) Method for measuring group delay by using wavelet transformation
CN109990713A (en) A kind of high-resolution phase detection method based on plane grating laser interferometer
CN111751834B (en) High-speed high-precision dynamic ranging method based on optical frequency modulation interference and single-frequency interference
CN112462380A (en) Dispersion compensation method based on laser frequency modulation continuous wave long-distance ranging
CN112798025B (en) Method for improving OFDR measurement spatial resolution and OFDR system
CN112034475A (en) FMCW laser radar frequency-sweep light source mode hopping compensation method
CN105716521A (en) Device and method for enlarging measurement range of frequency domain low coherent light interference
Wei et al. Frequency pair model for selection of signal spectral components to determine position of zero optical-path difference in a pulse-train interferometer
CN113253241B (en) Method for processing sweep frequency interference ranging signal
Liu et al. Optical-carried microwave interferometric ranging method based on frequency-scanning phase

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