Abstract

We consider the computationally efficient direction-of-arrival (DOA) and noncircular (NC) phase estimation problem of noncircular signal for uniform linear array. The key idea is to apply the noncircular propagator method (NC-PM) which does not require eigenvalue decomposition (EVD) of the covariance matrix or singular value decomposition (SVD) of the received data. Noncircular rotational invariance propagator method (NC-RI-PM) avoids spectral peak searching in PM and can obtain the closed-form solution of DOA, so it has lower computational complexity. An improved NC-RI-PM algorithm of noncircular signal for uniform linear array is proposed to estimate the elevation angles and noncircular phases with automatic pairing. We reconstruct the extended array output by combining the array output and its conjugated counterpart. Our algorithm fully uses the extended array elements in the improved propagator matrix to estimate the elevation angles and noncircular phases by utilizing the rotational invariance property between subarrays. Compared with NC-RI-PM, the proposed algorithm has better angle estimation performance and much lower computational load. The computational complexity of the proposed algorithm is analyzed. We also derive the variance of estimation error and Cramer-Rao bound (CRB) of noncircular signal for uniform linear array. Finally, simulation results are presented to demonstrate the effectiveness of our algorithm.

1. Introduction

Over the last several decades, the problem of estimating the direction-of-arrival (DOA) of multiple sources in the field of array signal processing has received considerable attention [13]. A variety of DOA estimation algorithms have been developed and applied in many fields, such as mobile communication system, radio astronomy, sonar, and radar [413]. Although the maximum likelihood estimator provides the optimum parameter estimation performance, its computational complexity is extremely demanding [810]. Simpler but suboptimal solutions can be achieved by the subspace-based approaches, which rely on the decomposition of observation space into signal subspace and noise subspace. For example, multiple signals classification (MUSIC) method [11] and estimation of signal parameters via rotational invariance technique (ESPRIT) algorithm [12, 13] are subspace-based DOA estimation algorithms both well known for their good angle estimation performance. However, conventional subspace techniques necessitate eigenvalue decomposition (EVD) of covariance matrix or singular value decomposition (SVD) of data matrix to estimate the signal and noise subspaces; thus, huge computation complexity will be involved, particularly when the number of sensors is large, such as the large towed arrays in sonar [14, 15]. It is well known that propagator method (PM) does not require EVD of covariance matrix or SVD of received data; thus, the computational load of PM algorithm can be significantly smaller [16]. But the spectral peak searching process is used in conventional PM algorithms [1720]; in order to save the complexity, rotational invariance PM (RI-PM) algorithms are proposed [2123], which avoid the spectral peak searching process and can obtain the closed-form solution of DOA.

The binary phase shift keying (BPSK), amplitude modulation (AM), and unbalanced quadrature phase shift keying (UQPSK) modulated signals frequently used in communication systems are noncircular (NC) signals [24]. The statistical parameters of noncircular signal, such as first and second moments, are rotational variant. The noncircularity of signal is investigated to enhance the performance of angle estimation algorithm [25]. Many DOA estimation methods of noncircular signals have been reported, which contain NC-MUSIC algorithms [2628], NC-ESPRIT algorithms [2931], and noncircular parallel factor (NC-PARAFAC) algorithm [32]. These algorithms we mention above can estimate more sources and have better angle estimation performance by introducing the noncircularity of signal into the conventional DOA estimation algorithms.

Many DOA estimation algorithms of noncircular signal based on PM have also been studied [3338]. These noncircular propagator method (NC-PM) algorithms can be divided into two kinds, that is, the noncircular spectral peak searching PM algorithms and noncircular rotational invariance PM (NC-RI-PM) algorithms. In fact, Liu et al. have proposed NC-PM algorithm [34] and real-valued NC-PM algorithm [35]; unfortunately, the spectral peak searching process is involved, which generally requires highly computational complexity. Subsequently, Liu et al. proposed NC-root-PM algorithm [36, 37] by exploiting the uniform distributed characteristic of the array elements, but its complexity is still high. The NC-RI-PM algorithm has been proposed by Sun and Zhou in [38], which has much lower computational complexity than Liu’s algorithm and can obtain the closed-form solution of DOA. However, we find that the reconstructed array output can be optimized, and some array elements in the extended propagator matrix have not been used; thus, the angle estimation performance of NC-RI-PM algorithm can be enhanced further. It is worth noting that the rotational invariance PM (RI-PM) algorithms have also been studied in [3942], but the noncircularity of signal has not been considered.

