Fault Reconstruction for a Giant Satellite Swarm Based on Hybrid Multi-Objective Optimization

: To perform indicator selection and veriﬁcation for the on-orbit fault reconstruction of a giant satellite swarm, a hybrid multi-objective fault reconstruction algorithm is proposed and then veriﬁed by Monte Carlo analysis. First, according to the on-orbit failure analysis of the satellite swarm, several optimization indicators, such as the health state of the satellite swarm, the total energy consumption of reconstruction, and the balance of fuel consumption, are proposed. Then, a hybrid multi-objective ﬁtness function is constructed, and a hybrid multi-objective genetic algorithm is used to optimize the objective function to obtain the optimal reconstruction strategy. Finally, the algorithm is statistically veriﬁed by Monte Carlo analysis. The simulation results not only show the algorithm’s validity but also reveal the relationship between the number of satellite faults and the health of the satellite swarm. From this, the maximum number of faulty satellites allowed in the giant satellite swarm is calculated, which is signiﬁcant for assessing the swarm’s health.


Introduction
In recent years, satellite systems have shown a trend of clustering [1]. For example, SpaceX plans to launch tens of thousands of low-Earth orbit satellites to achieve global high-speed internet coverage as part of its Starlink project [2]. The rapid increase in microand nanosatellites has made many tasks easier. Moreover, system complexity and fault risk have increased, which has led to higher requirements for the on-orbit maintenance and health management of ultra-large-scale satellite swarm systems [3]. The definition of an ultra-large-scale satellite cluster is not consistent throughout academia. To avoid confusion, this paper refers to satellite clusters consisting of more than 1000 satellites as ultra-large-scale satellite swarms.
To improve the ability of ultra-large-scale satellite swarms to handle and tolerate faults autonomously, scholars have carried out research on satellite swarm reconstruction technology. Yang H et al. [4] proposed a decentralized formation control scheme based on couplings within the satellite system and maintained the stable operation of the entire tethered satellite system by adjusting the controller of the faulty satellite itself. Zhang C et al. [5] considered the problem of spacecraft formation reconfiguration under directed and undirected communication. They proposed that the faulty satellite and the reconfigured healthy satellite cooperate to compensate for controller faults and synchronize the system attitude. Li J et al. [6] introduced a fuzzy logic system to solve the problem of tracking inconsistency among multiple satellites under fault conditions by cutting off the communication between faulty satellites and healthy satellites. From the perspective of fault-tolerant control, the above scheme focuses on the situation in which a single satellite fault has little impact on the entire system. To study the problem of mission failure caused by multi-satellite failure, Jiang B et al. [7] proposed a state-switching multibody satellite formation failure reconstruction method, which has faulty satellites leave the formation and healthy satellites compensate to achieve reconstruction.
The above fault reconstruction techniques have been developed over many years. However, for ultra-large-scale satellite clusters, factors such as the complexity of satellite coordination, task diversity, and uncertainty of reconstruction indicators introduce challenges to the satellite swarm reconstruction task. Since fault satellite swarm reconstruction is usually based on system redundancy and reconfiguration of the system structure [8], there may be multiple reconstruction control strategies. However, for different satellite swarm systems, there are differences in their reconstruction strategies in terms of energy consumption, health state, reliability, etc. Therefore, the problem of satellite swarm fault reconstruction is a kind of multi-objective optimization problem. Multi-objective optimization problems play an important role both in theory and in engineering applications. In the 1960s, R.S. Rosenberg [9] first tried to use genetic algorithms to solve multi-objective optimization problems, and gradually multi-objective optimization problems came to the public's attention. In the 1980s, J.D. Schaffer [10] designed the first multi-objective evolutionary algorithm, the vector-evaluated genetic algorithm. Since then, the multi-objective optimization algorithm has gone through two periods, from the initial development with simple rules to the current second stage with efficiency as the evolutionary goal. According to the different ways of solving multi-objective problems, refactoring methods can be divided into system-model-based and transformation-based methods. The former needs to establish an accurate system model and use inference diagnosis to solve the strategy, which cannot deal with large and complex systems and is not efficient. The latter converts the problem to other well-studied domain problems, such as the planning field [11][12][13].
Multi-objective optimization algorithms [14,15] are usually used to solve multi-objective problems. These algorithms have advantages and disadvantages. For instance, the multiobjective evolutionary algorithm based on decomposition (MOEA/D) [16,17] converges quickly, but the effect for high-dimensional multi-objective problems is poor. Multiobjective particle swarm optimization (MOPSO) algorithms [18][19][20] cannot solve discrete optimization problems. The nondominated sorting genetic algorithm-II (NSGA-II) [21][22][23] has low computational complexity and high population diversity and is widely used to solve multi-objective optimization problems. However, there is no unique solution to a given multi-objective optimization problem; that is, the goal of multi-objective algorithms is to find a global optimal solution on the Pareto front [24][25][26][27]. Therefore, for reconstruction problems that require a unique optimal solution, other solutions need to be considered. There are three common methods: (1) the generative method, i.e., generating the optimal solution set, which requires the decision maker to filter out the optimal solution according to the needs of the problem. (2) The interactive method finds the optimal solution through dialogue and interaction between the decision-maker and the analyst. (3) The transformation method requires the decision maker to judge the importance of the objective function, transforms the multi-objective problem into a single-objective problem for a solution and often uses an intelligent algorithm.
This paper takes the self-assembled space telescope [28] as the mission setting and introduces optimization indicators such as the health state of the satellite swarm, the total energy consumption of the reconstruction, and the balance of fuel consumption. Based on this framework, a reconstruction optimization model for each stage is constructed, and a hybrid multi-objective genetic algorithm [29,30] is used to perform the reconstruction. The core advantage of the method in this paper is that even if the satellite faults in the satellite swarm system are different, it still provides a multi-objective optimal reconstruction strategy.

