Computationally Efficient Direction Finding for a Mixture of Circular and Strictly Noncircular Sources with Uniform Rectangular Arrays

In this paper, a novel two-dimensional (2D) direction-of-arrival (DOA) estimation algorithm for the mixed circular and strictly noncircular sources is proposed. A general array model with a mixture of signals is firstly built based on uniform rectangular arrays (URAs), and then, the approach, which uses the rank-reduction-based ROOT-MUSIC, can solve 2D DOA estimation problem. Besides, the theoretical error of the proposed algorithm, a criterion of the performance for evaluation, is analyzed by the first-order Taylor expression using second-order statistics. As verified by the simulation results, a better DOA estimation performance and a lower computational complexity are achieved by the proposed algorithm than the existing methods resorting to the noncircularity of the incoming signals.


Introduction
Array signal processing has been widely used in the fields of sonar, radar, wireless communication, etc, and many excellent algorithms have been developed in the past few years. Especially the well-known multiple signal classification (MUSIC) algorithm has not only been applied in its time-reversal (TR) form to active location [1][2][3][4], but has also been widely used in direction-of-arrival (DOA) estimation (DOA is the research of passive location). In the last few years, DOA estimation became an important research branch of array signal processing [5][6][7]. Especially two-dimensional (2D) direction-of-arrival (DOA) estimation with different structured arrays, such as L-shaped uniform linear arrays (ULAs) [8][9][10], two-parallel ULAs [11][12][13], and uniform rectangular arrays (URAs) [14][15][16][17][18], has received much attention in past years. For URAs, the well-known multiple signal classification (MUSIC) algorithm can be used for 2D DOA estimation directly [16]; however, its computational complexity is very high. To overcome this problem, two efficient 2D DOA estimation methods have been proposed [17,18]. In [17], the reduced-dimension MUSIC algorithm was proposed, which reduced the computational complexity, and the 2D DOA estimation performance was very close to the 2D-MUSIC method. A preprocessing transformation matrix was introduced in [18], which transformed both the complex-valued covariance matrix and the complex-valued search vector into real-valued ones, then the 2D DOA estimation problem was decoupled into two successive real-valued one-dimensional (1D) DOA estimation problems with real-valued computations only.
However, only the covariance matrix was considered, which characterizes the circular Gaussian distribution in the above 2D DOA methods.
In recent years, the issue of utilizing noncircular information has been attracting increasing attention. Abeida et al. presented a theoretical analysis of the resolution of the conventional and noncircular MUSIC algorithms and proved that the noncircular MUSIC algorithm for the threshold array signal-to-noise ratios are very sensitive to the noncircularity phase separation of the sources [19]. Pascal et al. introduced the optimal widely linear (WL) minimum variance distortionless response (MVDR) beamformer for the reception of an unknown signal of interest (SOI) corrupted by potentially second-order (SO) noncircular background noise and interference [20]. Wan et al. [21] proposed an algorithm that utilized the noncircular characteristic to solve the DOA estimation of coherently distributed sources. Because the covariance matrix and the elliptic covariance matrix (which is also named the complementary covariance or pseudo-covariance) are used in non-circular signals direction-finding method simultaneously, the DOA estimation performance can be improved greatly by exploiting the noncircularity information. Gan et al. studied the non-circular characteristics of the signals and proposed an automatically-paired 2D DOAs estimation method based on non-circular signals [22]. Based on strictly non-circular signals, Steinwandt et al. proposed the high-resolution R − D NCstandard ESPRIT and the R − D NC unitary ESPRIT algorithms, which are applicable to arbitrary shift-invariant R − D antenna arrays and do not require a centro-symmetric array structure [23]. In [24], an extended rank reduction (ERARE) method was introduced for noncircular sources based on two-parallel ULAs, which made the estimation more accurate than that in [12].
However, in modern wireless communications, a usual situation is that some users send circular signals such as quadrature phase shift keying (QPSK) signals, but others send non-circular ones such as binary phase shift keying (BPSK) signals. Therefore, the signals impinging to the array may be mixed ones. As for the mixed signals situation, a key issue is how to estimate and distinguish the circular and non-circular signals. In [25,26], the mixed signals situation issue has been studied, and the corresponding methods based on 1D arrays have been proposed. Gao et al. [25] combined the observed data and their conjugate counterparts to construct two 1D DOA estimators for detecting the circular and non-circular signals, but the DOA estimation performance of the method degraded seriously in the condition of a small separating angle. Liu et al. proposed an improved algorithm in [26], which detected the circular and non-circular signals by using the difference between the circularity of the mixed sources. However, to the best of our knowledge, few research works have reported the 2D DOA estimation problem for mixed circular and non-circular signals. In [27], based on two-parallel uniform linear arrays (ULAs), Chen proposed an effective algorithm that combined the rank reduction method and MUSIC algorithm to form four 1D DOA estimators for solving the 2D mixed signals situation. Nevertheless, the computational complexity of the algorithm in [27] is still high, and the theoretical performance analysis is never mentioned in this research.
In this paper, we study the 2D DOA estimation problem based on uniform rectangular arrays (URAs) using the rank-reduction-based ROOT-MUSIC method and the theoretical performance analysis of the proposed algorithm. Firstly, we establish an array model with mixed circular and noncircular sources with URAs; secondly, to avoid seeking the peak of the spectrum and reduce the computation load, a novel algorithm based on ROOT-MUSIC and the rank-reduction method is proposed to solve the 2D DOA estimation issue; finally, the theoretical error of the proposed algorithm is derived as a benchmark. Particularly, the paper mainly discusses the uncorrelated signals impinging upon the array. If we utilize some decorrelation methods such as the spatial smoothing technologies [28,29] or four-order cumulants-based Toeplitz matrices reconstruction (FOC-TMR) method [13] to preprocess the correlated signals and obtain the de-correlated matrices, the proposed method can also be generalized to the case of correlated sources.
The rest of this paper is organized as follows. Section 2 presents the array signal model. The description of the proposed algorithm is introduced in Section 3. The theoretical error analysis of the proposed algorithm is derived in Section 4. Finally, the simulation results are given in Section 5, and conclusions are drawn in Section 6.
Notations: (·) * , (·) T and (·) H represent conjugation, transpose and conjugate transpose. E[·] is the expectation operation; diag(·) and blkdiag(·) stands for the diagonalization and block diagonalization operation, respectively; I ω denotes the ω × ω dimensional identity matrix; det[·] indicates the determinant of a matrix; arg(·) is the phase angle operator. Expect E n ,Ê n , E n1 ,Ê n1 , E n2 andÊ n2 representing the corresponding noise subspace; other variables have index n such as (·) n , which denote ones related to the non-circular signal. (·) c indicates a variable associated with the circular signal.
In practice, non-circularity and circularity are important properties of random variables; their concept directly comes from the geometrical interpretation of complex random variable. The signal would be called the circular source if its statistical characteristic has the rotational invariance characteristic, otherwise, it would be called the noncircular source. The work in [31] introduces the circularity and noncircularity in detail. Based on this, we only consider the rotational invariance characteristic of the first-and second-order statistical properties of the sources. For a complex random signal s k , k = 1, 2, ..., K, we defined E[s k ], E[s k s k * ] and E[s k 2 ] as the mean, the covariance and the elliptic covariance of the signal s k , respectively. If the source's first-and second-order statistical properties are rotational invariant for an arbitrary phase ϕ k as follows: E[s k e jϕ k (s k e jϕ k ) * ] = E[s k s k * ] The signal s k will be called the circular source. Conversely, the signal s k will be noncircular if the first-and second-order statistical properties are not rotational invariant. Beyond that, reference [32] proposes a model to describe signal sources with arbitrary second-order non-circularity. The source s k is defined as follows: where ϕ k is the rotation phase, η k (0 ≤ η k ≤ 1) denotes the non-circularity coefficient and s I k and s Q k represent the in-phase and quadrature components of the complex signal s k , respectively. Therefore, s k will represent a circular source if η k = 0, where the rotation phase ϕ k is irrelevant and undetermined. If η k = 1, s k will represent a strictly non-circular signal. Based on the above research, we defined the strictly non-circular signal as s n,k n (t) = b n,k ns n,k n (t), wheres n,k n (t) is a real signal and b n,k n = e jϕ kn (k n = 1, 2, . . . , K n ) is an arbitrary phase shift for the signal. Due to the phase information of circular source being irrelevant, it can be represented as s c,k c (t), k c = 1, 2, . . . , K c . Therefore, the mixed source signal vector can be modeled ass(t) = Bs(t). B can be expressed as: where

and s(t) is denoted as:
A new data vector f (t) is defined by concatenating the received data vectors x 1 (t), x 2 (t), ... , and x N (t) as follows: where C(θ, β) is the extend steering vector, and: In order to simplify the notation, the pair of angles (θ, β) and t is omitted. In Equation (11), a n,1 υ 1 n,1 a n,1 υ 2 n,1 . . . a n,1 υ N n,1 · · · a n,K n υ 1 n,Kn a n,K n υ 2 n,Kn . . .
C 1 is a N M × K n matrix of noncircular signals and C 2 is a N M × K c matrix of circular signals. As a classical procedure for existing noncircular DOA estimation algorithms [27,33,34], we can construct a new augmented data matrix f by combining the vector f and its conjugate counterpart f * as follows: The procedure can extend the array virtually and enlarge the aperture of the array antenna, and the estimating precision of DOA can be improved by utilizing the noncircularity of signals.
In Equation (14), C is a 2N M × (K n + 2K c ) matrix, which contains the new steering vectors of the impinging sources. s is a (K n + 2K c ) × 1 matrix of signals.
where c n,k n and c c,k c are the new steering vectors of non-circular signals and circular signals, respectively.
a n,k n υ 1 n,kn a n,k n υ 2 n,kn . . . a n,k n υ N n,kn which is a (K n + 2K c ) × 1 matrix of signals, and: which is a 2N M × 1 vector of noise.
The covariance matrix of f is calculated by: where R s = E[ s s H ] is the covariance matrix of s. R s is a full-rank matrix, since the incident signals are uncorrelated with each other. Then, the eigenvalue decomposition of R is: where the 2N M × (K n + 2K c ) matrix E s and the 2N M × (2NM − K n − 2K c ) matrix E n are the signal subspace and noise subspace, respectively.
Remark 1. In practice, the available observed data are finite. Thus, R can be approximated by: where L is the number of available data snapshots. Then, the eigenvalue decomposition ofˆ R is: therefore, the noise subspace E n and signal subspace E s can be approximated byÊ n andÊ s , respectively.

The Proposed Algorithm
In this section, we propose a 2D DOA estimation algorithm to solve the problem of estimating and distinguishing the mixed signals that are circular and non-circular in detail. Firstly, the method for estimating the mixed sources that are circular and noncircular is proposed; secondly, the algorithm, which only detects the circular signals, is proposed; and finally, we study how to distinguish these two kinds of sources.
Because both C and E s have the same signal subspace, orthogonal to the noise subspace spanned by the matrix E n , we devise estimators to obtain the 2D DOAs of noncircular and circular signals using the rank-reduction-based Root-MUSIC method.

2D DOA Estimation for Non-Circular Sources
Because the noise subspace E n is orthogonal to c n,k n (the steering vectors of non-circular signals), the following equation can be obtained directly: Then, together with Equations (16) and (24), we can get the following equation: a n,k n υ 1 n,kn a n,k n υ 2 n,kn . . . a n,k n υ N n,kn = E H n blkdiag[a n,k n , a n,k n , . . . , a n,k n , a * n,k n , a * n,k n , . . . , a * n,k n ].  Defining a vector p n (l) = [1, l, ..., l M−1 ] T , which is only related to l, where l = e −j(2π / λ)d x cos θ n,kn . Therefore, we can obtain the formula p n (l) = a n,k n and p n (l −1 ) = a * n,k n , then define a 2N M × 2N matrix Ω(l), which is only related to l, Ω(l) = blkdiag[p n (l), p n (l), · · · , p n (l), p n (l −1 ), p n (l −1 ), · · · , p n (l −1 )] = blkdiag[a n,k n , a n,k n , · · · , a n,k n , a * n,k n , a * n,k n , · · · , a * n,k n ] and define a 2N × 2N matrix: Note that if (2NM − K n − 2K c ) ≥ 2N and θ is not the true angle of the non-circular signal, Q n (l) is of full rank because in this case, the column rank of E n is not less than 2N. Then, Equation (25) holds true only when θ equals the true angle of the signal (Q n (l) drops rank). Since the covariance matrix of f is obtained from a finite number of samples, the reduction of the rank of Q n (l) can roughly be replaced by the minimum of the determinant of Q n (l). Therefore, we get the estimator of θ about non-circular signals as follows: Notice that f n (l) is a 2 × 3N(M − 1) order polynomial, showing that there are 3N(M − 1) pairs of conjugated roots. The estimates K n signal DOAs of θ can be obtained by finding the closest roots of the unit circle, which are given by: where l k n , k n = 1, ..., K n are the roots closest to the unit circle. Then, we take the estimatedθ k n of non-circular source into Equation (25) to get the estimator of β: where: u = e j(2π / λ)d y cos β n,kn ; it means that υ m n,kn = u m−1 , m = 1, · · · N and: To achieve the estimate of β, we need to solve the roots of f n (u). Due to f n (u) being a 2 × 3(N − 1)-order polynomial, it means there are 3(N − 1) pairs of conjugated roots. The estimate of β can be obtained by finding the closest roots of the unit circle; in this way, we can automatically pair the closest roots u 1 , u 2 , ..., u K n corresponding toθ 1 ,θ 2 , . . . ,θ K n , respectively, and the DOAs of β are given by:β where u k n , k n = 1, ..., K n are the roots closest to the unit circle, respectively.

2D DOA Estimation for Circular Sources
Because the noise subspace E n is also orthogonal to c c,k c (the steering vectors of circular signals), we can get the equation as follows: from Equation (34), we can get: Partitioning the noise subspace E n into E n = E n1 E n2 , where E n1 and E n2 are two submatrices of the same size N M × (2NM − K n − 2K c ), Equations (35) and (36) can be changed to: As proven in [27], Equations (37) and (38) are equivalent to each other; therefore, the estimator over θ, which corresponds to circular signals, can obtained based on Equation (37). Since: Defining a vector p c (l) = [1,l, ...,l M−1 ] T only related tol, wherel = e −j(2π / λ)d x cos θ c,kc . Therefore, we can get formula p c (l) = a c,k c , then defining a N M × N matrix Λ(l) only related tol, Λ(l)=blkdiag[p c (l), p c (l), · · · , p c (l)] = blkdiag[a c,k c , a c,k c , · · · , a c,k c ] and: Note that if (2NM − K n − 2K c ) ≥ N and when θ is not the true angle of circular signal, Q c (l) is of full rank because in this case, the column rank of E n1 is not less than N. Then, Equation (37) holds true only when θ equals the true angle of circular signal (Q c (l) drops rank). Since the covariance matrix of f is obtained from a finite number of samples, the reduction of the rank of Q c (l) can roughly be replaced by the minimum of the determinant of Q c (l). Therefore, we get the estimator of θ about circular signals as follows: Notice that f c (l) is a 2 × N(M − 1) order polynomial, which means that there are N(M − 1) pairs of conjugated roots, and the estimated K c signal DOAs of θ can be obtained by finding the closest roots of the unit circle, which are given by: wherel k c , k c = 1, ..., K c are the roots closest to the unit circle. Then, we take the estimatedθ k c of circular signals into Equation (37) to get the estimator of β: where: u = e j(2π / λ)d y cos β c, kc ; this means that υ m c,kc =ũ m−1 , m = 1, · · · N and: To obtain the estimate of β, we need to solve the roots of f c (ũ). According to f c (ũ) being a 2 × (N − 1)-order polynomial, this means there are N − 1 pairs of conjugated roots. The estimate of β can be obtained by finding the closest roots of the unit circle. Therefore, we can automatically pair the closest rootsũ 1 ,ũ 2 , ...,ũ K c corresponding toθ 1 ,θ 2 , . . . ,θ K c , respectively, and the DOAs of β are given by:β whereũ k c , k c = 1, ..., K c are the roots closest to the unit circle, respectively.

