CN109188387A - Distributed coherent radar target component estimation method based on Interpolation compensation - Google Patents
Distributed coherent radar target component estimation method based on Interpolation compensation Download PDFInfo
- Publication number
- CN109188387A CN109188387A CN201811014741.XA CN201811014741A CN109188387A CN 109188387 A CN109188387 A CN 109188387A CN 201811014741 A CN201811014741 A CN 201811014741A CN 109188387 A CN109188387 A CN 109188387A
- Authority
- CN
- China
- Prior art keywords
- value
- radar
- indicate
- matrix
- target
- 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.)
- Granted
Links
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
- 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
-
- 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/414—Discriminating targets with respect to background clutter
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 distributed coherent radar target component estimation method based on Interpolation compensation that the invention discloses a kind of mainly solves the problems, such as that the prior art is not high to the ground poor and real-time treatment effeciency of microinching target component estimated accuracy.Implementation step is: (1) generating base band echo matrix；(2) base band echo matrix is filtered；(3) estimate the rough position parameter of ground microinching target；(4) the fine location parameter of Interpolation compensation method estimation ground microinching target is utilized；(5) abscissa value and ordinate value of ground microinching target are found out.Compared with prior art, the present invention also improving the real-time treatment effeciency of distributed coherent system while improving distributed coherent system to ground microinching target component estimated accuracy.
Description
Technical field
The invention belongs to Radar Technology fields, further relate to one of motion platform Radar Technology field and are based on inserting
The distributed coherent radar target component estimation method that value complement is repaid.The present invention can be using distributed coherent system estimation ground at a slow speed
The position coordinates of moving target.
Background technique
Distributed coherent system clutter when making to marine worker or detecting to high-altitude target is weaker, at this moment can be by weak clutter
Distributed coherent system regards the estimation under noise background to target component as to the estimation of target component under background.Distributed coherent
System has the advantages that detection range is remote and coherent performance is high, thus is widely used in the demanding detection of parameters precision
In scene.
Paper " the distributed full phase parameter radar relevant parameters estimation property that Song Jing, Zhang Jianyun, Zheng Zhidong et al. are delivered at it
Can " (2014,8 electronics of Journal of Electronics&Information Technology and information journal) middle proposition
One kind being based on phase locked distributed radar target component estimation method.This method is in multiple-input and multiple-output (MIMO) mode
Under, derived mixing carat beauty-sieve circle (HCRB) closed solutions of time delay estimation first, then under full coherent mode, to point
After the received target echo of cloth system institute does energy accumulation, the analytic expression of the target output signal-to-noise ratio gain after providing accumulation,
And the configuration guidelines of dual-mode antenna number are studied, it finally obtains under the premise of system phase synchronization accuracy is sufficiently high, is based on phase
Synchronous processing mode can obtain the conclusion of higher time delay estimation precision and output signal-to-noise ratio gain.Existing for this method
Shortcoming is that, since energy accumulation is realized by registration phase, this can bring additional phase to target echo, lead
The target location coordinate estimated is caused to deviate true target location coordinate.
University of Electronic Science and Technology is in a kind of patent document " improved ML folded Clutter in Skywave Radars maneuvering target parameter Estimation side of its application
It is proposed in method " (application number: CN201610190528.9, application publication number: CN105676217A) a kind of improved based on most
The folded Clutter in Skywave Radars maneuvering target method for parameter estimation of maximum-likelihood function ML (Maximum Likelihood).This method is by sky wave thunder
The maneuvering target signal modeling reached is generalized phase multinomial, and it is non-thread that likelihood function maximization problems is then changed into ' overdetermination '
Property least-squares estimation optimization problem, finally propose that carrying out multi-dimensional search in airspace by maximum likelihood function realizes machine
The parameter Estimation of moving-target.Shortcoming existing for this method is that this method has introduced ' overdetermination ' Nonlinear least squares fitting
Optimization problem, this can bring least square fitting error to subsequent maneuvering target parameter Estimation, maneuvering target is caused to be joined
Number estimated accuracy is deteriorated.
BJ University of Aeronautics & Astronautics is in a kind of patent document " moving-target parameter Estimation side based on correlation function of its application
It proposes in method " (application number: CN201510256088.8 application publication number: CN104898119A) a kind of based on correlation function
Moving target parameter estimation method.This method including the following steps: step 1: original transient echo data and correlation are read in
Imaging parameters；Step 2: orientation Fourier transformation processing；Step 3: orientation is the same as compressed sensing CS (Compressive
Sensing) fac-tor compensates；Step 4: distance is handled to Fourier transformation；Step 5: distance is compensated to same distance
Fac-tor carries out Range compress processing；Step 6: distance is handled to inverse Fourier transform；Step 7: distance-Doppler domain
Carry out correlation function processing；Step 8: orientation inverse Fourier transform processing；Step 9: frequency domain zero padding, time domain increase sampling to phase
It closes processing result and carries out interpolation processing, seek maximum value；Step 10: target velocity is estimated by relevant treatment maximum value.This method
Existing shortcoming is, due in frequency domain zero padding step, time domain, which increases sampling and carry out interpolation processing to correlated results to increase, is
Calculation amount causes the real-time of system not high.
Summary of the invention
It is an object of the invention in view of the above shortcomings of the prior art, provide a kind of distributed phase based on Interpolation compensation
Join radar target method for parameter estimation.
Realizing the thinking of the object of the invention is: inside distributed coherent system, multiple transmitting unit radars emit letter simultaneously
Number, the echo after the microinching target scattering of ground is received by all receiving unit radars.In view of in real life can
Operability enables every unit radar inside distributed coherent system, all has transmitting unit radar and receiving unit radar
Effect.And in order to realize reception coherent, it should ensure that all receiving unit radar formation flights.In receiving coherent processing, first
To the three-dimensional echo matrix (range gate number * emits pulse process cycle number * and receives array number) of each receiving unit radar, it is
System error correction, to eliminate the inconsistency bring due to local oscillator phase and time synchronization phase between receiving unit radar
It influences.Down-converted is carried out to the echo matrix after Systematic Error Correction later, i.e., by the center frequency of matrix all elements,
Base band position is moved, the base band echo matrix of each receiving unit radar is obtained.In order in the receiver of each receiving unit radar
In, the echo through target scattering of different transmitting unit radars can be isolated by multiple baseband matched filters, it should be
Each transmitting unit radar selects mutually orthogonal transmitted waveform, to match each transmitting unit in each receiving unit radar
The echo component of radar contribution.In each receiving unit radar, after the corresponding filtering processing of different transmitting unit radars
Three-dimensional echo matrix finds out the rough position parameter of ground microinching target using peak extraction method, including ground is transported at a slow speed
Range gate, Doppler's channel number and the wave beam number of moving-target recycle Interpolation compensation method to obtain ground microinching target later
Fine location parameter, finally by each transmitting unit radar emission, each receiving unit radar receive under, transported at a slow speed with ground
Element value of the fine location parameter of moving-target as the matrix index of the three-dimensional matrice after Interpolation compensation, according to receiving unit thunder
Up to the sequence of serial number from small to large, target echo column vector is formed.Using the column vector, to distributed coherent system earth observation
Every bit in region carries out yield value matching, the corresponding location point of maxgain value is found out, by the corresponding horizontal seat of the location point
Scale value and ordinate value are realized respectively as the abscissa value and ordinate value of ground microinching target and utilize distributed phase
Join the position coordinates of system estimation ground microinching target.In order to improve the real-time treatment effeciency of system and reduce operand, benefit
With the rough position parameter for the target that peak extraction method has obtained, the part near rough position parameter, local area interpolating are only chosen
Three-dimensional echo matrix is compensated, system can be significantly improved under the premise of not reducing Parameter Estimation Precision in this way and handle effect in real time
Rate and reduction operand.
Specific steps of the invention include the following:
(1) base band echo matrix is generated:
(1a) to distributed coherent system inside, by each receiving unit radar receive ground microinching target scattering after
The echo matrix that echo data is constituted carries out Systematic Error Correction, the echo matrix after obtaining Systematic Error Correction；
The center frequency of element at each position of echo matrix after Systematic Error Correction is moved base by (1b)
Band frequency location, obtains the base band echo matrix of each receiving unit radar；
(2) base band echo matrix is filtered:
(2a) utilizes distance domain matched filtering formula, calculates each element in base band echo matrix and passes through distance domain matched filtering
After processing, corresponding to the echo values of each transmitting unit radar, by all numerical value according to range gate where it from small to large into
Row sequence, the three-dimensional echo matrix after forming distance domain matched filtering；
(2b) utilizes doppler filtering formula, and each element is passed through in the three-dimensional echo matrix after calculating distance domain matched filtering
All numerical value are ranked up from small to large according to Doppler's channel number where it, form by doppler filtering treated numerical value
Three-dimensional echo matrix after doppler filtering；
(2c) utilizes digital beam froming formula, and each element passes through number in the three-dimensional echo matrix after calculating doppler filtering
All numerical value are ranked up, composition number by word Wave beam forming treated numerical value from small to large according to wave beam serial number where it
Three-dimensional echo matrix after Wave beam forming；
(3) estimate the rough position parameter of ground microinching target；
(4) the fine location parameter of Interpolation compensation method estimation ground microinching target is utilized:
(4a) utilizes Interpolation compensation method, the three-dimensional echo matrix after obtaining Interpolation compensation；
(4b) utilizes peak extraction method, and ground microinching target is found out from the three-dimensional echo matrix after Interpolation compensation
Fine location parameter；
(5) abscissa value and ordinate value of ground microinching target are found out:
(5a) utilizes assignment method, obtains the abscissa value of ground microinching target and the region of search of ordinate value；
(5b) it is slow to find out ground from rectangular area composed by the region of search as abscissa value and ordinate value respectively
The abscissa value and ordinate value of fast moving target.
Compared with the prior art, the present invention has the following advantages:
First, since present invention utilizes Interpolation compensation methods, to the phase additional to target echo bring by energy accumulation
Position is compensated, and is overcome the prior art to since energy accumulation brings additional phase to target echo, is caused to estimate
Target location coordinate deviate the deficiency of true target location coordinate so that the present invention is not susceptible to additionally in engineering practice
Phase influence, improve distributed coherent system to the estimated accuracy of ground microinching target position parameter.
Second, since present invention utilizes peak extraction methods, it is slow to have obtained ground by the corresponding index of extraction maximum value
The location parameter of fast moving target improves the real-time of system, overcomes the prior art due to increasing sampling to phase using time domain
It closes result and carries out interpolation processing, increase system-computed amount, the deficiency for causing the real-time of system not high allows the invention to mention
The real-time treatment effeciency of high distribution coherent system.
Detailed description of the invention
Fig. 1 is flow chart of the invention；
Fig. 2 is analogous diagram of the invention.
Specific embodiment
The present invention will be further described with reference to the accompanying drawing.
Referring to Fig.1, step of the invention is further described.
Step 1, base band echo matrix is generated.
Echo inside distributed coherent system, after microinching target scattering in ground is received by each receiving unit radar
The echo matrix that data are constituted carries out Systematic Error Correction, the echo matrix after obtaining Systematic Error Correction.
Specific step is as follows for the Systematic Error Correction:
Step 1, using the systematic error estimation module inside distributed coherent system, estimate each receiving unit radar it
Between, due to the otherness bring systematic error matrix of local oscillator phase and time synchronization phase.
Step 2 is returned with what the echo data after receiving ground microinching target scattering by each receiving unit radar was constituted
Element in wave matrix at each position obtains systematic error divided by the element at the same position in systematic error matrix
Echo matrix after correction.
By the center frequency of the element at each position of the echo matrix after Systematic Error Correction, base band frequency is moved
Rate position obtains the base band echo matrix of each receiving unit radar.
Step 2, base band echo matrix is filtered.
Using distance domain matched filtering formula, calculates each element in base band echo matrix and handled by distance domain matched filtering
Afterwards, the echo values corresponding to each transmitting unit radar arrange all numerical value according to range gate where it from small to large
Sequence, the three-dimensional echo matrix after forming distance domain matched filtering.
The distance domain matched filtering formula is as follows:
Wherein, yp,q(l, k, n) indicates base band echo matrix by distance domain matched filtering treated p-th of transmitting singly
First radar emission, received first of the range gate of q-th of receiving unit radar, k-th of transmitting pulse process cycle, n-th of reception
The numerical value of array element, p=1,2 ..., mT, mTThe sum of transmitting unit radar in the distributed coherent system of expression, q=1,
2,...,mR, mRIndicate the sum of the receiving unit radar in distributed coherent system, l=1,2 ..., L, L indicate range gate
Sum, k=1,2 ..., K, K indicate the sum of all transmitting pulse process cycles, and n=1,2 ..., N, N indicate distributed phase
The sum of reception array element in each receiving unit radar of ginseng system, ∑ indicate sum operation, xq(r, k, n) indicates q-th of reception
Received r-th of unit radar apart from frequency domain sample point, k-th of transmitting pulse process cycle, the distance frequency for receiving array element for n-th
Domain echo data, sp(r) r-th of filtering at frequency domain sample point apart from frequency matching of p-th of transmitting unit radar is indicated
Wave number evidence, * indicate conjugate operation, and exp is indicated using natural logrithm as the index operation at bottom, and j indicates that imaginary unit's symbol, π indicate
Pi.
Using doppler filtering formula, each element is through excessive general in the three-dimensional echo matrix after calculating distance domain matched filtering
All numerical value are ranked up from small to large according to Doppler's channel number where it, how general form by the numerical value after strangling filtering processing
Strangle filtered three-dimensional echo matrix.
The doppler filtering formula is as follows:
Wherein, zp,q(l, a, n) indicates the three-dimensional echo matrix after distance domain matched filtering after doppler filtering is handled
P-th of transmitting unit radar emission, received first of the range gate of q-th of receiving unit radar, a-th of Doppler channel,
The numerical value of n reception array element, a=1,2 ..., K, b indicate the serial number of transmitting pulse process cycle.
Using digital beam froming formula, each element passes through digital wave in the three-dimensional echo matrix after calculating doppler filtering
Beam forms treated numerical value, and all numerical value are ranked up from small to large according to wave beam serial number where it, form digital beam
Three-dimensional echo matrix after formation.
The digital beam froming formula is as follows:
Wherein, fp,qIn three-dimensional echo matrix after (l, a, c) expression doppler filtering, handled by digital beam froming
Rear p-th of transmitting unit radar emission, received first of the range gate of q-th of receiving unit radar, a-th of Doppler channel,
The numerical value of c-th of wave beam, c=1,2 ..., N, g indicate to receive the serial number of array element.
Step 3, estimate the rough position parameter of ground microinching target.
Specific step is as follows for the rough position parameter of the estimation ground microinching target:
Step 1 finds out the maximum value of all elements from the three-dimensional echo matrix after digital beam froming.
Step 2 replaces the rough position parameter of ground microinching target with the matrix index of the maximum value.
Step 4, the fine location parameter of Interpolation compensation method estimation ground microinching target is utilized.
Three-dimensional echo matrix using Interpolation compensation method, after obtaining Interpolation compensation.
Specific step is as follows for the Interpolation compensation method:
Step 1 in the three-dimensional echo matrix after digital beam froming, does inverse fast Fourier transform to each dimension respectively,
Three-dimensional echo matrix after obtaining inverse fast Fourier transform.
Step 2 distinguishes zero padding to each dimension tail portion of the three-dimensional echo matrix after inverse fast Fourier transform, until three of them
Dimension respectively reach α L,Until α N, α indicates zero padding multiple, the positive integer in α ∈ [1,10], and ∈ expression belongs to symbol,
Three-dimensional echo matrix after obtaining zero padding.
Step 3, to each dimension of the three-dimensional echo matrix after zero padding, do respectively α L point,The fast Fourier of point and α N point
Transformation, the three-dimensional echo matrix after obtaining Interpolation compensation.
Using peak extraction method, the fine of ground microinching target is found out from the three-dimensional echo matrix after Interpolation compensation
Location parameter.
Specific step is as follows for the peak extraction method:
Step 1 finds out the maximum value of all elements from the three-dimensional echo matrix after Interpolation compensation.
Step 2 replaces the fine location parameter of ground microinching target with the matrix index of the maximum value.
Step 5, the abscissa value and ordinate value of ground microinching target are found out.
Using assignment method, the abscissa value of ground microinching target and the region of search of ordinate value are obtained.
Specific step is as follows for the assignment method:
Step 1 replaces ground microinching with the value interval of distributed coherent system earth observation region abscissa value
The region of search of target lateral coordinates value.
Step 2 replaces ground microinching with the value interval of distributed coherent system earth observation region ordinate value
The region of search of target ordinate value.
From rectangular area composed by the region of search as abscissa value and ordinate value, ground is found out respectively and is transported at a slow speed
The abscissa value and ordinate value of moving-target.
Specific step is as follows for the abscissa value for finding out ground microinching target respectively and ordinate value:
Step 1 finds out the fine location ginseng with ground microinching target from the three-dimensional echo matrix after Interpolation compensation
Number is the element value of matrix index, and the sequence by these element values according to receiving unit radar serial number from small to large forms target
Echo column vector.
Step 2 calculates the corresponding search steering vector of each location point in rectangular area according to the following formula:
Wherein, swIndicate that the corresponding search steering vector of w-th of location point in rectangular area, λ indicate launch wavelength, Tw1Table
Show distance of w-th of location point to the 1st transmitting unit radar, R in rectangular areaw1Indicate w-th of location point in rectangular area
To the distance of the 1st receiving unit radar, Tw2Indicate w-th of location point in rectangular area to the 2nd transmitting unit radar away from
From TwpIndicate distance of w-th of location point to p-th of transmitting unit radar in rectangular area, RwqIndicate w in rectangular area
Distance of a location point to q-th of receiving unit radar, T expression transposition operation.
Step 3 calculates the yield value in rectangular area at each location point according to the following formula:
Yw=sw Hz
Wherein, YwIndicate that the yield value in rectangular area at w-th of location point, H indicate that conjugate transposition operation, z indicate mesh
Mark echo column vector.
Step 4 finds out the corresponding location point of maxgain value from all location points in rectangular area, by the location point
Corresponding abscissa value and ordinate value, respectively as the abscissa value and ordinate value of ground microinching target.
Effect of the invention is described further below with reference to emulation experiment.
1. simulated conditions:
The environment of emulation experiment of the present invention are as follows: MATLAB 2017b, Intel (R) Xeon (R) CPU 2.20GHz, Window
7 professional versions.
2. emulation content and interpretation of result:
Emulation experiment of the invention is to be transported at a slow speed using method of the invention according to the received ground of distributed coherent system
The echo of moving-target estimates the abscissa value and ordinate value of ground microinching target.In distributed coherent system
Every unit radar have transceiver, i.e., each unit radar is both transmitting unit radar and receiving unit radar.
The sum of unit radar is 4, and the sum of range gate is 200, and the sum of all transmitting pulse process cycles is 128, respectively
The reception array number of receiving unit radar is 8, and zero padding multiple is 10 times, and the transmission power of each transmitting unit radar is
200kw, transmitting carrier frequency are 300MHz, and transmitted signal bandwidth is 1MHz, and pulse recurrence frequency is 5KHz, ground microinching target
Echo signal-to-noise ratio be 30dB.
Fig. 2 is analogous diagram of the invention.Fig. 2 (a) be for ground microinching target, between receiving unit radar between
When away from being 100 meters, ideally, in the case of Interpolation compensation and without the echo vector phase in the case of Interpolation compensation compare
Figure.Abscissa in Fig. 2 (a) indicates that send-receive pair, ordinate indicate the phase of echo vector.Without mark in Fig. 2 (a)
Curve indicates the simulation result curve of echo vector phase ideally, indicates Interpolation compensation with the curve of square mark
In the case of echo vector phase simulation result curve, in the case of indicating no Interpolation compensation with the curve that triangle indicates
The simulation result curve of echo vector phase.
When signal-to-noise ratio is that 30dB is compared when spacing is 100 meters between receiving unit radar it can be seen from Fig. 2 (a)
Echo vector phase curve in the case of no Interpolation compensation, echo vector phase curve in the case of Interpolation compensation and ideal
In the case of the goodness of fit largely put of echo vector phase curve it is more preferable, only at the position that abscissa is 11, due to sky
Between direction ambiguity cause the phase different from ideally corresponding points phase.It can be seen that between receiving unit radar
When spacing is uniform, using the echo vector in the case of Interpolation compensation, target lateral coordinates value and ordinate value can be more accurately estimated.
Echo vector phase ideally, only with transmitting range and reception distance dependent.Transmitting range refers to, each transmitting unit
Distance of the radar to ground microinching target.It receives distance to refer to, each receiving unit radar to ground microinching target
Distance.
Fig. 2 (b) be for ground microinching target, between four receiving unit radars spacing be respectively 100 meters,
At 200 meters and 300 meters, phase ideally, in the case of Interpolation compensation and without the echo vector in the case of Interpolation compensation
Compare figure.Abscissa indicates that send-receive pair, ordinate indicate echo vector phase in Fig. 2 (b).Without mark in Fig. 2 (b)
Curve indicates the simulation result curve of echo vector phase ideally, indicates Interpolation compensation feelings with the curve that diamond shape indicates
The simulation result curve of echo vector phase under condition indicates the echo in the case of no Interpolation compensation with the curve of circle mark
The simulation result curve of vector phase.
Receiving unit radar spacing parameter is adjusted it can be seen from Fig. 2 (b), so that four receiving unit radars
Between spacing when being respectively 100 meters, 200 meters and 300 meters, it is bent compared to the echo vector phase in the case of no Interpolation compensation
Line, the every bit of echo vector phase curve and echo vector phase curve ideally in the case of Interpolation compensation, kiss
It is right all more preferable.It can be seen that being sweared when spacing is non-homogeneous between receiving unit radar using the echo in the case of Interpolation compensation
Amount, can accurately estimate target lateral coordinates value and ordinate value.
Claims (10)
1. a kind of distributed coherent radar target component estimation method based on Interpolation compensation, which is characterized in that by each reception
Unit radar receives the echo matrix of the echo data composition after ground microinching target scattering, and Interpolation compensation method is utilized to estimate
The fine location parameter of ground microinching target finds out the abscissa value and ordinate value of ground microinching target；The party
The specific steps of method include the following:
(1) base band echo matrix is generated:
(1a) to distributed coherent system inside, by each receiving unit radar receive ground microinching target scattering after echo
The echo matrix that data are constituted carries out Systematic Error Correction, the echo matrix after obtaining Systematic Error Correction；
The center frequency of element at each position of echo matrix after Systematic Error Correction is moved base band frequency by (1b)
Rate position obtains the base band echo matrix of each receiving unit radar；
(2) base band echo matrix is filtered:
(2a) utilizes distance domain matched filtering formula, calculates each element in base band echo matrix and handles by distance domain matched filtering
Afterwards, the echo values corresponding to each transmitting unit radar arrange all numerical value according to range gate where it from small to large
Sequence, the three-dimensional echo matrix after forming distance domain matched filtering；
(2b) utilizes doppler filtering formula, and each element is through excessive general in the three-dimensional echo matrix after calculating distance domain matched filtering
All numerical value are ranked up from small to large according to Doppler's channel number where it, how general form by the numerical value after strangling filtering processing
Strangle filtered three-dimensional echo matrix；
(2c) utilizes digital beam froming formula, and each element passes through digital wave in the three-dimensional echo matrix after calculating doppler filtering
Beam forms treated numerical value, and all numerical value are ranked up from small to large according to wave beam serial number where it, form digital beam
Three-dimensional echo matrix after formation；
(3) estimate the rough position parameter of ground microinching target；
(4) the fine location parameter of Interpolation compensation method estimation ground microinching target is utilized:
(4a) utilizes Interpolation compensation method, the three-dimensional echo matrix after obtaining Interpolation compensation；
(4b) utilizes peak extraction method, and the fine of ground microinching target is found out from the three-dimensional echo matrix after Interpolation compensation
Location parameter；
(5) abscissa value and ordinate value of ground microinching target are found out:
(5a) utilizes assignment method, obtains the abscissa value of ground microinching target and the region of search of ordinate value；
(5b) finds out ground respectively and transports at a slow speed from rectangular area composed by the region of search as abscissa value and ordinate value
The abscissa value and ordinate value of moving-target.
2. the distributed coherent radar target component estimation method according to claim 1 based on Interpolation compensation, feature
It is, specific step is as follows for Systematic Error Correction described in step (1a):
The first step estimated between each receiving unit radar using the systematic error estimation module inside distributed coherent system, by
In the otherness bring systematic error matrix of local oscillator phase and time synchronization phase；
Second step, the echo square constituted with the echo data after receiving ground microinching target scattering by each receiving unit radar
Element in battle array at each position obtains Systematic Error Correction divided by the element at the same position in systematic error matrix
Echo matrix afterwards.
3. the distributed coherent radar target component estimation method according to claim 1 based on Interpolation compensation, feature
It is, distance domain matched filtering formula described in step (2a) is as follows:
Wherein, yp,q(l, k, n) indicates base band echo matrix by distance domain matched filtering treated p-th of transmitting unit thunder
Up to transmitting, received first of the range gate of q-th of receiving unit radar, k-th of transmitting pulse process cycle, n-th of reception array element
Numerical value, p=1,2 ..., mT, mTIndicate the sum of the transmitting unit radar in distributed coherent system, q=1,2 ..., mR,
mRIndicate the sum of the receiving unit radar in distributed coherent system, l=1,2 ..., L, L indicate the sum of range gate, k=
1,2 ..., K, K indicate that the sum of all transmitting pulse process cycles, n=1,2 ..., N, N indicate that distributed coherent system is each
The sum of reception array element in receiving unit radar, ∑ indicate sum operation, xq(r, k, n) indicates q-th of receiving unit radar
Received r-th apart from frequency domain sample point, k-th transmitting pulse process cycle, n-th reception array element apart from frequency domain echo number
According to sp(r) indicate r-th of p-th of transmitting unit radar at frequency domain sample point apart from frequency matching filtering data, *
Indicate conjugate operation, exp is indicated using natural logrithm as the index operation at bottom, and j indicates that imaginary unit's symbol, π indicate pi.
4. the distributed coherent radar target component estimation method according to claim 1 based on Interpolation compensation, feature
It is, doppler filtering formula described in step (2b) is as follows:
Wherein, zp,q(l, a, n) indicates the three-dimensional echo matrix after distance domain matched filtering, and by doppler filtering, treated the
P transmitting unit radar emission, received first of the range gate of q-th of receiving unit radar, a-th of Doppler channel, n-th
The numerical value of array element, a=1 are received, 2 ..., K, b indicate the serial number of transmitting pulse process cycle.
5. the distributed coherent radar target component estimation method according to claim 1 based on Interpolation compensation, feature
It is, digital beam froming formula described in step (2c) is as follows:
Wherein, fp,qIn three-dimensional echo matrix after (l, a, c) expression doppler filtering, by digital beam froming, treated
P-th of transmitting unit radar emission, received first of the range gate of q-th of receiving unit radar, a-th Doppler channel, c
The numerical value of a wave beam, c=1,2 ..., N, g indicate to receive the serial number of array element.
6. the distributed coherent radar target component estimation method according to claim 1 based on Interpolation compensation, feature
It is, specific step is as follows for the rough position parameter of estimation ground microinching target described in step (3):
The first step finds out the maximum value of all elements from the three-dimensional echo matrix after digital beam froming；
Second step replaces the rough position parameter of ground microinching target with the matrix index of the maximum value.
7. the distributed coherent radar target component estimation method according to claim 1 based on Interpolation compensation, feature
It is, specific step is as follows for Interpolation compensation method described in step (4a):
The first step in the three-dimensional echo matrix after digital beam froming, is done inverse fast Fourier transform to each dimension respectively, is obtained
Three-dimensional echo matrix after inverse fast Fourier transform；
Second step distinguishes zero padding to each dimension tail portion of the three-dimensional echo matrix after inverse fast Fourier transform, until three of them dimension
Number respectively reach α L,Until α N, α indicates zero padding multiple, the positive integer in α ∈ [1,10], and ∈ expression belongs to symbol, obtains
Three-dimensional echo matrix after to zero padding；
Third step, to each dimension of the three-dimensional echo matrix after zero padding, do respectively α L point,The fast Fourier of point and α N point becomes
It changes, the three-dimensional echo matrix after obtaining Interpolation compensation.
8. the distributed coherent radar target component estimation method according to claim 1 based on Interpolation compensation, feature
It is, specific step is as follows for peak extraction method described in step (4b):
The first step finds out the maximum value of all elements from the three-dimensional echo matrix after Interpolation compensation；
Second step replaces the fine location parameter of ground microinching target with the matrix index of the maximum value.
9. the distributed coherent radar target component estimation method according to claim 1 based on Interpolation compensation, feature
It is, specific step is as follows for assignment method described in step (5a):
The first step replaces ground microinching mesh with the value interval of distributed coherent system earth observation region abscissa value
Mark the region of search of abscissa value；
Second step replaces ground microinching mesh with the value interval of distributed coherent system earth observation region ordinate value
Mark the region of search of ordinate value.
10. the distributed coherent radar target component estimation method according to claim 1 based on Interpolation compensation, feature
It is, the specific steps of the abscissa value and ordinate value of finding out ground microinching target described in step (5b) respectively are such as
Under:
The first step is found out from the three-dimensional echo matrix after Interpolation compensation with the fine location parameter of ground microinching target
For the element value of matrix index, sequence by these element values according to receiving unit radar serial number from small to large, composition echo arrow
Amount；
Second step calculates the corresponding search steering vector of each location point in rectangular area according to the following formula:
Wherein, swIndicate that the corresponding search steering vector of w-th of location point in rectangular area, λ indicate launch wavelength, Tw1Indicate square
Distance of w-th of location point to the 1st transmitting unit radar, R in shape regionw1Indicate that w-th of location point is to the 1st in rectangular area
The distance of a receiving unit radar, Tw2Indicate distance of w-th of location point to the 2nd transmitting unit radar in rectangular area, Twp
Indicate distance of w-th of location point to p-th of transmitting unit radar in rectangular area, RwqIndicate w-th of position in rectangular area
O'clock to q-th of receiving unit radar distance, T indicate transposition operation；
Third step calculates the yield value in rectangular area at each location point according to the following formula:
Yw=sw Hz
Wherein, YwIndicate that the yield value in rectangular area at w-th of location point, H indicate that conjugate transposition operation, z indicate echo arrow
Amount；
4th step finds out the corresponding location point of maxgain value from all location points in rectangular area, by the location point pair
The abscissa value and ordinate value answered, respectively as the abscissa value and ordinate value of ground microinching target.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811014741.XA CN109188387B (en) | 2018-08-31 | 2018-08-31 | Target parameter estimation method for distributed coherent radar based on interpolation compensation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811014741.XA CN109188387B (en) | 2018-08-31 | 2018-08-31 | Target parameter estimation method for distributed coherent radar based on interpolation compensation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109188387A true CN109188387A (en) | 2019-01-11 |
| CN109188387B CN109188387B (en) | 2022-12-02 |
Family
ID=64917461
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811014741.XA Active CN109188387B (en) | 2018-08-31 | 2018-08-31 | Target parameter estimation method for distributed coherent radar based on interpolation compensation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109188387B (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110703211A (en) * | 2019-08-29 | 2020-01-17 | 西安电子科技大学 | Distributed coherent radar phase error estimation method and system based on polynomial iteration |
| CN110907910A (en) * | 2019-11-27 | 2020-03-24 | 中国船舶重工集团公司第七二四研究所 | Distributed coherent radar moving target echo coherent synthesis method |
| CN111884981A (en) * | 2020-08-04 | 2020-11-03 | 北京航空航天大学 | Signal processing method suitable for high-sensitivity satellite-borne ADS-B receiver |
| WO2020237663A1 (en) * | 2019-05-31 | 2020-12-03 | 深圳市大疆创新科技有限公司 | Multi-channel lidar point cloud interpolation method and ranging apparatus |
| CN112305526A (en) * | 2020-10-22 | 2021-02-02 | 电子科技大学 | Distributed array system synchronization method based on external calibration source |
| CN112505667A (en) * | 2020-11-19 | 2021-03-16 | 哈尔滨工程大学 | Two-dimensional sonar array motion attitude self-calibration method |
| CN112816945A (en) * | 2021-01-04 | 2021-05-18 | 北京无线电测量研究所 | Method and system for calculating out-of-plane distributed phased array radar coherent synthetic gain |
| CN113050053A (en) * | 2021-03-17 | 2021-06-29 | 中国人民解放军国防科技大学 | Method and system for acquiring coherent parameters of distributed coherent radar of moving platform |
| CN116203519A (en) * | 2023-05-05 | 2023-06-02 | 中国电子科技集团公司信息科学研究院 | Error calibration method for distributed radar system transceiver channel distance system |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100259442A1 (en) * | 2009-04-13 | 2010-10-14 | Abatzoglou Theagenis J | Fast implementation of a maximum likelihood algorithm for the estimation of target motion parameters |
| CN103424741A (en) * | 2013-08-29 | 2013-12-04 | 西安电子科技大学 | Smooth procession cone parameter estimation method based on high-resolution ISAR imaging |
| CN103472449A (en) * | 2013-09-16 | 2013-12-25 | 电子科技大学 | BP wideband synthesis method based on MIMO image domain |
-
2018
- 2018-08-31 CN CN201811014741.XA patent/CN109188387B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100259442A1 (en) * | 2009-04-13 | 2010-10-14 | Abatzoglou Theagenis J | Fast implementation of a maximum likelihood algorithm for the estimation of target motion parameters |
| CN103424741A (en) * | 2013-08-29 | 2013-12-04 | 西安电子科技大学 | Smooth procession cone parameter estimation method based on high-resolution ISAR imaging |
| CN103472449A (en) * | 2013-09-16 | 2013-12-25 | 电子科技大学 | BP wideband synthesis method based on MIMO image domain |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020237663A1 (en) * | 2019-05-31 | 2020-12-03 | 深圳市大疆创新科技有限公司 | Multi-channel lidar point cloud interpolation method and ranging apparatus |
| CN110703211A (en) * | 2019-08-29 | 2020-01-17 | 西安电子科技大学 | Distributed coherent radar phase error estimation method and system based on polynomial iteration |
| CN110703211B (en) * | 2019-08-29 | 2021-11-19 | 西安电子科技大学 | Distributed coherent radar phase error estimation method and system based on polynomial iteration |
| CN110907910A (en) * | 2019-11-27 | 2020-03-24 | 中国船舶重工集团公司第七二四研究所 | Distributed coherent radar moving target echo coherent synthesis method |
| CN110907910B (en) * | 2019-11-27 | 2022-06-03 | 中国船舶重工集团公司第七二四研究所 | Distributed coherent radar moving target echo coherent synthesis method |
| CN111884981B (en) * | 2020-08-04 | 2021-05-18 | 北京航空航天大学 | Signal processing method suitable for high-sensitivity satellite-borne ADS-B receiver |
| CN111884981A (en) * | 2020-08-04 | 2020-11-03 | 北京航空航天大学 | Signal processing method suitable for high-sensitivity satellite-borne ADS-B receiver |
| CN112305526B (en) * | 2020-10-22 | 2023-05-26 | 电子科技大学 | Distributed array system synchronization method based on external calibration source |
| CN112305526A (en) * | 2020-10-22 | 2021-02-02 | 电子科技大学 | Distributed array system synchronization method based on external calibration source |
| CN112505667A (en) * | 2020-11-19 | 2021-03-16 | 哈尔滨工程大学 | Two-dimensional sonar array motion attitude self-calibration method |
| CN112505667B (en) * | 2020-11-19 | 2022-07-15 | 哈尔滨工程大学 | Two-dimensional sonar array motion attitude self-calibration method |
| CN112816945A (en) * | 2021-01-04 | 2021-05-18 | 北京无线电测量研究所 | Method and system for calculating out-of-plane distributed phased array radar coherent synthetic gain |
| CN112816945B (en) * | 2021-01-04 | 2024-05-14 | 北京无线电测量研究所 | Method and system for calculating composite gain of different-surface distributed phased array Lei Daxiang |
| CN113050053A (en) * | 2021-03-17 | 2021-06-29 | 中国人民解放军国防科技大学 | Method and system for acquiring coherent parameters of distributed coherent radar of moving platform |
| CN113050053B (en) * | 2021-03-17 | 2023-07-28 | 中国人民解放军国防科技大学 | Method and system for acquiring phase parameters of distributed phase-coherent radar of mobile platform |
| CN116203519A (en) * | 2023-05-05 | 2023-06-02 | 中国电子科技集团公司信息科学研究院 | Error calibration method for distributed radar system transceiver channel distance system |
| CN116203519B (en) * | 2023-05-05 | 2023-06-27 | 中国电子科技集团公司信息科学研究院 | Error calibration method for distributed radar system transceiver channel distance system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109188387B (en) | 2022-12-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109188387A (en) | Distributed coherent radar target component estimation method based on Interpolation compensation | |
| CN106970371B (en) | A kind of object detection method based on Keystone and matched filtering | |
| CN106353744B (en) | Multi-parameter combined estimation method based on bistatic FDA-MIMO radars | |
| CN110412558A (en) | The vehicle-mounted fmcw radar velocity ambiguity method of solution based on TDM MIMO | |
| CN103901429B (en) | Maneuvering target inverse synthetic aperture radar imaging method based on sparse aperture | |
| CN108459307B (en) | Clutter-based MIMO radar transmit-receive array amplitude-phase error correction method | |
| CN110058202B (en) | High-performance hybrid positioning method based on angle and signal arrival time difference estimation and implementation device thereof | |
| CN103823217A (en) | Bistatic MIMO radar high-speed movement target parameter estimation method based on dual-frequency transmission | |
| CN108957387B (en) | Method and system for estimating two-dimensional arrival angle of satellite signal | |
| CN107255822A (en) | GNSS receiver modulated parameter estimating method under multi-path environment | |
| CN105717496A (en) | Realization method of FDA (Frequency Diverse Array) MIMO (Multiple-Input Multiple-Output) radar system based on matrix completion | |
| CN103383449A (en) | ESPRIT algorithm based short-range clutter suppression method for airborne radar | |
| CN108828504A (en) | MIMO radar target direction method for quick estimating based on part waveform correlation | |
| Yu et al. | Ground moving target motion parameter estimation using Radon modified Lv's distribution | |
| CN110109098B (en) | Scanning radar rapid super-resolution imaging method | |
| CN106125075B (en) | A kind of motion error extraction method of bistatic forward sight synthetic aperture radar | |
| CN109001671B (en) | Target detection and parameter estimation method and device for frequency hopping signal | |
| CN111198387A (en) | Space-time sampling navigation positioning method capable of resisting deception jamming | |
| CN109061599B (en) | STAP method based on cyclostationarity and symmetric prior knowledge | |
| CN111812644B (en) | MIMO radar imaging method based on sparse estimation | |
| US11428772B2 (en) | Compensating for crosstalk in determination of an angle of arrival of an electromagnetic wave at a receive antenna | |
| CN105785333A (en) | Airborne MIMO radar robust dimension-reduction space-time self-adaptive processing method | |
| CN111208511A (en) | Super-surface antenna three-dimensional imaging method based on linear frequency modulation signals | |
| CN106443623B (en) | A kind of sky-wave OTH radar target and Ionospheric Parameters combined estimation method | |
| CN108957394A (en) | A kind of day earthwave delay time estimation method applied to Loran |
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 |