CN109188346A - Macroscale homogenous cylindrical array list snap DOA estimation method - Google Patents
Macroscale homogenous cylindrical array list snap DOA estimation method Download PDFInfo
- Publication number
- CN109188346A CN109188346A CN201811006677.0A CN201811006677A CN109188346A CN 109188346 A CN109188346 A CN 109188346A CN 201811006677 A CN201811006677 A CN 201811006677A CN 109188346 A CN109188346 A CN 109188346A
- Authority
- CN
- China
- Prior art keywords
- array
- vector
- uniform
- signal
- even linear
- 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
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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a kind of Macroscale homogenous cylindrical array list snap DOA estimation method, mainly solve the problems, such as that prior art operand is high, hardware realization difficulty is big.Its implementation is: corresponding parameter is arranged according to the model of uniform cylindrical array, obtains guiding matrix A；The single snap signal of array received is lined up into matrix form X；FFT transform, no-coherence cumulating and maximum value detection are taken turns doing to each row of X, obtain pitch angle rough estimate evaluationIntroduce pitching phase compensation variable ηel, rightSection of closing on scan for related to matching, obtain the fine estimated value of pitch angleMatched transform is carried out to X and obtains imaginary circles battle array signal yUCA, then cyclic convolution and maximum value detection are carried out, obtain azimuth rough estimate evaluationIntroduce orientation phase compensation variable ηaz, rightSection of closing on scan for related to matching, obtain the fine estimated value in azimuthThe present invention can greatly reduce operand, can be used for target detection, target following or signal extraction under the premise of guaranteeing estimated accuracy.
Description
Technical field
The invention belongs to signal processing technology fields, in particular to direction of arrival DOA estimation method, can be used for target inspection
Survey, target following or signal extraction.
Background technique
Array signal processing is the different location composition sensor array that multiple sensors are arranged in space, and utilizes this
An array is received and is handled to spacing wave field, to extract array received signal and its characteristic information, is inhibited simultaneously
Interference and noise or uninterested information.
Direction of arrival DOA estimation is one of main research in array signal processing and radar, communication and sonar
It is research emphasis all the time etc. one of the vital task in many fields.Traditional DOA algorithm for estimating include MUSIC,
ESPRIT and related derivative algorithm, but these algorithms are all limited by operand, not can be used directly in large scale array,
Therefore the DOA algorithm for estimating based on conventional beamformer embodies advantage, while the low disadvantage of its angular resolution is advised due to array
The big factor of mould has obtained certain make up.
Most of DOA algorithm for estimating carries out research discussion both for even linear array at present, this is because uniform line
Battle array meets vandermonde matrix form, consequently facilitating processing mathematically.However the characteristics of even linear array is can only to provide 180 °
Nothing in range obscures azimuth information, and array normal directional resolution is high, and axis direction is poor.For disadvantage mentioned above, need
Study special construction array, including rectangular surfaces battle array, uniform circular array and conformal array etc..
It is main still to be formed based on phase pattern transform method currently for the DOA estimation method of uniform cylindrical array
Each uniform circular array of uniform cylindrical array is mapped as virtual line arrays, then applies UCA-RB-MUSIC, UCA-Root-MUSIC
Or the subsequent related algorithm such as UCA-ESPRIT.In the case of Macroscale homogenous cylindrical-array, operational precision needs phase mode
The guarantee of formula number, but phase pattern number increases the great growth that will lead to operand；Subsequent related algorithm needs to carry out
Eigenvalue Decomposition and polynomial rooting etc., operand is very huge；In addition, how soon subsequent related algorithm claps product if generally requiring
It is tired, also further cause the increase of operand.
Summary of the invention
It is an object of the invention to be directed to above-mentioned problem of the prior art, it is fast to provide a kind of Macroscale homogenous cylindrical array list
DOA estimation method is clapped, greatly to reduce operand under the premise of guaranteeing estimated accuracy.
In order to achieve the above objectives, technical solution of the present invention includes the following:
(1) corresponding parameter is arranged according to the model of uniform cylindrical array, obtains constituting the nicely rounded of uniform cylindrical array
Battle array steering vector aUCAWith the steering vector a of even linear arrayULA, and according to uniform circular array steering vector aUCAWith leading for even linear array
To vector aULAObtain the guiding matrix A of uniform cylindrical array；
(2) according to the guiding matrix A of uniform cylindrical array, the single snap signal of the array received is lined up into matrix formEach element respectively corresponds the array element accordingly numbered, wherein xmnIndicate the of the n-th circle uniform circular array
The reception signal of m array element, M indicate that the array number that each uniform circular array contains, N indicate the array element that each even linear array contains
Number；
(3) according to signal snap matrix X, estimate pitch angle of the signal source relative to uniform cylindrical array
Uniform cylindrical array is split into even linear array by (3a), i.e., every a line of X is individually carried out FFT transform, then will be each
No-coherence cumulating is done between row FFT transform result, obtains no-coherence cumulating result vector sel；
(3b) is to no-coherence cumulating result vector selMaximum value detection is done, and calculates signal source relative to uniform cylindrical array
Pitch angleRough estimate evaluation
(3c) introduces pitching phase compensation variable ηel, to pitch angle rough estimate evaluationClose on section scan for and
With correlation, and calculate pitch angleFine estimated value
(4) according to the fine estimated value of signal snap matrix X and pitch angleEstimate signal source relative to uniform cylindrical battle array
The azimuth angle theta of column:
A) according to the fine estimated value of pitch angleUniform cylindrical array is synthesized into a uniform circular array, i.e., by X process
Imaginary circles battle array signal y is got in distribution transforming in returnUCA；
(4b) is to imaginary circles battle array signal yUCACyclic convolution and maximum value detection are carried out, and calculates signal source relative to uniform
The rough estimate evaluation of the azimuth angle theta of cylindrical-array
(4c) introduces orientation phase compensation variable ηaz, azimuthal rough estimate evaluationSection of closing on scan for and match
Correlation, and the fine estimated value of computer azimuth angle θ
The invention has the benefit that
(1) compared to existing algorithm, under the premise of guaranteeing estimated accuracy, greatly reduce operand.
Existing to be directed to the signal processing algorithm comprising uniform circular array structure, being mostly based on phase pattern transform method will be equal
Nicely rounded battle array is mapped as virtual line arrays, then using subsequent phases such as UCA-RB-MUSIC, UCA-Root-MUSIC or UCA-ESPRIT
Close algorithm.In this process, the complexity of Eigenvalues Decomposition scheduling algorithm is Ο (M3), wherein M is the phase pattern of mapping
Number, and in order to guarantee that operational precision, phase pattern number are often suitable with the element number of array of uniform circular array, therefore for extensive
For array, the algorithm based on Eigenvalues Decomposition is almost not achievable.
The present invention is that the operation of rough estimate is reduced by FFT and cyclic convolution on the basis of conventional beamformer algorithm
Amount, be re-introduced into phase compensation and close on section search for it is related to matching, realize the raising of estimated accuracy.The present invention proposes algorithm
Computational complexity is Ο (Mlog2M+GM) magnitude, wherein M indicates the element number of array of even linear array or uniform circular array, and G expression is searched
Rope interval division number, G < < M usual for more massive array, this hour operation quantity can simplify as Ο (Mlog2M it) measures
Grade is much smaller than existing algorithm.
(2) algorithm structure is simple, is easily programmed or realizes on hardware configuration.
In existing processing method, include UCA-RB-MUSIC, UCA-Root-MUSIC or UCA-ESPRIT algorithm, need into
The operations such as row Eigenvalues Decomposition and polynomial solving, in common signal processor generally do not include such operation wrap or
Library.And the operation that the present invention needs is mainly FFT and searching algorithm, can in common signal processor operation packet or
Library fast programming is realized.
(3) by setting pitching or the division number in bearing search section, in arithmetic speed and can estimate according to the actual situation
Weigh between meter precision.
The present invention needs to set pitching or bearing search section when searching for optimal pitching or orientation phase compensation variable
Number is divided, division number is more, then estimated accuracy is higher, but arithmetic speed is slower, otherwise division number is fewer, then operation
Speed is faster, but estimated accuracy is lower, and therefore, user can be according to the actual situation between arithmetic speed and estimated accuracy
Tradeoff.
Detailed description of the invention
Fig. 1 is implementation flow chart of the invention；
Fig. 2 is the uniform cylindrical Array Model schematic diagram in the present invention；
Fig. 3 is the uniform circular array model schematic in the present invention；
Fig. 4 is the real part schematic diagram of the present invention emulation single snap signal matrix X of array received；
Fig. 5 is present invention emulation no-coherence cumulating result vector selSchematic diagram；
Fig. 6 is present invention emulation even linear array final output zULAWith pitching phase compensation variable ηelChange curve signal
Figure；
Fig. 7 is simulation cycles convolution modulus value vector s of the present inventionazSchematic diagram；
Fig. 8 is present invention emulation uniform circular array matching output modulus value zUCAWith orientation phase compensation variable ηazChange curve
Schematic diagram.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
Referring to Fig.1, of the invention to be accomplished by
Step 1, the guiding matrix A of uniform cylindrical array is calculated.
Corresponding parameter is arranged according to the model of uniform cylindrical array in (1a), obtains constituting the nicely rounded of uniform cylindrical array
Battle array steering vector aUCAWith the steering vector a of even linear arrayULA:
Referring to Fig. 2 and Fig. 3, uniform cylindrical array can both regard the uniform circular array group put by multiple equidistant straight lines as
At can also regard as and be made of the even linear array that multiple equidistant circumference are put, the model of uniform cylindrical array includes: array
Operation wavelength λ, the element number of array M of uniform circular array, the array element spacing d of uniform circular array, even linear array element number of array N, uniformly
Array element spacing h, pitch angle of the signal source relative to uniform cylindrical array of linear arraySignal source is relative to uniform cylindrical array
Azimuth angle theta；
In the above parameter, it includes λ, M, d, N and h that the present invention, which needs the parameter being arranged, and wherein d needs to meet condition
H needs to meet condition
In the above parameter, present invention parameter to be estimated includesAnd θ, whereinSignal source direction is defined as in cylinder method
Plane projection and reference direction, i.e., the radial angle of No. 0 array element, value interval be [0,2 π)；θ is defined as signal source direction and circle
The angle of mast axis positive direction, value interval are [0, π]；
According to the model parameter of uniform cylindrical array, the steering vector a of uniform circular array is obtainedUCA, it is expressed as follows:
Wherein j is imaginary unit,
According to the model parameter of uniform cylindrical array, the steering vector a of even linear array is obtainedULA, it is expressed as follows:
Wherein
(1b) is according to the steering vector a of uniform circular arrayUCAWith the steering vector a of even linear arrayULA, obtain uniform cylindrical array
Guiding matrix A:
WhereinIndicate Kronecker multiplication, T expression turns order.
Step 2, the single snap signal of array received is lined up into matrix form X.
According to the guiding matrix A of uniform cylindrical array, the single snap signal of the array received is lined up into matrix formEach element respectively corresponds the array element accordingly numbered in A in X, wherein xmnIndicate that the n-th circle is nicely rounded
The reception signal of m-th of array element of battle array.
Step 3, pitch angle of the estimation signal source relative to uniform cylindrical array
Uniform cylindrical array is split into even linear array by (3a), i.e., every a line of X is individually carried out FFT transform, then will be each
No-coherence cumulating is done between row FFT transform result, obtains no-coherence cumulating result vector sel, it is embodied as sel=Σ | XF |,
WhereinIt is FFT transform matrix, Σ=[1]1×NIt is column summation vector, p and q are namely for referring to
Show the variable of row matrix coordinate and column coordinate；
(3b) is to no-coherence cumulating result vector selMaximum value detection is done, and calculates signal source relative to uniform cylindrical array
Pitch angleRough estimate evaluation
If qelFor no-coherence cumulating result vector selMaximum value element position, qelValued space be 0,1 ...,
N-1 }, according to no-coherence cumulating result vector selMaximum value element position qel, calculate pitch angleRough estimate evaluation
WhenWhen,
WhenWhen,
(3c) introduces pitching phase compensation variable ηel, to pitch angle rough estimate evaluationSection of closing on scan for and match
Correlation, and calculate pitch angleFine estimated value
Set the division number G of the pitching region of searchel, define pitching phase compensation variable ηel, then ηelValued space be
Defining even linear array matching factor vector is
If even linear array matching output vector isWherein * indicates conjugation；
If even linear array final output is zULA=| | yULA||1, searching for makes the zULAObtain the pitching phase compensation of maximum value
VariableAnd it willIt is defined as optimal pitching phase compensation variable, wherein | | | |11 model of vector is sought in expression
Number；
According to optimal pitching phase compensation variableWith no-coherence cumulating result vector selMaximum value element position qel,
Calculate pitch angleFine estimated value
WhenWhen,
WhenWhen,
Step 4, azimuth angle theta of the estimation signal source relative to uniform cylindrical array.
(4a) is according to the fine estimated value of pitch angleUniform cylindrical array is synthesized into a uniform circular array, i.e., it is signal is fast
It claps matrix X and obtains imaginary circles battle array signal y by matched transformUCA, it is embodied asWhereinIt is even linear array Optimum Matching coefficient vector,For optimal pitching phase compensation variable, * table
Show conjugation；
(4b) is to imaginary circles battle array signal yUCACyclic convolution and maximum value detection are carried out, and calculates signal source relative to uniform
The rough estimate evaluation of the azimuth angle theta of cylindrical-array
Define uniform circular array initial matching coefficient vectorWhereinIt is that signal source is opposite
In the pitch angle of uniform cylindrical arrayFine estimated value；
If cyclic convolution result modulus value vector is saz=| hUCA_0⊙yUCA|, wherein ⊙ indicates cyclic convolution, | | it indicates
Modulus value；
If qazFor cyclic convolution result modulus value vector sazMaximum value element position, qazValued space be 0,
1,...,M-1}；
According to cyclic convolution result modulus value vector sazMaximum value element position qaz, the rough estimate evaluation of computer azimuth angle θ
For
(4c) introduces orientation phase compensation variable ηaz, azimuthal rough estimate evaluationSection of closing on scan for and match
Correlation, and the fine estimated value of computer azimuth angle θ
The division number in bearing search section is set as Gaz, define orientation phase compensation variable ηaz, then ηazValued space
For
Defining uniform circular array matching factor vector is
If uniform circular array matching output modulus value result isSearch makes the zUCAObtain the side of maximum value
Position phase compensation variableAnd it willIt is defined as top optimization direction phase compensation variable, wherein H indicates that conjugation turns
Order；
According to top optimization direction phase compensation variableWith cyclic convolution result modulus value vector sazMaximum value element position
qaz, the fine estimated value of computer azimuth angle θ is
Effect of the invention further verifies explanation by following Matlab l-G simulation test.
(1) simulated conditions:
This emulation directly gives the relevant parameter parameter of array, specifically: the array element of operation wavelength λ=1m, uniform circular array
The array element of number M=128, the array element spacing d=0.45m of uniform circular array, the element number of array N=64 of even linear array, even linear array
Spacing h=0.5m.Preset parameter to be estimated are as follows: pitch angle of the signal source relative to uniform cylindrical arraySignal source
Azimuth angle theta=70 ° relative to uniform cylindrical array.
Additional parameter are as follows: each array element receives Signal to Noise Ratio (SNR)=20dB of signal.
When carrying out DOA estimation, the setting pitching region of search divides number Gel=100, bearing search interval division
Number Gaz=100.
The error threshold of set angle estimation is 0.1 °.
(2) emulation content and result:
Emulation 1, constructs the single snap signal matrix of array received according to the method for the present invention using above-mentioned simulation parameter
X, and the real part of matrix X is recorded, as a result such as Fig. 4.It can be observed how the phase relation of each array element signals from Fig. 4.
Emulation 2, is calculated no-coherence cumulating result vector s according to the method for the present invention using above-mentioned simulation parameterel, such as
Fig. 5.The available vector s from Fig. 5elMaximum value be 5806, the position of maximum value is 6.
Utilize vector selMaximum value position, further calculate to obtain pitch angleRough estimate evaluationAnd
Calculating pitch angle rough estimate error isIt can be seen that pitch angle rough estimate error is greater than error threshold, no
Meet the required precision of angular estimation.
Emulation 3, using above-mentioned simulation parameter, according to the method for the present invention, calculates and draws even linear array final output zULAWith
Pitching phase compensation variable ηelChange curve, as shown in Figure 6.It is available from Fig. 6, η at that timeelValue be 0.0432,
zULA8203 are maximized, i.e., optimal pitching phase compensation variable
It utilizesIt further calculates to obtain pitch angleFine estimated valueAnd it is fine to calculate pitch angle
Evaluated error isIt can be seen that the fine evaluated error of pitch angle is less than error threshold, meet angular estimation
Required precision.
Emulation 4, is calculated cyclic convolution result modulus value vector s according to the method for the present invention using above-mentioned simulation parameteraz,
And by vector sazEach element opsition dependent draw point by point, as shown in Figure 7.It is available from Fig. 7, vector sazMaximum value is
8009, maximum value position 25.
Utilize vector sazMaximum value position, further calculate the rough estimate evaluation for obtaining azimuth angle theta
And computer azimuth angle rough estimate error isRough estimate error in visible bearing angle is greater than error threshold,
It is unsatisfactory for the required precision of angular estimation.
Emulation 5, according to the method for the present invention, calculates using above-mentioned simulation parameter and draws uniform circular array matching output modulus value
zUCAWith orientation phase compensation variable ηazChange curve, as shown in Figure 8.It is available from Fig. 8 to work as ηazValue be-
When 0.00589, zUCAIt is maximized 8201, i.e. top optimization direction phase compensation variable
It utilizesIt further calculates to obtain the fine estimated value of azimuth angle thetaAnd computer azimuth angle is finely estimated
Counting error isThe fine evaluated error in visible bearing angle is less than error threshold, meets the essence of angular estimation
Degree requires.
Those of ordinary skill in the art will appreciate that: realize that all or part of the steps of above method embodiment can pass through
The relevant hardware of program instruction is completed, and program above-mentioned can store in computer-readable storage medium, which exists
When execution, step including the steps of the foregoing method embodiments is executed；And storage medium above-mentioned includes: ROM, RAM, magnetic or disk
Etc. the various media that can store program code.
The above description is merely a specific embodiment, but scope of protection of the present invention is not limited thereto, any
Those familiar with the art in the technical scope disclosed by the present invention, can easily think of the change or the replacement, and should all contain
Lid is within protection scope of the present invention.Therefore, protection scope of the present invention should be based on the protection scope of the described claims.
Claims (10)
1. Macroscale homogenous cylindrical array list snap DOA estimation method, which is characterized in that described method includes following steps:
(1) corresponding parameter is arranged according to the model of uniform cylindrical array, the uniform circular array for obtaining constituting uniform cylindrical array is led
To vector aUCAWith the steering vector a of even linear arrayULA, and according to uniform circular array steering vector aUCAIt is sweared with the guiding of even linear array
Measure aULAObtain the guiding matrix A of uniform cylindrical array；
(2) according to the guiding matrix A of uniform cylindrical array, the single snap signal of the array received is lined up into matrix form X=
[xmn]m∈[0,M-1],n∈[0,N-1], each element respectively corresponds the array element accordingly numbered, wherein xmnIndicate the of the n-th circle uniform circular array
The reception signal of m array element, M indicate that the array number that each uniform circular array contains, N indicate the array element that each even linear array contains
Number；
(3) according to signal snap matrix X, estimate pitch angle of the signal source relative to uniform cylindrical array
Uniform cylindrical array is split into even linear array by (3a), i.e., every a line of X is individually carried out FFT transform, then by each row FFT
No-coherence cumulating is done between transformation results, obtains no-coherence cumulating result vector sel；
(3b) is to no-coherence cumulating result vector selMaximum value detection is done, and calculates signal source bowing relative to uniform cylindrical array
The elevation angleRough estimate evaluation
(3c) introduces pitching phase compensation variable ηel, to pitch angle rough estimate evaluationSection of closing on scan for and match phase
It closes, and calculates pitch angleFine estimated value
(4) according to the fine estimated value of signal snap matrix X and pitch angleEstimate signal source relative to uniform cylindrical array
Azimuth angle theta:
(4a) is according to the fine estimated value of pitch angleUniform cylindrical array is synthesized into a uniform circular array, i.e., is become X through overmatching
Get imaginary circles battle array signal y in returnUCA；
(4b) is to imaginary circles battle array signal yUCACyclic convolution and maximum value detection are carried out, and calculates signal source relative to uniform cylindrical
The rough estimate evaluation of the azimuth angle theta of array
(4c) introduces orientation phase compensation variable ηaz, azimuthal rough estimate evaluationSection of closing on scan for and match phase
It closes, and the fine estimated value of computer azimuth angle θ
2. the method according to claim 1, wherein in (1) uniform cylindrical array uniform circular array steering vector
aUCA, it is expressed as follows:
Whereinλ is the operation wavelength of array, and M is the element number of array of uniform circular array, and d is uniform
The array element spacing of circle battle array, meets
3. the method according to claim 1, wherein the guiding of the even linear array of uniform cylindrical array is sweared in (1)
Measure aULA, it is expressed as follows:
Whereinλ is the operation wavelength of array, and N is the element number of array of even linear array, and h is the array element spacing of even linear array,
Meet
4. the method according to claim 1, wherein in (1) uniform cylindrical array guiding matrix A, indicate such as
Under:
WhereinIndicate Kronecker multiplication, T expression turns order.
5. the method according to claim 1, wherein no-coherence cumulating result vector s in (3a)el, it is expressed as follows:
sel=Σ | XF |,
Wherein X is signal snap matrix,It is FFT transform matrix, N is the array element of even linear array
Number, Σ=[1]1×NIt is column summation vector.
6. the method according to claim 1, wherein (3b), is accomplished by
If qelFor no-coherence cumulating result vector selMaximum value element position, qelValued space be 0,1 ..., N-
1 }, wherein N is the element number of array of even linear array；
According to no-coherence cumulating result vector selMaximum value element position qel, calculate pitch angleRough estimate evaluation
WhenWhen,
WhenWhen,
Whereinλ is the operation wavelength of array, and h is the array element spacing of even linear array, is met
7. the method according to claim 1, wherein (3c), is accomplished by
Set the division number G of the pitching region of searchel, define pitching phase compensation variable ηel, then ηelValued space beWherein, N is the element number of array of even linear array；
Defining even linear array matching factor vector isWherein qelFor no-coherence cumulating result to
Measure selMaximum value element position；
If even linear array matching output vector isWherein X is signal snap matrix, and * indicates conjugation；
If even linear array final output is zULA=| | yULA||1, searching for makes the zULAObtain the pitching phase compensation variable of maximum valueAnd it willIt is defined as optimal pitching phase compensation variable, wherein | | | |11 norm of vector is sought in expression；
According to optimal pitching phase compensation variableWith no-coherence cumulating result vector selMaximum value element position qel, calculate
Pitch angleFine estimated value
WhenWhen,
WhenWhen,
Whereinλ is the operation wavelength of array, and h is the array element spacing of even linear array, is met
8. the method according to claim 1, wherein the imaginary circles battle array signal y in (4a)UCA, it is expressed as follows:
Wherein X is signal snap matrix, and * indicates conjugation,It is even linear array Optimum Matching coefficient
Vector, N are the element number of array of even linear array, qelFor no-coherence cumulating result vector selMaximum value element position,For most
Excellent pitching phase compensation variable.
9. the method according to claim 1, wherein (4b), is accomplished by
Define uniform circular array initial matching coefficient vectorWhereinIt is signal source relative to equal
The pitch angle of even cylindrical-arrayFine estimated value, λ is the operation wavelength of array, and M is equal
The element number of array of nicely rounded battle array, d are the array element spacing of uniform circular array, are met
If cyclic convolution result modulus value vector is saz=| hUCA_0⊙yUCA|, wherein yUCAIt is imaginary circles battle array signal, ⊙ indicates circulation
Convolution, | | indicate modulus value；
If qazFor cyclic convolution result modulus value vector sazMaximum value element position, qazValued space be 0,1 ..., M-
1 }, wherein M is the element number of array of uniform circular array；
According to cyclic convolution result modulus value vector sazMaximum value element position qaz, the rough estimate evaluation of computer azimuth angle θ is
10. the method according to claim 1, wherein (4c), is accomplished by
The division number in bearing search section is set as Gaz, define orientation phase compensation variable ηaz, then ηazValued space beWherein, M is the element number of array of uniform circular array；
Defining uniform circular array matching factor vector isWhereinIt is that signal source is opposite
In the pitch angle of uniform cylindrical arrayFine estimated value, qazFor cyclic convolution result modulus value vector sazMaximum value element
Position,λ is the operation wavelength of array, and d is the array element spacing of uniform circular array, is met
If uniform circular array matching output modulus value result isSearch makes the zUCAObtain the orientation phase of maximum value
Position compensation variableAnd it willIt is defined as top optimization direction phase compensation variable, wherein yUCAIt is imaginary circles battle array letter
Number, H indicates that conjugation turns order, | | indicate modulus value；
According to top optimization direction phase compensation variableWith cyclic convolution result modulus value vector sazMaximum value element position qaz,
The fine estimated value of computer azimuth angle θ is
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811006677.0A CN109188346B (en) | 2018-08-31 | 2018-08-31 | Single snapshot DOA estimation method for large-scale uniform cylindrical array |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811006677.0A CN109188346B (en) | 2018-08-31 | 2018-08-31 | Single snapshot DOA estimation method for large-scale uniform cylindrical array |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109188346A true CN109188346A (en) | 2019-01-11 |
| CN109188346B CN109188346B (en) | 2023-03-10 |
Family
ID=64917050
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811006677.0A Active CN109188346B (en) | 2018-08-31 | 2018-08-31 | Single snapshot DOA estimation method for large-scale uniform cylindrical array |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109188346B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111368256A (en) * | 2020-03-23 | 2020-07-03 | 电子科技大学 | Single snapshot direction finding method based on uniform circular array |
| CN112130111A (en) * | 2020-09-22 | 2020-12-25 | 南京航空航天大学 | Single-snapshot two-dimensional DOA estimation method for large-scale uniform cross array |
| CN112346005A (en) * | 2020-10-30 | 2021-02-09 | 哈尔滨工程大学 | Airspace rotation orientation estimation method applied to uniform circular hydrophone array |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101273389A (en) * | 2006-02-21 | 2008-09-24 | 中兴通讯股份有限公司 | Channel estimation method and device based on array antenna |
| CN102175989A (en) * | 2010-12-30 | 2011-09-07 | 电子科技大学 | Method for measuring incoherently distributed signal two-dimensional DOA (direction of arrival) |
| CN106507952B (en) * | 2011-08-04 | 2014-04-23 | 西南电子电信技术研究所 | A kind of quick spatial spectrum computational methods based on circle battle array |
| WO2015109869A1 (en) * | 2014-01-24 | 2015-07-30 | 深圳大学 | High resolution doa estimation method and system |
| CN104931923A (en) * | 2015-04-02 | 2015-09-23 | 刘松 | Grid iterative estimation of signal parameters via rotational invariance techniques (ESPRIT), namely, extensible rapid estimation algorithm capable of being used for uniform circular array 2-dimensional direction of arrival (2D DOA) |
| US20160219365A1 (en) * | 2013-07-24 | 2016-07-28 | Mh Acoustics, Llc | Adaptive Beamforming for Eigenbeamforming Microphone Arrays |
| CN107907853A (en) * | 2017-11-07 | 2018-04-13 | 中国人民解放军信息工程大学 | A kind of single distributed source DOA estimation method based on uniform circular array differential phase |
| CN107918108A (en) * | 2017-11-14 | 2018-04-17 | 重庆邮电大学 | A kind of uniform circular array 2-d direction finding method for quick estimating |
-
2018
- 2018-08-31 CN CN201811006677.0A patent/CN109188346B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101273389A (en) * | 2006-02-21 | 2008-09-24 | 中兴通讯股份有限公司 | Channel estimation method and device based on array antenna |
| CN102175989A (en) * | 2010-12-30 | 2011-09-07 | 电子科技大学 | Method for measuring incoherently distributed signal two-dimensional DOA (direction of arrival) |
| CN106507952B (en) * | 2011-08-04 | 2014-04-23 | 西南电子电信技术研究所 | A kind of quick spatial spectrum computational methods based on circle battle array |
| US20160219365A1 (en) * | 2013-07-24 | 2016-07-28 | Mh Acoustics, Llc | Adaptive Beamforming for Eigenbeamforming Microphone Arrays |
| WO2015109869A1 (en) * | 2014-01-24 | 2015-07-30 | 深圳大学 | High resolution doa estimation method and system |
| CN104931923A (en) * | 2015-04-02 | 2015-09-23 | 刘松 | Grid iterative estimation of signal parameters via rotational invariance techniques (ESPRIT), namely, extensible rapid estimation algorithm capable of being used for uniform circular array 2-dimensional direction of arrival (2D DOA) |
| CN107907853A (en) * | 2017-11-07 | 2018-04-13 | 中国人民解放军信息工程大学 | A kind of single distributed source DOA estimation method based on uniform circular array differential phase |
| CN107918108A (en) * | 2017-11-14 | 2018-04-17 | 重庆邮电大学 | A kind of uniform circular array 2-d direction finding method for quick estimating |
Non-Patent Citations (9)
| Title |
|---|
| ROALD GOOSSENS ET AL.: "UCA Root-MUSICWith Sparse Uniform Circular Arrays", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 * |
| YUGUAN HOU: "An iterative optimization for mutual coupling correction of UCA with single snapshot", 《IET INTERNATIONAL RADAR CONFERENCE 2015》 * |
| 刘艳等: "基于均匀圆阵的改进UCA-ESPRIT算法", 《计算机科学》 * |
| 施思等: "应用于波束域的改进TOFS算法", 《雷达科学与技术》 * |
| 杨永建等: "基于MUSIC算法的圆柱共形阵DOA 估计", 《空军工程大学学报( 自然科学版)》 * |
| 毛维平等: "均匀圆阵相干信源二维波达方向估计", 《***工程与电子技术》 * |
| 王洪雁等: "非均匀噪声条件下考虑互耦效应的DOA估计方法", 《吉林大学学报(工学版)》 * |
| 赵鸿图等: "循环卷积的时域与频域算法研究", 《计算机工程与设计》 * |
| 黄秀琼等: "基于循环特性的均匀圆阵快速空间谱计算", 《电讯技术》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111368256A (en) * | 2020-03-23 | 2020-07-03 | 电子科技大学 | Single snapshot direction finding method based on uniform circular array |
| CN111368256B (en) * | 2020-03-23 | 2023-03-03 | 电子科技大学 | Single snapshot direction finding method based on uniform circular array |
| CN112130111A (en) * | 2020-09-22 | 2020-12-25 | 南京航空航天大学 | Single-snapshot two-dimensional DOA estimation method for large-scale uniform cross array |
| CN112346005A (en) * | 2020-10-30 | 2021-02-09 | 哈尔滨工程大学 | Airspace rotation orientation estimation method applied to uniform circular hydrophone array |
| CN112346005B (en) * | 2020-10-30 | 2022-07-12 | 哈尔滨工程大学 | Airspace rotation orientation estimation method applied to uniform circular hydrophone array |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109188346B (en) | 2023-03-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109633525B (en) | Quick uniform circular array space spectrum direction finding realization method | |
| CN106019215B (en) | Nested array direction of arrival angle method of estimation based on fourth-order cumulant | |
| CN109188346A (en) | Macroscale homogenous cylindrical array list snap DOA estimation method | |
| CN108663653B (en) | Direction-of-arrival estimation method based on L-shaped electromagnetic vector sensor array | |
| CN104391183A (en) | Near-field-measurement-based rapid calculation method of antenna far-field region characteristic | |
| CN110197112B (en) | Beam domain Root-MUSIC method based on covariance correction | |
| CN110412514B (en) | Linear frequency modulation continuous wave waveform optimization method under MIMO system | |
| CN110161452B (en) | Direction-of-arrival estimation method based on cross-prime L-shaped electromagnetic vector sensor array | |
| CN105445718B (en) | A kind of DOA estimation method of the distributed overloading warship over-the-horizon radar based on array reconfiguration | |
| CN103323845B (en) | Image inversion method of non-uniform sampling comprehensive bore diameter radiometer | |
| CN109633520B (en) | Uniform circular array super-resolution spatial spectrum estimation method | |
| CN106970348B (en) | Electromagnetic Vector Sensor Array decorrelation LMS two dimension MUSIC method for parameter estimation | |
| CN108663654B (en) | 360-degree all-dimensional dynamic direction finding method based on continuous quantum pigeon group | |
| CN107576951B (en) | Direction-of-arrival estimation method based on nested electromagnetic vector sensor array | |
| Wang et al. | Linear frequency diverse array manifold geometry and ambiguity analysis | |
| CN106405487A (en) | General spatial spectrum estimation method based on extended ESPRIT | |
| CN105022026A (en) | Two-dimensional arrival angle estimation method of L-shaped array | |
| CN108535698A (en) | The low elevation estimate method of metre wave radar based on beam space | |
| CN104931923A (en) | Grid iterative estimation of signal parameters via rotational invariance techniques (ESPRIT), namely, extensible rapid estimation algorithm capable of being used for uniform circular array 2-dimensional direction of arrival (2D DOA) | |
| CN104502904A (en) | Torpedo homing beam sharpening method | |
| CN107037398A (en) | A kind of two-dimentional MUSIC algorithms estimate the parallel calculating method of direction of arrival | |
| CN106443571A (en) | Conformal mapping estimation of signal parameters via rotational invariance technique DOA estimation method | |
| CN109212466B (en) | Quantum dragonfly evolution mechanism-based broadband direction finding method | |
| CN102062851A (en) | Direction finding method based on improved L array star-carrying broadband multipurpose | |
| CN110445557A (en) | Deep space TT&C interferometry large aperture antenna is directed toward Calibration Method and device |
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 |