Identification of Circular and Noncircular Signals
In order to distinguish the 2D DOAs of circular and noncircular signals from the mixed signals, Equation (34) can be changed into: and: This means the noncircular method can also be applied to solve circular signals. Therefore, we can achieve the 2D DOAs of both non-circular and circular sources from Equations (28)- (30) and (33) and only obtain the 2D DOAs of circular sources from Equations (42)-(44) and (47). Then, the purpose of distinguishing the circular and noncircular signals from the mixtures can be accomplished.
The proposed algorithm can be outlined as: Step 1: Construct the new data vector f and calculate its approximate covarianceˆ R from Equations (14) and (22).
Step 2: Perform EVDtoˆ R and achieve the noise subspaceÊ n from Equation (23).
Step 6: Compare the estimate achieved by Step 3 and Step 5 to distinguish the K n 2D DOAs of the noncircular signals and the K c 2D DOAs of the circular ones.

Remark 2.
To calculateˆ R, a computational complexity of O((2N M) 2 L) is needed. The computational complexity of eigendecomposition operation is O ((2N M) 3 ). The proposed method employs four 1D estimators of polynomial-rooting; therefore, the complexity for the proposed method is While the algorithm in [27] employs several 1D spatial spectrum search procedures to obtain the 2D DOAs of signals, by defining the scanning interval of θ ∈ [0, π] with an interval of ∆θ and β ∈ [0, π] with an interval of ∆β, respectively, the complexity for the algorithm in [27] . The algorithm in [16] employs two direct 2D spatial spectrum search procedures, whose complexity for the method is O (NM) 2 L + (NM) 3 + π ∆θ π ∆β (2NM) 2 + π ∆θ π ∆β (NM) 2 . Therefore, the computational complexity of the proposed method has been reduced greatly.

Theoretical Performance Analysis
The theoretical DOA estimation error is caused by the finite data effect, the sensor errors and the unknown noise structure [31,35], and in this section, we will study the theoretical perturbation of the proposed algorithm as a criterion for evaluation.

