CN103760527B - Single base MIMO radar coherent source Wave arrival direction estimating method - Google Patents
Single base MIMO radar coherent source Wave arrival direction estimating method Download PDFInfo
- Publication number
- CN103760527B CN103760527B CN201410032063.5A CN201410032063A CN103760527B CN 103760527 B CN103760527 B CN 103760527B CN 201410032063 A CN201410032063 A CN 201410032063A CN 103760527 B CN103760527 B CN 103760527B
- Authority
- CN
- China
- Prior art keywords
- matrix
- theta
- mrow
- mimo radar
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000001427 coherent effect Effects 0.000 title claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims abstract description 79
- 230000009466 transformation Effects 0.000 claims abstract description 17
- 238000001914 filtration Methods 0.000 claims abstract description 11
- 238000006243 chemical reaction Methods 0.000 claims abstract description 6
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 5
- 238000009499 grossing Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 claims description 4
- 238000003491 array Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 2
- 230000003595 spectral effect Effects 0.000 abstract description 3
- 230000007704 transition Effects 0.000 abstract 2
- 238000004088 simulation Methods 0.000 description 11
- 238000004422 calculation algorithm Methods 0.000 description 10
- 238000004364 calculation method Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
- 230000004304 visual acuity Effects 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/288—Coherent receivers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/08—Systems for measuring distance only
- G01S13/10—Systems for measuring distance only using transmission of interrupted, pulse modulated waves
- G01S13/26—Systems for measuring distance only using transmission of interrupted, pulse modulated waves wherein the transmitted pulses use a frequency- or phase-modulated carrier wave
- G01S13/28—Systems for measuring distance only using transmission of interrupted, pulse modulated waves wherein the transmitted pulses use a frequency- or phase-modulated carrier wave with time compression of received pulses
- G01S13/284—Systems for measuring distance only using transmission of interrupted, pulse modulated waves wherein the transmitted pulses use a frequency- or phase-modulated carrier wave with time compression of received pulses using coded pulses
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention belongs to Radar Technology field, disclose a kind of single base MIMO radar coherent source Wave arrival direction estimating method, can be used for goal orientation and the tracking of single base MIMO radar of General Cell stream shape.Implementation step is: step 1, writes out the steering vector of single base MIMO radar according to array manifold; Step 2, arranges M Wei Fandemeng vector, obtains transition matrix G; Step 3, obtains the steering vector after converting by smooth transformation matrix F; Step 4, matched filtering is carried out to the reception data of MIMO radar, reception data after matched filtering are designated as X (k), obtain data Y (the k)=FGX (k) after conversion, then form autocorrelation matrix by transition matrix G and level and smooth transformation matrix F; Step 5, carries out feature decomposition to autocorrelation matrix, forms noise subspace; Step 6, utilizes noise subspace to form space zero spectral function; Step 7, adopts polynomial rooting method to obtain azimuth of target.
Description
Technical Field
The invention belongs to the technical field of radars, and particularly relates to a method for estimating the direction of arrival of coherent sources of any single-base MIMO radar, which can be used for target orientation and tracking of the single-base MIMO radar with any array manifold.
Background
The multi-input multi-output radar (namely the MIMO radar) fully utilizes the waveform diversity gain to improve the spatial resolution, improve the parameter estimation precision and increase the maximum number of the targets which can be positioned by the system. In recent years, MIMO radar has been developed very rapidly, wherein a direction of arrival estimation (i.e., DOA) method is a major issue in the research of MIMO radar. In various array structures, the uniform linear array is the foundation of the research of the current MIMO radar DOA algorithm due to the simple structure and easy realization, and various rapid DOA algorithms can be adopted. However, one-dimensional uniform lines can only provide 180 degrees of unambiguous azimuthal information, and the resolution in different directions is different, and sometimes the need for a particular system makes the array likely not to be uniform. For the MIMO radar, a space beam is not formed during transmission, an array element directional diagram covers the whole space domain and is used for detecting the whole space information of 360 degrees, and at the moment, the linear array is adopted, so that the ambiguity can occur in the range of 360 degrees, and the requirement cannot be met. Therefore, estimation of the direction of arrival for any array manifold is a problem that needs to be studied.
At present, the most basic angle super-resolution estimation method adopts a multi-signal classification MUSIC method, which can work under any array manifold, but needs to carry out full-space domain search on angles, and has very large computation amount. In addition, the subspace-like high-resolution DOA estimation algorithm represented by the MUSIC algorithm has good resolution performance and high estimation accuracy for incoherent spatial sources, but the resolution performance thereof gradually deteriorates as the degree of correlation between the spatial sources increases.
The traditional array space smoothing algorithm is a common preprocessing method for processing coherent sources, and then a sub-array covariance algorithm and a weighted space smoothing algorithm appear on the basis of a forward and backward space smoothing technology, and the improved algorithms reduce the aperture loss of the array to a certain degree, improve the resolving power of the coherent sources, but bring larger calculation amount. The various smoothing-technology-based solution coherence methods mentioned above are all based on a uniform linear array structure, and cannot process any array manifold; in addition, the algorithms obtain the decoherence capability by sacrificing the effective aperture of the array, and the resolution capability of the coherent source is greatly reduced after the decoherence.
Disclosure of Invention
In view of the above shortcomings in the prior art, the present invention aims to provide a method for estimating direction of arrival of coherent source of monostatic MIMO radar, which can be used for coherent source decorrelation processing of MIMO radar arrays of any array manifold, and can realize monostatic MIMO radar coherent source DOA estimation of any array manifold while reducing the amount of computation by using polynomial root-finding algorithm instead of spectral peak search algorithm.
The technical idea for realizing the invention is as follows: firstly, decomposing the MIMO radar steering vector into a structure of multiplying a matrix and a Van der Waals vector by using an array interpolation technology; then, the thought of space smoothing is used for reference, and the coherent source is subjected to decoherence; and finally, a polynomial root-finding method is adopted to replace a spectral peak searching mode in the traditional method to obtain target angle information, so that the quick estimation of the direction of arrival is realized.
The technical scheme of the invention is a coherent source direction-of-arrival estimation method of a single-base MIMO radar, which is characterized by comprising the following steps:
step 1, obtaining an array steering vector A (theta) according to an array manifold of the single-base MIMO radar, wherein the array has N array elements, and the polar coordinate of each array element is (r)n,βn) N is 1, 2,. N; n is the number of array elements;
step 2, setting an M-dimension vandermonde vector B (theta), and setting the M-dimension vandermonde vector B (theta) A according to a formula GH(θ)[A(θ)AH(θ)]-1Solving a conversion matrix G; wherein, it is madeSelecting the dimension M of the Van der Monte guide vector B (theta) to be approximately equal to 4kR, k to be 2 pi/lambda, wherein lambda is the wavelength of the emission signal;
step 3, performing similar space smoothing processing on the van der waals guide vector B (theta), wherein if F is a smooth transformation matrix, the smoothed guide vector C (theta) is FB (theta);
step 4, performing matched filtering on the received data of the MIMO radar, recording the matched filtered received data as X (k), converting the matrix G and the smooth transformation matrix F to obtain a transformed data matrix Y (k) ═ FGX (k), and forming an autocorrelation matrix R by using the data matrix Y (k)y;
Step 5, for the autocorrelation matrix RyPerforming feature decomposition to obtain noise subspace
Step 6, utilizing noise subspaceForming a MUSIC spatial zero spectrum function:
wherein | | | purple hair2The expression is given in the 2-norm,representation matrixThe conjugate transpose of (1);
step 7, let z equal to ejθThen f (θ) can be converted to:
and f (z) is equal to 0, and a polynomial root solving method is adopted to obtain a complex exponential form z of P azimuth angles which are closest to a unit circle and have amplitudes less than 1pP ═ 1, 2.. P; wherein P is the number of targets in the same range gate, thetapAzimuth for the p-th target;
step 8, using the complex exponential form z of the acceptance anglepObtaining the p-th acceptance azimuth angle thetap:
θp=angle(zp)
Wherein, angle () represents the phase angle;
p=1,2...P。
and solving the azimuth angles of the P targets, namely finishing the estimation of the direction of arrival of the MIMO radar coherent source.
Compared with the prior art, the invention has the following advantages: (1) the existing DOA estimation method independent of the array manifold is an MUSIC spectrum estimation method, but the method cannot process coherent source signals, the existing subspace type coherent source processing method has requirements on the array manifold, the array is required to be a linear uniform array or a two-dimensional linear uniform array, and the irregular array manifold cannot be processed; (2) the method capable of processing coherent source signals is a maximum likelihood method at present, but the method needs multidimensional search, the calculation amount is large, and particularly under the condition of multiple sources, the calculation amount is exponentially increased, and the engineering realization is difficult. The invention can change the MIMO radar steering vector into a Van der Mongolian form, and the DOA estimation adopts a polynomial root-solving form, thereby greatly reducing the operation amount.
Drawings
The present invention will be described in further detail with reference to the accompanying drawings and specific embodiments.
FIG. 1 is a flow chart of an implementation of a coherent source direction of arrival estimation method of a monostatic MIMO radar of the invention;
FIG. 2 is a graph of simulation results obtained using the present invention for 100 independent experiments on a target;
FIG. 3 is a graph of the mean square error of the target angle as a function of the SNR for locating a target 1 using the present invention;
fig. 4 is a graph of the mean square error of the target angle as a function of the signal-to-noise ratio SNR when the target 2 is positioned using the present invention.
Detailed Description
Referring to fig. 1, the method for estimating the direction of arrival of coherent source of single-base MIMO radar of the present invention comprises the following steps:
step 1, obtaining the array steering vector A (theta) according to the array manifold of the single-base MIMO radar.
The single-base MIMO radar array has N array elements (the N array elements are shared by transmitting and receiving), and the polar coordinate of each array element can be obtained according to the array element distribution as (r)n,βn) N is 1, 2,. N, N is an array element number;
according to the polar coordinates (r) of each array elementn,βn) Obtaining the transmitting guide vector a of the MIMO radar arrayt(theta) and the received steering vector ar(θ) is:
where k is 2 pi/λ, λ is the wavelength of the emitted signal, rnAnd betanThe distance between the nth array element and the origin of the coordinate and the angular position of the nth array element relative to the x axis in the polar coordinate are respectively, theta is an azimuth angle]TRepresenting a matrix transposition;
according to the transmission steering vector at(theta) and the received steering vector ar(theta) obtaining a MIMO radar arrayThe column steering vector A (θ) is:
in the formula,representing a Kronecker product operation;
step 2, setting an M-dimension vandermonde vector B (theta), and setting the M-dimension vandermonde vector B (theta) A according to a formula GH(θ)[A(θ)AH(θ)]-1And solving a conversion matrix G.
2a) Obtaining dimension M of Van der Monte vector B (theta)
Order toN is 1, 2, N is the number of array elements, the dimension M of the van der mond guide vector B (θ) is selected to be approximately equal to 4kR, k is 2 pi/λ, and λ is the wavelength of the transmitted signal;
neglecting the truncation error, the van der waals steering vector B (θ) is:
wherein, the [ alpha ], [ beta ]]TRepresenting a matrix transpose and j an imaginary number.
2b) Obtaining a transformation matrix G
Taking J angles at equal intervals in the observation range, namely: theta1,θ2,...,θJ;
According to theta1,θ2,...,θJRespectively form a steering vector matrixAnd vandermonde matrix
According to a steering vector matrixAnd vandermonde matrixObtaining a conversion matrix G by using a least square method:
wherein,representation matrixConjugate transpose of (c) ()-1Representing the matrix inversion.
And step 3, performing similar space smoothing processing on the van der waals guide vector B (theta), wherein if F is a smooth transformation matrix, the smoothed guide vector is represented as C (theta) being FB (theta).
Dividing a van der waals guide vector B (θ) (equivalent to a uniform linear array) into L mutually overlapped sub-arrays, wherein the dimension of each sub-array is Q-M-L + 1; wherein M is the dimension of B (θ);
constructing a smooth transformation matrix according to Q: fl=[0Q×(l-1)|IQ|0Q×(L-l)]Wherein L is 1, 2Q×(l-1)Denotes the zero matrix of Q × (l-1), IQDenotes an identity matrix of order Q, 0Q×(l-1)Represents zero matrix of Q × (L-L);
the smooth transformation matrix F can be obtained by:
the transformed steering vector C (θ) can be obtained by:
C(θ)=FB(θ)
step 4, performing matched filtering on the received data of the MIMO radar, recording the matched filtered received data as X (k), converting the matrix G and the smooth transformation matrix F to obtain a transformed data matrix Y (k) ═ FGX (k), and forming an autocorrelation matrix R by using the data matrix Y (k)y;
4a) Matched filtering of received data for radar
Assuming that the radar cross-sectional areas RCS of the targets are the same, the received signals are:
wherein P is the number of targets in the same range gate, thetapAzimuth of the p-th target, fdpIs the doppler frequency (P ═ 1, 2,. P) of the pth target, and V is the receive noise matrix of the MIMO radar array.
And performing matched filtering on the received data of the MIMO radar array, wherein the received data after matched filtering is recorded as X (k).
X(k)=E(SXr)
Where E () represents the mathematical expectation, S ═ S1 T,...Si T...,SN T]TRepresenting a matrix of N transmit signals,Siwhich represents the ith transmission signal, is transmitted,denotes SiTranspose of (i ═ 1.. N);
4b) obtaining a transformed data matrix Y (k) by converting the matrix G and the smooth transformation matrix F
Y(k)=FGX(k)
4c) Then, the transformed data matrix Y (k) is used to form an autocorrelation matrix Ry
Where K is the number of sampling points, YH(k) Denotes the conjugate transpose of Y (k).
Step 5, for the autocorrelation matrix RyPerforming feature decomposition to obtain noise subspace
For autocorrelation matrix RyPerforming characteristic decomposition to obtain a series of characteristic values lambdamAnd its corresponding feature vector em(m ═ 1, 2,. Q), the characteristic value λ is calculatedmCarry out the order from big to small lambda1≥λ2≥...≥λQ,e1,...,eQSelecting P large eigenvalues lambda for corresponding eigenvectors1,...,λpP is the number of targets, and the eigenvectors e corresponding to the eigenvalues are taken1,e2,...,epForming a signal subspace
Then the noise subspaceSatisfies the following conditions:
wherein I is an identity matrix with dimension Q.
Step 6, utilizing noise subspaceForming a MUSIC spatial zero spectrum function:
wherein | | | purple hair2The expression is given in the 2-norm,representation matrixThe conjugate transpose of (c).
Step 7, let z equal to ejθThen f (θ) can be converted to:
and f (z) is equal to 0, and a polynomial root solving method is adopted to obtain a complex exponential form z of P azimuth angles which are closest to a unit circle and have amplitudes less than 1p(p=1,2...P)。
Step 8, using the complex exponential form z of the acceptance anglepFinding the azimuth angle theta of the p-th targetp;
θp=angle(zp)
Wherein, angle () represents the phase angle;
p=1,2...P
and solving the azimuth angles of the P targets, namely finishing the estimation of the direction of arrival of the MIMO radar coherent source.
The effects of the present invention are further illustrated by the following simulation test.
(1) Simulation conditions are as follows:
according to the derivation steps of the invention, the method has no limitation of the array manifold and is suitable for any array manifold. Without loss of generality, assuming that the MIMO radar array is a uniform circular array, the wavelength of a transmitting signal is 0.75M, the array element number is 10, the array radius is 0.45M, the transmitting signal is in the form of a phase coding signal orthogonal to a carrier frequency, two coherent signal sources are arranged, the azimuth angle of a signal source 1 is 22 degrees, the azimuth angle of a signal source 2 is 45 degrees, the repetition period number of a receiving pulse is 100 degrees, the signal-to-noise ratio SNR is 20dB, the van der monster steering vector dimension M is 11 degrees, the length of a spatial smoothing sub-array is Q9 degrees, 100 independent monte carlo experiments are performed, and the mean square error calculation of a signal source p adopts a formulaWhereinAs an estimate of the azimuth angle of the signal source p, thetapIs the azimuth angle of the signal source p.
(2) Simulation content:
simulation 1, the target orientation angle is subjected to target positioning simulation by adopting the method, and the results of 100 independent simulation experiments on the target are shown in fig. 2. As can be seen from FIG. 2, the invention can realize the fast direction-of-arrival estimation of the single-base MIMO radar of any array manifold.
Simulation 2, when the signal source 1 is positioned by adopting the method, the simulation result of the variation of the target angle mean square error along with the signal-to-noise ratio SNR is shown in FIG. 3;
simulation 3, when the signal source 2 is positioned by adopting the method, the simulation result of the variation of the target angle mean square error along with the signal-to-noise ratio SNR is shown in figure 4.
It can be seen from fig. 3 and 4 that the mean square error of the target angle estimation decreases as the signal-to-noise ratio SNR increases, and the mean square error of the target angle estimation can reach below 0.1 ° when the signal-to-noise ratio is 5dB, and the positioning accuracy is high, indicating that the present invention is practical.
In summary, the invention can realize the estimation of the direction of arrival of the coherent source of the single-base multiple-input multiple-output (MIMO) radar with any array manifold, and has high positioning precision.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present invention without departing from the spirit and scope of the invention. Thus, it is intended that the present invention cover the modifications and variations of this invention provided they come within the scope of the appended claims and their equivalents.
Claims (4)
1. A coherent source direction-of-arrival estimation method for a single-base MIMO radar is characterized by comprising the following steps:
step 1, obtaining an array steering vector A (theta) according to an array manifold of the single-base MIMO radar, wherein the array has N array elements, and the polar coordinate of each array element is (r)n,βn) N is 1, 2, … N, N is array element number;
step 2, setting an M-dimension vandermonde vector B (theta), and setting the M-dimension vandermonde vector B (theta) A according to a formula GH(θ)[A(θ)AH(θ)]-1Solving a conversion matrix G; wherein, it is madeSelecting the dimension M of the Van der Monte guide vector B (theta) to be approximately equal to 4kR, k to be 2 pi/lambda, wherein lambda is the wavelength of the emission signal;
step 3, performing similar space smoothing processing on the van der mond guide vector B (theta), wherein F is a smooth transformation matrix, and the smoothed guide vector C (theta) is FB (theta);
step 4, performing matched filtering on the received data of the MIMO radar, recording the matched filtered received data as X (k), obtaining a transformed data matrix Y (k) ═ FGX (k) by converting the matrix G and the smooth transformation matrix F, and forming an autocorrelation matrix R by using the data matrix Y (k)y;
Step 5, for the autocorrelation matrix RyPerforming feature decomposition to obtain noise subspace
Step 6, utilizing noise subspaceForming a MUSIC spatial zero spectrum function:
wherein |)2The expression is given in the 2-norm,representation matrixThe conjugate transpose of (1);
step 7, let z equal to ejθThen f (θ) can be converted to:
and f (z) is equal to 0, and a polynomial root solving method is adopted to obtain a complex exponential form z of P azimuth angles which are closest to a unit circle and have amplitudes less than 1pP is 1, 2 … P; wherein P is the number of targets in the same range gate, thetapAzimuth for the p-th target;
step 8, using the complex exponential form z of the acceptance anglepFinding the azimuth angle theta of the p-th targetp:
θp=angle(zp)
Wherein, angle () represents the phase angle;
p=1,2…P
and solving the azimuth angles of the P targets, namely finishing the estimation of the direction of arrival of the MIMO radar coherent source.
2. The method of estimating the direction of arrival of a coherent source of a monostatic MIMO radar according to claim 1, wherein the M-dimensional vandermonde vector B (θ) and the transformation matrix G in step 2 are obtained by:
first, letrnThe distance between the nth array element and the origin of coordinates is defined, N is 1, 2, … N, N is the number of the array elements, the dimension M of the Van der Menu vector B (theta) is selected to be approximately equal to 4kR, k is 2 pi/lambda, and lambda is the wavelength of the emission signal;
neglecting truncation errors, the vandermonde vector B (θ) is:
wherein, the [ alpha ], [ beta ]]TRepresenting a matrix transpose, j represents an imaginary number,
secondly, taking J angles at equal intervals in the observation range, namely: theta1,θ2,…,θJ(ii) a According to theta1,θ2…,θJRespectively forming steering matricesAnd vandermonde matrix
According to a steering matrixAnd vandermonde matrixObtaining a conversion matrix G by using a least square method:
wherein,representation matrixConjugate transpose of (c) ()-1Representing the matrix inversion.
3. The method for estimating direction of arrival of coherent source of monostatic MIMO radar according to claim 1, wherein the smooth transformation matrix F in step 3 is obtained by the following steps:
firstly, dividing a van der waals guiding vector B (theta) into L mutually overlapped sub-arrays, wherein the dimension of each sub-array is Q-M-L + 1; wherein M is the dimension of B (θ);
next, a smooth transformation matrix is constructed from Q: fl=[OQ×(l-1)|IQ|OQ×(L-l)]Wherein L is 1, 2, … L, OQ×(l-1)Denotes the zero matrix of Q × (l-1), IQRepresenting an identity matrix of order Q, OQ×(L-1)Represents zero matrix of Q × (L-L);
finally, the smooth transformation matrix F can be obtained by:
。
4. the coherent source direction-of-arrival estimation method for monostatic MIMO radar according to claim 1, wherein the detailed substeps of step 4 are as follows:
4a) matched filtering of received data for radar
Assuming that the radar cross-sectional areas RCS of the targets are the same, the received signals are:
wherein P is the number of targets in the same range gate, thetapAzimuth of the p-th target, fdpIs the doppler frequency of the pth target, P1, 2, … P, V is the receive noise matrix of the MIMO radar array,
carrying out matched filtering on the received data of the MIMO radar array, and recording the received data after matched filtering as X (k);
X(k)=E(SXr)
wherein E () represents the mathematical expectation,representing a matrix, s, consisting of N transmitted signalsiWhich represents the ith transmission signal, is transmitted,denotes si1, i is 1 … N;
4b) obtaining a transformed data matrix Y (k) by converting the matrix G and the smooth transformation matrix F
Y(k)=FGX(k)
4c) Then, the transformed data matrix Y (k) is used to form an autocorrelation matrix Ry
Where K is the number of sampling points, YH(k) Denotes the conjugate transpose of Y (k).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410032063.5A CN103760527B (en) | 2014-01-23 | 2014-01-23 | Single base MIMO radar coherent source Wave arrival direction estimating method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410032063.5A CN103760527B (en) | 2014-01-23 | 2014-01-23 | Single base MIMO radar coherent source Wave arrival direction estimating method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103760527A CN103760527A (en) | 2014-04-30 |
CN103760527B true CN103760527B (en) | 2015-12-09 |
Family
ID=50527794
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410032063.5A Expired - Fee Related CN103760527B (en) | 2014-01-23 | 2014-01-23 | Single base MIMO radar coherent source Wave arrival direction estimating method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103760527B (en) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104991236B (en) * | 2015-06-18 | 2017-10-31 | 哈尔滨工程大学 | A kind of single base MIMO radar not rounded signal coherence source Wave arrival direction estimating method |
FR3045958B1 (en) * | 2015-12-16 | 2019-05-17 | Suez Groupe | METHOD FOR SPECTRALLY ESTIMATING THE SCALE OF A LIQUID ENVIRONMENT |
CN107255805B (en) * | 2017-08-01 | 2019-10-25 | 西安电子科技大学 | The prediction technique of radar target RCS based on weighted least-squares |
CN107576947B (en) * | 2017-08-08 | 2020-01-31 | 西安电子科技大学 | L-shaped array pair coherent information source two-dimensional direction of arrival estimation method based on time smoothing |
CN110389319B (en) * | 2019-07-22 | 2021-04-27 | 北京工业大学 | MIMO radar DOA estimation method based on low-altitude multipath condition |
CN110927661A (en) * | 2019-11-22 | 2020-03-27 | 重庆邮电大学 | Single-basis expansion co-prime array MIMO radar DOA estimation method based on MUSIC algorithm |
CN111781593A (en) * | 2020-07-09 | 2020-10-16 | 深圳大学 | Low data volume coherent signal DOA estimation method, device, equipment and medium |
CN113466796B (en) * | 2021-08-16 | 2023-05-26 | 电子科技大学 | Radar communication integrated method based on coherent phase modulation broadcast mode |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103364762A (en) * | 2013-07-13 | 2013-10-23 | 西安电子科技大学 | Estimation method for arriving direction of monostatic MIMO radar based on random array manifolds |
-
2014
- 2014-01-23 CN CN201410032063.5A patent/CN103760527B/en not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103364762A (en) * | 2013-07-13 | 2013-10-23 | 西安电子科技大学 | Estimation method for arriving direction of monostatic MIMO radar based on random array manifolds |
Non-Patent Citations (4)
Title |
---|
"Direction of arrival estimation for monostatic multiple-input multiple-output radar with arbitrary array structures";Y.Cao et al.;《Radar,Sonar & Navigation,IET》;20120831;第6卷(第7期);第679-686页 * |
"Fast DOA estimation for monostatic MIMO radar with arbitrary array configurations";Yunhe Cao et al.;《Radar (Radar), 2011 IEEE CIE International Conference on》;20111027;第1卷;第959-962页 * |
"一种单基地MIMO雷达多目标DOA估计算法研究";郭艺夺 等;《上海航天》;20120229;第29卷(第1期);第1-5,11页 * |
"基于多项式求根的双基地MIMO雷达多目标定位方法";谢荣 等;《电子与信息学报》;20100930;第32卷(第9期);第2197-2200页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103760527A (en) | 2014-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103760527B (en) | Single base MIMO radar coherent source Wave arrival direction estimating method | |
CN102707264B (en) | Estimating method of direction of arrival of bistatic MIMO (Multi-Input Multi-Output) radar based on circular array | |
CN108508423B (en) | Subarray digital sum and difference monopulse angle measurement method based on special-shaped array | |
CN106707257B (en) | MIMO radar Wave arrival direction estimating method based on nested array | |
CN111123192B (en) | Two-dimensional DOA positioning method based on circular array and virtual extension | |
CN103901417B (en) | Low-complexity space target two-dimensional angle estimation method of L-shaped array MIMO radar | |
CN103364762B (en) | Estimation method for arriving direction of monostatic MIMO radar based on random array manifolds | |
CN108303683B (en) | Single-base MIMO radar real-value ESPRIT non-circular signal angle estimation method | |
CN104991236B (en) | A kind of single base MIMO radar not rounded signal coherence source Wave arrival direction estimating method | |
Zhang et al. | Joint transmission and reception diversity smoothing for direction finding of coherent targets in MIMO radar | |
CN107576940B (en) | Low-complexity single-base MIMO radar non-circular signal angle estimation method | |
CN110927661A (en) | Single-basis expansion co-prime array MIMO radar DOA estimation method based on MUSIC algorithm | |
CN106483493B (en) | A kind of sparse double parallel linear array and estimating two-dimensional direction-of-arrival method | |
CN103852751B (en) | Based on the centralized MIMO radar waveform method for designing receiving Wave beam forming | |
CN106019234B (en) | The low computation complexity estimating two-dimensional direction-of-arrival method of inverted-L antenna battle array | |
CN106443610B (en) | MIMO radar receives and dispatches array mutual coupling error self-correcting method | |
CN107870315B (en) | Method for estimating direction of arrival of any array by using iterative phase compensation technology | |
CN109254272B (en) | Two-dimensional angle estimation method of concurrent polarization MIMO radar | |
CN108020812B (en) | Two-dimensional DOA estimation method based on special three-parallel line array structure | |
Zhang et al. | Two-dimensional direction of arrival estimation for coprime planar arrays via polynomial root finding technique | |
CN103353588B (en) | Two-dimensional DOA (direction of arrival) angle estimation method based on antenna uniform planar array | |
CN109765521B (en) | Beam domain imaging method based on subarray division | |
CN109061578B (en) | Notch directional diagram waveform comprehensive design method based on MIMO radar | |
CN104502904B (en) | Torpedo homing beam sharpening method | |
CN104020440B (en) | Interfere the two-dimentional direction of arrival estimation method of formula linear array based on L-type |
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 |
Granted publication date: 20151209 Termination date: 20200123 |
|
CF01 | Termination of patent right due to non-payment of annual fee |