敏捷对地观测卫星(agile earth observation satellite, AEOS)拥有良好的姿态机动能力，在民用和军用领域都有着广阔的应用前景，比如环境保护、国土普查、抗震减灾以及制空权、制天权、制海权等。在AEOS的运行过程中，对AEOS提前根据目标点特性进行任务规划至关重要，可有效发挥AEOS性能，获取更多的观测收益^[1]。而在任务规划过程中，AEOS需要对多个目标点进行规划，同时也要考虑到卫星在观测过程中所面临的多种状态、控制等复杂约束条件，因此，国内外越来越多的学者开始进行卫星任务规划方面的研究。

Michel等^[2]提出了对点目标以及区域目标的分割和观测方法，综合分析了时间约束、姿态机动约束、存储容量和能量的约束，同时考虑到了天气的不确定性，并提出了相应的优化目标函数。他们分别对比了贪婪算法、动态规划算法、局部搜索算法这几种算法的求解效率。赵琳等^[3]对敏捷卫星的单星单轨任务规划问题进行了研究，引入了任务—姿态协同规划思想，并根据任务—姿态协同规划数学模型设计了自适应伪谱遗传算法(APGA)，用以求解满足调整时间最优的任务规划问题。耿远卓等^[4]利用团划分算法对多个点目标进行了聚类，考虑了卫星的时间窗口约束以及最大加速度约束，在传统的蚁群算法的基础上引入了启发式寻优策略和新的信息素更新策略，加快了算法的收敛性。丁祎男等^[5]面向单星的任务规划问题，在综合分析遗传算法与禁忌搜索算法优缺点的基础上，提出了遗传禁忌混合算法，该算法能够解决算法的早熟问题，同时能够加快算法收敛速度。王海蛟等^[6]针对敏捷卫星的调度问题，采用了二进制与实数杂合的编码方式，将量子优化机制引入了遗传算法中，提出了改进的量子遗传算法，提高了搜索效率。Tangpattanakul等^[7]针对单星的任务规划问题，提出了一种基于知识的多目标局部搜索方法(IBMOLS)，该方法实质上是局部搜索算法与遗传算法的综合方法，能够比遗传算法更快地收敛。Chu等^[8]研究了由低分辨率和高分辨率的卫星组成的双星任务规划问题，对双星协同任务规划问题进行建模，提出了分支限定的求解方法。Lee等^[9]面向多星任务规划问题，主要考虑了卫星能源以及存储容量的约束，利用遗传算法进行求解，就不同的应用场景进行了仿真。Zheng等^[10]将多星的星上任务规划问题建模为约束优化问题，同时在现有的用于卫星任务规划的遗传算法的基础上，提出了求解速度更快的混合动态变异遗传算法(HDMGA)。黄生俊等^[11]对知识定义、知识更新规则和任务冲突处理策略做了详细描述，并综合蚁群算法的反馈特性和模拟退火算法的局部搜索特性, 设计了一种基于知识的改进模拟退火算法。何磊等^[12]考虑了光学成像敏捷卫星中的云层遮挡问题，采用了预判和二分法进行云层遮挡时间窗口的计算，并针对这一问题提出了相应的蚁群算法进行求解，提高了光学成像卫星的成像效率。Tipaldi等^[13]结合了规划系统和动态重规划的能力深入分析了马尔科夫决策过程(MDP)模型的自主规划方法，验证了该方法在实际应用中的性能。郝会成等^[14]将免疫遗传算法与蚁群算法等相结合，提出了基于混合遗传求解算法，同样提高了算法的求解速度。Xu等^[15]考虑了卫星的资源约束条件，在此基础上提出了基于优先权的结构算法，用于最大化观测收益。董轩鸿^[16]研究了对观测目标的条带划分策略，设计了改进的遗传算法。

综合来看，目前已有很多文献应用遗传算法等智能优化算法来解决卫星的任务规划问题^{[3, 5, 7, 10]}。虽然应用遗传算法可以求解出卫星的最优观测序列，但是传统的遗传算法在交叉、变异等寻优过程仍有许多值得改进的地方，因此有必要对遗传算法的变异过程做出改进，进而提升算法的求解效率。同时，如果规划问题的适应度函数过于简单^[17]，往往不能真实地反映卫星的实际观测情况，因此有必要对适应度函数做出改进。

基于上述考虑，本文提出了禁忌退火遗传混合算法，并将其应用于求解AEOS的任务规划问题。首先提出了本文的适应度函数.此函数综合考虑了卫星在观测过程中的多种约束条件，能够较为真实地反映卫星的实际观测情况。其次介绍了禁忌退火遗传混合算法。此算法对遗传算法的变异过程做出了改进，将禁忌搜索算法与模拟退火算法的寻优过程的优势引入到了遗传算法的变异过程，提升了遗传算法的邻域搜索能力，节省了算法的运行时间，适用于AEOS的任务规划问题。

1 任务规划问题建模 1.1 任务规划问题描述

AEOS在运行过程中，主要通过横滚轴和俯仰轴进行姿态机动。星上携带有摄像头，可以对地面实施观测任务。如图 1所示，卫星沿滚动轴的最大姿态机动角度为θ_max/2，沿俯仰轴的最大姿态机动角度为ξ_max/2。当卫星经过目标点上空时，需要提前对目标点的观测序列做出规划，得到满足观测约束下的最优/次优观测序列。

卫星观测目标点时，会获得相应的收益，在对单目标的持续观测和目标间的姿态机动过程中，卫星需要消耗一定能量。因此，构建的目标优化函数必须综合考虑观测收益及观测代价(即能量消耗)，在满足各种约束条件情况下，通过优化算法获得卫星对目标点的观测序列。

1.2 适应度函数的建立

定义如下变量：M为目标点集合，M={m₁, m₂, …, m_n}，其中m_i为第i个目标点，n为目标点总数；t_twsi为卫星对m_i可见的时间窗口的开始时间；t_twei为卫星对m_i可见的时间窗口的结束时间；t_si为卫星对m_i的实际开始观测时间；t_ei为卫星对m_i的实际结束观测时间；d_i为卫星对m_i的实际观测的持续时间，即d_i=t_ei-t_si；t_rij为卫星从观测m_i到观测m_j转过的姿态机动角度；p_i为m_i的任务优先级，p_i取值越大，意味着m_i的优先级越高；w_i为卫星观测m_i获得的收益，即w_i=d_i·p_i；pos_i=(lgt_i, lat_i)为m_i的经、纬度坐标，其中lgt_i为m_i的经度；lat_i为m_i的纬度，决策变量s_i和F_ij的定义如下：

$ s_{i}=\left\{\begin{array}{l} 1, m_{i} \text { 被观测 } \\ 0, m_{i} \text { 不被观测 } \end{array}\right. $

(1)

$ F_{i j}=\left\{\begin{array}{l} 1, m_{j} \text { 紧接着 } m_{i} \text { 被观测 } \\ 0, \text { 其他情况 } \end{array}\right. $

(2)

基于上述定义的变量，考虑如下适应度函数:

$ f=\sum\limits_{i=1}^{n} \sum\limits_{j=1}^{i}\left[s_{i} \cdot w_{i}-\eta\left(f_{1}+s_{i} \cdot d_{i}\right)\right] $

(3)

上式表示的函数由观测收益与卫星姿态机动的能量消耗两部分组成. 式中：s_i·w_i为卫星观测某一目标点所带来的收益，收益与目标点优先级成正比；f₁+s_i·d_i为卫星观测过程中消耗的能量，其中s_i·d_i为持续观测某一目标点的能耗，f₁为卫星姿态机动的能耗，与卫星姿态机动过程中匀加速与匀减速的时间Δt有关；η为卫星的能耗系数，是一个常数。Δt的计算方法如下。

本文假设卫星在目标点之间姿态机动时，每个轴均以最大力矩T_max进行机动，根据卫星需要转过的角度大小分为两种情况考虑，如图 2的两种情况所示。图 2(a)对应大角度机动情况，图 2(b)对应小角度的情况。

综上所述，考虑f₁为

$ f_{1}=F_{i j} \cdot \Delta t \cdot T_{\max }^{2} $

(4)

式中，Δt为卫星进行匀加速与匀减速过程的总时间。综合式(3)、(4)可得

$ f=\max \sum\limits_{i=1}^{n} \sum\limits_{j=1}^{i}\left[s_{i} \cdot w_{i}-\eta\left(F_{i j} \cdot \Delta t \cdot T_{\max }^{2}+s_{i} \cdot d_{i}\right)\right] $

(5)

从式(5)可以看出，f取得最大值是本文的优化目标，即在考虑到卫星的能量消耗的情况下尽可能多地获得观测收益。

1.3 AEOS任务规划的约束条件

下面给出AEOS任务规划的约束条件。

如果目标点m_i被观测，那么观测的持续时间区间必须被包含在卫星对m_i可见的时间窗口内，即

$ t_{t w s i} \leqslant t_{s i}, t_{t w e i} \geqslant t_{e i}, \text { 若 } s_{i}=1 $

(6)

卫星结束观测m_i以后，必须在m_j的时间窗口结束之前完成姿态机动才可以观测m_j，即

$ t_{s i}+t_{r i j}+d_{i} \leqslant t_{t w e j}, \text { 若 } F_{i j}=1 $

(7)

考虑如下卫星轨道动力学约束：

$ \ddot{\boldsymbol{r}}_{s}+\frac{\mu}{\boldsymbol{r}_{s}^{2}}=0 $

(8)

以及如下卫星姿态角约束、姿态角速度约束、姿态角加速度约束:

$ |\theta| \leqslant \theta_{\max },|\xi| \leqslant \xi_{\max } $

(9)

$ \omega_{\theta} \leqslant \omega_{\theta \max }, \omega_{\xi} \leqslant \omega_{\xi \max } $

(10)

$ \alpha_{\theta} \leqslant \alpha_{\theta \max }, \alpha_{\xi} \leqslant \alpha_{\xi \max } $

(11)

此外，还需要考虑如下卫星姿态转移时间约束:

$ t_{r i j}=f\left(r_{s}, \operatorname{pos}_{i}, \operatorname{pos}_{j}, T_{\max }\right) $

(12)

式中，时间t_rij与两个目标点的经纬度以及AEOS的最大机动力矩T_max相关，计算方法如图 2所示。其中AEOS需要机动的角度θ的求解采用了文献[17]的方法，如图 3所示。

在图 3中，R_ei和R_ej分别为m_i和m_j在地心惯性坐标系(ECI)中的坐标，具体公式为

$ \boldsymbol{R}_{e i}=\left[\begin{array}{ccc} \cos \omega t & -\sin \omega t & 0 \\ \sin \omega t & \cos \omega t & 0 \\ 0 & 0 & 1 \end{array}\right]\left[\begin{array}{c} R_{e} \cos \operatorname{lat}_{i} \cos \operatorname{lgt}_{i} \\ R_{e} \cos \operatorname{lat}_{i} \sin \lg t_{i} \\ R_{e} \sin \operatorname{lat}_{i} \end{array}\right] $

(13)

式中: ω为地球自转角速度，R_e为地球半径。

图 3中R_es为卫星在ECI中的坐标，可通过下式求得：

$ \boldsymbol{R}_{e s}=R_{s}\left[\begin{array}{l} \cos \Omega \cos (\psi+\alpha)-\sin \Omega \sin (\psi+\alpha) \cos i \\ \sin \Omega \cos (\psi+\alpha)+\cos \Omega \sin (\psi+\alpha) \cos i \\ \sin (\psi+\alpha) \sin i \end{array}\right] $

(14)

式中: R_s为卫星的轨道半径，Ω为升交点赤经，ψ为近地点幅角，α为卫星转过的角度，可由下式确定：

$ \alpha=E+\omega_{s} t $

(15)

式中: E为卫星开始观测时刻的真近点角，ω_s为卫星的角速度，t为当前的时刻。

由图 3中可得下式

$ \left\{\begin{array}{l} \boldsymbol{R}_{s i}=\boldsymbol{R}_{e i}-\boldsymbol{R}_{e s} \\ \boldsymbol{R}_{s j}=\boldsymbol{R}_{e j}-\boldsymbol{R}_{e s} \end{array}\right. $

(16)

因此，卫星从观测m_i姿态机动到m_j所需转动的角度θ可通过下式求得：

$ \theta=\arccos <\boldsymbol{R}_{s i}, \boldsymbol{R}_{s j}> $

(17)

综合式(17)与图 2，即可最终求得卫星的姿态机动时间t_rij。

2 禁忌退火遗传混合算法

对于大规模的卫星任务规划问题，传统的遗传算法很难在较短的运算时间内给出最优的观测序列。因此，有必要对遗传算法的寻优过程做出改进。相对于传统的遗传算法变异方法，本文提出的禁忌退火变异方法能够有效提高算法的求解搜索效率，节省算法的运行时间。

2.1 初始种群生成

本文采用整数编码来解决卫星的任务规划问题。整数编码与二进制编码等其他编码方式相比，能够更直观地表示卫星的观测序列，也便于算法搜寻邻域解。

采用整数编码方式的每个染色体(也称作解、个体)分别代表一种可行的观测序列。染色体上的每个基因所对应的数字即为目标点，如图 4所示。

图 4代表一条长度为6的染色体，表示卫星需要观测6个目标点，且最先5号目标点，最后观测6号目标点。

设地面上目标点总数为n，AEOS需要从n个目标点中选取N个目标点进行观测，种群规模为M。种群初始化的方法为，随机生成M个观测序列，其中每个观测序列都是N个不大于n且互不重复的正整数，即同一个目标点不会重复被观测。在此后算法的每一次迭代过程中，都会首先计算种群中每一条染色体所对应的适应度函数值。

2.2 适应度函数计算

适应度函数值的大小是衡量种群中每个个体优劣程度的重要指标。对于每一个给定的观测序列，为了计算此序列所对应的适应值，需要首先确定卫星能否依次观测此序列所有的目标点。过程如下：

Step1  比较当前的时间t与m_i的时间窗口结束时间t_twei。如果t>t_twei，则m_i无法被观测，令s_i=0，i=i+1，继续执行step1。如果t < t_twei，则执行step2。

Step2  比较t与t_twsi。如果t>t_twsi，则执行step3。如果t < t_twsi，则执行step4。

Step3  计算t_rij。如果t+t_rij>t_twei，则m_i无法被观测，令s_i=0，i=i+1，执行step1。如果t+t_rij < t_twei则m_i可以被观测，令s_i=1，t=t+t_rij+d_i，i=i+1执行step1。

Step4  计算t_rij。如果t+t_rij < t_twei，则m_i可以被观测，s_i=1，执行step5。如果t+t_rij>t_twei，则m_i无法被观测，令s_i=0，i=i+1，执行step1。

Step5  如果t+t_rij < t_twsi，令t=t_twsi+d_i，i=i+1，执行step1。如果t+t_rij>t_twsi，令t=t+t_rij+d_i，i=i+1，执行step1。

从i=1起，重复执行上述过程，直至i=N，即可得到每个观测序列所对应的s_i以及F_ij。将s_i、F_ij以及关于每个目标点的w_i等变量带入式(5)即可求得每个观测序列所对应的适应值的大小，每个个体的优劣程度也就可以依次确定。

2.3 选择、交叉操作

本文采用精英保留策略以及轮盘赌法对种群进行选择、交叉操作。精英保留策略将第g代种群中的若干优秀个体当做父本直接传递到第g+1代种群。随后根据轮盘赌法从父本中选择父本进行交叉运算，生成若干子代，使得第g+1代的种群规模为M。

轮盘赌法的计算过程如下所示：

$ Q_{i}=\frac{F_{i}}{\sum\limits_{i=1}^{n} F_{i}} $

(18)

$ P_{i}=\sum\limits_{j=1}^{i} Q_{j} $

(19)

式中: F_i为第i个父本对应的适应值，Q_i为第i个父本被选择的概率，P_i为从第1个父本直至第i父本的累加概率。其中，各父本按照适应度从大到小的顺序排列。

第i个父本被选中参与交叉过程的方法如下：产生一个介于0和1之间的随机数r，如果该随机数满足P_i-1 < r < P_i，则意味着第i个父本被选中，参与到交叉过程.轮盘赌法的操作过程也意味着越优秀的父本被选中进行交叉的几率越高。当两个父本被选中后，随机产生一个交叉点，这两个父本互相交换交叉点两边的染色体。交叉过程如图 5所示。

经过交叉操作得到的子代往往都含有相同的基因。比如图 5中，子代1的6和3都是重复出现的元素。因此，需要进行去重复操作，使每个子代的各个基因之间都两两互不重复。去重复操作是将每个子代中重复的基因随机用没有出现在该子代染色体中的基因替代，如图 6所示。

2.4 禁忌退火变异

本文给出了禁忌退火变异(tabu search-simulated annealing mutation，TSSAM)。首先给出了禁忌搜索算法和模拟退火算法的相关定义；其次介绍了TSSAM的具体过程；最后分析了TSSAM对于传统的变异方法的改进之处。

2.4.1 禁忌搜索算法

禁忌搜索算法(tabu search algorithm, TS)寻优策略为，从初始解的邻域解集中选取最优解或者是未被禁忌的最优解作为新的初始搜索解，直至算法收敛。禁忌搜索算法通过不断更新禁忌表来避免在一个邻域内重复搜索，进而提高算法的全局搜索能力，加快算法收敛速度。禁忌搜索算法的一些定义如下。

1) 初始解。参与禁忌搜索的最开始的解，在本文中即为参与变异的个体x₀。

2) 邻域解集。初始解x₀的邻域解构成的集合。对于采用整数编码方式的染色体来讲，邻域解采用2-opt形式，即通过互换两个基因位上的元素来产生邻域解，如图 7所示。邻域解集的规模与染色体的长度有关。

3) 候选解。x₀产生的邻域解集中的最优解或者是未被禁忌的最优解，用于同x₀比较，决定是否要用候选解替换x₀。

4) 禁忌表。用于记录被禁忌的前若干次的操作。对于采用整数编码的染色体来讲，禁忌表是一个N×N的矩阵。其中N为染色体长度。禁忌表中第i行第j列的数字a_{i, j}表示交换i, j两个位置上元素的操作当前被禁忌的次数。如果禁忌表中该位置上元素为0，则表示此操作当前未被禁忌。在算法最开始的时候，禁忌表初始化为零矩阵。

5) 禁忌长度。被禁忌的操作不允许被选取的最大次数。禁忌长度与染色体长度有关。禁忌长度记为L_tabu。

2.4.2 模拟退火算法

模拟退火算法(simulated annealing algorithm, SA)的寻优过程为，假定一个着火物体，其温度随迭代过程指数下降。每次迭代时，在当前解的邻域内随机搜寻可行解，利用Metropolis法则判断是否接受邻域内可行解，循环此过程直至算法收敛。

1) 退火温度。SA中物体的温度，记为T。物体的退火温度T从初始温度T₀开始，随迭代过程指数下降。

2) 退温率。物体的退火温度的衰减速度K。物体温度的退温策略为T(n+1)=K·T(n)，其中n为迭代次数。K满足0 < K < 1。

3) Metropolis法则。如果在当前解s₀的邻域内搜寻到的可行解s₁优于当前解，则以概率1接受s₁作为新的当前解。如果s₁劣于s₀，则以一定概率p接受s₁作为新的当前解。此概率p与当前的温度T以及s₁，s₀适应值的差有关，计算公式如下：

$ p=\exp (\Delta E / T) $

(20)

式中: ΔE为s₁与s₀适应度函数的差值，且ΔE < 0；T为当前系统的温度。

2.4.3 禁忌退火变异算法(TSSAM)

决定种群中某一个体是否参与变异的方法如下：生成随机数r，如果变异率大于r，那么对当前个体执行TSSAM操作。TSSAM步骤如下：

Step1  对于每一个参与变异的个体x₀，首先产生此个体的邻域解集。x₀也称作禁忌搜索的初始解。随后将x₀产生的邻域解按照从优到劣的顺序进行排列，记最优邻域解为x₁。先将x₁作为候选解。

Step2  将候选解x₁与初始解x₀进行比较。如果x₁的适应度函数值大于x₀，执行step3。如果x₁的适应值小于x₀，执行step4。

Step3  将x₁作为x₀变异后的个体传回到种群中，同时更新禁忌表。记由x₀得到x₁的方法为互换第i₁, j₁个基因位上的元素。更新禁忌表的方法为：对禁忌表中所有非零元素执行如下操作：

$ \left\{\begin{array}{l} a_{i, j}=a_{i, j}-1, \text { if } i \neq i_{1}, j \neq j_{1} \\ a_{i, j}=L_{\text {tabu }}, \text { if } i=i_{1}, j=j_{1} \end{array}\right. $

(21)

随后考虑下一个参加变异的个体，执行Step1。

Step4  搜索x₀的邻域解集中未被禁忌的最优解，将此最优解记为新的候选解x₂。将x₀与x₂的目标优化函数值的差记为ΔE，ΔE < 0。产生一个随机数r，如果满足exp(ΔE/T)>r，则接受x₂作为变异后的个体传回到种群中，同时更新禁忌表，随后考虑下一个参加变异的个体，执行step1。记由x₀得到x₂的方法为互换i₂, j₂基因位上的元素。更新禁忌表的方法为对禁忌表中所有非零元素执行如下操作：

$ \left\{\begin{array}{l} a_{i, j}=a_{i, j}-1, \text { if } i \neq i_{2}, j \neq j_{2} \\ a_{i, j}=L_{\text {tabu }}, \text { if } i=i_{2}, j=j_{2} \end{array}\right. $

(22)

生成随机数r，如果满足exp(ΔE/T) < r，则把初始解x₀传回种群中，即x₀不发生变异。同时更新禁忌表。随后考虑下一个参加变异的个体，执行step1。

重复执行上述过程，直至种群完成变异。

当种群中的所有参与变异的个体变异完成以后，更新退火温度，即令T(n+1)=K·T(n)，开始种群下一代的迭代进化过程。

传统遗传算法的变异过程本质是从参与变异的个体的邻域解集内随机选择一个解作为变异后的个体，而这一变异过程往往具有随机性。与传统的变异过程相比，TSSAM引入禁忌搜索算法中的邻域解集，在邻域解集中寻找可行解。如果变异个体的邻域解集中的最优解优于此个体，那么接受最优解为变异后的个体；反之，则通过Metropolis法则以一定概率接受未被禁忌的最优解作为变异后的个体。

综上分析，TSSAM兼具禁忌搜索算法与模拟退火算法的优点，既能够扩大解的搜索范围，又能够利用Metropolis法则提升算法的寻优能力，极大改善了遗传算法的变异过程，加快了算法的收敛过程。

本文设计的TSSAM算法的流程图如图 8所示。

3 典型AEOS任务规划问题验证

本文通过仿真验证了禁忌退火遗传混合算法提出的禁忌退火遗传混合算法在AEOS任务规划问题中的有效性。

3.1 仿真参数设置

考虑AEOS如下参数, 见表 1。

通过STK建立上述卫星的场景，在场景中随机生成若干个卫星可见的目标点，并为这些目标点随机指定1~10的任务优先级。这些目标点的范围为20°N~50°N，110°E~130°E。仿真时间从2020年3月24日04∶ 00∶ 00开始，到2020年3月25日04∶ 00∶ 00结束。通过STK可以解算出卫星对目标点可见的时间窗口。

为了充分验证算法的有效性，将禁忌退火遗传混合算法(TSSAGA)分别与禁忌遗传算法(TSGA)、退火遗传算法(SAGA)以及普通的遗传算法(GA)进行对比仿真实验。

遗传算法的求解性能与种群规模M、父本数量以及变异率有关。工程上通常取M=10N，父本数量为0.3M，变异率为0.09。禁忌搜索算法的寻优性能取决于禁忌长度以及邻域解集的规模，工程上通常取禁忌长度为$\sqrt{N(N-1) / 2}$，邻域解集的规模取为N/3。模拟退火算法的寻优性能取决于算法的初始温度与退火速率。初始温度与待规划的目标点数量有关，本文取退火速率为K=0.9。

综上所述，目前可确定的各算法参数见表 2。

3.2 仿真结果分析

本文设置了2组不同规模的对比仿真实验来验证TSSAGA在AEOS任务规划中的性能。第1组实验为从50个目标点中选20个目标点进行观测，即染色体长度N=20。上文中未给出的各个算法其他参数见表 3。

第2组实验为从100个目标点中选50个进行观测，染色体长度N=50。各算法其他参数见表 4。

为了尽可能详尽全面地比较算法的求解能力，每组对比实验中，每种算法进行50次蒙特卡洛仿真实验，随后比较算法的平均适应度函数值以及算法达到收敛的总耗时。算法收敛的判断标准为，连续15次迭代的最优解相同。对50个目标点进行规划的对比试验结果见表 5。

为了进一步验证算法的有效性，将每种算法分别运算10组，每组进行50次蒙特卡洛仿真实验并取平均值，对比试验结果如图 9所示。最优的观测序列见表 6。对100个点进行规划的对比实验结果见表 7。同样地，将每种算法分别运算10组，每组进行50次蒙特卡洛仿真实验并取平均值，对比试验结果如图 10所示。最优的观测序列见表 8。

综合表 5与表 7可知，求解相同规模的任务规划问题，TSSAGA收敛最快，且收益高于其他算法。相比于GA，TSAG和SAGA节省了20%~30%的运算时间，TSSAGA节省了约40%的运算时间。相比于传统遗传算法的变异过程，TSSAGA在变异过程中引入了禁忌搜索的方法以及Metropolis法则，极大地提升了算法的搜索效率以及寻优能力，节省了算法的收敛时间，有较高的工程应用价值。

4 结论

1) 面向敏捷对地观测卫星观测大规模地面点目标这一实际工程问题，考虑了卫星面临的多种约束条件，建立了对应的适应度函数。

2) 改进了传统遗传算法的变异过程，提出了禁忌退火变异方法，此方法兼具了禁忌搜索算法与模拟退火算法的有点，提高了整个算法的寻优能力，加快了算法的收敛速度。

3) 仿真结果表明，禁忌退火遗传混合算法与传统的遗传算法相比，节省了约40%的算法优化时间，大幅提高了智能优化算法求解敏捷对地观测卫星任务规划问题的求解效率。