CN113885003A - Target detection method based on non-coherent accumulation mode of coherent radar - Google Patents

Target detection method based on non-coherent accumulation mode of coherent radar Download PDF

Info

Publication number
CN113885003A
CN113885003A CN202111139490.XA CN202111139490A CN113885003A CN 113885003 A CN113885003 A CN 113885003A CN 202111139490 A CN202111139490 A CN 202111139490A CN 113885003 A CN113885003 A CN 113885003A
Authority
CN
China
Prior art keywords
pulse
data
coherent
coherent accumulation
distance
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.)
Pending
Application number
CN202111139490.XA
Other languages
Chinese (zh)
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.)
Shanghai Radio Equipment Research Institute
Original Assignee
Shanghai Radio Equipment Research Institute
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 Shanghai Radio Equipment Research Institute filed Critical Shanghai Radio Equipment Research Institute
Priority to CN202111139490.XA priority Critical patent/CN113885003A/en
Publication of CN113885003A publication Critical patent/CN113885003A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/415Identification of targets based on measurements of movement associated with the target

Landscapes

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

Abstract

The invention discloses a target detection method based on a non-coherent accumulation mode of coherent radar, which comprises the following steps: s1, sampling to obtain a frame of two-dimensional echo data; s2, performing pulse compression on each data of the S1; s3, carrying out distance walking calibration on the data of the S2; s4, carrying out coherent accumulation on the data of S3; s5, repeating S1-S4 to obtain multi-frame coherent accumulation data; s6, performing interframe alignment on the multiframe coherent accumulation data obtained in the S5; s7, performing non-coherent accumulation on the multi-frame coherent accumulation data after the inter-frame alignment; and S8, carrying out target detection based on the non-coherent accumulated data of S7. The advantages are that: the method realizes distance walking calibration by adopting an autocorrelation method between pulses, and realizes the detection of a moving target accumulated by long-time phase coherence; according to the method, after peak value target alignment is carried out on a two-dimensional data plane between adjacent multiframes, non-coherent accumulation is carried out, and then target detection is carried out, so that the purpose of improving the signal-to-noise ratio is achieved.

Description

Target detection method based on non-coherent accumulation mode of coherent radar
Technical Field
The invention relates to the technical field of target detection, in particular to a target detection method based on a non-coherent accumulation mode of a coherent radar.
Background
Space-based warning radar is an important application field for the development of space-based radar. The space-based warning radar has an important significance in realizing a moving target indication function, and a target detection method of the space-based warning radar becomes a hotspot of radar signal processing research. The accumulation of weak target energy is the main research direction for improving the signal to noise ratio of the echo and realizing the detection of the weak target. The method is limited by the design requirements of lightweight and low power consumption of the space-based radar, the transmitting power, the antenna gain and the like of the space-based radar are limited, and when a target is alarmed, capacity accumulation is often realized through long-time coherent accumulation, so that a target for remote detection is realized. However, as an alarm radar, the priori knowledge is often less, and the priori information such as the speed of the target cannot be obtained. On the one hand, under the condition of long-time coherent accumulation, for a moving object with unknown moving speed, the moving object has the problems of distance walking, phase walking and the like within an accumulation period. On the other hand, under the design principle of lightweight and low power consumption of space-based radar, the total resource of a radar system is limited, and when detection is performed, the detection signal-to-noise ratio is not high, so that the detection signal-to-noise ratio needs to be further improved through an information processing method, and the target detection probability of the radar system is improved.
Aiming at the first problem, namely the problem of distance and phase walk of a moving target with unknown moving speed under the condition of long-time coherent accumulation, the current commonly used algorithm comprises a keystone algorithm and the like. However, the algorithms are huge in operation amount and limited in application to real-time processing platforms of the space-based radar. Meanwhile, for space-based warning radars, if the motion speed is unknown and the problem of Doppler frequency folding exists, the difficulty is increased for distance walking compensation and phase supplement. Under the condition that the motion information of the target is unknown, if the speed information of the target is to be obtained by an accumulation method, the general method is to set a large speed search range and a certain speed search step length, then carry out envelope cyclic shift on pulse echoes according to the calculated number of echo envelope walking caused by each search speed, then carry out non-coherent accumulation on the pulse echoes subjected to the envelope cyclic shift, and finally compare the accumulated results of each possible search speed, wherein the speed corresponding to the maximum accumulation result is the speed estimation value of the target. This has the consequence of being at the expense of a high number of search operations.
On the other hand, the signal-to-noise ratio of the detection needs to be further improved by an information processing method. For coherent radar, pulse compression and coherent accumulation are conventional signal processing methods for improving the signal-to-noise ratio. How to maximize the utilization of echo data among multiple frames and further improve the target detection performance is a technical problem which can be further researched.
Disclosure of Invention
The invention aims to provide a target detection method based on a non-coherent accumulation mode of coherent radar, which is mainly used for solving the problem of detection of a moving target and a small target under a space-based alarm radar; meanwhile, the method can utilize the distance walking information to carry out preliminary estimation of the target speed; performing coherent accumulation on the calibrated data to obtain an accumulated two-dimensional data plane; after peak value target alignment is carried out on two-dimensional data planes between adjacent multiframes, non-coherent accumulation is carried out, and then target detection is carried out, so that the purpose of improving the signal-to-noise ratio is achieved, and the target detection probability of the warning radar is improved.
In order to achieve the purpose, the invention is realized by the following technical scheme:
a target detection method based on a non-coherent accumulation mode of coherent radar comprises the following steps:
s1, sampling to obtain a frame of two-dimensional echo data;
s2, performing pulse compression on the sampling data of each pulse of the two-dimensional echo data in the step S1;
s3, carrying out distance walking calibration on the data after pulse compression;
s4, carrying out coherent accumulation on the data after the distance walk calibration;
s5, repeating the steps S1-S4 to obtain multi-frame coherent accumulation data;
s6, performing interframe alignment on the multiframe coherent accumulation data obtained in the step S5;
s7, performing non-coherent accumulation on the multi-frame coherent accumulation data after the inter-frame alignment;
s8, target detection is performed based on the non-coherent accumulated data of the step S7.
Optionally, in step S1, the two-dimensional echo data is in the form of an LFM signal, the number of pulses of each frame of the LFM signal is N, each pulse samples M points, and the pulse repetition period is TrPulse width of TPEach sub-pulse has a frequency modulation slope of k and a sub-pulse bandwidth of B, and k is B/TP(ii) a Setting the initial radial distance of the uniform motion target as R0Radial velocity of motion v0The target distance of the ith pulse is delayed by taui=2(R0+v0iTr) C, then the target baseband echo signal sr(n, t) is:
Figure BDA0003282628490000031
wherein:
Figure BDA0003282628490000032
in the above formula, c represents the speed of light, and the length of t is M;
simultaneously, noise Noi (n, t) with a certain signal-to-noise ratio is superposed on a target baseband echo signal sr(n, t) obtaining two-dimensional echo data SNr(n,t):
SNr(n,t)=Noi(n,t)+sr(n,t) (3)。
Optionally, in step S2, the matching function of the pulse compression is:
Figure BDA0003282628490000033
the pulse compression is carried out in a frequency domain, and the two-dimensional echo data of the ith pulse is set as SNri(n, t), then the pulse compression algorithm is:
sci(n,t)=IFFT(FFT(SNri(n,t))·conj(FFT(sref(t)))) (5)
wherein, sci(N, t) is pulse pressure data of the ith pulse, i is 1,2, …, N.
Optionally, in step S3, the distance walking calibration is performed on the data after pulse compression by using an inter-pulse cross-correlation method.
Optionally, in step S3, the pulse pressure data sc of the 1 st pulse is used1(t) data sc obtained by compressing the pulse of the ith pulse with reference toi(t) pulse pressure data s of the 1 st pulsecm1(t) performing a cross-correlation process, a cross-correlation function R1i(t) is:
Figure BDA0003282628490000034
find R1i(t) position P where peak value is locatedXi,PXiI.e. represents sci(t) relative to sc1(t) the number of distance units, and pulse pressure data sc of the ith pulsei(t) Cyclic Displacement PXiPoint, obtain scxi(t), implementing distance walk calibration; distance walk calibration is carried out on the data after pulse compression of each pulse, distance walk calibration of one frame of two-dimensional echo data is realized, and a new two-dimensional data plane S is formedxc(n,t)。
sciCross-correlation function R of last pulse in (n, t) two-dimensional array1N(t) position P where peak value is locatedXNThe value is expressed in the Nth pulse and the 1 st pulseNumber of distance units P of distance traveled generated between pulsesXNThen, the corresponding walking distance dR is:
dR=C·PXN·ts (7)
the estimated target speed veComprises the following steps:
Figure BDA0003282628490000041
optionally, in step S4, a new two-dimensional data plane S is processedxc(n, t) performing coherent accumulation, and obtaining the column vector L corresponding to each distance unitxcm(N) coherent accumulation, where M is 1,2, …, M, N is 1,2, …, N, i.e. the column vector Lxcm(n) performing FFT to obtain Fxcm(n):
Fxcm(n)=FFT(Lxcm(n)) (9)
For each column vector Lxcm(n) all carry out coherent accumulation to form a new two-dimensional data plane Fxc1(n,t)。
Compared with the prior art, the invention has the following advantages:
the invention relates to a target detection method based on a non-coherent accumulation mode of a coherent radar, which is characterized in that under the condition of unknown target speed, distance walking calibration is realized by adopting an autocorrelation method based on pulses, and the detection of a moving target accumulated by long-time coherent is realized; meanwhile, the method can utilize the distance walking information to carry out preliminary estimation of the target speed; performing coherent accumulation on the calibrated data to obtain an accumulated two-dimensional data plane; after peak value target alignment is carried out on two-dimensional data planes between adjacent multiframes, non-coherent accumulation is carried out, and then target detection is carried out, so that the purpose of improving the signal-to-noise ratio is achieved, and the target detection probability of the warning radar is improved.
Furthermore, the method can realize the distance walk calibration of the moving target under the condition of long accumulation time without prior information such as the speed of the target and resolving the speed of the target, and can further meet the requirement of rapid detection of the alarm radar.
Furthermore, compared with other distance calibration algorithms, the distance walk calibration method provided by the method is small in calculation amount and more suitable for the signal processing capacity of the space-based platform.
Furthermore, the method can realize the preliminary estimation of the target speed by utilizing a single-frame signal, has short required time and can meet the requirement of the alarm radar for quick detection.
Furthermore, after coherent accumulation, non-coherent accumulation is carried out, so that the signal-to-noise ratio can be further improved, and the radar detection probability is improved.
Drawings
FIG. 1 is a schematic diagram of a target detection method based on a non-coherent accumulation mode of coherent radar according to the present invention;
FIG. 2(a) is a time domain simulation of an LFM echo signal;
FIG. 2(b) is a time domain diagram of an LFM echo signal superimposed with noise;
FIG. 3(a) is a two-dimensional data plane after pulse pressing a single pulse, i.e., the range dimension, in the echo signal of FIG. 2 (a);
FIG. 3(b) is a two-dimensional data plan view after a distance walk calibration of FIG. 3 (a);
FIG. 4(a) is a two-dimensional data plane obtained by direct coherent integration of FIG. 3 (a);
FIG. 4(b) is a two-dimensional data plane obtained by direct coherent integration of FIG. 3 (b);
FIG. 5(a) is data in a velocity dimension corresponding to a peak in the two-dimensional data plane of FIG. 4 (a);
FIG. 5(b) is data in a velocity dimension corresponding to a peak in the two-dimensional data plane of FIG. 4 (b);
FIG. 6(a) is a schematic diagram of a two-dimensional data plane after coherent accumulation of three consecutive frames of signals;
FIG. 6(b) is the distance dimension data corresponding to the peak of different frames;
FIG. 7(a) is a schematic diagram of a two-dimensional data plane after non-coherent accumulation of three consecutive frames of signals;
fig. 7(b) is a schematic diagram of distance dimension data corresponding to the peak of fig. 7 (a).
Detailed Description
The present invention will now be further described by way of the following detailed description of a preferred embodiment thereof, taken in conjunction with the accompanying drawings.
As shown in fig. 1, a method for detecting a target based on a coherent radar in a non-coherent accumulation mode according to the present invention includes:
and S1, sampling to obtain a frame of two-dimensional echo data.
In this embodiment, the signal form of the two-dimensional echo data is an LFM signal (linear frequency hopping signal), where N is the number of pulses of each frame of the LFM signal, M points are sampled per pulse, and the sampling interval is tsWith a pulse repetition period of TrPulse width of TPEach sub-pulse has a frequency modulation slope of k and a sub-pulse bandwidth of B, and k is B/TP(ii) a Setting the initial radial distance of the uniform motion target as R0Radial velocity of motion v0The target distance of the ith pulse is delayed by taui=2(R0+v0iTr) C, then the target baseband echo signal sr(n, t) is:
Figure BDA0003282628490000061
wherein:
Figure BDA0003282628490000062
in the above formula, c represents the speed of light, and t has a length of M.
Simultaneously, noise Noi (n, t) with a certain signal-to-noise ratio is superposed on a target baseband echo signal sr(n, t) obtaining two-dimensional echo data SNr(n,t):
SNr(n,t)=Noi(n,t)+sr(n,t) (3)。
In this embodiment, the target detection method based on the non-coherent accumulation mode of the coherent radar in fig. 1 is used to simulate the echo accumulated by long-time coherent accumulation, and the echo data is a two-dimensional array, where the X dimension represents the pulse accumulation number (slow time dimension) and the Y dimension represents the distance dimension (fast time dimension). Eye settingThe distance of the echo signal is R060km, 800m/s speed, 256 pulse accumulation, 1kHz repetition frequency, 200us sub-pulse width, 40MHz clock, t sampling periods25 ns. The number of sampling points in the range gate is 2000, the power ratio of the target echo signal to the noise is-40 dB, the echo time domain waveform without the noise superimposed is shown in fig. 2(a), the echo time domain waveform with the noise superimposed is shown in fig. 2(b), and as can be seen from fig. 2(a) and 2(b), the signal is completely submerged by the noise.
S2, the sampled data of each pulse of the two-dimensional echo data of step S1 is pulse compressed.
In step S2, the matching function of pulse compression is:
Figure BDA0003282628490000063
the pulse compression is carried out in a frequency domain, and the two-dimensional echo data of the ith pulse is set as SNri(n, t), then the pulse compression algorithm is:
sci(n,t)=IFFT(FFT(SNri(n,t))·conj(FFT(sref(t)))) (5)
wherein, sci(N, t) is pulse pressure data of the ith pulse, i is 1,2, …, N.
In this embodiment, the two-dimensional data plane after pulse compression is shown in fig. 3(a), which is a two-dimensional data plane after pulse compression of a single pulse, i.e., a distance dimension, in the echo data of fig. 2 (a). It can be seen that the pulse pressure peak of different pulses has a distance walk condition due to the existing speed of the target.
And S3, performing distance walking calibration on the data after pulse compression.
In this embodiment, the inter-pulse cross correlation method is used to perform the range walk calibration on the data after the pulse compression.
In the step S3, the pulse pressure data sc of the 1 st pulse is used in the step S31(t) data sc obtained by compressing the pulse of the ith pulse with reference toi(t) pulse pressure number of the 1 st pulseAccording to scm1(t) performing a cross-correlation process, a cross-correlation function R1i(t) is:
Figure BDA0003282628490000071
find R1i(t) position P where peak value is locatedXi,PXiI.e. represents sci(t) relative to sc1(t) the number of distance units, and pulse pressure data sc of the ith pulsei(t) Cyclic Displacement PXiPoint, obtain scxi(t), implementing distance walk calibration; distance walk calibration is carried out on the data after pulse compression of each pulse, distance walk calibration of one frame of two-dimensional echo data is realized, and a new two-dimensional data plane S is formedxc(n,t)。
The distance walk calibration was performed on fig. 3(a) using inter-pulse cross-correlation, and the two-dimensional data plane after calibration is shown in fig. 3 (b). As shown in fig. 3(b), a new two-dimensional data plane S is obtained after distance walk calibration by inter-pulse cross-correlationxc(n, t). It can be seen that the range walk has been significantly corrected, with the range peaks of the different pulses all being at the same location.
sciCross-correlation function R of last pulse in (n, t) two-dimensional array1N(t) position P where peak value is locatedXNThe value represents the number of distance units P traveled between the Nth pulse and the 1 st pulseXNThen, the corresponding walking distance dR is:
dR=C·PXN·ts (7)
the estimated target speed veComprises the following steps:
Figure BDA0003282628490000072
from the above, the motion velocity of the target can be preliminarily estimated through the formula (8), and the target velocity in one frame signal can be quickly estimated.
And S4, performing coherent accumulation on the data after the distance walk calibration.
In the step S4, a new two-dimensional data plane S is processedxc(n, t) performing coherent accumulation, and obtaining the column vector L corresponding to each distance unitxcm(n) coherent accumulation, where M is 1,2, …, M (M corresponds to the column vector L)xcmSubscript of (N), N is 1,2, …, N, i.e., a pair column vector Lxcm(n) performing FFT to obtain Fxcm(n):
Fxcm(n)=FFT(Lxcm(n)) (9)
For each column vector Lxcm(n) all carry out coherent accumulation to form a new two-dimensional data plane Fxc1(n, t) (see FIG. 4 (b)).
Fig. 4(a) is a two-dimensional data plane obtained by directly performing coherent integration on fig. 3(a), and fig. 4(b) is a two-dimensional data plane obtained by performing coherent integration on data subjected to distance walk calibration on fig. 3 (b). The comparison shows that the two-dimensional plane data subjected to coherent accumulation after calibration has higher signal peak value and more concentrated energy.
In this embodiment, fig. 5(a) is data of a velocity dimension corresponding to a peak value of the two-dimensional data plane of fig. 4(a), and fig. 5(b) is data of a velocity dimension corresponding to a peak value of the two-dimensional data plane of fig. 4(b), that is, data of a velocity dimension corresponding to a peak value after coherent accumulation, as can be seen by comparison, peak energy after distance walk calibration is more concentrated, and the signal-to-noise ratio of fig. 5(b) is improved by 10 dB.
And S5, repeating the steps S1-S4 to obtain multi-frame coherent accumulation data.
In this embodiment, taking three frames as an example, that is, the number of frames K is 3, the data planes after three coherent accumulation are respectively Fxc1(n,t)、Fxc2(n,t)、Fxc3(n,t)。
Fxc1(n,t)、Fxc2(n,t)、Fxc3The peak positions of (n, t) are moved from frame to frame in the distance dimension, as shown in fig. 6(a) and 6 (b).
Fig. 6(a) and 6(b) show the position intention of the target after coherent accumulation of signals of three continuous frames, and characterize that when the moving speed of the target exists, even after intra-frame distance walk calibration, distance walk still exists between frames. As shown in fig. 6(a), three points in the graph represent positions of objects in different frames, respectively, and fig. 6(b) is distance dimensional data corresponding to peaks in different frames.
And S6, performing interframe alignment on the multiframe coherent accumulation data obtained in the step S5.
Taking the number of frames K as 3 as an example, find F respectivelyxc1(n,t)、Fxc2(n,t)、Fxc3Position P where peak of (n, t) is located1(x1,y1)、P2(x2,y2)、P3(x3,y3) In a first frame Fxc1(n, t) is taken as a reference, the number dP of inter-frame distance walking units generated by the target in the 2 nd frame and the 3 rd frame is calculated, and the calculation method comprises the following steps:
dP2=y2-y1 (10)
dP3=y3-y1 (11)
respectively converting the 2 nd frame data Fxc2(n, t), frame No. 3Fxc3(n, t) cyclically shifting dP in the range dimension2、dP3Each distance unit obtains a new two-dimensional data plane and records the new two-dimensional data plane as Fnxc2(n,t)、Fnxc2And (n, t), realizing the alignment among the data of multiple frames.
And S7, performing non-coherent accumulation on the multi-frame coherent accumulation data after the inter-frame alignment.
Specifically, non-coherent accumulation is performed on the aligned data to obtain a new two-dimensional data plane Fnacc(n,t):
Fnacc(n,t)=Fxc1(n,t)+Fnxc2(n,t)+...+FnxK(n,t) (12)
FIG. 7(a) is a non-coherent accumulated two-dimensional data plane F obtained by performing modulo value calculation on three continuous frames of signals, performing data alignment with reference to the first frame, and accumulating the three frames of datanacc(n, t). Fig. 7(b) shows distance dimension data corresponding to the peak of fig. 7(a), and comparing fig. 7(b) and fig. 5(b), it can be seen that the signal-to-noise ratio is significantly improved.
S8, not based on step S7And carrying out target detection on the data obtained by coherent accumulation. Namely: data F based on non-coherent accumulationnaccAnd (n, t), calculating the signal-to-noise ratio and carrying out target detection.
In summary, the target detection method based on the non-coherent accumulation mode of the coherent radar is particularly suitable for the field of space-based target detection and alarm, typical long-time coherent accumulation scenes, and scenes such as ultrahigh-speed target motion, high-resolution radar and the like. The method is mainly used for solving the problem of detecting a moving target and a small target under a space-based warning radar, and under the condition that the target speed is unknown, the method realizes distance walking calibration by adopting an autocorrelation method based on pulses, and realizes the detection of the moving target accumulated by long-time phase coherence; meanwhile, the method utilizes the distance walking information to carry out preliminary estimation of the target speed; performing coherent accumulation on the calibrated data to obtain an accumulated two-dimensional data plane; after the two-dimensional data planes between the adjacent multiframes are subjected to target alignment, non-coherent accumulation is carried out, and then target detection is carried out, so that the aim of improving the signal-to-noise ratio is fulfilled, and the target detection probability of the warning radar is improved.
Furthermore, the method can solve the problems of distance walking, phase walking and the like caused by long-time coherent accumulation, can further improve the signal-to-noise ratio of target detection through multi-frame non-coherent accumulation, realizes stable detection and tracking and improves the detection probability of a radar system; the operation of the method has the characteristics of relatively small operation amount and high detection speed, and further meets the requirement of the space-based alarm radar for quick alarm.
Furthermore, the method realizes distance walking calibration by adopting a cross-correlation method based on pulses, realizes detection of a moving target accumulated by long-time phase coherence, and has the characteristic of small operand.
Furthermore, the method realizes the data alignment of a two-dimensional data plane through the peak value alignment among the multi-frame signals, and then carries out the non-coherent accumulation among the multi-frame signals, thereby being suitable for the static target detection and the target detection with speed, improving the signal-to-noise ratio of the detection, and having simple calculation process and small operand.
Furthermore, the method for realizing distance walking calibration and peak alignment among multi-frame signals based on the inter-pulse autocorrelation method has the advantages that target speed priori knowledge is not needed, the method can be well suitable for application scenes of space-based warning radars, the rapid estimation of the motion speed of the target can be realized through the algorithm provided by the patent, and the method can be more suitable for application scenes without target priori knowledge.
Furthermore, the non-coherent accumulation method provided by the method can further improve the detection signal-to-noise ratio and improve the target detection probability of the radar on the basis of the traditional coherent accumulation.
While the present invention has been described in detail with reference to the preferred embodiments, it should be understood that the above description should not be taken as limiting the invention. Various modifications and alterations to this invention will become apparent to those skilled in the art upon reading the foregoing description. Accordingly, the scope of the invention should be determined from the following claims.

Claims (7)

1. A target detection method based on a non-coherent accumulation mode of coherent radar is characterized by comprising the following steps:
s1, sampling to obtain a frame of two-dimensional echo data;
s2, performing pulse compression on the sampling data of each pulse of the two-dimensional echo data in the step S1;
s3, carrying out distance walking calibration on the data after pulse compression;
s4, carrying out coherent accumulation on the data after the distance walk calibration;
s5, repeating the steps S1-S4 to obtain multi-frame coherent accumulation data;
s6, performing interframe alignment on the multiframe coherent accumulation data obtained in the step S5;
s7, performing non-coherent accumulation on the multi-frame coherent accumulation data after the inter-frame alignment;
s8, target detection is performed based on the non-coherent accumulated data of the step S7.
2. The method for target detection based on non-coherent accumulation mode of coherent radar according to claim 1,
in step S1, the two-dimensional echo data is in the form of LFM signal, the number of pulses of each frame of LFM signal is N, each pulse samples M points, and the pulse repetition period is TrPulse width of TPEach sub-pulse has a frequency modulation slope of k and a sub-pulse bandwidth of B, and k is B/TP(ii) a Setting the initial radial distance of the uniform motion target as R0Radial velocity of motion v0The target distance of the ith pulse is delayed by taui=2(R0+v0iTr) C, then the target baseband echo signal sr(n, t) is:
Figure FDA0003282628480000011
wherein:
Figure FDA0003282628480000012
in the above formula, c represents the speed of light, and the length of t is M;
simultaneously, noise n (n, t) with a certain signal-to-noise ratio is superposed on the target baseband echo signal Noi (n, t) to obtain two-dimensional echo data SNr(n,t):
SNr(n,t)=Noi(n,t)+sr(n,t) (3)。
3. The method for target detection based on the non-coherent accumulation mode of coherent radar according to claim 2,
in step S2, the matching function of pulse compression is:
Figure FDA0003282628480000021
the pulse compression is carried out in a frequency domain, and the two-dimensional echo data of the ith pulse is set as SNri(n, t), the pulse pressure thereofThe reduction algorithm is as follows:
sci(n,t)=IFFT(FFT(SNri(n,t))·conj(FFT(sref(t)))) (5)
wherein, sci(N, t) is pulse pressure data of the ith pulse, i is 1,2, …, N.
4. The method for target detection based on non-coherent accumulation mode of coherent radar according to claim 1,
in step S3, the distance walk calibration is performed on the data after the pulse compression by using the inter-pulse cross correlation method.
5. The method for target detection based on the non-coherent accumulation mode of coherent radar according to claim 3,
in step S3, the pulse pressure data sc of the 1 st pulse is used1(t) data sc obtained by compressing the pulse of the ith pulse with reference toi(t) pulse pressure data s of the 1 st pulsecm1(t) performing a cross-correlation process, a cross-correlation function R1i(t) is:
Figure FDA0003282628480000022
find R1i(t) position P where peak value is locatedXi,PXiI.e. represents sci(t) relative to sc1(t) the number of distance units, and pulse pressure data sc of the ith pulsei(t) Cyclic Displacement PXiPoint, obtain scxi(t), implementing distance walk calibration; distance walk calibration is carried out on the data after pulse compression of each pulse, distance walk calibration of one frame of two-dimensional echo data is realized, and a new two-dimensional data plane S is formedxc(n,t)。
6. The method for target detection based on non-coherent accumulation mode of coherent radar according to claim 5,
sci(n, t) two-dimensional arrayIn the cross-correlation function R of the last pulse1N(t) position P where peak value is locatedXNWhich represents the number of distance units P traveled by the distance generated between the Nth pulse and the 1 st pulseXNThen, the corresponding walking distance dR is:
dR=C·PXN·ts (7)
the estimated target speed veComprises the following steps:
Figure FDA0003282628480000031
7. the method for target detection based on non-coherent accumulation mode of coherent radar according to claim 5,
in the step S4, a new two-dimensional data plane S is processedxc(n, t) performing coherent accumulation, and obtaining the column vector L corresponding to each distance unitxcm(N) coherent accumulation, where M is 1,2, …, M, N is 1,2, …, N, i.e. the column vector Lxcm(n) performing FFT to obtain Fxcm(n):
Fxcm(n)=FFT(Lxcm(n)) (9)
For each column vector Lxcm(n) all carry out coherent accumulation to form a new two-dimensional data plane Fxc1(n,t)。
CN202111139490.XA 2021-09-27 2021-09-27 Target detection method based on non-coherent accumulation mode of coherent radar Pending CN113885003A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111139490.XA CN113885003A (en) 2021-09-27 2021-09-27 Target detection method based on non-coherent accumulation mode of coherent radar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111139490.XA CN113885003A (en) 2021-09-27 2021-09-27 Target detection method based on non-coherent accumulation mode of coherent radar

Publications (1)

Publication Number Publication Date
CN113885003A true CN113885003A (en) 2022-01-04

Family

ID=79007179

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111139490.XA Pending CN113885003A (en) 2021-09-27 2021-09-27 Target detection method based on non-coherent accumulation mode of coherent radar

Country Status (1)

Country Link
CN (1) CN113885003A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115166648A (en) * 2022-09-08 2022-10-11 北京轩涌科技发展有限公司 Low signal-to-noise ratio radar signal processing method and device
CN116184346A (en) * 2022-12-27 2023-05-30 南京信息工程大学 5G downlink signal external radiation source radar coherent accumulation processing method

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115166648A (en) * 2022-09-08 2022-10-11 北京轩涌科技发展有限公司 Low signal-to-noise ratio radar signal processing method and device
CN116184346A (en) * 2022-12-27 2023-05-30 南京信息工程大学 5G downlink signal external radiation source radar coherent accumulation processing method
CN116184346B (en) * 2022-12-27 2023-11-03 南京信息工程大学 5G downlink signal external radiation source radar coherent accumulation processing method

Similar Documents

Publication Publication Date Title
CN103630888B (en) High-precision real-time microwave velocity and distance measuring device based on symmetrical triangle LFMCW (Linear Frequency Modulation Continuous Wave) radar
CN101825707B (en) Monopulse angular measurement method based on Keystone transformation and coherent integration
CN113885003A (en) Target detection method based on non-coherent accumulation mode of coherent radar
CN103616679B (en) Based on difference beam modulation and the PD radar range finding angle-measuring method of wave form analysis
Bar-Ilan et al. Sub-Nyquist radar via Doppler focusing
CN102628937B (en) Radar detection method based on generalized keystone transformation and non-coherent accumulation
CN104698459B (en) Stripe SAR (specific absorption resolution) compressed sensing and imaging method for missing data
CN110837081B (en) High-speed target detection method based on P/D (Peer-to-Peer) band radar signal fusion processing
CN107450055B (en) High-speed maneuvering target detection method based on discrete linear frequency modulation Fourier transform
CN107356922B (en) Multi-frequency radar target tracking method for ambiguity resolution by spread-table-aided method
CN104330791B (en) A kind of correlative accumulation method based on frequency domain shear
CN107356908B (en) Frequency agile signal coherent accumulation method
CN104849708A (en) High-speed maneuvering target parameter estimation method based on frequency domain polynomial phase transformation
CN111884984A (en) Fast carrier Doppler frequency shift capturing system
CN108809354B (en) Long-time accumulation method for high dynamic Doppler capture
CN110943758A (en) Secondary capturing method based on Doppler frequency change frequency and frequency joint search
CN114553260A (en) High-precision measurement system for DS/FH spread spectrum signal carrier frequency
CN106788584B (en) Improved PMF-FFT (pulse-modulated fast Fourier transform-fast Fourier transform) capturing method suitable for long code word application
CN111007473B (en) High-speed weak target detection method based on distance frequency domain autocorrelation function
CN109581350B (en) Radar distance and speed measurement method and device based on time-frequency integral interpolation
CN112968850A (en) Weak BPSK signal carrier capturing method
CN113820678B (en) Coherent accumulation method for high repetition frequency hopping signals of external radiation source radar
CN115685169A (en) Underwater sound weak moving target detection method based on broadband keystone transformation
CN108111452B (en) Method for estimating arrival time of broadband linear frequency modulation pulse signal
CN113030895B (en) Multi-frame coherent accumulation detection method for weak target

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