The aim of this paper is to develop a computationally efficient parameter estimation algorithm of noncircular signal. We propose an improved noncircular rotational invariance PM (improved NC-RI-PM) algorithm of noncircular signal for uniform linear array, which combine the noncircularity of signal and rotational invariance property between subarrays. Compared with NC-RI-PM algorithm, the proposed algorithm has the following advantages: (1) It has better angle estimation performance than that of NC-RI-PM algorithm. (2) The proposed algorithm has much lower computational complexity. (3) It can estimate elevation angles and noncircular phases with automatic pairing.

The reminder of this paper is organized as follows. In Section 2, we describe the data model for uniform linear array. The improved NC-RI-PM algorithm is presented in Section 3, as well as the computational complexity analysis and comparison. Section 4 derives the variance of estimation error and Cramer-Rao bound (CRB) of the proposed algorithm. Simulation results are presented in Section 5, while the conclusions are shown in Section 6.

Notations. , and denote inverse, conjugate, transpose, conjugate-transpose, and pseudoinverse operations, respectively. stands for a diagonal matrix, whose diagonal element is the vector . is identity matrix. is the expectation operator. means to get the phase. denotes the Hadamard product. means the estimated value of . denotes the Frobenius norm.

2. Data Model

We assume that there are uncorrelated narrowband signals impinging on a uniform linear array equipped with antennas as shown in Figure 1, and . We also assume that sources are far away from the array; thus, the incoming waves over the antennas are essentially planes. The noise is additive independent identically distributed Gaussian with zero mean and variance , independent of signals. We denote the elevation angles of sources as , where is the elevation angle of source, .

From [38], the array output at time can be modeled as where is the steering matrix, , and , is the distance between adjacent antennas, and is the wavelength. is the signal vector. is the noise vector, and , .

We just consider the signal of maximum noncircular rate in this paper. According to the noncircularity of signals, can be denoted as where . is a diagonal matrix; where is the noncircular phase of signal.

3. Angle Estimation Algorithm

In [38], Sun and Zhou have proposed the NC-RI-PM algorithm for uniform linear arrays which has much lower computational complexity than NC-PM algorithm based on spectral peak searching and can obtain the closed-form solution of DOA. However, we find that NC-RI-PM algorithm just uses part of the array elements in the extended propagator matrix, and the reconstructed array output data can be optimized. Thus, we propose an improved NC-RI-PM algorithm of noncircular signal for uniform linear array, which is presented as follows.

3.1. Improved NC-RI-PM Algorithm

According to (1), we construct the extended array output as where . is the permutation matrix, where is a diagonal matrix, and

Consider that channel state information is constant during sampling process. The sample data of extended array output can be expressed as where and are the sample signal matrix and sample noise matrix, respectively, , and is the number of snapshots.

We partition the extended steering matrix as where is a nonsingular matrix, . From [43], is a linear transformation of , where is the improved propagator matrix.

According to (4), the covariance matrix of the extended array output is [38]where is the covariance matrix.

Partition the covariance matrix as where , .

In the noiseless case, we can have

Actually, there is always noise, and the propagator matrix can be estimated by the following minimization problem [43]:where denotes the Frobenius norm. The estimate of is via

Then define where . In the noiseless case, according to (9)-(10), we can have

In the following, we estimate the elevation angles and noncircular phases via utilizing the rotational invariance property between subarrays. Partition the matrix into two parts as where , .

Define where and are the selective matrix.

Let

According to (6)-(7) and (21),

Define

According to (22)-(23), we can have

Perform the EVD of , which can be written as where , . The eigenvalues of in are corresponding to the diagonal elements of , and the eigenvectors of in are estimates of the column vectors in matrix .

Thus, the elevation angle of source can be estimated bywhere is the estimate of .

Next, we can also estimate the noncircular phases. According to (6) and (17)-(18), we can have where is a diagonal matrix.

Define

According to (27)-(28), we can get

Then can be written via EVD as where and . The eigenvalues of in corresponding to the diagonal elements of and the eigenvectors of in are estimates of the column vectors in matrix .

Note that the EVD of and are performed, respectively; we should consider the column and scale ambiguity between and before estimating the noncircular phases. According to (23) and (28), we can find that and are all computed from the matrix , so replace and in (28) with

Thus, the estimated noncircular phases can be expressed as where is the estimate of noncircular phase.

3.2. Algorithm Steps

The implementation of the proposed algorithm with finite array output data is summarized in this section. According to (4), the covariance matrix of sample extended array output data can be expressed as

We show the major steps of the improved NC-RI-PM algorithm of noncircular signals for uniform linear arrays as follows.(1)Construct the sample array output matrix and corresponding covariance matrix .(2)Partition into two parts as (12), and compute the propagator matrix .(3)Compute and , and perform the EVD of and as (25) and (30).(4)Estimate the elevation angles and noncircular phases via (26) and (32).

3.3. Differences between Our Works and NC-RI-PM Algorithm

NC-RI-PM algorithm [38] is simply reviewed for comparison in this section. According to (1), the extended subarray output in [38] is reconstructed aswhere , , , and . and are the selective matrices defined in (19) and (20), respectively.

Definewhere , , , and .

The sample array output of can be expressed as , where , . Suppose that is the sample covariance matrix, which is partitioned as , where and . Then the propagator matrix of NC-RI-PM algorithm can be computed by , where . Partition the propagator matrix as , where , , and .

Based on the rotational invariance property between and , we can get where , , , and , .

Thus, the elevation angles can be estimated by performing the EVD of . With the above analysis, we summarize the differences between the improved NC-RI-PM algorithm and NC-RI-PM algorithm in the following.(1)According to (34)-(35), , there are duplicated data in . These duplicated data have no benefit to improve the angle estimation performance but increase the computational complexity. The dimension of extended covariance matrix in NC-RI-PM algorithm is , whereas that of in our algorithm is just . So our algorithm has much lower computational complexity than NC-RI-PM algorithm.(2)According to (36), NC-RI-PM algorithm just uses part of the extended array elements in the propagator matrix , that is, the array elements in and , but the array elements in have not been used, whereas the improved NC-RI-PM algorithm fully uses all the elements in the extended propagator matrix. Thus, the angle estimation performance of our algorithm is better than that of NC-RI-PM algorithm.(3)The proposed algorithm can estimate elevation angles and noncircular phases with automatic pairing, while just the elevation angles are estimated in NC-RI-PM algorithm.

Remark 1. We assume that the number of sources is preknown, and it can be estimated by some methods shown in [4446].

Remark 2. According to (23) and (28), the elevation angles and noncircular phases are all estimated from the matrix , and the column and scale ambiguity can be solved by (31). Thus, the elevation angles and noncircular phases can be pairing automatically.

3.4. Complexity Analysis

Regarding the major computational complexity, we just consider matrix complex multiplication operations. It is known that the complexity of computing the covariance matrix with snapshots is in the order of and that of the EVD of dimension matrix is . In this paper, computational complexity of computing the matrices , , and and EVD of , are analyzed for the improved NC-RI-PM algorithm. Table 1 presents the total computational complexity of the improved NC-RI-PM algorithm, NC-ESPRIT algorithm in [29], and NC-RI-PM algorithm in [38]. Figure 2 shows the simulation results of computational complexity comparison versus and , respectively. From Figure 2, we can find that our algorithm has much lower computational complexity than that of NC-RI-PM algorithm and has approximate but still lower complexity than NC-ESPRIT algorithm.

4. Error Analysis and CRB

4.1. Error Analysis

This section analyses the variance of estimation error for the improved NC-RI-PM algorithm. We assume that where is the estimation error of covariance matrix.

According to (12), we can havewhere and are the estimation error matrix of and , respectively.

Combine (15) and (38); the estimate of the propagator is

Define where is the estimation error of the propagator matrix.

Thus, the estimation error of matrix is

According to (16)–(18) and (41), where and are the estimation error of matrices and , respectively. is the estimate of the steering matrix .

According to (21) and (42),

Define , . According to [43], we can obtain the estimate of matrix as

Let denote the eigenvalues of , which can be written as where is the estimation error of , , and is a unit column vector, .

Based on the first-order Taylor series expansion, the variance of estimation error of elevation angle can be expressed as

4.2. Cramer-Rao Bound

We derive the CRB of noncircular signal for uniform linear array of the proposed algorithm in this section. The parameters needed to be estimated can be denoted as [32]where and denote the real and imaginary parts of , respectively.

According to (8), the extended array output with snapshots can be rewritten as

The mean and covariance   of are

From [47], the element of the CRB matrix can be expressed as where and denote the first-order derivative of and with respect to the element of , respectively.

For the covariance matrix is just related to and the first term of (51) can be ignored. The element of CRB matrix can be simplified as

According to (48) and (52), where is the element of , , and is the column vector of extended direction matrix .

Define

Then we can have

According to (49) and (54)-(55), the first-order derivative of with respect to is

Combine (50) with (56); (52) can be rewritten as where

Let where and are the real and imaginary parts of , respectively.

According to (56) and (59), we can demonstrate thatwhere and . Consider

We just consider the elements related to the angles. According to (61), we can have aswhere denotes the parts unrelated to the elevation angles.

Until now, we obtain the CRB matrix as follows:

After further simplification, the CRB matrix can be rewritten as where , , , and is the power of noise. denotes the Hadamard product.

5. Simulation Results

The Monte Carlo simulations are adopted to evaluate the angle estimation performance of the proposed algorithm. We define the root-mean square error (RMSE) as where is the estimate of of the Monte Carlo trial and is the number of simulation loops. In the following simulations except Figures 6 and 8, we assume that there are sources, which are located at angles of , respectively. The corresponding noncircular phases are , respectively. We also assume that all sources have the same symbol duration and the same input signal-to-noise ratio (SNR), and the distance between adjacent antennas is equivalent to half of the wavelength, .

Figure 3 presents angle estimation results of elevation angles and noncircular phases of the proposed algorithm. and are used in the simulations, while SNR = 5 dB is used for elevation angles and SNR = 20 dB is used for noncircular phases. From Figure 3, we can see that the elevation angles and noncircular phases can be clearly observed.

Figure 4 shows the angle estimation performance comparison among the proposed algorithm, NC-RI-PM algorithm [38], NC-ESPRIT algorithm [29], and CRB of noncircular signals for uniform linear array. and are used. It is indicated in Figure 4 that the angle performance of the improved NC-RI-PM algorithm is better than NC-RI-PM algorithm, since the proposed algorithm fully uses all the elements in the extended propagator matrix. The angle estimation performance of our algorithm is close to that of NC-ESPRIT algorithm. Because the noncircular phases have not been estimated in NC-RI-PM algorithm and NC-ESPRIT algorithm, we just consider elevation angles in Figure 4 for angle performances comparison.

Figure 5 shows angle estimation performance of the proposed algorithm with and different values of . It is indicated in Figure 5 that the increasing of will lead to the improvement of angle estimation performance of the proposed algorithm.

Figure 6 depicts angle estimation performance of the proposed algorithm with and different values of . From Figure 6, the angle estimation performance of the proposed algorithm is enhanced with the number of snapshots increasing.

Figure 7 presents angle estimation performance of the proposed algorithm with , , and different values of . From Figure 7, we find that the angle performance of our algorithm degrades with the number of sources increasing.

Figure 8 displays the simulation result of the proposed algorithm with two closely spaced sources. We assume that the two close sources are located at angles of , and the corresponding noncircular phases are , respectively. , , and SNR = 15 dB are used in the simulation. It implies in Figure 8 that our algorithm works well when two sources are closely spaced.

6. Conclusions

In this paper, we have proposed an improved NC-RI-PM algorithm of noncircular signals for uniform linear array. The proposed algorithm has the following advantages: (1) It has much lower computational complexity, particularly when the ratio of the number of antennas to the number of sources is large. (2) It has better angle estimation performance than NC-RI-PM algorithm by constructing an improved propagator matrix, which fully uses all the array elements in the extended propagator matrix. (3) It can estimate elevation angles and noncircular phases with automatic pairing. (4) The maximum number of sources estimated by our algorithm is two times of that of PM algorithm via utilizing the noncircularity of signal.

We also analyze the computational complexity of the proposed algorithm. The error analysis of the proposed algorithm and CRB of noncircular signals for uniform linear array are also derived. It is well known that CRB expresses a lower bound on the variance of an unbiased estimator, which can be used to compare the angle performance of different algorithms. From Figure 4, we can see that the RMSE of the improved NC-RI-PM algorithm is much closer to CRB compared with that of NC-RI-PM algorithm. It is clearly indicated that the angle performance of the improved NC-RI-PM algorithm is better than that of NC-RI-PM algorithm, since the proposed algorithm fully uses all the elements in the extended propagator matrix. Finally, the angle estimation performance and computational complexity of the proposed algorithm are evaluated by computer simulations. Simulation results illustrate the effectiveness of the proposed algorithm in a variety of scenarios, particularly when the sources are closely spaced. Therefore, the proposed algorithm can be regarded not only as an improvement of the work in [38] but also as an improved DOA and noncircular phase estimation algorithm.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by the following research projects: China NSF Grants (61371169, 61471191, 61471192, and 61271327), Jiangsu Planned Project for Postdoctoral Research Funds (1201039C), China Postdoctoral Science Foundation (2012M521099, 2013M541661), Open Project of Key Laboratory of Modern Acoustic of Ministry of Education (Nanjing University), The Aeronautical Science Foundation of China (20120152001), Qing Lan Project and priority academic program development of Jiangsu high education institutions, The Fundamental Research Fund for the Central Universities, and the Colleges and Universities in Jiangsu Province Plans to Graduate Research and Innovation (CXLX11_0196). The authors are grateful to the anonymous referees for their constructive comments and suggestions in improving the quality of this paper.