# The Limited Aperture Sparse Array for DOA Estimation of Coherent Signals: A Mutual-Coupling-Optimized Array and Coherent DOA Estimation Algorithm

Authors: Yu Zhang; Meng Yang; Xiaofei Zhang; Shengxinlai Han; Jinke Cao
Journal: Mathematical Problems in Engineering
 (2022)

---

## Abstract

The paper investigates DOA estimation of coherent Signals with the limited aperture sparse array. Mutual coupling between the sensors of the array cannot be ignored in practical radar with a limited aperture of array sensors, which will result in a degradation in the performance of Direction of Arrival (DOA) estimation. This paper proposes a Mutual-coupling-optimized array (MCOA) with a limited aperture in this scenario to reduce the mutual coupling effect. Firstly, we prove the sparse uniform linear array (SULA) has the smallest mutual coupling leakage when the array aperture and the number of sensors is determined. Secondly, we modify the spacing of the array sensors in SULA to make sure that the spacing between all array sensors and the reference sensor are coprime aiming to estimate DOA without spatial aliasing. Thirdly, we give an expression for the array element spacing arrangement with reduced mutual coupling leakage. Finally, the coherent signals are well resolved by the Sparse Bayesian Learning (SBL) algorithm. Numerous simulations are conducted to validate the advantages of the proposed array compared to several sparse arrays for estimating coherent signals in the presence of mutual coupling.


---

## Body

## 1. Introduction

The problem of Direction of Arrival (DOA) estimation has attracted a lot of attention in the fields of radar, sonar, navigation, and astronomy [1–7], where the antenna arrays are utilized for collecting the spatial sampling of incident signals. Scholars propose many DOA estimation algorithms based on the uniform linear array (ULA), such as multiple signal classification (MUSIC) [5], estimation of signal parameters via rotational invariance techniques (ESPRIT) [6], propagator method (PM) [8], and parallel factor (PARAFAC) technique [9].However, the above-given algorithm is predicated on the assumption that the incident signals are uncorrelated. The received signals are typically coherent and the rank of the covariance matrix is insufficient due to the impact of the transmission environment in actual applications. The aforementioned DOA estimation algorithm will be invalid at this time [10]. To tackle this problem, some decoherence algorithms are proposed to deal with coherent signals. The most representative method is the spatial smoothing (SS) method [11], which regards ULA as many subarrays with the same array flow type and then averages the covariance matrix of each subarray to obtain the full rank covariance matrix. Later, people proposed the Forward/backward spatial smoothing techniques [12], improved spatial smoothing techniques [13] on this basis of SS and make a series of improvement on these algorithms [14]. In [15], the authors were devoted to establishing more accurate conditions by studying the positive definiteness of smoothed target covariance matrix. There are also algorithms that reconstruct the covariance matrix of the received signal, such as SVD algorithms and the Toeplitz decoherence method [16]. These methods estimate coherent signals at the expense of array aperture, which reduces DOA estimation performance. The compressed sensing (CS) algorithms [17–19] can estimate DOA by exploiting the sparsity of the target in the spatial domain without taking into account the coherence of the signals. In [20], an iterative adaptive approach (IAA) is given for the beamforming design based on the sparsity. In [21], the Orthogonal Matching Pursuit (OMP) is used to recover the spare signal with high probability, but the accuracy of OMP is lower than that of MUSIC. In [22], both the Sparse Bayesian Learning (SBL) and the relevance vector machine (RVM) are proposed. The weakness of CS algorithms is that they are more complex than the previously described DOA estimating techniques. The traditional DOA estimation algorithms are generally considered based on ULA, but sparse linear arrays are seldom utilized.Recently, sparse arrays such as Nested arrays (NA) [23] and coprime arrays (CA) [24] have attracted wide attention because such sparse arrays can achieve OM2 degrees of freedoms (DOFs) with only M antenna sensors. Though the DOA estimation performance of NA is better than that of CA, the mutual coupling leakage of NA is much greater than that of CA due to the influence of the dense ULA subarray, which reduces performance in the presence of mutual coupling. Despite the array positions of CA and NA can be expressed in closed-form, their continuous degrees of freedom are not the greatest. In comparison to CA and NA, the minimum redundant array (MRA) [25] has the most continuous degrees of freedom, allowing it to use more virtual arrays. However, MRA lacks a closed-form expression for the locations of its sensors, and its array design requires a significant amount of complicated calculations.The aperture of the array is usually limited in most applications, and the number of array elements is fixed. Because the spacing between the array elements is relatively close, the mutual coupling effect cannot be ignored. ULA and traditional sparse arrays will fail to estimate the DOA of coherent signals and the research of DOA estimation in this situation is relatively few. Though there are some methods [26–28] proposed to mitigate the mutual coupling effects by utilizing a complex mutual coupling model, these methods estimate mutual coupling coefficients at the cost of increased complexity and decreased degree of freedom (DOFs). Therefore, it is a good choice to consider how to reduce the mutual coupling effect when designing the array. Under the restrictions of a set array aperture and a number of array sensors, this paper determines the array design approach with the least mutual coupling leakage, and we further propose a Mutual-coupling-Optimized Array (MCOA) based on the theory of estimating DOA without spatial aliasing [29, 30]. To estimate the DOA of coherent signals, we use the sparse Bayesian learning-based (SBL) compressed sensing algorithm. In particular, we summarize our main contributions as follows:(1)
We propose a mutual-coupling-optimized array under the restriction of a fixed number of sensors and fixed array aperture. Then, we prove that the mutual coupling leakage of the suggested array is smaller than that of conventional sparse arrays and that it can estimate DOA without spatial aliasing.(2)
We employ the SBL algorithm to estimate coherent DOA for the proposed array and compare the SBL algorithm with other algorithms to demonstrate that the estimation performance of the SBL algorithm is better than other algorithms including OMP and IAA.The remainder of this paper is given as follows: we provide the mathematical model of the sparse array and the definition of mutual coupling matrix in Section2. In Section 3, we present how to design the mutual coupling optimized array with a limited aperture. Section 4 introduces the specific steps of the sparse Bayesian learning algorithm. Section 5 analyses the CRB of the array in this context. Section 6 verifies the theoretical performance of the proposed array through simulation analysis while Section 7 concludes this paper.Notations: scalars, vectors, matrices, and sets are denoted by lowercase lettersa, lowercase letters in boldface a, uppercase letters in boldface A, and letters in blackboard boldface A, respectively. AT, A∗, and AH are the transpose, complex conjugate, and complex conjugate transpose of A. Tr⋅ denotes the trace operator for a matrix. ⋅F represents the Frobenius norm and diag⋅ represents the matrix formed by the diagonal elements of the matrix. Ai,j denotes the i,j entry of A. gcdn1,n2,⋯,nM denotes the greatest common divisor of the elements.


## 2. Mathematical Model

Consider a sparse array consisting ofM sensors as shown in Figure 1, and the position of the ith sensor is denoted by zid with d=λ/2, zi represents the distance of the i- th sensor relative to the reference sensor and λ stands for the wavelength. Suppose that there are K far-field narrowband coherent signals from different directions θ=θ1,θ2,⋯,θKb2−4ac with powers p=σ12,σ22,⋯,σK2 impinge on this sparse array. The received signal of the array can be expressed as follows [2]:(1)X0=AS+N,where A=aθ1,aθ2,⋯,aθK represents the direction matrix and aθi=1,ej2πz2dsinθi/λ,⋯,ej2πzMdsinθi/λT denotes the direction vector of the ith signal. λ is the wavelength of the signal. S=α1,α2,⋯,αKs0∈ℂM×J means the narrowband coherent signals with J snapshots, where s0 is the generate signals and αi stands for the complex constant. N∈ℂM×J represents the additive white Gaussian noise vector with noise variance σ2.Figure 1
Array model.There is coupling between the array elements. The received signal model is expressed as follows [28]:(2)X=CAS+N,where C is the mutual coupling matrix. The mutual coupling matrix can be modelled as a B-banded symmetric Toeplitz matrix according to the assumption in the following equation [28]:(3)Cp,q=0,zp−zq>Bczp−zqzp−zq≤B,where 1>c1>⋯>cB>0, c1=c0ejπ/3, cs=c1e−js−1/8/s, s∈0,B. c0 is the mutual coupling constant and B denotes the maximum spacing of sensor pairs with mutual coupling. In addition, the mutual coupling is evaluated by the coupling leakage, i.e.,(4)γ=C−diagCFCF.


## 3. Mutual-Coupling-Optimized Array with Limited Aperture

In this section, we first show that the mutual coupling leakage of the sparse and uniform linear array is the smallest under the condition of finite aperture and number of elements. Then, we discussed how to change the position of the array elements, so that the mutual coupling leakage is still small, and there is no spatial aliasing in DOA estimation.


### 3.1. A Minimum Mutual Coupling Leakage Array

Considering theM- element array with an array aperture of N, we propose an array which has no spatial aliasing in DOA estimation and has the smaller mutual coupling leakage than most sparse arrays.Lemma 1.
For anM element array with an array aperture of N, the mutual coupling leakage is minimized if and only if the array elements are equally spaced.(5)zi+1d−zid=zid−zi−1d=NM−1,i=2,3,⋯,M−1.Proof.
According to equation (3), the mutual coupling matrix can be expressed as follows:(6)C=c0cz2−z1⋯czM−z1cz2−z1c0⋯czM−z2⋮⋮⋱⋮czM−z1czM−z2⋯c0.
Then, the expression of the coupling leakage can be calculated as follows:(7)γ=C−diagCFCF=2∑i=1M−1∑j=i+1Mczj−zi2Mc02+2∑i=1M−1∑j=i+1Mczj−zi2,where(8)∑i=1M−1∑j=i+1Mczj−zi2=c12∑i=1M−1∑j=i+1M1zjd−zid2=c12S,where S=∑i=1M−1∑j=i+1Mzjd−zid−2. The value of zi,i=1,2,⋯,M corresponding to the minimum value of S is the position of each array sensor when the mutual coupling leakage is minimum.
Denote theM−1 array spacings as x1,x2,⋯,xM−1 respectively. Then,(9)S=∑i=1M−1∑j=i+1Mzjd−zid−2=∑i=1M−1zi+1d−zid−2+∑i=1M−2zi+2d−zid−2+⋯+∑i=11zMd−zid−2=∑i=1M−1xi−2+∑i=1M−2xi+xi+1−2+⋯+∑i=11xi+xi+1+⋯+xi+M−2−2.
Calculate the minimum value ofS by Lagrange multiplier method.(10)gx1,⋯,xM−1,μ=S+μx1+x2+⋯+xM−1−N.
Take the partial derivative of each variable in the function and set the result equal to zero.(11)∂gx1,⋯,xM−1,μ∂x1=0⋮∂gx1,⋯,xM−1,μ∂xM−1=0∂gx1,⋯,xM−1,μ∂μ=0.
There is an extreme point in (10) when x1=x2=⋯=xM−1=N/M−1 and the cost function S has a minimum value at this time. Therefore, the value of zi,i=1,2,⋯,M corresponding to the minimum value of S is the position of each array element in the array when the mutual coupling leakage is minimum. From the above-given discussion, it can be seen that the M elements array with an aperture of the array N reach the minimum mutual coupling leakage when (5) holds. The array designed in (5) is a sparse uniform line array (SULA) when N>M−1λ.
However, SULA will cause spatial aliasing during DOA estimation [29] because the spacing of the adjacent sensors are larger than half-wavelength. Next, we discuss how to fine-tune the position of SULA's array elements to solve the angular ambiguity problem and maintain the advantage of low mutual coupling leakage.


### 3.2. B Mutual-Coupling-Optimized Array (MCOA)

