CN109581276A - It is a kind of based on the DOA estimation method for seeking poor nested battle array - Google Patents
It is a kind of based on the DOA estimation method for seeking poor nested battle array Download PDFInfo
- Publication number
- CN109581276A CN109581276A CN201811419344.0A CN201811419344A CN109581276A CN 109581276 A CN109581276 A CN 109581276A CN 201811419344 A CN201811419344 A CN 201811419344A CN 109581276 A CN109581276 A CN 109581276A
- Authority
- CN
- China
- Prior art keywords
- battle array
- poor
- seeking
- array
- covariance matrix
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- 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
Abstract
The present invention proposes a kind of based on the DOA estimation method of poor nested battle array is sought, and asks poor nested battle array to receive signal phasor method includes the following steps: summation is calculated；Poor nested battle array is asked to receive signal phasor three second-order statistics of building according to described；It is obtained seeking poor augmentation covariance matrix according to three second-order statistics；Poor augmentation covariance matrix is asked to acquire the virtual reception signal for asking difference to optimize battle array according to described；The DOA estimated value of incoming signal is obtained according to the virtual reception signal for seeking difference optimization battle array.The present invention asks difference optimization battle array to estimate using the summation of Sparse Array, and limitation of the signal number no more than submatrix number can be estimated by breaching；It is proposed asks poor nested battle array to give full play to the advantage that array seeks difference, maximises the continuous array element for virtually seeking difference optimization battle array, the Virtual array of redundancy is less, and freedom degree is bigger, and estimation performance is more preferable, while array aperture is bigger, and resolution ratio is higher.
Description
Technical field
The invention belongs to wirelessly communicate and Radar Signal Processing Technology field, and in particular to a kind of to use array antenna low
Direction of arrival (Direction-of-Arrival, DOA) estimation method of high-precision and low complex degree is realized under cost conditions.
Background technique
With the development of space division multiple access technique and intelligent antenna technology, signal is completed using the direction of arrival (DOA) of signal
Airspace capture and tracking attracted the research of a large amount of domestic and foreign scholars, especially in radar, sonar, navigation, communication, radio day
The numerous areas such as literature.
Many traditional DOA estimation methods are typically all to be based on even linear array (Uniform Linear Array, ULA),
But the estimable signal number of these methods will receive the limitation of array number, i.e., estimable signal number is no more than array number.
Moreover, limitation of the ULA due to array element spacing no more than half-wavelength, if wanting to increase array aperture, promoting DOA estimated accuracy and divide
Resolution must just increase array element number, therefore, will cause excessively complicated and system cost the increase of system.In view of existing for ULA
The above problem, in recent years, people have also been proposed Sparse Array.Wherein, nested battle array and relatively prime battle array are the main Sparse Arrays of two classes, and by
The extensive research of domestic and foreign scholars is arrived.As P.P.Vaidyanathan proposes a kind of space based on nested battle array with P.Pal
Smooth DOA estimation method.This method seeks the freedom degree and performance of poor Concept Promoting DOA estimation based on array.Meanwhile Piya
Pal and P.P.Vaidyanathan also proposed a kind of spatial-smoothing DOA estimation method based on relatively prime battle array, and this method is based on battle array
Column ask the concept of difference that can equally promote the freedom degree that DOA estimates, and mutual coupling is less than nested battle array.But relatively prime battle array
Ask difference optimization battle array there are hole, the freedom degree of virtual array could not be fully used.The Sparse Arrays such as nested battle array and relatively prime battle array with
Traditional ULA is compared, and in the identical situation of array element number, possesses bigger array aperture and smaller array element mutual coupling, is improved
DOA estimated accuracy, resolution ratio and maximum accessible signal number.On the other hand, array aperture under the same conditions, it is sparse
Array number needed for battle array is less, it means that more small-scale reception system and signal processing system etc. significantly reduce and be
System cost.
Although nested battle array and relatively prime battle array have used seeks difference optimization battle array, improves the freedom degree of DOA estimation to a certain extent,
But summation optimization battle array could not be used.
Summary of the invention
In view of the foregoing deficiencies of prior art, it is based on seeking poor nested battle array the purpose of the present invention is to provide one kind
DOA estimation method, with solve in the prior art DOA estimation freedom degree it is lower, estimation performance it is poor；Lower estimated accuracy
And resolution ratio；The disadvantages of virtually asking the redundancy array element of difference optimization battle array more.
In order to achieve the above objects and other related objects, the present invention provides a kind of based on asking the DOA of poor nested battle array to estimate
Meter method, method includes the following steps:
Summation, which is calculated, asks poor nested battle array to receive signal phasor x (t)；
Poor nested battle array is asked to receive signal phasor x (t) building three second-order statistics R1, R2, R3 according to described:
It is obtained seeking poor augmentation covariance matrix according to three second-order statistics
Poor augmentation covariance matrix is sought according to describedAcquire the virtual reception signal for seeking difference optimization battle array
According to the virtual reception signal for seeking difference optimization battle arrayObtain the DOA estimated value of incoming signal.
Optionally, described to be obtained seeking poor augmentation covariance matrix according to three second-order statisticsInclude:
T snap is carried out to second-order statistic R1, R2, R3 respectively, obtains the estimated value of second-order statistic R1, R2, R3Then seek poor augmentation covariance matrixAre as follows:
Optionally, poor augmentation covariance matrix is sought according to describedAcquire the virtual reception letter for seeking difference optimization battle array
NumberIt specifically includes:
Poor augmentation covariance matrix is sought in vectorizationObtain measurement vector
To measurement vectorCarry out the virtual reception signal that de-redundancy obtains asking difference optimization battle array
Optionally, according to the virtual reception signal for seeking difference optimization battle arrayThe DOA estimated value of incoming signal is obtained,
It specifically includes:
The smoothed covariance matrix R of difference optimization battle array reception signal is sought in buildingsmooth；
To smoothed covariance matrix RsmoothIt carries out feature decomposition and obtains signal subspace UsWith noise subspace Un；
Using MUSIC algorithm, spectral function is constructed；
The DOA estimated value of incoming signal is obtained by search spectrum peak.
In order to achieve the above objects and other related objects, the present invention also provides a kind of based on the DOA for seeking poor nested battle array
Estimating system, the estimating system include:
First computing module asks poor nested battle array to receive signal phasor x (t) for summation to be calculated；
Module is constructed, for asking poor nested battle array to receive three second-order statistics of signal phasor x (t) building according to described
R1, R2, R3:
Second computing module seeks poor augmentation covariance matrix for obtaining according to three second-order statistics
Third computing module, for seeking poor augmentation covariance matrix according to describedIt acquires and seeks difference optimization battle array
It is virtual to receive signal
Estimation module, for according to the virtual reception signal for seeking difference optimization battle arrayThe DOA for obtaining incoming signal estimates
Evaluation.
Optionally, described to be obtained seeking poor augmentation covariance matrix according to three second-order statisticsInclude:
T snap is carried out to second-order statistic R1, R2, R3 respectively, obtains the estimated value of second-order statistic R1, R2, R3Then seek poor augmentation covariance matrixAre as follows:
Optionally, poor augmentation covariance matrix is sought according to describedAcquire the virtual reception letter for seeking difference optimization battle array
NumberIt specifically includes:
Poor augmentation covariance matrix is sought in vectorizationObtain measurement vector
To measurement vectorCarry out the virtual reception signal that de-redundancy obtains asking difference optimization battle array
Optionally, according to the virtual reception signal for seeking difference optimization battle arrayThe DOA estimated value of incoming signal is obtained,
It specifically includes:
The smoothed covariance matrix R of difference optimization battle array reception signal is sought in buildingsmooth；
To smoothed covariance matrix RsmoothIt carries out feature decomposition and obtains signal subspace UsWith noise subspace Un；
Using MUSIC algorithm, spectral function is constructed；
The DOA estimated value of incoming signal is obtained by search spectrum peak.
As described above, of the invention is a kind of based on the DOA estimation method and system of seeking poor nested battle array, have with following
Beneficial effect:
The present invention asks difference optimization battle array to estimate using the summation of Sparse Array, and signal number can be estimated no more than by breaching
The limitation of submatrix number；It is proposed asks poor nested battle array to give full play to the advantage that array seeks difference, maximises and virtually asks
With the continuous array element for seeking difference optimization battle array, the Virtual array of redundancy is less, and freedom degree is bigger, and estimation performance is more preferable, while array hole
Diameter is bigger, and resolution ratio is higher.
Detailed description of the invention
In order to which the present invention is further explained, described content, with reference to the accompanying drawing makees a specific embodiment of the invention
Further details of explanation.It should be appreciated that these attached drawings are only used as typical case, and it is not to be taken as to the scope of the present invention
It limits.
Fig. 1 is that schematic diagram is arranged in array of the present invention；
Fig. 2 is the mentioned array of the present invention and the rooting mean square error of DOA estimation with SNR variation relation schematic diagram；
Fig. 3 is the mentioned array of the present invention and the rooting mean square error of DOA estimation with number of snapshots variation relation schematic diagram；
Fig. 4 is a kind of flow chart based on the DOA estimation method for seeking poor nested battle array described in the embodiment of the present invention.
Specific embodiment
Illustrate embodiments of the present invention below by way of specific specific example, those skilled in the art can be by this specification
Other advantages and efficacy of the present invention can be easily understood for disclosed content.The present invention can also pass through in addition different specific realities
The mode of applying is embodied or practiced, the various details in this specification can also based on different viewpoints and application, without departing from
Various modifications or alterations are carried out under spirit of the invention.It should be noted that in the absence of conflict, following embodiment and implementation
Feature in example can be combined with each other.
It should be noted that illustrating the basic structure that only the invention is illustrated in a schematic way provided in following embodiment
Think, only shown in schema then with related component in the present invention rather than component count, shape and size when according to actual implementation
Draw, when actual implementation kenel, quantity and the ratio of each component can arbitrarily change for one kind, and its assembly layout kenel
It is likely more complexity.
As shown in figure 4, it is a kind of based on the DOA estimation method for seeking poor nested battle array, method includes the following steps:
Step 1: summation is calculated, poor nested battle array is asked to receive signal phasor x (t)；
If poor nested battle array is asked to be made of two submatrixs, each submatrix is even linear array (ULA), and between two submatrixs
Every L, L is unknown parameter to be asked.Submatrix 1 is an intensive ULA, array number N1, array element spacing is unit spacing d, d=
λ/2, λ are signal wavelength.Then the element position of submatrix 1 is represented by set:
Submatrix 2 is a sparse ULA, array number N2, array element spacing is D=(2N1-1)d.D is also to be asked unknown
Parameter.Then the element position of submatrix 2 is represented by set:
(1a) determines submatrix interval L:
Submatrix 1 is an intensive ULA, and asking difference optimization battle array and summation optimization battle array is also intensive ULA.WithIndicate son
The summation of battle array 1 optimizes battle array, then has
Wherein,For array summation operator.Also referred to as submatrix 1 from and.Similarly, available submatrix 2 from and are as follows:
Use l2,iIndicate the position of i-th of the array element of submatrix 2,1≤i≤N2, then l2,i=(N1-1)d+L+(i-1)D.For benefit
With submatrix 1 from and filling submatrix interval caused by hole, then the position of the first array element of submatrix 2 should with submatrix 1 from and most
The latter element is adjacent, then has following equation establishment:
l2,1=2 (N1-1)d+d
Wherein, l2,1=(N1- 1) d+L, substitution solve: L=N1D, submatrix interval as to be asked.
(1b) determines the array element space D of submatrix 2:
The submatrix interval L=N obtained by (1a)1D indicates the element position of submatrix 2 again are as follows:
Certainly and similar with submatrix 1, the summation optimization battle array of submatrix 1 and submatrix 2 may be expressed as:
Submatrix 1 and submatrix 2 ask difference optimization battle array may be expressed as:
Wherein,Poor operator is sought for array.Also referred to as submatrix 1 and submatrix 2 mutual and,Also referred to as submatrix 1 and submatrix 2
Mutual deviation.
It submatrix 1 and submatrix 2 mutual and submatrix 1 can be decomposed into sum respectively with the array element in submatrix 2, then submatrix 1 and son
I-th of array element of battle array 2 sums to obtain virtual location set:
Similarly, the mutual deviation of submatrix 1 and submatrix 2 can be decomposed into submatrix 1 and ask poor with the array element in submatrix 2 respectively, then submatrix 1
Difference is asked to obtain location sets with i-th of array element of submatrix 2:
Then (i-1)-th array element of submatrix 1 and submatrix 2 sums to obtain location sets:
In order to using mutually andAnd mutual deviationHole present in Sparse Array is filled, thenThe last one element should be withHeader element it is adjacent, therefore have following equation establishment:
(3N1- 2) d+ (i-2) D+d=N1d+(i-1)D
It solves: D=(2N1- 1) d, the array element spacing of submatrix 2 as to be asked.
(1c) establishes array received model:
Utilize the submatrix interval L=N being previously obtained1Array element space D=(2N of d and submatrix 21- 1) d can be determined completely and be mentioned
The array structure for seeking poor nested battle array out.
As shown in Figure 1, poor nested battle array, array element sum N=N are asked in setting1+N2.Seek the element position of poor nested battle array
It is represented byVector d=[d also can be used1,d2,...,dN]TIndicate the element position of array, wherein
I=1,2 ..., N.
Assuming that there is K irrelevant far field narrowband real signal sk(t) from direction θkIt is incident on array, k=1,2 ..., K, θkTable
Show the incident angle of k-th signal.Noise is independent identically distributed additive white Gaussian noise, and independent with signal.Then ask
Poor nesting battle array, which receives signal phasor, may be expressed as:
Wherein, A=[a (θ1),a(θ2),…,a(θK)] it is the array manifold matrix for seeking poor nested battle array,Indicate the corresponding steering vector of k-th of signal.S (t)=[s1(t),s2
(t),...,sK(t)]TIndicate signal phasor,For noise vector, the multiple height of element independent same distribution and obedience
This distribution For noise power.
Step 2: poor augmentation covariance matrix is sought in building:
Three second-order statistics are constructed as follows using array received signal x (t):
R1=E [x (t) xH(t)]
R2=E [x (t) xT(t)]
R3=E [x*(t)xH(t)]
Wherein, xH(t) conjugate transposition for asking poor nested battle array to receive signal phasor x (t), x are indicatedT(t) it indicates to ask
Poor nesting battle array receives the transposition of signal phasor x (t).x*(t) conjugation for asking poor nested battle array to receive signal phasor x (t) is indicated
Matrix.
But R1It is unavailable ideal covariance matrix, in fact, estimating to obtain by T snap:
Wherein, T is number of snapshots.Similarly, can estimate to obtain:
It is further constructed followed by above-mentioned second-order statistic and seeks poor augmentation covariance matrix:
Step 3: acquire the virtual reception signal for seeking difference optimization battle array:
Summation in vectorization step 2 asks poor augmentation covariance matrix to obtain measurement vector:
Wherein,K=1,
2 ..., K,INFor N rank unit matrix, vec () is vectorization operator,For signal power.ThenIt can regard as
Poor nested battle array is asked virtually to optimize array manifold matrix corresponding to battle array, p, which can regard as, is incident on the single fast of the virtual optimization battle array
Clap signal phasor.But measurement vectorIt is middle there are the element of redundancy, i.e. the Virtual array of same position has corresponded to multiple observations
Data actually only need one.Therefore, it is necessary to rightIt goes to carry out de-redundancy, i.e. the Virtual array of same position is corresponding
Multiple observation data only take one or are averaged, and obtain:
Wherein,It is the irredundant measurement vector of array, γ=2N1N2+N1-N2.Then irredundant measurement vector
The virtual reception signal for seeking difference optimization battle array can be regarded as, which is the virtual ULA that an array number is 2 γ -1.
Step 4: using the DOA of smooth MUSIC algorithm estimation signal
The smoothed covariance matrix of difference optimization battle array reception signal is sought in (4a) building:
Difference optimization battle array is asked to be divided into γ overlapped and size identical submatrixs summation, each submatrix has γ battle array
Member.Then the virtual reception signal of i-th of submatrix is correspondingIn -1 element of the i-th to the i-th+γ, be denoted as zi, i=1,2 ..., γ.
The smoothed covariance matrix of difference optimization battle array reception signal is sought in following building:
Wherein,Indicate ziConjugate transposition.
(4b) building spectral function simultaneously obtains the estimated value of final DOA
To matrix RsmoothIt carries out feature decomposition and obtains signal subspace and noise subspace:
Wherein, ΛsIt is the diagonal matrix of K × K dimension, diagonal element RsmoothLarger K characteristic value；UsIt is (2 γ -1) × K
Dimensional signal subspace, by RsmoothThe corresponding feature vector of the biggish characteristic value of K at,Indicate UsConjugate transposition；
ΛnIt is the diagonal matrix of (2 γ -1-K) × (2 γ -1-K) dimension, diagonal element RsmoothLesser 2 γ -1-K characteristic values；Un
It is (2 γ -1) × (2 γ -1-K) dimension noise subspace, by RsmoothThe corresponding feature vector of lesser 2 γ -1-K characteristic values
At,Indicate UnConjugate transposition.
Then, using MUSIC algorithm, spectral function is constructed as follows
Wherein, avir(θ) is the corresponding steering vector of virtual ULA that array number is γ.It can then be obtained by search spectrum peak
The DOA estimated value of incoming signal.
Poor nested battle array and document [1] are asked in order to analyze summation proposed by the invention: " P.Pal,
P.P.Vaidyanathan.Nested arrays:a novel approach to array processing with
enhanced degrees of freedom.IEEE Transactions on Signal Processing,2010,58
(8): based on the nesting for seeking difference optimization battle array in 4167-4181 " (a kind of spatial-smoothing DOA estimation method based on nested battle array)
Battle array, document [2]: " P.Pal, P.P.Vaidyanathan.Coprime sampling and the MUSIC
algorithm.Proc.of IEEE Digital Signal Processing Workshop and IEEE Signal
A kind of Processing Education Workshop, Sedona, 2011,289-294 " (space smoothing based on relatively prime battle array
DOA estimation method) in based on ask difference optimization battle array relatively prime battle array and document [3]: " X.Wang, X.Wang, X.Lin.Co-prime
array processing with sum and difference co-array.49th Asilomar Conf.on
Signals, Systems, and Computers, 2015,380-384 " it (is proposed based on nested battle array and relatively prime battle array based on summation
Seek the DOA estimation method of poor concept) in based on ask difference optimization battle array nested battle array, the estimation performance of relatively prime battle array, devise two
Emulation experiment is organized to be compared.Therein it is proposed that summation to seek poor nested battle array array parameter be N1=N2=6, document [1] and text
The array parameter for the nested battle array that [3] use is offered as N1=N2=6, the array parameter for the relatively prime battle array that document [2] and document [3] use
Array element sum for M=4, N=5, all arrays is 12.Signal number is 16, and incident direction is equal in the range of -60 ° to 60 °
Even distribution.DOA estimation method is all made of space smoothing MUSIC algorithm.
Battery of tests number of snapshots are 3000, and carry out 1000 independent experiments, the rooting mean square error of DOA estimation
(RMSE) relationship changed with signal-to-noise ratio (SNR) is as shown in Figure 2.
Second group of test SNR is 0dB, equally carries out 1000 independent experiments, what the RMSE of DOA estimation changed with number of snapshots
Relationship is as shown in Figure 3.
It can be seen from the figure that the freedom degree for asking poor nested battle array that can greatly improve DOA estimation that the present invention is mentioned,
And significantly improve the estimation performance of DOA.
The present invention also provides a kind of based on the DOA estimating system for seeking poor nested battle array, which includes:
First computing module asks poor nested battle array to receive signal phasor x (t) for summation to be calculated；
Module is constructed, for asking poor nested battle array to receive three second-order statistics of signal phasor x (t) building according to described
R1, R2, R3:
Second computing module seeks poor augmentation covariance matrix for obtaining according to three second-order statistics
Third computing module, for seeking poor augmentation covariance matrix according to describedIt acquires and seeks difference optimization battle array
It is virtual to receive signal
Estimation module, for according to the virtual reception signal for seeking difference optimization battle arrayThe DOA for obtaining incoming signal estimates
Evaluation.
Preferably, described to be obtained seeking poor augmentation covariance matrix according to three second-order statisticsInclude:
T snap is carried out to second-order statistic R1, R2, R3 respectively, obtains the estimated value of second-order statistic R1, R2, R3Then seek poor augmentation covariance matrixAre as follows:
Preferably, poor augmentation covariance matrix is sought according to describedAcquire the virtual reception letter for seeking difference optimization battle array
NumberIt specifically includes:
Poor augmentation covariance matrix is sought in vectorizationObtain measurement vector
To measurement vectorCarry out the virtual reception signal that de-redundancy obtains asking difference optimization battle array
Preferably, according to the virtual reception signal for seeking difference optimization battle arrayThe DOA estimated value of incoming signal is obtained,
It specifically includes:
The smoothed covariance matrix R of difference optimization battle array reception signal is sought in buildingsmooth；
To smoothed covariance matrix RsmoothIt carries out feature decomposition and obtains signal subspace UsWith noise subspace Un；
Using MUSIC algorithm, spectral function is constructed；
The DOA estimated value of incoming signal is obtained by search spectrum peak.
Module/unit realization described in the present embodiment can be realized by method above-mentioned.
The above-described embodiments merely illustrate the principles and effects of the present invention, and is not intended to limit the present invention.It is any ripe
The personage for knowing this technology all without departing from the spirit and scope of the present invention, carries out modifications and changes to above-described embodiment.Cause
This, institute is complete without departing from the spirit and technical ideas disclosed in the present invention by those of ordinary skill in the art such as
At all equivalent modifications or change, should be covered by the claims of the present invention.
Claims (8)
1. a kind of based on the DOA estimation method for seeking poor nested battle array, which is characterized in that method includes the following steps:
Summation, which is calculated, asks poor nested battle array to receive signal phasor x (t)；
Poor nested battle array is asked to receive signal phasor x (t) building three second-order statistics R1, R2, R3 according to described；
It is obtained seeking poor augmentation covariance matrix according to three second-order statistics
Poor augmentation covariance matrix is sought according to describedAcquire the virtual reception signal for seeking difference optimization battle array
According to the virtual reception signal for seeking difference optimization battle arrayObtain the DOA estimated value of incoming signal.
2. according to claim 1 a kind of based on the DOA estimation method for seeking poor nested battle array, which is characterized in that described
It is obtained seeking poor augmentation covariance matrix according to three second-order statisticsInclude:
T snap is carried out to second-order statistic R1, R2, R3 respectively, obtains the estimated value of second-order statistic R1, R2, R3Then seek poor augmentation covariance matrixAre as follows:
3. according to claim 2 a kind of based on the DOA estimation method for seeking poor nested battle array, which is characterized in that according to
It is described to seek poor augmentation covariance matrixAcquire the virtual reception signal for seeking difference optimization battle arrayIt specifically includes:
Poor augmentation covariance matrix is sought in vectorizationObtain measurement vector
To measurement vectorCarry out the virtual reception signal that de-redundancy obtains asking difference optimization battle array
4. according to claim 3 a kind of based on the DOA estimation method for seeking poor nested battle array, which is characterized in that according to
The virtual reception signal for seeking difference optimization battle arrayThe DOA estimated value of incoming signal is obtained, is specifically included:
The smoothed covariance matrix R of difference optimization battle array reception signal is sought in buildingsmooth；
To smoothed covariance matrix RsmoothIt carries out feature decomposition and obtains signal subspace UsWith noise subspace Un；
Using MUSIC algorithm, spectral function is constructed；
The DOA estimated value of incoming signal is obtained by search spectrum peak.
5. a kind of based on the DOA estimating system for seeking poor nested battle array, which is characterized in that the estimating system includes:
First computing module asks poor nested battle array to receive signal phasor x (t) for summation to be calculated；
Construct module, for according to it is described ask poor nested battle array receive three second-order statistic R1 of signal phasor x (t) building,
R2, R3:
Second computing module seeks poor augmentation covariance matrix for obtaining according to three second-order statistics
Third computing module, for seeking poor augmentation covariance matrix according to describedIt acquires and asks the virtual of difference optimization battle array
Receive signal
Estimation module, for according to the virtual reception signal for seeking difference optimization battle arrayObtain the DOA estimation of incoming signal
Value.
6. according to claim 5 a kind of based on the DOA estimating system for seeking poor nested battle array, which is characterized in that described
It is obtained seeking poor augmentation covariance matrix according to three second-order statisticsInclude:
T snap is carried out to second-order statistic R1, R2, R3 respectively, obtains the estimated value of second-order statistic R1, R2, R3Then seek poor augmentation covariance matrixAre as follows:
7. according to claim 6 a kind of based on the DOA estimating system for seeking poor nested battle array, which is characterized in that according to
It is described to seek poor augmentation covariance matrixAcquire the virtual reception signal for seeking difference optimization battle arrayIt specifically includes:
Poor augmentation covariance matrix is sought in vectorizationObtain measurement vector
To measurement vectorCarry out the virtual reception signal that de-redundancy obtains asking difference optimization battle array
8. according to claim 7 a kind of based on the DOA estimating system for seeking poor nested battle array, which is characterized in that according to
The virtual reception signal for seeking difference optimization battle arrayThe DOA estimated value of incoming signal is obtained, is specifically included:
The smoothed covariance matrix R of difference optimization battle array reception signal is sought in buildingsmooth；
To smoothed covariance matrix RsmoothIt carries out feature decomposition and obtains signal subspace UsWith noise subspace Un；
Using MUSIC algorithm, spectral function is constructed；
The DOA estimated value of incoming signal is obtained by search spectrum peak.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811419344.0A CN109581276A (en) | 2018-11-26 | 2018-11-26 | It is a kind of based on the DOA estimation method for seeking poor nested battle array |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811419344.0A CN109581276A (en) | 2018-11-26 | 2018-11-26 | It is a kind of based on the DOA estimation method for seeking poor nested battle array |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109581276A true CN109581276A (en) | 2019-04-05 |
Family
ID=65924140
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811419344.0A Pending CN109581276A (en) | 2018-11-26 | 2018-11-26 | It is a kind of based on the DOA estimation method for seeking poor nested battle array |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109581276A (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110031794A (en) * | 2019-04-16 | 2019-07-19 | 中国人民解放军国防科技大学 | Coherent information source DOA estimation method based on difference common matrix reconstruction |
| CN110308417A (en) * | 2019-05-30 | 2019-10-08 | 电子科技大学 | Wave arrival direction estimating method and device under nested battle array element failure based on matrix fill-in |
| CN111474534A (en) * | 2020-04-16 | 2020-07-31 | 电子科技大学 | Two-dimensional DOA estimation method based on symmetric parallel nested array |
| CN111650552A (en) * | 2020-05-15 | 2020-09-11 | 北京理工大学 | Improved L-shaped array design method based on sum-difference cooperative array construction |
| CN113296049A (en) * | 2021-04-28 | 2021-08-24 | 南京航空航天大学 | Conjugate amplification DOA estimation method of non-circular signals in co-prime array pulse environment |
| CN113391261A (en) * | 2021-07-06 | 2021-09-14 | 南京航空航天大学 | Extended co-prime array rapid DOA estimation method based on Taylor expansion |
| CN113589224A (en) * | 2021-08-03 | 2021-11-02 | 宜宾电子科技大学研究院 | DOA estimation method based on enhanced nested array |
| CN113589223A (en) * | 2021-06-11 | 2021-11-02 | 南京邮电大学 | Direction finding method based on nested array under mutual coupling condition |
-
2018
- 2018-11-26 CN CN201811419344.0A patent/CN109581276A/en active Pending
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110031794A (en) * | 2019-04-16 | 2019-07-19 | 中国人民解放军国防科技大学 | Coherent information source DOA estimation method based on difference common matrix reconstruction |
| CN110308417A (en) * | 2019-05-30 | 2019-10-08 | 电子科技大学 | Wave arrival direction estimating method and device under nested battle array element failure based on matrix fill-in |
| CN110308417B (en) * | 2019-05-30 | 2021-07-02 | 电子科技大学 | Method and device for estimating direction of arrival under nested array element failure based on matrix filling |
| CN111474534A (en) * | 2020-04-16 | 2020-07-31 | 电子科技大学 | Two-dimensional DOA estimation method based on symmetric parallel nested array |
| CN111474534B (en) * | 2020-04-16 | 2023-04-07 | 电子科技大学 | Two-dimensional DOA estimation method based on symmetric parallel nested array |
| CN111650552A (en) * | 2020-05-15 | 2020-09-11 | 北京理工大学 | Improved L-shaped array design method based on sum-difference cooperative array construction |
| CN113296049A (en) * | 2021-04-28 | 2021-08-24 | 南京航空航天大学 | Conjugate amplification DOA estimation method of non-circular signals in co-prime array pulse environment |
| CN113296049B (en) * | 2021-04-28 | 2024-02-20 | 南京航空航天大学 | Conjugate augmented DOA estimation method for non-circular signals in intersubstance array pulse environment |
| CN113589223A (en) * | 2021-06-11 | 2021-11-02 | 南京邮电大学 | Direction finding method based on nested array under mutual coupling condition |
| CN113589223B (en) * | 2021-06-11 | 2023-05-05 | 南京邮电大学 | Direction finding method based on nested array under mutual coupling condition |
| CN113391261A (en) * | 2021-07-06 | 2021-09-14 | 南京航空航天大学 | Extended co-prime array rapid DOA estimation method based on Taylor expansion |
| CN113589224A (en) * | 2021-08-03 | 2021-11-02 | 宜宾电子科技大学研究院 | DOA estimation method based on enhanced nested array |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109581276A (en) | It is a kind of based on the DOA estimation method for seeking poor nested battle array | |
| CN109932680B (en) | Non-circular signal direction of arrival estimation method based on translational co-prime array | |
| CN108872929B (en) | Estimation method for direction of arrival of co-prime array based on rotation invariance of covariance matrix subspace of interpolated virtual array | |
| CN104749553B (en) | Direction of arrival angle method of estimation based on rapid sparse Bayesian learning | |
| CN104698433B (en) | DOA Estimation in Coherent Signal method based on single snapshot data | |
| US6311043B1 (en) | Method and measurement configuration for measuring the characteristics of radio channels | |
| CN109490820A (en) | A kind of arrival direction estimation method based on parallel nested battle array | |
| CN109655799A (en) | The non-homogeneous thinned array direction-finding method of covariance matrix vectorization based on IAA | |
| CN107290709B (en) | The relatively prime array Wave arrival direction estimating method decomposed based on vandermonde | |
| CN108710102B (en) | Direction-of-arrival estimation method based on second-order equivalent virtual signal inverse discrete Fourier transform of co-prime array | |
| CN108896954A (en) | A kind of direction of arrival estimation method based on joint real value subspace in relatively prime battle array | |
| CN102841344A (en) | Method for estimating parameters of near-field broadband signal resources by utilizing less array elements | |
| CN109683126A (en) | Direction of arrival measurement method, signal handling equipment and storage medium | |
| CN109471086A (en) | Relatively prime MIMO radar Wave arrival direction estimating method based on more sampling snap sum aggregate array signal discrete Fourier transforms | |
| CN111856402B (en) | Signal processing method and device, storage medium and electronic device | |
| CN110045323A (en) | A kind of relatively prime battle array robust adaptive beamforming algorithm based on matrix fill-in | |
| CN110286350A (en) | A kind of perfect match method and device of L-type Sparse Array DOA estimation | |
| CN112904272B (en) | Three-dimensional co-prime cubic array direction of arrival estimation method based on cross-correlation tensor | |
| CN106483493A (en) | A kind of sparse double parallel linear array and estimating two-dimensional direction-of-arrival method | |
| CN107907855A (en) | A kind of coprime array switching DOA estimation method and device for even linear array | |
| CN110018439A (en) | The direction-finding method that a kind of and difference beam is combined with wave beam MUSIC | |
| CN110286351A (en) | A kind of arrival direction estimation method and device based on L-type nesting battle array | |
| CN107302391A (en) | Adaptive beamforming method based on relatively prime array | |
| CN109782218A (en) | A kind of non-circular signal DOA estimation method of relevant distribution based on double parallel antenna array | |
| CN106680779B (en) | Beam-forming method and device under impulsive noise |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190405 |