２０１０矩 ６月 装备指挥技术学院学报 Ｊｏｕｒｎａｌ ｏｆ ｔｈｅ Ａｃａｄｅｍｙ ｏｆ Ｅｑｕｉｐｍｅｎｔ Ｃｏｍｍａｎｄ＆Ｔｅｃｈｎｏｌｏｇｙ Ｊｕｎｅ ２０１０ Ｖ０Ｉ．２１ Ｎｏ．３ 第２１卷 第３期 基于模拟退火算法的电子侦察卫星 任务规划问题研究 徐 欢， 祝江汉， 王慧林 （国防科技大学信息系统与管理学院，湖南长沙４１００７３） 摘 要：在分析电子侦察卫星的工作原理及其所担负的使命任务的基础 上，首先给出了电子侦察卫星任务规划的基本过程；然后在明确电子侦察卫星任务规 划问题的基本输入输出和一些基本假定的基础上，建立电子侦察卫星静态任务规划 的混合整数规划模型；最后提出了一种解决电子侦察卫星任务规划问题的模拟退火 算法（ｓｉｍｕｌａｔｅｄ ａｎｎｅａｌｉｎｇ，ＳＡ）。实验结果表明：该算法有效地解决了针对固定目标 的电子侦察卫星任务规划问题。 关 键 词：电子侦察卫星；任务规划；模拟退火算法 中图分类号：ＴＰ １８ 文献标识码：Ａ 文章编号：１６７３—０１２７（２０１０）０３－００６２－０５ ＤＯＩ：１０．３７８３／ｊ．ｉｓｓｎ．１６７３—０１２７．２０１０．０３．０１５ Ｒｅｓｅａｒｃｈ ｏｎ Ｍｉｓｓｉｏｎ Ｐｌａｎｎｉｎｇ ｆｏｒ ＥｌｅｃｔｒｏｎｉＣ ＲｅｃＯｎｎａｉｓｓａｎＣｅ Ｓａｔｅ…ｔｅｓ Ｂａｓｅｄ ０１３ Ｓｉｍｕｌａｔｅｄ ＡｎｎｅａＩ ｉｎｇ ＸＵ Ｈｕａｎ， ＺＨＵ Ｊｉａｎｇｈａｎ， ＷＡＮＧ Ｈｕｉｌｉｎ （Ｃｏｌｌｅｇｅ ｏｆ Ｉｎｆｏｒｍａｔｉｏｎ Ｓｙｓｔｅｍ ａｎｄ Ｍａｎａｇｅｍｅｎｔ，Ｎａｔｉｏｎａｌ Ｕｎｉｖｅｒｓｉｔｙ ｏｆ Ｄｅｆｅｎｓｅ Ｔｅｃｈｎｏｌｏｇｙ，Ｃｈａｎｇｓｈａ Ｈｕｎａｎ ４ １ ００７３，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｂａｓｅｄ ｏｎ ｔｈｅ ａｎａｌｙｓｉｓ ｏｆ ｔｈｅ ｗｏｒｋｉｎｇ ｐｒｉｎｃｉｐｌｅ ａｎｄ ｍｉｓｓｉｏｎｓ ｏｆ ｅｌｅｃｔｒｏｎｉｃ ｒｅｃｏｎｎａｉｓｓａｎｃｅ ｓａｔｅｌｌｉｔｅｓ，ｔｈｅ ｐａｐｅｒ，ｆｉｒｓｔｌｙ ｐｕｔｓ ｆｏｒｗａｒｄ ｔｈｅ ｂａｓｉｃ ｐｒｏｃｅｓｓ ｏｆ ｍｉｓｓｉｏｎ ｐｌａｎｎｉｎｇ ｆｏｒ ｅｌｅｃｔｒｏｎｉｃ ｒｅｃｏｎ－ ｎａｉｓｓａｎｃｅ ｓａｔｅｌｌｉｔｅｓ，ａｎｄ ｔｈｅｎ ｂａｓｅｄ ｏｎ ｔｈｅ ｄｅｆｉｎｉｔｉｏｎ ｏｆ ｔｈｅ ｂａｓｉｃ ｉｎｐｕｔ ａｎｄ ｏｕｔｐｕｔ ｏｆ ｔｈｅ ｐｒｏｂｌｅｍ ａｎｄ ｂａｓｉｃ ｈｙｐｏｔｈｅｓｉｓ，ｔｈｅ ｐａｐｅｒ ｅｓｔａｂｌｉｓｈｅｓ ｔｈｅ ｓｔａｔｉｃ ｍｉｓｓｉｏｎ ｐｌａｎｎｉｎｇ ｍｏｄｅ１．Ａ ｓｉｍｕｌａｔｅｄ ａｎｎｅａｌｉｎｇ（ＳＡ） ｉｓ ｄｅｓｉｇｎｅｄ ｆｏｒ ｔｈｅ ｍｉｓｓｉｏｎ ｐｌａｎｎｉｎｇ ｆｏｒ ｅｌｅｃｔｒｏｎｉｃ ｒｅｃｏｎｎａｉｓｓａｎｃｅ ｓａｔｅｌｌｉｔｅｓ ａｔ ｌａｓｔ．Ｔｈｅ ｒｅｓｕｌｔｓ ｓｈｏｗ ｔｈａｔ ｔｈｅ ｐｒｏｂｌｅｍ ｏｆ ｍｉｓｓｉｏｎ ｐｌａｎｎｉｎｇ ｆｏｒ ｅｌｅｃｔｒｏｎｉｃ ｒｅｃｏｎｎａｉｓｓａｎｃｅ ｓａｔｅｌｌｉｔｅｓ ａｉｍｅｄ ａｔ ｆｉｘｅｄ ｔａｒｇｅｔｓ ｃａｎ ｂｅ ｓｏｌｖｅｄ ｅｆｆｅｃｔｉｖｅｌｙ ｗｉｔｈ ｔｈｉｓ ｍｅｔｈｏｄ． Ｋｅｙ ｗｏｒｄｓ：ｅｌｅｃｔｒｏｎｉｃ ｒｅｃｏｎｎａｉｓｓａｎｃｅ ｓａｔｅｌｌｉｔｅｓ；ｍｉｓｓｉｏｎ ｐｌａｎｎｉｎｇ；ｓｉｍｕｌａｔｅｄ ａｎｎｅａｌｉｎｇ 电子侦察卫星是用于截获和监听敌方电子设 相关的防空系统、武器系统的配置情况；侦收、分 析敌方的遥控、遥测信号，掌握其战略武器系统的 性能、试验情况和发展动向；截获、监听敌方的无 备、重要设施和军事目标所辐射的电磁信号，以获 取战略或战术情报的人造地球卫星。它是随着电 子战和空间技术发展起来的一种现代电子侦察手 线电通信，分析信号特征并测定辐射源的位置，解 段，其主要任务是：侦收敌方的雷达信号，测定其 战术技术参数、位置，判明其类型、用途以及与之 收稿日期：２００９—０４—１７ 读破译其通信内容，从中分析掌握敌方潜在的军 事行动和作战计划；长期监视敌电磁辐射源的变 作者简介：徐欢，男，硕士研究生．主要研究方向：作战模型与模拟．ｘｈｓｔｒｏｎｇ＠ｙａｈｏｏ．ｅｏｒａ．ｃｎ． 祝江汉，男，副教授，硕士生导师． 第３期 徐欢，等：基于模拟退火算法的电子侦察卫星任务规划问题研究 ６３ 化情况，获取其电子设备发展水平以及部队配置 和活动规律等情报。 卫星将要执行的侦察任务，并确定了相应的数据 采集和数据下传活动的起止时间。航天侦察计划 电子侦察卫星按侦察任务的不同可分为：普 查型电子侦察卫星和详查型电子侦察卫星，一般 部署在高、中、低３类轨道上。 高轨（ｈｉｇｈ ｅａｒｔｈ ｏｒｂｉｔ，ＨＥＯ）包括地球同步 轨道和大椭圆轨道，卫星可以对地球的某一地区 的制订直接源于对用户的航天侦察需求进行任务 规划的结果。任务规划是指为了实现给定的目 标，对所有与特定任务相关的资源和活动进行规 划与调度的过程。电子侦察任务规划具体指：在 综合考虑卫星平台、有效载荷、地面接收处理系统 进行长时间的连续监视，所以主要用于窃听敌方 的能力和电子侦察任务要求的基础上，将资源分 配给相互竞争的多个任务，并确定任务中各具体 活动的起止时间，以排除不同任务之间的资源使 通信信号、监视导弹发射的测控信号和侦察远程 预警雷达信号，由于轨道很高，卫星对目标的定位 能力有限。 中轨（ｍｅｄｉｕｍ ｅａｒｔｈ ｏｒｂｉｔ，ＭＥｏ）是指６００～ ９００ ｋｍ的近圆轨道，轨道倾角在９Ｏ。左右，又称 “极轨”或太阳同步轨道，每天通过地球上某地的 时间不变。中轨卫星具有较长的寿命和中等定位 精度，一般用于实现普查和详查功能。 低轨（１ｏｗ ｅａｒｔｈ ｏｒｂｉｔ，ＬＥｏ）一般用于详查型 卫星，轨道选择高度为１７５～４００ ｋｍ的椭圆轨 道，轨道倾角为７０。左右，轨道周期约为９０ ｍｉｎ， 每周可以覆盖全球１次。低轨电子侦察卫星作用 距离较近，可以获取较低辐射功率的目标信号，同 时获得较好的定位精度，但由于轨道较低，卫星在 轨工作寿命较短。 在轨的１颗或多颗电子侦察卫星和地面上的 １个或多个数据接收站（可包含天上的数据中继 卫星）组成了１个卫星电子侦察系统。在执行侦 察任务时，电子侦察卫星将利用其各种星载遥感 器来采集敌重要目标的电子情报，并将捕获的情 报数据直接或经中继卫星传回地面接收站，以供 分析判读使用。 在现代战场上，电磁信号几乎覆盖了从短波 一直到微波、毫米波、红外、激光等电磁频谱。现 代通信系统基本上完成了从模拟传输体制向数字 传输的过渡。现代雷达在同一频段内，不仅有常 规的脉冲雷达、连续波雷达、相干雷达和相控阵雷 达，还衍生出多频率、捷变频、脉冲重复频率参差、 抖动与捷变、线性调频脉冲以及频谱扩展等技术 的新体制雷达信号。这样复杂的电磁环境要求侦 察卫星的天线必须具有很高的增益，星载侦察接 收机必须具有极宽的工作频率范围，极低的噪声 性能（提高接收灵敏度）和足够宽的动态范围，具 备在空域、频域和时域内稀释、分选信号的 能力 。 一般情况下，电子侦察卫星是按照预定的侦 察计划来实施侦察的。侦察计划中指定了具体的 用冲突，消除资源使用冗余，最大化资源使用效 率，并满足用户的需求。 卫星电子侦察系统的任务规划问题是一个非 常复杂的问题，涉及的资源和活动很多，且它们之 间的关系又十分复杂。为了科学合理地解决这一 问题，本文在分析卫星电子侦察系统的工作原理 和过程的基础上，进行合理的假设和简化，并建立 了问题的数学模型，然后采用模拟退火算法进行 求解。 １ 电子侦察卫星任务规划问题分析 １．１基本假设与简化 任何实际问题的建模求解过程都需要对客观 世界进行一定的简化，这里只能抓住我们感兴趣 的对象和它们之间的关系，而不可能完全真实地 再现客观世界的复杂内在关联。对于电子侦察卫 星任务规划问题来说，由于涉及的对象较多，对象 之间的关系复杂，因此不可能考虑到全部的细节 和实际约束。本文在研究电子侦察卫星任务规划 问题时，作了如下假设和简化方式：①侦察目标 为固定单电磁目标或地理上分布在比较小区域的 多个固定单电磁目标；②已经获取了固定目标的 概略电磁特征参数；③通过对平时侦察结果的分 析统计，将单个电子目标的侦察持续时间设定为 固定值；④星载电子侦察设备的校准或调整转换 时间可以设为固定值；⑤将每个侦察任务的存储 需求设定为固定值；⑥不考虑卫星的能量限制， 卫星具有足够的能量满足所有的观测任务和通信 任务。 １．２ 电子侦察卫星任务规划流程 卫星电子侦察系统是一个集信息截获、多重 信息传输和处理、涉及天上、地面、星际网的复杂 系统，侦察获取信息多，系统中的信息关系和数据 处理复杂，卫星系统和地面系统的相互依赖程度 高。所以，必须在地面管理中心的统一协调下，系 ６４ 装备指挥技术学院学报 ２０１０年 统才能同步、协调、互联和可靠地运行。卫星电子 侦察任务的组织实施实际是一个从各军兵种用户 提交需求到最终获得情报产品的循环过程。主要 ３）变量： 所占用的资源。 活动Ｊ的后续活动；Ｊ 活动Ｊ 如果活动ｉ与活动Ｊ在执行中占用同一个资 源，且活动Ｊ紧跟着活动ｉ执行，则Ｊ… ＝ｉ；相反 如果ｉ 一 ，则可推出ｉ 一 。 包含以下过程：①由各应用任务提交需求；②初 步检查需求的合法性并确定它们的优先级别； ③分析应用任务需求，通过任务细化产生电子侦 察任务订单；④根据时域、空域、频域覆盖等基本 约束，对电子侦察任务订单中的任务进行预处理， 删除任务订单中不满足执行基本条件的任务； ⑤对电子侦察任务订单中可行的任务需求进行 ｔ 为任务ｉ的开始时间，ｉ∈ＩＵＪ。 ２．２约束条件及优化目标表达 下面给出模型的约束条件及优化目标表达： ｆｌ ｍａｘ＞ｌｃＩ （１） 任务规划，消除对卫星电子侦察资源的使用冲突， 确定相关的数据采集执行时间、传感器工作模式、 工作参数以及数据下传活动的执行时间；⑥将任 务规划结果提交给测控部门，并生成对电子卫星 有效载荷的遥控指令；⑦由地面接收网络按计划 完成对侦察数据的接收（接收站与测控站可以是 合并的）；⑧对后续的数据处理以及数据分发等 地面支持活动进行规划调度；⑨根据调度结果完 成数据产品的存储整编和分发；⑩用户接收到订 购的产品。 电子侦察卫星任务规划过程如图１所示。 一信息 需求 （丫用户 ｆ ｆ ｆ ．塑 笪塑鱼．Ｉ电子侦察茬ｌ ｌ 翌Ｈ ｌ 侦察二 活动计划ｉ ４任务订单务豳ｌ 数传活动计划 图１ 电子侦察卫星任务规戈Ｉｌ过程 ２ 电子侦察卫星任务规划问题模型 ２．１相关参数定义 为了建立电子侦察卫星任务规划的数学模 型，首先定义相关符号如下。 １）集合：Ｒ为卫星资源集； 为侦察元任务 集；Ｊ为一个所有卫星所有可能的数据下传任务。 ２）参数：ｔ 为资源五上活动ｉ的持续时间； ｔｆ． 为资源ｋ上任务ｉ和任务 之问的转换时间； 为资源ｋ在执行活动ｉ之前已被占据的存储 空间；Ｃｋ为资源忌的存储容量；ｑ 为资源忌执行 任务ｉ时的存储需求； 为资源忌数据下传速率。 如果ｉ ＝＝： ∈Ｉ～ｉ。一是，则 ： ＋ ；如果 ｉ 一 ∈Ｊ，ｉ 一是，贝 Ｑ 一Ｑ 一 ×ｔ 。 ｓ？为活动ｉ的时间窗口开始时间；ｅ？为活动ｉ 的时间窗口结束时间。 ｓ ｔ Ｊ ｔ ≥ｔ ＋￡ ＋ （２） ｌ ０≤ ≤Ｃ （３） 【５７ ≤ｔ≤Ｐ？ —ｔ （４） 目标函数式（１）保证了安排活动的总权值最 大，其中Ｃ 为活动ｉ的权值，“ 为布尔变量。约束 条件式（２）保证了每个资源在任何时候只能满足 一个任务的需求。约束条件式（３）保证了所有卫 星在任何时刻存储的数据总量都不能超过其存储 容量的限制。约束条件式（４）保证了所有任务必 须在可选资源的时间窗口内执行，其中 为任务 ｉ的时间窗口总数。 ３ 电子侦察卫星任务规划模型求解 ３．１模拟退火算法 模拟退火算法是由Ｍｅｔｒｏｐｏｌｉｓ【２ 等人１９５３ 年提出的，并由ＫｉｒｋｐａｔｒｉｃｋＩ３ 在１９８３年首次成 功应用到组合优化问题中。ＳＡ算法是一种基于 Ｍｏｎｔｅ—Ｃａｒｌｏ迭代求解策略的随机搜索算法，它 源于对固体退火过程的模拟，采用Ｍｅｔｒｏｐｏｌｉｓ接 受准则，并用一组称为冷却进度表的参数控制算 法进程，使算法在多项式时间内给出近似最优解。 该算法具有描述简单、适用范围广、运行效率高、 受初始条件限制小等优点，已经被成功地应用于 各类组合优化问题。 ３．２模拟退火算法基本应用设计 ３．２．１ 邻域结构 本文采用两交换法实施邻域操作，该方法是 指随机选择解中的２个元素，并交换其值的邻域 操作方法。 ３．２．２候选解取舍准则 本文所建立的任务规划模型是一个最大化问 题，即最大化成功调度的侦察元任务之总权值，可 按照Ｍｅｔｒｏｐｏｌｉｓ准则对候选解进行取舍： ｆ１， ｆ（Ｓ ）≥厂（Ｓ） ，）一１． ｌｅｘｐ —————— —一・，ｆ （Ｓ，）Ｊ＜ ）＜ｒＩ （５） 第３期 徐欢，等：基于模拟退火算法的电子侦察卫星任务规划问题研究 ６５ ３．２．３温度参数控制 １）初始温度ｔ。。初始温度应选得足够高，才 能不使算法过早陷入局部最优解；但又不能选得 太高，以减少冗余的迭代。根据本文所研究问题 的特点，设置初始温度ｔ。为：ｔ。一＞：Ｗ ，即所有 ｉｆｆ 侦察元任务的总权值。 ２）控制参数衰减函数△（ ）。常用的控制参 数下降函数有２种形式：①等比例下降，ｔ ＋ 一 ａｔ ，尼≥Ｏ，其中，降温系数ａ∈（０，１）；②等步长下 Ｌ，一Ｌ 降，ｔ 一 』、 。，其中，ｔ。为起始温度，Ｋ为算法温 度下降的总次数。 为避免过早的陷入局部最优解，本文选用等 比例下降形式的控制参数下降函数，降温系数ａ 根据多次测试结果，设置为０．９。 ３）算法终止规则。模拟退火算法的终止规 则很多＿４ｊ，常用的有以下几种：①零度法，即当温 度ｔ 下降至小于某小正数ｅ时，算法停止；②循 环总数控制法，即总的温度下降次数为一定值 Ｋ，当温度迭代次数达到Ｋ时，算法停止；③不改 进规则控制法，在某一温度，给定次数内没有改进 当前的局部最优解，算法停止。 在模拟退火算法中，设置了同时使用上述３ 种搜索算法终止规则，只要满足其中任意一种终 止规则，即停止搜索过程。 ４）Ｍａｒｋｏｖ链长度Ｌ 。记温度ｔ 对应的 Ｍａｒｋｏｖ链长度为Ｌ ，定义了长度随温度变化的 规则：Ｌ ＋ｌ：＝＝１．０２×Ｌ ，ｋ＝＝＝０，１，…，Ｎ。 这里初始温度ｔ。所对应的Ｍａｒｋｏｖ链长度 Ｌ。按如下规则取值：Ｌ。一 ７ｌ。ｎ是候选观测活动 厶 的总数量。 ３．３模拟退火算法的基本流程 在算法基本应用设计的基础上，给出模拟退 火算法寻优阶段的基本流程如下。 Ｓｔｅｐ １：得到初始解Ｓ。，令最优解记录为 Ｓ 一Ｓ。。 Ｓｔｅｐ ２：设置初始温度ｔ。，最终温度ｔ，初始 Ｍａｒｋｏｖ链长度Ｌ。，最大连续不更新次数为 Ｋｍａ 一ｄ，令￡一ｔ。，Ｌ—Ｌ。，不更新次数ｎＲ一０。 Ｓｔｅｐ ３：若ｔ＞ｔ，转Ｓｔｅｐ ４，否则转Ｓｔｅｐ ｌ１。 Ｓｔｅｐ ４：若ｎＲ＜ｄ，转Ｓｔｅｐ ５，否则转Ｓｔｅｐ １１。 Ｓｔｅｐ ５：若‰ ≤Ｌ，转Ｓｔｅｐ ６，否则转Ｓｔｅｐ ８。 Ｓｔｅｐ ６：采用随机抽样方式从当前解Ｓ的邻 域空间内取出一个候选解ｓ ，并计算候选解Ｓ 和 当前解Ｓ的目标函数之差Ａｆ，Ａｆ一 厂（Ｓ）－－ｆ（Ｓ ）。 Ｓｔｅｐ ７：如果候选解５ 和当前解Ｓ的目标函 数之差Ａｆ＜０，接受该候选解；如果Ａｆ＞０，仅当 条件ｅｘｐ（Ａｆ／ｒ）＞ｒａｎ （０，１）被满足时接受候选 解Ｓ ，否则拒绝候选解Ｓ ，此时无论是否接受候 选解Ｓ 都令计数器 ＝＝＝ ＋１，并转Ｓｔｅｐ５。 Ｓｔｅｐ ８：如果当前温度所对应的最优解同之 前温度所对应的最优解相比没有提高，则令不更 新次数 Ｒ一，ｚＲ＋１，否则 Ｒ一０。 Ｓｔｅｐ ９：若当前解Ｓ大于最优解记录 Ｓ ＝＝＝Ｓ。。 Ｓｔｅｐ １０：根据温度衰减函数△（ｆ）降低温度， 使￡一△（￡），并更新当前温度ｔ所对应的Ｍａｒｋｏｖ 链长度Ｌ，令 一０，转Ｓｔｅｐ ３。 Ｓｔｅｐ １１：输出最优解Ｓ 。 ４计算实例分析 本文用若干个随机生成的不同规模的测试算 例对ＳＡ算法进行验证。测试算例４个，规划周 期均为１ ｄ，分别为１星１０目标、２星２０目标、３ 星３０目标、４星４０目标。本文所有算例都是在 Ｐｅｎｔｉｕｍ（Ｒ）Ⅳ２．４ ＧＨｚ、５１２ ＭＢ ＲＡＭ的台式 机上运行。 为了使本文设计的电子侦察卫星系统资源更 加切合实际，从国外已投入使用的电子侦察卫星 中选取了４颗具有代表性的卫星（俄罗斯Ｃｏｓｍｏｓ 系列电子侦察卫星：Ｃｏｓｍｏｓ一２０８２、Ｃｏｓｍｏｓ一２１５１、 Ｃｏｓｍｏｓ一１５４４、Ｃｏｓｍｏｓ一２０５８）构成一个电子侦察 卫星网，假设每颗卫星携带１部电子侦察载荷，传 感器采用简单圆锥类型。然后在某热点区域 （３２．５。Ｅ～３４．５。Ｅ，４３。Ｎ－－４６．５。Ｎ）内随机地设置 ４Ｏ个电子侦察目标。 设仿真时间为２００８—０１—０１Ｔ００：００：００／２６０８— ０１—０２Ｔ００：００：００，在此时间段根据如下假设生成 ４０个电子侦察元任务： １）假设电子侦察元任务的执行时间最小值 为１ ｒａｉｎ，最大值为３ ｒａｉｎ。每颗卫星执行任意２ 个任务之间的工作模式转换时间统一设为 ５ ｍｉｎ； ２）任务优先级随机设定为１～１Ｏ之间的１ 个整数，数值越大表示任务优先级越高。 当目标数为４０个时，从图２所示的算法的收 敛过程可以看出，算法具有较好的收敛性。 装备指挥技术学院学报 ２０１０年 — …， 一 。。 躜 一 。１ ｔ． ’ … ‘ ｌ 图２进化收敛不恿图 下面重点分析模拟退火算法在上述４种具有 不同问题规模计算实例中的表现。每个实例运算 １Ｏ次，取算法１ｏ次运算的均值作为评价基础，主 要统计算法的模型求解效果和运行时问，结果如 表１所示。 表１ ＳＡ算法结果 问题规模 磬 从图２可以看出，经过ＳＡ算法优化的侦察 计划执行总收益非常高；从表１可看出，当问题规 模较小、资源相对丰富的情况下，模拟退火算法表 现较为优异。但随着问题规模的增大，资源冲突 凸显出来，算法的运行时间快速增加，这是由算法 的复杂度决定的，但侦察计划的执行总收益还是 比较高的，这说明ＳＡ算法是具有相当的优化效 果的。 ５ 结束语 对于电子侦察卫星任务规划问题，本文研究 了一种基于模拟退火算法的任务规划算法。在算 法设计前分析了电子侦察卫星任务规划问题的关 键约束，总结出电子侦察卫星任务规划问题的特 点，并在合理假设的基础上建立了规划模型。该 算法针对问题特点，选用侦察元任务的总权值作 为退火的初始温度，并设计了相应的控制参数衰 减函数及终止条件。最后设计了一个多星多目标 场景算例，并生成了一批电子侦察任务，然后分别 在不同的问题规模下，采用模拟退火算法进行求 解。结果表明：在一定资源冲突下，该算法能够取 得较好的规划结果。 参考文献 （Ｒｅｆｅｒｅｎｃｅｓ） Ｅ１］林象平．雷达对抗原理［Ｍ］．西安：西北电讯工程学院出版 社，１９８５：１３２ １６１． ［２］ＭＥＴＲＯＰＯＩ ＩＳ Ｎ，ＲＯＳＥＮＢＩ ＵＴＨ Ａ，ＲＯＳＥＮＢＬＵＴＨ Ｍ． Ｅｑｕａｔｉｏｎ ｏｆ ｓｔａｔｅ ｃａｌｃｕｌａｔｉｏｎｓ ｂｙ ｆａｓｔ ｃｏｍｐｕｔｉｎｇ ｍａｃｈｉｎｅｓ ［Ｊ］．Ｊｏｕｒｎａｌ ｏｆ Ｃｈｅｍｉｃａｌ Ｐｈｙｓｉｃｓ，１９５３，２１（６）：１０８７—１０９２． ［３］ＫＩＲＫＰＡＴＲＩＣＫ Ｓ，ＧＥＩ ＡＴＴ ＪＲ Ｃ Ｄ，ＶＥＣＣＨＩ Ｍ Ｐ．Ｏｐｔｉ— ｍｉｚａｔｉｏｎ ｂｙ ｓｉｍｕｌａｔｅｄ ａｎｎｅａｌｉｎｇ［Ｊ］．Ｓｃｉｅｎｃｅ，１９８３，２２０（８）： ６７１ ６８０． ［４］邢文训，谢金星．现代优化计算方法［Ｍ］．北京：清华大学出 版社。２００５：１０３． （编辑：孙陆青）