Suppose that the first sensor is located at the origin without loss of generality, then the position of the array can be expressed as follows:(12)ℤ=0,z2,⋯,zMd.In order to facilitate the subsequent discussion,zi needs to be adjusted to integer by choosing an appropriate d.Theorem 1.
(see [29]). The array manifold aθ=1,ej2πz2sinθ/λ,⋯,ej2πzMsinθ/λT is invertible if and only if the sensor locations zi (assumed integers) are coprime.(13)gcdz2,z3,⋯,zM=1.
According to heorem 1, we design the array structure as follows:(14)d=x1,x2,…,xM−1=a−1,a,…,a︸M−2/2,a+1,a,…,a,N=aM−1,a,…,a+1,…,a+1︸b,…,a,N=aM−1+b,b<M−1,where d1×M−1 represent the spacing of adjacent sensors in array and a,b are integers.
In fact, the different arrangement order of elements ind will also cause the structural change of Z, which leads to different mutual coupling leakage. Due to the large number of repetitions of elements in d, there will be many repeated combinations of corresponding Z. Obviously, when the larger distance between adjacent sensors in the center of the array, the mutual coupling between the middle sensors and the sensors on both sides can be effectively reduced so that the mutual coupling leakage of the whole array degrades significantly. This is the reason why we make d as (14). The relationship between Z and d is zi+1−zi=di, then the expression of the array position Z can be written as follows:(15)Z=0,a−1,2a−1,⋯,Nd,N=aM−1,0,a,⋯,ka+1,⋯,Nd,N=aM−1+b,b<N+1.
We will prove that the position in (15) satisfies Theorem 1 in the following part.Proof.
WhenN=aM−1, then(16)gcd2a−1,a−1=gcda−1,2a−1moda−1=1=1,where amodb represents the remainder of dividing a by b.
WhenN=aM−1+b,b<N+1, then.(17)gcda,ka+1=gcda,ka+1moda=gcda,1=1.
In summarygcdz2,z3,⋯,zM=1, that is, the designed array structure satisfies Theorem 1. Suppose that an antenna array with array aperture N=8λ and the number of sensors M=8 needs to be designed. According to (15), since N is not divisible by M−1, we can calculate that N=2M−1+2 with a=2,b=2 and write the expression of d and Z.(18)d=2,2,3,3,2,2,2,ℤ=0,2,4,7,10,12,14,16d.
Figure2 shows the example of the above MCOA. The array structure is very similar to that of the SULA. There are five spacing between the sensors are 2d and the remaining two are 3d. We place the 3d in the middle of the array to make the mutual coupling between the middle sensors and sensor on both sides decreasing, which effectively reduces the coupling leakage of the whole array. The array locations also satisfy Theorem 1.(19)gcd2,4,7,10,12,14,16=1.Figure 2
MCOA array with an array aperture of 8 wavelengths and a sensor count of 8.


## 3.1. A Minimum Mutual Coupling Leakage Array

Considering theM- element array with an array aperture of N, we propose an array which has no spatial aliasing in DOA estimation and has the smaller mutual coupling leakage than most sparse arrays.Lemma 1.
For anM element array with an array aperture of N, the mutual coupling leakage is minimized if and only if the array elements are equally spaced.(5)zi+1d−zid=zid−zi−1d=NM−1,i=2,3,⋯,M−1.Proof.
According to equation (3), the mutual coupling matrix can be expressed as follows:(6)C=c0cz2−z1⋯czM−z1cz2−z1c0⋯czM−z2⋮⋮⋱⋮czM−z1czM−z2⋯c0.
Then, the expression of the coupling leakage can be calculated as follows:(7)γ=C−diagCFCF=2∑i=1M−1∑j=i+1Mczj−zi2Mc02+2∑i=1M−1∑j=i+1Mczj−zi2,where(8)∑i=1M−1∑j=i+1Mczj−zi2=c12∑i=1M−1∑j=i+1M1zjd−zid2=c12S,where S=∑i=1M−1∑j=i+1Mzjd−zid−2. The value of zi,i=1,2,⋯,M corresponding to the minimum value of S is the position of each array sensor when the mutual coupling leakage is minimum.
Denote theM−1 array spacings as x1,x2,⋯,xM−1 respectively. Then,(9)S=∑i=1M−1∑j=i+1Mzjd−zid−2=∑i=1M−1zi+1d−zid−2+∑i=1M−2zi+2d−zid−2+⋯+∑i=11zMd−zid−2=∑i=1M−1xi−2+∑i=1M−2xi+xi+1−2+⋯+∑i=11xi+xi+1+⋯+xi+M−2−2.
Calculate the minimum value ofS by Lagrange multiplier method.(10)gx1,⋯,xM−1,μ=S+μx1+x2+⋯+xM−1−N.
Take the partial derivative of each variable in the function and set the result equal to zero.(11)∂gx1,⋯,xM−1,μ∂x1=0⋮∂gx1,⋯,xM−1,μ∂xM−1=0∂gx1,⋯,xM−1,μ∂μ=0.
There is an extreme point in (10) when x1=x2=⋯=xM−1=N/M−1 and the cost function S has a minimum value at this time. Therefore, the value of zi,i=1,2,⋯,M corresponding to the minimum value of S is the position of each array element in the array when the mutual coupling leakage is minimum. From the above-given discussion, it can be seen that the M elements array with an aperture of the array N reach the minimum mutual coupling leakage when (5) holds. The array designed in (5) is a sparse uniform line array (SULA) when N>M−1λ.
However, SULA will cause spatial aliasing during DOA estimation [29] because the spacing of the adjacent sensors are larger than half-wavelength. Next, we discuss how to fine-tune the position of SULA's array elements to solve the angular ambiguity problem and maintain the advantage of low mutual coupling leakage.


## 3.2. B Mutual-Coupling-Optimized Array (MCOA)

Suppose that the first sensor is located at the origin without loss of generality, then the position of the array can be expressed as follows:(12)ℤ=0,z2,⋯,zMd.In order to facilitate the subsequent discussion,zi needs to be adjusted to integer by choosing an appropriate d.Theorem 1.
(see [29]). The array manifold aθ=1,ej2πz2sinθ/λ,⋯,ej2πzMsinθ/λT is invertible if and only if the sensor locations zi (assumed integers) are coprime.(13)gcdz2,z3,⋯,zM=1.
According to heorem 1, we design the array structure as follows:(14)d=x1,x2,…,xM−1=a−1,a,…,a︸M−2/2,a+1,a,…,a,N=aM−1,a,…,a+1,…,a+1︸b,…,a,N=aM−1+b,b<M−1,where d1×M−1 represent the spacing of adjacent sensors in array and a,b are integers.
In fact, the different arrangement order of elements ind will also cause the structural change of Z, which leads to different mutual coupling leakage. Due to the large number of repetitions of elements in d, there will be many repeated combinations of corresponding Z. Obviously, when the larger distance between adjacent sensors in the center of the array, the mutual coupling between the middle sensors and the sensors on both sides can be effectively reduced so that the mutual coupling leakage of the whole array degrades significantly. This is the reason why we make d as (14). The relationship between Z and d is zi+1−zi=di, then the expression of the array position Z can be written as follows:(15)Z=0,a−1,2a−1,⋯,Nd,N=aM−1,0,a,⋯,ka+1,⋯,Nd,N=aM−1+b,b<N+1.
We will prove that the position in (15) satisfies Theorem 1 in the following part.Proof.
WhenN=aM−1, then(16)gcd2a−1,a−1=gcda−1,2a−1moda−1=1=1,where amodb represents the remainder of dividing a by b.
WhenN=aM−1+b,b<N+1, then.(17)gcda,ka+1=gcda,ka+1moda=gcda,1=1.
In summarygcdz2,z3,⋯,zM=1, that is, the designed array structure satisfies Theorem 1. Suppose that an antenna array with array aperture N=8λ and the number of sensors M=8 needs to be designed. According to (15), since N is not divisible by M−1, we can calculate that N=2M−1+2 with a=2,b=2 and write the expression of d and Z.(18)d=2,2,3,3,2,2,2,ℤ=0,2,4,7,10,12,14,16d.
Figure2 shows the example of the above MCOA. The array structure is very similar to that of the SULA. There are five spacing between the sensors are 2d and the remaining two are 3d. We place the 3d in the middle of the array to make the mutual coupling between the middle sensors and sensor on both sides decreasing, which effectively reduces the coupling leakage of the whole array. The array locations also satisfy Theorem 1.(19)gcd2,4,7,10,12,14,16=1.Figure 2
MCOA array with an array aperture of 8 wavelengths and a sensor count of 8.


## 4. Sparse Bayesian-Based Compressed Sensing Method

In the sparse signal representation framework [23, 24], the direction matrix A in (2) should be replaced by a transfer matrix Ag, thus the signal model in (2) can be rewritten as follows:(20)X=CAgSg+N,where Sg=s1,⋯,sJ∈ℂG×J represent the complex signal amplitudes containing G DOAs and J snapshots. The transfer matrix Ag=a1,⋯,aG∈ℂM×G consists of all hypothetical DOAs. The likelihood function of the received signal X can be represented as follows:(21)pX|Sg;σ2=exp−1/σ2X−CAgSgF2πσ2NL.The SBL algorithm treatss as a zero mean complex Gaussian random vector with unknown diagonal covariance Γ=diagγ1,⋯,γM=diagγ. The prior model is given by the following equation:(22)pSg=∏j=1Jpsj=∏j=1JAC0,Γ.For Gaussian prior and likelihood, the evidencepX is Gaussian and represented as follows:(23)pX=∫pSgpX|SgdSg=∏j=1JACxj;0,Σx,where Σx=σ2I+CAgΓAgHCH and I stands for the identity matrix of order M×M. The SBL algorithm is to estimate the diagonal entries of Γ by maximizing the evidence(24)γ^1,⋯,γ^M=argmaxγ−∑j=1JxjHΣy−1xj−LlogΣx,the derivative of (24) is(25)∂−∑j=1JxjH∑y−1xj−LlogΣx∂γm=trXH∑x−1amamH∑x−1X−LamHΣx−1am=XH∑x−1am22−LamH∑x−1am=γmoldγmnew2XH∑x−1am22−LamH∑x−1am.The factorγmold/γmnew2 is introduced to obtain an iterative equation in γm. Equate the derivatives to zero and we can get the fixed-point update rule [1, 2, 4].(26)γmnew=γmold1LXH∑x−1am22amH∑x−1am=γmoldTrSx∑x−1amamHΣx−1amH∑x−1am.where Sx=1/JXXH is the sample covariance matrix.The main steps of the SBL algorithm are summarized as follows:Step 1: Set the parameters asε=10−3 and get the input data X,Ag,σ2,kmax;Step 2: Initialization the parametersγmold=1,∀m;Step 3: CalculateΣx=σ2I+CAΓoldAHCH;Step 4: Updateγmnew by using equation (27);Step 5: Setγold=γnew,Γold=diagγold,k=k+1;Step 6: Ifγnew−γold1/γold1>ε and k<kmax, go back to step 3;Step 7: The K largest peaks inγ are the required DOA values.


## 5. Performance Analysis

### 5.1. Comparison of Mutual Coupling Leakage between Different Arrays with the Same Aperture

We select some sparse arrays for comparison. In order to make the array aperture of all arrays consistent, we compress the element spacing of CA and ECA. The result is shown in Table1, it can be seen that the mutual coupling leakage of the proposed array is the smallest except SULA, and its mutual coupling leakage is very close to SULA.Table 1
Comparison of mutual coupling leakage of different arrays.
Array aperturedZMutual coupling leakageProposed8λ0.5λ0,2,4,7,10,12,14,160.0826CA8λ0.5λ0 4 5 8 10 12 15 160.1076ECA8λ0.4444λ0 2 4 5 6 8 13 180.1140NA8λ0.4211λ0 1 2 3 7 11 15 190.1194SULA8λ1.1429λ0 1 2 3 4 5 6 70.0806ULA3.5λ0.5λ0 1 2 3 4 5 6 70.1819


