International Journal of Antennas and Propagation

Volume 2014 (2014), Article ID 394707, 14 pages

http://dx.doi.org/10.1155/2014/394707

Research Article

High Accuracy 2D-DOA Estimation for Conformal Array Using PARAFAC

Liangtian Wan,¹ Weijian Si,¹ Lutao Liu,¹ Zuoxi Tian,² and Naixing Feng³

¹Department of Information and Communication Engineering, Harbin Engineering University, Harbin 150001, China

²Science and Technology on Underwater Test and Control Laboratory, Dalian 116013, China

³Institute of Electromagnetics and Acoustics, Xiamen University, Xiamen 361005, China

Received 26 October 2013; Revised 14 December 2013; Accepted 16 December 2013; Published 16 January 2014

Academic Editor: Hon Tat Hui

Copyright © 2014 Liangtian Wan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Due to the polarization diversity (PD) of element patterns caused by the varying curvature of the conformal carrier, the conventional direction-of-arrival (DOA) estimation algorithms could not be applied to the conformal array. In order to describe the PD of conformal array, the polarization parameter is considered in the snapshot data model. The paramount difficulty for DOA estimation is the coupling between the angle information and polarization parameter. Based on the characteristic of the cylindrical conformal array, the decoupling between the polarization parameter and DOA can be realized with a specially designed array structure. 2D-DOA estimation of the cylindrical conformal array is accomplished via parallel factor analysis (PARAFAC) theory. To avoid parameter pairing problem, the algorithm forms a PARAFAC model of the covariance matrices in the covariance domain. The proposed algorithm can also be generalized to other conformal array structures and nonuniform noise scenario. Cramer-Rao bound (CRB) is derived and simulation results with the cylindrical conformal array are presented to verify the performance of the proposed algorithm.

1. Introduction

Conformal arrays mounted on curved surfaces are commonly applied in various areas such as radar, sonar, and wireless communication [1]. The conformal antennas could fulfill specific aerodynamics, space-saving, elimination of random-induced bore-sight error, potential increase in available aperture, and so on [2, 3]. Most researches focus on the design of antenna configuration [4–6], the transform between the local and global coordinate [7–10], and the pattern synthesis of conformal array [11–15].

The conventional direction-of-arrival (DOA) estimation algorithms, for example, the multiple signal classification (MUSIC) [16] and the estimation of signal parameters via rotational invariance techniques (ESPRIT) [17], are not applicable due to the varying curvature of the conformal array. The “shadow effect” of the conformal array is caused by the metallic shelter, leading to the condition that not all elements have the ability to receive the signal. As a result of the incomplete steering vector, common DOA estimation algorithms cannot be used for conformal array. Recently, DOA estimation algorithms for conformal array have been proposed for high resolution [18–22]. With the help of fourth-order cumulant and ESPRIT algorithm, a blind DOA estimation algorithm is proposed in [18]. Based on the mathematical technique of geometric algebra, 2D-DOA and polarization parameter estimations are completed by exploiting iterative ESPRIT algorithm [19]. The state space and propagator method are used for joint frequency and 2D-DOA estimations for the cylindrical conformal array [20]. However, the main drawback of the algorithms in [18–20] is the need of parameter pairing.

Parallel factor (PARAFAC) analysis attracted the attention of researchers when it was originally introduced in array signal processing in 2000 [21, 22]. It has been widely used for low-rank decomposition of three and higher way array. In this paper, a high accuracy 2D-DOA estimation algorithm for the conformal array is proposed using PARAFAC theory. The proposed algorithm does not need parameter pairing and performs well when the signal-to-noise ratio (SNR) is low or in situations where high estimation accuracy is needed. Simulation results demonstrate the efficiency and accuracy of the proposed algorithm.

The rest of this paper is structured as follows. Section 2 introduces the snapshot data model. Section 3 elaborates the design of the conformal array structure. Section 4 contains the core contributions of this paper, in which the PARAFAC model is used for high accuracy 2D-DOA estimation. Section 5 extends the proposed algorithm to more general cases. Section 6 analyses the computational complexity of the algorithm. Section 7 derives the Cramer-Rao bound (CRB). Section 8 discusses the array design. Section 9 presents the simulation result. Conclusions are drawn in Section 10.

2. The Snapshot Data Model

In order to achieve accurate parameter estimation, the snapshot data model of the conformal array must be precisely established. The pattern of the elements is different due to the varying curvature of carriers [8]. The mathematic model of the conformal array constructed in [8, 18] is adopted in this paper. The far-field narrowband signal impinging on the conformal array with elevationand azimuthis shown in Figure 1(a). The steering vector takes the following form:

where represents the position vector of the th element, is the propagation vector, and is wavelength of the incident signal. As shown in Figure 1(b), and are orthogonal unit vectors and and are the polarization components of the incident signal. denotes the th element’s response of the incident signal, is the pattern of the th element, and is the electric field direction of the incident signal. is the angle between vector and . denotes the transpose of matrix .

(a)

(b)

Figure 1: (a) The direction vector impinges on the array. (b) The th sensors response.

Assuming that incident signals impinge on the designed conformal array, the snapshot data model can be represented as

where denotes the pattern matrix and denotes the manifold matrix. and are the diagonal matrices whose diagonal entries are and , respectively.andare the components of th signal’s polarization parameters along the orthogonal unit vector and , respectively. denotes the signal vector. is additive Gaussian white noise with the covariance matrix

denotes conjugate transpose of matrix . Collecting snapshots,

where is the signal matrix and is the noise matrix.

The definition of the pattern is in the local coordinate, meaning that transformation must be performed from the global coordinate to the local coordinate [8]. In addition, the polarization parameters couple with the signal parameters; thus the decoupling is definitely necessary for the 2D-DOA estimation of the conformal array.

3. The Structure of the Array

Since the signal’s amplitude of the sensors received will decrease because of the “shadow effect,” the whole conformal array is divided into three subarrays, and each subarray covers . For the single curvature and symmetry of the cylinder, each subarray possesses and remains the same structure as well as the parameter estimation mechanism. Thus, this paper only considers one subarray.

As shown in Figure 2, the distance between two sensors which are in the same intersecting surface is and the distance between two adjacent intersecting surfaces is .

Figure 2: The structure of cylindrical conformal array.

As shown in Figure 2, constitute array 1, constitute array 2, constitute array 3, and constitute array 4 ( are sensors of the array). The distance vector between array 1 and array 2 is and . The distance vector between array 1 and array 3 is and as shown in Figure 3. The sensors’ patterns of the same generatrix are identical. The sensors possess the same pattern , and the sensors possess the same pattern . The decoupling between polarization parameter and angle information can be realized by taking full advantage of the characteristic of the array structure.

(a)

(b)

Figure 3: (a) The distance vector . (b) The distance vector .

Assume that , , , and represent the data that array 1, array 2, array 3, and array 4 receive, respectively. The origin of the coordinate is the reference point. So, the received data can be expressed as

The covariance matrices among the received data are

where represents the signal covariance matrix and the noise covariance matrices are , respectively.

Consider

where and represent the elevation and azimuth of the distance vector in the global coordinate. The array structure is shown in Figures 2 and 3 in which the distance vector is parallel to the -axis. The elevation and azimuth of in the global coordinate are

is parallel to the -axis. The elevation and azimuth of in the global coordinate are

In real applications, the ideal covariance matrix has to be replaced by the sample covariance matrix which is obtained by a finite number of snapshots   ; that is,

4. The DOA Estimation Based on PARAFAC

Parallel factor (PARAFAC) analysis is a method of multiple data analysis, which has been first introduced in psychometrics. It has been used in many fields such as statistics, arithmetic complexity, and chemometrics. In this section, first, the PARAFAC theory is introduced briefly; second, based on the PARAFAC model, the multiple way array is constructed and the identifiability of inherently unresolvable source permutation and scaling ambiguities are analyzed; finally, the trilinear alternating least squares (TALS) regression algorithm is used to fit the PARAFAC model; thus, the unknown parameters in the PARAFAC model can be estimated.

4.1. Parallel Factor Analysis

In this subsection, the basic idea of PARAFAC theory is introduced. First, the trilinear decomposition of the element of the three-way array is elaborated. Second, based on the symmetrical characteristic of the trilinear decomposition, the three-way array is decomposed into three dimensions.

Considering a three-way array (the entry of the array denotes ), the trilinear decomposition of can be expressed as [21]

where , , and . The three-way array is represented as the sum of rank-1 factor. Here, the rank of is defined as the minimum number of rank-1 three-way components which are needed to decompose . The vectors , , and are called load vector, score vector, and factor profiles, respectively.

Assume ; an example is given as follows and thus we can understand the idea of trilinear decomposition intuitively. The decomposition of three-way array is shown in Figure 4; it can be seen that three dimensions can be decomposed.

Figure 4: The decomposition of three-way array.

The typical element of the matrix is defined as ; the typical element of the matrix is defined as ; the typical element of the matrix is defined as . In addition, a matrix , a matrix , and a matrix are defined. The typical element of each matrix is . The three-way matrix can be “sliced” along three different dimensions as shown in Figure 4.

Consider

where is the diagonal matrix constructed by the th row of the matrix . According to the definition of Khatri-Rao product, can be written in another form as

where “” is used to denote the Khatri-Rao (KR) product. stands for the three-way array which is constructed by stacking the matrices into a column. The definition of Khatri-Rao (KR) product can be found in the appendix. The three-way array can also be expressed in two other forms:

The reason of decomposition along three different dimensions is that the PARAFAC model can be solved by using TALS algorithm. The matrices , , and can be updated alternately.

4.2. The Identifiability of the Model

In this subsection, the PARAFAC model is constructed by the received data of the conformal array. The decomposition of the PARAFAC model must be unique; if not so, the PARAFAC model cannot be solved. Thus, some requirements of the matrices which construct the model have to be satisfied. Theorem 2 guarantees the uniqueness decomposition of the model. If the PARAFAC model is constructed by the received data of the conformal array, then Theorem 2 can be used to judge the uniqueness decomposition of the model.

On the basis of the PARAFAC theory [23], the three-way array of the cylindrical conformal array is constructed by (7)

where    are the sample covariance matrices, respectively, and represents the noise in real observation. Let ; based on the definition of Khatri-Rao product, (20) can be transformed as

where is the row vector constructed by the diagonal entries of the diagonal matrix . stands for the matrix.

Definition 1 (see [21]). Considering a given matrix , if and only if the matrix contains at least but not linearly independent columns, then the rank of matrix is . If any columns of matrix are linearly independent, then the Kruskal rank (-rank) is [24]. Generally speaking, .

In order to ensure the uniqueness of the PARAFAC model in (21)–(23), the sufficient condition must be met:

For any two different incident signals, and ,

ensures the uniqueness of the PARAFAC model, where

and are the elevation and azimuth of the sensor’s position in the global coordinate. and are the elevation and azimuth of the th incident signal. Thus, (26) is equivalent to

Theorem 2. The sufficient conditions of the PARAFAC model in (21)–(23) require that at least one of and and is not zero.

Proof. The -rank of matrices , , and is , , and , respectively. The -rank decomposition of the model is unique with for (see Theorem  1 in [21]).

4.3. The Trilinear Alternating Least Squares Regression Algorithm

The principle of the trilinear alternating least squares (TALS) regression algorithm is to fit the PARAFAC model in the noisy observation. The idea of TALS is very simple. In each step, only one matrix is updated. The updating algorithm is based on the estimated results of the last updating, and the remaining matrices are updated by least squares method. Repeat the step as mentioned above until the algorithm converges. The advantage of TALS algorithm is that the parameter adjustment is not necessary. A standard least square problem is solved in each step, and the performance of TALS is good [25]. Empirically, if the -rank condition (Theorem 2) is satisfied, the TALS algorithm can converge to the global minimum [26]. The accelerating convergence technique of TALS algorithm can be found in [27, 28]. A brief introduction about TALS algorithm can be found in [21]. The application details of the TALS algorithm which is used to fit the PARAFAC model are elaborated as follows.

On the basis of noisy observation, the problem (21) can be transformed into solving a least square problem

The principle of alternating least squares (ALS) can be used to fit the problem in (30). In the noiseless condition, ALS can be used to solve the matrices , , and which constructed the three-way array . The least square estimation of matrixcan be expressed as

Similarly, the matrices and can be expressed as

In the iterative procedure, given matrices and , the matrix can be represented as

The expression of matrices and is

where denotes the pseudoinverse of matrix .

Now, the ALS algorithm steps can be summarized as follows:

(1)initialize , ;(2)initialize , ;(3)if , calculate matrices , , and by (35)–(37). Update just one matrix at each time; then ;(4)else ; the iteration is terminated.

4.4. The 2D-DOA Estimation Algorithm

The estimators of matrices , , and are obtained by the TALS algorithm which is introduced in the last subsection. In this section, the 2D-DOA estimation can be obtained by the matrix , which is shown as follows.

The matrix can be acquired by the TALS algorithm. and can be calculated by matrix

where represents the th row of . stands for the absolute value operation. Because and are real numbers, and are squared to solve the ambiguity caused by the positive and negative values of and .

Consider

Then,

Take (12), (13) into (10), (11), and the elevation and azimuth of th incident signal can be expressed as

Based on Theorem 2, the estimators of , , and have the same column permutation matrix; that is, the th column of the steering matrix corresponds to the th of the matrix . Thus, the elevation and azimuth pair with each other automatically.

Combining PARAFAC theory with the TALS algorithm, the 2D-DOA estimation for the cylindrical conformal array can be summarized as follows:

(1)calculate the signals’ covariance matrices received by each subarray using (7);(2)construct the PARAFAC model by (20)–(23);(3)estimate the matrix using the TALS algorithm;(4)Calculate (36) and (37) using the estimator of matrix ; then, obtain the estimators of and ;(5)acquire the elevation and azimuth estimation from (40) and (41).

5. The Extendable Array Structure

The proposed algorithm can be extended to other array structures with little modification. First, some elements are added in the cylindrical conformal array, and then the elements can be arranged more flexibly. Second, the proposed algorithm is extended to conical conformal array.

5.1. The Extendable Cylindrical Conformal Array

First, the design of cylindrical conformal array is introduced. Second, the PARAFAC model is constructed by the received data. Finally, the 2D-DOA estimation is obtained.

As shown in Figures 2 and 3, the proposed algorithm is feasible when the distance vector between array 1 and array 3 is parallel to -axis. However, the proposed algorithm must be modified to be used in general cases. Arrays 5 and array 6 are added so that the arrays can be designed more flexibly. The extendable cylindrical conformal array is shown in Figure 5. is array 5, and is array 6. The sensors possess the same pattern . The distance vector between array 1 and array 5 is , and . Under the design in Figure 5, the only restriction is that arrays are parallel to -axis.

Figure 5: The extendable cylindrical conformal array structure.

The received data of array 5 and array 6 are

where

The cross-covariance matrices among array 1, array 5, and array 6 are

whereandare noise covariance matrices. Then, (20) is modified into anthree-way array PARAFAC model

where is the observed noise. The form of Khatri-Rao product is applied in (44)

where

As long as Theorem 2 holds, (45) is unique. can be obtained similarly with (37)

after calculating matrix by ALS algorithm.

Assume that , , and ; similar to and   , , , and ; solving (10), (11), (42) (36), (37), and (24), we can obtain

The solution of (48) is

, , and can be acquired by solving (49). Two methods can be used to obtain, . The first method can be described as

The second one is

5.2. The Extendable Conical Conformal Array

In this section, the proposed algorithm is extended to conical conformal array. First, the design of conical conformal array is introduced. Second, the PARAFAC model is constructed by the received data. Finally, the 2D-DOA estimation is obtained. The analytic solution of the conical conformal array does not exist. Thus, the iterative numerical approximation method is used to obtain the optimal solution.

Figure 6 applies the proposed algorithm to the conical conformal array. is array 1, is array 2, is array 3, and is array 4. The distance vector between array 1 and array 2 is , and . The distance vector between array 3 and array 4 is , and . possess the same pattern , and possess the same pattern . The decoupling between polarization parameter and DOA information can be completed by making full use of the geometry characteristic of the array structure.

Figure 6: The extendable conical conformal array structure.

The received data from array 1 to array 4 are

where are the noise matrices received by each array. The cross-covariance matrices among the arrays are

where is the manifold matrix of array 1 and array 2 and is the manifold matrix of array 3 and array 4. are the noise covariance matrices received by the arrays. and possess the same forms as (8) and (9). According to PARAFAC theory, the three-way array for conical conformal array can be constructed by (56)

and is the observed noise. Let , and (57) can be written in the following form of Khatri-Rao product:

The matrix can be obtained by the ALS algorithm if the uniqueness of (57) is guaranteed. By using the similar method, and can be represented as

The 2D-DOA estimation can be acquired by solving (61) and (62). The equations are nonlinear and the analytical solution does not exist. Thus, the iterative numerical approximation method is used to obtain the optimal solution. The nonlinear equations can be written as

where and . The equation with the same solution of (63) is expressed as

The iterative form can be represented as

The initial vector is selected as . Calculating (65) until the sequence converges, that is, , and regarding as the iterative solution of (63), we can achieve the 2D-DOA estimation of the incident signal.

Assuming that the noise is nonuniform, in theory through calculating the cross-covariance matrices among different arrays, nonuniform noise could be suppressed.

6. The Computational Complexity Analysis

We only focus on the major part, which is the number of multiplications involved in calculating covariance matrices and the ALS algorithm. As mentioned above, , , and represent the snapshot, source number, and sensor number, respectively. The ESPRIT algorithm which is similar to the algorithm proposed in [18] needs to calculate the eigendecomposition of covariance matrices and parameter pairing (the covariance matrix is used instead of the four-order cumulant in [18]). The eigendecomposition of three matrices requires about . The algorithm in this paper uses COMFAC algorithm to fit an three-way array. The computational complexity per iteration is . For the simulation in Section 9, the proposed algorithm converges in only two iterations. Therefore, the total computation complexity of the proposed algorithm is ( stands for the number of iterations). Since the initialization takes up a large proportion of computation time in each iteration, the proposed algorithm suffers from a higher computational complexity.

7. CRB of Joint Polarization and DOA Estimation

Cramer-Rao bound (CRB) gives a lower bound of unbiased parameter estimation. In this section, the multiparameter estimation CRB is derived. For the sake of simplicity, the signal covariance matrix is assumed to be known, and the noise is normalized as one. parameters are contained in the covariance matrix , that is, elevation parameters, azimuth parameters, and polarization parameters, respectively. In this paper, the polarization parameters of the incident signals are assumed to be known. Only parameters remain to be estimated. The vector parameters to be estimated are represented as

The CRB of joint polarization parameters and angle parameters is defined as

The Fisher information matrix (FIM) for the parameter is given by

where is the block matrix of elevation estimator and is the block matrix of azimuth estimator. The entry of th row and th column of FIM is represented as

where is the trace of matrix and is the partial derivative of matrix .

8. Discussion

The proposed algorithm is able to operate on very irregular shapes, which is not following a circular geometry. However, some requirements have to be satisfied. The design of the array is discussed in two categories: the algorithm has an analytic solution; the algorithm does not have an analytic solution.

( 1) The Algorithm Has an Analytic Solution. Based on the design of array in Figure 2 and (12), (13), it can be seen that the distance vector is perpendicular to . This is one requirement for obtaining the analytic solution. Another requirement is that the distance vector or is parallel or perpendicular to the coordinate axis. Then, the analytic solution can be obtained. For example, the design of the spherical conformal array is given as follows.

The structure of spherical conformal array is shown in Figure 7. The elements are arranged on the surface of the spherical conformal array. constitute array 1, constitute array 2, constitute array 3, and constitute array 4. The elements of the array are divided into two subarrays. Subarray 1 consists of array 1 and array 2, and subarray 2 consists of array 3 and array 4. The distance vector between array 1 and array 2 is . The distance vector between array 3 and array 4 is . As shown in Figures 7 and 8, the distance vector is perpendicular to . Also, the distance vector is parallel to -axis.

Figure 7: The structure of spherical conformal array.

(a) The distance vector

(b) The distance vector

Figure 8: The schematic diagram of the distance vector.

( 2) The Algorithm Does Not Have an Analytic Solution. The only requirement is that is not parallel to as shown in Figure 2. In other words, the subarray 1 is not parallel to subarray 2 as shown in Figure 6. According to (61) and (62), the analytic solution of the proposed algorithm does not exist. Thus, the iterative numerical approximation method is used to obtain the optimal solution.

For very irregular shapes array, if we can find two subarrays as mentioned above and the two distance vectors are not parallel to each other, then the proposed algorithm can be used for DOA estimation.

9. Simulation Results

In this section, we demonstrate the performance of the proposed algorithm via numerical simulation. Taking a cylindrical conformal array as an example, we use the array structure in Figure 2 for simulation. The number of sensors is 16; that is, . Without loss of generality, , ; , . The sensor patterns are and in the global coordinate. and are the elevation and azimuth of the th incident signal in the th local coordinate, respectively. The details regarding how the pattern is transformed from the local coordinate to the global coordinate can be found in [8, 20]. 500 independent trials are considered. Root mean square error (RMSE) which indicates the performance of the proposed algorithm is defined as

where and are the elevation and azimuth estimators of the th incident signal in the th trial. The ESPRIT algorithm proposed in [18] is simulated in the same scenarios (the ESPRIT algorithm is used for the covariance matrix rather than the four-order cumulant).

For the following experiments, the COMFAC algorithm is used to fit the three-way array. The initialization and fitting of the COMFAC are done in compressed space. The Tucker3 three-way model is used in data compression [29].

In the first experiment, we show how the success rate and RMSE changed with different SNR. Here, a successful experiment is defined as the experiment with estimation error of less than 1 degree. The elevation and azimuth angles of two far-field narrow incident signals are and , respectively. 500 snapshots are used in this experiment. The success rate and RMSE are shown in Figures 9(a) and 9(b), respectively. Figure 9(a) shows that the success rate of the proposed PARAFAC algorithm is higher than that of the ESPRIT algorithm when the SNR is low. In addition, the RMSE of the proposed algorithm is much smaller than that of the ESPRIT algorithm at high SNR. As the SNR increases, the RMSE of the proposed algorithm approximates to the CRB. Four covariance matrices are used to estimate DOA in the proposed algorithm. However, only two covariance matrices are used for the ESPRIT algorithm. The proposed algorithm utilizes more data information than the ESPRIT algorithm, which leads to better performance.

(a) The success rate versus SNR

(b) The RMSE versus SNR

Figure 9: The estimation performance versus SNR.

We plot the curves of success rate and RMSE versus number of snapshots in Figures 10(a) and 10(b). SNR is 0 dB and 10 dB in Figures 10(a) and 10(b), respectively. Other simulation conditions are the same as those in the first experiment. Figure 10(a) shows that the success rate of the proposed algorithm is higher than that of the ESPRIT algorithm. In particular, the success rate of source 2 estimated by the ESPRIT algorithm is lower than others. As shown in Figure 10(b), the RMSE of the proposed algorithm is much smaller than that of the ESPRIT algorithm at the same number of snapshots. The RMSE approaches to the CRB as the number of snapshots increases. The RMSE of source 2 with ESPRIT algorithm is relatively larger which is mainly caused by the fact that the ESPRIT algorithm only uses the autocovariance matrices of the array. However, the proposed algorithm achieves higher accuracy by using both autocovariance and cross-covariance matrices of the array.

(a) The success rate versus number of snapshots

(b) The RMSE versus number of snapshots

Figure 10: The estimation performance versus number of snapshots.

Figure 11 shows the resolution probability versus SNR for both the proposed PARAFAC algorithm and the ESPRIT algorithm when the sources are closely spaced (6 degrees separation). The elevation and azimuth of the incident signals are and , respectively. Other simulation conditions are same as those in the first experiment. The two incident signals are resolved if , are smaller than . Meanwhile, , are smaller than . and represent the estimated and real elevations for th incident signal, respectively. and represent the estimated and real azimuths for th incident signal, respectively [30]. Figure 11 shows that the resolution of ESPRIT algorithm outperforms the proposed algorithm in relatively low SNR (−2 dB to 2 dB). When SNR is greater than 2 dB, the proposed algorithm performs better. According to (25) and Theorem 2, to ensure the uniqueness of the PARAFAC model, must be guaranteed. In other words, the two incident signals could not be too close. When the SNR is low, the effect of noise is fairly remarkable. All the above reasons result in the relatively poor resolution performance of the proposed algorithm in low SNR. Nevertheless, it is a good compromise for having excellent estimation accuracy and robustness to noise.

Figure 11: The resolution probability versus SNR.

10. Conclusion

In this paper, a novel high accuracy 2D-DOA estimation algorithm for the conformal array is proposed. The ordinary DOA estimation algorithm cannot be used on conformal array because of the polarization diversity of the varying curvature. To avoid parameter pairing problem, the algorithm forms a PARAFAC model of covariance matrices in the covariance domain to estimate the 2D-DOA. The uniqueness condition of the PARAFAC model is derived. It can be generalized to other conformal array structure, that is, conical conformal array, and so forth. Computer simulation verifies the effectiveness of the proposed algorithm in terms of accuracy and robustness to noise.

Appendix

The KR product of two matrices and of an identical number of columns is given by

where “” represents the Kronecker product. Then, the definition of Kronecker product of two vectors and is given by

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported in part by the National Science Foundation of China under Grant 61201410 and in part by Fundamental Research Focused on Special Fund Project of the Central Universities (Program no. HEUCF130804).

References

K. M. Tsui and S. C. Chan, “Pattern synthesis of narrowband conformal arrays using iterative second-order cone programming,” IEEE Transactions on Antennas and Propagation, vol. 58, no. 6, pp. 1959–1970, 2010.

M. Comisso and R. Vescovo, “Fast co-polar and cross-polar 3D pattern synthesis with dynamic range ratio reduction for conformal antenna arrays,” IEEE Transactions on Antennas and Propagation, vol. 61, pp. 614–626, 2013.

W.-J. Zhao, L.-W. Li, E.-P. Li, and K. Xiao, “Analysis of radiation characteristics of conformal microstrip arrays using adaptive integral method,” IEEE Transactions on Antennas and Propagation, vol. 60, no. 2, pp. 1176–1181, 2012.

A. Elsherbini and K. Sarabandi, “ENVELOP antenna: a class of very low profile UWB directive antennas for radar and communication diversity applications,” IEEE Transactions on Antennas and Propagation, vol. 61, no. 3, pp. 1055–1062, 2013.

B. R. Piper and N. V. Shuley, “The design of spherical conformal antennas using customized techniques based on NURBS,” IEEE Antennas and Propagation Magazine, vol. 51, no. 2, pp. 48–60, 2009.

J. L. Gómez-Tornero, “Analysis and design of conformal tapered leaky-wave antennas,” IEEE Antennas and Wireless Propagation Letters, vol. 10, pp. 1068–1071, 2011.

T. Milligan, “More applications of euler rotation,” IEEE Antennas and Propagation Magazine, vol. 41, no. 4, pp. 78–83, 1999.

B.-H. Wang, Y. Guo, Y.-L. Wang, and Y.-Z. Lin, “Frequency-invariant pattern synthesis of conformal array antenna with low cross-polarisation,” IET Microwaves, Antennas and Propagation, vol. 2, no. 5, pp. 442–450, 2008.

L. Zou, J. Laseby, and Z. He, “Beamformer for cylindrical conformal array of non-isotropic antennas,” Advances in Electrical and Computer Engineering, vol. 11, no. 1, pp. 39–42, 2011.

L. Zou, J. Lasenby, and Z. He, “Beamforming with distortionless co-polarisation for conformal arrays based on geometric algebra,” IET Radar, Sonar and Navigation, vol. 5, no. 8, pp. 842–853, 2011.

D. W. Boeringer, D. H. Werner, and D. W. Machuga, “A simultaneous parameter adaptation scheme for genetic algorithms with application to phased array synthesis,” IEEE Transactions on Antennas and Propagation, vol. 53, no. 1, pp. 356–371, 2005.

Y. Y. Bai, S. Xiao, C. Liu, and B. Z. Wang, “Optimisation method on conformal array element positions for low sidelobe pattern synthesis,” IEEE Transactions on Antennas and Propagation, vol. 61, no. 4, pp. 2328–2332, 2013.

