Next Article in Journal
A Hydrogen Gas Sensor Based on TiO2 Nanoparticles on Alumina Substrate
Next Article in Special Issue
Biologically Inspired Hierarchical Contour Detection with Surround Modulation and Neural Connection
Previous Article in Journal
Efficient Sensor Placement Optimization for Shape Deformation Sensing of Antenna Structures with Fiber Bragg Grating Strain Sensors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Generalized L-Shaped Nested Array Concept Based on the Fourth-Order Difference Co-Array

1
School of Information and Electronics, Beijing Institute of Technology, 5 South Zhongguancun Street, Haidian District, Beijing 100081, China
2
Oil Production Technology Institute of Da Gang Oilfield Company, No. 1278, Xingfu Road, Haibin Street, Binhai New District, Tianjin 300280, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(8), 2482; https://doi.org/10.3390/s18082482
Submission received: 21 May 2018 / Revised: 18 July 2018 / Accepted: 25 July 2018 / Published: 1 August 2018

Abstract

:
In this paper, a generalized L-shaped nested array based on the fourth-order difference co-array is proposed for two-dimensional (2D) directions’ estimation. The new structure framework makes full use of the physical sensor locations to form a virtual uniform rectangular array (URA) as large as possible. As it utilizes the fourth-order difference instead of the traditional second-order difference result, this structure framework can acquire a much higher degree-of-freedom (DOF) than the existing 2D sparse arrays. The proposed structures have two advantages. One is that the subarrays can be chosen as any nested-class arrays, which makes the sparse array design more flexible. We can choose arbitrary subarray structures for DOF enhancement purposes. Another advantage is that the relative position of two subarrays can be set as any integral multiple of half wavelength. This means that two subarrays can be located as far as possible so that the relative influence between two physical subarrays can be ignored. The DOFs of several typical generalized L-shaped nested arrays (GLNAs) are compared in this paper. By setting the subarrays as different types and the relative position as a special value, a special GLNA is presented. Simulations show that GLNAs have obvious superiority in 2D direction-of-arrival estimation.

1. Introduction

In the past decades, two-dimensional (2D) planar arrays have been widely used in the field of beamforming, imaging, communication and radar [1]. The traditional configurations people usually use are the uniform rectangular arrays (URAs) and hexagonal arrays. Such kinds of arrays have physical sensors uniformly located on the rectangular grids, so that the 2D subspace estimation algorithms, such as the 2D ESPRIT method [2], can be directly applied on the received data. Although URA can obtain robust direction-of-arrival (DOA) estimation, the number of detectable sources is limited in the order of O ( N ) . The mutual coupling effect between the adjacent sensors influences the performance of URAs seriously, as well. Therefore, the sparse arrays with larger spacings between the adjacent sensor pairs have attracted much attention in the last five years.
The nested array [3] and coprime array [4,5] are important one-dimensional (1D) sparse arrays based on the concept of the second-order difference co-array (SODC), which can detect O ( N 2 ) sources with only N sensors. To further reduce the mutual coupling arising from the dense uniform linear array (ULA) in the nested array, a super nested array [6,7] is proposed to remove the dense physical sensors from the sparse ULA area. For the purpose of degree-of-freedom (DOF) enhancement and mutual coupling reduction, augmented nested arrays [8] are proposed, as well. Based on the coprime array, the DOF of which is much lower than that of the nested array, several optimized structures are proposed and verified to be effective in increasing DOF, such as the coprime array with compressed inter-element spacing (CACIS) [9], the coprime array with displaced subarrays (CADiS) and the coprime array with multi-period subarrays (CAMpS) [10].
Inspired by the 1D sparse array design, a few planar sparse arrays (PSAs) have been proposed recently for 2D sub-sampling and signal recovery. Some well-known PSA geometries include the L-shaped nested array [11,12], 2D coprime array, 2D nested array [13,14], billboard arrays, open box arrays [15] and hourglass arrays [16]. Among all these structures, the L-shaped nested array and 2D co-prime array are provided mainly to resolve the ambiguity problem and implement fast DOA estimation. The other PSA structures make use of the SODC concept to acquire the hole-free virtual planar array with DOF increasing from O ( M + N ) to O ( M N ) . Although the performance is improved significantly, there is still much space for DOF enhancement.
Recently, the fourth-order difference co-array (FODC) concept was proposed and utilized for 1D DOA estimation [17], which is verified to be effective in DOF enhancement. The new sparse arrays based on the 2qth-order cumulant [18] and especially the fourth-order cumulant [19] are the 2q-level nested array and sparse array extension with the fourth-order difference co-array enhancement based on the two-level nested array (SAFOE-NA). Another algorithm that can generate FODC makes use of the quasi-stationary (QS) signal properties to acquire the virtual FODC array. As many speech and audio signals have been proven to be quasi-stationary, such an FODC construction scheme has great practical significance. It is verified that by exploiting the FODC methods, the DOF of the 1D sparse arrays can be increased significantly. Nowadays, the four-order difference concept is only used in 1D sparse array design. For the 2D sparse array, there has been no study conducted so far.
In this paper, the generalized L-shaped nested arrays (GLNAs) based on the FODC concept are proposed for DOF enhancement. The new 2D configuration is a generalized framework including two nested-class subarrays positioned at an arbitrary distance. The nested-class arrays guarantee that a cross-shaped structure can be formed in SODC. Therefore, a virtual URA with the center at the origin is constructed in the FODC. As the FODC concept instead of SODC is utilized in the new GLNA construction, the number of detectable sources will be increased greatly. Since the subarrays can be chosen as any kind of existing nested arrays, more choices are provided in sparse array design. By setting the subarrays with a large 1D virtual consecutive range and sparse array structure, GLNA can achieve a high number of DOFs. Another advantage of GLNA is that the relative 2D position of the subarrays can be chosen as any integer times the minimum spacial interval d. Since relative distance can be set as large as possible, the influence between two physical subarrays can be greatly reduced. To make complete comparisons among all GLNAs, we provide the expressions of URA size and list the values of several different GLNAs. Furthermore, a special GLNA with the highest number of DOFs is presented and analyzed.
This paper is organized as follows. The QS signal model and FODC concept are presented in Section 2. The proposed generalized L-shaped nested array is proposed in Section 3. The DOF characteristics are summarized. Several GLNAs are compared to other SODC-based PSAs in the aspect of DOF. In Section 4, a special GLNA (sGLNA) is developed and analyzed, showing its great performance improvement in DOF enhancement. Simulation results are shown in Section 5, and the conclusions are drawn in Section 6. Notations: We denote vectors and matrices by using lower-case and upper-case bold characters, respectively. Specifically, ( . ) * denotes complex conjugate, whereas ( . ) T and ( . ) H respectively denote the transpose and conjugate transpose of a matrix or vector. v e c ( . ) implies the vectorization operator, which turns a matrix into a column vector. E ( . ) is the statistical expectation operator. | . | 2 is the Frobenius norm of the 2D vector in it.

2. Signal Model

Consider D uncorrelated narrowband QSsignals s i ( t ) , i = 1 , , D from direction ( θ i , ϕ i ) impinging on an N-sensor planar array with physical sensors located at S = [ n 1 ; n 2 ; ; n N ] d . n j = ( n j x , n j y ) is the j-th normalized sensor location for j = 1 , 2 , , N . d is the unit spacing between two sensors, which is also the length unit of the x and y Cartesian axes, where the arrays are placed. The value of d is set to be half the minimal wavelength λ . For convenience, d is omitted in the following discussions. Assume that the source signals s i ( t ) , i = 1 , , D are wide-sense quasi-stationary within the frame length P ˜ . Denote s i ( p , k ) as the i-th source signal at the k-th snapshot within the p-th frame. Then, the power σ p , i 2 = E { s i ( p , k ) s i * ( p , k ) } for k { p · P ˜ , p · P ˜ + 1 , , ( p + 1 ) · P ˜ 1 } can be approximated by [20]:
σ p , i 2 1 P ˜ k = p · P ˜ ( p + 1 ) · P ˜ 1 s i ( p , k ) s i * ( p , k )
Then, the received signals from N sensors at the k-th snapshot within the p-th frame can be represented as:
x S ( p , k ) = i = 1 D a S ( θ i , ϕ i ) s i ( p , k ) + w S ( p , k ) ,
where the element of a S ( θ i , ϕ i ) for the i-th source impinging on the grid ( n j x , n j y ) is e j 2 π d / λ ( s i n θ i c o s ϕ i n j x + s i n θ i s i n ϕ i n j y ) ; w S ( p , k ) is the wide-sense stationary noise vector.
By taking all the K snapshots in the p-th frame into consideration, the correlation matrix is given by:
R x x ( p ) = E { x S ( p , k ) x S ( p , k ) H } = A R s s ( p ) A H + σ n 2 I ,
where A = [ a S ( θ ¯ 1 , ϕ ¯ 1 ) , , a S ( θ ¯ D , ϕ ¯ D ) ] is the array manifold; R s s ( p ) = d i a g { σ p , 1 2 , σ p , 2 2 , , σ p , D 2 } is the source covariance matrix; σ n 2 is the noise power; I C N × N is the identity matrix. Vectorizing R x x ( p ) yields:
z ( p ) = vec { R xx ( p ) } = ( A * A ) s ˜ ( p ) + l N ( p ) ,
where s ˜ ( p ) = [ σ p , 1 2 , σ p , 2 2 , , σ p , D 2 ] T can be viewed as the virtual source vector; l N ( p ) = v e c { σ 2 ( p ) I } C N 2 × 1 is the noise vector. Assume two sensor location matrices as S 1 = [ n 11 ; n 12 ; ; n 1 N 1 ] and S 2 = [ n 21 ; n 22 ; ; n 2 N 2 ] , the numbers of sensors if which are N 1 and N 2 , respectively. Define the sensor location difference operation as D i f f ( S 1 S 2 ) = [ ( n 11 n 21 ) ; ( n 11 n 22 ) ; ; ( n 11 n 2 N 2 ) ; ( n 21 n 21 ) ; ( n 21 n 22 ) ; ; ( n 21 n 2 N 2 ) ; ; ( n 1 N 1 n 21 ) ; ( n 1 N 1 n 22 ) ; ; ( n 1 N 1 n 2 N 2 ) ] . Then, the Kronecker product ( A * A ) is an equivalent manifold matrix, whose virtual sensors are located at C ¯ = D i f f ( S S ) . It is actually an SODC with the virtual elements at the difference positions of the physical sensor locations.
As the assumption given by [20], σ p , i 2 , i = 1 , , D are wide-sense stationary and uncorrelated with each other. Denote the self-correlation of σ p , i 2 as σ ˜ i 2 = E { ( σ p , i 2 E { σ p , i 2 } ) 2 } . Then, SODC can be conducted again to form FODC after the zero-mean equivalent impinging sources have been generated.
Denote the expectation over P frame as E [ z ( p ) ] . We have:
E [ z ( p ) ] = ( A * A ) E [ s ˜ ( p ) ] + l N ( p ) ,
Subtracting the above result from z ( p ) , we can obtain:
z ¯ ( p ) = z ( p ) E [ z ( p ) ] A ˜ s ¯ ( p ) ,
where A ˜ = ( A * A ) with the i-th column as a ( θ ¯ i , ϕ ¯ i ) and s ¯ ( p ) = s ˜ ( p ) E [ s ˜ ( p ) ] .
As s ¯ ( p ) can be viewed as a zero-mean virtual signal vector, the difference co-array concept can be applied again. Similarly to (3), the correlation matrix of the virtual received signals z ¯ ( p ) over P frames can be represented as:
R z z = E { z ¯ ( p ) z ¯ H ( p ) } = A ˜ R s ¯ s ¯ A ˜ H ,
where R s ¯ s ¯ = d i a g { σ ˜ 1 2 , , σ ˜ D 2 } .
Performing the vectorization again, one has:
y = vec { R zz } = ( A ˜ * A ˜ ) s ^ ,
where ( A ˜ * A ˜ ) is the virtual manifold matrix with the virtual sensors located at C ^ = D i f f ( C ¯ C ¯ ); s ^ = v e c ( R s ¯ s ¯ ) = [ σ ˜ 1 2 , , σ ˜ D 2 ] T , is the equivalent source signal vector. Perform the 2D spatial smoothing algorithm [21] and the 2D subspace algorithms such as 2D unitary-ESPRIT [22] and 2D MUSIC [23] directly on y ; 2D DOA estimations can be obtained.

3. The Proposed L-Shaped Nested Arrays

Different from the existing L-shaped nested array [11,12] and other second-order difference co-array-based structures [13,14,16], the new 2D sparse arrays proposed in this paper are based on the fourth-order difference co-array concept. That means as long as a consecutive URA can be extracted from FODC, the 2D spatial smoothing algorithm can be made use of directly on it. Therefore, in this paper, we focus on designing a kind of 2D sparse array that can generate a URA area as large as possible in FODC.
Since the 2D sparse arrays we will develop can be grouped into a unified framework, we first introduce the generalized concept and then compare the DOF of several sub-structures in the framework. For convenience, in the following sections, we use the normalized sensor locations to represent the actual physical positions, i.e., d is omitted in the following discussions.

3.1. Generalized L-Shaped Nested Array

Note that d is omitted in the following discussions for convenience. Then, the generalized L-shaped nested array can be described by the definition shown below.
Definition 1.
(Generalized L-shaped nested array (GLNA)): Given two 1D nested-class arrays described by the matrix C h and C v , where C h = [ s h 1 ; s h 2 ; ; s h N x ] and C v = [ s v 1 ; s v 2 ; ; s v N y ] , the sensor numbers of these two arrays are N x and N y , respectively. Assume N = N x + N y , then a generalized L-shaped nested array is an N-sensor 2D planar sparse array located on the grids, which is denoted as follows:
S GLNA = [ H ; V ]
where H = [ ( s h 1 , 0 ) ; ( s h 2 , 0 ) ; ; ( s h N x , 0 ) ] and V = [ ( L x , L y + s v 1 ) ; ( L x , L y + s v 2 ) ; ; ( L x , L y + s v N y ) ] . L x and L y can be any integers.
It can be summarized from the definition that a GLNA consists of two 1D nested-class arrays located in the horizontal and vertical direction, respectively. H and V are the sensor locations of two subarrays, which can be any kind of aforementioned nested array, such as the classic nested arrays (CNAs), super nested arrays (SNAs) and augmented nested arrays (ANAs). The relative position parameter pair in this definition, ( L x , L y ) , of the two subarrays is arbitrary. To have a better idea of the GLNA configuration, we take an example by setting CNAs as subarrays. As shown in Figure 1, two 1D CNAs with C h = C v = [ 1 ; 2 ; 3 ; 4 ; 5 ; 10 ; 15 ; 20 ; 25 ] are distributed on the x-axis, and line x = 6 . The reference points of the two subarrays are ( 0 , 0 ) and ( 6 , 6 ) with the relative displacement as ( L x , L y ) = ( 6 , 6 ) . Therefore, the vertical sensors are located on V = [ ( 6 , 7 ) ; ( 6 , 8 ) ; ( 6 , 9 ) ; ( 6 , 10 ) ; ( 6 , 11 ) ; ( 6 , 16 ) ; ( 6 , 21 ) ; ( 6 , 26 ) ; ( 6 , 31 ) ] . The horizontal subarray lies on positions H = [ ( 1 , 0 ) ; ( 2 , 0 ) ; ( 3 , 0 ) ; ( 4 , 0 ) ; ( 5 , 0 ) ; ( 10 , 0 ) ; ( 15 , 0 ) ; ( 20 , 0 ) ; ( 25 , 0 ) ] . The total number of sensors in this GLNA is 18. Note that in the following structure, as length unit d is defined to be the half-wavelength, it is clear that the units of the x and y Cartesian axes of the array geometry are also in half wavelengths.
Denote the second-order and fourth-order self-difference of any sparse array S as C ¯ S = D i f f ( S S ) and C ^ S = D i f f ( C ¯ S C ¯ S ) , respectively. The SODC and FODC of GLNA have the following conclusions.
Property 1.
For two 1D subarrays C h and C v of GLNA defined by (9), their second-order self-differences are denoted as C ¯ C h and C ¯ C v . If there exists consecutive lags in their self-difference arrays, i.e., [ L h ; ( L h + 1 ) ; ; ( L h 1 ) ; L h ] is a part of C ¯ C h and [ L v ; ( L v + 1 ) ; ; ( L v 1 ) ; L v ] is a part of C ¯ C v , then the SODC of GLNA must include a cross-shaped structure distributed at S d = [ S h ; S v ] with S h = [ ( L h , 0 ) ; ( L h + 1 , 0 ) ; ; ( L h 1 , 0 ) ; ( L h , 0 ) ] and S v = [ ( 0 , L v ) ; ( 0 , L v + 1 ) ; ; ( 0 , L v 1 ) ; ( 0 , L v ) ] . The FODC of GLNA contains a central URA with a minimum aperture of ( 2 L h + 1 ) × ( 2 L v + 1 ) .
Proof. 
First, let us consider the SODC of GLNA. According to the expression of GLNA in (9), SODC can be written as:
C ¯ S = [ D h ; D v ; D h v ; D v h ] ,
with D h = D i f f ( H H ) and D v = D i f f ( V V ) as the self-differences of subarrays H and V . [ D h v , D v h ] T is the cross-difference. As H is a nested-class array located on the x-axis, its self-difference is still a linear array with the form of D h = [ C ¯ C h , 0 ] . Since the assumption that there exist consecutive lags in the self-difference arrays of each subarray, we conclude that S h = [ ( L h , 0 ) ; ( L h + 1 , 0 ) ; ; ( L h 1 , 0 ) ; ( L h , 0 ) ] is a part of matrix D h . Similarly, the self-difference result of V is a linear structure at the y-axis located at D v = [ 0 , C ¯ C v ] . Its array geometry also has a consecutive segment included in the center, i.e., the matrix S v = [ ( 0 , L v ) ; ( 0 , L v + 1 ) ; ; ( 0 , L v 1 ) ; ( 0 , L v ) ] is a part of D v . Then, we turn to D v h . From the definition of H and V , it is not difficult to obtain that D v h is a matrix, the elements of which can be denoted as ( s h , i + L x , s v , j + L y ) . Here, s h , i and s v , j are the i-th and the j-th elements of the matrix C h and C v , respectively. Such a matrix is the position of the cross-difference array, which is in fact several copies of the vertical subarray with the i-th reference points positioned at ( s h , i + L x , L y ) where s h , i is the i-th element of the matrix C h . Another cross-difference matrix D h v is the reverse of D v h . Therefore, it generates the same number of vertical subarrays symmetrically. In conclusion, no matter what type of subarray in GLNA is, as long as each of them can generate a consecutive segment by making the self-difference, it is definite that SODC will contain a cross-shaped structure at the x-axis and y-axis, i.e., matrix [ S v , S h ] is a part of matrix C ¯ S .
Then, we extract the cross structure out and make the difference again to generate part of FODC. Let S d = [ S v ; S h ] . Again, we have matrix D i f f ( S d S d ) = [ D i f f ( S v S v ) ; D i f f ( S h S h ) ; D i f f ( S v S h ) ; D i f f ( S h S v ) ] is a part of C ^ S = D i f f ( C ¯ S C ¯ S ) . From the form of S h and S v , we deduce that D i f f ( S h S h ) = [ ( 2 L h , 0 ) ; ( 2 L h + 1 , 0 ) ; ; ( 2 L h 1 , 0 ) ; ( 2 L h , 0 ) ] , D i f f ( S v S v ) = [ ( 0 , 2 L v ) ; ( 0 , 2 L v + 1 ) ; ; ( 0 , 2 L v 1 ) ; ( 0 , 2 L v ) ] and D i f f ( S v S h ) = D i f f ( S h S v ) = [ ( L h , L v ) ; ( L h , L v + 1 ) ; ; ( L h , L v ) ; ( L h + 1 , L v ) ; ; ( L h , L v ) ] . It contains a URA with the size of ( 2 L h + 1 ) × ( 2 L v + 1 ) . ☐
Figure 2 illustrates the SODC and FODC of GLNA in Figure 1. As the two subarrays in GLNA are the same classic nested arrays with the self-difference being hole-free, i.e., C ¯ C h = C ¯ C v = [ L h ; L h + 1 ; ; L h ] with L h = 24 , the SODC in Figure 2a possesses a cross-shaped structure with two 49-element ULAs at the x-axis and y-axis in the range of [ 24 , 24 ] . The top left and lower right parts are the cross-difference results, which are in fact copies of the original vertical subarray. The starting points of the top left copies are (−7,7), (−8,7), (−9,7), (−10,7), (−11,7), (−16,7), (−21,7), (−26,7) and (−31,7). The FODC configuration is shown in Figure 2b. It is obvious that a URA with the size of 49 × 49 is located in the center. The URA part can be extracted to perform the 2D spatial smoothing algorithm on it. In this case, at most 49 × 49 far field uncorrelated sources can be detected, which is much more than that of the existing PSAs based on the SODC concept.
Remark 1.
According to the above analysis, the new structures possess two advantages. First, any 1D subarray geometry having consecutive segments in its self-difference co-array can be adopted. This means that the types of the subarrays can be different, as long as each can generate an equivalent ULA in the self-difference result. This enables us to choose any existing nested-class arrays as needed for further DOF enhancement or mutual coupling reduction. Second, the relative position of the two subarrays can be flexibly set. This is because the cross-shaped section in SODC is generated from the self-difference of subarrays. No matter where the subarrays are located, the difference result would have no change. This inspires us to separate the two subarrays with the parameters ( L x , L y ) as large as possible, if lower mutual coupling is desired in the GLNA design. When the relative locations grow large enough, the influence between two subarrays can be ignored. Note that in this case, qualitative analysis instead of quantitative analysis for mutual coupling effect is considered. Therefore, we assume that the mutual coupling between subarray can be ignored because of the long relative distance and the mutual coupling inside each subarray is fully calibrated. In the following subsection, we would only compare the DOF of different GLNAs with different types of subarrays.

3.2. DOF Comparisons

For convenience, in this subsection, we set the two subarrays in GLNA as the same type, e.g., the classic nested array [3] and augmented nested array [8]. The minimum DOF that GLNA can obtain is calculated. As a contrast, the DOF of two PSAs, i.e., the hourglass array [16] and 2D nested array [13,14], are compared, as well. Definitions of these structures can be seen in the relative references.
First, assume that subarrays in GLNA-CNA are both classic nested arrays. The number of physical sensors in the dense part and sparse part of the horizontal CNA is N x d and N x s . The corresponding values for the vertical CNA are N y d and N y s . Take the horizontal CNA as an example. It is the union of the matrices C i n n e r = [ 1 ; 2 ; . . . ; N x d ] and C o u t e r = [ ( N x d + 1 ) ; 2 ( N x d + 1 ) ; ; N x s ( N x d + 1 ) ] . The CNA can be denoted as:
C C N A = [ C i n n e r ; C o u t e r ] = [ 1 ; 2 ; . . . ; N x d ; ( N x d + 1 ) ; 2 ( N x d + 1 ) ; ; N x s ( N x d + 1 ) ] .
Its SODC is a ULA with 2 N x d ( N x s + 1 ) 1 elements positioned at [ L h , L h ] = [ 1 N x d ( N x s + 1 ) , N x d ( N x s + 1 ) 1 ] . Therefore, we have L h = N x d ( N x s + 1 ) 1 . Similarly, for the vertical CNA, the parameter L v = N y d ( N y s + 1 ) 1 . According to Property 1, the central URA has the minimum aperture of ( 2 L h + 1 ) × ( 2 L v + 1 ) . Therefore, the minimum DOF of GLNA-CNA is:
DOF GLNA CNA = [ 2 N x d ( N x s + 1 ) 1 ] [ 2 N y d ( N y s + 1 ) 1 ] .
For the GLNA with augmented nested arrays as subarrays (GLNA-ANA), we use ANAI-2 [8] as the subarrays. For ANAI-2, its definition is:
C A N A I 2 = [ L 1 ; M 1 ; R 1 ]
with L 1 = [ 1 ; 2 ; 4 ; ; 2 k 1 + 2 ] , M 1 = [ 0 ; ( N d + 1 ) ; 2 ( N d + 1 ) ; , ( N s + 1 ) ( N d + 1 ) ] , R 1 = [ ( N d + 1 ) N s + N d 2 ; ( N d + 1 ) N s + N d 4 ; ; ( N d + 1 ) N 2 + N d 2 k 2 ; ; ( N d + 1 ) N s ] . Here, N d and N s are the sensor numbers of the dense part and sparse part of its original classic nested array, respectively. k 1 and k 2 are the parameters with k 1 [ 0 , R ( N 1 1 ) / 2 ] and k 2 [ 1 , R N 1 / 2 ] . R α round the value to the nearest integer, where R α α . In this case, ANAI-2 instead of other kinds of augmented nested array is chosen, as it has similar parameters as other nested arrays, so that it is easier for us to compare.
Assume the array parameters as ( N x d , N x s ) , which denote ( N d , N s ) on the x-axis, and ( N y d , N y s ) m which is ( N d , N s ) in the y-axis direction. Based on the properties in [8], it is not difficult to find that the consecutive parameters in SODC are L h = ( N x d + 1 ) N x s + N x d 3 and L v = ( N y d + 1 ) N y s + N y d 3 . Therefore, the minimum DOF of GLNA-ANA is:
DOF GLNA ANA = [ 2 ( N x d + 1 ) N x s + 2 N x d 5 ] [ 2 ( N y d + 1 ) N y s + 2 N y d 5 ] .
For fair comparisons, we consider the 2D nested array (2D-NA) and perform the fourth-order difference algorithm on it for DOF analysis. Denote N x d and N x s as the number of dense and sparse elements on the x-axis. N d y and N s y are the corresponding parameters on the y-axis. Then, the dense and sparse URAs in 2D-NA have N x d N y d and N x s N y s physical sensors, respectively. According to the results in [13,14], the SODC of 2D-NA is a URA with the size of N x d N x s × ( 2 N y d N y s 1 ) . Perform the difference again. Then, the URA region in FODC has the aperture of:
DOF 2 D NA = ( 2 N x d N x s 1 ) × ( 4 N y d N y s 3 ) .
For the hourglass array (HA), the related definition and properties are provided in [16]. Consider N x sensors located in the horizontal direction. 2 N y 2 sensors are distributed on the vertical lines. Then, its SODC is proven to be a URA with the size of ( 2 N x 1 ) ( 2 N y 1 ) . The FODC is definitely a URA with the DOF of:
DOF HA = ( 4 N x 3 ) ( 4 N y 3 ) .
For better understanding, we take the aforementioned PSAs with the same number of sensors which is 18 as an example. The minimum DOFs of GLNAs and other two PSAs are summarized in Table 1. The relative position of two subarrays in GLNAs is set as ( L x , L y ) = ( 6 , 6 ) . As can be seen, all GLNAs can acquire a much higher number of DOF than the other two PSAs. This is because 2D-NA and HA are proposed for second-order difference design. When the fourth-order difference concept is applied on the two PSAs, much difference redundancy exists. As for the two kinds of GLNAs, the GLNA-ANA has the largest URA in FODC, since ANAI-2 has more consecutive lags in SODC than CNA. Meanwhile, GLNA-CNA is a little worse than the GLNA-ANA. The better the 1D sparse array we choose, the higher the DOF we obtain.

4. A Special GLNA Configuration with Enhanced DOF

For convenient analysis, GLNAs discussed in the above subsection are chosen as PSAs with the same type of subarrays. However, for the generalized structure defined by (9), the two subarrays are allowed to be totally different types. In this subsection, we will provide a special GLNA with a nested array and an augmented nested array, ANAI-2, as subarrays. By setting their relative position ( L x , L y ) as specific values, the special GLNA can acquire a much higher DOF than the aforementioned GLNAs. We use sGLNA to represent this kind of GLNA in this paper.
First, observe the SODC structure of GLNA-CNA in Figure 2a. The virtual ULAs D h and D v on the x-axis and y-axis are generated from the self-difference of subarrays. As long as the nested-class arrays are chosen, the cross-shaped ULAs can be obtained. For the sparse elements in the second and fourth quadrants, they come from the cross-difference D v h and D h v . The location of this part is determined by the relative position ( L x , L y ) . Note that in the definition of GLNA in the above section, L x and L y are set to be arbitrary integers. However, if ( L x , L y ) is chosen as special values, the cross-difference parts will probably move to make one of the consecutive lags in this area connected to the virtual ULA at the y-axis. As a result, the length of the cross-shaped structure in the y-direction will be extended. The size of FODC will be enlarged, as well.
Motivated by the above strategy, we use one of the four augmented nested array structures called ANAI-2 [8] as the subarray in the x-direction and another classic nested array in the y-direction to form sGLNA. Here, ANAI-2 instead of other kinds of augmented nested arrays is chosen as it has similar parameters as other nested arrays, so that it is easier for us to compare. In order to define this special GLNA, ( L x , L y ) are chosen as special values. L y = N y s ( N y d + 1 ) 1 , and let the value of L x be equal to the abscissa of subarray H , i.e., L x is an element of S h Then, we have the following property.
Property 2.
For sGLNA with ANAI-2 and classic nested array being the horizontal and vertical subarray having parameters N x d , N x s and N y d , N y s , the cross-shaped structure in SODC is located at [ S v ; S h ] with S v = [ ( 0 , L v ) ; ( 0 , L v + 1 ) ; ; ( 0 , L v 1 ) ; ( 0 , L v ) ] where L v = N y s ( N y d + 1 ) + N y d 1 and S h = [ ( L h , 0 ) ; ( L h + 1 , 0 ) ; ; ( L h 1 , 0 ) ; ( L h , 0 ) ] where L h = ( N x d + 1 ) N x s + N x d 3 . The size of FODC is [ 2 ( N x d + 1 ) N x s + 2 N x d 5 ] × [ 2 N y s ( N y d + 1 ) + 2 N y d 1 ] .
Proof. 
Since the horizontal subarray is an ANAI-2, the 1D physical sensors’ locations are given according to the definition of augmented nested arrays in [8] and Equation (13). According to [8], its self-difference is a hole-free ULA. Therefore, D h = S h = [ ( L h , 0 ) ; ( L h + 1 , 0 ) ; ; ( L h 1 , 0 ) ; ( L h , 0 ) ] with L h = ( N x d + 1 ) N x s + N x d 3 .
For the vertical subarray, as it is a nested array, the 1D sensor location is C v = [ L 1 ; R 1 ] with L 1 = [ 1 ; 2 ; ; N y d ] and R 1 = [ ( N y d + 1 ) ; 2 ( N y d + 1 ) ; ; N y s ( N y d + 1 ) ] . Its self-difference is hole-free, as well, with the one-side consecutive length as L v 1 = N y s ( N y d + 1 ) 1 .
Now, let us see the second-order cross-difference D v h . From the analysis in the proof of Property 1, we have the D v h is a matrix, the elements of which can be denoted as ( s h , i + L x , s v , j + L y ) . Here, s h , i and s v , j are the i-th and the j-th elements of the matrices C h and C v , respectively. Since for sGLNA, L x is an element of C h and L y = N y s ( N y d + 1 ) 1 , it can be deduced that L x must be equal to one of the horizontal coordinates of ANAI-2. Therefore, there must be cross difference array in SODC, the j-th element of the matrix of which can be denoted as ( 0 , s v , j + L y ) , which is located on the y-axis. As a ULA part exists in the C v , the consecutive lags in the above matrix will become [ ( 0 , N y s ( N y d + 1 ) ) ; ( 0 , N y s ( N y d + 1 ) + 1 ) ; ( 0 , N y s ( N y d + 1 ) + 2 ) ; ; ( 0 , N y s ( N y d + 1 ) + N y d 1 ) ] . As a result, the consecutive length is L v 2 = N y d . Combined with the self-difference result, the consecutive length in the vertical direction is L v = L v 1 + L v 2 = N y s ( N y d + 1 ) 1 + N y d .
As for FODC, URA is in the part ( 2 L v + 1 ) × ( 2 L h + 1 ) . Therefore, the size of URA is [ 2 ( N x d + 1 ) N x s + 2 N x d 5 ] × [ 2 N y s ( N y d + 1 ) + 2 N y d 1 ] . ☐
Figure 3 illustrates an sGLNA with the horizontal subarray as an ANAI-2 having parameters as ( N x d , N x s ) = ( 4 , 5 ) . The vertical subarray is a nested array with parameters as ( N y d , N d s ) = ( 4 , 5 ) . The relative position is ( L x , L y ) = ( 15 , 24 ) . As can be seen, the subarray 2 is aligned with one of the sensors in Subarray 1. Such an arrangement guarantees that one copy of Subarray 2 is located at the y-axis in SODC, which can connect to the virtual ULA from 1 to 24, generated from the self-difference of Subarray 2. As such, the consecutive lags extend from the original range of [−24,24] to [−29,29]. Since Subarray 1 is an ANAI-2, its self-difference is hole-free, located in the range of [ 26 , 26 ] . We obtain that the URA in FODC has dimensions of ( 2 × 26 + 1 ) ( 2 × 29 + 1 ) = 3127 .
For comparison, the DOF of sGLNA is listed in Table 1. We conclude that sGLNA has the highest number of DOFs among all GLNAs. That is because one of subarrays in sGLNA is ANAI-2, which possesses the best performance in DOF enhancement. Besides, another subarray is a nested array. Although its DOF is less than that of ANAI-2, it has a dense ULA part, which can be used to enlarge the consecutive range in the cross-shaped structure of SODC. The values of L x and L y make sure that the copy of the nested array can be connected to the virtual ULA at the y-axis in SODC. The extended consecutive range leads to the great enhancement of DOF.

5. Numerical Experiments

In this section, simulation results are presented to compare the performance of GLNA-CNA, GLNA-ANA and sGLNA. The hourglass array and 2D nested array are also simulated as a contrast. In the experiments, we assume the number of sensors to be 18. The parameters for GLNA-CNA, GLNA-SNA, GLNA-ANA and sGLNA are ( N x d , N x s , N y d , N y s ) = ( 4 , 5 , 4 , 5 ) . The relative position is defined by L x and L y . ( L x , L y ) = ( 5 , 29 ) . For the 2D nested array and the hourglass array, the configurations are quite different from GLNAs. To make sure the total number of sensors is the same as that of GLNAs, we choose the parameters as ( N d x , N d y , N s x , N s y ) = ( 3 , 3 , 3 , 3 ) and ( N x , N y ) = ( 6 , 7 ) , respectively. The QS signals with P = 10 frames are used.
In the first experiment, D = 20 uncorrelated sources with the directions ( θ , ϕ ) are distributed at (−34 ,45 ), (−26 ,63 ), (26 ,−63 ), (34 ,−45 ), (−26 ,26 ), (−16 ,45 ), (16 ,−45 ), (26 ,−26 ), (0 ,−23 ), (0 ,−11 ), (0 ,11 ), (0 ,23 ), (−26 ,−26 ), (−16 ,−45 ), (16 ,45 ), (26 ,26 ), (−34 ,−45 ), (−26 ,−63 ), (26 ,63 ), (34 ,45 ). They are assumed to impinge on 2D-NA, HA, GLNA-CNA, GLNA-ANA and sGLNA. These source angles are selected randomly. The SNR is 20 dB. P ˜ = 100 snapshots in each frame, in total K = P ˜ × P = 1000 , are sampled for R x x ( p ) construction. The DOA estimation results are examined and compared in Figure 4. For each configuration, the root mean square error (RMSE) is defined as:
RMSE = [ ( 1 / ( D N ˜ ) n = 1 N ˜ i = 1 D ( θ i n ^ θ i n ) 2 + ( ϕ i n ^ ϕ i n ) 2 ] 1 / 2
where ( θ ^ , ϕ ^ ) and ( θ , ϕ ) are the estimated DOAs and the true DOAs; N ˜ is the number of trials. In this experiment, N ˜ = 30 . As can be concluded from Figure 4a, 2D-NA fails to estimate some directions with RMSE = 8.0128 . The situation becomes better for HA. However, there are still obvious estimation errors with RMSE = 2.5737 . In contrast, as Figure 4c–e illustrate, all GLNAs can estimate DOAs correctly. The sGLNA has the lowest RMSE = 0.5441 , compared to the GLNA-CNA with RMSE = 1.1827 and GLNA-ANA with RMSE = 0.5488 .
In the second experiment, we compare the RMSE results of the five kinds of PSAs mentioned above at different SNRs and snapshots. D = 20 sources distributed in the same normalized directions as the first experiment are used in this experiment. Figure 5a depicts the RMSE versus SNR results. N ˜ = 100 independent trials are conducted for each SNR. The snapshot number is set as K = 100 × 100 . The SNR varies from −10 dB to 30 dB. Figure 5b provides the RMSE comparisons with the snapshot number varying from 50 × 100 to 350 × 100 . The SNR is 10 dB in this case. As can be seen, for both figures, that sGLNA has the lowest estimation error among all five sparse arrays. This is because sGLNA can acquire the highest number of DOF. For other PSAs, GLNA-ANA has a slightly worse performance than sGLNA, followed by GLNA-CNA, 2D-NA and HA. Therefore, GLNAs can achieve better DOA performance than the second-order-based PSAs. This result coincides with the conclusion in Table 1. The reason why sGLNA is less robust than ALNA-ANA in Figure 5b is that the robustness performance of one PSA is determined by the distribution of its difference redundancy. Note that for sGLNA, it has 12 redundancy functions not less than two, and for GLNA-ANA, it has 14 redundancy functions not less than two. ALNA-ANA has more redundancy. Therefore, it is more robust than sGLNA.

6. Conclusions

In this paper, a generalized L-shaped nested array framework is proposed. The expressions of the dimensions of the cross-shaped structure in SODC and URA size in FODC are provided and verified.
Several GLNAs, such as GLNA-CNA and GLNA-ANA, are analyzed and compared in the DOF aspects. For DOF improvement, a special GLNA configuration is proposed and compared to other PSAs. In conclusion, GLNAs enjoy closed-form sensor locations and large URA sections in FODC. They can acquire much higher DOFs than any other SODC-based PSAs by choosing the subarrays and setting the distance between them properly. Meanwhile, note that for the GLNA framework, as a sparse array structure can be constructed and a low difference redundancy can be achieved, the mutual coupling effect is reduced and can be quantitatively analyzed in the future. Several numerical experiments indicate that these structures perform better than the existing PSAs.

Author Contributions

The main idea was proposed by L.Z., S.R. and X.W. conceived of and designed the experiments. X.L. and S.R. performed the experiments. L.Z. and S.R. analyzed the data. G.R. and X.W. contributed materials. L.Z., S.R. wrote the paper. All the authors participated in amending the original manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Trees, H.V. Optimum Array Processing : Part IV of Detection, Estimation And Modulation Theory; Wiley-Intersci.: New York, NY, USA, 2002; pp. 293–303. [Google Scholar]
  2. Roy, R.; Kailath, T. ESPRIT-estimation of signal parameters via rotational invariance techniques. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 984–995. [Google Scholar] [CrossRef]
  3. Pal, P.; Vaidyanathan, P.P. Nested Arrays: A Novel Approach to Array Processing With Enhanced Degrees of Freedom. IEEE Trans. Signal Process. 2010, 58, 4167–4181. [Google Scholar] [CrossRef] [Green Version]
  4. Vaidyanathan, P.P.; Pal, P. Sparse Sensing With Co-Prime Samplers and Arrays. IEEE Trans. Signal Process. 2011, 59, 573–586. [Google Scholar] [CrossRef]
  5. Zhang, Y.D.; Amin, M.G.; Himed, B. Sparsity-based DOA estimation using co-prime arrays. In Proceedings of the 2013 IEEE International Conference on Acoustics, Speech and Signal Processing, Vancouver, BC, Canada, 26–31 May 2013; pp. 3967–3971. [Google Scholar]
  6. Liu, C.L.; Vaidyanathan, P.P. Super Nested Arrays: Linear Sparse Arrays With Reduced Mutual Coupling Part I: Fundamentals. IEEE Trans. Signal Process. 2016, 64, 3997–4012. [Google Scholar] [CrossRef]
  7. Liu, C.L.; Vaidyanathan, P.P. Super Nested Arrays: Linear Sparse Arrays With Reduced Mutual Coupling Part II: High-Order Extensions. IEEE Trans. Signal Process. 2016, 64, 4203–4217. [Google Scholar] [CrossRef]
  8. Liu, J.; Zhang, Y.; Lu, Y.; Ren, S.; Cao, S. Augmented Nested Arrays With Enhanced DOF and Reduced Mutual Coupling. IEEE Trans. Signal Process. 2017, 65, 5549–5563. [Google Scholar] [CrossRef]
  9. Qin, S.; Zhang, Y.D.; Amin, M.G. Generalized Coprime Array Configurations for Direction-of-Arrival Estimation. IEEE Trans. Signal Process. 2015, 63, 1377–1390. [Google Scholar] [CrossRef]
  10. Wang, W.; Ren, S.; Chen, Z. Unified coprime array with multi-period subarrays for direction-of-arrival estimation. Dig. Signal Process. 2018, 74, 30–42. [Google Scholar] [CrossRef]
  11. Dong, Y.Y.; Dong, C.X.; Zhu, Y.T.; Zhao, G.Q.; Liu, S.Y. Two-dimensional DOA estimation for L-shaped array with nested subarrays without pair matching. IET Signal Process. 2016, 10, 1112–1117. [Google Scholar] [CrossRef]
  12. Niu, C.; Zhang, Y.; Guo, J. Interlaced Double-Precision 2-D Angle Estimation Algorithm Using L-Shaped Nested Arrays. IEEE Signal Process. Lett. 2016, 23, 522–526. [Google Scholar] [CrossRef]
  13. Pal, P.; Vaidyanathan, P.P. Nested Arrays in Two Dimensions, Part I: Geometrical Considerations. IEEE Trans. Signal Process. 2012, 60, 4694–4705. [Google Scholar] [CrossRef]
  14. Pal, P.; Vaidyanathan, P.P. Nested Arrays in Two Dimensions, Part II: Application in Two Dimensional Array Processing. IEEE Trans. Signal Process. 2012, 60, 4706–4718. [Google Scholar] [CrossRef]
  15. Liu, C.L.; Vaidyanathan, P.P. Two-dimensional sparse arrays with hole-free coarray and reduced mutual coupling. In Proceedings of the 2016 50th Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 6–9 Novemebr 2016; pp. 1508–1512. [Google Scholar]
  16. Liu, C.L.; Vaidyanathan, P.P. Hourglass Arrays and Other Novel 2-D Sparse Arrays With Reduced Mutual Coupling. IEEE Trans. Signal Process. 2017, 65, 3369–3383. [Google Scholar] [CrossRef]
  17. Chevalier, P.; Albera, L.; Ferreol, A.; Comon, P. On the virtual array concept for higher order array processing. IEEE Trans. Signal Process. 2005, 53, 1254–1271. [Google Scholar] [CrossRef] [Green Version]
  18. Pal, P.; Vaidyanathan, P.P. Multiple Level Nested Array: An Efficient Geometry for 2q-th Order Cumulant Based Array Processing. IEEE Trans. Signal Process. 2012, 60, 1253–1269. [Google Scholar] [CrossRef]
  19. Shen, Q.; Liu, W.; Cui, W.; Wu, S. Extension of nested arrays with the fourth-order difference co-array enhancement. In Proceedings of the 2016 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Shanghai, China, 20–25 March 2016; pp. 2991–2995. [Google Scholar]
  20. Shen, Q.; Liu, W.; Cui, W.; Wu, S. Extension of Co-Prime Arrays Based on the Fourth-Order Difference Co-Array Concept. IEEE Signal Process. Lett. 2016, 23, 615–619. [Google Scholar] [CrossRef] [Green Version]
  21. Shan, T.J.; Wax, M.; Kailath, T. On spatial smoothing for direction-of-arrival estimation of coherent signals. IEEE Trans. Acoust. Speech Signal Process. 1985, 33, 806–811. [Google Scholar] [CrossRef]
  22. Zoltowski, M.D.; Haardt, M.; Mathews, C.P. Closed-form 2-D angle estimation with rectangular arrays in element space or beamspace via unitary ESPRIT. IEEE Trans. Signal Process. 1996, 44, 316–328. [Google Scholar] [CrossRef]
  23. Schmidt, R. Multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar] [CrossRef]
Figure 1. A generalized L-shaped nested array (GLNA) with two classic nested arrays (CNAs) (GLNA-CNA) having the same parameters ( N d , N s ) = ( 4 , 5 ) distributed in the x, y directions with ( 0 , 0 ) and ( 6 , 6 ) as the reference. The circles represent physical sensors. The cross is a hole.
Figure 1. A generalized L-shaped nested array (GLNA) with two classic nested arrays (CNAs) (GLNA-CNA) having the same parameters ( N d , N s ) = ( 4 , 5 ) distributed in the x, y directions with ( 0 , 0 ) and ( 6 , 6 ) as the reference. The circles represent physical sensors. The cross is a hole.
Sensors 18 02482 g001
Figure 2. Difference co-arrays of the GLNA in Figure 1. (a) second-order difference co-array (SODC) (b) fourth-order difference co-array (FODC) (note that the unit of length in this figure is d). URA, uniform rectangular array.
Figure 2. Difference co-arrays of the GLNA in Figure 1. (a) second-order difference co-array (SODC) (b) fourth-order difference co-array (FODC) (note that the unit of length in this figure is d). URA, uniform rectangular array.
Sensors 18 02482 g002
Figure 3. An example of sGLNA with ( N x d , N x s , N y d , N y s ) = ( 4 , 5 , 4 , 5 ) and ( L x , L y ) = ( 15 , 24 ) . The circles represent physical sensors. The cross is a hole.
Figure 3. An example of sGLNA with ( N x d , N x s , N y d , N y s ) = ( 4 , 5 , 4 , 5 ) and ( L x , L y ) = ( 15 , 24 ) . The circles represent physical sensors. The cross is a hole.
Sensors 18 02482 g003
Figure 4. DOA estimations of different PSAs (a) 2D-NA, RMSE = 8.0128 ; (b) HA, RMSE = 2.5737 ; (c) GLNA-CNA, RMSE = 1.1827 ; (d) GLNA-ANA, RMSE = 0.5488 ; (e) sGLNA, RMSE = 0.5441 .
Figure 4. DOA estimations of different PSAs (a) 2D-NA, RMSE = 8.0128 ; (b) HA, RMSE = 2.5737 ; (c) GLNA-CNA, RMSE = 1.1827 ; (d) GLNA-ANA, RMSE = 0.5488 ; (e) sGLNA, RMSE = 0.5441 .
Sensors 18 02482 g004
Figure 5. DOA estimation comparisons. (a) RMSE as a function of SNR; (b) RMSE as a function of snapshot.
Figure 5. DOA estimation comparisons. (a) RMSE as a function of SNR; (b) RMSE as a function of snapshot.
Sensors 18 02482 g005
Table 1. Comparisons of DOF with the number of sensors as N = 18 . HA, hourglass array; ANA, augmented nested array; sGLNA, special GLNA.
Table 1. Comparisons of DOF with the number of sensors as N = 18 . HA, hourglass array; ANA, augmented nested array; sGLNA, special GLNA.
Array StructuresParametersDOF
HA ( N x , N y ) = ( 6 , 7 ) 525
2D-NA ( N x d , N x s , N y d , N y s ) = ( 3 , 3 , 3 , 3 ) 561
GLNA-CNA ( N x d , N x s , N y d , N y s ) = ( 4 , 5 , 4 , 5 ) 2401
GLNA-ANA ( N x d , N x s , N y d , N y s ) = ( 4 , 5 , 4 , 5 ) 2809
sGLNA ( N x d , N x s , N y d , N y s ) = ( 4 , 5 , 4 , 5 ) 3127

Share and Cite

MDPI and ACS Style

Zhang, L.; Ren, S.; Li, X.; Ren, G.; Wang, X. Generalized L-Shaped Nested Array Concept Based on the Fourth-Order Difference Co-Array. Sensors 2018, 18, 2482. https://doi.org/10.3390/s18082482

AMA Style

Zhang L, Ren S, Li X, Ren G, Wang X. Generalized L-Shaped Nested Array Concept Based on the Fourth-Order Difference Co-Array. Sensors. 2018; 18(8):2482. https://doi.org/10.3390/s18082482

Chicago/Turabian Style

Zhang, Lei, Shiwei Ren, Xiangnan Li, Guishan Ren, and Xiaohua Wang. 2018. "Generalized L-Shaped Nested Array Concept Based on the Fourth-Order Difference Co-Array" Sensors 18, no. 8: 2482. https://doi.org/10.3390/s18082482

APA Style

Zhang, L., Ren, S., Li, X., Ren, G., & Wang, X. (2018). Generalized L-Shaped Nested Array Concept Based on the Fourth-Order Difference Co-Array. Sensors, 18(8), 2482. https://doi.org/10.3390/s18082482

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop