CN112100701A - 基于遗传算法的二维分布式天线子阵位置优化方法 - Google Patents
基于遗传算法的二维分布式天线子阵位置优化方法 Download PDFInfo
- Publication number
- CN112100701A CN112100701A CN202010760077.4A CN202010760077A CN112100701A CN 112100701 A CN112100701 A CN 112100701A CN 202010760077 A CN202010760077 A CN 202010760077A CN 112100701 A CN112100701 A CN 112100701A
- Authority
- CN
- China
- Prior art keywords
- chromosome
- dimensional
- array
- antenna
- genetic algorithm
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000005457 optimization Methods 0.000 title claims abstract description 52
- 230000002068 genetic effect Effects 0.000 title claims abstract description 45
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 41
- 210000000349 chromosome Anatomy 0.000 claims abstract description 108
- 238000003491 array Methods 0.000 claims description 23
- 230000035772 mutation Effects 0.000 claims description 23
- 238000010586 diagram Methods 0.000 claims description 17
- 108090000623 proteins and genes Proteins 0.000 claims description 9
- 238000012163 sequencing technique Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 abstract description 8
- 238000004088 simulation Methods 0.000 description 5
- 239000002131 composite material Substances 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000012804 iterative process Methods 0.000 description 2
- 238000010187 selection method Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010353 genetic engineering Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 125000006850 spacer group Chemical group 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- 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
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- 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)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Biophysics (AREA)
- Geometry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Computer Hardware Design (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Software Systems (AREA)
- Mathematical Optimization (AREA)
- Genetics & Genomics (AREA)
- Pure & Applied Mathematics (AREA)
- Medical Informatics (AREA)
- Computational Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Mathematical Analysis (AREA)
- Biomedical Technology (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了一种基于遗传算法的二维分布式天线子阵位置优化方法,包括:S1:获取二维天线阵列并初始化二维天线阵列参数;S2:根据二维天线阵列参数创建种群,并初始化遗传算法参数;其中,种群包括子阵的方位向间隔染色体集合和俯仰向间隔染色体集合;S3:根据二维天线阵列参数和遗传算法参数构造适应度函数,并计算每个染色体的适应度值;S4:根据每个染色体的适应度值对种群进行遗传迭代,并更新每个染色体的适应度值;S5:重复步骤S4,直至满足预设终止条件,得到最优二维天线阵列子阵排布构型。本发明提供的方法降低了计算量和计算复杂度,且能够得到具有较低峰值旁瓣比的二维分布式天线阵列,有效满足分布式雷达对大孔径天线的需求。
Description
技术领域
本发明属于雷达技术领域,具体涉及一种基于遗传算法的二维分布式天线子阵位置优化方法。
背景技术
随着科技的发展,雷达系统已在军用民用等各领域得到越来越广泛的应用,相应的,对雷达系统的要求也越来越高。在未来战争中,复杂的战场环境和来自目标的威胁要求雷达系统必须具备对目标远距离探测和精确跟踪能力,这意味着雷达天线需要具有高增益和大物理孔径。由于卫星载体平台的物理结构限制,很难找到足够的空间来有效安装大尺寸天线阵列,同时高昂的发射成本也限制了大尺寸天线的应用。因此,采用多个分布式小孔径子阵并对信号进行相参处理的方式来扩大天线阵列孔径是一种有效技术途径,可以提高探测距离与探测精度,增强雷达的威力。
目前,已有很多学者对一维直线阵列的稀布优化和综合方法进行了广泛而深入的研究,但较少有涉及大规模的二维分布式阵列天线的研究。二维分布式阵列天线的由于在方位向和俯仰向均布设天线子阵,相比于一维阵列具有更大的阵列孔径,可以大大提高探测距离和探测精度。因此,二维分布式天线的布阵方式研究成为天线优化设计带来了新的挑战,大规模分布式阵列的长基线间隔会导致栅瓣或高旁瓣的出现,这会严重影响雷达的工作性能。为了更大程度地发挥二维分布式雷达的优势,必须对系统发射和接收天线的位置进行优化求解。
然而,针对二维大规模的分布式阵列天线的优化,是一个复杂的非线性优化问题,一是由于二维阵列规模较大,其阵元数量数十倍于一维阵列,优化求解的计算量急剧扩大;二是约束条件更加复杂,既要保证一定的优化自由度又要保证子阵之间在方位向和俯仰向同时满足最小间距要求,且不发生重叠交叉,而常规的解析方法几乎无法对上述问题进行求解。遗传算法作为一种适用于非线性优化问题的全局优化算法,在阵列天线优化设计等复杂电磁学问题中被广泛应用。然而标准遗传算法直接推广应用于二维分布式阵列优化会带来交叉变异过程中很大概率产生不满足约束条件的解,且优化计算量大,容易陷入局部最优,优化收敛速度慢,从而无法获得最优解,不能满足分布式雷达对大孔径天线的需求。
发明内容
为了解决现有技术中存在的上述问题,本发明提供了一种基于遗传算法的二维分布式天线子阵位置优化方法。本发明要解决的技术问题通过以下技术方案实现:
一种基于遗传算法的二维分布式天线子阵位置优化方法,包括:
S1:获取二维天线阵列并初始化所述二维天线阵列参数;其中,所述二维天线阵列包括若干子阵;
S2:根据所述二维天线阵列参数创建种群,并初始化遗传算法参数;其中,所述种群包括所述子阵的方位向间隔染色体集合和俯仰向间隔染色体集合;
S3:根据所述二维天线阵列参数和所述遗传算法参数构造适应度函数,并计算每个染色体的适应度值;
S4:根据所述每个染色体的适应度值对所述种群进行遗传迭代,并更新所述每个染色体的适应度值;
S5:重复步骤S4,直至满足预设终止条件,得到最优二维天线阵列子阵排布构型。
在本发明的一个实施例中,所述子阵的排布方式满足:
所述子阵在方位向上的任意位置进行排布,所述子阵在俯仰向上均匀的分布在n条线上,且各相邻子阵之间不出现交叉重叠。
在本发明的一个实施例中,所述方位向间隔染色体集合和俯仰向间隔染色体集合中的每个基因均采用实值编码。
在本发明的一个实施例中,所述子阵的方位向间隔染色体集合的表达式为:
其中,p11<p12<…<p1m,p21<p22<…<p2m,…,pn1<pn2<…<pnm,Ga表示二维阵列方位向子阵之间的间隔总和,rij表示[0,1]区间上的随机数,m表示天线方位向子阵数目,n表示天线俯仰向子阵数目;
所述子阵的俯仰向间隔染色体集合的表达式为:
Q=[q1 q2…qn]=Ge×[r1 r2…rn];
其中,q1<q2<…<qn,ri表示[0,1]区间上的随机数,Ge表示二维阵列俯仰向子阵之间的间隔总和。
在本发明的一个实施例中,在步骤S2之后,还包括:
将所述方位向间隔染色体集合转换成行向量,并与所述俯仰向间隔染色体集合拼接成一维染色体集合。
在本发明的一个实施例中,根据所述二维天线阵列参数和所述遗传算法参数构造适应度函数,包括:
构造子阵的天线方向图;
根据所述子阵的天线方向图得到二维分布式阵列天线方向图;
以降低分布式阵列的峰值旁瓣电平为优化目标,根据所述二维分布式阵列天线方向图构造适应度函数。
在本发明的一个实施例中,所述适应度函数的表达式为:
F(θ)=max{20log10[GT(u)/max(GT(θ))]};
其中,F(·)表示适应度值,GT(·)表示二维分布式阵列天线方向图,θ表示入射信号的方位角,u的取值范围是θ中除主瓣以外的区域。
在本发明的一个实施例中,步骤S4包括:
S41:根据所述每个染色体的适应度值,采用轮盘赌方法对所述一维染色体集合进行选择操作,得到新一代种群;
S42:采用优势二维交叉方法对所述新一代种群进行交叉操作;
S43:根据变异概率对所述交叉后的种群中的每个染色体进行变异操作;
S44:更新所述每个染色体的适应度值。
在本发明的一个实施例中,步骤S42包括:
将所述新一代种群中的个体按照适应度值由高到低进行排序;
根据交叉概率Pc随机选择出第i个和第i+1个个体作为交叉的目标染色体;
在所述第i个染色体的方位向间隔染色体部分随机选择一交叉位置,作为第一交叉位置,同时确定出俯仰向间隔染色体部分的交叉位置,作为第二交叉位置;相应的,在所述第i+1个染色体上选择同样的第一交叉位置和第二交叉位置;
将所述第i个染色体中位于所述第一交叉位置前的部分与所述第i+1个染色体中位于所述第一交叉位置后的部分进行拼接,并将所述第i+1个染色体中位于所述第一交叉位置前的部分与所述第i个染色体中位于所述第一交叉位置后的部分进行拼接;同时,
将所述第i个染色体中位于所述第二交叉位置前的部分与所述第i+1个染色体中位于所述第二交叉位置后的部分进行拼接,并将所述第i+1个染色体中位于所述第二交叉位置前的部分与所述第i个染色体中位于所述第二交叉位置后的部分进行拼,生成两个新的染色体,以完成交叉操作。
在本发明的一个实施例中,步骤S43包括:
根据变异概率选择需要变异的染色体,并在该染色体中随机选择一个位置作为变异位置,随机变异该位置的基因值;
若所述变异位置位于所述方位向间隔染色体部分,则变异值取值范围为[0,Ga];若所述变异位置位于所述俯仰向间隔染色体部分,则变异值取值范围为[0,Ge];其中,Ga表示二维阵列方位向子阵之间的间隔总和,Ge表示二维阵列俯仰向子阵之间的间隔总和;
对变异后的基因值进行排序处理,以完成变异操作。
本发明的有益效果:
1、本发明通过采用子阵在方位向和俯仰向的间隔作为优化变量,使得交叉变异生成的新种群更好的满足子阵间距约束条件,从而有效地降低了计算量和计算复杂度;
2、本发明采用优势二维交叉方法来保证产生符合约束条件的新种群,该方法操作简单,同时可提高收敛效率及寻优能力,通过应用以上方法可以获得具有较低峰值旁瓣比的二维分布式天线阵列,有效满足分布式雷达对大孔径天线的需求。
以下将结合附图及实施例对本发明做进一步详细说明。
附图说明
图1是本发明实施例提供的一种基于遗传算法的二维分布式天线子阵位置优化方法示意图;
图2是本发明实施例提供的另一种一种基于遗传算法的二维分布式天线子阵位置优化方法流程图示意图;
图3是本发明实施例提供的优势二维交叉方法操作示意图;
图4是本发明实施例提供的二维分布式天线阵列示意图;
图5是本发明实施例提供的优化获得的最优子阵位置示意图;
图6是本发明实施例提供的优化获得的阵列合成天线方向图;
图7是本发明实施例提供的优化获得的阵列合成天线俯仰向截面图;
图8是本发明实施例提供的优化获得的阵列合成天线方位向截面图;
图9是本发明实施例提供的应用遗传算法进行优化的迭代进程图。
具体实施方式
下面结合具体实施例对本发明做进一步详细的描述,但本发明的实施方式不限于此。
实施例一
请参见图1,图1是本发明实施例提供的基于遗传算法的二维分布式天线子阵位置优化方法示意图,包括:
S1:获取二维天线阵列并初始化二维天线阵列参数。
本实施例获取的二维天线阵列包括若干均匀子阵,其中,天线方位向子阵数目为m,天线俯仰向子阵数目为n,则二维天线阵列由n×m个子阵构成。其中,各子阵的大小均为L×L,天线阵列发射信号的载波波长为λ;载波频率为f;目标散射点的俯仰向角度为方位向角度为θ;天线阵列方位向长度为La,俯仰向长度为Le;子阵内的阵元均匀分布且阵元间距均等于雷达发射信号的半波长λ/2。
进一步地,子阵的排布方式满足:子阵在方位向上的任意位置进行排布,在俯仰向上均匀的分布在n条线上,其中,直线位置可以任意选取;且子阵间方位向和俯仰向最小间距均为dc,各相邻子阵之间不出现交叉重叠。
本实施例通过采用子阵在方位向和俯仰向的间隔作为优化变量,使得交叉变异生成的新种群更好的满足子阵间距约束条件,从而有效地降低了计算量和计算复杂度。
S2:根据二维天线阵列参数创建种群,并初始化遗传算法参数;其中,种群包括子阵的方位向间隔染色体集合和俯仰向间隔染色体集合。
具体地,设置种群规模为M,交叉概率为Pc,变异概率为Pm,终止迭代次数为g;以子阵方位向间隔集合P和俯仰向间隔集合Q为优化变量创建种群,子阵方位向间隔集合P的取值范围为[0,Ga],俯仰向间隔集合Q的取值范围为[0,Ge],其中,
Ga是二维阵列方位向子阵之间的间隔总和,其表达式为:
Ga=La-L×m-dc×(m-1);
Ge是二维阵列俯仰向子阵之间的间隔总和,其表达式为:
Ge=Le-L×n-dc×(n-1)。
进一步地,方位向间隔染色体集合和俯仰向间隔染色体集合中的每个基因均采用实值编码。
具体地,由于二维分布式阵列天线有n×m个子阵构成,所以表示方位向间隔信息的集合P有n×m个取值,即:
其中,p11<p12<…<p1m,p21<p22<…<p2m,…,pn1<pn2<…<pnm,rij表示[0,1]区间上的随机数。
表示子阵俯仰向间隔信息的集合Q有n个取值,即:
Q=[q1 q2…qn]=Ge×[r1 r2…rn];
其中,q1<q2<…<qn,ri表示[0,1]区间上的随机数。
此外,在步骤S2之后,还包括:
将方位向间隔染色体集合转换成行向量,并与俯仰向间隔染色体集合拼接成一维染色体集合。
具体地,先将二维矩阵P按行依次首尾相连,转换成长度为n×m的行向量[p11p12…p1m p21…pnm],然后与长度为n的行向量Q拼接成长度为N=n×m+n一维染色体:
[p11 p12…p1m p21…pnm q1 q2…qn]。
S3:根据二维天线阵列参数和遗传算法参数构造适应度函数,并计算每个染色体的适应度值。
具体地,首先需将染色体中子阵间隔变量转化为每个子阵对应位置信息,以便适应度值计算。其中,每个子阵对应位置信息由方位向坐标xij和俯仰向坐标yj表示,为:
然后,构造子阵的天线方向图。
进一步地,在n×m个子阵构成的二维分布式阵列天线中,以每个子阵左下角的顶点处的阵元坐标(xi,yj)来表示该子阵的位置,设参考子阵位置为(x0,y0)且天线方向图为G(θ)。对于在二维平面中任一子阵,当已知该子阵与参考子阵方位向相距为xL,俯仰向相距为yL时,其位置可唯一确定为(x0+xL,y0+yL),该子阵的天线方向图可以进一步表示为:
然后根据子阵的天线方向图得到二维分布式阵列天线方向图,表示为:
最后,以降低分布式阵列的峰值旁瓣电平为优化目标,根据二维分布式阵列天线方向图构造适应度函数,表示为:
F(θ)=max{20log10[GT(u)/max(GT(θ))]};
其中,F(·)表示适应度值,GT(·)表示二维分布式阵列天线方向图,θ表示入射信号的方位角,u的取值范围是θ中除主瓣以外的区域。
然后根据上述适应度函数表达式计算每个染色体的适应度值。
S4:根据每个染色体的适应度值对种群进行遗传迭代,并更新每个染色体的适应度值。
请参见图2,图2是本发明实施例提供的另一种基于遗传算法的二维分布式天线子阵位置优化方法流程图示意图。
S41:根据每个染色体的适应度值,采用轮盘赌方法对一维染色体集合进行选择操作,得到新一代种群。
轮盘赌算法是一种常用的选择方法,又称比例选择方法,其基本思想是:每个个体被选中的概率与其适应度大小成正比。本实施例采用轮盘赌方法对种群进进行选择操作,得到新一代种群。
S42:采用优势二维交叉方法对新一代种群进行交叉操作。
首先,将新一代种群中的个体按照适应度值由高到低进行排序;
然后,根据交叉概率Pc随机选择出第i个和第i+1个个体作为交叉的目标染色体。
接着进行交叉操作,具体地:
在第i个染色体的方位向间隔染色体部分随机选择一交叉位置,作为第一交叉位置;同时确定出俯仰向间隔染色体部分的交叉位置,作为第二交叉位置;相应的,在第i+1个染色体上选择同样的第一交叉位置和第二交叉位置,然后进行交叉操作。
例如,请参见图3,图3是本发明实施例提供的优势二维交叉方法操作示意图,其中,在第i个染色体的P段(也即方位向间隔染色体部分)上随机选择一交叉点pmn(也即第一交叉位置),同时唯一确定出Q段(也即俯仰向间隔染色体部分)的交叉点为qm(也即第二交叉位置)。相应的,在第i+1个染色体上同样选择交叉点pmn和qm作为交叉位置。
将第i个染色体中位于第一交叉位置pmn前的部分与第i+1个染色体中位于第一交叉位置pmn后的部分进行拼接,并将第i+1个染色体中位于第一交叉位置pmn前的部分与第i个染色体中位于第一交叉位置pmn后的部分进行拼接;同时,
将第i个染色体中位于第二交叉位置qm前的部分与第i+1个染色体中位于第二交叉位置qm后的部分进行拼接,并将第i+1个染色体中位于第二交叉位置qm前的部分与第i个染色体中位于第二交叉位置qm后的部分进行拼接,生成两个新的染色体,以完成交叉操作。
本实施例采用优势二维交叉方法来保证产生符合约束条件的新种群,该方法操作简单,同时可提高收敛效率及寻优能力,通过应用以上方法可以获得具有较低峰值旁瓣比的二维分布式天线阵列,有效满足分布式雷达对大孔径天线的需求。
S43:根据变异概率对交叉后的种群中的每个染色体进行变异操作。
首先,根据变异概率选择需要变异的染色体,并在该染色体中随机选择一个位置作为变异位置,随机变异该位置的基因值。
具体地,本实施例对于每一个染色体根据变异概率Pm决定是否进行变异,若需要变异,则在染色体中随机选择一个L位置为变异点,随机变异该位置的值。
若变异位置位于方位向间隔染色体P部分,则变异值取值范围为[0,Ga];若变异位置位于俯仰向间隔染色体Q部分,则变异值取值范围为[0,Ge]。
最后,对变异后的基因值进行排序处理,以完成变异操作。
S44:更新每个染色体的适应度值。
在完成上述遗传操作后,重新对染色体的适应度值进行计算。
至此,完成一次遗传迭代。
S5:重复步骤S4,直至满足预设终止条件,得到最优二维天线阵列子阵排布构型。
具体地,对种群重复进行遗传迭代,并更新每个染色体的适应度值,然后判断迭代次数是否达到预设的终止迭代次数g,若不满足,则重复迭代更新,若满足,则终止迭代更新,并输出最低峰值旁瓣比值及对应染色体,由此得到最佳二维天线阵列子阵排布构型。
本发明提供的基于遗传算法的二维分布式天线子阵位置优化方法优化过程中采用子阵级优化方法,以均匀子阵在方位向和俯仰向的合理布设来合成大规模的二维阵列,有效地减小了计算复杂度和计算量,并且分布式子阵构型在实际应用中易于工程实现。本发明利用特殊的编码方式和改进的交叉操作方法提升了遗传算法性能,解决了二维分布式阵列优化过程中最小子阵间隔、最大阵列长度及遗传算法易产生不可行解、易早熟及收敛速度慢的问题,最终获得了阵列中每个子阵的最优位置。
实施例二
下面通过仿真实验对上述实施例一提供的基于遗传算法的二维分布式天线子阵位置优化方法的效果进行进一步说明。
实验参数:
请参见图4,图4是本发明实施例提供的二维分布式天线阵列示意图;设二维分布式阵列由10×5个均匀子阵构成,每个子阵大小为10米×10米,要求子阵间最小间隔大于1米,阵列方位向和俯仰向最大长度分别为179米、74米;所采用发射信号频率为1300MHz,信号入射到阵列时角度为0°,目标角度也为0°;设初始种群规模为200,变异概率Pm=0.2,交叉概率Pc=0.8,对图4所示二维天线阵列进行优化,并进行了100次独立重复仿真实验。
实验内容及结果:
将上述参数代入程序,在MATLAB2017a平台上进行仿真,得到如图5~9所示的实验结果,其中,图5是本发明实施例提供的优化获得的最优子阵位置示意图;图6是本发明实施例提供的优化获得的阵列合成天线方向图;图7是本发明实施例提供的优化获得的阵列合成天线俯仰向截面图;图8是本发明实施例提供的优化获得的阵列合成天线方位向截面图;图9是本发明实施例提供的应用遗传算法进行优化的迭代进程图。
通过100次独立重复仿真优化所得的适应度值的最优值、平均值、最差值、方差以及平均收敛代数如下:
最优值(dB) | 平均值(dB) | 最差值(dB) | 方差 | 平均收敛代数 |
-29.8256 | -28.7741 | -25.9545 | 0.50164 | 166.69 |
由上述仿真实验可知,通过采用本发明的方法可以将二维分布式天线方向图的峰值旁瓣比从-25dB左右优化至接近-30dB,较大幅度提升了二维分布式天线的性能;同时可以看出100次优化结果的方差也较小,说明了本发明所述方法对二维天线阵列优化能力的较高可靠性和稳定性;通过平均收敛代数可以看出,优化进程可以稳定在200代以内收敛至最优值,并且在图9中体现了本方法在获得较优结果后,适应度值依旧能够保持一定程度的下降趋势,说明了本发明所述方法能够很好地解决传统遗传算法易早熟及收敛速度慢的问题。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (10)
1.一种基于遗传算法的二维分布式天线子阵位置优化方法,其特征在于,包括:
S1:获取二维天线阵列并初始化所述二维天线阵列参数;其中,所述二维天线阵列包括若干子阵;
S2:根据所述二维天线阵列参数创建种群,并初始化遗传算法参数;其中,所述种群包括所述子阵的方位向间隔染色体集合和俯仰向间隔染色体集合;
S3:根据所述二维天线阵列参数和所述遗传算法参数构造适应度函数,并计算每个染色体的适应度值;
S4:根据所述每个染色体的适应度值对所述种群进行遗传迭代,并更新所述每个染色体的适应度值;
S5:重复步骤S4,直至满足预设终止条件,得到最优二维天线阵列子阵排布构型。
2.根据权利要求1所述的基于遗传算法的二维分布式天线子阵位置优化方法,其特征在于,所述子阵的排布方式满足:
所述子阵在方位向上的任意位置进行排布,所述子阵在俯仰向上均匀的分布在n条线上,且各相邻子阵之间不出现交叉重叠。
3.根据权利要求1所述的基于遗传算法的二维分布式天线子阵位置优化方法,其特征在于,所述方位向间隔染色体集合和俯仰向间隔染色体集合中的每个基因均采用实值编码。
5.根据权利要求1所述的基于遗传算法的二维分布式天线子阵位置优化方法,其特征在于,在步骤S2之后,还包括:
将所述方位向间隔染色体集合转换成行向量,并与所述俯仰向间隔染色体集合拼接成一维染色体集合。
6.根据权利要求1所述的基于遗传算法的二维分布式天线子阵位置优化方法,其特征在于,根据所述二维天线阵列参数和所述遗传算法参数构造适应度函数,包括:
构造子阵的天线方向图;
根据所述子阵的天线方向图得到二维分布式阵列天线方向图;
以降低分布式阵列的峰值旁瓣电平为优化目标,根据所述二维分布式阵列天线方向图构造适应度函数。
7.根据权利要求6所述的基于遗传算法的二维分布式天线子阵位置优化方法,其特征在于,所述适应度函数的表达式为:
F(θ)=max{20log10[GT(u)/max(GT(θ))]};
其中,F(·)表示适应度值,GT(·)表示二维分布式阵列天线方向图,θ表示入射信号的方位角,u的取值范围是θ中除主瓣以外的区域。
8.根据权利要求5所述的基于遗传算法的二维分布式天线子阵位置优化方法,其特征在于,步骤S4包括:
S41:根据所述每个染色体的适应度值,采用轮盘赌方法对所述一维染色体集合进行选择操作,得到新一代种群;
S42:采用优势二维交叉方法对所述新一代种群进行交叉操作;
S43:根据变异概率对所述交叉后的种群中的每个染色体进行变异操作;
S44:更新所述每个染色体的适应度值。
9.根据权利要求8所述的基于遗传算法的二维分布式天线子阵位置优化方法,其特征在于,步骤S42包括:
将所述新一代种群中的个体按照适应度值由高到低进行排序;
根据交叉概率Pc随机选择出第i个和第i+1个个体作为交叉的目标染色体;
在所述第i个染色体的方位向间隔染色体部分随机选择一交叉位置,作为第一交叉位置,同时确定出俯仰向间隔染色体部分的交叉位置,作为第二交叉位置;相应的,在所述第i+1个染色体上选择同样的第一交叉位置和第二交叉位置;
将所述第i个染色体中位于所述第一交叉位置前的部分与所述第i+1个染色体中位于所述第一交叉位置后的部分进行拼接,并将所述第i+1个染色体中位于所述第一交叉位置前的部分与所述第i个染色体中位于所述第一交叉位置后的部分进行拼接;同时,
将所述第i个染色体中位于所述第二交叉位置前的部分与所述第i+1个染色体中位于所述第二交叉位置后的部分进行拼接,并将所述第i+1个染色体中位于所述第二交叉位置前的部分与所述第i个染色体中位于所述第二交叉位置后的部分进行拼,生成两个新的染色体,以完成交叉操作。
10.根据权利要求6所述的基于遗传算法的二维分布式天线子阵位置优化方法,其特征在于,步骤S43包括:
根据变异概率选择需要变异的染色体,并在该染色体中随机选择一个位置作为变异位置,随机变异该位置的基因值;
若所述变异位置位于所述方位向间隔染色体部分,则变异值取值范围为[0,Ga];若所述变异位置位于所述俯仰向间隔染色体部分,则变异值取值范围为[0,Ge];其中,Ga表示二维阵列方位向子阵之间的间隔总和,Ge表示二维阵列俯仰向子阵之间的间隔总和;
对变异后的基因值进行排序处理,以完成变异操作。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010760077.4A CN112100701B (zh) | 2020-07-31 | 2020-07-31 | 基于遗传算法的二维分布式天线子阵位置优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010760077.4A CN112100701B (zh) | 2020-07-31 | 2020-07-31 | 基于遗传算法的二维分布式天线子阵位置优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112100701A true CN112100701A (zh) | 2020-12-18 |
CN112100701B CN112100701B (zh) | 2024-02-09 |
Family
ID=73750258
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010760077.4A Active CN112100701B (zh) | 2020-07-31 | 2020-07-31 | 基于遗传算法的二维分布式天线子阵位置优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112100701B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112986989A (zh) * | 2021-02-01 | 2021-06-18 | 西安电子科技大学 | 基于遗传算法的正交相位编码信号抑制距离模糊的方法 |
CN113540791A (zh) * | 2021-06-04 | 2021-10-22 | 中山大学 | 一种孔径级收发同时阵列优化方法 |
CN117113459A (zh) * | 2023-10-23 | 2023-11-24 | 成都睿沿芯创科技有限公司 | 一种天线阵列的旁瓣优化方法、装置及电子设备 |
CN117829003A (zh) * | 2024-03-05 | 2024-04-05 | 西安电子科技大学 | 一种多策略优化的mimo阵列设计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090012768A1 (en) * | 2004-05-13 | 2009-01-08 | Seong-Ho Son | Method for deciding array spacing of array antenna by using genetic algorithm and array antenna having sofa structure with irregular array spacing |
CN109725294A (zh) * | 2018-12-12 | 2019-05-07 | 西安电子科技大学 | 基于改进遗传算法的雷达阵列稀疏优化方法 |
CN110376557A (zh) * | 2019-06-03 | 2019-10-25 | 西安电子科技大学 | 一种基于非均匀嵌套mimo雷达的栅瓣抑制方法 |
-
2020
- 2020-07-31 CN CN202010760077.4A patent/CN112100701B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090012768A1 (en) * | 2004-05-13 | 2009-01-08 | Seong-Ho Son | Method for deciding array spacing of array antenna by using genetic algorithm and array antenna having sofa structure with irregular array spacing |
CN109725294A (zh) * | 2018-12-12 | 2019-05-07 | 西安电子科技大学 | 基于改进遗传算法的雷达阵列稀疏优化方法 |
CN110376557A (zh) * | 2019-06-03 | 2019-10-25 | 西安电子科技大学 | 一种基于非均匀嵌套mimo雷达的栅瓣抑制方法 |
Non-Patent Citations (1)
Title |
---|
包子阳;陈客松;何子述;韩春林;: "基于改进遗传算法的圆阵稀布方法", 电子技术应用, no. 10 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112986989A (zh) * | 2021-02-01 | 2021-06-18 | 西安电子科技大学 | 基于遗传算法的正交相位编码信号抑制距离模糊的方法 |
CN112986989B (zh) * | 2021-02-01 | 2023-02-24 | 西安电子科技大学 | 基于遗传算法的正交相位编码信号抑制距离模糊的方法 |
CN113540791A (zh) * | 2021-06-04 | 2021-10-22 | 中山大学 | 一种孔径级收发同时阵列优化方法 |
CN113540791B (zh) * | 2021-06-04 | 2022-12-06 | 中山大学 | 一种孔径级收发同时阵列优化方法 |
CN117113459A (zh) * | 2023-10-23 | 2023-11-24 | 成都睿沿芯创科技有限公司 | 一种天线阵列的旁瓣优化方法、装置及电子设备 |
CN117113459B (zh) * | 2023-10-23 | 2024-01-05 | 成都睿沿芯创科技有限公司 | 一种天线阵列的旁瓣优化方法、装置及电子设备 |
CN117829003A (zh) * | 2024-03-05 | 2024-04-05 | 西安电子科技大学 | 一种多策略优化的mimo阵列设计方法 |
CN117829003B (zh) * | 2024-03-05 | 2024-06-04 | 西安电子科技大学 | 一种多策略优化的mimo阵列设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112100701B (zh) | 2024-02-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112100701B (zh) | 基于遗传算法的二维分布式天线子阵位置优化方法 | |
CN106407723B (zh) | 面向低副瓣的稀疏排布阵列天线激励电流幅度的确定方法 | |
CN109725294B (zh) | 基于改进遗传算法的雷达阵列稀疏优化方法 | |
CN113343588B (zh) | 一种多约束毫米波车载mimo雷达天线阵列设计方法与装置 | |
CN107944133B (zh) | 基于多目标量子蜘蛛群演化机制的环形天线阵列稀疏方法 | |
CN104993251B (zh) | 一种大型平面阵列天线方向图综合级联优化方法 | |
CN107302140B (zh) | 基于量子蜘蛛群演化机制的平面天线阵列稀疏方法 | |
CN111160556B (zh) | 一种基于自适应遗传算法的阵列稀疏优化方法 | |
CN114117565A (zh) | 一种基于深度学习的平面阵列天线辐射方向图综合方法 | |
CN107658573B (zh) | 一种基于量子搜寻者搜索机制的圆环阵方向图综合方法 | |
CN111353605A (zh) | 基于改进遗传算法的新型平面分子阵天线阵列综合布阵方法 | |
CN114399044A (zh) | 一种子阵级稀疏阵列发射波束旁瓣电平优化方法 | |
CN110276086B (zh) | 基于空间映射的隐身超表面雷达散射截面减缩方法 | |
CN103513225B (zh) | 一种基于空间增益的稀疏平面阵形优化方法 | |
CN114218849A (zh) | 基于深度强化学习的复杂阵列天线智能设计方法 | |
CN109271735B (zh) | 基于量子启发式引力搜索算法的阵列方向图综合方法 | |
CN112016662A (zh) | 基于混合差分进化算法与加权总体最小二乘法的阵列方向图综合法 | |
CN115146544A (zh) | 一种采用知识和数据混合驱动的阵列天线设计方法 | |
CN113691294B (zh) | 一种近场稀疏阵列天线波束建立方法及装置 | |
CN115566442A (zh) | 稀疏布阵方法、装置、电子设备和可读存储介质 | |
CN111914427B (zh) | 一种基于面积归一化策略的多约束矩形阵列稀布优化方法 | |
CN109031216B (zh) | 一种基于改进遗传算法的平面阵稀疏优化方法 | |
CN110705680B (zh) | 一种分布式阵列的子阵列位置优化方法及系统 | |
CN114239380B (zh) | 基于自调整映射法则的矩形稀布阵列优化方法 | |
CN112114289B (zh) | 基于框架理论的最小阵元数稀疏阵的优化方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |