Tri-polarized Sparse Array Design for Mutual Coupling Reduction in Direction Finding and Polarization Estimation

: Multi-polarized antenna arrays have the ability to provide both the direction and polarization information of the incident signals, which is important in radar, sonar, wireless communication, remote sensing, and so on. In this paper, a diversely polarized linear array of sparsely located but identically oriented tri-polarized vector antennas (VAs) is designed for estimating the direction-of-arrival (DOA) and polarization parameters of the incident signals in the presence of antenna mutual coupling (MC). In order to reduce the inter-VA MC, a new type of sparse array geometry is proposed, wherein the minimum inter-VA spacing is constrained to be no less than one signal wavelength. Considering the intra-VA MC e ﬀ ect, a full-wave electromagnetic simulation is introduced to ﬁt the manifold vector of an isolated VA. Based on the sparse VA array, a polarimetric subspace scheme is proposed for DOA and polarization estimation. When the knowledge about the intra-VA MC is a priori unavailable, an algebraic polarimetric blind scheme is also provided for DOA estimation. Computer simulations and real-world experiments (using an S-band 24-channel tri-polarized array system) validate the e ﬃ cacy of the designed array geometry along with the parameter estimation methods. type-2 polarimetric subspace scheme, direction-of-arrival).

A VA represents an antenna group which generally consists of several co-centered but diversely polarized dipoles/loops. An important problem in VA array signal parameter estimation is the alleviation of the intra-VA and inter-VA mutual coupling (MC) effects. For intra-VA MC reduction, one of the current measures is to develop a "distributed" ("stretched") VA whose constituent dipoles/loops are not spatially co-centered [10,[16][17][18]. Note that the approximate characteristics of an isolated VA (either co-located or distributed) involving MC can be also calculated via method of moments (MoM)-based full-wave electromagnetic simulations [21][22][23].
For inter-VA MC reduction, one possible strategy, in addition to array calibration, is to utilize sparsely located VAs. The array geometry may be minimum redundancy (MRA) [24], coprime [25], nested [26], super nested [27], or CADiS (coprime array with displaced subarrays) [28]. The minimum inter-VA spacing of the former four array geometries is half the signal wavelength. The MC between the corresponding VA pairs may not be neglected [21,29]. For CADiS, the inter-VA spacings are larger than half the signal wavelength. This is beneficial for alleviating inter-VA MC. However, the achievable consecutive coarray sensor number, and, thus, the resolvability of CADiS, is reduced compared with the other mentioned sparse arrays. There is still room to improve the design of sparse arrays.
In this paper, a new type of tri-polarized VA sparse array is designed, and two parameter estimation schemes are proposed. In the sparse array, the minimum inter-VA spacing is constrained to be no less than one signal wavelength so that the inter-VA MC can be neglected. Moreover, the aperture extension capability of the newly designed array is higher than that of CADiS. Considering the intra-VA MC, we utilize full-wave electromagnetic simulation to fit the manifold vectors of an isolated VA, based on which a polarimetric subspace scheme is proposed to estimate the direction and polarization parameters. For when the isolated VA manifold vectors are unavailable, a polarimetric blind scheme is presented for direction finding only, where the spatial diversity of the array is effectively used. The designed VA array geometry and the parameter estimation methods are validated by simulated data and real-world data obtained from a newly developed tri-polarized VA array system. Throughout the paper, superscripts "T", "*", and "H" represent the transpose, complex conjugate, and Hermitian transpose, respectively; j = √ −1 is the imaginary unit; "⊗" denotes the Kronecker product; "E{·}" is the expectation operator; "|S|" denotes the cardinality of set S; and "S 1 \S 2 ", termed as the set difference of S 1 and S 2 , represents the set of elements in S 1 but not in S 2 .

Tri-polarized Antenna Array Signal Model
Consider K uncorrelated far-field narrowband signals s k (t) K k=1 received by a linear array consisting of N tri-polarized VAs along the y-axis. The direction-of-arrival (DOA) of the kth signal is (θ k , ϕ k ) with θ k and ϕ k denoting the azimuth and elevation angles, respectively, as shown in Figure 1. The VA location set is S = {d 1 , d 2 , · · · d N } (in half signal wavelength). The output of the array can be modeled as a 3N × 1 vector: x(t) = K k=1 a(θ k , ϕ k , γ k , η k )s k (t) + n(t) (1) where a(θ k , ϕ k , γ k , η k ) = [a h (θ k , ϕ k ), a v (θ k , ϕ k )]h(γ k , η k ) is the steering vector of the kth signal, with a h (θ k , ϕ k ) and a v (θ k , ϕ k ) respectively characterizing the array responses to the horizontal and vertical polarimetric components of the kth signal, while h(γ k , η k ) = [cos γ k , sin γ k e jη k ] T represents the polarization vector of the kth signal. Here, γ k ∈ [0, π/2] and η k ∈ (−π, π] respectively denote the auxiliary polarization angle and the polarization phase difference of the kth signal. n(t) is the noise vector, which is assumed to be independent of the signals of interest. nested [26], super nested [27], or CADiS (coprime array with displaced subarrays) [28]. The minimum inter-VA spacing of the former four array geometries is half the signal wavelength. The MC between the corresponding VA pairs may not be neglected [21,29]. For CADiS, the inter-VA spacings are larger than half the signal wavelength. This is beneficial for alleviating inter-VA MC. However, the achievable consecutive coarray sensor number, and, thus, the resolvability of CADiS, is reduced compared with the other mentioned sparse arrays. There is still room to improve the design of sparse arrays.
In this paper, a new type of tri-polarized VA sparse array is designed, and two parameter estimation schemes are proposed. In the sparse array, the minimum inter-VA spacing is constrained to be no less than one signal wavelength so that the inter-VA MC can be neglected. Moreover, the aperture extension capability of the newly designed array is higher than that of CADiS. Considering the intra-VA MC, we utilize full-wave electromagnetic simulation to fit the manifold vectors of an isolated VA, based on which a polarimetric subspace scheme is proposed to estimate the direction and polarization parameters. For when the isolated VA manifold vectors are unavailable, a polarimetric blind scheme is presented for direction finding only, where the spatial diversity of the array is effectively used. The designed VA array geometry and the parameter estimation methods are validated by simulated data and real-world data obtained from a newly developed tri-polarized VA array system. Throughout the paper, superscripts "T", "*", and "H" represent the transpose, complex conjugate, and Hermitian transpose, respectively; j 1 = − is the imaginary unit; " ⊗ " denotes the Kronecker product; " {} E ⋅ " is the expectation operator; " | |  " denotes the cardinality of set  ; and " \ 1 2    ", termed as the set difference of 1  and 2  , represents the set of elements in 1  but not in 2  .

Tri-polarized Antenna Array Signal Model
Consider K uncorrelated far-field narrowband signals respectively denote the auxiliary polarization angle and the polarization phase difference of the kth signal. ( ) t n is the noise vector, which is assumed to be independent of the signals of interest.  If the inter-VA MC can be ignored (for example, if the minimum inter-VA spacing is no less than one signal wavelength), x(t) can be rewritten as where a s (θ k , ϕ k ) is the N × 1 spatial steering vector of the kth signal: a s (θ k , ϕ k ) = [e jπd 1 sin θ k sin ϕ k , e jπd 2 sin θ k sin ϕ k · · · , e jπd N sin θ k sin ϕ k ] T . (3) in which a iso,h (θ k , ϕ k ) and a iso,v (θ k , ϕ k ) respectively characterize the isolated VA responses to the horizontal and vertical polarimetric components of the kth signal. For a tripole VA, whose constituent antennas are three co-centered, orthogonally oriented short dipoles, a iso,h (θ k , ϕ k ) and a iso,v (θ k , ϕ k ) are given by a iso,h (θ k , ϕ k ) = [− sin θ k , cos θ k , 0] T and a iso,v (θ k , ϕ k ) = [cos ϕ k cos θ k , cos ϕ k sin θ k , − sin ϕ k ] T [1]. However, for the reasons of unknown anomalies in antenna design and inevitable intra-VA MC, closed-form expressions of a iso,h (θ, ϕ) and a iso,v (θ, ϕ) (also known as the horizontal and vertical isolated VA manifold vectors) are generally unavailable. In practice, full-wave electromagnetic simulations can be used to determine the parametric models of a iso,h (θ, ϕ) and a iso,v (θ, ϕ) from the gain and phase patterns of horizontal and vertical polarizations [21][22][23]. Figure 2 shows the topology and prototype of a tri-polarized antenna used in the S-band 24-channel array system. This antenna consists of two orthogonally polarized microstrips (y-and z-directions) and a collocated monopole (x-direction), with physical size larger than a half-wavelength. The gain patterns of horizontal and vertical polarizations for this tri-polarized antenna are shown in Figure 3. If the inter-VA MC can be ignored (for example, if the minimum inter-VA spacing is no less than one signal wavelength), ( ) t x can be rewritten as a a n (2) where ( , ) (3) is the 3 1 × isolated VA steering vector: , , respectively characterize the isolated VA responses to the horizontal and vertical polarimetric components of the kth signal. For a tripole VA, whose constituent antennas are three co-centered, orthogonally oriented short dipoles,  [1]. However, for the reasons of unknown anomalies in antenna design and inevitable intra-VA MC, closed-form expressions of ( , ) iso,h θ φ a and ( , ) iso,v θ φ a (also known as the horizontal and vertical isolated VA manifold vectors) are generally unavailable.
In practice, full-wave electromagnetic simulations can be used to determine the parametric models of ( , ) iso,h θ φ a and ( , ) iso,v θ φ a from the gain and phase patterns of horizontal and vertical polarizations [21][22][23]. Figure 2 shows the topology and prototype of a tri-polarized antenna used in the S-band 24-channel array system. This antenna consists of two orthogonally polarized microstrips (y-and z-directions) and a collocated monopole (x-direction), with physical size larger than a halfwavelength. The gain patterns of horizontal and vertical polarizations for this tri-polarized antenna are shown in Figure 3.

Inter-VA Spacing-Constrained Sparse Array
In this section, an inter-VA spacing-constrained (ISC) array geometry is designed to reduce the inter-VA MC.
Recall that the VA location set of the physical array is { , , , } Here, we In what follows, we show that a virtual uniform linear array (ULA) If the minimum inter-VA spacing of the physical array is constrained to be no less than one signal wavelength, then  The attempt here is to find a virtual ULA with as many as possible virtual VAs located at 2, 3, 4, 5, etc. (see also Figure 4 for the mentioned virtual ULA: the consecutive subarray of the virtual array). More precisely, given N and the array aperture L , the VA locations

Inter-VA Spacing-Constrained Sparse Array
In this section, an inter-VA spacing-constrained (ISC) array geometry is designed to reduce the inter-VA MC.
Recall that the VA location set of the physical array is S = {d 1 , d 2 , · · · d N }; the second-order difference coarray is Here, we set d 1 = 0. In what follows, we show that a virtual uniform linear array (ULA) V ⊂ D of no less than 2(N − 1) − 2 virtual VAs can be synthesized via a proper design of S. The inter-VA spacing of the virtual ULA is of half the signal wavelength. If N > 4, more degrees of freedom (DOFs) can be achieved for estimating the DOAs of signals. The processing aperture of the virtual ULA of VAs is also larger than that of a physical ULA of N VAs located at {0, 1, 2, · · · , N − 1} given that N > 4.
If the minimum inter-VA spacing of the physical array is constrained to be no less than one signal wavelength, then d m,n 1 and d m,n d N − 1; therefore, there will be at least two "gaps" at the positions 1 and d N − 1 in D (see, for example, the ISC arrays and their virtual counterparts as shown in Figure 4).

Inter-VA Spacing-Constrained Sparse Array
In this section, an inter-VA spacing-constrained (ISC) array geometry is designed to reduce the inter-VA MC.
Recall that the VA location set of the physical array is Here, we In what follows, we show that a virtual uniform linear array (ULA) VAs is also larger than that of a physical ULA of N VAs located at { , , , , If the minimum inter-VA spacing of the physical array is constrained to be no less than one signal wavelength, then  The attempt here is to find a virtual ULA with as many as possible virtual VAs located at 2, 3, 4, 5, etc. (see also Figure 4 for the mentioned virtual ULA: the consecutive subarray of the virtual array). More precisely, given N and the array aperture L , the VA locations  The attempt here is to find a virtual ULA with as many as possible virtual VAs located at 2, 3, 4, 5, etc. (see also Figure 4 for the mentioned virtual ULA: the consecutive subarray of the virtual array). More precisely, given N and the array aperture L, the VA locations {d n } N n=1 of the ISC array can be designed according to the following criterion: Additionally, in order to achieve a higher DOF than that of the physical array, it is desired that the aperture of the virtual ULA N V = N 0 − 1 be much larger than N. If N is small, the ISC array configurations satisfying criterion (6) can be found by exhaustive search through all the possible combinations of VA locations, or by using the computationally more efficient Scrunch recursive algorithm based on [30], which is summarized in Algorithm 1. If N is large, numerical algorithms such as simulated annealing [31] and ant colony optimization [32] can be used to find the suboptimal solution to (6).
Listed in Table 1 are some 8-VA ISC array geometries under different specifications of L and N 0 (for a certain value of L, the ISC array geometry may be nonunique; the mirror image is omitted for brevity). Analogous to [33], letting R and G be the numbers of redundant lags and gaps in D associated with the subaperture of the virtual array between 2 and d N − 2, respectively, we have the following remarks: 1.
Given N, the aperture parameter L is important for the design of the ISC array (L ≥ 2(N − 1) + 1 is required). If N ≥ 3, there exists a certain set of aperture parameters achieving N 0 = L − 2 (i.e., G = 0). The obtained "fully augmentable" ISC arrays involve different redundancies. Concerning aperture extension, the ISC array with "minimum redundancy" (where R is minimal or even zero) is of interest. For example, if N = 8, the use of any L belonging to the set {15, 16, · · · , 25} will lead to a fully augmentable ISC array. The ISC array with L = 25 is of minimal redundancy and the virtual aperture N V is larger than others. For different values of N, some minimum-redundancy ISC array configurations are listed in Table A1 of Appendix A.

2.
For the same setting of N, the virtual aperture of the ISC array can be further increased by using a larger L. The obtained ISC arrays now are "partially augmentable" (G 0). A partially augmentable ISC array is said to be "maximum contiguous" if the synthesized virtual ULA has the largest aperture. For example, if N = 8, N 0 = 24 (N V = 23) for the maximum-contiguous ISC array (L = 31), while N 0 = 23 (N V = 22) for the minimum-redundancy ISC array (L = 25).

Algorithm 1 Scrunch recursive algorithm for an ISC array
Input: the number of VAs N the array aperture L Output: the ISC array configuration S the position of the rightmost VA in the consecutive virtual array N 0 Procedure:

1.
Set two VAs at the endpoints of the ISC array, i.e., d 1 = 0, d N = L and S = {0, L}, with the number of remaining VAs N r = N − 2. Define two markers M L and M R (2 ≤ M L < M R ≤ L − 2) which indicate the positions the VAs attempted to set in the scrunching procedure from the ends of the array towards the center.

2.
Initialize N 0 with N 0 = L − 2 and proceed to steps 3-23. If no S is found, decrease N 0 and repeat the same steps until at least one S is found, then return S and N 0 . 3.
Denote the uncovered section of the consecutive virtual array {2, · · · , N 0 } in the scrunching procedure by V u = p L , p L + 1, · · · , p R , and initialize it with p L = 2 and p R = N 0 . The number of positions in the uncovered section N u is then initialized to N u = N 0 − 1.

4.
Repeat the following recursive function Scrunch (N r , M L , M R , V u , S) in steps 5-23.

5.
Calculate N u = |V u | and N m = N r (N r − 1)/2 + N r (N − N r ) which represents the maximum possible lags that can be generated by the N r VAs.
Determine the amount of scrunch ∆d: if M L ∈ S or M R ∈ S, ∆d = 2; else ∆d = 1.
Try to add a VA at position p t = M L + ∆d. Update M L by M L ← M L + ∆d . 14. else 15.
Try to add a VA at position Calculate the difference coarray D of the set S ∪ p t by using (5). 18. Call Return S, which satisfies the ISC criterion (6). 23. end if

ISC-Based DOA-Polarization Parameter Estimation
In this section, two subspace approaches to signal parameter estimation are presented by using the ISC array of tri-polarized VAs designed in Section 3.
Consider first the covariance matrix of the array output vector x(t), which is defined as R = E x(t)x H (t) . Under the assumptions made in the paper, for the ISC array, it can be verified that where a iso,k = a iso (θ k , ϕ k , γ k , η k ), a s,k = a s (θ k , ϕ k ), σ 2 k is the power of the kth signal, and σ 2 is the power of noise. In practice, R is replaced by its sample estimateR, while σ 2 is estimated by averaging the smaller eigenvalues ofR. By using the properties of the Kronecker product, We then partition the 3N × 3N matrix R as follows: Consider first the covariance matrix of the array output vector ( ) t x , which is defined as Under the assumptions made in the paper, for the ISC array, it can be verified that where , ( , , , ) k σ is the power of the kth signal, and 2 σ is the power of noise. In practice, R is replaced by its sample estimate R , while 2 σ is estimated by averaging the smaller eigenvalues of R . By using the properties of the Kronecker product, We then partition the 3 3 N N × matrix R as follows: in which uv R ( , , , sub-matrix given by (9) in which R uv (u, v = 1, 2, 3) is the N × N sub-matrix given by where r uv k = a iso,k (u)a * iso,k (v). Vectorizing R uv yields further the following coarray model: Following [24], an N V × N 2 selection matrix J (depending on S) can be constructed to synthesize the aforementioned virtual ULA of VAs associated with V, i.e., z uv V = Jr uv . We then divide the virtual ULA z uv V into N z = N V /2 overlapping subarrays of size N z , where " x " signifies the smallest integer not less than x. Letting the position of the N z th VA be the phase reference point, the nth subarray data can be written as where ρ n−1 k = exp −jπ(n − 1) sin θ k sin ϕ k and v(θ k , ϕ k ) = [1, e jπd sin θ k sin ϕ k , · · · , e jπ(N z −1) sin θ k sin ϕ k ] It follows that To fulfill DOA and polarization joint estimation, let us consider now the following polarimetric covariance matrix of the virtual ULA of VAs: where , , Vectorizing uv R yields further the following coarray model: , , Following [24], an It follows that [ , , , ] ( , ) ( , ) To fulfill DOA and polarization joint estimation, let us consider now the following polarimetric covariance matrix of the virtual ULA of VAs: where , , , , , is the virtual steering vector of the kth signal.
orthonormal matrix whose columns are the eigenvectors of R  associated with its 3 z N K − smallest (zero) eigenvalues. With U , DOA and polarization estimation then can be accomplished based on the rank deficient rationale, by applying, e.g., the two optional schemes presented below. (15) where where , , , is the virtual steering vector of the kth signal.
orthonormal matrix whose columns are the eigenvectors of R  associated with its 3 z N K − smallest (zero) eigenvalues. With U , DOA and polarization estimation then can be accomplished based on the rank deficient rationale, by applying, e.g., the two optional schemes presented below.

Using the Nominal Isolated VA Manifold Vectors
be the nominal horizontal and vertical isolated VA manifold vectors obtained via full-wave electromagnetic simulations; following [4], the DOA and polarization parameters can be obtained by solving the following optimization problem: where (16) is the virtual steering vector of the kth signal. Let U be the 3N z × (3N z − K) orthonormal matrix whose columns are the eigenvectors of R V associated with its 3N z − K smallest (zero) eigenvalues. With U, DOA and polarization estimation then can be accomplished based on the rank deficient rationale, by applying, e.g., the two optional schemes presented below.

Using the Nominal Isolated VA Manifold Vectors
Letâ iso,h (θ, ϕ) andâ iso,v (θ, ϕ) be the nominal horizontal and vertical isolated VA manifold vectors obtained via full-wave electromagnetic simulations; following [4], the DOA and polarization parameters can be obtained by solving the following optimization problem: The object function belongs to the standard generalized Rayleigh quotient, and the DOA parameters and polarization parameters can be estimated separately. To be specific, the DOA parameter can be found from the peaks of the following spatial spectrum: where "λ min {A, B}" denotes the minimal generalized eigenvalue of matrix pencil {A, B}. Letĥ k be the corresponding generalized eigenvector of D 1 (θ k ,φ k ), H H 1 (θ k ,φ k )H 1 (θ k ,φ k ) , withθ k andφ k denoting the estimates of θ k and ϕ k , respectively; the polarization parameters of the kth signal can then be estimated asγ k = arctan ĥ k (2)/ĥ k (1) andη k = ∠ ĥ k (2)/ĥ k (1) . We call this method the type-1 polarimetric subspace scheme (PS1).

Without Using the Nominal Isolated VA Manifold Vectors
If the polarization states of the signals are not of interest, or no knowledge about the isolated VA manifold vectors (or intra-VA MC) is available, an optional polarimetric blind scheme can be used to fulfill DOA estimation. Note that The object function belongs to the standard generalized Rayleigh quotient, and the DOA parameters and polarization parameters can be estimated separately. To be specific, the DOA parameter can be found from the peaks of the following spatial spectrum: Here, Since the VA manifold vectors are not available, we put the 3 × 1 vector a iso (θ, ϕ, γ, η) aside and only estimate the DOA parameters from the minimization of the generalized Rayleigh quotient. Similar to (21), we redefine the spatial spectrum as follows: The DOA parameters can be found from the peaks. We call this method the type-2 polarimetric subspace scheme (PS2).
A flow chart of the key steps of PS1 and PS2 is illustrated in Figure 5.

Analysis of Intra-VA MC Alleviation in PS1 and PS2
The vector antenna considered in this paper is composed of several co-centered but diversely polarized dipoles. The interaction between each close pair of dipoles is inevitable and difficult to model in closed-form expressions, especially when polarization information is also required in the subsequent algorithm. However, by virtue of full-wave electromagnetic simulations, the polarization responses of each dipole within the vector antenna can be calculated approximately, wherein the MC among the dipoles is also involved (see Figure 3 for the designed tri-polarized antenna). The manifold vectors of these dipoles can be obtained from the respective responses. In PS1, these manifold vectors, , are applied to estimate the DOA and polarization parameters. This method of alleviating intra-VA MC works better when the difference between the designed VA and practical VA is smaller.
On the other hand, if the manifolds for each VA are not available, we can just exploit the phase differences of VAs in the array to estimate the DOA. As shown in Equation (26) of PS2, the properties of vector antennas do not participate in the algorithm. So, the DOA estimation is not affected by the intra-VA MC. This method works better when the difference between the practical VAs is smaller.

Simulation and Experimental Results
In this Section, simulated and real-world data are used to validate the efficacy of the designed ISC array configuration and the performance of the parameter estimation schemes. In the following experiments, we consider only the one-dimensional azimuth angles of the incoming signals with identical elevation angle / 2 π for simplicity.

Capability of Intra-VA MC Alleviation
The performance of the proposed methods in reducing the intra-VA MC was studied first. The ISC array considered included 8 tri-polarized VAs with positions 0, 2, 4, 6, 11, 16, 21, and 24 (fully augmentable). Each VA was a triad of orthogonally collocated half-wavelength strip dipole antennas.

Analysis of Intra-VA MC Alleviation in PS1 and PS2
The vector antenna considered in this paper is composed of several co-centered but diversely polarized dipoles. The interaction between each close pair of dipoles is inevitable and difficult to model in closed-form expressions, especially when polarization information is also required in the subsequent algorithm. However, by virtue of full-wave electromagnetic simulations, the polarization responses of each dipole within the vector antenna can be calculated approximately, wherein the MC among the dipoles is also involved (see Figure 3 for the designed tri-polarized antenna). The manifold vectors of these dipoles can be obtained from the respective responses. In PS1, these manifold vectors, a iso,h (θ, ϕ) andâ iso,v (θ, ϕ), are applied to estimate the DOA and polarization parameters. This method of alleviating intra-VA MC works better when the difference between the designed VA and practical VA is smaller.
On the other hand, if the manifolds for each VA are not available, we can just exploit the phase differences of VAs in the array to estimate the DOA. As shown in Equation (26) of PS2, the properties of vector antennas do not participate in the algorithm. So, the DOA estimation is not affected by the intra-VA MC. This method works better when the difference between the practical VAs is smaller.

Simulation and Experimental Results
In this Section, simulated and real-world data are used to validate the efficacy of the designed ISC array configuration and the performance of the parameter estimation schemes. In the following experiments, we consider only the one-dimensional azimuth angles of the incoming signals with identical elevation angle π/2 for simplicity.

Capability of Intra-VA MC Alleviation
The performance of the proposed methods in reducing the intra-VA MC was studied first. The ISC array considered included 8 tri-polarized VAs with positions 0, 2, 4, 6, 11, 16, 21, and 24 (fully augmentable). Each VA was a triad of orthogonally collocated half-wavelength strip dipole antennas. Following [21,22], the manifold of a VA that includes intra-VA MC was generated by using the results of MoM-based full-wave electromagnetic simulation. As a comparison, three manifolds, each merely corresponding to a single type of strip dipole, were also considered, where intra-VA MC are not involved in these manifolds.
Consider four uncorrelated equal power signals from −45 • , −10 • , 10 • , and 30 • . All the signals have the same polarization state of (45 • , 90 • ). The snapshot number is 500. The signal-to-noise ratio (SNR) varies from −5 dB to 10 dB. The performance of PS1 and PS2 in the alleviation of intra-VA MC was compared with the estimation results without considering intra-VA MC. Shown in Figure 6 are the root-mean-square error (RMSE) curves of the DOA and polarization parameter estimates (obtained by averaging the results of 1000 independent trials). It can be observed that the three schemes provided almost the same DOA estimation, but PS1, which considered the intra-VA MC, had slightly better performance than the others. Moreover, the polarization parameter estimation by PS1 was clearly superior to that by the scheme in which intra-VA MC was not considered. These results demonstrate that the proposed schemes are effective in dealing with intra-VA MC, especially for polarization parameter estimation. MC, had slightly better performance than the others. Moreover, the polarization parameter estimation by PS1 was clearly superior to that by the scheme in which intra-VA MC was not considered. These results demonstrate that the proposed schemes are effective in dealing with intra-VA MC, especially for polarization parameter estimation.
Consider nine uncorrelated equal power signals from 55 −°, 41 −°, 28 −°, 14 −°, 0°, 14°, 28°, 41°, and 55°. The SNR is 10dB. The snapshot number is 500. All the signals have the same polarization state of (45 ,90 )°°. Shown in Figure 7 are the PS2 spatial spectra (100 independent runs) of the compared arrays in the presence of inter-VA MC. For comparison, the results assuming no inter-VA MC are also included in the figure. It is observed that, without inter-VA MC, all the six arrays can be used to obtain good DOA estimation, while CADiS appears to be the worst in terms of resolvability (its consecutive coarray aperture is shorter than that of the others). In the presence of inter-VA MC, the resolution capabilities and DOA estimation accuracies of the MRA, nested array, super nested array, and coprime array are observed to be degraded, possibly also exhibiting spurious peaks. The ISC and CADiS arrays are observed to be less deteriorated while the former is superior over the latter in terms of resolvability. Note that the compared arrays differ only in the VA spacing; the results demonstrate that, due to the VA spacing constraint of no less than one signal wavelength, RMSE (deg.)
Consider nine uncorrelated equal power signals from −55 • , −41 • , −28 • , −14 • , 0 • , 14 • , 28 • , 41 • , and 55 • . The SNR is 10dB. The snapshot number is 500. All the signals have the same polarization state of (45 • , 90 • ). Shown in Figure 7 are the PS2 spatial spectra (100 independent runs) of the compared arrays in the presence of inter-VA MC. For comparison, the results assuming no inter-VA MC are also included in the figure. It is observed that, without inter-VA MC, all the six arrays can be used to obtain good DOA estimation, while CADiS appears to be the worst in terms of resolvability (its consecutive coarray aperture is shorter than that of the others). In the presence of inter-VA MC, the resolution capabilities and DOA estimation accuracies of the MRA, nested array, super nested array, and coprime array are observed to be degraded, possibly also exhibiting spurious peaks. The ISC and CADiS arrays are observed to be less deteriorated while the former is superior over the latter in terms of resolvability. Note that the compared arrays differ only in the VA spacing; the results demonstrate that, due to the VA spacing constraint of no less than one signal wavelength, ISC and CADiS involve less serious inter-VA mutual coupling in comparison with the other arrays.

Robustness Against VA Misalignment
An S-band (2.4 GHz) 24-channel array system (including 8 tri-polarized VAs, see Figure 2) was developed, as shown in Figure 8. The VA positions of the system are adjustable. 0 90°° while the signal polarization phase difference η was 0° or 180°. The SNR was 0 dB. The snapshot number was 50. Shown in Figure 9 are the RMSE curves of the DOA and auxiliary polarization angle estimates versus polarization variation (obtained by averaging the results of 10,000 independent trials). It was observed that the array system exhibits a good robustness against VA misalignment for signal parameter estimation. The satisfactory performance of PS1 in DOA estimation is also observed, although the polarization estimation performance is not so good.

Robustness Against VA Misalignment
An S-band (2.4 GHz) 24-channel array system (including 8 tri-polarized VAs, see Figure 2) was developed, as shown in Figure 8. The VA positions of the system are adjustable.

Robustness Against VA Misalignment
An S-band (2.4 GHz) 24-channel array system (including 8 tri-polarized VAs, see Figure 2) was developed, as shown in Figure 8. The VA positions of the system are adjustable. 0 90°° while the signal polarization phase difference η was 0° or 180°. The SNR was 0 dB. The snapshot number was 50. Shown in Figure 9 are the RMSE curves of the DOA and auxiliary polarization angle estimates versus polarization variation (obtained by averaging the results of 10,000 independent trials). It was observed that the array system exhibits a good robustness against VA misalignment for signal parameter estimation. The satisfactory performance of PS1 in DOA estimation is also observed, although the polarization estimation performance is not so good. We investigated the robustness of the system against VA misalignment (causing different levels of signal-VA polarization mismatch) under a relatively simple condition: estimating the DOA and polarization parameter of a single linearly polarized signal with an 8-VA ISC array using the PS1 scheme. The VA positions were 0, 2, 4, 6, 11, 16, 21, and 24. The VA misalignment was characterized by the orientation angle of the monopole antenna measured from the x-o-y plane, denoted ψ. Here, 0 • ≤ ψ ≤ 45 • . The DOA of the signal was 45 • . The signal auxiliary polarization angle γ varied over [0 • , 90 • ] while the signal polarization phase difference η was 0 • or 180 • . The SNR was 0 dB. The snapshot number was 50. Shown in Figure 9 are the RMSE curves of the DOA and auxiliary polarization angle estimates versus polarization variation (obtained by averaging the results of 10,000 independent trials). It was observed that the array system exhibits a good robustness against VA misalignment for signal parameter estimation. The satisfactory performance of PS1 in DOA estimation is also observed, although the polarization estimation performance is not so good.

Parameter Estimation Using Real-World Data
Data were collected from the array system mentioned in Section 5.3. As shown in Figure 6, two signal sources were placed about 40 m away from the array. One was a linearly polarized horn antenna (polarization of about ( , ) 45 15°°); the other was an unmanned aerial vehicle (UAV) controller (polarization unknown). The DOAs of the radiated single tone signal (from the horn) and orthogonal quadrature phase shift keying modulated signal (from the UAV controller) were about 0° and 15°. The snapshot number was 1000. Displayed in Figure 10 are the spatial spectra of PS1 and PS2. For comparison, the results of the existing polarimetric rank reduction (PR) [4] and polarimetric average (PA) [9] methods are also included in the figure. The two signals can be observed to be well identified by all the methods, while PS1 and PS2 performed better than PR and PA. The polarization estimations obtained by PS1 were about ( . , . ) 44 1 21 8°° and ( . , . ) 68 3 110 0°°. Figure 10. Spatial spectra using real-world data (PR: polarimetric rank reduction method; PA: polarimetric average method).

Conclusions
A new sparse tri-polarized VA array called ISC was designed for the alleviation of inter-VA mutual coupling. Moreover, two unambiguous direction finding and polarization estimation methods, called PS1 and PS2, were presented for ISC. The behavior of ISC with parameter estimation using PS1 and PS2 was illustrated using simulated data and real-world data measured by a practical tri-polarized ISC array system. The results show that ISC is a good candidate for practical direction finding, regarding inter-VA MC reduction as well as aperture extension.   patial Spectra (dB) Figure 9. RMSE versus γ, η, and ψ: (a) DOA estimates; (b) auxiliary polarization angle estimates.

Parameter Estimation Using Real-World Data
Data were collected from the array system mentioned in Section 5.3. As shown in Figure 6, two signal sources were placed about 40 m away from the array. One was a linearly polarized horn antenna (polarization of about (45 • , 15 • )); the other was an unmanned aerial vehicle (UAV) controller (polarization unknown). The DOAs of the radiated single tone signal (from the horn) and orthogonal quadrature phase shift keying modulated signal (from the UAV controller) were about 0 • and 15 • . The snapshot number was 1000. Displayed in Figure 10 are the spatial spectra of PS1 and PS2. For comparison, the results of the existing polarimetric rank reduction (PR) [4] and polarimetric average (PA) [9] methods are also included in the figure. The two signals can be observed to be well identified by all the methods, while PS1 and PS2 performed better than PR and PA.

Parameter Estimation Using Real-World Data
Data were collected from the array system mentioned in Section 5.3. As shown in Figure 6, two signal sources were placed about 40 m away from the array. One was a linearly polarized horn antenna (polarization of about ( , ) 45 15°°); the other was an unmanned aerial vehicle (UAV) controller (polarization unknown). The DOAs of the radiated single tone signal (from the horn) and orthogonal quadrature phase shift keying modulated signal (from the UAV controller) were about 0° and 15°. The snapshot number was 1000. Displayed in Figure 10 are the spatial spectra of PS1 and PS2. For comparison, the results of the existing polarimetric rank reduction (PR) [4] and polarimetric average (PA) [9] methods are also included in the figure. The two signals can be observed to be well identified by all the methods, while PS1 and PS2 performed better than PR and PA. The polarization estimations obtained by PS1 were about ( . , . ) 44 1 21 8°° and ( . , . ) 68 3 110 0°°. Figure 10. Spatial spectra using real-world data (PR: polarimetric rank reduction method; PA: polarimetric average method).

Conclusions
A new sparse tri-polarized VA array called ISC was designed for the alleviation of inter-VA mutual coupling. Moreover, two unambiguous direction finding and polarization estimation methods, called PS1 and PS2, were presented for ISC. The behavior of ISC with parameter estimation using PS1 and PS2 was illustrated using simulated data and real-world data measured by a practical tri-polarized ISC array system. The results show that ISC is a good candidate for practical direction finding, regarding inter-VA MC reduction as well as aperture extension.   patial Spectra (dB) Figure 10. Spatial spectra using real-world data (PR: polarimetric rank reduction method; PA: polarimetric average method).

Conclusions
A new sparse tri-polarized VA array called ISC was designed for the alleviation of inter-VA mutual coupling. Moreover, two unambiguous direction finding and polarization estimation methods, called PS1 and PS2, were presented for ISC. The behavior of ISC with parameter estimation using PS1 and PS2 was illustrated using simulated data and real-world data measured by a practical tri-polarized ISC array system. The results show that ISC is a good candidate for practical direction finding, regarding inter-VA MC reduction as well as aperture extension.