2.5d integrated circuit test path optimization method based on decomposition multi-objective evolutionary algorithm

By optimizing the test path of 2.5D integrated circuits based on the decomposition multi-objective evolutionary algorithm, the contradiction between test time and hardware overhead is resolved, and the cost is effectively reduced.

CN116843024BActive Publication Date: 2026-01-09HARBIN INST OF TECH AT WEIHAI
View PDF 0 Cites 0 Cited by

Patent Information

Application Number
CN202310808109.7
Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
Filing Date
2023-07-03
Publication Date
2026-01-09
Estimated Expiration
2043-07-03

AI Technical Summary

Technical Problem

Existing technologies struggle to efficiently balance testing time and hardware overhead in 2.5D integrated circuit testing, resulting in high testing costs and unsatisfactory results from traditional optimization methods.

Method used

A test path optimization method based on decomposition multi-objective evolution algorithm is adopted. Through steps such as initialization, crossover and mutation operators, and environment selection, the test path configuration is optimized to balance test time and hardware overhead.

Benefits of technology

It effectively reduces the testing cost of 2.5D integrated circuits, achieves an efficient balance between testing time and hardware overhead, and reduces the overall testing cost.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN116843024B_ABST
    Figure CN116843024B_ABST
Patent Text Reader

Abstract

The 2.5D integrated circuit test path optimization method based on decomposition multi-objective evolutionary algorithm belongs to the field of integrated circuit test optimization, and aims to solve the problems that the traditional optimization method cannot efficiently solve the multi-objective optimization problem and the effect of balancing test time and test hardware overhead is not ideal because the gradient information of the problem is usually required. The method comprises the following steps: step one, initialization, generating an initial population for the test path configuration scheme of the 2.5D integrated circuit, and evaluating the initial population index according to the multi-objective evolutionary algorithm, and selecting the optimal fitness index as the ideal point of the initial population; step two, generating a new solution by using a crossover and mutation operator, and evaluating the new solution according to the multi-objective evolutionary algorithm and updating the ideal point; step three, environment selection, including updating the neighborhood solution and the elite solution set; step four, executing the iteration process of steps two and three for multiple times until the iteration is completed, and outputting the elite solution set.
Need to check novelty before this filing date? Find Prior Art

Description

TECHNICAL FIELD

[0001] The present application relates to a 2.5D integrated circuit test path optimization method, and belongs to the field of integrated circuit test optimization. BACKGROUND

[0002] At present, 2.5D integrated circuits, as the forerunner of 3D integrated circuits, have become one of the mainstream directions of the development of the integrated circuit industry. The core of the 2.5D integrated circuit is the passive silicon interposer (PSI) technology. Multiple wafers composed of intellectual property (IP) cores and related functional circuits are placed side by side above the PSI, and are interconnected between the wafers through horizontal redistribution layers and micro solder balls, and are connected to the outside through through-silicon-via (TSV) and controllable collapse chip connection solder balls.

[0003] The test problem of the 2.5D integrated circuit runs through the entire development process and is an important guarantee of product quality. The complex design structure, high integration degree and large number of IP cores make the test process of the 2.5D integrated circuit more complex and the test more difficult. The high-density internal interconnection and external pins result in a large increase in the number of test paths, and the test time and test cost required are also increased. In order to reduce the test time of the internal interconnection, the common practice at present is to divide a single long test path into multiple short test paths to improve the test efficiency. However, multiple test paths usually require additional hardware overhead such as through-silicon-via structures and micro solder balls, thereby increasing the test cost. It can be seen that the test time and the hardware overhead required for testing are two contradictory objectives, and therefore can be abstracted as a multi-objective optimization problem.

[0004] For any given 2.5D integrated circuit test architecture, there can be a large number of different test path configurations. For example, there are 13 test path configurations for 3 wafers on the PSI, as shown in Figure 1 For 4 wafers, there are 73 test path configurations, and for 12 wafers, there are 12470162233 test path configurations. It can be seen that the number of test paths increases exponentially with the number of wafers, and therefore the optimization problem is an NP-hard problem, and efficient optimization techniques are needed to find the optimal test path configuration.

[0005] In 2.5D integrated circuit testing, a single long test path (usually refers to a test scan chain connecting a pair of TDI interface and TDO interface, TDI represents test data input interface, and TDO represents test data output interface) consumes extremely long test time. In order to reduce time consumption, a long test path is usually divided into multiple shorter test paths for testing, but multiple test paths require additional TSV and micro-bumps and other hardware, which increases hardware overhead. Therefore, balancing test time and hardware overhead and reducing comprehensive test cost are problems that must be considered in integrated circuit test path design.

[0006] In the testing process, a test path containing one or more internal scan chains is usually referred to as an input test access mechanism (in-TAM), and a test path containing one or more external scan chains is referred to as an output test access mechanism (out-TAM). The connection between two wafers in the form of in-TAM and out-TAM is referred to as a basic test path connection, as shown in FIG. 1. In the basic test path connection, the internal scan chains of the two wafers are connected through test path 1 (P1) between the TDI and TDO interfaces, and the external scan chains are connected through test path 2 (P2) between the TDI_new and SO_new interfaces. When a long in-TAM or out-TAM is divided into multiple short in-TAMs or out-TAMs, each in-TAM requires a vertical path to change the test mode of the probe, and each out-TAM requires a vertical path to output the test response from the 2.5D integrated circuit, so TSV and micro-bumps are required to form the above-mentioned vertical test path. Figure 2

[0007] Using a single long test path alone results in long test time, and using multiple short test paths reduces test time but also increases hardware overhead, so optimization techniques are needed to balance test time and test hardware overhead and thus reduce comprehensive test cost. In solving multi-objective optimization problems, traditional optimization methods usually require gradient information of the problem and cannot efficiently solve the problem, and the effect of balancing test time and test hardware overhead is not ideal. SUMMARY

[0008] ​In view of the problems that the traditional optimization method cannot efficiently solve the multi-objective optimization problem, usually needs gradient information of the problem, and the effect of balancing test time and test hardware overhead is not ideal, the application provides a 2.5D integrated circuit test path optimization method based on a decomposition multi-objective evolutionary algorithm, which reduces test time and test hardware overhead, thereby reducing the comprehensive test cost. The application designs a multi-objective evolutionary algorithm based on a decomposition strategy to balance test time and hardware overhead, thereby reducing the comprehensive test cost and realizing efficient optimization of the 2.5D integrated circuit test path.

[0009] The 2.5D integrated circuit test path optimization method based on the decomposition multi-objective evolutionary algorithm comprises the following steps:

[0010] Step one, initialization, generating an initial population for the test path configuration scheme of the 2.5D integrated circuit, and evaluating the initial population index according to the multi-objective evolutionary algorithm, and selecting the optimal fitness index as the ideal point of the initial population;

[0011] Step two, generating a new solution by using a crossover and mutation operator, and evaluating the new solution and updating the ideal point according to the multi-objective evolutionary algorithm

[0012] Step three, environment selection, including updating the neighborhood solution and the elite solution set;

[0013] Step four, executing the iteration process of steps two and three for multiple times until the iteration ends, outputting the elite solution set, and the elite solution set is a non-dominated solution, which is a set of multiple optimal test path configuration schemes, and a decision maker selects an arbitrary optimal test path configuration scheme from the elite solution set for testing.

[0014] Preferably, the process of step one initialization is:

[0015] Step one, determining the maximum iteration number G of the multi-objective evolutionary algorithm max , the population size NP, the uniformly distributed weight vector set λ=[λ1, λ2, …, λ NP ], np=1, 2, …, NP, np is the serial number of the individual in the population, the number of neighbor weights T, and the elite solution set EP;

[0016] Step one, calculating the Euclidean distance between any two weight vectors, determining T nearest neighbor weight vectors for each weight vector, i.e., B(np)={np1, np2, …, np T}, nb=np1, np2, …, np T , and the corresponding respectively represent T neighbor weights λ np of the weight vector λ nb ;

[0017] Step three, randomly generate initial population in feasible region satisfying constraint condition: s1, s2, …, s NP Each individual represents a test path configuration scheme, and there are NP test path configuration schemes, wherein each individual is coded as:

[0018] s np ={x 11 ,x 12 ,…,x 1m ,…x dm ,y 11 ,y 12 ,…,y 1n ,…,y dn},

[0019] Wherein, x ij represents the inclusion relationship between the internal scan chain of the wafer and the in-TAM, if the input scan chain of the i-th wafer is contained in the j-th in-TAM, then the value of x ij is 1, otherwise 0;

[0020] y ik represents the inclusion relationship between the external scan chain of the wafer and the out-TAM, if the input scan chain of the i-th wafer is contained in the k-th out-TAM, then the value of y ik is 1, otherwise 0;

[0021] i represents the serial number of the wafer, i = 1, 2, …, d, d is the number of wafers; j represents the serial number of the in-TAM, j = 1, 2, …, m, m represents the number of in-TAMs, k represents the serial number of the out-TAM, k = 1, 2, …, n, n represents the number of out-TAMs, and there is a relationship: m ≤ d, n ≤ d;

[0022] Step four, evaluate the fitness value of the initial population: F(s np ) = [F1(s np ), F2(s np )]

[0023] Wherein, F1(s np ) is the test time cost of individual s np , which is optimized as follows:

[0024]

[0025] In the formula, c1 is the test cost of the automatic test equipment on the unit length test path, c ATE is the test cost per second, f is the test frequency, and N is the total number of interconnections to be tested;

[0026] I iNin(i) represents the number of micro-bumps connected to the input port of the ith chip; O i Nout(i) represents the number of micro-bumps connected to the output port of the ith chip; and there is a constraint:

[0027]

[0028]

[0029] L represents the length of the longest test path in all in-TAMs and out-TAMs;

[0030] F2(s np ) is the hardware overhead cost of individual s np , which is optimized as:

[0031]

[0032] where c2 is the hardware overhead caused by adding one extra in-TAM, c3 is the hardware overhead caused by adding one extra out-TAM, S TSV represents the area of through-silicon via, c interposer represents the hardware overhead cost of per unit area interconnection, S bump represents the area of micro-bump, c die represents the hardware overhead cost of per unit area chip;

[0033] Step one five, initialize ideal point z * = [z1, z2], where z1 and z2 are the minimum fitness function values on each optimization objective respectively, and are calculated as:

[0034]

[0035]

[0036]

[0037] Preferably, the specific process of step two is as follows:

[0038] Step two one, for each weight vector λ np , np = 1, 2, …, NP, randomly select two weight vector indexes a, b ∈ B(np) in its neighbor weight vector set B(np), the corresponding weight vectors are λ a and λ b , the corresponding individuals are s a and s b , and a new solution y is generated by using crossover and mutation operators;

[0039] ​Step two two, the continuous encoding new solution y generated in step two one is repaired and converted into discrete encoding, the elements less than 0.5 in the new solution y are set to 0, and the elements greater than or equal to 0.5 are set to 1;

[0040] Step two three, the new solution index is evaluated according to the multi-objective evolutionary algorithm, and the fitness index F of the new solution generated in this iteration is calculated l (y), l = 1, 2, and the ideal point z is updated * :

[0041] If F l (y) < z l , then z l is updated = F l (y), otherwise the ideal point of the last iteration is retained.

[0042] Preferably, the specific process of step three is as follows:

[0043] Step three one, the neighborhood solution is updated;

[0044] For each neighborhood solution s nb of the new solution y and the np th individual, nb ∈ B(np), the scalar function values of the neighborhood solution s nb and the new solution y under Chebyshev decomposition are compared, if g tch (y|λ nb ,z * ) ≤ g tch (s nb |λ nb ,z * ), then the neighborhood solution is replaced by the new solution, that is, s nb = y, otherwise the result of the last iteration is retained;

[0045]

[0046]

[0047] Step three two, the elite solution set is updated, first, all the individuals in EP dominated by the new solution y are removed from EP, and then the new solution y is added to EP.

[0048] Preferably, in step four, it is checked whether the maximum number of iterations G max allowed is reached, if it is reached, the elite solution set EP is output, otherwise, it is jumped to step two.

[0049] Advantages of the present application:

[0050] (1) The decomposition-based multi-objective evolutionary algorithm is an efficient heuristic optimization technique, which converts the multi-objective optimization problem into multiple simple optimization sub-problems by using decomposition technology, and performs collaborative optimization by using neighborhood relationship, so as to obtain the approximate solution of the entire Pareto front.

[0051] (2) The decomposition-based multi-objective evolutionary algorithm is used for optimizing the test path of the 2.5D integrated circuit, which can effectively balance the test time cost and hardware overhead cost, determine the optimal test path, and reduce the overall test cost. BRIEF DESCRIPTION OF DRAWINGS

[0052] Figure 1 is a test path configuration method of a 2.5D integrated circuit of three chips;

[0053] Figure 2 is a basic test path schematic diagram;

[0054] Figure 3 is a flowchart of the decomposition-based multi-objective evolutionary algorithm for optimizing the test path of the 2.5D integrated circuit. DETAILED DESCRIPTION

[0055] The technical solutions in the embodiments of the present application will be described clearly and completely below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor fall within the scope of protection of the present application.

