基于非接触式生命体征监测系统的一维雷达布阵优化方法
【专利摘要】本发明公开了一种基于非接触式生命体征监测系统的一维雷达布阵优化方法，步骤如下：在考虑心脏尺寸的条件下，根据雷达系统载频，确定阵列的孔径与阵元的数目，对雷达一维阵列孔径进行离散化处理并对离散孔径拓扑使用带约束条件的二进制向量来表示；对得到的二进制向量代入MIMO雷达方向图公式来确立与方向图的一一对应关系，由方向图可以求解其对应的峰值旁瓣电平PSL，从而将PSL表示为拓扑结构二进制向量的函数，建立对应的二进制优化问题；采用智能优化算法确定优化后的阵列拓扑结构，得到对应的最优PSL方向图。本发明通过智能优化算法能够有效降低非接触式生命体征监测系统雷达的PSL，从而实现非接触式生命体征监测系统的一维雷达布阵优化设计。
【专利说明】
基于非接触式生命体征监测系统的一维雷达布阵优化方法
技术领域
[0001] 本发明属于雷达技术领域，特别是一种基于非接触式生命体征监测系统的一维雷 达布阵优化方法。
【背景技术】
[0002] 人体的呼吸、心跳、脉搏、血压等是生命体征的重要指标，常被用来判断人体的生 理健康状况。这些指标如果出现异常往往表征着被测人人体可能出现的病变。因此对于临 床医学的诊断和治疗来说，实时监测这些生命体征信息具有非常重要的意义。目前监测系 统多为接触式的，但接触式监测系统有很多缺陷，因此非接触式生命体征监测技术的研究 便进入了科学家的视线，成为了当今社会一个研究的热点。
[0003]近年来国外学者提出一种新的非接触式生命体征监测技术一多普勒雷达相位调 制法，该新技术引出了有关雷达方面的各项研究，例如优化设计雷达阵列，目的是使得雷达 理论更加适用于非接触式生命体征监测领域。因为雷达收发阵列中的阵元位置和激励分布 影响到雷达的方向图，从而直接影响整个雷达系统的性能，目前尚没有针对具体参数来优 化设计阵元位置和激励分布的方法。
【发明内容】
[0004] 本发明的目的在于提供一种基于非接触式生命体征监测系统的一维雷达布阵优 化方法，以获得更低的旁瓣电平，从而提高雷达性能。
[0005] 实现本发明目的的技术解决方案为:一种基于非接触式生命体征监测系统的一维 雷达布阵优化方法，步骤如下：
[0006] 步骤1，在考虑心脏尺寸的条件下，根据具体的雷达系统载频，确定阵列的孔径与 阵元的数目，对雷达一维阵列孔径进行离散化处理并对离散孔径拓扑使用带约束条件的二 进制向量来表示；
[0007] 步骤2，将步骤1得到的二进制向量代入ΜΜ0雷达方向图公式，确立二进制向量与 ΜΙΜΟ雷达方向图的--对应关系，由ΜΙΜΟ雷达方向图求解峰值旁瓣电平PSL，并将PSL表示 为拓扑结构二进制向量的函数，从而建立二进制优化问题；
[0008] 步骤3,对步骤2获得的二进制优化问题采用智能优化算法进行求解，得到优化后 的阵列拓扑结构，并确定最优PSL方向图。
[0009] 进一步地，步骤1所述在考虑心脏尺寸的条件下，根据具体的雷达系统载频，确定 阵列的孔径与阵元的数目，对雷达一维阵列孔径进行离散化处理并对离散孔径拓扑使用带 约束条件的二进制向量来表示，具体步骤为：
[0010] 步骤1.1，假设雷达的工作频率为f，人体组织的介电常数为er，光速为c，得到雷达 工作波长
[0011] 步骤1.2,确定阵列的孔径大小，孔径的大小满足心脏完全被覆盖，假设人体心脏 的尺寸为则设计一个大于Μ的一维线阵，即构造孔径L = aX>M，a为常数系数， M、N的单位均为cm;
[0012]步骤1.3,确定一维线阵的阵元数沁，1'11 =队=队，队、队分另1」为发射阵元数和接收阵 元数；
[0013] 步骤1.4,一维线阵的阵元间距满足大于或者等于0.5λ，将线阵按0.5λ的间距离散 化为7^7 +1个端点，即2a+l个端点；
[0014] 步骤1.5,一维线阵的阵元离散化位置应遵循以下规则：一维线阵发射阵列、接收 阵列首尾端点都放置阵元，剩余的见-2个发射阵元与Λ-2个接收阵元放置在各自剩余的2a+ 1 -2个端点，即2a-1个端点上;使用二进制向量I；表示ΜΜ0雷达发射阵列拓扑结构，二进制 向量€来表示Μ頂0雷达接收阵列拓扑结构，即
[0015]
[0016]
[0017]步骤1.6,对步骤1.5得到的离散孔径拓扑使用带约束条件的二进制向量来表示， 如下所示：
[0018]
[0019]
[0020]
[0021]
[0022]
[0023] 式中，i = l，2，...，2a+l。
[0024] 进一步地，步骤2所述将步骤1得到的二进制向量代入ΜΙΜΟ雷达方向图公式，确立 二进制向量与MMO雷达方向图的--对应关系，由ΜΜ0雷达方向图求解峰值旁瓣电平PSL， 并将PSL表示为拓扑结构二进制向量的函数，从而建立二进制优化问题，具体步骤如下：
[0025] 步骤2.1，MBTO雷达的一维方向图公式为：
[0026]
[0027]其中λ表示ΜΠΚ)雷达工作波长，Θ表示平面波与一维阵列法线的夹角，ρ表示方向图 函数，Xtm表示发射阵列的第m个阵元位置，Xrn表示接收阵列的第η个阵元位置；
[0028] 步骤2.2，将带约束条件的二进制向量代入步骤2.1中的ΜΙΜΟ雷达方向图公式来表 示拓扑结构的二进制向量与方向图的一一对应关系；
[0029] 步骤2.3,求解峰值旁瓣电平PSL(^.)，将PSL表示为拓扑结构二进制向量的函 数，从而建立二进制优化问题，则一维线阵ΜΜ0雷达优化问题表述为：
[0030]
[0031]
[0032]
[0033]
[0034]
[0035]
[0036] 式中，l = l，2，."，2a+l。
[0037] 进一步地，步骤3所述对步骤2获得的二进制优化问题采用智能优化算法进行求 解，得到优化后的阵列拓扑结构，并确定最优PSL方向图，具体步骤如下：
[0038] 步骤3.1，为解决雷达布阵问题，需求解变量为0或1的二进制优化问题：
[0039]
[0040] subjectto xqG {〇，1}，f〇rKq^iQ
[0041 ]其中，目标函数f〇(Xl，X2, . . .，XQ)是Q个实数变量X1，X2, . . .，XQ的函数，并且每个变 量Xq只能取值〇或l，q=l，2,…，Q;
[0042]步骤3.2，给定种群规模后，初始化种群中粒子的位置向量& = (λ?2,. 对于BinPS0，xgq表示第g个粒子的第q个变量，旯中的元素只能取0或者1，通过产生随机的 〇、1整数来初始化粒子的位置向量；
[0043] 步骤3.3，对于BinPSO，速度、位置更新公式为：
[0044]
[0045] 其中，w为惯性系数，C1是粒子追踪自己历史最优解的权重，C2是粒子追踪群体历史 最优解的权重，ξ与η是均匀分布在[0，1]的随机数，表示第i个粒子第k+l次迭代的速度 的第q维，表示第i个粒子第k次迭代的速度的第q维，/<表示第i个粒子第k次迭代后搜索 到的历史最优位置的第q维，·<表示第i个粒子第k次迭代后位置的第q维，表示粒子第 k次迭代后搜索到的全局最优位置的第q维，粒子群算法的最大迭代次数设定为Κ，当迭代次 数超过K时，粒子群算法结束；
[0046] 步骤3.4，步骤3.3中速度的选定存在限定条件：
[0047]
[0048] 若更新后的速度大于设定的最大值便强制约束在速度最大值Vqmax，或者小于设定 的最小值便强制约束在速度最小值Vqmin ;
[0049] 步骤3.5，位置的迭代公式为：
[0050] X?,7
[0051] 其中，xf表示第i个粒子第k+Ι次迭代的位置的第q维，< 表示第i个粒子第k次迭 代的位置的第q维；
[0052]对于BinPSO算法位置的选定存在如下限定条件：
[0053]
[0054] 其中，rand表示均勾分布在[0,1]的随机数，S()表示sigmoid函数，即 S=L·.，为了防止sigmoid函数的饱和性，粒子的速度范围[Vqmin，Vqmax]为[_ 4，4];丄表 1 +€' 示放置阵元而〇表示不放置阵元；
[0055] 步骤3.6,按照步骤3.3给出的速度、位置更新公式以及步骤3.5给出的位置的迭代 公式更新粒子的速度与位置；
[0056]步骤3.7，确定最优解并且判断是否满足终止条件，如果满足即算法完成，否则计 算新一代种群的适应度函数值并且进入循环，即返回步骤3.1，直到最优解满足终止条件， 或者当迭代的次数超过了预先设置的最大迭代次数K则强制结束，通过BinPSO算法求解min PS/來人)，进而得到最优PSL方向图。
[0057]本发明与现有技术相比，其显著优点在于：（1)与非接触式生命体征监测技术理论 相结合，考虑人体心脏尺寸及具体雷达载频，选取了合适的阵列长度和阵元数目，完成了对 雷达系统的一维及二维阵列的布阵优化设计，获得了理想的PSL; (2)应用范围广泛，同样适 用于一般性虚拟孔径雷达、一般性智能优化算法、一般性二维阵列；（3)能够快速、高效的优 化设计阵元位置和激励分布，得到理想的方向图，满足雷达系统的工程设计需求。
[0058]下面结合附图对本发明做进一步详细说明。
【附图说明】
[0059] 图1为本发明基于非接触式生命体征监测系统的一维雷达布阵优化方法的流程 图。
[0060] 图2为本发明中心脏尺寸示意图。
[0061 ]图3为本发明中PS0类算法的优化流程图。
[0062]图4为本发明中经PS0类算法优化后的一维线阵ΜΙΜΟ雷达阵列拓扑结构图，其中 (a)为发射阵列图，（b)为接收阵列图。
[0063] 图5为本发明经PS0类算法优化后的一维线阵MBTO雷达方向图。
【具体实施方式】
[0064] 结合图1，本发明基于非接触式生命体征监测系统的一维雷达布阵优化方法，基于 非接触式生命体征理论，采用粒子群(Particle Swarm Optimization，PS0)类算法对ΜΙΜΟ 雷达进行一维布阵优化设计，实际上该布阵方法不仅适用于ΜΜ0雷达阵列，对其他虚拟孔 径雷达阵列也适用；该布阵方法不仅适用于PS0类算法，对其他智能优化算法都适用，具体 步骤如下：
[0065] 步骤1，在考虑心脏尺寸的条件下，根据具体的雷达系统载频，确定阵列的孔径与 阵元的数目，对雷达一维阵列孔径进行离散化处理并对离散孔径拓扑使用带约束条件的二 进制向量来表示；
[0066] 步骤1在考虑心脏尺寸的条件下，根据具体的雷达系统载频，确定阵列的孔径与阵 元的数目，对雷达一维阵列孔径进行离散化处理并对离散孔径拓扑使用带约束条件的二进 制向量来表示，具体方法为：
[0067] 步骤1.1，假设雷达的工作频率为f，且人体组织的介电常数为er，光速为c，可得到 雷达工作波长：
[0068] 步骤1.2,在优化雷达阵列前，确定阵列的孔径大小，孔径的大小满足心脏完全被 覆盖，假设人体心脏的尺寸为则设计一个大于Μ的一维线阵，即构造孔径L = aA >M，a为常数系数，M、N的单位均为cm;此处举例假设人体心脏的尺寸大约为12cmXllcm，则 对于一维线阵考虑设id 个大于12cm的一维线阵。即构造孔径L = aAcm>12cm;
[0069] 步骤1.3,对于一维阵列确定阵元数他，此处结合考虑了非接触式生命体征监测雷 达是一种进场天线阵元数目较少的雷达，在考虑实际应用且保证雷达性能的前提下，对他 应适当取值;见=队=队，队、队分别为发射阵元数和接收阵元数；
[0070] 步骤1.4,为方便本发明计算机仿真且保证雷达的一定性能，此处设置步骤1.2中a =7.5。为了减小天线阵元之间的互耦效应，一维线阵的阵元间距满足大于或者等于0.5λ， α2 将线阵按〇. 5λ的间距离散化为叾二V + l个端点，g卩2a+l个端点；此处则将线阵按0.5λ的间距 离散化为16个端点；
[0071] 步骤1.5,为了获得最大的物理孔径，一维线阵的阵元离散化位置应遵循以下规 则:一维线阵发射阵列、接收阵列首尾端点都放置阵元。剩余的他-2个发射阵元与Λ-2个接 收阵元放置在各自剩余的2a+l-2个端点，即2a-l个端点上；使用二进制向量ξ表示MMO雷 达发射阵列拓扑结构，二进制向量g来表示Μ頂0雷达接收阵列拓扑结构，即
[0072]
[0073]
[0074] 此实施例中，剩余的6个发射阵元与6个接收阵元合理放置在各自剩余的14个端点 上。使用二进制向量^表示MMO雷达发射阵列拓扑结构，二进制向量民来表示ΜΜΟ雷达接收 阵列拓朴结构，艮口
[0075]
[0076]
[0077]步骤1.6,对步骤1.5得到的离散孔径拓扑使用带约束条件的二进制向量来表示， 如下所示：
[0078]
[0079]
[0080]
[0081]
[0082]
[0083] 式中，i = l，2，...，2a+l。
[0084] 本实施例表示如下：
[0085]
[0086]
[0087]
[0088]
[0089] 7=1 Z--J
[0090] 步骤2，将步骤1得到的二进制向量代入ΜΜ0雷达方向图公式，确立二进制向量与 ΜΙΜΟ雷达方向图的--对应关系，由ΜΙΜΟ雷达方向图求解峰值旁瓣电平PSL，并将PSL表示 为拓扑结构二进制向量的函数，从而建立二进制优化问题，包括如下步骤：
[0091] 步骤2.1，MBTO雷达的一维方向图公式为：
[0092]
[0093]其中λ表示ΜΠΚ)雷达工作波长，Θ表示平面波与一维阵列法线的夹角，ρ表示方向图 函数，Xtm表示发射阵列的第m个阵元位置，Xrn表示接收阵列的第η个阵元位置；
[0094]步骤2.2，将步骤1.6中带约束条件的二进制向量代入步骤2.1中的ΜΜ0雷达方向 图公式来确立表示拓扑结构的二进制向量与方向图的一一对应关系；
[0095] 步骤2 · 3，可以求解其对应的峰值旁瓣电平PSL(/7,./7,)，将PSL表示为拓扑结构二进 制向量的函数，从而建立对应的二进制优化问题，则一维线阵ΜΜ0雷达优化问题可以表述 为：
[0108]
[0109] 步滕?，对步滕2犾侍的二进制优化问题（目标是降低PSL)采用智能优化算法进行 求解，得到优化后的阵列拓扑结构，由此可以得出对应的最优PSL方向图；
[0110] 步骤3,对步骤2获得的二进制优化问题采用智能优化算法进行求解，得到优化后 的阵列拓扑结构，并确定最优PSL方向图，包括如下步骤：
[0111] 步骤3.1，为解决雷达布阵问题，需求解变量为0或1的二进制优化问题：
[0112]
[0113] subjectto xqe{〇,l},f〇r
[0114]其中，目标函数f〇(Xl,X2, . . .，XQ)是Q个实数变量的函数，并且每个变量Xq只能取值 0或 1，q = l，2,…，Q〇
[0115] 步骤3.2，给定种群规模后，初始化种群中粒子的位置向量弋=1\|,\2，...,' (」，对 于BinPS0，xgq表示第g个粒子的第q个变量，士:中的元素只能取0或者1，通过产生随机的0，1 整数来初始化粒子的位置向量；
[0116] 步骤3.3，对于BinPSO，速度、位置更新公式为：
[0117]
[0118]其中，w为惯性系数，C1是粒子追踪自己历史最优解的权重，C2是粒子追踪群体历史 最优解的权重，ξ与η是均匀分布在[0，1]的随机数，<+1表示第i个粒子第k+l次迭代的速度 的第q维，< 表示第i个粒子第k次迭代的速度的第q维，Pt表示第i个粒子第k次迭代后搜索 到的历史最优位置的第q维，表示第i个粒子第k次迭代后位置的第q维，;表示粒子第 k次迭代后搜索到的全局最优位置的第q维，粒子群算法的最大迭代次数设定为K，当迭代次 数超过K时，粒子群算法结束。
[0119]步骤3.4，步骤3.3中速度的选定存在限定条件：
[0120]
[0121]表明若更新后的速度大于设定的最大值或者小于设定的最小值，那么便强制其约 束在最大值或者最小值；
[0122] 步骤3.5，位置的迭代公式为：.<+1 = .< +<+1，其中，<+1表示第i个粒子第k+Ι次迭 代的位置的第q维，·<表示第i个粒子第k次迭代的位置的第q维。对于BinPSO算法位置的选 定存在如下限定条件：
[0123]
[0124] 其中，rand表示均勾分布在[0,1]的随机数，S()表示sigmoid函数，即 -外）=，为了防止sigmoid函数的饱和性，粒子的速度范围[Vqmin，Vqmax]建议为[-4，4]; l + e 1表示放置阵元而0表示不放置阵元；
[0125] 步骤3.6,按照步骤3.3给出的速度、位置更新公式以及步骤3.5给出的位置的迭代 公式更新粒子的速度与位置；
[0126] 步骤3.7，确定最优解并且判断是否满足终止条件，如果满足即算法完成，否则计 算新一代种群的适应度函数值并且进入循环，即返回步骤3.1，直到最优解满足终止条件， 或者当迭代的次数超过了预先设置的最大迭代次数K则强制结束，通过BinPSO算法求解min PSZ,进而得到最优PSL方向图。其中，适应度函数为步骤2.3中的峰值旁瓣电平 PSL(/7,A ),：
[0127] BinPSO参数设置如下：
[0128] 种群规模:对于一维阵列见=30;对于二维阵列N2 = 100
[0129] 最大迭代次数:对于一维阵列1^ = 30;对于二维阵列K2 = 50
[0130] 粒子速度范围：[Vqmin，Vqmax ] = [ -4，4 ]
[01 31 ] 惯性系数:wmax = 0.95，Wmin = 0.4，按照迭代次数线性递减
[0132] 粒子追踪自己历史最优解的权重:cl = 1.4
[0133] 粒子追踪群体历史最优解的权重:c2 = 1.4
[0134] 通过以上步骤最后就能得到经过PS0类算法优化后的最优PSL。
[0135] 结合图4、图5:图4(a)是一维阵列ΜΜ0雷达BinPSO优化后的发射阵列拓扑结构，图 4(b)是一维阵列ΜΜ0雷达BinPSO优化后的接收阵列拓扑结构，图5是由BinPSO优化后一维 阵列Μ頂0雷达的方向图，其峰值旁瓣电平为-23.04dB。
[0136] 综上所述，本发明通过智能优化算法能够有效降低非接触式生命体征监测系统雷 达的PSL，从而实现非接触式生命体征监测系统的一维雷达布阵优化设计。
【主权项】
1. 一种基于非接触式生命体征监测系统的一维雷达布阵优化方法，其特征在于，步骤 如下： 步骤1，在考虑屯、脏尺寸的条件下，根据具体的雷达系统载频，确定阵列的孔径与阵元 的数目，对雷达一维阵列孔径进行离散化处理并对离散孔径拓扑使用带约束条件的二进制 向量来表示； 步骤2,将步骤1得到的二进制向量代入ΜΙΜΟ雷达方向图公式，确立二进制向量与ΜΙΜΟ 雷达方向图的一一对应关系，由ΜΙΜΟ雷达方向图求解峰值旁瓣电平PSL，并将Ρ化表示为拓 扑结构二进制向量的函数，从而建立二进制优化问题； 步骤3,对步骤2获得的二进制优化问题采用智能优化算法进行求解，得到优化后的阵 列拓扑结构，并确定最优PSL方向图。2. 根据权利要求书1所述的基于非接触式生命体征监测系统的一维雷达布阵优化方 法，其特征在于，步骤1所述在考虑屯、脏尺寸的条件下，根据具体的雷达系统载频，确定阵列 的孔径与阵元的数目，对雷达一维阵列孔径进行离散化处理并对离散孔径拓扑使用带约束 条件的二进制向量来表示，具体步骤为： 步骤1.1，假设雷达的工作频率为f，人体组织的介电常数为Er，光速为C，得到雷达工作 波去步骤1.2,确定阵列的孔径大小，孔径的大小满足屯、脏完全被覆盖，假设人体屯、脏的尺 寸为MXN且M>N，则设计一个大于Μ的一维线阵，即构造孔径L = aA>M，a为常数系数，M、N的 单位均为cm; 步骤1.3,确定一维线阵的阵元数化，化= Nt = Nr，Nt、Nr分别为发射阵元数和接收阵元 数； 步骤1.4,一维线阵的阵元间距满足大于或者等于0.5λ，将线阵按〇.5λ的间距离散化为 +1个端点，即2a+1个端点； 步骤1.5, 一维线阵的阵元离散化位置应遵循W下规则：一维线阵发射阵列、接收阵列 首尾端点都放置阵元，剩余的化-2个发射阵元与化-2个接收阵元放置在各自剩余的2a+l-2 个端点，即2a-1个端点上;使用二进制向量ζ表示ΜΙΜΟ雷达发射阵列拓扑结构，二进制向量 ζ来表示ΜΙΜΟ雷达接收阵列拓扑结构，即步骤1.6,对步骤1.5得到的离散孔径拓扑使用带约束条件的二进制向量来表示，如下 所示：bti , bri Ε { 0,11 btl 二 bt2a+l 二 brl 二 br2a+l 二 1式中，i = l,2,···, 2a+l。3. 根据权利要求书1所述的基于非接触式生命体征监测系统的一维雷达布阵优化方 法，其特征在于，步骤2所述将步骤1得到的二进制向量代入ΜΙΜΟ雷达方向图公式，确立二进 制向量与ΜΙΜΟ雷达方向图的一一对应关系，由ΜΙΜΟ雷达方向图求解峰值旁瓣电平PSL，并将 Ρ化表示为拓扑结构二进制向量的函数，从而建立二进制优化问题，具体步骤如下： 步骤2.1 ,ΜΙΜΟ雷达的一维方向图公式为：其中λ表示ΜΙΜΟ雷达工作波长，Θ表示平面波与一维阵列法线的夹角，Ρ表示方向图函 数，表示发射阵列的第m个阵元位置，Xrn表示接收阵列的第η个阵元位置； 步骤2.2，将带约束条件的二进制向量代入步骤2.1中的ΜΙΜΟ雷达方向图公式来表示拓 扑结构的二进制向量与方向图的一一对应关系； 步骤2.3，求解峰值旁瓣电平巧L话，衣)，将Ρ化表示为拓扑结构二进制向量的函数，从 而建立二进制优化问题，则一维线阵ΜΙΜΟ雷达优化问题表述为：式中，1 = 1,2,···, 2a+l。4. 根据权利要求书1所述的基于非接触式生命体征监测系统的一维雷达布阵优化方 法，其特征在于，步骤3所述对步骤2获得的二进制优化问题采用智能优化算法进行求解，得 到优化后的阵列拓扑结构，并确定最优PSL方向图，具体步骤如下： 步骤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，给定种群规模后，初始化种群中粒子的位置向量^=(^^1，乂^2，...，^0>对于 BinPSO，Xgq表示第g个粒子的第q个变量，馬中的元素只能取0或者1，通过产生随机的0、1整 数来初始化粒子的位置向量； 步骤3.3，对于BinPSO，速度、位置更新公式为：其中，w为惯性系数，Cl是粒子追踪自己历史最优解的权重，C2是粒子追踪群体历史最优 解的权重，ξ与η是均匀分布在[〇，1]的随机数，<+1表示第i个粒子第k+1次迭代的速度的第q 维，4表示第i个粒子第k次迭代的速度的第q维，表示第i个粒子第k次迭代后捜索到的 历史最优位置的第q维，4表示第i个粒子第k次迭代后位置的第q维，表示粒子第k次 迭代后捜索到的全局最优位置的第q维，粒子群算法的最大迭代次数设定为K，当迭代次数 超过即寸，粒子群算法结束； 步骤3.4，步骤3.3中速度的选定存在限定条件：若更新后的速度大于设定的最大值便强制约束在速度最大值Vqmax，或者小于设定的最 小值便强制约束在速度最小值Vqmin; 步骤3.5，位置的迭代公式为：其中，表示第i个粒子第k+1次迭代的位置的第q维，表示第i个粒子第k次迭代的 位置的第q维； 对于BinPSO算法位置的选定存在如下限定条件：其中，rand表示均匀分布在[0, υ的随机数，S〇表示sigmoid函数，印，为了 防止sigmoid函数的饱和性，粒子的速度范围[Vqmin，Vqmax]为[-4,4];1表示放置阵元而0表示 不放置阵元； 步骤3.6，按照步骤3.3给出的速度、位置更新公式W及步骤3.5给出的位置的迭代公式 更新粒子的速度与位置； 步骤3.7，确定最优解并且判断是否满足终止条件，如果满足即算法完成，否则计算新 一代种群的适应度函数值并且进入循环，即返回步骤3.1，直到最优解满足终止条件，或者 当迭代的次数超过了预先设置的最大迭代次数K则强制结束，通过BinPSO算法求解min 尸57. K，6.)，进而得到最优P化方向图。
【文档编号】G01S13/88GK105974404SQ201610270185
【公开日】2016年9月28日
【申请日】2016年4月27日
【发明人】顾陈, 王芳, 李聪羚, 杨宇鑫, 洪弘, 李彧晟, 朱晓华
【申请人】南京理工大学