基于蚁群优化算法的相控阵天线随机馈相方案的优化.pdf

第4 4卷 第3期2 0 2 3年9月制 导 与 引 信GU I D AN C E&F U Z E V o l.4 4 N o.3S e p.2 0 2 3文章编号:1 6 7 1-0 5 7 6(2 0 2 3)0 3-0 0 2 3-0 6基于蚁群优化算法的相控阵天线随机馈相方案的优化李秋生1,2,杨雨萱2(1.赣南师范大学智能控制工程技术研究中心,江西 赣州 3 4 1 0 0 0;2.赣南师范大学物理与电子信息学院,江西 赣州 3 4 1 0 0 0)摘 要:为了降低相控阵天线的峰值副瓣电平,提高相控阵天线的抗干扰性能,以泰勒幅度加权的线阵为例,分析了随机馈相的数学模型,并基于蚁群优化算法提出了针对相控阵天线随机馈相的优化算法,并进行了仿真实验验证。仿真结果表明,采用优化后的随机馈相方案,相控阵天线的峰值副瓣电平相比优化前降低了8 d B以上。此外,所提出的优化算法经过2 0 0次左右的迭代即接近于最终的优化结果,具有较好的收敛性能。关键词:相控阵天线;随机馈相;蚁群算法;天线方向图优化中图分类号:T N 8 2 1+.8文献标志码:AD O I:1 0.3 9 6 9/j.i s s n.1 6 7 1-0 5 7 6.2 0 2 3.0 3.0 0 4收稿日期:2 0 2 2-1 2-1 5基金项目:国家自然科学基金(6 1 5 6 1 0 0 4);江西省教育厅科学技术研究项目(G J J 2 0 1 4 0 8)作者简介:李秋生(1 9 7 6),男,博士,教授,主要从事智能信息处理、目标识别与跟踪等方面的研究。O p t i m i z i n g o f R a n d o m i z i n g P h a s e F e e d i n g S c h e m e f o r P h a s e dA r r a y A n t e n n a B a s e d o n A n t C o l o n y O p t i m i z a t i o n A l g o r i t h mL I Q i u s h e n g1,2,Y ANG Y u x u a n2(1.R e s e a r c h C e n t e r o f I n t e l l i g e n t C o n t r o l E n g i n e e r i n g T e c h n o l o g y,G a n n a n N o r m a lU n i v e r s i t y,G a n z h o u 3 4 1 0 0 0,J i a n g x i,C h i n a;2.S c h o o l o f P h y s i c s a n d E l e c t r o n i cI n f o r m a t i o n,G a n n a n N o r m a l U n i v e r s i t y,G a n z h o u 3 4 1 0 0 0,J i a n g x i,C h i n a)A b s t r a c t:I n o r d e r t o r e d u c e t h e p e a k s i d e l o b e l e v e l o f p h a s e d a r r a y a n t e n n a a n d i m p r o v e t h e a n t i-j a mm i n g p e r f o r m a n c e o f p h a s e d a r r a y a n t e n n a,T a y l o r a m p l i t u d e w e i g h t e d l i n e a r a r r a y w a s t a k e n a s a n e x a m p l e t o a n a l y z e t h e m a t h e m a t i c a l m o d e l o f r a n d o m i z i n g p h a s e f e e d i n g,a n d a n o p t i m i z a t i o n a l g o r i t h m f o r r a n d o m i z i n g p h a s e f e e d i n g o f p h a s e d a r r a y a n t e n n a w a s p r o p o s e d b a s e d o n a n t c o l o n y o p t i m i z a t i o n(A C O)a l g o r i t h m,a n d s i m u l a t i o n e x p e r i m e n t s w e r e c a r r i e d o u t.S i m u l a t i o n r e s u l t s s h o w t h a t t h e p e a k s i d e l o b e l e v e l o f t h e o p t i m i z e d p h a s e d a r r a y a n t e n n a i s r e d u c e d b y m o r e t h a n 8 d B c o m p a r e d w i t h t h a t b e f o r e o p t i m i z a t i o n.I n a d d i t i o n,t h e p r o p o s e d o p t i m i z a t i o n a l g o r i t h m i s c l o s e t o t h e f i n a l o p t i m i z a t i o n r e s u l t a f t e r a b o u t 2 0 0 i t e r a t i o n s,a n d h a s g o o d c o n v e r g e n c e p e r f o r m a n c e.制 导 与 引 信 第4 4卷K e y w o r d s:p h a s e d a r r a y a n t e n n a;r a n d o m i z i n g p h a s e f e e d i n g;a n t c o l o n y o p t i m i z a t i o n a l g o r i t h m;a n t e n n a p a t t e r n o p t i m i z a t i o n0 引言 相控阵天线具有波束扫描速度较快、多波束或自适应波束形成简便等突出优点,在雷达和通信等领域得到了广泛的应用。相控阵天线可通过改变天线阵元的个数、间距、馈电相位等控制天线波束指向或实现天线波束扫描,通过控制各阵元所连接的数字移相器来达到改变馈电相位的目的。在相控阵天线中,移相器和辐射单元是最重要的两个元件。移相器分为连续和数字两类,连续移相器的移相值允许在02 内连续取值,而数字移相器的移相值只能离散取值。在相控阵天线中,移相值需取连续值,但考虑到移相器的体积、质量、成本、移相速度、功率容量和稳定性等因素,通常采用数字移相器。另外,数字移相器还可以通过计算机进行控制,故其被广泛应用于相控阵天线中1-2。然而,使用数字移相器会不可避免地产生量化相位误差,导致天线产生寄生副瓣,降低天线的抗干扰能力和波束指向精度。由于数字移相器的移相取值受限,仅能取最小移相值的整数倍,即2/2,2/22,2/2P,其中P表示数字移相器的位数,通常对所需移相值作舍尾或进位处理,所得相位即为量化相位。使用数字移相器按量化相位进行移相,会产生量化相位误差。如果在此过程中,所有移相值均选择舍尾法或均选择进位法进行量化,则在天线方向图的某些角度上会产生高电平副瓣,称为量化副瓣或寄生副瓣。寄生副瓣将导致天线的抗干扰性能下降,并产生预定方向的波束指向角偏离现象3。通常采用随机馈相法来降低数字移相器所产生的高电平副瓣和波束指向角偏离的影响。所谓随机馈相法就是按给定的概率函数对移相值作进位或舍尾处理,使其变为最小移相值的整数倍4。文献5 对随机馈相进行了详细的理论分析,并指出,随机馈相引起的天线波束副瓣电平的均值的分贝数约为-1 2P。然而,对于随机馈相方案的优化,到目前为止仍没有从理论上提出一个最佳的准则。为此,本文基于蚁群优化算法提出一种相控阵天线随机馈相的优化方案,并对其进行仿真实验验证。1 随机馈相的数学模型 为方便分析,以泰勒(T a y l o r)幅度加权的线阵为例,线阵由N个间距为d的天线阵元组成,如图1所示。相控阵雷达发射机的输出信号被功率分配器分配到N个天线阵元中,每个天线阵元含有一个可变数字移相器,相邻天线阵元移相器的移相值相差,即从左到右的数字移相器移相值分别为0,(N-1)。为天线波束扫描角度。图1 线阵列天线示意图波程差引起的相邻天线阵元的相位差3=2 ds i n/(1)式中:为天线的波长。天线阵远场方向图的电场强度3E=Nn=1Ine x pj2(n-1)d(s i n-s i nB)/+n(2)式中:In表示第n个天线阵元的激励幅度;B表示波束的最大值指向角;n表示第n个天线阵元的精确移相值和实际移相值的差值,即移相误差。采用进位法和舍尾法时,n可以表示为n=-an=-(n-n),进位法 bn=n-n,舍尾法(3)42第3期李秋生,等:基于蚁群优化算法的相控阵天线随机馈相方案的优化式中:-an,bn分别表示采用进位法和舍尾法时第n个天线阵元的移相误差;n表示第n个天线阵元的精确移相值;n和n分别表示采用进位法和舍尾法时第n个天线阵元的实际移相值。移相误差具体含义可以用图2说明。图中,为数字移相器的最小相移值,x轴为天线阵元排布的坐标轴。图2 相位误差图随机馈相的本质就是按一定的概率函数选取天线阵元的移相误差-an或bn。第n个天线阵元的精确移相值n=m+r(4)式中:m=n/,其中 表示取整运算;r为不能被最小移相值整除的小数部分,0r1。给定馈相门限c,0 c 0.5,并对每个天线阵元移相的取值给定一个概率。选取的c不同,随机馈相的方式不同。随机馈相包括完全随机馈相(c=0)、部分随机馈相(0c0.5)、四舍五入法随机馈相(c=0.5)等馈相方式3。则采用随机馈相方式移相时,n可表示为n=-an,pn=1-an,pn1bn,qn=1-pnbn,qn=1 1-cr1,cr1-c,0rc (5)式中:pn,qn分别表示采用进位法和舍尾法时第n个天线阵元的移相取值概率。所谓副瓣电平L,即副瓣电场强度ES L与主瓣波束的电场强度最大值E0之比,其计算方法可表示为L=2 0 l g(ES L/E0)(6)选择完全随机馈相方式进行线阵天线的方向图仿真。仿真参数如下:阵元选取全向天线;阵元电流采用T a y l o r加权;阵元数N取2 8;相邻阵元间的距离d取波长的一半,即/2;波束的最大值指向角B取1 0;数字移相器的位数P取4。则当选取pn=qn=0.5时,优化前完全随机馈相天线仿真方向图如图3所示。可以看出,除俯仰角为1 0.0 附近的主瓣外,在俯仰角为6 9.8 处还产生了峰值电平为-1 9.8 0 4 d B的副瓣,这势必影响相控阵天线的副瓣性能。故为了满足设计要求、提高相控阵天线抗干扰能力,必须对随机馈相方案进行优化。图3 优化前完全随机馈相天线仿真方向图2 蚁群优化算法的基本原理 蚁群优化(a n t c o l o n y o p t i m i z a t i o n,A C O)算法是一种概率型智能优化算法,主要用来寻找最优路径。1 9 9 2年,意大利学者M a r c o D o r i g o发现蚂蚁在寻找食物过程中总是可以找到从起点到终点的最短路径,并基于此提出了A C O算法。经观察发现,在寻找食物的过程中,蚂蚁会持续不断地在走过的路径上留下信息素。起初,各条路径上没有信息素,蚂蚁可任意选择其中一条路径。随后,各路径的信息素浓度随着时间的推移逐渐降低,受信息素浓度的影响,蚂蚁选择信息素浓度高的路径的概率更大。如果蚂蚁移动的速度相同,那么在相同时间内,走较短路径的次数会比走较长路径的多,所以在较短路径上留下的信息素会高于在较长路径(即信息素浓度和路程成反比)上留下的,则该较短路径被之后的蚂蚁选中的可能性也就更大。蚂蚁在较短路径上的往复次数越多,则该路径的信息素浓度就会越高,选择较短路径的蚂蚁就会越多,这是一种正反馈机制。随着时间推移,各路径信息素浓度差距会越发明显,故蚂蚁可以找到最优路径。下面以A C O算法的经典问题 旅行商问题(t r a v e l i n g s a l e s m a n p r o b l e m,T S P)为例来阐52 制 导 与 引 信 第4 4卷述A C O算法的基本原理。首先,设最初城市i和城市j之间的信息素浓度i j相同,i j=e,e为很小的常数。K只蚂蚁随机分布在C座城市中,将这C座城市的编号存放到记录蚂蚁移动路线的禁忌表(t a b u)的第一列,对应城市即为蚂蚁在本次迭代中到达的第一座城市。每只蚂蚁根据转移概率公式选择下一座要转移到的城市,t时刻第k只蚂蚁从城市i转移到城市j的概率pki j(t)6可以表示为pki j(t)=(i j(t)(i j(t)sJk(i s(t)(i s(t),jJk(t)0,其他(7)式中:,分别为信息素启发因子和期望因子;i j为城市i和城市j之间距离的倒数;Jk为第k只蚂蚁还没有走过的城市的集合;i s为城市i和城市s之间的路径上的信息素浓度;i s为城市i和城市s之间距离的倒数。当C座城市都加入到 t a b u之后,第k只蚂蚁就结束了一次周游。K只蚂蚁都结束一次周游后,i j的更新公式6为i j(t+C)=(1-)i j(t)+i j,01(8)其中i j=Kk=1ki j(9)式中:为信息素蒸发系数;i j表示本次迭代中所有蚂蚁在城市i和城市j之间的路径上留下的信息素的增量;ki j表示本次迭代中第k只蚂蚁在城市i和城市j之间的路径上留下的的信息素的增量。M a r c o D o r i g o提出了蚁周(A n t-C y c l e)、蚁量(A n t-Q u a n t i t y)及蚁密(A n t-D e s t i n y)模型6,各模型的差别主要在于ki j的定义。A n t-C y c l e模型的ki j表达式7为ki j=Q/Lk,第k只蚂蚁在此次周游中经过城市i和城市j0,其他 (1 0)式中:Q为正数;Lk表示第k只蚂蚁此次周游的路程。A n t-Q u a n t i t y模型的ki j表达式为ki j=Q/di j,第k只蚂蚁在t和t+1时刻经过城市i和城市j0,其他 (1 1)式中:di j表示城市i和城市j之间的距离。A n t-D e s t i n y模型的ki j表达式为ki j=Q,第k只蚂蚁在t和t+1时刻经过城市i和j0,其他 (1 2)A n t-C y c l e模型在K只蚂蚁均完成周游后,才对路径上的i j进行更新。但A n t-Q u a n t i t y模型和A n t-D e s t i n y模型在蚂蚁每完成一步移动后,就立刻对该路径上的i j进行更新。3 基于蚁群优化算法的随机馈相方案的优化 设有N个天线阵元,K只蚂蚁需分别依次走过每个天线阵元,最终回到起点。相邻天线阵元间有两条路径,分别对应各个天线阵元实际移相值用进位法或舍尾法处理,蚂蚁从第i个天线阵元走到第i+1个天线阵元可选取其中任意一条路径。3.1 信息素更新 采用蚁群优化算法对相控阵随机馈相进行优化,这里i j表示第i个天线阵元到第i+1个天线阵元的路径j上的信息素浓度,i j为本次迭代中所有蚂蚁在第i个和i+1个天线阵元间路径j上留下的信息素增量之和,ki j为第k只蚂蚁在第i个和i+1个天线阵元间路径j上留下的信息素增量,j可取1或2(j=1表示采用进位法,j=2表示采用舍尾法)。该优化方案的信息素浓度迭代公式为i j(t+N)=(1-)i j(t)+i j,00为常数;Sk为第k只蚂蚁遍历N个天线阵元后所得馈相方案的峰值副瓣电平。3.2 状态转移 在相控阵随机馈相方案中,设Pi(Pi 0.5)为第i个和i+1个天线阵元间信息素浓度较高的路径上的信息素浓度占两路径总信息素浓度的比例。利用Pi 来确定第k只蚂蚁下一次的路径,即第i个天线阵元实际移相的取值方式(进位法或舍尾法)。若产生的随机数大于Pi(小概率事件),那么选取信息素浓度较低的路径;反之,则选取信息素浓度较高的路径。Pi的计算公式为Pi=m a xki1ki1+ki2,ki2ki1+ki2 (1 5)式中:m a x()为取最大值函数;ki1,ki2分别表示第k只蚂蚁在第i个天线阵元到第i+1个天线阵元的两条路径上的信息素浓度。3.3 优化算法 采用A C O算法对相控阵天线随机馈相方案进行优化的流程如图4所示。优化算法步骤阐述如下:a)变量初始化,对A C O算法及天线阵列相关的参数进行初始化,如波长、阵元间距d、阵元数N、波束的最大值指向角B、数字移相器位数P、蚁群中蚂蚁的数量K、信息素蒸发 系数等;b)初始路径设置,随机设置K只蚂蚁的初始路径,每只蚂蚁在相邻天线阵元间按Pi=0.5的比例随机选取一条路线,完成一次遍历,并绘制随机馈相的天线方向图,计算其峰值副瓣电平;c)信息表更新,记录本次遍历的最优解,以及最优解对应的各个天线阵元的相移值,并更新各路径上的信息素浓度;d)天线阵元遍历,每只蚂蚁按状态转移公式选取路径,遍历所有天线阵元,并计算相应的峰值副瓣电平;e)迭代判决,判断迭代次数是否小于设定的最大值,若达到最大迭代次数,则终止计算,输出优化后的结果,反之,重复步骤c)和步骤d)。图4 基于A C O的随机馈相方案的优化流程图4 算法仿真与分析 以上述T a y l o r幅度加权的线阵为例,进行优化算法仿真分析。设置仿真参数K=2 0,=0.1,=0.0 5。优化后的随机馈相相控阵天线仿真方向图如图5所示。图5 优化后的随机馈相相控阵天线仿真方向图72 制 导 与 引 信 第4 4卷由图5可以看出,优化后天线方向图的峰值副瓣电平有了较大程度的下降,约为-2 8.2 6 d B,相比于图3优化前的-1 9.8 0 d B,降低了8 d B以上,优化效果相当理想。图6给出了优化过程中群体最优个体对应的相控阵天线方向图的峰值副瓣电平变化曲线。图6 最优个体对应天线峰值副瓣电平变化曲线由图6可知,随着迭代次数的增加,最优个体所对应的天线峰值副瓣电平呈阶梯式下降趋势。曲线按迭代次数大体可以分成三段:第一段为第(11 0 0)次,第二段为第(1 0 12 5 0)次,第三段为第(2 5 1 5 0 0)次。在第一段,峰值副瓣电平快速下降,并逐渐趋近于-2 7.0 0 d B;在第二段,峰值副瓣电平仍在下降,但下降速度变缓,经过这个阶段的迭代优化,峰值副瓣电平约为-2 7.5 0 d B,已经趋近于最终的优化结果;在第三段,峰值副瓣电平在第二段的基础上又降低了大约0.3 5 d B,曲线基本趋于平缓,没有大的变化。从图6还可以看出,实际算法在第2 7 6次迭代后即已收敛,相控阵天线峰值副瓣电平已达到最终优化结果。可知,算法具有较好的收敛性能,在满足设计要求的前提下,可 减 少 迭 代 次 数,降 低 优 化 算 法 的 运算量。5 结束语 本文以T a y l o r幅度加权的线阵为例,对相控 阵天线随机馈相的数学模型进行了理论分析,提出了一种基于蚁群优化算法的相控阵天线随机馈相的优化方法,并进行了仿真验证。仿真结果表明,利用蚁群优化算法可以对相控阵天线随机馈相的最优解进行有效搜索,优化后相控阵天线的峰值副瓣电平显著降低。此外,该算法具有较好的收敛性能,在满足设计要求的前提下,可以减少迭代次数,降低优化算法的运算量。因此,采用该算法对相控阵天线的随机馈相方案进行优化设计,可以有效提高天线的设计质量。同时,需要指出的是,本文提出的优化算法收敛速度还不够理想,还需进一步对该算法进行优化。参考文献1 张胜辉,余鹏程,杨鹏,等.随机馈相法对相控阵天线增益 和指 向精 度 的影 响J.航 天 器 工 程,2 0 1 4,2 3(3):6 7-7 1.2 陈枭煜,孙武.随机馈相方案的仿真对比及阈值参数C的优化分析J.载人航天,2 0 1 8,2 4(2):2 6 7-2 7 2.3 李秋生.部分随机馈相对相控阵天线波瓣性能的影响分析J.赣南师范学院学报,2 0 0 6,2 7(6):7-1 0.4 张光义,赵玉洁.相控阵雷达技术M.北京:电子工业出版社,2 0 0 6.5 郭燕昌,钱继曾,黄富雄,等.相控阵与频率扫描天线原理M.北京:国防工业出版社,1 9 7 8.6 刘加奇,王泰华,董征.基于改进蚁群算法的移动机器 人 路 径 规 划 J.传 感 器 与 微 系 统,2 0 2 2,4 1(5):1 4 0-1 4 4.7 曾胜,戴贤君,胡徐胜,等.基于蚁群算法对调度车辆进行 路径 最优 化 设计 J.自 动 化 与 仪 表,2 0 2 2,3 7(4):8 9-9 3,9 8.8 李秋生.基于遗传算法的相控阵随机馈相方案的优化J.系统 工程 与 电子 技术,2 0 0 6,2 8(6):8 6 1-8 6 3.82