Construction of On-Orbit Reconstruction Indicators for a Giant Satellite Swarm
Fault reconfiguration refers to the process in which the satellite swarm maintains a healthy state by adjusting the positions of healthy satellites to replace faulty satellites or by isolating faulty satellites to limit their influence on the swarm. In essence, reconfiguration Appl. Sci. 2023, 13, 6674 3 of 13 involves a change in the satellite swarm topology; for example, position replacement results in a configuration change, and fault isolation results in a change in the information path. In some scenarios, the latter is also an independent satellite swarm fault isolation problem. The core issue of reconstruction is how to select reconstruction indicators. Starting from actual engineering, this paper studies the reconstruction of the swarm configuration. It constructs mathematical models for properties such as the satellite swarm health state, the total energy consumption of reconstruction, and the balance of fuel consumption, which lays the foundation for the selection of reconstruction indicators.

Model Assumptions of Satellite Swarm Failure
Satellite swarm fault reconstruction involves multi-satellite deployment. Without loss of generality, it is assumed that the satellites participating in the reconstruction are near the faulty satellite swarm, and the faulty satellites retain deorbiting capabilities. The entire satellite swarm failure reconstruction process is shown in Figure 1.

Construction of On-Orbit Reconstruction Indicators for a Giant Satellite Swarm
Fault reconfiguration refers to the process in which the satellite swarm mainta healthy state by adjusting the positions of healthy satellites to replace faulty satellit by isolating faulty satellites to limit their influence on the swarm. In essence, reconfig tion involves a change in the satellite swarm topology; for example, position replace results in a configuration change, and fault isolation results in a change in the inform path. In some scenarios, the latter is also an independent satellite swarm fault isol problem. The core issue of reconstruction is how to select reconstruction indicators. S ing from actual engineering, this paper studies the reconstruction of the swarm con ration. It constructs mathematical models for properties such as the satellite swarm h state, the total energy consumption of reconstruction, and the balance of fuel consu tion, which lays the foundation for the selection of reconstruction indicators.

Model Assumptions of Satellite Swarm Failure
Satellite swarm fault reconstruction involves multi-satellite deployment. Wit loss of generality, it is assumed that the satellites participating in the reconstructio near the faulty satellite swarm, and the faulty satellites retain deorbiting capabilities entire satellite swarm failure reconstruction process is shown in Figure 1. Self-assembled space telescopes are representative of mission-oriented multi-sat scenarios. Its health and stability will impact the space telescope's imaging mission. T fore, this paper studies the self-assembled mirror satellite swarm of the space telesc To simplify the decision-making problem of reconfiguration and allocation of sat clusters, a two-dimensional schematic diagram shows the status of satellite swarm re figuration. The three main states of the satellite swarm are shown in Figure 2  Self-assembled space telescopes are representative of mission-oriented multi-satellite scenarios. Its health and stability will impact the space telescope's imaging mission. Therefore, this paper studies the self-assembled mirror satellite swarm of the space telescope. To simplify the decision-making problem of reconfiguration and allocation of satellite clusters, a two-dimensional schematic diagram shows the status of satellite swarm reconfiguration. The three main states of the satellite swarm are shown in Figure 2: Figure  Considering the characteristics of the satellite swarm, the following assumptions are made for the model: 1. The mission-oriented satellite swarm is composed of individual star satellites that are connected via inter-satellite communication links. Considering the characteristics of the satellite swarm, the following assumptions are made for the model: The mission-oriented satellite swarm is composed of individual star satellites that are connected via inter-satellite communication links. 2.
The number of satellites in the satellite swarm is N, each satellite is of the same type, and each satellite is a regular hexagon [28] with a side length of L.

3.
The satellites in the satellite swarm are evenly distributed, and the distance between adjacent satellites is a, which is much smaller than the orbital radius r, i.e., a r.

4.
The satellites are numbered clockwise from the virtual center of the satellite swarm, and the inter-satellite distance correlation matrix is D n×n .

Relative Motion Modeling
To better describe the position and velocity of the satellite, a coordinate system R is established with the virtual center O t of the satellite swarm as the origin. The X t axis is perpendicular to the orbit where O t is located; the Y t axis is perpendicular to the X t axis and tangent to the O t orbit, and the direction of movement is positive. The Z t axis is perpendicular to the orbital plane and forms a right-handed coordinate system with the X t and Y t axes. In this coordinate system, the reconstructed satellite CS moves toward the faulty satellite TS, as shown in Figure 3. Considering the characteristics of the satellite swarm, the following assumptions are made for the model: 1. The mission-oriented satellite swarm is composed of individual star satellites that are connected via inter-satellite communication links.
2. The number of satellites in the satellite swarm is N , each satellite is of the same type, and each satellite is a regular hexagon [28] with a side length of L .
3. The satellites in the satellite swarm are evenly distributed, and the distance between adjacent satellites is a , which is much smaller than the orbital radius r , i.e., a r  .
4. The satellites are numbered clockwise from the virtual center of the satellite swarm, and the inter-satellite distance correlation matrix is n n × D .

Relative Motion Modeling
To better describe the position and velocity of the satellite, a coordinate system R is established with the virtual center t O of the satellite swarm as the origin. The t X axis is perpendicular to the orbit where t O is located; the t Y axis is perpendicular to the t X axis and tangent to the t O orbit, and the direction of movement is positive. The t Z axis is perpendicular to the orbital plane and forms a right-handed coordinate system with the t X and t Y axes. In this coordinate system, the reconstructed satellite CS moves toward the faulty satellite TS, as shown in Figure 3. Due to the small relative distance of the satellite swarms, the motion of the satellite swarms satisfies the C-W equation [31]. Then, the satellite reconstruction problem is: Due to the small relative distance of the satellite swarms, the motion of the satellite swarms satisfies the C-W equation [31]. Then, the satellite reconstruction problem is: where x, y and z are the position components of the reconstructed satellite in the R coordinate system; u x , u y , and u z are the orbit control accelerations of the reconstructed satellite; and n is the orbital angular velocity of the faulty satellite.
Let the state sequence be x(t) = x y z .
x . y . z T and the orbital acceleration be u(t) = u x u y u z T . The relative motion Equation (1) is written in matrix form to obtain Equation (2): where According to optimal control theory, considering the finite time, the optimal energy is: In Equation (3), t 0 and t f are the initial and final times, respectively. The energyoptimal position-velocity equation is as follows: In the formula, c 1 ∼ c 12 are constants, which can be obtained by substituting the known position and speed states of CS and TS into Equation (4); in addition, a = sin nt and b = cos nt .

The Introduction of the Reconstruction Indicators
Because satellites carry limited fuel, fuel consumption, and energy consumption need to be considered during on-orbit reconstruction. Meanwhile, the health state also needs to be considered, so the following parameters are used to establish a reconstruction indicator system: (1) Satellite swarm health state The satellite swarm health state refers to the durability of individual satellites and topological structures in the satellite swarm and the ability of the satellite swarm to meet mission requirements [32]. The satellite swarm health state is an important indicator that measures the overall health of the satellite swarm and detects faults effectively. The health state is used to measure the effect of satellite swarm reconstruction. Considering the nonlinear composition characteristics of the satellite swarm, such as the multiple dimensions of a single satellite, inter-satellite communication link, and mission effectiveness in the satellite swarm, the health state at time H swarm is expressed as: In the formula, H sat is the weighted sum of the health states of all individual satellites in the swarm, where the single-satellite reliability of satellite i is R sat_i and its weight is ω sat_i ; in the same way, H link is the weighted sum of the health status of the inter-satellite communication links, where the link reliability of link i is R link_i and its weight is ω link_i . k 1 , k 2 , k 3 are the influencing parameters, which represent the relative influence of each part of the satellite swarm system. E swarm is the satellite swarm performance, indicating the ability of the satellite swarm to achieve a task. Different tasks have different forms, which are determined by the specific circumstances of the task. This paper considers a satellite swarm that performs space telescope missions. According to the requirements of optical imaging, the fewer faulty satellites are located near the center of the mirror, the better the imaging performance, i.e., the number of effective satellites can be used to measure the satellite swarm performance. Therefore, the satellites are divided into different layers according to the distance from the satellite to the center of the surface. Satellite swarm effectiveness is determined by the number of working satellites in each layer and the weight values of each layer: In the formula, ω sat_i is the weight value of layer i; N layer_i indicates the number of nonfaulty satellites contained in the i-th layer and satisfies the number of nonfaulty satellites in the satellite swarm.
(2) Total energy consumption of reconstruction The total energy consumption of reconstruction is the sum of the energy consumed by all satellites participating in the reconstruction to complete orbital maneuvers. Due to the limited fuel carried by each satellite, to satisfy the satellite swarm's requirements for operating life and mission performance, the lower the total energy consumption is, the better. Therefore, the total energy consumption of reconstruction is expressed as: In this formula, J i indicates the energy consumption of the orbital maneuver of the i-th satellite participating in the reconstruction, and k indicates the number of satellites participating in the reconstruction.
(3) Fuel consumption balance [33] To reduce the reconstruction time, the principle of nearest replacement is adopted. However, the remaining fuel of adjacent satellites may be insufficient, resulting in excessive fuel consumption of the satellite, a shorter working life, and negative effects on the health and mission performance of the satellite swarm. The constraint of fuel consumption balance is introduced to avoid such a situation. The variance of the satellite speed increment involved in the reconstruction is used to characterize the fuel consumption balance: In the formula, J is the average energy consumption of all nonfaulty satellites participating in the reconstruction.

Fault Reconstruction Based on a Hybrid Multi-Objective Genetic Algorithm
After the reconstruction indicator system is established, the reconstruction strategy must be determined. Since the satellite numbers are discontinuous, the process requires optimizing a discrete system. Commonly used multi-objective optimization algorithms obtain a set of optimal solutions at the Pareto front for continuous variables, so they are unsuitable for this research scenario. Therefore, consistent with the optimization indicators established in the previous section, the linear weighted sum method is adopted, and a utility function is formed to represent the overall optimization goal after normalization; then, the reconstruction strategy is determined by a genetic algorithm. In summary, the decision tree of multi-objective fault reconstruction is shown in Figure 4.

Genetic Coding for Satellite Numbering
Suppose the number of satellites in a given swarm is N , and the set of numbers is S = {s 1 , s 2 , · · · , s N }, s i ∈ N. At a certain moment, N F faulty satellites are detected, the set of which is S F = s F1 , s F2 , · · · , s FN F , s Fi ∈ S. Then, the set of nonfaulty satellites is S H = s H1 , s H2 , · · · , s HN H , s Hi ∈ (S − S F ), and N H = N − N F .
Since it is necessary to determine the number of satellites participating in the reconstruction, integer coding is used. The decision variable X = [x 1 , x 2 , · · · , x N F ] is a vector composed of satellite numbers participating in the reconstruction. The genetic code is a one-to-one correspondence between the x i and s Fi , which is shown in Figure 5: obtain a set of optimal solutions at the Pareto front for continuous variables, so they ar unsuitable for this research scenario. Therefore, consistent with the optimization indica tors established in the previous section, the linear weighted sum method is adopted, an a utility function is formed to represent the overall optimization goal after normalization then, the reconstruction strategy is determined by a genetic algorithm. In summary, th decision tree of multi-objective fault reconstruction is shown in Figure 4.

Genetic Coding for Satellite Numbering
Since it is necessary to determine the number of satellites participating in the reconstruction, integer coding is used. The decision variable 1 2 [ , , , ] is a vector composed of satellite numbers participating in the reconstruction. The genetic code is a one-to-one correspondence between the i x and Fi s , which is shown in Figure 5:   (9) In the formula,

Construction of Hybrid Multi-Objective Functions
The total energy consumption and fuel consumption balance must be made as small as possible, while the satellite swarm health must be made as large as possible. According to the gray clustering [34] weight determination method, weight values are assigned to the optimization indicators, and the multi-objective problem is transformed into a singleobjective problem. Since where ( ) i f x is the i-th indicator function. Therefore, the evaluation function ( ) F X can be expressed as: The detailed flow of the hybrid multi-objective genetic algorithm is shown in Figure   Figure 5. Schematic diagram of the genetic code.
To prevent the repetition of variable x i in the process of cross-mutation, the process draws on the principle of solving the traveling salesman problem. During initialization, the domain of X is set to S H , which can be expressed as: In the formula, i = 1, 2, · · · , N F .

Construction of Hybrid Multi-Objective Functions
The total energy consumption and fuel consumption balance must be made as small as possible, while the satellite swarm health must be made as large as possible. According to the gray clustering [34] weight determination method, weight values are assigned to the optimization indicators, and the multi-objective problem is transformed into a singleobjective problem.
Since the units of the three indicators are different, normalization and standardization are used to obtain the normalized indicator function f i (x) * : where f i (x) is the i-th indicator function. Therefore, the evaluation function F(X) can be expressed as: The detailed flow of the hybrid multi-objective genetic algorithm is shown in Figure 6.

Simulation and Analysis
According to the model in Section 2.1, fault reconstruction simulation analysis is carried out first to verify the effectiveness of the reconstruction algorithm; then, Monte Carlo simulation is used to verify the correctness of the reconstruction algorithm.

Fault Reconstruction Simulation
According to the model assumptions in Section 2.1, the faulty satellite distribution is represented as a vector F S ; the orbital radius of the virtual center of the satellite swarm is s R , and Earth's gravitational constant is μ . The reconstruction model parameters are shown in Table 1. The three-axis coordinate parameters of the faulty satellite in the relative coordinate system are shown in Table 2. The simulated running orbit is nearly circular, so the angular velocities are approximately equal. The initial values of the two moving satellites in the relative system are shown in Table 3.

Simulation and Analysis
According to the model in Section 2.1, fault reconstruction simulation analysis is carried out first to verify the effectiveness of the reconstruction algorithm; then, Monte Carlo simulation is used to verify the correctness of the reconstruction algorithm.

Fault Reconstruction Simulation
According to the model assumptions in Section 2.1, the faulty satellite distribution is represented as a vector S F ; the orbital radius of the virtual center of the satellite swarm is R s , and Earth's gravitational constant is µ. The reconstruction model parameters are shown in Table 1. The three-axis coordinate parameters of the faulty satellite in the relative coordinate system are shown in Table 2. The simulated running orbit is nearly circular, so the angular velocities are approximately equal. The initial values of the two moving satellites in the relative system are shown in Table 3.   Table 3. Initial values of the two satellite reconstructions in the relative coordinate system.

Z/(m/s)
The genetic algorithm population size is set to 100, and the maximum number of generations is set to 100. The simulation of the optimal fault reconstruction strategy calculated using the algorithm in Chapter 3 is shown in Figure 7.  Table 3. Initial values of the two satellite reconstructions in the relative coordinate system.

Module SF
.9615 The genetic algorithm population size is set to 100, and the maximum num generations is set to 100. The simulation of the optimal fault reconstruction strat culated using the algorithm in Chapter 3 is shown in Figure 7. From the scatter plot in Figure 6, the overall fitness value fluctuates and de The fluctuation is due to the inconsistency of the direction of the optimization ind in the iterative process, resulting in local optima. However, as the number of gene increases, the global optimum is gradually located. Therefore, after 100 iteratio From the scatter plot in Figure 6, the overall fitness value fluctuates and decreases. The fluctuation is due to the inconsistency of the direction of the optimization indicators in the iterative process, resulting in local optima. However, as the number of generations increases, the global optimum is gradually located. Therefore, after 100 iterations, the curve shows convergence. The histogram indicates that the optimal reconstruction scheme has been obtained through the genetic algorithm: X = [22,59,6,23,43,29], i.e., the healthy satellites numbered 22, 59, 6, 23, 43, and 29 correspond to the reconstructed faulty satellites 10, 18, 33, 40, 44, and 53. This reconfiguration scheme optimizes the health state, energy consumption, and fuel consumption balance of the entire satellite swarm.

Monte Carlo Simulation Verification
To prove the effectiveness and reliability of the reconstruction decision-making algorithm for satellite swarms with different degrees of faults, this section uses the Monte Carlo method. The steps are as follows: 1.
Set the number of faulty satellites randomly distributed in the range [1,30], i.e.
Randomly select N F satellites in a satellite swarm with N satellites and generate fault type S F . 3.
Set the health state indicator of the objective function in Section 3.2 and the rest of the indicators as constraints, solving the optimal reconfiguration strategy. 4.
Calculate the optimal health state and average health state of the satellite swarm corresponding to the reconstruction strategy. 5.
Repeat the above steps 300 times.
After collecting health state data 300 times, the health state data are sorted according to the severity of the fault from low to high, producing the simulation results in Figure 8. The ratio of faulty satellites to the total number of satellites (N F /N) is calculated, and the corresponding Monte Carlo optimal health state value is selected, yielding the simulation results in Figure 9.

Monte Carlo Simulation Verification
To prove the effectiveness and reliability of the reconstruction decision-making algorithm for satellite swarms with different degrees of faults, this section uses the Monte Carlo method. The steps are as follows: 1. Set the number of faulty satellites randomly distributed in the range [1,30] 3. Set the health state indicator of the objective function in Section 3.2 and the rest of the indicators as constraints, solving the optimal reconfiguration strategy. 4. Calculate the optimal health state and average health state of the satellite swarm corresponding to the reconstruction strategy. 5. Repeat the above steps 300 times.
After collecting health state data 300 times, the health state data are sorted according to the severity of the fault from low to high, producing the simulation results in Figure 8.
The ratio of faulty satellites to the total number of satellites ( / F N N ) is calculated, and the corresponding Monte Carlo optimal health state value is selected, yielding the simulation results in Figure 9. From the simulation results in Figure 7: 1. For a satellite swarm with 61 satellites, when the number of faulty satellites is in the range [1,30], the reconstruction decision algorithm restores the satellite swarm health state to 0.95~0.4. The optimal health is very close to the average health, and the dif- From the simulation results in Figure 7: 1.
For a satellite swarm with 61 satellites, when the number of faulty satellites is in the range [1,30], the reconstruction decision algorithm restores the satellite swarm health state to 0.95~0.4. The optimal health is very close to the average health, and the difference is within 10 −3 , which shows that the algorithm has high precision.

2.
The health state curve is a typical ponytail curve, and nearly 1/3 of the data show that the health state value of the restored satellite swarm is below 0.7; when the health state is less than 0.7, the health state begins to decrease rapidly and fluctuate. Figure 8 is obtained from Figure 7 by using the fault proportion as the abscissa, which makes the Monte Carlo simulation results more intuitive. It is obvious that the higher the proportion of faulty satellites, the lower the health of the satellite swarm and the more severe the health state jitter, i.e., the health state after fault reconstruction is very uncertain. At the same time, when the health state is 0.7, the proportion of faulty stars is 38%. That is, for the mission satellite swarm when the number of faulty satellites exceeds 38%, the health state of the entire swarm becomes unstable. Therefore, it is possible to obtain the fault threshold of the satellite swarm, which can be used to manage it.

Conclusions
In this paper, the problem of reconstructing a satellite swarm after faults occur is analyzed, and a reconstruction optimization algorithm based on hybrid multi-objective optimization is proposed. Taking the health state of the satellite swarm, the total energy consumption of reconstruction, and the balance of the fuel consumption as the optimization indicators, the optimal strategy of the fault reconstruction is solved using mixed multiobjective weighting and a genetic algorithm; then, simulation verification is carried out on a typical case. Finally, the effectiveness of the reconstruction algorithm is verified by Monte Carlo simulation. The simulation results show that as the number of population iterations increases, the fitness value will gradually approach the global optimum, and the curve shows convergence. Therefore, the reconstruction algorithm proposed in this paper provides the optimal reconstruction scheme for satellite swarms with different fault degrees. At the same time, the Monte Carlo simulation shows that the iteration of the algorithm has high precision, and the change in the health state of the satellite swarm conforms to the ponytail curve of the faulty spacecraft. When the number of faulty satellites exceeds 38%, the health state of the entire swarm becomes unstable. Therefore, a fault threshold can be given for faulty satellite swarms, which can be used to effectively manage satellite swarms.
In the future, more reconstruction indicators will be introduced into the system to improve the overall efficiency of the satellite swarm. In addition, the efficiency of the multi-objective reconstruction planning algorithm will be improved, thereby increasing the practicability of on-orbit calculations.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the study's design; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.