[0056] It should be noted that the embodiments in the present application and the features in the embodiments can be combined with each other without conflict.

[0057] The present application will be further described below with reference to the drawings and specific embodiments, but is not limited by the present application.

[0058] Specific embodiment one: the present embodiment will be described below with reference to the drawings Figure 3 The decomposition-based multi-objective evolutionary algorithm for optimizing the test path of the 2.5D integrated circuit, which includes the following steps:

[0059] Step one, initialization, generating an initial population for the test path configuration scheme of the 2.5D integrated circuit, and evaluating the initial population index according to the multi-objective evolutionary algorithm, and selecting the optimal fitness index as the ideal point of the initial population;

[0060] Step two, generating a new solution by using the crossover and mutation operators, and evaluating the new solution and updating the ideal point according to the multi-objective evolutionary algorithm

[0061] Step 3: Environment selection, including updating the neighborhood solution and elite solution set;

[0062] Step 4: Repeat the iteration process of Step 2 and Step 3 multiple times until the iteration ends, and output the elite solution set. The elite solution set is a non-dominated solution, which is a set of multiple optimal test path configuration schemes. The decision-maker selects any one of the optimal test path configuration schemes from the elite solution set for testing.

[0063] The optimization problem of this implementation is defined as follows: Given a 2.5D integrated circuit with a total of d wafers, assume that the testing cost of the automatic test equipment on a test path of unit length is c1, the hardware overhead caused by adding an additional in-TAM is c2, and the hardware overhead caused by adding an additional out-TAM is c3. The final optimization objective is to determine the optimal test path scheme to minimize the testing cost C.

[0064] The initialization process in step one is as follows:

[0065] Step 11: Determine the maximum number of iterations G for the multi-objective evolutionary algorithm. max Population size NP, uniformly distributed weight vector set λ = [λ1, λ2, ..., λ NP ], np = 1, 2, ..., NP, where np is the index of an individual in the population, T is the number of neighbor weights, and EP is the elite solution set;

[0066] Steps 1 and 2: Calculate the Euclidean distance between any two weight vectors, and determine T nearest neighbor weight vectors for each weight vector, i.e., B(np) = {np1, np2, ..., np}. T}, nb = np1, np2, ..., np T Corresponding These represent the weight vectors λ and λ, respectively. np The weights λ of the T neighbors nb ;

[0067] Step 13: Randomly generate an initial population within the feasible region that satisfies the constraints: s1, s2, ..., s NP Each individual represents a test path configuration scheme, with a total of NP test path configuration schemes. Each individual is coded as follows:

[0068] s np ={x 11 ,x 12 ,…,x 1m ,…x dm ,y 11 ,y 12 ,…,y 1n ,…,y dn},

[0069] wherein x ij represents the inclusion relationship between the inner scan chain of the wafer and the in-TAM, the input scan chain of the i-th wafer is included in the j-th in-TAM, then x ij = 1, otherwise 0;

[0070] y ik represents the inclusion relationship between the outer scan chain of the wafer and the out-TAM, the input scan chain of the i-th wafer is included in the k-th out-TAM, then y ik = 1, otherwise 0;

[0071] i represents the serial number of the wafer, i = 1, 2,..., d, d is the number of wafers; j represents the serial number of the in-TAM, j = 1, 2,..., m, m represents the number of in-TAMs, k represents the serial number of the out-TAM, k = 1, 2,..., n, n represents the number of out-TAMs, and there is a relationship: m≤d,n≤d;

[0072] Step one four, fitness value evaluation is carried out on the initial population: F(s np ) = [F1(s np ), F2(s np )]

[0073] wherein F1(s np ) is the test time cost of the individual s np , which is optimized according to the following formula:

[0074]

[0075] In the formula, c1 is the test cost of the automatic test equipment on the unit length test path, c ATE is the test cost per second, f is the test frequency, and N is the total number of interconnections to be tested;

[0076] I i represents the number of micro-soldering balls connected to the input port of the i-th wafer; O i represents the number of micro-soldering balls connected to the output port of the i-th wafer; and there is a constraint condition:

[0077]

[0078]

[0079] L represents the length of the longest test path in all in-TAMs and out-TAMs;

[0080] F2(s np ) is the hardware overhead cost of the individual s np , which is optimized according to the following formula:

[0081]

[0082] where c2 is the hardware overhead of adding one extra in-TAM, c3 is the hardware overhead of adding one extra out-TAM, S TSV represents the area of through-silicon via, c interposer represents the hardware overhead cost per unit area of interconnection, S bump represents the area of micro-bump, c die represents the hardware overhead cost per unit area of wafer;

[0083] In summary, the optimization objective function is F = [F1, F2], and the present embodiment is a two-objective optimization problem. The balance point of the two objectives will be solved and optimized as follows.

[0084] Step one five, initialize ideal point z * = [z1, z2], where z1 and z2 are the minimum fitness function values on each optimization objective, and are calculated according to the formula

[0085]

[0086]

[0087] .

[0088] The specific process of step two is as follows:

[0089] Step two one, for each weight vector λ np , np = 1, 2, …, NP, randomly select two weight vector indexes a, b ∈ B(np) in its neighbor weight vector set B(np), and the corresponding weight vectors are λ a and λ b , and the corresponding individuals are s a and s b , and a new solution y is generated by using crossover and mutation operators;

[0090] Step two two, repair and convert the continuous coding new solution y generated in step two one into discrete coding, and set the elements less than 0.5 in the new solution y to 0 and the elements greater than or equal to 0.5 to 1;

[0091] Step two three, evaluate the new solution index according to the multi-objective evolutionary algorithm, calculate the fitness index F l (y), l = 1, 2, of the new solution generated in this iteration, and update the ideal point z * :

[0092] If F l (y) < z l , then update zl =F l (y), otherwise retain the ideal point from the previous iteration.

[0093] In each iteration, the new solution is compared with the individual corresponding to the ideal point in the previous iteration, and the one with the smaller fitness value is retained.

[0094] In each iteration, in addition to determining the new solution in step two, the set of domain solutions and elite solutions also needs to be updated, as detailed in step three.

[0095] The specific process of step three is as follows:

[0096] Step 31: Update the neighborhood solution;

[0097] For each neighborhood solution s of the new solution y and the np-th individual nb For nb∈B(np), compare the neighborhood solutions s nb And the scalar function value of the new solution y under Chebyshev decomposition, if g tch (y|λ nb ,z * )≤g tch (s nb |λ nb ,z * If the neighborhood solution is replaced by the new solution, then s nb =y, otherwise retain the result of the previous iteration without replacing;

[0098] The scalar function value of the new solution y under Chebyshev decomposition is:

[0099]

[0100] Neighborhood solution s nb The scalar function value under Chebyshev decomposition is:

[0101]

[0102] Step 3.2: Update the elite solution set. First, remove all individuals in the EP that are dominated by the new solution y from the EP, and then add the new solution y to the EP.

[0103] In step four, check whether the maximum allowed number of iterations G has been reached. max If the condition is met, output the elite solution set EP; otherwise, proceed to step two.

[0104] The elite solutions EP obtained through the above steps are all non-dominated solutions, which are the set of optimal test paths found. Decision-makers can select an optimal test path from this set to reduce testing costs. The above process can be used... Figure 3 The flowchart shown is an example.

[0105] While the application has been described with reference to particular embodiments, it is to be understood that the application is not limited to the particulars disclosed. Rather, it is a continuation of the principles and applications of the present application. It is therefore to be understood that numerous modifications, both as to the details and embodiments illustrated and the application, can be made by those skilled in the art without departing from the spirit and scope of the application as defined by the appended claims. It should be understood that all the features described in connection with the various embodiments can be combined in other combinations than those explicitly described. It should also be understood that features described in connection with one embodiment can be combined with features described in connection with another embodiment.

Claims

1. A 2.5D integrated circuit test path optimization method based on a decomposition multi-objective evolutionary algorithm, characterized in that, The method comprises the following steps: Step one, initialization, generating an initial population for the test path configuration scheme of the 2.5D integrated circuit, and evaluating the initial population index according to the multi-objective evolutionary algorithm, and selecting the optimal fitness index as the ideal point of the initial population; Step two, generating a new solution by using a crossover and mutation operator, and evaluating the new solution and updating the ideal point according to the multi-objective evolutionary algorithm; Step three, environment selection, including updating the neighborhood solution and the elite solution set; Step four, executing the iteration process of steps two and three for multiple times until the iteration ends, outputting the elite solution set, the elite solution set being a non-dominated solution and being a set comprising multiple optimal test path configuration schemes, and a decision maker selecting an arbitrary optimal test path configuration scheme from the elite solution set for testing; The process of step one initialization is as follows: Step one, determine the maximum number of iterations G of the multi-objective evolutionary algorithm max , population size NP, uniformly distributed weight vector set λ = [λ1, λ2, …, λ NP ]np = 1, 2, …, NP, np is the serial number of individuals in the population, the number of neighbor weights T, elite solution set EP; Step one two, calculate the Euclidean distance between any two weight vectors, determine T nearest neighbor weight vectors for each weight vector, that is, the neighbor weight vector set B(np) = {np1, np2, …, np T}, nb= np1, np2, …, np T , corresponding respectively represent the T neighbors of the weight vector λ np λ nb ; Step one three, randomly generate initial population in the feasible region that meets the constraint condition: s1, s2, …, s NP Each individual represents a test path configuration scheme, and there are NP test path configuration schemes in total, wherein each individual is encoded as: s np = {x 11 ,x 12 ,…,x 1m ,…x dm ,y 11 ,y 12 ,…,y 1n ,…,y dn} where x ij The inclusion relationship between the inner scan chain of a wafer and in-TAM, the input scan chain of the ith wafer is included in the jth in-TAM, then x ij is 1, otherwise 0; y ik The inclusion relationship of the outer scan chain of the wafer and the out-TAM, the input scan chain of the i-th wafer is included in the k-th out-TAM, then y ik The value is 1, otherwise 0. i represents the serial number of the wafer, i = 1, 2,..., d, d being the wafer quantity; j represents the in-TAM serial number, j = 1, 2,..., m, m representing the in-TAM quantity, k representing the out-TAM serial number, k = 1, 2,..., n, n representing the out-TAM quantity, and there being a relationship formula: m ≤ d, n ≤ d; Step one four, fitness value evaluation for initial population: F(s np ) = [F1(s np ), F2(s np )] where F1(s np ) is the testing time cost of individual s np , which is optimized as follows: wherein c1 is the test cost per unit length of test path for the automated test equipment, c ATE is the test cost per second, f is the test frequency, and N is the total number of interconnections to be tested. I i represents the number of micro-bumps connected to the input port of the i-th die; O i represents the number of micro-bumps connected to the output port of the i-th die; and there are the constraints: L represents the length of the longest test path in all in-TAMs and out-TAMs; F2(s np ) is the hardware overhead cost of individual s np is optimized as follows: where c2 is the hardware overhead due to adding one extra in-TAM, c3 is the hardware overhead due to adding one extra out-TAM, S TSV represents the area of a through-silicon via, c interposer represents the hardware overhead cost per unit area of interconnect, S bump represents the area of a micro-bump, c die represents the hardware overhead cost per unit area of wafer; Step one five, initializing ideal point z * = [zl, z2], where zl and z2 are the minimum fitness function values on each optimization objective, respectively, and are given by Step two is specifically as follows:

2. The method of claim 1, wherein the method is based on a decomposition-based multi-objective evolutionary algorithm. Step two one, generating a new solution y of continuous coding; Step two, for each weight vector λ np np = 1,2,..., NP, randomly pick two weight vector indexes a, b ∈ B(np) in its neighbor weight vector set B(np), the corresponding weight vectors are λ a and λ b , the corresponding individuals are s a and s b , generate new solution y by using crossover and mutation operators; Step two two, repairing and converting the new solution y of continuous coding into discrete coding, setting the elements less than 0.5 in the new solution y to 0, and setting the elements greater than or equal to 0.5 to 1; Step two three, according to the multi-objective evolutionary algorithm evaluation of new solution index, the new solution generated this iteration fitness index F l (y), l = 1, 2, and update the ideal point z * : If F l (y) < z l , then update z l = F l (y), otherwise keep the ideal point of the last iteration.

3. The method of claim 2, wherein the method is based on a decomposition-based multi-objective evolutionary algorithm. Step three is specifically as follows: Step three one, updating the neighborhood solution; for each neighborhood solution s of the np-th individual nb , nb e B(np), compare the neighborhood solution s nb and the scalar function value of the new solution y under Chebyshev decomposition, if g tch (y | l nb , z * ) < g tch (s nb | l nb , z * ), then replace the neighborhood solution by the new solution, i.e. s nb = y, otherwise keep the result of the last iteration; Step three two, updating the elite solution set, first removing all individuals dominated by the new solution y from the EP, and then adding the new solution y to the EP.

4. The method of claim 1, wherein the method is based on a decomposition-based multi-objective evolutionary algorithm. Step four checks if the maximum number of iterations G is reached max If yes, the elite set EP is output, otherwise go to step two.