# Generalized Sparse Polarization Array for DOA Estimation Using Compressive Measurements

Authors: Tao Chen; Jian Yang; Weitong Wang; Muran Guo
Journal: Wireless Communications and Mobile Computing
 (2021)

---

## Abstract

The compressive array method, where a compression matrix is designed to reduce the dimension of the received signal vector, is an effective solution to obtain high estimation performance with low system complexity. While sparse arrays are often used to obtain higher degrees of freedom (DOFs), in this paper, an orthogonal dipole sparse array structure exploiting compressive measurements is proposed to estimate the direction of arrival (DOA) and polarization signal parameters jointly. Based on the proposed structure, we also propose an estimation algorithm using the compressed sensing (CS) method, where the DOAs are accurately estimated by the CS algorithm and the polarization parameters are obtained via the least-square method exploiting the previously estimated DOAs. Furthermore, the performance of the estimation of DOA and polarization parameters is explicitly discussed through the Cramér-Rao bound (CRB). The CRB expression for elevation angle and auxiliary polarization angle is derived to reveal the limit of estimation performance mathematically. The difference between the results given in this paper and the CRB results of other polarized reception structures is mainly due to the use of the compression matrix. Simulation results verify that, compared with the uncompressed structure, the proposed structure can achieve higher estimated performance with a given number of channels.


---

## Body

## 1. Introduction

In a traditional scalar sensor array, the time delay of the phased array is used to estimate the direction of arrival (DOA). In the practical application environment, however, the signal to be detected usually has certain polarization characteristics. The polarization sensitive array [1–3] can be used to obtain and utilize the spatial and polarization domain information of the signal source comprehensively, which lays a physical foundation for improving the overall performance of the array signal processing. Therefore, the concept of polarization is extended to wireless communication [4], radar systems [5], and many other fields of space science [6]. At present, the research on polarization parameter estimation mainly focuses on how to improve the estimation performance, such as estimation accuracy and degree of freedom (DOF). As the polarized state of a signal varies with polarization diversity, only if the polarization direction of a single antenna matches with the incoming wave, all energy of the incoming wave can be received; otherwise, the loss of energy will occur [7]. Since a dual-polarized antenna [8, 9] can receive the signal energy along the horizontal and vertical branches, it can increase the receiving efficiency of signals.To obtain a higher number of DOFs under the premise of a given number of sensors, sparse array structures [10–12] have been proposed under the coarray framework, which generate equivalent virtual array elements via extracting the correlation information of received signals. On the one hand, the nested array [11] can generate a difference coarray with all continuous lags, which is very useful for spatial smoothing-based estimation algorithms. However, since one subarray has the sensors placed with a half wavelength, the mutual coupling effects will compromise the estimation performance. In [13], a sparse nested array has been proposed, where the interelement spacing of the dense subarray is extended to suppress mutual coupling effects. The diversely polarized dipoles are used to further improve estimation performance. Moreover, the spatially spread orthogonal dipoles are exploited in the sparse nested array in [14], where a passive direction finding structure with high accuracy has been proposed. Then, [15] further takes spatially spread square acoustic vector sensors into account, thus constructing a high-accuracy DOA estimation structure for an underdetermined case. On the other hand, the coprime array [12] can suppress the mutual coupling effect by using two sparse uniform linear arrays (ULA). In general, for a coprime array composed of an M-element subarray and an N-element subarray, up to 2M+N−2 uncorrelated sources can be distinguished. Therefore, the coprime array and its related improved structures [16] are also introduced into other radar structures [17]. In addition, [18] proposed an adaptive beamforming approach based on the coprime array, where the output performance is improved.The DOA parameters, including the azimuth and elevation angle, are of great importance in many applications, especially the unmanned driving technology [19] which is one of the current hot research issues. Thus, different catalogs of DOA estimation methods have been proposed, for instance, the subspace methods [7], deep learning methods [20], and sparse reconstruction methods [21]. The compressive sensing (CS) theory is the kernel of sparse reconstruction methods. However, CS is first used in time domain to break through the limitation of the Nyquist Sampling Theorem [22], since a high sampling rate is usually required for wideband signals, thus leading to a high cost of analog-to-digital converters (ADCs). To be specific, if the number of nonzero elements in a vector is much less than that of the zero elements, then this vector is regarded as a sparse vector and can be recovered from a small sample set. Thus, the sampling rate is dramatically reduced. Then, as the polarization sensitive array develops, several DOA and polarization joint estimation algorithms have been proposed [23–25]. Moreover, in addition to the estimation algorithms, the CS has also been used in system design. The DOA estimation system exploiting compressive measurements can dramatically reduce the complexity and the computational burden. For example, the compressive measurement method has found many applications in the coprime arrays [26], one-bit quantization [27], modulated wideband converter [28], and MIMO radar [29].Motivated by the above facts, in this paper, we mainly consider the orthogonal dipole antennas and propose a compressive measurement-based orthogonal dipole sparse array structure for the joint estimation of signal parameters. It is worth noting that there is no demand on the receive array structure, meaning that all kinds of dipole sparse array can be used. Based on the proposed structure, we first estimate the DOAs using a CS-based approach, and then the estimated DOAs are utilized to analyze the polarization parameters by a least-square estimation approach. In [24], the CS-based joint estimation of DOA and polarization parameters is proposed using a sparse array consisting of dual-polarized antenna elements. However, we would like to emphasize the contribution of this paper, as well as the difference with [24]:
(a)
In this paper, a two-dimensional polarization signal model is established, in which the azimuth angle and elevation angle are both taken into account, while [24] only considers the one-dimensional case, that is, only the azimuth angle is included in the signal model. In addition, it should be noted that a steering vector matrix ψθk,ϕk,γk,ηk is introduced into the signal model to describe the coherent structure in the polarizational and spatial domains(b)
The compression measurement method is applied to the proposed structure to compress the signal dimension by introducing a compression matrixΦ, therefore effectively reducing the number of channels required for subsequent digitization operations(c)
Considering that the Cramér-Rao bound (CRB) indicates the lower bound of the estimation error for an unbiased system, we first derive the CRB expression for the elevation angle and auxiliary polarization angle of the proposed structure. Then, theoretical performance analysis and simulation verification are made via the CRB expression in this paperNote that we use the CS theory twice in this paper. One is in the system design part, in order to reduce the dimension of received signal vector. The other is in the DOA estimation algorithm, where the group sparsity is used to obtain an improved number of DOF. Using the proposed structure, the number of channels is effectively controlled, thus reducing the hardware cost. In addition, although the compression leads to a degradation on the estimation performance, the proposed structure still outperforms the conventional dipole sparse array with the same number of channels, which can be clearly observed from both the theoretical derived CRB expression and the experimentally obtained root mean square error (RMSE). Therefore, the proposed structure also provides a flexible alternative option for low complexity polarization sensitive array with a relatively high estimation performance. Numerical simulations are conducted to examine the performance of the proposed structure.The following of this paper is organized as follows: In Section 2, we build the system model of the proposed structure. Then, a CS-based algorithm is proposed in Section 3 to jointly estimate the DOAs and polarization parameters. In Section 4, the CRB expression for the estimation of DOAs and polarization parameters of the proposed structure is derived. Numerical simulation results are shown in Section 5, where the corresponding analysis is given simultaneously. Finally, Section 6 concludes the whole paper.


## 2. System Model of the Proposed Structure

First, we would like to briefly review the receiving model of the polarization signals. Polarization sensitive receive array is composed ofL orthogonal dipoles, each of which is aligned with the y-axis in the Cartesian coordinate system. For the sake of convenience, the set S=dy1,dy2,⋯,dyL is used to represent the positions of the sensors arranged in ascending order, and the antenna at the origin is assumed to be the reference, i.e., dy1=0. It is noted that the sensors can be sparsely placed. For instance, the sensors can be arranged as an extended coprime array according to set S=Bad,0≤a≤2A−1∪Abd,0≤b≤B−1, where A and B are a pair of coprime integers. As shown in Figure 1, the signal from each dipole is processed separately; therefore, the array is divided into two subarrays according to the polarization receiving direction of the dipole:
(a)
All dipoles pointing in the direction of thex-axis constitute Subarray 1(b)
All dipoles pointing in the direction of they-axis constitute Subarray 2Figure 1
Schematic diagram of electromagnetic wave propagation.AssumeK narrow-band transverse electromagnetic (TEM) waves impinge upon the polarization sensitive array from the azimuth angle θk and elevation angle ϕk, where θk∈0,π and ϕk∈−2/π,2/π. It is assumed that each signal has an arbitrary elliptical electromagnetic polarization. The polarization of a TEM wave is often specified by two real parameters, namely, the auxiliary polarization angle γk (γk∈0,π/2) and the polarization phase difference ηk (ηk∈−π,π). The signal vector received by the polarization sensitive array, which has the dimension of Lo=2L, is expressed as
(1)xt=∑k=1Kψθk,ϕk,γk,ηk⊗uθk,ϕkskt+nt=Ast+nt, where we use ⊗ to denote the Kronecker product. Denote uθk,ϕk as the L-dimensional spatial steering vector of the kth signal, expressed as
(2)uθk,ϕk=ej2πdy1sinθksinϕkλ,⋯,ej2πdyLsinθksinϕkλT,in which λ is the signal wavelength. For simplicity of notation, we denote uθk,ϕk as uk. A vector matrix ψθk,ϕk,γk,ηk which describes the polarization information of incoming signals is defined as follows:
(3)ψθk,ϕk,γk,ηk=ψk=Ξθk,ϕkhγk,ηk,where Ξ and h are defined as
(4)Ξθk,ϕk=−sinθkcosϕkcosθkcosθkcosϕksinθk,hγk,ηk=cosγksinγkejηk,respectively. Thus, A=ψ1⊗u1,ψ2⊗u2,⋯,ψK⊗uK is the manifold matrix of polarization received signal, and st=s1t,s2t,⋯,sKtT is the complex envelope. Noise vector nt is assumed to be the zero mean complex Gaussian processes, where each of its entries is statistically independent.Then, letx1t be the signal received on Subarray 1 and x2t the signal received on Subarray 2. Then, replace ψk1 and ψk2 with ψxk and ψyk, respectively. The received signal vector of subarrays are expressed as
(5)x1t=∑k=1Kψxk⊗ukskt+n1t=A1st+n1t,x2t=∑k=1Kψyk⊗ukskt+n2t=A2st+n2t.The proposed orthogonal dipole sparse array structure exploiting compressive measurements method is shown in Figure2. The core principle of the compressive measurement method is to insert a combining network consisting of phase shifters and accumulators at the antenna outputs before subsequent digitization operations, which is equivalent to introducing a compression matrix Φ∈ℂM×L (M<L) for linear operations mathematically. In this way, the received signal vector in the L-dimension is compressed to M-dimension and then output for subsequent signal processing [30].Figure 2
System model of the proposed compressive measurement-based orthogonal dipole sparse array structure.It should be noted that the entries in the compression matrix are usually randomly selected from independent identically distributed parameters, and it is assumed that no additional noise is introduced during the compression process. In order to avoid information loss caused by data compression, the compression matrix can be optimized by various methods [27, 31, 32]. In this paper, Φ is selected to satisfy row-orthonormal, namely, ΦΦH=IM.Then, the received signal vector of each subarray after compression is expressed as(6)y1t=Φ∑k=1Kukψxkskt+n1t=A¯1st+n¯1t,y2t=Φ∑k=1Kukψykskt+n2t=A¯2st+n¯2t.Stacking the received vectors into a column vector yields(7)yt=A~st+n~t,where yt=y1Tt,y2TtT and A~=A¯1T,A¯2TT is the steering vector after compression with the kth term being
(8)a~k=ψk⊗u~k=ψxk,ψykT⊗Φuk.In addition,n~=n¯1Tt,n¯2TtT is the noise vector after compression.Therefore, for the proposed structure, the length of the compressive received signal vector is2M, while the length of conventional sparse dipole array with the same configuration is 2L. Assume that the number of snapshots is T. When the signal snapshots are used to compute the covariance matrix, the computational complexity of the proposed structure is O4TM2, while that of the conventional sparse dipole array is O4TL2.


