一种方向图可重构的稀布直线阵综合方法与流程

本发明涉及无线通信领域，特别涉及一种方向图可重构的稀布直线阵综合方法。

背景技术：

近年来，可重构阵列天线因其仅需改变阵元激励就可实现多个方向图而被广泛应用于雷达、遥感、航天器及无线通信等领域(参考文献[1]：a.f.morabito,a.massa,p.rocca,andt.isernia,"aneffectiveapproachtothesynthesisofphase-onlyreconfigurablelineararrays,"ieeetrans.antennaspropagat.,vol.60,no.8,pp.3622-3631,aug.2012)。众所周知，阵元位置的非均匀排布可以有效地节省阵元数目，从而有利于减轻阵列天线整体重量、简化馈电网络以及降低应用系统的成本和造价(参考文献[2]：x.zhao,y.zhang,q.yang,"ahybridalgorithmforsynthesizinglinearsparsearrays,"prog.electromagn.res.c,vol.63,pp.75-83,2016)。目前，非均匀阵列天线主要分为稀疏阵列天线和稀布阵列天线这两类，其中稀疏阵列天线是指从相邻元间隔为半波长的均匀满阵中，挑选出起作用的阵元而构成，与之不同，稀布阵列天线中阵元不再被约束在半波长的栅格上，而是可以更为随机地分布在给定的阵列口径中。因此，相比于稀疏阵列天线，稀布阵列天线具有更多的优化自由度来以更少的阵元数实现期待的辐射特性(参考文献[3]：x.zhao,q.yang,y.zhang,“ahybridmethodfortheoptimalsynthesisof3-dpatternsofsparseconcentricringarrays,”ieeetrans.antennaspropagat.,vol.64,no.2,pp.515-524,jan.2016)。

由压缩感知理论可知，稀布直线阵的综合问题可等效为高维稀疏信号的重建问题，在此基础上，本申请人已提出了基于迭代加权l1范数最小化的稀布直线阵综合方法(参考文献[4]：x.zhao,q.yang,y.zhang,"compressedsensingapproachforpatternsynthesisofmaximallysparsenon-uniformlineararray,"ietmicrowav.antennaspropag.,vol.8,pp.301-307,2014)，该方法可以利用凸优化算法对阵元数、阵元位置和相应的阵元激励这些阵列参数实现整体优化，从而最大程度的节省阵列所需的阵元数目，但是该方法仅涉及单一方向图要求的阵列综合问题。

此外，公开的研究表明，众多的综合方法已广泛地应用于具有单一方向图的稀布阵列天线的设计中，而关于方向图可重构的稀布阵列天线的研究相对较少，这是由于这类天线的综合更为复杂，是一个多约束、多参数非线性优化问题，需要根据多个方向图的辐射要求，确定所需的最少阵元数目、相应的阵元位置以及多个方向图对应的各组阵元激励。最近，扩展的矩阵束方法(参考文献[5]：y.liu,q.h.liu,andz.nie,“reducingthenumberofelementsinmultiple-patternlineararraysbytheextendedmatrixpencilmethods,”ieeetrans.antennaspropag.,vol.62,no.2,pp.652-660,feb,2014.)和多次测量矢量欠定系统局域解方法(参考文献[6]：f.yan,p.yang,f.yang,l.zhou,andm.gao,“synthesisofpatternreconfigurablesparsearrayswithmultiplemeasurementvectorsfocussmethod,”ieeetrans.antennaspropag.,vol.,no.,pp.,2016.)被应用于综合方向图可重构的稀布直线阵，尽管这两种方法可以有效地求解这类复杂的阵列综合问题，但其均需预先给定具体的参考方向图，即需要将已知的多个参考方向图作为算法的先验信息，而无法应用于仅根据期待的多个方向图包络即多个主波束赋形和共同的副瓣电平要求来综合可重构稀布直线阵，且在实际工程中，根据方向图包络要求来综合可重构稀布直线阵更具有普遍实用性。

技术实现要素：

本发明的目的在于克服现有技术需要将已知的多个参考方向图作为算法的先验信息，而无法应用于仅根据期待的多个方向图包络即多个主波束赋形和共同的副瓣电平要求来综合可重构稀布直线阵的缺陷，从而提供一种更具有普遍实用性的稀布直线阵综合方法。

为了实现上述目的，本发明提供了一种方向图可重构的稀布直线阵综合方法，包括：

步骤s01、根据给定的阵列口径长度，以微小步长对其进行均匀离散化，在每个离散点上都排布一个阵元，从而构造一个虚拟的过采样均匀间隔直线阵；确定所述虚拟的过采样均匀间隔直线阵的阵元位置矢量d；

步骤s02、根据所述虚拟的过采样均匀间隔直线阵的阵元位置矢量d，确定该虚拟的过采样均匀间隔直线阵可重构的第q个方向图对应的激励列矢量w^(q)，q＝1,…,q，以及由这q组激励列矢量组成的激励矩阵w；然后由位置矢量d、激励列矢量w^(q)和观测角度θ共同决定的方向图函数，确定所述虚拟的过采样均匀间隔直线阵的导向矢量a(θ)；

步骤s03、根据联合稀疏恢复理论，以最小化激励矩阵w的混合l2/l1范数为优化目标，以q个方向图设计要求为约束条件，建立综合q组方向图的凸优化模型；

步骤s04、利用凸优化算法求解步骤s03建立的凸优化模型，计算出激励矩阵w，由此确定阵元数目、阵元位置以及可重构的q个方向图对应的q组阵元激励。

上述技术方案中，步骤s01具体包括：

以离散步长δ均匀离散化给定的阵列口径l，形成n个足够密的均匀栅格点，在每个栅格点上放置一个天线单元，构成虚拟的n元过采样等间隔直线阵，其中δ<<λ/2，λ为工作波长；该虚拟阵的位置矢量d＝[d1,d2…,dn]，其中第n个栅格点的位置dn＝(n-1)δ，n＝1,2,…,n。

上述技术方案中，步骤s01中，所述离散步长δ根据给定的阵列口径l来取值，以保证离散点足够密，在相应的阵列口径上形成一个过采样均匀间隔阵，对于l≤50λ，δ＝0.001λ。

上述技术方案中，步骤s02具体包括：

根据所述虚拟的过采样均匀间隔直线阵的阵元位置矢量d中的元素值与激励矢量w^(q)中的元素值一一对应的关系，由所述阵元位置矢量d确定相应的第q个方向图对应的激励矢量w^(q)＝[w1^(q),w2^(q)…,wn^(q)]^t，其中，q＝1,…,q，t表示转置运算；

由q组激励矢量构成激励矩阵w，激励矩阵w的表达式为：

由阵元位置矢量d中的阵元位置和激励矢量w^(q)中的阵元激励确定虚拟阵方向图，所述虚拟阵方向图的表达式为：

其中θ为观测角，θ∈[-90°,+90°]；

虚拟阵方向图的矢量形式则为：

f^(q)(θ)＝a(θ)w^(q)(3)

由此确定导向矢量a(θ)，其表达式为：

上述技术方案中，所述激励矢量w^(q)为稀疏矢量，若第n个激励值wn^(q)＝0，则表明对应位置dn上不存在阵元。

上述技术方案中，步骤s02中，所述激励矩阵w是行稀疏矩阵。

上述技术方案中，步骤s03具体包括：

首先，定义混合l2/l1范数为：

其中w^n→表示激励矩阵w的第n行，且w^n→的l2范数对应为：

然后，根据给定的期望设计的q个赋形主波束fd^(q)(θ)和共同的副瓣电平上限要求ub，结合混合l2/l1范数以及步骤s02得到的激励列矢量w^(q)、导向矢量a(θ)，建立多方向图稀布直线阵综合的凸优化模型，该凸优化模型的表达式为：

min||w||2,1

|a(θ)w^(q)|≤ub,θ∈sideloberegion

q＝1,2,…q(7)

其中ε为赋形误差，取10^-6；shapedreigon表示赋形区域，sideloberegion表示副瓣区域。

上述技术方案中，步骤s04具体包括：

利用求解器cvx求解步骤s03建立的凸优化模型，并定义求解器的变量为复数，计算出激励矩阵w，它的非零行数目则为稀布直线阵的阵元数；激励矩阵w各列中非零元的数目和出现的位置相同，第q列中的非零元则为相应的第q个方向图对应的阵元激励；根据阵元激励与阵元位置的一一对应关系，由各列激励矢量相同的非零元索引，从阵元位置矢量d中确定相应的阵元位置。

本发明的优点在于：

1、本发明方法适用于根据期望的多个方向图包络来综合可重构稀布直线阵，避免了要将已知的参考方向图作为综合设计的先验信息；

2、通过联合稀疏恢复理论，本发明方法将多约束、多目标、多参数的方向图可重构稀布直线阵列综合问题转化为最小化混合l2/l1范数的凸优化问题，从而采用凸优化算法对这样的复杂阵列综合问题进行有效求解；

3、本发明方法可以实现对阵列参数的整体优化，通过同时优化阵元数目、阵元位置以及不同方向图对应的各组阵元激励，在节省阵元数的同时，仅更新阵元激励来实现多方向图的重构。

附图说明

图1是本发明的一种方向图可重构的稀布直线阵综合方法的流程图；

图2是本发明在阵列口径为9.5λ上设计的14元稀布直线阵综合的双方向图；

图3是本发明在阵列口径为9.5λ上设计的14元稀布直线阵的阵元位置与激励分布图。

具体实施方式

现结合附图对本发明作进一步的描述。

如图1所示，本发明的一种方向图可重构的稀布直线阵综合方法包括以下步骤：

s01：构造虚拟的过采样均匀阵；

根据给定的阵列口径长度，以微小步长对其进行均匀离散化，并假定每个离散点上都排布一个阵元，从而构成一个虚拟的过采样均匀间隔直线阵，据此确定相应的阵元位置矢量d。

步骤s01具体包括如下过程：

采用微小步长δ(δ<<λ/2，λ为工作波长)均匀离散化给定的阵列口径l，形成n个足够密的均匀栅格点，假定每个栅格点上放置一个天线单元，则构成虚拟的n元过采样等间隔直线阵，该虚拟阵的位置矢量d＝[d1,d2…,dn]，其中第n个栅格点的位置dn＝(n-1)δ(n＝1,2,…,n)。需要指出的是，需要根据给定的阵列口径l，确定离散步长δ的取值，从而保证口径上的离散点足够多，在相应的阵列口径上可以形成一个过采样均匀间隔阵，对于l≤50λ，通常δ＝0.001λ。

s02：确定行稀疏矩阵和导向矢量；

根据虚拟阵的位置矢量d，确定该虚拟阵可重构的第q(q＝1,…,q)个方向图对应的激励列矢量w^(q)和由这q组激励列矢量组成的激励矩阵w。由虚拟阵的位置矢量d、激励列矢量w^(q)和观测角度θ共同决定的多方向图函数，确定该虚拟阵的导向矢量a(θ)；

步骤s02具体包括如下步骤：

在虚拟阵位置矢量d确定的情况下，若该虚拟阵可以仅改变阵元激励来实现q个方向图，那么相应的第q(q＝1,…,q)个方向图对应的激励矢量w^(q)＝[w1^(q),w2^(q)…,wn^(q)]^t，其中t表示转置运算。由此可知，激励矢量w^(q)中的元素值与位置矢量d中的元素值是一一对应的，若第n个激励值wn^(q)＝0，则表明对应位置dn上不存在阵元。尽管假定了足够多的阵元位于阵列口径内，但实际起作用的阵元很少，从虚拟阵中移除那些激励值为零的阵元，则可得到稀布阵，这就要求激励矢量w^(q)是一个稀疏矢量，即非零元的数目远远小于矢量长度。对于可重构q组方向图的稀布直线阵而言，相应的q组激励矢量w^(q)是联合稀疏的，即各组稀疏矢量中的非零元数目和位置相同。

由这q组激励矢量组成的激励矩阵w为：

由于q组激励矢量w^(q)是联合稀疏的，因此激励矩阵w为行稀疏矩阵，起作用的阵元数则对应为矩阵w的非零行数目。

由阵元位置和阵元激励确定的虚拟阵方向图为：

其中θ为观测角，且θ∈[-90°,+90°]。

(2)式的矢量形式则为：

f^(q)(θ)＝a(θ)w^(q)(3)

由此确定的导向矢量a(θ)为：

s03：建立最小化混合l2/l1范数的凸优化模型；

根据联合稀疏恢复理论，以最小化激励矩阵w的混合l2/l1范数为优化目标，以q个赋形方向图设计要求为约束条件，建立多方向图综合的凸优化模型；

步骤s03具体过程如下：

为了加强激励矩阵w各列向量之间的联合稀疏性，定义混合l2/l1范数为：

其中w^n→表示激励矩阵w的第n行，且w^n→的l2范数对应为：

根据给定的q个赋形主波束fd^(q)(θ)(q＝1,…,q，θ∈赋形主波束区域)和共同的副瓣电平上限要求ub，结合混合l2/l1范数以及步骤s02得到的激励列矢量w^(q)、导向矢量a(θ)，建立多方向图稀布直线阵综合的凸优化模型，该凸优化模型为：

min||w||2,1

|a(θ)w^(q)|≤ub,θ∈sideloberegion

q＝1,2,…q(7)

其中，ε为赋形误差，通常取10^-6；shapedreigon表示赋形区域，sideloberegion表示副瓣区域。

s04：输出阵元数目、阵元位置、多方向图对应的各组阵元激励；

利用凸优化算法求解步骤s03建立的优化模型(7)，计算出激励矩阵w，由此确定阵元数目、阵元位置以及可重构的q个方向图对应的q组阵元激励。

步骤s04具体过程如下：

利用求解器cvx求解步骤s03建立的优化模型(7)，并定义求解器的变量为复数，计算出激励矩阵w，它的非零行数目则为稀布阵的阵元数；激励矩阵w各列中非零元的数目和出现的位置相同，第q列中的非零元则为相应的第q个方向图对应的阵元激励；根据阵元激励与阵元位置的一一对应关系，由各列阵元激励矢量相同的非零元索引，从阵元位置矢量d中确定相应的阵元位置。

本发明所提出的一种方向图可重构的稀布直线阵综合方法，可以进一步通过以下具体的仿真实例来验证并说明。

仿真实例：

本实例在给定的阵列口径l＝9.5λ内综合一个双方向图可重构的稀布直线阵，其中一个为笔形波束方向图，另一个为平顶波束方向图，要求笔形波束主波束区域为[-8.7°，+8.7°]，主波束方向θo＝0°，平顶波束的主波束区域为[-20.5°，+20.5°]，且其波纹抖动不超过0.5db，并要求这两个方向图的副瓣电平均小于-25db。已有研究表明，若采用相邻元间隔为0.5λ的均匀满阵，则需要20个阵元；若将该满阵实现的方向图作为参考方向图，利用扩展的矩阵束方法综合稀布直线阵，则需要15个阵元。采用本发明提出的一种方向图可重构的稀布直线阵综合方法进行设计，仿真结果表明，本发明所需的阵元数仅为14，综合的可重构方向图如图2所示，由此可知，本发明方法综合的稀布直线阵实现的方向图满足期望的辐射特性。相应的阵元位置与激励分布如图3所示，由于阵列结构的对称性，表1也仅列出了位于正半轴上的阵元位置和各个波束所对应的阵元激励，由此可知，该稀布直线阵的最小阵元间隔为0.697λ，相邻元之间的耦合很小，几乎可以忽略，这具有非常显著的工程意义和应用价值。

表1

最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。