CN110579737A - Sparse array-based MIMO radar broadband DOA calculation method in clutter environment - Google Patents
Sparse array-based MIMO radar broadband DOA calculation method in clutter environment Download PDFInfo
- Publication number
- CN110579737A CN110579737A CN201910646812.6A CN201910646812A CN110579737A CN 110579737 A CN110579737 A CN 110579737A CN 201910646812 A CN201910646812 A CN 201910646812A CN 110579737 A CN110579737 A CN 110579737A
- Authority
- CN
- China
- Prior art keywords
- doa
- sparse
- vector
- array
- sparse array
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 6
- 238000000034 method Methods 0.000 claims abstract description 31
- 238000005457 optimization Methods 0.000 claims abstract description 11
- 239000013598 vector Substances 0.000 claims description 43
- 239000011159 matrix material Substances 0.000 claims description 17
- 238000005070 sampling Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000012804 iterative process Methods 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 abstract description 12
- 238000012545 processing Methods 0.000 abstract description 5
- 238000011156 evaluation Methods 0.000 abstract description 2
- 238000003491 array Methods 0.000 description 10
- 230000009977 dual effect Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 2
- 230000001174 ascending effect Effects 0.000 description 2
- 239000003990 capacitor Substances 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000017105 transposition Effects 0.000 description 2
- 208000001992 Autosomal Dominant Optic Atrophy Diseases 0.000 description 1
- 206010011906 Death Diseases 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000000593 degrading effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
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
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
- G01S3/143—Systems for determining direction or deviation from predetermined direction by vectorial combination of signals derived from differently oriented antennae
-
- 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/023—Interference mitigation, e.g. reducing or avoiding non-intentional interference with other HF-transmitters, base station transmitters for mobile communication or other radar systems, e.g. using electro-magnetic interference [EMI] reduction techniques
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S2013/0236—Special technical features
- G01S2013/0245—Radar with phased array antenna
- G01S2013/0254—Active array antenna
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a sparse array-based MIMO radar broadband DOA calculation method in a clutter environment, and belongs to the field of signal processing. Specifically WCSAB, in the method, clutter interference is suppressed by Capon beam forming, and then target DOA is estimated by a CS method by jointly utilizing different narrowband signal information. Considering that the DOA estimation performance is not only related to the beam forming weight value but also related to the sparse array structure, the invention provides a joint optimization problem of the beam forming weight value and the sparse array, and provides a simple algorithm for solving the optimization problem. The method provided by the invention can improve the performance of target DOA estimation in a clutter environment, including high resolution and low sidelobe, the sparse array reduces the system cost and complexity, the Bayesian Mean Square Error (BMSE) of the target DOA estimation is taken as a performance evaluation index, and the sparse array structure designed by the algorithm is similar to the optimal sparse array performance obtained by an exhaustive method and has better performance than a nested array and a co-prime array.
Description
Technical Field
The invention belongs to the field of signal processing, and particularly relates to a sparse array-based MIMO radar broadband DOA estimation problem in a clutter environment.
Background
An mimo (Multiple Input Multiple output) radar is a novel radar system that synchronously transmits signals by using Multiple transmitting antennas, receives echo signals by using Multiple receiving antennas, and processes the signals in a centralized manner. Compared with the traditional phased array radar, the MIMO radar has obvious advantages such as higher resolution, better target detection, positioning and tracking performance and better target parameter estimation and identification capability. DOA estimation research is an important content in array signal processing, and the application of the DOA estimation research relates to the fields of radar, communication, sonar, radio astronomy, survey, earthquake, biomedicine and the like. At present, there are many classical DOA Estimation methods, such as Multiple Signal Classification (MUSIC), Signal parameter Estimation based on rotation invariant technology (ESPRIT), and so on. In recent years, the theory of Compressive Sensing (CS) has gained wide attention of scholars at home and abroad, and compared with the traditional method, the CS-based MIMO radar DOA estimation has better estimation performance under the conditions of less sampling data and low signal-to-noise ratio.
According to the traditional array signal processing theory, in order to ensure the uniqueness of DOA estimation, the distance between adjacent array elements in the array is less than or equal to half wavelength of an incident signal, and the array meeting the condition is called a full array. The spatial resolution of the array is related to the array aperture, and increasing the resolution requires increasing the array aperture, which means that more antennas are required in a full array. However, due to the practical constraints of software and hardware resources, the number of antennas is usually limited. Sparse arrays have attracted considerable attention in order to increase the array aperture without increasing the number of antennas. When the target is sparse in the observation space, the target DOA can be accurately estimated by the sparse array. However, in a clutter environment, the sparsity of the target in the observation space may be destroyed, thereby causing degradation of DOA estimation performance.
Considering that broadband signals have the advantages of large information amount, strong anti-interference capability, high resolution and the like, relatively representative broadband DOA estimation methods include an Incoherent Signal Subspace (ISSM), a Coherent Signal Subspace (CSSM) and the like. The ISSM divides the broadband signal into a plurality of narrow-band signals on a frequency band, then processes each narrow-band signal respectively, and finally averages the processing results of all the narrow-band signals to obtain a final estimation result. CSSM transforms the covariance matrix of narrowband signals of different frequencies to a reference frequency by focusing and then obtains the final result by using a narrowband estimation method. However, CSSM requires an estimate of the target DOA, and performance is greatly affected by the accuracy of the estimate.
Disclosure of Invention
The invention provides a sparse array-based MIMO radar broadband DOA estimation method in a clutter environment, in particular to WCSAB (wideband compressed sensing after beamforming). Considering that the DOA estimation performance is not only related to the beam forming weight value but also related to the sparse array structure, the invention provides a joint optimization problem of the beam forming weight value and the sparse array, and provides a simple algorithm for solving the optimization problem.
The technical scheme of the invention is a sparse array-based MIMO radar broadband DOA calculation method in a clutter environment, which comprises the following steps:
Step 1: let the position of the transmitting antenna be determined, and the feasible region for placing the receiving antenna be [0, Dr]For simplicity of analysis, the feasible fields are separated by an interval ΔrDiscretization to NrA plurality of grid points, and N receiving antennas disposed on some of the grid points, N < Nr;
Step 2: establishing an MIMO radar echo signal model to obtain echo signal time domain sampling data n=1,...,NrAnd p 1., L, where p represents a time domain snapshot and L is a snapshot number;
And step 3: for received signalPerforming an L-point discrete Fourier transform to obtain frequency domain data, i.e.
And N isrThe data of each lattice point is expressed in vector form, i.e. y [ l ]]=[y1[l],...,yNr[l]]TWherein p 1, L1;
And 4, step 4: discretizing the target angle observation area into G grid points theta1,...,θGK < G, where K represents the number of targets, the signal model is represented in sparse form:
y[l]=Φ[l]x+c[l]+u[l]
WhereinWhere a isr(θ,fl) Indicating a received steering vector, at(θ,fl) Representing the transmit steering vector, s [ l ]]Representing a frequency domain transmitted signal, x ═ x1,...,xG]TIs K sparse, i.e. x has only K non-zero elements, and the values and positions of the non-zero elements are the target reflection coefficient and DOA, c [ l [ ]]Represents clutter u [ l ]]Representing noise;
And 5: will beam form weight vector wg,lact on y [ l]The beam forming output result is obtained:
Will r isg,l(G1., G, L ═ 1., L) is expressed as a vector:
r=[r1,1,...,rG,1,...,r1,L,...,rG,L]T
=WrΦx+Wrc+Wru
wherein the weight matrix Wr=Diag{W1,...,Wl,...,WLIs a block diagonal matrix;
and is provided withΦ=[ΦT[1],...,ΦT[L]]T,c=[cT[1],...,cT[L]]TRepresents clutter, u ═ uT[1],...,uT[L]]TRepresenting noise;
Step 6: reconstructing a sparse vector x by basis pursuit desiccation based on the CS theory;
wherein eta is equal to or greater than 0 is a regularization parameter;
And 7: to the solution obtained in step 6The values of the elements in the table are sorted from big to small, and the corresponding lattice point of each sorted element is expressed as { theta [ [ theta ] ](1),...,θ(G)Then the DOA estimation result can be expressed as;
and 8: based on minimum Bayes mean square errorsolving for optimal WrThe following optimization problem is established
s.t.Wr=Diag{W1,...,WL}
||W||0=N
W=[w1,1,...,w1,L,...,wG,1,...,wG,L]
Wherein the DOA vector theta of the real targetTIs random in nature and is not only easy to be recognized,Is expressed in the pair thetaTin the hope of expectation,Denotes thetaTMean square error of DOA estimation, w, when determinedg,lRepresents a weight vector, where G1.., G, L1., L;
And step 9: the problem provided by the step 8 is solved in an optimized way to obtain the optimal Wr。
Further, the specific method of step 9 is as follows:
Step 1: initialization: the number of iterations j is 1,According to the formulaComputing beamforming weightsWhereinRc(fl) Is a clutter c [ l ]]The covariance matrix of (a); during each iteration, a set of lattice point selection vectors z is randomly generated1,...,zαfor a given z;
Step 2: repeating the iterative process from the step 3 to the step 6:
And step 3: randomly generating a set of lattice point selection vectors z1,...,zα};
And 4, step 4: according to the formula wg,l=z⊙ξg,lComputingAnd form
According to the formulaCalculating rg,lR is tog,lExpressed as a vector r;
Will be provided withSubstituting r into the formulaObtaining x reconstruction resultsAnd target DOA estimation result
According to the formulaobtaining BMSE
And 5: obtaining based on minimum BMSE
Step 6: according toBased onUpdatingGet the corresponding weight valueAnd let j equal j + 1; whereinIs Dc [ l ]]The covariance matrix of (a) is determined,
And 7: when in useThen, iteration stops and the optimal antenna selection is outpute0Is a threshold value set in advance.
the method provided by the invention can improve the performance of target DOA estimation in a clutter environment, including high resolution and low sidelobe, the sparse array reduces the system cost and complexity, the Bayesian Mean Square Error (BMSE) of the target DOA estimation is taken as a performance evaluation index, and the sparse array structure designed by the algorithm is similar to the optimal sparse array performance obtained by an exhaustive method and has better performance than a nested array and a co-prime array.
Drawings
FIG. 1 shows the results of BMSE ordering in ascending order for all possible sparse array configurations, and for comparison, FIG. 1 also shows the results for nested arrays (nested arrays) and co-prime arrays (co-prime arrays).
fig. 2(a) shows the optimal sparse array structure based on the minimum BMSE condition, and fig. 2(b) shows the sparse array structure obtained by the proposed algorithm according to the present invention.
Fig. 3 shows DOA estimation results for different sparse array structures when using the WCSAB method.
FIG. 4 shows the DOA estimation results of different array structures using WCSAB and WCT (wideband capacitor technique) methods for single target cases.
FIG. 5 shows the DOA estimation results of different array structures when WCSAB and WCT (wideband capacitor technique) methods are used respectively, considering the dual target case.
Detailed Description
For convenience of description, the following definitions are first made:
Bold capital letters represent matrices, bold lowercase letters represent vectors, (.)*For conjugation, (.)TIs a transposition ofHFor the conjugate transposition, | x | | non-conducting phosphor0And | | x | | non-conducting phosphor1Respectively representing l of the vector x0Norm sum l1Norm, | W | count0Representing the number of non-zero rows of the matrix W, Diag {. cndot.) representing a block diagonal matrix, Diagr{. denotes the diagonal matrix after removing zero rows,Denotes the expectation with respect to theta, INFor a unit array of order N, 1 is a full 1 vector, and a symbol &indicatesa Hadamard product.
Consider a co-located MIMO radar system with both transmit and receive antennas placed on the horizontal axis of a two-dimensional cartesian coordinate system. Suppose there are M transmit antennas and the position on the horizontal axis is known as dt,m(M ═ 1.., M). Assume that the feasible domain for placing the receive antenna is [0, Dr]for simplicity of analysis, the feasible fields are separated by an interval ΔrDiscretization to NrAnd grid points on which the receiving antennas are placed. Due to the constraint of the number of antennas, the radar system is assumed to have only N (N < N)r) And available receiving antennas. Order torepresenting a wideband signal transmitted by the mth transmitting antenna, having a frequency range of [ -B ]m/2,Bm/2]where p denotes a time-domain snapshot, TsIndicating the sampling period and L the number of fast beats. Let DOA of K far-field point targets be thetaT,k(K1.. K), then the signal received at the nth bin is the signal received at the nth bin
wherein f iscdenotes the carrier frequency, betakThe reflection coefficient of the kth target is represented and assumed to be unknown. Let the first transmitting antenna and the first lattice point doFor reference, then τTt,m,k=(dt,m-dt,1)sinθT,kC represents the time delay of the signal from the m-th transmitting antenna to the k-th target relative to the reference array element, tauTr,n,k=(n-1)ΔrsinθT,kRepresenting the time delay of the signal from the kth target to the nth grid point relative to the first grid point. Q denotes the number of clutter scatterers, gammaq(Q ═ 1., Q) denotes the reflection coefficients of clutter scatterers, and it is assumed that there are independent identically distributed (iid) gaussian random variables between them. Tau isCt,m,q=(dt,m-dt,1)sinθC,qC represents the time delay of the signal from the m-th transmitting antenna to the q-th clutter-scattering body relative to the reference array element, tauCr,n,q=(n-1)ΔrsinθC,qRepresenting the time delay of the signal from the qth clutter scatterer to the nth grid point, θC,qIndicating the orientation of the q-th clutter scatterer relative to the array.Is a variance of σ2White gaussian noise.
by performing an L-point Discrete Fourier Transform (DFT) on the time domain discrete signal, the frequency point f can be obtainedl=lfsfrequency domain data of (1...., L), where fsIs the frequency sampling interval, fl∈[-B/2,B/2]And isThe signal being at frequency flThe DFT result of (a) is
wherein s ism[l]And un[l]Respectively representing the transmitted signalsAnd noiseDFT of (2). Order toAndRespectively expressed at an angle theta and a frequency flA receive steering vector and a transmit steering vector. Will NrThe signal received by each grid point is expressed as a vector
Wherein
Under the CS framework, to estimate DOA θ of K targetsT,k(K1.. K.) the target angle observation region is discretized into G (K < G) grid points θ1,...,θGIt is assumed that the dispersion error is negligible, i.e. the target falls exactly on the grid point. Then the formula (3) can be expressed as
y[l]=Φ[l]x+c[l]+u[l] (4)
WhereinVector x ═ x1,...,xG]Tis K sparse, i.e. x has only K non-zero elements, and the values and positions of the non-zero elements are the target reflection coefficient and DOA, which can be expressed as
CS theory estimates the target DOA using the sparsity of x, however, this sparsity is destroyed in a cluttered environment, thereby degrading the performance of DOA estimation. Is composed ofThe interference of the spurious wave is suppressed, and a beam forming method is adopted at a receiving end. Order toIs shown in the direction thetagFrequency flAnd the position of the non-zero element represents the lattice point where the antenna is selected to be placed. Since there are only N available receive antennas, the weight vector is required to satisfy | | wg,l||0N. The beamformed output is given by
Will r isg,l(G1., G and L1., L) is expressed as a vector GL × 1
Where phi is [ ]T[1],...,ΦT[L]]T,c=[cT[1],...,cT[L]]T,u=[uT[1],...,uT[L]]T,Wr=Diag{W1,...,WL},according to equation (7), the DOA estimation problem can be converted into a sparse signal reconstruction problem, and based on the CS theory, the K sparse vector x can be reconstructed through base pursuit elimination (BPDN)
where η ≧ 0 is the regularization parameter, for the optimization problem of equation (8), the CVX toolkit can be used to solve. Order toThe solution of the above formula is shown,is the estimation result of the target DOA. Considering a target DOA vector θT=[θT,1,...,θT,K]Tis a random case, then the average estimation performance can be given by the Bayesian Mean Square Error (BMSE)
From the equation (9), the DOA estimated performance and the matrix WrIn connection with, in order to optimize performance, the following optimization problem is given
(10) The last two constraints in the equation are to ensure that w is true for different g and lg,lThe positions of the non-zero elements in (a) are the same. Due to wg,lthe position of the non-zero element in the (10) expression indicates that the corresponding lattice point is selected to place the antenna, so that the formula is a joint optimization problem of the weighted value and the sparse array structure.
Considering that equation (10) is an NP-hard problem, a simple algorithm is proposed to solve the optimization problem. The algorithm firstly optimizes the sparse array structure when the weight is given, and then updates the weight for the next iteration. First, how to optimize the sparse array structure when the weight is given is explained. First, a lattice point selection vector is definedWherein z isnBelongs to {0,1}, only if the element is 1, the corresponding lattice point is selected to place the antenna, and because only N available receiving antennas exist, the requirement of | z | calculation of the virtues0N. The weight of the first iteration is given by the formation of a full array condition Capon beam
WhereinRc(fl) Is a clutter c [ l ]]the covariance matrix of (2). During each iteration, a set of lattice point selection vectors z is randomly generated1,...,zαFor a given z, there is
wg,l=z⊙ξg,l (12)
For different z, different w can be obtainedg,lAnd Wr(WrFrom wg,lComposition) of formula (9), BMSE and WrIn this regard, it can thus be seen that BMSE is also related to z, denoted as e (z). Based on minimum BMSE, the optimal lattice point selection vector z can be obtainedop。
Based on zopUpdating xig,lin the corresponding weight value, orderXi isg,lIn the formula zopThe selected element value is updated by
WhereinIs Dc [ l ]]The covariance matrix of (a) is determined,When BMSE e (z)op) Less than a certain threshold e0When so, the iteration stops. The detailed algorithm is given in table 1.
TABLE 1 iterative algorithm for solving optimization problem
In order to suppress the interference of the spurious waves, a beam forming method is adopted at the receiving end.is shown in the direction thetagFrequency flThe position of the non-zero element represents the lattice point where the antenna is selected to be placed, and since there are only N available receiving antennas, the weight vector is required to satisfy | | wg,l||0n. The array structure of the antenna should be the same for different g and l, i.e. wg,lThe positions of the non-zero elements in (a) are the same. To represent this constraint, a matrix is constructed
W=[w1,1,...,w1,L,...,wG,1,...,wG,L] (14)
And satisfy | | W | count of the north0By this constraint, w can be satisfied for different g and l, i.e. the number of non-zero rows of the matrix is Ng,lthe positions of the non-zero elements in (a) are the same. Weight vector wg,lacting on received signal y [ l]Obtaining a beam forming output r according to the formula (6)g,lIt is expressed as a G × 1 vector
From the above formula, it can be found that for different frequencies flThe vectors x are equally sparse. In order to jointly utilize signal information of different frequencies, r isl(L ═ 1., L) is expressed as a vector of GL × 1, i.e., expression (7). By converting the DOA estimation problem into the sparse signal reconstruction problem, the reconstruction result of the sparse vector x can be obtained according to the formula (8)ThenThe position of the maximum K elements in the target DOA is the estimation result and is expressed asTo pairthe values of the elements in the table are sorted from big to small, and the corresponding lattice point of each sorted element is expressed as { theta [ [ theta ] ](1),...,θ(G)Then the DOA estimation result can be expressed as
Full array beam forming weight vector xig,lCan make the direction thetagFrequency flThe signal passes through without distortion while suppressing interference and noise in other directions, and is expressed as
The optimal solution of the above equation is equation (11).
Two simulation examples are given for sparse array-based MIMO radar broadband DOA estimation in a clutter environment, and the parameters are set as follows: suppose Dr11 λ/2, where λ represents the wavelength corresponding to the highest frequency of the signal. The feasible region in which the receiving antenna can be placed is divided by deltarλ/2 is the interval dispersion of 12 grid points. The number of transmitting and receiving antennas available for the MIMO radar system is assumed to be M-N-6, and the transmitting end array structure is determined to be known.
to simplify the analysis, it is assumed that the transmitted signal bandwidths are the same, i.e. Bm200MHz (M1.., M), the carrier frequency is 1 GHz.
the target angular field of view is discretized into 41 grid points-20 °, -19 °,20 °.
The clutter consists of 250 scatterers, distributed at angles of-90 °, -90 ° +180 °/250 °,90 °.
Defining a signal-to-noise ratioAnd signal to noise ratioWithout loss of generality, assuming a target reflection coefficient of 1, SNR and SCR are set to-5 dB and-30 dB, respectively.
in simulation 1, it is assumed that the targets are uniformly and randomly distributed on the grid points after the angular observation domain is discretized. To ensure that the aperture of the array does not change, letThen shareDifferent sparse array structures. Fig. 1 shows the results of the BMSE arrangement in ascending order for all possible sparse array structures, diamonds represent the optimal sparse array structure under the minimum BMSE condition, squares represent the sparse array structures obtained according to the algorithm given in table 1, and it can be seen that the two structures have similar performance, and the specific structures of the two sparse arrays are given in fig. 2. For comparison, fig. 1 also shows the results of a nested array (nested array) and a co-prime array (co-prime array), which are marked by asterisks and circles, respectively, and it can be seen that the performance of various sparse array structures is better than that of the nested array and the co-prime array. Assuming that there is only one target, DOA is-14 °, fig. 3 shows the DOA estimation results of the above four sparse array structures when using the WCSAB method. As can be seen from the figure, the optimal sparse array and the sparse array obtained by the algorithm can accurately estimate the target DOA, and the nested array and the co-prime array have errors in estimation.
In simulation 2, the performance of two wideband DOA estimation methods, WCSAB and wct (wideband Capon technique), were compared. The WCT method belongs to one of ISSM, and obtains a corresponding DOA estimation result by using a Capon method for each narrow-band signal, and then averages all the results to obtain a final estimation result. Fig. 4 considers the single target case with a target DOA of 10 °, and fig. 5 considers the dual target case with target DOAs of 6 ° and 10 °. In the figure, the solid line represents a full array structure at intervals of λ/2, the dotted line represents an optimal sparse array structure, and the dotted line represents a sparse array structure obtained according to an algorithm. As can be seen from fig. 4 and 5, compared to the full array structure, the two sparse arrays have narrower main lobe widths, i.e. the sparse arrays have higher resolution. But sparse arrays result in higher sidelobes and this problem is more severe in the dual target case. By comparison, it can be found that the side lobe of the WCSAB method is lower than the WCT. It can also be seen from fig. 5 that the WCSAB method can accurately estimate the target DOA in both sparse arrays, while the WCT estimates have errors. By contrast, the WCSAB method performs better.
Claims (2)
1. A MIMO radar broadband DOA calculation method based on sparse array in clutter environment comprises the following steps:
Step 1: let the position of the transmitting antenna be determined, and the feasible region for placing the receiving antenna be [0, Dr]For simplicity of analysis, the feasible fields are separated by an interval ΔrDiscretization to NrA plurality of grid points, and N receiving antennas disposed on some of the grid points, N < Nr;
Step 2: establishing an MIMO radar echo signal model to obtain echo signal time domain sampling datan=1,...,NrAnd p 1., L, where p represents a time domain snapshot and L is a snapshot number;
And step 3: for received signalPerforming an L-point discrete Fourier transform to obtain frequency domain data, i.e.
And N isrThe data of the grid points being represented in vector form, i.e.Wherein p 1, L1;
and 4, step 4: discretizing the target angle observation area intoG grid points theta1,...,θGK < G, where K represents the number of targets, the signal model is represented in sparse form:
y[l]=Φ[l]x+c[l]+u[l]
WhereinWhere a isr(θ,fl) Indicating a received steering vector, at(θ,fl) Representing the transmit steering vector, s [ l ]]Representing a frequency domain transmitted signal, x ═ x1,...,xG]TIs K sparse, i.e. x has only K non-zero elements, and the values and positions of the non-zero elements are the target reflection coefficient and DOA, c [ l [ ]]Represents clutter u [ l ]]Representing noise;
And 5: will beam form weight vector wg,lAct on y [ l]the beam forming output result is obtained:
Will r isg,l(G1., G, L ═ 1., L) is expressed as a vector:
r=[r1,1,...,rG,1,...,r1,L,...,rG,L]T
=WrΦx+Wrc+Wru
Wherein the weight matrix Wr=Diag{W1,...,Wl,...,WLIs a block diagonal matrix;
And is provided withΦ=[ΦT[1],...,ΦT[L]]T,c=[cT[1],...,cT[L]]TRepresents clutter, u ═ uT[1],...,uT[L]]TRepresenting noise;
Step 6: reconstructing a sparse vector x by basis pursuit desiccation based on the CS theory;
Wherein eta is equal to or greater than 0 is a regularization parameter;
And 7: to the solution obtained in step 6The values of the elements in the table are sorted from big to small, and the corresponding lattice point of each sorted element is expressed as { theta [ [ theta ] ](1),...,θ(G)Then the DOA estimation result can be expressed as;
and 8: based on minimum Bayes mean square errorSolving for optimal WrThe following optimization problem is established
s.t.Wr=Diag{W1,...,WL}
||W||0=N
W=[w1,1,...,w1,L,...,wG,1,...,wG,L]
Wherein the DOA vector theta of the real targetTis random in nature and is not only easy to be recognized,Is expressed in the pair thetaTIn the hope of expectation,Denotes thetaTMean square error of DOA estimation, w, when determinedg,lRepresents a weight vector, where G1.., G, L1., L;
And step 9: the problem provided by the step 8 is solved in an optimized way to obtain the optimal Wr。
2. The method according to claim 1, wherein the method for calculating the MIMO radar wideband DOA based on the sparse array in the clutter environment comprises the following specific steps:
Step 9.1: initialization: the number of iterations j is 1,According to the formulaComputing beamforming weightswhereinRc(fl) Is a clutter c [ l ]]The covariance matrix of (a); during each iteration, a set of lattice point selection vectors z is randomly generated1,...,zαFor a given z;
Step 9.2: repeating the iterative process from step 9.3 to step 9.6:
Step 9.3: randomly generating a set of lattice point selection vectors z1,...,zα};
step 9.4: according to the formula wg,l=z⊙ξg,lComputingAnd form
According to the formulaCalculating rg,lR is tog,lExpressed as a vector r;
Will be provided withSubstituting r into the formulaObtaining x reconstruction resultsAnd target DOA estimation result
According to the formulaTo obtain
Step 9.5: obtaining based on minimum BMSE
Step 9.6: according toBased onUpdatingGet the corresponding weight valueAnd let j equal j + 1; wherein Is Dc [ l ]]The covariance matrix of (a) is determined,
Step 9.7: when in useThen, iteration stops and the optimal antenna selection is outpute0Is a threshold value set in advance.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910646812.6A CN110579737B (en) | 2019-07-17 | 2019-07-17 | Sparse array-based MIMO radar broadband DOA calculation method in clutter environment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910646812.6A CN110579737B (en) | 2019-07-17 | 2019-07-17 | Sparse array-based MIMO radar broadband DOA calculation method in clutter environment |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110579737A true CN110579737A (en) | 2019-12-17 |
CN110579737B CN110579737B (en) | 2022-10-11 |
Family
ID=68811087
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910646812.6A Expired - Fee Related CN110579737B (en) | 2019-07-17 | 2019-07-17 | Sparse array-based MIMO radar broadband DOA calculation method in clutter environment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110579737B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111693975A (en) * | 2020-05-29 | 2020-09-22 | 电子科技大学 | MIMO radar sparse array design method based on deep neural network |
CN112016209A (en) * | 2020-08-28 | 2020-12-01 | 哈尔滨工业大学 | Distributed nested circular array comprehensive array arrangement method based on ant colony algorithm |
CN113673419A (en) * | 2021-08-19 | 2021-11-19 | 西北工业大学 | Beam domain fast sparse Bayesian direction estimation method suitable for strong interference environment |
Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100265799A1 (en) * | 2007-11-01 | 2010-10-21 | Volkan Cevher | Compressive sensing system and method for bearing estimation of sparse sources in the angle domain |
CN102967852A (en) * | 2012-11-29 | 2013-03-13 | 电子科技大学 | Method for generating multi-input multi-output over-horizon (MIMO-OTH) radar waveforms based on digital signal processor (DSP) sequences |
CN103353595A (en) * | 2013-06-18 | 2013-10-16 | 西安电子科技大学 | Meter wave radar height measurement method based on array interpolation compression perception |
US20150198713A1 (en) * | 2014-01-13 | 2015-07-16 | Mitsubishi Electric Research Laboratories, Inc. | Method and System for Through-the-Wall Imaging using Compressive Sensing and MIMO Antenna Arrays |
US20160238694A1 (en) * | 2015-02-16 | 2016-08-18 | Panasonic Intellectual Property Management Co., Ltd. | Radar device |
CN105974366A (en) * | 2016-04-29 | 2016-09-28 | 哈尔滨工程大学 | Four-order cumulant sparse representation-based MIMO (multiple-input-multiple-output) radar direction of arrival estimation method under mutual coupling condition |
CN106291540A (en) * | 2016-09-14 | 2017-01-04 | 河北省电力勘测设计研究院 | A kind of multiple-input and multiple-output GPR backwards projection target imaging method estimated based on DOA |
CN106501785A (en) * | 2016-09-13 | 2017-03-15 | 深圳大学 | A kind of sane sparse recovery STAP methods and its system based on alternating direction multiplier method |
CN106772225A (en) * | 2017-01-20 | 2017-05-31 | 大连大学 | Beam Domain DOA based on compressed sensing estimates |
US20170315221A1 (en) * | 2016-05-01 | 2017-11-02 | Technion Research & Development Foundation Ltd. | Target recovery in multiple input multiple output (mimo) radar system |
CN107479053A (en) * | 2017-09-21 | 2017-12-15 | 电子科技大学 | Carrier-borne MIMO radar based on STAP is steadily and surely launched with receiving co-design method |
US20180128892A1 (en) * | 2016-11-09 | 2018-05-10 | Raytheon Company | Systems and Methods for Direction Finding Using Compressive Sensing |
CN108802705A (en) * | 2018-04-24 | 2018-11-13 | 深圳大学 | It is a kind of based on sparse space-time adaptive processing method and system |
CN108828551A (en) * | 2018-08-28 | 2018-11-16 | 中国人民解放军空军工程大学 | A kind of compressed sensing based flexible MIMO radar compound target DOA estimation method |
CN108957388A (en) * | 2018-05-21 | 2018-12-07 | 南京信息工程大学 | A kind of MIMO radar coherent DOA estimation method based on covariance matching SL0 algorithm |
CN109061554A (en) * | 2018-06-26 | 2018-12-21 | 哈尔滨工程大学 | A kind of target angle of arrival estimation method updated based on spatial spreading grid dynamic |
CN109407045A (en) * | 2018-10-10 | 2019-03-01 | 苏州大学 | A kind of non-homogeneous sensor array broadband signal Wave arrival direction estimating method |
CN109901148A (en) * | 2019-03-21 | 2019-06-18 | 西安电子科技大学 | Broadband signal DOA estimation method based on covariance matrix rarefaction representation |
-
2019
- 2019-07-17 CN CN201910646812.6A patent/CN110579737B/en not_active Expired - Fee Related
Patent Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100265799A1 (en) * | 2007-11-01 | 2010-10-21 | Volkan Cevher | Compressive sensing system and method for bearing estimation of sparse sources in the angle domain |
CN102967852A (en) * | 2012-11-29 | 2013-03-13 | 电子科技大学 | Method for generating multi-input multi-output over-horizon (MIMO-OTH) radar waveforms based on digital signal processor (DSP) sequences |
CN103353595A (en) * | 2013-06-18 | 2013-10-16 | 西安电子科技大学 | Meter wave radar height measurement method based on array interpolation compression perception |
US20150198713A1 (en) * | 2014-01-13 | 2015-07-16 | Mitsubishi Electric Research Laboratories, Inc. | Method and System for Through-the-Wall Imaging using Compressive Sensing and MIMO Antenna Arrays |
US20160238694A1 (en) * | 2015-02-16 | 2016-08-18 | Panasonic Intellectual Property Management Co., Ltd. | Radar device |
CN105974366A (en) * | 2016-04-29 | 2016-09-28 | 哈尔滨工程大学 | Four-order cumulant sparse representation-based MIMO (multiple-input-multiple-output) radar direction of arrival estimation method under mutual coupling condition |
US20170315221A1 (en) * | 2016-05-01 | 2017-11-02 | Technion Research & Development Foundation Ltd. | Target recovery in multiple input multiple output (mimo) radar system |
CN106501785A (en) * | 2016-09-13 | 2017-03-15 | 深圳大学 | A kind of sane sparse recovery STAP methods and its system based on alternating direction multiplier method |
CN106291540A (en) * | 2016-09-14 | 2017-01-04 | 河北省电力勘测设计研究院 | A kind of multiple-input and multiple-output GPR backwards projection target imaging method estimated based on DOA |
US20180128892A1 (en) * | 2016-11-09 | 2018-05-10 | Raytheon Company | Systems and Methods for Direction Finding Using Compressive Sensing |
CN106772225A (en) * | 2017-01-20 | 2017-05-31 | 大连大学 | Beam Domain DOA based on compressed sensing estimates |
CN107479053A (en) * | 2017-09-21 | 2017-12-15 | 电子科技大学 | Carrier-borne MIMO radar based on STAP is steadily and surely launched with receiving co-design method |
CN108802705A (en) * | 2018-04-24 | 2018-11-13 | 深圳大学 | It is a kind of based on sparse space-time adaptive processing method and system |
CN108957388A (en) * | 2018-05-21 | 2018-12-07 | 南京信息工程大学 | A kind of MIMO radar coherent DOA estimation method based on covariance matching SL0 algorithm |
CN109061554A (en) * | 2018-06-26 | 2018-12-21 | 哈尔滨工程大学 | A kind of target angle of arrival estimation method updated based on spatial spreading grid dynamic |
CN108828551A (en) * | 2018-08-28 | 2018-11-16 | 中国人民解放军空军工程大学 | A kind of compressed sensing based flexible MIMO radar compound target DOA estimation method |
CN109407045A (en) * | 2018-10-10 | 2019-03-01 | 苏州大学 | A kind of non-homogeneous sensor array broadband signal Wave arrival direction estimating method |
CN109901148A (en) * | 2019-03-21 | 2019-06-18 | 西安电子科技大学 | Broadband signal DOA estimation method based on covariance matrix rarefaction representation |
Non-Patent Citations (6)
Title |
---|
AZIMI-SADJADI, MR; PEZESHKI, A; (...); HOHIL, M: "Wideband DOA estimation algorithms for multiple target detection and tracking using unattended acoustic sensors", 《CONFERENCE ON UNATTENDED/UNMANNED GROUND, OCEAN AND AIR SENSOR TECHNOLOGIES AND APPLICATIONS VI》 * |
JIAN-QI WU;WEI ZHU;BAIXIAO CHEN: "Compressed sensing techniques for altitude estimation in multipath conditions", 《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》 * |
MAHMOOD R. AZIMI-SADJADI;ALI PEZESHKI;NICHOLAS ROSEVEARE: "Wideband DOA estimation algorithms for multiple moving sources using unattended acoustic sensors", 《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》 * |
YU, Y; SUN, SQ AND PETROPULU, AP: "A capon beamforming method for clutter suppression in colocated compressive sensing based MIMO radars", 《CONFERENCE ON COMPRESSIVE SENSING II》 * |
何玄: "稀疏重构在天波雷达目标参数估计中的应用研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
薛会祥: "基于压缩感知理论的DOA估计算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111693975A (en) * | 2020-05-29 | 2020-09-22 | 电子科技大学 | MIMO radar sparse array design method based on deep neural network |
CN112016209A (en) * | 2020-08-28 | 2020-12-01 | 哈尔滨工业大学 | Distributed nested circular array comprehensive array arrangement method based on ant colony algorithm |
CN112016209B (en) * | 2020-08-28 | 2021-09-03 | 哈尔滨工业大学 | Distributed nested circular array comprehensive array arrangement method based on ant colony algorithm |
CN113673419A (en) * | 2021-08-19 | 2021-11-19 | 西北工业大学 | Beam domain fast sparse Bayesian direction estimation method suitable for strong interference environment |
CN113673419B (en) * | 2021-08-19 | 2024-05-28 | 西北工业大学 | Beam domain rapid sparse Bayesian azimuth estimation method suitable for strong interference environment |
Also Published As
Publication number | Publication date |
---|---|
CN110579737B (en) | 2022-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106353744B (en) | Multi-parameter combined estimation method based on bistatic FDA-MIMO radars | |
Häcker et al. | Single snapshot DOA estimation | |
Ender | A brief review of compressive sensing applied to radar | |
CN105652273B (en) | A kind of sparse imaging algorithm of MIMO radar based on mixing matching pursuit algorithm | |
CN106680815B (en) | MIMO radar imaging method based on tensor sparse representation | |
CN107576940B (en) | Low-complexity single-base MIMO radar non-circular signal angle estimation method | |
CN101813765B (en) | Noise suppression method based on inhomogeneous space solid array distributed SAR (Specific Absorption Rate) | |
CN108562866B (en) | Bistatic MIMO radar angle estimation method based on matrix filling | |
CN110579737B (en) | Sparse array-based MIMO radar broadband DOA calculation method in clutter environment | |
CN111693975A (en) | MIMO radar sparse array design method based on deep neural network | |
CN107064926B (en) | Bistatic MIMO radar angle estimation method under spatial color noise background | |
CN111610512B (en) | Frequency control array radar off-network target positioning method based on sparse Bayesian learning | |
CN113189592B (en) | Vehicle-mounted millimeter wave MIMO radar angle measurement method considering amplitude mutual coupling error | |
EP4050364A1 (en) | Radar detection using angle of arrival estimation based on scaling parameter with pruned sparse learning of support vector | |
Tao et al. | A knowledge aided SPICE space time adaptive processing method for airborne radar with conformal array | |
CN112763972A (en) | Sparse representation-based double parallel linear array two-dimensional DOA estimation method and computing equipment | |
EP4050373A1 (en) | Radar-based detection using angle of arrival estimation based on sparse array processing | |
CN115656957A (en) | FDA-MIMO target parameter estimation method for accelerating iterative convergence | |
Zhang et al. | Explicit Joint Resolution Limit of Range and Direction of Arrival Estimation for Phased-Array Radar | |
Hyder et al. | A new target localization method for bistatic FDA radar | |
Daniel et al. | MIMO radar waveform design for multiple extended targets using compressed sensing | |
CN114325560A (en) | Super-resolution target direction finding method for beam scanning radar | |
Liu et al. | Compressive sensing for very high frequency radar with application to low-angle target tracking under multipath interference | |
CN109752688B (en) | Method for calculating angle difference of adjacent information sources for sensor array system | |
CN112327264B (en) | Bistatic FDA-MIMO radar multidimensional parameter joint estimation method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20221011 |