## 3. Signal Parameter Estimation Approach

Since two-dimensional DOA estimation based on linear array cannot be realized, it is generally assumed that the signals and the linear array are in theyz plane, that is, the azimuth angle θk=π/2. To avoid three-dimensional parameters, we search for ϕk, γk, and ηk; [24] proposed a different reformulation, where a CS-based approach is first used to estimate the DOAs, and then the polarization parameters are estimated utilizing the estimated DOAs and a least-square estimation approach. In this paper, the above method is improved and applied to the proposed structure.We start with the array output covariance matrixRyy=EytyHt. The self-lag covariance matrix for the data vector yit and the cross-lag covariance matrix between yit and yjt can be obtained as
(9)Ryyi=ΦUψiRssψiHUH+pnILΦH,Ryyi,j=ΦUψiRssψjHUHΦH,respectively (1≤i≠j≤2), with Rss representing the source covariance matrix and pn indicating the noise power. In addition, U=u1,u2,⋯,uK is a spatial phase matrix with L×K dimension.We denote the vectorized form ofRyy as ryy, which can be regarded as a received data vector at a virtual array with an extended coarray aperture. On the basis of the matrix algorithm, the vectorization covariance matrix of different subarrays are calculated as
(10)ryyi=Φ0U∗⊙UvecψiRssψiH+pni~L2,ryyi,j=Φ0U∗⊙UvecψiRssψjH,in which we denote Φ∗⊗Φ∈ℂM2×L2 as Φ0 for notational simplicity. A matrix U∗⊙U=u1∗⊗u1,⋯,uK∗⊗uK that leads to a series of virtual array elements is defined, and the Khatri-Rao product is denoted by ⊙. Define a set of integers D=dyα−dyβdyα,dyβ∈S to represent the locations of virtual sensors and arrange them in ascending order. Then, the corresponding array manifold can be represented as UD=uD1,uD2,⋯,uDK.The vectorized covariance matrices are stacked and simplified as(11)ryy1ryy1,2ryy2,1ryy2=Φ0U∗⊙Up1Φ0U∗⊙Up1,2Φ0U∗⊙Up2,1Φ0U∗⊙Up2+pnΦ0i~L200pnΦ0i~L2,where i~L2=vecIL and the k-th item containing signal power is represented as
(12)pk1=pkcosγk2,(13)pk2=pkcosϕksinγk2,(14)pk1,2=−pkcosϕksinγkcosγke−jηk,(15)pk2,1=pkcosϕksinγkcosγkejηk.Discretizing the spatial domainΩP,Q (P,Q≫K) into a grid, let dgridθp,ϕq (1≤p≤P, 1≤q≤Q) represent the steering vector. Thus, the discretized array manifold corresponding to this grid can be obtained as
(16)Ggrid,i=Φ0dgrid∗θ1,ϕ1⊗dgridθ1,ϕ1,⋯,dgrid∗θP,ϕQ⊗dgridθP,ϕQ.It can be known that there is aγk (k∈1,K) such that sinγk⟶0 or cosγk⟶0. In this case, p1,2 and p2,1 approach zero simultaneously due to the item sinγkcosγk. Thus, using p1,2 and p2,1 has no improvement on the estimation performance. Meanwhile, the computational complexity is increased. However, it is impossible for cosγk2 and sinγk2 to be equal to zero at the same time. Thus, both p1 and p2 can be utilized for DOA estimation, and we have
(17)r=G~gridp+pnI~=G~p~,where r=ryy1T,ryy2TT, p=p1T,p2TT, I~=diagΦ0vecIL,Φ0vecIL, and G~grid=diagGgrid,1,Ggrid,2. The nonzero entries of p1 and p2 share the same support corresponding to the same grid. Thus, to utilize the group sparsity, hereby we define a sparse vector ξp~ as the ℓ2-norm of each row in the matrix p1,p2. The group LASSO algorithm [21] is utilized to solve the group sparsity problem, and the minimization problem is as follows:
(18)p~^=argminp~^12∥r−G~gridp~∥2+μ0∥ξp~∥1.The nonzero items inp~^ at its respective positions are the estimated DOAs.The estimated elevation angle is denoted asϕ^=ϕ^1,⋯,ϕ^KT, and the estimation of the array manifold is defined as G^=Ggrid∈ℂM2×K, which is used for the next polarization parameter estimation. Then, the vectors p1, p1,2, p2,1, and p2 which contain the polarization parameters can be obtained using the following least-square estimation:
(19)p^i,j=G^HG^−1G^Hri,j,i,j=1,2,where the kth item is expressed in (12)–(15). Simplifying the above equations, we have the following equations:
(20)cos2γk=p^k1−p^k2/cosϕ^k2pk,cosηksin2γk=p^k2,1−p^k1,2pkcosϕ^k.The estimated polarization phase differenceη^k can be obtained by expressing sin2γk in terms of cos2γk, whereas the estimated auxiliary polarization angle can be calculated by the following equation:
(21)γ^k=12atan2p^k2/p^k11/2,cosϕ^k.Thus, the joint estimation of DOA and polarization parameters is completed. The procedure of the proposed algorithm is summarized in Algorithm1.Algorithm 1: Procedure of the proposed algorithm.
procedure polarization signal parameter estimationInitializeryyi⟵vecRyyi,ryyi,j⟵vecRyyi,jforθp,ϕq, 1≤p≤P, 1≤q≤QdoG~grid⟵Φ0dgrid∗θ1,ϕ1⊗dgridθ1,ϕ1,⋯,dgrid∗θP,ϕQ⊗dgridθP,ϕQp~^⟵argminp~1/2∥r−G~gridp~∥2+μ0∥ξp~∥1endifp~^i≠0, i∈1,Qthenϕ^k⟵thepositionofp~^i, k=1,⋯,Kendforϕ^k∈ϕ^1,⋯,ϕ^K, pk∈p1,⋯,pKdocos2γk⟵p^k1−p^k2/cosϕ^k2/pkcosηksin2γk⟵p^k2,1−p^k1,2/pkcosϕ^kγ^k⟵1/2atan2p^k2/p^k11/2,cosϕ^kendreturnϕ^⟵ϕ^1,⋯,ϕ^KT, η^⟵η^1,⋯,η^KT, γ^⟵γ^1,⋯,γ^KTend procedure


## 4. The Cramér-Rao Bound Analysis

Since the linear array is exploited, without loss of generality, the azimuth angleθk and polarization phase difference ηk are set as π/2, thus limiting the DOAs to the yz plane and the polarization state to the same great circle orbit of Poincare sphere. The CRB [33] provides a lower bound on the covariance matrix of any unbiased estimator. In fact, under mild regularity conditions, the maximum likelihood estimator achieves the CRB asymptotically, as the number of snapshots tends to infinity. Given T independent samples of a zero mean Gaussian process yt whose statistics depend on a parameter vector α(22)α=ΩT,YT,pT,pn1,pn2T,with Ω=ϕ1,ϕ2,⋯,ϕKT, Y=γ1,γ2,⋯,γKT, and p=p1,p2,⋯,pKT, the CRB is obtained by the inverse of the Fisher informative matrix (FIM) [34].
(23)CRBα=1TGHΠΔ⊥G−1,in which
(24)G=M∂ryy∂ϕ1,⋯,∂ryy∂ϕK,∂ryy∂γ1,⋯,∂ryy∂γK,(25)Δ=M∂ryy∂p1,⋯,∂ryy∂pK,∂ryy∂pn1,∂ryy∂pn2,and M=RyyT⊗Ryy−1/2.Then, denotediagΦ,Φ as Ψ∈ℂ2M×2L, and vectorizing Ryy yields
(26)ryy=A~∗⊗A~vecRss+vecRnn.ExpandA~∗⊗A~ as ∑k=1Kψk⊗u~k∗⊗ψk⊗u~k, so that the ryy is finally simplified to the following form due to space limitation:
(27)ryy=∑k=1KV0HkΓkJuDkpk+V0pnvecILo.Several definitions are given to illustrate the vectorized results in (27):
(28)V0=Ψ∗⊗Ψ,Hk=hk∗1ILhk∗2IL⊗hk1ILhk2IL,Γk=−ILcosϕkIL⊗−ILcosϕkIL.Besides, the binary matrixJ has the definition as follows:
(29)<J>:,d=vecId,d∈D,where Id satisfies
(30)<Id>a,b=1,ifa−b=d,0,otherwise.Thus, by utilizing the relationshipuk∗⊗uk=JuDk, uk and uDk can be bridged.Taking the derivatives ofryy with respect to the DOA, polarization, signal power, and noise power, we have following results:
(31)∂ryy∂ϕk=2V0HkΓ′kJukpk+jπV0HkΓkJdiagDcosϕkukpk,(32)∂ryy∂γk=V0H′k+H′k∗ΓkJukpk,(33)∂ryy∂pk=V0HkΓkJuk,(34)∂ryy∂pn1=V0vecIL000,∂ryy∂pn2=V0vec000IL,where
(35)Γ′k=0−sinϕkIL⊗−ILcosϕkIL+−ILcosϕkIL⊗0−sinϕkIL,H′k=−sinγkILjcosγkIL⊗cosγkILjsinγkIL.Substituting (31)–(34) into (23), (24), and (25) leads to the CRB for the proposed structure.


