基于免疫遗传算法的希布阵天线的优化设计方法与流程

1.本发明涉及天线阵列技术领域，尤其涉及一种基于免疫遗传算法的希布阵天线的优化设计方法。


背景技术：

2.当前，天线布阵的设计方法有很多种，传统的方法如穷举法、阵元法、动态规划法、智能算法中的遗传算法、粒子群算法、凸优化算法等，都得到了广泛的应用。近几年来，随着大数据技术的发展，由于其简单、通用性强、鲁棒性强以及适合于并行处理等特点，其应用领域主要有函数优化、组合优化、作业调度、智能控制、机器人学、图像处理、人工生命、机器学习等。该方法采用群体搜索技术，将种群表示一组问题解，通过一系列的遗传操作，如选择、交叉和变异，对当前种群进行选择、交叉和变异，从而产生新一代种群，并逐步使种群处于包含近似最佳解的状态。
3.稀布天线阵列的优化布阵技术就是天线阵的辐射方向图的优化。天线阵的辐射方向图组成要素是阵列单元数目、分布形式、单元间距、激励幅度和相位，因此我们通常是通过优化算法来优化这些要素的某些方面以得到最逼近的期望方向图。对稀布天线阵列而言，如何设计一组阵元间距的最优解，使得稀布天线阵列的峰值旁瓣电平在整个可见区最小，是天线布阵技术要解决的主要问题。
4.由于稀布线阵没有稀疏线阵那样具有规则的栅格约束阵元的位置，所以要减小阵元间的互藕又必须约束阵元间距不小于某个值。目前解决的方法主要有两种，第一种方法是约束阵元间距的最小和最大。第二种方法限制最小阵元间距和孔径。阵列孔径决定主瓣的宽度，若给定阵元数之后，方法一布阵结束后才能确定阵列孔径因此相当于不受主瓣宽度的限制。另一种方法已包含在优化目标中对主瓣宽度要求。


技术实现要素：

5.本发明所要解决的技术问题是如何提供一种能够使天线阵列峰值旁瓣电平显著降低，计算量小、收敛速度快，效率高的基于免疫遗传算法的希布阵天线的优化设计方法。
6.为解决上述技术问题，本发明所采取的技术方案是：一种基于免疫遗传算法的希布阵天线的优化设计方法，其特征在于包括如下步骤：
7.随机生成初始群体：希布阵天线中阵元数为n，阵列孔径为l，利用随机法对随机赋值一个约束范围内的实数，组成列向量即为一个个体，考虑到每个个体的适应度计算量大，采用取相邻两个阵元距离的差值为个体作为种群个体，d为阵元间距；
8.抽取疫苗：将公式作为疫苗，进行归一化方向图计算，其中λ为雷达工作波长，d为阵元间距，θ为波束指向角；θ0为方位与天线法线夹角；
9.适应度的计算：假设阵列孔径l固定不变，方向图的函数表达式为
其中d
n
是与第一个阵元的距离；归一化方向图旁瓣电压表达式为将归一化方向图作为适应度函数；
10.若满足上述的条件，结束算法，不满足进行如下遗传操作：
11.遗传操作包括选择、交叉和变异：选择是通过个体的适应度大小决定是否被选择，被选择的概率为交叉或变异引入的新元素阵元很可能使得新个体不满足最小阵元间距约束，因此在个体作经过一定变化后再进行交叉和变异；遗传操作结束之后，加入免疫算子，提高适应度和防止记忆退化现象的出现；对新生成的种群的个体也就是相邻两个阵元距离的差值重新进行排序。
12.子代种群完成更新，继续循环此算法过程，直到选出最优的停止。
13.采用上述技术方案所产生的有益效果在于：所述方法将免疫算法和遗传算法结合进行直线阵列天线稀布阵，不但可以减少天线阵元，大幅降低成本，达到抑制阵列天线方向图旁瓣电平的目的，而且能够克服计算量大，收敛速度较慢的问题，收敛速度显著提升。
附图说明
14.下面结合附图和具体实施方式对本发明作进一步详细的说明。
15.图1是本发明实施例所述方法的流程图；
16.图2是本发明实施例所述天线阵列整体结构几何模型图；
具体实施方式
17.下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
18.在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是本发明还可以采用其他不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类似推广，因此本发明不受下面公开的具体实施例的限制。
19.本发明实施例公开了一种于免疫遗传算法的希布阵天线的优化设计方法，所述方法设计思想是在传统遗传算法的加入免疫算子，使稀布阵列天线形成过程中不仅抑制阵列天线方向图旁瓣电平，还提高了运算效率。
20.总体的，首先随机产生初始群体，相邻两个阵元距离的差值为个体作为种群个体；通过雷达工作波长λ、阵元间距d等之间的约束关系作为疫苗进行限定；之后计算适应度函数，如果当前种群中有最佳个体，则结束算法，否则对当前种群继续进行选择交叉变异操作，得到子代种群；对子代种群进行接种疫苗操作，使得种群中具有更高的适应度的个体留下来，得到接种疫苗后种群；为防止种群个体出现退化现象，再种群进行免疫选操作，得到新一代种群，然后返回到适应度函数的计算这一步骤寻找最优个体，直到找到最佳时结束
算法。
21.进一步的，如图1所示，本技术所述方法包括如下步骤：
22.步骤一：种群初始化，建立一个均匀分布的直线阵列模型。如图2所示，阵元数为n，阵列孔径为l，d为阵元间距，采用取相邻两个阵元距离的差值为个体作为种群个体，对每个个体基因进行排序，阵元间距为n
‑
1个，在[0,l]上组成一个(n
‑
1)
×
1为列向量。表示为：
[0023][0024]
其中x
i
满足x1＜x2＜x3…
x
n
‑1∈[0,l
‑
(n
‑
2)d]，搜索空间由[0,l]缩小到[0,l
‑
(n
‑
2)d]。
[0025]
为了防止混叠现象的发生，在孔径一定的情况下，相邻两个阵元距离的差值即阵元间距应该不小于波长的一半，当阵元间距等于波长一半时最优，该步骤可以减小搜索区域，提高算法的效率。
[0026]
步骤二：该模型n为阵元个数，λ为雷达工作波长，d为阵元间距，θ为波束指向角；θ0为方位与天线法线夹角：
[0027]
则归一化方向图函数表示为：
[0028]
当上式中(n为整数)方向图函数值为1，此时存在无数个最大值点，其中(n为整数)时出现栅瓣，这将会与doa估计产生严重误差，因此为了保证不出现栅瓣，要满足因此在进行归一化方向图计算时，将公式作为疫苗。
[0029]
步骤三：适应度的计算，假设阵列孔径为l固定不变，方向图的函数表达式为其中d
n
是与第一个阵元的距离，归一化方向图旁瓣电压表达式为把归一化方向图作为适应度函数。
[0030]
步骤四：若满足上述的条件，结束算法，若不满足则进行下面遗传操作。
[0031]
步骤五：遗传操作包括选择、交叉、变异。采用轮盘赌的选择法，选择是通过个体的
适应度大小决定是否被选择，若某个个体的适应度为m，种群大小p，被选择的概率为
[0032]
上述遗传操作过程中对个体进行适当的变化，遗传操作结束之后，接种疫苗和选则疫苗实际上是加入免疫算子，免疫算子的加入目的是提高适应度和防止记忆退化现象的出现；为了保证个体满足不出现栅瓣和最小约束阵元距离。对新生成的种群的个体也就是相邻两个阵元距离的差值重新进行排序。
[0033]
步骤六：子代种群完成更新，继续循环步骤一至步骤五过程，直到选出最优停止。
[0034]
所述算法把免疫算法和遗传算法结合进行直线阵列天线稀布阵，理论层面来说，不但可以减少天线阵元，大幅降低成本，还可以达到抑制阵列天线方向图旁瓣电平的目的，在仿真小规模的天线阵模型中效果有效率明显的提升。