Theoretical Perturbation for Non-Circular Sources
According to Equation (25), defining a n,k n = [1, l, ..., l M−1 ] T , υ m n,kn = u m−1 , m = 1, · · · , N and b n,k n = z, where l = e −j(2π / λ)d x cos θ n,kn , u = e j(2π / λ)d y cos β n,kn , z = e jϕ n,kn , ϕ n,k n is the non-circular phase of non-circular sources. Therefore, Equation (25) can be changed into: E H n blkdiag[a n,k n , a n,k n , · · · , a n,k n , a * n,k n , a * n,k n , · · · , a * n,k n ] M(l, u, z) = blkdiag[a n,k n , a n,k n , · · · , a n,k n , a * n,k n , a * n,k n , · · · , a * n,k n ].  Then, we can define the equation as follows: In practical situations, due to the influence of noise, the sensor errors and the finite data effect, we can only achieveÊ n , the approximation of E n , according to Equation (53), we can get the following formula:f (l, u, z) = M H (l, u, z)Ê nÊ H n M(l, u, z) The locations of the local minima of either the noise-free or perturbed spectral polynomials are obtained when the first derivative equal to zero. In other words, the roots of the derivative of f (l, u, z) give the locations of the relative minima off (l, u, z). Thus, to find the perturbations of the DOAs (which correspond to minima off (l, u, z)), we must calculate the perturbations of the roots of the derivative off (l, u, z). Define (l k n ,û k n ,ẑ k n ) as the estimate of (l, u, z) when the true angle is θ = θ n,k n , β = β n,k n , and ϕ = ϕ n,k n . Therefore, the first-order partial derivatives off (l, u, z) at (l k n ,û k n ,ẑ k n ) are zeros. Approximating the perturbation of DOAs uses the first two terms in the Taylor series expansion of the first-order partial derivative off (l, u, z) about the true angles of arrival.
f l (k n ) + ∆l k nf l,l (k n ) + ∆u k nf l,u (k n ) + ∆z k nf l,z (k n ) = 0 (55) where ∆l k n =l k n − l k n , ∆u k n =û k n − u k n and ∆z k n =ẑ k n − z k n are the perturbations of l k n , u k n and z k n , respectively. Together with Equations (55)-(57), solve the set of equations as follows: where: According to the expressions for Equation (54), the first derivative,f η (k n )(η = l, u, z), and the second derivative,f τ,κ (k n )(τ, κ = l, u, z), can be achieved, respectively. Using the value of (l k n , u k n , z k n ) and Equations (58) and (59), we can obtain the perturbation ∆l and ∆u. Then, the perturbation of θ and β can be achieved from ∆l and ∆u, respectively. The relation between an arrival angle and a signal root is given in [36]: ∆θ n,k n = C θ,k n Im ∆l k n l k n (61) ∆β n,k n = C β,k n Im ∆u k n u k n where C θ,k n = λ 2πd x sin θ n,kn and C β,k n = − λ 2πd y sin β n,kn . Therefore, the theoretical perturbation of θ n,k n and β n,k n can be obtained from Equations (61) and (62).

Theoretical Perturbation for Circular Sources
Similar to the noncircular method, defining a c,k c = [1,l, ...,l M−1 ] T and υ m c,kc =ũ m−1 , m = 1, · · · N, wherel = e −j(2π / λ)d x cos θ c,kc andũ = e j(2π / λ)d y cos β c,kc . Equation (35) can be changed to: defining: therefore, the equation can be defined as follows: In practical application, we can only achieveÊ n1 , the approximation of E n1 ; according to Equation (65), we can get the following equation: the roots of the derivative off (l,ũ) give the locations of the relative minima off (l,ũ). Define (l k c ,û k c ) as the estimate of (l,ũ) when the true angle is θ = θ c,k c and β = β c,k c . Therefore, the first-order partial derivatives off (l,ũ) at (l k c ,û k c ) are zeros. Approximating the perturbation of DOAs uses the first two terms in the Taylor series expansion of the first-order partial derivative off (l,ũ) at the true angles of arrival.f l (k c ) + ∆l k cf l ,l (k c ) + ∆ũ k cf l ,ũ (k c ) = 0 (67) where ∆l k c =l k c −l k c and ∆ũ k c =û k c −ũ k c are the perturbations ofl k c andũ k c , respectively. Together with Equations (67) and (68) to solve the set of equations as follows: where:δ According to the expressions for Equation (66), the first derivative,f η (k c )(η =l,ũ) and the second derivative,f τ,κ (k c ), (τ,κ =l,ũ) can be achieved, respectively. Using the value of (l k c ,ũ k c ) and Equations (69) and (70), we can obtain the perturbation ∆l and ∆ũ. Then, the relation between an arrival angle and a signal root is as follows: where C θ,k c = λ 2πd x sin θ c,kc and C β,k c = − λ 2πd y sin β c,kc ; therefore, the theoretical perturbation of θ c,k c and β c,k c can be obtained.

Simulation Results
In this section, simulation results are provided to demonstrate the performance of the proposed algorithm. For Simulations (1)-(3), the URAs have N = 3 rows and M = 8 columns; both d x and d y are half wavelength. For all simulations, we utilize BPSK and QPSK to represent the strictly non-circular signal and the circular signal, respectively, and the sources can be realized base on (7). For instance, if we take η k = 0 and ϕ k = 0, the QPSK signal s k ∈ 2 j can be obtained. If we take η k = 1 and ϕ k = 0, the BPSK signal s k ∈ (−1, +1) can be acquired. The power of additive white Gaussian noise is σ 2 n , and the signal-to-noise (SNR) is defined as SNR = 10log 10 (σ 2 s /σ 2 n ). We use the root mean square error (RMSE) to evaluate the estimation performance, which is defined as: where M c is the number of Monte Carlo simulations, K is the number of signals,ε m,k is the estimated θ k or β k in the m-th Monte Carlo simulation and ε k is the true value for either θ k or β k of the k-th signal.
The proposed algorithm in theoretical analysis and experimental results, the algorithm of 2D-MUSIC, the algorithm in [18], the algorithm in [27] and the deterministic CRB (Cramer-Rao bound) in [37], are compared in terms of RMSE. Figure 3a,b shows that the proposed method is steadily better than the other three algorithms, and the experimental values of the proposed algorithm overlap together with the theoretical error ones.    2D-MUSIC algorithm in [18] algorithm in [27] proposed algorithm(experimental) proposed algorithm(theoretical) CRB 2D-MUSIC algorithm in [18] algorithm in [27] proposed algorithm(experimental) proposed algorithm(theoretical) CRB

Performance versus Snapshots
We study the performance of the proposed algorithm with a varying snapshot number from 30-850. The SNR is fixed at 10 dB; the number of Monte Carlo simulations and incident signals are the same as in the second experiment. The RMSE results for angle estimation are shown in Figure 4a,b. The proposed method is steadily better than the other three algorithms, and the experimental values of the proposed algorithm overlap together with the theoretical error ones. 2D-MUSIC algorithm in [18] algorithm in [27] proposed algorithm(experimental) proposed algorithm(theoretical) CRB 2D-MUSIC algorithm in [18] algorithm in [27] proposed algorithm(experimental) proposed algorithm(theoretical) CRB

Performance versus the Number of Noncircular Mixed Signals
Finally, we consider the performance of the number of noncircular signals based on the proposed algorithm. The URAs have N = 3 rows and M = 6 columns; the SNRs vary form 0 dB-30 dB; the number of snapshots is 1200; and the number of Monte Carlo simulations is 1000. There are four uncorrelated signals from directions (θ 1 = 65 • , β 1 = 50 • ), (θ 2 = 105 • , β 2 = 95 • ), (θ 3 = 75 • , β 3 = 40 • ) and (θ 4 = 115 • , β 4 = 70 • ), and the total number of sources remains unchanged. We consider cases about one, two, three and four BPSK signals, respectively. The results are shown in Figure 5a,b. We can see that the RMSE of the experimental values of the proposed algorithm will overlap together with the theoretical error ones in a relatively high SNR, and the 2D DOA estimation performance of the proposed method improves from Case 1 to Case 4 because the dimension of the noise subspace has been extended by the increasing number of BPSK signals.

Conclusions
In this paper, a novel low-complexity 2D DOA estimation algorithm for mixed circular and non-circular signals has been proposed based on the URAs; besides, the theoretical error of the proposed algorithm is analyzed. As verified by the simulation results, the proposed method has a lower computational complexity and a better DOA estimation performance than the existing methods due to the usage of the noncircularity of the incoming signals.