## 5. Simulation Results

Throughout our simulations, the 10-element coprime array withA=3 and B=5 is considered. Without additional instructions, we assume the channel number after compression is M=8, and Φ∈ℂ8×10 is generated from the standard complex Gaussian distribution. The incident signal is generated uniformly in the range of ϕk∈−40∘,50∘ with θk=π/2, and the auxiliary polarization angle and polarization phase difference are evenly distributed in γk∈5∘,85∘ and ηk∈−120∘,150∘. Under the condition that SNR=5 dB and 1,000 snapshots, the elevation angle is estimated using the CS-based approach. The angle range from −60∘ to 60∘ is uniformly divided into grids with a 0.1∘ searching step.The signal parameter estimation results of 17 sources using the proposed structure are shown in Figure3, where the actual parameter values are marked in red. Compared with the results in [24], the proposed structure with channel number compression can still correctly analyze 17 signals with 10 physical elements. Then, the variation of CRB derived with the number of sources is simulated to verify the DOF that the proposed structure can obtain. To further analyze the estimation performance of the proposed structure, the following coprime array configurations are considered as the comparison structure: (a) the proposed structure with L=10,M=8; (b) the 10-element coprime array with the compression matrix Φ is equal to a unit matrix, that is, without compression, denoted as dipole coprime array 1; and (c) the 8-element coprime array with A=2 and B=5, namely, L=M=8, denoted as dipole coprime array 2, in which the idea of the CACIS configuration idea proposed in [35] is used for reference.Figure 3
Signal parameter estimation of the proposed structure (SNR=5dB and 1,000 snapshots).
(a)
Normalized spatial spectra(b)
Polarization parameter estimationFigure4 describes the CRB curves of ϕ and γ versus the number of sources when taking 1,000 snapshots at 10 dB SNR. When the compression ratio is L/M=1.25, the DOFs obtained by the proposed structure is not much less than that obtained by an uncompressed structure with the same number of physical elements, namely, dipole coprime array 1. However, compared with dipole coprime array 2 with the same number of channels, the proposed structure can resolve more uncorrelated signal sources, showing the superiority of the proposed structure for the number of DOFs.Figure 4
CRB versus the number of sources (SNR=10dB and 1,000 snapshots).
(a)
CRB (ϕ) versus no. of sources(b)
CRB (γ) versus no. of sourcesThe RMSE versus SNR of the comparison structure is shown in Figure5. The CRB (ϕ) and CRB (γ) of the three array configurations are presented by dash lines. The order of RMSE depicted by solid lines is consistent with that of CRB shown in Figure 4. It can be observed in Figure 4 that the downward trend of the RMSE curves can also fit the CRB curves well, while in Figure 4, the RMSE curves tend to be flat as the SNR increases since the search step is limited by the RIP criterion. Due to the largest number of channels, the estimation accuracy of dipole coprime array 1 is the highest among the three configurations, while the computational complexity is also the highest. In the case of the same number of channels, the proposed structure has a lower RMSE than dipole coprime array 2. In general, from the perspective of the number of inequalities, dimensionality reduction can inevitably lead to a decrease in estimation performance, such as DOF and estimation accuracy. In order to achieve the purpose of avoiding excessive system complexity, better estimation performance can be obtained by using the proposed structure for DOA estimation, which also verifies our previous analysis. On the other hand, in Figure 4, it can be observed that, in a low SNR region, the performance of the proposed structure is almost the same as that of the dipole coprime array 2, indicating that the compressive measurement-based structure has no significant improvement on the estimation of polarization parameters. However, in a large SNR region, the accuracy of the proposed structure approaches the dipole coprime array 1. We must note that this improvement is mainly caused by limitation of the step of the searching grid. To be specific, as the SNR increases, the estimation performance achieves the ceiling of current searching step, thus leading to the phenomenon that the proposed structure has almost the same performance as dipole coprime array 1. Theoretically speaking, by observing the CRB curves shown in Figures 3 and 4, the improvement on the estimation of polarization parameters is negligible.Figure 5
RMSE versus SNR (1,000 snapshots).
(a)
RMSE (ϕ) versus SNR(b)
RMSE (γ) versus SNR


## 6. Conclusion

In this paper, we proposed a compressive measurement-based orthogonal dipole sparse array structure, which can be used for high-performance signal parameter estimation with a small number of given elements. In the joint estimation algorithm of DOA and polarization parameters, the CS-based algorithm and the least-square estimation method were adopted. Then, based on compressive measurements, we derived the CRB expression for the elevation angle and auxiliary polarization angle. By comparing to the array configurations with the unit matrix as the compression matrix, we considered the CRB curves versus the number of independent signal sources and SNR. Thus, we come to the conclusion that, under the condition that we control the system complexity by reducing the number of channels, better parameter estimation performance can be obtained by the proposed orthogonal dipole array structure, especially for the estimation of DOAs. Simulation results also verified the theoretical analysis.



---

*Source: 5539709-2021-03-31.xml*