Next Article in Journal
Authentication Technology in Internet of Things and Privacy Security Issues in Typical Application Scenarios
Next Article in Special Issue
Stochastic Maximum Likelihood Direction Finding in the Presence of Nonuniform Noise Fields
Previous Article in Journal
Anti-Jamming Low-Latency Channel Hopping Protocol for Cognitive Radio Networks
Previous Article in Special Issue
Wideband Direction-of-Arrival Estimation Based on Hierarchical Sparse Bayesian Learning for Signals with the Same or Different Frequency Bands
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hybrid T-Shaped Sensor Array Composed of Acoustic Vector Sensors and Scalar Sensors

1
School of Information Engineering, Nanchang Institute of Technology, Nanchang 330099, China
2
Key Laboratory for Information Science of Electromagnetic Waves and the Research Center of Smart Networks and Systems, School of Information Science and Technology, Fudan University, Shanghai 200433, China
*
Authors to whom correspondence should be addressed.
Electronics 2023, 12(8), 1813; https://doi.org/10.3390/electronics12081813
Submission received: 28 February 2023 / Revised: 4 April 2023 / Accepted: 7 April 2023 / Published: 11 April 2023
(This article belongs to the Special Issue Advances in Array Signal Processing)

Abstract

:
Through the more available acoustic information or the polarization information provided, vector sensor arrays outperform the scalar sensor arrays in accuracy of localization. However, the cost of a vector sensor array is higher than that of a scalar sensor array. To reduce the cost of a two-dimensional (2-D) vector sensor array, a hybrid T-shaped sensor array consisting of two orthogonal uniform linear arrays (ULAs) is proposed, where one ULA is composed of acoustic vector sensors and the other is composed of scalar sensors. By utilizing the cross-correlation tensor between the received signals from the two ULAs, two virtual uniform rectangular arrays (URAs) of acoustic vector sensors are obtained, and they can be combined into a larger URA. It is shown that a larger acoustic vector sensor URA with M2+1 degrees of freedom (DOFs) can be obtained from the specially designed T-shaped array with M acoustic vector sensors and 2M scalar sensors. Furthermore, by means of the proposed tensor model for the larger URA, the inter-sensor spacing can be allowed to exceed greatly a half-wavelength. Accordingly, the proposed method can achieve both a high DOF and a large array aperture. Simulation results show that the proposed method has a better performance in 2-D direction-of-arrival estimation than some existing methods under the same array cost.

1. Introduction

The two-dimensional (2-D) direction-of-arrival (DOA) estimation is known to be a fundamental problem in many fields such as wireless communications, radar, and sonar [1,2,3]. It is well known that the DOA estimation performance mainly depends on the degree of freedom (DOF) and array aperture [4,5]. Accordingly, achieving a higher DOF or/and a larger array aperture under the same array cost has received extensive attention.
To increase the achievable DOF with a limited number of sensors, the methods reported in the literature can be divided into two categories. One is to exploit new array configurations [6,7,8,9,10,11]. An effective method for doing that is to construct a virtual array, i.e., the difference coarray from the physical array covariance, with a higher degree of freedom (DOF) than that of the physical array. Two of the typical schemes reported in the literature are the nested array [6] and the coprime array [7]. In [8], a parallel coprime array structure and a novel algorithm for 2-D DOA estimation were proposed. By vectorizing the cross-covariance matrix of subarray data, the resulting virtual difference coarray enables re-solving more signals than the number of antennas. In [9], a generalized coprime planar array (GCPA) geometry for 2-D DOA estimation was proposed, where two rectangular uniform planar subarrays are used. GCPA geometry allows a more flexible array layout and extends the array aperture to achieve a great performance improvement. Zheng and Mu [10] proposed a method based on an augmented covariance matrix which is constructed using the output signals of two parallel difference coarrays. A method for two-dimensional (2-D) direction-of-arrival (DOA) estimation using two parallel nested arrays was proposed in [10], which utilized an increased array aperture and an enhanced DOF by forming an augmented covariance matrix based on the TPDC output signals. A novel sparse planar array consisting of multiple coprime and nested subarrays is designed and a corresponding 2-D DOA estimation method is developed in [11], and by vectorizing two covariance matrices, two virtual coprime planar subarrays made are available, which have many more virtual elements than the physical ones. Another approach is to develop new signal processing methods for conventional array configurations. The method in [12] utilizes the conjugate symmetry property of the ULA manifold matrix to increase the effective array aperture and virtual snapshots. To further increase the degree of freedom of the array, a tensor-based approach [13] divides each ULA in the L-shaped array into the optimum number of subarrays, and using the conjugate symmetry property of the ULA manifold matrix, a virtual URA with a much higher DOF is constructed. To make full use of the inherent spatial relevance among these coarray tensors, a coupled coarray tensor CPD approach is proposed to jointly decompose them for high-accuracy DOA estimation in a closed-form manner [14]. Although the maximum number of identifiable sources (i.e., DOF) for these methods can exceed the total actual number of the physical sensors, the arrays (e.g., coarray) used directly for DOA estimation in these methods require at least one array element spacing no larger than half a wavelength, which avoids the ambiguous angle estimation but limits the array aperture extension.
To extend the array aperture to be much higher than a half-wavelength with a limited number of sensors, vector sensors are widely used in sensing systems [15,16,17,18,19,20,21,22,23]. This is because a complete acoustic vector sensor (AVS) consists of four components, three orthogonal velocity sensors and another pressure sensor, co-located in space. An acoustic vector sensor can therefore measure all three particle-velocity field components plus the acoustic pressure induced by any acoustic incidence [24]. An AVS is capable of acquiring both acoustic pressure and three-dimensional particle velocities at any point in space using one omnidirectional sensor and three orthogonally co-located directional sensors, respectively. The 4D information obtained by a single AVS could help improve the signal processing performance [24]. Similarly, an electromagnetic vector sensor usually consists of three orthogonally oriented dipoles to measure the electric field, plus three orthogonally oriented loops to measure the magnetic field of the source. It can not only provide the DOA of the signal but also give polarization information [25]. Regarding the aforementioned vector arrays, [15,16] are two classic studies on the use of the ESPRIT algorithm for vector sensor URA, which also means that the two-dimensional aperture was extended. However, these algorithms require nontrivial pair-matching computations between two independent sets of direction estimates. To extend the array aperture, the sensor/element-spacing of methods in [17,18,19,20] can be intentionally extended to be much higher than a half-wavelength and hence would provide an enhanced spatial resolution. To increase the achievable DOF, the idea of associating electromagnetic vector sensors and acoustic vector sensors with nested arrays were proposed in [22,23]. However, they were confined to the usage of a linear array configuration. The array aperture in another direction was limited, and therefore, the performance of DOA estimation cannot be significantly improved. More importantly, the elements in these arrays are all individual complete vector sensors, i.e., four or six components. Under the condition of the same array aperture, the cost of this array is high. In order to reduce array costs while maintaining a large aperture in two directions, a coprime L-shaped array composed of a triangular SS-EMVS plus multiple components of SS-EMVS was proposed in [21]. As in [15,16], the method in [21] can extend the array aperture effectively in both directions, but its DOF is also less than the total number of physical sensors.
To simultaneously improve DOF and array aperture under the same array cost, a hybrid T-shaped sensor array for 2-D DOA estimation is proposed. This array configuration is composed of acoustic vector sensors and scalar sensors. From the second-order statistics of the received signal from the array, we can obtain two third-order data tensors. According to the structural characteristics of the T-shaped array, we reorder the slices of these tensors, and then we can obtain two new third-order tensors which correspond to two virtual uniform rectangular arrays (URAs) of acoustic vector sensors. Furthermore, we can directly combine them into another third-order tensor which corresponds to a larger URA. To maximize the DOF of the larger URA and simply calculate DOA, a fourth-order tensor model for the larger URA is derived. It is shown that a larger acoustic vector sensor URA with M2+1 DOFs can be obtained from the proposed array with M acoustic vector sensors and 2M scalar sensors. Furthermore, the virtual URA aperture can be further enlarged by extending the inter-sensor spacing to multiple half-wavelengths. Although cyclic ambiguity is introduced, it can be simply resolved by the proposed tensor model. Accordingly, the proposed method can provide not only a much higher DOF than the number of physical array elements, but also a large array aperture. Under the same array cost, the proposed method has a higher 2-D DOA estimation performance than the methods reported in the literature.
(·)*, (·)T, , denote conjugate, transpose, outer product, and Kronecker product, respectively.

2. Hybrid T-Shaped Sensor Array

As shown in Figure 1, in order to extend the same virtual array aperture in two directions, the proposed hybrid T-shaped sensor array consists of two orthogonal uniform linear arrays (ULAs) (denoted by X-ULA and Z-ULA, respectively) in the XZ plane, where the X-ULA has M acoustic vector sensors with spacing d, and the Z-ULA has 2M scalar sensors with spacing d. More precisely, the element positions in the X-ULA are given by {0.5d, 1.5d, …, 0.5d + (M − 1)d}, and the ones in the Z-ULA are given by {−0.5d − (M − 1)d, …, −1.5d, −0.5d, 0.5d, 1.5d, …, 0.5d + (M − 1)d}. The K narrow-band far-field uncorrelated signals with the power { σ k 2 } k = 1 K impinge on the array. The presence of kth source sk with elevation angle θ k [ 0 , π / 2 ) and azimuth angle φ k [ 0 , π ) is shown in Figure 1. Let Θ = e j 2 π cos ( θ k ) d / λ and Φ = e j 2 π cos ( φ k ) d / λ , where λ denotes the wavelength of the incident signal, and then the spatial steering vectors of X-ULA and Z-ULA can be represented as a x k = [ Φ 0.5 , , Φ M 0.5 ] T and a z k = [ Θ M + 0.5 , , Θ 0.5 , Θ 0.5 , , Θ M 0.5 ] T , respectively.
The received signal vector of Z-ULA and matrix of X-ULA from K sources at the nth snapshot can be represented as
z ( n ) = k = 1 K a z k s k ( n ) + w z ( n ) 2 M
X ( n ) = k = 1 K ( a x k p k ) s k ( n ) + W x ( n ) M × 4
where w z ( n ) and W x ( n ) are temporally and spatially Gaussian white with zero-mean additive noise vector and matrix corresponding to the Z-ULA and X-ULA, respectively, and p k is the spatial response vector of the acoustic vector sensor located at the origin. An acoustic vector sensor consists of an omnidirectional pressure sensor and up to three orthogonal particle velocity sensors, which can be expressed as [26]
p k = [ 1 , cos ( φ k ) , sin 2 ( θ k ) cos 2 ( φ k ) , cos ( θ k ) ] T

3. Tensor Model for Hybrid T-Shaped Sensor Array

3.1. Virtual URA with Acoustic Vector Sensors

We begin by constructing a virtual URA with acoustic vector sensors from the cross-correlation tensor [27] of z ( n ) and X ( n ) , i.e.,
Z X = E [ z ( n ) X * ( n ) ] = k = 1 K ( a z k a x k * p k * ) = k = 1 K ( a z k a x k * p k * ) σ k 2 2 M × M × 4
where a x k * = [ Φ 0.5 , , Φ M + 0.5 ] T . Since the noise is assumed to be spatially independent, and X-ULA and Z-ULA have no common array elements, the cross-correlation tensor Z X no longer contains the terms related to noise.
According to a x k = [ Φ 0.5 , , Φ M 0.5 ] T , we reverse the order of the lateral slices [27] of the tensor Z X , and we obtain a new tensor
1 = k = 1 K ( a z k a x k 1 p k ) σ k 2 2 M × M × 4
where a x k 1 = [ Φ M + 0.5 , , Φ 0.5 ] T . Moreover, since p k * and p k are the same in data Equation (3), the p k * can be directly written as p k . Obviously, 1 can be considered as the received data tensor of a virtual URA with acoustic vector sensors, where a x k 1 and a z k are the steering vectors for the kth source along the z and x axes, respectively, and p k is the spatial response vector of the acoustic vector sensor located at the origin of the kth source.
Therefore, tensor 1 can correspond to a (virtual) URA of 2M2 acoustic vector sensors which lies in the left half of the XZ plane, as depicted in Figure 2.
To obtain another virtual URA, we calculate the conjugate tensor of Z X firstly, i.e.,
Z X * = E [ z * ( n ) X ( n ) ] = k = 1 K ( a z k * a x k p k ) σ k 2 2 M × M × 4
According to a z k = [ Θ M + 0.5 , , Θ 0.5 , Θ 0.5 , , Θ M 0.5 ] T , we reverse the order of the horizontal slices [21] of Z X * , and then we can obtain a new tensor
2 = k = 1 K ( a z k a x k p k ) σ k 2 2 M × M × 4
Similarly, 2 can be considered as the received data tensor of another virtual URA of 2M2 acoustic vector sensors, which lies in the right half of the XZ plane, as shown in Figure 3.
Note that 1 and 2 have the same dimensions (i.e., a z k and p k ), so we can combine them into a larger one. Let V be the tensor of size 2 M × 2 M × 4 , which is given by
V ( : , 1 : M , : ) = 1
V ( : , M + 1 : 2 M , : ) = 2
and then the tensor V can be described as
V = k = 1 K ( a z k h x k p k ) σ k 2 2 M × 2 M × 4
where a z k = [ Θ k - M + 0.5 , , Θ k 0.5 , Θ k 0.5 , , Θ k M 0.5 ] T and h x k = [ a x k 1 , a x k ] = [ Φ k - M + 0.5 , , Φ k 0.5 , Φ k 0.5 , , Φ k M 0.5 ] T .
Similar to Equations (5) and (7), tensor V corresponds to a larger virtual URA, which contains all the acoustic vector sensors in the XZ plane as depicted in Figure 4. In addition, the element spacing of this larger virtual URA is d, which depends on the inter-sensor spacing of the physical T-shaped array.

3.2. Maximization of DOF

Consider V as a tensor of size 2 M × 2 M × 4 × 1 ; then, V behaves like the equivalent single-snapshot signals of the virtual URA. To increase the number of resolvable sources of the virtual URA, i.e., the DOF of the virtual URA, we increase the number of equivalent snapshots of V by applying the spatial smoothing method [6]. It is worth noting that each sensor in our virtual URA is an acoustic vector sensor, while the one in [6] is a scalar sensor. Therefore, we utilize tensor algebra to achieve spatial smoothing.
Let us divide the virtual URA into N Z N X overlapping identical subarrays of size L Z × L X . Note that L Z + N Z = 2 M + 1 and L X + N X = 2 M + 1 . Then, the received signal at the (nz, nx)th ( n z [ 1 , N Z ] , n x [ 1 , N X ] ) subarray in the virtual URA can be given by
V ( n z : n z + L Z 1 , n x : n x + L X 1 , : ) = k = 1 K a z k ( 1 ) Θ k n z 1 h x k ( 1 ) Φ k n x 1 p k σ k 2
where a z k ( 1 ) = [ Θ k M + 0.5 , Θ k M + 1.5 , , Θ k M + L Z 0.5 ] T and h x k ( 1 ) = [ Φ k M + 0.5 , Φ k M + 1.5 , , Φ k M + L X 0.5 ] T .
From Equation (11), we can see that V has a Vandermonde structure. Therefore, let Q be the tensor of size L Z × N Z × L X × N X × 4 , which is given by
Q ( : , n z , : , n x , : ) = V ( n z : n z + L Z 1 , n x : n x + L X 1 , : )
and then Q can be expressed as
Q = k = 1 K a z k ( 1 ) b z k h x k ( 1 ) b x k p k σ k 2
where b z k = [ 1 , Θ k , , Θ k N Z 1 ] T and b x k = [ 1 , Φ k , , Φ k N X 1 ] T .
Via Definition 1 in [13], one can combine the second and fourth dimensions of Q to construct the dimension corresponding to equivalent snapshots, i.e.,
  N = Q { 1 } { 3 } { 5 } { 2 , 4 } = k = 1 K a z k ( 1 ) h x k ( 1 ) p k n k σ k 2 L Z × L X × 4 × N Z N X
where n k = b z k b x k .
To further increase the DOF, we combine the first and second dimensions of N , i.e.,
  N = N { 1 , 2 } { 3 } { 4 } = k = 1 K f k p k n k σ k 2 L Z L X × 4 × N Z N X
where f k = a z k ( 1 ) h x k ( 1 ) is the steering vector of our virtual URA.
From Equation (15), we can see that N corresponds to a URA with L Z L X acoustic vector sensors and contains N Z N X equivalent snapshots.
From N , the 2-D DOA estimation of the URA can be accomplished by the Vandermonde recovery of the factor matrix. Even so, the 2-D DOA estimation via a simpler method is recommended. Therefore, once the 2-D DOA estimation is achieved by the methods mentioned in the following text, the complex calculations of Vandermonde recovery can be avoided. Similar to [13,28], in order to easily calculate the DOA of the incident signals, we use the Vandermonde structure of Q to rearrange Q as a seventh-order tensor:
N = k = 1 K a z k ( 1 ) b z k ( 1 ) c k h x k ( 1 ) b x k ( 1 ) e k p k σ k 2
where c k = [ 1 , Θ k ] T and e k = [ 1 , Φ k ] T ; b z k ( 1 ) = [ 1 , Θ k , , Θ k N Z 2 ] T and b x k ( 1 ) = [ 1 , Φ k , , Φ k N X 2 ] T .
We combine the dimensions of the seventh-order tensor N as follows:
T = N { 1 , 4 } { 2 , 5 } { 3 , 6 } { 7 } = k = 1 K f k p k r k n k ( 1 ) σ k 2
where r k = c k e k and n k ( 1 ) = b z k ( 1 ) b x k ( 1 ) .
Finally, T is the fourth-order canonical polyadic (CP) tensor model for our virtual URA.
Let k ( A ) represent Kruskal’s rank [27] of the matrix A , and then our tensor model T is sufficiently unique if [27]
k ( F ) + k ( P ) + k ( R ) + k ( N ( 1 ) ) 2 K + 3
where F = [ f 1 , , f K ] , P = [ p 1 , , p K ] , R = [ r 1 , , r K ] , and N ( 1 ) = [ n K ( 1 ) , , n K ( 1 ) ] .
We assume that the DOA pairs are chosen by the condition given by [6], which guarantees k ( F ) = min ( L Z L X , K ) , k ( P ) = min ( 2 , K ) , k ( R ) = min ( 4 , K ) , and k ( N ( 1 ) ) = min ( ( N Z 1 ) ( N X 1 ) , K ) . Then, applying the Lagrange multiplier method to Equation (18), we can find that when L Z = L X , the maximum number of identifiable sources for the proposed method is obtained, i.e.,
K M 2 + 1
Remark 1. 
For the proposed hybrid T-shaped sensor array composed of M acoustic vector sensors and 2M scalar sensors, we have established a fourth-order tensor T to process its received signals, which corresponds to a URA with M2 acoustic vector sensors and contains (M + 1)2 equivalent snapshots. Thus, this virtual URA has M2 + 1 DOFs.

3.3. Two-Dimensional Angle Estimation via an Extended-Aperture Hybrid T-Shaped Sensor Array

According to the proposed tensor model T , DOA information exists in different dimensions of the tensor. Utilizing the cp4_alsls MATLAB function [29] to carry out the CP tensor decomposition for T , one can obtain the estimations of R and P , i.e., R ^ and P ^ .
According to r k = c k e k , DOA estimates can be obtained from R ^ as follows:
u ¯ k = λ 4 π d { [ R ^ ( 2 , k ) R ^ ( 1 , k ) ] + [ R ^ ( 4 , k ) R ^ ( 3 , k ) ] } k = 1 , 2 , , K
v ¯ k = λ 4 π d { [ R ^ ( 3 , k ) R ^ ( 1 , k ) ] + [ R ^ ( 4 , k ) R ^ ( 2 , k ) ] } k = 1 , 2 , , K
where denotes returning the phase angle in the interval [−π,π] for each element of a complex array z. u ¯ k = cos ( φ ^ k ) and v ¯ k = cos ( θ ^ k ) are the cosine estimations of azimuth and elevation, respectively.
According to Equation (3), DOA estimates can also be obtained from P ^ as follows:
u ˜ k = P ^ ( 2 , k ) / P ^ ( 1 , k ) ,   k = 1 , 2 , , K
v ˜ k = P ^ ( 4 , k ) / P ^ ( 1 , k ) ,   k = 1 , 2 , , K
On the one hand, compared with u ˜ k and v ˜ k inherently extracted from information based on a single vector hydrophone which has no effective geometric aperture, u ˜ k and v ˜ k are extracted from information that encompasses the entire virtual array aperture and elements. Therefore, the estimates from R ^ have higher accuracy [15].
On the other hand, when the inter-sensor spacing of our array is set to be multiple half-wavelengths, there is a set of DOA estimates with cyclic ambiguity [30]. Assume a single source is impinging on the array from 2-D DOA ( φ , θ ) . The phase differences between the received signals at two adjacent sensors in the x-direction and z-direction are respectively denoted as
{ Δ φ x = mod ( 2 π λ d cos φ , 2 π ) Δ φ z = mod ( 2 π λ d cos θ , 2 π )
where the mod operation is based on the principle that the phase of a signal rotates 2 π for every λ distance the signal travels. Hence the relationship between the phase difference and inter-element spacing is that
{ Δ φ x + 2 k x π = 2 π λ d cos φ , 2 k x [ n , n ] Δ φ z + 2 k z π = 2 π λ d cos θ , 2 k z [ n , n ]
Since φ [ 0 , π ] and θ [ 0 , π / 2 ] , the above ranges for the two integers k x and k z are simplistically given according to their constraints 1 cos φ 1 and 0 cos θ 1 . For particular phase differences Δ φ x and Δ φ z , there exists one or a set of angles that satisfies Equation (25).
Specifically, when d 0.5 λ , both k x and k z can only take the value 0, which means that the azimuth angle and elevation angle have a one-to-one correspondence with cosine values within the range of 0 to 2 π . As d increases, the numbers of possible k x and k z values increase. Thus, all possible values of k x and k z exist in two sets u ¯ i , k and v ¯ j , k in the azimuth and elevation angles of the kth signal that satisfy Equation (25).
Compared with u ¯ k and v ¯ k , u ˜ k and v ˜ k do not suffer any similar extended-aperture ambiguity regardless of the inter-sensor spacing [15]. It should be emphasized that based on the fourth-order model T , we can simultaneously obtain these two sets of DOA estimates from different dimensions of the tensor. Thus, these two sets of estimates can be used for mutual disambiguation to yield a set of fine and unambiguous estimates without additional angle pairing, which is given by
I = min i , j ( | u ¯ i , k u ˜ k | + | v ¯ j , k v ˜ k | )
where u ^ i , k and v ^ j , k are taken as the true estimates if I obtains the minimum value. From the direction-cosine estimates derived above, the kth signal’s azimuth and elevation arrival angles may be estimated as
φ ^ = cos 1 ( u ^ i , k )
θ ^ k = cos 1 ( v ^ j , k )
The overall procedure of the proposed algorithm is summarized in Algorithm 1.
Remark 2. 
As in [15,16,21], the inter-sensor spacing in our array can also be allowed to exceed a half-wavelength greatly to extend the array aperture, while our method has a much higher DOF than others. As in [10,12,13], the proposed method can also provide a much higher DOF than the number of physical sensors, while only our array is able to enlarge the inter-sensor spacing to extend the array aperture. Accordingly, the proposed method can simultaneously offer high DOF and large array aperture.
Algorithm 1. Summary of the proposed algorithm.
Input: z and X of the Equations (1) and (2)
Calculate   Z X   and   Z X * in Equations (4) and (6)
Transform   Z X   and   Z X * ,   and   build   1   with   2
Combine   the   1   with   2 , and build V
Maximize DOF for V from Equations (11) to (17) and build T
Conduct CP decomposition on T and obtain two sets of DOA estimation from R ^ and P ^
Resolve cyclic ambiguity to obtain accurate estimates.
Output :   { θ ^ k }   and   { φ ^ k }

4. Numerical Simulation

To prove that our array can achieve better performance under the same array cost, we compare it with the vector sensor arrays of large array aperture [15,16,21]. In addition, since our array contains scalar sensors, some scalar arrays with high DOFs [10,12,13] are used as benchmarks. All the simulation results are obtained via 100 Monte Carlo trials. The root mean square error (RMSE) of parameter estimation is defined as
R M S E = 1 K J j = 1 J k = 1 K ( ( θ ^ k j θ k ) 2 + ( φ ^ k j φ k ) 2 )
where θ ^ k j and φ ^ k j are the estimations of elevation angle and azimuth angle in the jth experiment for the kth signal and J is the number of Monte Carlo trials.

4.1. Identifiability of the Proposed Method

We consider the hybrid T-shaped sensor array with 6 acoustic vector sensors and 12 scalar sensors. The snapshot number N, and signal-to-noise ratio (SNR) are set to be 500 and 10 dB, respectively. The sensor spacing of the proposed method d = 5 × ( λ / 2 ) . The 36 narrow-band waves impinge on the array, whose DOAs ( θ k , φ k ) are chosen by 36 pairs as {(52 + m1(74 − 52)/7, 52 + m2(74 − 52)/7), m1 = 1, …,6, m2 = 1, …,6} according to the condition given by [6]. It can be seen from the estimated results shown in Figure 5 that the proposed method can effectively handle the 36 sources which are more than 18 array elements.

4.2. RMSE vs. SNR

This example investigates the processing performance of the proposed method concerning SNR. In the proposed method, we utilize 6 acoustic vector sensors (containing 24 sensor components) and 12 scalar sensors. For the purpose of comparison, we consider scalar sensors as a component of vector sensors. Now, this proposed array contains 36 sensor components. Therefore, 36 vector sensor components are used in [15,16,21], under the same array cost. The SNR is varied from −5 dB to 20 dB. The azimuth and elevation angles of the two sources are set to (46°, 55°) and (48°, 62°), respectively. The snapshot number is set to 500. References [15,16,21] can allow 100× half-wavelength spacing, but the graph is constructed setting this to only 5×, i.e., d = 5 × (λ/2). The inter-sensor spacing of other methods is set to be (λ/2). From the results in Figure 6, we can see that [15,16,21] have better performance than [10,12,31] due to the array aperture advantage. Note that although the method in [13] utilizes the scalar sensor array needing at least one inter-sensor spacing not larger than half a wavelength, it can achieve a much higher DOF than [15,16,21]. Therefore [13] also shows better performance. The proposed method can simultaneously offer O ( M 2 ) DOF as in [13] and large array aperture as in [15,16,21]. Accordingly, the proposed method has the smallest RMSE among all considered methods.

4.3. RMSE vs. N

In this example, we examine the performance of the proposed method against the number of snapshots N. N is varied from 100 to 800, and SNR is fixed at 10 dB, while the other parameters are the same as in the second example. From the results shown in Figure 7, we can see that for all cases, the proposed method yields the best estimation results.

4.4. Runtime vs. the Number of Sensor Components

In this experiment, we compared the runtime required for each algorithm to run once in MATLAB 2022b on the same computer. The results are shown in Figure 8. For the purpose of comparison, we consider scalar sensors as a component of vector sensors. Other conditions are the same as Experiment 2, except that the total number of sensor components changes from 24 to 72. From Figure 6, it can be seen that the running time of the proposed method increases with an increase in the total number of sensor components. This is because the larger the number of sensor components is, the larger the size of T is, resulting in more time being needed to achieve the tensor decomposition. However, it is important to note that the larger size of T also means a higher DOF. Compared with other faster algorithms, although our algorithm may not be the most computationally efficient, it nonetheless belongs to the class of algorithms with relatively low computational complexity.

5. Conclusions

To reduce the cost of the 2-D vector sensor array, a new hybrid T-shaped sensor array composed of acoustic vector sensors and scalar sensors has been proposed. The tensor-based approach to this array processing with enhanced DOF and extended array aperture has been provided. Specifically, the proposed array contains M acoustic vector sensors and 2M scalar sensors, which are respectively placed along the x-axis and z-axis. Utilizing the cross-correlation tensor of the array received signals, a data tensor corresponding to a virtual URA with acoustic vector sensors is constructed. Applying the conjugate symmetry property of the ULA manifold matrix to the cross-correlation tensor, the data tensor of another virtual URA with acoustic vector sensors is also constructed. The analysis shows that these two virtual URAs can be combined into a larger virtual URA with 4M2 acoustic vector sensors. It is shown that a virtual URA with approximately M2 + 1 DOFs can be obtained. Since all the array elements in the virtual URA are acoustic vector sensors, the inter-sensor spacing can be extended with the help of the acoustic vector sensor characteristics. Thus, the proposed method has both a higher DOF and a larger array aperture. As demonstrated by simulation results, the proposed method can achieve superior 2-D DOA estimation performance to many existing methods under the same array cost. Additionally, 2-D DOA estimation is achieved by the CP decomposition, which means the multidimensional search can be avoided.
The hybrid concept proposed herein may be applied to other types of sensors, such as electromagnetic vector sensors, for highly accurate angle estimation with increased DOF and array aperture. Naturally, it can also be applied to other types of nonuniform arrays, such as nested arrays, coprime arrays, and minimum redundancy arrays.

Author Contributions

Conceptualization, W.R.; methodology, W.R.; software, Y.L.; validation, W.R.; formal analysis, W.R.; writing—original draft preparation, Y.L.; writing—review and editing, W.R.; supervision, W.R. and D.L.; All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China under Grant 61961025 and in part by the Jiangxi Provincial Natural Science Foundation under Grant 20202BABL202001.

Data Availability Statement

Data sharing is not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Krim, H.; Viberg, M. Two decades of array signal processing research: The parametric approach. IEEE Signal Process. Mag. 1996, 13, 67–94. [Google Scholar] [CrossRef]
  2. Zheng, Z.; Fu, M.; Wang, W.-Q.; Zhang, S.; Liao, Y. Localization of mixed near-field and far-field sources using symmetric double-nested arrays. IEEE Trans. Antennas Propag. 2019, 67, 7059–7070. [Google Scholar] [CrossRef]
  3. Wang, X.; Wan, L.; Huang, M.; Shen, C.; Han, Z.; Zhu, T. Low-complexity channel estimation for circular and noncircular signals in virtual MIMO vehicle communication systems. IEEE Trans. Veh. Technol. 2020, 69, 3916–3928. [Google Scholar] [CrossRef]
  4. He, W.; Yang, X.; Wang, Y. A high-resolution and low-complexity DOA estimation method with unfolded coprime linear arrays. Sensors 2019, 20, 218. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Aboumahmoud, I.; Muqaibel, A.H.; Alhassoun, M.; Alawsh, S.A. A Review of Sparse Sensor Arrays for Two-Dimensional Direction-of-Arrival Estimation. IEEE Access 2021, 9, 92999–93017. [Google Scholar] [CrossRef]
  6. Pal, P.; Vaidyanathan, 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]
  7. Vaidyanathan, P.P.; Pal, P. Sparse sensing with co-prime samplers and arrays. IEEE Trans. Signal Process. 2010, 59, 573–586. [Google Scholar] [CrossRef]
  8. Qin, S.; Zhang, Y.D.; Amin, M.G. Two-dimensional DOA estimation using parallel coprime subarrays. In Proceedings of the 2016 IEEE Sensor Array and Multichannel Signal Processing Workshop (SAM), Rio de Janeiro, Brazil, 10–13 July 2016; pp. 1–4. [Google Scholar]
  9. Zheng, W.; Zhang, X.; Zhai, H. Generalized coprime planar array geometry for 2-D DOA estimation. IEEE Commun. Lett. 2017, 21, 1075–1078. [Google Scholar] [CrossRef]
  10. Zheng, Z.; Mu, S. Two-dimensional DOA estimation using two parallel nested arrays. IEEE Commun. Lett. 2019, 24, 568–571. [Google Scholar] [CrossRef]
  11. Si, W.; Zeng, F.; Qu, Z.; Peng, Z. Two-dimensional DOA estimation via a novel sparse array consisting of coprime and nested subarrays. IEEE Commun. Lett. 2020, 24, 1266–1270. [Google Scholar] [CrossRef]
  12. Dong, Y.-Y.; Dong, C.-X.; Liu, W.; Chen, H.; Zhao, G.-q. 2-D DOA estimation for L-shaped array with array aperture and snapshots extension techniques. IEEE Signal Process. Lett. 2017, 24, 495–499. [Google Scholar] [CrossRef] [Green Version]
  13. Rao, W.; Li, D.; Zhang, J.Q. A tensor-based approach to L-shaped arrays processing with enhanced degrees of freedom. IEEE Signal Process. Lett. 2017, 25, 1–5. [Google Scholar] [CrossRef]
  14. Zheng, H.; Shi, Z.; Zhou, C.; Haardt, M.; Chen, J. Coupled Coarray Tensor CPD for DOA Estimation With Coprime L-Shaped Array. IEEE Signal Process. Lett. 2021, 28, 1545–1549. [Google Scholar] [CrossRef]
  15. Wong, K.T.; Zoltowski, M.D. Extended-aperture underwater acoustic multisource azimuth/elevation direction-finding using uniformly but sparsely spaced vector hydrophones. IEEE J. Ocean. Eng. 1997, 22, 659–672. [Google Scholar] [CrossRef] [Green Version]
  16. Zoltowski, M.D.; Wong, K.T. ESPRIT-based 2-D direction finding with a sparse uniform array of electromagnetic vector sensors. IEEE Trans. Signal Process. 2000, 48, 2195–2204. [Google Scholar] [CrossRef]
  17. Yang, M.; Ding, J.; Chen, B.; Yuan, X. A multiscale sparse array of spatially spread electromagnetic-vector-sensors for direction finding and polarization estimation. IEEE Access 2018, 6, 9807–9818. [Google Scholar] [CrossRef]
  18. He, J.; Zhang, Z.; Shu, T.; Yu, W. Sparse nested array with aperture extension for high accuracy angle estimation. Signal Process. 2020, 176, 107700. [Google Scholar] [CrossRef]
  19. Liu, S.; Zhao, J.; Zhang, Y.; Wu, D. 2D DOA Estimation Algorithm by Nested Acoustic Vector-Sensor Array. Circuits Syst. Signal Process. 2021, 41, 1115–1130. [Google Scholar] [CrossRef]
  20. He, J.; Li, L.; Shu, T. Sparse nested arrays with spatially spread square acoustic vector sensors for high-accuracy underdetermined direction finding. IEEE Trans. Aerosp. Electron. Syst. 2021, 57, 2324–2336. [Google Scholar] [CrossRef]
  21. Yang, M.; Ding, J.; Chen, B.; Yuan, X. Coprime L-shaped array connected by a triangular spatially-spread electromagnetic-vector-sensor for two-dimensional direction of arrival estimation. IET Radar Sonar Navig. 2019, 13, 1609–1615. [Google Scholar] [CrossRef]
  22. Fu, M.; Zheng, Z.; Wang, W.-Q.; So, H.C. Coarray interpolation for DOA estimation using coprime EMVS array. IEEE Signal Process. Lett. 2021, 28, 548–552. [Google Scholar] [CrossRef]
  23. Qu, X.; Lou, Y.; Zhao, Y.; Lu, Y.; Qiao, G. Augmented Tensor MUSIC for DOA Estimation Using Nested Acoustic Vector-Sensor Array. IEEE Signal Process. Lett. 2022, 29, 1624–1628. [Google Scholar] [CrossRef]
  24. Cao, J.; Liu, J.; Wang, J.; Lai, X. Acoustic vector sensor: Reviews and future perspectives. IET Signal Process. 2017, 11, 1–9. [Google Scholar] [CrossRef]
  25. Rahamim, D.; Tabrikian, J.; Shavit, R. Source localization using vector sensor array in a multipath environment. IEEE Trans. Signal Process. 2004, 52, 3096–3103. [Google Scholar] [CrossRef]
  26. Nehorai, A.; Paldi, E. Acoustic vector-sensor array processing. IEEE Trans. Signal Process. 1994, 42, 2481–2491. [Google Scholar] [CrossRef]
  27. Kolda, T.G.; Bader, B.W. Tensor decompositions and applications. SIAM Rev. 2009, 51, 455–500. [Google Scholar] [CrossRef]
  28. Sidiropoulos, N.D. Generalizing Caratheodory’s uniqueness of harmonic parameterization to N dimensions. IEEE Trans. Inf. Theory 2001, 47, 1687–1690. [Google Scholar] [CrossRef] [Green Version]
  29. Vervliet, N.; Debals, O.; De Lathauwer, L. Tensorlab 3.0—Numerical optimization strategies for large-scale constrained and coupled matrix/tensor factorization. In Proceedings of the 2016 50th Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 6–9 November 2016; pp. 1733–1738. [Google Scholar]
  30. Wu, Q.; Sun, F.; Lan, P.; Ding, G.; Zhang, X. Two-dimensional direction-of-arrival estimation for co-prime planar arrays: A partial spectral search approach. IEEE Sens. J. 2016, 16, 5660–5670. [Google Scholar] [CrossRef]
  31. Kawitkar, R. Performance of different types of array structures based on multiple signal classification (MUSIC) algorithm. In Proceedings of the 2009 5th International Conference on MEMS NANO, and Smart Systems, Dubai, United Arab Emirates, 28–30 December 2009; pp. 159–161. [Google Scholar]
Figure 1. Proposed hybrid T-shaped sensor array configuration.
Figure 1. Proposed hybrid T-shaped sensor array configuration.
Electronics 12 01813 g001
Figure 2. Virtual URA with Acoustic Vector Sensors in the left half of the XZ plane.
Figure 2. Virtual URA with Acoustic Vector Sensors in the left half of the XZ plane.
Electronics 12 01813 g002
Figure 3. Virtual URA with Acoustic Vector Sensors in the right half of the XZ plane.
Figure 3. Virtual URA with Acoustic Vector Sensors in the right half of the XZ plane.
Electronics 12 01813 g003
Figure 4. Virtual URA with Acoustic Vector Sensors in XZ plane.
Figure 4. Virtual URA with Acoustic Vector Sensors in XZ plane.
Electronics 12 01813 g004
Figure 5. Estimations of 36 sources.
Figure 5. Estimations of 36 sources.
Electronics 12 01813 g005
Figure 6. RMSE versus SNR.
Figure 6. RMSE versus SNR.
Electronics 12 01813 g006
Figure 7. RMSE versus N.
Figure 7. RMSE versus N.
Electronics 12 01813 g007
Figure 8. Runtime versus the number of sensor components.
Figure 8. Runtime versus the number of sensor components.
Electronics 12 01813 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rao, W.; Li, Y.; Li, D. Hybrid T-Shaped Sensor Array Composed of Acoustic Vector Sensors and Scalar Sensors. Electronics 2023, 12, 1813. https://doi.org/10.3390/electronics12081813

AMA Style

Rao W, Li Y, Li D. Hybrid T-Shaped Sensor Array Composed of Acoustic Vector Sensors and Scalar Sensors. Electronics. 2023; 12(8):1813. https://doi.org/10.3390/electronics12081813

Chicago/Turabian Style

Rao, Wei, Yuanqing Li, and Dan Li. 2023. "Hybrid T-Shaped Sensor Array Composed of Acoustic Vector Sensors and Scalar Sensors" Electronics 12, no. 8: 1813. https://doi.org/10.3390/electronics12081813

APA Style

Rao, W., Li, Y., & Li, D. (2023). Hybrid T-Shaped Sensor Array Composed of Acoustic Vector Sensors and Scalar Sensors. Electronics, 12(8), 1813. https://doi.org/10.3390/electronics12081813

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