第19卷 第3期太赫兹科学与电子信息学报Vo1.19，No.3 2021年6月 Journal of Terahertz Science and Electronic Information Technology Jun.，2021 文章编号：2095-4980(2021)03-0442-06
秦自立1a,1b,2，杨冠1a,1b,2，王方力1a,1b,2，李超1a,1b，纪奕才1a,1b,2
(1.中国科学院 a.空天信息创新研究院；b.电磁辐射与探测技术重点实验室，北京 100190；
2.中国科学院大学电子电气与通信工程学院，北京 100049)
摘 要：为了解决多输入多输出(MIMO)天线阵列由于阵元间隔过大造成的阵列方向图出现栅瓣，在雷达回波成像时出现影响目标识别的虚假目标的问题，提出了一种改进的遗传算法对阵列 排布进行优化。在传统标准遗传算法上进行改进，用多个矩阵组合表示MIMO阵列，针对在矩形 平面随机分布的稀疏阵列的方向图旁瓣问题进行优化设计，并采用基于Logistic混沌序列的方法
产生种扰动，避免优化过程进入局部最优状态。通过实例对比22发射天线22接收天线的均匀
规则排布MIMO阵列和经改进遗传算法优化的稀布MIMO阵列，结果表明，改进遗传算法可以有
效解决规则排布阵列方向图中出现的栅瓣，并且降低方向图旁瓣，提高雷达成像性能。该优化算
法变量可控，具有很强的实用性，为MIMO雷达的阵列排布提供了解决方法。
关键词：平面稀布阵列；多输入多输出；改进遗传算法；栅瓣；旁瓣
中图分类号：TN951文献标志码：A doi：10.11805/TKYDA2019552
Modified genetic algorithm for optimizing MIMO sparse array
QIN Zili1a,1b,2，YANG Guan1a,1b,2，WANG Fangli1a,1b,2，LI Chao1a,1b，JI Yicai1a,1b,2 (1a.Aerospace Information Research Institute；1b.Key Laboratory of Electromagnetic Radiation and Sensing Technology，
Chinese Academy of Sciences，Beijing 100190，China；
2.School of Electrical and Communication Engineering，University of Chinese Academy of Sciences，Beijing 100049，China)
Abstract：The array pattern of Multiple Input Multiple Output(MIMO) antenna array has grating lobes due to the large array element spacing, and the false targets that affect the target recognition appear
in the radar echo imaging. An improved genetic algorithm is proposed to optimize the array arrangement.
The traditional standard genetic algorithm is improved to represent MIMO array with multiple matrix
combinations, and the pattern sidelobes of sparse array with random distribution in rectangular plane are
optimized. The method based on Logistic chaotic sequence is adopted to generate population disturbance
and avoid the optimization process entering into the local optimal state. A uniform regular arrangement of
MIMO array with 22 transmitting antennas and 22 receiving antennas and a modified MIMO array optimized
by modified genetic algorithm are compared by examples. The results show that the modified genetic
algorithm can effectively avoid the grating lobes in the regular array pattern, reduce the side lobes of the
pattern, and improve the radar imaging performance. The optimization algorithm has strong practicability
for its variable is controllable, which provides a solution for array arrangement of MIMO radar.
Keywords：planar sparse array；Multiple Input Multiple Output；modified genetic algorithm；
grating lobe；side lobe
多输入多输出(MIMO)雷达近年来得到了越来越多研究学者的关注[1−2]，由于其独特的多个发射天线独立发射正交信号并且多个接收天线同时接收信号的特点，其获得了远超过实际阵元数目的观测通道，具有很高的数据获取效率，相比同样孔径的单站雷达系统，远远节约了硬件成本，以空间积累代替时间积累，从而具有缩短积累时间，实现单次快拍成像的优势。
收稿日期：2019-12-18；修回日期：2020-02-18
基金项目：国家自然科学基金资助项目(61671432)
作者简介：秦自立(1995–)，男，在读硕士研究生，主要研究方向为新型天线设计、毫米波雷达设计。email:****************
第3期 秦自立等：改进遗传算法优化MIMO 稀布阵列 443
稀布MIMO 雷达成像紧随着MIMO 成像技术在各领域的发展而发展，因为进一步降低了成本，减小系统复杂度，以尽可能小的成本实现尽可能大的目标，也是各类阵列成像技术走向实用、获得推广的必经之路[3−4]。稀布平面MIMO 阵列在实际应用中存在的问题是由于阵列单元间隔大于半波长，阵列方向图出现栅瓣，并且对于阵列的副瓣电平难以控制。本文采用遗传算法对阵列排布进行优化，目的在于去除阵因子方向图栅瓣，降低副瓣。遗传算法在对天线阵列进行优化方面早有研究。文献[5−6]是较早使用遗传算法对稀疏阵列进行优化的实例，文献[7]中提出采用改进遗传算法并结合降维算法优化椭圆环阵列的旁瓣特性，文献[8]提出了一种应用于线性天线阵综合中多目标优化遗传算法，文献[9]针对传统遗传算法在迭代后期容易陷入局部最优的问题，提出采用混沌理论和遗传算法相结合的方法对T 形MIMO 线阵的低旁瓣进行优化。文献[10]引入距离摄动，同时对发射天线和接收天线进行几何修正，优化稀疏MIMO 直线阵列方向图，实现了具有较低的相对旁瓣电平和最大自由度的模式合成。遗传算法对天线阵列性能优化设计的研究已有较长历史，对MIMO 线阵的研究也有较多的文献，但是对于平面稀疏MIMO 阵列的研究并不多见。本文介绍一种针对平面稀布MIMO 天线阵列方向图进行优化的遗传算法，将传统平面阵列收发一体改进为收发分离的MIMO 阵列，用两个二维矩阵来对天线单元进行分类以及定位，可优化任意排布的面阵方向图。将此算法进行改进还可用于对一维线性稀布MIMO 阵列以及更为复杂的三维空间稀布MIMO 阵列方向图进行优化。
1 改进遗传算法
遗传算法是模拟生物种进化过程中优胜劣汰、适者生存法则的全局随机优化算法。类比生物繁衍过程，种中的大量个体进行基因交叉和基因突变，选择性能优异的后代个体，对保留的个体进行繁殖后代，直到满足需求的个体出现。作为全局性的概率优化算法之一，遗传算法适用范围广，简单易操作，对于解决传统搜索方法中难以处理的复杂非线性问题有着很强的适用性[5]。 本文中的改进遗传算法是在传统一维标准遗传算法的编码方式，种的生成方法，选择、交叉、变异算子的设计基础上进行改进的二维遗传算法，针对在矩形平面随机分布的稀疏阵列的方向图旁瓣问题进行优化设计，约束条件包括阵元数、阵列孔径和阵元间距。采用实值编码优化变量，将二维平面上的点的横纵坐标对改写成复数的实部和虚部，摒弃了传统的二进制编码方式，这些代表阵元位置的复数构成了个体的二维基因序列。在进行遗传操作前对种进行预处理，遗传操作结束后还需对种进行后处理[11]。 1.1 初始种的生成
假定初始种由M 个个体组成，个体的有用基因为阵元数N ，其中发射单元数为t N ，接收单元数为r N ，满足条件t r N N N =+。基因用复数(j )n n x y +表示，代表着阵元的平面坐标，每个个体均为Q 行P 列的复数矩阵，对于稀布的平面阵列，N <P Q ⨯。每个基因需要满足的空间约束条件如下：
<
/2/2, /2/2 ,, 1,, y m x
n n m n n L x L H d y H x x y y m n Z m n P Q n d m ⎧⎪-⎪⎨
-⎪
⨯⎩
-∈-⎪≠≤≤≤≤≥≥≤≤ (1) 式中：x d 为x 方向的最小间距；y d 为y 方向上的最小间距；L 和H 为阵列在x 和y 方向的最大口径。对于x 方向的P 个阵元，满足最小阵元间隔之后x 方向总长度L 还剩余的区间长度为：
(1)x x S L P d =--⨯ (2)
在这个区间中，P 个阵元随机分布。同理，y 方向的Q 个阵元随机分布的区间为：
(1)y y S H Q d =--⨯ (3)
将个体的产生约束在区间，相比较原始阵列排布空间L H ⨯，算法搜索空间减小，更容易寻到更优解。在
[]0,x S 区间上生成
Q 行P 列的实数矩阵并将每行的元素按从小到大重新排序即可得到约束后的阵元x 坐标矩阵
R ，同理在0,y S ⎡⎤⎣⎦区间上生成Q 行P 列的实数矩阵并将每列的元素按从小到大重新排序即可得到约束后的阵元y 坐标矩阵I 。定义Q 行P 列的个体平面约束矩阵C 为：
(),[(1)/2]j [(1)/2]x y m n n d L m d H =-⨯-+⨯-⨯-C (4)
因此，一个满基因个体就可以表示为：
j =+G R I +C (5)
444
太赫兹科学与电子信息学报 第
19卷
矩阵G 为未稀疏的普通平面阵列，为了得到稀布平面MIMO 阵列，还需要对上述个体矩阵G 进行两步处理，即阵元稀疏化以及分类。
定义个体的特性矩阵S 以及个体稀疏矩阵U ，矩阵S 是包含3种元素(0,1,-1)的Q 行P 列矩阵，0代表满基因个体复数矩阵G 中需要去掉的阵元，即稀疏化，对应G 中相应位置元素取0得到U 中元素，1和−1分别代表发射天线和接收天线阵元，其个数分别为发射天线和接收天线的个数t N 和r N ，对应G 中相应位置元素不变得到U 中元素。个体稀疏矩阵U 与特性矩阵S 共同组成了一个个体所需的所有基因。
例如，若个体特性矩阵S 为：
10011011-⎡⎤⎢⎥-⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦
S (6)
则对应的个体稀疏矩阵U 为：
111111212111j 00j j 0j j P P Q Q QP QP x y x y x y x y x y ++⎡⎤⎢⎥
+⎢⎥
⎢
⎥=⎢⎥
阵列天线
⎢⎥⎢⎥++⎣⎦
U (7) 多个随机个体组成了种，个体数量决定了种的复杂度。与个体相对应，分别定义种特性矩阵M S 、种稀疏矩阵M U 、种平面约束矩阵M C ，均为M 页Q 行P 列矩阵，第m 页则代表第m 个个体。
1.2 适应度函数
本文以去除平面稀布MIMO 阵列阵因子方向图栅瓣、最小化方向图旁瓣为目标，稀布平面MIMO 阵列的优
化模型为：
{}{}
MIMO min max [(,)]s.t. /2/2/2/2
,, 1, , ,n n m n m x
y
n SLL E L x L H y H x x y y m n Z m n N m n d d
θϕ⎧⎪--⎪
-⎨⎪-⎪
∈≠⎩
≤≤≤≤≥≥≤≤ (8) 式中：MIMO E 表示MIMO 阵列的阵因子方向图；SLL 表示方向图旁瓣。
计算平面MIMO 阵列方向图时先计算出MIMO 阵列的虚拟阵列阵元位置，再用计算普通平面阵列方向图的方法得到原始MIMO 阵列的方向图。对于发射单元和接收单元数量分别为t N 和r N 的MIMO 阵列，虚拟阵元数为t r N N ⨯，可由发射单元和接收单元的位置计算得到虚拟阵列[3]：
t r
()/2()/21,2,,i j i j n t r n t r x x x y y y n N N =+⎧⎪
=+⎨⎪=⨯⎩ (9) 可以得到MIMO 阵列方向图为
MIMO (,)exp[j ()]n n n
E k x u y v θϕ=+∑ (10)
式中：2π
k λ
=
；sin cos u θϕ=；sin sin v θϕ=。
个体的适应度为最大旁瓣电平的绝对值，适应度函数值越大则代表性能更优，适应度函数表达式为：
{}MIMO max |[(,)]|SLL E fit θϕ= (11)
1.3 遗传操作
遗传算法与大多数全局优化算法都存在相同的问题，即优化过程容易陷入局部最优并且后期收敛速度变
慢。为解决此问题，本文采用一种基于Logistic 混沌序列的扰动方法[12]。若优化迭代过程进入长时间的局部最优状态(在此以最优状态经多次迭代后不发生变化为判定准则)，则对体产生一个大的扰动，其表达式如下：
第3期
秦自立等：改进遗传算法优化MIMO 稀布阵列 445
11
()
()n x n x n n y n x n x x S x y y S y μμ++=-⎧⎨=-⎩ (12)
式中：4/x x S μ=,4/y y S μ=；n x ,n y 为当前个体阵元的位置坐标，对应的1n x +,1n y +为产生扰动后阵元的位置。 容易证明，扰动之后的阵元依旧满足空间约束条件式(1)。对扰动后的体进行阵元重新排序，排序只针对未被稀疏的元素，将矩阵元素的实部按行从小到大重新排序，将矩阵元素的虚部按列从小到大重新排序后重新组成种。
对种进行的遗传操作都在约束区间x y S S ⨯上进行，因此在进行遗传操作之前需要对种稀疏矩阵进行预处理得到可操作矩阵P ，预处理表达式为：
M M =-P U C (13)
遗传操作后需将得到的矩阵'P 进行后处理得到子代种'M U ，后处理表达式为：
M
M '''=+U P C (14) 式中：'P 是P 经过遗传操作后得到的矩阵；M
'C 是'P 对应的约束矩阵。 遗传操作的选择过程实质是按照某种方法有概率地选择一部分父代个体进入到子代中。赌法是较为常
用的方法，是一种回放式的随机抽样方法，个体被选择从而可以进入下一代继续繁殖的概率等于其适应度函数值占种适应度函数值的总和的比例。此种方法
在一定程度上可以避免优化过程进入局部最优状态，但是选择过程的误差较大。在此采用锦标赛
选择方法[13]，每次从种中选出一定数量的个
体，从中选择适应度函数值最高的个体进入子代种之中。重复该操作即可得到新的种。
在对P 进行交叉操作时，依据交叉概率p c
选择需要进行操作的个体，在其中选择任意两个(即矩阵P 的某两页以及对应的特性矩阵S )，随机选择个体的一个或多个非稀疏元素进行交换，
即单点交叉或者多点交叉。将得到的两个个体矩阵元素进行排序，排序方法与产生种扰动后的排序方法相同。
对交叉操作后的子代个体，依据变异概率p m 选择需要进行操作的个体，可进行单点变异
或者多点变异。在个体中随机选择非稀疏元素，用区间x y S S ⨯上的随机复数代替。将变异后的个体矩阵元素进行排序，排序方法与上述方法相同，最终得到遗传操作后的子代矩阵'P ，对'
P 进行后处理即可得到子代种。
多次重复上述的遗传操作，即可得到满足需求的后代个体。由于所有的遗传操作都是在进行约束后的种中进行，得到的子代体中的基因均满足间隔约束条件，都是可行解。算法流程图如图1所示。
1.4 算法的延展
鉴于稀布MIMO 雷达并不局限于平面阵列，此方法亦可应用于一维直线阵列稀布MIMO 和三维空间稀布MIMO 阵列。对于直线稀布MIMO 阵列，其空间自由度为1，在生成种时位置坐标用实数即可表示，与平面阵列相比，代表种的矩阵对应地就减少一个维度。三维稀布MIMO 阵列与之类似，只是代表种的矩阵相应要增加一个维度。此时天线单元的位置坐标不再用复数表示，可用一个包含三个元素的数组表示其空间位置。对应的算法复杂度增加，运行时间加长。
Fig.1 Flow chart of algorithm
图1 算法流程图
446
太赫兹科学与电子信息学报 第
19卷
2 阵列优化实例
本节以22发22收的平面稀布MIMO 阵列为例，P =Q =10，稀布率为(22+22)/(10×10)=44%，交叉概率p c =0.8，变异概率p m =0.1，随机五点交叉和三点变异，初始种个体数为100，工作频率为77 GHz ，
单元间隔1cm 2.5x y d d λ===。阵元分布空间为L =11cm,H =11cm ，单元天线等幅同相，主波束指向法向。
当阵列规则排布时，发射阵元、接收阵元分布以及虚拟阵列如图2(a)所示，阵因子方向图如图2(b)所示，φ=0°,φ=90°两个主平面方向图如图3所示。可以看出阵列方向图存在大量栅瓣，旁瓣电平为−13.29 dB 。
采用上述改进遗传算法对22发22收平面稀布MIMO 阵列进行阵因子方向图优化，表1中列出采用改进遗传算法对阵列排布进行5次优化分别得到的阵列方向图旁瓣电平最优值。优化后发射阵元、接收阵元分布以及虚拟阵列如图4(a)所示，阵因子方向图如图4(b)所示，φ=0°,φ=90°两个主平面方向图如图5所示。从图中可以看出，优化后阵列排布呈无序排布，阵因子方向图旁瓣消除，方向图旁瓣电平为−16.62 dB ，φ=0°平面的方向图旁瓣为−16.87 dB ，φ=90°平面的方向图旁瓣为−16.76 dB 。同均匀矩形阵列相比，旁瓣电平降低了3.33 dB 。
将两种稀布MIMO 阵列进行成像仿真对比，对距离天线阵列5 m 处的点目标进行成像，仿真结果如图6所示。从图6(b)可以看出，矩形排布MIMO 阵列成像中出现虚假目标，且虚假目标幅度与真实点目标幅度相当，使得无法分辨真假目标。而图6(c)所示的经过优化后的平面稀布MIMO 阵列成像结果，在准确的位置处识别出目标，且不存在虚假目标影响，验证了本文方法的有效性。
Fig.5 Patterns of principal planes
图5 主平面方向图:
(a) φ=0° plane -80 -60 -40 -20 0 20 40 60 80
θ/(°)
0 -10 -20 -30 -40
-50
r e l a t i v e p a t t e r n /d B
-80 -60 -40 -20 0 20 40 60 80
θ/(°)
0 -10 -20 -30 -40
-50
r e l a t i v e p a t t e r n /d B
(b) φ=90° plane x =-5 y =-16.37
x =-32 y =-16.76
Fig.4 Random arrangement of MIMO array
图4 随机排布MIMO 阵列
-0.10 -0.05 0 0.05 0.10 x /m
(a) array element distribution and virtual array arrangement 0.10
0.05 0
-0.05
-0.10 y /m
(b) array pattern 1 u =sin θ cos φ v =sin θ sin φ -1 -1 1
0 0 r e l a t i v e p a t t e r n /d B
0 -10
-20
-30 -40 TX RX virtual array
Fig.2 Rectangular MIMO array
图2 矩形MIMO 阵列
-0.10 -0.05 0 0.05 0.10
x /m
(a) array element distribution
0.10
0.05 0
-
0.05 -0.10 y /m
(b) array pattern 1
u =sin θ cos φ v =sin θ sin φ -1 -1 0 0 r e l a t i v e p a t t e r n /d B
-10
-20
-30 -40 TX RX virtual array
1 Fig.3 Principal planes pattern
图3 主平面方向图
(a) φ=0°plane
(b) φ=90°plane
-80 -60 -40 -20 0 20 40 60 80
θ/(°) 0
-10 -20 -30 -40 -50
r e l a t i v e p a t t e r n /d B
0 -10 -20 -30 -40 -50
r e l a t i v e p a t t e r n /d B
x =3
y =-13.29
-80 -60 -40 -20 0 20 40 60 80
θ/(°)
x =3
y =-13.29
表1 多次优化结果 Table1 Multiple optimization results
no. optimal value 1 16.14 2 16.44 3 16.22 4 16.62 5 16.04
Fig. 6 Comparison of point target imaging simulation 图6 点目标成像仿真对比 (a) point target model (b) simulation results of rectangular
MIMO array (c) simulation results of sparse MIMO array 1 u =sin θ cos φ v =sin θ sin φ -1 -1 1 0 0 z /m
6.0 5.5
5.0 4.5
4.0 3.5 3.0 - -3 -2 -1 0 1 2 3 x /m 3 2 1 0 -1 -2 -3 y /m
-3 -2 -1 0 1 2 3 x /m 3 2
1 0 -1
-2 -3 y /m
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议