Low-orbit satellite deorbit control method and system based on particle swarm algorithm

The disclosure describes a low-orbit satellite deorbit control method and system based on a particle swarm algorithm. The method includes: constructing an objective function according to a position parameter and a perturbation acceleration of each of particles in a particle swarm as well as a preset Lagrangian multiplier and a penalty factor; determining an objective fitness of each particle and a population fitness of the particle swarm based on the objective function, updating the position parameter and the velocity of each particle, and obtaining a population optimal fitness at a maximum number of iterations; updating the Lagrangian multiplier and the penalty factor according to the population optimal fitness, comparing the objective deviation value with a preset convergence condition, and determining an objective population fitness and an objective position parameter according to a comparison result to obtain an objective deorbit mode.

CROSS REFERENCE TO RELATED APPLICATIONS

The present application claims the benefit of Chinese Patent Application No. 202010629155.7 filed on Jul. 1, 2020, the contents of which are incorporated herein by reference in their entirety.

FIELD OF THE INVENTION

The application relates to the technical field of constellation operation management, and in particular to a low-orbit satellite deorbit control method and system based on a particle swarm algorithm.

BACKGROUND OF THE INVENTION

With continuous increase of manned aerospace launching missions, low-orbital and medium-high orbital space is gradually congested, result in a progressive increase of the collision probability of satellite in-orbit operation. According to the “IADC Space Debris Mitigation Guidelines” compiled by Inter-Agency Space Debris Coordination Committee (IADC), when a natural orbit will serve more than 25 years, a low-orbit satellite at end of its lifetime needs to enter a disposal orbit and then burn down in the atmosphere within 25 years. Moreover, development of small satellite technology and reduction of launching costs may cause a sharp increase in the number of satellites in a low-orbit satellite constellation. OneWeb plans to launch a constellation composed of 2,620 satellites, Samgsung plans to launch a constellation composed of 4,600 satellites, Boeing plans to launch a constellation composed of 2956 satellites and Starlink eventually plans to launch a constellation composed of 42,000 satellites. Such constellations are huge and their satellites will affect the low-orbit space environment after expiration of their service life. Therefore, it is particularly important to take a reasonable approach to control the satellite at end of its lifetime in the low-orbit satellite constellation.

In related technologies, for low earth orbit (LEO) satellites, deorbit methods are generally divided into electric power cable deorbit, atmospheric drag deorbit, and satellite's own thruster deorbit. The satellite's own thruster deorbit can exhibit characteristics of short deorbit duration and flexible control when compared to other deorbit methods. As a typical small thruster, electric thrusters have been widely used in deep space exploration, orbit maintenance and orbit maneuvering. All satellites in the Starlink constellation may be equipped with electric thrusters. Fromm et al. simulated the deorbit duration of a satellite under the action of different small thrusters, and obtained a thrust range of the small thruster suitable for the deorbit of the simulated objective satellite. Huang et al. simulated and analyzed a deorbit process of a single satellite and OneWeb constellation, and obtained an optimal deorbit strategy of the single satellite and the constellation. However, the research on deorbit control in related technologies mainly focuses on the derivation of control rate, and has rarely mentioned algorithm implementation. Particle swarm algorithm has achieved excellent effects in various continuous space optimization problems, neural network training and other fields, but little on deorbit control, which may limit further research and accurate control of satellite space movement to a certain extent.

SUMMARY OF THE INVENTION

The problem to be solved by embodiments of the present invention is to provide a low-orbit satellite deorbit control method and system based on a particle swarm algorithm.

In order to solve above technical problem, the invention may be implemented as follows.

According to one aspect of the disclosure, it is provided a low-orbit satellite deorbit control method based on a particle swarm algorithm, including:

at step S1, constructing an objective function according to a position parameter and a perturbation acceleration of each of particles in a particle swarm as well as a preset Lagrangian multiplier and a preset penalty factor;

at step S2, calculating an objective function value of each particle in the particle swarm based on the objective function, determining an objective fitness of each particle and a population fitness of the particle swarm according to a calculation result, updating the position parameter and the velocity of each particle, and obtaining a population optimal fitness at a maximum number of iterations and the position parameters of particles corresponding to the population optimal fitness; and

at step S3, updating the Lagrangian multiplier and the penalty factor according to the population optimal fitness, obtaining a current number of iterations and an objective deviation value based on the population optimal fitness, comparing the current number of iterations and the objective deviation value with a preset iteration condition and a preset convergence condition, and determining an objective population fitness and an objective position parameter according to a comparison result to obtain an objective deorbit mode.

In an exemplary embodiment of the disclosure, the step of constructing the objective function according to the position parameter and the perturbation acceleration of each of particles in the particle swarm as well as the preset Lagrangian multiplier and the preset penalty factor further includes:

at step S11, obtaining a rate of change of each orbital element over time under the thrust applied by a thruster, and determining a control rate containing current costate variables according to the rate of change;

at step S12, calculating a thrust vector, a fuel consumption rate, and a perturbation acceleration of the thrust applied by the thruster;

at step S13, determining a particle position vector according to the current costate variables, and using the particle position vector as the position parameter of each of the particles in the particle swarm;

at step S14, performing orbit integration of the control rate based on the position parameter and the perturbation acceleration, and constructing the objective function according to a integration result, the preset Lagrangian multiplier and the preset penalty factor.

In an exemplary embodiment of the disclosure, the steps of comparing the current number of iterations and the objective deviation value with the preset iteration condition and the preset convergence condition respectively and determining the objective population fitness and the objective position parameter according to the comparison result to obtain the objective deorbit mode further includes:

if the current number of iterations meets the preset iteration condition and the objective deviation value meets the preset convergence condition, determining the population optimal fitness as the objective population fitness, and determining the position parameter corresponding to the objective population fitness as the objective position parameter;

if the current number of iterations does not meet the preset iteration condition or the objective deviation value does not meet the preset convergence condition, updating the objective function based on the updated Lagrangian multiplier and the updated penalty factor, and continuing to perform step S14and subsequent steps, until the current number of iterations meets the preset iteration condition and the objective deviation value meets the preset convergence condition.

In an exemplary embodiment of the disclosure, in step S11, the step of determining the control rate containing the current costate variables according to the rate of change further includes:

based on the rate of change and the current costate variables, obtaining a Hamiltonian function by equation (7):

where H is a Hamiltonian function value, λa, λeand λiare the current costate variables, a, e and i are semi-major axis, eccentricity and orbital inclination of the operational orbit respectively, and t is time;

deriving a preset angle corresponding to the thrust applied by the thruster according to the Hamilton function (7) to obtain the control rate containing the current costate variables.

In an exemplary embodiment of the disclosure, in step S12, the thrust vector is obtained by equation (11), the fuel consumption rate is obtained by equation (12), and the perturbation acceleration is obtained by equation (14);

where η is small thruster efficiency, P is small thruster power, m is spacecraft mass, g is acceleration of gravity, Ispis small thruster specific impulse, GMeis gravitational coefficient of the earth, having a value of 3.986×1014, Reis radius of the earth, having a size of 6378137 m,

KD=CD⁢Am⁢kr,CD,Am
is the satellite drag coefficient and area-to-mass ratio,

kr=1-2⁢σ⁢h⁢cos⁢iv2,
α=7.292×10−5is rotation velocity of the earth, ρ is atmospheric density at the altitude of the orbit on which the satellite is located, v is tangential velocity, r is radial distance, θ is true anomaly, ω is argument of perigee, i is orbital inclination.

In an exemplary embodiment of the disclosure, in step S14, the objective function J shown in equation (16) is constructed:
J=f(xik)+Σj=6λjhi(xik)+Σj=16rj[hj(xik)]2(16);

where hj(xik)=|σj−σj*|, j=1, 2, . . . , 6 (17); hj(xik) is an absolute value of the difference between the current orbital elements and the objective orbital elements generated by the control rate of a particle i in the k-th iteration, λj, rjare the Lagrangian multiplier and the penalty factor respectively, f(xik) is deorbit duration generated by the control rate of the particle i in the k-th iteration, which is obtained by orbital integration of the control rate based on the position parameter and the perturbation acceleration.

In an exemplary embodiment of the disclosure, in the process of orbit integration, if the orbit altitude increases and/or the eccentricity deviates from a threshold, the obtained integration result is deleted.

In an exemplary embodiment of the disclosure, in step S2, the position parameter and the velocity of each particle are updated by equation (18):

wherein the position and the velocity of a certain particle in the k-th iteration are xikand vikrespectively, ω is inertia weight, pikis the best historical position of an individual in the past k iterations, pgkis the best population historical position of the particle swarm after k iterations, m1is a random number between [0,1], ρkis a random direction parameter, ns, nfrepresents the relationship between current best population historical position pgkand that of the previous iteration pgk-1, and if both are the same, nfis accumulated by one, otherwise nsis accumulated by one, s, f represent thresholds of ns, nfrespectively, ωmaxis maximum value of inertia weight, ωminis minimum value of inertia weight, kmaxis maximum number of iterations;

wherein in each previous iteration, an individual particle is compared with an individual historical value, and the particle swarm is compared with a population historical value, so that respective particle positions pikand pgkcorresponding to the minimum fitness are obtained.

In an exemplary embodiment of the disclosure, in step S14, the Lagrangian multiplier and the penalty factor are updated by equation (20):

wherein the value of xbestnis equal to the population optimal fitness pgkobtained in step S2, εfis a constraint threshold, λjn, rjnare the Lagrangian multiplier and the penalty factor respectively, where n is sequence number, hj(xik) is an absolute value of the difference between the current orbital elements and the objective orbital elements generated by the control rate of a particle i in the k-th iteration.

In an exemplary embodiment of the disclosure, the orbital elements include semi-major axis a, eccentricity e, orbital inclination i, right ascension of the ascending node Ω, argument of perigee ω and true anomaly θ.

In an exemplary embodiment of the disclosure, the objective fitness and the population fitness are fitnesses corresponding to the minimum value of the objective function.

According to another aspect of the disclosure, it is provided a low-orbit satellite deorbit control system based on a particle swarm algorithm, including:

a construction unit configured to construct an objective function according to a position parameter and a perturbation acceleration of each of particles in a particle swarm as well as a preset Lagrangian multiplier and a preset penalty factor;

an updating unit configured to calculate an objective function value of each particle in the particle swarm based on the objective function, determine an objective fitness of each particle and a population fitness of the particle swarm according to a calculation result, update the position parameter and the velocity of each particle, and obtain a population optimal fitness at a maximum number of iterations and the position parameters of particles corresponding to the population optimal fitness; and

a processing unit configured to update the Lagrangian multiplier and the penalty factor according to the population optimal fitness, obtain a current number of iterations and an objective deviation value based on the population optimal fitness, compare the current number of iterations and the objective deviation value with a preset iteration condition and a preset convergence condition, and determine an objective population fitness and an objective position parameter according to a comparison result to obtain an objective deorbit mode.

According to yet another aspect of the disclosure, it is provided a readable storage medium storing one or more readable programs adapted to cause a processor is executable by one or more processors to implement the steps of a low-orbit satellite deorbit control method based on a particle swarm algorithm according to the first aspect of the disclosure.

According to still yet another aspect of the disclosure, it is provided a computer device comprising a memory configured to store computer programs thereon, and a processor configured to execute the computer programs stored on the memory to implement the steps of a low-orbit satellite deorbit control method based on a particle swarm algorithm the first aspect of the disclosure.

Advantageous effects obtained by embodiments of the invention disclosed herein will be briefly explained.

(1) Compared with the traditional shooting method, the low-orbit satellite deorbit control method based on the particle swarm algorithm provided by the disclosure may be less sensitive to the initial value, and the search can be performed upon a given search range, which will greatly reduce difficulties in initial value determination.

(2) The low-orbit satellite deorbit control method based on the particle swarm algorithm provided by the disclosure has been improved for the position, velocity iteration formula and inertia weight iteration formula of the particle swarm algorithm, reducing the possibility of the algorithm falling into a local optimal solution, and improving the overall search ability and global convergence of the algorithm.

(3) The low-orbit satellite deorbit control method based on the particle swarm algorithm provided by the disclosure adopts the augmented Lagrangian method to introduce a secondary penalty term, and then corrects the Lagrangian multiplier by iteration to avoid ill-conditioned problems due to excessively large parameter values by the penalty function method, as well as transforms some constraints into the fitness function through the orbit integration program, which may reduce the constraints of the objective function and improve the calculation speed of the program. The present method can be utilized to obtain an optimal control rate that minimizes the low-thrust deorbit duration of low-orbit satellites, which introduces the advantages of insensitivity to initial values and fast convergence speed, and provides an effective solution for the engineering application of low-orbit satellite deorbit with small thrust.

DETAILED DESCRIPTION OF THE EMBODIMENTS

The technical problems to be solved and technical solutions and beneficial effects of the present disclosure will become apparent from the following detailed description with reference to the accompanying drawings and embodiments. It should be noted that the exemplary embodiments described herein are merely intended to explain the present invention, rather than limit the scope of the present invention.

The exemplary embodiment of the disclosure provides a low-orbit satellite deorbit control method based on a particle swarm algorithm. Firstly, three dimensions of particles are used to represent three costate variables, and all particles in the population are randomly initialized within a feasible region and a Lagrangian multiplier and a penalty factor are initialized respectively. Subsequently, an objective function is used as a fitness function of the particle swarm algorithm, a fitness function value is calculated for each particle, an optimal fitness of each particle and a population optimal fitness are updated, and the position and velocity of each of the particles are updated until a maximum number of iterations is reached; Finally, augmented Lagrangian function parameters are update, a new objective function is generate, the previous step is repeat until the number of iterations and constraint threshold requirements are met, and an optimization result is updated. This process combines the Hamiltonian function and the augmented Lagrangian particle swarm algorithm. Under the augmented Lagrangian function, the particle swarm algorithm can obtain costate variables corresponding to the shortest deorbit duration after multiple iterations, thereby achieving the optimal control rate.

Illustratively,FIG.1shows a flowchart of a low-orbit satellite deorbit control method based on a particle swarm algorithm according to an exemplary embodiment of the disclosure; andFIG.2shows an outer block diagram of an augmented Lagrangian particle swarm algorithm according to an exemplary embodiment of the disclosure. The method of embodiments of the disclosure will be described below with reference toFIGS.1and2.

At S1, an objective function is constructed according to a position parameter and a perturbation acceleration of each of particles in a particle swarm as well as a preset Lagrangian multiplier and a penalty factor.

The method further includes the following steps.

At step S11, a rate of change of each orbital element over time under the thrust applied by a thruster is obtained, and a control rate containing current costate variables is determined according to the rate of change.

In this embodiment, the orbital elements include semi-major axis a, eccentricity e, orbital inclination i, right ascension of the ascending node Ω, argument of perigee ω and true anomaly θ, and the control rate containing current costate variables is obtained by the following steps S111, S112and S113.

At step S111, the rate of change of each orbital element over time under the thrust applied by the thruster is obtained. In this case, the influence of the thrust by the thruster on the orbital element is calculated by:

At S112, based on the rate of change and the current costate variables, a Hamiltonian function is obtained by equation (7).

In this embodiment, the Hamiltonian function is calculated by:

where H is a Hamiltonian function value, λa, λeand λiare the current costate variables, a, e and i are semi-major axis, eccentricity and orbital inclination of the operational orbit respectively, and t is time.

At step S113, a preset angle corresponding to the thrust applied by the thruster is derived according to the Hamilton function (7) to obtain the control rate containing the current costate variables.

Herein,FIG.4shows a schematic diagram of the thrust direction of a small thruster according to an exemplary embodiment of the disclosure. The process of determining the control rate in the embodiment will be described below with reference toFIG.4. During the process, a control rate containing current costate variables at the angles α and β as shown inFIG.4is calculated, where the partial derivatives of the Hamilton function may be taken with respect to α and β respectively, and a plurality of simultaneous equations may be set as follows:

At step S12, a thrust vector, a fuel consumption rate, and a perturbation acceleration of the thrust applied by the thruster are calculated.

In this embodiment, the vector of the small thrust and the fuel consumption rate shown inFIG.4are calculated, and the non-spherical perturbation and the perturbation acceleration of the atmospheric drag perturbation are calculated. J2

At step S121, the thrust vector of the thruster is calculated by:

The fuel consumption rate is calculated by:

where η is small thruster efficiency, P is small thruster power, m is spacecraft mass, g is acceleration of gravity, and Ispis small thruster specific impulse.

At step S122, the nonspherical J2perturbation and the perturbation acceleration of the atmospheric drag perturbation are calculated.

The perturbation acceleration is calculated by:

where, GMeis gravitational coefficient of the earth, having a value of 3.986×1014, Reis radius of the earth, having a size of 6378137 m,

KD=CD⁢Am⁢kr,CD,Am
is the satellite drag coefficient and area-to-mass ratio,

kr=1-2⁢σ⁢h⁢⁢⁢cos⁢⁢iv2,
σ=7.292×10−5is rotation velocity of the earth, ρ is atmospheric density at the altitude of the orbit on which the satellite is located, v is tangential velocity, r is radial distance, θ is true anomaly, ω is argument of perigee, and i is orbital inclination.

At step S13, a particle position vector is determined according to the current costate variables, and the particle position vector is used as the position parameter of each of the particles in the particle swarm.

In this embodiment, the costate variable is used as the particle position vector, and the particle position vector is calculated by:
λi=[λaλeλi]  (15).

At step S14, orbit integration of the control rate is performed based on the position parameter and the perturbation acceleration, and the objective function is constructed according to a integration result, the preset Lagrangian multiplier and the penalty factor;

In this embodiment, the objective function J shown in equation (16) is constructed:
J=f(xik)+Σj=16λjhj(xik)+Σj=16rj[hj(xik)]2(16);

where hj(xik)=|σj−σj*|, j=1, 2, . . . , 6 (17); hj(xik) is an absolute value of the difference between the current orbital elements and the objective orbital elements generated by the control rate of a particle i in the k-th iteration, λj,rjare the Lagrangian multiplier and the penalty factor respectively, f(xik) is deorbit duration generated by the control rate of the particle i in the k-th iteration, which is obtained by orbital integration of the control rate based on the position parameter and the perturbation acceleration (implemented by integration of equations (1) to (6)). λj,rjremains unchanged during the sub-problem solving operation.

It should be noted that some constraints may be implemented in the orbital integration process through corresponding processing. Since the initial position value of each of the particles is random, some particles may cause abnormal conditions such as an increase in orbital altitude and abnormal eccentricity during orbital integration. To solve such a problem, an interruption condition can be set during the orbital integration process and the deorbit duration is returned as an infinite value. The fitness value of each of these particles also becomes infinite, which is eliminated in fitness screening.

At step S2, an objective function value of each particle in the particle swarm is calculated based on the objective function, an objective fitness of each particle and a population fitness of the particle swarm are determined according to a calculation result, the position parameter and the velocity of each particle are updated, and a population optimal fitness at a maximum number of iterations and the position parameters of particles corresponding to the population optimal fitness is obtained, wherein the population optimal fitness is a population fitness that is output when the preset maximum number of iterations is met.

In this embodiment, the objective fitness and the population fitness are fitnesses corresponding to the minimum value of the objective function.

The optimal fitness of each particle pikand the population optimal fitness pgkare obtained by comparing respective objective function values of respective particles in the population (the smallest is the optimum), the position and the velocity of each particle are updated, and an optimization result of the subroutine is output when the preset number of iterations is met (seeFIG.3).

In this case, the position parameter and the velocity of each particle are updated by equation (18):

wherein the position and velocity of a certain particle in the k-th iteration are xikand vikrespectively, ω is inertia weight, pikis the best historical position of the individual in the past k iterations, pgkis the best population historical position of the particle swarm after k iterations, m1is a random number between [0,1], ρkis a random direction parameter, ns,nfrepresents the relationship between current best population historical position pgkand that of the previous iteration pgk-1, and if both are the same, nfis accumulated by one, otherwise nsis accumulated by one, s, f represent thresholds of ns,nfrespectively, ωmaxis maximum value of inertia weight, ωminis minimum value of inertia weight, kmaxis maximum number of iterations;

wherein in each previous iteration, an individual particle is compared with an individual historical value, and the particle swarm is compared with a population historical value, so that respective particle positions pikand pgkcorresponding to the minimum fitness are obtained.

At step S3, the Lagrangian multiplier and the penalty factor are updated according to the population optimal fitness, a current number of iterations and an objective deviation value are obtained based on the population optimal fitness, the current number of iterations and the objective deviation value are compared with a preset iteration condition and a preset convergence condition, and an objective population fitness and an objective position parameter are determined according to a comparison result to obtain an objective deorbit mode. In this case, the objective deviation value based on the population optimal fitness is the deviation value between the semi-major axis or the perigee generated by the position parameter integration corresponding to the population optimal fitness obtained in step S2and the target.

In this embodiment, if the current number of iterations meets the preset iteration condition and the objective deviation value meets the convergence condition, the population optimal fitness is determined as the objective population fitness, and the position parameter corresponding to the objective population fitness is determined as the objective position parameter.

If the current number of iterations does not meet the preset iteration condition or the objective deviation value does not meet the convergence condition, the objective function is updated based on the updated Lagrangian multiplier and the penalty factor, and step S14and subsequent steps are continued to be performed, until the current number of iterations meets the preset iteration condition and the objective deviation value meets the convergence condition.

In this case, the Lagrange multiplier λjnand the penalty factor rjnis updated by calculation of:

In this case, if the optimal solution xbestndoes not satisfy the constraint condition, the procedure returns to step S14, a new objective function is calculated by the updated Lagrangian multiplier and penalty factor, and then the procedure proceeds to subsequent steps until the constraint condition is satisfied εf.

In the low-orbit satellite deorbit control method based on the particle swarm algorithm of the disclosure, firstly a control rate with costate variables is derived from a Hamiltonian function; subsequently, costate variable vectors are used as position parameters of particles in a particle swarm to operate a particle swarm algorithm subroutine with a current objective function (FIG.3), an optimal fitness pikof each particle, a population optimal fitness pgkand the position and velocity of each of the particles are updated in each iteration, and an optimization result is output when a maximum number of iterations is met; finally, according to the optimization result obtained by the subroutine, parameters in the augmented Lagrangian function are updated, and then a new objective function is obtained for the next subroutine iteration, and an optimization result is output when the number of iterations and constraint threshold requirements are met. The method can be utilized to obtain an optimal control rate that minimizes the low-thrust deorbit duration of low-orbit satellites, which introduces the advantages of insensitivity to initial values and fast convergence speed, and provides an effective solution for the engineering application of low-orbit satellite deorbit with small thrust.

The low-orbit satellite deorbit control method based on the particle swarm algorithm of the disclosure will be described below in conjunction with this embodiment. The embodiment may be achieved by the following steps: at the first step, parameters of a spacecraft itself and parameters of a small thruster are given; at the second step, semi-major axis, eccentricity and orbital inclination of an operational orbit and two disposal orbits of the spacecraft are given; at the third step, parameters of an augmented Lagrangian particle swarm algorithm are given; at the fourth step, deorbit duration without thrust is calculate; at the fifth step, an optimal control rate, deorbit duration and fuel consumption of the two disposal orbits are calculated.

Herein the parameters of the spacecraft itself and the parameters of the small thruster are shown in Table 1, the semi-major axis, eccentricity and orbit inclination on an operational orbit and two disposal orbits of the spacecraft are shown in Table 2, the parameters of the augmented Lagrangian particle swarm algorithm are shown in Table 3, and the optimal control rate, deorbit duration and fuel consumption of the two disposal orbits are shown in Table 4.

