Abstract

For tacking and localizing sources in the mobile wireless sensor network, underdetermined direction of arrival (DOA) estimation with high-accuracy is a crucial issue. In this paper, a novel sparse array configuration is developed for accurate DOA estimation from the perspective of sum-difference coarray (SDCA). As compared with most of the existing sparse array configurations, the proposed array can effectively reduce the overlap between difference coarray (DCA) and sum coarray (SCA) and can achieve more consecutive degrees of freedom (DOF), more sources can be resolved accordingly. Additionally, the proposed array has hole-free DCA and SDCA. Then, the concept of coarray redundancy ratio (CRR) is introduced for evaluating the coarray overlap quantitatively and the closed-form CRR expressions of the proposed array are derived in detail. Based on the good properties of the proposed array, vectorized conjugate augmented MUSIC (VCAM) is adopted for underdetermined DOA estimation. The theoretical propositions and numerical simulations demonstrate the superior performance of the proposed array in terms of CRR, consecutive DOF, and DOA estimation accuracy.

1. Introduction

In the large-scale mobile wireless communication network [1, 2], underdetermined direction of arrival (DOA) estimation [36] with high-accuracy is a key issue to solve for tacking and localizing sources. The optimized array configuration design is a prerequisite for accurate DOA estimation. Recently, the emerging sparse arrays break through the constraint of spatial sampling theorem by sparsely arranging the array elements and have remarkable advantages in array element layout flexibility, degrees of freedom (DOF), and virtual array aperture [7, 8]. Benefiting from the sparse arrays, more sources than the number of the physical sensors can be resolved with high-accuracy.

From the view of array configuration design, coprime array (CPA) [9] and nested array (NA) [10] are the two most typical sparse arrays. In comparison to the minimum redundancy array and the minimum hole array, CPA and NA are easy to construct since their sensor locations are analytically tractable. Specifically, CPA is obtained from a pair of coprime uniform linear arrays (ULA), which can offer DOF with sensors, but has holes in the virtual coarray. By contrast, NA contains two “nesting” ULAs with increasing element spacing, which has a consecutive virtual coarray with DOF from sensors. On the basis of these prototypes, several modifications are developed, such as generalized CPA [11], thinned CPA [12], super NA [13], and augmented NA [14], aiming at increasing consecutive DOF and expanding effective array aperture. Nevertheless, note that these arrays construct synthetic virtual arrays from the view of difference coarray (DCA), the number of consecutive DOF in the virtual array cannot be more than twice of the physical aperture. Accordingly, the number of the resolvable sources and the DOA estimation accuracy are also affected.

Alternatively, several array configurations based on sum-difference coarray (SDCA) have attracted considerable interests, wherein the introduction of sum coarray (SCA) originating from multiple-input multiple-output (MIMO) system [15, 16] can bring more DOF. In [17], a novel CPA based on SDCA is presented, from which the consecutive DOF more than twice of array aperture can be obtained. Following this, a modified nested array configuration named as sum-diff NA (SdNA) is proposed in [18]. To further improve the available DOF, two kinds of improved NAs [19], termed as INAwSDCA-I and INAwSDCA-II, are developed by rearranging the translational NA. In [20], the unfold CPA configuration with SDCA is proposed. However, for the above array configurations based on SDCA, lots of overlapping virtual sensors between their DCA and SCA lead to unsatisfactory coarray utilization. In addition, there also exist holes in the SDCA. Toward this end, transformed NA (TNA) is designed in [21], which can reduce the overlap between DCA and SCA to some extent, but the degree of overlap under different number of array elements is only evaluated via simulations, lacking quantificational analysis.

To tackle these problems, a novel sparse array configuration with low coarray overlap is proposed in this paper. The proposed array is constructed with two concatenated subarrays with different interelement spacings. From the perspective of virtual array, both the DCA and SDCA of the proposed array are hole-free and can provide more consecutive DOF accordingly. Moreover, the coarray redundancy ratio (CRR) of the proposed array is derived in detail to evaluate the coarray overlap quantitatively. Based on the proposed array, the vectorized conjugate augmented MUSIC (VCAM) approach can be employed for DOA estimation. Now, we briefly summarize the contributions of this work as follows:(1)A novel array configuration is proposed to reduce the overlap between DCA and SCA(2)The proposed array has hole-free DCA and SDCA, which has been proved via theoretical propositions(3)More sources than twice of the number of sensors can be resolved

The mathematical notations used throughout this paper are denoted as follows. Vectors and matrices are denoted by lowercase and uppercase bold-face letters, respectively. , , and denote transpose, conjugate transpose, and conjugate. and represent the statistical expectation and modulus of the internal entity, denotes the derivation of with respect to a, denotes the nonnegative integer set, denotes the positive integer set, denotes a set of consecutive positive integers ranging from 1 to n in the case of , and holds in the case of ; otherwise, denotes a set of consecutive negative integers ranging from −n to −1. The symbol denotes the Kronecker product and denotes the Khatri–Rao product.

The remainder of this paper is organized as follows. The DOA estimation model based on SDCA is formulated in Section 2. Section 3 presents the proposed low coarray redundancy array configuration in detail. Section 4 exhibits the simulation results in terms of CRR, consecutive DOF, and DOA estimation accuracy. Conclusions are drawn in Section 5. The proofs of Propositions 1 and 2 are derived in Appendix A and Appendix B.

2. Problem Formulation

Assume that K far-field narrowband sources from directions impinge on an N-element nonuniform sparse antenna array. Then, the array output at time t can be represented aswhere the value of K is known in advance which can be estimated via source number estimation methods [22]. is the source vector with the element being , is the deterministic complex amplitude, and is the frequency offset. denotes the array manifold matrix. After normalizing by the unit element spacing d, the steering vector can be expressed as with being the sensor location vector. is the Gaussian white noise vector with zero mean and variance being . By collecting samples, the time average function can be expressed aswhere is the time lag, and are respectively the ith and the jth () row of , , and . For the pseudosnapshots, the pseudodata matrix is defined aswhere is the pseudosampling period that satisfies the Nyquist sampling theorem. For convenience, let the 1st sensor be the reference, . Then, the covariance matrix of is calculated as

Vectorizing yieldswhere and the kth column of can be denoted aswhere the term in (6) behaves like a virtual steering vector of a longer virtual ULA. More specifically, the union of and is respectively corresponding to DCA and its mirror version; and are respectively corresponding to SCA and its mirror version. Then, the spatial smoothing MUSIC or sparse construction approach can be performed for underdetermined DOA estimation, and the above procedure is termed as the vectorized conjugate augmented MUSIC (VCAM) [17].

3. Low Coarray Redundancy Array Configuration

Definition 1. Assume , , and , after normalizing by d, the configuration of the proposed array with antenna elements is defined byIntuitively, an example of is given in Figure 1, where white circles and black circles denote the sensor locations of subarray and , respectively.
From the perspective of coarray equivalence, the proposed array, NA, and TNA are compared in Figure 2 with the number of antenna elements being N = 6, where the red circles and blue triangles indicate the virtual sensors in DCA and SCA, respectively, and black crosses indicate holes. Note that the DCA, SCA, and SDCA are origin-symmetric, here the nonnegative parts of the abovementioned coarrays are selected for performance comparison, which are respectively abbreviated as n-DCA, n-SCA, and n-SDCA. The sensors in the rectangular dashed box indicate the overlapping virtual sensors between the n-DCA and the consecutive n-SCA. We can see from Figure 2 that the n-SDCA is continuous in the range of (0, 24) for the proposed array, while (0, 16) for NA and (0, 22) for TNA. Moreover, there exists only 1 overlapping virtual sensor for the proposed array, while 10 for NA and 4 for TNA.

Proposition 1. Some properties of the proposed array configuration in the virtual coarray are listed as follows:(i)The DCA of the proposed array is continuous in set of with (ii)The SCA of the proposed array is continuous in set of with , let , and the value of is given as follows:(a)if , (b)if , (c)if ,

Proof. See Appendix A.

Proposition 2. The proposed array has hole-free DCA and SDCA, and the maximum consecutive DOF is .

Proof. See Appendix B. Note that for a sparse array configuration with the fixed number of antenna elements, the higher consecutive DOF means larger virtual array aperture and more accurate DOA estimation. Therefore, it is essential to look for appropriate and to maximize the consecutive DOF.

Corollary 1. After optimizing and , the optimal consecutive DOF and the corresponding solutions are given as

Proof. The optimal DOF under the constraint of can be cast as the following optimization problem:Using arithmetic mean-geometric mean (AM-GM) inequalities, the solutions are obtained from (8). According to the above analysis, the proposed array exhibits the advantages in terms of consecutive DOF and virtual array aperture. From the view of coarray efficiency, the overlapping virtual sensors between the DCA and SCA lead to coarray redundancy, which, to some extent, reduce the number of available virtual coarrays. In view of this, the concept of CRR is introduced to evaluate the coarray overlap.

Definition 2. CRR denotes the ratio of the overlapping virtual sensors between the DCA and SCA, which is given byMore specifically, CRR of the proposed array configuration can be calculated as

Proposition 3. For the optimal DOF given in Corollary 1, if , the corresponding CRR is no more than 0.0496; if , the corresponding CRR is no more than 0.1053.

Proof. If , the corresponding CRR can be simplified to . Then, the first-order derivative of is taken as . Since is even, the maximum of is calculated as . If , the corresponding CRR can be simplified to . Similarly, the first-order derivative of is taken as . Note that is odd; the maximum of is calculated as .

4. Numerical Simulations

In this section, numerical simulations have been carried out to evaluate the performance of the proposed array, where SdNA [18] and TNA [21] are selected for performance comparison. In the first simulation, the CRR scatters versus are plotted in Figure 3, and the blue asterisks indicate the case where N is even and the red squares indicate the case where N is odd. It can be observed that the CRR of the proposed array tends to increase with the increase of , thus the optimal CRR can be obtained from when N is even and when N is odd.

In the second simulation, the CRR of three array configurations with and are respectively compared in Figures 4 and 5. The results show that the proposed array can achieve the smallest CRR compared to SdNA and TNA with the same number of antenna elements. Also, it has been concluded that the CRR of the proposed array is no more than 0.0496 when N is even and no more than 0.1053 when N is odd, which has been verified in Proposition 3.

In the third simulation, we compare the consecutive DOF of three array configurations and the simulation results are depicted in Figure 6. As can be seen from Figure 6, the proposed array can obtain more consecutive DOF than the other two arrays benefitting from the lower CRR. This implies that there exist few overlapping sensors between the DCA and SCA of the proposed array, thus more effective virtual sensors can be utilized for the construction of SDCA.

Then, the DOA estimation performance is investigated in the last simulation. We consider K = 27 narrowband sources uniformly distributed between and impinging on the proposed sparse array with 8 sensors, where SNR being 0 dB and . Figure 7 depicts the MUSIC spectrum of the proposed array where the blue solid lines denote the angle estimates and the red dotted lines denote the incident sources. It can be seen that the proposed array can detect 27 impinging sources with only 8 sensors, and the estimated spectrum has sharp and clearly discernible peaks in the vicinity of the true impinging sources.

Figures 8 and 9 respectively depict the root mean square error (RMSE) of DOA estimation versus SNR and snapshots via 200 Monte Carlo trials, where RMSE is defined aswith being the estimate of the kth impinging source for the nth Monte Carlo trial. All the simulation conditions are the same as Figure 7, except that the snapshot number is fixed at 200 in Figure 8 and the SNR is fixed at −6 dB in Figure 9. The results show that the proposed array is superior to the other two array configurations in terms of DOA estimation accuracy, which mainly attributed to the larger consecutive DOF and the lower CRR of the proposed array.

5. Conclusions

This paper presents a novel sparse array configuration for accurate DOA estimation, which has fewer coarray overlapping sensors and more consecutive DOF as compared with most of the existing sparse array configurations. From the perspective of coarray equivalence, the proposed array has hole-free DCA and SDCA. The closed-form CRR expressions of the proposed array are derived in order to evaluate the coarray overlap quantitatively. With the utilization of the proposed array, the VCAM approach is adopted for DOA estimation. The results of the theoretical analysis and the numerical simulations demonstrate the effectiveness and favorable performance of the proposed array configuration and DOA estimation performance.

Appendix

A. Proof of Proposition 1

(i)The positive part of the difference set between and is , which distributes continuously in the set of with . Similarly, the positive part of is continuous in the set of . Thus, the whole DCA (containing both the positive part and the negative part) of the proposed array is hole-free, which is continuous in the set of .(ii)The positive part of the sum set is continuous in the set of and is also continuous in the set with . Thus, the positive set of is continuous in the set . By contrast, the set is a discontinuous set composed of discrete points which are directly determined by and . More specifically, if , the largest element in is , where when and when . Then, the positive SCA can be written as with in the case of . If , the largest element in is , which means the positive SCA is with . If and , the largest element in is and , then the positive SCA is with . If and , the largest element in is , the corresponding positive SCA is with . Since the negative SCA and the positive SCA are symmetrical about the zero, the whole SCA of the proposed array is continuous in set of .

B. Proof of Proposition 2

Based on Proposition 1, the whole DCA is continuous in the set of with , which means the DCA of the proposed array is hole-free. Also, holds in the case of , then we have since (, ). This implies that falls within the continuous interval of DCA; as a result, the SDCA is hole-free. Similar analyses can be utilized for the cases of , , i.e., if or , ; if and , ; if and , . Consequently, the SDCA of the proposed array is continuous in the set of and can be simplified to , which means the proposed array has hole-free SDCA. Accordingly, the maximum consecutive DOF is .

Data Availability

No data were used in this study.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grant no. 62101223 and Major Basic Research Project of the Natural Science Foundation of the Jiangsu Higher Education Institutions of China under Grant nos. 20KJB510027 and 20KJA510008.