### 5.2. Crarmer–Rao Bound (CRB)

According to the knowledge of the literature [32], the Crarmer–Rao Bound (CRB) matrix can be represented as follows:(27)CRB=σn22JReDH∏A⊥DP^T−1,where Re· stands for the operation of taking the real part, A represents the manifold matrix of the array, ΠA⊥=I−AAHA−1AH is the orthogonal projection of A, and I stands for the identity matrix of order M×M, P^=1/J∑t=1JstsHt, σn2 denotes the average power of signal source, D can be written as follows:(28)D=∂aθ1∂θ1,∂aθ2∂θ2,⋯,∂aθK∂θK,where aθK denotes steering vector.


### 5.3. Computational Complexity

In this section, we provide the complexity of the SBL method compared with OMP and MUSIC. Assuming that the number of array elements isM, there are G grid points and J snapshots, the maximum number of iterations is kmax, then the computational complexity of main operations are as follows: (a) calculate the covariance matrix: OM2J; (b) update the Σx in Step 3: O2M2G+G2M+M3; (C) update γmnew in Step 4: OG2M3+3M2+M. The computational complexity of SBL is OM2J+kmax2GM3+5GM2+G2M+GM+M3.


## 5.1. Comparison of Mutual Coupling Leakage between Different Arrays with the Same Aperture

We select some sparse arrays for comparison. In order to make the array aperture of all arrays consistent, we compress the element spacing of CA and ECA. The result is shown in Table1, it can be seen that the mutual coupling leakage of the proposed array is the smallest except SULA, and its mutual coupling leakage is very close to SULA.Table 1
Comparison of mutual coupling leakage of different arrays.
Array aperturedZMutual coupling leakageProposed8λ0.5λ0,2,4,7,10,12,14,160.0826CA8λ0.5λ0 4 5 8 10 12 15 160.1076ECA8λ0.4444λ0 2 4 5 6 8 13 180.1140NA8λ0.4211λ0 1 2 3 7 11 15 190.1194SULA8λ1.1429λ0 1 2 3 4 5 6 70.0806ULA3.5λ0.5λ0 1 2 3 4 5 6 70.1819


## 5.2. Crarmer–Rao Bound (CRB)

According to the knowledge of the literature [32], the Crarmer–Rao Bound (CRB) matrix can be represented as follows:(27)CRB=σn22JReDH∏A⊥DP^T−1,where Re· stands for the operation of taking the real part, A represents the manifold matrix of the array, ΠA⊥=I−AAHA−1AH is the orthogonal projection of A, and I stands for the identity matrix of order M×M, P^=1/J∑t=1JstsHt, σn2 denotes the average power of signal source, D can be written as follows:(28)D=∂aθ1∂θ1,∂aθ2∂θ2,⋯,∂aθK∂θK,where aθK denotes steering vector.


## 5.3. Computational Complexity

In this section, we provide the complexity of the SBL method compared with OMP and MUSIC. Assuming that the number of array elements isM, there are G grid points and J snapshots, the maximum number of iterations is kmax, then the computational complexity of main operations are as follows: (a) calculate the covariance matrix: OM2J; (b) update the Σx in Step 3: O2M2G+G2M+M3; (C) update γmnew in Step 4: OG2M3+3M2+M. The computational complexity of SBL is OM2J+kmax2GM3+5GM2+G2M+GM+M3.


## 6. Simulation Results

In this part, we provide numerical simulations of the performance of the proposed MCOA as well as a comparison with the other sparse arrays and the CRB. The array aperture is limited to 8 wavelengths and the number of array sensors is 8. Two coherent signals with equal power impinge on the array with directionsθ1=10∘, θ2=40∘, and the correlation coefficient is set to α1,α2=1,ejπ/4. Define the Root Mean Square Error (RMSE) of the DOA estimates as follows:(29)RMSE=1K1Q∑k=1K∑q=1Qθ^q,k−θk2,where Q and K are the number of Monte Carlo trials and the total number of coherent signals, respectively. θ^q,k means the qth estimate of the real angle θk . Unless other stated, we assume that the mutual coupling constant is c0=0.12 and the maximum spacing of sensor pairs is B=100. For each Figure, 1000 Monte Carlo simulations were run to estimate the Root Mean Square Error (RMSE).First, we compare the CRB of different arrays in Figure3. The result shows that the CRB of the coprime array is very close to the CRB of the proposed array, but the CRB of the proposed array is the smallest among all arrays, which indicates that its performance is optimal.Figure 3
CRB comparison of different arrays.Figure4 depicts the spatial spectrum of the SBL method with the proposed array. The signal-to-noise ratio (SNR) is 10 dB and the snapshot is J=200. There are many spectral peaks in the estimation result, but the peak of the incident signal is the highest, which is 40 dB higher than other spectral peaks. It shows that the estimation result of this method is accurate enough to be used for the coherent signal estimation when there is mutual coupling between array sensors.Figure 4
The spatial spectrum for DOA estimation.The proposed array is also compared with other arrays with different SNRs in Figure5. The SNR varies from −5 dB to 20 dB and the number of snapshots is J=200. All five kinds of arrays can accurately estimate the incident angle of the relevant signals. Due to the influence of mutual coupling, the angle estimation of the CA and ULA decreases greatly, and their RMSE value is larger than the proposed array. Because the influence of mutual coupling leakage is less than that of other arrays, the proposed array can better estimate the DOA of coherent signals.Figure 5
The DOA estimation performance with different SNRs.Figure6 shows the performance of different arrays with snapshots changing, the coherent signal is consistent with the previous simulation and the SNR is 5 dB. The snapshot varies from 10 to 600. When the snapshot is less than 100, all five arrays cannot work well and the performance is not good enough. However, the RMSE of all the arrays reduces with the snapshot increasing. When the snapshot is larger than 100, the curve of RMSE of the CA and ULA is almost a straight line, because the impact of the mutual coupling leakage at this time is greater than the performance improvement brought about by the increase of snapshots. On the other hand, the performance of the proposed array and NA is better than the other three arrays, and the proposed array has the best performance because the mutual coupling leakage of the proposed array is the smallest among these arrays.Figure 6
The DOA estimation performance with different snapshots.Finally, we compared the SBL algorithm with two other compressed sensing algorithms including the IAA [20] and the OMP [21]. In Figure 7, the SNR changes from −5 dB to 20 dB and the snapshot is 200. The other simulation conditions were the same as before. At low SNR, the SBL has the same performance as IAA and OMP, but when the SNR increases, the performance of OMP hardly improves, and the performance of IAA is better than OMP. The curve of SBL is smoother than the other two algorithms, which means that in this case, its performance is more stable than other algorithms.Figure 7
Comparison of RMSE of different algorithms with SNR.


## 7. Conclusions

In this paper, the mutual coupling optimization array with a given number of sensors under the condition of the finite aperture is studied. Compared with the sparse arrays including CA, NA, and ECA, the mutual coupling leakage of the proposed array is smaller. When there is mutual coupling between array sensors, we apply the SBL to DOA estimation of coherent signals and compare its performance with other compressed sensing methods including IAA and OMP. Finally, various simulations are carried out to prove the superior performance of the proposed array for estimating coherent signals in the condition of mutual coupling. In fact, this paper mainly focuses on the coherent signal estimation problem of the sparse array in the 1D-DOA case. By using the previous related research, this result can be extended to an L-shaped array to realize DOA estimation of coherent signal in the 2D-DOA case. Reference [31].



---

*Source: 8800833-2022-10-05.xml*