基于pso类算法的mimo雷达方阵布阵优化方法
【专利摘要】本发明提供一种基于PSO类算法的MIMO雷达方阵布阵优化方法，包括以下步骤：步骤1，对方阵进行离散化处理，并对离散孔径拓扑使用带约束条件的二进制向量来表示；步骤2，对步骤1得到的二进制向量代入MIMO雷达方向图公式来确立与方向图的一一对应关系。由方向图可以求解其对应的峰值旁瓣电平PSL，从而将PSL表示为拓扑结构二进制向量的函数，建立对应的二进制优化问题；步骤3，对步骤2获得的二进制优化问题(目标是降低PSL)采用PSO类算法进行求解，得到优化后的阵列拓扑结构，由此可以得出对应的最优PSL方向图。本发明通过PSO类算法能够有效降低MIMO雷达PSL，从而实现MIMO雷达方阵的优化设计。
【专利说明】
基于PSO类算法的ΜI MO雷达方阵布阵优化方法
技术领域
[0001]本发明属于雷达技术领域，特别是一种基于PS0类算法的ΜΙΜΟ雷达方阵布阵优化 方法。
【背景技术】
[0002] 2008年，Yarovoy研究小组于微波会议上发表了《基于近场成像的超宽带天线阵列 PS0算法布阵优化设计》一文。文中为了实现近场成像给出了两种特殊的MMO超宽带阵列， 并运用PS0算法进行优化设计使得可以在足够的横向距离分辨力和合适的主旁瓣电平比之 间得到一个权衡。研究的仿真结果显示通过对适应度函数适当的优化，可在不同场景得到 合适的非均匀线阵结构且达到了期望的权衡目标。
[0003] 国内，2013年，中南大学信息科学与工程学院的施荣华，朱炫滋研究小组给出了将 粒子群与遗传算法相混合的ΜΙΜ0雷达线阵优化设计，该小组提出了一种改进的自适应操 作，将加入自适应操作的遗传算法与粒子群算法进行混合，在混合过程中完成了 2次信息传 递:GA的初始种群由PS0中最优个体产生;经过遗传操作后再由PS0更新所有个体的速度、位 置。从仿真结果看，该算法能够有效解决非线性、全局寻优问题，无论是收敛速度还是算法 稳定性，以及寻优能力均有明显的改进。
[0004] 同年，电子科技大学电子工程学院的张伟研究小组发现传统的ΜΜ0稀疏布阵虽然 可以获得最大的连续虚拟孔径，但在机载应用中却无法进行规则的稀疏布阵，于是他们采 用模拟退火算法以Μ頂0接收端的虚拟收发联合波束为优化对象，对ΜΜ0雷达的稀疏阵列进 行了布阵优化。仿真结果证明通过该方法可以在保持主瓣不展宽的情况下获得更好的旁瓣 水平。
[0005] 但是现有技术中，没有将PS0算法与DE算法应用于ΜΜ0雷达布阵的优化，因此旁瓣 电平较高，且目前研究仅针对一维线阵进行优化，没有探究二维线阵优化的工作。
【发明内容】
[0006] 本发明的目的在于提供一种基于PS0类算法的ΜΜ0雷达方阵布阵优化方法，以获 得更低的旁瓣电平，从而提高雷达性能。
[0007] 实现本发明目的的技术解决方案为:一种基于PS0类算法的ΜΜ0雷达方阵布阵优 化方法，步骤如下：
[0008] 步骤1，对方阵进行离散化处理，并对离散孔径拓扑使用带约束条件的二进制向量 表不；
[0009] 步骤2，将步骤1得到的二进制向量代入ΜΜ0雷达方向图公式，确立该二进制向量 与ΜΜ0雷达方向图的--对应关系，并求解ΜΙΜ0雷达方向图对应的峰值旁瓣电平PSL，将 PSL表示为拓扑结构二进制向量的函数，从而建立二进制优化问题；
[0010] 步骤3,对步骤2建立的二进制优化问题采用PS0类算法进行求解，得到优化后的阵 列拓扑结构，最终得到最优PSL方向图。
[0011] 进一步地，步骤1所述对方阵进行离散化处理，并对离散孔径拓扑使用带约束条件 的二进制向量表示，具体步骤如下：
[0012] 步骤1.1，设正方形边长L为:L = aA，a为常数，λ为雷达工作波长;ΜΙΜΟ雷达的发射 阵元数目Nt和接收阵元数目Nr满足:Nt = Nr = N，N为正整数；
[0013] 步骤1.2,阵元的间距满足大于或者等于0.5λ，将正方形按0.5λ的间距离散化为 /7 2 个端点，即8a个端点； 0.5 Λ
[0014] 步骤1.3,阵元离散化位置应遵循以下规则:发射、接收阵列中，正方形的四个端点 处都放置阵元，其他的Ν-4个阵元需要放置在剩余的8a-4个端点上，用g表示ΜΜ0雷达发射 阵列拓扑结构，/;,来表示Μ頂0雷达接收阵列拓扑结构，即：
[0015]
[0016] 其中，btk，brke {〇，1}，1表示第k个端点放置阵元，0表示第k个端点不放置阵元，k =1，2,3."8a;
[0017] 步骤1.4,方阵的离散孔径拓扑使用带约束条件的二进制向量表示如下： k~x li-i :〇
[0023] 进一步地，步骤2所述将步骤1得到的二进制向量代入ΜΙΜΟ雷达方向图公式，确立 该二进制向量与ΜΙΜΟ雷达方向图的一一对应关系，并求解ΜΙΜΟ雷达方向图对应的峰值旁瓣 电平PSL，将PSL表示为拓扑结构二进制向量的函数，从而建立二进制优化问题，具体步骤如 下：
[0024] 步骤2.1，二维MBTO雷达方向图/7(0,供)的公式为：
[0025]
[0026] 其中，λ表示为雷达工作波长，Θ表示入射波方向的仰角，表示入射波方向的方位 角，Nt和Nr分别为ΜΙΜΟ雷达的发射阵列与接收阵列的阵兀数量，&为发射阵列的阵兀位置， &为接收阵列的阵兀位置，i?为入射波的入射方向单位矢量；
[0027] 步骤2.2，将带约束条件的二进制向量代入步骤2.1的Μ頂0雷达方向图公式来确立 表示拓扑结构的二进制向量与方向图的一一对应关系；
[0028] 步骤2.3，由步骤2.2求解峰值旁瓣电平PSL(/〗,,/i,.):，将PSL表示为拓扑结构二进制 向量的函数，从而建立二进制优化问题如下：
[0029] min PSL^k.b^
[0035] 进一步地，步骤3所述对步骤2建立的二进制优化问题采用PS0类算法进行求解，得 到优化后的阵列拓扑结构，最终得到最优PSL方向图，具体步骤如下：
[0036] 步骤3.1，为解决雷达布阵问题，需求解变量为0或1的二进制优化问题：
[0037]
[0038] subjectto xqe{〇,l},f〇r
[0039] 其中，目标函数fo(X1，x2, . . .，xQ)是Q个实数变量Χ1，χ2, . . .，xQ的函数，并且每个变 量Xq只能取值〇或l，q=l，2,…，Q;
[0040] 步骤3.2,给定种群规模和最大迭代次数后，初始化种群中第m个粒子的位置向量 毛=^,1八2，...,\^，对于祀1^0^?1为第111个粒子在〇维空间中每个维度上的位置坐标，4: 中的元素只能取0或者1，通过产生随机的0、1整数来初始化粒子的位置向量；
[0041 ] 步骤3.3，对于BinPSO，速度、位置更新公式为：
[0042]
[0043]其中，w为惯性系数，C1是粒子追踪自己历史最优解的权重，c2是粒子追踪群体历史 最优的权重，和分别为第k次迭代中第m个粒子在第q维的速度和位置坐标，为第k 次迭代中第m个粒子在第q维自己搜索到历史最优位置坐标，为第k+1次迭代中第m个粒 子在第q维的速度，为第m个粒子在第k次迭代中第q维搜索到的全局最优位置；
[0044] 步骤3.4，步骤3.3中速度的选定存在限定条件：
[0045]
[0046] 若更新后的速度大于设定的最大值便强制约束在速度最大值Vqmax，或者小于设定 的最小值便强制约束在速度最小值Vqmin ;
[0047] 步骤3.5，位置的迭代公式为：
[0048]
[0049] 其中，为第m个粒子在第k+Ι次迭代中第m个粒子在第q维的位置坐标，为第m 个粒子在第k次迭代中第m个粒子在第q维的位置坐标；
[0050] 对于BinPSO算法位置的选定存在如下限定条件：
[0051]
[0052] 其中，rand表示均勾分布在[0,1]的随机数，S()表示sigmoid函数，即 为了防止sigmoid函数的饱和性，粒子的速度范围[Vqmin，Vqmax]为[-4,4];
[0053] 步骤3.6,按照步骤3.3给出的速度、位置更新公式以及步骤3.5给出的位置的迭代 公式更新粒子的速度与位置；
[0054]步骤3.7，确定最优解并且判断是否满足终止条件，如果满足即算法完成，否则返 回步骤3.1计算新一代种群的适应度函数值并进入循环，直到最优解满足终止条件，或者当 迭代的次数超过了步骤3.2预先设置的最大迭代次数，则强制结束，得到最优PSL方向图。
[0055] 本发明与现有技术相比，其显著优点为：（1)采用PS0类算法，选取了合适的阵列长 度和阵元数目，完成了对ΜΜ0雷达方阵的布阵优化设计，获得了理想的PSL; (2)适用于一般 性虚拟孔径雷达、一般性智能优化算法、一般性二维阵列，获得更低的旁瓣电平，提高了雷 达性能；（3)有效可靠，优化效果明显。
【附图说明】
[0056] 图1为本发明基于PS0类算法的MBTO雷达方阵布阵优化方法的流程图。
[0057] 图2为本发明中PS0类算法优化流程图。
[0058]图3为本发明中经PS0类算法优化后的方阵ΜΜ0雷达阵列拓扑结构图，其中（a)为 发射阵列图，（b)为接收阵列图。
[0059] 图4为本发明经PS0类算法优化后的方阵MBTO雷达方向图。
【具体实施方式】
[0060] 结合图1，本发明基于粒子群（Particle Swarm 0ptimization，PS0)PS0类算法的 Μπω雷达方阵布阵优化方法，该布阵方法不仅适用于本发明所述的方阵，对任何可做离散 化处理的二维阵型皆适用，具体步骤如下：
[0061] 步骤1，对方阵进行离散化处理，并对离散孔径拓扑使用带约束条件的二进制向量 表示，具体步骤如下：
[0062] 步骤1. 1，设正方形边长L为:L = aA，a为常数，λ为雷达工作波长(约3cm) ;ΜΙΜ0雷达 的发射阵元数目Nt和接收阵元数目Nr满足:Nt = Nr = N，N为正整数；
[0063]步骤1.2,为了减小天线阵元之间的互耦效应，阵元的间距满足大于或者等于0.5 (1 / 入，将正方形按0.5λ的间距离散化为4X&-个端点，g卩8a个端点；
[0064] 步骤1.3，为了获得最大的物理孔径，阵元离散化位置应遵循以下规则:发射、接收 阵列中，正方形的四个端点处都放置阵元，其他的N-4个阵元需要放置在剩余的8a_4个端点 上，用￡表示Μ頂0雷达发射阵列拓扑结构，反来表示Μ頂0雷达接收阵列拓扑结构，即：
[0065]
[0066] 其中，btk，brk e {〇，1}，1表示第k个端点放置阵元，0表示第k个端点不放置阵元，k =1，2,3."8a;
[0067] 步骤1.4,方阵的离散孔径拓扑使用带约束条件的二进制向量表示如下：
[0068]
[0069]
[0070]
[0071]
[0072] .A：.=丄 妒=丄.
[0073] 步骤2,将步骤1得到的二进制向量代入MM0雷达方向图公式，确立该二进制向量 与ΜΜ0雷达方向图的--对应关系，并求解ΜΙΜΟ雷达方向图对应的峰值旁瓣电平PSL，将 PSL表示为拓扑结构二进制向量的函数，从而建立二进制优化问题，具体步骤如下：
[0074] 步骤2 · 1，二维MBTO雷达方向图舛武命)的公式为：
[0075]
[0076] 其中，λ表示为雷达工作波长(约3cm)，Θ表示入射波方向的仰角，史表示入射波方向 的方位角，Nt和N r分别为ΜΠΚ)雷达的发射阵列与接收阵列的阵元数量，ξ:为阵元分布在Χ0Υ 平面上发射阵列的阵元位置，< 为阵元分布在Χ0Υ平面上接收阵列的阵元位置，δ为入射波 的入射方向单位矢量；
[0077] 步骤2.2，将带约束条件的二进制向量代入步骤2.1的Μ頂0雷达方向图公式来确立 表示拓扑结构的二进制向量与方向图的一一对应关系；
[0078] 步骤2.3,由步骤2.2求解峰值旁瓣电平PSL(/^.)，将PSL表示为拓扑结构二进制 向量的函数，从而建立二进制优化问题如下：
[0079]
[0080]
[0081]
[0082]
[0083]
[0084]
[0085] 步骤3,对步骤2建立的二进制优化问题（目标是降低PSL)采用PS0类算法进行求 解，得到优化后的阵列拓扑结构，最终得到最优PSL方向图，具体步骤如下：
[0086] 步骤3.1，为解决雷达布阵问题，需求解变量为0或1的二进制优化问题：
[0087] minimize ./；, (.vr A2,...,.ry) ·ΥΙ，·γ:'···，.、.{?
[0088] subjeetto xqe{〇,l},f〇r
[0089] 其中，目标函数fo(X1，X2, . . .，xQ)是Q个实数变量Χ1，χ2, . . .，xQ的函数，并且每个变 量Xq只能取值〇或l，q=l，2,…，Q;
[0090] 步骤3.2,给定种群规模和最大迭代次数后，初始化种群中第m个粒子的位置向量 无》=(:^，^,-.,"^〇，对于扮1^〇，、为第 111个粒子在〇维空间中每个维度上的位置坐标，4 中的元素只能取0或者1，通过产生随机的0、1整数来初始化粒子的位置向量；
[0091] 步骤3.3，对于BinPSO，速度、位置更新公式为：
[0092]
[0093] 其中，w为惯性糸数，(^是敉于退踪目己)力史最优觯的权重，c2是粒子追踪群体历史 最优的权重，it和/4分别为第k次迭代中第m个粒子在第q维的速度和位置坐标，为第k 次迭代中第m个粒子在第q维自己搜索到历史最优位置坐标，为第k+1次迭代中第m个粒 子在第q维的速度，为第m个粒子在第k次迭代中第q维搜索到的全局最优位置；
[0094] 步骤3.4，步骤3.3中速度的选定存在限定条件：
[0095]
[0096] 若更新后的速度大于设定的最大值便强制约束在速度最大值Vqmax，或者小于设定 的最小值便强制约束在速度最小值Vqmin ;
[0097] 步骤3.5，位置的迭代公式为：
[0098]
[0099] 其中，为第m个粒子在第k+Ι次迭代中第m个粒子在第q维的位置坐标，为第m 个粒子在第k次迭代中第m个粒子在第q维的位置坐标；
[0?00]对于BinPSO算法位置的选定存在如下限定条件：
[0101]
[0102] 其中，rand表示均勾分布在[0，1 ]的随机数，S ()表示s i gmo i d函数，即 5(^) = '!~^~^:，为了防止8181]1〇1(1函数的饱和性，粒子的速度范围[￥-，'\^111￡?]为[-4,4] ;该限 7 1 -f e ' 定条件对应步骤1.3的1表示放置阵元而0表示不放置阵元；
[0103] 步骤3.6,按照步骤3.3给出的速度、位置更新公式以及步骤3.5给出的位置的迭代 公式更新粒子的速度与位置；
[0104] 步骤3.7，确定最优解并且判断是否满足终止条件，如果满足即算法完成，否则返 回步骤3.1计算新一代种群的适应度函数值并进入循环，直到最优解满足终止条件，或者当 迭代的次数超过了步骤3.2预先设置的最大迭代次数，则强制结束，根据步骤2.3得到最优 PSL方向图。
[0105] 实施例1
[0106] 本发明基于PS0类算法的MBTO雷达方阵布阵优化方法，具体步骤如下：
[0107] 步骤1，对方阵进行离散化处理，并对离散孔径拓扑使用带约束条件的二进制向量 来表示；
[0108] 步骤2，对步骤1得到的二进制向量代入ΜΜ0雷达方向图公式来确立与方向图的一 一对应关系。由方向图可以求解其对应的峰值旁瓣电平PSL，从而将PSL表示为拓扑结构二 进制向量的函数，建立对应的二进制优化问题；
[0109] 步骤3,对步骤2获得的二进制优化问题（目标是降低PSL)采用PS0类算法进行求 解，得到优化后的阵列拓扑结构，由此可以得出对应的最优PSL方向图；
[0110]步骤1选择对方阵进行离散化处理，并对离散孔径拓扑使用带约束条件的二进制 向量来表示的具体方法为：
[0111] 步骤1.1，此处举例假设正方形边长为? = 4λ，ΜΜ0雷达的发射、接收阵元数目Nt = Nr=16;
[0112] 步骤1.2,为了减小天线阵元之间的互耦效应，阵元的间距需要满足大于或者等于 0.5λ，可将正方形按〇. 5λ的间距离散化为32个端点；
[0113] 步骤1.3,为了获得最大的物理孔径，阵元离散化位置应遵循:发射、接收阵列中， 正方形的四个端点处都放置阵元，其他的12个阵元需要合理的放置在剩余的28个端点上。 用ζ表示ΜΙΜΟ雷达发射阵列拓扑结构，民来表示ΜΙΜΟ雷达接收阵列拓扑结构，即 4=1? 夂2 …t3J，其中bti，brie {〇，1}，1表示第i个端点放置 阵元而0表示不放置阵元；
[0114] 步骤1.4,方阵的离散孔径拓扑使用带约束条件的二进制向量来表示为：
[0115]
[0116]
[0117]
[0118]
[0119] 1--1. 1 - 1
[0120] 步骤2对步骤1得到的二进制向量代入ΜΜ0雷达方向图公式来确立与方向图的一 一对应关系。由方向图可以求解其对应的峰值旁瓣电平PSL，从而将PSL表示为拓扑结构二 进制向量的函数，建立对应的二进制优化问题主要包括如下步骤：
[0121] 步骤2.1，二维MBTO雷达方向图公式为：
[0122]
[0123] 其中λ表示波长，Θ表示入射波方向的仰角，@表示入射波方向的方位角；
[0124] 步骤2.2，将步骤1.4中带约束条件的二进制向量代入步骤2.1的ΜΜ0雷达方向图 公式来确立表示拓扑结构的二进制向量与方向图的一一对应关系；
[0125] 步骤2 · 3，由步骤2 · 2可以求解其对应的峰值旁瓣电平PSL(/^.)，将PSL表示为拓 扑结构二进制向量的函数，从而建立对应的二进制优化问题：
[0126]
[0127]
[0128]
[0129] bti，brie {〇，ι}
[0130] btl = bt9 = btl7 = bt25 = brl = br9 = brl7 = br32= 1
[0131]
i~l /=1
[0132] 结合图2,步骤3对步骤2获得的二进制优化问题采用BinPSO算法进行求解，得到优 化后的阵列拓扑结构，由此可以得出对应的最优PSL方向图主要包括如下步骤：
[0133] 步骤3.1，为解决雷达布阵问题，需求解变量为0或1的二进制优化问题：
[0134]
[0135] subjectto xqe{〇,l},f〇r
[0136] 其中目标函数fo(xi,X2, . . .，xq)是Q个实数变量的函数，并且每个变量xq只能取值0 或1。
[0137] 步骤3.2,给定种群规模后，初始化种群中粒子的位置向量萬对于 B inPSO，焉中的元素只能取0或者1，通过产生随机的0,1整数来初始化粒子的位置向量；
[0138] 步骤3.3，对于BinPSO，速度、位置更新公式为：
[0139]
[0140]其中，w为惯性系数，C1是粒子追踪自己历史最优解的权重，c2是粒子追踪群体历史 最优解的权重，ξ与η是均匀分布在[0，1 ]的随机数。
[0141] 步骤3.4，步骤3.3中速度的选定存在限定条件：
[0142]
[0143] 表明若更新后的速度大于设定的最大值或者小于设定的最小值，那么便强制其约 束在最大值或者最小值；
[0144] 步骤3.5，位置的迭代公式为：xf二為X1。对于BinPSO算法位置的选定存在如 下限定条件：
[0145]
[0146] 其中rand表示均勾分布在[0,1]的随机数，S()表示sigmoid函数，即5'(?) = ~丁:， 为了防止sigmoid函数的饱和性，粒子的速度范围[Vqmin，Vqmax]建议为[_4,4];该限定条件对 应步骤1.3的1表示放置阵元而0表示不放置阵元；
[0147] 步骤3.6,按照步骤3.3给出的速度、位置更新公式以及步骤3.5给出的位置的迭代 公式更新粒子的速度与位置；
[0148] 步骤3.7，确定最优解并且判断是否满足终止条件，如果满足即算法完成，否则计 算新一代种群的适应度函数值并且进入循环，直到某一代的最优解是满足终止条件的，或 者当迭代的次数超过了预先设置的最大迭代次数，则算法也强制完成。其中，适应度函数为 步骤2.4中的峰值旁瓣电平PSL( Hj。
[0149] BinPSO参数设置如下：
[0150] 种群规模:N=100
[0151] 最大迭代次数:K = 50
[0152] 粒子速度范围：[Vqmin，Vqmax] = [-4，4]
[0153] 惯性系数:wmax = 0.95，Wmin = 0.4，按照迭代次数线性递减
[0154] 粒子追踪自己历史最优解的权重:cl = 1.4
[0155] 粒子追踪群体历史最优解的权重:c2 = l .4
[0156] 通过以上步骤最后就能得到经过PS0类算法优化后的最优PSL。
[0157] 结合图3、图4:图3(a)是ΜΗ?雷达方阵采用BinPSO优化后的发射阵列拓扑结构，图 3(b)是方阵ΜΠΚ)雷达BinPSO优化后的接收阵列拓扑结构，图4是由BinPSO优化后方阵ΜΜ0 雷达的方向图，其峰值旁瓣电平为-9.5915dB。
[0158] 综上所述，本发明通过PS0类算法能够有效降低ΜΜ0雷达PSL，实现了 ΜΜ0雷达方 阵的优化设计。
【主权项】
1. 一种基于PSO类算法的ΜΙΜΟ雷达方阵布阵优化方法，其特征在于，步骤如下： 步骤1，对方阵进行离散化处理，并对离散孔径拓扑使用带约束条件的二进制向量表 示； 步骤2,将步骤1得到的二进制向量代入ΜΙΜΟ雷达方向图公式，确立该二进制向量与 ΜΙΜΟ雷达方向图的一一对应关系，并求解ΜΙΜΟ雷达方向图对应的峰值旁瓣电平PSL，将PSL 表示为拓扑结构二进制向量的函数，从而建立二进制优化问题； 步骤3,对步骤2建立的二进制优化问题采用PS0类算法进行求解，得到优化后的阵列拓 扑结构，最终得到最优Ρ化方向图。2. 根据权利要求书1所述的基于PS0类算法的ΜΙΜΟ雷达方阵布阵优化方法，其特征在 于，步骤1所述对方阵进行离散化处理，并对离散孔径拓扑使用带约束条件的二进制向量表 示，具体步骤如下： 步骤1.1，设正方形边长L为:L = aA，a为常数，λ为雷达工作波长;ΜΙΜΟ雷达的发射阵元 数目Nt和接收阵元数目Nr满足:Nt = Nr = N，N为正整数； 步骤1.2,阵元的间距满足大于或者等于0.5λ，将正方形按〇.5λ的间距离散化为4χ^ 个端点，即8a个端点； 步骤1.3,阵元离散化位置应遵循W下规则:发射、接收阵列中，正方形的四个端点处都 放置阵元，其他的N-4个阵元需要放置在剩余的8a-4个端点上，用4表示ΜΙΜΟ雷达发射阵列 拓扑结构，4来表示ΜΙΜΟ雷达接收阵列拓扑结构，即：其中，6*1<，13*￡{0，1}，1表示第4个端点放置阵元，0表示第4个端点不放置阵元^=1，2，3...8a; 步骤1.4,方阵的离散孔径拓扑使用带约束条件的二进制向量表示如下：3. 根据权利要求书1所述的基于PS0类算法的ΜΙΜΟ雷达方阵布阵优化方法，其特征在 于，步骤2所述将步骤1得到的二进制向量代入ΜΙΜΟ雷达方向图公式，确立该二进制向量与 ΜΙΜΟ雷达方向图的一一对应关系，并求解ΜΙΜΟ雷达方向图对应的峰值旁瓣电平PSL，将PSL 表示为拓扑结构二进制向量的函数，从而建立二进制优化问题，具体步骤如下： 步骤2.1，二维ΜΙΜΟ雷达方向图ρ( A刮的公式为：其中，λ表示为雷达工作波长，θ表示入射波方向的仰角，g表示入射波方向的方位角，Nt 和Nr分别为ΜΙΜΟ雷达的发射阵列与接收阵列的阵元数量，克为发射阵列的阵元位置，是,为 接收阵列的阵元位置，友为入射波的入射方向单位矢量； 步骤2.2，将带约束条件的二进制向量代入步骤2.1的ΜΙΜΟ雷达方向图公式来确立表示 拓扑结构的二进制向量与方向图的一一对应关系； 步骤2.3，由步骤2.2求解峰值旁瓣电平PSL(6，&).将PSL表示为拓扑结构二进制向量 的函数，从而建立二进制优化问题如下：4.根据权利要求书1所述的基于PS0类算法的ΜΙΜΟ雷达方阵布阵优化方法，其特征在 于，步骤3所述对步骤2建立的二进制优化问题采用PS0类算法进行求解，得到优化后的阵列 拓扑结构，最终得到最优Ρ化方向图，具体步骤如下： 步骤3.1，为解决雷达布阵问题，需求解变量为0或1的二进制优化问题：subjectto XqG {〇,!},for l^q^Q 其中，目标函数f〇(Xl，X2, . . .，XQ)是Q个实数变量X1，X2, . . .，XQ的函数，并且每个变量Xq 只能取值0或l，q = l，2，'-，Q; 步骤3.2,给定种群规模和最大迭代次数后，初始化种群中第m个粒子的位置向量 4二(而,1，爲2，..，，而2)，对于81沾50，义。。为第111个粒子在9维空间中每个维度上的位置坐标，4 中的元素只能取0或者1，通过产生随机的〇、1整数来初始化粒子的位置向量； 步骤3.3，对于BinPSO，速度、位置更新公式为：其中，W为惯性系数，C1是粒子追踪自己历史最优解的权重，C2是粒子追踪群体历史最优 的权重，乂W和分别为第k次迭代中第m个粒子在第q维的速度和位置坐标，诚g为第k次迭 代中第m个粒子在第q维自己捜索到历史最优位置坐标，1^+/为第k+1次迭代中第m个粒子在 第q维的速度，；4。,,/为第m个粒子在第k次迭代中第q维捜索到的全局最优位置； 步骤3.4，步骤3.3中速度的选定存在限定条件：若更新后的速度大于设定的最大值便强制约束在速度最大值Vq max,或者小于设定的最 小值便强制约束在速度最小值Vq min; 步骤3.5,位置的迭代公式为：其中，为第m个粒子在第k+1次迭代中第m个粒子在第q维的位置坐标，苗g为第m个粒 子在第k次迭代中第m个粒子在第q维的位置坐标； 对于BinPSO算法位置的选定存在如下限定条件：其中，rand表示均匀分布在[0，：[]的随机数，S〇表示sigmoid函数，即，为了 防止sigmoid函数的饱和性，粒子的速度范围[Vq min, Vq max]为[-4,4]; 步骤3.6，按照步骤3.3给出的速度、位置更新公式w及步骤3.5给出的位置的迭代公式 更新粒子的速度与位置； 步骤3.7,确定最优解并且判断是否满足终止条件，如果满足即算法完成，否则返回步 骤3.1计算新一代种群的适应度函数值并进入循环，直到最优解满足终止条件，或者当迭代 的次数超过了步骤3.2预先设置的最大迭代次数，则强制结束，得到最优P化方向图。
【文档编号】G06F19/00GK105975746SQ201610270932
【公开日】2016年9月28日
【申请日】2016年4月27日
【发明人】顾陈, 侯云飞, 李聪羚, 杨宇鑫, 洪弘, 李彧晟, 朱晓华
【申请人】南京理工大学