# 2D-DOA and Mutual Coupling Estimation in Vehicle Communication System via Conformal Array

Authors: Yan Zou; Hong Xie; Liangtian Wan; Guangjie Han; Wei Li
Journal: Mobile Information Systems
 (2015)

---

## Abstract

Many direction-of-arrival (DOA) estimation algorithms have been proposed recently. However, the effect of mutual coupling among antenna elements has not been taken into consideration. In this paper, a novel DOA and mutual coupling coefficient estimation algorithm is proposed in intelligent transportation systems (ITS) via conformal array. By constructing the spectial mutual coupling matrix (MCM), the effect of mutual coupling can be eliminated via instrumental element method. Then the DOA of incident signals can be estimated based on parallel factor (PARAFAC) theory. The PARAFAC model is constructed in cumulant domain using covariance matrices. The mutual coupling coefficients are estimated based on the former DOA estimation and the matrix transformation between MCM and the steering vector. Finally, due to the drawback of the parameter pairing method in Wan et al., 2014, a novel method is given to improve the performance of parameter pairing. The computer simulation verifies the effectiveness of the proposed algorithm.


---

## Body

## 1. Introduction

Intelligent transportation systems (ITS) have emerged as an effective way of improving the performance of transportation systems and enhancing travel security. The term connected vehicles refers to applications, services, and technologies that connect a vehicle to its surroundings. With the rapid development of 5G, the peak data rate will likely be in the range of tens of Gbps, which is suitable for real-time communication between vehicles and BS (base station). In order to minimize the cost of the vehicles and reduce the friction resistance between the surface of the vehicles and atmosphere, the conformal array antennas installed in the vehicles should be a good choice.Direction-of-arrival (DOA) estimation plays an important role in array signal processing, which has been widely used in radar, sonar, and smart antenna [1–3]. Multiple signal classification- (MUSIC-) based algorithm [4] and estimation of signal parameters via rotational invariance techniques- (ESPRIT-) based algorithm [5] are two types of DOA estimation algorithms which can achieve superresolution. However, the effect of mutual coupling will reduce the performance of DOA estimation severely.In order to solve this problem, many algorithms have been proposed to deal with it. An effective compensation method for the effect of mutual coupling in uniform circular arrays (UCAs) employed for two-dimensional (2D) DOA estimations was proposed [6]. A DOA estimation algorithm for mixed signals with unknown mutual coupling was proposed in [7]. The noncoherent (uncorrelated or partially correlated) signals were firstly estimated with unknown mutual coupling; then the mutual coupling coefficients can be obtained by these estimation. Finally, the noncoherent signals were eliminated and the effect of mutual coupling was compensated. The bias in uniform linear array (ULA) and general linear array caused by effect of mutual coupling was studied in [8]. The DOA estimation for noncoherent and coherent signals has been proposed in [9], which is used for patient’s localization.Parallel factor (PARAFAC) analysis attracted the attention of researchers when it was originally introduced in array signal processing in 2000 [10, 11]. It has been widely used for low-rank decomposition of three-way and higher way array. Based on PARAFAC and cumulant, a new array structure was designed for 2D-DOA estimation [12]. The DOA and polarization estimation for a single electromagnetic vector-sensor were acquired by using PARAFAC [13], and the proposed algorithm performs better than the ESPRIT algorithm.Conformal array is the array integrated on surface of object [14, 15]. Some DOA estimation algorithms have been proposed for conformal array. The propagator method and three rotational invariance relationships have been used for fast DOA estimation for cylindrical conformal array [16]. A ultrawideband DOA estimation algorithm has been proposed for DOA estimation based on spatial baseline method and multiple subarray technique [17]. The manifold separation technique has been used for DOA estimation of wireless sensor network with arbitrary array configuration [18, 19]. A real-time DOA estimation algorithm has been proposed in [20]. The DOA estimation is accomplished via ESPRIT and the elimination of the effect of mutual coupling [21]. Based on spatial baseline technique, the DOA can be estimated with low computational complexity. However, only one source can be estimated [22]. In [23], the authors have extended the PARAFAC theory for frequency and DOA estimation. By using space-time matrix technique, the frequency and DOA estimation have been completed as well. However, the positions of four guiding elements have to be known as a prior [24].In this paper, a novel DOA and mutual coupling coefficient estimation algorithm for conformal array is proposed in vehicle communication system. In realistic application, the effect of mutual coupling should not be ignored. Based on the selection matrix, the spectial mutual coupling matrix (MCM) is constructed. The effect of mutual coupling can be eliminated by using the spectial MCM. Then two rotational invariant matrices are constructed in order to estimate the DOA of incident signals based on conformal array. Based on PARAFAC theory, the model is constructed in cumulant domain using covariance matrices. The mutual coupling coefficient is estimated finally via matrix transformation.


