CN108303683A - Single not rounded signal angle methods of estimation of base MIMO radar real value ESPRIT - Google Patents
Single not rounded signal angle methods of estimation of base MIMO radar real value ESPRIT Download PDFInfo
- Publication number
- CN108303683A CN108303683A CN201810081868.7A CN201810081868A CN108303683A CN 108303683 A CN108303683 A CN 108303683A CN 201810081868 A CN201810081868 A CN 201810081868A CN 108303683 A CN108303683 A CN 108303683A
- Authority
- CN
- China
- Prior art keywords
- signal
- esprit
- estimation
- real
- rounded
- 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
- 238000000034 method Methods 0.000 title claims abstract description 34
- 239000011159 matrix material Substances 0.000 claims abstract description 53
- 230000009467 reduction Effects 0.000 claims abstract description 21
- 238000010276 construction Methods 0.000 claims abstract description 13
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 13
- 238000001914 filtration Methods 0.000 claims abstract description 8
- 230000009466 transformation Effects 0.000 claims description 8
- 238000003491 array Methods 0.000 claims description 3
- 230000003750 conditioning effect Effects 0.000 claims description 3
- 230000010363 phase shift Effects 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 abstract description 43
- 238000005516 engineering process Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000000686 essence Substances 0.000 description 1
- 230000005713 exacerbation Effects 0.000 description 1
- JEIPFZHSYJVQDO-UHFFFAOYSA-N ferric oxide Chemical compound O=[Fe]O[Fe]=O JEIPFZHSYJVQDO-UHFFFAOYSA-N 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 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/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
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention belongs to Radar Technology fields, disclose a kind of list not rounded signal angle methods of estimation of base MIMO radar real value ESPRIT, and array element is received data carries out matched filtering with transmitting signal respectively, obtains observation data vector;Dimensionality reduction pretreatment is carried out to observation data, lower dimensional space is obtained and receives data vector;Data vector is received using the real value that the not rounded characteristic and Euler formula construction array apertures of signal double;Construct the invariable rotary relationship of the virtual array of aperture extension;The covariance matrix that extension receives data is calculated, Eigenvalues Decomposition is carried out to it, estimation obtains real-valued signal subspace;New real-valued signal subspace is defined, the invariable rotary equation of new real-valued signal subspace is solved;The DOA estimated values of target are calculated.The computation complexity of ESPRIT algorithms can be greatly lowered in the present invention while significantly improving DOA estimated accuracies, be suitable for low signal-to-noise ratio and low number of snapshots occasion.
Description
Technical field
The invention belongs to Radar Technology field more particularly to a kind of single base not rounded signal angles MIMO radar real value ESPRIT
Spend method of estimation.
Background technology
Currently, the prior art commonly used in the trade is such:Multiple-input and multiple-output (multiple input multiple
Output, MIMO) radar is a kind of new system radar to be got up based on MIMO development communication technologies.MIMO radar utilizes waveform point
The thought of collection emits mutually orthogonal waveform using multiple transmitting antennas, while receiving target using multiple reception antennas simultaneously
Reflect signal.Compared with traditional phased-array radar, MIMO radar has higher angular resolution and more degree of freedom, tool
There is better angle estimation performance.Direction of arrival (direction of arrival, DOA) estimation is MIMO radar parameter Estimation
An important research content.Invariable rotary sub-space technique (Estimate signal parameters via
Rotational invariance technique, ESPRIT) it is a kind of subspace class high resolution DOA estimation algorithm of classics.
By being utilized respectively the invariable rotary structure of MIMO radar emission array and receiving array, ESPRIT algorithms can apply to MIMO
In radar target DOA estimations.Studies have shown that using the not rounded characteristic of signal, the precision of radar parameter estimation can be significantly improved,
Improve estimation performance.Nearly ten years, numerous scholars are around the expansion further investigation of ESPRIT algorithms, it is proposed that various to be suitable for MIMO
The ESPRIT innovatory algorithms of radar.U-ESPRIT algorithms (Electronics Letters, 2012,48 (3):179-181) exist
Data covariance matrix real value will be received by unitary transformation on the basis of ESPRIT, reduce the computational complexity of algorithm, and
And improve the angle estimation performance under the conditions of low signal-to-noise ratio and low number of snapshots.But this method does not utilize the spy of transmitting signal
Point, therefore its asymptotic estimates performance is identical as ESPRIT algorithms.C-ESPRIT algorithms (Electronics Letters, 2010,
46(25):The virtual array that the characteristics of not rounded signal 1692-1694) is utilized construction array aperture doubles, can significantly improve
Angle estimation precision, but computation complexity is dramatically increased also with doubling for matrix dimensionality, is unfavorable for the real-time of algorithm for estimating
It realizes.RV-ESPRIT algorithms (Journal of Applied Remote Sensing, 2016,10 (2):025003) it is a kind of
The real value ESPRIT algorithms for the characteristics of not rounded signal is utilized, although it uses real value processing means, it calculates complicated
Degree is as the increase of array number is in a cube time growth, and when MIMO port numbers are larger, calculation amount is still considerable.
In conclusion problem of the existing technology is:
(1) existing majority ESPRIT algorithms do not make full use of the not rounded characteristic of transmitting signal, in low signal-to-noise ratio and low fast
Under the conditions of umber of beats, since the inaccuracy of subspace estimation can cause angle estimation precision low or even fail;(2) emit to utilize
The not rounded characteristic of signal, existing ESPRIT algorithms typically directly construct the virtual array that aperture doubles to improve target angle estimation
Precision, this certainly will cause the computation complexity of algorithm to steeply rise, and be unfavorable for the real-time implementation of algorithm.
Solve the meaning of above-mentioned technical problem:The present invention makes full use of the not rounded characteristic of transmitting signal, can improve
The angle estimation precision of ESPRIT algorithms solves ESPRIT algorithms performance severe exacerbation under the conditions of low signal-to-noise ratio and low number of snapshots
The problem of, it is its practical application based theoretical.The present invention can reduce the complexity of existing not rounded signal ESPRIT algorithms,
It provides efficient MIMO radar not rounded angle estimating method, accelerates the speed of target direction estimation, be conducive to ESPRIT algorithms
Real-time implementation promotes the practical application of DOA algorithm for estimating.
Invention content
In view of the problems of the existing technology, the present invention provides a kind of list not rounded letters of base MIMO radar real value ESPRIT
Number angle estimating method.
The invention is realized in this way a kind of list not rounded signal angle methods of estimation of base MIMO radar real value ESPRIT,
Array element is received data and is carried out with transmitted waveform by the list not rounded signal angle methods of estimation of base MIMO radar real value ESPRIT
Matched filtering obtains observation data vector, carries out dimensionality reduction pretreatment to observation data, obtains lower dimensional space and receive data vector;Profit
Data vector is received with the real value that the not rounded characteristic of signal and Euler formula construction array apertures double, the extension of construction aperture
The invariable rotary relationship of virtual array calculates the covariance matrix for receiving data, carries out Eigenvalues Decomposition, and estimation obtains signal subspace
Space;New real-valued signal subspace is defined, the invariable rotary equation of the real-valued signal subspace is solved, target is calculated
DOA。
Further, the list not rounded signal direction of arrival methods of base MIMO radar real value ESPRIT include the following steps:
Each reception data for receiving array element are carried out matched filtering with transmitted waveform respectively, obtain matching filter by step 1
Observation data vector x (t) after wave;
Step 2 selects dimensionality reduction transformation matrix U, carries out dimensionality reduction to observation data vector x (t), obtains the observation after dimensionality reduction
Data vector y (t)=Ux (t);
Observation data vector y (t) is decomposed into real part y by step 3 using Euler formulac(t) and imaginary part ys(t), it utilizes
The not rounded characteristic of signal connects the real and imaginary parts for observing data, and the real value that construction array aperture doubles receives data vector
yr(t);
Step 4 defines two selection matrix J1And J2, the invariable rotary relationship J of construction extension virtual array2Gr=J1Gr
Ω;
Step 5 calculates yr(t) data covariance matrix Ry, Eigenvalues Decomposition is carried out to it, estimation obtains real-valued signal
Subspace Us;
Step 6 defines new signal subspaceNew reality is solved using total least square method
The invariable rotary equation of value signal subspaceReal value matrix Ψ is calculated;
Step 7 carries out Eigenvalues Decomposition to real value matrix Ψ, obtains P characteristic value, and then the wave for obtaining P target reaches
Direction estimation.
Further, data vector is observed in the step 1 is:
X (t)=As (t)+n (t);
Wherein,Combine steering vector matrix for send-receive,
θ1,θ2,…,θPFor the direction of arrival of P target,It is oriented to and swears for emission array
Amount, M are transmitting antenna number,For receiving array steering vector, N is to receive
Antenna number,Operator is accumulated for Kronecker;S (t)=[s1(t),s2(t),…,sP(t)]TFor signal phasor;n(t)∈
CMN×1It is zero-mean, covariance matrix σ2The white complex gaussian noise vector of I.Under not rounded signal conditioning, s (t) can be indicated
For:
S (t)=Λ r (t);
Wherein, r (t) is not rounded signal and meets r (t)=r*(t),Table
Show the additional phase shift of the P signal.
Further, the observation data vector in the step 2 after dimensionality reduction is:
Y (t)=V1/2GΛr(t)+nT(t);
Wherein, G=[g (θ1),g(θ2),…,g(θP)],Ne=M+
N-1 is the effective array number of virtual line arrays;nT=V1/2FHN (t) is the noise vector after dimensionality reduction;For diagonal matrix, diag () is indicated
Element diagonal matrixization operates, and transformation matrix F is defined as:
Further, the reception data vector y that the aperture constructed in the step 3 doublesr(t) it is:
Wherein, yc(t) and ys(t) be respectively y (t) real and imaginary parts;
Gc=[gc(θ1),...,gc(θP)], gc(θp)=[cos βp,...,cos((Ne-1)πsinθp+βp)]T, Gs=[gs(θ1),...,
gs(θP)], gs(θp)=[sin βp,...,sin((Ne-1)πsinθp+βp)]T;For expanded noise vector, ns(t)
=Im [n (t)], nc(t)=Re [n (t)].
Further, the invariable rotary equation J of virtual array steering vector matrix is extended in the step 42Gr=J1GrΩ
In, selection matrix J1And J2It is defined as;
Wherein, T1And T2It is defined as:
The invariable rotary equation J2Gr=J1GrIn Ω,It is right
Angular moment battle array, diagonal element include the DOA information of target.
Further, the reception data vector y extended in step 5r(t) covariance matrix is Ry=E { y (t)ryr(t)H,
Its Eigenvalues Decomposition is:
Ry=UsΣsUs H+UnΣnUn H;
Wherein, ΣsFor by RyThe diagonal matrix that constitute of the big characteristic value of P, UsFor corresponding signal subspace;Σn
For by residue (2Ne- 1-P) diagonal matrix that constitutes of a small characteristic value, UnFor corresponding noise subspace.
Further, real-valued signal subspace new in step 6 is defined asNew real-valued signal subspace
Invariable rotary equation be
Further, the DOA estimated values of P target described in step 7 can be calculated by following formula:
Wherein, λ1,λ2,......,λpFor the P characteristic value of real value matrix Ψ,For the DOA of P target
Estimated value.
Another object of the present invention is to provide a kind of application list not rounded signals of base MIMO radar real value ESPRIT
The MIMO radar of angle estimating method.
In conclusion advantages of the present invention and good effect are:The present invention is converted using dimensionality reduction drops observation data
Dimension pretreatment, can substantially reduce operational data dimension on the whole, ensure that subsequent calculate of algorithm carries out in lower dimensional space, so
Euler formula construction real values are utilized to receive data vector afterwards so that it is subsequent to be calculated as real-valued calculation, therefore there is lower meter
Complexity is calculated, the real-time implementation of algorithm is conducive to.The present invention utilizes the real value of the characteristics of not rounded signal construction extension to receive data
Array aperture is expanded to original twice by vector, and the DOA that target is then carried out using the invariable rotary structure of array extending is estimated
Meter, can significantly improve the angle estimation precision of ESPRIT algorithms, be suitable for low signal-to-noise ratio and low number of snapshots occasion.Therefore, originally
Invention can provide the angle estimation precision significantly improved with lower computational complexity.
Description of the drawings
Fig. 1 is the list not rounded signal angle method of estimation streams of base MIMO radar real value ESPRIT provided in an embodiment of the present invention
Cheng Tu.
Fig. 2 be it is provided in an embodiment of the present invention carried out under the conditions of M=8, N=6, L=100, SNR=10dB 100 times it is imitative
The target DOA estimation value schematic diagram that true experiment obtains.
Fig. 3 is the present invention under the conditions of M=8, N=6, L=100, the root-mean-square error of angle estimation with signal-to-noise ratio change
Change relation curve schematic diagram.
Fig. 4 be the present invention under the conditions of M=8, N=6, SNR=10dB, the root-mean-square error of angle estimation is with number of snapshots
Variation relation curve synoptic diagram.
Fig. 5 be the present invention under the conditions of M=N, L=200, P=3, computational complexity shows with array number variation relation curve
It is intended to.
Specific implementation mode
In order to make the purpose , technical scheme and advantage of the present invention be clearer, with reference to embodiments, to the present invention
It is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, it is not used to
Limit the present invention.
As shown in Figure 1, the not rounded signal angle estimations of list base MIMO radar real value ESPRIT provided in an embodiment of the present invention
Method includes the following steps:
S101:Array element is received into data and obtains observation data vector with transmitted waveform progress matched filtering;
S102:Dimensionality reduction pretreatment is carried out to observation data, lower dimensional space is obtained and receives data vector;
S103:Data vector is received using the real value that the not rounded characteristic and Euler formula construction array apertures of signal double;
S104:Construct the invariable rotary relationship of the virtual array of aperture extension;
S105:The covariance matrix for receiving data is calculated, Eigenvalues Decomposition is carried out, estimation obtains signal subspace;
S106:New real-valued signal subspace is defined, the invariable rotary equation of the real-valued signal subspace is solved;
S107:The DOA of target is calculated.
The list not rounded signal angle methods of estimation of base MIMO radar real value ESPRIT provided in an embodiment of the present invention are specifically wrapped
Include following steps:
(1) MIMO radar emits mutually orthogonal pulse code signal using M transmitting antenna, it is assumed that far field space exists
P incoherent narrowband targets receive target echo using N number of reception antenna in receiving terminal, matched filter pair are used in combination
Each reception data carry out matched filtering, obtain observation data vector;
Involved observation data vector is:
X (t)=As (t)+n (t);
Wherein,Combine steering vector matrix for send-receive,
θ1,θ2,…,θPFor the direction of arrival of P target,It is oriented to and swears for emission array
Amount,For receiving array steering vector,Operator is accumulated for Kronecker;s
(t)=[s1(t),s2(t),…,sP(t)]TFor signal phasor;n(t)∈CMN×1It is zero-mean, covariance matrix σ2The multiple height of I
This white noise vector.Under not rounded signal conditioning, s (t) can be expressed as:
S (t)=Λ r (t);
Wherein, r (t) is not rounded signal and meets r (t)=r*(t),Table
Show the additional phase shift of the P signal.Therefore, observation data vector can be expressed as:
X (t)=A Λ r (t)+n (t);
(2) dimensionality reduction transformation matrix U=V is chosen-1/2FH, dimensionality reduction is carried out to observation data vector x (t), obtains the number after dimensionality reduction
According to matrix y (t)=Ux (t);
Observation data vector after involved dimensionality reduction is:
Y (t)=V1/2GΛr(t)+nT(t);
Wherein, G=[g (θ1),g(θ2),…,g(θP)],Ne=M+
N-1 is the effective array number of virtual line arrays;
For diagonal matrix, diag () indicates that the operation of element diagonal matrixization, transformation matrix F are defined as:
nT=V1/2FHN (t) is the noise vector after dimensionality reduction.
(3) Euler formula are utilized, observation data vector y (t) is decomposed into real part yc(t) and imaginary part ys(t), signal is utilized
Not rounded characteristic, the reception data vector y that two parts series configuration aperture is doubledr(t);
The real part y of involved data vector y (t)c(t) and imaginary part ys(t) it is respectively:
The reception data vector y of involved extensionr(t) it is:
Wherein,To extend the steering vector matrix of virtual array, Gc=
[gc(θ1),...,gc(θP)], gc(θp)=[cos βp,...,cos((Ne-1)πsinθp+βp)]T,
Gs=[gs(θ1),...,gs(θP)], gs(θp)=[sin βp,...,sin((Ne-1)πsinθp+βp)]T;For expanded noise vector, ns(t)=Im [n (t)], nc(t)=Re [n (t)].
(4) two selection matrix J are defined1And J2, the invariable rotary equation of construction extension virtual array steering vector matrix;
Involved selection matrix J1And J2For:
Wherein, T1And T2It is defined as:
The invariable rotary equation of involved extension virtual array steering vector matrix is:
J2Gr=J1GrΩ;
Wherein,Diagonal matrix, diagonal element include target
DOA information.
(5) it calculates extension and receives data vector yr(t) autocorrelation matrix Ry, Eigenvalues Decomposition is carried out to it, estimation obtains
Signal subspace Us;
Involved yr(t) autocorrelation matrix is Ry=E { y (t)ryr(t)H, Eigenvalues Decomposition can be expressed as:
Ry=UsΣsUs H+UnΣnUn H;
Wherein, ΣsFor by RyThe diagonal matrix that constitute of the big characteristic value of P, UsFor corresponding signal subspace;Σn
For by residue (2Ne- 1-P) diagonal matrix that constitutes of a small characteristic value, UnFor corresponding noise subspace.
(6) new real-valued signal subspace is calculatedNew letter is solved using least square method or total least square method
Real value matrix Ψ is calculated in the invariable rotary equation in work song space;
Involved new signal subspace isThe invariable rotary equation of involved signal subspace
For
(7) Eigenvalues Decomposition is carried out to Ψ, obtains P eigenvalue λ1,λ2,......,λp, and then target is calculated
DOA estimates.
The DOA of P involved target can be estimated to obtain by following formula:
Wherein,For the DOA estimated values of P target.
The application effect of the present invention is explained in detail with reference to emulation.
(1) simulated conditions and content
Consider the single base MIMO radar system being made of even linear array, emit array number M=8, receives array number N=6,
Each array element spacing is half-wavelength.Assuming that far field space, there are 3 incoherent narrowband targets, the azimuth of each target is respectively
θ1=100,θ2=150,θ3=200.To verify effectiveness of the invention, the present invention and U-ESPRIT algorithms, RV-ESPRIT are calculated
Method and C-ESPRIT algorithms are compared.The root-mean-square error of angle estimation is defined as:
Wherein, K is total Monte-Carlo experiment numbers,Indicate p-th of mesh in kth time Monte-Carlo experiments
Target DOA estimated values, θpFor the angle actual value of p-th of target.
(2) simulation result
1, MIMO radar target positioning performance
Fig. 2 is to carry out 100 emulation experiments under the conditions of M=8, N=6, L=100, SNR=10dB using the present invention and obtain
The target DOA estimation value arrived.From figure 2 it can be seen that multiple targets can accurately determine simultaneously using inventive algorithm
Position.
2, the root-mean-square error of MIMO radar angle estimation with signal-to-noise ratio variation relation
Fig. 3 is that the present invention carries out the angle that 500 Monte-Carlo are tested under the conditions of M=8, N=6, L=100
The root-mean-square error of estimation with signal-to-noise ratio variation relation curve.From figure 2 it can be seen that under Low SNR, this hair
The bright estimated accuracy with C-ESPRIT algorithms, RV-ESPRIT algorithms is superior to U-ESPRIT algorithms, wherein the present invention and C-
ESPRIT algorithms are apparent to the improvement of angle estimated accuracy.The estimated accuracy of each algorithm is improved with the increase of signal-to-noise ratio,
Inventive algorithm and C-ESPRIT algorithms asymptotic estimates performance having the same.
3, the root-mean-square error of MIMO radar angle estimation with number of snapshots variation relation
Fig. 4 is that the present invention carries out 500 Monte-Carlo under the conditions of M=8, N=6, SNR=10dB and tests
The root-mean-square error of angle estimation with number of snapshots variation relation curve.From figure 3, it can be seen that the present invention and C-ESPRIT are calculated
Method, the estimated accuracy of RV-ESPRIT algorithms are superior to U-ESPRIT algorithms.Since the present invention and C-ESPRIT algorithms are utilized
The angle estimation precision of the reception data that virtual aperture doubles, the two obtains large increase, and angle estimation precision is basic
It is close.
4, the computational complexity of MIMO radar angle estimation with dual-mode antenna number variation relation
From figure 5 it can be seen that the computational complexity (real value multiplication number) of each algorithm with the increase of array number and
Increase.Wherein, the computational complexity of U-ESPRIT algorithms, C-ESPRIT algorithms and RV-ESPRIT algorithms is with array number
Increase and steeply rises, and the computational complexity of inventive algorithm is slower with array number variation, and computational complexity is minimum.This
It is due to the present invention while to use dimensionality reduction transformation and real valueization operation, the operation for greatly reducing ESPRIT algorithms is complicated
Degree.
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, all essences in the present invention
All any modification, equivalent and improvement etc., should all be included in the protection scope of the present invention made by within refreshing and principle.
Claims (10)
1. a kind of list not rounded signal angle methods of estimation of base MIMO radar real value ESPRIT, which is characterized in that the list base
Array element is received data and is obtained with transmitted waveform progress matched filtering by the not rounded signal angle methods of estimation of MIMO radar real value ESPRIT
To observation data vector, dimensionality reduction pretreatment is carried out to observation data, lower dimensional space is obtained and receives data vector;Utilize the non-of signal
The real value that circle characteristic and Euler formula construction array apertures double receives data vector, the virtual array of structure aperture extension
Invariable rotary relationship calculates the covariance matrix for receiving data, carries out Eigenvalues Decomposition, and estimation obtains signal subspace;Definition
New real-valued signal subspace solves the invariable rotary equation of the real-valued signal subspace, the DOA of target is calculated.
2. the list not rounded signal angle methods of estimation of base MIMO radar real value ESPRIT as described in claim 1, feature exist
In the list not rounded method for estimating signal wave direction of base MIMO radar real value ESPRIT includes the following steps:
Each reception data for receiving array element are carried out matched filtering, after obtaining matched filtering by step 1 with transmitted waveform respectively
Observation data vector x (t);
Step 2 selects dimensionality reduction transformation matrix U, carries out dimensionality reduction to observation data vector x (t), obtains the observation data after dimensionality reduction
Vector y (t)=Ux (t);
Observation data vector y (t) is decomposed into real part y by step 3 using Euler formulac(t) and imaginary part ys(t), signal is utilized
Not rounded characteristic, the real and imaginary parts series connection of data will be observed, the real value that doubles of construction array aperture receives data vector yr
(t);
Step 4 defines two selection matrix J1And J2, the invariable rotary relationship J of construction extension virtual array2Gr=J1GrΩ;
Step 5 calculates yr(t) data covariance matrix Ry, Eigenvalues Decomposition is carried out to it, it is empty that estimation obtains real-valued signal
Between Us;
Step 6 defines new real-valued signal subspaceNew real-valued signal subspace is solved using total least square method
Invariable rotary equationReal value matrix Ψ is calculated;
Step 7 carries out Eigenvalues Decomposition to real value matrix Ψ, obtains its P characteristic value, and then the wave for estimating P target reaches
Direction.
3. the list not rounded signal angle methods of estimation of base MIMO radar real value ESPRIT as claimed in claim 2, feature exist
In observing data vector in the step 1 is:
X (t)=As (t)+n (t);
Wherein,Combine steering vector matrix, θ for send-receive1,
θ2,…,θPFor the direction of arrival of P target,For emission array steering vector, M
For transmitting antenna number,For receiving array steering vector, N is reception antenna
Number,Operator is accumulated for Kronecker;S (t)=[s1(t),s2(t),…,sP(t)]TFor signal phasor;n(t)∈CMN×1
It is zero-mean, covariance matrix σ2The white complex gaussian noise vector of I.Under not rounded signal conditioning, s (t) can be expressed as:
S (t)=Λ r (t);
Wherein, r (t) is not rounded signal and satisfaction Indicate P
The additional phase shift of a reflection signal.
4. the list not rounded signal angle methods of estimation of base MIMO radar real value ESPRIT as claimed in claim 2, feature exist
In the observation data vector in the step 2 after dimensionality reduction is:
Y (t)=V1/2GΛr(t)+nT(t);
Wherein, G=[g (θ1),g(θ2),…,g(θP)],Ne=M+N-1 is
The effective array number of virtual line arrays;nT=V1/2FHN (t) is the noise vector after dimensionality reduction;For diagonal matrix, diag () indicates member
Plain diagonal matrixization operation, transformation matrix F are defined as:
5. the list not rounded signal angle methods of estimation of base MIMO radar real value ESPRIT as claimed in claim 2, feature exist
In the reception data vector y that the aperture constructed in the step 3 doublesr(t) it is:
Wherein, yc(t) and ys(t) be respectively y (t) real and imaginary parts; Gc=
[gc(θ1),...,gc(θP)], gc(θp)=[cos βp,...,cos((Ne-1)πsinθp+βp)]T, Gs=[gs(θ1),...,gs
(θP)], gs(θp)=[sin βp,...,sin((Ne-1)πsinθp+βp)]T;For expanded noise vector, nc(t)=
Re [n (t)], ns(t)=Im [n (t)].
6. the list not rounded signal angle methods of estimation of base MIMO radar real value ESPRIT as claimed in claim 2, feature exist
In the invariable rotary equation J of extension virtual array steering vector matrix in the step 42Gr=J1GrIn Ω, selection matrix J1
And J2It is defined as:
Wherein, T1And T2It is defined as:
The invariable rotary equation J2Gr=J1GrIn Ω,For to angular moment
Battle array, diagonal element include the DOA information of target.
7. the list not rounded signal angle methods of estimation of base MIMO radar real value ESPRIT as claimed in claim 2, feature exist
In the reception data vector y extended in step 5r(t) covariance matrix is Ry=E { y (t)ryr(t)H, Eigenvalues Decomposition
For:
Ry=UsΣsUs H+UnΣnUn H;
Wherein, ΣsFor by RyThe diagonal matrix that constitute of the big characteristic value of P, UsFor corresponding signal subspace;ΣnIt serves as reasons
Residue (2Ne- 1-P) diagonal matrix that constitutes of a small characteristic value, UnFor corresponding noise subspace.
8. the list not rounded signal angle methods of estimation of base MIMO radar real value ESPRIT as claimed in claim 2, feature exist
In new real-valued signal subspace is defined as in step 6The invariable rotary side of new real-valued signal subspace
Cheng Wei
9. the list not rounded signal angle methods of estimation of base MIMO radar real value ESPRIT as claimed in claim 2, feature exist
In the DOA estimated values of P target described in step 7 can be calculated by following formula:
Wherein, λ1,λ2,......,λpFor the P characteristic value of real value matrix Ψ,For the DOA estimations of P target
Value.
10. the list not rounded signal angles of base MIMO radar real value ESPRIT described in a kind of application claim 1~9 any one are estimated
The MIMO radar of meter method.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810081868.7A CN108303683B (en) | 2018-01-29 | 2018-01-29 | Single-base MIMO radar real-value ESPRIT non-circular signal angle estimation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810081868.7A CN108303683B (en) | 2018-01-29 | 2018-01-29 | Single-base MIMO radar real-value ESPRIT non-circular signal angle estimation method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108303683A true CN108303683A (en) | 2018-07-20 |
CN108303683B CN108303683B (en) | 2021-10-22 |
Family
ID=62866621
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810081868.7A Expired - Fee Related CN108303683B (en) | 2018-01-29 | 2018-01-29 | Single-base MIMO radar real-value ESPRIT non-circular signal angle estimation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108303683B (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109407048A (en) * | 2018-11-12 | 2019-03-01 | 华南理工大学 | Underwater DOA estimation method and device based on non-circular signal and the adjustable battle array of angle |
CN109521393A (en) * | 2018-11-05 | 2019-03-26 | 昆明理工大学 | A kind of DOA estimation algorithm based on signal subspace revolving property |
CN109521392A (en) * | 2018-10-24 | 2019-03-26 | 华南理工大学 | Underwater one-dimensional DOA estimation method and device based on non-circular signal and L-type linear array |
CN109975775A (en) * | 2019-04-03 | 2019-07-05 | 西安邮电大学 | A kind of half measured data emulation mode of radar return |
CN110929371A (en) * | 2019-09-18 | 2020-03-27 | 中国电子科技集团公司第三十八研究所 | Virtual interpolation subarray calculation method and system based on least square estimation |
CN112327264A (en) * | 2020-10-14 | 2021-02-05 | 海南大学 | Bistatic FDA-MIMO radar multi-dimensional parameter joint estimation method |
CN112731303A (en) * | 2020-12-21 | 2021-04-30 | 南昌工程学院 | Interference array meter wave radar under non-Gaussian noise, and steady height measurement method and application |
CN112965059A (en) * | 2021-02-06 | 2021-06-15 | 西安邮电大学 | Real-value root-finding multi-signal classification method, system, computer equipment and application |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8238496B1 (en) * | 2010-05-18 | 2012-08-07 | Applied Micro Circuits Corporation | Multiuser multiple-input multiple-output (MU-MIMO) channel estimation for multicarrier communications |
CN103901417A (en) * | 2014-04-02 | 2014-07-02 | 哈尔滨工程大学 | Low-complexity space target two-dimensional angle estimation method of L-shaped array MIMO radar |
CN104730517A (en) * | 2015-03-23 | 2015-06-24 | 西安电子科技大学 | Bistatic MIMO radar multi-target tracking method |
CN105093185A (en) * | 2015-08-23 | 2015-11-25 | 哈尔滨工程大学 | Sparse representation-based single-base multi-output multi-input radar target direction of arrival estimation method |
CN107576940A (en) * | 2017-07-10 | 2018-01-12 | 西安邮电大学 | A kind of not rounded signal angle method of estimation of low complex degree list base MIMO radar |
-
2018
- 2018-01-29 CN CN201810081868.7A patent/CN108303683B/en not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8238496B1 (en) * | 2010-05-18 | 2012-08-07 | Applied Micro Circuits Corporation | Multiuser multiple-input multiple-output (MU-MIMO) channel estimation for multicarrier communications |
CN103901417A (en) * | 2014-04-02 | 2014-07-02 | 哈尔滨工程大学 | Low-complexity space target two-dimensional angle estimation method of L-shaped array MIMO radar |
CN104730517A (en) * | 2015-03-23 | 2015-06-24 | 西安电子科技大学 | Bistatic MIMO radar multi-target tracking method |
CN105093185A (en) * | 2015-08-23 | 2015-11-25 | 哈尔滨工程大学 | Sparse representation-based single-base multi-output multi-input radar target direction of arrival estimation method |
CN107576940A (en) * | 2017-07-10 | 2018-01-12 | 西安邮电大学 | A kind of not rounded signal angle method of estimation of low complex degree list base MIMO radar |
Non-Patent Citations (1)
Title |
---|
JIAN GONG,SHUNTIAN LOU,YIDUO GUO: "ESPRIT-like algorithm for computational-efficient angle estimation in bistatic multiple-input multiple-output radar", 《JOURNAL OF APPLIED REMOTE SENSING》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109521392B (en) * | 2018-10-24 | 2022-12-16 | 华南理工大学 | Underwater one-dimensional DOA estimation method and device based on non-circular signal and L-shaped linear array |
CN109521392A (en) * | 2018-10-24 | 2019-03-26 | 华南理工大学 | Underwater one-dimensional DOA estimation method and device based on non-circular signal and L-type linear array |
CN109521393A (en) * | 2018-11-05 | 2019-03-26 | 昆明理工大学 | A kind of DOA estimation algorithm based on signal subspace revolving property |
CN109407048B (en) * | 2018-11-12 | 2023-09-29 | 华南理工大学 | Underwater DOA estimation method and device based on non-circular signal and included angle adjustable array |
CN109407048A (en) * | 2018-11-12 | 2019-03-01 | 华南理工大学 | Underwater DOA estimation method and device based on non-circular signal and the adjustable battle array of angle |
CN109975775B (en) * | 2019-04-03 | 2022-12-23 | 西安邮电大学 | Radar echo semi-actual measurement data simulation method |
CN109975775A (en) * | 2019-04-03 | 2019-07-05 | 西安邮电大学 | A kind of half measured data emulation mode of radar return |
CN110929371A (en) * | 2019-09-18 | 2020-03-27 | 中国电子科技集团公司第三十八研究所 | Virtual interpolation subarray calculation method and system based on least square estimation |
CN112327264A (en) * | 2020-10-14 | 2021-02-05 | 海南大学 | Bistatic FDA-MIMO radar multi-dimensional parameter joint estimation method |
CN112327264B (en) * | 2020-10-14 | 2023-10-13 | 海南大学 | Bistatic FDA-MIMO radar multidimensional parameter joint estimation method |
CN112731303A (en) * | 2020-12-21 | 2021-04-30 | 南昌工程学院 | Interference array meter wave radar under non-Gaussian noise, and steady height measurement method and application |
CN112731303B (en) * | 2020-12-21 | 2024-04-09 | 南昌工程学院 | Interference array meter wave radar under non-Gaussian noise and steady height measurement method and application |
CN112965059A (en) * | 2021-02-06 | 2021-06-15 | 西安邮电大学 | Real-value root-finding multi-signal classification method, system, computer equipment and application |
Also Published As
Publication number | Publication date |
---|---|
CN108303683B (en) | 2021-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108303683A (en) | Single not rounded signal angle methods of estimation of base MIMO radar real value ESPRIT | |
CN103901417B (en) | Low-complexity space target two-dimensional angle estimation method of L-shaped array MIMO radar | |
CN102213761B (en) | Multi-target location method of bistatic common-address multi-input-multi-output radar | |
CN106707257B (en) | MIMO radar Wave arrival direction estimating method based on nested array | |
JP2999266B2 (en) | High-resolution signal evaluation method for one-dimensional or two-dimensional direction estimation or frequency estimation | |
CN107576940A (en) | A kind of not rounded signal angle method of estimation of low complex degree list base MIMO radar | |
CN109471082A (en) | Array element defect MIMO radar angle estimating method based on signal subspace reconstruct | |
CN108732549A (en) | A kind of array element defect MIMO radar DOA estimation method based on covariance matrix reconstruct | |
CN109143197B (en) | 2D-DOA and polarization parameter estimation method of polarization MIMO radar based on auxiliary array element | |
CN103983952A (en) | Low-complexity receiving and transmitting angle joint estimation method for non-circular signal double-base MIMO radar | |
CN107436429A (en) | The bistatic MIMO radar method for parameter estimation of polarization based on sparse reconstruct under impulsive noise environment | |
CN102981152A (en) | Multiple-target and send-receive angle estimation method of double-base multiple-input and multiple-output radar | |
CN104251989A (en) | Compression spatial spectrum-based single base MIMO (Multiple Input Multiple Output) radar target DOA (Direction of Arrival) estimation method | |
CN106054121A (en) | Method and device of determining radar target angle | |
CN103926555B (en) | A kind of method that utilization not rounded signal measuring antenna array receiver machine width is mutually responded | |
CN109143228A (en) | The angle estimating method of the bistatic non-circular target of MIMO radar under array mutual-coupling condition | |
CN109375154A (en) | Coherent signal method for parameter estimation based on uniform circular array under a kind of impulsive noise environment | |
CN107783078A (en) | A kind of wave beam Doppler tenth of the twelve Earthly Branches ESPRIT multiple target angle estimating method | |
CN112130111A (en) | Single-snapshot two-dimensional DOA estimation method for large-scale uniform cross array | |
CN109683151A (en) | Tenth of the twelve Earthly Branches rooting MUSIC angle estimating method under non-uniform noise environment based on matrix completion | |
CN103217671B (en) | Multi-input and multi-output fast estimation method for radar receiving and transmitting angles under color-noise environment | |
CN108828504A (en) | MIMO radar target direction method for quick estimating based on part waveform correlation | |
CN108363048A (en) | A kind of angle estimating method of the polarization MIMO radar sparse based on block | |
CN106872936A (en) | Near field sources L-type acoustic vector-sensor array row ambiguity solution Multiple Parameter Estimation Methods | |
Hassanien et al. | Subspace-based direction finding using transmit energy focusing in MIMO radar with colocated antennas |
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: 20211022 |
|
CF01 | Termination of patent right due to non-payment of annual fee |