# A Knowledge-Based Simulated Annealing Algorithm to Multiple Satellites Mission Planning Problems

Authors: Da-Wei Jin; Li-Ning Xing
Journal: Journal of Applied Mathematics
 (2013)

---

## Abstract

The multiple satellites mission planning is a complex combination optimization problem. A knowledge-based simulated annealing algorithm is proposed to the multiple satellites mission planning problems. The experimental results suggest that the proposed algorithm is effective to the given problem. The knowledge-based simulated annealing method will provide a useful reference for the improvement of existing optimization approaches.


---

## Body

## 1. Introduction

The imaging satellite plays a significant role in various fields such as the disaster prevention and environmental protection and attracts great attention of many scholars in the world [1, 2]. Although the number of satellites in orbit is increasing, the imaging satellite resources are still limited compared with the rapid growing imaging requirements [3, 4]. To improve the performance of satellite resources, the appropriate mathematical models and software tools must be employed to realize a better management and allocation of satellite resources [5, 6].With the development of aerospace industry, the number and classification of imaging satellites are increasing. Also, both the number and the category of imaging requirements are largely increased [7]. It is very necessary to make a comprehensive scheduling of multiple imaging satellites [8, 9]. Therefore, this paper proposes a novel method to deal with the multiple satellites mission planning problem [10, 11].The motivation of this work can be summarized as follows. At first, the multiple satellites mission planning is a classical combination optimization problem and it is very urgent to develop some new methods to deal with this problem. At second, the interaction between evolution and learning is more popular than before, so a knowledge-based simulated annealing algorithm is designed and implemented according to this fashion.


## 2. Multiple Satellites Mission Planning Problems

The multiple satellites mission planning problem can be summarized as follows. Under the conditions of satisfying the resource constraints and imaging requirements, make a comprehensive scheduling of spot targets and regional targets, allocate satellite resources and working time for each target, and develop an optimized observation plan to maximize the summation of priority of completed targets. There are two subproblems to be solved in this problem. One is how to make a comprehensive scheduling of spot targets and regional targets, based on which the other one is how to optimize the observation plan of satellite resources. The multiple satellites mission planning is a complex combination optimization problem.


### 2.1. Hypothesis

(1)
There is only one sensor in each satellite. In fact, if there are multiple sensors in one satellite, then this satellite can be divided into multiple satellites with one sensor. For this reason, both the sensorsi and the satellite si have the same meaning.
(2)
The regional target can be divided into multiple atomic tasks. In fact, a spot target just is one atomic task. In this paper, the atomic taskoijkv means the vth subtask of task tj which can be observed in the kth time window of satellite si. The regional target tj is divided into the following sets:
(1)Oj=⋃i=1NS⋃k=1Nij⋃v=1Nijkoijkv,
whereNS denotes the number of satellites, Nij denotes the number of visible time windows about sensor si and target tj, and Nijk denotes the number of atomic tasks which can be observed in the kth time window of satellite si.
(3)
In order to improve the performance of imaging resources, multiple atomic tasks are integrated into one composite task, which can be observed at one time.


### 2.2. Inputs

(1)
The number of satellites isNS. Spani denotes the longest working time of satellite si within one single booting.
(2)
M
i and Pi denote the maximum of storage capacity and energy of satellite si,mi, and pi refer to the storage and energy consumption of satellite si for observation in unit time.
(3)
R
i denotes the maximum number of side-view imaging in each circle, ϑi is the swaying speed of satellite si (degree per second), ρi refers to the energy consumption of swaying in each degree, di is the settling time of swaying of satellite si.
(4)
The number of tasks isNT. Suppose the first h targets are spot ones, and the last NT-h targets are regional ones. LPi means the priority of target i.
(5)
V
W
i
j
=
{
v
w
i
j
1
,
v
w
i
j
2
,
…
,
v
w
i
j
N
i
j
} denotes the set of time windows about sensor si and target tj,and Nij is the total number of visible time windows about sensor si and target tj.
(6)
w
s
i
j
k
v
,
w
e
i
j
k
v, and gijkv refer to the starting time, ending time, and observation angle of atomic task oijkv, respectively.
(7)
b
i
l denotes the lth composite task of satellite si, while wsil,weil, and gil refer to the starting time, ending time, and observation angle of composite task bil.
(8)
B
i
=
{
b
i
1
,
b
i
2
,
…
,
b
i
N
B
i
} denotes the set of composite tasks of satellite si, in which NBi is the number of composite tasks.


### 2.3. Outputs

In the multiple satellites mission planning problem, the outputs (decision variables) can be summarized as follows:(2)xijkvl={1,Theatomictaskoijkvwasexecutedinthelthcompositetasksofsatellitesi0,Otherwise.


### 2.4. Objectives

The objective of this problem is to maximize the priority summation of fulfilled tasks. That is,(3)maxJ=CSpot+CPolygon,
where CSpot denotes the priority summation of spot targets, and it is
(4)CSpot=∑j=1hCSpot(tj)=∑j=1h∑i=1NS∑k=1Nij∑l=1NBjxijk1l×LPj.CPolygon means the priority summation of regional targets. This computation process is a little complex. At first, we should compute the overlap area between the multiple observed atomic tasks and regional target tj. That is,
(5)SPolygon=(⋃i=1NS⋃k=1Nij⋃v=1Nijk⋃l=1NBixijkvlψ(oijkv))⋂ψ(tj).
Here, ψ(oijkv) denotes the polygon area of atomic task oijkv, ∩ denotes the intersection operation of multiple sets, ψ(tj) denotes the polygon area of regional target tj, and ∪ denotes the union operation of different sets. At second, we should compute the coverage rate of each regional target. It is
(6)Cover(tj)=f(SPolygon)f(tj),
where f(SPolygon) denotes the proportion of polygon area SPolygon, and f(tj) denotes the proportion of regional target tj. In total, the priority summation of regional targets can be computed as follows:
(7)CPolygon=∑j=h+1NT[Cover(tj)×LPj].


### 2.5. Constraints

(1)
Uniqueness constraint on spot targets; that is, each spot target can only be observed once:(8)∑i=1NS∑k=1Nij∑l=1NBjxijk1l≤1,∀j=1,…,h.
(2)
Conversion time constraints between composite tasks, it is(9)weil+|gil-gi(l+1)|ϑi+di≤wsi(l+1),0∀i=1,…,NS;l=1,…,NBi-1.
Each composite task corresponds to an observation activity of the satellite, and every two composite tasks must satisfy the corresponding conversion time constraint. The conversion time include the rotation time and settling time of the sway.
(3)
Storage constraints of the satellite:(10)∑l=1NBi(weil-wsil)mi≤Mi,∀i=1,…,NS.
(4)
Energy constraints of the satellite:(11)∑l=1NBi(weil-wsil)pi+∑l=1NBi-1|gi(l+1)-gil|ρi≤Pi,00000000000000000000000000∀i=1,…,NS.The satellites usually consume some energy during the imaging and swaying process, so the corresponding energy consumption can be represented in a function about the observation time of satellites and sensor’s swaying. The energy consumption should not exceed the largest energy limitation.
(5)
Constraint on the maximum number of side-view imaging in a single circle. We can obtain the identifier of orbit circles according to the starting time and ending time of observation activities and then verify this constraint. Suppose the orbit set of satellitesi is
(12)Orbiti={orbiti1,orbiti2,…,orbitiNOi}
in whichNOi is the number of orbits, and orbitir refers to orbit r of satellite si. Given that orbitir contains Nir composite tasks {bir1,bir2,…,birNir},
(13)Nir≤Ri,∀i=1,…,NS,r=1,…,NOi.


## 2.1. Hypothesis

(1)
There is only one sensor in each satellite. In fact, if there are multiple sensors in one satellite, then this satellite can be divided into multiple satellites with one sensor. For this reason, both the sensorsi and the satellite si have the same meaning.
(2)
The regional target can be divided into multiple atomic tasks. In fact, a spot target just is one atomic task. In this paper, the atomic taskoijkv means the vth subtask of task tj which can be observed in the kth time window of satellite si. The regional target tj is divided into the following sets:
(1)Oj=⋃i=1NS⋃k=1Nij⋃v=1Nijkoijkv,
whereNS denotes the number of satellites, Nij denotes the number of visible time windows about sensor si and target tj, and Nijk denotes the number of atomic tasks which can be observed in the kth time window of satellite si.
(3)
In order to improve the performance of imaging resources, multiple atomic tasks are integrated into one composite task, which can be observed at one time.


## 2.2. Inputs

(1)
The number of satellites isNS. Spani denotes the longest working time of satellite si within one single booting.
(2)
M
i and Pi denote the maximum of storage capacity and energy of satellite si,mi, and pi refer to the storage and energy consumption of satellite si for observation in unit time.
(3)
R
i denotes the maximum number of side-view imaging in each circle, ϑi is the swaying speed of satellite si (degree per second), ρi refers to the energy consumption of swaying in each degree, di is the settling time of swaying of satellite si.
(4)
The number of tasks isNT. Suppose the first h targets are spot ones, and the last NT-h targets are regional ones. LPi means the priority of target i.
(5)
V
W
i
j
=
{
v
w
i
j
1
,
v
w
i
j
2
,
…
,
v
w
i
j
N
i
j
} denotes the set of time windows about sensor si and target tj,and Nij is the total number of visible time windows about sensor si and target tj.
(6)
w
s
i
j
k
v
,
w
e
i
j
k
v, and gijkv refer to the starting time, ending time, and observation angle of atomic task oijkv, respectively.
(7)
b
i
l denotes the lth composite task of satellite si, while wsil,weil, and gil refer to the starting time, ending time, and observation angle of composite task bil.
(8)
B
i
=
{
b
i
1
,
b
i
2
,
…
,
b
i
N
B
i
} denotes the set of composite tasks of satellite si, in which NBi is the number of composite tasks.


## 2.3. Outputs

In the multiple satellites mission planning problem, the outputs (decision variables) can be summarized as follows:(2)xijkvl={1,Theatomictaskoijkvwasexecutedinthelthcompositetasksofsatellitesi0,Otherwise.


## 2.4. Objectives

The objective of this problem is to maximize the priority summation of fulfilled tasks. That is,(3)maxJ=CSpot+CPolygon,
where CSpot denotes the priority summation of spot targets, and it is
(4)CSpot=∑j=1hCSpot(tj)=∑j=1h∑i=1NS∑k=1Nij∑l=1NBjxijk1l×LPj.CPolygon means the priority summation of regional targets. This computation process is a little complex. At first, we should compute the overlap area between the multiple observed atomic tasks and regional target tj. That is,
(5)SPolygon=(⋃i=1NS⋃k=1Nij⋃v=1Nijk⋃l=1NBixijkvlψ(oijkv))⋂ψ(tj).
Here, ψ(oijkv) denotes the polygon area of atomic task oijkv, ∩ denotes the intersection operation of multiple sets, ψ(tj) denotes the polygon area of regional target tj, and ∪ denotes the union operation of different sets. At second, we should compute the coverage rate of each regional target. It is
(6)Cover(tj)=f(SPolygon)f(tj),
where f(SPolygon) denotes the proportion of polygon area SPolygon, and f(tj) denotes the proportion of regional target tj. In total, the priority summation of regional targets can be computed as follows:
(7)CPolygon=∑j=h+1NT[Cover(tj)×LPj].


## 2.5. Constraints

(1)
Uniqueness constraint on spot targets; that is, each spot target can only be observed once:(8)∑i=1NS∑k=1Nij∑l=1NBjxijk1l≤1,∀j=1,…,h.
(2)
Conversion time constraints between composite tasks, it is(9)weil+|gil-gi(l+1)|ϑi+di≤wsi(l+1),0∀i=1,…,NS;l=1,…,NBi-1.
Each composite task corresponds to an observation activity of the satellite, and every two composite tasks must satisfy the corresponding conversion time constraint. The conversion time include the rotation time and settling time of the sway.
(3)
Storage constraints of the satellite:(10)∑l=1NBi(weil-wsil)mi≤Mi,∀i=1,…,NS.
(4)
Energy constraints of the satellite:(11)∑l=1NBi(weil-wsil)pi+∑l=1NBi-1|gi(l+1)-gil|ρi≤Pi,00000000000000000000000000∀i=1,…,NS.The satellites usually consume some energy during the imaging and swaying process, so the corresponding energy consumption can be represented in a function about the observation time of satellites and sensor’s swaying. The energy consumption should not exceed the largest energy limitation.
(5)
Constraint on the maximum number of side-view imaging in a single circle. We can obtain the identifier of orbit circles according to the starting time and ending time of observation activities and then verify this constraint. Suppose the orbit set of satellitesi is
(12)Orbiti={orbiti1,orbiti2,…,orbitiNOi}
in whichNOi is the number of orbits, and orbitir refers to orbit r of satellite si. Given that orbitir contains Nir composite tasks {bir1,bir2,…,birNir},
(13)Nir≤Ri,∀i=1,…,NS,r=1,…,NOi.


## 3. The Proposed Approach

In the existing intelligent optimization methods, they have not involved the effective management (i.e., learning, storage, and utilization) to domain knowledge of practical problems, and thus they cannot effectively obtain the optimal solution of optimization problems. Based on our knowledge, the complex optimization problem can be solved through learning some related knowledge from the optimization process and then employing it to guide the subsequent optimization process. In this paper, we mainly focus on the explicit knowledge, namely, the knowledge that can be definitely expressed. All the explicit knowledge can be expressed and stored through computer language and updated and adopted by other models.In this paper, the knowledge-based simulated annealing algorithm is defined as a hybrid approach which effectively combines the simulated annealing model and knowledge model, in which the simulated annealing model searches the solution space of optimization problems by “neighborhood search” strategy, and the knowledge model learns useful knowledge from the optimization process and then uses it to guide the subsequent optimization process. The knowledge-based simulated annealing algorithm employs the modeling idea of integrating the knowledge model and the simulated annealing model. Its basic framework is shown in Figure1.Figure 1
The basic framework of knowledge-based simulated annealing algorithm.The working mechanism of knowledge-based simulated annealing algorithm is shown in Figure2. The top illustrates the optimization process of simulated annealing algorithm: search the feasible space of optimization problems through the “neighborhood search” strategy and then converge to the near-optimal solution through continually iterations. And the bottom shows the effect of knowledge model: learn the useful knowledge from the early optimization process and then apply it to guide the subsequent optimization process.Figure 2
The operation mechanism of knowledge-based simulated annealing algorithm.During the early operation process of knowledge-based simulated annealing algorithm, due to that the available samples are not sufficient, so the reliability of learned knowledge is not high, and the guidance to optimization process is not distinct. However, with the advance of iterations, the reliability of learned knowledge is higher than before, and the guidance to optimization process is more distinct. Compared with traditional simulated annealing algorithm, with the help of knowledge model, the knowledge-based simulated annealing algorithm can either converge to a satisfactory solution more quickly or converge to a satisfactory solution with higher quality.


### 3.1. Knowledge Definition, Learning, and Application

(
1) Operator Knowledge. Usually, different operators have different application scopes. When solving practical problems, a single operator may be more effective to some problems than to others. Therefore, it is very difficult to find out a general operator that can solve all kinds of practical problems effectively. In order to improve the performance of simulated annealing algorithm, we can adopt several kinds of operators for different operations in simulated annealing algorithm, with the hope of trying to mine some operators that can effectively solve the present problem.When an operator is employed to one operation, suppose the original individual set isCB and the newly generated individual set is CA, and if the optimal individual in CA is better than that in CB, then the present operation is considered to be successful. The optimization performance of one operator can be defined as the number of successful operations realized by the given operator. Operator knowledge mainly refers to the accumulated knowledge for the optimization performance of various operators.The simple case for extraction and application of operator knowledge is displayed as follows. Suppose there are three operators in the simulated annealing algorithm, if the current operation implemented by thekth operator is successful, then the optimization performance N(k) of the kth operator will be added by one. During the next operation, the selection probability for each operator will be calculated according to the following formula:
(14)P(i)=N(i)∑i=13N(i),
where P(i) denotes the probability that the ith operator is selected for the operation.(
2) Parameter Knowledge. One simple and easy method to improve the optimization performance of simulated annealing algorithm is to adjust the parameter values dynamically during the optimization process. To lower the sensitivity of parameters to experimental results, we employ several different parameter combinations to develop the evolution process and decide the parameter combination selected for the next iteration according to their optimization performance. After a single iteration, if the global optimal solution has been improved, then it is a successful iteration. In the process of solving current cases, the number of successful iterations obtained by the given parameter combination is regarded as its optimization performance. Parameter knowledge mainly refers to the accumulated knowledge of the optimization performance for various parameter combinations.During the initialization stage, the orthogonal method can be adopted to generate different parameter combinations and all their optimization performance is initialized to 1. Before any iteration of simulated annealing algorithm, a parameter combination will be selected stochastically from several parameter combinations through roulette selection (based on the optimization performance of parameter combinations) as the parameter of current iteration. If the global optimal solution has been improved in this iteration, then the optimization performance value of the parameter combination used in this iteration will be increased.


### 3.2. The Available Neighborhood Structures

The following neighborhood structures have been designed based on the basic ones and passed the constraint checking of the basic neighborhood structure, thus greatly simplifying the optimization calculation process.(
1) Arrange. The scheduling of imaging satellites for composing tasks usually involves the task arrangement which is performed through composing or inserting the contained atomic tasks. Neighborhood composition is preferred in the arrangement of atomic tasks, and the insertion operation is adopted only when atomic tasks cannot be composed or the task has not been finished, which contributes to the observation composition of much more tasks and helps to arrange more tasks.(
2) Relocate. This kind of neighborhood can reallocate the satellite resource and time windows of tasks. There are several time windows about satellites and tasks, so we can reallocate time windows for the task as to the same satellite. Two kinds of operations are involved here: one is to reallocate the task to different time windows about the same satellite and the other is to reallocate the task to another satellite. The relocating aims to increase the chance of arranging other tasks through adjusting the satellite resource and time windows of tasks.(
3) Swap. This neighborhood aims to exchange two tasks. About parallel scheduling problems, three kinds of neighborhood structures for task exchange are mentioned: the exchange between two tasks with the same resource (exchanging the execution order); the exchange between two tasks with different resources; the exchange between the scheduled tasks and unscheduled ones. To the multiple satellites mission planning problem, the relevant experiments indicate that the first two exchanges are unsuccessful. Therefore, this paper adopts the third kind of neighborhood structure for task exchange, namely, to exchange tasks under the real satellite resource and the virtual one.(
4) Sample. Random sampling neighborhood is a special measure adopted to handle large-scale neighborhood effectively. It is not a new neighborhood structure itself, but a method to generate new neighbourhood structures through sampling the existing ones randomly. To be specific, before each move, random sampling neighbourhood will sample all neighbors generated by the given neighbourhood structure, so as to limit the move to part of the sampled neighbourhood, and the proportion of sampling is usually specified in advance upon requirements.Among the designs of neighbourhood structure, the arrangement of neighbourhood is designed to implement the insertion and composition operations among atomic tasks in which the task composition is preferred, and it is a kind of improved neighbourhood; the neighbourhood relocation and swap disturb the solutions, respectively, through adjusting the resources and time windows of the task as well as exchanging the scheduled tasks with unscheduled ones, and they belong to a regulative neighbourhood which can increase the chance to arrange other tasks after the relocation and swap of tasks, so they should be tried on all the unfinished tasks.


## 3.1. Knowledge Definition, Learning, and Application

(
1) Operator Knowledge. Usually, different operators have different application scopes. When solving practical problems, a single operator may be more effective to some problems than to others. Therefore, it is very difficult to find out a general operator that can solve all kinds of practical problems effectively. In order to improve the performance of simulated annealing algorithm, we can adopt several kinds of operators for different operations in simulated annealing algorithm, with the hope of trying to mine some operators that can effectively solve the present problem.When an operator is employed to one operation, suppose the original individual set isCB and the newly generated individual set is CA, and if the optimal individual in CA is better than that in CB, then the present operation is considered to be successful. The optimization performance of one operator can be defined as the number of successful operations realized by the given operator. Operator knowledge mainly refers to the accumulated knowledge for the optimization performance of various operators.The simple case for extraction and application of operator knowledge is displayed as follows. Suppose there are three operators in the simulated annealing algorithm, if the current operation implemented by thekth operator is successful, then the optimization performance N(k) of the kth operator will be added by one. During the next operation, the selection probability for each operator will be calculated according to the following formula:
(14)P(i)=N(i)∑i=13N(i),
where P(i) denotes the probability that the ith operator is selected for the operation.(
2) Parameter Knowledge. One simple and easy method to improve the optimization performance of simulated annealing algorithm is to adjust the parameter values dynamically during the optimization process. To lower the sensitivity of parameters to experimental results, we employ several different parameter combinations to develop the evolution process and decide the parameter combination selected for the next iteration according to their optimization performance. After a single iteration, if the global optimal solution has been improved, then it is a successful iteration. In the process of solving current cases, the number of successful iterations obtained by the given parameter combination is regarded as its optimization performance. Parameter knowledge mainly refers to the accumulated knowledge of the optimization performance for various parameter combinations.During the initialization stage, the orthogonal method can be adopted to generate different parameter combinations and all their optimization performance is initialized to 1. Before any iteration of simulated annealing algorithm, a parameter combination will be selected stochastically from several parameter combinations through roulette selection (based on the optimization performance of parameter combinations) as the parameter of current iteration. If the global optimal solution has been improved in this iteration, then the optimization performance value of the parameter combination used in this iteration will be increased.


## 3.2. The Available Neighborhood Structures

The following neighborhood structures have been designed based on the basic ones and passed the constraint checking of the basic neighborhood structure, thus greatly simplifying the optimization calculation process.(
1) Arrange. The scheduling of imaging satellites for composing tasks usually involves the task arrangement which is performed through composing or inserting the contained atomic tasks. Neighborhood composition is preferred in the arrangement of atomic tasks, and the insertion operation is adopted only when atomic tasks cannot be composed or the task has not been finished, which contributes to the observation composition of much more tasks and helps to arrange more tasks.(
2) Relocate. This kind of neighborhood can reallocate the satellite resource and time windows of tasks. There are several time windows about satellites and tasks, so we can reallocate time windows for the task as to the same satellite. Two kinds of operations are involved here: one is to reallocate the task to different time windows about the same satellite and the other is to reallocate the task to another satellite. The relocating aims to increase the chance of arranging other tasks through adjusting the satellite resource and time windows of tasks.(
3) Swap. This neighborhood aims to exchange two tasks. About parallel scheduling problems, three kinds of neighborhood structures for task exchange are mentioned: the exchange between two tasks with the same resource (exchanging the execution order); the exchange between two tasks with different resources; the exchange between the scheduled tasks and unscheduled ones. To the multiple satellites mission planning problem, the relevant experiments indicate that the first two exchanges are unsuccessful. Therefore, this paper adopts the third kind of neighborhood structure for task exchange, namely, to exchange tasks under the real satellite resource and the virtual one.(
4) Sample. Random sampling neighborhood is a special measure adopted to handle large-scale neighborhood effectively. It is not a new neighborhood structure itself, but a method to generate new neighbourhood structures through sampling the existing ones randomly. To be specific, before each move, random sampling neighbourhood will sample all neighbors generated by the given neighbourhood structure, so as to limit the move to part of the sampled neighbourhood, and the proportion of sampling is usually specified in advance upon requirements.Among the designs of neighbourhood structure, the arrangement of neighbourhood is designed to implement the insertion and composition operations among atomic tasks in which the task composition is preferred, and it is a kind of improved neighbourhood; the neighbourhood relocation and swap disturb the solutions, respectively, through adjusting the resources and time windows of the task as well as exchanging the scheduled tasks with unscheduled ones, and they belong to a regulative neighbourhood which can increase the chance to arrange other tasks after the relocation and swap of tasks, so they should be tried on all the unfinished tasks.


## 4. Experimental Results

Six sets of even-distributed targets are generated, respectively, with the number of 100, 200, 300, 400, 500, and 600. Based on a different number of satellites, several test cases are then constructed for each set of targets. Each case will be calculated for fifty times and the comparison will be made on the averages of the results to avoid the influence of randomness. To validate the performance of this algorithm, the Heuristic Method (HM) and the Standard Simulated Annealing Algorithm (SSAA) are employed to be compared with the Knowledge-Based Simulated Annealing Algorithm (KSAA).


### 4.1. Experimental Results of Small Scale Testing Examples

Six small scale testing examples(100×2,100×4,100×6,200×2,200×4,200×6) are applied to test the performance of different algorithms. Experimental results of small scale testing examples are displayed in Figures 3 and 4. Experimental results in Figure 3 suggest that there is no remarkable difference in optimization performance among HM, SSAA, and KSSA. Experimental results in Figure 4 suggest that the computational time of SSAA and KSSA is little larger than that of HM. Also, there is no distinct difference in computational time between SSAA and KSSA. Since SSAA and KSSA are all intelligent optimization approaches, so the computational time of SSAA and KSSA is larger than that of HM.Figure 3
Experimental results in optimization performance of small scale testing examples.Figure 4
Experimental results in computational time of small scale testing examples.


### 4.2. Experimental Results of Medium Scale Testing Examples

Six medium scale testing examples(300×2,300×4,300×6,400×2,400×4,400×6) are applied to test the performance of different algorithms. Experimental results of medium scale testing examples are displayed in Figures 5 and 6. Experimental results in Figure 5 suggest that there is no significant difference in optimization performance among SSAA and KSSA, but the optimization performance of SSAA and KSSA is better than that of HM. Since SSAA and KSSA are all intelligent optimization approaches, so the optimization performance of SSAA and KSSA should be better than that of HM, and the experimental result indicates this viewpoint. Experimental results in Figure 6 suggest that the computational time of SSAA and KSSA is little larger than that of HM. Also, there is no distinct difference in computational time between SSAA and KSSA.Figure 5
Experimental results in optimization performance of medium scale testing examples.Figure 6
Experimental results in computational time of medium scale testing examples.


### 4.3. Experimental Results of Large Scale Testing Examples

Six large scale testing examples(500×4,500×6,500×8,600×4,600×6,600×8) are applied to test the performance of different algorithms. The experimental results of large scale testing examples are displayed in Figures 7 and 8. Experimental results in Figure 7 suggest that the optimization performance of SSAA is better than that of HM, and the optimization performance of KSAA is better than that of SSAA. Since SSAA is an intelligent optimization approach, so the optimization performance of SSAA should be better than that of HM, and experimental results indicate this viewpoint. At the same time, KSAA extended the learning and guidance function based on the SSAA, so its performance should be better than that of SSAA, and experimental results validate this opinion. Experimental results in Figure 8 suggest that the computational time of SSAA is larger than that of HM, and the computational time of KSAA is larger than that of SSAA. Since SSAA is an intelligent optimization approach, so it needs more computation cost than that of HM to deal with large scale testing examples. Also, the number of modules in KSSA is more than that in SSAA, so it needs more computation cost than that of SSAA to deal with large scale testing examples.Figure 7
Experimental results in optimization performance of large scale testing examples.Figure 8
Experimental results in computational time of large scale testing examples.


## 4.1. Experimental Results of Small Scale Testing Examples

Six small scale testing examples(100×2,100×4,100×6,200×2,200×4,200×6) are applied to test the performance of different algorithms. Experimental results of small scale testing examples are displayed in Figures 3 and 4. Experimental results in Figure 3 suggest that there is no remarkable difference in optimization performance among HM, SSAA, and KSSA. Experimental results in Figure 4 suggest that the computational time of SSAA and KSSA is little larger than that of HM. Also, there is no distinct difference in computational time between SSAA and KSSA. Since SSAA and KSSA are all intelligent optimization approaches, so the computational time of SSAA and KSSA is larger than that of HM.Figure 3
Experimental results in optimization performance of small scale testing examples.Figure 4
Experimental results in computational time of small scale testing examples.


## 4.2. Experimental Results of Medium Scale Testing Examples

Six medium scale testing examples(300×2,300×4,300×6,400×2,400×4,400×6) are applied to test the performance of different algorithms. Experimental results of medium scale testing examples are displayed in Figures 5 and 6. Experimental results in Figure 5 suggest that there is no significant difference in optimization performance among SSAA and KSSA, but the optimization performance of SSAA and KSSA is better than that of HM. Since SSAA and KSSA are all intelligent optimization approaches, so the optimization performance of SSAA and KSSA should be better than that of HM, and the experimental result indicates this viewpoint. Experimental results in Figure 6 suggest that the computational time of SSAA and KSSA is little larger than that of HM. Also, there is no distinct difference in computational time between SSAA and KSSA.Figure 5
Experimental results in optimization performance of medium scale testing examples.Figure 6
Experimental results in computational time of medium scale testing examples.


## 4.3. Experimental Results of Large Scale Testing Examples

Six large scale testing examples(500×4,500×6,500×8,600×4,600×6,600×8) are applied to test the performance of different algorithms. The experimental results of large scale testing examples are displayed in Figures 7 and 8. Experimental results in Figure 7 suggest that the optimization performance of SSAA is better than that of HM, and the optimization performance of KSAA is better than that of SSAA. Since SSAA is an intelligent optimization approach, so the optimization performance of SSAA should be better than that of HM, and experimental results indicate this viewpoint. At the same time, KSAA extended the learning and guidance function based on the SSAA, so its performance should be better than that of SSAA, and experimental results validate this opinion. Experimental results in Figure 8 suggest that the computational time of SSAA is larger than that of HM, and the computational time of KSAA is larger than that of SSAA. Since SSAA is an intelligent optimization approach, so it needs more computation cost than that of HM to deal with large scale testing examples. Also, the number of modules in KSSA is more than that in SSAA, so it needs more computation cost than that of SSAA to deal with large scale testing examples.Figure 7
Experimental results in optimization performance of large scale testing examples.Figure 8
Experimental results in computational time of large scale testing examples.


## 5. Conclusions

The innovation of this paper lies in that it has established a multisatellite mission planning model based on the constraint satisfaction and proposed a knowledge-based simulated annealing algorithm to solve multisatellite mission planning problems, which has been proved effective in the relevant experiment.The future research directions can be summarized as follows. (1) Extend the types of knowledge. Extracting the experiential knowledge of experts to multiple satellites mission planning problems, we could employ this knowledge to guide the optimization process effectively so as to enhance optimization efficiency as much as possible. (2) Adopt new knowledge mining techniques. We should consider employing advanced knowledge mining through machine learning or data mining to obtain useful knowledge from the optimization process.



---

*Source: 913450-2013-11-27.xml*