## 2. The Snapshot Data Model

The algorithm performance would deteriorate significantly because of the mutual coupling among different elements. In order to estimate DOA estimation accurately, the effect of mutual coupling should be considered in the snapshot data model of the conformal array. Due to the varying curvature of carriers, each element possesses different patterns [18]. The incident signals are shown in Figure 1(a) with elevation θ and azimuth φ. We can give the form of steering vector as follows:(1)aθ,φ=h1e-j2πp1·u/λ,h2e-j2πp2·u/λ,…,h2Me-j2πp2M·u/λT,(2)hm=gmθ2+gmφ21/2kθ2+kφ21/2cos⁡θmgk=gmplcos⁡θmgk=gm·pl=gmθkθ+gmφkφ.Figure 1
(a) The incident signalu. (b) The ith element’s response.
(a)
(b)The notation of alphabet mentioned above can be found in [25]. Taking (2) into (1), we have (3)aθ,φ=aθθ,φ+aφθ,φ,where (4)aθθ,φ=g1θkθe-j2πp1·u/λ,g2θkθe-j2πp2·u/λ,…,g2Mθkθe-j2πp2M·u/λT,aφθ,φ=g1φkφe-j2πp1·u/λ,g2φkφe-j2πp2·u/λ,…,g2Mφkφe-j2πp2M·u/λT.WhenP incident signals arrive in the array, the data matrix of array output can be expressed as(5)Xn=G⊙CASn+Nn=CGθ⊙AθKθ+Gφ⊙AφKφSn+Nn=CBSn+Nn,(6)Gθ=gθθ1,φ1,gθθ2,φ2,…,gθθP,φP,(7)Gφ=gφθ1,φ1,gφθ2,φ2,…,gφθP,φP,(8)Aθ=aθθ1,φ1,aθθ2,φ2,…,aθθP,φP,(9)Aφ=aφθ1,φ1,aφθ2,φ2,…,aφθP,φP,(10)Kθ=diag⁡k1θ,k2θ,…,kPθ,(11)Kφ=diag⁡k1φ,k2φ,…,kPφ,(12)Sn=s1n,s2n,…,sPnT,(13)Wn=w1n,w2n,…,w2MnT,where G denotes the pattern matrix and A denotes the manifold matrix. ⊙ stands for the Kronecker product. k1θ,k2θ,…,krθ and k1φ,k2φ,…,krφ construct Kθ and Kφ, respectively. The ith signal’s polarization parameters are kiθ and kiφ, respectively. The signal vector is Sn. The covariance matrix of noise Wn is(14)EWnWnH=Q=σ2I.·H denotes conjugate transpose of matrix ·. I is the identical matrix.For simplicity, only the effect of mutual coupling among the same uniform linear array (ULA) is taken into consideration. The mutual coupling matrix (MCM)C is modeled as a banded symmetric Toeplitz matrix. c=c1,c2,…,cq,0,…,0 is the first row of C. The element of vector c satisfying 0<cq<cq-1<⋯<c1=1. This matrix is formulated as(15)C=Toeplitzc=Toeplitzc1,c2,…,cq,0,…,0=1c1⋯cq⋯0c11c1⋯⋱0⋮c11⋱⋯cqcq⋯⋱⋱c1⋮0⋱⋯c11c10⋯cq⋯c11.CollectingN snapshots, we have the matrix form of (5) as(16)X=CBS+W.There is a characteristic called “shadow effect” which most conformal arrays possess. It means that when the incident signal impinges on the array, not all elements could receive it. In order to solve this problem, the subarray divided technique (SDT) proposed in [18, 19] is adopted in this paper. The whole conformal array consists of several subarrays, and each subarray covers a certain range of angle. For the single curvature and symmetry of the cylinder, the parameter estimation mechanism and array design of each subarray are identical. Thus only one subarray is considered in this paper.


## 3. The Mechanism of Mutual Coupling Elimination

The design of cylindrical conformal array is shown in Figure2, the distance between two elements in the identical plane is λ/2, and the distance between two neighbouring planes is λ/2. The radius of the cylinder is 5λ.Figure 2
The structure of cylindrical conformal array.According the form of the MCMC, a mutual coupling elimination mechanism is introduced based on instrumental elements [26]. To show this mechanism in an intuitive way, we give an example based on 1~M elements. The position vector pm of the mth element can be expressed as(17)pm=xme→x+yme→y+zme→z.Since M elements are arranged on the same generatrix, the x-coordinate has the relationship x1=x2=⋯=xM=x as well as y1=y2=⋯=yM=y and h1=h2=⋯=hM. The manifold matrix of this subarray is B1=b1θ1,b1θ2,…,b1θP, and the steering vector is (18)b1θp=Hp,Hpvp1,…,HpvpM-1T,where vp=exp⁡-j2πdcos⁡θp/λ and(19)Hp=h1e-j2πxsin⁡θpcos⁡φp+ysin⁡θpsin⁡φp/λ.The instrumental elements are set as the first and lastq elements. Although the array aperture is reduced, the effect of mutual coupling is eliminated as well. We construct a selection matrix that has the following form: F=0M-2q×qIM-2q0M-2q×q, and the received data X- of center array can be given by(20)X-=PCB1S+PW=C-B1S+PW,where the new M-×MC- is expressed as(21)C-=PC=cq⋯1⋯cq0⋯00cq⋯1⋯cq⋯0⋮⋱⋱⋯⋱⋯⋱⋮0⋯0cq⋯1⋯cq.By using instrumental elements, we can get an important equation as follows:(22)C-b1θp=Hpcqvp-q+⋯+c1vp-q+1+c1vp+⋯+cqvpqHpcqvp1-q+⋯+c1+vp+c1vp2+⋯+cqvpq+1⋮HpcqvpM-q-2+⋯+c1vpM-2+vpM-1+c1vpM+⋯+cqvpM+q-2HpcqvpM-q-1+⋯+c1vpM-1+vpM+c1vpM+1+⋯+cqvpM+q-1=Hpcqvp-q+⋯+c1vp-1+1+c1vp+⋯+cqvpq1vp⋮vpM--1=2Hp∑q′=1qcq′cos⁡2q′πsin⁡θpdλ+1b-1θp=cθpb-1θp,where b-1θp is the ideal steering vector of the center array (q+1~M-q elements), M-=M-2p, and cθp is a scalar which can be written as (23)cθp=2∑q′=1qcq′cos⁡2q′πsin⁡θpdλ+1.Based on the (21), the new covariance matrix can be expressed as (24)R-x=EX-X-H=C-B1R-sB1HC-H+σ2I=B-1DR-sDHB-1H+σ2I,where(25)B-1=b-1θ1,b-1θ2,…,b-1θP,D=diag⁡cθ1,cθ2,…,cθP.Then (24) can be written as(26)R-x=B-1RB-1H+σ2I,where (27)R=diag⁡σ12cθ12,σ22cθ22,…,σP2cθP2and σp2 stands for the power of pth incident signal. A very interesting result is shown in (27) where the mutual coupling coefficients are coming into the novel signal covariance matrix completely. When cθp≠0, p=1,2,…,P, rank⁡R=P, the elimination of effect of mutual coupling can be accomplished. When cθp=0, p=1,2,…,P, we have(28)cθpb-θp=0.Equation (33) implies that the center array is blind at some particular angles. It could not receive incident signal from certain directions. Thus, these angles are called “blind angles.” Fortunately, for given mutual coefficients, cθ-l is a continuous function, the probability of cθ-l=0 is approximate to zero, which means that “blind angles” phenomenon rarely happens in practice. More details can be found in [27].As shown in Figure3, there are two rotational invariance relations in the designed array. The ESPRIT algorithm which is similar to the algorithm described in [17] can be used for DOA estimation. The parameter pairing between elevation and azimuth has to be considered. However, the algorithm is not suitable for estimating conformal array. Based on PARAFAC theory, a novel high-accuracy 2D-DOA estimation algorithm is proposed in this paper for conformal array.Figure 3
The distance vector. (a) The distance vectorΔP1 between array 1 and array 2. (b) The distance vector ΔP2 between array 3 and array 4.
(a)
(b)


## 4. DOA Estimation Based on PARAFAC

