CN114065896A - Multi-objective Decomposition Evolutionary Algorithm Based on Neighbor Adjustment and Angle Selection Strategy - Google Patents

Multi-objective Decomposition Evolutionary Algorithm Based on Neighbor Adjustment and Angle Selection Strategy Download PDF

Info

Publication number
CN114065896A
CN114065896A CN202111239833.XA CN202111239833A CN114065896A CN 114065896 A CN114065896 A CN 114065896A CN 202111239833 A CN202111239833 A CN 202111239833A CN 114065896 A CN114065896 A CN 114065896A
Authority
CN
China
Prior art keywords
individual
neighborhood
parent
offspring
population
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.)
Pending
Application number
CN202111239833.XA
Other languages
Chinese (zh)
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.)
Xian University of Technology
Original Assignee
Xian University of Technology
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 Xian University of Technology filed Critical Xian University of Technology
Priority to CN202111239833.XA priority Critical patent/CN114065896A/en
Publication of CN114065896A publication Critical patent/CN114065896A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • 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]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/12Computing arrangements based on biological models using genetic models
    • G06N3/126Evolutionary algorithms, e.g. genetic algorithms or genetic programming

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Linguistics (AREA)
  • General Health & Medical Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Software Systems (AREA)
  • Artificial Intelligence (AREA)
  • Biomedical Technology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Evolutionary Biology (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Genetics & Genomics (AREA)
  • Physiology (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供了一种基于邻域调整和角度选择策略的多目标分解进化算法,包括:初始化获得初始种群P(t)、与初始种群P(t)各父代个体关联的权重向量、历史记忆邻域存档器的所有值ST初始化为定值;基于自适应邻域调整策略,分配初始种群P(t)中每个父代个体关联的权重向量邻域大小Ti;繁殖产生子代个体形成新种群Q(t);基于角度选择机制,遍历新种群Q(t)中所有子代个体选择优秀子代个体替换其父代个体;更新ST值为ST,pos,t+1并判断种群进化是否完成。本发明的多目标分解进化算法比现有的MOEA/D及其改进算法更能够有效地解决连续的复杂多目标问题,具有更好的收敛性和准确性。

Figure 202111239833

The invention provides a multi-objective decomposition evolution algorithm based on neighborhood adjustment and angle selection strategy, including: initializing to obtain an initial population P(t), a weight vector associated with each parental individual of the initial population P(t), a historical memory All values S T of the neighborhood archiver are initialized to fixed values; based on the adaptive neighborhood adjustment strategy, the weight vector neighborhood size T i associated with each parent individual in the initial population P(t) is assigned; reproduction produces offspring individuals Form a new population Q(t); based on the angle selection mechanism, traverse all the offspring individuals in the new population Q(t) to select excellent offspring individuals to replace their parent individuals; update the value of S T to S T, pos, t+1 and Determine whether the population evolution is complete. Compared with the existing MOEA/D and its improved algorithm, the multi-objective decomposition evolutionary algorithm of the invention can effectively solve the continuous complex multi-objective problem, and has better convergence and accuracy.

Figure 202111239833

Description

Multi-target decomposition evolution algorithm based on neighborhood adjustment and angle selection strategy
Technical Field
The invention belongs to the technical field of multi-objective optimization research in evolutionary computation, and particularly relates to a multi-objective decomposition evolutionary algorithm based on a neighborhood adjustment and angle selection strategy.
Background
In scientific research and engineering practice, it is often involved in simultaneously solving multiple conflicting objectives, which are collectively referred to in the field of optimization as Multi-Objective optimization Problems (Multi-Objective issues). Solving the multi-target problem has very important practical significance and scientific research value. The multi-objective optimization theory is urgently needed to guide the production and the life of people.
In general, all objective functions in a multi-objective optimization are conflicting, so there is no point in the decision space or objective space where all objective function values are optimal at the same time. Due to this special property, unlike the single objective optimization problem, the optimal solution of the multi-objective optimization problem is not a single solution but a set of compromise solutions. In practice, the trade-off between goals can be defined in pareto optima by giving the following definitions: 1. let u, n ∈ RmIf, if
Figure BDA0003319021050000011
Satisfy ui≤υiAnd is
Figure BDA0003319021050000012
uj<vjIf u dominates v, it is noted that u < v. 2. If there is no solution x ∈ Ω, so that F (x) < F (x)*) Then the variable x is decided*Referred to as Pareto optimal solution. Pareto optimal solution set (PS) is a set of Pareto optimal solutions, namely:
Figure BDA0003319021050000013
Figure BDA0003319021050000014
pareto Frontier (PF) is the set of target vectors corresponding to the Pareto optimal solution, i.e.: PF ═ f (x) x ∈ PS }. Currently, commonly used multi-objective evolutionary algorithms can be divided into three categories: 1. guiding the population to search and evolve by utilizing the pareto advantage relation among population solutions based on the pareto advantage algorithm; 2. an index-based algorithm by buildingA performance index to guide the search and evolution of the population; 3. a multi-objective optimization problem is decomposed into a plurality of single-objective subproblems by using weight vectors based on a decomposition algorithm.
The multi-objective evolutionary algorithm solves each sub-problem to finally obtain the PF of the whole multi-objective optimization problem, for example: firstly, converting a multi-objective optimization problem into a series of single-objective optimization subproblems by utilizing a weight vector preset in a target space based on a decomposed multi-objective evolutionary algorithm-MOEA/D; then, a group of subproblems is collaboratively evolved through a neighborhood concept defined by the relation between the subproblems, and a solution which is wanted by us is finally found, a more complex calculation index is not introduced, only less calculation resources are needed in the iteration process of the population, and the time complexity is lower; the existence of the weight vector enables the whole algorithm to have a hidden diversity protection mechanism, which is beneficial to the diversity distribution of the population and can generate a uniformly distributed solution by using a smaller population; the neighborhood in the MOEA/D realizes information sharing among individuals, and an aggregation function is introduced into the algorithm to optimize corresponding targets without considering conflicts among the targets, so that the method has the advantages of being very beneficial to adaptability distribution, rapid convergence of populations and the like. However, the decomposition-based multi-objective evolutionary algorithm-MOEA/D also has the following disadvantages:
the neighborhood used in MOEA/D is completely fixed, which not only can not meet the demand of the species group on prospectiveness and development, but also can cause certain resource waste.
2. The weight vector in the decomposition-based algorithm implies a population diversity protection mechanism, but MOEA/D adopting a TCH aggregation method does not fully utilize the characteristic.
Therefore, there is a need for improvement of the existing multi-objective evolutionary algorithm, MOEA/D.
Disclosure of Invention
In order to solve the defects of the multi-objective evolutionary algorithm MOEA/D, the invention discloses a multi-objective decomposition evolutionary algorithm based on neighborhood adjustment and angle selection strategies, and the adaptive neighborhood adjustment strategy is adopted to distribute the neighborhood size of the weight vector associated with each parent individual, so that the weight vectors of the parent individuals/child individuals have neighborhoods with different sizes in different evolution stages, and the requirements of different conditions on exploration and development are met. Meanwhile, the multi-objective decomposition evolution algorithm of the invention also improves the environment selection process, and selects excellent offspring individuals close to the weight vector as much as possible according to an angle selection mechanism by fully utilizing the distribution information of the individuals in the population so as to increase the diversity of the population and make the finally obtained solution distribution more uniform.
The technical scheme for realizing the purpose of the invention is as follows: a multi-objective decomposition evolution algorithm based on neighborhood adjustment and angle selection strategies comprises the following steps:
s1, acquiring a population in the problem search space, and initializing to obtain an initial population P (t) containing N parent individuals; carrying out grid design on a target space in a search space, and obtaining weight vectors associated with each grid and each parent individual of an initial population P (t) through WS-conversion; establishing a history memory neighborhood archiver with H entries, and storing all values S of the history memory neighborhood archiverTInitializing to a fixed value, and initializing an ideal point according to a population objective function value;
s2, distributing neighborhood size T to weight vector associated with each parent individual in initial population P (T) based on adaptive neighborhood adjustment strategyi
S3, reproducing each parent individual in the initial population P (t) to generate an offspring individual, and forming a new population Q (t);
s4, traversing all the descendant individuals in the new population Q (t) based on an angle selection mechanism, and selecting excellent descendant individuals to replace parent individuals;
s5, S for updating history memory neighborhood archiverTA value of ST,pos,t+1
And S6, judging whether the population evolution is finished or not, outputting the population after the evolution if the population evolution is finished, and repeating the steps S2 to S5 until the evolution is finished if the population evolution is not finished.
Further, in step S2, a neighborhood size T is assigned to the weight vector associated with each parent individual in the initial population p (T)iThe method comprises:
S201, randomly selecting an index number k in a history memory neighborhood archive device, wherein the k belongs to [1, H ];
s202, obtaining a numerical value S corresponding to the index number k according to the index number kT,kAnd according to the formula Ti=randci(ST,k0.1) calculating the neighborhood size TiWherein randc (μ, σ) represents a Cauchy distributed random value with μ as a mean and σ as a variance.
Furthermore, in step S202, if T is detectediIf the set threshold range is exceeded, the T is again determinedi=randci(ST,k0.1) calculating until T satisfying the threshold range is generatediThe value is obtained.
Further, in the step S4, based on the angle selection mechanism, the method of traversing all the offspring individuals in the new population q (t) and selecting the superior offspring individual to replace the parent individual comprises:
s401, calculating an aggregation function value of each offspring individual in the new population Q (t), and calculating an angle theta between each weight vector in a neighborhood where the offspring individual is located and the associated parent individualiCalculating the included angle theta between the child individual and all the weight vectors in the neighborhoodi
S402, updating the offspring individuals in the neighborhood, and selecting excellent offspring individuals to replace parent individuals;
s4021, when t is less than or equal to 0.8 × maxgen, maxgen is the maximum evolution algebra, and for a certain weight vector K in the neighborhood of a descendant individual, K belongs to [1, i ];
if the aggregation function value of the offspring is superior to that of the parent
Figure BDA0003319021050000041
The child individual is an excellent child individual, and the parent individual associated with the weight vector K is replaced;
if the aggregation function value of the offspring individual is superior to that of the parent individual, but
Figure BDA0003319021050000042
Then by the formula
Figure BDA0003319021050000043
Determining the probability that the descendant replaces the parent associated with the weight vector K,
Figure BDA0003319021050000044
an improved value of an aggregation function for the offspring individual relative to the parent individual;
s4022, when 0.8 × maxgen is less than or equal to t, only comparing the aggregation function value of the filial generation individual with that of the parent generation individual, and if the filial generation individual has the aggregation function value
Figure BDA0003319021050000051
If the number of the child individuals is smaller than that of the parent individuals, the child individuals replace the parent individuals associated with the weight vectors K for excellent child individuals;
s403, adding the excellent offspring individuals replacing the parent individuals into the GoodT matrix, and recording the excellent offspring individuals
Figure BDA0003319021050000052
Into dif _ val.
Further, in the above step S5,
Figure BDA0003319021050000053
wherein G represents an algebraic generation of evolution; k is an element of [1, H ]]Representing the index value in the archive, and determining the updated position in the history memory neighborhood archiver; mean is a measure ofWL(GoodT) A weighted lehmer average representing the size of the neighborhood.
Further, the mean mentioned aboveWL(GoodT) The calculation formula of (2) is as follows:
Figure BDA0003319021050000054
wherein
Figure BDA0003319021050000055
Is the improved value of the aggregation function for the offspring individual relative to the parent individual.
Compared with the prior art, the invention has the beneficial effects that:
1. after the neighborhood self-adaptive strategy based on historical performance is adopted, the weight vectors associated with the parent individuals/the child individuals can have neighborhoods with different sizes in different evolution stages, on one hand, the waste of computing resources can be reduced, on the other hand, the balance between population development and exploration is favorably maintained, the global search capability and the local search capability of the algorithm are balanced, and the population is prevented from falling into local optimum.
2. The new selection mechanism (angle selection mechanism) can effectively select the child individuals which are close to the weight vector and have better distribution, improves the diversity of the population, enables the obtained solution to be more uniformly distributed, and optimizes the performance of the algorithm.
3. Compared with the existing MOEA/D and partial improved algorithm thereof, the multi-target decomposition evolution algorithm (MOEA/D-ADN) provided by the invention can effectively solve the continuous complex multi-target problem.
Drawings
In order to more clearly illustrate the technical solution of the embodiment of the present invention, the drawings used in the description of the embodiment will be briefly introduced below. It should be apparent that the drawings in the following description are only for illustrating the embodiments of the present invention or technical solutions in the prior art more clearly, and that other drawings can be obtained by those skilled in the art without any inventive work.
FIG. 1 is a block diagram of a flow of a multi-objective decomposition-evolution algorithm based on neighborhood adjustment and angle selection strategies according to the present invention;
FIG. 2 is a flowchart of an embodiment of an algorithm of a multi-objective decomposition evolution algorithm based on a neighborhood adjustment and angle selection strategy;
FIG. 3 is a solution set obtained by running the multi-objective decomposition evolution algorithm based on the neighborhood adjustment and angle selection strategy on the test function WFG8 according to the present invention;
FIG. 4 is a solution set obtained by running the prior art algorithm MOEA/D on the test function WFG 8;
FIG. 5 is a solution set resulting from the prior art algorithm MOEA/D-CMA running on test function WFG 8.
FIG. 6 is a solution set resulting from the prior art algorithm MOEA/D-DE running on test function WFG 8.
FIG. 7 is a solution set resulting from running the prior art algorithm MOEA/D-M2M on test function WFG 8.
FIG. 8 is a solution set resulting from the prior art algorithm MOEA/D-Pas running on test function WFG 8.
Fig. 9 is a standard Pareto front of test function WFG 8.
Detailed Description
The invention will be further described with reference to specific embodiments, and the advantages and features of the invention will become apparent as the description proceeds. These examples are illustrative only and do not limit the scope of the present invention in any way. It will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention, and that such changes and modifications may be made without departing from the spirit and scope of the invention.
The specific embodiment provides a multi-objective decomposition evolution algorithm based on neighborhood adjustment and angle selection strategies, as shown in fig. 1 and 2, comprising the following steps:
step 1: the initialization process of the population, the weight vector and the history memory neighborhood archiver comprises the steps of S1:
s1, acquiring a population in the problem search space, and initializing to obtain an initial population P (t) containing N parent individuals; carrying out grid design on a target space in a search space, and obtaining weight vectors associated with each grid and each parent individual of an initial population P (t) through WS-conversion; establishing a history memory neighborhood archiver with H entries, and storing all values S of the history memory neighborhood archiverTInitializing to a fixed value, and initializing an ideal point according to a population objective function value.
Specifically, a simplex mesh design method may be selected in the target space to divide the target space into N meshes, and WS-transform initialization is performed on the N meshes to obtain N weight vectors, where a more specific WS-transform method is as follows:
Figure BDA0003319021050000071
specifically, the size H of the neighborhood archive in the history memory neighborhood archiver is set to be a fixed value, and the value S in the neighborhood archive is initializedT,i(i is more than or equal to 1 and less than or equal to H) is a constant value, and the filing update position pos is initialized to 1;
specifically, the population is initialized to an initial population p (t) containing N parent individuals, where t is an evolution algebra and an initial value thereof is 0, and the population is initialized in the following manner: first, a continuous, uniformly distributed decision variable is generated in the decision space, i.e.
Figure BDA0003319021050000081
i is the dimension of the decision variable and satisfies 1 ≦ i ≦ N, k is the population size and satisfies 1 ≦ k ≦ N, xi_up、xi_lowRepresenting the upper and lower bounds of a certain dimension of the decision variables; secondly, calculating and generating an objective function value of the parent individual; thirdly, initializing the ideal point according to the target function value of the individual in P (t). It should be noted that the objective function value and the objective function value initialization ideal point are calculated by using the existing method, and are not described again.
Step 2: the evolutionary propagation process, namely selecting the parents to perform operations such as mutation, intersection and the like to generate the next generation until the population evolution is completed, comprises the following steps S2 to S3.
S2, distributing neighborhood size T to weight vector associated with each parent individual in initial population P (T) based on adaptive neighborhood adjustment strategyi
S201, randomly selecting an index number k in a history memory neighborhood archive device, wherein the k belongs to [1, H ];
s202, obtaining a numerical value S corresponding to the index number k according to the index number kT,kAnd according to the formula Ti=randci(ST,k0.1) calculating the neighborhood size TiWherein randc (μ, σ) represents a Cauchy distributed random value with μ as a mean and σ as a variance. Preferably, if TiIf the set threshold range is exceeded, the T is again determinedi=randci(ST,k,01) calculating until T satisfying a threshold range is generatediThe value is obtained.
S3, breeding each parent in the initial population p (t) to generate an offspring individual, and forming a new population q (t), specifically, breeding each parent in the initial population p (t) generates an offspring individual by a variant crossing manner.
Step 3: the environment selecting process, which selects the excellent descendant to enter the next generation by traversing all descendant individuals in q (t), includes the steps S4:
s4, traversing all the descendant individuals in the new population Q (t) based on an angle selection mechanism, and selecting excellent descendant individuals to replace parent individuals;
s401, calculating an aggregation function value of each offspring individual in the new population Q (t), and calculating an angle theta between each weight vector in a neighborhood where the offspring individual is located and the associated parent individualiCalculating the included angle theta between the child individual and all the weight vectors in the neighborhoodi
S402, updating the offspring individuals in the neighborhood, and selecting excellent offspring individuals to replace parent individuals;
s4021, when t is less than or equal to 0.8 × maxgen, maxgen is the maximum evolution algebra, and for a certain weight vector K in the neighborhood of a descendant individual, K belongs to [1, i ];
if the aggregation function value of the offspring is superior to that of the parent
Figure BDA0003319021050000091
The child individual is an excellent child individual, and the parent individual associated with the weight vector K is replaced;
if the aggregation function value of the offspring individual is superior to that of the parent individual, but
Figure BDA0003319021050000092
Then by the formula
Figure BDA0003319021050000093
Determining the probability that the descendant replaces the parent associated with the weight vector K,
Figure BDA0003319021050000094
an improved value of an aggregation function for the offspring individual relative to the parent individual;
s4022, when 0.8 × maxgen is less than or equal to t, only comparing the aggregation function value of the filial generation individual with that of the parent generation individual, and if the filial generation individual has the aggregation function value
Figure BDA0003319021050000095
If the number of the child individuals is smaller than that of the parent individuals, the child individuals replace the parent individuals associated with the weight vectors K for excellent child individuals;
s403, adding the excellent offspring individuals replacing the parent individuals into the GoodT matrix, and recording the excellent offspring individuals
Figure BDA0003319021050000096
In dif _ val, in this step, the superior offspring of the parent is successfully replaced, and the neighborhood size T of the weight vector associated with the parent that produces the superior offspring is considered as the neighborhood size T of the weight vector associated with the parentiIs excellent.
Step 4: updating the file, i.e. updating the S of the archive archiver in the history memory neighborhood after a generation of individuals in the population has evolvedTValue including step S5, S for updating history memory neighborhood archiverTA value of ST,pos,t+1
In particular, the method comprises the following steps of,
Figure BDA0003319021050000101
wherein G represents an algebraic generation of evolution; k is an element of [1, H ]]Representing the index value in the archive, and determining the updated position in the history memory neighborhood archiver; mean is a measure ofWL(GoodT) A weighted lehmer average representing the size of the neighborhood. Mean aboveWL(GoodT) The calculation formula of (2) is as follows:
Figure BDA0003319021050000102
Figure BDA0003319021050000103
wherein
Figure BDA0003319021050000104
Is the improved value of the aggregation function for the offspring individual relative to the parent individual.
It should be noted that, at the beginning, the archive update position pos is initialized to 1, and S is updated every timeTThe value of pos is incremented by one, and if pos > H, the value of pos needs to be reset to 1.
Step 5: judging whether population evolution is finished or not and whether termination conditions are met or not, wherein the method comprises the following steps of S6:
and S6, judging whether the population evolution is finished or not, outputting the population after the evolution if the population evolution is finished, and repeating the steps S2 to S5 until the evolution is finished if the population evolution is not finished.
Specifically, if the function evaluation frequency FES is less than or equal to MAXFES, it indicates that the evolution is not completed, and the steps S2 to S5 need to be repeated before judgment; if FES > MAXFES, the population evolution is completed, and the current population P (t) is output, wherein MAXFES represents the maximum function evaluation times.
The effect of the multi-objective decomposition evolution algorithm based on the neighborhood adjustment and angle selection strategy is further illustrated by the following simulation experiment:
first, given test functions, the performance of the proposed algorithm was tested by selecting multi-mode multi-target test functions with various problem characteristics (e.g., spoofing, deflection, multi-mode) from the WFG test suite, WFG 1-9, the test set of which is specifically shown in Table 1 below.
Next, algorithm parameters are initialized, and in this experiment, the target number is set to m-3. The number of decision variables is set to n-2 (m-1) + 20. The population size N of the test function is 91, the maximum evaluation times MAXFES is 36400, and the test function is independently operated for 20 times. During initialization, H is set to N/10, ST,i(1. ltoreq. i.ltoreq.H) is set to N/10.
Table 1: WFG 1-9 test function set
Figure BDA0003319021050000111
Figure BDA0003319021050000121
In order to better explain the multi-objective decomposition evolution algorithm machine, the existing MOEA/D algorithm and the variant algorithms thereof MOEA/D-DE, MOEA/D-M2M, MOEA/D-CMA and MOEA/D-Pas are adopted to carry out experiments on the same test functions, and the performances of the multi-objective decomposition evolution algorithm machine are compared.
The selected test performance index is an Inverted Generated Distance (IGD), the IGD represents an average distance between a real Pareto front edge of the multi-target problem and an individual solution set obtained by the algorithm, the IGD index is adopted to evaluate the performance of the algorithm, and reliable information about diversity and convergence of the obtained solution can be provided.
Suppose P*For a set evenly distributed on the real Pareto frontier, P is the individual solution set obtained by the algorithm, IGD measures P*And P, the calculation formula is as follows:
Figure BDA0003319021050000131
Figure BDA0003319021050000132
wherein d (v, P) represents from P*The euclidean distance between the starting point v and the nearest member of P. As can be known from the definition of IGD, the IGD index can measure the diversity and convergence of the solution finally obtained by the algorithm at the same time, and the smaller the IGD value is, the more uniformly the obtained solution set is distributed in the target space, and the smaller the distance from the real PF is.
Comparing the operation result of the multi-objective decomposition evolution algorithm with IGD indexes obtained by the operation results of the algorithms such as MOEA/D, MOEA/D-DE, MOEA/D-M2M, MOEA/D-CMA, MOEA/D-Pas and the like: all tests were run 20 times independently and the mean and standard deviation of the metric values were recorded and the results were statistically analyzed using a Wilcoxon rank sum test with a significance level of 0.05, as shown in table 2 below, where the symbols "+", "-" and "═ respectively indicate that the current algorithm gave significantly better, significantly worse and statistically similar results to the MOEA/D.
TABLE 2 results of experiments with IGD index obtained by running five algorithms 20 times on WFG 1-WFG 9
Figure BDA0003319021050000133
Figure BDA0003319021050000141
The data in the table is divided into two rows, the first row being the mean of the indices found in 20 independent runs, the second row being the standard deviation, the bold data representing the best results for this example of the problem in this experiment. From the data in the table, one can derive: in 9 examples, the multi-objective decomposition evolution algorithm (MOEA/D-AND) of the invention is optimal on 6 example problems, AND has better convergence AND accuracy in 6 algorithms.
In order to better compare the solution sets obtained by the operation results of the multi-objective decomposition evolution algorithm (MOEA/D-ADN) of the invention with the solution sets obtained by the operation results of the algorithms such as MOEA/D, MOEA/D-DE, MOEA/D-CMA, MOEA/D-Pas and MOEA/D-M2M, the operation results of the algorithms on a three-objective test function WFG8 are compared, and the operation results are shown in FIGS. 3 to 9, wherein 9 represents a standard Pareto front edge of the three-objective test function WFG8, gray points in the graphs 3 to 9 represent the distribution of the solution sets finally obtained by the algorithms on a target space, and the distribution can be seen from the graphs: on the three-target WFG8 problem, the solution sets obtained by the two algorithms in FIG. 3 and FIG. 5 are distributed well, the coincidence rate with the standard Pareto front edge is the highest among several algorithms, and the rest algorithms are inferior to the two algorithms and have disordered solution distribution.
Finally, it should be noted that the above embodiments are only used for illustrating the technical solutions of the present invention, and are preferred embodiments of the present invention, not limiting the same; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art can modify the technical solutions described in the foregoing embodiments, or make various changes and modifications to the present invention without departing from the spirit and scope of the present invention, and therefore, the present invention should be covered by the protection scope of the present invention.

Claims (6)

1.一种基于邻域调整和角度选择策略的多目标分解进化算法,其特征在于,包括以下步骤:1. a multi-objective decomposition evolutionary algorithm based on neighborhood adjustment and angle selection strategy, is characterized in that, comprises the following steps: S1、获取问题搜索空间内的种群,初始化后得到含有N个父代个体的初始种群P(t);对搜索空间内目标空间进行网格设计,经WS-转换得到各网格与初始种群P(t)各个父代个体关联的权重向量;建立带有H个条目的历史记忆邻域存档器,将历史记忆邻域存档器的所有值ST初始化为定值,根据种群目标函数值初始化理想点;S1. Obtain the population in the problem search space, and obtain the initial population P(t) containing N parent individuals after initialization; design the grid for the target space in the search space, and obtain each grid and the initial population P through WS-transformation (t) The weight vector associated with each parent individual; establish a historical memory neighborhood archiver with H entries, initialize all the values S T of the historical memory neighborhood archiver to a fixed value, and initialize the ideal value according to the population objective function value point; S2、基于自适应邻域调整策略,为初始种群P(t)中每个父代个体关联的权重向量分配邻域大小TiS2. Based on the adaptive neighborhood adjustment strategy, assign a neighborhood size T i to the weight vector associated with each parent individual in the initial population P(t); S3、对初始种群P(t)中各个父代个体进行繁殖产生子代个体,形成新种群Q(t);S3. Propagating each parent individual in the initial population P(t) to generate offspring individuals to form a new population Q(t); S4、基于角度选择机制,遍历新种群Q(t)中所有子代个体,选择优秀子代个体替换其父代个体;S4. Based on the angle selection mechanism, traverse all the offspring individuals in the new population Q(t), and select excellent offspring individuals to replace their parent individuals; S5、更新历史记忆邻域存档器的ST值为ST,pos,t+1S5. The S T value of the updated historical memory neighborhood archiver is S T,pos,t+1 ; S6、判断种群进化是否完成,若完成则输出进化后的种群,若未完成则重复步骤S2至S5直至进化完成。S6. Determine whether the evolution of the population is completed. If completed, output the evolved population. If not, repeat steps S2 to S5 until the evolution is completed. 2.根据权利要求1所述的基于邻域调整和角度选择策略的多目标分解进化算法,其特征在于,步骤S2中,为初始种群P(t)中每个父代个体关联的权重向量分配邻域大小Ti的方法包括:2. The multi-objective decomposition evolutionary algorithm based on neighborhood adjustment and angle selection strategy according to claim 1, is characterized in that, in step S2, it is the weight vector associated with each parent individual in the initial population P(t) assigned Methods of neighborhood size T i include: S201、在历史记忆邻域存档器随机选取一个索引号k,k∈[1,H];S201, randomly select an index number k, k∈[1,H] in the historical memory neighborhood archiver; S202、依据索引号k获取与索引号k对应的数值ST,k,并依据公式Ti=randci(ST,k,0.1)计算邻域大小Ti的值,其中randc(μ,σ)代表的是以μ为均值,σ为方差的柯西分布随机值。S202, obtain the value S T, k corresponding to the index number k according to the index number k, and calculate the value of the neighborhood size T i according to the formula T i =randc i (S T,k ,0.1), where randc(μ,σ ) represents a random value from a Cauchy distribution with μ as the mean and σ as the variance. 3.根据权利要求2所述的基于邻域调整和角度选择策略的多目标分解进化算法,其特征在于:步骤S202中,若Ti超出设定阈值范围,则重新根据Ti=randci(ST,k,0.1)计算直至产生满足阈值范围的Ti值。3. the multi-objective decomposition evolutionary algorithm based on neighborhood adjustment and angle selection strategy according to claim 2, is characterized in that: in step S202, if T i exceeds the setting threshold range, then again according to T i =randc i ( S T,k ,0.1) is calculated until a T i value that satisfies the threshold range is generated. 4.根据权利要求1所述的基于邻域调整和角度选择策略的多目标分解进化算法,其特征在于,步骤S4中基于角度选择机制,遍历新种群Q(t)中所有子代个体,选择优秀子代个体替换其父代个体的方法为:4. The multi-objective decomposition evolutionary algorithm based on neighborhood adjustment and angle selection strategy according to claim 1, it is characterized in that, based on angle selection mechanism in step S4, traverse all the offspring individuals in the new population Q(t), select The method for the excellent offspring to replace its parent is: S401、计算新种群Q(t)中每个子代个体的聚合函数值,计算子代个体所在邻域中每个权重向量与其相关联的父代个体之间的角度θi,计算子代个体与其邻域中所有权重向量的夹角θiS401. Calculate the aggregation function value of each child individual in the new population Q(t), calculate the angle θ i between each weight vector and its associated parent individual in the neighborhood where the child individual is located, and calculate the child individual and its associated parent. the included angle θ i of all weight vectors in the neighborhood; S402、在邻域中进行子代个体更新,选择优秀子代个体替换父代个体;S402, update the offspring individuals in the neighborhood, and select excellent offspring individuals to replace the parent individuals; S4021、在t≤0.8*maxgen时,maxgen为最大进化代数,对于子代个体的邻域中某个权重向量K,K∈[1,i];S4021. When t≤0.8*maxgen, maxgen is the maximum evolutionary algebra, and for a certain weight vector K in the neighborhood of the offspring individual, K∈[1,i]; 若子代个体的聚合函数值比父代个体的聚合函数值优越,且θk≤θk,则该子代个体为优秀子代个体,替换权重向量K所关联的父代个体;If the aggregate function value of the offspring individual is superior to that of the parent individual, and θ k ≤ θ k , the offspring individual is an excellent offspring individual, and the parent individual associated with the weight vector K is replaced; 若子代个体的聚合函数值比父代个体的聚合函数值优越,但θkk,则以公式
Figure FDA0003319021040000021
确定子代个体替换权重向量K所关联的父代个体的概率,
Figure FDA0003319021040000022
为子代个体相对于父代个体的聚合函数的改进值;
If the aggregate function value of the offspring individual is superior to that of the parent individual, but θ kk , then the formula
Figure FDA0003319021040000021
Determine the probability that the offspring individual replaces the parent individual associated with the weight vector K,
Figure FDA0003319021040000022
is the improvement value of the aggregate function of the offspring individual relative to the parent individual;
S4022、在0.8*maxgen≤t时,则只比较子代个体的聚合函数值与父代个体的聚合函数值,若子代个体的
Figure FDA0003319021040000023
小于父代个体,则该子代个体为优秀子代个体替换权重向量K所关联的父代个体;
S4022. When 0.8*maxgen≤t, only the aggregate function value of the offspring individual is compared with the aggregate function value of the parent individual. If the offspring individual’s aggregate function value is
Figure FDA0003319021040000023
is smaller than the parent individual, then the child individual is an excellent child individual to replace the parent individual associated with the weight vector K;
S403、将替换父代个体的优秀子代个体加入矩阵GoodT中,并记录该优秀子代个体的
Figure FDA0003319021040000031
到dif_val中。
S403. Add the outstanding offspring individual that replaces the parent individual into the matrix GoodT, and record the outstanding offspring individual
Figure FDA0003319021040000031
into dif_val.
5.根据权利要求1所述的基于邻域调整和角度选择策略的多目标分解进化算法,其特征在于:步骤S5中,
Figure FDA0003319021040000032
其中G代表的是进化的代数;k∈[1,H]代表的是存档中的索引值,决定了历史记忆邻域存档器中更新的位置;meanWL(GoodT)代表邻域大小的加权莱默平均值。
5. the multi-objective decomposition evolutionary algorithm based on neighborhood adjustment and angle selection strategy according to claim 1, is characterized in that: in step S5,
Figure FDA0003319021040000032
where G represents the algebra of evolution; k∈[1,H] represents the index value in the archive, which determines the updated position in the historical memory neighborhood archiver; mean WL (Good T ) represents the weighted neighborhood size Laimer average.
6.根据权利要求5所述的基于邻域调整和角度选择策略的多目标分解进化算法,其特征在于:meanWL(GoodT)的计算公式为:6. the multi-objective decomposition evolutionary algorithm based on neighborhood adjustment and angle selection strategy according to claim 5, is characterized in that: the calculation formula of mean WL (Good T ) is:
Figure FDA0003319021040000033
其中
Figure FDA0003319021040000034
为子代个体相对于父代个体的聚合函数的改进值。
Figure FDA0003319021040000033
in
Figure FDA0003319021040000034
is the improvement value of the aggregate function of the offspring individuals relative to the parent individuals.
CN202111239833.XA 2021-10-25 2021-10-25 Multi-objective Decomposition Evolutionary Algorithm Based on Neighbor Adjustment and Angle Selection Strategy Pending CN114065896A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111239833.XA CN114065896A (en) 2021-10-25 2021-10-25 Multi-objective Decomposition Evolutionary Algorithm Based on Neighbor Adjustment and Angle Selection Strategy

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111239833.XA CN114065896A (en) 2021-10-25 2021-10-25 Multi-objective Decomposition Evolutionary Algorithm Based on Neighbor Adjustment and Angle Selection Strategy

Publications (1)

Publication Number Publication Date
CN114065896A true CN114065896A (en) 2022-02-18

Family

ID=80235393

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111239833.XA Pending CN114065896A (en) 2021-10-25 2021-10-25 Multi-objective Decomposition Evolutionary Algorithm Based on Neighbor Adjustment and Angle Selection Strategy

Country Status (1)

Country Link
CN (1) CN114065896A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117035373A (en) * 2023-10-09 2023-11-10 中国电建集团山东电力管道工程有限公司 Intelligent management method and system for pipeline prefabrication production line
CN117077981A (en) * 2023-10-16 2023-11-17 四川大学 Method and device for distributing stand by fusing neighborhood search variation and differential evolution
CN118430835A (en) * 2024-07-02 2024-08-02 中国人民解放军总医院 A multi-objective clinical decision-making method and system based on collaboration of large and small models

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117035373A (en) * 2023-10-09 2023-11-10 中国电建集团山东电力管道工程有限公司 Intelligent management method and system for pipeline prefabrication production line
CN117035373B (en) * 2023-10-09 2024-01-23 中国电建集团山东电力管道工程有限公司 Intelligent management method and system for pipeline prefabrication production line
CN117077981A (en) * 2023-10-16 2023-11-17 四川大学 Method and device for distributing stand by fusing neighborhood search variation and differential evolution
CN117077981B (en) * 2023-10-16 2024-02-02 四川大学 Method and device for distributing stand by fusing neighborhood search variation and differential evolution
CN118430835A (en) * 2024-07-02 2024-08-02 中国人民解放军总医院 A multi-objective clinical decision-making method and system based on collaboration of large and small models
CN118430835B (en) * 2024-07-02 2024-09-27 中国人民解放军总医院 A multi-objective clinical decision-making method and system based on collaboration of large and small models

Similar Documents

Publication Publication Date Title
CN114065896A (en) Multi-objective Decomposition Evolutionary Algorithm Based on Neighbor Adjustment and Angle Selection Strategy
Liu et al. Many-objective evolutionary optimization based on reference points
Drugan et al. Designing multi-objective multi-armed bandits algorithms: A study
US8825573B2 (en) Controlling quarantining and biasing in cataclysms for optimization simulations
CN109271320B (en) A method for prioritization of upper-level multi-objective test cases
CN108153918B (en) A Personalized Movie Recommendation Method Based on Adaptive Orthogonal Crossover Multi-objective Optimization Algorithm
CN114862090A (en) Workshop scheduling method and system based on improved multi-target genetic algorithm
CN115034557A (en) Agile satellite emergency task planning method
JP5019744B2 (en) Multi-objective optimization apparatus, multi-objective optimization method, and multi-objective optimization program
CN118095070A (en) Multi-scene multi-objective problem optimization method and device, electronic equipment and storage medium
Song et al. A hybrid differential evolution for multi-objective optimisation problems
Kuk et al. A new adaptive operator selection for NSGA-III applied to CEC 2018 many-objective benchmark
JP7137074B2 (en) Optimization calculation method, optimization calculation device, and optimization calculation program
CN118297353B (en) Industrial production process multi-objective optimization method based on branch non-dominant sorting algorithm
Saha et al. Exploiting linear interpolation of variational autoencoders for satisfying preferences in evolutionary design optimization
CN118708323A (en) A workflow scheduling method and device based on multi-objectives
CN109934344B (en) An Improved Rule-Based Model-Based Multi-objective Distribution Estimation Method
Jameel et al. A Reference Point‐Based Evolutionary Algorithm Solves Multi and Many‐Objective Optimization Problems: Method and Validation
CN113094973B (en) Civil aircraft demand optimization method based on multi-objective optimization algorithm
Starkey et al. Fuzzy dominance rules for real-world many objective optimization
Alkafaween et al. CSGA: A DUAL POPULATION GENETIC ALGORITHM BASED ON MEXICAN CAVEFISH GENETIC DIVERSITY.
Ruane et al. Simulation and Genetic Algorithms to Improve the Performance of an Automated Manufacturing Line
Kong et al. Large-dimensional multi-objective evolutionary algorithms based on improved average ranking
CN119989540B (en) Target design method for lightweight and modal optimization of rear subframe structure based on ML
Yang et al. An Many-Objective Optimization Algorithm with Reference Point and Angle based on Non-Dominated Sorting Approach

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20220218

RJ01 Rejection of invention patent application after publication