# Thermos Array: Two-Dimensional Sparse Array with Reduced Mutual Coupling

Authors: Lei Sun; Minglei Yang; Baixiao Chen
Journal: International Journal of Antennas and Propagation
 (2018)

---

## Abstract

Sparse planar arrays, such as the billboard array, the open box array, and the two-dimensional nested array, have drawn lots of interest owing to their ability of two-dimensional angle estimation. Unfortunately, these arrays often suffer from mutual-coupling problems due to the large number of sensor pairs with small spacingd (usually equal to a half wavelength), which will degrade the performance of direction of arrival (DOA) estimation. Recently, the two-dimensional half-open box array and the hourglass array are proposed to reduce the mutual coupling. But both of them still have many sensor pairs with small spacing d, which implies that the reduction of mutual coupling is still limited. In this paper, we propose a new sparse planar array which has fewer number of sensor pairs with small spacing d. It is named as the thermos array because its shape seems like a thermos. Although the resulting difference coarray (DCA) of the thermos array is not hole-free, a large filled rectangular part in the DCA can be facilitated to perform spatial-smoothing-based DOA estimation. Moreover, it enjoys closed-form expressions for the sensor locations and the number of available degrees of freedom. Simulations show that the thermos array can achieve better DOA estimation performance than the hourglass array in the presence of mutual coupling, which indicates that our thermos array is more robust to the mutual-coupling array.


---

## Body

## 1. Introduction

Many planar arrays, such as the billboard array [1], open box array [2], and the two-dimensional nested array [3], can estimate the azimuth and elevation angles of sources simultaneously. And these arrays have been demonstrated to be able to resolve more sources than number of sensors by exploiting the second-order statistic information of data received by sensors. However, in real-word application, the mutual coupling among the adjacent sensors [4] with small spacing cannot be ignored. And the angle estimation performance of these aforementioned arrays often degrade significantly due to the mutual coupling.In order to solve the problem of mutual coupling, two different approaches have been studied. The first one is developing some estimation methodologies based on mutual-coupling models [5–7]. Unfortunately, as mentioned in [8], these estimation methodologies often fail when there is model mismatch. Based on the principle that large interelement spacing can reduce the effect of mutual coupling, the second approach is to design new arrays in which the number of sensor pairs with small spacing d (usually equal to a half wavelength) is greatly reduced. The coprime array [9, 10] and the super nested array [8] belong to this type for linear array design. However, these linear arrays can only estimate the azimuth angles or the elevation angles of the sources, that is, they cannot estimate the two-dimensional angles simultaneously. Recently, the two-dimensional half-open box array and hourglass array were proposed [11], which are capable of reducing mutual coupling by placing the sensors properly. In the meantime, the difference coarrays (DCAs) of these arrays can also be kept as filled rectangular arrays. Among the planar arrays mentioned above, the hourglass arrays show the best performance when mutual coupling exists.Figure1(a) shows an example of the hourglass array. We can observe that the four corners in the hourglass array are dense parts with three or four sensors with small spacing, which will bring in the mutual-coupling problem among these sensors. As reported in [11], the number of sensor pairs with crucial small separation of (0, d) is four or five times larger than those with separation of (d, 0), where d is the smallest interelement spacing. This inspires us that if we reduce the number of sensor pairs with separation of (0, d) to the same level as those with separation of (d, 0), we can further reduce the effect of mutual coupling.Figure 1
(a) Array geometries of the hourglass array [11] and (c) thermos array (normalized by d) in which Nx=15 and Ny=27, and the total number of sensors is K=67. (b) The resulting DCA of the hourglass array and (d) resulting DCA of the thermos array.
(a)
(b)
(c)
(d)In this paper, we propose a new array configuration called thethermos array, which has fewer number of sensor pairs with small separations than the hourglass array, and hence, the thermos array is more robust to mutual coupling. More quantitative analyses are listed in Table 1, wherein Nx≥2L+1 and Ny≥9 (the size of the array), where
(1)L=Ny8+1,ifNyiseven,Ny+14,ifNyisodd,Table 1
Crucial weight comparison of two array geometries.
Array geometries
Nx, Ny
V1,0
V0,1
V1,0+V0,1
Hourglass array
Ny is even
2
10
12
Ny is odd
2
8
10
Thermos array
Nx is even
3
3
6
Nx is odd
2
2
4and it is defined in [11]. Comparing with the hourglass array, the thermos array can provide more degrees of freedom (DOFs) in the available difference coarray when Nx and Ny (the size of the hourglass array) satisfy certain conditions. Furthermore, the thermos array enjoys closed-form expressions for both the sensor locations and the number of available DOF.The rest of the paper is organized as follows. The signal model is presented in Section2. The proposed array geometry and array characteristic comparison with the hourglass array are developed in Section 3, and Section 4 shows the DOA estimation algorithm. Section 5 presents numerical examples. Section 6 concludes the paper.


## 2. Signal Model

Consider aK-element planar array with sensors located at
(2)v⋅d=v1,v2,…,vKd,where vi=xi,yii=1K denotes the position vector of sensors in the plane and xi,yii=1,…,K are integers. The d is the unit interelement spacing, usually equal to a half wavelength. For simplicity and intuition, we set d as 1 in the array position analysis hereafter. Let Q serve as the uncorrelated narrowband sources that impinge on the array from directions ϕq,θq,q=1,2,…,Q, where ϕq and θq indicate the azimuth and elevation of the qth source, respectively. A steering vector is the array response to a unit strength source at direction ϕq,θq
(3)aϕq,θq=ejκx1θ¯+y1ϕ¯,ejκx2θ¯+y2ϕ¯,…,ejκxKθ¯+yKϕ¯T=ejκv1θ¯,ϕ¯T,ejκv2θ¯,ϕ¯T,…,ejκvKθ¯,ϕ¯TT,where κ=2πd/λ,θ¯=cosϕqsinθq,ϕ¯=sinϕqsinθq, and λ is the signal wavelength.The data received at the array can be expressed as(4)xt=Ast+nt,t=1,2,…,T,where A=aϕ1,θ1,aϕ2,θ2,…,aϕQ,θQ is the array manifold matrix and T is the number of snapshots. The source signals st=s1t,s2t,…,sQtT are assumed unknown, but each source sqt is assumed to follow a complex Gaussian distribution CN0,σq2. The source signals are assumed mutually independent and uncorrelated over time. The components of the noise vector nt are assumed independent and identically distributed (i.i.d.) additive white Gaussian noise with zero mean and variance σn2 and are independent from the sources. The covariance matrix of the received data xt can be expressed as
(5)Rxx=ExtxHt=ARssAH+σn2IK=∑q=1Qσq2aϕq,θqaHϕq,θq+σn2IK,where Rss is a Q×Q diagonal matrix with diagonal elements σ12,σ22,…,σQ2 and IK is a K×K identity matrix.Vectorizing matrixRxx, we obtain a K2×1 vector
(6)z=vecRxx=Bc+σn21n,where vec· is the vectorization operation that stacks all columns of a matrix on top of one another into a vector.
(7)B=a∗ϕ1,θ1⊗aϕ1,θ1,…,a∗ϕQ,θQ⊗aϕQ,θQT=bϕ1,θ1,bϕ2,θ2,…,bϕQ,θQ∈ℂK2×Q,and c=σ12,σ22,⋯,σQ2T is interpreted as an equivalent source signal vector. The symbol ⊗ in (7) stands for the Kronecker product, and 1n=e1T,e2T,…,eKTT, with ei denoting a column vector of all zeros, except the ith element, which is equal to one. The vector z is equivalent to the data received by a virtual array with elements located at the location set vi−vj,i,j=1,2,…,K. By averaging the multiple data associated with the virtual sensors at the same location in vi−vj, we obtain another virtual array that has sensors with averaged data z~ located at the distinct location of vi−vj, which are named as difference coarray and defined as D.When considering the mutual coupling among sensors, we can represent the received data as(8)xt=CAst+nt,where C is the mutual-coupling matrix with the (i,j)th element written as [11]
(9)Cij=cvi−vj2,ifvi−vj2≤B,0,others,where c· represents the mutual-coupling coefficient, ·2 denotes the l2-norm operation, and B means the maximum spacing between two sensors in the planar array.As we all know, sensor pairs with larger separation can reduce the effect of mutual coupling. We use a crucial weight to analyze the mutual coupling among sensors, and it is defined as follows.Definition 1
(crucial weight). The crucial weightVL is the number of sensor pairs with separation L:
(10)VL=cardvi−vj=cardxi−xj,yi−yj=cardL∣L=1,0or1,0,where card· and ⋅ are the operation of cardinality and absolute value, respectively. Take the four sensors located at the lower left corner in Figure 1(a), for example, their positions are (0, 0), (0, 1), (0, 2), and (1, 0), respectively. Then there is only one sensor pair with separation L=1,0 and two sensor pairs with separation L=0,1. Thus, V1,0=1 and V0,1=2.This definition gives the number of sensor pairs with the smallest spacingd in a planar array, and thus, the smaller number of crucial weight implies lower mutual coupling among sensors. Therefore, when designing an array which is robust to the unknown mutual coupling, a smaller value of the crucial weight is desired. It should be noted that the definition of crucial weight here is just a small part of the weight function defined in [11]. But our definition is more convenient for the analysis of mutual coupling.In the following section, we will design a new planar array with less mutual coupling and analyze its performance with the crucial weight.


