CN109087023B - 一种基于分治策略的多星观测分层调度方法及*** - Google Patents
一种基于分治策略的多星观测分层调度方法及*** Download PDFInfo
- Publication number
- CN109087023B CN109087023B CN201810967096.7A CN201810967096A CN109087023B CN 109087023 B CN109087023 B CN 109087023B CN 201810967096 A CN201810967096 A CN 201810967096A CN 109087023 B CN109087023 B CN 109087023B
- Authority
- CN
- China
- Prior art keywords
- task
- scheduling
- tasks
- scheduled
- scheduling sequence
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Economics (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Marketing (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种基于分治策略的多星观测分层调度方法及***,具体包括以下步骤:S1、利用分配算法把任务分配至各轨道圈次上,形成各轨道圈次的任务集;S2、利用近似优化算法求解各轨道圈次上任务集的调度序列;S3、分配算法根据各轨道圈次调度序列反馈的未生成调度序列的任务重新更新分配算法的分配方案,进而形成新的各轨道圈次的新的任务集;S4、重复步骤S1、S2、S3直至达到多星观测分层调度的终止条件。通过把一个复杂的组合优化问题进行有效的分解简化,变成一个双层规划,有效降低了问题求解的复杂度,尤其在求解大规模多星观测调度问题时表现出卓越的性能。本发明应用于卫星调度技术领域。
Description
技术领域
本发明涉及卫星调度技术领域,尤其涉及一种基于分治策略的多星观测分层调度方法及***。
背景技术
对地观测卫星(Earth Observing Satellites,EOSs)绕地球飞行,能够获取地球表面指定区域的图像信息以满足用户的观测需求,在环境监控、情报侦察等领域扮演了关键角色。为了充分利用稀缺的卫星资源,研究高效的观测调度方法以提高卫星应用水平,具有重要的意义。由于用户的对地观测需求通常超过现有卫星资源的观测能力,因此卫星观测调度问题(Satellite Observation Scheduling Problem,SOSP)是一类过度订阅(oversubscribed)问题。
早期与SOSP相关的研究侧重于单颗卫星的调度。Vasquez和Hao将卫星观测调度问题转化为经典背包模型,提出了一种禁忌搜索算法求解该模型。在另一篇论文中,Vasquez和Hao提出了基于分解的方法来获得优化问题最小上界,以评估禁忌搜索算法求解所得卫星观测调度方案的质量。Bensana等将Spot5卫星的调度转化为约束满足问题,并利用精确或近似算法分别对该问题进行求解。Gabrel和Vanderpooten提出了一种无环图模型来描述Spot5卫星的调度问题,首先生成多条有效路径,然后从中选定一条最好的路径。Lin等采用数学规划方法来获得近似最优卫星观测调度方案。Baek等采用一种新的遗传算法模拟卫星实际观测调度问题。Miguel等基于在节点打包和3-正则独立***的有效不等式,为Spot5卫星的日常成像调度提供规划方案。
卫星观测调度是要在满足一系列观测任务需求和约束条件下,确定每颗卫星在每个轨道圈次上的观测活动,实现观测收益的最大化,是提高卫星使用效率的关键。传统的多星调度算法以每颗卫星为调度单位,但是如果调度时间区间较大,一颗卫星会绕地球飞行多个轨道圈次(例如成像侦察卫星一般每天绕地飞行14轨道圈次左右),不同的卫星可能观测到相同的目标,且同一颗卫星在其不同的轨道圈次上可能可以观测到相同的目标,造成观测冲突更加复杂(一般假定一个目标只需观测一次),尤其在大规模多星观测调度问题时的求解程度更加复杂。
发明内容
针对现有技术存在的不足,本发明的目的是提供一种基于分治策略的多星观测分层调度方法及***,该方法有效降低了问题求解的复杂度,尤其在求解大规模多星观测调度问题时表现出卓越的性能。
本发明在进行多星观测分层调度时,首先将卫星分解到轨道圈次,把每个轨道圈次看成是具备对地观测能力的资源,即将多星多轨道圈次调度问题转化为多轨道圈次调度问题。
设集合O中包含一般性的具有对地观测能力轨道圈次资源,O={oj|j=1,2,...,M},M为集合的基,即是调度问题中轨道圈次资源的数量。设集合T是待观测的目标集合,包括点目标和分解后的条带目标,T={ti|i=1,2,...,N},N表示任务数量。0-1变量xij用来表示在调度方案中任务ti是否在轨道圈次oj上完成:任务ti在轨道圈次oj上完成时xij=1,否则xij=0。多星观测调度的目标是最大化观测收益,为所有完成任务的收益之和,即:
式中,pi是完成任务ti可获得的收益。
调度过程中需要满足复杂的约束条件,包括唯一性约束、能量约束、存储容量约束和测摆时间约束,具体如下:
唯一性约束:每个任务只需最多完成一次,即:
能量约束:轨道圈次oj上,观测活动和传感器侧摆活动消耗的能量不能多于该该轨道圈次上允许消耗的最大能量:
式中,eoj表示在轨道圈次oj上卫星和传感器观测活动消耗能量的速率;tei表示任务ti的结束时间;tsi表示任务ti的开始时间;yih是1-0变量,表示ti是否先于th完成,先于则取1,否则取0;esj表示在轨道圈次oj上传感器侧摆活动消耗能量的速率;θi表示在轨道oj上观测任务ti时传感器所需测摆的角度;θh表示在轨道oj上观测任务th时传感器所需测摆的角度;Ej表示轨道圈次oj的允许消耗的最大储能量。
存储容量约束:在轨道圈次oj上,观测活动获取对地观测数据并进行星上存储,但是消耗的存储资源不能超过该轨道圈次上所能允许的最大存储容量:
式中,wj表示轨道圈次oj上观测活动消耗存储资源的速率;Wj表示轨道圈次oj的允许消耗最大存储容量。
侧摆时间约束:卫星在观测不同目标时,需要重新开机,同时由于相对的空间几何关系发生变化,需要调整传感器的侧摆角度以对准观测目标。因此同一个轨道圈次上,相邻的两个观测任务之间,需要充足的时间让卫星传感器进行开机和校准(调整侧摆角度):
tsh-tei≥aj+|θh-θi|/vj,i,h=1,2,...,N,yih=1,j=1,2,...M
式中,tsh表示任务th的开始时间;aj表示在轨道oj上在对某个目标进行观测之前的传感器开机时间。
假定一个多星调度问题考虑N个任务,M个轨道圈次,则变量数为N·M,约束条件的数量为N+3M。真实场景的卫星调度问题一般都问题规模很大,考虑一般规模的调度问题:例如待观测目标数(包括分解后的观测条带)为500,考虑5颗卫星,调度周期为一天,每颗卫星绕地飞行14圈,则共有大约70个轨道圈次,则变量的规模为500x70=35000,约束条件的规模为500+3x70=710。不难发现多星观测调度问题是一个复杂的离散优化问题,如何快速高效的求解这个问题,需要设计有效的调度框架和算法。
综上,本发明采用的技术方案是:一种基于分治策略的多星观测分层调度方法,包括以下步骤:
S1、利用分配算法把任务分配至各轨道圈次上,形成各轨道圈次的任务集;
S2、利用近似优化算法求解各轨道圈次上任务集的调度序列;
S3、分配算法根据各轨道圈次调度序列反馈的未生成调度序列的任务重新更新分配算法的分配方案,进而形成新的各轨道圈次的新的任务集;
S4、重复步骤S1、S2、S3直至达到多星观测分层调度的终止条件。
作为上述技术方案的进一步改进,步骤S2中,所述近似优化算法为模拟退火算法。
作为上述技术方案的进一步改进,所述模拟退火算法获取当前轨道圈次的初始调度序列的步骤包括:
S211、将单轨道圈次任务集中的待规划调度任务根据权重从大到小依次排列,形成单轨道圈次未被调度的初始任务集U;
S212、从U中选择权重最大的任务tk,应用***和删除邻域结构变换,把tk加入到局部解S中生成新的局部解S′,若S′的适应度高于S,则对S进行更新,令S=S′,否则不对S进行更新,同时将tk从U中移除至Uj中,其中Uj表示在轨道圈次上未被调度的初始任务集;
S213、重复步骤S212直至U中的任务全部遍历完成,获取当前轨道上的调度序列S与当前轨道上未被调度的任务集Uj。
作为上述技术方案的进一步改进,所述模拟退火算法的领域搜索结构包括基于贪婪准则的邻域搜索结构与基于概率搜索的邻域搜索结构,所述模拟退火算法包括邻域搜索结构动态选择策略:根据两种邻域搜索结构在之前分配方案内的优化性能,决定在后续分配方案中领域搜索结构的选择。
作为上述技术方案的进一步改进,基于贪婪准则的邻域搜索结构的模拟退火算法具体包括:
S221、获取当前轨道圈次的调度序列S与当前轨道圈次未被调度的任务集Uj;
S222、在Uj中取出一个尚未调度的权重最大的任务tk,令Uj=Uj-tk;
S223、将tk***到当前轨道圈次的调度序列S中去,判断是否满足侧摆时间约束,若满足则直接形成新的调度序列S″,若不满足则删除与tk冲突的任务后形成新的调度序列S″,其中删除的任务存入Su中,Su是所有轨道圈次未被调度的任务集合;
S224、判断新的调度序列是否满足能量约束和存储容量约束:
若满足,则重复步骤S221、S222、S223、S224直至Uj中所有的任务都遍历完成,获取当前轨道上的最终调度序列S″并将当前轨道未被调度的任务存入集合Su中;
若不满足,则依次删除权重最小的任务直到满足能量约束和存储容量约束后重复步骤S221、S222、S223、S224直至Uj中所有的任务都遍历完成,获取当前轨道上的最终调度序列S″并将当前轨道未被调度的任务存入集合Su中,其中Su是所有轨道圈次未被调度的任务集合。
作为上述技术方案的进一步改进,基于概率搜索的邻域搜索结构的模拟退火算法具体包括:
S231、获取当前轨道圈次的调度序列S与当前轨道圈次未被调度的任务集Uj;
S233、将tk***到当前轨道圈次的调度序列S中去,判断是否满足侧摆时间约束,若满足则直接形成新的调度序列S″,若不满足则删除与tk冲突的任务后形成新的调度序列S″,其中删除的任务存入Su中,Su是所有轨道圈次未被调度的任务集合;
S234、判断新的调度序列是否满足能量约束和存储容量约束:
若满足,则重复步骤S231、S232、S233、S234直至Uj中所有的任务都遍历完成,获取当前轨道上的最终调度序列S″并将当前轨道未被调度的任务存入集合Su中;
若不满足,则删除任务tj的概率为其中indexj为任务tj优先级指标,直到满足能量约束和存储容量约束后重复步骤S231、S232、S233、S234直至Uj中所有的任务都遍历完成,获取当前轨道上的最终调度序列S″并将当前轨道未被调度的任务存入集合Su中,其中Su是所有轨道圈次未被调度的任务集合。
作为上述技术方案的进一步改进,所述分配算法为蚁群优化算法。
作为上述技术方案的进一步改进,步骤S3中具体包括:
S31、根据各轨道圈次调度序列获取所有轨道圈次未被调度的任务集合Su;
S32、从Su中选择一个权重最高的任务,将其***到冲突权重最小的轨道圈次上,若该任务的冲突权重小于该任务自身的权重,则***成功并更新分配方案,否则不更新分配方案;
S33、重复步骤S32直至Su中所有的未被调度的任务都遍历完成。
作为上述技术方案的进一步改进,步骤S4中,所述多星观测分层调度的终止条件为计算时间达到时间预设值或迭代次数达到次数预设值。
本发明还公开了一种基于分治策略的多星观测分层调度***,其采用的技术方案是:
一种基于分治策略的多星观测分层调度***,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现上述方法的步骤。
本发明的有益技术效果:
本发明通过分配算法把任务分到不同的轨道圈次上去,每个轨道圈次对应一个任务集,此为本方法的任务分配层,随后通过近似优化算法得到每个轨道圈次上卫星的观测方案,此为本方法的任务调度层,最终根据各轨道圈次调度序列的反馈重新更新分配算法的分配方案,进而形成新的各轨道圈次的新的任务集。反复迭代直至满足算法终止条件,把一个复杂的组合优化问题进行有效的分解简化,变成一个双层规划,有效降低了问题求解的复杂度,尤其在求解大规模多星观测调度问题时表现出卓越的性能。
附图说明
图1是本实施例的框架结构图。
具体实施方式
为了便于本发明的实施,下面结合具体实例作进一步的说明。
如图1所示的一种基于分治策略的多星观测分层调度方法,包括以下步骤:
S1、利用分配算法把任务分配至各轨道圈次上,形成各轨道圈次的任务集;
S2、利用近似优化算法求解各轨道圈次上任务集的调度序列;
S3、分配算法根据各轨道圈次调度序列反馈的未生成调度序列的任务重新更新分配算法的分配方案,进而形成新的各轨道圈次的新的任务集;
S4、重复步骤S1、S2、S3直至达到多星观测分层调度的终止条件,多星观测分层调度的终止条件为计算时间达到时间预设值或迭代次数达到次数预设值。
步骤S2中,近似优化算法为模拟退火算法,模拟退火算法获取当前轨道圈次的初始调度序列的步骤包括:
S211、将单轨道圈次任务集中的待规划调度任务根据权重从大到小依次排列,形成单轨道圈次未被调度的初始任务集U;
S212、从U中选择权重最大的任务tk,应用***和删除邻域结构变换,把tk加入到局部解S中生成新的局部解S′,若S′的适应度高于S,则对S进行更新,令S=S′,否则不对S进行更新,同时将tk从U中移除至Uj中,其中Uj表示在轨道圈次上未被调度的初始任务集,其中***和删除邻域变换表示先***任务,再删除与新***任务冲突任务的一种邻域搜索方式,局部解是指的有潜力继续加任务的解,适应度衡量解的好坏,这里的解即所需求取的初始调度序列;
S213、重复步骤S212直至U中的任务全部遍历完成,获取当前轨道圈次的调度序列S与当前轨道圈次未被调度的任务集Uj。
本实施例中模拟退火算法的领域搜索结构包括基于贪婪准则的邻域搜索结构与基于概率搜索的邻域搜索结构,模拟退火算法包括邻域搜索结构动态选择策略:根据两种邻域搜索结构在之前分配方案内的优化性能,决定在后续分配方案中领域搜索结构的选择:
在对问题优化求解过程中,不同的邻域结构搜索策略可能具有不同的优势,因此本实施例提出一种自适应机制实现两种不同邻域搜索结构的动态调整,根据每种邻域搜索结构在之前分配方案内的优化性能,即之前给定迭代次数内的优化性能,决定他们在后续迭代过程中被选择的概率,这是一种基于强化学习的思想。假定pro1和pro2分别表示为选择基于贪婪准则的邻域搜索结构和基于概率搜索的邻域搜索结构的概率,在算法初始化时,设定proi=0.5,i=1,2,每隔一定的迭代次数Itr,通过以下规则更新每个邻域结构的选择概率:
proi=proi′/∑i=1,2proi′
式中,proi′表示中间变量;η为惯性权重因子,表示邻域结构之前的选择概率所占的比重;(1-η)表示当前最新历史搜索经验对更新选择概率的权重;Seli为最近Itr次迭代过程中,第i个邻域结构被选中的次数;Suci表示利用第i个邻域结构产生了更高质量的解的次数。
最后对不同邻域结构的选择概率进行标准化。上述邻域结构选择概率的更新规则表明,能够产生更好解的邻域结构,将会得到更大的选择概率。从而实现在优化算法运行过程中邻域结构的动态自适应选择。
在邻域结构搜索中,***/删除邻域是一种经典有效的邻域结构,又通常被称为交换邻域。每次从未被调度且尚未遍历的任务集合中选择权重最大的任务,***到当前的观测调度序列中,如果在调度序列中存在与***任务冲突的任务,则把冲突任务全部删除,因此基于贪婪准则的邻域搜索结构的模拟退火算法具体包括:
S221、获取当前轨道圈次的调度序列S与当前轨道圈次未被调度的任务集Uj;
S222、在Uj中取出一个尚未调度的权重最大的任务tk,tk未包含在禁忌表中,令Uj=Uj-tk,其中tk未包含在禁忌表中表示一种算法机制,即放在禁忌表中的任务不参与邻域结构搜索;
S223、将tk***到当前轨道圈次的调度序列S中去,判断是否满足侧摆时间约束,若满足则直接形成新的调度序列S″,若不满足则删除与tk冲突的任务后形成新的调度序列S″,其中删除的任务存入Su中,Su是所有轨道圈次未被调度的任务集合;
S224、判断新的调度序列是否满足能量约束和存储容量约束:
若满足,则重复步骤S221、S222、S223、S224直至Uj中所有的任务都遍历完成,获取当前轨道上的最终调度序列S″并将当前轨道未被调度的任务存入集合Su中;
若不满足,则依次删除权重最小的任务直到满足能量约束和存储容量约束后重复步骤S221、S222、S223、S224直至Uj中所有的任务都遍历完成,获取当前轨道上的最终调度序列S″并将当前轨道未被调度的任务存入集合Su中,其中Su是所有轨道圈次未被调度的任务集合。
基于贪婪准则的邻域搜索结构有利于局部搜索。为了提高解的多样性,还设计了一种基于概率搜索的邻域搜索结构,在该邻域结构中,不但考虑了任务的权重,还考虑完成每个任务需潜在消耗的侧摆资源和时间窗口资源。综合权重、侧摆角度和时间窗口,在***和替换邻域结构变换过程中,为每个任务tk计算一个优先级指标indexk,首先分别计算权重指标,时间窗口指标和侧摆角度指标:
权重指标为:
式中,iWk为权重指标;wk是第k个任务的权重;wi是第i个任务的权重;nj是Uj任务数量。
时间窗口指标为:
式中,iTk为时间窗口指标;spank为tk的时间窗口长度;spani为ti的时间窗口长度。
侧摆角度指标为:
式中,iTk为侧摆角度指标;θi在轨道oj上观测任务ti时传感器所需测摆的角度;θk在轨道oj上观测任务tk时传感器所需测摆的角度。
优先级指标indexk最终表示为:
indexk=iWk α·iTk β·iθk γ
式中,α,β,γ分别表示不同要素的影响因子,由用户针对具体问题具体设置。
基于概率搜索的邻域搜索结构的模拟退火算法具体包括:
S231、获取当前轨道圈次的调度序列S与当前轨道圈次未被调度的任务集Uj;
S233、将tk***到当前轨道圈次的调度序列S中去,判断是否满足侧摆时间约束,若满足则直接形成新的调度序列S″,若不满足则删除与tk冲突的任务后形成新的调度序列S″,其中删除的任务存入Su中,Su是所有轨道圈次未被调度的任务集合;
S234、判断新的调度序列是否满足能量约束和存储容量约束:
若满足,则重复步骤S231、S232、S233、S234直至Uj中所有的任务都遍历完成,获取当前轨道上的最终调度序列S″并将当前轨道未被调度的任务存入集合Su中;
若不满足,则删除任务tj的概率为其中indexj为任务tj优先级指标,直到满足能量约束和存储容量约束后重复步骤S231、S232、S233、S234直至Uj中所有的任务都遍历完成,获取当前轨道上的最终调度序列S″并将当前轨道未被调度的任务存入集合Su中,其中Su是所有轨道圈次未被调度的任务集合。
分配算法为蚁群优化算法,步骤S3中具体包括:
S31、根据各轨道圈次调度序列获取所有轨道圈次未被调度的任务集合Su;
S32、从Su中选择一个权重最高的任务,将其***到冲突权重最小的轨道圈次上,若该任务的冲突权重小于该任务自身的权重,则***成功并更新分配方案,否则不更新分配方案;
S33、重复步骤S32直至Su中所有的未被调度的任务都遍历完成。
其中,在步骤S1中采用蚁群优化算法第一次把任务分配至各轨道圈次上时,蚁群算法会根据参数设置进行初始分配,直接形成第一次的分配方案。
作为上述技术方案的进一步改进,步骤S4中,算法终止条件为计算时间达到时间预设值或迭代次数达到次数预设值。
本实施例还公开了一种基于分治策略的多星观测分层调度***,其采用的技术方案是:
一种基于分治策略的多星观测分层调度***,包括存储器和处理器,存储器存储有计算机程序,处理器执行计算机程序时实现上述方法的步骤。
为了评估基于本实施例中分治策略的ACO-SA算法的性能,将其与禁忌搜索(TS),遗传算法(GA),模拟退火(SA)最高优先级优先调度算法(HPFS)和带局部搜索的蚁群优化算法进行对比(ACO-LS)。ACO-LS利用蚁群优化算法与局部搜索策略结合,生成多星规划调度问题的解,但是在该算法中没有用到分而治之的策略,而是把调度问题看成一个整体进行求解。
为了全面更加全面的评估算法的性能,本实施例提出了六个仿真场景。所使用卫星包含两组,第一组是8颗资源和尖兵系列侦察卫星,第二组是16颗侦察卫星。每颗卫星绕地球运行一圈的时间约100分钟,一天运行14个轨道圈次。传感器的横向侧摆角度范围为[-330,330]。调度周期是24小时。运用专业软件STK计算得到卫星与目标的可见时间窗口和侧摆角度。六个场景分别是:第一个场景8颗卫星对大兴安岭林区的观测;第二个场景是8颗卫星对长白山林区的观测;第三个场景是8颗卫星对大兴安岭林区和长白山林区的观测;第四个场景16颗卫星对大兴安岭林区的观测;第五个场景是16颗卫星对长白山林区的观测;第六个场景是16颗卫星对大兴安岭林区和长白山林区的观测。
使用C++语言编程实现该算法,并在Intel Core(TM)[email protected] CPU、16.0GB RAM的Windows 7***的电脑上运行实验。
第一组仿真实验场景至第六组仿真实验场景的计算结果如表1到表6中的数据所示。通过对实验数据的观测,不难发现如下几个现象:
1.ACO-SA在每一组实验中都稳定地表现出最佳性能,显示出采用分治策略集成蚁群算法和模拟退火算法可以有效的解决多星观测调度问题。尤其是ACO-SA算法总是比ACO-LS生成更好的观测调度方案,表明分治策略对于求解大规模卫星观测调度问题具有显著优势;
2.在8颗卫星的仿真场景下,对森林区域的观测覆盖率显著小于16颗卫星仿真场景下对森林区域的观测覆盖率,表明适当增加卫星观测资源是提高森林观测覆盖能力的重要手段;
3.通过第三组仿真实验的数据可以发现,8颗卫星24小时内难以有效完成对大兴安岭和长白山地区的有效观测覆盖,采用最高效的ACO-SA算法亦只能达到大概50%的覆盖率。16颗卫星24小时内则基本可以有效观测覆盖长白山地区的森林资源,观测完成率达到97%以上;
4.多星观测规划调度的时效性看,ACO-SA的时效性处在中间水平,ACO-SA生成卫星观测方案的时间随着任务数和卫星数量规模的增长呈近似线性增长,在最大规模(第六组仿真场景)下,能够在900秒内得到高质量的规划调度方案,表明ACO-SA算法的适用性。
表1不同卫星观测调度算法在第一组试验场景的计算结果
算法 | CT | NT | Cov(%) | Time(s) |
TS | 1862 | 1215 | 65.22% | 321.52 |
GA | 1862 | 1244 | 66.52% | 425.62 |
SA | 1862 | 1198 | 62.12% | 564.24 |
ACO-LS | 1862 | 1254 | 67.25% | 368.45 |
ACO-SA | 1862 | 1351 | 72.65% | 388.78 |
表2不同卫星观测调度算法在第二组试验场景的计算结果
算法 | CT | NT | Cov(%) | Time(s) |
TS | 1456 | 1025 | 70.40% | 280.25 |
GA | 1456 | 1085 | 74.52% | 355.32 |
SA | 1456 | 987 | 67.79% | 451.52 |
ACO-LS | 1456 | 1102 | 75.69% | 298.12 |
ACO-SA | 1456 | 1203 | 83.24% | 326.45 |
表3不同卫星观测调度算法在第三组试验场景的计算结果
算法 | CT | NT | Cov(%) | Time(s) |
TS | 3318 | 1425 | 42.95% | 582.35 |
GA | 3318 | 1524 | 45.93% | 565.15 |
SA | 3318 | 1324 | 39.90% | 684.60 |
ACO-LS | 3318 | 1588 | 46.85% | 535.40 |
ACO-SA | 3318 | 1658 | 50.02% | 606.66 |
表4不同卫星观测调度算法在第四组试验场景的计算结果
算法 | CT | NT | Cov(%) | Time(s) |
TS | 1862 | 1625 | 87.27% | 541.12 |
GA | 1862 | 1678 | 90.12% | 655.45 |
SA | 1862 | 1598 | 85.82% | 743.15 |
ACO-LS | 1862 | 1699 | 91.25% | 524.36 |
ACO-SA | 1862 | 1742 | 93.56% | 560.65 |
表5不同卫星观测调度算法在第五组试验场景的计算结果
算法 | CT | NT | Cov(%) | Time(s) |
TS | 1456 | 1355 | 93.06% | 482.32 |
GA | 1456 | 1382 | 94.92% | 505.24 |
SA | 1456 | 1275 | 87.57% | 661.32 |
ACO-LS | 1456 | 1398 | 96.02% | 458.55 |
ACO-SA | 1456 | 1418 | 97.39% | 500.48 |
表6不同卫星观测调度算法在第六组试验场景的计算结果
上表中技术术语ACO-SA为本发明算法,TS为禁忌搜索,GA为遗传算法,SA为模拟退火,ACO-LS为利用蚁群优化算法与局部搜索策略结合,生成多星规划调度问题的解的算法。
本实施例提出的分而治之框架具有通用性,适用于很多类型的规划调度,例如无人机调度和车辆路径规划等。后续的工作是进一步研究敏捷卫星和超敏卫星的调度问题,以及成像和数传的一体化调度问题。
以上包含了本发明优选实施例的说明,这是为了详细说明本发明的技术特征,并不是想要将发明内容限制在实施例所描述的具体形式中,依据本发明内容主旨进行的其他修改和变型也受本专利保护。本发明内容的主旨是由权利要求书所界定,而非由实施例的具体描述所界定。
Claims (5)
1.一种基于分治策略的多星观测分层调度方法,其特征在于,包括以下步骤:
S1、利用分配算法把任务分配至各轨道圈次上,形成各轨道圈次的任务集;
S2、利用近似优化算法求解各轨道圈次上任务集的调度序列;
S3、根据各轨道圈次调度序列反馈的未生成调度序列的任务重新更新分配算法的分配方案,进而形成新的各轨道圈次的新的任务集;
S4、重复步骤S1、S2、S3直至达到多星观测分层调度的终止条件;
步骤S2中,所述近似优化算法为模拟退火算法;
所述模拟退火算法获取当前轨道圈次的初始调度序列的步骤包括:
S211、将单轨道圈次任务集中的待规划调度任务根据权重从大到小依次排列,形成单轨道圈次未被调度的初始任务集U;
S212、从U中选择权重最大的任务tk,应用***和删除邻域结构变换,把tk加入到局部解S中生成新的局部解S′,若S′的适应度高于S,则对S进行更新,令S=S′,否则不对S进行更新,同时将tk从U中移至Uj中,其中Uj表示在轨道圈次上未被调度的任务集;
S213、重复步骤S212直至U中的任务全部遍历完成,获取当前轨道圈次的调度序列S与当前轨道圈次未被调度的任务集Uj;
所述模拟退火算法的领域搜索结构包括基于贪婪准则的邻域搜索结构与基于概率搜索的邻域搜索结构,所述模拟退火算法包括邻域搜索结构动态选择策略:根据两种邻域搜索结构在之前分配方案内的优化性能,决定在后续分配方案中领域搜索结构的选择,具体为:
假定pro1和pro2分别表示为选择基于贪婪准则的邻域搜索结构和基于概率搜索的邻域搜索结构的概率,在算法初始化时,设定proi=0.5,i=1,2,每隔一定的迭代次数Itr,通过以下规则更新每个邻域结构的选择概率:
式中,表示中间变量;η为惯性权重因子,表示邻域结构之前的选择概率所占的比重;(1-η)表示当前最新历史搜索经验对更新选择概率的权重;Seli为最近Itr次迭代过程中,第i个邻域结构被选中的次数;Suci表示利用第i个邻域结构产生了更高质量的解的次数;
基于概率搜索的邻域搜索结构的模拟退火算法具体包括:
S231、获取当前轨道圈次的调度序列S与当前轨道圈次未被调度的任务集Uj;
S233、将tk***到当前轨道圈次的调度序列S中去,判断是否满足侧摆时间约束,若满足则直接形成新的调度序列S″,若不满足则删除与tk冲突的任务后形成新的调度序列S″,其中删除的任务存入Su中,Su是所有轨道圈次未被调度的任务集合;
S234、判断新的调度序列是否满足能量约束和存储容量约束:
若满足,则重复步骤S231、S232、S233、S234直至Uj中所有的任务都遍历完成,获取当前轨道上的最终调度序列S″并将当前轨道未被调度的任务存入集合Su中;
若不满足,则删除任务tj的概率为其中indexj为任务tj优先级指标,直到满足能量约束和存储容量约束后重复步骤S231、S232、S233、S234直至Uj中所有的任务都遍历完成,获取当前轨道上的最终调度序列S″并将当前轨道未被调度的任务存入集合Su中,其中Su是所有轨道圈次未被调度的任务集合;
在基于概率搜索的邻域搜索结构中,不但考虑了任务的权重,还考虑完成每个任务需潜在消耗的侧摆资源和时间窗口资源,综合权重、侧摆角度和时间窗口,在***和替换邻域结构变换过程中,为每个任务tk计算一个优先级指标indexk,首先分别计算权重指标,时间窗口指标和侧摆角度指标:
权重指标为:
式中,iWk为权重指标;wk是第k个任务的权重;wi是第i个任务的权重;nj是Uj任务数量;
时间窗口指标为:
式中,iTk为时间窗口指标;spank为tk的时间窗口长度;spani为ti的时间窗口长度;
侧摆角度指标为:
式中,iTk为侧摆角度指标;θi在轨道oj上观测任务ti时传感器所需测摆的角度;θk在轨道oj上观测任务tk时传感器所需测摆的角度;
优先级指标indexk最终表示为:
indexk=iWk α·iTk β·iθk γ
式中,α,β,γ分别表示不同要素的影响因子,由用户针对具体问题具体设置;
基于贪婪准则的邻域搜索结构的模拟退火算法具体包括:
S221、获取当前轨道圈次的调度序列S与当前轨道圈次未被调度的任务集Uj;
S222、在Uj中取出一个尚未调度的权重最大的任务tk,令Uj=Uj-tk;
S223、将tk***到当前轨道圈次的调度序列S中去,判断是否满足侧摆时间约束,若满足则直接形成新的调度序列S″,若不满足则删除与tk冲突的任务后形成新的调度序列S″,其中删除的任务存入Su中,Su是所有轨道圈次未被调度的任务集合;
S224、判断新的调度序列S″是否满足能量约束和存储容量约束:
若满足,则重复步骤S221、S222、S223、S224直至Uj中所有的任务都遍历完成,获取当前轨道上的最终调度序列S″并将当前轨道未被调度的任务存入集合Su中;
若不满足,则依次删除权重最小的任务直到满足能量约束和存储容量约束后重复步骤S221、S222、S223、S224直至Uj中所有的任务都遍历完成,获取当前轨道上的最终调度序列S″并将当前轨道未被调度的任务存入集合Su中,其中Su是所有轨道圈次未被调度的任务集合。
2.根据权利要求1所述基于分治策略的多星观测分层调度方法,其特征在于,所述分配算法为蚁群优化算法。
3.根据权利要求2所述基于分治策略的多星观测分层调度方法,其特征在于,步骤S3中具体包括:
S31、根据各轨道圈次调度序列获取所有轨道圈次未被调度的任务集合Su;
S32、从Su中选择一个权重最高的任务,将其***到冲突权重最小的轨道圈次上,若该任务的冲突权重小于该任务自身的权重,则***成功并更新分配方案,否则不更新分配方案;
S33、重复步骤S32直至Su中所有的未被调度的任务都遍历完成。
4.根据权利要求1所述基于分治策略的多星观测分层调度方法,其特征在于,步骤S4中,所述多星观测分层调度的终止条件为计算时间达到时间预设值或迭代次数达到次数预设值。
5.一种基于分治策略的多星观测分层调度***,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至4中任一项所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810967096.7A CN109087023B (zh) | 2018-08-23 | 2018-08-23 | 一种基于分治策略的多星观测分层调度方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810967096.7A CN109087023B (zh) | 2018-08-23 | 2018-08-23 | 一种基于分治策略的多星观测分层调度方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109087023A CN109087023A (zh) | 2018-12-25 |
CN109087023B true CN109087023B (zh) | 2021-03-02 |
Family
ID=64794304
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810967096.7A Active CN109087023B (zh) | 2018-08-23 | 2018-08-23 | 一种基于分治策略的多星观测分层调度方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109087023B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109767128B (zh) * | 2019-01-15 | 2021-06-11 | 中国人民解放军国防科技大学 | 基于机器学习的成像卫星自主任务规划方法 |
CN110766284B (zh) * | 2019-09-24 | 2022-04-08 | 合肥工业大学 | 多星任务的合成方法和*** |
CN110728447B (zh) * | 2019-10-10 | 2021-03-09 | 合肥工业大学 | 一种最早完成区域目标覆盖的分区卫星任务规划方法 |
CN111176807B (zh) * | 2019-12-09 | 2023-12-12 | 中国空间技术研究院 | 一种多星协同任务规划方法 |
CN112016812B (zh) * | 2020-08-06 | 2022-07-12 | 中南大学 | 多无人机任务调度方法、***及存储介质 |
CN112070370B (zh) * | 2020-08-21 | 2024-03-08 | 中南大学 | 中继卫星任务规划方法、***及存储介质 |
CN112613644B (zh) * | 2020-12-09 | 2022-11-04 | 合肥工业大学 | 多星多站集成调度方法和*** |
CN112737660B (zh) * | 2020-12-09 | 2022-06-10 | 合肥工业大学 | 多星多站数据下传调度方法和*** |
CN113158384B (zh) * | 2021-03-03 | 2021-10-08 | 东北石油大学 | 一种油气管道路由规划方法及*** |
CN113408063B (zh) * | 2021-07-20 | 2022-06-21 | 中南大学 | 面向多星多侦察目标的轨道机动优化方法 |
CN113487220B (zh) * | 2021-07-23 | 2022-07-15 | 中南大学 | 面向静态目标观测的空天异构对地观测资源协同调度方法 |
CN113919122A (zh) * | 2021-08-05 | 2022-01-11 | 合肥工业大学 | 基于模拟退火算法的多星任务调度方法和*** |
CN114604443B (zh) * | 2022-02-17 | 2023-06-06 | 中国电子科技集团公司第十研究所 | 一种真实复杂约束条件下的多星测控调度方法及*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104361234A (zh) * | 2014-11-15 | 2015-02-18 | 北京理工大学 | 一种复杂约束条件下多星多任务观测调度优化方法 |
CN106951312A (zh) * | 2016-08-31 | 2017-07-14 | 合肥工业大学 | 多星任务的调度方法及*** |
CN108335012A (zh) * | 2017-12-26 | 2018-07-27 | 佛山科学技术学院 | 一种智能遥感卫星层次化分布式自主协同任务规划*** |
-
2018
- 2018-08-23 CN CN201810967096.7A patent/CN109087023B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104361234A (zh) * | 2014-11-15 | 2015-02-18 | 北京理工大学 | 一种复杂约束条件下多星多任务观测调度优化方法 |
CN106951312A (zh) * | 2016-08-31 | 2017-07-14 | 合肥工业大学 | 多星任务的调度方法及*** |
CN108335012A (zh) * | 2017-12-26 | 2018-07-27 | 佛山科学技术学院 | 一种智能遥感卫星层次化分布式自主协同任务规划*** |
Non-Patent Citations (1)
Title |
---|
"An adaptive Simulated Annealing-based satellite observation scheduling method combined with a dynamic task clustering strategy";WuGuoHua等;《Computers&Industrial Engineering》;20171130;摘要,第3、5、6节 * |
Also Published As
Publication number | Publication date |
---|---|
CN109087023A (zh) | 2018-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109087023B (zh) | 一种基于分治策略的多星观测分层调度方法及*** | |
CN109829636B (zh) | 基于动态优先级的卫星应急任务调度规划方法 | |
CN109948944B (zh) | 一种卫星任务调度方法及*** | |
US10726368B2 (en) | Innovative satellite scheduling method based on genetic algorithms and simulated annealing and related mission planner | |
CN113128828B (zh) | 一种基于多智能体强化学习的卫星观测分布式在线规划方法 | |
He et al. | Scheduling multiple agile earth observation satellites with an edge computing framework and a constructive heuristic algorithm | |
CN109039428B (zh) | 基于冲突消解的中继卫星单址天线调度随机搜索方法 | |
CN105225003A (zh) | 一种布谷鸟搜索算法解决uav多任务侦察决策问题的方法 | |
Wu et al. | Ensemble of metaheuristic and exact algorithm based on the divide-and-conquer framework for multisatellite observation scheduling | |
CN109933842A (zh) | 一种基于约束满足遗传算法的移动目标单星任务规划方法 | |
Song et al. | Improved genetic algorithm with local search for satellite range scheduling system and its application in environmental monitoring | |
CN103345504A (zh) | 单星调度的算子构建方法 | |
CN111612384B (zh) | 一种具有时间分辨率约束的多星接力任务规划方法 | |
CN109034670B (zh) | 一种卫星在轨活动规划方法及*** | |
Eddy et al. | Markov decision processes for multi-objective satellite task planning | |
CN112990478B (zh) | 联邦学习数据处理*** | |
CN113051815A (zh) | 一种基于独立指针网络的敏捷成像卫星任务规划方法 | |
CN117829277B (zh) | 面向多类型卫星的应急观测任务调度方法、***及装置 | |
CN113919122A (zh) | 基于模拟退火算法的多星任务调度方法和*** | |
CN109299490A (zh) | 基于带禁忌列表蚁群算法视频卫星对地凝视任务规划方法 | |
CN116579513A (zh) | 基于改进鲸鱼优化算法的物资配送路径优化方法和*** | |
Long et al. | Improved particle swarm optimization with reverse learning and neighbor adjustment for space surveillance network task scheduling | |
CN116629511A (zh) | 一种不确定环境下基于两阶段混合调度的多星动态任务规划方法及装置 | |
CN113657718B (zh) | 一种多机器人动态联盟任务分配方法及相关装置 | |
Veldkamp et al. | Designing item pools for adaptive testing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |