The interelement spacing of a coprime array breaks through the half-wavelength limitation, so that a larger array aperture can be obtained with a fixed number of array elements. In this paper, the symmetry of the noncircular signal is used to virtually expand the L-shaped array into an orthogonal cross array. Furthermore, the virtual received signal of the augmented array is obtained by the second-order statistic of the received data. Decoupling and dimension reduction of elevation and azimuth are realized by a z-axis subarray and y-axis subarray. Finally, the sparse reconstruction of the signal is realized by the minimum absolute convergence and selection operator method. This method can enlarge the aperture and freedom of array and has higher accuracy and resolution of DOA estimation. It has the advantages of automatic parameter pairing without additional pairing operation and is effective for coherent and incoherent signals. The final numerical simulation results prove the effectiveness of the method in this paper.
National Natural Science Foundation of China6177239861972239Key Research and Development Program Projects of Shaanxi Province2019SF-2572020GY-0242021GY-182Education Department of Shaanxi Province19JK0169Science and Technology Program of Hantai District2019KX-211. Introduction
Direction of arrival (DOA) estimation is an important research content of array signal processing, which is widely used in navigation, communication, and radar [1, 2]. After decades of development, a complete set of processing theory and method for the estimation of DOA has been formed. Subspace-like high-resolution algorithms such as MUSIC and ESPRIT break through the restriction of Rayleigh limits and have good high-resolution performance [3, 4]. In the actual environment, due to the existence of multipath effect and electromagnetic interference, the signals received by the array are often coherent signals. The existence of coherent signals leads to the performance degradation or even complete ineffectiveness of the traditional subspace-like high-resolution algorithms. The existence of the coherent signal results in the rank loss of the data covariance matrix. The rank of the data covariance matrix must be restored before the subspace algorithm can be applied to the DOA estimation of the coherent signal, which also leads to the increase of the calculation amount.
The application of a coprime array makes the degree of freedom of DOA estimation no longer limited by the number of physical array elements, thus greatly reducing the hardware overhead [5, 6]. The interelement spacing of the coprime array based on the sparse model is larger than half a wavelength, so the aperture of the array is effectively expanded and the resolution is improved. In recent years, the coprime array has become a research hotspot in the field of array signal processing. In [7, 8], a virtual array model based coprime array is presented [9]. In [10], a spatially smooth DOA estimation algorithm based on the continuous virtual array element is presented. However, this method leads to the abandonment of partial array elements due to the need of continuous virtual arrays. In order to compensate for this deficiency, a DOA algorithm based on the covariance matrix reconstruction of the coprime array is proposed in [11]. In [12], a two-dimensional angle estimation based on the L-shaped coprime array is proposed. In [13], the azimuth of the signals’ source is estimated by constructing the cumulant matrix, which reduces the aperture loss of the array. In [14], a two-dimensional direction-finding algorithm for sparsely L-shaped coprime arrays is represented.
Because of its unique statistical characteristics, the noncircular signals can be processed by using both the covariance matrix and the elliptic covariance matrix of the received data at the receiving end, which can effectively improve the performance of signal detection and estimation. With the wide application of noncircular signals in communication and other fields, such as common BPSK and AM signals, noncircular signal direction finding has become a new hot spot of research [15, 16].
The compressed sensing technology can achieve high-efficiency signal processing under the condition of undersampling by compressing redundant information. The compressed sensing algorithm does not require the eigen-decomposition of the data correlation matrix, so it can directly distinguish the coherent signals without additional decoherence processing. The degree of freedom of the nonuniform array can be fully utilized, unlike the spatial smoothing method and the Toeplitz method, which utilize only half of the degrees of freedom [17, 18].
In this paper, the advantages of the sparse coprime array and the characteristics of the noncircular signal are combined to receive the data from the actual array to obtain the received data of the symmetric virtual array, and the actual and virtual array data constitute the whole array data, so as to realize the further expansion of the array data. By vectorizing the covariance matrix of the whole array data and by deduplication, a new array receiving data is obtained. Combined with the structural advantages of the L-shaped array, the two-dimensional DOA estimation of the incoherent noncircular signal is realized by the compressed sensing algorithm, which enlarges the aperture of the array and the degree of freedom of the array to a certain extent and improves the resolution and estimation accuracy of the DOA.
2. Array Signal Model2.1. Array Architecture
As shown in Figure 1, the coprime array is a kind of nonuniform array, which consists of two subarrays.
The structure of the coprime array.
Here, the number of elements of the two subarrays satisfies the relatively prime relationship, and the interelement spacing also satisfies the relatively prime relationship, and the interelement spacing is no longer limited to half wavelength, so it is a sparse array. The structure of the coprime array and its subarrays is shown in Figure 2.
The structure of the coprime array and two subarrays.
By contrast, as shown in Figures 3 and 4, the array element position of the symmetric uniform array and symmetric coprime array is 0,d,2d,3d,…,M−1d and0,Md,Nd,2Md,2Nd,…,M−1Nd,N−1Md. The received data of the negative half-axis are also obtained by virtual expansion based on the symmetry of the circular signal, and virtual expansion is carried out by vectorization of the data correlation matrix. The array element position of the negative half-axis of the uniform array and coprime array is −d,−2d,−3d,…,−M−1d and −d,−Nd,−2Md,−2Nd,…,−M−1Nd,−N−1Md.
The structure of the symmetric uniform array, where the black array elements are real elements and the blue array elements are the corresponding virtual array elements.
The structure of the symmetric coprime array, where the white and black array elements are the real elements of the two subarrays and the blue and purple array elements are the virtual array elements corresponding to the real array elements.
2.2. Receiving Data Model
The receiving array is an L-shaped array composed of a coprime subarray Za distributed in the z-axis and a coprime subarray Ya distributed in the y-axis. The coprime array is composed of two sparse uniform arrays with the first element overlapped; one array is N array elements with interelement spacing Md, and the other array is M array elements with interelement spacing Nd. M and N are two mutually prime numbers. d is a fixed length, which is less than or equal to the minimum half wavelength. The schematic diagram of the array structure is shown in Figure 5.
The schematic diagram of the L-shaped array structure.
The array element numbers of the subarray Za and Ya are both Q = N + M−1, and the corresponding array element coordinates are, respectively, z˜=z˜1,z˜2,…,z˜q,…,z˜Qd and y˜=y˜1,y˜2,…,y˜q,…,y˜Qd, where z˜q and y˜q are integers, 1≤q≤Q.
K far-field, coherent, narrowband noncircular signals θk,ϕk are incident on the receiving array from different directions. θk∈0,π is the elevation angle of the kth signal, that is, the angle between the incident signal and the positive direction of the z-axis, and ϕk∈−π/2,π/2 is the azimuth angle of the kth incident signal, that is, the angle between the projection of the incident signal on the xoy plane and the positive direction of the x-axis. The received signal of the z-axis subarray Za is(1)Z1t=z1t,…,zqt,…,zQtT=AzSt+Nzt,where zqt=∑k=1Ksktejμkz˜q+nzqt is the received signal of the q-th array element on the z-axis, St=s1t,…,skt,…,sKtT is the incident signal matrix, Nzt is the noise of the z-axis subarray, and Az=azθ1,…,azθk,…,azθK is the array steering vector of the z-axis subarray. azθk=ejμkz˜1,ejμkz˜2,…,ejμkz˜q,…,ejμkz˜QT is the array steering vector corresponding to the kth signal, nzqt is the real array noise received by the q-th array element on the z-axis subarray, μk=−2πdcosθk/λ, λ is the wavelength of the incident signal, and skt is the kth incident signal, 1≤k≤K.
The received signal of the y-axis subarray Ya is(2)Y1t=y1t,…,yqt,…,yQtT=AySt+Nyt,where yqt=∑k=1Ksktejφky˜q+nyqt is the received signal of the q-th array element on the y-axis,Ay=ayθ1,ϕ1,…,ayθk,ϕk,…,ayθK,ϕK is the array steering vector of the y-axis subarray, and Nyt is the real array noise of the y-axis subarray. ayθk,ϕk=ejφky˜1,ejφky˜2,…,ejφky˜q,…,ejφky˜QT is the y-axis subarray steering vector corresponding to the kth signal, nyqt is noise received by the q-th array element on the y-axis subarray, and φk=−2πdsinθksinϕk/λ.
2.3. Virtual Expansion Array
The array elements of the original array are arranged in the positive half-axis of the coordinate axis. Using the characteristics of the noncircular signal, the received data of the negative half-axis is obtained by taking the conjugate operation to the received data of the positive half-axis array element.
In Figure 6, the array elements of the virtual subarrays Za′ and Ya′ are represented by hollow circles, and those of the actual subarrays Za and Ya are represented by solid circles.
The schematic diagram after virtual expansion.
The noncircular signal is actually the real signal, which satisfies the following relationship:(3)St=S∗t,where ·∗ denotes the conjugate operation.
The received signal Z1t of the z-axis subarray Za and the received signal Y1t of the y-axis subarray Ya are, respectively, conjugated to obtain the received data of the virtual array:(4)Z2t=zQ∗t,zQ−1∗t,…,zq∗t,…,z1∗tT,Y2t=yQ∗t,yQ−1∗t,…,yq∗t,…,y1∗tT,where ·T denotes the transpose of the matrix, zq∗t=z−qt, and yq∗t=y−qt.
The real and virtual array received signal construct the whole array received data:(5)Zt=Z2tZ1t=z−Qt,…,z−1t,z1t,…,zQT=A˜zSt+N˜zt,Yt=Y2tY1ty−Qt,…,y−1t,y1t,…,yQT=A˜ySt+N˜yt.
The steering vector in z-axis direction A˜z=a˜zθ1,…,a˜zθk,…,a˜zθK is only related to the elevation angle, the steering vector in the y-axis direction A˜y=a˜yθ1,ϕ1,…,a˜yθk,ϕk,…,a˜yθK,ϕK is related to the elevation and azimuth angles, N˜zt and N˜yt are, respectively, and the whole array noises of y-axis and z-axis subarray, where a˜zθk=e−jμkz˜Q…,e−jμkz˜q,…,e−jμkz˜1,ejμkz˜1,…,ejμkz˜q,…,ejμkz˜QT and a˜yθk,ϕk=e−jφky˜Q…,e−jφky˜q,…,ejφky˜2,e−jφky˜1,ejφky˜1,ejφky˜2,…,ejφky˜q,…,ejφky˜QT; then, the position coordinate of the array elements of the whole array is as follows:(6)z˜a=y˜a=−NM,−M−1N,…,−N,−M,M,N,…M−1N,NMd,where M and N are two coprime numbers, M < N, and d is less than or equal to the half wavelength of the incident signal.
3. Virtual Extended DOA Estimate Based on Vectorizaiton of the Covariance Matrix
The covariance matrix of data received by the whole array is(7)Rzz=EZZH=A˜zRsA˜zH+σnz2I2Q,Ryy=EYYH=A˜yRsA˜yH+σny2I2Q,where σnz2 and σny2 denote the white noise power of z-axis and y-axis subarrays, respectively, I2Q represents the identity matrix of 2Q×2Q, Rs=EStSHt is the signal correlation matrix, and ·H denotes the complex conjugate transpose.
The data covariance matrix Zv and Yv are vectorized to obtain the equivalent received data of the array under the augmented array, that is,(8)Zv=vecRzz=A˜z∗⊙A˜zP+σn2I4Q2=BzP+σnz2I4Q2,Yv=vecRyy=A˜y∗⊙A˜yP+σn2I4Q2=ByP+σny2I4Q2,where Zv and Yv are the column vectors of 4Q2×1, the rank of the two matrix is 1, vec· denotes the vectorization of the matrix, ⊙ represents the Khatri–Rao product, I4Q2 represents the identity column vector of 4Q2×1, and P=vecRs.(9)Bz=A˜z∗⊙A˜z=az∗θ1⊗azθ1,⋯,az∗θK⊗azθK,(10)By=A˜y∗⊙A˜y=ay∗θ1,ϕ1⊗ayθ1,ϕ1,…,ay∗θK,ϕK⊗ayθK,ϕK,where ⊗ represents the Kronecker product. Bz and By are often not the steering vectors of the uniform array. Using the method of transforming the coprime array into the uniform array, ZV and YV are transformed into the data Z˜V and Y˜V of the uniform augmented array, and the steering vectors Bz and By are transformed into the steering vectors B˜z and B˜y of the uniform augmented array. If the MUSIC method or ESPRIT method is used, the rank of the matrix needs to be restored by means of spatial smoothing and other methods, and the array will have half of the degree-of-freedom loss. Considering that the compressed sensing method does not require the decomposition of the data covariance matrix and is not affected by rank loss, the compressed sensing method is used to estimate the DOA of the signal.
3.1. Estimation of Elevation Angle
The estimation of elevation angle is obtained by using the data received from coprime subarrays arranged along the z-axis. From equation (9), it can be seen that the steering vector of the z-axis subarray is a single-valued function of the elevation angle. The overcomplete elevation angle sparse dictionary B˜zθ¯∈CF×Nθ is constructed according to the position of the array elements of the augmented array, where θ¯ is the potential signal of elevation angle and F is the number of elements of the uniform augmented array. The equivalent received data Z˜v under the z-axis augmented array can be expressed as(11)Z˜v=B˜zθ¯x′+σnz2IF,where x′∈CNθ×1 represents the variance of the potential signals, x′ has a sparse structure and K nonzero lines, each nonzero position corresponds to the DOA of the sound source signal, Nθ is the number of the potential elevation angle signal, Nθ≫K, Nθ≫F, and IF represents the identity matrix of F×F.
The 1 norm x^=argminx′x′1 satisfying the constraintZ˜v−B˜zθ¯x′−σnz2IF2≤ε1 obtains x′. x′ is a column vector, each of its elements corresponds to a column in dictionary B˜zθ¯, a column in dictionary corresponds to an incidence angle, and x′ has K nonzero elements, which correspond to K incidence angles. With the consideration of noise, the elevation angle θ^k can be solved by the minimum absolute convergence method and selection operator method, and the optimization problem of formula (12) can be solved by the OMP optimization package in MATLAB.(12)x^=argminx′x′1s.t.Z˜v−B˜zθ¯x′−σnz2IF2≤ε1,where ·1 and ·2 denote the 1-norm and 2-norm of matrices and ε1 is the noise constraint, and cross validation is used to select the appropriate value.
3.2. Estimation of Azimuth Angle
The estimates of azimuth angle can be achieved by using the data received from coprime subarrays arranged along the y-axis. It can be seen from formula (10) that the steering vector of the y-axis subarray is a function of elevation and azimuth angle. Substituting the estimated value of the elevation angle, the overcomplete azimuth angle sparse dictionary B˜yϕ¯∈CF×Nϕ is constructed according to the array element position of the y-axis augmented array, where ϕ¯ is the potential signal of azimuth angle. The equivalent received data Y˜v under the z-axis augmented array can be represented as(13)Y˜v=B˜yϕ¯y′+σny2IF,where y′∈CNϕ×1 indicates the variance of the potential signals, y′ has a sparse structure and K nonzero lines, and each nonzero position corresponds to the DOA of the sound signal. Nϕ is the number of potential azimuth angle, Nϕ≫K, Nϕ≫F.
In the case of considering noise, the azimuth angle can be estimated by solving the ℓ1 norm of equation (14); the optimization procedure of formula (14) is similar to that of formula (12):(14)y^=argminy′y′1s.t.Y˜v−B˜yϕ¯y′−σny2IF2≤ε2,where ε2 is the noise constraint, and cross validation is used to select the appropriate value.
The parameter estimation algorithm of compressed sensing dimension reduction based on the L-shaped array is proposed. Thanks to the proposed B˜zθ¯ in (12) and B˜yϕ¯ in (14) being one-dimensional dictionaries, we achieve decoupling of the elevation and azimuth. The automatic matching of azimuth and elevation angles does not need an additional pairing operation, which reduces the amount of calculation. If MUSIC or ESPRIT algorithm is used because Z˜v in (11) and Y˜v in (13) are column vectors with rank 1, it is necessary to reconstruct the data matrix with the idea of spatial smoothing to restore the rank of the matrix, resulting in half loss of the array aperture. In this paper, the optimization equations of formulas (12) and (14) are used to solve the problem by using compressed sensing, which does not need signal subspace, does not need to restore the rank of the matrix, and does not have the loss of the array aperture. The compressed sensing algorithm can process coherent and incoherent signals.
4. Simulation and Performance Analysis
The receiving array is an L-shaped array composed of a coprime subarray Za distributed in the z-axis and a coprime subarray Ya distributed in the y-axis. The parameters of two coprime arrays take M = 2, N = 3, and the positions of the array elements are (0, 2d, 3d, 4d), where d=λ/2, λ is the incident signal wavelength. Considering that two far-field narrowband signals are incident on the array, the signals are independent of each other. The incident angles of the two signals are, respectively, 25°,20° and 70°,50°, and the noise of the array is Gaussian white noise. The number of Monte Carlo experiments corresponding to each SNR is 50 times. For the sake of convenience, the array structure used in this paper is called the symmetric coprime array, and the corresponding DOA estimation method is called the proposed method.
Simulation 1.
The elevation and azimuth angles estimated by the proposed method are simulated. The SNR is 10 dB. The number of snapshots is 200, that is to say, all array sensors are sampled 200 times synchronously, and 200 synchronous sampling data are, herein, obtained. The more the snapshots, the more the information accumulated and the more accurate the estimation results. Especially for the subspace method, the number of snapshots is too small, the performance is degraded, and even the algorithm is invalid.
It can be seen from Figure 7 that the estimated values are in good agreement with the actual values, and the estimated disturbance is nearly within the 0.5-degree range of the real values, which shows that the proposed method can correctly estimate the elevation angle and azimuth angles of the incident signals.
The DOA estimation versus the number of experiments.
Simulation 2.
The variation of RMSE of elevation and azimuth with SNR is studied. The snapshot number is set as 200, and the SNR varies from 0 dB to 40 dB. The array element position of the symmetric coprime array is (0, 2d, 3d, 4d). By contrast, the array element position of the symmetric uniform array is set as (0, d, 2d, 3d), the received data of negative half-axis (−d, −2d, −3d) is also obtained by virtual expansion based on the symmetry of circular signal, and virtual expansion is carried out by vectorization of the data correlation matrix. Furthermore, the array element position of the coprime array is set as (0, 2d, 3d, 4d), and it does not make use of circular signal for virtual expansion or vectorization of the data correlation matrix for virtual expansion. The following is a comparison of the corresponding estimation methods of these three arrays.
Figures 8 and 9 are the performance analysis of the DOA estimation, estimated by the symmetric uniform array method, coprime array method, and the proposed method. The RMSE of the three array methods is compared with the change of the SNR. It can be seen that the parameter estimation performance of the symmetric coprime array method and the symmetric uniform array method has been improved significantly; especially, in the low SNR, the estimation accuracy of the symmetric coprime array method is significantly higher than that of the method of the symmetric uniform array and coprime array. The RMSE of DOA estimation of the coprime array does not decrease with the increase of SNR, which is biased estimation.
The simulation results also show that there is little difference in the calculation time of the three methods. The average time of the symmetric uniform array is 2.346 s, coprime array is 2.458 s, and symmetric coprime array is 2.509 s. DOA estimation-time increase of the coprime and symmetric coprime array is mainly caused by vectorization and deredundancy.
The RMSE of elevation versus SNR.
The RMSE of azimuth versus SNR.
Simulation 3.
The variation of RMSE of DOA with the number of snapshots is studied, the SNR is set as 10 dB, and the range of the snapshots is from 100 to 1000. It can be seen from Figures 10 and 11 that when the number of snapshots reaches 200, the estimated performance almost no longer changes with the increase of the number of snapshots. The estimation error of the coprime array does not decrease with the increase of the number of snapshots, which is biased estimation.
The RMSE of elevation versus snapshot number.
The RMSE of azimuth versus snapshot number.
Simulation 4.
The variation of success probability of DOA (elevation and azimuth angles) with SNR is studied, the snapshot number is set as 200, and the SNR varies from 0 dB to 40 dB.
From Figures 12 and 13, it can be seen that the DOA estimation performance of the proposed symmetric coprime array method is better than that of the other two array methods. This is because the coprime symmetric array adopted in this paper effectively expands the array aperture and increases the resolution under the same number of array elements.
The probability of success of elevation.
The probability of success of azimuth.
5. Conclusions
In this paper, a two-dimensional DOA estimation method of the coprime array is studied. An L-shaped coprime array arranged in the z-axis and y-axis is used as a receiving array. The second-order statistic of the coprime array is used to construct the virtual augmented array to receive signals. Two-dimensional DOA estimation is realized by the sparse reconstruction method and optimal solution by minimum absolute convergence and selection operator. The decoupling of elevation and azimuth angles is realized by using the advantage that the z-axis subarray is only related to elevation angle. The two-dimensional parameter estimation is realized by two one-dimensional sparse reconstructions, and the parameters are automatically paired, and no additional pairing operation is needed. It can be seen that if the coprime array does not use the vectorization of the covariance matrix for the virtual expansion of the array, there will be ambiguity in the DOA estimation, resulting in the biased estimation. However, the proposed method greatly improves the freedom of the array on the premise of maintaining the same number of elements as that of the two other array methods. The proposed method also improves the accuracy and resolution of DOA estimation under low SNR.
Appendix
Due to the great amount of different variables and notations in this manuscript, some points are hard to be followed by the reader, so Table 1 is introduced to facilitate the reader to read.
The explanation of variables and notations.
Variable
Meaning
Variable
Meaning
Za
The coprime subarray distributed in the z-axis
Za′
The virtual subarrays distributed in the z-axis
Ya
The coprime subarray distributed in the y-axis
Ya′
The virtual subarrays distributed in the y-axis
Y1t
The received signal of the y-axis real subarray
Z1t
The received signal of the z-axis real subarray
Y2t
The received signal of the y-axis virtual subarray
Z2t
The received signal of the z-axis virtual subarray
Yt
The whole array received signal of the y-axis subarray
Zt
The whole array received signal of the y-axis subarray
Ay
The y-axis real subarray steering vector
Az
The z-axis real subarray steering vector
A˜y
The whole array steering vector of the y-axis subarray
A˜z
The whole array steering vector of the z-axis subarray
Nyt
The real y-axis subarray noise
Nzt
The noise of the z-axis subarray
N˜yt
The whole array noises of the y-axis subarray
N˜zt
The whole array noises of the z-axis subarray
Ryy
The covariance matrix of the whole received data of the y-axis subarray
Rzz
The covariance matrix of the whole received data of the y-axis subarray
Yv
The vectorized data covariance matrix of the y-axis subarray
Zv
The vectorized data covariance matrix of the y-axis subarray
Y˜v
The equivalent received data under the z-axis augmented array
Z˜v
The equivalent received data under the z-axis augmented array
By
The Kronecker product operation of two y-axis steering vectors
Bz
The Kronecker product operation of two z-axis steering vectors
B˜y
The y-axis steering vectors of the uniform augmented array
B˜z
The z-axis steering vectors of the uniform augmented array
B˜y
The overcomplete azimuth angle sparse dictionary
B˜z
The overcomplete elevation angle sparse dictionary
y′
The variance of the potential azimuth angle signals
x′
The variance of the potential azimuth angle signals
x^
The estimated elevation angle
y^
The estimated azimuth angle
Data Availability
All data included in this study are available upon request by contact with the corresponding author.
Conflicts of Interest
The authors declare no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Acknowledgments
This work was supported by the National Natural Science Foundation of China under Contract nos. 61772398 and 61972239, the Key Research and Development Program Projects of Shaanxi Province (2019SF-257, 2020GY-024, and 2021GY-182), the Special Scientific Research Project of the Shaanxi Provincial Education Department (19JK0169), and the Science and Technology Program of Hantai District (2019KX-21).
WongK. T.Linshan LiL. S.ZoltowskiM. D.Root-MUSIC-based direction-finding and polarization estimation using diversely polarized possibly collocated antennas2004312913210.1109/lawp.2004.8310832-s2.0-11944258087WangL.WangG.ChenZ.Joint DOA-polarization estimation based on uniform concentric circular array201327131702171410.1080/09205071.2013.8231222-s2.0-84882328960YuanX.WongK. T.AgrawalK.Polarization estimation with a dipole-dipole pair, a dipole-loop pair, or a loop-loop pair of various orientations20126052442245210.1109/tap.2012.21897402-s2.0-84860684573WongK. T.LaiA. K.-Y.Inexpensive upgrade of base-station dumb antennas by two magnetic loops for “blind” adaptive downlink beamforming200547118919310.1109/map.2005.14362742-s2.0-20544469440SunF. G.WuQ. H.SunY. M.DingG. R.LanP.An iterative approach for sparse direction-of- arrival estimation in coprime arrays with off-grid targets2016613542LiuC.-L.VaidyanathanP. P.Cramér-Rao bounds for coprime and other sparse arrays, which find more sources than sensors201761436110.1016/j.dsp.2016.04.0112-s2.0-84971656929PalP.VaidyanathanP. P.Nested arrays: a novel approach to array processing with enhanced degrees of freedom20105884167418110.1109/tsp.2010.20492642-s2.0-77954583774SunF.LanP.GaoB.Partial spectral search‐based DOA estimation method for co‐prime linear arrays201551242053205510.1049/el.2015.22612-s2.0-84948416091VaidyanathanP. P.PalP.Sparse sensing with coprime arraysProceedings of the IEEE 44th Asilomar Conference on Signals Systems and ComputersNovember 2010Pacific Grove, CA, USA14051409PalP.VaidyanathanP. P.Coprime Sampling and the Music AlgorithmProceedings of the 2011 Digital Signal Processing and Signal Processing Education MeetingJanuary 2011Sedona, AZ, USA289294ZhangH.WuC. X.DOA estimation based on coprime array of covariance matrix reconstruction20184056468LiuQ. H.YiX. D.JinL. N.ChenW.Two dimensional direction of arrival estimation for co-prime L-shaped array using sparse reconstructionProceedings of the IEEE 8th International Congress on Image and Signal Processing (CISP)October 2015Shenyang, China14991503LiangJ.LiuD.Passive localization of near-field sources using cumulant20099895396010.1109/jsen.2009.20255802-s2.0-67749084306ShaoH.SuW. M.GuH.WangC.Two-dimensional direction finding using a sparse coprime L-sharped array [J]2012275886891YueY.XuY.ShenL.LiuZ.Parameter estimation of coexisted circular and strictly non‐circular signals2017531386486610.1049/el.2017.07302-s2.0-85021295449YangX.ZhengZ.HuB.Off‐grid DOA estimation of incoherently distributed non‐circular sources via generalised approximate message passing201652426226410.1049/el.2015.19732-s2.0-84958544433SebastianP.MartinH.CarlosM.DOA Estimation using random linear arrays via compressive sensing [J]2014125859863NorthardtE. T.BilikI.AbramovichY. I.Spatial compressive sensing for direction-of-arrival estimation with bias mitigation via expected likelihood20136151183119510.1109/tsp.2012.22326542-s2.0-84873934241