## 3. Thermos Array

The hourglass array [11] is a generalized open box array, and its name comes from the array geometry resembling an hourglass with two pillars on both sides. By reducing the number of sensor pairs with separation (1, 0) and (0, 1), the hourglass array shows better DOA estimation performance than the half-open box array (HOBA) and the half-open box array with two layers (HOBA-2) when the mutual coupling is presented. Unfortunately, as mentioned above, the value of the crucial weight of V0,1 is still larger than that of V1,0 as listed in Table 1. If we can reduce the number of sensor pairs with separation (0, 1), the array should have less mutual coupling. Following this idea, we propose a new planar array dubbed thermos array.


### 3.1. Thermos Array

Our proposed array consists of six uniform linear arrays (ULAs) with the interelement spacing of2d and several certain sensors at the four corners. For intuition, we first give an example of our proposed array as shown in Figure 1(c), where Nx=15, Ny=27, and the total number of sensors K=2Ny+Nx−2=67. Nx and Ny are two positive integers used for array designing. The filled circles denote the positions of sensors. As we can see, the geometry of the proposed array resembles a thermos with double-deck glass on both sides. There are six ULAs, like T,B,L1,L2,ℝ1, and ℝ2, with the same interelement spacing of 2d. Apparently, L1 and ℝ2 have the same structure, the same as L2 and ℝ1. Next, we give the construction method of our proposed thermos array using the following definition.Definition 2
(thermos array). Given two positive integersNx and Ny, the total number of sensors is K=2Ny+Nx−2. Then the thermos array is given by set v.
(11)v=B∪T∪L1∪L2∪ℝ1∪ℝ2∪F,where the position sets of six ULAs can be expressed as
(12)B=xi,0xi∈b1,T=xi,Ny−2xi∈b2,L1=−2,yiyi∈11,L2=0,yiyi∈12,ℝ1=Nx−1,yiyi∈12,ℝ2=Nx+1,yiyi∈11.
(1)
WhenNy is even, the position set of certain sensors is
(13)F=−1,Ny−2,0,Ny−2,Nx−1,Ny−2,Nx,Ny−2.
Herein, setb1, b2, l1, and l2 as
(14)Nxiseven,b1=−1+2l0≤l≤Nx2,b2=1+2l0≤l≤Nx2−2,l1=2l1≤l≤Ny2−2,l2=−1+2l1≤l≤Ny2−1,Nxisodd,b1=−1+2l0≤l≤Nx+12,b2=2+2l0≤l≤Nx−12−2,l1=2l1≤l≤Ny2−2,l2=−1+2l1≤l≤Ny2−1.
(2)
WhenNy is odd,
(15)F=−1,2,Nx,2,−1,Ny−4,0,Ny−3,Nx−1,Ny−3,Nx,Ny−4,−1,Ny−2,0,Ny−2,Nx−1,Ny−2,Nx,Ny−2.
Setb1,b2,l1, and l2 satisfying
(16)Nxiseven,b1=−1+2l∣0≤l≤Nx2,b2=1+2l∣0≤l≤Nx2−2,l1=2l∣2≤l≤Ny−12−2,l2=−1+2l∣1≤l≤Ny−12−2,Nxisodd,b1=−1+2l∣0≤l≤Nx+12,b2=2+2l∣0≤l≤Nx−12−2,l1=2l∣2≤l≤Ny−12−2,l2=−1+2l∣1≤l≤Ny−12−2.


### 3.2. Crucial Weight

Since the six ULAs have the same interelement spacing of2d, instead of d, sensor pairs with the smallest spacing d only appear at the joints of the six ULAs and the region F of certain sensors (which are the four corners of the thermos array). When the positions of certain sensors are fixed, the value of the crucial weight will be determined once Nx and Ny are given. To evaluate the values of the crucial weights, we just need to analyze the sensors at the four corners in addition. Let us take the case when Nx and Ny are both odd numbers as an example. There exist two sensor pairs contributing to V1,0, which are
(17)vi−1,Ny−2,vi,Ny−2vi=0,Nx.ForV0,1, there exist two sensor pairs
(18)vi,Ny−3,vi,Ny−2vi=0,Nx−1.Table1 gives the crucial weight comparison of the hourglass array and the thermos array. It can be observed that the value of V0,1 is reduced to the same level of V1,0. For the reason that small separations of (1, 0) and (0, 1) have the same mutual coupling according to (9), we also give the summation of small separations V1,0+V0,1 in Table 1. It can be seen that the crucial weight of the thermos array is about half of the hourglass array and thus has less mutual coupling.


### 3.3. DOF in the Difference Array

In order to utilize the second-order statistic information for DOA estimation, we analyze the DCA of the thermos array and its number of DOFs. Figure1(d) plots the DCA of the thermos array. We observe that although the resulting DCA is not a filled rectangular array, the central part of the DCA, circled out in the rectangular line box, is hole-free and its size is 2Nx+1by2Ny−3. Owing to the fact that the DOA estimation algorithm can only use the filled rectangular part (see Section 4), we care more about this rectangular part and name it as the available difference array. Therefore, we only analyze the number of DOFs in the available DCA, and it is dubbed as the available DOF.The number of available DOFs in the available difference array of the thermos array is(19)fthermos=2Nx+1⋅2Ny−3=4NxNy−6Nx+2Ny−3.The DCA of the hourglass array [11] is hole-free, and its number of DOFs is fhourglass=2Nx−1×2Ny−1.The difference betweenfthermos and fhourglass is
(20)Δf=fthermos−fhourglass=−4Nx+4Ny−4.Assuming thatΔf>0, we obtain
(21)Ny>Nx+1.Therefore, whenNx and Ny satisfy (21), the thermos array will provide more DOFs than the hourglass array. Figures 1(a) and 1(b) give one example of the hourglass array and its resulting DCA, where Nx=15 and Ny=27. From Figures 1(b) and 1(d), we can see that the number of DOFs in the DCA of the hourglass array is fhourglass=2Nx−1×2Ny−1=29×53=1537, while the available DOFs in the available DCA of our proposed array is fthermos=4×15×27−6×15+2×27−3=1581>fhourglass. If we take the virtual elements outside the rectangle line box in the DCA into consideration, the number of DOFs for our proposed array will be larger.


## 3.1. Thermos Array

Our proposed array consists of six uniform linear arrays (ULAs) with the interelement spacing of2d and several certain sensors at the four corners. For intuition, we first give an example of our proposed array as shown in Figure 1(c), where Nx=15, Ny=27, and the total number of sensors K=2Ny+Nx−2=67. Nx and Ny are two positive integers used for array designing. The filled circles denote the positions of sensors. As we can see, the geometry of the proposed array resembles a thermos with double-deck glass on both sides. There are six ULAs, like T,B,L1,L2,ℝ1, and ℝ2, with the same interelement spacing of 2d. Apparently, L1 and ℝ2 have the same structure, the same as L2 and ℝ1. Next, we give the construction method of our proposed thermos array using the following definition.Definition 2
(thermos array). Given two positive integersNx and Ny, the total number of sensors is K=2Ny+Nx−2. Then the thermos array is given by set v.
(11)v=B∪T∪L1∪L2∪ℝ1∪ℝ2∪F,where the position sets of six ULAs can be expressed as
(12)B=xi,0xi∈b1,T=xi,Ny−2xi∈b2,L1=−2,yiyi∈11,L2=0,yiyi∈12,ℝ1=Nx−1,yiyi∈12,ℝ2=Nx+1,yiyi∈11.
(1)
WhenNy is even, the position set of certain sensors is
(13)F=−1,Ny−2,0,Ny−2,Nx−1,Ny−2,Nx,Ny−2.
Herein, setb1, b2, l1, and l2 as
(14)Nxiseven,b1=−1+2l0≤l≤Nx2,b2=1+2l0≤l≤Nx2−2,l1=2l1≤l≤Ny2−2,l2=−1+2l1≤l≤Ny2−1,Nxisodd,b1=−1+2l0≤l≤Nx+12,b2=2+2l0≤l≤Nx−12−2,l1=2l1≤l≤Ny2−2,l2=−1+2l1≤l≤Ny2−1.
(2)
WhenNy is odd,
(15)F=−1,2,Nx,2,−1,Ny−4,0,Ny−3,Nx−1,Ny−3,Nx,Ny−4,−1,Ny−2,0,Ny−2,Nx−1,Ny−2,Nx,Ny−2.
Setb1,b2,l1, and l2 satisfying
(16)Nxiseven,b1=−1+2l∣0≤l≤Nx2,b2=1+2l∣0≤l≤Nx2−2,l1=2l∣2≤l≤Ny−12−2,l2=−1+2l∣1≤l≤Ny−12−2,Nxisodd,b1=−1+2l∣0≤l≤Nx+12,b2=2+2l∣0≤l≤Nx−12−2,l1=2l∣2≤l≤Ny−12−2,l2=−1+2l∣1≤l≤Ny−12−2.


## 3.2. Crucial Weight

Since the six ULAs have the same interelement spacing of2d, instead of d, sensor pairs with the smallest spacing d only appear at the joints of the six ULAs and the region F of certain sensors (which are the four corners of the thermos array). When the positions of certain sensors are fixed, the value of the crucial weight will be determined once Nx and Ny are given. To evaluate the values of the crucial weights, we just need to analyze the sensors at the four corners in addition. Let us take the case when Nx and Ny are both odd numbers as an example. There exist two sensor pairs contributing to V1,0, which are
(17)vi−1,Ny−2,vi,Ny−2vi=0,Nx.ForV0,1, there exist two sensor pairs
(18)vi,Ny−3,vi,Ny−2vi=0,Nx−1.Table1 gives the crucial weight comparison of the hourglass array and the thermos array. It can be observed that the value of V0,1 is reduced to the same level of V1,0. For the reason that small separations of (1, 0) and (0, 1) have the same mutual coupling according to (9), we also give the summation of small separations V1,0+V0,1 in Table 1. It can be seen that the crucial weight of the thermos array is about half of the hourglass array and thus has less mutual coupling.


## 3.3. DOF in the Difference Array

In order to utilize the second-order statistic information for DOA estimation, we analyze the DCA of the thermos array and its number of DOFs. Figure1(d) plots the DCA of the thermos array. We observe that although the resulting DCA is not a filled rectangular array, the central part of the DCA, circled out in the rectangular line box, is hole-free and its size is 2Nx+1by2Ny−3. Owing to the fact that the DOA estimation algorithm can only use the filled rectangular part (see Section 4), we care more about this rectangular part and name it as the available difference array. Therefore, we only analyze the number of DOFs in the available DCA, and it is dubbed as the available DOF.The number of available DOFs in the available difference array of the thermos array is(19)fthermos=2Nx+1⋅2Ny−3=4NxNy−6Nx+2Ny−3.The DCA of the hourglass array [11] is hole-free, and its number of DOFs is fhourglass=2Nx−1×2Ny−1.The difference betweenfthermos and fhourglass is
(20)Δf=fthermos−fhourglass=−4Nx+4Ny−4.Assuming thatΔf>0, we obtain
(21)Ny>Nx+1.Therefore, whenNx and Ny satisfy (21), the thermos array will provide more DOFs than the hourglass array. Figures 1(a) and 1(b) give one example of the hourglass array and its resulting DCA, where Nx=15 and Ny=27. From Figures 1(b) and 1(d), we can see that the number of DOFs in the DCA of the hourglass array is fhourglass=2Nx−1×2Ny−1=29×53=1537, while the available DOFs in the available DCA of our proposed array is fthermos=4×15×27−6×15+2×27−3=1581>fhourglass. If we take the virtual elements outside the rectangle line box in the DCA into consideration, the number of DOFs for our proposed array will be larger.


## 4. DOA Estimation Algorithm

Recalling (6), since the powers of the sources in the equivalent source vector c are real values, the rank of matrix ccH is one. Hence, these equivalent sources behave like fully coherent sources in the DCA. In order to decorrelate the equivalent sources, there are mainly two types of algorithms. One is the spatial-smoothing operation [9, 12] and the other is the sparse signal recovery [13, 14] technique. Here we employ the two-dimensional (2D) spatial-smoothing algorithm to enhance the rank.Owing to the fact that the spatial-smoothing operation needs several rotational invariance structures, we can only use the central filled rectangular part (i.e., the available DCA defined in Section3.3) for DOA estimation. We denote that the size of the virtual rectangular array is M by N, where M=2Nx+1 and N=2Ny−3. We extract the available DCA in the rectangular line box of Figure 1(d) and plot it in Figure 2. In order to apply the spatial smoothing to decorrelate the equivalent sources, we divide the available DCA into K^=M+1N+1/4 overlapping subarrays as shown in Figure 2. Each subarray has M+1/2×N+1/2 virtual sensors. We use the upper-right corner element of the subarray as its reference element, and the reference element of the kth subarray is located at position k−1−M+1/22k−1/M+1,2k−1/M+1, where ⋅·stands for round down operation.Figure 2
Subarray division in the available DCA, andM=2Nx+1 and N=2Ny−3.The equivalent data of thekth subarray can be expressed as follows:
(22)z~k=B~kc+σn2eK^−k+1,k=1,2,…,K^,where B~k=b~kϕ1,θ1,b~kϕ2,θ2,…,b~kϕQ,θQ is the array manifold matrix of the kth subarray and the corresponding steering vector is
(23)b~ϕq,θq=ejκD1θ¯,ϕ¯T,ejκD2θ¯,ϕ¯T,…,ejκDK^θ¯,ϕ¯TT,where Dk′ denotes the position of the k′th virtual sensor in the kth subarray k′=1,2,…,K^. The covariance matrix of the equivalent data z~k is given by
(24)Rk=z~kz~kH.Then taking the average ofRk over k, we can obtain the spatially smoothed matrix
(25)Rave=1K^∑k=1K^Rk.The eigenvalue decomposition ofRave yields
(26)Rave=UsΛsUsH+UnΛnUnH,where Λs=Diagλ1,λ2,…,λQ with λii=1,2,…,Q being the first Q-largest signal eigenvalues, and Us∈ℂK^×Q and Un∈ℂK^×K^−Q are the signal and noise subspace, respectively.Then we employ the MUSIC [15] algorithm to obtain estimation of the angles:
(27)ϕ^,θ^=argmaxϕ,θ1b~Hϕ,θUnUnHb~ϕ,θ.Although we only use the filled rectangular part of the DCA for DOA estimation, the following simulations will show that it can provide good performance when mutual coupling exists. We can also adopt the sparse signal recovery technique [14] to estimate the DOA, then all the virtual elements in the DCA including the parts with holes can be used. The DOA estimation performance can be improved furthermore.


## 5. Numerical Examples

In this section, we conduct simulation experiments to show the characteristics of our proposed array with the hourglass array when mutual coupling exists.


### 5.1. DOA Estimation

Firstly, we show the capability of our proposed array to estimate the DOAs of the sources when there is mutual coupling and compare its result with the hourglass array. We use the thermos array and the hourglass array depicted in Figures1(a) and 1(c), both of which have 67 physical sensors. The mutual coefficients in (9) are c1=0.19,cl=c1e−j2πl−1/5/l,l=2,2,5,22,3,…,B,B=4. We consider 72 uncorrelated narrowband sources impinging on the array with equal power, whose true DOAs are shown in Figure 3 with empty circles. We use the 2D spatial-smoothing MUSIC algorithm in Section 4 to estimate the angles, and the estimated results are also shown in Figure 3, denoted by the ×-shaped symbols. The snapshots are 1000 and the SNR is 0. It can be observed that our proposed array can identify all the sources, but the hourglass array fails for several sources.Figure 3
DOA estimations in the presence of mutual coupling of (a) the hourglass array and (b) the thermos array.
(a)
ϕ (degrees)
(b)
ϕ (degrees)


### 5.2. Estimation Performance

Next, we use Monte Carlo simulations to evaluate the average root-mean-square error (RMSE) of the estimated DOAs as a function of signal to noise ratio (SNR). The RMSE is defined asRMSE=1/Q∑q=1QEϕ^q−ϕq2+θ^q−θq2, where ϕ^q and θ^q are the estimated values of ϕq and θq, respectively. We also use both array geometries in Figures 1(a) and 1(c) and consider Q=36 narrowband uncorrelated sources with their true DOAs plotted in Figure 4. We assume that c1=0.1, the SNR is -10dB, and the snapshot number is 1000. Figure 3 gives the estimation results of one trial for the two array geometries, which shows that both arrays can identify all these 36 sources with slight differences. In order to evaluate their DOA estimation performance, we use 200 independent Monte Carlo trials to calculate their RMSEs, and the results as a function of SNR are shown in Figure 5. It is evident that the DOA estimation performance is improved with the increase of the input SNR, and our proposed array outperforms the hourglass array.Figure 4
DOA estimation for a single trial of two array geometries.Figure 5
RMSE versus SNR of two array geometries with 1000 snapshots.


## 5.1. DOA Estimation

Firstly, we show the capability of our proposed array to estimate the DOAs of the sources when there is mutual coupling and compare its result with the hourglass array. We use the thermos array and the hourglass array depicted in Figures1(a) and 1(c), both of which have 67 physical sensors. The mutual coefficients in (9) are c1=0.19,cl=c1e−j2πl−1/5/l,l=2,2,5,22,3,…,B,B=4. We consider 72 uncorrelated narrowband sources impinging on the array with equal power, whose true DOAs are shown in Figure 3 with empty circles. We use the 2D spatial-smoothing MUSIC algorithm in Section 4 to estimate the angles, and the estimated results are also shown in Figure 3, denoted by the ×-shaped symbols. The snapshots are 1000 and the SNR is 0. It can be observed that our proposed array can identify all the sources, but the hourglass array fails for several sources.Figure 3
DOA estimations in the presence of mutual coupling of (a) the hourglass array and (b) the thermos array.
(a)
ϕ (degrees)
(b)
ϕ (degrees)


## 5.2. Estimation Performance

Next, we use Monte Carlo simulations to evaluate the average root-mean-square error (RMSE) of the estimated DOAs as a function of signal to noise ratio (SNR). The RMSE is defined asRMSE=1/Q∑q=1QEϕ^q−ϕq2+θ^q−θq2, where ϕ^q and θ^q are the estimated values of ϕq and θq, respectively. We also use both array geometries in Figures 1(a) and 1(c) and consider Q=36 narrowband uncorrelated sources with their true DOAs plotted in Figure 4. We assume that c1=0.1, the SNR is -10dB, and the snapshot number is 1000. Figure 3 gives the estimation results of one trial for the two array geometries, which shows that both arrays can identify all these 36 sources with slight differences. In order to evaluate their DOA estimation performance, we use 200 independent Monte Carlo trials to calculate their RMSEs, and the results as a function of SNR are shown in Figure 5. It is evident that the DOA estimation performance is improved with the increase of the input SNR, and our proposed array outperforms the hourglass array.Figure 4
DOA estimation for a single trial of two array geometries.Figure 5
RMSE versus SNR of two array geometries with 1000 snapshots.


## 6. Conclusion

In this paper, we proposed a new two-dimensional sparse array design approach, named as the thermos array. It features having fewer sensor pairs with the smallest spacing and being more robust to mutual coupling over the hourglass array. Our proposed array has closed-form expressions for the sensor locations and the number of available DOFs. Although the resulting difference coarray of our proposed array is not hole-free, the available difference coarray can also facilitate us to perform spatial-smoothing-based DOA estimation. And the simulations have shown that our proposed array can achieve better DOA estimation performance than the hourglass array when mutual coupling exists.It is worth mentioning that the 2D spatial-smoothing algorithm in Section4 only employs a part, instead of the whole, of the difference coarray. If we adopt the sparse signal reconstruction-based DOA estimation algorithm, the whole difference coarray can be utilized. In this way, the thermos array can provide more DOFs and better DOA estimation performance.



---

*Source: 3624514-2018-04-29.xml*