CN112287564A - Electrode array optimization method based on goblet sea squirt group algorithm - Google Patents

Electrode array optimization method based on goblet sea squirt group algorithm Download PDF

Info

Publication number
CN112287564A
CN112287564A CN202011311093.1A CN202011311093A CN112287564A CN 112287564 A CN112287564 A CN 112287564A CN 202011311093 A CN202011311093 A CN 202011311093A CN 112287564 A CN112287564 A CN 112287564A
Authority
CN
China
Prior art keywords
electrode array
electrode
salps
salp
space
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
Application number
CN202011311093.1A
Other languages
Chinese (zh)
Other versions
CN112287564B (en
Inventor
蒋毅
许源
梁运华
刘功能
郭奉仁
李泽文
王梓糠
尹骏刚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
State Grid Hunan Electric Power Co Ltd
Maintenance Co of State Grid Hunan Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Hunan Electric Power Co Ltd
Maintenance Co of State Grid Hunan Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, State Grid Hunan Electric Power Co Ltd, Maintenance Co of State Grid Hunan Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202011311093.1A priority Critical patent/CN112287564B/en
Publication of CN112287564A publication Critical patent/CN112287564A/en
Application granted granted Critical
Publication of CN112287564B publication Critical patent/CN112287564B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial 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]

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Data Mining & Analysis (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Computational Linguistics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses an electrode array optimization method based on a goblet sea squirt group algorithm, which comprises the steps of taking the set number of electrodes as the dimension of a search space of the goblet sea squirt group, taking the arrangement space of the electrode array as the set search space of the goblet sea squirt group, utilizing iteration of the goblet sea squirt group algorithm until the set target is met or the iteration reaches the set times, outputting the final food position of the goblet sea squirt group, and taking the final food position as the arrangement position of the electrode array to complete optimization. The electrode array optimization method based on the goblet sea squirt group algorithm has the advantages of simplicity, practicability, small calculated amount, high optimization speed, high efficiency and the like.

Description

基于樽海鞘群算法的电极阵列优化方法Electrode array optimization method based on salps swarm algorithm

技术领域technical field

本发明涉及电极检测技术领域,尤其涉及一种基于樽海鞘群算法的电极阵列优化方法。The invention relates to the technical field of electrode detection, in particular to an electrode array optimization method based on a salps swarm algorithm.

背景技术Background technique

电力行业的检修、预试等人员在高压场地内工作,经常遇到作业范围周边设备带电运行的情形,存在安全距离不足的隐患。普通近电告警装置难以在复杂电磁环境中准确报警。基于多电极的近电报警装置具有电场信号空间分辨率高和抗干扰能力强等优势,可通过电极阵列定向检测设备是否带电,从而解决普通近电报警装置误报率较高等问题。The maintenance, pre-test and other personnel in the power industry work in high-voltage fields, and often encounter the situation that the surrounding equipment is running under power, and there is a hidden danger of insufficient safety distance. It is difficult for ordinary near-electrical alarm devices to accurately alarm in complex electromagnetic environments. The multi-electrode-based near-electrical alarm device has the advantages of high spatial resolution of electric field signals and strong anti-interference ability. It can detect whether the equipment is charged through the electrode array, so as to solve the problem of high false alarm rate of ordinary near-electrical alarm devices.

多电极空间分布和激励方式是该近电报警装置的技术关键,宜采用智能算法对其进行优化设计。智能算法被广泛应用在工程领域中,用于最优化问题特别是非线性优化问题的求解。近年来,智能算法因其在解决非线性优化问题上的效率优势,被应用于各类传感器阵列设计中,比如遗传算法、粒子群算法等都已得到了广泛的应用。然而,一个复杂电极阵列模型的优化过程可能需要几小时甚至更长时间,计算代价极其巨大。The multi-electrode spatial distribution and excitation mode are the technical keys of the near-electrical alarm device, and an intelligent algorithm should be used to optimize the design. Intelligent algorithms are widely used in engineering to solve optimization problems, especially nonlinear optimization problems. In recent years, intelligent algorithms have been widely used in the design of various sensor arrays due to their efficiency advantages in solving nonlinear optimization problems, such as genetic algorithms and particle swarm optimization. However, the optimization process of a complex electrode array model can take hours or even longer, and the computational cost is extremely high.

发明内容SUMMARY OF THE INVENTION

本发明要解决的技术问题是克服现有技术的不足,提供一种简单实用、计算量较小、优化速度快和效率高的基于樽海鞘群算法的电极阵列优化方法。The technical problem to be solved by the present invention is to overcome the deficiencies of the prior art, and to provide a simple and practical electrode array optimization method based on the salps swarm algorithm, which is simple and practical, has a small amount of calculation, and has fast optimization speed and high efficiency.

为解决上述技术问题,本发明提出的技术方案为:In order to solve the above-mentioned technical problems, the technical scheme proposed by the present invention is:

一种基于樽海鞘群算法的电极阵列优化方法,以设定电极数作为樽海鞘群搜索空间维度,以电极阵列排布空间作为樽海鞘群的设定搜索空间,利用樽海鞘群算法迭代直至满足设定目标或迭代达到设定次数,输出樽海鞘群最终食物位置,并将其作为电极阵列的排列位置,完成优化。An electrode array optimization method based on the salp swarm algorithm, the set number of electrodes is used as the dimension of the salp group search space, the electrode array arrangement space is used as the set search space of the salp group, and the salp swarm algorithm is used to iterate until it satisfies Set a target or iterate to a set number of times, output the final food position of the salps group, and use it as the arrangement position of the electrode array to complete the optimization.

利用樽海鞘群算法迭代直至满足设定目标的过程包括以下步骤:The process of iterating using the salps swarm algorithm until the set goal is met includes the following steps:

S1:初始化樽海鞘群的个体个数N=3n+1,n为自然数,将电极数作为搜索空间的维度数D,设定迭代次数Tmax;初始化种群为

Figure BDA0002789862080000011
S1: Initialize the number of individuals of the salps group N=3n+1, n is a natural number, take the number of electrodes as the dimension D of the search space, and set the number of iterations T max ; the initialization population is
Figure BDA0002789862080000011

S2:计算樽海鞘群各个体的适应度,将适应度最大的个体作为食物,其余的个体均分为n个领航者、n个追踪者以及n个后援者;S2: Calculate the fitness of each individual in the sea squirt group, take the individual with the greatest fitness as the food, and divide the rest of the individuals into n leaders, n trackers, and n supporters;

S3:更新领航者位置为:S3: Update the navigator position to:

Figure BDA0002789862080000021
Figure BDA0002789862080000021

其中,

Figure BDA0002789862080000022
为第i只领航者樽海鞘个体在第j维空间的位置,ubj与lbj分别为在第j维搜索空间的搜索上限与搜索下限;Fj为食物在第j维空间的位置,c2与c3为区间[0,1]内产生的随机数;参数c1的定义为:
Figure BDA0002789862080000023
t为当前迭代次数,Tmax为最大迭代次数;in,
Figure BDA0002789862080000022
is the position of the i-th pilot salp individual in the j-th dimension space, ub j and lb j are the search upper limit and search limit in the j-th dimension search space, respectively; F j is the position of the food in the j-th dimension space, c 2 and c 3 are random numbers generated in the interval [0,1]; the definition of parameter c 1 is:
Figure BDA0002789862080000023
t is the current number of iterations, and Tmax is the maximum number of iterations;

S4:更新追踪者位置为:S4: Update tracker location to:

Figure BDA0002789862080000024
Figure BDA0002789862080000024

式中,

Figure BDA0002789862080000025
为第i只追踪者樽海鞘个体在第j维空间的位置,
Figure BDA0002789862080000026
为共生量;B为受益参数,随机选取1或2;r为[0,1]区间的随机数;In the formula,
Figure BDA0002789862080000025
is the position of the i-th tracker salp individual in the j-th dimension space,
Figure BDA0002789862080000026
is the symbiotic amount; B is the benefit parameter, 1 or 2 is randomly selected; r is a random number in the [0,1] interval;

S5:后援者位置更新为:S5: Supporter location updated to:

Figure BDA0002789862080000027
Figure BDA0002789862080000027

Figure BDA0002789862080000028
Figure BDA0002789862080000028

式中,

Figure BDA0002789862080000029
为第i只后援者樽海鞘个体在第j维空间的位置,t为当前迭代次数;Tmax为最大迭代次数;In the formula,
Figure BDA0002789862080000029
is the position of the i-th supporter salps in the j-th dimension space, t is the current number of iterations; T max is the maximum number of iterations;

S6:重复步骤S2~S5,直至迭代次数达到Tmax或食物适应度值达到设定终止阈值。S6: Repeat steps S2 to S5 until the number of iterations reaches T max or the food fitness value reaches the set termination threshold.

更新第i只追踪者樽海鞘个体在第j维空间的位置时,存储更新前的第i只追踪者樽海鞘个体在第j维空间的位置为

Figure BDA00027898620800000210
并判定
Figure BDA00027898620800000211
Figure BDA00027898620800000212
的适应度优劣,当
Figure BDA00027898620800000213
优于
Figure BDA00027898620800000214
时,再次更新第i只追踪者樽海鞘个体在第j维空间的位置为
Figure BDA00027898620800000215
When updating the position of the i-th tracker salp individual in the j-th dimension space, the position of the i-th tracker-salp individual in the j-th dimension space before the update is stored as
Figure BDA00027898620800000210
and judge
Figure BDA00027898620800000211
and
Figure BDA00027898620800000212
The fitness is good or bad, when
Figure BDA00027898620800000213
better than
Figure BDA00027898620800000214
When , the position of the i-th tracker salps in the j-th dimension space is updated again as
Figure BDA00027898620800000215

樽海鞘群个体还包括多个后援者,后援者位置更新为:The individual salps group also includes multiple supporters, and the positions of the supporters are updated to:

Figure BDA00027898620800000216
Figure BDA00027898620800000216

Figure BDA00027898620800000217
Figure BDA00027898620800000217

式中,

Figure BDA0002789862080000031
为第i只后援者樽海鞘个体在第j维空间的位置,t为当前迭代次数;Tmax为最大迭代次数。In the formula,
Figure BDA0002789862080000031
is the position of the i-th supporter salps in the j-th dimensional space, t is the current number of iterations; T max is the maximum number of iterations.

目标电极阵列排布空间为同心圆电极阵列,设定目标为最小化峰值旁瓣电平,各所述樽海鞘个体的适应度如下运算得到:The arrangement space of the target electrode array is a concentric electrode array, and the target is set to minimize the peak side lobe level. The fitness of each individual salps is calculated as follows:

Figure BDA0002789862080000032
Figure BDA0002789862080000032

Figure BDA0002789862080000033
Figure BDA0002789862080000033

其中,θ为偏移法线的角度,

Figure BDA0002789862080000034
为电极子的相移量,j是虚数符号,Fmax为同心圆电极阵列的主瓣峰值,Ny为电极排布的圆环数量,Nx为环上电极子数量,rm指第m个圆环上第n个电极子(n,m)的半径,
Figure BDA0002789862080000035
馈电幅度是I(n,m)各电极子同激励,Ψ(n,m)是电极子与x轴的夹角。where θ is the angle offset from the normal,
Figure BDA0002789862080000034
is the phase shift of the electrode, j is the sign of the imaginary number, F max is the peak value of the main lobe of the concentric electrode array, N y is the number of electrode rings, N x is the number of electrodes on the ring, r m refers to the mth The radius of the nth electrode element (n,m) on the ring,
Figure BDA0002789862080000035
The feeding amplitude is I(n,m) the co-excitation of each electrode, and Ψ(n,m) is the angle between the electrode and the x-axis.

目标电极阵列排布空间为矩形电极阵列,设定目标为最小化峰值旁瓣电平,各所述樽海鞘个体的适应度如下运算得到:The arrangement space of the target electrode array is a rectangular electrode array, and the target is set to minimize the peak side lobe level. The fitness of each individual salps is calculated as follows:

Figure BDA0002789862080000036
Figure BDA0002789862080000036

Figure BDA0002789862080000037
Figure BDA0002789862080000037

其中,θ为偏移法线的角度,

Figure BDA0002789862080000038
为电极子的相移量,j是虚数符号,Fmax为同心圆电极阵列的主瓣峰值,Ny为电极排布的圆环数量,Nx为环上电极子数量,rm指第m个圆环上第n个电极子(n,m)的半径,
Figure BDA0002789862080000039
pmx为电极在x方向上的位置,pmy是电极子在y方向上的位置,馈电幅度是I(n,m)各电极子同激励。where θ is the angle offset from the normal,
Figure BDA0002789862080000038
is the phase shift of the electrode, j is the sign of the imaginary number, F max is the peak value of the main lobe of the concentric electrode array, N y is the number of electrode rings, N x is the number of electrodes on the ring, r m refers to the mth The radius of the nth electrode element (n,m) on the ring,
Figure BDA0002789862080000039
p mx is the position of the electrode in the x direction, p my is the position of the electrode element in the y direction, and the feeding amplitude is I(n, m) where each electrode element is co-excited.

与现有技术相比,本发明的优点在于:Compared with the prior art, the advantages of the present invention are:

本发明的基于樽海鞘群算法的电极阵列优化方法,以设定电极数作为樽海鞘群搜索空间维度,以电极阵列排布空间作为樽海鞘群的设定搜索空间,利用樽海鞘群算法迭代直至满足设定目标或迭代达到设定次数,输出樽海鞘群最终食物位置,并将其作为电极阵列的排列位置,完成优化。这种优化方式步骤简单,容易实现,并且求解精度能够满足需求,收敛速度较快,探索范围不会陷入局部最优解中,因此大大降低了所需的计算量,缩短了优化过程所需时间,提升电极阵列性能的同时节约了时间与成本。In the electrode array optimization method based on the salp group algorithm of the present invention, the set number of electrodes is used as the dimension of the search space of the salp group, the electrode array arrangement space is used as the set search space of the salp group, and the salp group algorithm is used to iterate until When the set target is met or the iteration reaches the set number of times, the final food position of the salps group is output, and it is used as the arrangement position of the electrode array to complete the optimization. This optimization method is simple and easy to implement, and the solution accuracy can meet the requirements, the convergence speed is fast, and the exploration range will not fall into the local optimal solution, so the required calculation amount is greatly reduced, and the time required for the optimization process is shortened. , which saves time and cost while improving the performance of the electrode array.

附图说明Description of drawings

图1为本发明基于樽海鞘群算法的电极阵列优化方法得到的同心圆电极阵列图;1 is a diagram of a concentric electrode array obtained by an electrode array optimization method based on the salp group algorithm of the present invention;

图2为图1所示的电极阵列得到的检测图。FIG. 2 is a detection diagram obtained by the electrode array shown in FIG. 1 .

具体实施方式Detailed ways

为了便于理解本发明,下文将结合说明书附图和较佳的实施例对本文发明做更全面、细致地描述,但本发明的保护范围并不限于以下具体实施例。In order to facilitate the understanding of the present invention, the present invention will be described more comprehensively and in detail below with reference to the accompanying drawings and preferred embodiments of the specification, but the protection scope of the present invention is not limited to the following specific embodiments.

实施例:Example:

本实施例的基于樽海鞘群算法的电极阵列优化方法,以设定电极数作为樽海鞘群搜索空间维度,以电极阵列排布空间作为樽海鞘群的设定搜索空间,利用樽海鞘群算法迭代直至满足设定目标或迭代达到设定次数,输出樽海鞘群最终食物位置,并将其作为电极阵列的排列位置,完成优化。这种优化方式步骤简单,容易实现,并且求解精度能够满足需求,收敛速度较快,探索范围不会陷入局部最优解中,因此大大降低了所需的计算量,缩短了优化过程所需时间,提升电极阵列性能的同时节约了时间与成本。In the electrode array optimization method based on the salp swarm algorithm in this embodiment, the number of electrodes is set as the dimension of the salp group search space, the electrode array arrangement space is used as the set search space of the salp group, and the salp swarm algorithm is used to iterate Until the set target is met or the iteration reaches the set number of times, the final food position of the salps group is output, and it is used as the arrangement position of the electrode array to complete the optimization. This optimization method is simple and easy to implement, and the solution accuracy can meet the requirements, the convergence speed is fast, and the exploration range will not fall into the local optimal solution, so the required calculation amount is greatly reduced, and the time required for the optimization process is shortened. , which saves time and cost while improving the performance of the electrode array.

本实施例中,对同心圆电极阵列进行优化,约束方向图主瓣区域为

Figure BDA0002789862080000041
设定目标为最小化峰值旁瓣电平,初始目标值设定为-15dB,若在达到最大迭代数之前满足设定目标值,则目标值降低0.5dB。In this embodiment, the concentric circular electrode array is optimized, and the main lobe region of the constraint pattern is
Figure BDA0002789862080000041
The set goal is to minimize the peak sidelobe level. The initial target value is set to -15dB. If the set target value is met before the maximum number of iterations is reached, the target value is reduced by 0.5dB.

利用樽海鞘群算法迭代直至满足设定目标的过程包括以下步骤:The process of iterating using the salps swarm algorithm until the set goal is met includes the following steps:

S1:初始化樽海鞘群的个体个数N=3n+1,n为自然数,N一般取值范围在50~100之间,本实施例为61,将电极数作为搜索空间的维度数D,本实施例为24,设定迭代次数Tmax=500;初始化种群为

Figure BDA0002789862080000042
初始化种群个体中各电极坐标
Figure BDA0002789862080000043
其中θ为[0,2π]之间的随机数,多个电极沿三个同心圆环阵列排布,且同一环上,相邻电极子距离应大于设定距离。S1: Initialize the number of individuals in the sea squirt group N=3n+1, n is a natural number, and N generally ranges from 50 to 100. In this embodiment, it is 61. The number of electrodes is used as the dimension D of the search space. The embodiment is 24, and the number of iterations T max =500 is set; the initialization population is
Figure BDA0002789862080000042
Initialize the coordinates of each electrode in the individual population
Figure BDA0002789862080000043
Among them, θ is a random number between [0, 2π]. Multiple electrodes are arranged in an array of three concentric rings, and on the same ring, the distance between adjacent electrodes should be greater than the set distance.

S2:计算樽海鞘群各个体的适应度:S2: Calculate the fitness of each individual in the salps group:

Figure BDA0002789862080000044
Figure BDA0002789862080000044

Figure BDA0002789862080000045
Figure BDA0002789862080000045

其中,θ为偏移法线的角度,

Figure BDA0002789862080000051
为电极子的相移量,j是虚数符号,Fmax为同心圆电极阵列的主瓣峰值,Ny为电极排布的圆环数量,Nx为环上电极子数量,rm指第m个圆环上第n个电极子(n,m)的半径,
Figure BDA0002789862080000052
馈电幅度是I(n,m)各电极子同激励,Ψ(n,m)是电极子与x轴的夹角;where θ is the angle offset from the normal,
Figure BDA0002789862080000051
is the phase shift of the electrode, j is the sign of the imaginary number, F max is the peak value of the main lobe of the concentric electrode array, N y is the number of electrode rings, N x is the number of electrodes on the ring, r m refers to the mth The radius of the nth electrode element (n,m) on the ring,
Figure BDA0002789862080000052
The feeding amplitude is I(n,m) the co-excitation of each electrode, and Ψ(n,m) is the angle between the electrode and the x-axis;

将其中适应度最大的个体作为食物,其余的个体均分为n个领航者、n个追踪者以及n个后援者;The individual with the greatest fitness is used as food, and the rest are divided into n leaders, n trackers, and n supporters;

S3:更新领航者位置为:S3: Update the navigator position to:

Figure BDA0002789862080000053
Figure BDA0002789862080000053

其中,

Figure BDA0002789862080000054
为第i只领航者樽海鞘个体在第j维空间的位置,ubj与lbj分别为在第j维搜索空间的搜索上限与搜索下限,搜索空间坐标可依据
Figure BDA0002789862080000055
求得,求取上限时,m=3,求取下限时m=1;Fj为食物在第j维空间的位置,c2与c3为区间[0,1]内产生的随机数;参数c1的定义为:
Figure BDA0002789862080000056
t为当前迭代次数,Tmax为最大迭代次数;in,
Figure BDA0002789862080000054
is the position of the i-th pilot salp individual in the j-th dimension space, ub j and lb j are the search upper limit and search limit in the j-th dimension search space, respectively, and the search space coordinates can be based on
Figure BDA0002789862080000055
Obtain, when obtaining the upper limit, m=3, when obtaining the lower limit, m=1; F j is the position of the food in the jth dimension space, and c 2 and c 3 are random numbers generated in the interval [0,1]; The parameter c 1 is defined as:
Figure BDA0002789862080000056
t is the current number of iterations, and Tmax is the maximum number of iterations;

S4:更新追踪者位置为:S4: Update tracker location to:

Figure BDA0002789862080000057
Figure BDA0002789862080000057

式中,

Figure BDA0002789862080000058
为第i只追踪者樽海鞘个体在第j维空间的位置,
Figure BDA0002789862080000059
为共生量;B为受益参数,随机选取1或2;r为[0,1]区间的随机数;更新第i只追踪者樽海鞘个体在第j维空间的位置时,存储更新前的第i只追踪者樽海鞘个体在第j维空间的位置为
Figure BDA00027898620800000510
并判定
Figure BDA00027898620800000511
Figure BDA00027898620800000512
的适应度优劣,当
Figure BDA00027898620800000513
优于
Figure BDA00027898620800000514
时,再次更新第i只追踪者樽海鞘个体在第j维空间的位置为
Figure BDA00027898620800000515
In the formula,
Figure BDA0002789862080000058
is the position of the i-th tracker salp individual in the j-th dimension space,
Figure BDA0002789862080000059
is the symbiotic amount; B is the benefit parameter, 1 or 2 is randomly selected; r is a random number in the interval [0,1]; when updating the position of the i-th tracker salps in the j-th dimension space, store the pre-update The position of the individual i tracker salps in the jth dimension is:
Figure BDA00027898620800000510
and judge
Figure BDA00027898620800000511
and
Figure BDA00027898620800000512
The fitness is good or bad, when
Figure BDA00027898620800000513
better than
Figure BDA00027898620800000514
When , the position of the i-th tracker salps in the j-th dimension space is updated again as
Figure BDA00027898620800000515

S5:后援者位置更新为:S5: Supporter location updated to:

Figure BDA00027898620800000516
Figure BDA00027898620800000516

Figure BDA0002789862080000061
Figure BDA0002789862080000061

式中,

Figure BDA0002789862080000062
为第i只后援者樽海鞘个体在第j维空间的位置,t为当前迭代次数;Tmax为最大迭代次数。通过后援者引入余弦函数使后援者更新方式更具灵活性,加上1使位置更新距离总体提高,又由于余弦函数的最大值为1的缘故,不致使距离过于大,从而,使后援者群体全局搜索和局部搜索得到平衡。领航者、跟随者、后援者三个子种群有利于在一定程度上提升算法的探索性能以及后期跳出局部最优的能力。In the formula,
Figure BDA0002789862080000062
is the position of the i-th supporter salps in the j-th dimensional space, t is the current number of iterations; T max is the maximum number of iterations. The introduction of the cosine function by the supporters makes the update method of the supporters more flexible, and the addition of 1 increases the overall position update distance, and because the maximum value of the cosine function is 1, the distance is not too large, so that the supporters group Global search and local search are balanced. The three sub-populations of leader, follower, and supporter are beneficial to improve the exploration performance of the algorithm to a certain extent and the ability to jump out of the local optimum in the later stage.

S6:重复步骤S2~S5,直至迭代次数达到Tmax或食物适应度值达到设定终止阈值。S6: Repeat steps S2 to S5 until the number of iterations reaches T max or the food fitness value reaches the set termination threshold.

之后即可以输出食物位置作为电极排布位置,按此位置排布电极可以使电极阵列达到目标要求。图1即为根据上述优化得到的最优电极阵列方向图,根据其得到图2中所示的检测图,峰值旁瓣电平为-17.3dB,且具有较窄主瓣。After that, the food position can be output as the electrode arrangement position, and arranging electrodes according to this position can make the electrode array meet the target requirements. Figure 1 is the optimal electrode array pattern obtained according to the above optimization. According to it, the detection diagram shown in Figure 2 is obtained. The peak sidelobe level is -17.3dB and has a narrow mainlobe.

本实施例中,当目标电极阵列排布空间为矩形电极阵列时,设定目标为最小化峰值旁瓣电平,各所述樽海鞘个体的适应度如下运算得到:In this embodiment, when the target electrode array arrangement space is a rectangular electrode array, the target is set to minimize the peak side lobe level, and the fitness of each individual salps is calculated as follows:

Figure BDA0002789862080000063
Figure BDA0002789862080000063

Figure BDA0002789862080000064
Figure BDA0002789862080000064

其中,θ为偏移法线的角度,

Figure BDA0002789862080000065
为电极子的相移量,j是虚数符号,Fmax为同心圆电极阵列的主瓣峰值,Ny为电极排布的圆环数量,Nx为环上电极子数量,rm指第m个圆环上第n个电极子(n,m)的半径,
Figure BDA0002789862080000066
pmx为电极在x方向上的位置,pmy是电极子在y方向上的位置,馈电幅度是I(n,m)各电极子同激励。where θ is the angle offset from the normal,
Figure BDA0002789862080000065
is the phase shift of the electrodes, j is the sign of the imaginary number, F max is the peak value of the main lobe of the concentric electrode array, N y is the number of electrode rings, N x is the number of electrodes on the ring, r m refers to the mth The radius of the nth electrode element (n,m) on the ring,
Figure BDA0002789862080000066
p mx is the position of the electrode in the x direction, p my is the position of the electrode element in the y direction, and the feeding amplitude is I(n, m) where each electrode element is co-excited.

以上所述仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例。对于本领域的技术人员来说,在不脱离本发明的技术构思前提下所得到的改进和变换也应视为本发明的保护范围。The above descriptions are only preferred embodiments of the present invention, and the protection scope of the present invention is not limited to the above-mentioned embodiments. For those skilled in the art, improvements and transformations obtained without departing from the technical concept of the present invention should also be regarded as the protection scope of the present invention.

Claims (5)

1.一种基于樽海鞘群算法的电极阵列优化方法,其特征在于:以设定电极数作为樽海鞘群搜索空间维度,以电极阵列排布空间作为樽海鞘群的设定搜索空间,利用樽海鞘群算法迭代直至满足设定目标或迭代达到设定次数,输出樽海鞘群最终食物位置,并将其作为电极阵列的排列位置,完成优化。1. an electrode array optimization method based on the salp group algorithm, is characterized in that: with the set electrode number as the salp group search space dimension, with the electrode array arrangement space as the set search space of the salp group, using the salp group The ascidian swarm algorithm iterates until the set target is met or the iteration reaches the set number of times, and the final food position of the salp swarm is output, which is used as the arrangement position of the electrode array to complete the optimization. 2.根据权利要求1所述的基于樽海鞘群算法的电极阵列优化方法,其特征在于:利用樽海鞘群算法迭代直至满足设定目标的过程包括以下步骤:2. The electrode array optimization method based on the salps swarm algorithm according to claim 1, is characterized in that: the process that utilizes salps swarm algorithm iteration until meeting the set target comprises the following steps: S1:初始化樽海鞘群的个体个数N=3n+1,n为自然数,将电极数作为搜索空间的维度数D,设定迭代次数Tmax;初始化种群为
Figure FDA0002789862070000011
S1: Initialize the number of individuals of the salps group N=3n+1, n is a natural number, take the number of electrodes as the dimension D of the search space, and set the number of iterations T max ; the initialization population is
Figure FDA0002789862070000011
S2:计算樽海鞘群各个体的适应度,将适应度最大的个体作为食物,其余的个体均分为n个领航者、n个追踪者以及n个后援者;S2: Calculate the fitness of each individual in the sea squirt group, take the individual with the greatest fitness as the food, and divide the rest of the individuals into n leaders, n trackers, and n supporters; S3:更新领航者位置为:S3: Update the navigator position to:
Figure FDA0002789862070000012
Figure FDA0002789862070000012
其中,
Figure FDA0002789862070000013
为第i只领航者樽海鞘个体在第j维空间的位置,ubj与lbj分别为在第j维搜索空间的搜索上限与搜索下限;Fj为食物在第j维空间的位置,c2与c3为区间[0,1]内产生的随机数;参数c1的定义为:
Figure FDA0002789862070000014
t为当前迭代次数,Tmax为最大迭代次数;
in,
Figure FDA0002789862070000013
is the position of the i-th pilot salp individual in the j-th dimension space, ub j and lb j are the search upper limit and search limit in the j-th dimension search space, respectively; F j is the position of the food in the j-th dimension space, c 2 and c 3 are random numbers generated in the interval [0,1]; the definition of parameter c 1 is:
Figure FDA0002789862070000014
t is the current number of iterations, and Tmax is the maximum number of iterations;
S4:更新追踪者位置为:S4: Update tracker location to:
Figure FDA0002789862070000015
Figure FDA0002789862070000015
式中,
Figure FDA0002789862070000016
为第i只追踪者樽海鞘个体在第j维空间的位置,
Figure FDA0002789862070000017
为共生量;B为受益参数,随机选取1或2;r为[0,1]区间的随机数;
In the formula,
Figure FDA0002789862070000016
is the position of the i-th tracker salp individual in the j-th dimension space,
Figure FDA0002789862070000017
is the symbiotic amount; B is the benefit parameter, 1 or 2 is randomly selected; r is a random number in the [0,1] interval;
S5:后援者位置更新为:S5: Supporter location updated to:
Figure FDA0002789862070000018
Figure FDA0002789862070000018
Figure FDA0002789862070000019
Figure FDA0002789862070000019
式中,
Figure FDA00027898620700000110
为第i只后援者樽海鞘个体在第j维空间的位置,t为当前迭代次数;Tmax为最大迭代次数;
In the formula,
Figure FDA00027898620700000110
is the position of the i-th supporter salps in the j-th dimension space, t is the current number of iterations; T max is the maximum number of iterations;
S6:重复步骤S2~S5,直至迭代次数达到Tmax或食物适应度值达到设定终止阈值。S6: Repeat steps S2 to S5 until the number of iterations reaches T max or the food fitness value reaches the set termination threshold.
3.根据权利要求2所述的基于樽海鞘群算法的电极阵列优化方法,其特征在于:更新第i只追踪者樽海鞘个体在第j维空间的位置时,存储更新前的第i只追踪者樽海鞘个体在第j维空间的位置为
Figure FDA0002789862070000021
并判定
Figure FDA0002789862070000022
Figure FDA0002789862070000023
的适应度优劣,当
Figure FDA0002789862070000024
优于
Figure FDA0002789862070000025
时,再次更新第i只追踪者樽海鞘个体在第j维空间的位置为
Figure FDA0002789862070000026
3. The electrode array optimization method based on the salp group algorithm according to claim 2, characterized in that: when updating the position of the i-th tracker salp individual in the j-th dimension space, the i-th tracker before the update is stored. The position of the individual salps in the jth dimension space is
Figure FDA0002789862070000021
and judge
Figure FDA0002789862070000022
and
Figure FDA0002789862070000023
The fitness is good or bad, when
Figure FDA0002789862070000024
better than
Figure FDA0002789862070000025
When , the position of the i-th tracker salps in the j-th dimension space is updated again as
Figure FDA0002789862070000026
4.根据权利要求2所述的基于樽海鞘群算法的电极阵列优化方法,其特征在于:目标电极阵列排布空间为同心圆电极阵列,设定目标为最小化峰值旁瓣电平,各所述樽海鞘个体的适应度如下运算得到:4. The electrode array optimization method based on the salps group algorithm according to claim 2, characterized in that: the target electrode array arrangement space is a concentric electrode array, and the set target is to minimize the peak side lobe level. The fitness of Ascidian individuals is calculated as follows:
Figure FDA0002789862070000027
Figure FDA0002789862070000027
Figure FDA0002789862070000028
Figure FDA0002789862070000028
其中,θ为偏移法线的角度,
Figure FDA0002789862070000029
为电极子的相移量,j是虚数符号,Fmax为同心圆电极阵列的主瓣峰值,Ny为电极排布的圆环数量,Nx为环上电极子数量,rm指第m个圆环上第n个电极子(n,m)的半径,
Figure FDA00027898620700000210
馈电幅度是I(n,m)各电极子同激励,Ψ(n,m)是电极子与x轴的夹角。
where θ is the angle offset from the normal,
Figure FDA0002789862070000029
is the phase shift of the electrode, j is the sign of the imaginary number, F max is the peak value of the main lobe of the concentric electrode array, N y is the number of electrode rings, N x is the number of electrodes on the ring, r m refers to the mth The radius of the nth electrode element (n,m) on the ring,
Figure FDA00027898620700000210
The feeding amplitude is I(n,m) the co-excitation of each electrode, and Ψ(n,m) is the angle between the electrode and the x-axis.
5.根据权利要求2所述的基于樽海鞘群算法的电极阵列优化方法,其特征在于:目标电极阵列排布空间为矩形电极阵列,设定目标为最小化峰值旁瓣电平,各所述樽海鞘个体的适应度如下运算得到:5. The electrode array optimization method based on the salps group algorithm according to claim 2, wherein the target electrode array arrangement space is a rectangular electrode array, and the set target is to minimize the peak side lobe level, each of the described The fitness of individual salps is calculated as follows:
Figure FDA00027898620700000211
Figure FDA00027898620700000211
Figure FDA00027898620700000212
Figure FDA00027898620700000212
其中,θ为偏移法线的角度,
Figure FDA00027898620700000213
为电极子的相移量,j是虚数符号,Fmax为同心圆电极阵列的主瓣峰值,Ny为电极排布的圆环数量,Nx为环上电极子数量,rm指第m个圆环上第n个电极子(n,m)的半径,
Figure FDA00027898620700000214
pmx为电极在x方向上的位置,pmy是电极子在y方向上的位置,馈电幅度是I(n,m)各电极子同激励。
where θ is the angle offset from the normal,
Figure FDA00027898620700000213
is the phase shift of the electrodes, j is the sign of the imaginary number, F max is the peak value of the main lobe of the concentric electrode array, N y is the number of electrode rings, N x is the number of electrodes on the ring, r m refers to the mth The radius of the nth electrode element (n,m) on the ring,
Figure FDA00027898620700000214
p mx is the position of the electrode in the x direction, p my is the position of the electrode element in the y direction, and the feeding amplitude is I(n, m) where each electrode element is co-excited.
CN202011311093.1A 2020-11-20 2020-11-20 Electrode array optimization method based on goblet sea squirt group algorithm Active CN112287564B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011311093.1A CN112287564B (en) 2020-11-20 2020-11-20 Electrode array optimization method based on goblet sea squirt group algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011311093.1A CN112287564B (en) 2020-11-20 2020-11-20 Electrode array optimization method based on goblet sea squirt group algorithm

Publications (2)

Publication Number Publication Date
CN112287564A true CN112287564A (en) 2021-01-29
CN112287564B CN112287564B (en) 2023-04-07

Family

ID=74399322

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011311093.1A Active CN112287564B (en) 2020-11-20 2020-11-20 Electrode array optimization method based on goblet sea squirt group algorithm

Country Status (1)

Country Link
CN (1) CN112287564B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114584921A (en) * 2022-03-18 2022-06-03 中山大学 A wireless sensor network node location method based on the ring salps algorithm

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003105080A2 (en) * 2002-06-10 2003-12-18 Universite Libre De Bruxelles Metapopulation genetic algorithm for combinatorial optimisation problems
KR20180065472A (en) * 2016-12-08 2018-06-18 한국항공우주연구원 Apparatus and method for pattern synthesis of antenna array
CN109873810A (en) * 2019-01-14 2019-06-11 湖北工业大学 A phishing detection method based on salps swarm algorithm support vector machine
CN111027663A (en) * 2019-11-12 2020-04-17 天津大学 A method to improve the salp swarm algorithm

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003105080A2 (en) * 2002-06-10 2003-12-18 Universite Libre De Bruxelles Metapopulation genetic algorithm for combinatorial optimisation problems
KR20180065472A (en) * 2016-12-08 2018-06-18 한국항공우주연구원 Apparatus and method for pattern synthesis of antenna array
CN109873810A (en) * 2019-01-14 2019-06-11 湖北工业大学 A phishing detection method based on salps swarm algorithm support vector machine
CN111027663A (en) * 2019-11-12 2020-04-17 天津大学 A method to improve the salp swarm algorithm

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
A. DURMUŞ: "Circular Antenna Array Synthesis Using Salp Swarm Optimization", 《BALKAN JOURNAL OF ELECTRICAL & COMPUTER ENGINEERING》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114584921A (en) * 2022-03-18 2022-06-03 中山大学 A wireless sensor network node location method based on the ring salps algorithm

Also Published As

Publication number Publication date
CN112287564B (en) 2023-04-07

Similar Documents

Publication Publication Date Title
CN107657098B (en) Circular antenna array sparse method based on quantum chicken swarm evolution mechanism
CN103675799B (en) A kind of phased array sonar system transducer sparse battle array optimization method
CN106342323B (en) The submatrix weighted value of phased-array radar difference beam Sidelobe Suppression is determined method
CN105842702A (en) Multi-beam imaging sonar sidelobe suppression method and array sparsity method
CN110940970B (en) MIMO radar target detection method and system for floating oil sea surface
CN107657070A (en) Decoupling method and antenna based on the conformal array beamses of genetic algorithm frequency diversity
CN104408295B (en) A kind of Loads of Long-span Bridges substructure stormy waves coupling load method for numerical simulation
CN107390194B (en) A kind of radar target detection method based on the graceful divergence of full Donald Bragg
CN106569171A (en) Dual-layer-hybrid-array-based estimation method for direction angle of arrival
CN113791379B (en) Orthogonal matching pursuit DOA estimation method in nested array non-Gaussian environment
CN105954731A (en) Direction finding method of bistatic MIMO radar based on immune cuckold search under impulsive noise
CN112287564B (en) Electrode array optimization method based on goblet sea squirt group algorithm
CN109212466A (en) A kind of broadband direction-finding method based on quantum dragonfly mechanism of Evolution
CN109343006B (en) NFLM signal optimization method and device based on augmented Lagrange genetic algorithm
Jagannath et al. Grid matching for sparse signal recovery in compressive sensing
Liang et al. Identification of ship and corner reflector in sea clutter environment
CN109635486A (en) A kind of high resolution three-dimensional imaging sonar transducer array sparse optimization method
CN109212465A (en) A kind of particular array dynamic direction-finding method based on cultural ant lion mechanism
CN103323810A (en) L-array azimuthal angle and pitch angle paired signal processing method
CN113420912B (en) A method for identifying users with abnormal low voltage in distribution network
CN110113815A (en) A kind of improved wireless sensor network locating method based on IWO
CN114706084A (en) Underwater sound passive tracking multi-target method based on YOLO algorithm
CN104050686A (en) Novel intensive space target tracking method
CN118509804A (en) DV-Hop positioning method based on improved gold jackal optimization
CN104166795B (en) A kind of multiple sine wave frequency estimating methods based on many observation vector rarefaction representations

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