一种带激励幅度约束的低副瓣波束赋形方法

摘要

本发明提供一种带激励幅度约束的低副瓣波束赋形方法，在预设动态范围或者预设激励幅度的前提下，实现副瓣电平最小化的约束波束赋形。本发明直接对各阵元激励幅度进行约束，通过引入缩放因子，得到满足任意赋形波束下的阵列激励权重，是带幅度约束波束赋形的一般化方法，适用于激励动态范围约束或幅度预设问题，避免波形约束和激励幅度约束的能量失配问题，相较于传统的迭代方法，在相同动态范围约束或给定幅度约束的条件下，有效降低了阵列天线的副瓣电平，适用范围更广。

说明书

技术领域

本发明涉及阵列天线波束赋形技术，特别涉及带激励幅度约束的波束赋形技术。

背景技术

阵列天线波束赋形是一类非凸、高维、非线性问题，其涉及电磁学和数值优化相关学科，在雷达、无线通信和遥感等领域具有广泛的应用前景。

阵列单元激励幅度约束主要包括：激励动态范围约束和预给定幅度约束。针对动态范围约束，可以有效降低阵元互耦和功分器设计复杂度；针对预给定幅度值约束，设计满足相应波形的激励幅度。

传统波束赋形方法主要包含三类：经典解析方法、智能优化算法和数值优化方法。对于经典解析方法，包括切比雪夫综合法、泰勒综合法和Woodward-Lawson综合法等方法，不能解决任意波形和激励幅度有效约束；对于智能优化算法，包括遗传算法、粒子群算法和免疫算法等方法，不能解决大规模阵列问题，并且难以施加复杂波形和激励幅度约束。数值优化方法由于其理论体系较为完整，在信号处理和机器学习领域引起关注，通过构建目标函数和波形约束，实现相关波束赋形。

幅度约束问题通常分为两类，动态范围约束和约束激励幅度约束。在归一化方向图约束的前提下，直接约束激励动态范围或预设激励幅度值时，方向图和激励幅度存在能量失配，即在该激励幅度约束的前提下，不能优化得到满足该赋形波束的激励权重。

数值优化方法分为两个步骤：1.以阵列天线方向图为基础，施加激励幅度约束和和主瓣方向图约束；2.在此基础上，设计迭代数值算法求解，获得各阵元激励，从而实现波束赋形。然而，对于现有波束赋形方法，激励幅度约束通常只考虑激励动态范围约束或者预设幅度值约束，适用范围受限。

传统数值优化方法中带幅度约束的波束赋形方法主要有两类：

一类是B.J,Liu Y,Cheng J,et al.在2016IEEE的电磁研究进展研讨会上公开的减小动态距离比的赋形功率方向图天线阵综合方法Shaped power pattern antennaarray synthesis with reduction of dynamic range ratio[C]//2016Progress inElectromagnetic Research Symposium(PIERS).IEEE,2016以及Fuchs B在2014年天线与传播期刊上公开的凸松弛在阵列综合问题中的应用Application of Convex Relaxationto Array Synthesis Problems[J].IEEE Transactions on Antennas&Propagation,2014,62(2):634-640.中给出了对于半正定规划的波束赋形综合方法，其目标寻找激励权重w满足以下约束条件：

find w

s.t.L

|a

其中，

引入辅助变量ξ，使得各阵元激励幅度满足如下条件：

引入新的未知量

find W

