Next Article in Journal
Skill-Learning-Based Trajectory Planning for Robotic Vertebral Plate Cutting: Personalization Through Surgeon Technique Integration and Neural Network Prediction
Previous Article in Journal
A Multi-Objective Optimization Framework That Incorporates Interpretable CatBoost and Modified Slime Mould Algorithm to Resolve Boiler Combustion Optimization Problem
Previous Article in Special Issue
Multi-Level Thresholding Color Image Segmentation Using Modified Gray Wolf Optimizer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Space Telescope Scheduling Approach Combining Observation Priority Coding with Problem Decomposition Strategies

1
School of Information Science and Engineering, Jiaxing University, Jiaxing 314001, China
2
School of Science, Jiangxi University of Science and Technology, Ganzhou 341000, China
3
School of Computer Science and Cyber Engineering, Guangzhou University, Guangzhou 510006, China
4
Institute of Information Network and Artificial Intelligence, Jiaxing University, Jiaxing 314001, China
5
National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China
*
Authors to whom correspondence should be addressed.
Biomimetics 2024, 9(12), 718; https://doi.org/10.3390/biomimetics9120718
Submission received: 17 October 2024 / Revised: 14 November 2024 / Accepted: 17 November 2024 / Published: 21 November 2024

Abstract

:
With the increasing number of space debris, the demand for telescopes to observe space debris is also constantly increasing. The telescope observation scheduling problem requires algorithms to schedule telescopes to maximize observation value within the visible time constraints of space debris, especially when dealing with large-scale problems. This paper proposes a practical heuristic algorithm to solve the telescope observation of space debris scheduling problem. In order to accelerate the solving speed of algorithms on large-scale problems, this paper combines the characteristics of the problem and partitions the large-scale problem into multiple sub-problems according to the observation time. In each sub-problem, a coding method based on the priority of the target going into the queue is proposed in combination with the actual observation data, and a decoding method matching the coding method is designed. In the solution process for each sub-problem, an adaptive variable neighborhood search is used to solve the space debris observation plan. When solving all sub-problems is completed, the observation plans obtained on all sub-problems are combined to obtain the observation plan of the original problem.

1. Introduction

The space debris monitoring network regularly tracks and catalogs over 28,000 pieces of debris, and the total mass of artificial objects in low Earth orbit exceeds 9200 tons [1,2,3]. Given the harm of space debris to the operational safety of spacecraft and the potential impact threat of near-Earth objects on Earth, the demand for space debris detection and early warning in various countries is constantly increasing [4,5]. Space observation equipment is being constructed in various countries worldwide to meet space observation needs [6,7,8]. In 2016, the Five-hundred-meter Aperture Spherical Radio Telescope (FAST), the world’s largest and most sensitive single-aperture radio telescope built independently by China, was fully completed and put into operation. The scheduling of FAST requires maximizing the number of observable proposals and the overall scientific priority and minimizing the slew cost generated by telescope shifting while considering the constraints. The researchers modeled the problem as a minimum-cost maximum-flow problem, designed a method based on maximum matching edge detection to reduce the problem size, and proposed a backtracking algorithm to minimize the transformation cost of the optimal scheduling [9,10,11,12]. Compared with individual telescopes, telescope arrays allow for much greater observational accuracy and range by synchronizing multiple telescopes.
Researchers have done a great deal of research on the telescope scheduling algorithm to meet the needs of the observing equipment. The Gravitational Wave Optical Transient Observer prototype instrument (GOTO) was inaugurated in July 2017 in La Palma, Canary Islands [13,14]. GOTO’s telescope scheduling system consists of several independent control processes, with a master named ‘pilot’ overseeing the other processes. Observations are determined by an instantaneous scheduler instructing the ‘pilot’ on the targets to observe in real time and providing fast follow-up of transient events [13]. The scheduler for the GLORIA telescope network was designed with three algorithms to maximize the total network acceptance rate and minimize the time elapsed between observation submissions and results [15]. The first algorithm is based on the weather forecasting of the telescope position, the second algorithm is based on fuzzy logic using different input parameters, and the third algorithm is based on predicting the conditional probability that each telescope will receive an observation. After that, GLORIA researchers explore new machine learning methods, such as neural networks, support vector machines, etc., comparing these methods with the three algorithms mentioned above [15]. A prediction-type scheduler is used in Algeria’s National Aures Observatory (NAO), which solves the scheduling problem using a multi-objective genetic algorithm (named NSGA-II) based on Pareto optimality [16]. The INO340 telescope of the Iranian National Observatory aims to minimize the idle time of the telescope and reduce its mechanical motion costs while obtaining the best quality image results [17]. INO340 uses genetic algorithms to consider predictable factors that affect observation conditions and obtain the optimal scheduling plan [17]. In 2023, Zhang et al. proposed a multilevel scheduling model for the time-domain survey telescope array scheduling problem. They encapsulated the functionality in software with a hierarchical architecture, developing a flexible framework and proposing an optimal metric to maintain uniform coverage and efficient time utilization from a global perspective [18].
The current research has fewer concerns about reducing the dimensionality and improving the solution efficiency for large-scale telescope scheduling problems. However, the above problems often exist in the scheduling problems of the actual telescope observation of space debris. In this paper, we design a problem decomposition strategy to reduce the dimensionality and solution time of the problem according to the problem characteristics. By drawing on similar job shop scheduling problems, mutually cooperative coding and decoding methods are proposed to enhance problem-solving efficiency. The performance of the algorithm and the proposed strategy are verified using adaptive variable domain search in conjunction with the above strategy on real-world arithmetic cases provided by the Observatory.

2. Telescope Observation Scheduling Problem

Telescope observation of space debris is a typical practical scheduling problem, where it is necessary to determine which telescope will observe which space target at what time before telescope observation. The role of algorithmic scheduling in the process of telescope observation of space debris is crucial, and this has attracted a large number of researchers to carry out related research work. The model proposed in this paper is as follows.
We suppose there are N targets in space. For any target n { 1 , , N } , each target is assigned a priority level l n , where l n 1 , , R and R is a constant (in this article, it is 9) representing its importance. These levels are used to rank and allocate telescope resources among the N independent targets. The bigger the number of observation levels of the target the higher the observation value. Each target can be observed only within the time window, denoted as [ O n , D n ] , and D n O n T n . O n is the initial time when the target can be observed, D n is the deadline, and T n is the time taken by the telescope to observe the target n. This paper assumes that the time T n required for a telescope to observe any target is 90 s.
Let there be a total of M telescopes on the ground and any telescope m { 1 , , M } . A telescope can observe only one target at a time, and each target can be assigned to only one telescope for observation. When a telescope starts observing a target, it needs to observe the target for 90 s before it can observe the next target, with the switching time between telescope targets being ignored. When a target has already been observed, all telescopes will not observe it again. We define a Boolean variable x n , m , where x n , m = 1 to indicate that target n is in the observation queue of telescope m, and x n , m = 0 to indicate that target n is not in the observation queue of telescope m. S n and C n are the start and end observation times of target n in the telescope observation queue. We assume that q is a target of { 1 , , N } and n q . C q is the observation end time of q. The calculation formulas for S n and C n are as follows:
S n = O n , i f i = 1 N z i , n m = 0 , x n , m = 1 max O n , C q , i f z q , n m = 1 , x n , m = 1 , n u m ( n , m ) n u m ( q , m ) = 1 max O n , T m , i f j = 1 N z n , j m = 0 , x n , m = 1
C n = S n + T n
where z q , n m is a Boolean variable, and z q , n m = 1 is denoted by the fact that both target q and target n are in the observation queue of telescope m, and q precedes n. In all other cases z q , n m = 0 . n u m n , m denotes the observing order when target n is inserted into the current observing queue of telescope m. The meaning of n u m ( n , m ) n u m ( q , m ) = 1 is that both target q and target n are in the current observation queue of telescope m, and target q is the previous neighboring target of target n at the position to be inserted in the current observation queue of telescope m. z q , n m and n u m ( n , m ) are shown in Figure 1. The horizontal axis in the figure is the time axis, representing the starting and ending observation times of the target in the telescope observation queue. The vertical axis represents the telescope number, indicating that targets q, n, i, and j are in the observation queue of telescope m. S q represents the start observation time of target q, and C q represents the end observation time of target q. T j is the observation time of target j in the telescope queue, and T m represents the end observation time of the observation queue for telescope m. Due to q preceding n, z q , n m = 1 and z n , q m = 0 . In addition, q is the first target observed in the observation queue of telescope m, so n u m q , m = 1 , and, similarly, n u m n , m = 2 . We define T m as the end observation time of the last target in the observation queue of telescope m, which updates with the change of the observation queue. When there is no target in telescope m’s queue, we set the T m value to 0.
As shown in Equation (1), the values of S n are divided into three cases:
(1)
The target is inserted at the head of the observing queue. Target n is added to the current observation queue of telescope m as the head observation target (see Figure 6 in Section 3.3.1 for details). At this point, target n is at the head of the queue, so, for any target i (regardless of whether it is in the queue of m), i { 1 , , N } , i n , there is z i , n m = 0 , i.e., i = 1 N z i , n m = 0 ,   x n , m = 1 . To simplify the experiment and select S n as early as possible, let S n = O n .
(2)
The target is inserted into the middle of the observation queue. Target n is inserted into the middle of the observation queue of telescope m. At this point, target q is the previous target adjacent to the insertion position of target n (see Figure 7 in Section 3.3.1 for details), i.e., z q , n m = 1 , x n , m = 1 and n u m ( n , m ) n u m ( q , m ) = 1 . We simplify the experiment to obtain the value of S n as early as possible while satisfying the constraints; then, the value of S n is max ( O n , C q ) .
(3)
The target is inserted at the end of the observing queue. The target n is added to the current observing queue of telescope m as the tail of the queue (see Figure 8 in Section 3.3.1 for details); so, for any target j (regardless of whether it is in the queue of m), j { 1 , , N } , j n , there is z n , j m = 0 , i.e., j = 1 N z n , j m = 0 , x n , m = 1 . The value of S n is max ( O n , T m ) .
Equation (2) represents the calculation method for the end observation time C n . After determining S n , T n is set to 90s, i.e., C n = S n + T n . If the target is inserted into the tail of the observation queue, T n also needs to be updated, that is, T m = C n .
It should be noted that, in the insertion method designed in this article, every time a new target is inserted into the telescope observation queue, a free time window that meets the constraints in the observation queue will be selected (without affecting the start and end observation times of other targets that have already entered the observation queue). The start observation time should be set as early as possible. If the target does not have a free window time in the queue that meets the constraints, the target will not be inserted into the observation queue of the current telescope.
In order to observe as many high-observation-value targets as possible, the telescope scheduling goal is the total target successful observation value L, denoted as
L = n = 1 N m = 1 M x n , m · l n
The telescope scheduling model is established as follows:
max L = n = 1 N m = 1 M x n , m · l n
C n S q + z q , n m · Y
C q S n + z n , q m · Y
z n , q m + z q , n m = 1 · Y
m = 1 M x n , m 1
C n = S n + T n
S n O n
D n C n
x n , m = { 0 , 1 }
z n , q m = { 0 , 1 } , z q , n m = { 0 , 1 }
n { 1 , , N } , q { 1 , , N } , n q , l n { 1 , , R } , and m { 1 , , M } . Constraints (5), (6), and (7) indicate that a telescope can observe only one target at a time, as shown in Figure 2. Targets q, j, and i are in the observation queue of telescope m. Since the observation time of targets i and j in the queue coincide with the observable period of n, and the telescope can only observe one target at a time, n cannot join the observation queue of m.
Among them, Y is a large positive number. For the constraint (5), if z q , n m = 1 , the constraint loses its restrictiveness; for the constraint (6), if z n , q m = 1 , the constraint loses its restrictiveness.
The constraint (8) means that each target can only be assigned at most one telescope.
The constraint (9) represents the value of the observation end time C n of target n.
The constraint (10) indicates that the observation start time must be greater than or equal to the initial time when the target can be observed. Constraint (11) states that the observation end time must be less than the observable deadline of the target.
The constraints (12) and (13) represent x n , m , and z n , q are Boolean variables.
The established model aims to find the optimal scheduling policy σ * = x n , m * , S n * , ( n { 1 , , N } , m { 1 , , M } ) to maximize the total target observation value.
The simplified data of some observation targets are shown in Table 1.
In Table 1, each target has its unique corresponding target number n. In an instance, each target has at least one observable period (from initial time to deadline) and may have multiple observable periods. For example, Target 16 has three observable periods. Not every target can ultimately be observed. For each target’s observable period, the target will not be observed if no available time slot matching this period can be found in the observation schedules of all telescopes. The target level l n represents the target’s observation level; the bigger the target’s observation level value, the higher the observation value. O n is the initial time of the observable period for target n, and D n is the deadline of the observable period for target n, measured in seconds. To simplify the experiment, the observation window for each target in the article is the same for each telescope.

3. A Heuristic Approach Combining Problem Decomposition

3.1. Algorithm’s Overall Flow

According to the particularity of the telescope observation scheduling problem, the original problem can be divided into multiple sub-problems for solution. The observation queues obtained by solving the divided sub-problems do not overlap in time. After solving the sub-problems separately, the observation plans obtained can be combined to obtain the observation plan of the original problem.
The original problem can be divided into sub-problems recursively. After using the problem decomposition strategy, the algorithm flow chart is shown in Figure 3.
The algorithm’s overall execution process is as follows:
1.
A sub-problem is solved to obtain an observation plan for that sub-problem. If the space debris in that observation plan appears again in other unsolved sub-problems, it is removed from the unsolved sub-problems.
2.
Repeat the above steps until all sub-problems are solved and corresponding observation plans are obtained.
3.
The observation plan obtained by solving all sub-problems is combined to obtain the observation plan for the original problem.

3.2. Problem Decomposition Strategy

