Augmented Nested Arrays With Enhanced DOF and Reduced Mutual Coupling
摘要——最近,非均匀线性阵列(例如,互质/嵌套阵列)因其能够生成虚拟差分协同阵列而在阵列信号处理领域引起了研究人员的极大关注。在阵列设计中,一个关键问题是如何放置传感器以获得最佳性能,旨在同时实现最大自由度容量和最小互耦比。本文提出了一种增广嵌套阵列的概念,通过将嵌套阵列的密集子阵列分成几个部分,这些部分可以重新排列在嵌套阵列稀疏子阵列的两侧。具体来说,对于任意给定的阵元数量,我们推导出了物理传感器位置和虚拟传感器位置的四个闭式表达式。与具有相同阵元数量的(超级)嵌套阵列相比,新形成的增广嵌套阵列拥有更高的自由度容量和更小的互耦。最后,数值仿真结果验证了所提出阵列的有效性。
索引术语——增广嵌套阵列,自由度,互耦,DOA估计,稀疏线性阵列。
文章目录
- I. 引言
- II. REVIEW OF DIFFERENCE CO-ARRAY SIGNAL MODEL
- A. Difference Co-Array Signal Model
- B. Mutual Coupling
- D. Properties of Nested Array
- III. AUGMENTED NESTED ARRAYS
- A.Two Level Augmented Nested Arrays
- B. ANAI-1
- C. ANAI-2
- IV. NUMERICAL EXPERIMENTS
- **A. 自由度**
- **B. 互耦**
- **C. 使用ESPRIT算法进行DOA估计**
I. 引言
阵列信号处理在许多应用中扮演着至关重要的角色,例如雷达[1],通信[2]和导航[3]。传统上,研究人员主要关注均匀线性阵列(ULAs),其中相邻传感器之间的阵元间距小于λ/2\lambda/2λ/2以避免空间混叠。众所周知,阵列孔径大小在基于阵列的系统中起着基础性作用。对于此类ULA,孔径的增加通常需要增加天线阵元的数量,从而导致更高的硬件成本和计算复杂性。为了获得最大的空间分辨率/自由度(DOF)和最小的互耦,非均匀线性阵列(NLA)比ULA吸引了更多的关注,这可以追溯到最小冗余阵列(MRAs)[4]或最小最小孔洞阵列(MHAs)[5]。许多重要的优点在文献[6]–[8]中通过直接或间接形成增广协方差矩阵得到了验证。通过利用增广协方差矩阵,这些非均匀阵列可以成功地避免空间混叠,并能够探测到比传感器数量更多的信源。
最近,嵌套阵列[9]和互质阵列[10]的引入和分析重新激发了人们对此类非均匀线性阵列(NLA)几何结构的兴趣。互质/嵌套阵列由两个具有不同阵元间距的均匀线性子阵列组成,可以生成一组均匀分布的虚拟差分协同阵列。因此,许多阵列信号处理算法,例如波束形成[9], [10]、波达方向(DOA)估计[9]–[12]和稀疏DOA恢复[13]–[15],都可以基于虚拟阵列(而非物理阵列)来执行。由于差分协同阵列可以将阵列孔径和自由度(DOF)从O(N)O(N)O(N)增加到O(N2)O(N^2)O(N2),因此不仅可以获得更好的空间分辨率和更高的自由度容量,而且还能成功避免空间混叠。此外,通过结合多频率操作[16], [17]或联合时频变换[18],这类NLA的最大自由度可以得到适度增加。遵循这些协同阵列的概念,一些更有效的NLA几何结构被进一步开发出来,例如超级嵌套阵列(SNA)[19]–[21]、嵌套最小冗余阵列(NMRA)[22]以及一些优化和广义的互质阵列[23]–[25]。
这些NLA几何结构仍然存在一些局限性。从没有闭式表达式的MRA/MHA/NMRA开始,只有一些结果是通过穷举搜索[26]列出的。(广义)互质阵列的协同阵列中存在孔洞(见图1(a))。嵌套阵列包含一个密集的子均匀线性阵列(sub-ULA),导致了高的互耦。相对而言,SNA获得了最佳性能。它保留了嵌套阵列的所有优良特性,同时通过重新分布密集子ULA的阵元来降低互耦。
理论上,对于任何具有无孔洞差分协同阵列的 MMM 元NLA,其冗余度 RRR 和单边均匀自由度 LuL_uLu 将满足以下条件:R+Lu=M(M−1)/2R+L_u = M(M-1)/2R+Lu=M(M−1)/2,其中R=∑ℓ=1Lu(w(ℓ)−1)R=\sum_{\ell=1}^{L_u}(w(\ell)-1)R=∑ℓ=1Lu(w(ℓ)−1),w(ℓ)w(\ell)w(ℓ)表示间距为 ℓλ/2\ell\lambda/2ℓλ/2 的传感器对的数量。因此,互耦比主要由重复传感器对的数量(即RRR)及其分布决定。根据[20], [21],嵌套阵列和SNA拥有相同的自由度 LuL_uLu 和相同的冗余度 RRR。对于SNA而言,其较低的互耦比很大程度上是由于对RRR个重复传感器对的分布进行了优化(即减少了小间距传感器对的数量)。
在本文中,我们关注一种双赢策略:单边均匀自由度LuL_uLu的增加可以带来冗余度RRR的减少,这可能会降低互耦。在最佳情况下,如果 Lu=M(M−1)/2L_u = M(M-1)/2Lu=M(M−1)/2,冗余度 RRR 将变为0。基于此策略,我们提出了一种增广嵌套阵列(augmented nested array,ANA)的概念,通过将嵌套阵列的密集ULA(均匀线性阵列)拆分成若干个左/右子阵列,这些子阵列可以重新排列在嵌套阵列的稀疏ULA的两侧。此外,我们建立了一些充分条件来保证相应差分协同阵列的无孔洞特性。结果,所形成的ANA具有以下优点:
- i) ANA具有闭式解的传感器位置和无孔洞的协同阵列(hole-free co-arrays)。
- ii) 在相同阵元数下,ANA比(generalized,广义)互质阵列和(super,超级)嵌套阵列具有更高的 uniform DOF capacity容量。
- iii) ANA比嵌套阵列和前几阶的超级嵌套阵列具有更低的互耦。
II. REVIEW OF DIFFERENCE CO-ARRAY SIGNAL MODEL
我们首先介绍一些集合运算,以便在下文中描述物理阵列的几何形状和相应的虚拟阵列配置。
定义1:对于两个给定的整数集合 A\mathbb AA 和 B\mathbb BB,可以定义以下运算:
- 交集: A∩B={a∣∀a∈Aand a∈B}A \cap B = \{a|∀a \in A \text{ and } a \in B\}A∩B={a∣∀a∈A and a∈B}
- 并集: A∪B={a∣∀a∈Aor a∈B}A \cup B = \{a|∀a \in A \text{ or } a \in B\}A∪B={a∣∀a∈A or a∈B}
- 和集: A+B={a+b∣∀a∈A,∀b∈B}A + B = \{a + b|∀a \in A, ∀b \in B\}A+B={a+b∣∀a∈A,∀b∈B}
- 差集: A−B={a−b∣∀a∈A,∀b∈B}A - B = \{a - b|∀a \in A, ∀b \in B\}A−B={a−b∣∀a∈A,∀b∈B}
- 平移集: λ+A={λ+a∣∀a∈A,λ∈Z}λ + A = \{λ + a|∀a \in A, λ \in Z\}λ+A={λ+a∣∀a∈A,λ∈Z}
- 伸缩集: λA={λa∣∀a∈A,λ∈Z}λA = \{λa|∀a \in A, λ \in Z\}λA={λa∣∀a∈A,λ∈Z}
A. Difference Co-Array Signal Model
考虑一个非均匀线性阵列(NLA),其在位置A={a1,a2,...,aM}1\mathbb{A} = \{a_1, a_2, ..., a_M\}^1A={a1,a2,...,aM}1处有MMM个传感器,并且有DDD个窄带远场不相关信源从方向θi\theta_iθi入射到阵列上,其功率为δi2\delta_i^2δi2,其中i=1,2,3,...,Di = 1, 2, 3, ..., Di=1,2,3,...,D。因此,第nnn个接收快照可以表示为
x(n)=As(n)+w(n),(1)x(n) = As(n) + w(n), \tag{1} x(n)=As(n)+w(n),(1)
其中s(n)s(n)s(n)和w(n)w(n)w(n)分别表示输入信源向量和噪声向量;AAA表示一个M×DM \times DM×D维的阵列流形矩阵,其第(l,i)(l,i)(l,i)个元素Al,iA_{l,i}Al,i为e−jπalsin(θi)e^{-j\pi a_l \sin(\theta_i)}e−jπalsin(θi),其中l=1,2,...,Ml=1,2,...,Ml=1,2,...,M;i=1,2,...,Di=1,2,...,Di=1,2,...,D。然后通过采集NNN个快照,可以估计出协方差矩阵为
Rxx=E{xxH}=ARssAH+σn2IM≈1N∑n=1Nx(n)x(n)H,(2)R_{xx} = E\{xx^H\} = AR_{ss}A^H + \sigma_n^2I_M \approx \frac{1}{N}\sum_{n=1}^{N}x(n)x(n)^H, \tag{2} Rxx=E{xxH}=ARssAH+σn2IM≈N1n=1∑Nx(n)x(n)H,(2)
其中Rss=E{ssH}=diag{p}R_{ss} = E\{ss^H\} = diag\{p\}Rss=E{ssH}=diag{p}, p=[δ12,δ22,...,δD2]Tp = [\delta_1^2, \delta_2^2, ..., \delta_D^2]^Tp=[δ12,δ22,...,δD2]T,σn2\sigma_n^2σn2表示噪声功率。根据[9], [10], [20],RxxR_{xx}Rxx可以被向量化为一个长向量
z=vec(Rxx)=(A∗⊙A)p+σn2vec(IM)∈CM2×1,(3)z = vec(R_{xx}) = (A^* \odot A)p + \sigma_n^2vec(I_M) \in \mathbb{C}^{M^2 \times 1}, \tag{3} z=vec(Rxx)=(A∗⊙A)p+σn2vec(IM)∈CM2×1,(3)
¹本文中所有的阵元位置都以半个最小波长为单位。
其中 iM=vec{IM}i_M = vec\{I_M\}iM=vec{IM}。在上述向量化模型中,项 (A∗⊙A)\left(\mathbf{A}^{*} \odot \mathbf{A}\right)(A∗⊙A) 的行为类似于一个虚拟阵列的流形矩阵(the manifold matrix of a virtual array),该虚拟阵列的传感器位置由物理阵元位置的差分值给出(即A−A\mathbb{A} - \mathbb{A}A−A)。因此,这个虚拟阵列被称为差分协同阵列(difference co-array,DCA)。根据[9], [15],可以从z中选择并平滑所需的行,以构建来自虚拟ULA(均匀线性阵列)的单快照数据,其传感器位于位置 U=[−Lu,Lu]\mathbb{U} = [-L_u, L_u]U=[−Lu,Lu]。因此,单快照模型可以表示为
z1=Jz=A1p+σn2eLu∈C(2Lu+1)×1(4)\mathbf{z}_{1}=\mathbf{J} \mathbf{z}=\mathbf{A}_{1} \mathbf{p}+\sigma_{n}^{2} \mathbf{e}_{L_{u}} \in \mathbb{C}^{\left(2 L_{u}+1\right) \times 1} \tag{4} z1=Jz=A1p+σn2eLu∈C(2Lu+1)×1(4)
其中eLu=[01×Lu101×Lu]Te_{L_u} = [0_{1 \times L_u} \ 1 \ 0_{1 \times L_u}]^TeLu=[01×Lu 1 01×Lu]T;A1A_1A1表示虚拟ULA的一个(2Lu+1)×D(2L_u+1) \times D(2Lu+1)×D维流形矩阵;J\bf JJ 表示一个(2Lu+1)×M2(2L_u+1) \times M^2(2Lu+1)×M2 的选择矩阵。在其每一行中只有一个非零元素,其位置对应于从 M2M^2M2 个虚拟传感器中选出的虚拟阵元的索引。许多阵列信号处理算法可以基于虚拟信号(4)而不是(1)来执行。由于自由度(DOF)的增加,虚拟阵列信号模型可以分辨多达 O(M2)O(M^2)O(M2) 个不相关信源。
B. Mutual Coupling
公式(1)中的阵列信号模型是在不考虑物理传感器之间互耦的情况下建立的。在实践中,应考虑小间距传感器对之间的耦合效应。因此,通过引入一个耦合矩阵C,公式(1)可以修改为
x(n)=CAs(n)+w(n).(5)\mathbf{x}(n)=\mathbf{C A s}(n)+\mathbf{w}(n) . \tag{5} x(n)=CAs(n)+w(n).(5)
矩阵 C\mathbf{C}C 的元素取决于多种因素,如工作频率、天线类型、相控阵的工作模式、传感器之间的距离、阻抗等。C\mathbf{C}C 的闭式表达式很复杂。在本文中,由于我们只关注线性阵列的几何结构,根据[20], [27]–[29],C可以由一个 a B-banded mode 来近似
Ci,j={c∣ai−aj∣∣ai−aj∣≤B0∣ai−aj∣>B,(6)\mathbf{C}_{i, j}=\left\{\begin{array}{cl} c_{\left|a_{i}-a_{j}\right|} & \left|a_{i}-a_{j}\right| \leq B \\ 0 & \left|a_{i}-a_{j}\right|>B \end{array},\right. \tag{6} Ci,j={c∣ai−aj∣0∣ai−aj∣≤B∣ai−aj∣>B,(6)
其中 Ci,j\mathbf{C}_{i, j}Ci,j 的大小与传感器间距成反比,并满足1=c0>∣c1∣>∣c2∣>⋯>∣cB∣>∣cB+1∣=01 = c_0 > |c_1| > |c_2| > \cdots > |c_B| > |c_{B+1}| = 01=c0>∣c1∣>∣c2∣>⋯>∣cB∣>∣cB+1∣=0。例如,文献[20]中使用的耦合系数如下:c0=1,c1=0.3ejπ/3c_0 = 1, c_1 = 0.3e^{j\pi/3}c0=1,c1=0.3ejπ/3以及cl=c1e−j(l−1)π/8/l,2≤l≤Bc_l = c_1e^{-j(l-1)\pi/8}/l, 2 \le l \le Bcl=c1e−j(l−1)π/8/l,2≤l≤B。此外,对于一个给定的NLA,总互耦可以通过耦合泄漏[20] L(M)L(\mathbf M)L(M) 来评估,即
L(M)=∥C−C~∥F∥C~∥F,(7)L(M) = \frac{\|\mathbf{C}-\widetilde{\mathbf{C}}\|_{F}}{\|\tilde{C}\|_F}, \tag{7} L(M)=∥C~∥F∥C−C∥F,(7)
其中 C~i,j={0i≠jCi,ji=j\widetilde{\mathbf{C}}_{i, j}=\left\{\begin{array}{cc} 0 & i \neq j \\ \mathbf{C}_{i, j} & i=j \end{array}\right.Ci,j={0Ci,ji=ji=j;∥C−C~∥F\|\mathbf{C}-\widetilde{\mathbf{C}}\|_{F}∥C−C∥F是所有非对角线分量的能量,它表征了互耦的量。因此 L(M)L(\mathbf M)L(M) 是这两个分量之间的能量比,在本文中也称为互耦比。从概念上讲,L(M)L(\mathbf M)L(M) 越小,互耦就越小。
C. 具有无孔洞协同阵列的NLA的性质
在介绍差分协同阵列模型之后,我们从建立一个用于评估NLA性能的定义开始,重点关注物理NLA的构建以及相应虚拟阵列的性质。
定义2:对于一个阵列A\mathbb{A}A,令D\mathbb{D}D表示其对应的DCA;U\mathbb UU 表示具有最大孔径的中心ULA;M(ℓ)\mathbb M(\ell)M(ℓ)表示在延迟 ℓ\ellℓ 处对协同阵列支撑有贡献的所有对。
A={a1,a2,…,aM}D={ai−aj∣i,j∈[1,M]}U={−Lu,Lu}⊆DM(ℓ)={(ai,aj)∣ai−aj=ℓ;i,j∈[1,M]}(8)\begin{aligned} \mathbb{A} &= \{a_1, a_2, \dots, a_M\} \\ \mathbb{D} &= \{a_i - a_j | i, j \in [1, M]\} \\ \mathbb{U} &= \{-L_u, L_u\} \subseteq \mathbb{D} \\ \mathbb M(\ell) &= \{(a_i, a_j) | a_i - a_j = \ell; i, j \in [1, M]\} \end{aligned} \tag{8} ADUM(ℓ)={a1,a2,…,aM}={ai−aj∣i,j∈[1,M]}={−Lu,Lu}⊆D={(ai,aj)∣ai−aj=ℓ;i,j∈[1,M]}(8)
那么,D\mathbb{D}D, U\mathbb{U}U 和 M(ℓ)\mathbb{M}(\ell)M(ℓ) 的基数分别被称为自由度DOF、均匀自由度(uDOF)和权重函数(表示为d(M),u(M)d(M), u(M)d(M),u(M) 和 w(ℓ)w(\ell)w(ℓ))。此外,LuL_uLu被称为单边均匀自由度(one-side uniform DOF)。
为了空间采样的目的,完美的NLA应满足以下标准:
- 对于任意给定的阵元数M,具有闭式解的传感器位置,即 A=F(M)={a1,a2,…,aM}\mathbb{A} = F(M) = \{a_1, a_2, \dots, a_M\}A=F(M)={a1,a2,…,aM}。
- 最大化均匀自由度uDOF u(M)u(M)u(M),其值至少不小于嵌套阵列的uDOF。
- 最小化小间距的权重函数w(ℓ)w(\ell)w(ℓ),即w(1),w(2)w(1), w(2)w(1),w(2) 和 w(3)w(3)w(3)。
在提出这三个标准之后,我们首先将分析此类具有无孔洞协同阵列的NLA的 LuL_uLu 和 RRR 的边界,这可以为设计最优NLA提供理论基础。
定理1:考虑一个具有无孔洞协同阵列的 MMM 元NLA。假设其单边均匀自由度和冗余度分别表示为LuL_uLu和RRR,其中 R=∑ℓ=1Lu(w(ℓ)−1)R = \sum_{\ell=1}^{L_u}(w(\ell)-1)R=∑ℓ=1Lu(w(ℓ)−1) 且 Lu+R=M(M−1)/2L_u+R = M(M-1)/2Lu+R=M(M−1)/2。那么,LuL_uLu和RRR的理论边界可以表示为
Lu≤M2/(2+4/3π)R≥M2/(3π+2)(9)\begin{aligned} L_u &\le M^2/(2+4/3\pi) \\ R &\ge M^2/(3\pi+2) \end{aligned} \tag{9} LuR≤M2/(2+4/3π)≥M2/(3π+2)(9)
定理1的证明 在附录A-A中给出。由于嵌套阵列是一种著名的具有无孔洞协同阵列的NLA,我们将从数论的角度对其性质进行深入分析。所提出的ANA几何结构的灵感便来源于此。
D. Properties of Nested Array
嵌套阵列(The nested array)由一个由 N1+1N_1+1N1+1 个传感器的密集子ULA(均匀线性阵列)和一个有 N2N_2N2 个传感器的稀疏子ULA组成,其中密集子ULA可以生成从1到 N1N_1N1 的所有连续延迟,而稀疏ULA可以生成 (N1+1)(N_1+1)(N1+1) 倍的延迟,如图2所示。这种几何结构可以用一个有趣的数论来解释:任何非负整数vvv都可以被 (N1+1)(N_1+1)(N1+1) 整除,并表示为:v=v1(N1+1)+v2v = v_1(N_1+1) + v_2v=v1(N1+1)+v2,
²由于具有少量阵元的最优NLA已经通过计算机搜索找到。因此,在本文中,我们主要关注阵元数不小于10的情况。
其中 v1v_1v1 是非负商,v1∈[1,N2]v_1 \in [1, N_2]v1∈[1,N2];v2v_2v2 是非负余数,v2∈[0,N1]v_2 \in [0, N_1]v2∈[0,N1]。实际上,嵌套阵列的两个子阵列(一个密集,一个稀疏)分别对应于 v1v_1v1 和 v2v_2v2。密集子阵列由位于位置 As2=(N1+1)−ℓ\mathbb{A}_{s_2} = (N_1+1)-\ellAs2=(N1+1)−ℓ(其中 ℓ∈[0,N1]\ell \in [0, N_1]ℓ∈[0,N1])的 (N1+1)(N_1+1)(N1+1) 个传感器组成。稀疏子阵列由位于位置 As1=(N1+1)ℓ\mathbb{A}_{s_1} = (N_1+1)\ellAs1=(N1+1)ℓ(其中 ℓ∈[1,N2]\ell \in [1, N_2]ℓ∈[1,N2])的N2N_2N2个传感器组成。基于以上描述和解释,嵌套阵列可以表示为
ANested(M)=As1∪As2,(10)\mathbb{A}_{\text{Nested}}(M) = \mathbb{A}_{s_1} \cup \mathbb{A}_{s_2}, \tag{10} ANested(M)=As1∪As2,(10)
其中 MMM 表示物理传感器的数量,两个子阵列As1\mathbb{A}_{s1}As1和As2\mathbb{A}_{s2}As2由下式给出
{As1=(N1+1)ℓ,ℓ∈V1(N2)=[1,N2]As2=(N1+1)−ℓ,ℓ∈V2(N1)=[0,N1].(11)\begin{cases} \mathbb{A}_{s_1} = (N_1+1)\ell, & \ell \in V_1(N_2) = [1, N_2] \\ \mathbb{A}_{s_2} = (N_1+1)-\ell, & \ell \in V_2(N_1) = [0, N_1] \end{cases}. \tag{11} {As1=(N1+1)ℓ,As2=(N1+1)−ℓ,ℓ∈V1(N2)=[1,N2]ℓ∈V2(N1)=[0,N1].(11)
[23]中的CACIS-嵌套阵列与[9]中的嵌套阵列略有不同。其密集子ULA的特征是V2(N1)=[0,N1+1]V_2(N_1)=[0, N_1+1]V2(N1)=[0,N1+1]。实际上,任何在(10)中定义的阵列只有当[0,N1][0, N_1][0,N1]是 V2(N1)\mathbb{V}_2(N_1)V2(N1) 的子集时,才能保持嵌套阵列的基本特性(即无孔洞)。由于在(10)中定义的稀疏子阵列的阵元间距为(N1+1)(N_1+1)(N1+1),我们称此类阵列为基-(N1+1)(N_1+1)(N1+1)嵌套阵列。
性质1:在(10)和(11)中定义的嵌套阵列ANested(M)\mathbb{A}_{\text{Nested}}(M)ANested(M)的协同阵列是无孔洞的。此外,小间距传感器对的数量(the number of sensor pairs with small separation)由包含N1+1N_1+1N1+1个传感器的密集子ULA决定。
{Lu=(N1+1)N2−1w(ℓ)=ψN1+1(ℓ)=N1+1−ℓ,if 1≤ℓ≤N1.(12)\begin{cases} L_u = (N_1+1)N_2 - 1 \\ w(\ell) = \psi_{N_1+1}(\ell) = N_1+1-\ell, & \text{if } 1 \le \ell \le N_1 \end{cases}. \tag{12} {Lu=(N1+1)N2−1w(ℓ)=ψN1+1(ℓ)=N1+1−ℓ,if 1≤ℓ≤N1.(12)
性质1的证明在 附录A-B 中给出。由于密集子阵列的影响,通常嵌套阵列具有很高的互耦比。对于超级嵌套阵列,密集子阵列被分成许多部分,其中一些部分被插入到稀疏子阵列中,如图1(c)(c)(c)和(d)所示。
在本文中,提出了一种称为增广嵌套阵列(augmented nested array,ANA)的新几何结构。与超级嵌套阵列不同,密集子阵列被分裂成左/右子阵列。密集子阵列的一些传感器将被重新排列在稀疏子阵列的两侧(而不是插入到稀疏子阵列中),这可以同时带来自由度(DOF)容量的增加和互耦的降低。
III. AUGMENTED NESTED ARRAYS
A.Two Level Augmented Nested Arrays
如前所述,嵌套阵列可以由v=(N1+1)ℓ1+ℓ2v=(N_1+1)\ell_1+\ell_2v=(N1+1)ℓ1+ℓ2来表征,它也可以被重写为
v=((N1+1)N2+ℓ2)−(N1+1)(N2−ℓ1)=((N1+1)N2+ℓ2)−(N1+1)ℓ3,其中 ℓ3=N2−ℓ1(13)\begin{aligned} v &= ((N_1+1)N_2 + \ell_2) - (N_1+1)(N_2-\ell_1) \\ &= ((N_1+1)N_2+\ell_2) - (N_1+1)\ell_3, \quad \text{其中 } \ell_3=N_2-\ell_1 \end{aligned} \tag{13} v=((N1+1)N2+ℓ2)−(N1+1)(N2−ℓ1)=((N1+1)N2+ℓ2)−(N1+1)ℓ3,其中 ℓ3=N2−ℓ1(13)
这意味着位于左侧位置(N1+1)−ℓ2(N_1+1)-\ell_2(N1+1)−ℓ2的传感器可以被移动到右侧位置(N1+1)N2+ℓ2(N_1+1)N_2+\ell_2(N1+1)N2+ℓ2。受此概念启发,可以优化嵌套阵列的几何结构;非负余数V2V_2V2可以被分裂成两个子集:V21V_{21}V21和V22V_{22}V22。由V21V_{21}V21定义的传感器可以列在稀疏子阵列的左侧,而由V22V_{22}V22定义的传感器被移动到右侧。因此,这种两级ANA可以表示为
AANAI(M)=L1∪M1∪R1,(14)\mathbb{A}_{\text{ANAI}}(M) = \mathbb{L}_1 \cup \mathbb{M}_1 \cup \mathbb{R}_1, \tag{14} AANAI(M)=L1∪M1∪R1,(14)
其中L1\mathbb{L}_1L1, M1\mathbb{M}_1M1和R1\mathbb{R}_1R1分别表示左、中、右子阵列,可以表示为
{L1=−(N1+1)−ℓ,ℓ∈V21M1=(N1+1)ℓ,ℓ∈V1(N2)R1=(N1+1)N2+ℓ,ℓ∈V22.(15)\begin{cases} \mathbb{L}_1 = -(N_1+1)-\ell, & \ell \in V_{21} \\ \mathbb{M}_1 = (N_1+1)\ell, & \ell \in V_1(N_2) \\ \mathbb{R}_1 = (N_1+1)N_2+\ell, & \ell \in V_{22} \end{cases}. \tag{15} ⎩⎨⎧L1=−(N1+1)−ℓ,M1=(N1+1)ℓ,R1=(N1+1)N2+ℓ,ℓ∈V21ℓ∈V1(N2)ℓ∈V22.(15)
需要注意的是,默认情况下,0∈V210 \in V_{21}0∈V21且0∉V220 \notin V_{22}0∈/V22。
定理2 (几乎无孔洞):考虑在(14)和(15)中定义的两级ANA AANAI(M)\mathbb{A}_{\text{ANAI}}(M)AANAI(M)。如果V21∪V22=V2(N1)V_{21} \cup V_{22} = V_2(N_1)V21∪V22=V2(N1),那么该ANA的单边均匀自由度不小于((N1+1)N2−1)((N_1+1)N_2-1)((N1+1)N2−1)。此外,具有小间距的传感器对仅存在于左、右子NLA中,并且对于任何小的延迟ℓ\ellℓ,w(ℓ)w(\ell)w(ℓ)不大于(N1+1−ℓ)(N_1+1-\ell)(N1+1−ℓ)。
{Lu≥(N1+1)N2−1w(ℓ)=wL(ℓ)+wR(ℓ)≤N1+1−ℓ,if 1≤ℓ≤N1.(16)\begin{cases} L_u \ge (N_1+1)N_2-1 \\ w(\ell)=w_L(\ell)+w_R(\ell) \le N_1+1-\ell, & \text{if } 1 \le \ell \le N_1 \end{cases}. \tag{16} {Lu≥(N1+1)N2−1w(ℓ)=wL(ℓ)+wR(ℓ)≤N1+1−ℓ,if 1≤ℓ≤N1.(16)
定理2的证明在附录A-C中给出。根据该定理,与具有相同阵元数的父嵌套阵列相比,ANA可以同时增加(或至少保持)单边均匀自由度并减少互耦。其协同阵列可能在延迟±ℓ\pm \ell±ℓ处有孔洞,但幸运的是ℓ>(N1+1)N2−1\ell > (N_1+1)N_2-1ℓ>(N1+1)N2−1。因此,这种协同阵列是几乎无孔洞的。例如,考虑图3中的父嵌套阵列和ANA,嵌套阵列的协同阵列
B. ANAI-1
C. ANAI-2
IV. NUMERICAL EXPERIMENTS
A. 自由度
在第一个实验中,我们评估了五种具有相同阵元数MMM的嵌套阵列的增强自由度(DOF)。因此,定义自由度比率为
γ(M)=M2/Lu(M),(34)\gamma(M) = M^2/L_u(M), \tag{34} γ(M)=M2/Lu(M),(34)
其中MMM是物理阵元数,Lu(M)L_u(M)Lu(M)是(8)中 U\mathbb{U}U 的最大单边孔径。γ(M)\gamma(M)γ(M) 越小,自由度容量就越高。图9展示了当MMM从11变化到100时对应的γ(M)\gamma(M)γ(M)。四种增广嵌套阵列(ANA)比(超级)嵌套阵列拥有更高的自由度。特别是,当 MMM 很大时,两种四级ANA比其他阵列实现了显著的性能提升。根据定理1,γ(M)\gamma(M)γ(M) 的理论最大下界约为2.4。四级ANA γ(M)\gamma(M)γ(M) 值接近3.0。
B. 互耦
在第二个实验中,我们评估了五种嵌套阵列的互耦。图10展示了当MMM从11变化到100时,公式(7)中对应的L(M)L(M)L(M)。嵌套阵列和ANAI-1具有较高的L(M)L(M)L(M)值,这意味着它们遭受着严重的互耦效应。相对而言,ANAII-1可以适度地减少互耦。然后,超级嵌套阵列、ANAI-2和ANAII-2则具有非常低的L(M)L(M)L(M)值。总的来说,ANAII-2能够将互耦比保持在与4阶超级嵌套阵列相同的水平。毫无疑问,通过增加超级嵌套阵列的阶数,互耦可以被进一步降低,结果是其权重函数可以低于ANAII-2。因此,ANA(增广嵌套阵列)的重点在于自由度与互耦之间的综合平衡的重要性。
C. 使用ESPRIT算法进行DOA估计
在第三个实验中,基于互耦模型(5),我们展示了使用空间平滑[31]和ESPRIT算法,用于由14个传感器组成的九种非均匀线性阵列(NLA)
(注:AMRA={0,1,2,8,15,16,26,36,46,56,59,63,65,68}\mathbb{A}_{\text{MRA}} = \{0, 1, 2, 8, 15, 16, 26, 36, 46, 56, 59, 63, 65, 68\}AMRA={0,1,2,8,15,16,26,36,46,56,59,63,65,68}, w(1)=3w(1)=3w(1)=3, w(2)=2w(2)=2w(2)=2, w(3)=2w(3)=2w(3)=2)。
对于10种NLA的γ(14)\gamma(14)γ(14)和L(14)L(14)L(14)值列于表II中。使用信噪比SNR=15 dB的2000个快照来估计协方差矩阵。21个不相关的窄带信源在−30∘-30^\circ−30∘和30∘30^\circ30∘之间均匀分布,如垂直网格线所示。在所有后续仿真中,信源数量均设置为其真实值。图11列出了所有估计的波达方向(DOA)。由于其有限的自由度(DOF)容量,使用互质阵列的估计结果中存在明显的虚假DOA。 (CADiS)嵌套阵列和ANAI-1拥有更高的自由度,因此,它们的性能得到了适度的改善。结果中仍存在一些带有小误差的虚假DOA和漏检的DOA。MRA, SNA, ANAI-2和ANAII-2拥有增加的自由度和降低的互耦,因此,这四种NLA可以清晰且正确地检测到所有信源。
L(14)L(14)L(14) 是互耦比,γ(14)\gamma(14)γ(14) 是自由度比例,这两个值都是越低越好。