CN103257344B - Iteration-adaptive-algorithm-based method for detecting coherent MIMO radar target - Google Patents
Iteration-adaptive-algorithm-based method for detecting coherent MIMO radar target Download PDFInfo
- Publication number
- CN103257344B CN103257344B CN201310162479.4A CN201310162479A CN103257344B CN 103257344 B CN103257344 B CN 103257344B CN 201310162479 A CN201310162479 A CN 201310162479A CN 103257344 B CN103257344 B CN 103257344B
- Authority
- CN
- China
- Prior art keywords
- target
- centerdot
- represent
- reflection coefficient
- iteration
- 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
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses an iteration-adaptive-algorithm-based method for detecting a coherent MIMO radar target. First, a relation model among a first echo signal, a transmitting signal reflection coefficients and noise interference of the target is confirmed, wherein the relation model is observed on the n<th> observation point, secondly, the relation model is linearized, further the reflection coefficients of the target are initialized through a delayed superposition DAS algorithm, and finally, a bayesian model order selecting tool is used for optimizing the reflection coefficients of the target. An iteration adaptive algorithm can effectively solve existing encountered problems of moving target parameter detection and greatly improves accuracy of detecting of the coherent MIMO radar moving target.
Description
Technical field
The present invention relates to a kind of MIMO radar object detection method, is more particularly a kind of relevant MIMO radar object detection method based on iteration self-adapting algorithm.
Background technology
Multiple-input and multiple-output (MIMO) radar is compared phased-array radar and is had detection perform, good anti-stealthy effect and the good anti-lethality higher to weak signal target, has attracted increasing scholar to study it.MIMO radar generally can be divided into two classes: the first kind is statistics MIMO radar or incoherent MIMO radar, the i.e. bare cloth of bay space, thus can obtain space diversity gain, effectively overcomes the impact that target RCS rises and falls on Studies of Radar Detection performance; Equations of The Second Kind is relevant MIMO, adopts less aerial array just can form very large virtual array aperture, improves radar angular resolution and interference rejection capability.
Space-time adaptive (STAP) is a kind of moving target detection method, has clutter suppression capability strong, the advantage that robustness is good, interference rejection capability is strong.Clutter suppression capability is by force because STAP technology has space-time two-dimensional filtering characteristic.Robustness is well because STAP technology, by adaptive characteristic, can realize mating with the effective of complicated external environment; Meanwhile, it also can compensate multiple inevitable systematic error to a certain extent.Interference rejection capability is by force because STAP can realize the suppression to the multiple interference under complex electromagnetic environment.STAP technology is combined with relevant MIMO radar, not only system can be made to have the superiority of MIMO performance but also have the superiority of STAP, better raising is to clutter recognition performance, thus the radar system that can make obtains the ability of higher target detection, the tracking performance of moving-target, the performance of Faint target detection and stronger signal interception.
No matter for moving-target or quiet target, as long as the reflection coefficient obtaining target then can know the location parameter (complex magnitude, distance, angle) of the location parameter (complex magnitude, distance, angle, Doppler) of moving-target, quiet target.
In prior art, relevant MIMO radar space-time adaptive processing method sharply increases because suppression jamming causes Clutter Degrees of Freedom under the electromagnetic environment of complexity, and the degree of freedom of traditional self-adaptive processing algorithm (STAP) cannot meet clutter recognition requirement; Duplicity interference causes the seriously non-homogeneous of the distribution of clutter, makes STAP algorithm cannot obtain enough training samples meeting independent same distribution condition.In addition, existing STAP needs auxiliary data to go to estimate the master data of clutter and noise harmonics matrix, but the auxiliary data with high accuracy is very un-come-at-able, especially for clutter environment heterogeneous.In order to reduce the dependence for auxiliary data, propose the space-time adaptive algorithm not needing auxiliary data, such as DAS method, but the problem of the low high secondary lobe of resolution is led in this class methods face.Therefore this patent is intended proposing innovatory algorithm on the basis of existing space-time adaptive method, overcomes the shortcoming of existing STAP and DAS, better plays the advantage of relevant MIMO radar and space-time adaptive method.
Summary of the invention
The object of the invention is to propose a kind of relevant MIMO radar moving target detection method calculated based on iteration self-adapting, the method obtains its location parameter by being optimized the reflection coefficient of moving-target and quiet target respectively, is specially:
A, when observed object is quiet target, do not consider Doppler's effect
A1, determine first echoed signal that the n-th observation station observes and transmit, between the reflection coefficient of target and noise, relational model is as follows:
In formula, Y (n) is first echoed signal that the n-th observation station observes, wherein
for the number of observation station; Y
hn Y (n) is carried out conjugate transpose operation by () expression; P represents along the minimum range unit i.e. number of distance quantifying unit that directions of rays is divided into by distance in radar echo signal process, r=1,2 ... P; K represents that the angle intervals that there is target quantizes number of unit, a=1 ... ..K; α
r,arepresent the reflection coefficient of target region internal object; b
an () is that the receiving array of the n-th observation station guides vector; a
an () is that the emission array of the n-th observation station guides vector;
represent a
an () does transpose operation; E (n) represents the additive noise of the n-th observation station Received signal strength, E
hn E (n) is carried out conjugate transpose operation by () expression;
representing will
carry out conjugate transpose operation,
For the waveform matrix that transmits obtained after zero padding process,
For the waveform matrix that transmits, wherein, s
mrepresent the orthogonal waveforms of m article of transmission antennas transmit, m=1 ... M
t; M
tfor the number of emitting antenna, L is the number of subpulse of transmitting;
J
r∈ C
(L+P-1) × (L+P-1)for transformation matrix, be used for describing different distance and quantize the signal that receives of section, its expression formula is:
A2, to above-mentioned relation model linearization;
A3, time delay is utilized to superpose the reflection coefficient of DAS algorithm initialization target, as follows:
Wherein:
r=1,...P,a=1,...K;
Wherein
represent initialized target reflection factor,
represent V
r,acarry out conjugate transpose operation; M
rrepresent the number of receiving antenna; Vec represents matrix is carried out vectorization operation, and arranging successively by the order of rectangular array by matrix is row;
A4, Bayesian model exponent number selection tool is utilized to carry out the reflection coefficient of optimization aim, be specially: time delay is superposed the target reflection factor that DAS is initial, (1) formula of bringing into obtains the Bayes BIC value of initial value, the reflection coefficient of iterative target also upgrades Bayes BIC value, twice Bayes BIC value relatively, when Bayes BIC value no longer reduces, then iteration is stopped to obtain the reflection coefficient of target;
In formula, || .||
2represent Euclid norm,
represent Euclid norm square, η represents the target number of current selection; 4 represent the number of parameters that each target will estimate; J (η) represents the target sequence number set selected.When carrying out first time iteration
now hypothesis does not have target;
represent a undetermined target sequence number, wherein
represent the position of undetermined target in distance quantifying unit,
represent and do not determine the position of target in angular spacing quantifying unit,
after iteration terminates, beyond J (η), the reflection coefficient of part is set as 0, is considered to do not have target; Wherein, the iterative formula of reflection coefficient is as follows:
In formula,
wherein in R
represent the reflection coefficient of the target that last iteration obtains, when first time iteration, in R
represent initial target reflection factor;
B, when target is moving-target, consider Doppler's effect
B1, determine first echoed signal that the n-th observation station observes and transmit, the reflection of target
Between coefficient and noise, relational model is as follows:
In formula, Y (n) is first echoed signal that the n-th observation station receives, wherein
p represents along the minimum range unit i.e. number of distance quantifying unit that directions of rays is divided into by distance in radar echo signal process, r=1 ... P; K represents that the angle intervals that there is target quantizes number of unit, a=1 ... ..K; G represents the interval quantizing number of unit of the Doppler that there is target, d=1 ... G; α
r, a, drepresent the reflection coefficient of target region internal object; b
an () is that the receiving array of the n-th observation station guides vector; a
an () is that the emission array of the n-th observation station guides vector;
represent a
an () carries out transpose operation,
will
conjugate transpose operation, wherein
For the waveform matrix that transmits of introducing Doppler obtained after zero padding process,
Wherein:
m=1 ... M
t, d=1 ... G;
S
mrepresent the waveform of m transmission antennas transmit orthogonal signal; ω
dit is the angle Doppler frequency of d Doppler's quantifying unit; M
tfor the number of emitting antenna; L is the number of subpulse of transmitting;
the Hadamard product of computing representing matrix; E (n) represents the additive noise in the n-th observation station Received signal strength, E
hn E (n) is carried out conjugate transpose operation by () expression;
B2, to above-mentioned relation model linearization;
B3, time delay is utilized to superpose the reflection coefficient of DAS algorithm initialization target, as follows:
Wherein:
Wherein
represent initialized target reflection factor,
represent v
r, a, dcarry out conjugate transpose operation; Vec represents matrix is carried out vectorization operation, and arranging successively by the order of rectangular array by matrix is row;
B4, Bayesian model exponent number selection tool is utilized to carry out the reflection coefficient of optimization aim, be specially: time delay is superposed the target reflection factor that DAS is initial, (2) formula of bringing into obtains the Bayes BIC value of initial value, the reflection coefficient of iterative target also upgrades Bayes BIC value, twice Bayes BIC value relatively, when Bayes BIC value no longer reduces, then iteration is stopped to obtain the reflection coefficient of target;
In formula, || ||
2represent Euclid norm,
represent Euclid norm square, η represents the target number of current selection; 5 represent the number of parameters that each target will estimate; J (η) represents the target sequence number set selected.When carrying out first time iteration
now hypothesis does not have target;
represent a undetermined target sequence number, wherein
represent the position of undetermined target in distance quantifying unit,
represent and do not determine the position of target in angular spacing quantifying unit,
represent and do not determine the position of target in Doppler's quantifying unit,
after iteration terminates, beyond J (η), the reflection coefficient of part is set as 0, is considered to do not have target;
Wherein, the iterative formula of reflection coefficient is as follows:
In formula,
r=1,...P,a=1,...K,d=1,...G
Wherein in R
represent the reflection coefficient of the target that last iteration obtains; When first time iteration, in R
represent initial target reflection factor.
The present invention compared with prior art has following significant progress:
(1) method of the present invention does not need auxiliary data, the v obtained after only needing the y by receiving and conversion goes to obtain the target component of target, thus solve existing space-time adaptive algorithm and need auxiliary data to go to estimate the master data of clutter and noise harmonics matrix, under non-homogeneous clutter environment, the problem of the auxiliary data of condition cannot be met.
(2) existing adaptive algorithm, such as Capon method and amplitude-phase algorithm for estimating (APES) can reduce the impact of clutter.But when the signal received is in low signal-to-noise ratio, if when the number of snapshot is insufficient, MIMO target detection performance will weaken greatly, but iteration self-adapting algorithm is customer parameter freedom and non-parametric adaptive algorithm, in several situation even only having a snapshot, still can obtain the parameter of target comparatively accurately.
(3) this method further improves the performance that DAS estimates the low and high secondary lobe of resolution.By the initial method that DAS method is estimated as target component, BIC model order selection tool is adopted to improve the accuracy of target component estimation.
(4) existing space-time adaptive algorithm operation quantity is high, but this algorithm has fast convergence rate, feature that stable performance operand is low.
Below in conjunction with the drawings and specific embodiments, the present invention is described in further detail;
Accompanying drawing explanation
Fig. 1 is the process flow diagram of relevant MIMO object detection system of the present invention;
Fig. 2 is DAS method moving target detect analogous diagram;
Fig. 3 is iteration self-adapting algorithm (IAA) moving target detect analogous diagram.
Embodiment
As shown in Figure 1, a kind of relevant MIMO radar object detection method based on iteration self-adapting algorithm of the present invention, the method obtains its location parameter by being optimized the reflection coefficient of moving-target and quiet target respectively, is specially:
A, when observed object is quiet target, does not consider Doppler's effect
A1, determine first echoed signal that the n-th observation station observes and transmit, between the reflection coefficient of target and noise, relational model is as follows:
In formula, Y (n) is first echoed signal that the n-th observation station receives, wherein
for the number of observation station; Y
hn Y (n) is carried out conjugate transpose operation by () expression; P represents along the minimum range unit i.e. number of distance quantifying unit that directions of rays is divided into by distance in radar echo signal process, r=1,2 ... P; K represents that the angle intervals that there is target quantizes number of unit, a=1 ... ..K; α
r,arepresent the reflection coefficient of target region internal object; b
an () is that the receiving array of the n-th observation station guides vector; a
an () is that the emission array of the n-th observation station guides vector;
represent a
an () does transpose operation; E (n) represents the additive noise of the n-th observation station Received signal strength, E
hn E (n) is carried out conjugate transpose operation by () expression;
representing will
carry out conjugate transpose operation,
For the waveform matrix that transmits obtained after zero padding process,
For the waveform matrix that transmits, wherein, s
mrepresent the orthogonal waveforms of m article of transmission antennas transmit, m=1 ... M
t; M
tfor the number of emitting antenna, L is the number of subpulse of transmitting; J
r∈ C
(L+P-1) × (L+P-1)for transformation matrix, be used for describing different distance and quantize the signal that receives of section, its expression formula is:
A2, to above-mentioned relation model linearization,
wherein
x=[α
1,1?α
1,2?…α
r,a…?α
P,K]
T∈C
PK×1 r=1,...P,a=1,...K;
M
rrepresent the number of receiving antenna; α
r,arepresent the reflection coefficient of target region internal object; Vec represents matrix is carried out vectorization operation, and arranging successively by the order of rectangular array by matrix is row;
A3, time delay is utilized to superpose the reflection coefficient of DAS algorithm initialization target, as follows:
Wherein
represent initialized target reflection factor,
represent V
r,acarry out conjugate transpose operation;
A4, Bayesian model exponent number selection tool is utilized to carry out the reflection coefficient of optimization aim, be specially: time delay is superposed the target reflection factor that DAS is initial, (1) formula of bringing into obtains the Bayes BIC value of initial value, the reflection coefficient of iterative target also upgrades Bayes BIC value, twice Bayes BIC value relatively, when Bayes BIC value no longer reduces, then iteration is stopped to obtain the reflection coefficient of target;
In formula, || .||
2represent Euclid norm,
represent Euclid norm square, η represents the target number of current selection; 4 represent the number of parameters that each target will estimate; J (η) represents the target sequence number set selected.When carrying out first time iteration
now hypothesis does not have target;
represent a undetermined target sequence number, wherein
represent the position of undetermined target in distance quantifying unit,
represent and do not determine the position of target in angular spacing quantifying unit,
after iteration terminates, beyond J (η), the reflection coefficient of part is set as 0, is considered to do not have target;
Wherein, the iterative formula of reflection coefficient is as follows:
In formula,
wherein in R
represent the reflection coefficient of the target that last iteration obtains, when first time iteration, in R
represent initial target reflection factor; Pass through
calculate new
upgrade Bayes BIC value, compare twice, front and back Bayes BIC value, be met condition when Bayes BIC value no longer reduces
B, when target is moving-target, consider Doppler's effect
B1, determine first echoed signal that the n-th observation station receives and transmit, between the reflection coefficient of target and noise, relational model is as follows:
In formula, Y (n) is first echoed signal that the n-th observation station receives, wherein
p represents along the minimum range unit i.e. number of distance quantifying unit that directions of rays is divided into by distance in radar echo signal process, r=1 ... P; K represents that the angle intervals that there is target quantizes number of unit, a=1 ... ..K; G represents the interval quantizing number of unit of the Doppler that there is target, d=1 ... G; α
r, a, drepresent the reflection coefficient of target region internal object; b
an () is that the receiving array of the n-th observation station guides vector; a
an () is that the emission array of the n-th observation station guides vector;
represent a
an () carries out transpose operation,
will
conjugate transpose operation, wherein
For the waveform matrix that transmits of introducing Doppler obtained after zero padding process,
Wherein:
m=1 ... M
t, d=1 ... G;
S
mrepresent the waveform of m transmission antennas transmit orthogonal signal; ω
dit is the angle Doppler frequency of d Doppler's quantifying unit; M
tfor the number of emitting antenna; L is the number of subpulse of transmitting;
the Hadamard product of computing representing matrix; E (n) represents the additive noise in the n-th observation station Received signal strength, E
hn E (n) is carried out conjugate transpose operation by () expression;
B2, to above-mentioned relation model linearization, as follows:
y=Ax+e
Wherein, A=[v
1,1,1v
1,1,2v
r, a, dv
p, K, G] x=[α
1,1,1α
1,1,2α
r, a, dα
p, K, G]
t
M
rrepresent the number of receiving antenna; α
r, a, drepresent the reflection coefficient of target region internal object; Vec represents matrix is carried out vectorization operation, and arranging successively by the order of rectangular array by matrix is row;
B3, time delay is utilized to superpose the reflection coefficient of DAS algorithm initialization target, as follows:
Wherein
represent initialized target reflection factor,
represent v
r, a, dcarry out conjugate transpose operation;
B4, Bayesian model exponent number selection tool is utilized to carry out the reflection coefficient of optimization aim, be specially: time delay is superposed the target reflection factor that DAS is initial, (2) formula of bringing into obtains the Bayes BIC value of initial value, the reflection coefficient of iterative target also upgrades Bayes BIC value, twice Bayes BIC value relatively, when Bayes BIC value no longer reduces, then iteration is stopped to obtain the reflection coefficient of target;
In formula, || .||
2represent Euclid norm,
represent Euclid norm square, η represents the target number of current selection; 5 represent the number of parameters that each target will estimate; J (η) represents the target sequence number set selected.When carrying out first time iteration
now hypothesis does not have target;
represent a undetermined target sequence number, wherein
represent the position of undetermined target in distance quantifying unit,
represent and do not determine the position of target in angular spacing quantifying unit,
represent and do not determine the position of target in Doppler's quantifying unit,
after iteration terminates, beyond J (η), the reflection coefficient of part is set as 0, is considered to do not have target;
Wherein, the iterative formula of reflection coefficient is as follows:
In formula,
R=1 ... P, a=1 ... K, d=1 ... G is wherein in R
represent the reflection coefficient of the target that last iteration obtains; When first time iteration, in R
for initial target reflection factor.
Effect of the present invention can be illustrated by following emulation:
1, simulated conditions and content
And the moving target detect of the to the effect that relevant MIMO radar of our research comparatively complicated due to the situation of can not ignore Doppler, therefore the present invention carries out the detection simulation of the parameter when target is moving-target, the target of random generation 4 motion, suppose that all targets are all positioned at the position that Doppler's number is 1, the region that we are concerned about is angle is from-30 degree to 30 degree, and doppler angle is from-30 degree to 30 degree.Emitting antenna Mt=5, accept antenna Mr=5, the number of subpulse is L=32, the frequency of operation of carrier wave is 1.24GHz, the wavelength of the distance dt=2.5* carrier wave between emitting antenna, accept the wavelength of the distance dr=0.5* carrier wave between antenna, to-noise ratio SNA=20, wherein SNA=10log (tr (S
hs)/L σ
2), the mark of tr representing matrix, noise
for M
rwhite Gaussian noise, σ that row L+P-1 arranges
2for the variance of Gaussian noise; Distance is divided into 12 frequency bands and P=12, angular divisions is 7 frequency bands and K=7, doppler angle are divided into 5 frequency bands and H=5.Observation station
namely an observation station is only had.Iteration self-adapting algorithm is adopted to emulate target localization.
2, the result of emulation as shown in Figures 2 and 3, when adopting DAS method to carry out parameter estimation to target, there is lower resolution and more weak AF panel, be difficult to estimate target accurately, a lot of mistakes can be produced as the place that will be target, being estimated as target, when adopting iteration self-adapting algorithm (IAA) to carry out moving target detect to it, the parameter of target can being oriented under similar circumstances comparatively accurately.Moreover due to this method be customer parameter freely and non-parametric adaptive algorithm, so still target can be oriented very accurately when the situation that snapshot number lacks even only has a snapshot.In addition, the method does not need auxiliary data.The v obtained after only needing the y by receiving and conversion goes to obtain the target component of target, thus solve existing space-time adaptive algorithm and need auxiliary data to go to estimate the master data of clutter and noise harmonics matrix, but under non-homogeneous clutter environment, the problem of the auxiliary data of condition cannot be met.And this algorithm is compared with DAS method, when improving resolution and reducing secondary lobe, time and the DAS of this method computing are almost equal, compared to the space-time adaptive algorithm needing auxiliary data, because this method does not need aided algorithm therefore to improve arithmetic speed to a certain extent.And the method has the feature of fast convergence rate, stable performance.
To sum up, the problem that the present invention adopts iteration self-adapting algorithm effectively can overcome existing moving-target parameter detecting to run into, greatly improves the accuracy of relevant MIMO radar moving target detect.
Claims (1)
1., based on a relevant MIMO radar object detection method for iteration self-adapting algorithm, obtain its location parameter by being optimized the reflection coefficient of moving-target and quiet target respectively, the method is specially:
A, when observed object is quiet target, do not consider Doppler's effect
A1, determine first echoed signal that the n-th observation station observes and transmit, the reflection of target
Between coefficient and noise, relational model is as follows:
In formula, Y (n) is first echoed signal that the n-th observation station observes, wherein
for the number of observation station; Y
hn Y (n) is carried out conjugate transpose operation by () expression; P represents along the minimum range unit i.e. number of distance quantifying unit that directions of rays is divided into by distance in radar echo signal process, r=1,2 ... P; K represents that the angle intervals that there is target quantizes number of unit, a=1 ... ..K; α
r,arepresent the reflection coefficient of target region internal object; b
an () is that the receiving array of the n-th observation station guides vector; a
an () is that the emission array of the n-th observation station guides vector;
represent a
an () does transpose operation; E (n) represents the additive noise of the n-th observation station Received signal strength, E
hn E (n) is carried out conjugate transpose operation by () expression;
representing will
carry out conjugate transpose operation,
For the waveform matrix that transmits obtained after zero padding process,
For the waveform matrix that transmits, wherein, s
mrepresent the orthogonal waveforms of m article of transmission antennas transmit, m=1 ... M
t; M
tfor the number of emitting antenna, L is the number of subpulse of transmitting;
J
r∈ C
(L+P-1) × (L+P-1)for transformation matrix, be used for describing different distance and quantize the signal that receives of section, its expression formula is:
A2, to above-mentioned relation model linearization;
A3, time delay is utilized to superpose the reflection coefficient of DAS algorithm initialization target, as follows:
Wherein:
r=1,...P,a=1,...K;
Wherein
represent initialized target reflection factor,
represent V
r,acarry out conjugate transpose operation; M
rrepresent the number of receiving antenna; Vec represents matrix is carried out vectorization operation, and arranging successively by the order of rectangular array by matrix is row;
A4, Bayesian model exponent number selection tool is utilized to carry out the reflection coefficient of optimization aim, be specially: time delay is superposed the target reflection factor that DAS is initial, (1) formula of bringing into obtains the Bayes BIC value of initial value, the reflection coefficient of iterative target also upgrades Bayes BIC value, twice Bayes BIC value relatively, when Bayes BIC value no longer reduces, then iteration is stopped to obtain the reflection coefficient of target;
In formula, || .||
2represent Euclid norm,
represent Euclid norm square, η represents the target number of current selection; 4 represent the number of parameters that each target will estimate; J (η) represents the target sequence number set selected when this iteration, when carrying out first time iteration
now hypothesis does not have target;
represent a undetermined target sequence number, wherein
dong represents the position of undetermined target in distance quantifying unit,
represent and do not determine the position of target in angular spacing quantifying unit,
after iteration terminates, beyond J (η), the reflection coefficient of part is set as 0, is considered to do not have target;
Wherein, the iterative formula of reflection coefficient is as follows:
In formula,
wherein in R
represent the reflection coefficient of the target that last iteration obtains, when first time iteration, in R
represent initial target reflection factor;
B, when target is moving-target, consider Doppler's effect
B1, determine first echoed signal that the n-th observation station observes and transmit, between the reflection coefficient of target and noise, relational model is as follows:
In formula, Y (n) is first echoed signal that the n-th observation station observes, wherein
p represents along the minimum range unit i.e. number of distance quantifying unit that directions of rays is divided into by distance in radar echo signal process, r=1 ... P; K represents that the angle intervals that there is target quantizes number of unit, a=1 ... ..K; G represents the interval quantizing number of unit of the Doppler that there is target, d=1 ... G; α
r, a, drepresent the reflection coefficient of target region internal object; b
an () is that the receiving array of the n-th observation station guides vector; a
an () is that the emission array of the n-th observation station guides vector;
represent a
an () carries out transpose operation,
will
conjugate transpose operation, wherein
For the waveform matrix that transmits of introducing Doppler obtained after zero padding process,
Wherein:
m=1 ... M
t, d=1 ... G;
S
mrepresent the waveform of m transmission antennas transmit orthogonal signal; ω
dit is the angle Doppler frequency of d Doppler's quantifying unit; M
tfor the number of emitting antenna; L is the number of subpulse of transmitting;
the Hadamard product of computing representing matrix; E (n) represents the additive noise in the n-th observation station Received signal strength, E
hn E (n) is carried out conjugate transpose operation by () expression;
B2, to above-mentioned relation model linearization;
B3, time delay is utilized to superpose the reflection coefficient of DAS algorithm initialization target, as follows:
Wherein:
Wherein
represent initialized target reflection factor,
represent v
r, a, dcarry out conjugate transpose operation; Vec represents matrix is carried out vectorization operation, and arranging successively by the order of rectangular array by matrix is row;
B4, Bayesian model exponent number selection tool is utilized to carry out the reflection coefficient of optimization aim, be specially: time delay is superposed the target reflection factor that DAS is initial, (2) formula of bringing into obtains the Bayes BIC value of initial value, the reflection coefficient of iterative target also upgrades Bayes BIC value, twice Bayes BIC value relatively, when Bayes BIC value no longer reduces, then iteration is stopped to obtain the reflection coefficient of target;
In formula, || .||
2represent Euclid norm,
represent Euclid norm square, η represents the target number of current selection; 5 represent the number of parameters that each target will estimate; J (η) represents the target sequence number set selected; When carrying out first time iteration
now hypothesis does not have target;
represent a undetermined target sequence number, wherein
represent the position of undetermined target in distance quantifying unit,
represent and do not determine the position of target in angular spacing quantifying unit,
represent and do not determine the position of target in Doppler's quantifying unit,
after iteration terminates, beyond J (η), the reflection coefficient of part is set as 0, is considered to do not have target;
Wherein, the iterative formula of reflection coefficient is as follows:
In formula,
r=1,...P,a=1,...K,d=1,...G
Wherein in R
represent the reflection coefficient of the target that last iteration obtains; When first time iteration, in R
represent initial target reflection factor.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310162479.4A CN103257344B (en) | 2013-05-03 | 2013-05-03 | Iteration-adaptive-algorithm-based method for detecting coherent MIMO radar target |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310162479.4A CN103257344B (en) | 2013-05-03 | 2013-05-03 | Iteration-adaptive-algorithm-based method for detecting coherent MIMO radar target |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103257344A CN103257344A (en) | 2013-08-21 |
CN103257344B true CN103257344B (en) | 2015-01-14 |
Family
ID=48961368
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310162479.4A Expired - Fee Related CN103257344B (en) | 2013-05-03 | 2013-05-03 | Iteration-adaptive-algorithm-based method for detecting coherent MIMO radar target |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103257344B (en) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103744076B (en) * | 2013-12-25 | 2016-07-06 | 河海大学 | MIMO radar moving target detection method based on non-convex optimization |
CN103810325B (en) * | 2014-01-08 | 2016-08-17 | 西安电子科技大学 | The linear thinned array antenna optimization method of low sidelobe based on SQP |
CN103777190B (en) * | 2014-02-26 | 2016-03-30 | 南京信息工程大学 | A kind of angle estimating method of bistatic MIMO radar high speed highly maneuvering target |
CN104515975B (en) * | 2014-12-12 | 2017-02-22 | 中国电子科技集团公司电子科学研究院 | Coherent MIMO (multiple input multiple output) radar waveform design method facing clutter suppression |
CN104635229B (en) * | 2014-12-12 | 2017-02-22 | 中国电子科技集团公司电子科学研究院 | Coherent MIMO (multiple-input multiple-output) radar waveform design method oriented to unambiguous zone area expansion |
CN104868946B (en) * | 2015-06-12 | 2017-11-17 | 哈尔滨工业大学 | The disturbance restraining method of adaptive weighted Subarray mixing MIMO phased array systems |
TWI588507B (en) * | 2015-10-14 | 2017-06-21 | 國立成功大學 | Radio frequency energy-transmitting apparatus with location detection function and radio frequency energy-harvesting apparatus and radio frequency energy-transmitting method with location detection function |
TWI633710B (en) * | 2017-01-17 | 2018-08-21 | 台達電子工業股份有限公司 | A polarization-tracking rf power transmitting apparatus, a positioning and polarization-tracking rf energy-transmitting apparatus, an rf energy hunting apparatus and method of rf power transfer the same |
US11592548B2 (en) * | 2019-09-27 | 2023-02-28 | Intel Corporation | Methods and apparatus to improve doppler velocity estimation |
CN111487609A (en) * | 2020-04-20 | 2020-08-04 | 中国人民解放军海军航空大学 | Multi-frequency continuous wave MIMO array radar system and target parameter estimation method thereof |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3113148B2 (en) * | 1994-06-01 | 2000-11-27 | 三菱電機株式会社 | Radar equipment |
CN101915911B (en) * | 2010-08-26 | 2013-05-01 | 哈尔滨工业大学 | Space arbitrary configuration distributed SAR moving target parameter estimation method based on cancellation accumulation space-time spectrum |
-
2013
- 2013-05-03 CN CN201310162479.4A patent/CN103257344B/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN103257344A (en) | 2013-08-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103257344B (en) | Iteration-adaptive-algorithm-based method for detecting coherent MIMO radar target | |
CN102279387B (en) | Method for estimating target arrival angle of multiple input multiple output (MIMO) radar | |
CN105699945B (en) | Waveform optimization design method in frequency control battle array MIMO radar system | |
CN102156279B (en) | Method for detecting moving target on ground by utilizing bistatic radar based on MIMO (Multiple Input Multiple Output) | |
CN103901417B (en) | Low-complexity space target two-dimensional angle estimation method of L-shaped array MIMO radar | |
CN102520395B (en) | Clutter suppression method based on bistatic multiple-input and multiple-output radar | |
CN105738879B (en) | Radar clutter space-time adaptive pre-filtering method based on sparse recovery | |
CN103353591B (en) | Bistatic radar localization dimension reduction clutter suppression method based on MIMO | |
CN103901410B (en) | Airborne bistatic MIMO radar clutter suppression method based on sparse recovery | |
CN103744076B (en) | MIMO radar moving target detection method based on non-convex optimization | |
CN106353744A (en) | Multi-parameter combined estimation method based on bi-static FDA-MIMO radars | |
CN107561486B (en) | Shallow sea target direction-of-arrival estimation method based on active time reversal | |
CN103353592B (en) | Bistatic radar multichannel combination dimension reduction clutter suppression method based on MIMO | |
CN102135617A (en) | Multi-target positioning method of bistatic multi-input multi-output radar | |
CN102608595B (en) | Target location method based on distributed coherent process for meterwave MIMO (multiple-input multiple-output) radars | |
CN103728606A (en) | Doppler channel correlation two-stage dimension reduction method for onboard multiple input multiple output (MIMO) radar | |
CN102866388B (en) | Iterative computation method for self-adaptive weight number in space time adaptive processing (STAP) | |
CN106291540A (en) | A kind of multiple-input and multiple-output GPR backwards projection target imaging method estimated based on DOA | |
CN103760546B (en) | A kind of radar low target Wave arrival direction estimating method | |
CN103983958A (en) | MIMO radar continuous target angle estimation method based on multi-measurement-vector sparse representation | |
CN104155633B (en) | Clutter suppression method of non-positive side-looking bistatic MIMO radar | |
CN105182313A (en) | MIMO-STAP steady waveform design method based on incomplete clutter prior knowledge | |
CN103969640A (en) | Bistatic MIMO radar target sparsity imaging method | |
CN105807275A (en) | MIMO-OFDM-STAP steady waveform design method based on partial clutter priori knowledge | |
CN105182325A (en) | Rank-1 constraint based metrewave MIMO radar low-elevation-angle object height measuring method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
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: 20150114 Termination date: 20170503 |