When the elimination of mutual coupling is done, the decoupling between polarization and direction can be completed based on well array design. Define twoM--1×M- selection matrices:(29)P1=1⋯00⋮⋱⋮⋮0⋯10,P2=01⋯0⋮⋮⋱⋮00⋯1.q
+
1
~
M
-
q
-
1, q+2~M-q, q+1~M-q, and M+q+1~2M-q elements constitute array 1, array 2, array 3, and array 4, respectively. The distance vector between arrays 1 and 2 is ΔP1 and d1=ΔP1=λ/2. The distance vector between arrays 1 and 3 is ΔP2 and d2=ΔP2=λ/2, which is shown in Figure 3. Since the direction pattern is identical through one generatrix, the effect of polarization parameter can be ignored.The distance among different subarrays is shown in Figures2 and 3. In the global coordinate, the elevation and azimuth of ΔP1 and ΔP2 can be represented as θΔP1=0, φΔP1=π/2, and θΔP2=π/2, φΔP2=0, respectively.Assume thatX1, X2, X3, X4 represent the constructed data of array 1, array 2, array 3, and array 4, respectively. The origin of the coordinate is the reference point. So X1, X2, X3, X4 can be expressed as(30)X1=BS+N1,(31)X2=BΨ1S+N2,(32)X3=BΨ2S+N3,(33)X4=BΨ1Ψ2S+N4.The covariance matrices among the received data are(34)R1=EX1X1H=BRsBH+Q1,R2=EX2X1H=BΨ1RsBH+Q2,R3=EX3X1H=BΨ2RsBH+Q3,R4=EX4X1H=BΨ1Ψ2RsBH+Q4,where B=P1B-1, Rs=diag⁡s12,…,sr2 represents the signal covariance matrix, and the noise covariance matrices are Q1~Q4, respectively:(35)Φ1=diag⁡exp⁡-jω11,exp⁡-jω12,…,exp⁡-jω1r,Φ2=diag⁡h3θ1,φ1h1θ1,φ1exp⁡-jω21,h3θ2,φ2h1θ2,φ2exp⁡-jω22,…,h3θr,φrh1θr,φrexp⁡-jω2r,ω1i=2πλΔP1·ui=2πd1λsin⁡θΔP1cos⁡φΔP1sin⁡θicos⁡φi+sin⁡θΔP1sin⁡φΔP1sin⁡θisin⁡φi+cos⁡θΔP1cos⁡θi,ω2i=2πλΔP2·ui=2πd2λsin⁡θΔP2cos⁡φΔP2sin⁡θicos⁡φi+sin⁡θΔP2sin⁡φΔP2sin⁡θisin⁡φi+cos⁡θΔP2cos⁡θi,where θΔPi and φΔPi (i=1,2) represent the elevation and azimuth of the distance vector in the global coordinate, respectively.On the basis of the PARAFAC theory [11], the m×m×4 three-way array of the cylindrical conformal array is constructed by (34):(36)R:,:,1R:,:,2R:,:,3R:,:,4=R1R2R3R4=BRsBHBΨ1RsBHBΨ2RsBHBΨ1Ψ2RsBH+Q~,where Ri (i=1,…,4) are the sample covariance matrices, respectively, and Q~ represents the noise in real observation. Let C=BH, and, based on the definition of Khatri-Rao product, (36) can be transformed as (37)R=D⊙BC+Q~,(38)RX=CT⊙DBT+Q~X,(39)RY=B⊙CTDT+Q~Y,(40)D=Λ-1RsΛ-1Ψ1RsΛ-1Ψ2RsΛ-1Ψ1Ψ2Rs,where Λ-1Rs is the row vector constructed by the diagonal entries of the diagonal matrix Rs.On the basis of noisy observation, problem (37) can be transformed into solving a least square problem:(41)minD,B,CR-D⊙BCF2.The principle of alternating least squares (ALS) can be used to fit the problem in (41). In the noiseless condition, ALS can be used to solve the matrices B, C, and D which constructed the three-way array R. The least square estimation of matrix C can be expressed as (42)C=arg⁡minCR-D⊙BCF2.Similarly, the matrices B and D can be expressed as (43)BT=arg⁡minBRX-CT⊙DBTF2,DT=arg⁡minDRY-B⊙CTDTF2.In the iterative procedure, given matrices B and D, the matrix C can be represented as (44)C=D⊙B†R.The expression of matrices BT and DT is (45)BT=CT⊙D†RX,DT=B⊙CT†RY,where ·† denotes the pseudoinverse of matrix ·.The matrixD can be acquired by the TALS algorithm. ω1i and ω2i can be calculated by matrix D:(46)ω1i=-12angleD2iD1i+angleD4iD3i· stands for the matrix and · stands for the absolute value operation, where Dji represents the jth row of D. Because h1, h3, and h4 are real numbers, D3i/D1i and D4i/D2i are squared to solve the ambiguity caused by the positive and negative values of h1, h3, and h4:(47)ω2i=-12angleh3θi,φih1θi,φiexp⁡-jω2i2=-12angleexp⁡-j2ω2i=-12angleD3iD1i2=-12angleD4iD2i2,and then (48)ω2i=-14angleD3iD1i2+angleD4iD2i2.The elevation θi and azimuth φi of ith incident signal can be expressed as(49)θi=arccosλω1i2πd2=arccos2ω1iπ,(50)φi=arccosλω2i2πd2sin⁡θi=arccos2ω2iπsin⁡θi.


## 5. Mutual Coupling Coefficient Estimation

Based on the 2D-DOA estimation of incident signals, the estimations of mutual coupling coefficients can be obtained. The covariance matrix ofX can be expressed as(51)RX=EXXH=CARsAHCH+σn2IN.The rank of CARsAHCH is P. Take the eigendecomposition of RX; the P big eigenvalues λ1,λ2,…,λP and M-P small eigenvalues λP+1=λP+2=⋯=λK=σn2 can be obtained, respectively. Their corresponding eigenvectors are u1,u2,…,uM. The signal subspace is spanned by CA, which is orthogonal to the noise subspace spanned by uKu+P+1,uKu+P+2,…,uM:(52)UnHCaθi=0,i=1,2,…,P.By considering the complex symmetric Toeplitz form of C, Caθ can be expressed as(53)Caθ=1c1⋯cN-1c11⋯cN-2⋮⋮⋮⋮cN-1cN-2⋯11v⋮vM-1=1v⋯vM-2vM-1v1+v2⋯vM-10v2v+v3⋯00⋮⋮⋱⋮0vM-2vM-3+vM-1⋯10vM-1vM-2⋯v11c1⋮cN-1=Tθc,where Tθ is the sum of two M×M0+1 matrices [28].Based on (52) and the estimations of θ1,θ2,…,θP, we have(54)UnHTθic=0,i=1,2,…,P.It can be seen that (54) is the linear equations of mutual coupling coefficients c. The coefficient matrix can be defined as(55)Q=UnHTθ1⋮UnHTθP.Then (50) can be written as(56)Qc=0,where Q is a PM-P×M0+1 matrix and Q=q1,q2,…,qP. Due to c(1)=1, we have(57)Qc=q1,q2,…,qM0+11c1⋮cM0=0.When PM-P≥M0-1, the least square solution is obtained as(58)c1,c2,…,cM0=-q2,…,qM0+1†q1.The mutual coupling coefficients estimation is completed.


## 6. Discussion

The estimators ofB, C, and D have the same column permutation matrix; that is, the ith column of the steering matrix B corresponds to the ith of the matrix D. Thus, the elevation and azimuth pair with each other automatically.WhenP≥2, the parameter pairing between ω1i  i=1,2,…,P and ω2i  i=1,2,…,P is a big problem in [21]. Here, we give another method to solve the parameter pairing problem in [21]. Theoretically, the eigenvector of Φ1, Φ2, and Φ3 is identical. However, the EVD of them are done independently. The order of eigenvector may be distinct. In order to obtain accurate DOA estimation, this problem must be solved. We construct another subarray: array 1 and array 5 constructed by M+q+2~2M-q elements. The distance vector is ΔP3. Based on the relationship between ΔP1, ΔP2, and ω1i, ω2i, the relationship between ΔP3 and ω3kk=1,2,…,P can be obtained similarly. As shown in Figure 4, ΔP3=ΔP1+ΔP2.Figure 4
The distance vector between distinct subarrays.For the identical source, we have(59)angleh2θp,φph1θp,φpexp-jω3p2=anglet3k2=angleexp-j2ω1i+angleexp-j2ω2i=anglet1i2+anglet2j2,where t3k is the eigenvalue of Φ3. Φ3 can be solved easily according to the solution procedure of Φ1 and Φ2. Thus, the parameter pairing problem can be transformed as a problem which minimizes(60)min⁡anglet3k2-anglet1i2+anglet2j2.Based on (60), the parameter pairing is completed finally.The steps of DOA estimation for cylindrical conformal array in the presence of mutual coupling can be summarized as follows:(1)
Use the selection matrixF=0M-2q×qIM-2q0M-2q×q to eliminate the effect of mutual coupling.(2)
Set array 1~array 4 to construct three subarrays.(3)
Construct the PARAFAC model based on (36).(4)
Update (44), (45), and the matrix D can be obtained by TALS algorithm.(5)
Calculateω1i and ω2i based on (46) and (47).(6)
Estimate the 2D-DOA for cylindrical conformal array based on (49) and (50).(7)
Calculate the covariance matrix of the data received by element1~M based on (51).(8)
Obtain the mutual coupling coefficients based on (58).


## 7. Simulation Results

In this section, we present numerical simulation results to illustrate the performance of the proposed algorithm. We compare the proposed algorithm with ESPRIT without mutual coupling (ESPRIT-NMC) [5] and the ESPRIT-E algorithm proposed in [21]. In this paper, root mean square error (RMSE) is calculated by the formula(61a)RMSE=1200∑η=1200θ^k,η-θ2+φ^k,η-φ2,(61b)RMSEc=1200c∑η=1200c^η-c×100%,where θ^k,η and φ^k,η are the estimate of the elevation and azimuth of the kth incident signal in the ηth trial. For the sake of simplicity, only the cylindrical conformal array is utilized for simulation, and the array structure is shown in Figure 2. The number of elements is 16; that is, m=8. Without loss of generality, k1θ=0.5, k1φ=0.5; k2θ=0.3, k2φ=0.7. The pattern transformation can be found in [25]. The azimuth and elevation of two incident signals are (100°, 60°) and (95°, 50°), respectively. The freedom of degree (FOD) of mutual coupling is 2, and c1=0.4500+i0.5362 and c2=0.2598-i0.1500. Since the instrumental elements are used, the elements actually used are 12; that is, m=6. For ESPRIT-NMC, the elements actually used are 12; that is, m=6. 200 Monte Carlo trials are considered in the simulations.The RMSE and successful probability of different algorithms against SNR are shown in Figures5(a) and 5(b), respectively. The snapshot number is 200. It can be seen from Figure 5(a) that the RMSE of the proposed algorithm is smaller than that of ESPRIT-E. Since the mutual coupling coefficient has been estimated, this result can be used to recover the original covariance matrix. The array aperture is extended compared with ESPRIT-E and ESPRIT-NMC. It can be seen from Figure 5(b) that the successful probability increase as SNR increases. The successful probability of ESPRIT-E is higher than that of ESPRIT-NMC. The proposed algorithm outperforms ESPRIT-E and ESPRIT-NMC.Figure 5
The performance against SNR (a). The RMSE against SNR (b). The successful probability against SNR.
(a)
(b)The RMSE and successful probability of different algorithms against snapshot number are shown in Figures5(a) and 5(b), respectively. SNR is fixed at 0 dB. It can be seen from Figure 6(a) that the RMSE of proposed algorithm is smaller than that of ESPRIT-E. When snapshot number is larger than 700, the RMSE of proposed algorithm approximates ESPRIT-NMC. The successful probability of proposed algorithm is higher than that of ESPRIT-E. However, the proposed algorithm does not reach 100% because of the low SNR (0 dB). Thus, the proposed algorithm has better performance compared with ESPRIT-E.Figure 6
The performance against snapshot number (a). The RMSE against snapshot number (b). The successful probability against snapshot number.
(a)
(b)The mutual coupling estimation versus input SNR is shown in Figure7. The results show that the mutual coupling has little influence on the proposed algorithm. Although some array apertures are used to eliminate the effect of mutual coupling, the proposed algorithm still has high estimation accuracy and approaches CRB at high SNR.Figure 7
The RMSE against SNR.


## 8. Conclusions

In this paper, a novel DOA and mutual coupling coefficient estimation algorithm for conformal array is proposed in ITS. By constructing the spectical mutual coupling matrix (MCM), the effect of mutual coupling can be eliminated via instrumental element method. Then the DOA of incident signals can be estimated based on parallel factor (PARAFAC) theory. The proposed algorithm can also be extended to other conformal arrays. However, only the mutual coupling in the identical linear array is considered in this paper. The mutual coupling between different generatrices should be studied. The Cramér-Rao bound is also an interesting problem to research, and we will focus on it in the future work.



---

*Source: 841341-2015-12-28.xml*