CN109932710B - Long-distance target speed correction method based on sawtooth LFMCW waveform system speed measurement radar - Google Patents

Long-distance target speed correction method based on sawtooth LFMCW waveform system speed measurement radar Download PDF

Info

Publication number
CN109932710B
CN109932710B CN201810898984.8A CN201810898984A CN109932710B CN 109932710 B CN109932710 B CN 109932710B CN 201810898984 A CN201810898984 A CN 201810898984A CN 109932710 B CN109932710 B CN 109932710B
Authority
CN
China
Prior art keywords
target
distance
frequency
speed
lfmcw
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810898984.8A
Other languages
Chinese (zh)
Other versions
CN109932710A (en
Inventor
王勋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Leihua Electronic Technology Research Institute Aviation Industry Corp of China
Original Assignee
Leihua Electronic Technology Research Institute Aviation Industry Corp of China
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 Leihua Electronic Technology Research Institute Aviation Industry Corp of China filed Critical Leihua Electronic Technology Research Institute Aviation Industry Corp of China
Priority to CN201810898984.8A priority Critical patent/CN109932710B/en
Publication of CN109932710A publication Critical patent/CN109932710A/en
Application granted granted Critical
Publication of CN109932710B publication Critical patent/CN109932710B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention relates to a remote target speed correction method based on a sawtooth LFMCW waveform system speed measurement radar, which comprises the following steps: fast time dimension FFT measures the target distance; obtaining Doppler frequency shift and long-distance echo frequency shift by slow time dimension FFT; and correcting the carrier frequency to obtain the speed value of the remote target. The correction method of the invention can fundamentally solve the error generated when the radar slow time dimension FFT is calculated by the linear frequency modulation continuous wave system on the long-distance target and the large frequency modulation slope, has small resource consumption in signal processing, is easy to realize, and can be widely used in the phase-coherent system LFMCW radar.

Description

Long-distance target speed correction method based on sawtooth LFMCW waveform system speed measurement radar
Technical Field
The invention belongs to the technical field of radar signal processing, and particularly relates to a remote target speed correction method based on a sawtooth LFMCW waveform system speed measurement radar.
Background
Linear Frequency Modulated Continuous Wave (LFMCW) radar has high range resolution and range accuracy, and has more obvious advantages in low range blind area than pulse system radar.
The saw-toothed LFMCW waveform is shown in fig. 1 below.
The LFMCW radar transmit signal model may be expressed as:
s 1 (t)=cos(2πf 0 t+πμt 2 ),0≤t≤T e (1)
wherein μ = B/T e Mu is LFM slope, B is bandwidth, T e Frequency modulation period is the sweep frequency period; f. of 0 Is the start frequency (carrier frequency) of the chirp.
Assuming that there is a point scatterer at a distance R, the radar receives its echo signal as:
s r (t)=acos[2πf 0 (t-Δτ)+πμ(t-Δτ) 2 ],0≤t-Δτ≤T e (2)
in the formula, a is a signal amplitude, which is related to a target radar scattering cross section (RCS), a distance and an antenna gain; Δ τ =2R/c is the time delay and c is the speed of light.
The reference signal input by the mixer is:
s ref (t)=cos(2πf r0 +πμt 2 ),0≤t≤T e (3)
in the formula, f r0 For the start frequency of the reference signal LFM, let f be r0 =f 0
The complex signal model after the frequency mixing and the low-pass filtering of the received signal and the reference signal is as follows:
Figure BDA0001758905970000021
the instantaneous frequency of the signal is:
Figure BDA0001758905970000022
the above equation shows that the range of the target is proportional to the instantaneous frequency. So that the received signal is sampled and the sample sequence is FFT-ed at frequency f i The peak position of (a) corresponds to a target distance of
Figure BDA0001758905970000023
In the actual LFMCW system radar, fast time dimension FFT is carried out in a frequency sweep section, and the target distance is calculated by the frequency peak value according to the formula (6). And performing slow time dimension FFT between every two frequency scanning bands, wherein the frequency peak value corresponds to the target speed.
Figure BDA0001758905970000024
However, when the sawtooth LFMCW waveform is used for speed measurement, the remote target causes carrier frequency change, which further affects the calculation result of the remote target, so that the method cannot be ignored, and the detected speed value needs to be corrected to improve the calculation accuracy.
Disclosure of Invention
The invention aims to provide a method for correcting the speed of a long-distance target based on a sawtooth LFMCW waveform system speed measurement radar, which is used for correcting the problem that when the sawtooth LFMCW waveform is used for speed measurement, the long-distance target causes carrier frequency change, so that the calculation result of the long-distance target is influenced.
In order to achieve the purpose, the invention adopts the technical scheme that: a long-distance target speed correction method based on a sawtooth LFMCW waveform system speed measurement radar comprises the following steps:
fast time dimension FFT measures the target distance;
obtaining Doppler frequency shift and long-distance echo frequency shift by slow time dimension FFT;
and correcting the carrier frequency to obtain the speed value of the remote target.
Further, the target distance is:
Figure BDA0001758905970000031
in the formula: f. of i For frequency, c is the speed of light, μ is the LFM slope, μ = B/T e B is bandwidth of frequency modulation, T e Is a frequency modulation period.
Further, the corrected target speed value is
Figure BDA0001758905970000032
In the formula: f. of i Is a frequency f o Is the original carrier frequency, c is the speed of light, μ is the LFM slope, μ = B/T e B is bandwidth of frequency modulation, T e For frequency-modulated period, R 0 Is the target distance.
The correction method of the invention can fundamentally solve the error generated when the radar slow time dimension FFT is calculated by the linear frequency modulation continuous wave system on the long-distance target and the large frequency modulation slope, has small resource consumption in signal processing, is easy to realize, and can be widely used in the phase-coherent system LFMCW radar.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the invention and together with the description, serve to explain the principles of the invention.
Fig. 1 is a diagram illustrating the relationship between the instantaneous frequency of a transmitted signal and time in the prior art.
FIG. 2 is a flow chart of a method for correcting a speed of a remote target according to the present invention.
Fig. 3 is a three-dimensional spectrum diagram of an echo signal obtained according to simulation parameters in an embodiment of the present invention.
Detailed Description
In order to make the implementation objects, technical solutions and advantages of the present invention clearer, the technical solutions in the embodiments of the present invention will be described in more detail below with reference to the accompanying drawings in the embodiments of the present invention.
The actual near and far target echo simulation processes are given below, and the parameters of the simulation processes are described in table 1 below:
table 1 simulation program input parameters
Description of parameters Unit State of state Simulation parameter
Frequency modulation period T e S Input the method 1e-3
Bandwidth B of frequency modulation Hz Input the method 1e8
Target distance R 0 m Input the method [1.1e6,5e3,2.5e4]
Target velocity v m/s Input the method [23,7,33]
Target signal to noise ratio dB Input device [20,15,10]
Carrier frequency f o Hz Input device 1e10
Fast time dimension sampling rate f s Hz Input device 2.5e7
Fast time dimension FFT point number Nfft 2 Input device 4096
FFT point number Nfft in slow time dimension 1 Input device 64
As shown in FIG. 2, the simulation results of the three target range gates can be judged to be [2732,3551,1367] in sequence according to the echo signal-to-noise ratio, and the real distances are [1.1e6,5e3,2.5e4].
Firstly, whether the simulation is correct is verified: unambiguous distance R u Namely:
R u =c*T e /2 (9)
wherein, T e The frequency modulation period, c is the speed of light.
Distance after blur R 0 Namely:
Figure BDA0001758905970000041
wherein f is s For fast time dimension sampling rate, nfft2 is the fast time dimension FFT point number. The three target deblurring back range gates are [2731.8,3551,1366.2 ]]。
Therefore, the simulation result of the range gate is consistent with the calculation result, the simulation result is correct, and the simulation data is available.
The three target speed dimension simulation results are [28,31,14] in sequence, and the real speed is [23,7,33].
Unambiguous velocity V u Namely:
Figure BDA0001758905970000042
wherein f is 0 As carrier frequency, velocity after blurring V 0 Namely:
Figure BDA0001758905970000051
wherein, nfft1 is the number of FFT points in the slow time dimension. Carrier frequency in accordance with f 0 The results of the three targets after removing the non-fuzzy speed are calculated to be 35.13,30.86 and 13.8]. It can be seen that the target is in close range (about kilometer), and the speed is according to the carrier frequency f 0 The calculation result is basically the same as the simulation value, but the far-distance echo speed and the simulation value are greatly different.
The method of the invention gives out the reason and the solution that the difference between the remote echo speed and the simulation value is large in principle. As can be seen from the previous formula, the slow time dimension echo signal is actually linear frequency modulation, the frequency obtained by actual FFT calculation is broadband, and the center frequency is
Figure BDA0001758905970000052
The carrier frequency after correction is f 0 ', i.e.:
Figure BDA0001758905970000053
the carrier frequency is calculated by adopting the formula, and the results of three targets without fuzzy speed are obtained as 27.4274,30.86 and 13.8, which are basically consistent with the simulation result.
The method can correct the error of the LFMCW speed measurement algorithm of the long-distance echo caused by carrier frequency change, and has important application prospect in practical engineering.
The method can solve the error generated in the calculation of the slow time dimension FFT of the radar of the linear frequency modulation continuous wave system on the long-distance target and the large frequency modulation slope in principle, and compared with the prior art, the method has the following advantages that:
(1) The method is easy to realize in engineering, and the signal processing consumes less resources;
(2) The method of the invention solves the error caused by LFMCW waveform speed measurement in principle;
(3) The method can be widely applied to the phase-coherent system LFMCW radar.
The above description is only for the best mode of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the appended claims.

Claims (2)

1. The method for correcting the speed of the long-distance target based on the sawtooth LFMCW waveform system speed measurement radar is characterized by comprising the following steps:
fast time dimension FFT measures the target distance;
obtaining Doppler frequency shift and long-distance echo frequency shift by slow time dimension FFT;
and (3) correcting the carrier frequency to obtain a speed value of the remote target, wherein the relation between the corrected target speed v and the frequency is as follows:
Figure FDA0004027363740000011
in the formula: f. of i Is a frequency f o Is the original carrier frequency, c is the speed of light, μ is the LFM slope, μ = B/T e B is bandwidth of frequency modulation, T e For frequency-modulated period, R 0 Is the target distance.
2. The method for correcting the speed of the long-distance target based on the sawtooth LFMCW waveform system speed radar according to claim 1, wherein the target distance is as follows:
Figure FDA0004027363740000012
CN201810898984.8A 2018-08-08 2018-08-08 Long-distance target speed correction method based on sawtooth LFMCW waveform system speed measurement radar Active CN109932710B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810898984.8A CN109932710B (en) 2018-08-08 2018-08-08 Long-distance target speed correction method based on sawtooth LFMCW waveform system speed measurement radar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810898984.8A CN109932710B (en) 2018-08-08 2018-08-08 Long-distance target speed correction method based on sawtooth LFMCW waveform system speed measurement radar