The observation target dimension of the telescope observation scheduling problem is enormous, and, on average, the observation scheduling of thousands of space targets needs to be processed daily. In order to reduce the dimension of the problem, a problem decomposition method is proposed based on the characteristics of the problem in this article:
1.
Sort all target data in ascending order according to O n , as shown in Table 2.
2.
After sorting, add the minimum value of O n and the maximum value of D n in all data and divide by 2 as the segmented time point, as shown in Figure 4.
3.
The data with O n and D n less than the segmented time point are added to the set of the first sub-problem.
4.
The data with O n and D n greater than the segmented time point are added to the set of the second sub-problem.
5.
If the segmented time point is greater than O n and less than D n , the observable period of the target is divided into two periods. One is from O n to the segment time point, and the other is from the segment time point to D n . If the length of the time period after separation is less than T n , it will be discarded. Otherwise, the previous time segment is added to the first sub-problem. The following data segment is added to the set of the second sub-problem.
Taking the data in Table 2 as an example, the minimum value of O n is 77,539 in the data of target 2378, and the maximum value of D n is 80,267 in the data of target 2121. We add 77,539 to 80,267 and divide by 2 to obtain the segmented time point 78,903.
The observation period for target 2378 is divided into two segments, 77,539 to 78,903 and 78,903 to 79,309. Both segments have a range greater than 90 s. The data for the first segment are added to sub-problem 1, and the data for the second segment are added to sub-problem 2.
The O n and D n of targets 6848, 2377, and 4259 are all less than 78,903, and the data are directly added to the set of sub-problem 1. The O n and D n of targets 4617, 6266, and 2121 are all greater than 78,903, and the data are directly added to the set of sub-problem 2. The data of the two sub-problems obtained after problem decomposition are shown in Table 3 and Table 4.
The observation plan based on the data in Table 2 on a telescope is shown in Table 5. The observation plan obtained from Table 3 and Table 4 is also the same as the observation plan in Table 5.
It is demonstrated through experiments in the subsequent chapters that the use of the problem decomposition strategy significantly reduces the algorithm solution time when compared to the case where the solution is done directly without using the problem decomposition strategy.
It is worth mentioning that the experimental results in subsequent chapters show that using the problem decomposition strategy to decompose the problem and reduce its dimensionality each time continuously saves less and less computational time. At the same time, the observation plan obtained by using a small number of times of problem decomposition is better than the one obtained by using multiple times of problem decomposition. In the case where the 1000 pieces of data are taken from each of the two two-telescope examples, the observation plan obtained by decomposing the original problem into four sub-problems is better than that obtained by decomposing the original problem into two sub-problems. Taking instance 1 as an example, the period of the original problem is 74,143 to 113,650. After decomposing into four sub-problems, the period of sub-problem 1 is 74,143 to 84,019, the period of sub-problem 2 is 84,019 to 93,896, the period of sub-problem 3 is 93,896 to 103,565, and the period of sub-problem 4 is 103,565 to 113,650.

3.3. Coding and Decoding

The main factors affecting the scheduling of telescope observations in the problem studied in this paper are as follows [19]:
1.
Level of observation targets.
2.
Target observable time window.
3.
Observation time of the telescope on the target.
The problem has many constraints, and the dimension of the solution space of the problem is vast, necessitating the minimization of the solution space as much as possible, and, at the same time, requires the efficiency of the coding operation, which is the basis for the subsequent algorithms to be able to solve the problem efficiently. Designing a reasonable decoding method with appropriate coding to reduce the solution space size is the focus of the research problem.

3.3.1. Coding Method

In order to solve the above problems, this paper uses the target observation data provided by the observatory to propose a coding method based on the priority of the target entering the telescope observation queue.
This paper takes the random initial solution 58 , 12 , 29 , 16 , 13 of the problem obtained from the data in Table 1 as an example. Each variable in the solution indicates the number of the observation target, and the earlier the target appears in the solution, the earlier it is added to telescope observation by the decoding method. According to the target entry priority in the above solution, the sorted data to be added to the observation queue are obtained, as shown in Table 6.
Taking Table 6 as an example, the coding method in this paper reduces the problem size of 8 to the solution coding of five decision variables, which reduces the solution space and also facilitates the operators to operate efficiently on the solution, which lays the foundation for the subsequent algorithms to solve the problem efficiently. Which telescope the subsequent target is added to and when it is added are handled by a decoding method that matches the coding method. The coding part is only responsible for generating a data table sorted according to priority based on the target priority in the solution, and the decoding method reads the data table to decide which telescope the target joins and when it joins the telescope.

3.3.2. Decoding Method

Applying the decoding approach has achieved significant results in similar scheduling problems (job shop scheduling) and, together with the appropriate coding approach, it can significantly reduce the search space [20,21]. The essence of the algorithm operation solution in the telescope scheduling problem is to allow as many targets as possible to join the observation queue of the telescope while selecting high observation value targets as much as possible to join the observation queue. Based on this property of the telescope scheduling problem, this paper proposes a decoding method that matches the abovementioned coding method. The flowchart of the decoding method is shown in Figure 5.
The judgment conditions for the four cases are given first:
1.
Queue head insertion condition: the target was not observed, and S q O n T n , where q is the first target in the observation queue of the current telescope m, and the observation start time is S q . The queue head insertion is shown in Figure 6.
2.
Insertion condition in the middle of the queue: S k C j T n and S k O n T n are satisfied simultaneously, and the target is not observed, where j is the last neighboring observed target at the location to be inserted, and k is the next neighboring observed target at the location to be inserted. The team insertion is shown in Figure 7.
3.
Insertion conditions at the end of the queue: the target is not observed, and D n T m T n . The end of the queue is inserted, as shown in Figure 8.
4.
Replacement condition: the target is not observed while satisfying l n < l p , O n S p and C p D n , where p is the observed target in the current queue. The replacement insertion is shown in Figure 9.
It is worth mentioning that, after a large number of repeated experiments, it was found that the observation queue obtained by the decoding method is only related to the order in which the target enters the queue, and is not related to the order in which the three insertion methods are used.
The execution steps of the decoding method are as follows:
1.
Obtain the first data in the list of observed target data.
2.
Determine whether the acquired target data can be added to the first telescope observation queue. Judge whether the target can be added to a telescope queue in order to perform the following four steps: (A) Determine whether the target can be inserted into the tail of the observation queue, and if the insertion conditions are met, insert it into the tail of the observation queue. (B) Determine whether the target can be inserted into the observation queue head, and if the insertion conditions are met, insert it into the observation queue head. (C) Determine whether the target can be inserted into the middle of the observation queue, and if the insertion criteria are met, insert it into the middle of the observation queue. (D) If the target is inserted into the current telescope observation queue, the target number will be added to the collection of observed queues; otherwise, continue to judge whether the target can be added to the observation queue of the second telescope until the target is inserted into the telescope queue or each telescope queue can not be inserted, the end of the insertion process.
3.
After completing step 2 above, regardless of whether the target is added to the observation queue or cannot be added to the observation queue, read the target data in the following data table and continue with step 2 until all data have been executed in step 2.
4.
After all target data in the data table have been inserted and judged, unobserved targets are added to the unobserved target set.
5.
Determine whether each unobserved target can replace the observed low-value target. If the conditions are met, the low-value target will be removed from the observation queue, and the high-value target will occupy the position of the removed low-value target in the observation queue. Finally, the scheduling ends after all unobserved targets have executed the replacement judgment.
Multiple replacement rounds were attempted in the experiment, testing the strategy mentioned in Step 5 until no low-value targets could be replaced in the telescope queues. This improved the total observation value but significantly increased the algorithm’s computation time. Therefore, fewer replacement rounds strike a better balance between total observation value and algorithm-solving time.
For the convenience of setting up only one telescope, the solution 58 , 12 , 29 , 16 , 13 in Section 3.3.1 uses the decoding method proposed in this article to obtain the observation queue formation process, as shown in Table 7 and  Table 8.
Table 7 shows the observation queue after the insertion of all data tables. The process by which targets enter the observation queue is as follows:
1.
Target 58 is the first to enter the observation queue, S n is 76,925, and C n is 77,015.
2.
Target 12 does not fulfill the conditions to enter the observation queue.
3.
Target 29 enters the end of the observation queue with an S n of 106,755 and a C n of 106,845.
4.
Target 16 enters the observation queue header with an S n of 75,833 and a C n of 75,923.
5.
Target 13 enters the observation queue queue with an S n of 77,015 and a C n of 77,105.
Table 8 shows the observation queue after performing the replacement operation. The observed value of the unobserved target 12 is higher than the replaceable observed target 58, executing the replacement operation.

3.4. Adaptive Variable Neighborhood Search Algorithm

Adaptive variable neighborhood search (AVNS) has achieved better results in many huge neighborhood problems [22,23,24,25,26,27,28]. By combining several operators to dynamically compete for neighborhood search opportunities, adaptive variable neighborhood search has good adaptability and scalability [29,30,31]. This paper’s adaptive variable neighborhood search uses the insertion, commutative, and two-opt operators. The operator dynamic competition strategy from the literature [24] is adopted, and the dynamic competition strategy is as follows: the insertion operator, the commutative operator, and the two-opt operator are set with weights W 1 , W 2 , and W 3 , respectively, and the initial value of W 1 , W 2 , and W 3 is 1.0. When an operator searches for a better solution, its weight increases by 1.0. If an operator fails to find a better solution, the weight of this operator remains unchanged, while the weight of other operators increases by 0.5. The W update formula is as follows:
W k = W k + 1.0 W i = W i , i 1 , 2 , 3 k i f L ( n e w _ s o l u t i o n ) > L ( o l d _ s o l u t i o n ) W k = W k W i = W i + 0.5 , i 1 , 2 , 3 k i f L ( n e w _ s o l u t i o n ) L ( o l d _ s o l u t i o n )
where L ( s o l u t i o n ) is the total value of successful observations of the solution (Formula (3)). In this formula, the following conditions apply:
1.
If an operator k finds an improved solution (i.e., L ( n e w _ s o l u t i o n ) > L ( o l d _ s o l u t i o n ) ), then W k is increased by 1.0, and other operator weights remain unchanged.
2.
If the operator does not yield a better solution, then W k remains the same, and each other operator’s weight increases by 0.5.
Each neighborhood search uses a roulette wheel method to select operators, and the probabilities of the three operators being selected are P 1 , P 2 , and P 3 . The calculation formula is as follows [24]:
P j = W j j = 1 3 W j , j = 1 , 2 , 3
Among them, W j is the weight corresponding to operator j, and j = 1 3 W j is the sum of the weights of all operators.

3.4.1. AVNS Algorithm Process

The flowchart of the AVNS in this paper is shown in Figure 10.
The algorithm flow is as follows:
1.
Initialize weights and generate a certain number of initial solutions.
2.
Total observed values are calculated using the decoding method proposed in Section 3.3.2.
3.
Determine if the algorithm termination condition is met.
4.
Select a roulette wheel operator for each individual in the population and search for neighborhoods using the selected operator.
5.
Total observed values are calculated using the decoding method proposed in Section 3.3.2.
6.
Update the population and save the optimal solution.
7.
Repeat steps 3 to 6 until the termination condition is met and the algorithm terminates.
During the initialization process, we first set the initial weights of all three operators to 1.0. Subsequently, we generate initial solutions based on population size. When the algorithm generates an initial solution, it first forms a set using all the target numbers in the instance. The set obtained from the example in Table 1 is 12 , 13 , 16 , 29 , 58 . Subsequently, randomly sorting the target numbers in the set to generate a sequence of target numbers yields an initial solution, such as 29 , 58 , 12 , 16 , 13 .
The pseudo-code of the algorithm in this paper is shown in Algorithm 1.
Algorithm 1 AVNS
Require: Best solution
Ensure: Target data, pop_size, telescope_number, gen, i=0
 1: Initial_population()
 2: Initialize_operator_weights()
 3: for Each solution do
 4:     Section 3.3.2 decoding_method()
 5: end for
 6: while i<gen do
 7:     for Each solution do
 8:         Adaptive_variable_neighborhood_search()
 9:         Section 3.3.2 decoding_method()
 10:     end for
 11:     Update_population()
 12:     Update_operator_weights()
 13:     i+=1
 14: end while
pop_size is the population size, telescope_number is the number of telescopes, and gen is the maximum number of iterations. Section 3.3.2’s decoding_method() is used to calculate the total observation value of each solution.

3.4.2. Insert Operator

The insertion operator randomly selects a target number in the solution, removes the selected target number from its original position, and inserts it into other positions in the solution to generate a new neighborhood solution. The process of the insertion operator to generate a new neighborhood solution is shown in Figure 11.
Taking the solution in Section 3.3.1 as an example, the insertion operator selects target number 12 in the second position in the solution and inserts target 12 into target number 16 to produce a new solution.

3.4.3. Commutative Operator

The commutative operator generates a new neighborhood solution by exchanging the target numbers at two different positions in the solution. The process of generating new neighborhood solutions through the commutative operator is shown in Figure 12.
Taking the solution in Section 3.3.1 as an example, the commutative operator swaps the positions of the target number 12 and the target number 16 in the solution, producing a new solution.

3.4.4. Two-Opt Operator

The two-opt operator selects two different positions in the solution and flips the sequence of target numbers between the two positions to generate a new neighborhood solution. The process of generating a new neighborhood solution by the two-opt operator is shown in Figure 13.
Taking the solution in Section 3.3.1 as an example, the two-opt operator selects the queue from target number 12 to target number 13 and inverts the target numbers to form a new solution.

3.4.5. Complexity Analysis

We suppose that there are m telescopes and n space debris in the telescope observation of space debris scheduling problem, and the population size of the algorithm is N. The computational complexity of the total successful observation value is O n m . The time complexity of generating N random initial solutions is O N n . Therefore, the time complexity of generating initial solutions and calculating the total successful observation value for each solution is O N m n 2 .
In the adaptive variable neighborhood search stage, the time complexity of all three operators is O m n 2 , so the time complexity of this stage is O N m n 2 . Accordingly, the computation complexity of the AVNS is equal to O N m n 2 .

4. Experimental Studies

In order to verify the performance of the algorithm in this paper, the experiment used ten spatial target observation instances provided by relevant units, each containing about 9000 target observation information units. According to the problem’s difficulty, instances with more than 800 data points are considered large-scale problems, instances with 500–800 data points are considered medium-scale problems, and instances with less than 500 data points are considered small-scale problems. The computer’s operating system was Windows 11, the CPU was i7-12700h, and Python was the programming language. Due to the simple and efficient characteristic of the adaptive variable neighborhood search algorithm, the algorithm had only one parameter, P O P _ S I Z E , and population size was set to 10. To obtain the algorithm’s source code in this article, please visit https://github.com/shixin63/AVNS (accessed on 10 November 2024).
In order to verify the validity of the decoding approach proposed in this paper, the decoding approach in this paper was compared with the greedy decoding approach. The greedy decoding method process is as follows:
1.
Sort the data in ascending order according to when the target started to be observed.
2.
Determine whether the target can be inserted into the end of the observation queue. If it meets the insertion conditions, insert the end of the observation queue.
3.
If the target is inserted into the end of the observation queue of the current telescope, then the target number will be added to the set of observed queues; otherwise, continue to determine whether the target can be added to the end of the observation queue of the second telescope until the target is inserted into the end of the telescope queue or each telescope queue cannot be inserted into and the end of the insertion process is over.
4.
Read the target data in the following data table and continue with steps 2 and 3 until all data have been executed through steps 2 and 3, ending the decoding.
We took the 2000 target observations from each of the ten instances, set the number of telescopes to two, and used greedy decoding and the decoding method proposed in this paper for the data in each instance, respectively. The experimental results are shown in Table 9.
The experimental results in Table 9 show that the decoding method proposed in this paper achieves significantly better decoding results than the greedy decoding method in all instances.
In order to verify the effectiveness of the problem decomposition strategy proposed in this article, an ablation experiment was conducted while ensuring that the algorithm’s operating environment and parameters remained unchanged. The experiment was divided into three situations: (1) Decompose the original problem into four sub-problems. (2) Decompose the original problem into two sub-problems. (3) Solve the original problem directly without using the problem decomposition strategy. Each case was calculated independently ten times on ten instances. The number of telescopes was set to two, each calculation instance took 1000 pieces of data, and the algorithm termination condition was to run for 50 generations. The experimental results are shown in Table 10. The time in the table represents the average time (in seconds) for the algorithm to run independently ten times on each instance from initialization to obtaining all observation plans.
The data in Table 10 show that when solving four sub-problems, the solution speed is the fastest, and the solution quality is higher than when directly solving the original problem, accounting for about 29% of the original problem-solving time. The best solution is achieved when solving two sub-problems, with a solution time of approximately 47% of the original problem. From this, it can be seen that the problem decomposition strategy proposed in this article can effectively improve the algorithm’s solving efficiency.
The experimental results show that, if each sub-problem can be solved in parallel, the algorithm’s execution time can be significantly shortened. However, as the problem is continually decomposed into smaller sub-problems, the reduction in solving time becomes progressively more minor, and some local optimality issues may arise. If parallel computing is to be used, after selecting the appropriate number of sub-problems, the communication between different sub-problem-solving processes must also be addressed. For example, when a target n is assigned to the observation queue in sub-problem 1, there is no need to add target n to the observation queue while solving other sub-problems. Some studies have explored communication methods in the parallel solving of scheduling problems but, so far, suitable methods for communication between sub-problem solutions in this paper have yet to be identified.
The algorithm in this article independently calculates ten times on ten instances, with the number of telescopes set to five. Each instance takes 2000 pieces of data. Due to the larger instance scale compared to the ablation experiments in Table 10, the original problem was divided into eight sub-problems on an instance with 2000 data points to speed up the algorithm’s solution. The termination condition of the algorithm was to run for 50 generations. The experimental results are shown in Table 11, where target number represents the number of observation targets included in 2000 data points, the average count represents the average number of observation targets after ten calculations, the average represents the average observation value in 10 calculations, and the maximum represents the maximum observation value obtained in 10 calculations.
From the data in Table 11, it can be seen that the algorithm in this paper can ensure that most of the targets are observed in the case of five telescopes and 2000 pieces of data, and only a tiny part of the targets that cannot be observed due to the observation time constraints are missed. In the problem-solving process, the replacement strategy of the decoding method can try to observe high-value targets first to ensure that high-value targets enter the observation queue to the maximum extent.
The effectiveness of the algorithm strategy proposed in the article has been demonstrated through previous experiments. We continued to compare the algorithm presented in this article with Genetic Algorithm (GA) [32,33] and Simulated Annexing Algorithm (SA) [34,35]. It is worth mentioning that the fewer dispatchable telescopes there are, the less observation queue space is available for the target, and the more complex the solution to the corresponding problem is. Let the number of telescopes be two; all algorithms are computed independently ten times on each of the ten instances, and 1000 data points are taken for each instance.
The algorithm in this article divided the problem into four sub-problems, with a population size of 10 and an iteration count of 50. GA used the POX crossover operator during the crossover process and the commutative operator during the mutation process. In each iteration, another solution was selected for each solution in the population based on the fitness roulette. Subsequently, the similarity between the observation queues of the two solutions was assessed. If they were dissimilar, crossover occurred; if they were similar, no crossover was performed. If a better solution was generated after crossover and added to the population, the population remained unchanged. In each round, a swapping operation was performed on each solution, where only the positions of two target numbers within the solution were exchanged, and only the better solution was accepted. The population size of GA was set to 10. The initial temperature of the SA was set to 100, with an initial temperature decrease rate of 0.98. The number of iterations for the isothermal process was set to 20, and the neighborhood search operator was selected as the swapping operator. The running time of GA and SA was the same as the CPU time required for 50 iterations of the algorithm in this article. Both comparison algorithms and AVNS used the decoding method proposed in this paper, and the experimental results are shown in Table 12 and Table 13.
As seen from the data in Table 12 and Table 13, the present algorithm ensures the scheduling of most of the targets into the observation queue, even in the more complex case of two telescopes. At the same time, the maximum value, average value, and average number of observed targets obtained by the proposed algorithm are superior to the comparison algorithms. AVNS uses three combined operators to have more robust local and global search capabilities than GA and SA. When the search falls into a local optimal state, it can adjust the weight to select the two-opt operator with a more extensive search range. When the operator needs local search, it can choose the insertion and commutative operators for search.
The fewer telescopes there are, the less observation time they provide, making it more difficult to obtain an observation plan with a higher total observation value. In order to demonstrate this phenomenon more intuitively, experiments were conducted with 10 telescopes. A total of 1000 data points were taken from each instance, and the algorithm independently calculated 10 times for each instance. The termination condition for the algorithm was that it stopped running after ten or five iterations without producing a better solution. The experimental results are shown in Table 14.
According to the experimental results in Table 14, it can be seen that, with ten telescopes compared to two telescopes, better solution results can be obtained with less solving time, and all targets in each instance are successfully observed. However, an excessive number of telescopes can cause some telescopes to remain idle for a long time, resulting in resource waste. Therefore, further research is needed to set the number of telescopes in the algorithm in future studies.
In summary, the performance of the scheduling methods and algorithms in this paper was tested by simulation experiments on 2000 data points taken in each of the ten arithmetic cases used in real work. The algorithm’s performance was also tested by taking 1000 data points in each example in the more complex case of two telescopes. The experimental results all show that the algorithms in this paper perform well in the data sets used in the actual work and can find a better solution within an acceptable time.

5. Conclusions

In this paper, a telescope scheduling model was established according to the actual needs of space observation and telescope array scheduling. A coding method based on target priority was proposed relying on the actual data, and a decoding method matching the coding method was designed. In order to solve the problem that the algorithm has difficulty searching on huge dimensional arithmetic cases, a problem decomposition strategy was designed. On this basis, an adaptive variable neighborhood search algorithm was used to solve the telescope scheduling problem, and the effectiveness of the algorithm proposed in this paper was verified.
The combination of coding and decoding methods proposed in this paper effectively reduced the size of the solution space and facilitated the algorithm solving the problem efficiently. Three insertion methods and one substitution method for the target to enter the telescope observation queue were proposed, and the effectiveness was proved in simulation experiments using actual data. The proposed problem decomposition method can split a problem into several small-scale sub-problems, and the split sub-problems can be solved separately and combined to shorten the problem-solving time. Based on the above strategy, excellent results have been achieved in solving the telescope scheduling problem using the adaptive variable neighborhood search algorithm.
The effectiveness of the model, coding method, decoding method, and problem decomposition strategy established in this paper were verified by simulation experiments using actual data, which proved that the algorithm used in this paper can meet the requirements of telescope observation scheduling. The problem decomposition strategy proposed in this article abandoned some observable periods of the target when decomposing the original problem into sub-problems, as they were less than T n after decomposition. These periods are observable and will be sought to solve this problem in future research. The next step is to conduct more in-depth research on the coding and decoding methods to further improve the search speed of the algorithm in such a large problem dimension as 9000 data points. The application of deep reinforcement learning methods in combinatorial optimization problems is a hot topic in current research, and future research work can explore the use of deep reinforcement learning to solve telescope scheduling problems [36,37,38,39]. In practice, in addition to maximizing the scientific value of observations, the scheduling of telescope observations usually considers factors such as operating costs [11,15]. Therefore, studying multi-objective algorithms for the telescope observation scheduling problem is also of great value [40,41,42,43,44,45,46].

Author Contributions

K.Z.: Conceptualization, methodology, formal analysis, investigation, and writing—original draft. B.-L.Y.: conceptualization, methodology, supervision, and writing—review and editing. X.X.: conceptualization, methodology, and writing—review and editing. Z.W.: conceptualization, methodology, supervision, and writing—review and editing. X.Z.: validation and writing—review and editing. H.J.: validation and writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Space Debris and Near-Earth Asteroid Defense project of China (KJSP2020020202), Teaching Reform Project of Zhejiang Higher Education “14th Five-Year Plan” (jg20220434), National Natural Science Foundation of China (62106055, 61703183), National Key Research and Development Program of China (2023YFC3305900, 2023YFC3305903), Natural Science Foundation of Zhejiang Province (LTGS23F030002, LGG19F030010), Guangdong Natural Science Foundation (2022A1515011825), Guangzhou Science and Technology Planning Project (2023A04J0388, 2023A03J0662), and the Qin Shen Scholar Program of Jiaxing University.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon request. There are no restrictions on data availability.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Sen, T.; Yan, W. Research on data processing method for detection of small and weak targets in space. In Proceedings of the AOPC 2022: Optical Sensing, Imaging, and Display Technology, Online, 18–19 December 2022; SPIE: Bellingham, DC, USA, 2023; Volume 12557, pp. 489–493. [Google Scholar]
  2. Zhao, M.; Li, G.; Li, H.; Li, S. Reliable Scheduling Algorithm for Space Debris Monitoring Resources Using Adaptive Multipopulation Differential Evolutionary Optimization With Reinforcement Learning. IEEE Trans. Reliab. 2022, 71, 687–697. [Google Scholar] [CrossRef]
  3. Mark, C.P.; Kamath, S. Review of Active Space Debris Removal Methods. Space Policy 2019, 47, 194–206. [Google Scholar] [CrossRef]
  4. Tang, J.; Cheng, H. The origin, status and future of space debris. Physics 2021, 50, 317–323. [Google Scholar] [CrossRef]
  5. Li, G.; Jiang, H.; Cheng, H.; Liu, J.; Yang, C.; Jiang, P. Space debris observation performance research of CSTAR at Kunlun Station in Antarctica. Adv. Space Res. 2019, 64, 1527–1536. [Google Scholar] [CrossRef]
  6. Labate, M.G.; Waterson, M.; Alachkar, B.; Hendre, A.; Lewis, P.; Bartolini, M.; Dewdney, P. Highlights of the Square Kilometre Array Low Frequency (SKA-LOW) Telescope. J. Astron. Telesc. Instruments Syst. 2022, 8, 011024. [Google Scholar] [CrossRef]
  7. Suchodolski, T. Active Control Loop of the BOROWIEC SLR Space Debris Tracking System. Sensors 2022, 22, 2231. [Google Scholar] [CrossRef]
  8. Zhang, H.; Liu, H.; Xu, W.; Lu, Z. Large aperture diffractive optical telescope: A review. Opt. Laser Technol. 2020, 130, 106356. [Google Scholar] [CrossRef]
  9. Li, W.M.; Zhou, L.J. A preliminary study on photogrammetry for the FAST main reflector measurement. Res. Astron. Astrophys. 2021, 21, 156. [Google Scholar] [CrossRef]
  10. Wang, Z.J.; Zhan, Z.H.; Yu, W.J.; Lin, Y.; Zhang, J.; Gu, T.L.; Zhang, J. Dynamic Group Learning Distributed Particle Swarm Optimization for Large-Scale Optimization and Its Application in Cloud Workflow Scheduling. IEEE Trans. Cybern. 2020, 50, 2715–2729. [Google Scholar] [CrossRef]
  11. Luo, Q.; Zhao, L.; Yu, C.; Xiao, J.; Sun, J.; Zhu, M.; Zhong, Y. Cost-efficient scheduling of FAST observations. Exp. Astron. 2018, 45, 107–126. [Google Scholar] [CrossRef]
  12. Yu, N.P.; Qian, L.; Zhang, C.P.; Jiang, P.; Xu, J.L.; Wang, J.J. HI detection of J030417.78+002827.4 by the Five-hundred-meter Aperture Spherical Radio Telescope. Res. Astron. Astrophys. 2021, 21, 100. [Google Scholar] [CrossRef]
  13. Dyer, M.J.; Dhillon, V.S.; Littlefair, S.; Steeghs, D.; Ulaczyk, K.; Chote, P.; Galloway, D.; Rol, E. A telescope control and scheduling system for the Gravitational-wave Optical Transient Observer (GOTO). In Proceedings of the Observatory Operations: Strategies, Processes, and Systems VII, Montreal, QC, Canada, 25–27 June 2014; Peck, A.B., Seaman, R.L., Benn, C.R., Eds.; International Society for Optics and Photonics, SPIE: Bellingham, DC, USA, 2018; Volume 10704, p. 107040. [Google Scholar] [CrossRef]
  14. Mong, Y.L.; Ackley, K.; Galloway, D.K.; Dyer, M.; Cutter, R.; Brown, M.J.I.; Lyman, J.; Ulaczyk, K.; Steeghs, D.; Dhillon, V.; et al. Searching for Fermi GRB optical counterparts with the prototype Gravitational-wave Optical Transient Observer (GOTO). Mon. Not. R. Astron. Soc. 2021, 507, 5463–5476. [Google Scholar] [CrossRef]
  15. López-Casado, C.; del Pulgar, C.P.; Muñoz, V.F.; Castro-Tirado, A.J. Observation scheduling and simulation in a global telescope network. Future Gener. Comput. Syst. 2019, 95, 116–125. [Google Scholar] [CrossRef]
  16. Bekhti, Y.A.; Lamraoui, K.; Seghouani, N.; Benatchba, K.; Bouzidi-Hassini, S. A scheduler for the National Aures Observatory. J. Phys. Conf. Ser. 2019, 1269, 012004. [Google Scholar] [CrossRef]
  17. Shalchian, H.; Sotoudeh, M.H.; Khosroshahi, H.G.; Ravanmehr, R.; Fatemi, S.; Altafi, H. A metaheuristic approach for INO340 telescope flexible scheduling. In Software and Cyberinfrastructure for Astronomy VII; Ibsen, J., Chiozzi, G., Eds.; International Society for Optics and Photonics, SPIE: Bellingham, DC, USA, 2022; Volume 12189. [Google Scholar] [CrossRef]
  18. Zhang, Y.; Yu, C.; Sun, C.; Shang, Z.; Hu, Y.; Zhi, H.; Yang, J.; Tang, S. A Multilevel Scheduling Framework for Distributed Time-domain Large-area Sky Survey Telescope Array. Astron. J. 2023, 165, 77. [Google Scholar] [CrossRef]
  19. Brandt, P. Designing an ngVLA Dynamic Scheduler ngVLA Computing Memo# 6. 2021. Available online: https://library.nrao.edu/public/memos/ngvla/NGVLAC_06.pdf (accessed on 15 October 2024).
  20. Zhang, C.; Guan, Z.; Liu, Q.; Shao, X.; Li, P. New scheduling type applied to solving job-shop scheduling problem. J. Mech. Eng. 2008, 44, 24–31. [Google Scholar] [CrossRef]
  21. Zhu, H.; Chen, M.; Zhang, Z.; Tang, D. An Adaptive Real-Time Scheduling Method for Flexible Job Shop Scheduling Problem With Combined Processing Constraint. IEEE Access 2019, 7, 125113–125121. [Google Scholar] [CrossRef]
  22. He, M.; Wei, Z.; Wu, X.; Peng, Y. An Adaptive Variable Neighborhood Search Ant Colony Algorithm for Vehicle Routing Problem With Soft Time Windows. IEEE Access 2021, 9, 21258–21266. [Google Scholar] [CrossRef]
  23. Xue, Z.F.; Wang, Z.J.; Zhan, Z.H.; Kwong, S.; Zhang, J. Neural Network-Based Knowledge Transfer for Multitask Optimization. IEEE Trans. Cybern. 2024, 2024, 1–14. [Google Scholar] [CrossRef]
  24. Pan, Z.X.; Wang, L.; Chen, J.F.; Wu, Y.T. A Novel Evolutionary Algorithm with Adaptation Mechanism for Fuzzy Permutation Flow-Shop Scheduling. In Proceedings of the 2021 IEEE Congress on Evolutionary Computation (CEC), Krakow, Poland, 28 June–1 July 2021; pp. 367–374. [Google Scholar] [CrossRef]
  25. Wang, Z.J.; Jian, J.R.; Zhan, Z.H.; Li, Y.; Kwong, S.; Zhang, J. Gene Targeting Differential Evolution: A Simple and Efficient Method for Large-Scale Optimization. IEEE Trans. Evol. Comput. 2023, 27, 964–979. [Google Scholar] [CrossRef]
  26. Zheng, X.L.; Wang, L. A Collaborative Multiobjective Fruit Fly Optimization Algorithm for the Resource Constrained Unrelated Parallel Machine Green Scheduling Problem. IEEE Trans. Syst. Man, Cybern. Syst. 2018, 48, 790–800. [Google Scholar] [CrossRef]
  27. Wang, Z.J.; Zhou, Y.R.; Zhang, J. Adaptive Estimation Distribution Distributed Differential Evolution for Multimodal Optimization Problems. IEEE Trans. Cybern. 2022, 52, 6059–6070. [Google Scholar] [CrossRef] [PubMed]
  28. Karakostas, P.; Sifaleras, A.; Georgiadis, M.C. Adaptive variable neighborhood search solution methods for the fleet size and mix pollution location-inventory-routing problem. Expert Syst. Appl. 2020, 153, 113444. [Google Scholar] [CrossRef]
  29. Sabar, N.R.; Kendall, G. An iterated local search with multiple perturbation operators and time varying perturbation strength for the aircraft landing problem. Omega 2015, 56, 88–98. [Google Scholar] [CrossRef]
  30. Wang, Z.J.; Zhan, Z.H.; Li, Y.; Kwong, S.; Jeon, S.W.; Zhang, J. Fitness and Distance Based Local Search With Adaptive Differential Evolution for Multimodal Optimization Problems. IEEE Trans. Emerg. Top. Comput. Intell. 2023, 7, 684–699. [Google Scholar] [CrossRef]
  31. Marinho Diana, R.O.; de Souza, S.R. Analysis of variable neighborhood descent as a local search operator for total weighted tardiness problem on unrelated parallel machines. Comput. Oper. Res. 2020, 117, 104886. [Google Scholar] [CrossRef]
  32. Hu, K.; Wang, L.; Cai, J.; Cheng, L. An improved genetic algorithm with dynamic neighborhood search for job shop scheduling problem. Math. Biosci. Eng. 2023, 20, 17407–17427. [Google Scholar] [CrossRef]
  33. Tan, S.; Yao, Q.; Li, J.; Shi, S.C. Modified genetic algorithm for China Space Station Telescope high-sensitivity terahertz detection module observation schedule. J. Astron. Telesc. Instruments Syst. 2024, 10, 037002. [Google Scholar] [CrossRef]
  34. Teramoto, K.; Morinaga, E.; Wakamatsu, H.; Arai, E. A neighborhood limitation method for job-shop scheduling based on simulated annealing. Trans. Inst. Syst. Control Inf. Eng. 2020, 33, 171–181. [Google Scholar] [CrossRef]
  35. Liu, Y.; Zhang, S.; Hu, H. A conflict avoidance algorithm for space-based collaborative stereo observation mission scheduling of space debris. Adv. Space Res. 2022, 70, 2302–2314. [Google Scholar] [CrossRef]
  36. Mazyavkina, N.; Sviridov, S.; Ivanov, S.; Burnaev, E. Reinforcement learning for combinatorial optimization: A survey. Comput. Oper. Res. 2021, 134, 105400. [Google Scholar] [CrossRef]
  37. Wang, X.; Wang, S.; Liang, X.; Zhao, D.; Huang, J.; Xu, X.; Dai, B.; Miao, Q. Deep Reinforcement Learning: A Survey. IEEE Trans. Neural Netw. Learn. Syst. 2022, 35, 5064–5078. [Google Scholar] [CrossRef] [PubMed]
  38. Li, K.; Zhang, T.; Wang, R. Deep Reinforcement Learning for Multiobjective Optimization. IEEE Trans. Cybern. 2021, 51, 3103–3114. [Google Scholar] [CrossRef] [PubMed]
  39. Li, K.; Ni, W.; Tovar, E.; Guizani, M. Joint Flight Cruise Control and Data Collection in UAV-Aided Internet of Things: An Onboard Deep Reinforcement Learning Approach. IEEE Internet Things J. 2021, 8, 9787–9799. [Google Scholar] [CrossRef]
  40. Liu, W.; Wang, R.; Zhang, T.; Li, K.; Li, W.; Ishibuchi, H. Hybridization of evolutionary algorithm and deep reinforcement learning for multi-objective orienteering optimization. IEEE Trans. Evol. Comput. 2022, 27, 1260–1274. [Google Scholar] [CrossRef]
  41. Wang, Z.J.; Zhan, Z.H.; Kwong, S.; Jin, H.; Zhang, J. Adaptive Granularity Learning Distributed Particle Swarm Optimization for Large-Scale Optimization. IEEE Trans. Cybern. 2021, 51, 1175–1188. [Google Scholar] [CrossRef]
  42. Li, K.; Zhang, T.; Wang, R.; Wang, L.; Ishibuchi, H. An Evolutionary Multiobjective Knee-Based Lower Upper Bound Estimation Method for Wind Speed Interval Forecast. IEEE Trans. Evol. Comput. 2022, 26, 1030–1042. [Google Scholar] [CrossRef]
  43. Wang, Z.J.; Zhan, Z.H.; Lin, Y.; Yu, W.J.; Wang, H.; Kwong, S.; Zhang, J. Automatic Niching Differential Evolution With Contour Prediction Approach for Multimodal Optimization Problems. IEEE Trans. Evol. Comput. 2020, 24, 114–128. [Google Scholar] [CrossRef]
  44. Liang, Z.; Dong, H.; Liu, C.; Liang, W.; Zhu, Z. Evolutionary Multitasking for Multiobjective Optimization With Subspace Alignment and Adaptive Differential Evolution. IEEE Trans. Cybern. 2022, 52, 2096–2109. [Google Scholar] [CrossRef]
  45. Wang, Z.J.; Zhan, Z.H.; Lin, Y.; Yu, W.J.; Yuan, H.Q.; Gu, T.L.; Kwong, S.; Zhang, J. Dual-Strategy Differential Evolution With Affinity Propagation Clustering for Multimodal Optimization Problems. IEEE Trans. Evol. Comput. 2018, 22, 894–908. [Google Scholar] [CrossRef]
  46. Cai, X.; Geng, S.; Wu, D.; Cai, J.; Chen, J. A Multicloud-Model-Based Many-Objective Intelligent Algorithm for Efficient Task Scheduling in Internet of Things. IEEE Internet Things J. 2021, 8, 9645–9653. [Google Scholar] [CrossRef]
Figure 1. Telescope m’s observation queue Gantt chart.
Figure 1. Telescope m’s observation queue Gantt chart.
Biomimetics 09 00718 g001
Figure 2. Schematic of observation queue time constraints.
Figure 2. Schematic of observation queue time constraints.
Biomimetics 09 00718 g002
Figure 3. Algorithm’s overall flow chart.
Figure 3. Algorithm’s overall flow chart.
Biomimetics 09 00718 g003
Figure 4. Schematic diagram of problem decomposition.
Figure 4. Schematic diagram of problem decomposition.
Biomimetics 09 00718 g004
Figure 5. Flow chart of proposed decoding method.
Figure 5. Flow chart of proposed decoding method.
Biomimetics 09 00718 g005
Figure 6. Queue head insertion.
Figure 6. Queue head insertion.
Biomimetics 09 00718 g006
Figure 7. Insertion in the middle of the queue.
Figure 7. Insertion in the middle of the queue.
Biomimetics 09 00718 g007
Figure 8. Insertion at the end of the queue.
Figure 8. Insertion at the end of the queue.
Biomimetics 09 00718 g008
Figure 9. Replacement.
Figure 9. Replacement.
Biomimetics 09 00718 g009
Figure 10. AVNS algorithm flowchart.
Figure 10. AVNS algorithm flowchart.
Biomimetics 09 00718 g010
Figure 11. Insertion operator.
Figure 11. Insertion operator.
Biomimetics 09 00718 g011
Figure 12. Commutative operator.
Figure 12. Commutative operator.
Biomimetics 09 00718 g012
Figure 13. Two-opt Operator.
Figure 13. Two-opt Operator.
Biomimetics 09 00718 g013
Table 1. An example of target data.
Table 1. An example of target data.
Target NumberTarget LevelInitial TimeDeadline
n l n O n D n
12376,92577,025
13377,01577,115
16175,83376,121
16176,60977,109
16184,58585,345
291106,755106,921
291112,615113,022
58176,92577,020
Table 2. Sorted target observation data.
Table 2. Sorted target observation data.
Target NumberTarget LevelInitial TimeDeadline
n l n O n D n
2378177,53979,309
6848977,56278,349
2377177,56877,876
4259877,61478,082
4617779,25579,448
6266179,90880,176
2121379,99480,267
Table 3. First set of target observation data.
Table 3. First set of target observation data.
Target NumberTarget LevelInitial TimeDeadline
n l n O n D n
2378177,53978,903
6848977,56278,349
2377177,56877,876
4259877,61478,082
Table 4. Second set of target observation data.
Table 4. Second set of target observation data.
Target NumberTarget LevelInitial TimeDeadline
n l n O n D n
2378178,90379,309
4617779,25579,448
6266179,90880,176
2121379,99480,267
Table 5. The observation plan obtained from the original problem.
Table 5. The observation plan obtained from the original problem.
Target NumberTarget LevelStart TimeEnd Time
n l n S n C n
2378177,53977,629
6848977,56277,652
2377177,65277,742
4259877,74277,832
4617779,25579,345
6266179,90879,998
2121379,99880,088
Table 6. Prioritized data table.
Table 6. Prioritized data table.
Target NumberTarget LevelInitial TimeDeadline
n l n O n D n
58176,92577,020
12376,92577,025
291106,755106,921
291112,615113,022
16175,83376,121
16176,60977,109
16184,58585,345
13377,01577,115
Table 7. Observation queue after insertion.
Table 7. Observation queue after insertion.
Target NumberTarget LevelInitial TimeDeadline
n l n S n C n
16175,83375,923
58176,92577,015
13377,01577,105
291106,755106,845
Table 8. Replaced observation queue.
Table 8. Replaced observation queue.
Target NumberTarget LevelInitial TimeDeadline
n l n S n C n
16175,83375,923
12376,92577,015
13377,01577,105
291106,755106,845
Table 9. Comparison of decoding methods.
Table 9. Comparison of decoding methods.
InstanceTotal Successful Observation Value
Greedy Decoding MethodDecoding Method of This Article
111481366
211141372
311581368
411431383
511611407
611591418
711751405
811571375
912011410
1012001374
Average1161.61387.8
Table 10. Ablation experiment.
Table 10. Ablation experiment.
Instance4 Sub-Problems2 Sub-ProblemsOriginal Problem
TimeValueTimeValueTimeValue
12083.87943028.7801.37663763
21847.3767.92938.6776.47243.5741.4
32020.1765.53105.3767.76578.4732.6
41976.6774.53218781.17028.5740.4
52044786.93443.6798.36091.2751
62085.2768.13394774.36690.5735.3
71943.1781.83352.2787.36575.6745.9
81890.1757.13287.1766.36766.6723.3
91890.17523200.6765.46873.6720.4
101995.9727.53042.37366352.3697.4
Average1977.62767.53201.04775.46786.32735.1
Table 11. Experimental results.
Table 11. Experimental results.
InstanceTarget NumberMaximumAverageStd. DeviationAvg. Count
1117517101700.64.781130.3
21161168616803.461122
3117116881680.33.881125.3
4116916931685.64.101129.9
51195174917453.701151
6119617561744.43.891150.3
7119617491752.42.501150.4
81191174017372.271155.9
9118917561752.43.961145.3
10115116861680.62.561103.6
Table 12. Comparative experiment.
Table 12. Comparative experiment.
InstanceTarget NumberTimeAVNSGA
MaximumAverageStd. DeviationAvg. CountMaximumAverageStd. DeviationAvg. Count
15922118.581479411.80536.5765762.76.24518.6
25731905788767.98.12525.625736731.45.71496.4
35782213.3778765.58.65522.75730728.34.74500
45901999.7778774.53.00539737733.79.21504.4
56012021.3791786.93.02538.375750747.15.29499.9
65921971.4770768.191.96538.625733729.95.83502
75912000.6784781.82.54539.5744742.45.39502.9
85812013.7761757.12.42558.1257217198.70497.1
95822068.97587523.32528.5720716.35.14494
105692007.6730727.52.55519.5698691.95.82489.4
Average584.92032775.2767.5394.74534.65733.4730.36.21500.5
Table 13. Comparative experiment.
Table 13. Comparative experiment.
InstanceTarget NumberTimeAVNSSA
MaximumAverageStd. DeviationAvg. CountMaximumAverageStd. DeviationAvg. Count
15922118.581479411.80536.5779771.85.42529.5
25731905788767.98.12525.625746733.36.28504.3
35782213.3778765.58.65522.75748740.46.92513.8
45901999.7778774.53.00539755744.65.28518.4
56012021.3791786.93.02538.375771763.45.72519.2
65921971.4770768.191.96538.625753729.96.72516.2
75912000.6784781.82.54539.57447404.76517
85812013.7761757.12.42558.125758752.85.27509.2
95822068.97587523.32528.57417307.47505
105692007.6730727.52.55519.5712706.33.69503.3
Average584.92032775.2767.5394.74534.65750.5740.975.75513.59
Table 14. 10 Telescope experiment results.
Table 14. 10 Telescope experiment results.
InstanceTarget NumberMaximumAverageAvg. CountAvg. Time
1592835835592145.4
2573806806573189.6
3578807807578128.3
4590819819590199.7
5601848848601221
6592819819592204.2
7591830830591156.6
8581804804581148
9582805805582209.4
10569774774569180.6
Average584.9814.7814.7584.9178.3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, K.; Ye, B.-L.; Xia, X.; Wang, Z.; Zhang, X.; Jiang, H. A Space Telescope Scheduling Approach Combining Observation Priority Coding with Problem Decomposition Strategies. Biomimetics 2024, 9, 718. https://doi.org/10.3390/biomimetics9120718

AMA Style

Zhang K, Ye B-L, Xia X, Wang Z, Zhang X, Jiang H. A Space Telescope Scheduling Approach Combining Observation Priority Coding with Problem Decomposition Strategies. Biomimetics. 2024; 9(12):718. https://doi.org/10.3390/biomimetics9120718

Chicago/Turabian Style

Zhang, Kaiyuan, Bao-Lin Ye, Xiaoyun Xia, Zijia Wang, Xianchao Zhang, and Hai Jiang. 2024. "A Space Telescope Scheduling Approach Combining Observation Priority Coding with Problem Decomposition Strategies" Biomimetics 9, no. 12: 718. https://doi.org/10.3390/biomimetics9120718

APA Style

Zhang, K., Ye, B. -L., Xia, X., Wang, Z., Zhang, X., & Jiang, H. (2024). A Space Telescope Scheduling Approach Combining Observation Priority Coding with Problem Decomposition Strategies. Biomimetics, 9(12), 718. https://doi.org/10.3390/biomimetics9120718

Article Metrics

Back to TopTop