WO2021155567A1 - Frequency modulation nonlinear correction-based ranging method and related device - Google Patents

Frequency modulation nonlinear correction-based ranging method and related device Download PDF

Info

Publication number
WO2021155567A1
WO2021155567A1 PCT/CN2020/074495 CN2020074495W WO2021155567A1 WO 2021155567 A1 WO2021155567 A1 WO 2021155567A1 CN 2020074495 W CN2020074495 W CN 2020074495W WO 2021155567 A1 WO2021155567 A1 WO 2021155567A1
Authority
WO
WIPO (PCT)
Prior art keywords
signal
frequency
beat
flight time
term
Prior art date
Application number
PCT/CN2020/074495
Other languages
French (fr)
Chinese (zh)
Inventor
巫红英
李强
Original Assignee
华为技术有限公司
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 华为技术有限公司 filed Critical 华为技术有限公司
Priority to PCT/CN2020/074495 priority Critical patent/WO2021155567A1/en
Priority to CN202080004674.2A priority patent/CN112639528B/en
Publication of WO2021155567A1 publication Critical patent/WO2021155567A1/en

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
    • 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
    • G01S17/32Systems determining position data of a target for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
    • G01S17/34Systems determining position data of a target for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of continuous, frequency-modulated waves while heterodyning the received signal, or a signal derived therefrom, with a locally-generated signal related to the contemporaneously transmitted signal

Definitions

  • This application relates to the field of lidar, and in particular to a ranging method and related devices based on frequency modulation nonlinear correction.
  • Frequency Modulated Continuous Wave Frequency Modulated Continuous, FMCW
  • Lidar Light Detection and Ranging, LiDAR
  • ranging ranging
  • speed measurement target detection
  • tracking and imaging recognition Light Detection and Ranging
  • autonomous driving atmospheric environment monitoring
  • geographic surveying and mapping unmanned aerial vehicles and other fields.
  • FMCW LiDAR obtains the position information of the measured target by measuring the coherence of the emitted light signal and the echo signal to obtain the frequency domain response of the beat signal.
  • the basic process is shown in Figure 1.
  • the frequency difference f b between the transmitted signal and the echo signal is estimated, and the flight time ⁇ is calculated according to the frequency difference f b , and finally the distance of the measured target is calculated according to the flight time.
  • FIG. 2 is a block diagram of a typical FMCW LiDAR system architecture.
  • the laser By injecting appropriate driving current into the tunable laser, the laser is driven to emit a laser signal. After the beam passes through the isolator, it is divided into two laser beams by the coupler. The measured target is called the measuring light. The measuring light is reflected by the target surface and collected by the lens and passed through the circulator. The reference light passes through the coupler.
  • Coherence occurs on the photosensitive surface of the balanced detector to generate a beat signal whose frequency is proportional to the flight time ⁇ ; through the beat signal Perform frequency domain analysis to obtain the beat signal frequency f b , that is, the frequency difference between the transmitted signal and the echo signal; then according to known parameters such as the frequency modulation period and the frequency modulation bandwidth, the flight time ⁇ can be calculated, and then the target distance can be estimated .
  • the frequency modulation curves of the transmitted signal and the echo signal should be consistent and the frequency changes linearly with time. There is only a flight delay between the two, so the frequency modulation linearity of the transmitted signal will greatly affect the FMCWLiDAR The accuracy of ranging.
  • the output frequency of the laser is controlled by changing the drive current of the laser. Due to the inherent frequency modulation effect of the semiconductor laser, the output frequency of the laser exhibits obvious nonlinearity to the modulation current, as shown in Figure 3.
  • the beat signal is no longer a single-frequency signal, its beat frequency will change with time, and the frequency spectrum will be broadened, resulting in a decrease in the signal-to-noise ratio of the beat signal, which will seriously affect To the ranging accuracy and ranging range.
  • the embodiments of the present application provide a ranging method and related devices based on frequency modulation nonlinear correction, so as to improve the ranging accuracy and ranging range of the lidar.
  • an embodiment of the present application provides a ranging method based on frequency modulation nonlinear correction, including:
  • the reference road beat signal and the measurement road beat signal are respectively obtained based on the laser emission signal; the frequency modulation nonlinearity is calculated according to the reference road beat signal; and according to the measurement road
  • the beat frequency signal is calculated to obtain the initial flight time ⁇ ; it is determined that the frequency modulation nonlinearity is not less than the first preset threshold, and the frequency modulation nonlinearity term ⁇ (t) of the laser emission signal is calculated according to the reference path beat frequency signal.
  • the frequency modulation nonlinear term ⁇ (t) of the laser emission signal performs nonlinear iterative correction on the measurement road beat signal to obtain the target flight time; the target distance is calculated according to the target flight time.
  • the non-linear iterative correction of the beat signal of the measurement circuit is performed based on the flight time ⁇ and the frequency modulation non-linear term ⁇ (t) of the laser emission signal, which eliminates the influence of the non-linear frequency modulation on the beat signal of the measurement circuit and makes the final estimation
  • the deviation between the target flight time and the actual flight time of the target is the smallest, and the target distance can be determined according to the target flight time, which improves the ranging accuracy of the lidar.
  • the non-linear iterative correction of the measurement path beat frequency signal according to the initial flight time ⁇ and the frequency modulation nonlinear term ⁇ (t) of the laser emission signal to obtain the target flight time includes:
  • the frequency modulation nonlinearity correction is performed on the measurement road beat signal, which eliminates the influence of the frequency modulation nonlinearity on the measurement road beat signal, and improves the flight time based on the corrected beat signal. Accuracy, thereby improving the accuracy of ranging.
  • determining the target flight time according to the corrected beat signal includes:
  • peak i-1 is the initial frequency spectrum peak
  • flight time ⁇ i-1 is the initial flight time ⁇ .
  • the predetermined condition peak i -peak i-1) / peak i-1 ⁇ Thr rate, Thr rate as the third predetermined threshold value.
  • the non-linear iterative correction of the measurement path beat frequency signal according to the initial flight time ⁇ and the frequency modulation nonlinear term ⁇ (t) of the laser emission signal to obtain the target flight time includes:
  • the frequency modulation nonlinear term ⁇ (t) of the laser emission signal is delayed according to the delay flight time ⁇ 'i-1 to obtain the frequency modulation corresponding to the echo signal of the measurement path
  • the non-linear term ⁇ (t- ⁇ 'i -1 ); the delayed flight time ⁇ 'i -1 is obtained based on the flight time ⁇ i -1 ;
  • the influence of FM nonlinearity on the measurement road beat signal can be eliminated to a large extent, and the accuracy of the flight time obtained based on the corrected beat signal is improved. , Thereby improving the accuracy of ranging.
  • obtaining the flight time ⁇ i according to the corrected beat signal includes:
  • the correction signal comprises a beat frequency signal a first correction or the second correction beat beat signal; an echo signal corresponding to the measured channel FM nonlinear term ⁇ (t- ⁇ 'i- 1) of Measure the beat signal of the road and perform FM nonlinear correction to obtain the corrected beat signal, including:
  • the frequency modulation nonlinearity term ⁇ (t) of the laser emission signal and the frequency modulation nonlinearity term ⁇ (t- ⁇ ) corresponding to the echo signal of the measurement path ' I-1 ) Perform FM nonlinear term correction on the beat signal of the measurement road to obtain the first corrected beat signal; when the signal-to-noise ratio of the beat signal of the measurement road is lower than the second preset threshold, according to the measurement road
  • the FM nonlinear term ⁇ (t- ⁇ ′i -1 ) corresponding to the wave signal performs receiving end FM nonlinear term correction on the compensated beat signal to obtain a second corrected beat signal;
  • the compensated beat signal is obtained based on the FM nonlinear term ⁇ (t) of the laser emitting signal to compensate the FM nonlinear term of the transmitter of the measurement path beat signal.
  • the frequency modulation nonlinear term ⁇ (t) of the laser emission signal the beat signal of the measurement path is compensated, and the compensated beat signal is obtained. Therefore, when correcting the beat signal of the measurement path, it is only necessary to calibrate the beat signal of the measurement path.
  • the FM nonlinear term of the echo signal performs FM nonlinear correction on the compensated beat signal, and also eliminates the influence of FM nonlinearity on the beat signal of the measurement path.
  • calculating the initial flight time ⁇ according to the measured road beat frequency signal includes:
  • the signal-to-noise ratio of the measured road beat signal is not lower than the second preset threshold, perform FFT on the measured road beat signal to obtain the frequency domain signal corresponding to the measured road beat signal;
  • the frequency domain signal determines the frequency of the measured road beat signal, and calculates the initial flight time ⁇ according to the frequency of the measured road beat signal;
  • the frequency of the measurement road beat signal is the frequency corresponding to the peak point position of the frequency domain signal corresponding to the measurement road beat signal, or the frequency corresponding to the midpoint position of the half-height width of the frequency domain signal.
  • calculating the initial flight time ⁇ according to the measured road beat frequency signal includes:
  • the measurement road beat signal is compensated by the FM nonlinear term at the transmitting end.
  • the frequency modulation nonlinear term of the laser emission signal is calculated based on the beat frequency signal of the reference path;
  • the frequency of the compensated beat signal is the frequency corresponding to the peak point position of the frequency domain signal corresponding to the compensated beat signal, or the frequency corresponding to the midpoint position of the half-height width of the frequency domain signal.
  • the initial flight time calculated based on the measured road beat signal will have a larger estimation error. Therefore, in order to reduce the calculation error of the initial flight time and reduce the number of iterations, Use the FM nonlinear term ⁇ (t) of the laser emission signal to compensate the FM nonlinear term of the transmitter on the measuring path beat signal to obtain the beat signal after FM nonlinear term compensation at the transmitter. Only the beat signal remains The frequency modulation nonlinear term corresponding to the echo signal is calculated based on the beat signal to obtain the initial flight time ⁇ with a smaller estimation error.
  • the frequency modulation nonlinear term ⁇ (t) of the laser emission signal is calculated according to the reference path beat signal, including:
  • an embodiment of the present application provides a ranging device based on frequency modulation nonlinear correction, including:
  • the acquisition unit is used to acquire the reference road beat signal and the measurement road beat signal, the reference road beat signal and the measurement road beat signal are respectively obtained based on the laser emission signal;
  • the calculation unit is used to calculate the frequency modulation nonlinearity according to the reference road beat frequency signal; and calculate the initial flight time ⁇ according to the measured road beat frequency signal;
  • the calculation unit is further configured to calculate the frequency modulation nonlinearity term ⁇ (t) of the laser emission signal according to the reference path beat signal when the determining unit determines that the frequency modulation nonlinearity is not less than the first preset threshold,
  • the correction unit is used to perform nonlinear iterative correction of the measurement path beat frequency signal according to the initial flight time ⁇ and the frequency modulation nonlinear term ⁇ (t) of the laser emission signal to obtain the target flight time;
  • the calculation unit is used to calculate the target distance according to the target flight time.
  • the influence of the frequency modulation nonlinearity on the beat signal of the measurement circuit can be eliminated to a large extent.
  • the influence minimizes the deviation between the last estimated target flight time and the target actual flight time, and the target distance can be determined according to the target flight time, which improves the ranging accuracy of the lidar.
  • the correction unit is specifically used for:
  • the correction unit in terms of determining the target flight time according to the corrected beat signal, is specifically configured to:
  • peak i-1 is the initial frequency spectrum peak
  • flight time ⁇ i-1 is the initial flight time ⁇ .
  • the predetermined condition peak i -peak i-1) / peak i-1 ⁇ Thr rate, Thr rate as the third predetermined threshold value.
  • the correction unit is specifically configured to include:
  • the frequency modulation nonlinear term ⁇ (t) of the laser emission signal is delayed according to the delay flight time ⁇ 'i-1 to obtain the frequency modulation corresponding to the echo signal of the measurement path
  • the non-linear term ⁇ (t- ⁇ 'i -1 ); the delayed flight time ⁇ 'i -1 is obtained based on the flight time ⁇ i -1 ;
  • the calculation unit is specifically configured to:
  • the corrected beat frequency signal includes the first corrected beat frequency signal or the second corrected beat frequency signal; according to the frequency modulation nonlinearity term ⁇ (t- ⁇ ' i-1 ) corresponding to the echo signal of the measurement path Perform FM nonlinear correction on the beat frequency signal of the measurement path to obtain the aspect of correcting the beat frequency signal.
  • the correction unit is specifically used for:
  • the frequency modulation nonlinearity term ⁇ (t) of the laser emission signal and the frequency modulation nonlinearity term ⁇ (t- ⁇ ) corresponding to the echo signal of the measurement path ' I-1 ) Perform FM nonlinear term correction on the beat signal of the measurement road to obtain the first corrected beat signal; when the signal-to-noise ratio of the beat signal of the measurement road is lower than the second preset threshold, according to the measurement road
  • the FM nonlinear term ⁇ (t- ⁇ ′i -1 ) corresponding to the wave signal performs receiving end FM nonlinear term correction on the compensated beat signal to obtain a second corrected beat signal;
  • the compensated beat signal is obtained based on the FM nonlinear term ⁇ (t) of the laser emitting signal to compensate the FM nonlinear term of the transmitter of the measurement path beat signal.
  • the calculation unit is specifically configured to:
  • the signal-to-noise ratio of the measured road beat signal is not lower than the second preset threshold, perform FFT on the measured road beat signal to obtain the frequency domain signal corresponding to the measured road beat signal;
  • the frequency domain signal determines the frequency of the measured road beat signal, and calculates the initial flight time ⁇ according to the frequency of the measured road beat signal;
  • the frequency of the measurement road beat signal is the frequency corresponding to the peak point position of the frequency domain signal corresponding to the measurement road beat signal, or the frequency corresponding to the midpoint position of the half-height width of the frequency domain signal.
  • the calculation unit is specifically configured to:
  • the measurement road beat signal is compensated by the FM nonlinear term at the transmitting end.
  • the frequency modulation nonlinear term of the laser emission signal is calculated based on the beat frequency signal of the reference path;
  • the frequency of the compensated beat signal is the frequency corresponding to the peak point position of the frequency domain signal corresponding to the compensated beat signal, or the frequency corresponding to the midpoint position of the half-height width of the frequency domain signal.
  • the calculation unit in terms of calculating the frequency modulation nonlinear term ⁇ (t) of the laser emission signal according to the reference path beat signal, is specifically configured to:
  • an embodiment of the present application provides a ranging device based on frequency modulation nonlinear correction, including:
  • Memory for storing instructions
  • At least one processor coupled with the memory
  • the instruction when the at least one processor executes the instruction, the instruction causes the processor to execute all or part of the method shown in the first aspect.
  • an embodiment of the present application provides a computer storage medium that stores a computer program.
  • the computer program includes program instructions that, when executed by a processor, cause the processor to execute All or part of the method shown on the one hand.
  • an embodiment of the present application provides a laser radar, including a laser and a processor coupled with the laser, and the processor executes all or part of the method shown in the first aspect.
  • Figure 1 is a schematic diagram of the FMCW lidar ranging principle
  • FIG. 2 is a block diagram of the system architecture of FMCW lidar
  • Figure 3 shows the beat signal under non-linear frequency modulation
  • FIG. 4 is a schematic diagram of a processing method based on modulated laser drive current predistortion in the prior art
  • FIG. 5 is a schematic flowchart of a ranging method based on frequency modulation nonlinear correction according to an embodiment of the application
  • FIG. 6 is a schematic flowchart of a non-linear iterative correction of a measurement road beat frequency signal provided by an embodiment of the application;
  • FIG. 7 is a schematic diagram of another flow chart for performing nonlinear iterative correction on measurement road beat frequency signals according to an embodiment of the application.
  • FIG. 8 is a schematic flowchart of another ranging method based on frequency modulation nonlinear correction provided by an embodiment of the application.
  • Figure 9a is a comparison diagram of ranging errors before and after nonlinear correction
  • Figure 9b is the frequency spectrum of the beat signal after iterative correction of the nonlinear term
  • FIG. 10 is a schematic structural diagram of a ranging device based on frequency modulation nonlinear correction provided by an embodiment of the application;
  • FIG. 11 is a schematic structural diagram of another ranging device based on frequency modulation nonlinear correction provided by an embodiment of the application.
  • the main problem of this scheme is that by adjusting the drive current to control the output frequency of the laser, the nonlinearity can only be reduced to a relatively small range, but it cannot be completely eliminated in principle.
  • the results of the current nonlinear correction can only be used for the next measurement. If due to the interference of external factors such as temperature and vibration, the output frequency of the laser in the next measurement is controlled by the drive current estimated by the predistortion in the last measurement, and the frequency value does not match the expected value, which leads to the previous predistortion The estimation result is invalid, and it fails to achieve the purpose of iterative correction.
  • the resampling algorithm is a method of FM nonlinear correction at the receiving end.
  • the core idea is to resample the beat frequency signal through the estimated sampling time sequence, so that the signal becomes a single frequency signal again.
  • the specific principle is as follows:
  • the premise of the re-sampling algorithm is that the delay ⁇ of the beat signal is a small value, and the higher-order term of ⁇ can be ignored.
  • the delay ⁇ corresponding to the beat signal does not meet the premise assumptions. , Resulting in the corresponding ranging error after resampling is still very large. Therefore, the performance of resampling is related to the distance to be measured, which has great practical limitations.
  • this application proposes a schematic flow chart of a ranging method based on frequency modulation nonlinear correction. As shown in Figure 5, the method includes:
  • S501 Obtain a reference road beat signal and a measurement road beat signal.
  • the laser signal emitted by the laser is divided into two paths, one path is the optical signal of the reference path, and the other path is the optical signal of the measurement path.
  • the frequency and phase of the optical signal of the reference path and the optical signal of the measurement path are the same.
  • the optical signal of the reference path is divided into two paths, and the reference path beat signal s ref (t) can be obtained by a Mach-Zehnder interferometer (MZI) with a preset fixed arm length difference of ⁇ D. Specifically, one optical signal passes through the short arm of the MZI, and the other optical signal passes through the long arm of the MZI. Because the two optical signals pass through different optical paths, the two optical signals are coherent on the photosensitive surface of the balanced detector to obtain a reference Road beat signal s ref (t).
  • MZI Mach-Zehnder interferometer
  • the optical signal of the measurement path is divided into two paths, one of which is a local oscillator signal, and the other is emitted to the measured target, and the echo signal is obtained by reflection on the surface of the target.
  • the echo signal and the local oscillator signal are coherent on the photosensitive surface of the balanced detector to obtain the measurement road beat signal s b (t).
  • S502 Calculate the frequency modulation nonlinearity according to the reference road beat frequency signal; and calculate the initial flight time ⁇ according to the measured road beat frequency signal.
  • calculating the frequency modulation nonlinearity according to the reference road beat signal includes:
  • H b (t) hilbert(s ref (t))
  • real(H b (t) ) Is the real part of the transformed signal H b (t)
  • imag(H b (t)) is the imaginary part of the transformed signal H b (t).
  • the non-linearity of the frequency modulation f ideal (t) is the ideal chirp frequency
  • f ideal (t) f 0 t+0.5 ⁇ t 2
  • f 0 is the carrier frequency
  • is the corresponding FM slope when ideal chirp
  • B is the FM bandwidth.
  • the frequency modulation nonlinearity can also be based on the root mean square sum of the frequency modulation nonlinearity term of the laser emission signal (that is, the difference between the frequency modulation frequency f(t) of the laser emission signal and the ideal chirp frequency f ideal (t)).
  • the frequency modulation bandwidth (B) is calculated, and the frequency modulation nonlinearity can be expressed as:
  • the degree of frequency modulation nonlinearity of the laser can be measured according to the degree of deviation between the frequency modulation frequency f(t) of the laser emission signal and the ideal chirp frequency f ideal (t), that is, f(t) and f ideal (t) Root Mean Squard Error (RMSE) between
  • both the FM frequency f(t) and the ideal chirp frequency f ideal (t) can be regarded as discrete signals, so f k is the frequency corresponding to the k-th discrete point on the FM frequency f(t), f ideal, k is the frequency corresponding to the k-th discrete point on the ideal chirp frequency f ideal (t).
  • the frequency modulation nonlinearity term ⁇ (t) of the laser emission signal is calculated according to the reference path beat frequency signal, and the frequency modulation nonlinearity term ⁇ (t) of the laser emission signal is calculated according to the initial flight time ⁇ and the frequency modulation nonlinearity of the laser emission signal
  • the term ⁇ (t) performs nonlinear iterative correction on the measured road beat signal to obtain the target flight time.
  • the first preset threshold Thr nonRate is the maximum non-linearity of the laser that can be received by the FMCW ranging system.
  • the frequency modulation nonlinear term ⁇ (t) of the laser emission signal is obtained based on the following formula:
  • calculating the initial flight time ⁇ according to the measured road beat signal s b (t) includes:
  • beat signal s compensated 'b (t) of frequency f' b is a beat signal s compensated 'peak point b (t) corresponding to the frequency domain signal corresponding to the position of the frequency or the frequency The frequency corresponding to the center of the half-height width of the domain signal.
  • the non-linear iterative correction is performed according to the initial flight time ⁇ and the frequency modulation non-linear term ⁇ (t) of the laser emission signal to obtain the target flight time, including:
  • determining the target flight time according to the corrected beat frequency signal includes:
  • peak i-1 is the initial spectral peak
  • flight time ⁇ i-1 is the initial flight time ⁇ .
  • the preset condition is (peak i -peak i-1 )/peak i-1 ⁇ Thr rate , and Thr rate is the third preset threshold.
  • performing nonlinear iterative correction on the measurement path beat frequency signal according to the initial flight time ⁇ and the frequency modulation nonlinear term ⁇ (t) of the laser emission signal to obtain the target flight time includes:
  • the frequency modulation nonlinear term ⁇ (t) of the laser emission signal is delayed according to the delay flight time ⁇ ′ i-1 to obtain the measurement path echo
  • obtaining the flight time ⁇ i according to the corrected beat signal includes:
  • the correction signal comprises a beat frequency signal a first correction or the second correction beat beat signal; an echo signal corresponding to the measured channel FM nonlinear term ⁇ (t- ⁇ 'i- 1) of Measure the beat signal of the road and perform FM nonlinear correction to obtain the corrected beat signal, including:
  • the frequency modulation nonlinearity term ⁇ (t) of the laser emission signal and the frequency modulation nonlinearity term ⁇ (t- ⁇ ) corresponding to the echo signal of the measurement path ' I-1 ) Perform frequency modulation nonlinear term correction on the beat frequency signal of the measurement path to obtain the first corrected beat frequency signal;
  • the compensated beat frequency signal ⁇ (t- ⁇ 'i-1) corresponding to the frequency modulation nonlinearity of the echo signal of the measurement path Performing correction of the non-linear term of the FM at the receiving end to obtain the second corrected beat signal;
  • the compensated beat signal is based on the FM nonlinear term ⁇ (t) of the laser emission signal to compensate for the FM nonlinear term at the transmitting end of the beat signal of the measuring path, so the compensated beat signal can also be called It is a beat signal that has been compensated for FM nonlinearity at the transmitter.
  • the signal-to-noise ratio of the measured road beat signal s b (t) is not lower than the second preset threshold, according to the frequency modulation nonlinear term ⁇ (t) of the laser emission signal and the frequency modulation nonlinear term of the echo signal ⁇ (t- ⁇ i ) Perform frequency modulation nonlinear term correction on the measurement path beat signal s b (t) to obtain the first corrected beat signal
  • the FM nonlinearity term ⁇ (t- ⁇ i ) of the measured road echo signal is performed on the transmitter.
  • the delayed flight time ⁇ ′ i-1 ⁇ i-1 +p, where p is the preset step length.
  • the delayed flight time ⁇ ′ i-1 can also be obtained based on other methods.
  • the FM nonlinear term ⁇ (t) of the laser emission signal is delayed by ⁇ ′ i-1 to obtain the FM nonlinear term ⁇ (t- ⁇ ′ i-1 ) of the echo signal of the measurement path.
  • the signal-to-noise ratio of the measured road beat signal s b (t) is not lower than the second preset threshold, according to the frequency modulation non-linear term ⁇ (t) of the laser emission signal and the frequency modulation non-linearity of the echo signal of the measurement road
  • the term ⁇ (t- ⁇ 'i -1 ) performs FM nonlinear term correction on the measurement path beat signal s b (t) to obtain the first corrected beat signal
  • the first corrected beat signal is: s b (t) is the initial beat signal of the measurement path, in other words, s b (t) is the corresponding beat signal before the iterative correction of the nonlinear term of the measurement path.
  • the transmitted frequency modulation is performed according to the frequency modulation non-linear term ⁇ (t- ⁇ ' i-1 ) of the measured road echo signal the beat signal s nonlinearity compensation 'b (t) FM reception terminal nonlinear term correction, the second correction to obtain a beat signal Among them, the second corrected beat signal is: The first corrected beat signal Or the second corrected beat signal Perform FFT to obtain the first corrected beat signal Or the second corrected beat signal The corresponding frequency domain signal is used to obtain the flight time ⁇ i and the spectral peak value peak i based on the frequency domain signal.
  • the frequency of the first corrected beat signal is determined according to the frequency domain signal corresponding to the first corrected beat signal, and the flight time ⁇ i is determined according to the frequency of the first corrected beat signal; or corresponding to the second corrected beat signal Determine the frequency of the second corrected beat signal from the frequency domain signal, and determine the flight time ⁇ i according to the frequency of the second corrected beat signal.
  • the frequency of the first corrected beat signal is the frequency corresponding to the peak i position of the frequency domain signal corresponding to the first corrected beat signal, or the frequency corresponding to the center position of the half-height width of the frequency domain signal;
  • the second corrected beat frequency The frequency of the signal is the frequency corresponding to the peak i position of the frequency domain signal corresponding to the second corrected beat frequency signal, or the frequency corresponding to the center position of the half-height width of the frequency domain signal.
  • Thr rate is a third preset threshold
  • the third preset threshold is a peak growth rate threshold.
  • peak i-1 is the initial peak of the spectrum
  • flight time ⁇ i-1 is the initial flight time ⁇ .
  • the initial spectral peak value is the peak value of the frequency domain signal corresponding to the measured road beat signal s b (t); when the signal-to-noise ratio of the measured road beat signal s b (t) is lower than the second preset threshold, Sign initial spectrum peak frequency compensated signal s peak 'b (t) corresponding to a frequency domain signal.
  • the delayed flight time ⁇ 'i-1 according to the flight time ⁇ i-1 .
  • the delayed flight time ⁇ ′ i-1 ⁇ i-1 +p, where p is the preset step length.
  • the delayed flight time ⁇ ′ i-1 can also be obtained based on other methods.
  • the first corrected beat signal is: s b (t) is the measurement road beat frequency signal.
  • the second corrected beat signal is:
  • the frequency of the first corrected beat signal is determined according to the frequency domain signal corresponding to the first corrected beat signal, and the flight time ⁇ i is determined according to the frequency of the first corrected beat signal; or corresponding to the second corrected beat signal Determine the frequency of the second corrected beat signal from the frequency domain signal, and determine the flight time ⁇ i according to the frequency of the second corrected beat signal.
  • the frequency of the first corrected beat signal is the frequency corresponding to the peak i position of the frequency domain signal corresponding to the first corrected beat signal, or the frequency corresponding to the center position of the half-height width of the frequency domain signal;
  • the second corrected beat frequency The frequency of the signal is the frequency corresponding to the peak i position of the frequency domain signal corresponding to the second corrected beat frequency signal, or the frequency corresponding to the center position of the half-height width of the frequency domain signal.
  • Target flight time flight time ⁇ i .
  • S504 Calculate the target distance according to the target flight time.
  • the target distance R cT/2
  • c is the speed of light
  • T is the target flight time
  • the beat signal of the reference path is obtained, and the optical signal of the reference path is specifically divided into two paths.
  • the reference path can be obtained by presetting the MZI with the fixed arm length difference ⁇ D.
  • the beat signal s ref (t) One of the optical signals passes through the short arm of the MZI, and the other optical signal passes through the long arm of the MZI. Because the two optical signals pass through different optical paths, the two optical signals are coherent on the photosensitive surface of the balanced detector to obtain the reference path beat frequency
  • the frequency modulation frequency f(t) of the laser emission signal is calculated according to the reference path beat frequency signal s ref (t).
  • H b, imag is the imaginary part of H b (t)
  • H b, real is the real part of H b (t).
  • the frequency modulation nonlinearity term ⁇ (t) and the frequency modulation nonlinearity non rate of the laser emission signal are calculated. Specifically, according to the instantaneous phase of the frequency modulation frequency f(t) of the laser emission signal and the beat frequency signal of the reference path Calculate the FM frequency f(t) from the relationship, where, Then calculate the frequency modulation nonlinearity term ⁇ (t) and the frequency modulation nonlinearity non rate of the laser emission signal according to the frequency modulation frequency f(t).
  • f ideal (t) is the ideal chirp frequency
  • f ideal (t) f 0 t+0.5 ⁇ t 2
  • f 0 is the carrier frequency
  • is the corresponding FM slope when ideal chirp
  • B is the FM bandwidth.
  • the frequency modulation nonlinearity can also be based on the root mean square sum of the frequency modulation nonlinearity term of the laser emission signal (that is, the difference between the frequency modulation frequency f(t) of the laser emission signal and the ideal chirp frequency f ideal (t)).
  • the frequency modulation bandwidth (B) is calculated, and the frequency modulation nonlinearity can be expressed as:
  • the degree of frequency modulation nonlinearity of the laser can be measured according to the degree of deviation between the frequency modulation frequency f(t) of the laser emission signal and the ideal chirp frequency f ideal (t), that is, f(t) and f ideal (t) Root Mean Squard Error (RMSE) between
  • both the FM frequency f(t) and the ideal chirp frequency f ideal (t) can be regarded as discrete signals, so f k is the frequency corresponding to the k-th discrete point on the FM frequency f(t), f ideal, k is the frequency corresponding to the k-th discrete point on the ideal chirp frequency f ideal (t).
  • the beat signal of the measurement path specifically, divide the optical signal of the measurement path into two paths, where one signal is used as a local oscillator signal, and the other signal is transmitted to the measured target, and the echo signal is obtained by reflection on the target surface.
  • the echo signal and the local oscillator signal are coherent on the photosensitive surface of the balanced detector to obtain the measurement road beat signal s b (t).
  • the FM nonlinear term is corrected according to the FM nonlinear term ⁇ (t) and the FM nonlinear term ⁇ (t- ⁇ ) on the measurement path beat signal s b (t) to obtain the first corrected beat signal
  • the first corrected beat signal Or the second corrected beat signal Perform FFT to get the first corrected beat signal Or the second corrected beat signal Corresponding frequency domain signal, calculate the first corrected beat frequency signal based on the frequency domain signal Or the second corrected beat signal The frequency and spectral peak peakcurr, according to the first corrected beat signal Or the second corrected beat signal The frequency is calculated to get the flight time ⁇ .
  • Figure 9a shows that the FM nonlinearity correction method of the present application effectively improves the accuracy of distance estimation in the case of FM nonlinearity. It can be seen from Fig. 9b that after the iterative correction of the nonlinear term of this application, the spectral peak of the beat signal gradually increases and is close to the ideal value. Therefore, the nonlinear correction method of this application can effectively increase the beat frequency in the case of frequency modulation nonlinearity. The signal-to-noise ratio in the frequency domain of the signal can effectively increase the farthest ranging range.
  • the frequency modulation nonlinearity of the measurement path signal can be corrected according to the frequency modulation nonlinearity term of the laser emission signal and the frequency modulation nonlinearity term of the measurement path echo signal. Eliminate the influence of non-linear frequency modulation on the measured road beat signal, improve the accuracy of the flight time based on the corrected beat signal, and further improve the ranging accuracy.
  • the distance measuring device 1000 includes:
  • the acquiring unit 1001 is used to acquire the reference road beat signal and the measurement road beat signal, the reference road beat signal and the measurement road beat signal are obtained based on the laser emission signal respectively;
  • the calculation unit 1002 is configured to calculate the frequency modulation nonlinearity according to the reference road beat signal; and calculate the initial flight time ⁇ according to the measured road beat signal;
  • the calculation unit 1002 is further configured to calculate the frequency modulation nonlinearity term ⁇ (t) of the laser emission signal according to the reference path beat signal when the determining unit 1004 determines that the frequency modulation nonlinearity is not less than the first preset threshold,
  • the correction unit 1003 is configured to perform nonlinear iterative correction on the measurement path beat frequency signal according to the initial flight time ⁇ and the frequency modulation nonlinear term ⁇ (t) of the laser emission signal to obtain the target flight time;
  • the calculation unit 1002 is configured to calculate the target distance according to the target flight time.
  • correction unit 1003 is specifically configured to:
  • the correction unit 1003 is specifically configured to:
  • peak i-1 is the initial frequency spectrum peak
  • flight time ⁇ i-1 is the initial flight time ⁇ .
  • the predetermined condition peak i -peak i-1) / peak i-1 ⁇ Thr rate, Thr rate as the third predetermined threshold value.
  • the correction unit 1003 is specifically configured to include:
  • the frequency modulation nonlinear term ⁇ (t) of the laser emission signal is delayed according to the delay flight time ⁇ 'i-1 to obtain the frequency modulation corresponding to the echo signal of the measurement path
  • the non-linear term ⁇ (t- ⁇ 'i -1 ); the delayed flight time ⁇ 'i -1 is obtained based on the flight time ⁇ i -1 ;
  • the calculation unit 1002 is specifically configured to:
  • the corrected beat frequency signal includes the first corrected beat frequency signal or the second corrected beat frequency signal; according to the frequency modulation nonlinearity term ⁇ (t- ⁇ ' i-1 ) corresponding to the echo signal of the measurement path Performing FM nonlinear correction on the beat signal of the measurement path to obtain the aspect of correcting the beat signal, the correction unit 1003 is specifically used for:
  • the frequency modulation nonlinearity term ⁇ (t) of the laser emission signal and the frequency modulation nonlinearity term ⁇ (t- ⁇ ) corresponding to the echo signal of the measurement path ' I-1 ) Perform FM nonlinear term correction on the beat signal of the measurement road to obtain the first corrected beat signal; when the signal-to-noise ratio of the beat signal of the measurement road is lower than the second preset threshold, according to the measurement road
  • the FM nonlinear term ⁇ (t- ⁇ ′i -1 ) corresponding to the wave signal performs receiving end FM nonlinear term correction on the compensated beat signal to obtain a second corrected beat signal;
  • the compensated beat signal is obtained based on the FM nonlinear term ⁇ (t) of the laser emitting signal to compensate the FM nonlinear term of the transmitter of the measurement path beat signal.
  • the calculation unit 1002 is specifically configured to:
  • the signal-to-noise ratio of the measured road beat signal is not lower than the second preset threshold, perform FFT on the measured road beat signal to obtain the frequency domain signal corresponding to the measured road beat signal;
  • the frequency domain signal determines the frequency of the measured road beat signal, and calculates the initial flight time ⁇ according to the frequency of the measured road beat signal;
  • the frequency of the measurement road beat signal is the frequency corresponding to the peak point position of the frequency domain signal corresponding to the measurement road beat signal, or the frequency corresponding to the midpoint position of the half-height width of the frequency domain signal.
  • the calculation unit 1002 is specifically configured to:
  • the measurement road beat signal is compensated by the FM nonlinear term at the transmitting end.
  • the frequency modulation nonlinear term of the laser emission signal is calculated based on the beat frequency signal of the reference path;
  • the frequency of the compensated beat signal is the frequency corresponding to the peak point position of the frequency domain signal corresponding to the compensated beat signal, or the frequency corresponding to the midpoint position of the half-height width of the frequency domain signal.
  • calculation unit 1002 is specifically configured to:
  • the aforementioned units are used to execute the relevant steps of the aforementioned method.
  • the acquiring unit 1001 is used to execute the related content of step S501
  • the calculating unit 1002 is used to execute the related content of steps S502 and S504
  • the correction unit 1003 and the determining unit 1004 are used to execute the related content of step S503.
  • the ranging device 1000 based on frequency modulation nonlinear correction is presented in the form of a unit.
  • the "unit” here can refer to an application-specific integrated circuit (ASIC), a processor and memory that executes one or more software or firmware programs, an integrated logic circuit, and/or other devices that can provide the above-mentioned functions .
  • ASIC application-specific integrated circuit
  • the above acquisition unit 1001, calculation unit 1002, correction unit 1003, and determination unit 1004 may be implemented by the processor 1101 of the distance measurement device based on frequency modulation nonlinear correction shown in FIG. 11.
  • the distance measuring device 1100 shown in FIG. 11 can be implemented with the structure in FIG. 11, and the distance measuring device 1100 includes at least one processor 1101, at least one memory 1102 and at least one communication interface 1103.
  • the processor 1101, the memory 1102, and the communication interface 1103 are connected through the communication bus and complete mutual communication.
  • the processor 1101 may be a general-purpose central processing unit (CPU), a microprocessor, an application-specific integrated circuit (ASIC), or one or more integrated circuits for controlling the execution of the above program programs.
  • CPU central processing unit
  • ASIC application-specific integrated circuit
  • the communication interface 1103 is used to communicate with other devices or communication networks, such as Ethernet, radio access network (RAN), wireless local area network (Wireless Local Area Networks, WLAN), etc.
  • RAN radio access network
  • WLAN Wireless Local Area Networks
  • the memory 1102 may be a read-only memory (ROM) or other types of static storage devices that can store static information and instructions, random access memory (RAM), or other types that can store information and instructions
  • the dynamic storage device can also be electrically erasable programmable read-only memory (Electrically Erasable Programmable Read-Only Memory, EEPROM), CD-ROM (Compact Disc Read-Only Memory, CD-ROM) or other optical disc storage, optical disc storage (Including compact discs, laser discs, optical discs, digital versatile discs, Blu-ray discs, etc.), magnetic disk storage media or other magnetic storage devices, or can be used to carry or store desired program codes in the form of instructions or data structures and can be used by a computer Any other media accessed, but not limited to this.
  • the memory can exist independently and is connected to the processor through a bus.
  • the memory can also be integrated with the processor.
  • the memory 1102 is used to store application program codes for executing the above solutions, and the processor 1101 controls the execution.
  • the processor 1101 is configured to execute application program codes stored in the memory 1102.
  • the code stored in the memory 1102 can execute the above-provided ranging method based on FM nonlinear correction, such as obtaining the reference road beat signal and measuring road beat signal; calculating the FM nonlinearity based on the reference road beat signal; and according to the measurement
  • the initial flight time ⁇ is calculated from the road beat frequency signal; when the frequency modulation nonlinearity is not less than the first preset threshold, the frequency modulation nonlinear term ⁇ (t) of the laser emission signal is calculated according to the reference road beat frequency signal, according to the initial flight time ⁇ and the frequency modulation non-linear term ⁇ (t) of the laser emission signal perform non-linear iterative correction on the measurement road beat signal to obtain the target flight time; the target distance is calculated according to the target flight time.
  • An embodiment of the present application further provides a computer storage medium, wherein the computer storage medium may store a program, and when the program is executed, it includes any part or All steps.
  • the disclosed device may be implemented in other ways.
  • the device embodiments described above are merely illustrative, for example, the division of the units is only a logical function division, and there may be other divisions in actual implementation, for example, multiple units or components may be combined or may be Integrate into another system, or some features can be ignored or not implemented.
  • the displayed or discussed mutual coupling or direct coupling or communication connection may be indirect coupling or communication connection through some interfaces, devices or units, and may be in electrical or other forms.
  • the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the objectives of the solutions of the embodiments.
  • the functional units in the various embodiments of the present application may be integrated into one processing unit, or each unit may exist alone physically, or two or more units may be integrated into one unit.
  • the above-mentioned integrated unit can be implemented in the form of hardware or software functional unit.
  • the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer readable memory.
  • the technical solution of the present application essentially or the part that contributes to the existing technology or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a memory.
  • a number of instructions are included to enable a computer device (which may be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods described in the various embodiments of the present application.
  • the aforementioned memory includes: U disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), mobile hard disk, magnetic disk or optical disk and other media that can store program codes.
  • the program can be stored in a computer-readable memory, and the memory can include: a flash disk , Read-only memory (English: Read-Only Memory, abbreviation: ROM), random access device (English: Random Access Memory, abbreviation: RAM), magnetic disk or optical disc, etc.

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Optical Radar Systems And Details Thereof (AREA)

Abstract

A frequency modulation nonlinear correction-based ranging method, comprising: acquiring a reference path beat frequency signal and a measurement path beat frequency signal (S501); calculating frequency modulation nonlinearity according to the reference path beat frequency signal, and calculating an initial flight time τ according to the measurement path beat frequency signal (S502); when the frequency modulation nonlinearity is not less than a first preset threshold, calculating a frequency modulation nonlinearity term ε(t) of a laser emission signal according to the reference path beat frequency signal, and performing nonlinear iterative correction on the measurement path beat frequency signal according to the initial flight time τ and the frequency modulation nonlinear term ε(t) of the laser emission signal so as to obtain target flight time (S503); and calculating a target distance according to the target flight time (S504). The ranging accuracy and range of lidars may be improved.

Description

基于调频非线性校正的测距方法及相关装置Ranging method and related device based on frequency modulation nonlinear correction 技术领域Technical field
本申请涉及激光雷达领域,具体涉及一种基于调频非线性校正的测距方法及相关装置。This application relates to the field of lidar, and in particular to a ranging method and related devices based on frequency modulation nonlinear correction.
背景技术Background technique
调频连续波(Frequency Modulated Continuous,FMCW)激光雷达(Light Detection and Ranging,LiDAR)是一种光学遥感技术,它能完成测距、测速、目标探测、跟踪和成像识别等功能,可应用于智能交通、自动驾驶、大气环境监测、地理测绘、无人机等领域。Frequency Modulated Continuous Wave (Frequency Modulated Continuous, FMCW) Lidar (Light Detection and Ranging, LiDAR) is an optical remote sensing technology that can complete functions such as ranging, speed measurement, target detection, tracking and imaging recognition, and can be applied to intelligent transportation , Autonomous driving, atmospheric environment monitoring, geographic surveying and mapping, unmanned aerial vehicles and other fields.
FMCW LiDAR是通过测量发射光信号与回波信号相干得到拍频信号的频域响应来获取被测目标的位置信息。基本过程如图1所示。通过对拍频信号进行频域分析,估计发射信号与回波信号之间的频率差f b,进而根据频率差f b计算飞行时间τ,最后根据飞行时间计算得到被测目标的距离。 FMCW LiDAR obtains the position information of the measured target by measuring the coherence of the emitted light signal and the echo signal to obtain the frequency domain response of the beat signal. The basic process is shown in Figure 1. Through frequency domain analysis of the beat signal, the frequency difference f b between the transmitted signal and the echo signal is estimated, and the flight time τ is calculated according to the frequency difference f b , and finally the distance of the measured target is calculated according to the flight time.
图2为典型的FMCW LiDAR系统架构框图。通过向可调谐激光器注入适当驱动电流,驱动激光器发射激光信号,光束经过隔离器后被耦合器分为两束激光,其中一束光束作为参考光,另外一束光束先后经过环形器和镜头出射至被测目标物,称为测量光。测量光经目标表面反射后被镜头收集并经过环形器与参考光经过耦合器,在平衡探测器的光敏面发生相干,产生频率与飞行时间τ成正比的拍频信号;通过对该拍频信号进行频域分析,得到拍频信号频率值f b,也即发射信号与回波信号的频率差;再根据调频周期和调频带宽等已知参数,便可计算出飞行时间τ,进而估计目标距离。 Figure 2 is a block diagram of a typical FMCW LiDAR system architecture. By injecting appropriate driving current into the tunable laser, the laser is driven to emit a laser signal. After the beam passes through the isolator, it is divided into two laser beams by the coupler. The measured target is called the measuring light. The measuring light is reflected by the target surface and collected by the lens and passed through the circulator. The reference light passes through the coupler. Coherence occurs on the photosensitive surface of the balanced detector to generate a beat signal whose frequency is proportional to the flight time τ; through the beat signal Perform frequency domain analysis to obtain the beat signal frequency f b , that is, the frequency difference between the transmitted signal and the echo signal; then according to known parameters such as the frequency modulation period and the frequency modulation bandwidth, the flight time τ can be calculated, and then the target distance can be estimated .
在理想线性调频情况下,发射信号与回波信号的调频曲线应该是一致的且频率随着时间线性变化,二者只存在飞行延时,故发射信号的调频线性度会很大程度上影响FMCWLiDAR的测距精度。实际情况中,激光器的输出频率是通过改变激光器的驱动电流来控制的,由于半导体激光器本身固有的频率调制效应,激光器的输出频率对于调制电流呈现明显的非线性,如图3所示。对比于线性调频下,当存在调频非线性时,拍频信号不再是单频信号,其拍频频率会随时间变化,频谱会被展宽,导致拍频信号的信噪比下降,进而严重影响到测距精度和测距范围。In an ideal chirp condition, the frequency modulation curves of the transmitted signal and the echo signal should be consistent and the frequency changes linearly with time. There is only a flight delay between the two, so the frequency modulation linearity of the transmitted signal will greatly affect the FMCWLiDAR The accuracy of ranging. In actual situations, the output frequency of the laser is controlled by changing the drive current of the laser. Due to the inherent frequency modulation effect of the semiconductor laser, the output frequency of the laser exhibits obvious nonlinearity to the modulation current, as shown in Figure 3. Compared with linear frequency modulation, when there is FM nonlinearity, the beat signal is no longer a single-frequency signal, its beat frequency will change with time, and the frequency spectrum will be broadened, resulting in a decrease in the signal-to-noise ratio of the beat signal, which will seriously affect To the ranging accuracy and ranging range.
发明内容Summary of the invention
本申请实施例提供了一种基于调频非线性校正的测距方法及相关装置,以提高激光雷达的测距精度与测距范围。The embodiments of the present application provide a ranging method and related devices based on frequency modulation nonlinear correction, so as to improve the ranging accuracy and ranging range of the lidar.
第一方面,本申请实施例提供了一种基于调频非线性校正的测距方法,包括:In the first aspect, an embodiment of the present application provides a ranging method based on frequency modulation nonlinear correction, including:
获取参考路拍频信号和测量路拍频信号,参考路拍频信号和测量路拍频信号是分别基于激光器发射信号得到的;根据参考路拍频信号计算得到调频非线性度;并根据测量路拍频信号计算得到初始飞行时间τ;确定调频非线性度不小于第一预设阈值,根据参考路拍频信号计算得到激光器发射信号的调频非线性项ε(t),根据初始飞行时间τ和激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行非线性迭代校正,以得到目标飞行时间;根据目标飞行时间计算得到目标距离。Obtain the reference road beat signal and the measurement road beat signal. The reference road beat signal and the measurement road beat signal are respectively obtained based on the laser emission signal; the frequency modulation nonlinearity is calculated according to the reference road beat signal; and according to the measurement road The beat frequency signal is calculated to obtain the initial flight time τ; it is determined that the frequency modulation nonlinearity is not less than the first preset threshold, and the frequency modulation nonlinearity term ε(t) of the laser emission signal is calculated according to the reference path beat frequency signal. According to the initial flight time τ and The frequency modulation nonlinear term ε(t) of the laser emission signal performs nonlinear iterative correction on the measurement road beat signal to obtain the target flight time; the target distance is calculated according to the target flight time.
通过基于飞行时间τ和所述激光器发射信号的调频非线性项ε(t)对测量路的拍频信号进行非线性迭代校正,消除了调频非线性对测量路拍频信号的影响,使得最后估计的目标飞行时间与目标实际飞行时间偏差最小,进而可根据目标飞行时间确定目标距离,提高了激光雷达的测距精度。The non-linear iterative correction of the beat signal of the measurement circuit is performed based on the flight time τ and the frequency modulation non-linear term ε(t) of the laser emission signal, which eliminates the influence of the non-linear frequency modulation on the beat signal of the measurement circuit and makes the final estimation The deviation between the target flight time and the actual flight time of the target is the smallest, and the target distance can be determined according to the target flight time, which improves the ranging accuracy of the lidar.
在一个可行的实施例中,根据初始飞行时间τ和激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行非线性迭代校正,以得到目标飞行时间,包括:In a feasible embodiment, the non-linear iterative correction of the measurement path beat frequency signal according to the initial flight time τ and the frequency modulation nonlinear term ε(t) of the laser emission signal to obtain the target flight time includes:
S1:根据延时飞行时间τ' i-1对激光器发射信号的调频非线性项ε(t)进行延时,以得到测量路回波信号对应的调频非线性项ε(t-τ' i-1),延时飞行时间τ' i-1是基于飞行时间τ i-1得到的; S1: Delay the FM nonlinear term ε(t) of the laser emission signal according to the delay flight time τ′ i-1 to obtain the FM nonlinear term ε(t-τ′ i-) corresponding to the echo signal of the measurement path. 1 ), the delayed flight time τ'i -1 is obtained based on the flight time τ i-1 ;
S2:根据测量路回波信号对应的调频非线性项ε(t-τ' i-1)对测量路拍频信号进行调频非线性校正,以得到校正拍频信号; S2: Perform FM nonlinear correction on the beat frequency signal of the measurement path according to the FM nonlinear term ε(t-τ'i -1) corresponding to the echo signal of the measurement path to obtain a corrected beat frequency signal;
S3、根据校正拍频信号确定目标飞行时间。S3. Determine the target flight time according to the corrected beat frequency signal.
通过根据激光器发射信号的调频非线性项对测量路拍频信号进行调频非线性校正,消除了调频非线性对测量路拍频信号的影响,提高了基于校正后的拍频信号得到的飞行时间的精确性,进而提高了测距精度。According to the frequency modulation nonlinear term of the laser emission signal, the frequency modulation nonlinearity correction is performed on the measurement road beat signal, which eliminates the influence of the frequency modulation nonlinearity on the measurement road beat signal, and improves the flight time based on the corrected beat signal. Accuracy, thereby improving the accuracy of ranging.
在一个可行的实施例中,根据校正拍频信号确定目标飞行时间,包括:In a feasible embodiment, determining the target flight time according to the corrected beat signal includes:
根据校正拍频信号获取飞行时间τ i;当频谱峰值peak i满足预设条件时,将飞行时间τ i确定为目标飞行时间;当频谱峰值peak i不满足预设条件时,令i=i+1,并重复执行上述步骤S1-S3;频谱峰值peak i为校正拍频信号对应的频域信号的峰值; Obtain the flight time τ i according to the corrected beat signal; when the peak frequency peak i meets the preset condition, the flight time τ i is determined as the target flight time; when the peak frequency peak i does not meet the preset condition, let i=i+ 1. Repeat the above steps S1-S3; the peak value of the spectrum peak i is the peak value of the frequency domain signal corresponding to the corrected beat signal;
当i=1时,peak i-1为所述初始频谱峰值,飞行时间τ i-1为初始飞行时间τ。 When i=1, peak i-1 is the initial frequency spectrum peak, and the flight time τ i-1 is the initial flight time τ.
通过对测量路拍频信号进行迭代校正,可在很大程度上消除调频非线性对测量路拍频信号的影响,提高了基于校正后的拍频信号得到的飞行时间的精确性,进而提高了测距精度。By iteratively correcting the beat signal of the measurement road, the influence of FM nonlinearity on the beat signal of the measurement road can be eliminated to a large extent, and the accuracy of the flight time obtained based on the corrected beat signal is improved, thereby improving Ranging accuracy.
在一个可行的实施例中,预设条件为(peak i-peak i-1)/peak i-1<Thr rate,该Thr rate为第三预设阈值。 In one possible embodiment, the predetermined condition (peak i -peak i-1) / peak i-1 <Thr rate, Thr rate as the third predetermined threshold value.
在一个可行的实施例中,根据初始飞行时间τ和激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行非线性迭代校正,以得到目标飞行时间,包括:In a feasible embodiment, the non-linear iterative correction of the measurement path beat frequency signal according to the initial flight time τ and the frequency modulation nonlinear term ε(t) of the laser emission signal to obtain the target flight time includes:
根据初始飞行时间τ和激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行M次非线性迭代计算,以得到目标飞行时间,M为大于1的整数,According to the initial flight time τ and the frequency modulation non-linear term ε(t) of the laser emission signal, perform M non-linear iterative calculations on the measurement path beat frequency signal to obtain the target flight time, M is an integer greater than 1,
其中,在进行第i次非线性迭代计算时,根据延时飞行时间τ' i-1对激光器发射信号的调频非线性项ε(t)进行延时,以得到测量路回波信号对应的调频非线性项ε(t-τ' i-1);延时飞行时间τ' i-1是基于飞行时间τ i-1得到的; Among them, when the i-th nonlinear iterative calculation is performed, the frequency modulation nonlinear term ε(t) of the laser emission signal is delayed according to the delay flight time τ'i-1 to obtain the frequency modulation corresponding to the echo signal of the measurement path The non-linear term ε(t-τ'i -1 ); the delayed flight time τ'i -1 is obtained based on the flight time τi -1 ;
根据测量路回波信号对应的调频非线性项ε(t-τ' i-1)对测量路拍频信号进行调频非线性校正,以得到校正拍频信号;根据校正拍频信号获取飞行时间τ iPerform FM nonlinear correction on the measurement path beat signal according to the FM nonlinear term ε(t-τ' i-1 ) corresponding to the echo signal of the measurement path to obtain the corrected beat signal; obtain the flight time τ according to the corrected beat signal i ;
其中,当i=1时,飞行时间τ i-1为初始飞行时间τ;当i=M时,目标飞行时间为飞行时间τ iAmong them, when i=1, the flight time τ i-1 is the initial flight time τ; when i=M, the target flight time is the flight time τ i .
通过对测量路拍频信号进行多次调频非线性校正,可在很大程度上消除调频非线性对测量路拍频信号的影响,提高了基于校正后的拍频信号得到的飞行时间的精确性,进而提 高了测距精度。By performing multiple FM nonlinear corrections on the measurement road beat signal, the influence of FM nonlinearity on the measurement road beat signal can be eliminated to a large extent, and the accuracy of the flight time obtained based on the corrected beat signal is improved. , Thereby improving the accuracy of ranging.
在一个可行的实施例中,根据校正拍频信号获取飞行时间τ i,包括: In a feasible embodiment, obtaining the flight time τ i according to the corrected beat signal includes:
对校正拍频信号进行FFT,以得到校正拍频信号对应的频域信号;根据校正拍频信号对应的频域信号获取飞行时间τ iPerform FFT on the corrected beat frequency signal to obtain the frequency domain signal corresponding to the corrected beat frequency signal; obtain the flight time τ i according to the frequency domain signal corresponding to the corrected beat frequency signal.
在一个可行的实施例中,校正拍频信号包括第一校正拍频信号或第二校正拍频信号;根据测量路回波信号对应的调频非线性项ε(t-τ' i-1)对测量路拍频信号进行调频非线性校正,以得到校正拍频信号,包括: In one possible embodiment, the correction signal comprises a beat frequency signal a first correction or the second correction beat beat signal; an echo signal corresponding to the measured channel FM nonlinear term ε (t-τ 'i- 1) of Measure the beat signal of the road and perform FM nonlinear correction to obtain the corrected beat signal, including:
当测量路拍频信号的信噪比不低于第二预设阈值时,根据激光器发射信号的调频非线性项ε(t)和测量路回波信号对应的调频非线性项ε(t-τ' i-1)对测量路拍频信号进行调频非线性项校正,以得到第一校正拍频信号;当测量路拍频信号的信噪比低于第二预设阈值时,根据测量路回波信号对应的调频非线性项ε(t-τ' i-1)对补偿后的拍频信号进行接收端调频非线性项校正,以得到第二校正拍频信号; When the signal-to-noise ratio of the beat frequency signal of the measurement path is not lower than the second preset threshold, the frequency modulation nonlinearity term ε(t) of the laser emission signal and the frequency modulation nonlinearity term ε(t-τ) corresponding to the echo signal of the measurement path ' I-1 ) Perform FM nonlinear term correction on the beat signal of the measurement road to obtain the first corrected beat signal; when the signal-to-noise ratio of the beat signal of the measurement road is lower than the second preset threshold, according to the measurement road The FM nonlinear term ε(t-τ′i -1 ) corresponding to the wave signal performs receiving end FM nonlinear term correction on the compensated beat signal to obtain a second corrected beat signal;
其中,补偿后的拍频信号是基于根据激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行发射端的调频非线性项补偿得到的。Among them, the compensated beat signal is obtained based on the FM nonlinear term ε(t) of the laser emitting signal to compensate the FM nonlinear term of the transmitter of the measurement path beat signal.
当测量路拍频信号的信噪比很低时,由于基于该测量路拍频信号计算得到的初始飞行时间会有较大的估计误差,因此为了降低初始飞行时间的计算误差和减少迭代次数,根据激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行了补偿,并得到了补偿后的拍频信号,从而在对测量路拍频信号进行校正时,只需要根据测量路回波信号的调频非线性项对补偿后的拍频信号进行调频非线性校正,也消除调频非线性对测量路拍频信号的影响。When the signal-to-noise ratio of the measured road beat signal is very low, since the initial flight time calculated based on the measured road beat signal will have a large estimation error, in order to reduce the calculation error of the initial flight time and reduce the number of iterations, According to the frequency modulation nonlinear term ε(t) of the laser emission signal, the beat signal of the measurement path is compensated, and the compensated beat signal is obtained. Therefore, when correcting the beat signal of the measurement path, it is only necessary to calibrate the beat signal of the measurement path. The FM nonlinear term of the echo signal performs FM nonlinear correction on the compensated beat signal, and also eliminates the influence of FM nonlinearity on the beat signal of the measurement path.
在一个可行的实施例中,根据测量路拍频信号计算得到初始飞行时间τ,包括:In a feasible embodiment, calculating the initial flight time τ according to the measured road beat frequency signal includes:
当测量路拍频信号的信噪比不低于第二预设阈值时,对测量路拍频信号进行FFT,以得到测量路拍频信号对应的频域信号;根据测量路拍频信号对应的频域信号确定测量路拍频信号的频率,并根据测量路拍频信号的频率计算得到初始飞行时间τ;When the signal-to-noise ratio of the measured road beat signal is not lower than the second preset threshold, perform FFT on the measured road beat signal to obtain the frequency domain signal corresponding to the measured road beat signal; The frequency domain signal determines the frequency of the measured road beat signal, and calculates the initial flight time τ according to the frequency of the measured road beat signal;
其中,测量路拍频信号的频率为测量路拍频信号对应的频域信号的峰值点位置对应的频率,或者频域信号的半高宽度中点位置对应的频率。Wherein, the frequency of the measurement road beat signal is the frequency corresponding to the peak point position of the frequency domain signal corresponding to the measurement road beat signal, or the frequency corresponding to the midpoint position of the half-height width of the frequency domain signal.
在一个可行的实施例中,根据测量路拍频信号计算得到初始飞行时间τ,包括:In a feasible embodiment, calculating the initial flight time τ according to the measured road beat frequency signal includes:
当测量路拍频信号的信噪比低于第二预设阈值时;根据激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行发射端的调频非线性项补偿,以得到补偿后的拍频信号;激光器发射信号的调频非线性项是基于参考路拍频信号计算得到的;When the signal-to-noise ratio of the measurement road beat signal is lower than the second preset threshold; according to the FM nonlinear term ε(t) of the laser emitting signal, the measurement road beat signal is compensated by the FM nonlinear term at the transmitting end. After the beat frequency signal; the frequency modulation nonlinear term of the laser emission signal is calculated based on the beat frequency signal of the reference path;
对补偿后的拍频信号进行FFT,以得到补偿后的拍频信号对应的频域信号;根据补偿后的拍频信号对应的频域信号确定补偿后的拍频信号的频率,并根据补偿后的拍频信号的频率计算得到初始飞行时间τ;Perform FFT on the compensated beat signal to obtain the frequency domain signal corresponding to the compensated beat signal; determine the frequency of the compensated beat signal according to the frequency domain signal corresponding to the compensated beat signal, and according to the compensated beat signal Calculate the frequency of the beat signal to get the initial flight time τ;
其中,补偿后的拍频信号的频率为补偿后的拍频信号对应的频域信号的峰值点位置对应的频率,或者频域信号的半高宽度中点位置对应的频率。The frequency of the compensated beat signal is the frequency corresponding to the peak point position of the frequency domain signal corresponding to the compensated beat signal, or the frequency corresponding to the midpoint position of the half-height width of the frequency domain signal.
当测量路拍频信号的信噪比很低时,基于该测量路拍频信号计算得到的初始飞行时间会有较大的估计误差,因此为了降低初始飞行时间的计算误差,并减少迭代次数,利用激光器发射信号的调频非线项ε(t)对测量路拍频信号进行发射端的调频非线性项补偿,以得到发射端调频非线性项补偿后的拍频信号,该拍频信号只剩下回波信号对应的调频非线性项, 基于该拍频信号计算得到估计误差较小的初始飞行时间τ。When the signal-to-noise ratio of the measured road beat signal is very low, the initial flight time calculated based on the measured road beat signal will have a larger estimation error. Therefore, in order to reduce the calculation error of the initial flight time and reduce the number of iterations, Use the FM nonlinear term ε(t) of the laser emission signal to compensate the FM nonlinear term of the transmitter on the measuring path beat signal to obtain the beat signal after FM nonlinear term compensation at the transmitter. Only the beat signal remains The frequency modulation nonlinear term corresponding to the echo signal is calculated based on the beat signal to obtain the initial flight time τ with a smaller estimation error.
在一个可行的实施例中,根据参考路拍频信号计算得到激光器发射信号的调频非线性项ε(t),包括:In a feasible embodiment, the frequency modulation nonlinear term ε(t) of the laser emission signal is calculated according to the reference path beat signal, including:
对参考路拍频信号进行希尔伯特变换,以得到变换后的拍频信号;根据变换后的拍频信号计算得到激光器发射信号的调频频率;根据激光器发射信号的调频频率和理想线性调频频率计算得到激光器发射信号的调频非线性项ε(t)。Perform Hilbert transform on the beat signal of the reference path to obtain the transformed beat signal; calculate the FM frequency of the laser emission signal according to the transformed beat signal; Calculate the FM frequency of the laser emission signal according to the FM frequency and ideal chirp frequency of the laser emission signal The frequency modulation nonlinear term ε(t) of the laser emission signal is calculated.
第二方面,本申请实施例提供一种基于调频非线性校正的测距装置,包括:In a second aspect, an embodiment of the present application provides a ranging device based on frequency modulation nonlinear correction, including:
获取单元,用于获取参考路拍频信号和测量路拍频信号,参考路拍频信号和测量路拍频信号是分别基于激光器发射信号得到的;The acquisition unit is used to acquire the reference road beat signal and the measurement road beat signal, the reference road beat signal and the measurement road beat signal are respectively obtained based on the laser emission signal;
计算单元,用于根据参考路拍频信号计算得到调频非线性度;并根据测量路拍频信号计算得到初始飞行时间τ;The calculation unit is used to calculate the frequency modulation nonlinearity according to the reference road beat frequency signal; and calculate the initial flight time τ according to the measured road beat frequency signal;
计算单元,还用于在确定单元确定调频非线性度不小于第一预设阈值时,根据参考路拍频信号计算得到激光器发射信号的调频非线性项ε(t),The calculation unit is further configured to calculate the frequency modulation nonlinearity term ε(t) of the laser emission signal according to the reference path beat signal when the determining unit determines that the frequency modulation nonlinearity is not less than the first preset threshold,
校正单元,用于根据初始飞行时间τ和激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行非线性迭代校正,以得到目标飞行时间;The correction unit is used to perform nonlinear iterative correction of the measurement path beat frequency signal according to the initial flight time τ and the frequency modulation nonlinear term ε(t) of the laser emission signal to obtain the target flight time;
计算单元,用于根据目标飞行时间计算得到目标距离。The calculation unit is used to calculate the target distance according to the target flight time.
通过基于飞行时间τ和所述激光器发射信号的调频非线性项ε(t)对测量路的拍频信号进行非线性迭代校正,可在很大程度上消除调频非线性对测量路拍频信号的影响,使得最后估计的目标飞行时间与目标实际飞行时间偏差最小,进而可根据目标飞行时间确定目标距离,提高了激光雷达的测距精度。By performing nonlinear iterative correction on the beat signal of the measurement circuit based on the flight time τ and the frequency modulation nonlinear term ε(t) of the laser emission signal, the influence of the frequency modulation nonlinearity on the beat signal of the measurement circuit can be eliminated to a large extent. The influence minimizes the deviation between the last estimated target flight time and the target actual flight time, and the target distance can be determined according to the target flight time, which improves the ranging accuracy of the lidar.
在一个可行的实施例中,校正单元具体用于:In a feasible embodiment, the correction unit is specifically used for:
S1:根据延时飞行时间τ' i-1对激光器发射信号的调频非线性项ε(t)进行延时,以得到测量路回波信号对应的调频非线性项ε(t-τ' i-1),延时飞行时间τ' i-1是基于飞行时间τ i-1得到的; S1: Delay the FM nonlinear term ε(t) of the laser emission signal according to the delay flight time τ′ i-1 to obtain the FM nonlinear term ε(t-τ′ i-) corresponding to the echo signal of the measurement path. 1 ), the delayed flight time τ'i -1 is obtained based on the flight time τ i-1 ;
S2:根据测量路回波信号对应的调频非线性项ε(t-τ' i-1)对测量路拍频信号进行调频非线性校正,以得到校正拍频信号; S2: Perform FM nonlinear correction on the beat frequency signal of the measurement path according to the FM nonlinear term ε(t-τ'i -1) corresponding to the echo signal of the measurement path to obtain a corrected beat frequency signal;
S3、根据校正拍频信号确定目标飞行时间。S3. Determine the target flight time according to the corrected beat frequency signal.
在一个可行的实施例中,在根据校正拍频信号确定目标飞行时间的方面,校正单元具体用于:In a feasible embodiment, in terms of determining the target flight time according to the corrected beat signal, the correction unit is specifically configured to:
根据校正拍频信号获取飞行时间τ i;当频谱峰值peak i满足预设条件时,将飞行时间τ i确定为目标飞行时间;当频谱峰值peak i不满足预设条件时,令i=i+1,并重复执行上述步骤S1-S3;频谱峰值peak i为校正拍频信号对应的频域信号的峰值; Obtain the flight time τ i according to the corrected beat signal; when the peak frequency peak i meets the preset condition, the flight time τ i is determined as the target flight time; when the peak frequency peak i does not meet the preset condition, let i=i+ 1. Repeat the above steps S1-S3; the peak value of the spectrum peak i is the peak value of the frequency domain signal corresponding to the corrected beat signal;
当i=1时,peak i-1为所述初始频谱峰值,飞行时间τ i-1为初始飞行时间τ。 When i=1, peak i-1 is the initial frequency spectrum peak, and the flight time τ i-1 is the initial flight time τ.
在一个可行的实施例中,预设条件为(peak i-peak i-1)/peak i-1<Thr rate,该Thr rate为第三预设阈值。 In one possible embodiment, the predetermined condition (peak i -peak i-1) / peak i-1 <Thr rate, Thr rate as the third predetermined threshold value.
在一个可行的实施例中,校正单元具体用于,包括:In a feasible embodiment, the correction unit is specifically configured to include:
根据初始飞行时间τ和激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行M次非线性迭代计算,以得到目标飞行时间,M为大于1的整数,According to the initial flight time τ and the frequency modulation non-linear term ε(t) of the laser emission signal, perform M non-linear iterative calculations on the measurement road beat frequency signal to obtain the target flight time, M is an integer greater than 1,
其中,在进行第i次非线性迭代计算时,根据延时飞行时间τ' i-1对激光器发射信号的调频非线性项ε(t)进行延时,以得到测量路回波信号对应的调频非线性项ε(t-τ' i-1);延时飞行时间τ' i-1是基于飞行时间τ i-1得到的; Among them, when the i-th nonlinear iterative calculation is performed, the frequency modulation nonlinear term ε(t) of the laser emission signal is delayed according to the delay flight time τ'i-1 to obtain the frequency modulation corresponding to the echo signal of the measurement path The non-linear term ε(t-τ'i -1 ); the delayed flight time τ'i -1 is obtained based on the flight time τi -1 ;
根据测量路回波信号对应的调频非线性项ε(t-τ' i-1)对测量路拍频信号进行调频非线性校正,以得到校正拍频信号;根据校正拍频信号获取飞行时间τ iPerform FM nonlinear correction on the measurement path beat signal according to the FM nonlinear term ε(t-τ' i-1 ) corresponding to the echo signal of the measurement path to obtain the corrected beat signal; obtain the flight time τ according to the corrected beat signal i ;
其中,当i=1时,飞行时间τ i-1为初始飞行时间τ;当i=M时,目标飞行时间为飞行时间τ iAmong them, when i=1, the flight time τ i-1 is the initial flight time τ; when i=M, the target flight time is the flight time τ i .
在一个可行的实施例中,在根据校正拍频信号获取飞行时间τ i的方面,计算单元具体用于: In a feasible embodiment, in terms of obtaining the flight time τ i according to the corrected beat signal, the calculation unit is specifically configured to:
对校正拍频信号进行FFT,以得到校正拍频信号对应的频域信号;根据校正拍频信号对应的频域信号获取飞行时间τ iPerform FFT on the corrected beat frequency signal to obtain the frequency domain signal corresponding to the corrected beat frequency signal; obtain the flight time τ i according to the frequency domain signal corresponding to the corrected beat frequency signal.
在一个可行的实施例中,校正拍频信号包括第一校正拍频信号或第二校正拍频信号;在根据测量路回波信号对应的调频非线性项ε(t-τ' i-1)对测量路拍频信号进行调频非线性校正,以得到校正拍频信号的方面,校正单元具体用于: In a feasible embodiment, the corrected beat frequency signal includes the first corrected beat frequency signal or the second corrected beat frequency signal; according to the frequency modulation nonlinearity term ε(t-τ' i-1 ) corresponding to the echo signal of the measurement path Perform FM nonlinear correction on the beat frequency signal of the measurement path to obtain the aspect of correcting the beat frequency signal. The correction unit is specifically used for:
当测量路拍频信号的信噪比不低于第二预设阈值时,根据激光器发射信号的调频非线性项ε(t)和测量路回波信号对应的调频非线性项ε(t-τ' i-1)对测量路拍频信号进行调频非线性项校正,以得到第一校正拍频信号;当测量路拍频信号的信噪比低于第二预设阈值时,根据测量路回波信号对应的调频非线性项ε(t-τ' i-1)对补偿后的拍频信号进行接收端调频非线性项校正,以得到第二校正拍频信号; When the signal-to-noise ratio of the beat frequency signal of the measurement path is not lower than the second preset threshold, the frequency modulation nonlinearity term ε(t) of the laser emission signal and the frequency modulation nonlinearity term ε(t-τ) corresponding to the echo signal of the measurement path ' I-1 ) Perform FM nonlinear term correction on the beat signal of the measurement road to obtain the first corrected beat signal; when the signal-to-noise ratio of the beat signal of the measurement road is lower than the second preset threshold, according to the measurement road The FM nonlinear term ε(t-τ′i -1 ) corresponding to the wave signal performs receiving end FM nonlinear term correction on the compensated beat signal to obtain a second corrected beat signal;
其中,补偿后的拍频信号是基于根据激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行发射端的调频非线性项补偿得到的。Among them, the compensated beat signal is obtained based on the FM nonlinear term ε(t) of the laser emitting signal to compensate the FM nonlinear term of the transmitter of the measurement path beat signal.
在一个可行的实施例中,在根据测量路拍频信号计算得到初始飞行时间τ的方面,计算单元具体用于:In a feasible embodiment, in terms of calculating the initial flight time τ according to the measured road beat signal, the calculation unit is specifically configured to:
当测量路拍频信号的信噪比不低于第二预设阈值时,对测量路拍频信号进行FFT,以得到测量路拍频信号对应的频域信号;根据测量路拍频信号对应的频域信号确定测量路拍频信号的频率,并根据测量路拍频信号的频率计算得到初始飞行时间τ;When the signal-to-noise ratio of the measured road beat signal is not lower than the second preset threshold, perform FFT on the measured road beat signal to obtain the frequency domain signal corresponding to the measured road beat signal; The frequency domain signal determines the frequency of the measured road beat signal, and calculates the initial flight time τ according to the frequency of the measured road beat signal;
其中,测量路拍频信号的频率为测量路拍频信号对应的频域信号的峰值点位置对应的频率,或者频域信号的半高宽度中点位置对应的频率。Wherein, the frequency of the measurement road beat signal is the frequency corresponding to the peak point position of the frequency domain signal corresponding to the measurement road beat signal, or the frequency corresponding to the midpoint position of the half-height width of the frequency domain signal.
在一个可行的实施例中,在根据测量路拍频信号计算得到初始飞行时间τ的方面,计算单元具体用于:In a feasible embodiment, in terms of calculating the initial flight time τ according to the measured road beat signal, the calculation unit is specifically configured to:
当测量路拍频信号的信噪比低于第二预设阈值时;根据激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行发射端的调频非线性项补偿,以得到补偿后的拍频信号;激光器发射信号的调频非线性项是基于参考路拍频信号计算得到的;When the signal-to-noise ratio of the measurement road beat signal is lower than the second preset threshold; according to the FM nonlinear term ε(t) of the laser emitting signal, the measurement road beat signal is compensated by the FM nonlinear term at the transmitting end. After the beat frequency signal; the frequency modulation nonlinear term of the laser emission signal is calculated based on the beat frequency signal of the reference path;
对补偿后的拍频信号进行FFT,以得到补偿后的拍频信号对应的频域信号;根据补偿后的拍频信号对应的频域信号确定补偿后的拍频信号的频率,并根据补偿后的拍频信号的频率计算得到初始飞行时间τ;Perform FFT on the compensated beat signal to obtain the frequency domain signal corresponding to the compensated beat signal; determine the frequency of the compensated beat signal according to the frequency domain signal corresponding to the compensated beat signal, and according to the compensated beat signal Calculate the frequency of the beat signal to get the initial flight time τ;
其中,补偿后的拍频信号的频率为补偿后的拍频信号对应的频域信号的峰值点位置对应的频率,或者频域信号的半高宽度中点位置对应的频率。The frequency of the compensated beat signal is the frequency corresponding to the peak point position of the frequency domain signal corresponding to the compensated beat signal, or the frequency corresponding to the midpoint position of the half-height width of the frequency domain signal.
在一个可行的实施例中,在根据参考路拍频信号计算得到激光器发射信号的调频非线性项ε(t)的方面,计算单元具体用于:In a feasible embodiment, in terms of calculating the frequency modulation nonlinear term ε(t) of the laser emission signal according to the reference path beat signal, the calculation unit is specifically configured to:
对参考路拍频信号进行希尔伯特变换,以得到变换后的拍频信号;根据该变换后的拍频信号计算得到激光器发射信号的调频频率;根据激光器发射信号的调频频率和理想线性调频频率计算得到激光器发射信号的调频非线性项ε(t)。Perform Hilbert transform on the beat signal of the reference path to obtain the transformed beat signal; calculate the FM frequency of the laser emission signal according to the transformed beat signal; According to the FM frequency and ideal chirp of the laser emission signal Frequency calculation obtains the frequency modulation nonlinear term ε(t) of the laser emission signal.
第三方面,本申请实施例提供一种基于调频非线性校正的测距装置,包括:In a third aspect, an embodiment of the present application provides a ranging device based on frequency modulation nonlinear correction, including:
存储器,用于存储指令;以及Memory for storing instructions; and
至少一台处理器,与所述存储器耦合;At least one processor coupled with the memory;
其中,当所述至少一台处理器执行所述指令时,所述指令致使所述处理器执行如第一方面所示的全部或者部分方法。Wherein, when the at least one processor executes the instruction, the instruction causes the processor to execute all or part of the method shown in the first aspect.
第四方面,本申请实施例提供一种计算机存储介质,该计算机存储介质存储有计算机程序,所述计算机程序包括程序指令,所述程序指令当被处理器执行时使所述处理器执行如第一方面所示的全部或者部分方法。In a fourth aspect, an embodiment of the present application provides a computer storage medium that stores a computer program. The computer program includes program instructions that, when executed by a processor, cause the processor to execute All or part of the method shown on the one hand.
第五方面,本申请实施例提供一种激光雷达,包括激光器以及与该激光器耦合的处理器,该处理器执行如第一方面所示的全部或部分方法。In a fifth aspect, an embodiment of the present application provides a laser radar, including a laser and a processor coupled with the laser, and the processor executes all or part of the method shown in the first aspect.
附图说明Description of the drawings
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍。In order to more clearly describe the technical solutions in the embodiments of the present application or the prior art, the following will briefly introduce the drawings that need to be used in the description of the embodiments or the prior art.
图1为FMCW激光雷达测距原理示意图;Figure 1 is a schematic diagram of the FMCW lidar ranging principle;
图2为FMCW激光雷达的系统架构框图;Figure 2 is a block diagram of the system architecture of FMCW lidar;
图3为非线性调频下的拍频信号;Figure 3 shows the beat signal under non-linear frequency modulation;
图4为现有技术中基于调制激光器驱动电流预失真的处理方法示意图;4 is a schematic diagram of a processing method based on modulated laser drive current predistortion in the prior art;
图5为本申请实施例提供的一种基于调频非线性校正的测距方法的流程示意图;FIG. 5 is a schematic flowchart of a ranging method based on frequency modulation nonlinear correction according to an embodiment of the application;
图6为本申请实施例提供的一种对测量路拍频信号进行非线性迭代校正的流程示意图;FIG. 6 is a schematic flowchart of a non-linear iterative correction of a measurement road beat frequency signal provided by an embodiment of the application;
图7为本申请实施例提供的另一种对测量路拍频信号进行非线性迭代校正的流程示意图;FIG. 7 is a schematic diagram of another flow chart for performing nonlinear iterative correction on measurement road beat frequency signals according to an embodiment of the application;
图8为本申请实施例提供的另一种基于调频非线性校正的测距方法的流程示意图;FIG. 8 is a schematic flowchart of another ranging method based on frequency modulation nonlinear correction provided by an embodiment of the application;
图9a为非线性校正前后测距误差对比图;Figure 9a is a comparison diagram of ranging errors before and after nonlinear correction;
图9b为非线性项迭代校正后拍频信号的频谱;Figure 9b is the frequency spectrum of the beat signal after iterative correction of the nonlinear term;
图10为本申请实施例提供的一种基于调频非线性校正的测距装置的结构示意图;10 is a schematic structural diagram of a ranging device based on frequency modulation nonlinear correction provided by an embodiment of the application;
图11为本申请实施例提供的另一种基于调频非线性校正的测距装置的结构示意图。FIG. 11 is a schematic structural diagram of another ranging device based on frequency modulation nonlinear correction provided by an embodiment of the application.
具体实施方式Detailed ways
为了使本技术领域的人员更好地理解本申请方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述。In order to enable those skilled in the art to better understand the solutions of the present application, the technical solutions in the embodiments of the present application will be described clearly and completely in conjunction with the accompanying drawings in the embodiments of the present application.
为了提高测距精度和增大测距范围,人们提出了校正调频非线性的方法。在现有的校 正调频非线性的方法中,较为普遍的一种是基于调制激光器驱动电流预失真的处理方法。如图4所示,通过估计当前测量激光器实际输出频率v k(t)与理想线性频率v d(t)之间的频率差e k(t),来调整下一次测量激光器驱动电流u k+1(t),从而改变下一次测量时激光器的输出频率。经过多轮这样的迭代校正,最终使得激光输出频率趋于线性。 In order to improve the accuracy of the ranging and increase the range of the ranging, a method of correcting the non-linearity of the frequency modulation has been proposed. Among the existing methods for correcting frequency modulation nonlinearity, a more common one is a processing method based on the predistortion of the modulated laser drive current. As shown in Figure 4, by estimating the frequency difference e k (t) between the actual output frequency v k (t) of the current measurement laser and the ideal linear frequency v d (t), the next measurement laser driving current u k+ is adjusted 1 (t) to change the output frequency of the laser in the next measurement. After multiple rounds of such iterative corrections, the laser output frequency tends to be linear in the end.
这个方案的主要问题在于,通过调整驱动电流的方式来控制激光器的输出频率,只能将非线性减小到一个比较小的范围内,并不能从原理上完全消除。此外,当前非线性校正的结果只能用于下一次测量。倘若由于温度、振动等外界因素的干扰,通过上一次测量时预失真估计出的驱动电流来控制下一次测量时激光器的输出频率,其频率值与预期值不匹配,这就导致上一次预失真估计结果失效,而没有达迭代校正的目的。The main problem of this scheme is that by adjusting the drive current to control the output frequency of the laser, the nonlinearity can only be reduced to a relatively small range, but it cannot be completely eliminated in principle. In addition, the results of the current nonlinear correction can only be used for the next measurement. If due to the interference of external factors such as temperature and vibration, the output frequency of the laser in the next measurement is controlled by the drive current estimated by the predistortion in the last measurement, and the frequency value does not match the expected value, which leads to the previous predistortion The estimation result is invalid, and it fails to achieve the purpose of iterative correction.
重采样算法是一种在接收端进行调频非线性校正的方法。核心思想是,通过估计的采样时间序列对拍频信号进行重采样,使得信号重新变为单频信号。具体原理如下:The resampling algorithm is a method of FM nonlinear correction at the receiving end. The core idea is to resample the beat frequency signal through the estimated sampling time sequence, so that the signal becomes a single frequency signal again. The specific principle is as follows:
(1)发射信号:
Figure PCTCN2020074495-appb-000001
(1) Transmission signal:
Figure PCTCN2020074495-appb-000001
(2)回波信号:s r(t)=s t(t-τ) (2) Echo signal: s r (t) = s t (t-τ)
(3)拍频信号:
Figure PCTCN2020074495-appb-000002
(忽略τ的高阶项)
(3) Beat signal:
Figure PCTCN2020074495-appb-000002
(Ignore the higher-order terms of τ)
(4)令
Figure PCTCN2020074495-appb-000003
则s′ b(t)=exp{j2π[f 0+αt′]τ}
(4) Order
Figure PCTCN2020074495-appb-000003
Then s′ b (t)=exp{j2π[f 0 +αt′]τ}
将t′看作新的采样时间序列,对拍频信号s′ b(t)进行重采样,使得重采样后的拍频信号变回为单频信号。
Figure PCTCN2020074495-appb-000004
是调频非线性项。
The t 'treated as new sample time series of the beat signal s' b (t) for resampling, so that the beat frequency signal is resampled back single frequency signal.
Figure PCTCN2020074495-appb-000004
Is the FM nonlinear term.
重采样算法的前提假设是拍频信号的延时τ为一个较小的值,则可以忽略τ的高阶项,对于较远的目标,其拍频信号对应的延时τ由于不满足前提假设,导致重采样处理后其对应的测距误差仍然很大。所以,重采样的性能与待测距离有关,具有很大的实用局限性。The premise of the re-sampling algorithm is that the delay τ of the beat signal is a small value, and the higher-order term of τ can be ignored. For a far target, the delay τ corresponding to the beat signal does not meet the premise assumptions. , Resulting in the corresponding ranging error after resampling is still very large. Therefore, the performance of resampling is related to the distance to be measured, which has great practical limitations.
基于上述原因,本申请提出一种基于调频非线性校正的测距方法的流程示意图。如图5所示,该方法包括:Based on the foregoing reasons, this application proposes a schematic flow chart of a ranging method based on frequency modulation nonlinear correction. As shown in Figure 5, the method includes:
S501、获取参考路拍频信号和测量路拍频信号。S501: Obtain a reference road beat signal and a measurement road beat signal.
具体地,激光器发射的激光信号分成两路,一路为参考路的光信号,另一路为测量路的光信号。其中,参考路光信号和测量路光信号的频率和相位相同。Specifically, the laser signal emitted by the laser is divided into two paths, one path is the optical signal of the reference path, and the other path is the optical signal of the measurement path. Among them, the frequency and phase of the optical signal of the reference path and the optical signal of the measurement path are the same.
参考路的光信号分为两路,可通过预设固定臂长差为τ D的马赫-曾德尔干涉仪(Mach-Zehnder interferometer,MZI)得到参考路拍频信号s ref(t)。具体地,一路光信号经过MZI的短臂,另一路光信号经过MZI的长臂,因为两路光信号经过光程不同,所以两路光信号在平衡探测器的光敏面上进行相干,得到参考路拍频信号s ref(t)。 The optical signal of the reference path is divided into two paths, and the reference path beat signal s ref (t) can be obtained by a Mach-Zehnder interferometer (MZI) with a preset fixed arm length difference of τ D. Specifically, one optical signal passes through the short arm of the MZI, and the other optical signal passes through the long arm of the MZI. Because the two optical signals pass through different optical paths, the two optical signals are coherent on the photosensitive surface of the balanced detector to obtain a reference Road beat signal s ref (t).
测量路的光信号分为两路,其中,一路为本振信号,另一路发射至被测目标物,经该目标物表面反射得到回波信号。回波信号与本振信号在平衡探测器的光敏面进行相干,以得到测量路拍频信号s b(t)。 The optical signal of the measurement path is divided into two paths, one of which is a local oscillator signal, and the other is emitted to the measured target, and the echo signal is obtained by reflection on the surface of the target. The echo signal and the local oscillator signal are coherent on the photosensitive surface of the balanced detector to obtain the measurement road beat signal s b (t).
S502、根据参考路拍频信号计算得到调频非线性度;并根据测量路拍频信号计算得到初始飞行时间τ。S502: Calculate the frequency modulation nonlinearity according to the reference road beat frequency signal; and calculate the initial flight time τ according to the measured road beat frequency signal.
具体地,根据参考路拍频信号计算得到调频非线性度,包括:Specifically, calculating the frequency modulation nonlinearity according to the reference road beat signal includes:
对参考路拍频信号s ref(t)进行希尔伯特变换,以得到变换后的信号H b(t),根据变换后的信号得到参考路拍频信号的瞬时相位
Figure PCTCN2020074495-appb-000005
Perform Hilbert transform on the reference road beat signal s ref (t) to obtain the transformed signal H b (t), and obtain the instantaneous phase of the reference road beat signal according to the transformed signal
Figure PCTCN2020074495-appb-000005
其中,
Figure PCTCN2020074495-appb-000006
H b,real=real(H b(t)),H b,imag=imag(H b(t)),H b(t)=hilbert(s ref(t)),real(H b(t))为变换后的信号H b(t)的实部,imag(H b(t))为变换后的信号H b(t)的虚部。
in,
Figure PCTCN2020074495-appb-000006
H b, real = real(H b (t)), H b, imag = imag(H b (t)), H b (t)=hilbert(s ref (t)), real(H b (t) ) Is the real part of the transformed signal H b (t), and imag(H b (t)) is the imaginary part of the transformed signal H b (t).
根据激光器发射信号的调频频率f(t)与参考路拍频信号瞬时相位
Figure PCTCN2020074495-appb-000007
的对应关系,计算激光器发射信号的调频频率f(t),其中,
Figure PCTCN2020074495-appb-000008
According to the instantaneous phase of the frequency modulation frequency f(t) of the laser emission signal and the beat frequency signal of the reference path
Figure PCTCN2020074495-appb-000007
Calculate the frequency modulation frequency f(t) of the laser emission signal, where,
Figure PCTCN2020074495-appb-000008
根据调频频率计算调频非线性度,该调频非线性度
Figure PCTCN2020074495-appb-000009
f ideal(t)为理想线性调频频率,f ideal(t)=f 0t+0.5αt 2,f 0为载波频率,α是理想线性调频时对应的调频斜率,B为调频带宽。
Calculate the non-linearity of the frequency modulation according to the frequency of the frequency modulation, the non-linearity of the frequency modulation
Figure PCTCN2020074495-appb-000009
f ideal (t) is the ideal chirp frequency, f ideal (t)=f 0 t+0.5αt 2 , f 0 is the carrier frequency, α is the corresponding FM slope when ideal chirp, and B is the FM bandwidth.
可选地,调频非线性度还可以是根据激光器发射信号的调频非线性项(即激光器发射信号的调频频率f(t)与理想线性调频频率f ideal(t)的差)的均方根和调频带宽(B)计算得到,调频非线性度可表示为: Optionally, the frequency modulation nonlinearity can also be based on the root mean square sum of the frequency modulation nonlinearity term of the laser emission signal (that is, the difference between the frequency modulation frequency f(t) of the laser emission signal and the ideal chirp frequency f ideal (t)). The frequency modulation bandwidth (B) is calculated, and the frequency modulation nonlinearity can be expressed as:
Figure PCTCN2020074495-appb-000010
Figure PCTCN2020074495-appb-000010
Figure PCTCN2020074495-appb-000011
Figure PCTCN2020074495-appb-000011
可选地,可根据激光器发射信号的调频频率f(t)与理想线性调频频率f ideal(t)的偏离程度来衡量激光器的调频非线性程度,即用f(t)与f ideal(t)之间的均方根误差(Root Mean Squard Error,RMSE)来衡量,即 Optionally, the degree of frequency modulation nonlinearity of the laser can be measured according to the degree of deviation between the frequency modulation frequency f(t) of the laser emission signal and the ideal chirp frequency f ideal (t), that is, f(t) and f ideal (t) Root Mean Squard Error (RMSE) between
Figure PCTCN2020074495-appb-000012
Figure PCTCN2020074495-appb-000012
在此需要说明的是,调频频率f(t)和理想线性调频频率f ideal(t)均可看成离散信号,因此f k为调频频率f(t)上第k个离散点对应的频率,f ideal,k为理想线性调频频率f ideal(t)上第k个离散点对应的频率。 It should be noted here that both the FM frequency f(t) and the ideal chirp frequency f ideal (t) can be regarded as discrete signals, so f k is the frequency corresponding to the k-th discrete point on the FM frequency f(t), f ideal, k is the frequency corresponding to the k-th discrete point on the ideal chirp frequency f ideal (t).
S503、当调频非线性度不小于第一预设阈值时,根据参考路拍频信号计算得到激光器发射信号的调频非线性项ε(t),根据初始飞行时间τ和激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行非线性迭代校正,以得到目标飞行时间。S503. When the frequency modulation nonlinearity is not less than the first preset threshold, the frequency modulation nonlinearity term ε(t) of the laser emission signal is calculated according to the reference path beat frequency signal, and the frequency modulation nonlinearity term ε(t) of the laser emission signal is calculated according to the initial flight time τ and the frequency modulation nonlinearity of the laser emission signal The term ε(t) performs nonlinear iterative correction on the measured road beat signal to obtain the target flight time.
其中,第一预设阈值Thr nonRate为FMCW测距系统所能接收的激光器的最大非线性度。 Wherein, the first preset threshold Thr nonRate is the maximum non-linearity of the laser that can be received by the FMCW ranging system.
具体地,激光器发射信号的调频非线性项ε(t)基于以下公式得到:Specifically, the frequency modulation nonlinear term ε(t) of the laser emission signal is obtained based on the following formula:
Figure PCTCN2020074495-appb-000013
Figure PCTCN2020074495-appb-000013
Figure PCTCN2020074495-appb-000014
Figure PCTCN2020074495-appb-000014
在一个可行的实施例中,根据测量路拍频信号s b(t)计算得到初始飞行时间τ,包括: In a feasible embodiment, calculating the initial flight time τ according to the measured road beat signal s b (t) includes:
当测量路拍频信号s b(t)的信噪比不小于第二预设阈值时,对测量路拍频信号s b(t)进行快速FFT,得到该测量路拍频信号s b(t)对应的频域信号,基于该频域信号确定测量路拍频信号s b(t)的频率f b;当测量路拍频信号s b(t)的信噪比小于第二预设阈值时,根据激光器发射信号的调频非线性项ε(t)对测量路拍频信号s b(t)进行发射端的调频非线性项补偿,以得到补偿后的拍频信号s' b(t);对补偿后的拍频信号s' b(t)进行FFT,以得到补偿后的拍频信号对应的频域信号;根据补偿后的拍频信号对应的频域信号确定补偿后的拍频信号s' b(t)的频率f' b;根据拍频信号的频率f b或补偿后的拍频信号s' b(t)的频率f' b计算得到初始飞行时间τ;初始飞行时间τ为f b/α或f' b/α,α为理想线性调频时对应的调频斜率; When the signal-to-noise ratio of the measured road beat signal s b (t) is not less than the second preset threshold, fast FFT is performed on the measured road beat signal s b (t) to obtain the measured road beat signal s b (t ) corresponding to frequency domain signal based on the frequency-domain channel measurement signal to determine the beat signal s b (t) of frequency f b; when measured channel beat signal s b (t) is less than a second predetermined threshold SNR , According to the frequency modulation non-linear term ε(t) of the laser emission signal, perform the frequency modulation non-linear term compensation at the transmitting end of the beat frequency signal s b (t) of the measuring path to obtain the compensated beat frequency signal s b (t); beat signal s compensated 'b (t) for FFT, to obtain a frequency-domain signal of the beat signal corresponding to the compensated; according shot compensated pilot signal corresponding to a frequency domain signal to determine the beat signal s compensated' b (t) of frequency f 'b; the beat signal s of frequency f b or compensation of the beat signal' b (t) of frequency f 'b calculated from the initial flight time [tau]; an initial flight time [tau] is f b when the corresponding ideal chirp / α or f 'b / α, α chirp rate;
其中,测量路拍频信号s b(t)的频率f b为测量路拍频信号s b(t)对应的频域信号的峰值点位置对应的频率,或者该频域信号半高宽度中心位置对应的频率;补偿后的拍频信号s' b(t)的频率f' b为补偿后的拍频信号s' b(t)对应的频域信号的峰值点位置对应的频率,或者该频域信号半高宽度中心位置对应的频率。 Wherein the measuring channel beat signal s b (t) of frequency f b for the measurement channel beat signal s b (t) corresponding to the peak position frequency corresponding to the frequency domain signal, or a signal half-width of the center position of the frequency domain corresponding to a frequency; beat signal s compensated 'b (t) of frequency f' b is a beat signal s compensated 'peak point b (t) corresponding to the frequency domain signal corresponding to the position of the frequency or the frequency The frequency corresponding to the center of the half-height width of the domain signal.
在一个可行的实施例中,根据初始飞行时间τ和激光器发射信号的调频非线性项ε(t)进行非线性迭代校正,以得到目标飞行时间,包括:In a feasible embodiment, the non-linear iterative correction is performed according to the initial flight time τ and the frequency modulation non-linear term ε(t) of the laser emission signal to obtain the target flight time, including:
S1:根据延时飞行时间τ' i-1对激光器发射信号的调频非线性项ε(t)进行延时,以得到测量路回波信号的调频非线性项ε(t-τ' i-1),其中,所述延时飞行时间τ' i-1是基于飞行时间τ i-1得到的; S1: Delay the FM nonlinear term ε(t) of the laser emission signal according to the delay flight time τ′ i-1 to obtain the FM nonlinear term ε(t-τ′ i-1) of the echo signal of the measurement path ), wherein the delayed flight time τ'i -1 is obtained based on the flight time τ i-1 ;
S2、根据测量路回波信号对应的调频非线性项对测量路拍频信号进行调频非线性项校正,以得到校正拍频信号;S2, according to the FM nonlinear term corresponding to the echo signal of the measurement path, perform FM nonlinear term correction on the beat frequency signal of the measurement path to obtain a corrected beat frequency signal;
S3、根据校正拍频信号确定目标飞行时间。S3. Determine the target flight time according to the corrected beat frequency signal.
进一步地,根据校正拍频信号确定目标飞行时间,包括:Further, determining the target flight time according to the corrected beat frequency signal includes:
根据校正拍频信号获取飞行时间τ i;当频谱峰值peak i满足预设条件时,将飞行时间τ i确定为目标飞行时间;当频谱峰值peak i不满足预设条件时,令i=i+1,并重复执行步骤S1-S3;频谱峰值peak i为校正拍频信号对应的频域信号的峰值; Obtain the flight time τ i according to the corrected beat signal; when the peak frequency peak i meets the preset condition, the flight time τ i is determined as the target flight time; when the peak frequency peak i does not meet the preset condition, let i=i+ 1. Repeat steps S1-S3; the peak value of the spectrum peak i is the peak value of the frequency domain signal corresponding to the corrected beat signal;
当i=1时,peak i-1为初始频谱峰值,飞行时间τ i-1为初始飞行时间τ。 When i=1, peak i-1 is the initial spectral peak, and the flight time τ i-1 is the initial flight time τ.
其中,预设条件为(peak i-peak i-1)/peak i-1<Thr rate,Thr rate为第三预设阈值。 The preset condition is (peak i -peak i-1 )/peak i-1 <Thr rate , and Thr rate is the third preset threshold.
在一个可行的实施例中,根据初始飞行时间τ和所述激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行非线性迭代校正,以得到目标飞行时间,包括:In a feasible embodiment, performing nonlinear iterative correction on the measurement path beat frequency signal according to the initial flight time τ and the frequency modulation nonlinear term ε(t) of the laser emission signal to obtain the target flight time includes:
根据初始飞行时间τ和激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行M次非线性迭代计算,以得到目标飞行时间,M为大于1的整数,According to the initial flight time τ and the frequency modulation non-linear term ε(t) of the laser emission signal, perform M non-linear iterative calculations on the measurement path beat frequency signal to obtain the target flight time, M is an integer greater than 1,
其中,在进行第i次非线性迭代计算时,根据所述延时飞行时间τ' i-1对所述激光器发射信号的调频非线性项ε(t)进行延时,以得到测量路回波信号对应的调频非线性项ε(t-τ' i-1);延时飞行时间τ' i-1是基于飞行时间τ i-1得到的; Wherein, when the i-th nonlinear iterative calculation is performed, the frequency modulation nonlinear term ε(t) of the laser emission signal is delayed according to the delay flight time τ′ i-1 to obtain the measurement path echo The non-linear term ε(t-τ'i -1 ) corresponding to the frequency modulation of the signal; the delayed flight time τ'i -1 is obtained based on the flight time τi -1 ;
根据测量路回波信号对应的调频非线性项ε(t-τ' i-1)对测量路拍频信号进行调频非线性校正,以得到校正拍频信号;根据校正拍频信号获取飞行时间τ iPerform FM nonlinear correction on the measurement path beat signal according to the FM nonlinear term ε(t-τ' i-1 ) corresponding to the echo signal of the measurement path to obtain the corrected beat signal; obtain the flight time τ according to the corrected beat signal i ;
其中,当i=1时,飞行时间τ i-1为初始飞行时间τ;当i=M时,目标飞行时间为飞行时间τ iAmong them, when i=1, the flight time τ i-1 is the initial flight time τ; when i=M, the target flight time is the flight time τ i .
在一个可行的实施例中,根据校正拍频信号获取飞行时间τ i,包括: In a feasible embodiment, obtaining the flight time τ i according to the corrected beat signal includes:
对校正拍频信号进行FFT以得到校正拍频信号对应的频域信号,根据校正拍频信号对应的频域信号获取飞行时间τ iPerform FFT on the corrected beat frequency signal to obtain a frequency domain signal corresponding to the corrected beat frequency signal, and obtain the flight time τ i according to the frequency domain signal corresponding to the corrected beat frequency signal.
在一个可行的实施例中,校正拍频信号包括第一校正拍频信号或第二校正拍频信号;根据测量路回波信号对应的调频非线性项ε(t-τ' i-1)对测量路拍频信号进行调频非线性校正,以得到校正拍频信号,包括: In one possible embodiment, the correction signal comprises a beat frequency signal a first correction or the second correction beat beat signal; an echo signal corresponding to the measured channel FM nonlinear term ε (t-τ 'i- 1) of Measure the beat signal of the road and perform FM nonlinear correction to obtain the corrected beat signal, including:
当测量路拍频信号的信噪比不低于第二预设阈值时,根据激光器发射信号的调频非线性项ε(t)和测量路回波信号对应的调频非线性项ε(t-τ' i-1)对测量路拍频信号进行调频非线性项校正,以得到第一校正拍频信号; When the signal-to-noise ratio of the beat frequency signal of the measurement path is not lower than the second preset threshold, the frequency modulation nonlinearity term ε(t) of the laser emission signal and the frequency modulation nonlinearity term ε(t-τ) corresponding to the echo signal of the measurement path ' I-1 ) Perform frequency modulation nonlinear term correction on the beat frequency signal of the measurement path to obtain the first corrected beat frequency signal;
当测量路拍频信号的信噪比低于所述第二预设阈值时,根据测量路回波信号对应的调频非线性项ε(t-τ' i-1)对补偿后的拍频信号进行接收端调频非线性项校正,以得到第二校正拍频信号; When the signal-to-noise ratio of the beat frequency signal of the measurement path is lower than the second preset threshold, the compensated beat frequency signal ε(t-τ'i-1) corresponding to the frequency modulation nonlinearity of the echo signal of the measurement path Performing correction of the non-linear term of the FM at the receiving end to obtain the second corrected beat signal;
其中,补偿后的拍频信号是基于根据激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行发射端的调频非线性项补偿得到的,因此补偿后的拍频信号又可称为已做发射端调频非线性补偿的拍频信号。Among them, the compensated beat signal is based on the FM nonlinear term ε(t) of the laser emission signal to compensate for the FM nonlinear term at the transmitting end of the beat signal of the measuring path, so the compensated beat signal can also be called It is a beat signal that has been compensated for FM nonlinearity at the transmitter.
具体地,当测量路拍频信号s b(t)的信噪比不低于第二预设阈值时,根据激光器发射信号的调频非线性项ε(t)和回波信号的调频非线性项ε(t-τ i)对测量路拍频信号s b(t)进行调频非线性项校正,以得到第一校正拍频信号
Figure PCTCN2020074495-appb-000015
当测量路拍频信号s b(t)的信噪比低于第二预设阈值时,根据测量路回波信号的调频非线性项ε(t-τ i)对已做发射端调频非线性补偿后的拍频信号s' b(t)进行接收端调频非线性项校正,以得到第二校正拍频信号
Figure PCTCN2020074495-appb-000016
Specifically, when the signal-to-noise ratio of the measured road beat signal s b (t) is not lower than the second preset threshold, according to the frequency modulation nonlinear term ε(t) of the laser emission signal and the frequency modulation nonlinear term of the echo signal ε(t-τ i ) Perform frequency modulation nonlinear term correction on the measurement path beat signal s b (t) to obtain the first corrected beat signal
Figure PCTCN2020074495-appb-000015
When the signal-to-noise ratio of the measured road beat signal s b (t) is lower than the second preset threshold, the FM nonlinearity term ε(t-τ i ) of the measured road echo signal is performed on the transmitter. beat signal s' b (t) FM reception terminal nonlinear term compensation corrected to obtain a second corrected beat signal
Figure PCTCN2020074495-appb-000016
具体地,如图6所示:Specifically, as shown in Figure 6:
S10A、根据飞行时间τ i-1获取延时飞行时间τ' i-1S10A. Obtain the delayed flight time τ'i-1 according to the flight time τ i-1 .
可选地,延时飞行时间τ' i-1=τ i-1+p,p为预设步长。当然还可以基于其他方式得到延时飞行时间τ' i-1Optionally, the delayed flight time τ′ i-1i-1 +p, where p is the preset step length. Of course, the delayed flight time τ′ i-1 can also be obtained based on other methods.
S20A、根据延时飞行时间τ' i-1对激光器发射信号的调频非线性项ε(t)进行延时,以得到测量路回波信号的调频非线性项ε(t-τ' i-1)。 S20A. Delay the FM nonlinear term ε(t) of the laser emission signal according to the delay flight time τ′ i-1 to obtain the FM nonlinear term ε(t-τ′ i-1) of the echo signal of the measurement path ).
换句话说,对激光器发射信号的调频非线性项ε(t)进行一个τ' i-1的延时,得到测量路回波信号的调频非线性项ε(t-τ' i-1)。 In other words, the FM nonlinear term ε(t) of the laser emission signal is delayed by τ′ i-1 to obtain the FM nonlinear term ε(t-τ′ i-1 ) of the echo signal of the measurement path.
S30A、当测量路拍频信号s b(t)的信噪比不低于第二预设阈值时,根据激光器发射信号的调频非线性项ε(t)和测量路回波信号的调频非线性项ε(t-τ' i-1)对测量路拍频信号s b(t)进行调频非线性项校正,以得到第一校正拍频信号
Figure PCTCN2020074495-appb-000017
其中,第一校正拍频信号为:
Figure PCTCN2020074495-appb-000018
s b(t)为测量路的初始拍频信号,换句话说,s b(t)为测量路非线性项迭代校正前对应的拍频信号。当测量路拍频信号s b(t)的信噪比低于第二预设阈值时,根据测量路回波信号的调频非线性项ε(t-τ' i-1)对已做发端调频非线性项补偿后的拍频信号s' b(t)进行接收端调频非线性项校正,以得到第二校正拍频信号
Figure PCTCN2020074495-appb-000019
其中,第二校正拍频信号为:
Figure PCTCN2020074495-appb-000020
对第一校正拍频信号
Figure PCTCN2020074495-appb-000021
或第二校正拍频信号
Figure PCTCN2020074495-appb-000022
进行FFT,以得到该第一校正拍频信号
Figure PCTCN2020074495-appb-000023
或第二校正拍频信号
Figure PCTCN2020074495-appb-000024
对应的频域信号,基于该频域信号获取飞行时间τ i和频谱峰值peak i
S30A. When the signal-to-noise ratio of the measured road beat signal s b (t) is not lower than the second preset threshold, according to the frequency modulation non-linear term ε(t) of the laser emission signal and the frequency modulation non-linearity of the echo signal of the measurement road The term ε(t-τ'i -1 ) performs FM nonlinear term correction on the measurement path beat signal s b (t) to obtain the first corrected beat signal
Figure PCTCN2020074495-appb-000017
Among them, the first corrected beat signal is:
Figure PCTCN2020074495-appb-000018
s b (t) is the initial beat signal of the measurement path, in other words, s b (t) is the corresponding beat signal before the iterative correction of the nonlinear term of the measurement path. When the signal-to-noise ratio of the measured road beat signal s b (t) is lower than the second preset threshold, the transmitted frequency modulation is performed according to the frequency modulation non-linear term ε(t-τ' i-1 ) of the measured road echo signal the beat signal s nonlinearity compensation 'b (t) FM reception terminal nonlinear term correction, the second correction to obtain a beat signal
Figure PCTCN2020074495-appb-000019
Among them, the second corrected beat signal is:
Figure PCTCN2020074495-appb-000020
The first corrected beat signal
Figure PCTCN2020074495-appb-000021
Or the second corrected beat signal
Figure PCTCN2020074495-appb-000022
Perform FFT to obtain the first corrected beat signal
Figure PCTCN2020074495-appb-000023
Or the second corrected beat signal
Figure PCTCN2020074495-appb-000024
The corresponding frequency domain signal is used to obtain the flight time τ i and the spectral peak value peak i based on the frequency domain signal.
具体地,根据第一校正拍频信号对应的频域信号确定第一校正拍频信号的频率,并根据第一校正拍频信号的频率确定飞行时间τ i;或者根据第二校正拍频信号对应的频域信号确定第二校正拍频信号的频率,并根据第二校正拍频信号的频率确定飞行时间τ iSpecifically, the frequency of the first corrected beat signal is determined according to the frequency domain signal corresponding to the first corrected beat signal, and the flight time τ i is determined according to the frequency of the first corrected beat signal; or corresponding to the second corrected beat signal Determine the frequency of the second corrected beat signal from the frequency domain signal, and determine the flight time τ i according to the frequency of the second corrected beat signal.
其中,第一校正拍频信号的频率为第一校正拍频信号对应的频域信号的峰值peak i位置对应的频率,或者该频域信号半高宽度中心位置对应的频率;第二校正拍频信号的频率为第二校正拍频信号对应的频域信号的峰值peak i位置对应的频率,或者该频域信号半高宽度中心位置对应的频率。 The frequency of the first corrected beat signal is the frequency corresponding to the peak i position of the frequency domain signal corresponding to the first corrected beat signal, or the frequency corresponding to the center position of the half-height width of the frequency domain signal; the second corrected beat frequency The frequency of the signal is the frequency corresponding to the peak i position of the frequency domain signal corresponding to the second corrected beat frequency signal, or the frequency corresponding to the center position of the half-height width of the frequency domain signal.
S40A、当(peak i-peak i-1)/peak i-1>=Thr rate时,令i=1+1,并重复执行步骤S10A-S40A;当(peak i-peak i-1)/peak i-1<Thr rate时,执行步骤S50A。其中,Thr rate为第三预设阈值,该第三预设阈值为峰值增长速度门限值。 S40A. When (peak i -peak i-1 )/peak i-1 >= Thr rate , set i=1+1, and repeat steps S10A-S40A; when (peak i -peak i-1 )/peak When i-1 &lt; Thr rate , step S50A is executed. Wherein, Thr rate is a third preset threshold, and the third preset threshold is a peak growth rate threshold.
其中,当i=1时,peak i-1为初始频谱峰值,飞行时间τ i-1为初始飞行时间τ,当测量路拍频信号s b(t)的信噪比不低于第二预设阈值时,初始频谱峰值为测量路拍频信号s b(t)对应的频域信号的峰值;当测量路拍频信号s b(t)的信噪比低于第二预设阈值时,初始频谱峰值为补偿后的拍频信号s' b(t)对应的频域信号的峰值。 Among them, when i=1, peak i-1 is the initial peak of the spectrum, and the flight time τ i-1 is the initial flight time τ. When the signal-to-noise ratio of the measured road beat signal s b (t) is not lower than the second predetermined When the threshold is set, the initial spectral peak value is the peak value of the frequency domain signal corresponding to the measured road beat signal s b (t); when the signal-to-noise ratio of the measured road beat signal s b (t) is lower than the second preset threshold, Sign initial spectrum peak frequency compensated signal s peak 'b (t) corresponding to a frequency domain signal.
S50A、将目标飞行时间确定为飞行时间τ iS50A. Determine the target flight time as the flight time τ i .
再具体地,如图7所示:More specifically, as shown in Figure 7:
S10B、根据飞行时间τ i-1获取延时飞行时间τ' i-1。可选地,延时飞行时间τ' i-1=τ i-1+p,p为预设步长。当然还可以基于其他方式得到延时飞行时间τ' i-1S10B. Obtain the delayed flight time τ'i-1 according to the flight time τ i-1 . Optionally, the delayed flight time τ′ i-1i-1 +p, where p is the preset step length. Of course, the delayed flight time τ′ i-1 can also be obtained based on other methods.
S20B、根据延时飞行时间τ' i-1对激光器发射信号的调频非线性项ε(t)进行延时,以得到测量路回波信号的调频非线性项ε(t-τ' i-1)。换句话说,对激光器发射信号的调频非线性项ε(t)进行一个τ' i-1的延时,得到测量路回波信号的调频非线性项ε(t-τ' i-1)。 S20B. Delay the FM nonlinear term ε(t) of the laser emission signal according to the delay flight time τ′ i-1 to obtain the FM nonlinear term ε(t-τ′ i-1) of the echo signal of the measurement path ). In other words, the FM nonlinear term ε(t) of the laser emission signal is delayed by τ′ i-1 to obtain the FM nonlinear term ε(t-τ′ i-1 ) of the echo signal of the measurement path.
S30B、当测量路拍频信号的信噪比不低于第二预设阈值时,根据激光器发射信号的调频非线性项ε(t)和测量路回波信号的调频非线性项ε(t-τ' i-1)对测量路拍频信号s b(t)进行调频非线性项校正,以得到第一校正拍频信号
Figure PCTCN2020074495-appb-000025
并根据第一校正拍频信号
Figure PCTCN2020074495-appb-000026
获取飞行时间τ i和频谱峰值peak i。其中,第一校正拍频信号为:
Figure PCTCN2020074495-appb-000027
Figure PCTCN2020074495-appb-000028
s b(t)为测量路拍频信号。当测量路拍频信号的信噪比低于第二预设阈值时,根据测量路回波信号的调频非线性项ε(t-τ' i-1)对已做发端调频非线性项补偿后的拍频信号s' b(t)进行接收端调频非线性项校正,以得到第二校正拍频信号
Figure PCTCN2020074495-appb-000029
并根据第二校正拍频信号
Figure PCTCN2020074495-appb-000030
获取飞行时间τ i和频谱峰值peak i。其中,第二校正拍频信号为:
Figure PCTCN2020074495-appb-000031
Figure PCTCN2020074495-appb-000032
S30B. When the signal-to-noise ratio of the measured road beat signal is not lower than the second preset threshold, according to the frequency modulation nonlinear term ε(t) of the laser emission signal and the frequency modulation nonlinear term ε(t- τ' i-1 ) Perform FM nonlinear term correction on the measurement path beat signal s b (t) to obtain the first corrected beat signal
Figure PCTCN2020074495-appb-000025
And according to the first corrected beat signal
Figure PCTCN2020074495-appb-000026
Obtain the flight time τ i and the spectral peak peak i . Among them, the first corrected beat signal is:
Figure PCTCN2020074495-appb-000027
Figure PCTCN2020074495-appb-000028
s b (t) is the measurement road beat frequency signal. When the signal-to-noise ratio of the beat frequency signal of the measurement path is lower than the second preset threshold, according to the FM nonlinearity term ε(t-τ' i-1 ) of the echo signal of the measurement path, after the non-linear term of the transmitter FM is compensated beat signal s' b (t) FM reception terminal nonlinear term correction, the second correction to obtain a beat signal
Figure PCTCN2020074495-appb-000029
And according to the second correction beat signal
Figure PCTCN2020074495-appb-000030
Obtain the flight time τ i and the spectral peak peak i . Among them, the second corrected beat signal is:
Figure PCTCN2020074495-appb-000031
Figure PCTCN2020074495-appb-000032
具体地,根据第一校正拍频信号对应的频域信号确定第一校正拍频信号的频率,并根据第一校正拍频信号的频率确定飞行时间τ i;或者根据第二校正拍频信号对应的频域信号确定第二校正拍频信号的频率,并根据第二校正拍频信号的频率确定飞行时间τ iSpecifically, the frequency of the first corrected beat signal is determined according to the frequency domain signal corresponding to the first corrected beat signal, and the flight time τ i is determined according to the frequency of the first corrected beat signal; or corresponding to the second corrected beat signal Determine the frequency of the second corrected beat signal from the frequency domain signal, and determine the flight time τ i according to the frequency of the second corrected beat signal.
其中,第一校正拍频信号的频率为第一校正拍频信号对应的频域信号的峰值peak i位置对应的频率,或者该频域信号半高宽度中心位置对应的频率;第二校正拍频信号的频率为第二校正拍频信号对应的频域信号的峰值peak i位置对应的频率,或者该频域信号半高宽度 中心位置对应的频率。 The frequency of the first corrected beat signal is the frequency corresponding to the peak i position of the frequency domain signal corresponding to the first corrected beat signal, or the frequency corresponding to the center position of the half-height width of the frequency domain signal; the second corrected beat frequency The frequency of the signal is the frequency corresponding to the peak i position of the frequency domain signal corresponding to the second corrected beat frequency signal, or the frequency corresponding to the center position of the half-height width of the frequency domain signal.
S40B、当i<M时,令i=1+1,并重复执行步骤S10B-S40B;当i=M时,执行步骤S50B。S40B. When i<M, set i=1+1, and repeat steps S10B-S40B; when i=M, perform step S50B.
S50B、目标飞行时间=飞行时间τ iS50B. Target flight time = flight time τ i .
其中,当i=1时,飞行时间τ i-1为初始飞行时间;当i=M时,目标飞行时间为飞行时间τ iAmong them, when i=1, the flight time τ i-1 is the initial flight time; when i=M, the target flight time is the flight time τ i .
在此需要说明的是,在执行如图7所示的方法之前,对i进行初始化,以使i=1。It should be noted here that before the method shown in FIG. 7 is executed, i is initialized so that i=1.
S504、根据目标飞行时间计算得到目标距离。S504: Calculate the target distance according to the target flight time.
其中,目标距离R=cT/2,c为光速,T为目标飞行时间。Among them, the target distance R=cT/2, c is the speed of light, and T is the target flight time.
在一个具体地实施例中,如图8所示,获取参考路拍频信号,具体地将参考路的光信号分为两路,可通过预设固定臂长差为τ D的MZI得到参考路拍频信号s ref(t)。其中一路光信号经过MZI的短臂,另一路光信号经过MZI的长臂,因为两路光信号经过光程不同,所以两路光信在平衡探测器的光敏面上进行相干,得到参考路拍频信号s ref(t)。 In a specific embodiment, as shown in FIG. 8, the beat signal of the reference path is obtained, and the optical signal of the reference path is specifically divided into two paths. The reference path can be obtained by presetting the MZI with the fixed arm length difference τ D. The beat signal s ref (t). One of the optical signals passes through the short arm of the MZI, and the other optical signal passes through the long arm of the MZI. Because the two optical signals pass through different optical paths, the two optical signals are coherent on the photosensitive surface of the balanced detector to obtain the reference path beat frequency The signal s ref (t).
在获取参考路拍频信号s ref(t)后,根据参考路拍频信号s ref(t)计算激光器发射信号的调频频率f(t)。具体地,对该参考路拍频信号s ref(t)进行希尔伯特变换,得到变换后的拍频信号H b(t),该拍频信号可表示为H b(t)=hilbert(s ref(t)),再根据变换后的拍频信号得到参考路拍频信号的瞬时相位
Figure PCTCN2020074495-appb-000033
该瞬时相位
Figure PCTCN2020074495-appb-000034
可表示为:
Figure PCTCN2020074495-appb-000035
其中,H b,imag为H b(t)的虚部,H b,real为H b(t)的实部。
After obtaining the reference path beat frequency signal s ref (t), the frequency modulation frequency f(t) of the laser emission signal is calculated according to the reference path beat frequency signal s ref (t). Specifically, the Hilbert transform is performed on the reference road beat frequency signal s ref (t) to obtain the transformed beat frequency signal H b (t), which can be expressed as H b (t)=hilbert( s ref (t)), and then obtain the instantaneous phase of the beat frequency signal of the reference path according to the transformed beat frequency signal
Figure PCTCN2020074495-appb-000033
The instantaneous phase
Figure PCTCN2020074495-appb-000034
Can be expressed as:
Figure PCTCN2020074495-appb-000035
Among them, H b, imag is the imaginary part of H b (t), and H b, real is the real part of H b (t).
根据激光器发射信号的调频频率f(t)计算得到激光器发射信号的调频非线性项ε(t)和调频非线性度non rate。具体地,根据激光器发射信号的调频频率f(t)与参考路拍频信号的瞬时相位
Figure PCTCN2020074495-appb-000036
的关系计算调频频率f(t),其中,
Figure PCTCN2020074495-appb-000037
再根据调频频率f(t)计算得到激光器发射信号的调频非线性项ε(t)和调频非线性度non rate。激光器发射信号的调频非线性项
Figure PCTCN2020074495-appb-000038
调频非线性度
Figure PCTCN2020074495-appb-000039
其中,f ideal(t)为理想线性调频频率,f ideal(t)=f 0t+0.5αt 2,f 0为载波频率,α为理想线性调频时对应的调频斜率,B为调频带宽。
According to the frequency modulation frequency f(t) of the laser emission signal, the frequency modulation nonlinearity term ε(t) and the frequency modulation nonlinearity non rate of the laser emission signal are calculated. Specifically, according to the instantaneous phase of the frequency modulation frequency f(t) of the laser emission signal and the beat frequency signal of the reference path
Figure PCTCN2020074495-appb-000036
Calculate the FM frequency f(t) from the relationship, where,
Figure PCTCN2020074495-appb-000037
Then calculate the frequency modulation nonlinearity term ε(t) and the frequency modulation nonlinearity non rate of the laser emission signal according to the frequency modulation frequency f(t). Frequency modulation nonlinear term of laser emission signal
Figure PCTCN2020074495-appb-000038
FM nonlinearity
Figure PCTCN2020074495-appb-000039
Among them, f ideal (t) is the ideal chirp frequency, f ideal (t)=f 0 t+0.5αt 2 , f 0 is the carrier frequency, α is the corresponding FM slope when ideal chirp, and B is the FM bandwidth.
可选地,调频非线性度还可以是根据激光器发射信号的调频非线性项(即激光器发射信号的调频频率f(t)与理想线性调频频率f ideal(t)的差)的均方根和调频带宽(B)计算得到,调频非线性度可表示为: Optionally, the frequency modulation nonlinearity can also be based on the root mean square sum of the frequency modulation nonlinearity term of the laser emission signal (that is, the difference between the frequency modulation frequency f(t) of the laser emission signal and the ideal chirp frequency f ideal (t)). The frequency modulation bandwidth (B) is calculated, and the frequency modulation nonlinearity can be expressed as:
Figure PCTCN2020074495-appb-000040
Figure PCTCN2020074495-appb-000040
Figure PCTCN2020074495-appb-000041
Figure PCTCN2020074495-appb-000041
可选地,可根据激光器发射信号的调频频率f(t)与理想线性调频频率f ideal(t)的偏离程度来衡量激光器的调频非线性程度,即用f(t)与f ideal(t)之间的均方根误差(Root Mean  Squard Error,RMSE)来衡量,即 Optionally, the degree of frequency modulation nonlinearity of the laser can be measured according to the degree of deviation between the frequency modulation frequency f(t) of the laser emission signal and the ideal chirp frequency f ideal (t), that is, f(t) and f ideal (t) Root Mean Squard Error (RMSE) between
Figure PCTCN2020074495-appb-000042
Figure PCTCN2020074495-appb-000042
在此需要说明的是,调频频率f(t)和理想线性调频频率f ideal(t)均可看成离散信号,因此f k为调频频率f(t)上第k个离散点对应的频率,f ideal,k为理想线性调频频率f ideal(t)上第k个离散点对应的频率。 It should be noted here that both the FM frequency f(t) and the ideal chirp frequency f ideal (t) can be regarded as discrete signals, so f k is the frequency corresponding to the k-th discrete point on the FM frequency f(t), f ideal, k is the frequency corresponding to the k-th discrete point on the ideal chirp frequency f ideal (t).
获取测量路拍频信号,具体地,将测量路的光信号分为两路,其中,一路信号作为本振信号,另一路信号发射至被测目标物,经目标物表面反射得到回波信号。回波信号与本振信号在平衡探测器的光敏面进行相干,以得到测量路拍频信号s b(t)。 Obtain the beat signal of the measurement path, specifically, divide the optical signal of the measurement path into two paths, where one signal is used as a local oscillator signal, and the other signal is transmitted to the measured target, and the echo signal is obtained by reflection on the target surface. The echo signal and the local oscillator signal are coherent on the photosensitive surface of the balanced detector to obtain the measurement road beat signal s b (t).
对测量路拍频信号进行FFT,计算拍频信号的频率f b和频谱峰值peak last。具体地,当测量路拍频信号的信噪比不低于第二预设阈值时,对测量路拍频信号s b(t)进行FFT,得到测量路拍频信号对应的频域信号,基于该频域信号计算得到测量路拍频信号的频率f b和频谱峰值peak last,并根据测量路拍频信号的频率f b计算得到飞行时间τ;当测量路拍频信号的信噪比低于第二预设阈值时,根据激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行发射端的调频非线性项补偿,以得到补偿后的拍频信号s' b(t);对补偿后的拍频信号s' b(t)进行FFT,得到补偿后的拍频信号对应的频域信号,基于该频域信号计算得到测量路补偿后拍频信号s' b(t)的频率f' b和频谱峰值peak last,此时频谱峰值peak last为补偿后拍频信号对应的频域信号的峰值,并根据测量路补偿后的拍频信号的频率f' b计算得到飞行时间τ。 Perform FFT on the beat frequency signal of the measurement path, and calculate the frequency f b and the spectral peak peak last of the beat frequency signal. Specifically, when the signal-to-noise ratio of the measured road beat signal is not lower than the second preset threshold, perform FFT on the measured road beat signal s b (t) to obtain the frequency domain signal corresponding to the measured road beat signal, based on The frequency domain signal is calculated to obtain the frequency f b of the measured road beat signal and the peak value of the spectrum peak last , and the flight time τ is calculated according to the frequency f b of the measured road beat signal; when the signal to noise ratio of the measured road beat signal is lower than when a second preset threshold, compensating for the nonlinear term FM transmitter according to the FM nonlinear term ε (t) of the signal emitted by the laser measurement channel beat signal, a beat signal s to obtain the compensated 'b (t); of the beat signal s compensated 'b (t) for an FFT, a beat signal corresponding to the frequency-domain signal is compensated for, based on the beat signal s from the frequency domain signal calculated measurement path compensating' b (t) of frequency f 'b and peak last spectral peaks, the spectrum peak at this time is compensated peak last beat frequency peak frequency domain signal corresponding to the signal, and the frequency f in accordance with the beat frequency signal after the measurement path compensating' b is calculated flight time τ .
当调频非线性度non rate大于第一预设阈值Thr nonRate时,根据步长p更新飞行时间τ,更新后的飞行时间为τ=τ+p,对激光器发射信号的调频非线性项ε(t)进行一个更新后的飞行时间τ的延时,得到测量路对应的回波信号的调频非线性项ε(t-τ);当测量路拍频信号的信噪比不低于第二预设阈值时,根据调频非线性项ε(t)和调频非线性项ε(t-τ)对测量路拍频信号s b(t)进行调频非线性项校正,得到第一校正拍频信号
Figure PCTCN2020074495-appb-000043
当测量路拍频信号的信噪比不低于第二预设阈值时,根据测量路回波信号对应的调频非线性项ε(t-τ)对已做发端调频非线性项补偿后的拍频信号s' b(t)进行接收端调频非线性项校正,得到第二校正拍频信号
Figure PCTCN2020074495-appb-000044
When the frequency modulation nonlinearity non rate is greater than the first preset threshold Thr nonRate , the flight time τ is updated according to the step length p, and the updated flight time is τ=τ+p, and the frequency modulation nonlinearity term ε(t ) Perform an updated flight time τ delay to obtain the FM nonlinear term ε(t-τ) of the echo signal corresponding to the measurement path; when the signal-to-noise ratio of the beat frequency signal of the measurement path is not lower than the second preset When the threshold is used, the FM nonlinear term is corrected according to the FM nonlinear term ε(t) and the FM nonlinear term ε(t-τ) on the measurement path beat signal s b (t) to obtain the first corrected beat signal
Figure PCTCN2020074495-appb-000043
When the signal-to-noise ratio of the measured road beat signal is not lower than the second preset threshold, the beat frequency that has been compensated for the originating FM nonlinearity term ε(t-τ) corresponding to the measured road echo signal pilot signal s' b (t) for correcting the nonlinear term FM receiving end, the beat signal to obtain a second corrected
Figure PCTCN2020074495-appb-000044
对第一校正拍频信号
Figure PCTCN2020074495-appb-000045
或第二校正拍频信号
Figure PCTCN2020074495-appb-000046
进行FFT,以得到第一校正拍频信号
Figure PCTCN2020074495-appb-000047
或第二校正拍频信号
Figure PCTCN2020074495-appb-000048
对应的频域信号,基于该频域信号计算第一校正拍频信号
Figure PCTCN2020074495-appb-000049
或第二校正拍频信号
Figure PCTCN2020074495-appb-000050
的频率和频谱峰值peakcurr,根据第一校正拍频信号
Figure PCTCN2020074495-appb-000051
或第二校正拍频信号
Figure PCTCN2020074495-appb-000052
的频率计算得到飞行时间τ。
The first corrected beat signal
Figure PCTCN2020074495-appb-000045
Or the second corrected beat signal
Figure PCTCN2020074495-appb-000046
Perform FFT to get the first corrected beat signal
Figure PCTCN2020074495-appb-000047
Or the second corrected beat signal
Figure PCTCN2020074495-appb-000048
Corresponding frequency domain signal, calculate the first corrected beat frequency signal based on the frequency domain signal
Figure PCTCN2020074495-appb-000049
Or the second corrected beat signal
Figure PCTCN2020074495-appb-000050
The frequency and spectral peak peakcurr, according to the first corrected beat signal
Figure PCTCN2020074495-appb-000051
Or the second corrected beat signal
Figure PCTCN2020074495-appb-000052
The frequency is calculated to get the flight time τ.
当(peak curr-peak last)/peak last>=Thr rate时,对根据第一校正拍频信号
Figure PCTCN2020074495-appb-000053
或第二校正拍频信号
Figure PCTCN2020074495-appb-000054
的频率计算得到的飞行时间τ进行更新,得到更新后的飞行时间,并重复执行上述步骤,直至(peak curr-peak last)/peak last<Thr rate
When (peak curr -peak last )/peak last >= Thr rate , to correct the beat signal according to the first correction
Figure PCTCN2020074495-appb-000053
Or the second corrected beat signal
Figure PCTCN2020074495-appb-000054
Update the flight time τ calculated by the frequency of, to obtain the updated flight time, and repeat the above steps until (peak curr -peak last )/peak last <Thr rate ;
当(peak curr-peak last)/peak last<Thr rate,基于根据校正后的拍频信号的频率计算得到的飞行时间τ计算得到目标距离R=cτ/2,c为光速,即被测目标物与雷达之间的距离。 When (peak curr -peak last )/peak last <Thr rate , the target distance R=cτ/2 is calculated based on the flight time τ calculated according to the frequency of the corrected beat signal, and c is the speed of light, that is, the measured target The distance from the radar.
图9a表明,本申请调频非线性校正方法在调频非线性情况下有效提升了距离估计的精度。从图9b中可以看出,经过本申请非线项迭代校正后,拍频信号的频谱峰值逐渐增大且 靠近理想值,故本申请非线性校正方法在调频非线性情况下可以有效提升拍频信号频域信噪比,也即可以有效提升最远测距范围。Figure 9a shows that the FM nonlinearity correction method of the present application effectively improves the accuracy of distance estimation in the case of FM nonlinearity. It can be seen from Fig. 9b that after the iterative correction of the nonlinear term of this application, the spectral peak of the beat signal gradually increases and is close to the ideal value. Therefore, the nonlinear correction method of this application can effectively increase the beat frequency in the case of frequency modulation nonlinearity. The signal-to-noise ratio in the frequency domain of the signal can effectively increase the farthest ranging range.
可以看出,在本申请实施例中,通过根据激光器发射信号的调频非线性项和测量路回波信号的调频非线性项对测量路拍频信号进行调频非线性校正,可在很大程度上消除调频非线性对测量路拍频信号的影响,提高了基于校正后的拍频信号得到的飞行时间的精确性,进而提高了测距精度。It can be seen that, in the embodiment of the present application, the frequency modulation nonlinearity of the measurement path signal can be corrected according to the frequency modulation nonlinearity term of the laser emission signal and the frequency modulation nonlinearity term of the measurement path echo signal. Eliminate the influence of non-linear frequency modulation on the measured road beat signal, improve the accuracy of the flight time based on the corrected beat signal, and further improve the ranging accuracy.
参见图10,图10为本申请实施例提供的一种基于调频非线性校正的测距装置的结构示意图。如图10所示,该测距装置1000,包括:Refer to FIG. 10, which is a schematic structural diagram of a distance measurement device based on frequency modulation nonlinear correction provided by an embodiment of the application. As shown in FIG. 10, the distance measuring device 1000 includes:
获取单元1001,用于获取参考路拍频信号和测量路拍频信号,参考路拍频信号和测量路拍频信号是分别基于激光器发射信号得到的;The acquiring unit 1001 is used to acquire the reference road beat signal and the measurement road beat signal, the reference road beat signal and the measurement road beat signal are obtained based on the laser emission signal respectively;
计算单元1002,用于根据参考路拍频信号计算得到调频非线性度;并根据测量路拍频信号计算得到初始飞行时间τ;The calculation unit 1002 is configured to calculate the frequency modulation nonlinearity according to the reference road beat signal; and calculate the initial flight time τ according to the measured road beat signal;
计算单元1002,还用于在确定单元1004确定调频非线性度不小于第一预设阈值时,根据参考路拍频信号计算得到激光器发射信号的调频非线性项ε(t),The calculation unit 1002 is further configured to calculate the frequency modulation nonlinearity term ε(t) of the laser emission signal according to the reference path beat signal when the determining unit 1004 determines that the frequency modulation nonlinearity is not less than the first preset threshold,
校正单元1003,用于根据初始飞行时间τ和激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行非线性迭代校正,以得到目标飞行时间;The correction unit 1003 is configured to perform nonlinear iterative correction on the measurement path beat frequency signal according to the initial flight time τ and the frequency modulation nonlinear term ε(t) of the laser emission signal to obtain the target flight time;
计算单元1002,用于根据目标飞行时间计算得到目标距离。The calculation unit 1002 is configured to calculate the target distance according to the target flight time.
在一个可行的实施例中,校正单元1003具体用于:In a feasible embodiment, the correction unit 1003 is specifically configured to:
S1:根据延时飞行时间τ' i-1对激光器发射信号的调频非线性项ε(t)进行延时,以得到测量路回波信号对应的调频非线性项ε(t-τ' i-1),延时飞行时间τ' i-1是基于飞行时间τ i-1得到的; S1: Delay the FM nonlinear term ε(t) of the laser emission signal according to the delay flight time τ′ i-1 to obtain the FM nonlinear term ε(t-τ′ i-) corresponding to the echo signal of the measurement path. 1 ), the delayed flight time τ'i -1 is obtained based on the flight time τ i-1 ;
S2:根据测量路回波信号对应的调频非线性项ε(t-τ' i-1)对测量路拍频信号进行调频非线性校正,以得到校正拍频信号; S2: Perform FM nonlinear correction on the beat frequency signal of the measurement path according to the FM nonlinear term ε(t-τ'i -1) corresponding to the echo signal of the measurement path to obtain a corrected beat frequency signal;
S3、根据校正拍频信号确定目标飞行时间。S3. Determine the target flight time according to the corrected beat frequency signal.
在一个可行的实施例中,在根据校正拍频信号确定目标飞行时间的方面,校正单元1003具体用于:In a feasible embodiment, in terms of determining the target flight time according to the corrected beat signal, the correction unit 1003 is specifically configured to:
根据校正拍频信号获取飞行时间τ i;当频谱峰值peak i满足预设条件时,将飞行时间τ i确定为目标飞行时间;当频谱峰值peak i不满足预设条件时,令i=i+1,并重复执行上述步骤S1-S3;频谱峰值peak i为校正拍频信号对应的频域信号的峰值; Obtain the flight time τ i according to the corrected beat signal; when the peak frequency peak i meets the preset condition, the flight time τ i is determined as the target flight time; when the peak frequency peak i does not meet the preset condition, let i=i+ 1. Repeat the above steps S1-S3; the peak value of the spectrum peak i is the peak value of the frequency domain signal corresponding to the corrected beat signal;
当i=1时,peak i-1为所述初始频谱峰值,飞行时间τ i-1为初始飞行时间τ。 When i=1, peak i-1 is the initial frequency spectrum peak, and the flight time τ i-1 is the initial flight time τ.
在一个可行的实施例中,预设条件为(peak i-peak i-1)/peak i-1<Thr rate,该Thr rate为第三预设阈值。 In one possible embodiment, the predetermined condition (peak i -peak i-1) / peak i-1 <Thr rate, Thr rate as the third predetermined threshold value.
在一个可行的实施例中,校正单元1003具体用于,包括:In a feasible embodiment, the correction unit 1003 is specifically configured to include:
根据初始飞行时间τ和激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行M次非线性迭代计算,以得到目标飞行时间,M为大于1的整数,According to the initial flight time τ and the frequency modulation non-linear term ε(t) of the laser emission signal, perform M non-linear iterative calculations on the measurement path beat frequency signal to obtain the target flight time, M is an integer greater than 1,
其中,在进行第i次非线性迭代计算时,根据延时飞行时间τ' i-1对激光器发射信号的调频非线性项ε(t)进行延时,以得到测量路回波信号对应的调频非线性项ε(t-τ' i-1);延时 飞行时间τ' i-1是基于飞行时间τ i-1得到的; Among them, when the i-th nonlinear iterative calculation is performed, the frequency modulation nonlinear term ε(t) of the laser emission signal is delayed according to the delay flight time τ'i-1 to obtain the frequency modulation corresponding to the echo signal of the measurement path The non-linear term ε(t-τ'i -1 ); the delayed flight time τ'i -1 is obtained based on the flight time τi -1 ;
根据测量路回波信号对应的调频非线性项ε(t-τ' i-1)对测量路拍频信号进行调频非线性校正,以得到校正拍频信号;根据校正拍频信号获取飞行时间τ iPerform FM nonlinear correction on the measurement path beat signal according to the FM nonlinear term ε(t-τ' i-1 ) corresponding to the echo signal of the measurement path to obtain the corrected beat signal; obtain the flight time τ according to the corrected beat signal i ;
其中,当i=1时,飞行时间τ i-1为初始飞行时间τ;当i=M时,目标飞行时间为飞行时间τ iAmong them, when i=1, the flight time τ i-1 is the initial flight time τ; when i=M, the target flight time is the flight time τ i .
在一个可行的实施例中,在根据校正拍频信号获取飞行时间τ i的方面,计算单元1002具体用于: In a feasible embodiment, in terms of obtaining the flight time τ i according to the corrected beat signal, the calculation unit 1002 is specifically configured to:
对校正拍频信号进行FFT,以得到校正拍频信号对应的频域信号;根据校正拍频信号对应的频域信号获取飞行时间τ iPerform FFT on the corrected beat frequency signal to obtain the frequency domain signal corresponding to the corrected beat frequency signal; obtain the flight time τ i according to the frequency domain signal corresponding to the corrected beat frequency signal.
在一个可行的实施例中,校正拍频信号包括第一校正拍频信号或第二校正拍频信号;在根据测量路回波信号对应的调频非线性项ε(t-τ' i-1)对测量路拍频信号进行调频非线性校正,以得到校正拍频信号的方面,校正单元1003具体用于: In a feasible embodiment, the corrected beat frequency signal includes the first corrected beat frequency signal or the second corrected beat frequency signal; according to the frequency modulation nonlinearity term ε(t-τ' i-1 ) corresponding to the echo signal of the measurement path Performing FM nonlinear correction on the beat signal of the measurement path to obtain the aspect of correcting the beat signal, the correction unit 1003 is specifically used for:
当测量路拍频信号的信噪比不低于第二预设阈值时,根据激光器发射信号的调频非线性项ε(t)和测量路回波信号对应的调频非线性项ε(t-τ' i-1)对测量路拍频信号进行调频非线性项校正,以得到第一校正拍频信号;当测量路拍频信号的信噪比低于第二预设阈值时,根据测量路回波信号对应的调频非线性项ε(t-τ' i-1)对补偿后的拍频信号进行接收端调频非线性项校正,以得到第二校正拍频信号; When the signal-to-noise ratio of the beat frequency signal of the measurement path is not lower than the second preset threshold, the frequency modulation nonlinearity term ε(t) of the laser emission signal and the frequency modulation nonlinearity term ε(t-τ) corresponding to the echo signal of the measurement path ' I-1 ) Perform FM nonlinear term correction on the beat signal of the measurement road to obtain the first corrected beat signal; when the signal-to-noise ratio of the beat signal of the measurement road is lower than the second preset threshold, according to the measurement road The FM nonlinear term ε(t-τ′i -1 ) corresponding to the wave signal performs receiving end FM nonlinear term correction on the compensated beat signal to obtain a second corrected beat signal;
其中,补偿后的拍频信号是基于根据激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行发射端的调频非线性项补偿得到的。Among them, the compensated beat signal is obtained based on the FM nonlinear term ε(t) of the laser emitting signal to compensate the FM nonlinear term of the transmitter of the measurement path beat signal.
在一个可行的实施例中,在根据测量路拍频信号计算得到初始飞行时间τ的方面,计算单元1002具体用于:In a feasible embodiment, in terms of calculating the initial flight time τ according to the measured road beat signal, the calculation unit 1002 is specifically configured to:
当测量路拍频信号的信噪比不低于第二预设阈值时,对测量路拍频信号进行FFT,以得到测量路拍频信号对应的频域信号;根据测量路拍频信号对应的频域信号确定测量路拍频信号的频率,并根据测量路拍频信号的频率计算得到初始飞行时间τ;When the signal-to-noise ratio of the measured road beat signal is not lower than the second preset threshold, perform FFT on the measured road beat signal to obtain the frequency domain signal corresponding to the measured road beat signal; The frequency domain signal determines the frequency of the measured road beat signal, and calculates the initial flight time τ according to the frequency of the measured road beat signal;
其中,测量路拍频信号的频率为测量路拍频信号对应的频域信号的峰值点位置对应的频率,或者频域信号的半高宽度中点位置对应的频率。Wherein, the frequency of the measurement road beat signal is the frequency corresponding to the peak point position of the frequency domain signal corresponding to the measurement road beat signal, or the frequency corresponding to the midpoint position of the half-height width of the frequency domain signal.
在一个可行的实施例中,在根据测量路拍频信号计算得到初始飞行时间τ的方面,计算单元1002具体用于:In a feasible embodiment, in terms of calculating the initial flight time τ according to the measured road beat signal, the calculation unit 1002 is specifically configured to:
当测量路拍频信号的信噪比低于第二预设阈值时;根据激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行发射端的调频非线性项补偿,以得到补偿后的拍频信号;激光器发射信号的调频非线性项是基于参考路拍频信号计算得到的;When the signal-to-noise ratio of the measurement road beat signal is lower than the second preset threshold; according to the FM nonlinear term ε(t) of the laser emitting signal, the measurement road beat signal is compensated by the FM nonlinear term at the transmitting end. After the beat frequency signal; the frequency modulation nonlinear term of the laser emission signal is calculated based on the beat frequency signal of the reference path;
对补偿后的拍频信号进行FFT,以得到补偿后的拍频信号对应的频域信号;根据补偿后的拍频信号对应的频域信号确定补偿后的拍频信号的频率,并根据补偿后的拍频信号的频率计算得到初始飞行时间τ;Perform FFT on the compensated beat signal to obtain the frequency domain signal corresponding to the compensated beat signal; determine the frequency of the compensated beat signal according to the frequency domain signal corresponding to the compensated beat signal, and according to the compensated beat signal Calculate the frequency of the beat signal to get the initial flight time τ;
其中,补偿后的拍频信号的频率为补偿后的拍频信号对应的频域信号的峰值点位置对应的频率,或者频域信号的半高宽度中点位置对应的频率。The frequency of the compensated beat signal is the frequency corresponding to the peak point position of the frequency domain signal corresponding to the compensated beat signal, or the frequency corresponding to the midpoint position of the half-height width of the frequency domain signal.
在一个可行的实施例中,在根据参考路拍频信号计算得到激光器发射信号的调频非线性项ε(t)的方面,计算单元1002具体用于:In a feasible embodiment, the calculation unit 1002 is specifically configured to:
对参考路拍频信号进行希尔伯特变换,以得到变换后的拍频信号;根据该变换后的拍频信号计算得到激光器发射信号的调频频率;根据激光器发射信号的调频频率和理想线性调频频率计算得到激光器发射信号的调频非线性项ε(t)。Perform Hilbert transform on the beat signal of the reference path to obtain the transformed beat signal; calculate the FM frequency of the laser emission signal according to the transformed beat signal; According to the FM frequency and ideal chirp of the laser emission signal Frequency calculation obtains the frequency modulation nonlinear term ε(t) of the laser emission signal.
需要说明的是,上述各单元(获取单元1001、计算单元1002、校正单元1003和确定单元1004)用于执行上述方法的相关步骤。比如获取单元1001用于执行步骤S501的相关内容,计算单元1002用于执行步骤S502和S504的相关内容,校正单元1003和确定单元1004用于执行步骤S503的相关内容。It should be noted that the aforementioned units (the acquisition unit 1001, the calculation unit 1002, the correction unit 1003, and the determination unit 1004) are used to execute the relevant steps of the aforementioned method. For example, the acquiring unit 1001 is used to execute the related content of step S501, the calculating unit 1002 is used to execute the related content of steps S502 and S504, and the correction unit 1003 and the determining unit 1004 are used to execute the related content of step S503.
在本实施例中,基于调频非线性校正的测距装置1000是以单元的形式来呈现。这里的“单元”可以指特定应用集成电路(application-specific integrated circuit,ASIC),执行一个或多个软件或固件程序的处理器和存储器,集成逻辑电路,和/或其他可以提供上述功能的器件。此外,以上获取单元1001、计算单元1002、校正单元1003和确定单元1004可通过图11所示的基于调频非线性校正的测距装置的处理器1101来实现。In this embodiment, the ranging device 1000 based on frequency modulation nonlinear correction is presented in the form of a unit. The "unit" here can refer to an application-specific integrated circuit (ASIC), a processor and memory that executes one or more software or firmware programs, an integrated logic circuit, and/or other devices that can provide the above-mentioned functions . In addition, the above acquisition unit 1001, calculation unit 1002, correction unit 1003, and determination unit 1004 may be implemented by the processor 1101 of the distance measurement device based on frequency modulation nonlinear correction shown in FIG. 11.
如图11所示测距装置1100可以以图11中的结构来实现,该测距装置1100包括至少一个处理器1101,至少一个存储器1102以及至少一个通信接口1103。所述处理器1101、所述存储器1102和所述通信接口1103通过所述通信总线连接并完成相互问的通信。The distance measuring device 1100 shown in FIG. 11 can be implemented with the structure in FIG. 11, and the distance measuring device 1100 includes at least one processor 1101, at least one memory 1102 and at least one communication interface 1103. The processor 1101, the memory 1102, and the communication interface 1103 are connected through the communication bus and complete mutual communication.
处理器1101可以是通用中央处理器(CPU),微处理器,特定应用集成电路(application-specific integrated circuit,ASIC),或一个或多个用于控制以上方案程序执行的集成电路。The processor 1101 may be a general-purpose central processing unit (CPU), a microprocessor, an application-specific integrated circuit (ASIC), or one or more integrated circuits for controlling the execution of the above program programs.
通信接口1103,用于与其他设备或通信网络通信,如以太网,无线接入网(RAN),无线局域网(Wireless Local Area Networks,WLAN)等。The communication interface 1103 is used to communicate with other devices or communication networks, such as Ethernet, radio access network (RAN), wireless local area network (Wireless Local Area Networks, WLAN), etc.
存储器1102可以是只读存储器(read-only memory,ROM)或可存储静态信息和指令的其他类型的静态存储设备,随机存取存储器(random access memory,RAM)或者可存储信息和指令的其他类型的动态存储设备,也可以是电可擦可编程只读存储器(Electrically Erasable Programmable Read-Only Memory,EEPROM)、只读光盘(Compact Disc Read-Only Memory,CD-ROM)或其他光盘存储、光碟存储(包括压缩光碟、激光碟、光碟、数字通用光碟、蓝光光碟等)、磁盘存储介质或者其他磁存储设备、或者能够用于携带或存储具有指令或数据结构形式的期望的程序代码并能够由计算机存取的任何其他介质,但不限于此。存储器可以是独立存在,通过总线与处理器相连接。存储器也可以和处理器集成在一起。The memory 1102 may be a read-only memory (ROM) or other types of static storage devices that can store static information and instructions, random access memory (RAM), or other types that can store information and instructions The dynamic storage device can also be electrically erasable programmable read-only memory (Electrically Erasable Programmable Read-Only Memory, EEPROM), CD-ROM (Compact Disc Read-Only Memory, CD-ROM) or other optical disc storage, optical disc storage (Including compact discs, laser discs, optical discs, digital versatile discs, Blu-ray discs, etc.), magnetic disk storage media or other magnetic storage devices, or can be used to carry or store desired program codes in the form of instructions or data structures and can be used by a computer Any other media accessed, but not limited to this. The memory can exist independently and is connected to the processor through a bus. The memory can also be integrated with the processor.
其中,所述存储器1102用于存储执行以上方案的应用程序代码,并由处理器1101来控制执行。所述处理器1101用于执行所述存储器1102中存储的应用程序代码。Wherein, the memory 1102 is used to store application program codes for executing the above solutions, and the processor 1101 controls the execution. The processor 1101 is configured to execute application program codes stored in the memory 1102.
存储器1102存储的代码可执行以上提供的基于调频非线性校正的测距方法,比如获取参考路拍频信号和测量路拍频信号;根据参考路拍频信号计算得到调频非线性度;并根据测量路拍频信号计算得到初始飞行时间τ;当调频非线性度不小于第一预设阈值时,根据参考路拍频信号计算得到激光器发射信号的调频非线性项ε(t),根据初始飞行时间τ和激光器发射信号的调频非线性项ε(t)对测量路拍频信号进行非线性迭代校正,以得到目标飞行时间;根据目标飞行时间计算得到目标距离。The code stored in the memory 1102 can execute the above-provided ranging method based on FM nonlinear correction, such as obtaining the reference road beat signal and measuring road beat signal; calculating the FM nonlinearity based on the reference road beat signal; and according to the measurement The initial flight time τ is calculated from the road beat frequency signal; when the frequency modulation nonlinearity is not less than the first preset threshold, the frequency modulation nonlinear term ε(t) of the laser emission signal is calculated according to the reference road beat frequency signal, according to the initial flight time τ and the frequency modulation non-linear term ε(t) of the laser emission signal perform non-linear iterative correction on the measurement road beat signal to obtain the target flight time; the target distance is calculated according to the target flight time.
本申请实施例还提供一种计算机存储介质,其中,该计算机存储介质可存储有程序, 该程序执行时包括上述方法实施例中记载的任何一种基于调频非线性校正的测距方法的部分或全部步骤。An embodiment of the present application further provides a computer storage medium, wherein the computer storage medium may store a program, and when the program is executed, it includes any part or All steps.
需要说明的是,对于前述的各方法实施例,为了简单描述,故将其都表述为一系列的动作组合,但是本领域技术人员应该知悉,本申请并不受所描述的动作顺序的限制,因为依据本申请,某些步骤可以采用其他顺序或者同时进行。其次,本领域技术人员也应该知悉,说明书中所描述的实施例均属于优选实施例,所涉及的动作和模块并不一定是本申请所必须的。It should be noted that for the foregoing method embodiments, for the sake of simple description, they are all expressed as a series of action combinations, but those skilled in the art should know that this application is not limited by the described sequence of actions. Because according to this application, some steps can be performed in other order or at the same time. Secondly, those skilled in the art should also know that the embodiments described in the specification are all preferred embodiments, and the actions and modules involved are not necessarily required by this application.
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述的部分,可以参见其他实施例的相关描述。In the above-mentioned embodiments, the description of each embodiment has its own emphasis. For parts that are not described in detail in an embodiment, reference may be made to related descriptions of other embodiments.
在本申请所提供的几个实施例中,应该理解到,所揭露的装置,可通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性或其它的形式。In the several embodiments provided in this application, it should be understood that the disclosed device may be implemented in other ways. For example, the device embodiments described above are merely illustrative, for example, the division of the units is only a logical function division, and there may be other divisions in actual implementation, for example, multiple units or components may be combined or may be Integrate into another system, or some features can be ignored or not implemented. In addition, the displayed or discussed mutual coupling or direct coupling or communication connection may be indirect coupling or communication connection through some interfaces, devices or units, and may be in electrical or other forms.
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。The units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the objectives of the solutions of the embodiments.
另外,在本申请各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。In addition, the functional units in the various embodiments of the present application may be integrated into one processing unit, or each unit may exist alone physically, or two or more units may be integrated into one unit. The above-mentioned integrated unit can be implemented in the form of hardware or software functional unit.
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储器中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储器中,包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储器包括:U盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer readable memory. Based on this understanding, the technical solution of the present application essentially or the part that contributes to the existing technology or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a memory. A number of instructions are included to enable a computer device (which may be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods described in the various embodiments of the present application. The aforementioned memory includes: U disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), mobile hard disk, magnetic disk or optical disk and other media that can store program codes.
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,该程序可以存储于一计算机可读存储器中,存储器可以包括:闪存盘、只读存储器(英文:Read-Only Memory,简称:ROM)、随机存取器(英文:Random Access Memory,简称:RAM)、磁盘或光盘等。Those of ordinary skill in the art can understand that all or part of the steps in the various methods of the above-mentioned embodiments can be completed by a program instructing relevant hardware. The program can be stored in a computer-readable memory, and the memory can include: a flash disk , Read-only memory (English: Read-Only Memory, abbreviation: ROM), random access device (English: Random Access Memory, abbreviation: RAM), magnetic disk or optical disc, etc.
以上对本申请实施例进行了详细介绍,本文中应用了具体个例对本申请的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本申请的方法及其核心思想;同时,对于本领域的一般技术人员,依据本申请的思想,在具体实施方式及应用范围上均会有改变之处,综上上述,本说明书内容不应理解为对本申请的限制。The embodiments of the application are described in detail above, and specific examples are used in this article to illustrate the principles and implementation of the application. The descriptions of the above embodiments are only used to help understand the methods and core ideas of the application; at the same time, for Those of ordinary skill in the art, based on the ideas of the application, will have changes in the specific implementation and the scope of application. In summary, the content of this specification should not be construed as limiting the application.

Claims (23)

  1. 一种基于调频非线性校正的测距方法,其特征在于,包括:A ranging method based on frequency modulation nonlinear correction, which is characterized in that it comprises:
    获取参考路拍频信号和测量路拍频信号;所述参考路拍频信号和测量路拍频信号是分别基于激光器发射信号得到的;Acquiring a reference road beat signal and a measurement road beat signal; the reference road beat signal and the measurement road beat signal are respectively obtained based on the laser emission signal;
    根据所述参考路拍频信号计算得到调频非线性度;并根据所述测量路拍频信号计算得到初始飞行时间τ;Calculate the frequency modulation nonlinearity according to the reference road beat signal; and calculate the initial flight time τ according to the measured road beat signal;
    确定所述调频非线性度不小于第一预设阈值,根据所述参考路拍频信号计算得到所述激光器发射信号的调频非线性项ε(t),根据所述初始飞行时间τ和所述激光器发射信号的调频非线性项ε(t)对所述测量路拍频信号进行非线性迭代校正,以得到目标飞行时间;It is determined that the frequency modulation nonlinearity is not less than a first preset threshold, the frequency modulation nonlinearity term ε(t) of the laser emission signal is calculated according to the reference path beat frequency signal, and the frequency modulation nonlinearity term ε(t) is calculated according to the initial flight time τ and the The frequency modulation nonlinear term ε(t) of the laser emission signal performs nonlinear iterative correction on the measurement path beat frequency signal to obtain the target flight time;
    根据所述目标飞行时间计算得到目标距离。The target distance is calculated according to the target flight time.
  2. 根据权利要求1所述的方法,其特征在于,所述根据初始飞行时间τ和所述激光器发射信号的调频非线性项ε(t)对所述测量路拍频信号进行非线性迭代校正,以得到目标飞行时间,包括:The method according to claim 1, characterized in that the non-linear iterative correction is performed on the measurement path beat frequency signal according to the initial flight time τ and the frequency modulation non-linear term ε(t) of the laser emission signal to Get the target flight time, including:
    S1:根据所述延时飞行时间τ′ i-1对所述激光器发射信号的调频非线性项ε(t)进行延时,以得到所述测量路回波信号对应的调频非线性项ε(t-τ′ i-1),所述延时飞行时间τ′ i-1是基于所述飞行时间τ i-1得到的; S1: Delay the FM nonlinear term ε(t) of the laser emission signal according to the delay flight time τ′i -1 to obtain the FM nonlinear term ε(t) corresponding to the echo signal of the measurement path t-τ′ i-1 ), the delayed flight time τ′ i-1 is obtained based on the flight time τ i-1 ;
    S2:根据所述测量路回波信号对应的调频非线性项ε(t-τ′ i-1)对所述测量路拍频信号进行调频非线性校正,以得到校正拍频信号; S2: Perform FM nonlinear correction on the measurement path beat signal according to the FM nonlinear term ε(t-τ′ i-1) corresponding to the measurement path echo signal to obtain a corrected beat signal;
    S3、根据所述校正拍频信号确定所述目标飞行时间。S3. Determine the target flight time according to the corrected beat frequency signal.
  3. 根据权利要求2所述的方法,其特征在于,所述根据所述校正拍频信号确定所述目标飞行时间,包括:The method according to claim 2, wherein the determining the target flight time according to the corrected beat frequency signal comprises:
    根据所述校正拍频信号获取飞行时间τ i Obtain the flight time τ i according to the corrected beat frequency signal;
    当所述频谱峰值peak i满足预设条件时,将所述飞行时间τ i确定为所述目标飞行时间;当所述频谱峰值peak i不满足所述预设条件时,令i=i+1,并重复执行步骤S1-S3;所述频谱峰值peak i为所述校正拍频信号对应的频域信号的峰值; When the spectral peak i meets a preset condition, the flight time τ i is determined as the target flight time; when the spectral peak i does not meet the preset condition, let i=i+1 , And repeat steps S1-S3; the spectrum peak peak i is the peak value of the frequency domain signal corresponding to the corrected beat signal;
    当所述i=1时,所述peak i-1为所述初始频谱峰值,所述飞行时间τ i-1为所述初始飞行时间τ。 When i=1, the peak i-1 is the initial frequency spectrum peak, and the flight time τ i-1 is the initial flight time τ.
  4. 根据权利要求3所述的方法,其特征在于,所述预设条件为(peak i-peak i-1)/peak i-1<Thr rate,所述Thr rate为第三预设阈值。 The method according to claim 3, wherein said predetermined condition is a (peak i -peak i-1) / peak i-1 <Thr rate, the Thr rate is a third predetermined threshold value.
  5. 根据权利要求1所述的方法,其特征在于,所述根据初始飞行时间τ和所述激光器发射信号的调频非线性项ε(t)对所述测量路拍频信号进行非线性迭代校正,以得到目标飞行时间,包括:The method according to claim 1, characterized in that the non-linear iterative correction is performed on the measurement path beat frequency signal according to the initial flight time τ and the frequency modulation non-linear term ε(t) of the laser emission signal to Get the target flight time, including:
    根据所述初始飞行时间τ和所述激光器发射信号的调频非线性项ε(t)对所述测量路拍频信号进行M次非线性迭代计算,以得到所述目标飞行时间,所述M为大于1的整数,According to the initial flight time τ and the frequency modulation nonlinear term ε(t) of the laser emission signal, perform M non-linear iterative calculations on the measurement path beat frequency signal to obtain the target flight time, where M is An integer greater than 1,
    其中,在进行第i次非线性迭代计算时,根据所述延时飞行时间τ′ i-1对所述激光器发射信号的调频非线性项ε(t)进行延时,以得到测量路回波信号对应的调频非线性项ε(t-τ′ i-1);所述延时飞行时间τ′ i-1是基于所述飞行时间τ i-1得到的; Wherein, when the i-th nonlinear iterative calculation is performed, the frequency modulation nonlinear term ε(t) of the laser emission signal is delayed according to the delay flight time τ′ i-1 to obtain the measurement path echo The FM nonlinear term ε(t-τ′ i-1 ) corresponding to the signal; the delayed flight time τ′ i-1 is obtained based on the flight time τ i-1 ;
    根据所述测量路回波信号对应的调频非线性项ε(t-τ′ i-1)对所述测量路拍频信号进行调频非线性校正,以得到校正拍频信号;根据所述校正拍频信号获取飞行时间τ iPerform FM nonlinear correction on the measurement path beat signal according to the FM nonlinear term ε(t-τ′ i-1 ) corresponding to the echo signal of the measurement path to obtain a corrected beat signal; Frequency signal to obtain flight time τ i ;
    其中,当i=1时,所述飞行时间τ i-1为所述初始飞行时间τ;当所述i=M时,所述目标飞行时间为飞行时间τ iWherein, when i=1, the flight time τ i-1 is the initial flight time τ; when i=M, the target flight time is the flight time τ i .
  6. 根据权利要求3-5任一项所述的方法,其特征在于,其特征在于,所述根据所述校正拍频信号获取飞行时间τ i,包括: The method according to any one of claims 3 to 5, characterized in that the obtaining the flight time τ i according to the corrected beat frequency signal comprises:
    对所述校正拍频信号进行快速傅里叶变换FFT,以得到所述校正拍频信号对应的频域信号;Performing fast Fourier transform FFT on the corrected beat frequency signal to obtain a frequency domain signal corresponding to the corrected beat frequency signal;
    根据所述校正拍频信号对应的频域信号获取飞行时间τ i Obtain the flight time τ i according to the frequency domain signal corresponding to the corrected beat frequency signal.
  7. 根据权利要求2-6任一项所述的方法,其特征在于,所述校正拍频信号包括第一校正拍频信号或第二校正拍频信号;所述根据所述测量路回波信号对应的调频非线性项ε(t-τ′ i-1)对所述测量路拍频信号进行调频非线性校正,以得到校正拍频信号,包括: The method according to any one of claims 2-6, wherein the corrected beat frequency signal comprises a first corrected beat frequency signal or a second corrected beat frequency signal; the echo signal corresponding to the measurement path corresponds to The FM nonlinear term ε(t-τ′ i-1 ) performs FM nonlinear correction on the beat frequency signal of the measurement path to obtain a corrected beat frequency signal, including:
    当所述测量路拍频信号的信噪比不低于第二预设阈值时,根据所述激光器发射信号的调频非线性项ε(t)和所述测量路回波信号对应的调频非线性项ε(t-τ′ i-1)对所述测量路拍频信号进行调频非线性项校正,以得到所述第一校正拍频信号; When the signal-to-noise ratio of the measurement path beat signal is not lower than the second preset threshold, the frequency modulation nonlinearity term ε(t) of the laser emission signal and the frequency modulation nonlinearity corresponding to the measurement path echo signal The term ε(t-τ′ i-1 ) performs frequency modulation nonlinear term correction on the measurement path beat signal to obtain the first corrected beat signal;
    当所述测量路拍频信号的信噪比低于所述第二预设阈值时,根据所述测量路回波信号对应的调频非线性项ε(t-τ′ i-1)对所述补偿后的拍频信号进行接收端调频非线性项校正,以得到所述第二校正拍频信号; When the signal-to-noise ratio of the measurement path beat frequency signal is lower than the second preset threshold value, the frequency modulation nonlinearity term ε(t-τ′ i-1 ) corresponding to the measurement path echo signal is compared to the The compensated beat frequency signal is corrected by a frequency modulation nonlinear term at the receiving end to obtain the second corrected beat frequency signal;
    其中,所述补偿后的拍频信号是基于根据所述激光器发射信号的调频非线性项ε(t)对所述测量路拍频信号进行发射端的调频非线性项补偿得到的。Wherein, the compensated beat frequency signal is obtained based on the frequency modulation nonlinearity term ε(t) of the measurement path beat frequency signal at the transmitting end according to the frequency modulation nonlinearity term ε(t) of the laser transmitting signal.
  8. 根据权利要求7所述的方法,其特征在于,所述根据所述测量路拍频信号计算得到初始飞行时间τ,包括:The method according to claim 7, wherein the calculating the initial flight time τ according to the measured road beat frequency signal comprises:
    当所述测量路拍频信号的信噪比不低于第二预设阈值时,对所述测量路拍频信号进行FFT,以得到所述测量路拍频信号对应的频域信号;When the signal-to-noise ratio of the measurement road beat signal is not lower than a second preset threshold, perform an FFT on the measurement road beat signal to obtain a frequency domain signal corresponding to the measurement road beat signal;
    根据所述测量路拍频信号对应的频域信号确定所述测量路拍频信号的频率,并根据所述测量路拍频信号的频率计算得到所述初始飞行时间τ;Determine the frequency of the measurement road beat signal according to the frequency domain signal corresponding to the measurement road beat signal, and calculate the initial flight time τ according to the frequency of the measurement road beat signal;
    其中,所述测量路拍频信号的频率为所述测量路拍频信号对应的频域信号的峰值点位置对应的频率,或者所述频域信号的半高宽度中点位置对应的频率。Wherein, the frequency of the measurement path beat signal is a frequency corresponding to a peak point position of a frequency domain signal corresponding to the measurement path beat signal, or a frequency corresponding to a midpoint position of a half-height width of the frequency domain signal.
  9. 根据权利要求7所述的方法,其特征在于,所述根据所述测量路拍频信号计算得到初始飞行时间τ,包括:The method according to claim 7, wherein the calculating the initial flight time τ according to the measured road beat frequency signal comprises:
    当所述测量路拍频信号的信噪比低于第二预设阈值时;根据所述激光器发射信号的调 频非线性项ε(t)对所述测量路拍频信号进行发射端的调频非线性项补偿,以得到补偿后的拍频信号;激光器发射信号的调频非线性项是基于所述参考路拍频信号计算得到的;When the signal-to-noise ratio of the measurement path beat frequency signal is lower than the second preset threshold; perform the frequency modulation nonlinearity at the transmitting end of the measurement path beat frequency signal according to the frequency modulation nonlinearity term ε(t) of the laser emission signal Term compensation to obtain the compensated beat frequency signal; the frequency modulation nonlinear term of the laser emission signal is calculated based on the reference path beat frequency signal;
    对所述补偿后的拍频信号进行FFT,以得到所述补偿后的拍频信号对应的频域信号;Performing FFT on the compensated beat signal to obtain a frequency domain signal corresponding to the compensated beat signal;
    根据所述补偿后的拍频信号对应的频域信号确定所述补偿后的拍频信号的频率,并根据所述补偿后的拍频信号的频率计算得到所述初始飞行时间τ;Determine the frequency of the compensated beat signal according to the frequency domain signal corresponding to the compensated beat signal, and calculate the initial flight time τ according to the frequency of the compensated beat signal;
    其中,所述补偿后的拍频信号的频率为所述补偿后的拍频信号对应的频域信号的峰值点位置对应的频率,或者所述频域信号的半高宽度中点位置对应的频率。Wherein, the frequency of the compensated beat frequency signal is the frequency corresponding to the peak point position of the frequency domain signal corresponding to the compensated beat frequency signal, or the frequency corresponding to the midpoint position of the half-height width of the frequency domain signal .
  10. 根据权利要求1-9任一项所述的方法,其特征在于,所述根据所述参考路拍频信号计算得到激光器发射信号的调频非线性项ε(t),包括:The method according to any one of claims 1-9, wherein the calculation of the frequency modulation nonlinear term ε(t) of the laser emission signal according to the reference path beat frequency signal comprises:
    对所述参考路拍频信号进行希尔伯特变换,以得到变换后的拍频信号;Performing Hilbert transform on the reference road beat signal to obtain a transformed beat signal;
    根据所述变换后的拍频信号计算得到所述激光器发射信号的调频频率;Calculating the frequency modulation frequency of the laser emission signal according to the converted beat frequency signal;
    根据所述激光器发射信号的调频频率和理想线性调频频率计算得到所述激光器发射信号的调频非线性项ε(t)。According to the frequency modulation frequency of the laser emission signal and the ideal chirp frequency, the frequency modulation nonlinear term ε(t) of the laser emission signal is calculated.
  11. 一种基于调频非线性校正的测距装置,其特征在于,包括:A distance measuring device based on frequency modulation nonlinear correction, which is characterized in that it comprises:
    获取单元,用于获取参考路拍频信号和测量路拍频信号;The acquisition unit is used to acquire the reference road beat frequency signal and the measurement road beat frequency signal;
    计算单元,用于根据所述参考路拍频信号计算得到调频非线性度;并根据所述测量路拍频信号计算得到初始飞行时间τ;The calculation unit is configured to calculate the frequency modulation nonlinearity according to the reference road beat signal; and calculate the initial flight time τ according to the measured road beat signal;
    所述确定单元,用于确定所述调频非线性度不小于第一预设阈值,所述计算单元还用于根据所述参考路拍频信号计算得到激光器发射信号的调频非线性项ε(t),The determining unit is configured to determine that the frequency modulation nonlinearity is not less than a first preset threshold, and the calculation unit is further configured to calculate the frequency modulation nonlinearity term ε(t ),
    校正单元,用于根据所述初始飞行时间τ和所述激光器发射信号的调频非线性项ε(t)对所述测量路拍频信号进行非线性迭代校正,以得到目标飞行时间;A correction unit, configured to perform nonlinear iterative correction on the measurement path beat frequency signal according to the initial flight time τ and the frequency modulation nonlinear term ε(t) of the laser emission signal to obtain the target flight time;
    所述计算单元,还用于根据所述目标飞行时间计算得到目标距离。The calculation unit is further configured to calculate the target distance according to the target flight time.
  12. 根据权利要求11所述的装置,其特征在于,所述校正单元具体用于:The device according to claim 11, wherein the correction unit is specifically configured to:
    S1:根据所述延时飞行时间τ′ i-1对所述激光器发射信号的调频非线性项ε(t)进行延时,以得到所述测量路回波信号对应的调频非线性项ε(t-τ′ i-1),所述延时飞行时间τ′ i-1是基于所述飞行时间τ i-1得到的; S1: Delay the FM nonlinear term ε(t) of the laser emission signal according to the delay flight time τ′i -1 to obtain the FM nonlinear term ε(t) corresponding to the echo signal of the measurement path t-τ′ i-1 ), the delayed flight time τ′ i-1 is obtained based on the flight time τ i-1 ;
    S2:根据所述测量路回波信号对应的调频非线性项ε(t-τ′ i-1)对所述测量路拍频信号进行调频非线性校正,以得到校正拍频信号; S2: Perform FM nonlinear correction on the measurement path beat signal according to the FM nonlinear term ε(t-τ′ i-1) corresponding to the measurement path echo signal to obtain a corrected beat signal;
    S3、根据所述校正拍频信号确定所述目标飞行时间。S3. Determine the target flight time according to the corrected beat frequency signal.
  13. 根据权利要求12所述的装置,其特征在于,在所述根据所述校正拍频信号确定所述目标飞行时间的方面,所述校正单元具体用于:The device according to claim 12, wherein, in the aspect of determining the target flight time according to the corrected beat signal, the correction unit is specifically configured to:
    根据所述校正拍频信号获取飞行时间τ i Obtain the flight time τ i according to the corrected beat frequency signal;
    当所述频谱峰值peak i满足预设条件时,将所述飞行时间τ i确定为所述目标飞行时间;当所述频谱峰值peak i不满足所述预设条件时,令i=i+1,并重复执行步骤S1-S3;所述频谱 峰值peak i为所述校正拍频信号对应的频域信号的峰值; When the spectral peak i meets a preset condition, the flight time τ i is determined as the target flight time; when the spectral peak i does not meet the preset condition, let i=i+1 , And repeat steps S1-S3; the spectrum peak peak i is the peak value of the frequency domain signal corresponding to the corrected beat signal;
    当所述i=1时,所述peak i-1为所述初始频谱峰值,所述飞行时间τ i-1为所述初始飞行时间τ。 When i=1, the peak i-1 is the initial frequency spectrum peak, and the flight time τ i-1 is the initial flight time τ.
  14. 根据权利要求13所述的装置,其特征在于,所述预设条件为(peak i-peak i-1)/peak i-1<Thr rate,所述Thr rate为第三预设阈值。 The apparatus according to claim 13, wherein said predetermined condition is a (peak i -peak i-1) / peak i-1 <Thr rate, the Thr rate is a third predetermined threshold value.
  15. 根据权利要求11所述的装置,其特征在于,所述校正单元具体用于:The device according to claim 11, wherein the correction unit is specifically configured to:
    根据所述初始飞行时间τ和所述激光器发射信号的调频非线性项ε(t)对所述测量路拍频信号进行M次非线性迭代计算,以得到所述目标飞行时间,所述M为大于1的整数,According to the initial flight time τ and the frequency modulation nonlinear term ε(t) of the laser emission signal, perform M non-linear iterative calculations on the measurement path beat frequency signal to obtain the target flight time, where M is An integer greater than 1,
    其中,在进行第i次非线性迭代计算时,根据所述延时飞行时间τ′ i-1对所述激光器发射信号的调频非线性项ε(t)进行延时,以得到测量路回波信号对应的调频非线性项ε(t-τ′ i-1);所述延时飞行时间τ′ i-1是基于所述飞行时间τ i-1得到的; Wherein, when the i-th nonlinear iterative calculation is performed, the frequency modulation nonlinear term ε(t) of the laser emission signal is delayed according to the delay flight time τ′ i-1 to obtain the measurement path echo The FM nonlinear term ε(t-τ′ i-1 ) corresponding to the signal; the delayed flight time τ′ i-1 is obtained based on the flight time τ i-1 ;
    根据所述测量路回波信号对应的调频非线性项ε(t-τ′ i-1)对所述测量路拍频信号进行调频非线性校正,以得到校正拍频信号;根据所述校正拍频信号获取飞行时间τ iPerform FM nonlinear correction on the measurement path beat signal according to the FM nonlinear term ε(t-τ′ i-1 ) corresponding to the echo signal of the measurement path to obtain a corrected beat signal; Frequency signal to obtain flight time τ i ;
    其中,当i=1时,所述飞行时间τ i-1为所述初始飞行时间τ;当所述i=M时,所述目标飞行时间为飞行时间τ iWherein, when i=1, the flight time τ i-1 is the initial flight time τ; when i=M, the target flight time is the flight time τ i .
  16. 根据权利要求13-15任一项所述的装置,其特征在于,在所述根据所述校正拍频信号获取飞行时间τ i的方面,所述校正单元具体用于: The device according to any one of claims 13-15, wherein, in the aspect of acquiring the flight time τ i according to the corrected beat frequency signal, the correction unit is specifically configured to:
    对所述校正拍频信号进行快速傅里叶变换FFT,以得到所述校正拍频信号对应的频域信号;Performing fast Fourier transform FFT on the corrected beat frequency signal to obtain a frequency domain signal corresponding to the corrected beat frequency signal;
    根据所述校正拍频信号对应的频域信号获取飞行时间τ i Obtain the flight time τ i according to the frequency domain signal corresponding to the corrected beat frequency signal.
  17. 根据权利要求12-16任一项所述的装置,其特征在于,在所述校正拍频信号包括第一校正拍频信号或第二校正拍频信号;所述根据所述测量路回波信号对应的调频非线性项ε(t-τ′ i-1)对所述测量路拍频信号进行调频非线性校正,以得到校正拍频信号的方面,所述校正单元具体用于: The device according to any one of claims 12-16, wherein the corrected beat frequency signal comprises a first corrected beat frequency signal or a second corrected beat frequency signal; the echo signal according to the measurement path The corresponding FM nonlinear term ε(t-τ′ i-1 ) performs FM nonlinear correction on the measurement path beat signal to obtain the aspect of correcting the beat signal, and the correction unit is specifically configured to:
    当所述测量路拍频信号的信噪比不低于第二预设阈值时,根据所述激光器发射信号的调频非线性项ε(t)和所述测量路回波信号对应的调频非线性项ε(t-τ′ i-1)对所述测量路拍频信号进行调频非线性项校正,以得到所述第一校正拍频信号; When the signal-to-noise ratio of the measurement path beat signal is not lower than the second preset threshold, the frequency modulation nonlinearity term ε(t) of the laser emission signal and the frequency modulation nonlinearity corresponding to the measurement path echo signal The term ε(t-τ′ i-1 ) performs frequency modulation nonlinear term correction on the measurement path beat signal to obtain the first corrected beat signal;
    当所述测量路拍频信号的信噪比低于所述第二预设阈值时,根据所述测量路回波信号对应的调频非线性项ε(t-τ′ i-1)对所述补偿后的拍频信号进行接收端调频非线性项校正,以得到所述第二校正拍频信号; When the signal-to-noise ratio of the measurement path beat frequency signal is lower than the second preset threshold value, the frequency modulation nonlinearity term ε(t-τ′ i-1 ) corresponding to the measurement path echo signal is compared to the The compensated beat frequency signal is corrected by a frequency modulation nonlinear term at the receiving end to obtain the second corrected beat frequency signal;
    其中,所述补偿后的拍频信号是基于根据所述激光器发射信号的调频非线性项ε(t)对所述测量路拍频信号进行发射端的调频非线性项补偿得到的。Wherein, the compensated beat frequency signal is obtained based on the frequency modulation nonlinearity term ε(t) of the measurement path beat frequency signal at the transmitting end according to the frequency modulation nonlinearity term ε(t) of the laser transmitting signal.
  18. 根据权利要求17所述的装置,其特征在于,在所述根据所述测量路拍频信号计算 得到初始飞行时间τ的方面,所述计算单元具体用于:The device according to claim 17, wherein, in terms of calculating the initial flight time τ according to the measured road beat frequency signal, the calculating unit is specifically configured to:
    当所述测量路拍频信号的信噪比不低于第二预设阈值时,对所述测量路拍频信号进行FFT,以得到所述测量路拍频信号对应的频域信号;When the signal-to-noise ratio of the measurement road beat signal is not lower than a second preset threshold, perform an FFT on the measurement road beat signal to obtain a frequency domain signal corresponding to the measurement road beat signal;
    根据所述测量路拍频信号对应的频域信号确定所述测量路拍频信号的频率,并根据所述测量路拍频信号的频率计算得到所述初始飞行时间τ;Determine the frequency of the measurement road beat signal according to the frequency domain signal corresponding to the measurement road beat signal, and calculate the initial flight time τ according to the frequency of the measurement road beat signal;
    其中,所述测量路拍频信号的频率为所述测量路拍频信号对应的频域信号的峰值点位置对应的频率,或者所述频域信号的半高宽度中点位置对应的频率。Wherein, the frequency of the measurement path beat signal is a frequency corresponding to a peak point position of a frequency domain signal corresponding to the measurement path beat signal, or a frequency corresponding to a midpoint position of a half-height width of the frequency domain signal.
  19. 根据权利要求17所述的装置,其特征在于,在所述根据所述测量路拍频信号计算得到初始飞行时间τ的方面,所述计算单元具体用于:The device according to claim 17, wherein in the aspect of calculating the initial flight time τ according to the measured road beat frequency signal, the calculating unit is specifically configured to:
    当所述测量路拍频信号的信噪比低于第二预设阈值时;根据所述激光器发射信号的调频非线性项ε(t)对所述测量路拍频信号进行发射端的调频非线性项补偿,以得到补偿后的拍频信号;激光器发射信号的调频非线性项是基于所述参考路拍频信号计算得到的;When the signal-to-noise ratio of the measurement path beat frequency signal is lower than the second preset threshold; perform the frequency modulation nonlinearity at the transmitting end of the measurement path beat frequency signal according to the frequency modulation nonlinearity term ε(t) of the laser emission signal Term compensation to obtain the compensated beat frequency signal; the frequency modulation nonlinear term of the laser emission signal is calculated based on the reference path beat frequency signal;
    对所述补偿后的拍频信号进行FFT,以得到所述补偿后的拍频信号对应的频域信号;Performing FFT on the compensated beat signal to obtain a frequency domain signal corresponding to the compensated beat signal;
    根据所述补偿后的拍频信号对应的频域信号确定所述补偿后的拍频信号的频率,并根据所述补偿后的拍频信号的频率计算得到所述初始飞行时间τ;Determine the frequency of the compensated beat signal according to the frequency domain signal corresponding to the compensated beat signal, and calculate the initial flight time τ according to the frequency of the compensated beat signal;
    其中,所述补偿后的拍频信号的频率为所述补偿后的拍频信号对应的频域信号的峰值点位置对应的频率,或者所述频域信号的半高宽度中点位置对应的频率。Wherein, the frequency of the compensated beat frequency signal is the frequency corresponding to the peak point position of the frequency domain signal corresponding to the compensated beat frequency signal, or the frequency corresponding to the midpoint position of the half-height width of the frequency domain signal .
  20. 根据权利要求11-19任一项所述的装置,其特征在于,在根据所述参考路拍频信号计算得到激光器发射信号的调频非线性项ε(t)的方面,所述计算单元具体用于:The device according to any one of claims 11-19, wherein, in terms of calculating the frequency modulation nonlinear term ε(t) of the laser emission signal according to the reference path beat frequency signal, the calculation unit specifically uses At:
    对所述参考路拍频信号进行希尔伯特变换,以得到变换后的拍频信号;Performing Hilbert transform on the reference road beat signal to obtain a transformed beat signal;
    根据所述变换后的拍频信号计算得到所述激光器发射信号的调频频率;Calculating the frequency modulation frequency of the laser emission signal according to the converted beat frequency signal;
    根据所述激光器发射信号的调频频率和理想线性调频频率计算得到所述激光器发射信号的调频非线性项ε(t)。According to the frequency modulation frequency of the laser emission signal and the ideal chirp frequency, the frequency modulation nonlinear term ε(t) of the laser emission signal is calculated.
  21. 一种基于调频非线性校正的测距装置,其特征在于,包括:A distance measuring device based on frequency modulation nonlinear correction, which is characterized in that it comprises:
    存储器,用于存储指令;以及Memory for storing instructions; and
    至少一台处理器,与所述存储器耦合;At least one processor coupled with the memory;
    其中,当所至少一台处理器执行所述指令时,所述指令致使所述处理器执行如权利要求1-10任一项所述的方法。Wherein, when at least one processor executes the instruction, the instruction causes the processor to execute the method according to any one of claims 1-10.
  22. 一种计算机存储介质,其特征在于,所述计算机存储介质存储有计算机程序,所述计算机程序包括程序指令,所述程序指令当被处理器执行时使所述处理器执行如权利要求1-10任一项所述的方法。A computer storage medium, wherein the computer storage medium stores a computer program, the computer program includes program instructions, and when executed by a processor, the program instructions cause the processor to execute as claimed in claims 1-10. Any of the methods.
  23. 一种雷达,其特征在于,包括:A radar, characterized in that it includes:
    激光器、以及与所述激光器耦合的处理器;A laser, and a processor coupled with the laser;
    所述处理器用于执行如权利要求1至10中任一项所述的方法。The processor is configured to execute the method according to any one of claims 1 to 10.
PCT/CN2020/074495 2020-02-07 2020-02-07 Frequency modulation nonlinear correction-based ranging method and related device WO2021155567A1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
PCT/CN2020/074495 WO2021155567A1 (en) 2020-02-07 2020-02-07 Frequency modulation nonlinear correction-based ranging method and related device
CN202080004674.2A CN112639528B (en) 2020-02-07 2020-02-07 Ranging method based on frequency modulation nonlinear correction and related device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2020/074495 WO2021155567A1 (en) 2020-02-07 2020-02-07 Frequency modulation nonlinear correction-based ranging method and related device

Publications (1)

Publication Number Publication Date
WO2021155567A1 true WO2021155567A1 (en) 2021-08-12

Family

ID=75291197

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2020/074495 WO2021155567A1 (en) 2020-02-07 2020-02-07 Frequency modulation nonlinear correction-based ranging method and related device

Country Status (2)

Country Link
CN (1) CN112639528B (en)
WO (1) WO2021155567A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023151706A1 (en) * 2022-02-14 2023-08-17 北京万集科技股份有限公司 Photoelectric phase-locked loop linear correction system, and ranging apparatus

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110096630A1 (en) * 2009-10-28 2011-04-28 Agency For Defense Development Method for estimating target range error and sonar system thereof
CN103163513A (en) * 2013-03-13 2013-06-19 哈尔滨工业大学 Frequency modulated continuous wave (FMCW) laser radar high-accuracy signal measurement method based on phase demodulation method
CN108663684A (en) * 2018-06-08 2018-10-16 天津大学 A kind of phase difference ranging method based on equal optical frequency intervals resampling
CN109521436A (en) * 2018-10-16 2019-03-26 天津大学 A kind of moving object dynamic distance measurement method based on double light path CW with frequency modulation
CN109633670A (en) * 2018-10-25 2019-04-16 上海无线电设备研究所 It is a kind of to utilize the laser pulse ranging method for receiving signal width amendment measurement error
CN109917368A (en) * 2019-04-18 2019-06-21 重庆大学 A method of realizing active refelction type micro-wave radar absolute distance measurement

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5289252A (en) * 1992-12-08 1994-02-22 Hughes Aircraft Company Linear frequency modulation control for FM laser radar
DE10018553A1 (en) * 2000-04-14 2001-10-18 Bosch Gmbh Robert Detecting and correcting nonlinearities in microwave radar system e.g. for motor vehicle by updating voltage-frequency characteristic of fm generating component based on reference signal
JP2008079165A (en) * 2006-09-22 2008-04-03 Toyota Industries Corp Fm detection circuit
US8026843B2 (en) * 2008-01-31 2011-09-27 Infineon Technologies Ag Radar methods and systems using ramp sequences
US10191154B2 (en) * 2016-02-11 2019-01-29 Massachusetts Institute Of Technology Methods and apparatus for time-of-flight imaging
CN107907878B (en) * 2017-11-08 2021-05-11 零八一电子集团有限公司 Method for obtaining FMCW radar distance measurement value with high precision
US10969465B2 (en) * 2018-05-01 2021-04-06 Mitsubishi Electric Research Laboratories, Inc. Reference-free nonlinearity correction for FMCW-based sensing systems
CN108761429A (en) * 2018-05-17 2018-11-06 北京航天计量测试技术研究所 A kind of CW with frequency modulation laser radar Nonlinear frequency modulation high-precision real-time correction method
CN109031340B (en) * 2018-07-26 2022-07-08 天津大学 Continuous frequency modulation laser radar device for measuring object movement speed
CN109490857B (en) * 2018-10-29 2023-01-06 北京遥感设备研究所 Method and system for determining frequency modulation nonlinearity of LFM pulse signal of radar equipment
CN110146012B (en) * 2019-06-11 2021-07-02 西安工业大学 Semiconductor laser nonlinear frequency modulation pre-correction method
CN112987025B (en) * 2019-08-23 2023-10-13 长沙天仪空间科技研究院有限公司 Laser imaging device and method based on synthetic aperture

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110096630A1 (en) * 2009-10-28 2011-04-28 Agency For Defense Development Method for estimating target range error and sonar system thereof
CN103163513A (en) * 2013-03-13 2013-06-19 哈尔滨工业大学 Frequency modulated continuous wave (FMCW) laser radar high-accuracy signal measurement method based on phase demodulation method
CN108663684A (en) * 2018-06-08 2018-10-16 天津大学 A kind of phase difference ranging method based on equal optical frequency intervals resampling
CN109521436A (en) * 2018-10-16 2019-03-26 天津大学 A kind of moving object dynamic distance measurement method based on double light path CW with frequency modulation
CN109633670A (en) * 2018-10-25 2019-04-16 上海无线电设备研究所 It is a kind of to utilize the laser pulse ranging method for receiving signal width amendment measurement error
CN109917368A (en) * 2019-04-18 2019-06-21 重庆大学 A method of realizing active refelction type micro-wave radar absolute distance measurement

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023151706A1 (en) * 2022-02-14 2023-08-17 北京万集科技股份有限公司 Photoelectric phase-locked loop linear correction system, and ranging apparatus

Also Published As

Publication number Publication date
CN112639528A (en) 2021-04-09
CN112639528B (en) 2021-12-14

Similar Documents

Publication Publication Date Title
JP6830311B2 (en) Phase calibration of stepped chirp signals for synthetic aperture radar
US11119209B2 (en) Robust dual-radar-beam systems and methods for traffic monitoring
US8344942B2 (en) System and method for resolving ambiguity in radar, lidar, and acoustic systems
CN103323822B (en) Method and device for estimating channel errors
CN103399307B (en) A kind of method of correction channel error and device
CN106802416B (en) Fast factorization back projection SAR self-focusing method
CN109471095A (en) Fmcw radar distance estimating algorithm based on iteratively faster interpolation
US11269046B2 (en) Phase-comparison of multi-frequency transmissions for assisting the determination of position or time
WO2021155567A1 (en) Frequency modulation nonlinear correction-based ranging method and related device
CN108535719A (en) CW with frequency modulation landing radar speed-measuring method based on the correction of Doppler frequency spectrum center of gravity
CN116299551A (en) Terahertz SAR two-dimensional self-focusing imaging algorithm
CN115034993A (en) Method and device for defocusing removal processing of moving target image and storage medium
CN110596659B (en) Single-pulse angle measurement internal calibration source loop self-calibration method based on secondary radar
CN109358330B (en) Bistatic forward-looking SAR linear time synchronization error estimation method
CN108445477A (en) The precision distance measurement method of airdrome scene foreign bodies detection radar
CN112415512B (en) SAR moving target focusing method based on advance and retreat method and golden section method
KR102172085B1 (en) Radar apparatus and method for measuring distance of target using the same
CN110988860B (en) Frequency modulation continuous wave SAR motion compensation method based on triangular wave modulation
WO2021184206A1 (en) Linear swept frequency correction method and device, storage medium, and system
KR102053203B1 (en) Multiple altitude operating interferometer radar altimeter apparatus based on valid angle identification and method thereof
Mei et al. Research on amplitude-phase error for LFMCW radar
Wang et al. Robust TDOA/FDOA estimation from emitter signals for hybrid localization using UAVs
JP3835035B2 (en) Radar cross section measuring device
Samarasekera et al. Iterative minimum-entropy based algorithm for phase noise removal in fmcw radars
Yang et al. Autofocus Algorithm for SAR

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 20917475

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 20917475

Country of ref document: EP

Kind code of ref document: A1