Publications (2)

Publication Number Publication Date
CN109932710A CN109932710A (en) 2019-06-25
CN109932710B true CN109932710B (en) 2023-03-03

Family

ID=66984467

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810898984.8A Active CN109932710B (en) 2018-08-08 2018-08-08 Long-distance target speed correction method based on sawtooth LFMCW waveform system speed measurement radar

Country Status (1)

Country Link
CN (1) CN109932710B (en)

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7474255B2 (en) * 2006-12-05 2009-01-06 Chung Shan Institute Of Science And Technology, Armaments Bureau, M.N.D. Target tracking method of radar with frequency modulated continuous wave
US9739881B1 (en) * 2016-03-24 2017-08-22 RFNAV, Inc. Low cost 3D radar imaging and 3D association method from low count linear arrays for all weather autonomous vehicle navigation
CN106842166B (en) * 2017-03-31 2019-04-30 东南大学 A kind of solution velocity ambiguity method suitable for LFMCW radar system
CN108279403B (en) * 2018-01-04 2020-03-17 电子科技大学 Keystone transformation parallel implementation method based on distance framing
CN108287335A (en) * 2018-04-19 2018-07-17 广州合智瑞达科技有限公司 A method of ranging and range rate being carried out to multiple target using the frequency modulated signal of LFMCW radars

Also Published As

Publication number Publication date
CN109932710A (en) 2019-06-25

Similar Documents

Publication Publication Date Title
CN107966688B (en) Broadband radar target speed ambiguity resolving method based on phase interference technology
CN108415010B (en) Radar multi-target detection method based on trapezoidal LFMCW modulation
CN107861117B (en) Multi-target parameter measuring method suitable for continuous wave perimeter surveillance radar
CN107561508B (en) Coherent accumulation detection method for uniformly accelerated moving target
CN109521417B (en) Multi-target detection calculation method based on FMCW radar waveform
WO2018163677A1 (en) Radar device
CN108535719B (en) Frequency modulation continuous wave landing radar speed measurement method based on Doppler frequency spectrum gravity center correction
CN110632587A (en) Weak moving object monitoring method based on rapid FMCW radar
CN104237866A (en) Improved constant false alarm rate detecting method for ship-borne linear frequency modulation continuous wave radar
CN103454637A (en) Terahertz inverse synthetic aperture radar imaging method based on frequency modulation step frequency
CN108828602B (en) Signal processing method for eliminating velocity ambiguity in pulse phase dry method velocity measurement
CN111580063B (en) Radar target detection method based on generalized solution frequency modulation-wedge transform
CN116087908B (en) Radar high-precision level meter measuring method based on cooperative operation
CN115508820A (en) Target detection method of linear frequency modulation pulse radar
JP5699405B2 (en) Radar received signal processing apparatus and method
Molchanov et al. On micro-Doppler period estimation
CN109932710B (en) Long-distance target speed correction method based on sawtooth LFMCW waveform system speed measurement radar
CN110441749A (en) Frequency stepping radar target motion parameter estimation method
JP2002131421A (en) Radar signal processing device and method thereof
CN113835077B (en) Search radar target detection method and system based on variable pulse repetition frequency
Willetts et al. Optimal time-frequency distribution selection for LPI radar pulse classification
KR101185480B1 (en) Improvement method of the accuracy of beat frequency in fmcw radar for detecting and tracking a high speed moving target
CN107678002B (en) Method for efficiently measuring ground clutter improvement factor of MTD radar
Cao et al. High resolution range profile imaging of high speed moving targets based on fractional Fourier transform
JP2009109457A (en) Radar device

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant