Multi-Adaptive Strategies-Based Higher-Order Quantum Genetic Algorithm for Agile Remote Sensing Satellite Scheduling Problem

The agile remote sensing satellite scheduling problem (ARSSSP) for large-scale tasks needs to simultaneously address the difficulties of complex constraints and a huge solution space. Taking inspiration from the quantum genetic algorithm (QGA), a multi-adaptive strategies-based higher-order quantum genetic algorithm (MAS-HOQGA) is proposed for solving the agile remote sensing satellites scheduling problem in this paper. In order to adapt to the requirements of engineering applications, this study combines the total task number and the total task priority as the optimization goal of the scheduling scheme. Firstly, we comprehensively considered the time-dependent characteristics of agile remote sensing satellites, attitude maneuverability, energy balance, and data storage constraints and established a satellite scheduling model that integrates multiple constraints. Then, quantum register operators, adaptive evolution operations, and adaptive mutation transfer operations were introduced to ensure global optimization while reducing time consumption. Finally, this paper demonstrated, through computational experiments, that the MAS-HOQGA exhibits high computational efficiency and excellent global optimization ability in the scheduling process of agile remote sensing satellites for large-scale tasks, while effectively avoiding the problem that the traditional QGA has, namely low solution efficiency and the tendency to easily fall into local optima. This method can be considered for application to the engineering practice of agile remote sensing satellite scheduling for large-scale tasks.


Introduction
Remote sensing satellites are particularly effective tools for observing and analyzing earth's resources and environment and play an irreplaceable role in earth resource exploration, natural disaster prevention, and environmental protection [1].Compared to traditional non-agile remote sensing satellites, as shown in Figure 1, the advantage of agile remote sensing satellites lies in their flexible and fast maneuverability, which not only greatly enhances the observation capabilities of satellites but also leads to a significant increase in the complexity of solving satellite execution plans [2].Before the rapid development of agile remote sensing satellites, many scholars had conducted research over a long period on the scheduling problem of non-agile remote sensing satellites.Previous research mainly focused on explaining the concept of satellite scheduling.Due to the lack of large-scale task sets, satellite scheduling models were often simplified [3][4][5][6].Abramson [7] considered satellite constraints, including satellite lateral sway and storage and introducing an integer linear programming method for finding a solution.Bianchesi [8] takes the COSMO SkyMed constellation as the research object, with the goal of maximizing the number of images obtained, and adopts a heuristic construction algorithm to solve the satellite scheduling problem, but does not introduce enough real constraints.Florio [9] considered the mobility constraints of satellites in his research, while also considering the energy and storage capabilities of satellites as constraints.Lemaître [10] took into account the strong maneuverability of agile remote sensing satellites and considered the satellite's maneuverability model in detail, but simplified the processing by incorporating additional constraints for the satellite.Reference [11] treated the decomposition of regional targets as a type of constraint for satellites.Reference [12] proposed and to some extent solved the resource matching problem based on conventional observation tasks and emergency observation tasks, with task efficiency as the global optimization objective in the scheduling process.Peng's research also focuses on this aspect [13].Xu [14] focused on satellite observation time window constraints and energy constraints, simplifying other constraints.Reference [15] focused on addressing the observation time dependency in the scheduling process of agile remote sensing satellites in his research, but simplified other constraints.The algorithm proposed by Reference [16] can effectively improve the real-time performance of computation, but its feasibility for largescale tasks has not been verified.Du's research focuses on multi-star mission scheduling, and his proposed method is beneficial for the efficiency of large-scale task allocation, but does not consider the single-star constraint problem in detail [17].
Starting from the end of the twentieth century, intelligent heuristic optimization algorithms have been widely used in the field of satellite mission planning and scheduling due to their good optimization capabilities, such as the ant colony algorithm [18,19], the particle swarm optimization algorithm [20], the genetic algorithm [21,22], the improved genetic algorithm [23][24][25][26][27], the simulated annealing algorithm [28][29][30], the taboo search algorithm [31][32][33][34][35], and evolutionary algorithms [36,37].The complexity of agile remote sensing satellite mission scheduling based on large-scale tasks has increased Before the rapid development of agile remote sensing satellites, many scholars had conducted research over a long period on the scheduling problem of non-agile remote sensing satellites.Previous research mainly focused on explaining the concept of satellite scheduling.Due to the lack of large-scale task sets, satellite scheduling models were often simplified [3][4][5][6].Abramson [7] considered satellite constraints, including satellite lateral sway and storage and introducing an integer linear programming method for finding a solution.Bianchesi [8] takes the COSMO SkyMed constellation as the research object, with the goal of maximizing the number of images obtained, and adopts a heuristic construction algorithm to solve the satellite scheduling problem, but does not introduce enough real constraints.Florio [9] considered the mobility constraints of satellites in his research, while also considering the energy and storage capabilities of satellites as constraints.Lemaître [10] took into account the strong maneuverability of agile remote sensing satellites and considered the satellite's maneuverability model in detail, but simplified the processing by incorporating additional constraints for the satellite.Reference [11] treated the decomposition of regional targets as a type of constraint for satellites.Reference [12] proposed and to some extent solved the resource matching problem based on conventional observation tasks and emergency observation tasks, with task efficiency as the global optimization objective in the scheduling process.Peng's research also focuses on this aspect [13].Xu [14] focused on satellite observation time window constraints and energy constraints, simplifying other constraints.Reference [15] focused on addressing the observation time dependency in the scheduling process of agile remote sensing satellites in his research, but simplified other constraints.The algorithm proposed by Reference [16] can effectively improve the real-time performance of computation, but its feasibility for large-scale tasks has not been verified.Du's research focuses on multi-star mission scheduling, and his proposed method is beneficial for the efficiency of large-scale task allocation, but does not consider the single-star constraint problem in detail [17].
Starting from the end of the twentieth century, intelligent heuristic optimization algorithms have been widely used in the field of satellite mission planning and scheduling due to their good optimization capabilities, such as the ant colony algorithm [18,19], the particle swarm optimization algorithm [20], the genetic algorithm [21,22], the improved genetic algorithm [23][24][25][26][27], the simulated annealing algorithm [28][29][30], the taboo search algorithm [31][32][33][34][35], and evolutionary algorithms [36,37].The complexity of agile remote sensing satellite mission scheduling based on large-scale tasks has increased exponentially, and intelligent heuristic optimization algorithms have shown the disadvantage of low solution efficiency during the solution process.Therefore, traditional algorithms are in urgent need of further innovation and development.In recent years, artificial intelligence has been applied in various fields.Many scholars have combined artificial intelligence algorithms with satellite scheduling problems and conducted a series of studies [38][39][40][41][42][43].In addition, many scholars have applied physics-informed deep learning approaches to many fields and proved the algorithm's superior optimization ability [44].Similarly, it is feasible to apply this method to the field of satellite mission scheduling.The biggest feature of this type of method is that it can improve the timeliness of satellite execution plans while achieving global optimization.However, these methods are limited by the performance of on-board computers and have certain limitations in the current field of engineering applications.
Currently, in the engineering practice process of satellite mission scheduling, how to efficiently and reliably generate satellite task execution plans is a key issue that urgently needs to be solved.Narayanan [45] integrated the ideas of genetic algorithms and quantum theory and demonstrated that this algorithm has a significant effect on improving search capability and computational speed.Silveira [46] innovatively proposed a quantum-inspired evolutionary algorithm and demonstrated its effectiveness through sorting problems.Nowotniak [47] introduced the concept of high-order on the basis of a quantum genetic algorithm and proved that this method has significant effects on improving computational speed.
In our study, we propose a novel MAS-HOQGA to solve the agile remote sensing satellite scheduling problem.First, taking into account the constraints, including the time-dependent characteristics of agile remote sensing satellites, attitude maneuverability, payload observation capability, data transmission resources, storage capability, and energy balance, a refined satellite constraint model with comprehensive multiple constraints was established.The scheduling objective function was established based on the comprehensive value derived from the number of tasks and their priorities.Then, on the basis of the traditional QGA, the quantum register operator, the adaptive evolution strategy, and the adaptive mutation transfer strategy are introduced.The proposed MAS-HOQGA is used to solve the agile satellite scheduling problem.
Compared to traditional QGA, the quantum chromosome in this paper is composed of quantum register operators, ensuring that the proposed algorithm has advantages in both individual measurement and updates.Quantum registers can reduce the time complexity of individual measurements and updates, thereby improving algorithm performance.In addition, the adaptive evolution strategy and adaptive mutation transfer strategy introduced in this paper are beneficial for improving the convergence speed of the algorithm and preventing it from falling into local optima.
The experimental results show that the MAS-HOQGA has achieved significant improvements in both comprehensive revenue and algorithm running time for scheduling results compared with the QGA. Figure 2 shows the framework of the research content of this paper.
Specifically, the specific work and innovation of the paper are as follows: (1) Considered the agile remote sensing satellite mission scheduling problem for large-scale tasks scenarios.
(2) Established a comprehensive multi-constraint satellite scheduling model and described the optimization objective function based on the total task number and the total task priorities.
(3) A higher-order QGA based on multi-adaptive strategies was proposed.
(4) The proposed algorithm shows excellent computing power and global optimization capabilities in the field of agile remote sensing satellite mission scheduling.egy introduced in this paper are beneficial for improving the convergence speed of the algorithm and preventing it from falling into local optima.
The experimental results show that the MAS-HOQGA has achieved significant improvements in both comprehensive revenue and algorithm running time for scheduling results compared with the QGA. Figure 2 shows the framework of the research content of this paper.This paper is divided into five parts.After the introduction, Section 2 introduces the establishment process of the multi-constraint refined satellite scheduling model.Section 3 introduces MAS-HOQGA, including quantum individual encoding and quantum register initialization, the measurement and adaptive evolution strategies of quantum register operators, and the adaptive mutation transfer strategy of quantum registers.In Section 4, the comprehensive revenue and algorithm running time of the MAS-HOQGA and QGA for scheduling tasks of different scales are compared, and the impact of the probability amplitude adjustment parameter setting on the results is analyzed.The final section presents the conclusions of this paper.

Refined Satellite Scheduling Model with Multi-Constraints
The establishment of a satellite mission scheduling model needs to consider the usage constraints of the satellite and the requirements of its various subsystems.In this section, the relevant parameters in the satellite scheduling model are defined, and the constraints of each subsystem in the satellite scheduling model are given.

Parameter Definition
The satellite properties are defined as follows: M max : satellite storage maximum; m ut : satellite payload data rate; E max : satellite energy maximum; e p ut : the energy consumption per unit of time during satellite payload operation; e d ut : the energy consumption per unit of time during satellite maneuvering; A ut : the angle at which a satellite can maneuver its orientation per unit of time; Fov: the satellite payload field of view angle.
The task properties are defined as follows: X = X 1 , X 2 , . . ., X i , . . ., X N : the collection of tasks, where N is the total number included in X, and X i is the i-th observation task.
X i = id i , prio i , lat i , lon i , alt i , dur i , obs i,bgn , obs i,end : id i is the task ID, prio i is the task priority, lat i is the task latitude, lon i is the task longitude, alt i is the task altitude, dur i is the task observation duration, obs i,bgn is the task observation began time, and obs i,end is the task observation end time.
The task window properties are defined as follows: the observation windows set of task X i , m is the total number of observation windows corresponding to the mission X i , where win m,i = Orb m.i , win Bgn m,i , win End m,i , Orb m.i is the number of cycles in which the observation window is located, win Bgn m,i is the beginning time of the window, and win End m,i is the end time of the window.

Payload Constraints
Based on the working capacity of the payload, specify the maximum working time within one orbit of the payload and require that the working time of each orbit of the payload be less than the maximum working time.where k is the number of observation tasks arranged within the j-th orbital cycle, Payload orb j dur i is the observation duration of the i-th payload mission arranged within the orbital cycle j, and Payload_T orb max is the maximum working duration of the payload within one orbital cycle.

Maneuverability Constraints
The time interval between two observation tasks needs to be greater than the fastest maneuvering time of the satellite, as shown in Figure 3.

Payload Constraints
Based on the working capacity of the payload, specify the maximum working time within one orbit of the payload and require that the working time of each orbit of the payload be less than the maximum working time.
where k is the number of observation tasks arranged within the j-th orbital cycle,      is the observation duration of the i-th payload mission arranged within the orbital cycle j, and _   is the maximum working duration of the payload within one orbital cycle.

Maneuverability Constraints
The time interval between two observation tasks needs to be greater than the fastest maneuvering time of the satellite, as shown in Figure 3.

Energy Constraints
During the operation of the satellite, it is necessary to always ensure the satellite energy balance; that is, the energy consumed by the satellite during operation should be less than the maximum energy consumption.
N is the count of successfully arranged observation tasks, N+1 is the number of attitude maneuvers arranged,  0 is the initial attitude of the satellite, and  +1 is the attitude that the satellite needs to maintain after completing all planned tasks.

Energy Constraints
During the operation of the satellite, it is necessary to always ensure the satellite energy balance; that is, the energy consumed by the satellite during operation should be less than the maximum energy consumption.
N is the count of successfully arranged observation tasks, N + 1 is the number of attitude maneuvers arranged, A 0 is the initial attitude of the satellite, and A N+1 is the attitude that the satellite needs to maintain after completing all planned tasks.

Data Storage Constraints
According to the design requirements for satellite data storage, it is required that the data generated by the payload be promptly arranged to the ground, and they cannot exceed the satellite storage maximum.
2.2.5.Comprehensive Revenue Scheduling Function In the actual engineering scheduling process, satellite users hope that the satellite can observe more tasks on one hand, and on the other hand, they also hope that the satellite can prioritize completing higher priority tasks.This paper considers both the number of observation tasks and their corresponding priorities.A comprehensive scheduling revenue function is established as the objective function for the scheduling solution.
where ω 1 represents the revenue weight corresponding to the number of observation missions, ω 2 represents the revenue weight corresponding to the priority of the observation missions, and y i serves as a decision variable.

Quantum Chromosome Encoding and Register Initialization
In quantum computing, quantum genes are composed of qubits, with | 0⟩ and | 1⟩ being the basic elements that make up a qubit.The qubits | q⟩ are represented as follows: In the equation above, the probability amplitudes of | 0⟩ and | 1⟩ being α and β, respec- tively, the probability of obtaining 0 for this gene locus after measurement is |α| 2 , and the probability of obtaining 1 for this gene locus after measurement is |β| 2 .For a scheduling scheme consisting of N observation tasks, the encoding of quantum chromosome q [X] can be represented as In the above-mentioned quantum chromosome encoding, each qubit of the quantum chromosome is represented by a set of independent binary genes, and there is no connection between qubits.Such a quantum chromosome is called a quantum chromosome of order-1.
The difference between higher-order quantum chromosomes and quantum chromosomes of order-1 is that the qubits of higher-order quantum chromosomes will affect each other, and the scope of influence depends on the order of the higher-order quantum chromosomes.
Taking quantum chromosome of order-2 as an example, two consecutive qubits constitute an independent quantum register, and each quantum register is a superposition state composed of any linear combination of four basic states | 00⟩ , | 01⟩ , | 10⟩, and | 11⟩ .
In the above equation, the probability amplitudes of | 00⟩ , | 01⟩ , | 10⟩, and | 11⟩ are α 0 , α 1 , α 2 , and α 3 , respectively.The probability of obtaining | 00⟩ for this register after measurement is |α 0 | 2 , the probability of obtaining | 01⟩ for this register after measurement is |α 1 | 2 , the probability of obtaining | 10⟩ for this register after measurement is |α 2 | 2 , the probability of obtaining | 11⟩ for this register after measurement is |α 3 | 2 .The probability amplitudes meet the normalization requirements and are all initialized to 1  2 .This operation can ensure that the initial state of each individual in the initial population is consistent, thereby ensuring that each individual in the initial solution has the same probability of being selected.
For a scheduling scheme consisting of N observation tasks, the encoding of a quantum chromosome of order-2 R [X] can be represented as follows: The schematic diagram of a quantum chromosome composed of independent qubits and quantum registers is shown in Figure 4.
| 1 | 2 , the probability of obtaining |10⟩ for this register after measurement is | probability of obtaining |11⟩ for this register after measurement is | 3 | 2 .
The probability amplitudes meet the normalization requirements and are all ized to 1 2 .This operation can ensure that the initial state of each individual in the population is consistent, thereby ensuring that each individual in the initial solut the same probability of being selected.
For a scheduling scheme consisting of N observation tasks, the encoding of a tum chromosome of order-2 | [] ⟩ can be represented as follows: The schematic diagram of a quantum chromosome composed of independent and quantum registers is shown in Figure 4.

Quantum Registers Measurement
Quantum chromosomes require a step involving a measurement operator to c each quantum register from a superposition state into a definite state.For a quantum mosome of order-2, it is necessary to collapse each quantum register composed consecutive qubits, from a superposition state composed of any linear combination four basic states, into a definite measurement value.The measurement function the measurement value, while recording the index corresponding to the basic st

q q
Rgt 1 =R 1 3 4 q q chromosome containing N qubits in order-2

Quantum Registers Measurement
Quantum chromosomes require a step involving a measurement operator to collapse each quantum register from a superposition state into a definite state.For a quantum chromosome of order-2, it is necessary to collapse each quantum register composed of two consecutive qubits, from a superposition state composed of any linear combination of the four basic states, into a definite measurement value.The measurement function returns the measurement value, while recording the index corresponding to the basic states of | 00⟩ , | 01⟩ , | 10⟩ , and | 11⟩ .Algorithm 1 is the quantum registers measurement strategy pseudocode.

Adaptive Evolution for Quantum Registers
In quantum genetic algorithms, quantum gate rotation is a way for each generation of individuals to evolve.The basic principle for quantum gate rotation is to ensure that the current quantum individual evolves towards the optimal quantum individual.The evolution of higher-order quantum chromosomes also ensures that each quantum individual evolves toward the optimal quantum individual.The difference lies in that the evolution method of quantum registers is not the same as that of traditional qubits.specifically, the basic evolution principle of quantum registers is to reduce the probability amplitude of the basic state that does not corresponding to the probability amplitude index.
Let quantum register j (j ∈ (1, 2, . . . ,N/2)) of the quantum individual i (i ∈ (1, 2, . . . ,Q)) currently be represented as R i,j .Figure 5 shows the schematic diagram of the quantum register evolution operator.R b j is the best individual from the previous generation.The length of the vertical line segment represents the probability of obtaining the corresponding state of the quantum register after measurement of R i,j .According to the logic of quantum register evolution, for the four states of the j-th register of the current individual R i,j , which needs to evolve, except for the probability amplitude corresponding to α 1 i,j , which needs to be increased, the probability amplitude of the other three states needs to be reduced, where µ ∈ (0,1) is the probability amplitude adjustment parameter.probability amplitude index ← 3 16: end if

Adaptive Evolution for Quantum Registers
In quantum genetic algorithms, quantum gate rotation is a way for each generatio of individuals to evolve.The basic principle for quantum gate rotation is to ensure th the current quantum individual evolves towards the optimal quantum individual.Th evolution of higher-order quantum chromosomes also ensures that each quantum ind vidual evolves toward the optimal quantum individual.The difference lies in that the ev lution method of quantum registers is not the same as that of traditional qubits.speci cally, the basic evolution principle of quantum registers is to reduce the probability am plitude of the basic state that does not corresponding to the probability amplitude inde Let quantum register j ( ∈ (1,2, … , /2)) of the quantum individual i ( ∈ (1,2, … ,  currently be represented as | , ⟩. Figure 5 shows the schematic diagram of the quantu register evolution operator.|   ⟩ is the best individual from the previous generation.Th length of the vertical line segment represents the probability of obtaining the correspon ing state of the quantum register after measurement of | , ⟩.According to the logic quantum register evolution, for the four states of the j-th register of the current individu | , ⟩, which needs to evolve, except for the probability amplitude corresponding to  which needs to be increased, the probability amplitude of the other three states needs be reduced, where  ∈ (0,1) is the probability amplitude adjustment parameter.In the progress of quantum evolution, the probability amplitude adjustment param eter will influence the converge result.When the value of  is too large, premature co vergence may occur, while when the value of  is too small, the convergence rate is to slow.Accelerating the speed of quantum individual evolution can effectively improv computational efficiency.This paper sets the probability amplitude adjustment paramet to a dynamic adaptive value that changes with the evolution state.The probability amp tude adjustment parameter is determined by the fitness of the currently quantum regist to be evolved, the fitness of the optimal quantum register.
Algorithm 2 presents the pseudocode for the adaptive evolution strategy of quantu registers.In the progress of quantum evolution, the probability amplitude adjustment parameter will influence the converge result.When the value of µ is too large, premature convergence may occur, while when the value of µ is too small, the convergence rate is too slow.Accelerating the speed of quantum individual evolution can effectively improve computational This paper sets the probability amplitude adjustment parameter to a dynamic adaptive value that changes with the evolution state.The probability amplitude adjustment parameter is determined by the fitness of the currently quantum register to be evolved, the fitness of the optimal quantum register.

µ =
Fitness opt − Fitness cur Fitness opt (12) Algorithm 2 presents the pseudocode for the adaptive evolution strategy of quantum registers.

Adaptive Mutation Transfer in Quantum Registers
Quantum mutation operations are usually not considered in traditional quantum genetic algorithms because, compared to genetic algorithms, QGA contains richer encoding information, ensuring population diversity.However, the evolution principle states that each quantum individual evolves toward the optimal quantum individual, which may lead to algorithm premature convergence and falling into local optima.This paper innovatively proposes an adaptive quantum mutation transfer strategy which increases gene diversity and improves search ability through this operation, preventing the process from falling into local optima.
The adaptive mutation transfer strategy of quantum registers requires the introduction of mutation transfer probability parameters, can affect gene diversity and algorithm efficiency.Algorithm 3 presents the pseudocode for the adaptive mutation transfer strategy applied to quantum registers.

Computational Experiment
This section first introduces the experimental environment, scenario parameters, and algorithm parameter settings.Then, the experimental results are presented and some analyses are given.

Experimental Environment
The environment for the computational experiment was implemented on an Intel Core I7 CPU (2.80 GHz) with 16 GB of RAM, running the Win10 operating system on a 64-bit architecture.MATLAB R2021b was used for coding.

Computational Experiment
This section first introduces the experimental environment, scenario parameters, and algorithm parameter settings.Then, the experimental results are presented and some analyses are given.

Experimental Environment
The environment for the computational experiment was implemented on an Intel Core I7 CPU (2.80 GHz) with 16 GB of RAM, running the Win10 operating system on a 64-bit architecture.MATLAB R2021b was used for coding.

Parameter Settings
Table 1 shows the orbit parameter information.The satellite attributes are as follows: the payload field of the view angle is set to 45 • , the minimum solar altitude angle is set to 5 • , the maximum maneuvering angular velocity is set to 1.3 • /s, the maximum maneuvering angular acceleration is set to 0.01 • /s 2 .The satellite attributes mentioned above are typical parameters in realistic satellite operations.
The order of the quantum chromosome was set to 2. The satellite attributes of individuals' population scale number Q was set to 50.The satellite attributes of the number of iterative evolutions were set to 100.The mutation transfer probability parameter P m = 0.05.The revenue weight corresponding to the number of observation missions ω 1 = 0.35.The revenue weight corresponding to the priority of observation missions ω 2 = 0.65.
Target information includes the following: within 80 • north and south latitude, the number of 226, 426, and 626 observation targets are randomly generated.The distribution of observation targets is shown in Figure 7, where the red line represents the satellite trajectory.This paper uses the Monte Carlo simulation method for the complete experimental verification across tasks of different scales, and the results are the average values of the algorithm after running 10 times.In addition to comprehensive revenue, the algorithm's running time is also used as an evaluation indicator for the proposed algorithm.In this paper, the best revenue of each generation will gradually increase and stabilize with the iteration of the algorithm until the revenue no longer changes.When the best revenue no longer changes, it means that the evolution has converged.The shorter the algorithm's running time to reach the convergence state for the best revenue, the faster the algorithm converges.
By comparing Figure 8 and Table 2, it can be seen that under the same number of iterations, the MAS-HOQGA can achieve higher revenue, and the convergence speed of the MAS-HOQGA is faster, with fewer iterations corresponding to achieving this optimal revenue.When the task scale is 226, compared with the QGA, the MAS-HOQGA has a 6.44% increase in comprehensive revenue and a 46.9% decrease in algorithm running time.When the task scale is 426, compared with the QGA, the MAS-HOQGA has a comprehensive revenue increase of 10.23% and an algorithm running time decrease of 24.5%.When the task scale is 626, the MAS-HOQGA has a 7.21% increase in comprehensive revenue and a 32.8% decrease in algorithm running time compared to QGA.By comparing Figure 8 and Table 2, it can be seen that under the same number of iterations, the MAS-HOQGA can achieve higher revenue, and the convergence speed of the MAS-HOQGA is faster, with fewer iterations corresponding to achieving this optimal revenue.When the task scale is 226, compared with the QGA, the MAS-HOQGA has a 6.44% increase in comprehensive revenue and a 46.9% decrease in algorithm running time.When the task scale is 426, compared with the QGA, the MAS-HOQGA has a comprehensive revenue increase of 10.23% and an algorithm running time decrease of 24.5%.When the task scale is 626, the MAS-HOQGA has a 7.21% increase in comprehensive revenue and a 32.8% decrease in algorithm running time compared to QGA.The impact of probability amplitude adjustment parameter on the MAS-HOQGA The iterations number was set to 300.The experimental results with different probability amplitude adjustment parameter settings are statistically analyzed and shown in Table 3.Based on the maximum revenue in the experimental results and the iteration number corresponding to the maximum revenue, it can be concluded that, within a certain range, setting a large probability amplitude adjustment parameter can achieve convergence faster.However, if it is too large, it can easily lead to the algorithm falling into premature convergence, causing local optima.If the parameter is set too small, it can lead to a slow convergence speed.The probability amplitude adjustment parameter adaptive strategy proposed in this paper can reduce time consumption and avoid falling into local optima.

Computational complexity analysis of MAS-HOQGA
The calculation speed of the scheduling algorithm is an important factor in evaluating the satellite mission scheduling problems.The faster the algorithm's calculation speed, the higher the timeliness of the mission scheduling becomes.The speed advantage of quantum computing comes from quantum superposition and quantum entanglement.In this paper, by introducing quantum register operators, the characteristics of quantum superposition and quantum entanglement are utilized to improve the speed of quantum individual measurement and update, effectively reducing the computational complexity of MAS-HOQGA.In this paper, the order of the MAS-HOQGA is set to 2, and theoretically, the time complexity of quantum individual measurement and update can be increased by 50% compared to the traditional QGA.At the same time, the encoding complexity of quantum individuals has not increased, so the running time of the MAS-HOQGA has been greatly improved.However, as the order increases, the encoding complexity of quantum individuals will greatly increase, which will also lead to an increase in the computational complexity of the algorithm.This problem needs to be solved by improving computer performance.

3.
The Impact of Task Scale and Constraints on MAS-HOQGA Through theoretical analysis, it can be determined that the MAS-HOQGA has good scalability in terms of task scale and scheduling model constraints.Firstly, the task scale determines the length of the quantum chromosome, which will only affect the computational complexity of the MAS-HOQGA, but it will not affect the solution methodology proposed in this paper.Secondly, increasing the constraints of the scheduling model will affect both the running time of the scheduling algorithm and the number of executable tasks in the final solution, but it will not affect the solution methodology of the MAS-HOQGA.

Figure 1 .
Figure 1.Imaging schematic diagram of agile remote sensing satellite.

Figure 1 .
Figure 1.Imaging schematic diagram of agile remote sensing satellite.

Figure 2 .
Figure 2. Framework diagram of the paper.

Figure 3 .
Figure 3. Schematic diagram showing constraints for agile earth observation satellite maneuverability.

Figure 3 .
Figure 3. Schematic diagram showing constraints for agile earth observation satellite maneuverability.

Figure 4 .
Figure 4.Quantum chromosomes composed of independent qubits and quantum registers.

Figure 4 .
Figure 4.Quantum chromosomes composed of independent qubits and quantum registers.

Figure 6 Figure 6 .
Figure6shows the implementation flow chart for the MAS-HOQGA.

Table 2 .
The experimental results obtained by QGA and MAS-HOQGA.

Table 2 .
The experimental results obtained by QGA and MAS-HOQGA.

Table 3 .
The experimental results for different probability amplitude adjustment parameters.