CN116882147A - 基于混合麻雀搜索算法的阵列天线合成方法 - Google Patents
基于混合麻雀搜索算法的阵列天线合成方法 Download PDFInfo
- Publication number
- CN116882147A CN116882147A CN202310799861.XA CN202310799861A CN116882147A CN 116882147 A CN116882147 A CN 116882147A CN 202310799861 A CN202310799861 A CN 202310799861A CN 116882147 A CN116882147 A CN 116882147A
- Authority
- CN
- China
- Prior art keywords
- array
- value
- individuals
- algorithm
- individual
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 241000287127 Passeridae Species 0.000 title claims abstract description 51
- 238000010845 search algorithm Methods 0.000 title claims abstract description 27
- 238000001308 synthesis method Methods 0.000 title claims abstract description 7
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 48
- 230000003044 adaptive effect Effects 0.000 claims abstract description 13
- 230000008859 change Effects 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 41
- 238000000034 method Methods 0.000 claims description 21
- 230000035772 mutation Effects 0.000 claims description 16
- 238000013461 design Methods 0.000 claims description 7
- 238000009826 distribution Methods 0.000 claims description 7
- 238000009827 uniform distribution Methods 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 238000003786 synthesis reaction Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000012163 sequencing technique Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 2
- 230000002194 synthesizing effect Effects 0.000 claims 1
- 238000005457 optimization Methods 0.000 description 12
- 238000003491 array Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 230000007547 defect Effects 0.000 description 4
- 238000013507 mapping Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 241000255581 Drosophila <fruit fly, genus> Species 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000019637 foraging behavior Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Evolutionary Computation (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Software Systems (AREA)
- Computing Systems (AREA)
- Artificial Intelligence (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computational Linguistics (AREA)
- Computer Hardware Design (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Geometry (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了一种基于混合麻雀搜索算法的阵列天线合成方法,属于射频天线技术领域。本发明根据阵列类型及预期指标,建立合适的阵列模型和适应度函数,利用混合麻雀搜索算法对模型的参数进行优化,获得最优的解决方案。本发明的混合麻雀搜索算法中混合的策略包括自适应搜索步长策略和差分变异策略,前者引入自适应参数控制搜索步长动态变化,以提升麻雀搜索算法的搜索范围、精度并加快算法收敛;后者使用差分进化算法的变异操作对种群进行扰动,以增加种群多样性,避免算法陷入局部最优,进而提高算法在高维问题中的搜索效率。
Description
技术领域
本发明属于射频天线技术领域,具体提出了一种基于混合策略的麻雀搜索算法,应用于天线阵列方向图合成的技术领域,可实现稀疏阵列和稀布阵列的设计。
背景技术
随着社会的发展,无线通信技术有了更广泛的应用场景,比如雷达、导航、卫星等,传统的单天线已经无法满足实际应用场景中对方向图的要求,还存在着衰减大、传输效率低等问题。因此由一定数量天线单元组成的阵列天线逐渐有了广泛的应用。根据电磁波在空间相互干涉的原理,阵列天线可以通过对天线单元的电流、相位以及间距等的控制,产生各种各样的波束特性,灵活满足不同应用场景。
以阵元间距作为优化变量,可将天线阵列分为均匀间距阵列和不等间距阵列。相对于均匀间距阵列,不等间距阵列使用较少的阵元数目获得了较窄的波束和较低的峰值旁瓣电平(PSLL),从而具有成本低、体积小、重量轻等优点,被广泛应用在卫星通信、雷达和远程无线传感系统等领域中。不等间距阵列通常被分为两类:具有随机间距的稀布阵列和通过对初始等间距数组的某些元素进行选择性归零得到的稀疏阵列。其中,稀布阵列具有更大的自由度去优化阵元位置使PSLL更低,近年来受到了更多的关注。然而,阵列响应是阵元位置的复指数函数。考虑到互耦效应和栅瓣等因素,还要对阵列孔径、阵元间距、阵元数目、波束宽度等进行约束,因此不均匀阵列的合成是一个较难的多约束、多维、非线性优化问题。传统的确定性解析算法,比如矩阵铅笔法和贝叶斯压缩感知方法等,是通过对已知的参考方向图进行匹配来重建阵元位置的,并且可能产生不可行解,应用范围较小。而随机搜索算法是使用智能优化算法,比如遗传算法、粒子群算法等,对阵列模型进行迭代优化,并且在迭代过程中通过约束条件不断维护解的可行性,最终得到满足设计目标的阵列布局。但有着陷于局部最优的风险,而且面对多维问题时计算时间较长。
麻雀搜索算法(SSA)是2020年提出的一种群体智能优化算法,模仿自然界麻雀的觅食行为和反捕食行为,将麻雀群体分为发现者、加入者和预警者三个角色。发现者负责向外探索,加入者会在最优个体附近进行探索,预警者负责突变以增加种群多样性,三个角色相互配合,使得麻雀算法具有调整参数少,搜索效率高、通用性强等优点。自提出以来,该算法已被广泛应用于各个领域,并被证明可以有效地优化高难度多维问题。有学者将SSA与A-star算法结合,应用于车辆的全局路径规划;也有学者提出了一种改进麻雀搜索算法(ISSA)解决光伏微电网系统中部分遮阳下的功率失配损失问题。但SSA因为发现者的搜索步长变化单一、预警者的突变方向固定等缺点,所以也存在着收敛精度低、易陷入局部最优等常见问题。
发明内容
针对现有技术的缺陷,本发明的目的在于提出一种基于混合麻雀搜索算法(MSSA)的阵列天线合成方法。该方法使用自适应搜索步长策略和差分变异策略分别改进了SSA的发现者和预警者模型,自适应搜索步长通过控制迭代前后期及搜索能力,提升了算法的搜索范围和精度,而差分变异增加了种群多样性,避免算法陷入局部最优,使其在阵列设计问题中具有更高的搜索效率和稳定性。
为了实现上述目的,本发明采用的技术方案如下:
基于混合麻雀搜索算法的阵列天线合成方法,包括以下步骤:
步骤一,建立阵列天线模型,采用稀疏阵列或稀布阵列;阵列的初始化参数包括阵元数目、阵列口径、阵元位置,阵列的约束条件包括最小阵元间距、是否固定阵列孔径、阵列是否对称;
步骤二,确定混合麻雀搜索算法的参数;初始化参数包括麻雀种群大小NP、最大迭代次数tmax、优化问题维度d、第t次迭代中第i个麻雀个体的第j维位置信息发现者比例PD、预警者比例SD、安全值ST、变异算子F;
步骤三,生成初始种群;根据选择的阵列模型,随机生成初始种群,种群中每个个体都是一个矩阵形式的完整阵列布局;如果阵列模型是稀疏阵列,则每个个体为由0和1组成的二维矩阵,矩阵的维度与阵列规模保持一致;如果阵列模型是稀布阵列,则每个个体由X、Y轴的相对间距矩阵及阵元的选择矩阵联合组成;
步骤四,设计适应度函数,计算初始适应度值;对步骤三中的每个个体根据对应模型,转化为具体的平面坐标并计算适应度值;
步骤五,排序;根据适应度值从小到大的顺序对个体进行排序,记录排序第一和末尾的个体,作为全局最优个体Xbest和全局最差个体Xworst;按照发现者比例PD将种群分为发现者和加入者;
步骤六,发现者更新;对发现者个体使用引入自适应搜索步长策略的发现者更新公式进行更新,生成新个体,并计算适应度值;
步骤七,加入者更新;对加入者个体使用加入者更新公式进行更新,生成新个体,并计算适应度值;
步骤八,预警者更新;按照预警者比例SD在原种群中随机挑选部分个体记为预警者,对预警者使用差分变异操作进行更新,生成新个体,并计算适应度值;
步骤九,约束条件控制;对于步骤六、步骤七和步骤八生成的新个体,如果个体值中某个数值超出设定的变量边界值,则将越界的数值纠正至合理数值,然后再计算适应度值;
步骤十,选择;依次对比每个个体新旧两代的适应度值,取较优的那一代个体作为最终的新一代个体,得到新一代种群;
步骤十一,迭代;重复执行步骤五至步骤十,直至满足波束宽度或峰值旁瓣电平的设计要求,或者达到最大迭代次数,此时的Xbest就是该阵列天线合成中的最优阵列布局。
进一步地,步骤九中,将越界的数值纠正至合理数值的方式如下:
Xnew=rand(1)×Xold
其中,Xold表示旧值,Xnew表示纠正后的新值,rand(1)是在0~1之间符合均匀分布的随机数。
进一步地,步骤四中的适应度函数定义为和/>平面内峰值旁瓣电平之和:
其中,S1和S2分别为两个平面的旁瓣角度区域,AF(·,·)为阵列因子,max表示取最大值,M、BWc、BWd、μ分别代表权重系数、计算波束宽度、预设波束宽度、允许误差。
进一步地,步骤六中的发现者更新公式为:
其中,Q是服从标准正态分布的随机数;L表示一个1×d的单位矩阵;ST∈[0.5,1]是设定的安全值;R是预警值,为0~1的随机数;w是自适应参数,控制搜索步长的变化,其计算公式如下:
取wmax=10,wmin=0.1,rand(1)是在0~1之间符合均匀分布的随机数。
进一步地,步骤七中的加入者更新公式为:
其中,是目前发现者所占据的最优位置,/>是当前全局最差的位置;A表示一个1×d的矩阵,每个元素随机赋值为1或-1,并且A+=AT(AAT)-1;L表示一个1×d的单位矩阵。
进一步地,步骤八中的预警者更新公式为:
其中,F是变异算子,rand(1)是在0~1之间符合均匀分布的随机数;r1、r2、r3是种群内互不相同的三个个体。
本发明较现有技术具有如下优点:
1、本发明的麻雀搜索算法混合了自适应搜索步长策略。引入自适应参数w,由个体序号i和迭代次数t联合控制发现者向外搜索的步长,使得发现者会在迭代初期、高个体序号(较差个体)时进行宽泛搜索;高迭代次数、低个体序号(较优个体)时进行精确搜索,兼顾了搜索范围和精度,加快了算法的收敛,提高了算法的搜索效率。使用正态分布的随机数Q替换0~1的随机数α,使得搜索步长可以大于1,进一步扩大了搜索范围。
2、本发明的麻雀搜索算法混合了差分进化算法的变异操作。采用“current-to-rand/1”变异操作,用三个不同的个体对现有个体进行干扰,以产生新个体,避免算法过早收敛。相对于原预警者在全局最优和最差个体间移动进行更新的方式,变异操作克服了新个体大多仍处在种群内部的缺陷,产生了更强、更稳定的扰动能力,从而提高种群的多样性,进一步提升算法的搜索能力和稳定性。
3、本发明将混合麻雀搜索算法应用于稀疏阵列和稀布阵列设计中,在约束条件下,获得满足预期波束宽度和最大旁瓣电平的阵列设计结果。具有更快的收敛速度和搜索能力,有效地解决了智能算法稳定性较差的弊端,通用性较强。
附图说明
图1为本发明混合麻雀搜索算法流程图。
图2为本发明的仿真实例一稀布阵列天线单元分布图。
图3为本发明的仿真实例一稀布阵列天线方向图。
图4为本发明的仿真实例二稀布阵列天线单元分布图。
图5为本发明的仿真实例二稀布阵列天线方向图。
具体实施方式
以下结合仿真实例与附图,对本发明进行详细描述。该仿真实例及其说明仅用于解释本发明,并不作为对本发明的限定。
一种基于混合麻雀搜索算法的阵列天线合成方法,包括以下步骤:
(1)建立天线阵列模型。可选择的天线阵列模型为稀疏阵列与稀布阵列。阵列的初始化参数包括阵元数目、阵列口径、阵元位置等。阵列的约束条件包括最小阵元间距、是否固定阵列孔径、阵列是否对称等。
假如阵元数目为N,对于稀疏阵列模型,将阵列按照最小阵列间距划分为网格,每个网格点代表一个阵元。则初始种群为数字0与1组成的阵列,其中0表示该位置天线关闭,1表示该位置天线开启,需要保证1的个数正好是N;对于稀布阵列,将阵元随机分布在阵列中,每个阵元的坐标都用一对实数(xn,yn)表示,但是需要使用非对称映射方法约束阵元间距。首先设置P和Q来定义阵元矩阵的行数和列数,需要保证P×Q≥N,使用选择矩阵W∈RP×Q来打开N个阵元,关闭P×Q-N个阵元。再设置矩阵A和B分别表示每个阵元与上一个阵元的X轴相对距离和Y轴相对距离。最后使用映射函数T1和T2将A和B矩阵映射成坐标矩阵X=T1(A)和Y=T2(B)。
对阵元数目为4N并且关于原点对称的平面阵列进行等幅度激励(In=1)和零相位激励则该阵列模型的阵列因子可表示为:
其中,θ和分别代表俯仰角和方位角;k=2π/λ,是空间波数,λ是波长。
(2)确定算法参数。初始化参数包括麻雀种群大小NP、最大迭代次数tmin、优化问题维度为d、第t次迭代中第i个麻雀个体的第j维位置信息发现者比例PD、预警者比例SD、安全值ST、变异算子F等。
(3)生成初始种群。根据选择的阵列模型,随机生成初始种群,种群中每个个体都是一个矩阵形式的完整阵列布局。如果是稀疏阵列,则每个个体为0和1组成的选择矩阵,维度与阵列规模保持一致。如果是稀布阵列,则每个个体由X、Y轴的相对间距矩阵及阵元的选择矩阵联合组成。
(4)计算初始适应度值。对步骤(3)中的每个个体根据对应模型,转化为具体的平面坐标并计算适应度值。为了得到满足波束宽度条件的最低的最大旁瓣电平,适应度函数定义为和/>平面内峰值旁瓣电平之和:
其中,S1和S2分别为两个平面的旁瓣角度区域。公式后部分是对两个平面波束宽度的约束,M、BWc、BWd、μ分别代表权重系数、计算波束宽度、预设波束宽度、允许误差。可以通过较大的M,将波束宽度限制在预设区间中。
(5)排序。根据适应度值连同个体进行从小到大排序,记录排序第一和末尾的个体,记为全局最优个体Xbest和全局最差个体Xworst。种群的前1:NP×PD个体记为发现者,剩余个体NP×PD+1:NP记为加入者。
(6)发现者更新。对发现者使用如下公式进行更新,并计算适应度值:
其中,Q是服从标准正态分布的随机数;L表示一个1×d的单位矩阵;ST∈[0.5,1]是设定的安全值;R是预警值,为0~1的随机数;w是自适应参数,控制搜索步长的变化,其计算公式如下:
取wmax=10,wmin=0.1,使w随迭代次数的增加而自适应增大。
当随机生成的预警值小于安全值时,发现者会以自己为中心向四周探索,搜索步长在1上下波动,变化趋势由个体序号i和自适应参数w联合控制。使得发现者在迭代初期、高个体序号(较差个体)时进行宽泛搜索;高迭代次数、低个体序号(较优个体)时进行精确搜索。兼顾了搜索范围和搜索精度,提高了算法的搜索效率。反之则说明当前不安全,麻雀个体会向四周随机跳跃,增加了探索方式。
(7)加入者更新。对加入者使用如下公式进行更新,并计算适应度值:
其中,是目前发现者所占据的最优位置,/>是当前全局最差的位置。A表示一个1×d的矩阵,每个元素随机赋值为1或-1,并且A+=AT(AAT)-1。当i≤NP/2时,加入者会飞到发现者中的最优个体附近觅食;反之表明较差的加入者没有获得食物,此时需要飞往远离最差个体的地方觅食。
(8)预警者更新。在原种群中随机挑选NP×SD个体记为预警者,对预警者使用如下公式进行更新,并计算适应度值:
其中,F是变异算子;r1、r2、r3是种群内互不相同的三个个体。“current-to-rand/1”变异操作以三个个体对当前个体进行干扰,能够提高种群的多样性,进一步提升算法的搜索能力和稳定性。
(9)约束条件控制。对于步骤(6)、步骤(7)和步骤(8)生成的新个体,如果个体值中某个数值超出设定的变量边界值,则利用边界控制算法,将越界的数值纠正至合理数值,然后再计算适应度值。所选择的实际的边界控制算法是利用随机值将该数值拉回变量区间中,即:
Xnew=rand(1)×Xold
(10)选择。对比新旧个体的适应度值,取较优的个体组成新一代种群。
此处,将当前准备更新的种群被称为旧种群,也就是上一代种群,经过更新后的种群称为新种群,最后会一一对比旧种群和新种群中相同序号的两个个体,保留适应度值较低(因为本方法的目标是最低的峰值旁瓣电平,所以是保留适应度值较低的个体,换种说法也是更优秀的个体)的那个个体,全部比较完毕后,保留的全部个体组成真正的下一代种群,用它们参与下一次迭代,下一次迭代中它们又被当做为旧种群。
(11)迭代。重复执行步骤(5)至步骤(110),直至满足设定的指标,或者达到最大迭代次数,此时的Xbest就是最终的最优解。
该方法采用了混合麻雀搜索算法,首先设定算法参数麻雀种群大小NP、最大迭代次数tmax、发现者比例PD、预警者比例SD、安全值ST、变异算子F,再根据阵列的参数如阵元数目、阵列口径、阵元最小间距等确定优化问题维度d,最后根据选择的阵列天线类型随机生成麻雀个体组成初代麻雀种群,其中,稀疏阵列天线的麻雀个体是0和1组成的二维矩阵,稀布阵列天线的麻雀个体是由X、Y轴的相对间距矩阵及阵元的选择矩阵联合组成的矩阵。
结合图1,以20×10平面稀布阵列的合成为示例,对该方法进行具体阐述。具体有两种情况:
实例一:最小间距0.5λ、固定阵列孔径。
实例二:最小间距0.2λ。
下面,结合实例一进行说明。该方法包含以下步骤:
步骤一:建立天线阵列模型。选择平面稀布阵列模型。初始化参数包括天线阵列的阵元数目4N=108(1/4区域有27个阵元),阵列孔径2L=9.5λ,2H=4.5λ。阵列的约束条件为最小阵元间距dc=0.5λ、固定阵列孔径(即阵列的四个角上固定有阵元)、阵列上下左右对称。对该阵列进行等幅度激励(In=1)和零相位激励则阵列因子可表示为:
其中,θ和分别代表俯仰角和方位角;k=2π/λ,是空间波数,λ是波长。
步骤二:确定算法参数。初始化参数包括麻雀种群大小NP=200、最大迭代次数tmin=3000、优化问题维度为d=10×20、第t次迭代中第i个麻雀个体的第j维位置信息发现者比例PD=0.2、预警者比例SD=0.2、安全值ST=0.8、变异算子F=0.8(侧重提高种群多样性)。
步骤三:生成初始种群。由于1/4区域有27个阵元,设置阵元形式为4×9,即X轴上9列阵元,Y轴上4行阵元,由阵列孔径和阵元最小间距可以计算出X轴、Y轴的剩余最大可分配空间。定义优化变量矩阵A、B分别代表X、Y轴的位置信息,其中A矩阵的维度为4×10,每个元素的范围都是[0,1],每一行的第10列代表当前行的实际可分配空间占最大可分配空间的比例,而每一行的前9列分别代表相邻阵元的相对间距占实际可分配空间的比例。同理B矩阵的维度为5×9,每个元素的范围都是[0,1],每一列的第5行代表当前列的实际可分配空间占最大可分配空间的比例,而每一列的前4行分别代表相邻阵元的相对间距占实际可分配空间的比例。最后定义4×9矩阵C,由0和1组成,总共27个1代表对应的27个阵元是开启状态,剩余的9个0代表对应的9个阵元是关闭状态。由A矩阵可以计算出相邻阵元的X轴实际间距,累加得到每个阵元的实际X轴坐标。由B矩阵可以计算出相邻阵元的Y轴实际间距,但是阵元Y轴实际坐标是前一行阵元的最大Y轴坐标加上自身的相对间距得到的,以约束阵元间距。由于设置了固定阵列孔径的约束,还需要对最后一个阵元的坐标进行修正。由A、B矩阵得到阵元的实际坐标后,再使用C矩阵就可以获得该个体的实际阵列布局。因此初始化的每个个体都由A、B、C三个矩阵构成,矩阵的每个参数都是0~1的随机数。
步骤四:计算初始适应度值。对步骤(3)中的每个个体转化为具体阵列布局并计算适应度值。为了得到满足波束宽度条件的最低的最大旁瓣电平,适应度函数定义为和/>平面内最大旁瓣电平之和:
其中,S1和S2分别为两个平面的旁瓣角度区域。公式后部分是对两个平面波束宽度的约束,M1=M2=100,BWd1=BWd2=15,μ1=μ2=10分别代表权重系数、计算波束宽度、预设波束宽度、允许误差。可以通过较大的M,将波束宽度限制在预设区间[5,25]之间,由此让算法获得波束宽度在该范围内的最低峰值旁瓣电平。
步骤五:排序。根据适应度值连同个体进行从小到大排序,记录排序第一和末尾的个体,记为全局最优个体Xbest和全局最差个体Xworst。种群的前1:40(NP×PD)个体记为发现者,剩余个体41:200记为加入者。
步骤六:发现者更新。对发现者使用如下公式进行更新,并计算适应度值:
其中,Q是服从标准正态分布的随机数;L表示一个1×d的单位矩阵;ST=0.8是设定的安全值;R是预警值,为0~1的随机数;w是自适应参数,控制搜索步长的变化,其计算公式如下:
取wmax=10,wmin=0.1,使w随迭代次数的增加而自适应增大。
当随机生成的预警值小于安全值时,发现者会以自己为中心向四周探索,反之则说明当前不安全,麻雀个体会向四周随机跳跃。
步骤七:加入者更新。对加入者使用如下公式进行更新,并计算适应度值:
其中,是目前发现者所占据的最优位置,/>是当前全局最差的位置。A表示一个1×d的矩阵,每个元素随机赋值为1或-1,并且A+=AT(AAT)-1。当i≤NP/2时,加入者会飞到发现者中的最优个体附近觅食;反之表明较差的加入者没有获得食物,此时需要飞往远离最差个体的地方觅食。
步骤八:预警者更新。在原种群中随机挑选NP×SD个体记为预警者,对预警者使用如下公式进行更新,并计算适应度值:
其中,F是变异算子,为了侧重提高种群多样性,F取0.8;r1、r2、r3是种群内互不相同的三个个体。
步骤九:约束条件控制。对于步骤六、步骤七和步骤八生成的新个体,如果个体值中某个数值超出设定的变量边界值,则利用边界控制算法,将越界的数值纠正至合理数值。所选择的实际的边界控制算法是利用随机值将该数值拉回变量区间中,即:
Xnew=rand(1)×Xold
步骤十:选择。对比新旧个体是适应度值,取较低适应度值的个体组成新一代种群。
步骤十一:迭代。重复执行步骤五至步骤十,直至满足设定的指标,或者达到最大迭代次数,此时的Xbest就是最终的最优解。
实例一:混合麻雀搜索算法和原麻雀搜索算法、遗传算法、差分进化算法的结果如表1所示,可以看出混合麻雀搜索算法的最佳结果和平均结果都由于其他算法,具有更好的优化效果。
表1不同算法在实例一中结果比较
最佳结果的阵列布局如图2所示、阵元坐标如表2所示。最佳结果的阵列方向图如图3所示。
表2实例一中的阵列布局的阵元坐标
表中,第一、第二个数字分别代表X、Y轴坐标,单位为波长。
实例二:混合麻雀搜索算法取得的最佳结果的阵列布局如图4所示、阵元坐标如表3所示、阵列方向图如图5示。混合麻雀搜索算法取得最佳结果为-89.707dB(-46.170dB在平面、-43.537dB在/>平面),相同实验中,果蝇优化算法、麻雀搜索算法、和差分进化算法取得的最佳结果分别为-72.523dB,-74.76dB和-80.427dB。混合麻雀搜索算法的平均结果为-85.781dB,同样优于上述算法取得的最佳结果,进一步突出了该算法的搜索效果。
表3实例二中的阵列布局的阵元坐标
表中,第一、第二个数字分别代表X、Y轴坐标,单位为波长。
总之,本发明根据阵列类型及预期指标,建立合适的阵列模型和适应度函数,利用混合麻雀搜索算法对模型的参数进行优化,获得最优的解决方案。本发明的混合麻雀搜索算法中混合的策略包括自适应搜索步长策略和差分变异策略,前者引入自适应参数控制搜索步长动态变化,以提升麻雀搜索算法的搜索范围、精度并加快算法收敛;后者使用差分进化算法的变异操作对种群进行扰动,以增加种群多样性,避免算法陷入局部最优,进而提高算法在高维问题中的搜索效率。
Claims (6)
1.基于混合麻雀搜索算法的阵列天线合成方法,其特征在于,包括以下步骤:
步骤一,建立阵列天线模型,采用稀疏阵列或稀布阵列;阵列的初始化参数包括阵元数目、阵列口径、阵元位置,阵列的约束条件包括最小阵元间距、是否固定阵列孔径、阵列是否对称;
步骤二,确定混合麻雀搜索算法的参数;初始化参数包括麻雀种群大小NP、最大迭代次数tmax、优化问题维度d、第t次迭代中第i个麻雀个体的第j维位置信息发现者比例PD、预警者比例SD、安全值ST、变异算子F;
步骤三,生成初始种群;根据选择的阵列模型,随机生成初始种群,种群中每个个体都是一个矩阵形式的完整阵列布局;如果阵列模型是稀疏阵列,则每个个体为由0和1组成的二维矩阵,矩阵的维度与阵列规模保持一致;如果阵列模型是稀布阵列,则每个个体由X、Y轴的相对间距矩阵及阵元的选择矩阵联合组成;
步骤四,设计适应度函数,计算初始适应度值;对步骤三中的每个个体根据对应模型,转化为具体的平面坐标并计算适应度值;
步骤五,排序;根据适应度值从小到大的顺序对个体进行排序,记录排序第一和末尾的个体,作为全局最优个体Xbest和全局最差个体Xworst;按照发现者比例PD将种群分为发现者和加入者;
步骤六,发现者更新;对发现者个体使用引入自适应搜索步长策略的发现者更新公式进行更新,生成新个体,并计算适应度值;
步骤七,加入者更新;对加入者个体使用加入者更新公式进行更新,生成新个体,并计算适应度值;
步骤八,预警者更新;按照预警者比例SD在原种群中随机挑选部分个体记为预警者,对预警者使用差分变异操作进行更新,生成新个体,并计算适应度值;
步骤九,约束条件控制;对于步骤六、步骤七和步骤八生成的新个体,如果个体值中某个数值超出设定的变量边界值,则将越界的数值纠正至合理数值,然后再计算适应度值;
步骤十,选择;依次对比每个个体新旧两代的适应度值,取较优的那一代个体作为最终的新一代个体,得到新一代种群;
步骤十一,迭代;重复执行步骤五至步骤十,直至满足波束宽度或峰值旁瓣电平的设计要求,或者达到最大迭代次数,此时的Xbest就是该阵列天线合成中的最优阵列布局。
2.根据权利要求1所述的基于混合麻雀搜索算法的阵列天线合成方法,其特征在于,步骤九中,将越界的数值纠正至合理数值的方式如下:
Xnew=rand(1)×Xold
其中,Xold表示旧值,Xnew表示纠正后的新值,rand(1)是在0~1之间符合均匀分布的随机数。
3.根据权利要求1所述的基于混合麻雀搜索算法的阵列天线合成方法,其特征在于,步骤四中的适应度函数定义为和/>平面内峰值旁瓣电平之和:
其中,S1和S2分别为两个平面的旁瓣角度区域,AF(·,·)为阵列因子,max表示取最大值,M、BWc、BWd、μ分别代表权重系数、计算波束宽度、预设波束宽度、允许误差。
4.根据权利要求1所述的基于混合麻雀搜索算法的阵列天线合成方法,其特征在于,步骤六中的发现者更新公式为:
其中,Q是服从标准正态分布的随机数;L表示一个1×d的单位矩阵;ST∈[0.5,1]是设定的安全值;R是预警值,为0~1的随机数;w是自适应参数,控制搜索步长的变化,其计算公式如下:
取wmax=10,wmin=0.1,rand(1)是在0~1之间符合均匀分布的随机数。
5.根据权利要求1所述的基于混合麻雀搜索算法的阵列天线合成方法,其特征在于,步骤七中的加入者更新公式为:
其中,是目前发现者所占据的最优位置,/>是当前全局最差的位置;A表示一个1×d的矩阵,每个元素随机赋值为1或-1,并且A+=AT(AAT)-1;L表示一个1×d的单位矩阵。
6.根据权利要求1所述的基于混合麻雀搜索算法的阵列天线合成方法,其特征在于,步骤八中的预警者更新公式为:
其中,F是变异算子,rand(1)是在0~1之间符合均匀分布的随机数;r1、r2、r3是种群内互不相同的三个个体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310799861.XA CN116882147A (zh) | 2023-07-03 | 2023-07-03 | 基于混合麻雀搜索算法的阵列天线合成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310799861.XA CN116882147A (zh) | 2023-07-03 | 2023-07-03 | 基于混合麻雀搜索算法的阵列天线合成方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116882147A true CN116882147A (zh) | 2023-10-13 |
Family
ID=88256081
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310799861.XA Pending CN116882147A (zh) | 2023-07-03 | 2023-07-03 | 基于混合麻雀搜索算法的阵列天线合成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116882147A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117892097A (zh) * | 2024-03-14 | 2024-04-16 | 陕西省杂交油菜研究中心 | 基于大数据的油菜化杀效果评估方法及系统 |
-
2023
- 2023-07-03 CN CN202310799861.XA patent/CN116882147A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117892097A (zh) * | 2024-03-14 | 2024-04-16 | 陕西省杂交油菜研究中心 | 基于大数据的油菜化杀效果评估方法及系统 |
CN117892097B (zh) * | 2024-03-14 | 2024-05-28 | 陕西省杂交油菜研究中心 | 基于大数据的油菜化杀效果评估方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Guney et al. | Amplitude-only pattern nulling of linear antenna arrays with the use of bees algorithm | |
CN107944133B (zh) | 基于多目标量子蜘蛛群演化机制的环形天线阵列稀疏方法 | |
CN111160556B (zh) | 一种基于自适应遗传算法的阵列稀疏优化方法 | |
CN107657098B (zh) | 基于量子鸡群演化机制的环形天线阵列稀疏方法 | |
CN116882147A (zh) | 基于混合麻雀搜索算法的阵列天线合成方法 | |
CN109992845B (zh) | 一种基于matlab和hfss的波束控制和阵列方向图优化方法 | |
CN111432368A (zh) | 一种适用于稀疏锚节点wsn的测距定位方法 | |
CN104020448A (zh) | 等阵元约束的雷达子阵级和波束/差波束形成优化方法 | |
CN111276822B (zh) | 一种天线方向图主瓣可控的天线阵列设计方法 | |
CN109271735B (zh) | 基于量子启发式引力搜索算法的阵列方向图综合方法 | |
CN116882149A (zh) | 一种基于混合差分果蝇优化算法的天线阵列综合方法 | |
CN116205143A (zh) | 一种基于物理信息神经网络实现天线方向图的设计方法 | |
CN114239439B (zh) | 一种基于部落算法的自动化滤波器的设计方法 | |
CN115133291A (zh) | 非规则天线子阵、相控阵天线及相控阵天线的设计方法 | |
CN115374695A (zh) | 基于麻雀搜索算法与阵列加权的稀布阵列天线优化方法 | |
CN114218849A (zh) | 基于深度强化学习的复杂阵列天线智能设计方法 | |
CN114861446A (zh) | 天线阵列稀疏稀布优化方法、装置、设备及可读存储介质 | |
GB2398429A (en) | Partitioning an antenna array | |
CN113067615B (zh) | 一种近场阵列波束控制方法 | |
CN114239380B (zh) | 基于自调整映射法则的矩形稀布阵列优化方法 | |
GB2398428A (en) | Partitioning process for antenna or sensor arrays | |
Zhang et al. | A Fast and Efficient Beamforming Algorithm Imitating Plant Growth Gene for Phased Array Antenna | |
CN110705680B (zh) | 一种分布式阵列的子阵列位置优化方法及系统 | |
CN116632561B (zh) | 一种稀布高增益天线阵列设计方法、存储介质 | |
Zhao et al. | A hybrid algorithm for synthesizing linear sparse arrays |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |