CN109143155A - Coherent signal Wave arrival direction estimating method and system based on mutual pixel array - Google Patents
Coherent signal Wave arrival direction estimating method and system based on mutual pixel array Download PDFInfo
- Publication number
- CN109143155A CN109143155A CN201810847429.2A CN201810847429A CN109143155A CN 109143155 A CN109143155 A CN 109143155A CN 201810847429 A CN201810847429 A CN 201810847429A CN 109143155 A CN109143155 A CN 109143155A
- Authority
- CN
- China
- Prior art keywords
- array
- signal
- target aerial
- subarray
- aerial array
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
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
A kind of coherent signal Wave arrival direction estimating method and system based on mutual pixel array provided by the invention, the signal received to mutual pixel array carry out sampling and obtain time domain discrete signal;Again by the array number of setting subarray, the fourth order cumulant of sub-array signal is obtained according to time domain discrete signal, and the fourth order cumulant of sub-array signal is subjected to front and back to obtaining smoothed out fourth order cumulant after smoothing processing;Eigenvalues Decomposition is carried out to smoothed out fourth order cumulant and obtains projection operator, the final direction of arrival for estimating signal according to all projection operators using Root-MUSIC algorithm.This method and system, which pass through, combines mutual pixel array and fourth order cumulant to carry out signal Mutual coupling, and the precision of estimated result can be effectively ensured in the case where not increasing array element number, has preferable effect when signal number is more than array element number;And can effectively ring off signal it is related, enhance the spatial resolution to signal, various complex environments can be suitable for.
Description
Technical field
The present invention relates to signal processing technology fields, more particularly, to a kind of coherent signal wave based on mutual pixel array
Arrival direction estimating method and system.
Background technique
In array signal processing, direction of arrival (Direction-of-Arrival, DOA) estimation is a kind of determining signal
The effective method of information source position is estimated in direction, has application in fields such as radar, communication, medicine.In DOA estimation field,
It is obtained with the subspace method that multiple signal classification (MUltiple Signal Classification, MUSIC) algorithm is representative
To being widely applied.MUSIC algorithm is decomposed by the correlation matrix to signal, is obtained orthogonal signal subspace and is made an uproar
Phonon space, the characteristic with airspace super-resolution, therefore be widely used.
However, in traditional MUSIC algorithm, it usually needs a large amount of sampled data, to guarantee the estimated accuracy of algorithm;
Simultaneously in order to obtain a large amount of sampled data, it need to usually increase the quantity of array element.Although can be mentioned by the quantity for increasing array element
The estimated accuracy of high MUSIC algorithm, but this mode not only increases reception data volume, and to data transmission, storage and place
Reason brings huge pressure, so that the realization difficulty of hardware increases, it is difficult to meet actual engineering demand.
In view of this, it is urgent to provide the signals that one kind can ensure estimated accuracy in the case where not increasing array element quantity
Wave arrival direction estimating method and system.
Summary of the invention
The present invention is usual using traditional DOA estimation method progress signal Mutual coupling in the prior art in order to overcome
The quantity of array element need to be increased to ensure estimated accuracy, cause to be difficult to the problem of being suitable for complicated applications environment, one kind is provided and is based on
The coherent signal Wave arrival direction estimating method and system of mutual pixel array.
On the one hand, the present invention provides a kind of coherent signal Wave arrival direction estimating method based on mutual pixel array, comprising:
It is sampled according to signal received by Target Aerial Array of the time domain to preset quantity, obtains each Target Aerial Array pair
The time domain discrete signal answered, wherein every two Target Aerial Array mutual pixel array each other;
For any one Target Aerial Array, the array number of the subarray of the Target Aerial Array is set, according to the array element of subarray
Number obtains the corresponding sub-array signal of all subarrays of the Target Aerial Array from the corresponding time domain discrete signal of the Target Aerial Array;
For any two Target Aerial Array, obtained according to the corresponding sub-array signal of all subarrays of two Target Aerial Arrays
Obtain the corresponding smoothed out fourth order cumulant matrix of two Target Aerial Arrays;
Eigenvalues Decomposition is carried out to each smoothed out fourth order cumulant matrix and obtains each smoothed out fourth order cumulant
The corresponding projection operator of matrix;
The direction of arrival of signal is estimated according to all projection operators using Root-MUSIC algorithm.
Preferably, which is obtained from the corresponding time domain discrete signal of the Target Aerial Array according to the array number of subarray
The corresponding sub-array signal of all subarrays of column, specifically:
For any one Target Aerial Array, using the array number of the corresponding subarray of the Target Aerial Array as destination number,
In all array elements of the Target Aerial Array, the adjacent array element of destination number is combined, obtains all submatrixs of the Target Aerial Array
Column;
The corresponding sub-array signal of each subarray is isolated from the corresponding time domain discrete signal of the Target Aerial Array, is obtained
The corresponding sub-array signal of all subarrays of the Target Aerial Array.
Preferably, which is obtained according to the corresponding sub-array signal of all subarrays of two Target Aerial Arrays
Corresponding smoothed out fourth order cumulant matrix is arranged, specifically:
For any two Target Aerial Array, a subarray is randomly selected from two respective subarrays of Target Aerial Array
Combination of two is carried out, multiple subarray combinations is obtained, any one subarray is combined, according to the institute in subarray combination
There is the corresponding sub-array signal of subarray to obtain the subarray and combines corresponding fourth order cumulant matrix;
It all subarrays combines corresponding fourth order cumulant matrix and carries out front and back to smoothing processing, obtain two targets
The corresponding smoothed out fourth order cumulant matrix of array.
Preferably, subarray combination is obtained according to the corresponding sub-array signal of all subarrays in subarray combination
Corresponding fourth order cumulant matrix, specific formula for calculation are as follows:
Wherein, x1aIt (t) is the corresponding sub-array signal of subarray a of Target Aerial Array 1;x2bIt (t) is the son of Target Aerial Array 2
The corresponding sub-array signal of array b;RabCorresponding fourth order cumulant square is combined for the subarray that subarray a and subarray b is constituted
Battle array.
Preferably, it all subarrays combines corresponding fourth order cumulant matrix and carries out front and back to smoothing processing, be somebody's turn to do
The corresponding smoothed out fourth order cumulant matrix of two Target Aerial Arrays, specific formula for calculation are as follows:
Wherein, R12For Target Aerial Array 1 and the corresponding smoothed out fourth order cumulant of Target Aerial Array 2;A is Target Aerial Array 1
Any one subarray;B is any one subarray of Target Aerial Array 2;N1For the array number of Target Aerial Array 1;M1For Target Aerial Array
The array number of 1 each subarray;N2For the array number of Target Aerial Array 2;M2For the array number of each subarray of Target Aerial Array 2;
J is the permutation matrix that opposition angle array element element is 1.
Preferably, Eigenvalues Decomposition is carried out to each smoothed out fourth order cumulant matrix and obtains each smoothed out quadravalence
The corresponding projection operator of cumulant matrices, specifically:
For any one smoothed out fourth order cumulant matrix, feature is carried out to the smoothed out fourth order cumulant matrix
Value is decomposed, and multiple feature vectors are obtained;
The feature vector for meeting preset condition is filtered out from all feature vectors according to the characteristic value of all feature vectors,
As target feature vector;
All target feature vectors are combined into objective matrix, which is obtained according to objective matrix
The corresponding projection operator of matrix.
Preferably, the direction of arrival of signal is estimated according to all projection operators using Root-MUSIC algorithm, specifically:
For any one projection operator, the corresponding multinomial letter of the projection operator is constructed using preset function building method
Number;
Target equation is established according to the corresponding polynomial function of all projection operators, target equation is solved, is distributed in
Equattion root on unit circle;
The direction of arrival of signal is estimated according to equattion root.
On the one hand, the present invention provides a kind of coherent signal Mutual coupling system based on mutual pixel array, comprising:
Signal acquisition module, for being sampled according to signal received by Target Aerial Array of the time domain to preset quantity,
The corresponding time domain discrete signal of each Target Aerial Array is obtained, wherein every two Target Aerial Array mutual pixel array each other;
Signal pre-processing module, for setting the array element of the subarray of the Target Aerial Array for any one Target Aerial Array
Number, all submatrixs of the Target Aerial Array are obtained according to the array number of subarray from the corresponding time domain discrete signal of the Target Aerial Array
Arrange corresponding sub-array signal;
Smoothing computation module is used for for any two Target Aerial Array, according to all subarrays of two Target Aerial Arrays
Corresponding sub-array signal obtains the corresponding smoothed out fourth order cumulant matrix of two Target Aerial Arrays;
Eigenvalues Decomposition module obtains each for carrying out Eigenvalues Decomposition to each smoothed out fourth order cumulant matrix
The corresponding projection operator of smoothed out fourth order cumulant matrix;
Mutual coupling module, for estimating that the wave of signal reaches according to all projection operators using Root-MUSIC algorithm
Direction.
On the one hand, the present invention provides a kind of electronic equipment, comprising:
At least one processor;And
At least one processor being connect with the processor communication, in which:
The memory is stored with the program instruction that can be executed by the processor, and the processor calls described program to refer to
It enables and is able to carry out any of the above-described method.
On the one hand, the present invention provides a kind of non-transient computer readable storage medium, and the non-transient computer is readable to deposit
Storage media stores computer instruction, and the computer instruction makes the computer execute any of the above-described method.
A kind of coherent signal Wave arrival direction estimating method and system based on mutual pixel array provided by the invention, to mutual primitive matrix
It arranges received signal and carries out sampling acquisition time domain discrete signal;Again by setting subarray array number, according to time domain from
Scattered signal obtains the fourth order cumulant of sub-array signal, and the fourth order cumulant of sub-array signal is carried out front and back to after smoothing processing
Obtain smoothed out fourth order cumulant;Eigenvalues Decomposition is carried out to smoothed out fourth order cumulant and obtains projection operator, final benefit
The direction of arrival of signal is estimated according to all projection operators with Root-MUSIC algorithm.This method and system, which pass through, combines mutual primitive matrix
Column and fourth order cumulant carry out signal Mutual coupling, can effectively ensure to estimate knot in the case where not increasing array element number
The precision of fruit has preferable effect when signal number is more than array element number;And can effectively ring off signal it is related, enhance
To the spatial resolution of signal, various complex environments can be suitable for, can effectively meet various Practical Project demands.
Detailed description of the invention
Fig. 1 is a kind of bulk flow of coherent signal Wave arrival direction estimating method based on mutual pixel array of the embodiment of the present invention
Journey schematic diagram;
Fig. 2 is a kind of whole knot of coherent signal Mutual coupling system based on mutual pixel array of the embodiment of the present invention
Structure schematic diagram;
Fig. 3 is the structural framing schematic diagram of a kind of electronic equipment of the embodiment of the present invention.
Specific embodiment
With reference to the accompanying drawings and examples, specific embodiments of the present invention will be described in further detail.Implement below
Example is not intended to limit the scope of the invention for illustrating the present invention.
Fig. 1 is a kind of bulk flow of coherent signal Wave arrival direction estimating method based on mutual pixel array of the embodiment of the present invention
Journey schematic diagram, as shown in Figure 1, the present invention provides a kind of coherent signal Wave arrival direction estimating method based on mutual pixel array, comprising:
S1 is sampled according to signal received by Target Aerial Array of the time domain to preset quantity, obtains each target array
Corresponding time domain discrete signal is arranged, wherein every two Target Aerial Array mutual pixel array each other;
Specifically, the Target Aerial Array of preset quantity is arranged first, and wherein Target Aerial Array can be the array day of reception signal
Line, and every two Target Aerial Array mutual pixel array each other.Mutual pixel array is a kind of thinned array with specific arrangement mode, it can
To construct virtual uniform array, there is preferable effect when signal number is more than array element number.For example, a pair of mutually pixel array
Array number is respectively A and B, then this pair of coprime array can construct the virtual array that an array number is A*B, therefore this is a pair of
The signal number that mutual pixel array is capable of handling is A*B, and A*B is greater than A+B, and the signal number that can be handled is more than array element sum.
Therefore, the coherent signal Wave arrival direction estimating method phase based on mutual pixel array in the case where identical array element, in the present embodiment
More signals are capable of handling compared with traditional DOA estimation method.
In the present embodiment, 3 Target Aerial Arrays, respectively Target Aerial Array 1, Target Aerial Array 2 and Target Aerial Array are arranged in total
3, and the array element quantity of Target Aerial Array 1, Target Aerial Array 2 and Target Aerial Array 3 is respectively N1、N2And N3.In order to enable Target Aerial Array 1,
The location arrangements of Target Aerial Array 1 can be by every two Target Aerial Array mutual pixel array each other in Target Aerial Array 2 and Target Aerial Array 3
{An1d|n1=0 ..., N1-1};It is { Bn by the location arrangements of Target Aerial Array 22d|n2=0 ..., N2-1};By Target Aerial Array 3
Location arrangements be { Cn3d|n3=0 ..., N3-1}.The array element spacing of Target Aerial Array 1 is A*d as a result,;The battle array of Target Aerial Array 2
First spacing is B*d;The array element spacing of Target Aerial Array 3 is C*d.Wherein, A, B and C are positive integer coprime two-by-two;D is between array element
Away from a basic unit, be usually no more than half-wavelength.In other embodiments, the preset quantity of Target Aerial Array can basis
Actual demand is configured, and is not specifically limited herein.
After the Target Aerial Array for arranging preset quantity, signal is received using the Target Aerial Array of preset quantity.Then it presses again
It is sampled according to signal received by Target Aerial Array of the time domain to preset quantity, certain pretreatment is carried out to sampled signal
Afterwards, it can be obtained the corresponding time domain discrete signal of each Target Aerial Array.Wherein, the number of snapshots of sampling can according to actual needs into
Row setting, is not specifically limited herein.
S2 sets the array number of the subarray of the Target Aerial Array for any one Target Aerial Array, according to the battle array of subarray
First number obtains the corresponding subarray letter of all subarrays of the Target Aerial Array from the corresponding time domain discrete signal of the Target Aerial Array
Number;
Specifically, based on the above technical solution, for any one target in the Target Aerial Array of preset quantity
Array sets the array number of the subarray of the Target Aerial Array, and the array number of each subarray of same Target Aerial Array is identical, submatrix
The array number of column should be less than the array element total quantity of the Target Aerial Array, can specifically be configured, not do according to actual needs herein
It is specific to limit.It is all subarrays that can determine in the Target Aerial Array according to the array number of subarray, finally again from the target array
Arrange the corresponding sub-array signal of all subarrays (institute i.e. in subarray that the Target Aerial Array is obtained in corresponding time domain discrete signal
There is the signal that array element receives).Similarly, according to the method described above step can get other Target Aerial Arrays all subarrays it is corresponding
Sub-array signal.
S3 believes any two Target Aerial Array according to the corresponding subarray of all subarrays of two Target Aerial Arrays
Number obtain the corresponding smoothed out fourth order cumulant matrix of two Target Aerial Arrays;
Specifically, after the corresponding sub-array signal of all subarrays for obtaining each Target Aerial Array, for any two
Target Aerial Array randomly selects a subarray from two respective subarrays of Target Aerial Array and carries out combination of two, thus needle
Multiple subarray combinations are obtained to two Target Aerial Arrays, and include two subarrays in the combination of each subarray.For multiple
Any one subarray combination in subarray combination, obtains two corresponding submatrixs of subarray in subarray combination
The corresponding sub-array signal of two subarrays is carried out fourth order cumulant and calculates acquisition subarray combination correspondence by column signal
Fourth order cumulant matrix.On this basis, it combines all subarrays before and after corresponding fourth order cumulant matrix carries out to flat
Sliding processing, can be obtained the corresponding smoothed out fourth order cumulant matrix of two Target Aerial Arrays.Similarly, it is walked according to the above method
Suddenly the corresponding smoothed out fourth order cumulant matrix of any other two Target Aerial Arrays can be obtained.Smoothed out fourth order cumulant
Matrix effectively relieves signal correlation.
It should be noted that fourth order cumulant illustrates the statistical relationship between the signal of four array element, compared to than two
Rank statistic, fourth order cumulant contain more information, are conducive to the effective aperture for increasing array, increase spatial resolution.
Also, fourth order cumulant can also inhibit Gaussian noise, have preferable effect in the coloured noise environment of Gaussian Profile.
S4 carries out Eigenvalues Decomposition to each smoothed out fourth order cumulant matrix and obtains each smoothed out fourth-order cumulant
The corresponding projection operator of moment matrix;
Specifically, based on the above technical solution, it for any two Target Aerial Array, can get corresponding smooth
Fourth order cumulant matrix afterwards can be obtained multiple smoothed out fourth order cumulant matrix as a result,.On this basis, for appointing
A smoothed out fourth order cumulant matrix of anticipating obtains more smoothed out fourth order cumulant matrix progress Eigenvalues Decomposition
A feature vector, then all feature vectors are arranged according to the order of magnitude of characteristic value, isolate biggish characteristic value
With lesser characteristic value, wherein biggish characteristic value is generally 10 times or more of smaller characteristic value.Finally, according to smaller characteristic value
The matrix that character pair vector is constituted can be obtained the corresponding projection operator of the smoothed out fourth order cumulant matrix.Similarly, it presses
It can get the corresponding projection operator of other smoothed out fourth order cumulant matrix according to above method step.
S5 estimates the direction of arrival of signal using Root-MUSIC algorithm according to all projection operators.
Specifically, after obtaining the corresponding projection operator of each smoothed out fourth order cumulant matrix, Root- is utilized
MUSIC algorithm estimates the direction of arrival of signal according to all projection operators.Root-MUSIC algorithm specifically: construction is each first
The corresponding polynomial function of projection operator establishes equation further according to the corresponding polynomial function of all projection operators, and to the party
Cheng Jinhang is solved, and obtains the equattion root being distributed on unit circle, the direction of arrival of signal can be finally estimated according to equattion root.
A kind of coherent signal Wave arrival direction estimating method based on mutual pixel array provided by the invention, connects mutual pixel array
The signal being subject to carries out sampling and obtains time domain discrete signal;Again by the array number of setting subarray, according to time domain discrete signal
The fourth order cumulant of sub-array signal is carried out front and back to being put down after smoothing processing by the fourth order cumulant for obtaining sub-array signal
Fourth order cumulant after cunning;Eigenvalues Decomposition is carried out to smoothed out fourth order cumulant and obtains projection operator, finally utilizes Root-
MUSIC algorithm estimates the direction of arrival of signal according to all projection operators.This method, which passes through, combines mutual pixel array and fourth-order cumulant
Amount carries out signal Mutual coupling, and the precision of estimated result can be effectively ensured in the case where not increasing array element number,
Signal number has preferable effect when being more than array element number;And can effectively ring off signal it is related, enhance the sky to signal
Between resolving power, can be suitable for various complex environments, can effectively meet various Practical Project demands.
Based on any of the above-described embodiment, a kind of coherent signal Wave arrival direction estimating method based on mutual pixel array, root are provided
All subarrays pair of the Target Aerial Array are obtained from the corresponding time domain discrete signal of the Target Aerial Array according to the array number of subarray
The sub-array signal answered, specifically: for any one Target Aerial Array, the array number of the corresponding subarray of the Target Aerial Array is made
The adjacent array element of destination number is combined, the target array is obtained in all array elements of the Target Aerial Array for destination number
All subarrays of column;The corresponding subarray letter of each subarray is isolated from the corresponding time domain discrete signal of the Target Aerial Array
Number, obtain the corresponding sub-array signal of all subarrays of the Target Aerial Array.
Specifically, in the present embodiment, for any one Target Aerial Array, the array element of the subarray of the Target Aerial Array is set
Number, then using the array number of subarray as destination number, in all array elements of the Target Aerial Array, by the adjacent battle array of destination number
Member is combined, and can be obtained all subarrays of the Target Aerial Array.On this basis, then from the corresponding time domain of the Target Aerial Array
The corresponding sub-array signal of each subarray is isolated in discrete signal, all subarrays that can be obtained the Target Aerial Array are corresponding
Sub-array signal.
For example, for Target Aerial Array 1, Target Aerial Array 2 and Target Aerial Array 3, to Target Aerial Array 1, Target Aerial Array 2 and mesh
It marks signal received by array 3 and carries out TsThe sampling of a snap can obtain Target Aerial Array 1, Target Aerial Array 2 and target array respectively
The corresponding time domain discrete signal x of column 31(t)、x2(t) and x3(t), x1(t)、x2(t) and x3(t) it may particularly denote are as follows:
x1(t)=[x1(0, t) ..., x1(N1-1,t)]T, t=1 ..., Ts;
x2(t)=[x2(0, t) ..., x2(N2- 1, t)]T, t=1 ..., Ts;
x3(t)=[x3(0, t) ..., x3(N3-1,t)]T, t=1 ..., Ts;
Wherein, N1For the array element quantity of Target Aerial Array 1;N2For the array element quantity of Target Aerial Array 2;N3For the battle array of Target Aerial Array 3
First quantity.
Further, the array number of the subarray of Target Aerial Array 1 is set as M1, by the battle array of the subarray of Target Aerial Array 2
First number is set as M2, the array number of the subarray of Target Aerial Array 3 is set as M3.On this basis, for Target Aerial Array 1
Obtain N1-M1+ 1 subarray;It can be obtained N for Target Aerial Array 22-M2+ 1 subarray;It can be obtained for Target Aerial Array 3
N3-M3+ 1 subarray.Finally, from Target Aerial Array 1, Target Aerial Array 2 and the corresponding time domain discrete signal x of Target Aerial Array 31(t)、
x2(t) and x3In 9t), it respectively is isolated by the corresponding submatrix of each subarray of Target Aerial Array 1, Target Aerial Array 2 and Target Aerial Array 3
Column signal x1a(t)、x2b(t) and x3c(t), it may particularly denote are as follows:
x1a(t)=[x1(a,t),...,x1(a+M1-1,t)]T, a=0 ..., N1-M1;
x2b(t)=[x2(b,t),...,x2(b+M2-1,t)]T, b=0 ..., N2-M2;
x3c(t)=[x3(c,t),...,x3(c+M3-1,t)]T, c=0 ..., N3-M3;
Wherein, a is any one subarray in Target Aerial Array 1;B is any one subarray in Target Aerial Array 2;c
For any one subarray in Target Aerial Array 3.
A kind of coherent signal Wave arrival direction estimating method based on mutual pixel array provided by the invention, for any one mesh
Array is marked, using the array number of the corresponding subarray of the Target Aerial Array as destination number, in all array elements of the Target Aerial Array,
The adjacent array element of destination number is combined, all subarrays of the Target Aerial Array are obtained;When corresponding from the Target Aerial Array
The corresponding sub-array signal of each subarray is isolated in the discrete signal of domain, all subarrays for obtaining the Target Aerial Array are corresponding
Sub-array signal is conducive to the fourth order cumulant for combining all sub-array signals to obtain sub-array signal, and then is conducive to enhance
To the spatial resolution of signal.
Based on any of the above-described embodiment, a kind of coherent signal Wave arrival direction estimating method based on mutual pixel array, root are provided
It is corresponding smoothed out that two Target Aerial Arrays are obtained according to the corresponding sub-array signal of all subarrays of two Target Aerial Arrays
Fourth order cumulant matrix, specifically: it is random from two respective subarrays of Target Aerial Array for any two Target Aerial Array
It chooses a subarray and carries out combination of two, obtain multiple subarray combinations, any one subarray is combined, according to the son
The corresponding sub-array signal of all subarrays in array combination obtains the subarray and combines corresponding fourth order cumulant matrix;It will
All subarrays combine corresponding fourth order cumulant matrix and carry out front and back to smoothing processing, and it is corresponding to obtain two Target Aerial Arrays
Smoothed out fourth order cumulant matrix.
Specifically, in the present embodiment, for any two Target Aerial Array, from two respective subarrays of Target Aerial Array
It randomly selects a subarray and carries out combination of two, can get multiple subarrays thus directed towards two Target Aerial Arrays and combine, and
It include two subarrays in each subarray combination.For any one subarray combination in the combination of multiple subarrays, obtain
Two corresponding sub-array signals of subarray in subarray combination believe the corresponding subarray of two subarrays
It number carries out fourth order cumulant and calculates obtaining the subarray and combining corresponding fourth order cumulant matrix.Specific formula for calculation are as follows:
Wherein, x1aIt (t) is the corresponding sub-array signal of subarray a of Target Aerial Array 1;x2bIt (t) is the son of Target Aerial Array 2
The corresponding sub-array signal of array b;RabCorresponding fourth order cumulant square is combined for the subarray that subarray a and subarray b is constituted
Battle array;E indicates expectation.
It should be noted that the preset quantity of Target Aerial Array is 3 in the present embodiment, that is, there are 3 Target Aerial Arrays, respectively
For Target Aerial Array 1, Target Aerial Array 2 and Target Aerial Array 3;A, b and c is respectively in Target Aerial Array 1, Target Aerial Array 2 and Target Aerial Array 3
Any one subarray.Any one subarray a and target obtained in Target Aerial Array 1 can be calculated by above-mentioned calculation formula
The subarray that any one subarray b in array 2 is constituted combines corresponding fourth order cumulant Rab.Similarly, according to above-mentioned calculating
Formula can also calculate the composition of any one subarray c in any one subarray a and Target Aerial Array 3 in Target Aerial Array 1
Subarray combines corresponding fourth order cumulant RacAnd in any one subarray b and Target Aerial Array 3 in Target Aerial Array 2
The subarray that any one subarray c is constituted combines corresponding fourth order cumulant Rbc。
On this basis, it all subarrays combines corresponding fourth order cumulant matrix and carries out front and back to smoothing processing, i.e.,
It can get the corresponding smoothed out fourth order cumulant matrix of two Target Aerial Arrays.Specific formula for calculation are as follows:
Wherein, R12For Target Aerial Array 1 and the corresponding smoothed out fourth order cumulant of Target Aerial Array 2;A is Target Aerial Array 1
Any one subarray;B is any one subarray of Target Aerial Array 2;N1For the array number of Target Aerial Array 1;M1For Target Aerial Array
The array number of 1 each subarray;N2For the array number of Target Aerial Array 2;M2For the array number of each subarray of Target Aerial Array 2;
J is the permutation matrix that opposition angle array element element is 1.
It should be noted that since the array number of Target Aerial Array 1 is N1, and the array element of each subarray of Target Aerial Array 1
Number M1, then the submatrix number of columns of Target Aerial Array 1 is N1-M1+1;Therefore the value of a is (0, N in above-mentioned calculation formula1-M1).Due to
The array number of Target Aerial Array 2 is N2, and the array number M of each subarray of Target Aerial Array 22, then the submatrix columns of Target Aerial Array 2
Amount is N2-M2+1;Therefore the value of b is (0, N in above-mentioned calculation formula2-M2)。
Similarly, the corresponding smoothed out quadravalence of any other two Target Aerial Arrays can be obtained according to above method step to tire out
Accumulated amount matrix R13And R23.It can effective ring off signal phase by the smoothed out fourth order cumulant matrix that above method step obtains
It closes.
A kind of coherent signal Wave arrival direction estimating method based on mutual pixel array provided by the invention, for any two mesh
Array is marked, a subarray is chosen respectively from all subarrays of two Target Aerial Arrays and is combined, obtain multiple submatrixs
Column combination, combines any one subarray, according to the corresponding sub-array signal of all subarrays in subarray combination
It obtains the subarray and combines corresponding fourth order cumulant matrix;Corresponding fourth order cumulant matrix all subarrays is combined to carry out
Front and back obtains the corresponding smoothed out fourth order cumulant matrix of two Target Aerial Arrays to smoothing processing.This method passes through antithetical phrase
Array signal carries out fourth order cumulant and calculates the fourth order cumulant matrix for obtaining sub-array signal, then the quadravalence of sub-array signal
Cumulant matrices are smoothed, can effectively ring off signal it is related, be conducive to enhance the spatial resolution to signal, in turn
Be conducive to improve the accuracy of signal estimated result.
Based on any of the above-described embodiment, a kind of coherent signal Wave arrival direction estimating method based on mutual pixel array is provided, it is right
It is corresponding that each smoothed out fourth order cumulant matrix carries out each smoothed out fourth order cumulant matrix of Eigenvalues Decomposition acquisition
Projection operator, specifically: for any one smoothed out fourth order cumulant matrix, to the smoothed out fourth order cumulant matrix
Eigenvalues Decomposition is carried out, multiple feature vectors are obtained;It is screened from all feature vectors according to the characteristic value of all feature vectors
The feature vector for meeting preset condition out, as target feature vector;All target feature vectors are combined into objective matrix, root
The corresponding projection operator of the smoothed out fourth order cumulant matrix is obtained according to objective matrix.
Specifically, for any one smoothed out fourth order cumulant matrix, to the smoothed out fourth order cumulant matrix
Eigenvalues Decomposition is carried out, multiple feature vectors is obtained, is screened from all feature vectors according to the characteristic value of all feature vectors
The feature vector for meeting preset condition out, as target feature vector.In the present embodiment, by all feature vectors according to characteristic value
Order of magnitude arranged, biggish characteristic value and lesser characteristic value are isolated, wherein biggish characteristic value is generally
10 times or more of smaller characteristic value, will the corresponding feature vector of smaller characteristic value as target feature vector.Finally, will smaller spy
Value indicative character pair vector is combined into objective matrix, and it is corresponding to obtain the smoothed out fourth order cumulant matrix according to objective matrix
Projection operator.Similarly, step can get the corresponding projection operator of other smoothed out fourth order cumulant matrix according to the method described above.
For example, fourth order cumulant matrix R smoothed out for 3 obtained in above method embodiment12、R13With
R23, with R12For, to R12After carrying out feature decomposition, the objective matrix that the smaller characteristic value of acquisition is constituted is V12n, then corresponding
Projection operatorSimilarly, it can get R13And R23Corresponding projection operator T13And T23, wherein H representing matrix is total to
Yoke transposition.
A kind of coherent signal Wave arrival direction estimating method based on mutual pixel array provided by the invention is flat for any one
Fourth order cumulant matrix after cunning carries out Eigenvalues Decomposition to the smoothed out fourth order cumulant matrix, obtain multiple features to
Amount;The feature vector for meeting preset condition is filtered out from all feature vectors according to the characteristic value of all feature vectors, as
Target feature vector;All target feature vectors are combined into objective matrix, which is obtained according to objective matrix
The corresponding projection operator of cumulant matrices.This method is by carrying out Eigenvalues Decomposition acquisition to smoothed out fourth order cumulant matrix
Projection operator is conducive to the direction of arrival that signal is estimated according to projection operator.
Based on any of the above-described embodiment, a kind of coherent signal Wave arrival direction estimating method based on mutual pixel array is provided, benefit
The direction of arrival of signal is estimated according to all projection operators with Root-MUSIC algorithm, specifically: any one projection is calculated
Son constructs the corresponding polynomial function of the projection operator using preset function building method;It is corresponding according to all projection operators
Polynomial function establishes target equation, solves target equation, obtains the equattion root being distributed on unit circle;Estimated according to equattion root
The direction of arrival of signal.
Specifically, on the basis of above-mentioned acquisition multiple projection operators, for any one projection operator, default letter is utilized
Number building method constructs the corresponding polynomial function of the projection operator.The projection operator T in the above way obtained in embodiment12For
Example corresponds to the polynomial function f of construction12(z) are as follows:
Z=exp (j ω);
Wherein, j is imaginary unit, i.e. radical sign -1;ω is arbitrary number;M1For the array element of each subarray of Target Aerial Array 1
Number;M2For the array number of each subarray of Target Aerial Array 2.Similarly, it can get projection operator T13And T23Corresponding multinomial letter
Number f13(z) and f23(z)。
Further, target equation is established according to the corresponding polynomial function of all projection operators, solves target equation, obtains
The equattion root that must be distributed on unit circle.In above-mentioned acquisition polynomial function f12(z)、f13(z) and f23(z) it on the basis of, establishes
Target equation are as follows:
F (z)=f12(z)+f13(z)+f23(z)=0
Above-mentioned equation is solved, obtaining i-th of equattion root being distributed on unit circle is
Finally, estimate the direction of arrival of signal, according to equattion root with equattion rootFor, the signal wave finally estimated reaches side
To may be expressed as:
A kind of coherent signal Wave arrival direction estimating method based on mutual pixel array provided by the invention, for any one throwing
Shadow operator constructs the corresponding polynomial function of the projection operator using preset function building method;According to all projection operators pair
The polynomial function answered establishes target equation, solves target equation, obtains the equattion root being distributed on unit circle;According to equattion root
Estimate the direction of arrival of signal.This method is based on mutual prime matrix and fourth order cumulant obtains projection operator, further according to projection operator
The direction of arrival for estimating signal, the precision of estimated result can be effectively ensured in the case where not increasing array element number.
Fig. 2 is a kind of whole knot of coherent signal Mutual coupling system based on mutual pixel array of the embodiment of the present invention
Structure schematic diagram provides a kind of coherent signal direction of arrival based on mutual pixel array as shown in Fig. 2, being based on any of the above-described embodiment
Estimating system, comprising:
Signal acquisition module 1, for being sampled according to signal received by Target Aerial Array of the time domain to preset quantity,
The corresponding time domain discrete signal of each Target Aerial Array is obtained, wherein every two Target Aerial Array mutual pixel array each other;
Signal pre-processing module 2, for setting the array element of the subarray of the Target Aerial Array for any one Target Aerial Array
Number, all submatrixs of the Target Aerial Array are obtained according to the array number of subarray from the corresponding time domain discrete signal of the Target Aerial Array
Arrange corresponding sub-array signal;
Smoothing computation module 3 is used for for any two Target Aerial Array, according to all subarrays of two Target Aerial Arrays
Corresponding sub-array signal obtains the corresponding smoothed out fourth order cumulant matrix of two Target Aerial Arrays;
Eigenvalues Decomposition module 4 obtains often for carrying out Eigenvalues Decomposition to each smoothed out fourth order cumulant matrix
The corresponding projection operator of a smoothed out fourth order cumulant matrix;
Mutual coupling module 5, for estimating the wave of signal according to all projection operators using Root-MUSIC algorithm
Up to direction.
Specifically, the present invention provides a kind of coherent signal Mutual coupling system based on mutual pixel array, including signal
Acquisition module 1, signal pre-processing module 2, smoothing computation module 3, Eigenvalues Decomposition module 4 and Mutual coupling module 5 are led to
The cooperation crossed between each module realizes that the method in any of the above-described embodiment of the method, specific implementation process can be found in above method reality
Example is applied, details are not described herein again.
A kind of coherent signal Mutual coupling system based on mutual pixel array provided by the invention, connects mutual pixel array
The signal being subject to carries out sampling and obtains time domain discrete signal;Again by the array number of setting subarray, according to time domain discrete signal
The fourth order cumulant of sub-array signal is carried out front and back to being put down after smoothing processing by the fourth order cumulant for obtaining sub-array signal
Fourth order cumulant after cunning;Eigenvalues Decomposition is carried out to smoothed out fourth order cumulant and obtains projection operator, finally utilizes Root-
MUSIC algorithm estimates the direction of arrival of signal according to all projection operators.The system, which passes through, combines mutual pixel array and fourth-order cumulant
Amount carries out signal Mutual coupling, and the precision of estimated result can be effectively ensured in the case where not increasing array element number,
Signal number has preferable effect when being more than array element number;And can effectively ring off signal it is related, enhance the sky to signal
Between resolving power, can be suitable for various complex environments, can effectively meet various Practical Project demands.
Fig. 3 shows the structural block diagram of a kind of electronic equipment of the embodiment of the present invention.Referring to Fig. 3, the electronic equipment, packet
It includes: processor (processor) 31, memory (memory) 32 and bus 33;Wherein, the processor 31 and memory 32 are logical
It crosses the bus 33 and completes mutual communication;The processor 31 is used to call the program instruction in the memory 32, with
Execute method provided by above-mentioned each method embodiment, for example, received according to Target Aerial Array of the time domain to preset quantity
To signal sampled, obtain the corresponding time domain discrete signal of each Target Aerial Array, wherein every two Target Aerial Array is mutual each other
Pixel array;For any one Target Aerial Array, the array number of the subarray of the Target Aerial Array is set, according to the array number of subarray
The corresponding sub-array signal of all subarrays of the Target Aerial Array is obtained from the corresponding time domain discrete signal of the Target Aerial Array;It is right
In any two Target Aerial Array, which is obtained according to the corresponding sub-array signal of all subarrays of two Target Aerial Arrays
Mark the corresponding smoothed out fourth order cumulant matrix of array;Eigenvalues Decomposition is carried out to each smoothed out fourth order cumulant matrix
Obtain the corresponding projection operator of each smoothed out fourth order cumulant matrix;It is calculated using Root-MUSIC algorithm according to all projections
The direction of arrival of son estimation signal.
The present embodiment discloses a kind of computer program product, and the computer program product includes being stored in non-transient calculating
Computer program on machine readable storage medium storing program for executing, the computer program include program instruction, when described program instruction is calculated
When machine executes, computer is able to carry out method provided by above-mentioned each method embodiment, for example, according to time domain to present count
Signal received by the Target Aerial Array of amount is sampled, and the corresponding time domain discrete signal of each Target Aerial Array is obtained, wherein often
Two Target Aerial Arrays mutual pixel array each other;For any one Target Aerial Array, the array number of the subarray of the Target Aerial Array is set,
All subarrays of the Target Aerial Array are obtained from the corresponding time domain discrete signal of the Target Aerial Array according to the array number of subarray
Corresponding sub-array signal;For any two Target Aerial Array, according to the corresponding son of all subarrays of two Target Aerial Arrays
Array signal obtains the corresponding smoothed out fourth order cumulant matrix of two Target Aerial Arrays;To each smoothed out fourth-order cumulant
Moment matrix carries out Eigenvalues Decomposition and obtains the corresponding projection operator of each smoothed out fourth order cumulant matrix;Utilize Root-
MUSIC algorithm estimates the direction of arrival of signal according to all projection operators.
The present embodiment provides a kind of non-transient computer readable storage medium, the non-transient computer readable storage medium
Computer instruction is stored, the computer instruction makes the computer execute method provided by above-mentioned each method embodiment, example
As included: that signal received by the Target Aerial Array according to time domain to preset quantity samples, each Target Aerial Array pair is obtained
The time domain discrete signal answered, wherein every two Target Aerial Array mutual pixel array each other;For any one Target Aerial Array, the mesh is set
The array number for marking the subarray of array, obtains from the corresponding time domain discrete signal of the Target Aerial Array according to the array number of subarray
The corresponding sub-array signal of all subarrays of the Target Aerial Array;For any two Target Aerial Array, according to two target arrays
The corresponding sub-array signal of all subarrays of column obtains the corresponding smoothed out fourth order cumulant matrix of two Target Aerial Arrays;
It is corresponding that each smoothed out fourth order cumulant matrix of Eigenvalues Decomposition acquisition is carried out to each smoothed out fourth order cumulant matrix
Projection operator;The direction of arrival of signal is estimated according to all projection operators using Root-MUSIC algorithm.
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 be stored in a computer readable storage medium, the program
When being executed, step including the steps of the foregoing method embodiments is executed;And storage medium above-mentioned includes: ROM, RAM, magnetic disk or light
The various media that can store program code such as disk.
The embodiments such as electronic equipment described above are only schematical, wherein it is described as illustrated by the separation member
Unit may or may not be physically separated, and component shown as a unit may or may not be object
Manage unit, it can it is in one place, or may be distributed over multiple network units.It can select according to the actual needs
Some or all of the modules therein is selected to achieve the purpose of the solution of this embodiment.Those of ordinary skill in the art are not paying wound
In the case where the labour for the property made, it can understand and implement.
Through the above description of the embodiments, those skilled in the art can be understood that each embodiment can
It realizes by means of software and necessary general hardware platform, naturally it is also possible to pass through hardware.Based on this understanding, on
Stating technical solution, substantially the part that contributes to existing technology can be embodied in the form of software products in other words, should
Computer software product may be stored in a computer readable storage medium, such as ROM/RAM, magnetic disk, CD, including several fingers
It enables and using so that a computer equipment (can be personal computer, server or the network equipment etc.) executes each implementation
Method described in certain parts of example or embodiment.
Finally, the present processes are only preferable embodiment, it is not intended to limit the scope of the present invention.It is all
Within the spirit and principles in the present invention, any modification, equivalent replacement, improvement and so on should be included in protection of the invention
Within the scope of.
Claims (10)
1. a kind of coherent signal Wave arrival direction estimating method based on mutual pixel array characterized by comprising
It is sampled according to signal received by Target Aerial Array of the time domain to preset quantity, it is corresponding to obtain each Target Aerial Array
Time domain discrete signal, wherein every two Target Aerial Array mutual pixel array each other;
For any one Target Aerial Array, set the array number of the subarray of the Target Aerial Array, according to the array number of subarray from
The corresponding sub-array signal of all subarrays of the Target Aerial Array is obtained in the corresponding time domain discrete signal of the Target Aerial Array;
For any two Target Aerial Array, it is somebody's turn to do according to the corresponding sub-array signal of all subarrays of two Target Aerial Arrays
The corresponding smoothed out fourth order cumulant matrix of two Target Aerial Arrays;
Eigenvalues Decomposition is carried out to each smoothed out fourth order cumulant matrix and obtains each smoothed out fourth order cumulant matrix
Corresponding projection operator;
The direction of arrival of signal is estimated according to all projection operators using Root-MUSIC algorithm.
2. the method according to claim 1, wherein corresponding from the Target Aerial Array according to the array number of subarray
The corresponding sub-array signal of all subarrays of the Target Aerial Array is obtained in time domain discrete signal, specifically:
For any one Target Aerial Array, using the array number of the corresponding subarray of the Target Aerial Array as destination number, in the mesh
In all array elements for marking array, the adjacent array element of destination number is combined, all subarrays of the Target Aerial Array are obtained;
The corresponding sub-array signal of each subarray is isolated from the corresponding time domain discrete signal of the Target Aerial Array, obtains the mesh
Mark the corresponding sub-array signal of all subarrays of array.
3. the method according to claim 1, wherein corresponding according to all subarrays of two Target Aerial Arrays
Sub-array signal obtains the corresponding smoothed out fourth order cumulant matrix of two Target Aerial Arrays, specifically:
For any two Target Aerial Array, a subarray is randomly selected from two respective subarrays of Target Aerial Array and is carried out
Combination of two obtains multiple subarray combinations, any one subarray is combined, according to all sons in subarray combination
The corresponding sub-array signal of array obtains the subarray and combines corresponding fourth order cumulant matrix;
It all subarrays combines corresponding fourth order cumulant matrix and carries out front and back to smoothing processing, obtain two Target Aerial Arrays
Corresponding smoothed out fourth order cumulant matrix.
4. according to the method described in claim 3, it is characterized in that, corresponding according to all subarrays in subarray combination
Sub-array signal obtains the subarray and combines corresponding fourth order cumulant matrix, specific formula for calculation are as follows:
Wherein, x1aIt (t) is the corresponding sub-array signal of subarray a of Target Aerial Array 1;x2bIt (t) is the subarray b of Target Aerial Array 2
Corresponding sub-array signal;RabCorresponding fourth order cumulant matrix is combined for the subarray that subarray a and subarray b is constituted.
5. according to the method described in claim 3, it is characterized in that, all subarrays combining corresponding fourth order cumulant matrix
Front and back is carried out to smoothing processing, obtains the corresponding smoothed out fourth order cumulant matrix of two Target Aerial Arrays, it is specific to calculate public affairs
Formula are as follows:
Wherein, R12For Target Aerial Array 1 and the corresponding smoothed out fourth order cumulant of Target Aerial Array 2;A is any of Target Aerial Array 1
One subarray;B is any one subarray of Target Aerial Array 2;N1For the array number of Target Aerial Array 1;M1For Target Aerial Array 1
The array number of each subarray;N2For the array number of Target Aerial Array 2;M2For the array number of each subarray of Target Aerial Array 2;J is
Oppose that array element element in angle is 1 permutation matrix.
6. the method according to claim 1, wherein carrying out feature to each smoothed out fourth order cumulant matrix
Value, which is decomposed, obtains the corresponding projection operator of each smoothed out fourth order cumulant matrix, specifically:
For any one smoothed out fourth order cumulant matrix, characteristic value point is carried out to the smoothed out fourth order cumulant matrix
Solution, obtains multiple feature vectors;
The feature vector for meeting preset condition is filtered out from all feature vectors according to the characteristic value of all feature vectors, as
Target feature vector;
All target feature vectors are combined into objective matrix, which is obtained according to objective matrix
Corresponding projection operator.
7. the method according to claim 1, wherein being estimated using Root-MUSIC algorithm according to all projection operators
The direction of arrival of signal is counted, specifically:
For any one projection operator, the corresponding polynomial function of the projection operator is constructed using preset function building method;
Target equation is established according to the corresponding polynomial function of all projection operators, solves target equation, acquisition is distributed in unit
Equattion root on circle;
The direction of arrival of signal is estimated according to equattion root.
8. a kind of coherent signal Mutual coupling system based on mutual pixel array characterized by comprising
Signal acquisition module is obtained for sampling according to signal received by Target Aerial Array of the time domain to preset quantity
The corresponding time domain discrete signal of each Target Aerial Array, wherein every two Target Aerial Array mutual pixel array each other;
Signal pre-processing module, for setting the array number of the subarray of the Target Aerial Array, root for any one Target Aerial Array
All subarrays pair of the Target Aerial Array are obtained from the corresponding time domain discrete signal of the Target Aerial Array according to the array number of subarray
The sub-array signal answered;
Smoothing computation module, for being corresponded to according to all subarrays of two Target Aerial Arrays for any two Target Aerial Array
Sub-array signal obtain the corresponding smoothed out fourth order cumulant matrix of two Target Aerial Arrays;
Eigenvalues Decomposition module, it is each smooth for carrying out Eigenvalues Decomposition acquisition to each smoothed out fourth order cumulant matrix
The corresponding projection operator of fourth order cumulant matrix afterwards;
Mutual coupling module, for estimating that the wave of signal reaches side according to all projection operators using Root-MUSIC algorithm
To.
9. a kind of electronic equipment characterized by comprising
At least one processor;And
At least one processor being connect with the processor communication, in which:
The memory is stored with the program instruction that can be executed by the processor, and the processor calls described program to instruct energy
Enough methods executed as described in claim 1 to 7 is any.
10. a kind of non-transient computer readable storage medium, which is characterized in that the non-transient computer readable storage medium is deposited
Computer instruction is stored up, the computer instruction makes the computer execute the method as described in claim 1 to 7 is any.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810847429.2A CN109143155B (en) | 2018-07-27 | 2018-07-27 | Correlated signal direction of arrival estimation method and system based on mutual prime array |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810847429.2A CN109143155B (en) | 2018-07-27 | 2018-07-27 | Correlated signal direction of arrival estimation method and system based on mutual prime array |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109143155A true CN109143155A (en) | 2019-01-04 |
CN109143155B CN109143155B (en) | 2020-06-02 |
Family
ID=64799238
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810847429.2A Active CN109143155B (en) | 2018-07-27 | 2018-07-27 | Correlated signal direction of arrival estimation method and system based on mutual prime array |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109143155B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110208736A (en) * | 2019-07-11 | 2019-09-06 | 西安电子科技大学 | Non-circular signal uniform array direction of arrival angle estimation method based on fourth-order cumulant |
CN110531312A (en) * | 2019-08-29 | 2019-12-03 | 深圳市远翰科技有限公司 | A kind of DOA estimation method and system based on sparse symmetric matrix column |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105182293A (en) * | 2015-08-25 | 2015-12-23 | 西安电子科技大学 | Method for estimating DOA and DOD of MIMO radar based on co-prime array |
CN106226729A (en) * | 2016-07-15 | 2016-12-14 | 西安电子科技大学 | Relatively prime array direction of arrival angular estimation method based on fourth-order cumulant |
CN106443574A (en) * | 2016-11-08 | 2017-02-22 | 西安电子科技大学 | Direction-of-arrival (DOA) angle estimation method based on dual-layer nested array |
CN107015190A (en) * | 2017-03-01 | 2017-08-04 | 浙江大学 | Relatively prime array Wave arrival direction estimating method based on the sparse reconstruction of virtual array covariance matrix |
CN107290709A (en) * | 2017-05-05 | 2017-10-24 | 浙江大学 | The relatively prime array Wave arrival direction estimating method decomposed based on vandermonde |
CN108267712A (en) * | 2018-01-02 | 2018-07-10 | 天津大学 | A kind of DOA estimation method and device that mutual pixel array is translated based on compression |
-
2018
- 2018-07-27 CN CN201810847429.2A patent/CN109143155B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105182293A (en) * | 2015-08-25 | 2015-12-23 | 西安电子科技大学 | Method for estimating DOA and DOD of MIMO radar based on co-prime array |
CN106226729A (en) * | 2016-07-15 | 2016-12-14 | 西安电子科技大学 | Relatively prime array direction of arrival angular estimation method based on fourth-order cumulant |
CN106443574A (en) * | 2016-11-08 | 2017-02-22 | 西安电子科技大学 | Direction-of-arrival (DOA) angle estimation method based on dual-layer nested array |
CN107015190A (en) * | 2017-03-01 | 2017-08-04 | 浙江大学 | Relatively prime array Wave arrival direction estimating method based on the sparse reconstruction of virtual array covariance matrix |
CN107290709A (en) * | 2017-05-05 | 2017-10-24 | 浙江大学 | The relatively prime array Wave arrival direction estimating method decomposed based on vandermonde |
CN108267712A (en) * | 2018-01-02 | 2018-07-10 | 天津大学 | A kind of DOA estimation method and device that mutual pixel array is translated based on compression |
Non-Patent Citations (1)
Title |
---|
YANG HU等: "DOA Estimation of Coherent Signals on Coprime Arrays Exploiting Fourth-Order Cumulants", 《SENSORS》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110208736A (en) * | 2019-07-11 | 2019-09-06 | 西安电子科技大学 | Non-circular signal uniform array direction of arrival angle estimation method based on fourth-order cumulant |
CN110208736B (en) * | 2019-07-11 | 2023-02-10 | 西安电子科技大学 | Non-circular signal uniform array direction-of-arrival angle estimation method based on fourth-order cumulant |
CN110531312A (en) * | 2019-08-29 | 2019-12-03 | 深圳市远翰科技有限公司 | A kind of DOA estimation method and system based on sparse symmetric matrix column |
CN110531312B (en) * | 2019-08-29 | 2021-09-17 | 深圳市远翰科技有限公司 | DOA estimation method and system based on sparse symmetric array |
Also Published As
Publication number | Publication date |
---|---|
CN109143155B (en) | 2020-06-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108872929B (en) | Estimation method for direction of arrival of co-prime array based on rotation invariance of covariance matrix subspace of interpolated virtual array | |
Rangarao et al. | gold-MUSIC: A Variation on MUSIC to Accurately Determine Peaks of the Spectrum | |
CN107561484B (en) | Direction-of-arrival estimation method based on interpolation co-prime array covariance matrix reconstruction | |
CN110463147A (en) | The method of solution code sign and reception and the receiver for solving code sign | |
CN108020812B (en) | Two-dimensional DOA estimation method based on special three-parallel line array structure | |
CN107870315B (en) | Method for estimating direction of arrival of any array by using iterative phase compensation technology | |
CN107037392A (en) | A kind of relatively prime array Wave arrival direction estimating method of free degree increase type based on compressed sensing | |
CN107589399A (en) | Based on the relatively prime array Wave arrival direction estimating method for sampling virtual signal singular values decomposition more | |
CN108802669B (en) | Two-dimensional direction of arrival estimation method, two-dimensional direction of arrival estimation device and terminal | |
CN109143155A (en) | Coherent signal Wave arrival direction estimating method and system based on mutual pixel array | |
CN109245814A (en) | Adaptive beamformer method based on maximum likelihood resampling | |
KR101958337B1 (en) | The method and apparatus for estimating the direction of arrival of a signal | |
CN110244273A (en) | It is a kind of based on the target angle estimation method for being uniformly distributed formula array | |
CN106680779B (en) | Beam-forming method and device under impulsive noise | |
Yuri et al. | Probability of false peaks occurring via circular and concentric antenna arrays DOA estimation | |
CN112731280B (en) | ESPRIT-DOA estimation method in inter-mass array mixed noise environment | |
CN108614234B (en) | Direction-of-arrival estimation method based on multi-sampling snapshot co-prime array received signal fast Fourier inverse transformation | |
Magdy et al. | Direction of arrival estimation based on maximum likelihood criteria using gravitational search algorithm | |
CN113314138A (en) | Sound source monitoring and separating method and device based on microphone array and storage medium | |
CN110245706B (en) | Lightweight target detection method for embedded application | |
CN110427669B (en) | Neural network model calculation method for phased array scanning radiation beams | |
CN108872930B (en) | Extended aperture two-dimensional joint diagonalization DOA estimation method | |
CN110208736A (en) | Non-circular signal uniform array direction of arrival angle estimation method based on fourth-order cumulant | |
CN114185002B (en) | Three-dimensional parameter estimation method based on beam space matrix beam | |
Al-Sadoon et al. | New and less complex approach to estimate angles of arrival |
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 |