TABLE 3Augmented Lagrangian Particle Swarm Algorithm ParametersAlgorithm ParameterNumerical ValueNumber of Populations N500Maximum Number of Iterations kmax50Maximum Inertia Weight ωmax0.9Minimum Inertia Weight ωmin0.4Threshold of nss15Threshold of nff5Initial Value of Random Direction Parameter ρ01Initial Penalty Factor r01000Constraint Threshold εf0.0001

FIG.5is a graph showing variation of semi-major axis over time without thrust according to Embodiment 1 of the disclosure. As seen fromFIG.5, the deorbit duration without thrust under natural conditions will be 140 years, far beyond the 25 years deorbit requirement.FIGS.6and7show the results of iterations in the two disposal orbits respectively.FIG.6illustrates a total of 60 iterations in the process, and the semi-major axis deviation convergence condition will be met at the 57th iteration.FIG.7illustrates a total of 25 iterations in the process, and the perigee deviation convergence condition will be met at the 12th iteration.FIGS.8and9show graphs of variation of semi-major axis, eccentricity and orbital inclination over time on the two disposal orbit respectively. As seen fromFIGS.8and9, the spacecraft can reach the corresponding objective orbit (disposal orbit) under the optimal control rate obtained by the disclosure. Based on this effect, the effectiveness of the optimal control rate obtained by the disclosed algorithm is proved.

Another exemplary embodiment of the disclosure provides a low-orbit satellite deorbit control system based on a particle swarm algorithm, including: a construction unit configured to construct an objective function according to a position parameter and a perturbation acceleration of each of particles in a particle swarm as well as a preset Lagrangian multiplier and a penalty factor; an updating unit configured to calculate an objective function value of each particle in the particle swarm based on the objective function, determine an objective fitness of each particle and a population fitness of the particle swarm according to a calculation result, update the position parameter and the velocity of each particle, and obtain a population optimal fitness at a maximum number of iterations and the position parameters of particles corresponding to the population optimal fitness; and a processing unit configured to update the Lagrangian multiplier and the penalty factor according to the population optimal fitness, obtain a current number of iterations and an objective deviation value based on the population optimal fitness, compare the current number of iterations and the objective deviation value with a preset iteration condition and a preset convergence condition, and determine an objective population fitness and an objective position parameter according to a comparison result to obtain an objective deorbit mode.

The detailed description of each unit in the aforementioned control system can refer to the description of each step in the control method, which will not be repeatedly described. The aforementioned control system can achieve the same function as the control method.

On the basis of the above-mentioned embodiments, an embodiment of the disclosure provides a computer readable storage medium storing one or more readable programs adapted to cause a processor is executable by one or more processors to implement the steps of a low-orbit satellite deorbit control method based on a particle swarm algorithm as described in the first embodiment, for example.

Optionally, the aforementioned storage medium may be a non-transitory computer-readable storage medium. For example, the non-transitory computer-readable storage medium may be ROM, random access memory (RAM), CD-ROM, magnetic tape, floppy disk, and optical data storage device, etc.

Additionally, another embodiment of the disclosure further provides a computer device, including a processor and a memory.

The memory is configured to store computer programs therein.

The processor is configured to execute the computer programs stored in the memory to implement a low-orbit satellite deorbit control method based on a particle swarm algorithm as described in the first embodiment, for example.

The aforementioned processor may be a general-purpose processor, including central processing unit (CPU), network processor (NP), etc.; it may also be digital signal processor (DSP), application specific integrated circuit (ASIC), field-programmable gate array (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components.

Accordingly, the above description is only considered as one embodiment of the disclosure, rather than limitation of the invention in any way. Although the invention has been described in the above with the preferred embodiments, the invention is not limited to those. Those skilled in the art, without departing from the scope of the invention, can make some equivalent changes or modifications using the content disclosed above, which will fall within the scope of the invention.