A simulation method for the time-saving orbit change mission of "one-to-many" approach to space targets

Abstract

The invention relates to a method for simulating a time-saving orbit change task for "one-to-many" approach of space targets. A service aircraft adopts a double-pulse orbit change strategy to sequentially approach each target satellite, and determines an objective function for the time-saving orbit change; an optimization level is set, wherein, in the first level, a total orbit change time optimization variable is set; in the second level, an approach sequence optimization variable is set, and an approach sequence of the service aircraft arriving at the target satellite is determined by the variable size; in the third and fourth levels, within the total time limit range for the service aircraft to complete the approach of all target satellites, a value range of an approach orbit change end time optimization variable for each target satellite that satisfies both a measurement and control window and an illumination window is determined; and a value range of an approach orbit change start time optimization variable for each target satellite is determined; in the fifth level, under a given approach sequence, given approach orbit change end time and approach orbit change start time are optimized according to the objective function to obtain an optimization result; and in accordance with the set optimization level order, optimization is started by level, and the optimal result is selected from the optimization results to obtain the time-saving orbit change strategy.

Description

Technical Field

[0001] The invention belongs to the technical field of on-orbit service and maintenance of spacecraft. Background Art

[0002] On-orbit service and maintenance is one of the hot spots at the forefront of aerospace technology. After the service vehicle approaches the space target satellite, it can carry out close observation, module replacement, fuel filling, auxiliary orbit change and other service operations, which can realize the on-orbit fault recovery and extend the life of the target satellite, so as to continue to play the value of on-orbit application and generate great economic benefits. In order to carry out service and maintenance operations on the space target satellite, the prerequisite is to approach the target satellite and dock at zero distance.

[0003] The service vehicle itself has a cost, and operations such as approaching, docking, and repairing target satellites also require a certain amount of cost. A service vehicle approaching and repairing multiple target satellites is a necessary way to achieve the economic benefits of in-orbit servicing. Therefore, the research on the simulation method of the most time-saving orbit change task for the "one-to-many" approach of a service vehicle to multiple target satellites has very important practical significance.

[0004] At present, there are no public literature reports or patents at home and abroad on the same spatial "one-to-many" approach and the most time-saving trajectory change task simulation method as this patent. Summary of the invention

[0005] The technical problem solved by the present invention is to provide a method for simulating orbit change tasks for "one-to-many" approach to space targets, which is suitable for optimizing the orbit change control strategy of a high-orbit service aircraft to approach and perform maintenance tasks on multiple space targets in the same high orbit within a given time during an on-orbit service mission, so as to achieve the purpose of saving the approach process time as much as possible.

[0006] The technical solution of the present invention is: a method for simulating a time-saving orbit change task for a "one-to-many" approach to a space target, comprising:

[0007] The service spacecraft uses a double-pulse orbit change strategy to approach each target satellite in sequence and determine the objective function of the orbit change with the least time.

[0008] Set the optimization level, where

[0009] The first layer sets the total orbit change time optimization variable, which is the time T that the service spacecraft takes to complete the orbit change mission of approaching all n target satellites.

[0010] The second layer: setting the approach sequence optimization variables, and determining the approach sequence of the service spacecraft to the target satellite by the variable size;

[0011] The third and fourth layers: within the total time limit for the service spacecraft to complete the approach to all target satellites, determine the value range of the optimization variable of the approach change end time for each target satellite that satisfies both the measurement and control window and the illumination window; determine the value range of the optimization variable of the approach change start time for each target satellite;

[0012] The fifth layer: Under a given approach sequence, the given approach change end time and approach change start time are optimized according to the above objective function to obtain the optimization result;

[0013] According to the set optimization level order, start optimization by level, and select the best result from the optimization results to obtain the most time-saving track change strategy.

[0014] Preferably, the objective function is to ensure that the fuel consumed by the service aircraft to complete all target approach missions is less than a given fuel constraint and the absolute value of the difference is minimal.

[0015] Preferably, the objective function is:

[0016]

[0017] Where Δv u The velocity increment consumed by the service vehicle to complete the mission of approaching and changing orbits to n target satellites, Δv g is a given consumable speed increment constraint; σ takes a large positive number, at least on the order of 1e8.

[0018] Preferably, the approach sequence is set by:

[0019] Number the n target satellites as 1, 2, ..., n, and set the approach sequence S 1 , S 2 ,…,S n The value range is [0,1], and the service aircraft is based on S 1 , S 2 ,…,S n The target satellites with corresponding numbers are approached from first to last in the order from small to large.

[0020] Preferably, the value range of the optimization variable of the approaching orbit change end time for each target satellite that satisfies both the measurement and control window and the illumination window is determined by the following method:

[0021] Within the total time limit for the service spacecraft to complete the approach to all target satellites, the tracking and control window interval of the target satellite is calculated based on the position of the ground tracking and control station, the trackable azimuth and pitch angle, and the orbital parameters of the target satellite;

[0022] The illumination window interval of the target satellite is calculated according to the orbital parameters of the target satellite;

[0023] The intersection of the measurement and control window interval and the illumination window interval is taken to obtain the window interval in which the service aircraft can approach the target satellite;

[0024] According to the order in which the target satellites are approached by the serviced spacecraft, the approach and orbit change end time optimization variable is set to a floating point number with a value greater than 1. The value range of the approach and orbit change end time optimization variable is determined according to the principle that the integer part represents which approach window the approach and orbit change end time belongs to the target satellite, and the decimal part represents the proportion of the approach and orbit change end time in the window.

[0025] Preferably, the value range of the optimization variable of the end time of approaching orbit change of the i-th target satellite is:

[0026]

[0027] Among them, the number of approachable windows of the i-th target satellite is n wi , n represents the total number of target satellites; b 1 The value range is [0.5,1), b 2 The value range is (1,1.5).

[0028] Preferably, the value range of the optimization variable for the approaching orbit change start time of each target satellite is [τ i ,λ i ];

[0029] τ i =λ i-1 +h i-1

[0030] In the formula, λ i-1 The time when the service vehicle approaches the i-1th target satellite and ends the orbit change; i h is the end time of the service spacecraft approaching the i-th target satellite; i-1 The operating time used by the service vehicle to provide maintenance services to the i-1th target.

[0031] Preferably, in the fifth layer of optimization iteration, an optimization algorithm is used to iteratively simulate the double-pulse velocity increment of the service spacecraft approaching the orbit change. After the simulation, the angle between the sunlight to the service spacecraft vector and the target satellite to the service spacecraft vector are calculated when the service spacecraft approaches. If the angle is less than the front-light constraint angle, it is determined that the front-light observation condition does not meet the requirements, and the sum of the double-pulse velocity increment moduli is set to a very large number, the order of magnitude of which is greater than the order of magnitude of the result after optimization by the optimization algorithm.

[0032] Preferably, optimization is started in the following manner by level:

[0033] The first layer of optimization variables obtain fixed values ​​within the value range, the second layer iterates and optimizes according to the fixed values ​​of the first layer, during the second layer optimization process, the third layer needs to iterate according to the fixed values ​​of the optimization variables of the second layer, during the third layer optimization process, the fourth layer needs to iterate according to the fixed values ​​of the optimization variables of the third layer, during the fourth layer optimization process, the fifth layer needs to iterate according to the fixed values ​​of the optimization variables of the fourth layer.

[0034] Preferably, the optimization sequence is as follows: when the fifth layer completes the iteration, the fourth layer changes the value of the optimization variable, and the fifth layer iterates again according to the value of the fourth layer, until the fourth layer completes the iteration of all values; when the fourth layer completes the iteration, the third layer changes the value of the optimization variable, and the fourth layer iterates again according to the value of the third layer, until the third layer completes the iteration of all values; when the third layer completes the iteration, the second layer changes the value of the optimization variable, and the third layer iterates again according to the value of the second layer, until the second layer completes the iteration of all values; when the second layer completes the iteration, the first layer changes the value of the optimization variable, and the second layer iterates again according to the value of the first layer, until the first layer completes the iteration of all values, at which point the optimization simulation is completed, and the time-saving and optimal orbit change strategy for the service vehicle approaching multiple target satellites is obtained.

[0035] The beneficial effects of the present invention compared with the prior art are:

[0036] (1) The present invention provides a trajectory change strategy for a high-orbit service vehicle to reach and rendezvous with multiple space targets in the same high orbit, and provides a simulation calculation method for the most time-saving trajectory change sequence for the service vehicle to sequentially approach n targets, which can achieve the most time-saving "one-to-many" arrival trajectory change rendezvous under a given fuel constraint;