s.t.Tr(Q

Tr(Q

Tr(Q

Tr(Q

Tr(Q

W≥0

rank(W)＝1

其中，

通过将最后一个非凸约束秩约束rank(W)＝1进行松弛，从而通过转化为凸问题进行求解。

缺点：该方法将向量问题转为矩阵求解，增加了计算复杂度，并且不能保证宽波束赋形得到的解rank(W)＝1；同时，该方法没有考虑副瓣电平最小化，且仅仅适用于动态范围约束问题，不能解决一般化的幅度约束问题。

另一类是Fan,Xuhui,Liang,et al.在2019年天线与传播期刊上公开的基于ADMM(交替方向乘子法)的最小动态范围赋形波束综合方法Shaped Power Pattern SynthesisWith Minimization of Dynamic Range Ratio[J].IEEE Transactions on Antennas andPropagation,2019，其中构造问题如下：

L

|w

其中，

再通过引入附加变量y

L

|v

y

z

v

其中，

该方法是在给定副瓣电平约束前提下，通过ADMM框架方法，迭代优化使动态范围达到最小化。

缺点：基于ADMM的动态范围最下化方法不能解决针对副瓣电平优化和预设激励幅度下的方向图优化问题。

在工程应用中，通常需要在预设动态范围或者预设激励幅度的前提下，实现副瓣电平最小化。

发明内容

本发明所要解决的技术问题是，在预设动态范围或者预设激励幅度的前提下，实现副瓣电平最小化的约束波束赋形方法。

本发明为解决上述技术问题所采用的技术方案是，一种带激励幅度约束的低副瓣波束赋形方法，如图1所示，包括以下步骤：

步骤1确定阵列天线中的总阵元数为N，阵列天线的远场辐射方向图为a(θ)

设置初始值s

步骤2固定{u

其中，A为向导矢量共轭转置矩阵

步骤3固定{s

步骤4固定{s

步骤5固定{s

其中，B为对称矩阵，B＝A

再根据情况利用辅助变量

5a)激励动态范围约束：

5b)预设幅度约束：

步骤6计算

步骤7计算

||·||

步骤8判断是否满足

步骤9根据确定的阵列激励权重w完成波束赋形。

本发明直接对各阵元激励幅度进行约束，是带幅度约束波束赋形的一般化方法，适用于激励动态范围或幅度问题，相较于传统的迭代方法，在相同动态范围约束或给定幅度约束的条件下，有效降低了阵列天线的副瓣电平。

本发明的有益效果是，直接针对每个阵元的激励幅度进行灵活约束，通过引入缩放因子，避免波形约束和激励幅度约束的能量失配问题，得到满足任意赋形波束下的阵列激励权重，进而实现低副瓣波束要求，适用范围更广。

附图说明

图1为本发明实施流程示意图；

图2不同动态范围约束的平顶方向图；

图3副瓣电平随动态范围变化曲线；

图4不同动态范围的激励幅度；

图5笔形波束方向图；

图6阵元激励幅度。

具体实施方式

本发明的实施流程如1所示。

以线阵列天线为例，平面阵列天线或更高维阵列天线的结论可以以此类推。假定天线具有任意分布特性的N个阵元(均匀的或非均匀的)，则阵列天线的远场辐射方向图可描述为：

F(θ)＝a(θ)

其中

低副瓣波束赋形问题表示为：

s.t.L

|a

其中，L

本发明引入缩放因子，在任意赋形波束约束和激励幅度约束的前提下，通过优化缩放因子，实现波形约束和激励幅度约束的能量完美匹配。

将激励幅度约束的低副瓣波束赋形优化问题表示为：

其中,η表示待优化的目标副瓣电平值，D

将动态范围约束低副瓣波束赋形问题和预设幅度约束的低副瓣波束赋形问题进行统一，通过定义缩放因子s，避免波形约束和激励幅度约束的能量失配问题，实现任意赋形波束和幅度约束匹配。

引入辅助变量u＝[u

s.t.L

|v

1≤|w

u

v

该问题是一个非凸非线性问题，不能通过凸优化方法和ADMM方法求解。由于存在变量耦合约束，采用惩罚对偶分解框架，定义主瓣区域和副瓣区域的对偶变量ξ＝[ξ

其中，A

算法的具体流程如下：

步骤1设置初始值s

步骤2固定变量{u

其中

因此，根据一阶最优化条件可得：

其中||·||

步骤3固定变量{s

s.t.L

其中辅助变量

因此，将

步骤4固定变量{s

其中

步骤5固定变量{s

其中

构造函数r(w)的上界函数

其中，B＝A

因此可以下面两种情况的最优解：

1)激励动态范围约束：

2)预设幅度约束：

步骤6计算

步骤7计算

步骤8若

步骤9根据确定的阵列激励权重w用于波束赋形。

实验验证

采用如下30阵元均匀分布的线阵列天线，假设每个天线单元都是全向的，即g

表1 30阵元阵列天线的位置信息

Case 1：不同动态范围约束下的波束赋形

本发明实验考虑平顶波束赋形，主瓣区域为[-10°,10°]，副瓣区域为[-90°,-15°]和[15°,90°],主瓣纹波为0.6dB的平顶波束。动态范围区间为[1,6],通过上述算法流程迭代技术可求得阵列激励权重(各阵元激励权重)，并计算相应的阵列方向图，如图2所示，在不同的激励动态范围约束前提下，主瓣内的纹波小于等于0.6dB，满足预设主瓣约束要求，同时，在动态范围增加时，副瓣电平降低；图3表示不同动态范围下的副瓣电平变化情况，在主瓣区域和纹波大小不变的前提下，副瓣电平随动态范围增加而降低；图4表示不同动态范围下各阵元激励幅度变化情况。

Case2：幅度预设约束下的波束赋形

本发明实验为笔形波束赋形，波束指向为5°，副瓣区域为[-90°，0°]和[10°，90°]，激励幅度预先给定。通过上述算法流程，迭代计算可得各阵元激励权重，并计算相应阵列方向图。图5表示笔形波束方向图；图6表示预给定的各阵元激励幅度和提出的激励幅度，表明各个阵元能够满足预设幅度要求。