K. Yang, Z. Zhao, J. Ouyang, Z. Nie, and Q. H. Liu, “Optimisation method on conformal array element positions for low sidelobe pattern synthesis,” IET Microwave, vol. 6, no. 6, pp. 646–652, 2012.

R. Karimizadeh, M. Hakkak, A. Haddadi, and K. Forooraghi, “Conformal array pattern synthesis using the weighted alternating reverse projection method considering mutual coupling and embedded-element pattern effects,” IET Microwave, vol. 6, no. 6, pp. 621–626, 2012.

L. I. Vaskelainen, “Constrained least-squares optimization in conformal array antenna synthesis,” IEEE Transactions on Antennas and Propagation, vol. 55, no. 3, pp. 859–867, 2007.

J. He, M. O. Ahmad, and M. N. S. Swamy, “Near-field localization of partially polarized sources with a cross-dipole array,” IEEE Transactions on Aerospace and Electronic Systems, vol. 49, no. 2, pp. 859–867, 2013.

X. Yuan, “Estimating the DOA and the polarization of a polynomial-phase signal using a single polarized vector-sensor,” IEEE Transactions on Signal Processing, vol. 60, no. 3, pp. 1270–1282, 2012.

Z.-S. Qi, Y. Guo, and B.-H. Wang, “Blind direction-of-arrival estimation algorithm for conformal array antenna with respect to polarisation diversity,” IET Microwaves, Antennas and Propagation, vol. 5, no. 4, pp. 433–442, 2011.

L. Zou, J. Lasenby, and Z. S. He, “Direction and polarization estimation using polarized cylindrical conformal arrays,” IET Signal Processing, vol. 6, no. 5, pp. 395–403, 2012.

W. J. Si, L. T. Wan, L. T. Liu, and Z. X. Tian, “Fast estimation of frequency and 2-D DOAs for cylindrical conformal array antenna using state-space and propagator method,” Progress in Electromagnetics Research, vol. 137, pp. 51–71, 2013.

N. D. Sidiropoulos, G. B. Giannakis, and R. Bro, “Blind PARAFAC receivers for DS-CDMA systems,” IEEE Transactions on Signal Processing, vol. 48, no. 3, pp. 810–823, 2000.

N. D. Sidiropoulos, R. Bro, and G. B. Giannakis, “Parallel factor analysis in sensor array processing,” IEEE Transactions on Signal Processing, vol. 48, no. 8, pp. 2377–2388, 2000.

J. B. Kruskal, “Rank decomposition, and uniqueness for 3-way and Nway arrays,” in Multiway Data Analysis, pp. 8–18, 1988.

J. B. Kruskal, “Three-way arrays: rank and uniqueness of trilinear decompositions, with application to arithmetic complexity and statistics,” Linear Algebra and Its Applications, vol. 18, no. 2, pp. 95–138, 1977.

R. Bro, “PARAFAC. Tutorial and applications,” Chemometrics and Intelligent Laboratory Systems, vol. 38, no. 2, pp. 149–171, 1997.

P. K. Hopke, P. Paatero, H. Jia, R. T. Ross, and R. A. Harshman, “Three-way (PARAFAC) factor analysis: examination and comparison of alternative computational methods as applied to ill-conditioned data,” Chemometrics and Intelligent Laboratory Systems, vol. 43, no. 1-2, pp. 25–42, 1998.

R. Bro, Multi-Way Analysis in the Food Industry: Models, Algorithms, and Applications, University of Amsterdam (NL) & Royal Veterinary and Agricultural University (DK), Amsterdam, The Netherlands, 1998.

H. A. L. Kiers and W. P. Krijnen, “An efficient algorithm for PARAFAC of three-way data with large numbers of observation units,” Psychometrika, vol. 56, no. 1, pp. 147–152, 1991.

N. D. Sidiropoulos, “COMFAC: Matlab Code for LS Fitting of the Complex PARAFAC Model in 3-D,” 1998, http://www.ece.umn.edu/~nikos/comfac.m

P. Stoica and A. B. Gershman, “Maximum-likelihood DOA estimation by data-supported grid search,” IEEE Signal Processing Letters, vol. 6, no. 10, pp. 273–275, 1999.