CN110398730B - Maneuvering target coherent detection method based on coordinate rotation and non-uniform Fourier transform - Google Patents
Maneuvering target coherent detection method based on coordinate rotation and non-uniform Fourier transform Download PDFInfo
- Publication number
- CN110398730B CN110398730B CN201910562876.8A CN201910562876A CN110398730B CN 110398730 B CN110398730 B CN 110398730B CN 201910562876 A CN201910562876 A CN 201910562876A CN 110398730 B CN110398730 B CN 110398730B
- Authority
- CN
- China
- Prior art keywords
- target
- fourier transform
- maneuvering
- maneuvering target
- signal
- 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.)
- Expired - Fee Related
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 42
- 230000001427 coherent effect Effects 0.000 title claims abstract description 33
- 230000005012 migration Effects 0.000 claims abstract description 38
- 238000013508 migration Methods 0.000 claims abstract description 38
- 238000000034 method Methods 0.000 claims abstract description 32
- 230000033001 locomotion Effects 0.000 claims abstract description 19
- 238000012545 processing Methods 0.000 claims abstract description 14
- 230000001133 acceleration Effects 0.000 claims description 33
- 238000009825 accumulation Methods 0.000 claims description 29
- 230000006870 function Effects 0.000 claims description 22
- 230000000694 effects Effects 0.000 claims description 21
- 230000009021 linear effect Effects 0.000 claims description 21
- 230000009466 transformation Effects 0.000 claims description 14
- 230000003595 spectral effect Effects 0.000 claims description 10
- 230000008878 coupling Effects 0.000 claims description 7
- 238000010168 coupling process Methods 0.000 claims description 7
- 238000005859 coupling reaction Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 5
- 230000008030 elimination Effects 0.000 claims description 2
- 238000003379 elimination reaction Methods 0.000 claims description 2
- 238000005316 response function Methods 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 abstract description 46
- 238000004364 calculation method Methods 0.000 abstract description 10
- 238000004088 simulation Methods 0.000 abstract description 10
- 238000010586 diagram Methods 0.000 description 8
- 230000006835 compression Effects 0.000 description 5
- 238000007906 compression Methods 0.000 description 5
- 230000014509 gene expression Effects 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 230000036039 immunity Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000012292 cell migration Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035485 pulse pressure Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
- G01S13/581—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of interrupted pulse modulated waves and based upon the Doppler effect resulting from movement of targets
- G01S13/582—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of interrupted pulse modulated waves and based upon the Doppler effect resulting from movement of targets adapted for simultaneous range and velocity measurements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
- G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/295—Means for transforming co-ordinates or for evaluating data, e.g. using computers
- G01S7/2955—Means for determining the position of the radar coordinate system for evaluating the position data of the target in another coordinate system
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details 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/415—Identification of targets based on measurements of movement associated with the target
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention belongs to the technical field of radar target signal estimation, and particularly relates to a maneuvering target coherent detection method based on coordinate rotation and non-uniform Fourier transform, which comprises the following steps: aiming at radar signals, establishing a maneuvering target signal model; and performing range migration processing on the maneuvering target signal model, and estimating target motion parameters through coordinate rotation and non-uniform Fourier transform. The method can relieve the contradiction between the calculation complexity and the detection performance, eliminates QRM by utilizing second-order Keystone transform (SoKT), and constructs a phase compensation function by combining the relation between the rotation angle and the Doppler frequency; and simulation experiments verify that compared with the existing representative algorithm, the method can obtain nearly ideal detection performance with lower operation complexity and has stronger application prospect.
Description
Technical Field
The invention belongs to the technical field of radar target signal estimation, and particularly relates to a maneuvering target coherent detection method based on coordinate rotation and non-uniform Fourier transform.
Background
In recent years, with the rapid development of stealth airplanes and Unmanned Aerial Vehicles (UAVs), the radar maneuvering weak target detection is required to be higher and higher. To detect such low radar cross-sectional area (RCS) targets, long-term coherent accumulation is an indispensable tool. Unfortunately, due to the high mobility of the above targets, Linear Range Migration (LRM), Quadratic Range Migration (QRM) and Doppler Frequency Migration (DFM) are difficult to avoid for the targets within the accumulation time. These adverse factors will severely affect the detection performance of traditional accumulation algorithms (e.g. MTD). Therefore, how to robustly detect a maneuvering weak target becomes a hot topic in the field of radar signal processing.
The LRM will cause a cross-range cell migration of the target envelope. To eliminate the LRM, many successful algorithms are proposed, such as inverse variable-scale Fourier transform (SCIFT), Keystone Transform (KT), modified coordinate rotation transform (MLRT), and Radon-Fourier transform. These algorithms enable satisfactory detection capabilities with moderate computational effort. However, they suffer from severe cumulative performance loss due to neglect of QRM and DFM effects caused by target acceleration. To solve this problem, researchers have expended a great deal of effort and developed many approaches, which can roughly be divided into two categories in general: (a) search-type algorithms, representative methods include generalized Radon-Fourier transform (GRFT), KT and Lu-distribution algorithms (KT-LVD), KT and linear canonical transform (KT-LCT), modified axial rotation transform and Lu-distribution (MART-LVT), modified axial rotation and fractional Fourier transform (IAR-FRFT), and Radon-Lu-distribution (RLVD). These algorithms exhibit superior detection performance at low signal-to-noise ratios (SNR) through parameter searching. However, the great computational complexity makes it unacceptable in practical applications. (b) Non-search class algorithms, typical methods include symmetric autocorrelation function-variable scale fourier transform (SAF-SFT), three-dimensional variable scale transform (TDST), frequency second order phase difference (FD-SoPD), and Adjacent Cross Correlation Function (ACCF). These algorithms reduce the coupling order by correlation operations, greatly reducing the computational burden. However, the correlation belongs to a nonlinear operation, which causes cross term interference under multiple targets and reduces the noise immunity.
Disclosure of Invention
Therefore, the invention provides the maneuvering target coherent detection method and device based on coordinate rotation and non-uniform Fourier transform, which can relieve the contradiction between the calculation complexity and the detection performance, obtain nearly ideal detection performance with lower computation amount and have stronger application prospect.
According to the design scheme provided by the invention, the maneuvering target coherent detection method based on coordinate rotation and non-uniform Fourier transform comprises the following contents:
A) aiming at radar signals, establishing a maneuvering target signal model;
B) and performing range migration processing on the maneuvering target signal model, and estimating target motion parameters through coordinate rotation and non-uniform Fourier transform.
In the above, a) is to establish a maneuvering target signal model, which includes the following contents: assuming that the radar transmits a linear frequency modulation signal, the maneuvering target moves at constant acceleration, and the instantaneous slope distance of the target from the radar is obtained according to the slow time variable, the initial slope distance of the target, the radial speed and the acceleration; obtaining radar receiving signals according to the instantaneous slope distance, the target radial speed, the signal amplitude and the bandwidth; and carrying out Fourier transform on the radar receiving signals along a fast time axis to obtain a maneuvering target signal model.
And B) eliminating the effects of linear range migration, secondary range migration and Doppler frequency migration.
Preferably, in B), second-order Keystone transformation is adopted to eliminate second-order range migration caused by the acceleration of the maneuvering target.
Preferably, the second-order Keystone transform eliminates the second-order range migration process, and comprises the following contents: scaling the slow time of each distance frequency unit, and bringing the scaled slow time variable into a maneuvering target signal model; and simplifying and representing the maneuvering target signal model obtained after scaling through Taylor expansion, performing Fourier inverse transformation along the range frequency, completing secondary coupling between Doppler frequency and slow time variable, and correcting a secondary range migration effect caused by the maneuvering target.
Preferably, in B), the jointly estimated target motion parameters include the following: representing the maneuvering target signal model as a discrete form under a coordinate system; eliminating residual linear distance migration effect by introducing coordinate rotation transformation, extracting signal energy along the azimuth direction, and estimating the equivalent speed, the unambiguous speed and the real speed of the maneuvering target by using the rotation angle; constructing a phase compensation function according to the estimated non-fuzzy speed, and multiplying the phase compensation function by the extracted signal energy to compensate a linear phase term; and carrying out parameter estimation by utilizing non-uniform Fourier transform.
Preferably, in the parameter estimation by using the non-uniform fourier transform, the non-uniform fourier transform is performed on the multiplication result, and the maneuvering target energy accumulation is determined to be a single peak condition according to the non-uniform fourier transform, and the maneuvering target acceleration is estimated from the single peak condition.
Preferably, a phase compensation function is constructed according to the estimated target speed and acceleration; carrying out slow time Fourier transform and range frequency Fourier inverse transform on the maneuvering target signal model to complete coherent accumulation; obtaining a single spectral peak of a distance-Doppler domain after the mobile target receiving signal is accumulated; target detection is accomplished by a single spectral peak.
Preferably, the target detection is finished by using a constant false alarm aiming at a single spectrum peak value, if the peak value exceeds a self-adaptive threshold, a moving target is judged to exist, and a moving target parameter is obtained; and if the peak value is smaller than the threshold, judging that the moving target is not detected.
Furthermore, the invention also provides a coherent detection device for a maneuvering target based on coordinate rotation and non-uniform Fourier transform, which comprises: a model building module and a joint estimation module, wherein,
the model establishing module is used for establishing a maneuvering target signal model aiming at the radar signal;
and the joint estimation module is used for carrying out range migration processing on the maneuvering target signal model and estimating target motion parameters through coordinate rotation and non-uniform Fourier transform.
The invention has the beneficial effects that:
aiming at radar signals, a maneuvering target signal model is established; the method comprises the steps of carrying out range migration processing on a maneuvering target signal model, and estimating target motion parameters through coordinate rotation and non-uniform Fourier transform in a combined manner, so that the contradiction between the calculation complexity and the detection performance is effectively relieved; and second-order Keystone transform (SoKT) is utilized to eliminate QRM, and a phase compensation function is constructed by combining the relation between the rotation angle and the Doppler frequency, so that almost ideal detection performance can be obtained with lower operation complexity compared with the conventional representative algorithm. Finally, the effectiveness of the technical scheme of the invention is further verified through simulation and actual measurement radar data processing results.
Description of the drawings:
FIG. 1 is a flow chart of a coherent detection method in an embodiment;
FIG. 2 is a schematic diagram of an embodiment of a coherent detection device;
FIG. 3 is a flow chart of the LRT-NuFFT algorithm in an embodiment;
FIG. 4 is a comparison of algorithm computation complexity in the examples;
FIG. 5 is a schematic diagram of coherent accumulation simulation results of a single maneuvering target in an embodiment;
FIG. 6 is a diagram illustrating a multi-target coherent accumulation simulation result of the LRT-NuFFT algorithm in the embodiment;
FIG. 7 is a graph showing the variation of the target detection probability with SNR in the example;
fig. 8 is a schematic diagram of the processing result of the measured UAV data in the embodiment.
The specific implementation mode is as follows:
in order to make the objects, technical solutions and advantages of the present invention clearer and more obvious, the present invention is further described in detail below with reference to the accompanying drawings and technical solutions.
The long-time coherent accumulation can obviously improve the detection and motion parameter estimation capability of the radar on the maneuvering target. However, Linear Range Migration (LRM), Quadratic Range Migration (QRM), and Doppler Frequency Migration (DFM) during the coherent processing severely degrade radar detection and estimation performance. In view of this, in the embodiment of the present invention, referring to fig. 1, a method for coherent detection of a maneuvering target based on coordinate rotation and non-uniform fourier transform is provided, which includes the following steps:
s101) aiming at radar signals, establishing a maneuvering target signal model;
s102) carrying out range migration processing on the maneuvering target signal model, and jointly estimating target motion parameters through coordinate rotation and non-uniform Fourier transform.
Further, in the embodiment of the invention, assuming that the radar transmits a linear frequency modulation signal, the maneuvering target moves at a constant acceleration, and the instantaneous slope distance of the target from the radar is obtained according to the slow time variable, the initial slope distance of the target, the radial speed and the acceleration; obtaining radar receiving signals according to the instantaneous slope distance, the target radial speed, the signal amplitude and the bandwidth; and carrying out Fourier transform on the radar receiving signals along a fast time axis to obtain a maneuvering target signal model.
The radar transmits a Linear Frequency Modulated (LFM) signal,
wherein,
Tpand fcFor the pulse duration and the carrier frequency of the signal,and gamma denotes the fast time variation and chirp rate, respectively.
The maneuvering target moves with constant acceleration, its instantaneous slope distance R (t) from the radarm) Satisfy the requirement of
Wherein R is0V and a are the initial slope, radial velocity and acceleration of the target, respectively. t is tm=mT(m=1,2,…,Na) Representing a slow time variable, T being the pulse repetition time, NaTo accumulate the number of pulses.
Neglecting the effects of noise, the received signal after pulse pressure can be expressed as:
wherein A iscAnd B signal amplitude and bandwidth, respectively.
The formula (2) can be substituted into the formula (3)
Wherein λ ═ c/fcIs the signal wavelength.
Doppler ambiguity often occurs due to the low pulse repetition frequency of the radar and the high velocity of the target. Thus, the radial velocity of the target may be expressed as
v=Nbvb+v0 (5)
Wherein v isb=λfpA/2 and fpRespectively representing radar blind speed and pulse repetition frequency. N is a radical ofbIs a target Doppler ambiguity number, v0Is free of blurring speed and satisfies-vb/2≤v0<vb/2。
The formula (5) is inserted into the formula (4) to obtain:
in the above equation, the equation exp (-j2 π f is usedcNbvbtm/c)=1。
It can be observed from equation (6) that the signal envelope varies non-linearly with slow time. The LRM effect occurs when the offset exceeds a distance resolution unit Δ r ═ c/2B during the integration time. The QRM effect can also be observed if the target has high mobility (i.e. possesses large acceleration). Fourier Transform (FT) is performed on equation (6) along the fast time axis, and the Fourier Transform (FT) can be obtained
Wherein f isrIs relative toThe doppler frequency of (d). Formula (7) shows thatrAnd tmThe coupling relationship between them is the root cause for LRM and QRM. The doppler frequency brought by the target is defined as:
due to the acceleration of the target, a linear Doppler shift occurs, with a frequency shift of 2atmAnd/lambda. If the value exceeds the Doppler resolution cell 1/NaT, then the DFM effect occurs and causes the doppler domain energy to diverge. Therefore, in order to detect a maneuvering target at low signal-to-noise ratio, the effects of LRM, QRM, and DFM must be effectively eliminated. Therefore, in the embodiment of the invention, the range migration processing comprises the elimination processing of the effects of linear range migration, quadratic range migration and Doppler frequency migration.
Further, the second-order Keystone transformation is adopted to eliminate the second-order range migration caused by the acceleration of the maneuvering target in the embodiment of the invention. Preferably, the second-order Keystone transform eliminates the second-order range migration process, and comprises the following contents: scaling the slow time of each distance frequency unit, and bringing the scaled slow time variable into a maneuvering target signal model; and simplifying and representing the maneuvering target signal model obtained after scaling through Taylor expansion, performing Fourier inverse transformation along the range frequency, completing secondary coupling between Doppler frequency and slow time variable, and correcting a secondary range migration effect caused by the maneuvering target.
Second order coupling caused by target acceleration is eliminated by using a second order Keystone transform (SoKT). SoKT scales the slow time of each range frequency unit, which can be expressed as
Wherein, taIs a scaled slow time variable.
The formula (9) may be substituted for the formula (7):
for narrow band radar systems, one-order taylor spreads are typically satisfied:
therefore, the formula (10) can be simplified to
Wherein, Ve=Nbvb+v0The/2 is defined as the equivalent speed of the target.
Distance frequency f along equation (12)rBy performing an Inverse Fourier Transform (IFT), we can obtain
As can be observed from formulae (12) and (13), SoKT eliminates frAnd tmA second order coupling between them. Thus, the QRM effect caused by the target is effectively corrected.
Further, in the embodiment of the present invention, the jointly estimating the target motion parameter includes the following contents: representing the maneuvering target signal model as a discrete form under a coordinate system; eliminating residual linear distance migration effect by introducing coordinate rotation transformation, extracting signal energy along the azimuth direction, and estimating the equivalent speed, the unambiguous speed and the real speed of the maneuvering target by using the rotation angle; constructing a phase compensation function according to the estimated non-fuzzy speed, and multiplying the phase compensation function by the extracted signal energy to compensate a linear phase term; and carrying out parameter estimation by utilizing non-uniform Fourier transform. Preferably, in the parameter estimation by using the non-uniform fourier transform, the non-uniform fourier transform is performed on the multiplication result, and the maneuvering target energy accumulation is determined to be a single peak condition according to the non-uniform fourier transform, and the maneuvering target acceleration is estimated from the single peak condition.
The residual LRM and QFM in equation (13) pose a significant difficulty for coherent accumulation. The embodiment of the invention provides the method for eliminating LRM and QFM by comprehensively utilizing LRT and NuFFT. Notice ta=mT,fsKB, where fsIs the signal sampling rate and K is the over-sampling rate. Therefore, R ═ ρ n, R0=ρnR0Where ρ ═ c/2fsDistance sampling unit, n and nR0Are respectively R and R0The distance unit identification bit. Equation (13) can be expressed in discrete form in the (n, m) coordinate system, i.e.
In order to compensate for the LRM effect caused by the equivalent velocity, an LRT is introduced, which is defined as a coordinate rotation operation as shown below,
Formula (15) may be substituted for formula (14):
as shown in equation (17), the target energies are concentrated in the same range bin. The signal energy is then extracted in the azimuth direction. For a particular range bin, the extracted signal can be expressed as:
obviously, sazi(m') is a chirp signal. The LVT is adopted to perform coherent accumulation and motion parameter estimation, and although satisfactory accumulation performance and noise immunity are exhibited, huge calculation amount is daunting.
In the embodiment of the present invention, it is noted from the expressions (17) and (18) that the target equivalent speed VeNo blurring velocity v0And the true velocity v can be estimated simultaneously by means of the angle of rotation, i.e.
According to equation (12), the unambiguous velocity can be solved by
Wherein ROUND (-) is an integer function,is the doppler ambiguity number of the target. So the true velocity of the target is estimated as
the linear phase term can be compensated by multiplying equation (22) by equation (18). Then, energy accumulation and parameter estimation can be completed by utilizing efficient NuFFT:
As can be seen from equation (23), the target energy accumulates as a single peak, and from its position, the acceleration of the target can be estimated as:
further, in the embodiment of the invention, a phase compensation function is constructed according to the estimated target speed and acceleration; carrying out slow time Fourier transform and range frequency Fourier inverse transform on the maneuvering target signal model to complete coherent accumulation; obtaining a single spectral peak of a distance-Doppler domain after the mobile target receiving signal is accumulated; target detection is accomplished by a single spectral peak. Preferably, the target detection is finished by using a constant false alarm aiming at a single spectrum peak value, if the peak value exceeds a self-adaptive threshold, a moving target is judged to exist, and a moving target parameter is obtained; and if the peak value is smaller than the threshold, judging that the moving target is not detected.
And constructing a phase compensation function according to the estimated target speed and acceleration to eliminate the influence of LRM, QRM and DFM:
finally, the coherent accumulation can be completed by carrying out the slow time FT and the distance frequency IFT on the formula (7):
wherein A isCIIs a complex amplitude, fdTo correspond to a slow time tmThe doppler frequency of (d).
In equation (26), the received signal of the maneuvering target is accumulated as a single spectral peak in the range-doppler domain, whose peak value is | SCI(2R0/c,-2v0Lambda) |. Target detection can then be performed using Constant False Alarm (CFAR) techniques:
if the peak value exceeds the self-adaptive threshold eta, the moving object is indicated, and the moving parameter isAndotherwise, if | SCI(2R0/c,-2v0And/lambda) is less than the threshold eta, no moving object is detected.
Suppose Q targets are observed within the coherent integration time. Similar to equation (14), the signal after passing through SoKT can be expressed as:
LRT of formula (28) gives:
when in useThen the LRM of the ith target is corrected, while the LRMs of the remaining q-1 targets still exist, i.e.
Then, along the distance unit nR0,iExtracting the direction signal can obtain:
wherein,
similar to equations (19) to (22), the phase compensation function is constructed from the rotation angle
Multiplying equation (31) by equation (33) and performing NuFFT to focus the target energy and estimate acceleration yields:
note that only the ith target energy is concentrated to a single peak, while the other target energies are divergent, mainly due to two reasons: (a) the envelopes of the other q-1 objects are distributed among different range bins; (b) the linear phase term exp (-j4 π v) in equation (32)0,qm' T/λ) cannot be compensated by the phase compensation function H1(m') so that NuFFT cannot gather its energy.
Finally, coherent accumulation and parameter estimation of the ith target can be completed, however, the rest q-1 targets cannot realize focusing by using the parameters of the ith target.
LRT-NuFFT correction of LRM is realized through rotation angle search. However, the LRT process requires a large number of interpolation operations, which inevitably results in numerical errors, and in addition, the rotation angle is not a natural representation of the target motion parameters. Because, in order to solve this problem, the embodiment of the present invention further provides an interpolation-free implementation form of LRT-NuFFT, i.e. speed search-NuFFT (VS-NuFFT).
Considering the coordinate rotation transformation in equation (16), it can be rewritten as:
FT along the fast time for equation (35) can be obtained:
as can be seen from equation (36), LRT operation can be achieved by phase compensation from the frequency domain, where the compensation function is
When V iss=Ve,frAnd m' are eliminated. IFT is performed along the range frequency for equation (36) to obtain:
it is noted that the velocity search interval should be less than Δ v ═ λ/2N, depending on doppler frequency resolutionaT。
Further, based on the above method, an embodiment of the present invention further provides a coherent detection apparatus for a maneuvering target based on coordinate rotation and non-uniform fourier transform, as shown in fig. 2, including: a model building module 101 and a joint estimation module 102, wherein,
the model establishing module 101 is used for establishing a maneuvering target signal model aiming at radar signals;
and the joint estimation module 102 is used for performing range migration processing on the maneuvering target signal model and estimating target motion parameters through coordinate rotation and non-uniform Fourier transform.
In order to further verify the effectiveness of the invention, the following explains the technical scheme of the invention by combining a specific implementation algorithm and a simulation experiment:
referring to fig. 3, the specific steps of the LRT-NuFFT algorithm can be designed as follows:
Step 2: the distance direction FT. SoKT eliminates QRM. Distance IFT is given as sSoKT(n,m);
And step 3: according to the detectionDetermining a rotation angle search rangeSearch intervalShould be less than arctan (lambda/2N)aρ);
And 4, step 4: LRT is carried out on the formula (17) under a certain rotation angle, and then a rotated signal s is obtainedrot(n′,m′);
And 5: constructing a phase compensation function H according to the rotation angle1(m') as shown in formulas (19) to (22). Extracting the slow time signal to obtain sazi(m′);
Step 6: to sazi(m′)·H1(m') performing NuFFT to realize signal energy accumulation and estimate acceleration, as shown in equation (23);
and 7: when the initial distance and the search angle match with the target true value, the output of NuFFT reaches a maximum value. Estimating a target speed and acceleration from equations (21) to (24);
and 8: construction of a further phase compensation function H using the estimated motion parameters2(fr,tm) And realizing coherent accumulation as shown in formulas (25) and (26).
Compared with the existing MART-LVT and LRT-NuFFT algorithms, the method has the following improvements and advantages: (1) in MART-LVT, QRM effects due to target acceleration are ignored, which can cause cumulative performance loss in certain scenarios. Whereas in LRT-NuFFT the QRM is first eliminated with SoKT, which makes the algorithm more suitable for highly mobile object detection. (2) Unlike MART-LVT, LRT-NuFFT utilizes the rotation angleAnd a no-blur velocity v0The relationship between them, the phase compensation is innovatively constructed. Therefore, the subsequent energy accumulation and acceleration estimation can be quickly achieved by NuFFT, which greatly eases the computational burden of computing the LVT algorithm. Compared with MART-LVT and KT-LVD. TDST, SAF-SFT and ACCF algorithm, LRT-NuFFT algorithm proposed in the embodiment of the present invention does not involve nonlinear operation, therefore, the algorithm does not have cross terms, and satisfies linear property; although LVD, TDST, SAF-SFT and ACCF have satisfactory cross term suppression capability, weak targets are still easily submerged in cross terms caused by strong targets. The rotation angle search is equivalent to the speed search, so that the interpolation process can be realized by utilizing FT, phase compensation and IFT, and the introduction of numerical errors is avoided; in addition, the speed search is more intuitive than the angle search, and the result is also the natural expression of the target motion parameter; in terms of search interval, a uniform search of angles corresponds to a non-uniform sampling of velocity, which is approximately equal at small angles but exhibits significant differences at large angles.
The LRT-NuFFT algorithm's calculated quantities were analyzed and compared to representative MART-LVT, KT-LVD, TDST, SAF-SFT, and MLRT algorithms as follows:
suppose thatNF、NrAnd NaThe angle search number, folding factor search number, distance cell number, and pulse number are respectively represented.
The main steps of the MART-LVT algorithm include the MART operationAnd LVT operationTherefore, the total calculation amount isMagnitude.
For the KT-LVD algorithm, the main step involves a folding factor search (O (N)FNr) And LVD operationTherefore, the total calculation amount is aboutMagnitude.
For the TDST algorithm, the need arisesAnda two-step scale-variable fourier transform (SFT) is calculated. Therefore, the total calculation amount is aboutMagnitude.
The main steps of the SAF-SFT algorithm include two SFT operations with a computational complexity of O (3N) respectivelyaNrlog2Na) Andtherefore, the total calculation amount is O (3N)a(Na+Nr)log2Na) Magnitude.
Unlike the MART-LVT algorithm, the algorithm proposed in the embodiment of the present invention uses phase compensation and NuFFT to reduce the computational complexity, so the total computation amount is also about
The calculated amount of the above algorithm is given in table 1. Suppose thatNKTFig. 4 shows a graphical representation of the calculated amount of the algorithm, 10. It is clear that the LRT-NuFFT algorithm has a lower computational load than MART-LVT, KT-LVD and TDST.
TABLE 1 different algorithm computational complexity
The algorithm of the invention is analyzed in combination with radar simulation data, and the simulation radar parameters are shown in table 2:
TABLE 2 simulated Radar parameters
First, the coherent accumulation capability of the proposed algorithm is evaluated with a single target. The target motion parameters are: r0=150km,v=150m/s,a=8m/s2. Figure 5(a) shows the pulse compression result, where the signal-to-noise ratio is-10 dB. Obviously, the target trajectory is overwhelmed by strong noise and cannot be observed. In addition, fig. 5(b) shows the accumulation result of MLRT. MLRT cannot detect maneuvering targets because QRM and DFM are not considered. Fig. 5(c) and 5(d) show the accumulation results of the SAF-SFT and TDST algorithms, respectively. Because the non-linear transformation loses some of the signal energy, the output results of both algorithms remain buried in noise. In contrast, MART-LVT and LRT-NuFFT are able to accumulate targets as a single peak at corresponding locations in the range-Doppler domain, as shown in FIGS. 5(e) and 5 (f). The simulation experiment shows the coherent accumulation capacity of the proposed algorithm in a low signal-to-noise ratio environment.
Secondly, the accumulation performance under multiple targets is analyzed. The motion parameter settings for the two motorized targets (Tr1 and Tr2) are given in table 3. Fig. 6(a) shows the pulse compression results. Fig. 6(b) shows the velocity search results, where two distinct spectral peaks give the equivalent velocities of Tr1 and Tr 2. From this, the true speeds of Tr1 and Tr2 are estimated as:the LRMs of the two targets can be corrected separately using this parameter, as shown in fig. 6(c) and 6 (d). The accelerations through NuFFT, Tr1 and Tr2 can be estimated in the corresponding distance units, resulting inAs shown in fig. 6(e) and 6 (f). Finally, FIG. 6(g) and FIG. 6(h) show the coherent integration results for the two targets, respectively. The simulation experiment verifies the accumulation capacity of the algorithm to multiple targets.
TABLE 3 two maneuvering target simulation parameters
It is noted that many spurious peaks can be observed in fig. 6(b), which are spaced at a half-blind speed. Therefore, these spurious peaks are referred to as semi-blind velocity side lobes (HBSSLs). In practical situations, a weak target may be submerged in the HBSSLs of a strong target and may be processed using the CLEAN algorithm.
The target detection performance of the proposed algorithm was evaluated by Monte Carlo experiments, where the signal-to-noise ratio after pulse compression varied from-25 dB to 5 dB. At each signal-to-noise ratio, 500 independent Monte Carlo experiments were performed. False alarm rate is set to Pfa=10-6. Five representative algorithms (MART-LVT, KT-LVD, TDST, SAF-SFT, and MLRT) were used as comparisons. The target detection probability curve is shown in fig. 7. Obviously, the LRT-NuFFT algorithm can obtain nearly ideal detection performance. This is because the LRT-NuFFT effectively eliminates the effects of LRM, QRM and DFM effects. MART-LVT and KT-LVD suffer some performance loss due to ignoring QRM caused by target acceleration. The detection thresholds of TDST and SAF-SFT algorithms are respectively 3dB and 6dB higher than the ideal case, which indicates that the bilinear transformation can lose a large amount of signal energy and deteriorate the anti-noise performance. MLRT does not take into account QRM and DFM effects and thus the detection performance is the worst.
The algorithm in the embodiment of the invention is further explained below by combining a commercial unmanned aerial vehicle 3 Xinjiang spirit, and data is collected in a certain campus. Fig. 8(a) and 8(b) show an experimental scenario and the Frequency Modulated Continuous Wave (FMCW) radar system used. The radar parameters are listed in table 4. To obtain doppler ambiguity of the target, the PRF of the radar is artificially reduced.
TABLE 4 FMCW Radar System parameters
Fig. 8(c) shows the target motion trajectory after pulse compression. Within the coherent accumulation time of 0.92s, the drone moves more than 7 range cells, causing severe range migration phenomena. Fig. 8(d) shows the equivalent velocity estimation result. The speed search range is set to [ -10,10 [)]m/s, the search interval is 0.01 m/s. From the peak position, the equivalent velocity of the UAV can be derived asTrue velocity ofThe LRM is effectively corrected by phase compensation, and the acceleration of the UAV is obtained through NuFFT, namelyAs shown in fig. 8(e) and 8(f), respectively. Fig. 8(g) shows the accumulation of LRT-NuFFT, and it can be seen that a well-focused spectral peak is formed in the range-doppler domain. Meanwhile, the accumulation results of the MTD and MLRT are given in fig. 8(h) and 8(i) for reference. However, it can be seen that the energy of the target is spread among multiple range and doppler cells, making target detection difficult; further verifies that the LRT-NuFFT algorithm can obtain nearly ideal detection performance with lower calculation amount in the embodiment of the invention, and verifies the effectiveness of the technical scheme in the embodiment of the invention through UAV experimental results.
Unless specifically stated otherwise, the relative steps, numerical expressions, and values of the components and steps set forth in these embodiments do not limit the scope of the present invention.
Based on the foregoing method, an embodiment of the present invention further provides a server, including: one or more processors; a storage device for storing one or more programs which, when executed by the one or more processors, cause the one or more processors to implement the method described above.
Based on the above method, the embodiment of the present invention further provides a computer readable medium, on which a computer program is stored, wherein the program, when executed by a processor, implements the above method.
The device provided by the embodiment of the present invention has the same implementation principle and technical effect as the method embodiments, and for the sake of brief description, reference may be made to the corresponding contents in the method embodiments without reference to the device embodiments.
It is clear to those skilled in the art that, for convenience and brevity of description, the specific working processes of the system and the apparatus described above may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
In all examples shown and described herein, any particular value should be construed as merely exemplary, and not as a limitation, and thus other examples of example embodiments may have different values.
It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures.
The flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present invention. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a non-volatile computer-readable storage medium executable by a processor. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
Finally, it should be noted that: the above-mentioned embodiments are only specific embodiments of the present invention, which are used for illustrating the technical solutions of the present invention and not for limiting the same, and the protection scope of the present invention is not limited thereto, although the present invention is described in detail with reference to the foregoing embodiments, those skilled in the art should understand that: any person skilled in the art can modify or easily conceive the technical solutions described in the foregoing embodiments or equivalent substitutes for some technical features within the technical scope of the present disclosure; such modifications, changes or substitutions do not depart from the spirit and scope of the embodiments of the present invention, and they should be construed as being included therein. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (4)
1. A maneuvering target coherent detection method based on coordinate rotation and non-uniform Fourier transform is characterized in that,
A) aiming at radar signals, establishing a maneuvering target signal model;
B) carrying out range migration processing on the maneuvering target signal model, and jointly estimating target motion parameters through coordinate rotation and non-uniform Fourier transform;
A) the method for establishing the maneuvering target signal model comprises the following steps: assuming that the radar transmits a linear frequency modulation signal, the maneuvering target moves at constant acceleration, and the instantaneous slope distance of the target from the radar is obtained according to the slow time variable, the initial slope distance of the target, the radial speed and the acceleration; obtaining radar receiving signals according to the instantaneous slope distance, the target radial speed, the signal amplitude and the bandwidth;
fourier transformation is carried out on the radar receiving signals along a fast time axis to obtain a maneuvering target signal model, and the maneuvering target signal model is expressed as follows:wherein f isrIs relative toDoppler frequency of fcIs the carrier frequency of the signal and,representing a fast time variable, tmDenotes a slow time variable, AfAnd B is the signal amplitude and bandwidth, v, respectivelyb=λfpA/2 and fpRespectively representing radar blind speed and pulse repetition frequency, NbIs a target Doppler ambiguity number, v0For no blurring speed, R0And a is the initial slope and acceleration of the target, respectively, and λ ═ c/fcIs the signal wavelength;
B) the intermediate range migration treatment comprises elimination treatment of linear range migration, secondary range migration and Doppler frequency migration;
B) in the method, second-order Keystone transformation is adopted to eliminate secondary range migration caused by maneuvering target acceleration; the second-order Keystone transformation eliminates the secondary range migration process, and comprises the following contents: scaling the slow time of each distance frequency unit, and bringing the scaled slow time variable into a maneuvering target signal model; simplifying and expressing a maneuvering target signal model obtained after scaling through Taylor expansion, performing Fourier inverse transformation along range frequency, completing secondary coupling between Doppler frequency and slow time variable, and correcting a secondary range migration effect caused by a maneuvering target;
B) in the method, the jointly estimated target motion parameters include the following: representing the maneuvering target signal model as a discrete form under a coordinate system; eliminating residual linear distance migration effect by introducing coordinate rotation transformation, extracting signal energy along the azimuth direction, and estimating the equivalent speed, the unambiguous speed and the real speed of the maneuvering target by using the rotation angle; constructing a phase compensation function according to the estimated non-fuzzy speed, and multiplying the phase compensation function by the extracted signal energy to compensate a linear phase term; carrying out parameter estimation by utilizing non-uniform Fourier transform;
in the parameter estimation by utilizing the non-uniform Fourier transform, firstly, the non-uniform Fourier transform is carried out on the multiplication result, the situation that the energy of the maneuvering target is accumulated into a single peak value is judged according to the non-uniform Fourier transform, and the acceleration of the maneuvering target is estimated from the single peak value situation, wherein the non-uniform Fourier transform is expressed as:the maneuvering target acceleration estimate is expressed as:p (-) represents the impulse response function of the non-uniform fourier transform,to correspond to (m' T)2Is the rotated coordinate, nR0Is R0T is the pulse repetition time and K is the over-sampling rate.
2. The method for coherent detection of maneuvering targets based on coordinate rotation and non-uniform Fourier transform according to claim 1, characterized in that a phase compensation function is constructed according to the estimated target speed and acceleration; carrying out slow time Fourier transform and range frequency Fourier inverse transform on the maneuvering target signal model to complete coherent accumulation; obtaining a single spectral peak of a distance-Doppler domain after the mobile target receiving signal is accumulated; target detection is accomplished by a single spectral peak.
3. The method for detecting the coherent moving target based on the coordinate rotation and the non-uniform Fourier transform as claimed in claim 1 or 2, wherein the target detection is completed by using a constant false alarm for a single spectral peak value, and if the peak value exceeds a self-adaptive threshold, a moving target is determined to exist, and a moving target parameter is obtained; and if the peak value is smaller than the threshold, judging that the moving target is not detected.
4. A mobile object coherent detection device based on coordinate rotation and non-uniform Fourier transform, which is realized based on the method of claim 1 and comprises the following steps: a model building module and a joint estimation module, wherein,
the model establishing module is used for establishing a maneuvering target signal model aiming at the radar signal;
and the joint estimation module is used for carrying out range migration processing on the maneuvering target signal model and estimating target motion parameters through coordinate rotation and non-uniform Fourier transform.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910562876.8A CN110398730B (en) | 2019-06-26 | 2019-06-26 | Maneuvering target coherent detection method based on coordinate rotation and non-uniform Fourier transform |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910562876.8A CN110398730B (en) | 2019-06-26 | 2019-06-26 | Maneuvering target coherent detection method based on coordinate rotation and non-uniform Fourier transform |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110398730A CN110398730A (en) | 2019-11-01 |
CN110398730B true CN110398730B (en) | 2021-07-06 |
Family
ID=68323539
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910562876.8A Expired - Fee Related CN110398730B (en) | 2019-06-26 | 2019-06-26 | Maneuvering target coherent detection method based on coordinate rotation and non-uniform Fourier transform |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110398730B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113109779B (en) * | 2020-02-28 | 2024-03-22 | 加特兰微电子科技(上海)有限公司 | Method for improving target detection precision, integrated circuit, radio device and electronic equipment |
CN112748427B (en) * | 2020-12-09 | 2024-04-23 | 中国人民解放军战略支援部队信息工程大学 | Joint estimation method for distance difference and distance difference change rate of pulse signals |
CN113708860B (en) * | 2021-10-26 | 2022-01-11 | 南京天朗防务科技有限公司 | Method and device for estimating LFM signal multipath time delay |
CN115828074B (en) * | 2022-09-08 | 2023-07-04 | 中国人民解放军军事科学院系统工程研究院 | Combined estimation method for target positioning parameters of radiation source |
CN117421576B (en) * | 2023-09-13 | 2024-05-17 | 中国人民解放军军事科学院系统工程研究院 | High-speed maneuvering target positioning parameter estimation method based on adjacent cross correlation |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104330791A (en) * | 2014-10-24 | 2015-02-04 | 上海无线电设备研究所 | Phase-coherent accumulation method based on frequency domain shear |
CN106896358A (en) * | 2017-04-27 | 2017-06-27 | 电子科技大学 | A kind of high-speed target phase-coherent accumulation detection method based on position rotation transformation |
CN108549067A (en) * | 2018-07-27 | 2018-09-18 | 电子科技大学 | A kind of phase-coherent accumulation detection method being applied to three rank maneuvering targets |
CN109521410A (en) * | 2018-11-16 | 2019-03-26 | 西安电子科技大学 | High-speed maneuver target phase-coherent accumulation detection method based on time reversal transformation |
CN109541568A (en) * | 2019-01-24 | 2019-03-29 | 中国人民解放军海军航空大学 | A kind of radar maneuvering target span from the quick phase-coherent accumulation detection method of doppler cells |
CN109581318A (en) * | 2019-01-10 | 2019-04-05 | 中国人民解放军海军航空大学 | Radar highly maneuvering target phase-coherent accumulation detection method based on time reversal nonuniform sampling |
CN109613507A (en) * | 2018-12-21 | 2019-04-12 | 北京理工大学 | A kind of detection method for high-order maneuvering target radar return |
-
2019
- 2019-06-26 CN CN201910562876.8A patent/CN110398730B/en not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104330791A (en) * | 2014-10-24 | 2015-02-04 | 上海无线电设备研究所 | Phase-coherent accumulation method based on frequency domain shear |
CN106896358A (en) * | 2017-04-27 | 2017-06-27 | 电子科技大学 | A kind of high-speed target phase-coherent accumulation detection method based on position rotation transformation |
CN108549067A (en) * | 2018-07-27 | 2018-09-18 | 电子科技大学 | A kind of phase-coherent accumulation detection method being applied to three rank maneuvering targets |
CN109521410A (en) * | 2018-11-16 | 2019-03-26 | 西安电子科技大学 | High-speed maneuver target phase-coherent accumulation detection method based on time reversal transformation |
CN109613507A (en) * | 2018-12-21 | 2019-04-12 | 北京理工大学 | A kind of detection method for high-order maneuvering target radar return |
CN109581318A (en) * | 2019-01-10 | 2019-04-05 | 中国人民解放军海军航空大学 | Radar highly maneuvering target phase-coherent accumulation detection method based on time reversal nonuniform sampling |
CN109541568A (en) * | 2019-01-24 | 2019-03-29 | 中国人民解放军海军航空大学 | A kind of radar maneuvering target span from the quick phase-coherent accumulation detection method of doppler cells |
Non-Patent Citations (3)
Title |
---|
A Coherent Detection and Velocity Estimation Algorithm for the High-Speed Target Based on the Modified Location Rotation Transform;Zhi Sun 等;《IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING》;20180731;第11卷(第7期);第2346-2361页 * |
高机动小RCS目标长时间相参积累检测新方法;战立晓 等;《系统工程与电子技术》;20130331;第35卷(第3期);全文 * |
高速高机动雷达目标检测与参数估计算法研究;周莹;《中国优秀硕士学位论文全文数据库 信息科技辑(月刊)》;20190215(第02期);正文第10-19、60-62页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110398730A (en) | 2019-11-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110398730B (en) | Maneuvering target coherent detection method based on coordinate rotation and non-uniform Fourier transform | |
CN108549067B (en) | Coherent accumulation detection method applied to third-order maneuvering target | |
CN106772309B (en) | Radar target detection method based on telecommunication satellite radiation source | |
CN107561508B (en) | Coherent accumulation detection method for uniformly accelerated moving target | |
CN108089171B (en) | A kind of radar rapid detection method for unmanned plane target | |
CN107450055B (en) | High-speed maneuvering target detection method based on discrete linear frequency modulation Fourier transform | |
CN111736128A (en) | Coherent accumulation method based on SKT-SIAF-MSCFT | |
CN110646774B (en) | Maneuvering target coherent detection method and device based on product variable-scale periodic Lu distribution | |
CN111551922B (en) | Three-dimensional space double/multi-base radar high-speed target detection method | |
CN109613507B (en) | Detection method for high-order maneuvering target radar echo | |
CN113504522B (en) | Space-time decoupling and super-resolution angle measurement method based on random switching of transmitting antennas | |
CN111045002A (en) | Maneuvering target coherent accumulation method based on TRT and SNuFFT | |
CN104849708A (en) | High-speed maneuvering target parameter estimation method based on frequency domain polynomial phase transformation | |
Jin et al. | Radar coherent detection for Doppler-ambiguous maneuvering target based on product scaled periodic Lv's distribution | |
Huang et al. | Robust detection and motion parameter estimation for weak maneuvering target in the alpha-stable noise environment | |
CN116449320A (en) | Long-time accumulation and parameter estimation method under frequency agile radar system | |
CN110221264B (en) | Doppler fuzzy target coherent detection method and device based on periodic Keystone transformation | |
CN114545351A (en) | Maneuvering target coherent detection method and system based on range frequency axis inversion transformation and second-order WVD (WVD) | |
Zhou et al. | Coherent integration for targets with constant Cartesian velocities based on accurate range model | |
CN112230220A (en) | Method for detecting dynamic target of Deramp-STAP based on Radon transformation | |
CN107064927A (en) | A kind of high speed missile-borne radar target correlative accumulation method | |
CN115808661A (en) | Distance fuzzy high-speed target accumulation detection method based on remainder processing | |
CN115685169A (en) | Underwater sound weak moving target detection method based on broadband keystone transformation | |
CN113900088A (en) | Long-time coherent accumulation method and system for uniform acceleration maneuvering target | |
Yang et al. | An Efficient Coherent Integration Method for High-Speed Maneuvering Target with Nonlinear Motion |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210706 |