[0037] (2) The present invention comprehensively considers the constraints of measurement and control, lighting, time, etc. during the process of the service aircraft changing orbits to multiple targets, and uses genetic, ant colony, particle swarm and other optimization algorithms to optimize the design of the orbital maneuvering strategy of the service aircraft changing orbits to multiple targets in succession. There are many optimization design variables, including the order of approach of the service aircraft to multiple targets, the start approach time, the end approach time, and the vector value of the orbit change speed increment. If the number of target satellites is set to n, the number of optimization variables reaches 1+9n. The increase in the number of variables leads to a large scale of optimization design and difficulty in convergence. The process of the service aircraft approaching multiple targets usually needs to be calculated within the set mission time, which is convenient for providing time range constraints for the design of the start approach time and the end approach time. In order to solve the approach change with the most time saving, the time range has changed from a design variable to a design target, which increases the difficulty of problem solving. To address this problem, the present invention proposes a five-level optimization strategy. The first level guesses the mission time of the "one-to-many" orbit change of the service aircraft, the second level optimizes the approach order, the third level optimizes the end approach time, the fourth level optimizes the start approach time, and the fifth level optimizes the vector value of the orbit change speed increment. The optimization goal is that the fuel consumed by the service aircraft to complete the approach mission to all targets is less than the given fuel constraint and the absolute value of the difference is the smallest. During the optimization search process, the objective function is simulated and calculated under the five-level optimization strategy, which effectively improves the convergence speed. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] Figure 1 Service vehicle begins orbital maneuvers

[0039] Figure 2 Service vehicle approaches target satellite 3

[0040] Figure 3 Service spacecraft approaches target satellite 4

[0041] Figure 4 Service vehicle approaches target satellite 2

[0042] Figure 5 Service vehicle approaches target satellite 1 DETAILED DESCRIPTION

[0043] The present invention will be further described below in conjunction with the embodiments.

[0044] The present invention relates to a method for simulating an orbit change mission for a "one-to-many" approach to a space target, and is suitable for optimizing the orbit change control strategy for an on-orbit service mission in which a high-orbit service vehicle approaches a plurality of space targets in the same high orbit within a given fuel constraint to perform a maintenance mission, so as to achieve the purpose of minimizing the approach process time.

[0045] A method for simulating a time-saving orbit change mission for a "one-to-many" approach to a space target, the implementation steps of which are as follows:

[0046] (1) Design objective function

[0047] The purpose of the present invention is to find a simulation method for the most time-saving orbit change task for a space "one-to-many" approach task. First, the objective function needs to be designed.

[0048] Fuel is a key constraint for space orbit change missions. Generally, the greater the orbit change fuel consumption, the shorter the orbit change time; the smaller the orbit change fuel consumption, the longer the orbit change time. For an orbit change mission with a given fuel consumption constraint, the orbit change mission can be completed as quickly as possible by making full use of the given fuel.

[0049] Therefore, the optimization design objective function for the time-saving trajectory change task of the "one-to-many" approach of space targets is to ensure that the fuel consumed by the service vehicle to complete the approach mission to all targets is less than the given fuel constraint and the absolute value of the difference is the smallest. In this way, the time-saving problem is converted into making the most effective use of the given amount of fuel, so that the problem can be solved under the condition of taking time as a variable.

[0050] The number of target satellites for the “one-to-many” approach mission is set to n. The service spacecraft uses a double-pulse orbit change strategy to approach each target satellite in sequence. The velocity increment vector modulus required to start the orbit change for approaching the i-th target satellite is Δv i1 , the velocity increment vector modulus required to end the track change is Δv i2 , then the speed increment consumed by the service spacecraft to complete the approach and orbit change mission to n target satellites is expressed as The consumable velocity increment constraint is given by Δv g , then the objective function of the “one-to-many” approach task is:

[0051]

[0052] σ is taken as a large positive number, such as 10 16 .

[0053] (2) Design optimization variables

[0054] The optimization variables include:

[0055] 1) The time T taken by the service vehicle to complete all the approach and orbit change tasks;

[0056] 2) The approach sequence S of the service spacecraft to n target satellites 1 , S 2 ,…,S n , the number of variables is n;

[0057] 3) The starting orbit change time t of the service spacecraft approaching the i-th target satellite i1 , start the track change speed increment Δvix1 , Δv iy1 , Δv iz1 , end track change time t i2 , the speed increment Δv at the end of the track change ix2 , Δv iy2 , Δv iz2 , the number of variables for each target satellite is 8, and the number of variables for n target satellites is 8n.

[0058] The total number of optimization variables is 1+9n. It can be seen that as the number of target satellites increases, the number of optimization variables increases sharply, resulting in great difficulties in the optimization convergence of the "one-to-many" approach change strategy.

[0059] (3) Five-layer optimized design

[0060] In order to solve the problem of difficulty in optimization convergence due to the large number of optimization variables, a five-layer optimization design method is adopted:

[0061] 1) First layer: Generate total track change time

[0062] The total orbit change time is the time T taken by the service vehicle to complete the orbit change mission for all n target satellites.

[0063] 2) Second layer: Generate approach sequence

[0064] Number the n target satellites as 1, 2, ..., n, and set the approach sequence S 1 , S 2 ,…,S n The value range is [0,1]. During the optimization iteration, the service aircraft is based on S 1 , S 2 ,…,S n The target satellites with corresponding numbers are approached from first to last in the order from small to large.

[0065] 3) The third layer: Generate the end time of approaching orbit change for each target satellite

[0066] When the second-layer approach sequence is given, the approach change end time for each target satellite is generated according to the sequence.

[0067] It is a very complex operation for a service spacecraft to implement close docking with a target satellite, which usually needs to be carried out under the support of ground measurement and control. At the same time, the service spacecraft needs to perform precise measurements of the target satellite during the approach process, which requires the support of direct sunlight observation conditions. Therefore, the end time of the approach and orbit change of the target satellite must meet the requirements of being in both the measurement and control area and the direct sunlight observation area.

[0068] Within the total time limit T for the service aircraft to complete the approach to all target satellites, the target satellite's measurement and control window interval can be calculated based on the position of the ground measurement and control station, the measurable azimuth and pitch angle, and the orbital parameters of the target satellite. The target satellite's illumination window interval can be calculated based on the orbital parameters of the target satellite. The measurement and control window interval and the illumination window interval are intersected to obtain the service aircraft's approach window interval for the target satellite. Let the number of approach windows for the i-th target satellite be n wi , then the approachable window interval can be expressed as:

[0069]

[0070] Where t is1 Indicates the start time of the first window, t is2 Indicates the end time of the second window, and so on.

[0071] According to the order in which the target satellites are approached by the serviced spacecraft, the value range of the optimization variable for the end time of the approach of the i-th target satellite is:

[0072]

[0073] Where b 1 It can be a decimal close to 1, such as 0.75; b 2 It can be a decimal greater than 1, such as 1.25.

[0074] It can be seen from the above formula that in the optimization iteration, the optimization variable of the approach change end time will take the value of a floating point number greater than 1. Let its integer part be k, which means that the approach change end time belongs to the kth approach window of the target satellite, and the decimal part is p, which means the proportion of the approach change end time in the window. For example, the optimization variable value of 2.6 means that the approach change end time is t is2 +0.6×(t ie2 -t is2 ).

[0075] 4) The fourth layer: Generate the start time of approaching orbit change for each target satellite

[0076] Assume that the time it takes for the service vehicle to complete the service operation on the i-th target satellite is τ i , which is calculated as:

[0077] τ i =λ i-1 +h i-1

[0078] Where λ i-1 The time when the service vehicle approaches the i-1th target satellite and the orbit change ends. i-1The operation time used by the service vehicle to provide maintenance services to the i-1th target is used as the input value of the trajectory change mission simulation.

[0079] The value range of the optimization variable for the start time of approaching orbit change for each target satellite is:

[0080] [τ i ,λ i ]

[0081] 5) Fifth layer: Simulate and generate double-pulse velocity increments for approaching orbit changes of each target satellite

[0082] According to the conditions given by the previous two to four levels of optimization, the start time and end time of the service spacecraft approaching the orbit change can be known. The orbit parameters at the beginning of the service spacecraft's orbit change can be extrapolated from the orbit parameters of the service spacecraft's previous state, and the orbit parameters at the end of the service spacecraft's orbit change can be extrapolated from the orbit parameters of the target satellite at the initial moment. Then, the time and orbit parameters of the start of the orbit change, and the time and orbit parameters of the end of the orbit change are used as the initial point and end point of the orbit change. The Lambert algorithm can be used to iteratively simulate the double pulse velocity increment of the service spacecraft approaching the orbit change. The optimization goal is to minimize the sum of the moduli of the two pulse velocity increments.

[0083] When the spacecraft approaches the target satellite, it needs to conduct optical detection of the target. This requires that the spacecraft should be in a state of direct sunlight observation when approaching the target satellite. When the service spacecraft approaches, the angle between the vector from sunlight to the service spacecraft and the vector from the target satellite to the service spacecraft is calculated. If the angle is less than the direct sunlight constraint angle (for example, 90°), it is determined that the direct sunlight observation condition does not hold. The sum of the two pulse velocity increment moduli is set to a large number (for example, 1×10 16 ).

[0084] 6) Five-level optimization method

[0085] After the first-layer optimization variables obtain fixed values ​​within the value range, the second layer iterates and optimizes according to the fixed values ​​of the first layer. During the second-layer optimization process, the third layer needs to iterate according to the fixed values ​​of the optimization variables of the second layer. During the third-layer optimization process, the fourth layer needs to iterate according to the fixed values ​​of the optimization variables of the third layer. During the fourth-layer optimization process, the fifth layer needs to iterate according to the fixed values ​​of the optimization variables of the fourth layer.

[0086] When the fifth layer completes the iteration, the fourth layer changes the value of the optimization variable, and the fifth layer iterates again according to the value of the fourth layer, until the fourth layer completes the iteration of all values; when the fourth layer completes the iteration, the third layer changes the value of the optimization variable, and the fourth layer iterates again according to the value of the third layer, until the third layer completes the iteration of all values; when the third layer completes the iteration, the second layer changes the value of the optimization variable, and the third layer iterates again according to the value of the second layer, until the second layer completes the iteration of all values; when the second layer completes the iteration, the first layer changes the value of the optimization variable, and the second layer iterates again according to the value of the first layer, until the first layer completes the iteration of all values. At this time, the optimization simulation is completed, and the time-saving and optimal orbit change strategy for the service spacecraft approaching multiple target satellites can be obtained.

[0087] Example

[0088] According to this method, an optimization simulation example of the track change control strategy is given, such as Figure 1-5 shown.

[0089] The “one-to-four” approach change of the service vehicle to the four target satellites is simulated and optimized under the constraint of 1200m / s orbit change speed increment.

[0090] At the initial time (recorded as time 0), the orbital elements of the service spacecraft and the target satellite are shown in the following table:

[0091] Table 1 Orbital elements of the service vehicle and target satellite

[0092]

[0093] After optimization calculation, the approach sequence of the service spacecraft to the four target satellites is: target satellite 3, target satellite 4, target satellite 2, and target satellite 1.

[0094] The service vehicle faced four target satellites and made eight orbit changes, consuming a total speed increment of 1181.68 m / s in 111.4 hr.

[0095] Table 2 Orbit change time and orbit change speed increment of service aircraft

[0096]

[0097] Although the present invention has been disclosed as above in the form of a preferred embodiment, it is not intended to limit the present invention. Any person skilled in the art may make possible changes and modifications to the technical solution of the present invention by using the methods and technical contents disclosed above without departing from the spirit and scope of the present invention. Therefore, any simple modifications, equivalent changes and modifications made to the above embodiments based on the technical essence of the present invention without departing from the content of the technical solution of the present invention shall fall within the protection scope of the technical solution of the present invention.

Claims

1. A method for simulating a time-saving orbit change mission for a "one-to-many" approach to a space target, characterized in that include: The service spacecraft uses a double-pulse orbit change strategy to approach each target satellite in sequence and determine the objective function of the orbit change with the least time. Set the optimization level, where The first layer sets the total orbit change time optimization variable, which is the time T that the service spacecraft takes to complete the orbit change mission of approaching all n target satellites. The second layer: setting the approach sequence optimization variables, and determining the approach sequence of the service spacecraft to the target satellite by the variable size; The third and fourth layers: within the total time limit for the service spacecraft to complete the approach to all target satellites, determine the value range of the optimization variable of the approach change end time for each target satellite that satisfies both the measurement and control window and the illumination window; determine the value range of the optimization variable of the approach change start time for each target satellite; The fifth layer: Under a given approach sequence, the given approach change end time and approach change start time are optimized according to the above objective function to obtain the optimization result; According to the set optimization level order, start optimization by level, and select the best result from the optimization results to obtain the most time-saving track change strategy; In the fifth level of optimization iteration, the optimization algorithm is used to iteratively simulate the double pulse velocity increment of the service spacecraft approaching the orbit change. After the simulation, the angle between the vector from the sunlight to the service spacecraft and the vector from the target satellite to the service spacecraft is calculated when the service spacecraft approaches. If the angle is less than the front light constraint angle, it is determined that the front light observation condition does not hold, and the sum of the double pulse velocity increment modulus is set to a very large number, the order of magnitude of which is greater than the order of magnitude of the result after the optimization algorithm is optimized. The objective function is: Where Δv u The velocity increment consumed by the service vehicle to complete the mission of approaching and changing orbits to n target satellites, Δv g is a given consumable speed increment constraint; σ is a large positive number, at least 1e8; The value range of the optimization variable of the approaching orbit change end time for each target satellite that satisfies both the measurement and control window and the illumination window is determined by the following method: Within the total time limit for the service spacecraft to complete the approach to all target satellites, the tracking and control window interval of the target satellite is calculated based on the position of the ground tracking and control station, the trackable azimuth and pitch angle, and the orbital parameters of the target satellite; The illumination window interval of the target satellite is calculated according to the orbital parameters of the target satellite; The intersection of the measurement and control window interval and the illumination window interval is taken to obtain the window interval in which the service aircraft can approach the target satellite; According to the order in which the target satellites are approached by the serviced spacecraft, the approach and orbit change end time optimization variable is set to a floating point number with a value greater than 1, and the value range of the approach and orbit change end time optimization variable is determined according to the principle that the integer part represents which approach window the approach and orbit change end time belongs to the target satellite, and the decimal part represents the proportion of the approach and orbit change end time in the window; The value range of the optimization variable for the end time of approaching orbit change of the i-th target satellite is: Among them, the number of approachable windows of the i-th target satellite is n wi , n represents the total number of target satellites; b1 ranges from [0.5,1], and b2 ranges from (1,1.5]; The value range of the optimization variable for the start time of approaching orbit change of each target satellite is [τ i ,λ i ]; t i =λ i-1 +h i-1 In the formula, λ i-1 The time when the service vehicle approaches the i-1th target satellite and ends the orbit change; i h is the end time of the service spacecraft approaching the i-th target satellite; i-1 The operating time used by the service vehicle to provide maintenance services to the i-1th target.

2. The simulation method according to claim 1, characterized in that: Start optimizing by level as follows: The first layer of optimization variables obtain fixed values ​​within the value range, the second layer iterates and optimizes according to the fixed values ​​of the first layer, during the second layer optimization process, the third layer needs to iterate according to the fixed values ​​of the optimization variables of the second layer, during the third layer optimization process, the fourth layer needs to iterate according to the fixed values ​​of the optimization variables of the third layer, during the fourth layer optimization process, the fifth layer needs to iterate according to the fixed values ​​of the optimization variables of the fourth layer.

3. The simulation method according to claim 2, characterized in that: The optimization order is as follows: When the fifth layer completes the iteration, the fourth layer changes the value of the optimization variable, and the fifth layer iterates again according to the value of the fourth layer, until the fourth layer completes the iteration of all values; when the fourth layer completes the iteration, the third layer changes the value of the optimization variable, and the fourth layer iterates again according to the value of the third layer, until the third layer completes the iteration of all values; when the third layer completes the iteration, the second layer changes the value of the optimization variable, and the third layer iterates again according to the value of the second layer, until the second layer completes the iteration of all values; when the second layer completes the iteration, the first layer changes the value of the optimization variable, and the second layer iterates again according to the value of the first layer, until the first layer completes the iteration of all values, at which point the optimization simulation is completed, and the time-saving and optimal orbit change strategy for the service vehicle approaching multiple target satellites is obtained.

4. The simulation method according to claim 1, characterized in that: The objective function is to ensure that the fuel consumed by the service aircraft to complete all target approach tasks is less than the given fuel constraint and the absolute value of the difference is the smallest.

5. The simulation method according to claim 1, characterized in that: The approach sequence is set up in the following way: Number the n target satellites as 1, 2, ..., n, and set the approach sequence S1, S2, ..., S n The value range is [0,1], and the service aircraft is divided into S1, S2, ..., S n The target satellites with corresponding numbers are approached from first to last in the order from small to large.

Images