CN108023637B - 一种同构多星在线协同方法 - Google Patents
一种同构多星在线协同方法 Download PDFInfo
- Publication number
- CN108023637B CN108023637B CN201711279412.3A CN201711279412A CN108023637B CN 108023637 B CN108023637 B CN 108023637B CN 201711279412 A CN201711279412 A CN 201711279412A CN 108023637 B CN108023637 B CN 108023637B
- Authority
- CN
- China
- Prior art keywords
- observation
- satellite
- task
- communication
- earth
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000004891 communication Methods 0.000 claims abstract description 246
- 230000006854 communication Effects 0.000 claims abstract description 239
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 124
- 230000007246 mechanism Effects 0.000 claims abstract description 14
- 230000001360 synchronised effect Effects 0.000 claims description 46
- 239000013598 vector Substances 0.000 claims description 27
- 239000003999 initiator Substances 0.000 claims description 22
- 238000010276 construction Methods 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 16
- 238000009826 distribution Methods 0.000 claims description 10
- 238000003384 imaging method Methods 0.000 claims description 10
- 238000012163 sequencing technique Methods 0.000 claims description 6
- 230000008520 organization Effects 0.000 claims description 3
- 230000008439 repair process Effects 0.000 claims description 3
- 230000007175 bidirectional communication Effects 0.000 claims description 2
- 230000001186 cumulative effect Effects 0.000 claims description 2
- 230000000750 progressive effect Effects 0.000 claims description 2
- 230000003993 interaction Effects 0.000 description 16
- 238000005096 rolling process Methods 0.000 description 12
- 239000008186 active pharmaceutical agent Substances 0.000 description 8
- 230000004044 response Effects 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000001427 coherent effect Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 239000003795 chemical substances by application Substances 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008447 perception Effects 0.000 description 2
- 206010037844 rash Diseases 0.000 description 2
- 230000001960 triggered effect Effects 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000012782 phase change material Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004043 responsiveness Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000011273 social behavior Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002459 sustained effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/1851—Systems using a satellite or space-based relay
- H04B7/18513—Transmission in a satellite or space-based system
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/18521—Systems of inter linked satellites, i.e. inter satellite service
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/19—Earth-synchronous stations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Aviation & Aerospace Engineering (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Radio Relay Systems (AREA)
Abstract
本发明公开了一种同构多星在线协同方法。在所述同构多星在线协同方法中,采用集中‑分布式协同架构,协同决策功能仅配置在地球静止轨道通信卫星上,而对地观测卫星进行分布式计算,所述同构多星在线协同方法采用基于拍卖机制的任务协同分配算法,拍卖由作为拍卖商的地球静止轨道通信卫星组织进行,每个对地观测卫星根据自身的局部信息来计算与应急观测任务相对应的投标出价值,并以此进行投标,每个对地观测卫星将其出价通知拍卖商,然后地球静止轨道通信卫星决定每个应急观测任务的胜出者,如果有多个通信中继节点相互关联,则拍卖商的角色能够在这些通信中继节点之间进行转移。
Description
技术领域
本发明涉及卫星技术领域,特别是涉及一种同构多星在线协同方法,更具体地,涉及一种采用集中-分布式协同架构的同构多星在线协同方法。
背景技术
随着星上载荷探测能力和图像处理能力的提升,对地观测卫星可在星上发现有价值的观测目标并生成进一步的观测任务请求(通常是具体有较高收益或优先等级的应急观测任务),也可接收其他卫星通信传递过来的观测任务请求(应急观测任务请求),而这些观测任务请求是动态随机到达的且时效性要求高,需要在线解决,在多星(尤其是同构多星)的情况下,从而需要有多星在线协同调度机制来进行协同调度。
发明内容
本发明的目的在于提供一种同构多星在线协同方法来实现同构多星的在线协同,尤其是更好地应对应急观测任务。
为实现上述目的,本发明提供一种同构多星在线协同方法,在所述方法中,对地观测卫星作为资源观测卫星,地球静止轨道通信卫星作为通信中继节点,对地观测卫星与地球静止轨道通信卫星通信,对地观测卫星之间不能直接进行卫星通信,所述同构多星在线协同方法采用集中-分布式协同架构,协同决策功能仅配置在地球静止轨道通信卫星上,而对地观测卫星进行分布式计算,所述同构多星在线协同方法采用基于拍卖机制的任务协同分配算法,拍卖由作为拍卖商的地球静止轨道通信卫星组织进行,每个对地观测卫星根据自身的局部信息来计算与应急观测任务相对应的投标出价值,并以此进行投标,每个对地观测卫星将其出价通知拍卖商,然后地球静止轨道通信卫星决定每个应急观测任务的胜出者,如果有多个通信中继节点相互关联,则拍卖商的角色能够在这些通信中继节点之间进行转移。
优选地,所述地球静止轨道通信卫星采用下述的中标决策方式:
1).如果有多个对地观测卫星能够在截止时间之前完成应急观测任务,则选择使***总收益增量最大的对地观测卫星作为中标者;如果有多个卫星的投标出价值对应于最高的总收益增量,则以实际开始观测时间最早的对地观测卫星具有中标优先权;
2).如果只有一个对地观测卫星可以在截止日期之前完成应急观测任务,同时***总收益的增量为正值,则选择该卫星作为中标者;
3).如果当对地观测卫星执行该应急观测任务产生的***总收益增量均是负值,或者没有对地观测卫星可以在截止时间之前完成应急观测任务,则地球静止轨道通信卫星舍弃该应急观测任务。
优选地,以下式计算应急观测任务的时间可用性wioj,计算如下:
其中,olsioj为第o批次中第j个应急观测任务的最晚观测开始时间,且olsioj=max{olfioj-poj-(2·maxθ)/sproll,oesioj}。
优选地,每颗对地观测卫星具有六个信息向量,具体如下:
1.)任务束,1≤n≤|BUio|,该任务束表示对地观测卫星i从第o批应急观测任务中选取并调度成功的任务集合,同时按照各任务添加到束中的时间先后进行排序。当前任务束的长度为|BUio|,且要小于一批次内的应急观测任务数量u,|BUio|≤u;当任务束为空时,用BUio=φ和|BUio|=0表示,
4.)胜出者列表WAio@{waio1,...,waiou},长度为u,其中waioj表示对地观测卫星i当前所认为的针对第o批应急观测任务中第j个任务的胜出者,具体取值与胜出者编号相对应,当waioj=φ时,对地观测卫星i认为当前对于该任务无胜出者,
5.)中标出价列表WBio@{wbio1,...,wbiou},长度为u,其中wbioj表示对应的胜出者所给出的出价,当取值为0时,表示该任务当前无胜出者,
优选地,所述任务协同分配算法包括两个阶段的迭代:束构建阶段和一致性构建阶段,其中前者对应于每个对地观测卫星以贪婪搜索的形式生成排序后的任务束,后者对应于识别任务分配冲突,并通过与相邻对地观测卫星之间的局部通信进行冲突消解,上述的两个阶段不断地迭代重复,直至达到收敛。
优选地,对地观测卫星与地球静止轨道通信卫星之间的通信是间歇性的,所述任务协同分配算法进一步包括一个新的阶段:同步通信环路预测阶段,其中,同步通信环路是用于两个或多个对地观测卫星(卫星)的双向通信链路,以向其他对地观测卫星共享信息向量并从环路中接收来自其他对地观测卫星的反馈信息,一个同步通信环路构建于发起者和一个或多个响应者的通信时间窗口之间均存在重叠时间的基础上,
其描述形式是一个多元组<ID,I,R,ES,LF>SCL:
—ID是一个标识符;
—I是通信发起者编号,即第一个对其信息进行分享;
—R是接收并分享信息的通信响应者编号;
—ES是通信环路的最早开始时间;
—LF是通信环路的最晚结束时间,
其中通信发起者I和响应者R都是对地观测卫星,并且ES和LF是由发起者和所有响应者的通信时间窗口上的重叠时间区间确定的,同步通信环路预测阶段用于预测同步通信环路的存在,且提供所述同步通信环路的多元组<ID,I,R,ES,LF>SCL,同步通信环路预测阶段是束构建阶段的前提。
附图说明
图1是同步通信环路示意图。
图2是异步通信环路的两种情况示意图。
图3是整个调度区间上在线协同调度求解的总体概念图。
图4是中继节点与对地观测卫星之间的通信示意图。
图5是合同网协议中的信息流示意图。
图6示出空间和地面段之间可用的通信。
图7示出针对GEO之间不互联与互联情况的m-CBBA算法结果。
图8示出各算法在任务空间分布和不同滚动调度周期下的结果。
具体实施方式
下面结合附图对本发明的实施例进行详细说明。
本发明中的卫星(对地观测卫星LEO)为自主智能体(自主卫星),具有下述特征:自主感知,自主规划调度和自主执行。
自主的体现方面包括精密姿态确定、轨道确定、姿态敏感器/驱动器与载荷校准、姿态控制、轨道机动、状态监视与预测分析、故障检测、诊断、隔离和修复,预报式建模(在计划调度、状态监视分析和轨道确定保持等中,需要预报式模型的支持,如星历、太阳强度等信息),任务计划与调度,载荷控制与配置,数据存储与通信和图像数据处理等方面。例如,能够通过感知来自主生成应急调度任务。
分布式对地观测卫星***具有独特的轨道特性和间隙性通信特性,通信约束较多,***内的通信连通性并不能实时保证。因此,在动态任务环境中,面对有限的星上计算能力和时变的通信约束,如何设计出高效适用的多星在线协同调度机制与算法,在***效能和通信成本代价之间取得平衡,是个非常具有挑战性的问题。
分布式对地观测卫星***是位于不同轨道上的卫星集合,其任务是对地球表面的各个热点区域进行观测成像,比如对应于火山活动或森林火灾的热点。各个对地观测卫星配备有单一的成像载荷,其相机中心视场线可以左右摆动。侧摆活动需要有持续时长限制,因此观测活动太靠近的任务请求不能由同一颗卫星实现。此外,卫星配备有异常探测设备,能够发现热点事件,并在星上生成观测任务请求。
观测任务请求的在线到达主要来自三个主要来源:星座内卫星自身生成的任务,星座内其他卫星上生成并传递过来的任务以及地面站的应急上注。
对地观测卫星处于低轨道上,星间链路并不总是可用的。例如,RapidEye是德国勃兰登堡RapidEye AG公司的商业多光谱对地观测卫星项目,包括由五颗小卫星组成的卫星星座。五颗RapidEye对地观测卫星部署在高度为630公里的同一太阳同步轨道上,并且等距放置以确保一致的成像条件和短的重访时间间隔。这些卫星在其轨道平面上彼此时间间隔约为19分钟。可见五颗卫星之间并不存在星间直接通信链路,因为卫星的数量和高度不允许星间链路的存在,因此我们只能在地面站中离线操作和协同这些卫星。本发明假设存在数量有限的通信中继节点,可提供比以前更多的通信机会。而这些中继节点可以位于地面上,如地面通信站或移动通信车辆,也可以是部署在太空的通信卫星等等。
在卫星星座的在线调度问题中,有两种类型的任务,常规观测任务和应急观测任务。地面控制中心接收用户的观测请求并对其进行安排,然后通过星地链路上传已生成的观测计划。在本发明中,将这些作为预先收集的观测请求称为常规观测任务。相比之下,从星座内外的卫星上随机到达的观测请求被称为应急观测任务,而不可预见的应急观测任务需要在无地面干预的情况下进行在线调度。应急观测任务对应的收益值通常高于常规观测任务,因为这类任务具有高价值并且亟需快速响应。对于更一般的情况,应急观测任务往往是以批次的形式随机到达。
从更广泛的多对地观测卫星***协同角度,对地观测卫星一次只能对地面上的一个目标进行成像。其次,各任务针对不同的观测目标,并且相互之间是独立的,每个任务可以由一颗卫星完成。第三,卫星调度问题是任务请求数量明显大于可用的对地观测卫星数量的过度订购问题,因此为每个卫星构建任务计划。
为了表述的统一简便,对本发明后续中所用到的字符进行统一定义:
下标
i,i′,i″:低轨LEO(Low Earth Orbit)卫星编号,i=1,2,K,n1
g:通信中继节点编号,h=1,2,K,n2
j,j′:应急观测任务编号,j=1,2,K,u
k:常规观测任务编号,k=1,2,...,vi
q:通信时间窗口编号,q=1,2,K,mig
o,o′,o″:任务批次编号,o=1,2,K,l
参数量
H:整个调度区间
tl:整个调度区间H的总时长
n1:***中对地观测卫星的总数量
n2:可用的通信中继节点数量
u:一批次任务中的应急观测任务数量
l:整个调度区间H内应急观测任务总批次数
vi:对地观测卫星i上已上传的常规任务数量
ro:中继节点向对地观测卫星发布第o批次应急观测任务的时间
trio:第o批次应急观测任务的时间到达对地观测卫星i的时间
sproll:对地观测卫星的侧摆角速度
maxθ:对地观测卫星的最大侧摆角
oesioj:针对对地观测卫星i,第o批次中第j个应急观测任务的最早观测开始时间
olfioj:针对对地观测卫星i,第o批次中第j个应急观测任务的最晚观测结束时间
olsioj:针对对地观测卫星i,第o批次中第j个应急观测任务的最晚观测开始时间
otwioj:针对对地观测卫星i,第o批次中第j个应急观测任务的观测时间窗口
poj:第o批次中第j个应急观测任务的成像时长
eoj:第o批次中第j个应急观测任务的收益,由管理者或星上决策给出
θioj:针对对地观测卫星i,第o批次中第j个应急观测任务的侧摆角度
bioj:针对对地观测卫星i,第o批次中第j个应急观测任务的实际开始时间
cioj:针对对地观测卫星i,第o批次中第j个应急观测任务的实际结束时间
oesik:对地观测卫星i上第k个常规任务的最早开始时间
olfik:对地观测卫星i上第k个常规任务的最晚结束时间
pik:对地观测卫星i上第k个常规任务的成像时长
eik:对地观测卫星i上第k个常规任务的收益,由管理者给出
bik:对地观测卫星i上第k个常规任务的实际开始时间
cik:对地观测卫星i上第k个常规任务的实际结束时间
siojk:对地观测卫星i上,当第k个常规任务紧接着第o批次中第j个应急观测任务执行时的次序依赖转换时间
sikoj:对地观测卫星i上,当第o批次中第j个应急观测任务紧接着第k个常规任务执行时的次序依赖转换时间
siojo′j′:对地观测卫星i上,当第o′批次中第j′个应急观测任务紧接着第o批次中第j个应急观测任务执行时的次序依赖转换时间,其中(o≠o′)||(j≠j′)=1
ctwigq:对地观测卫星i与中继节点g之间的第q个通信时间窗口
mig:对地观测卫星i与中继节点g之间的通信时间窗口数量
utioj:针对对地观测卫星i,第o批次中第j个应急观测任务的表示形式
wioj:当第o批次中第j个应急观测任务可在对地观测卫星i上调度执行时,等于1,否则,等于0。
通信中继节点与***内所有对地观测卫星之间的通信时间窗口集合,定义如下:
具体而言,通信中继节点g与对地观测卫星i之间的通信时间窗口集合,定义如下:
其中通信时间窗口ctwigq定义为:
ctwigq=[cesigq,clfigq]
其中cesigq是最早通信开始时间,clfigq是最晚通信结束时间。
当且仅当对地观测卫星i与中继节点之间存在一个或多个通信时间窗口时,第o批次应急观测任务才可以从中继节点发布给对地观测卫星i。对批次任务的时间可用性wio,计算如下:
当且仅当第o批次应急观测任务在最晚观测开始时间olsioj之前到达对地观测卫星i,第o批次中第j个应急观测任务可在对地观测卫星i上进行调度执行。对应急观测任务的时间可用性wioj,计算如下:
其中olsioj=max{olfioj-poj-(2·maxθ)/sproll,oesioj}。
根据实际应用情况面向整个***的假设如下:
1.不同来源的应急观测任务在通信中继节点进行汇总,并由其发布;
2.星上载荷是光学成像载荷,需考虑地影的影响;
3.对地观测卫星与中继节点之间的通信延迟时间较短,可忽略不计;
4.在整个调度区间内,星上固存空间和持续工作时间是足够的。
当第o″批次应急观测任务到达卫星星座时,在线协同调度问题可用混合整数线性规划模型MILP(mixed integer linear programming)进行表示:
决策变量
cioj+(siojk+pik)yiojk+olfioj(yiojk-1)≤cik (4.2)
cik+(sikoj+poj)yikoj+olfik(yikoj-1)≤cioj (4.3)
cio'j'+(sio'j'oj+pioj)yio'j'oj+olfio'j'(yio'j'oj-1)≤cioj (4.4)
(trio+poj)xioj+sikojyikoj+sio'j'ojyio'j'oj≤cioj (4.5)
pikzik+siojkyiojk≤cik (4.6)
(oesioj+poj)xioj≤cioj (4.7)
(oesik+pik)zik≤cik (4.8)
cioj≤olfiojxioj (4.9)
cik≤olfikzik (4.10)
bioj+poj=cioj (4.11)
bik+pik=cik (4.12)
xioj={0,1},yiojk={0,1},yikoj={0,1},yio'j'oj={0,1},zik={0,1},wioj={0,1} (4.15)
在集中-分布式协同架构中,协同决策功能仅配置在中心节点上,而其他节点可进行分布式计算,而与之相对应的是基于市场机制的任务协同分配方法,包括拍卖机制。拍卖是由拍卖商组织进行,其中每个对地观测卫星根据自身的局部信息计算与任务相对应的投标出价值,并以此进行投标。每个对地观测卫星将其出价通知拍卖商,然后拍卖商再决定每个任务的胜出者。该算法机制可保证生成无冲突的解决方案,因为拍卖商只选择一个对地观测卫星作为胜出者。由于该算法机制面向集中-分布式架构的,则投标计算是分布式进行的。
在多星协同应用场景中,通信中继节点扮演拍卖商的角色,负责对应急观测任务进行招标,评估各投标卫星的出价并对其进行中标结果反馈。如果有多个通信中继节点相互关联,则拍卖商的角色可在这些通信中继节点之间进行转移,以提高***的可靠性。作为具体的观测资源对地观测卫星,对地观测卫星之间不能直接进行卫星通信,同时拍卖商和资源对地观测卫星之间的通信链路是间歇性的。
针对该算法机制中的中标决策阶段,提出以下决策规则:
1.如果有多个对地观测卫星对地观测卫星可以在截止日期之前完成应急观测任务,则选择使***总收益增量最大的对地观测卫星作为中标者。在这种情况下,如果有多个卫星的投标出价值对应于最高的总收益增量,则实际开始最早的对地观测卫星具有中标优先权。
2.如果只有一个对地观测卫星可以在截止日期之前完成应急观测任务,同时***总收益的增量为正值,则选择该卫星作为中标者。
3.如果当对地观测卫星执行该应急观测任务产生的***总收益增量均是负值,或者没有对地观测卫星可以在截止日期之前完成应急观测任务,则拍卖商舍弃该应急观测任务。
在图5中,中继节点1作为拍卖商,中继节点2充当通信中继的角色,而LEO_3中标。中继节点1通过可用的通信链路进行拍卖发布,投标收集和中标反馈。
集中-分布式架构下基于合同网协议的在线协同算法,其伪代码如下:
输入:
应急观测任务集合Ut={ut1,ut2,...,ut(u·l)}
通信时间窗口集合CTW
每批次中应急观测任务的总数量u
输出:
针对具体合同网协议的实现,本发明提出两种方式:单项任务下的招投标和批次任务下的招投标,其中前者是指每次迭代只对单个任务进行招投标并决定该任务的胜出者,后者则是指每次迭代针对批量任务进行招投标并决定多个任务的胜出者。
单项任务下的合同网协议算法SingleItem-CNP-based Online CoordinationScheduling Algorithm
每批次应急观测任务的冒泡排序计算复杂度是O(u2),采用二分法对给定的批次任务和给定对地观测卫星进行时间可用性检查,对应的计算复杂度是 针对给定的对地观测卫星,每个应急观测任务的在线调度计算复杂度为则单项任务下的合同网协议算法全局计算复杂度为
批次任务下的合同网协议算法Algorithm 2.Batch-CNP-based OnlineCoordination Scheduling Algorithm
%%批次任务下的投标阶段
%%批次任务下的选标与中标阶段
面向分散式架构,采用分散式算法,特别是拍卖算法和一致性协议,对任务进行无冲突分配,适用于很难达到全局一致性态势感知的通信环境不可靠情况。
基于同步通信的改进一致性束算法
针对通信约束下的多星在线协同任务调度这一复杂组合优化问题,并面向分散式架构,本发明提出基于同步通信的改进一致性束算法m-CBBA算法。
采用同步通信的算法对何时进行通信有严格的规则定义,同时,同步策略会对通信情况进行排队,只能在特定的事件触发下才能执行,这保证了在程序执行过程中算法具有预计的状态。在大多数迭代算法中,均采用同步算法。在同步迭代算法中,各模块会执行并行计算,分享状态变量,然后等待下一个触发事件,以进行算法的下一次迭代过程。在这个过程中,对每个对地观测卫星的状态进行保证。
采用同步算法的成本代价在于需要维持同步行为,特别是在分布式和分散式算法中,要实现同步协同行为,则需要满足协同各方的同时通信可见,才能进行协同同步化。
一致性束算法CBBA是一分布式拍卖算法,以实时分布式的方式处理复杂的约束,为多对地观测卫星多任务分配问题给出任务分配方案,则在实际应用中非常高效。
一致性束算法具有多个特征,可用于实现异构***协同任务调度,假定Nt为待分配的总任务数,Lt为每束的最大任务数量,Na为可用的对地观测卫星数量,D为该对地观测卫星网络中的对地观测卫星数量,具体如下:
(1)一致性束算法是面向分散式架构的;
(2)一致性束算法属于多项式时间算法,则束构建阶段的计算复杂度是O(NtLt);
(3)一致性束算法会在max{Nt,LtNa}D次迭代内收敛,从而一致束算法能够很好地适应网络规模和任务数量的变化;
(4)不同的求解目标,对地观测卫星内部模型以及相关约束都可以通过设计合适的评分函数进行体现。如果评分函数满足边际收益衰减特性,可保证得到好的可行解。
针对分布式多星***,要执行一致性束算法,则需要每颗对地观测卫星具有六个信息向量,具体如下:
任务束,1≤n≤|BUio|。该任务束表示对地观测卫星i从第o批应急观测任务中选取并调度成功的任务集合,同时按照各任务添加到束中的时间先后进行排序。当前任务束的长度为|BUio|,且要小于一批次内的应急观测任务数量u,|BUio|≤u;当任务束为空时,用BUio=φ和|BUio|=0表示。
胜出者列表WAio@{waio1,...,waiou},长度为u,其中waioj表示对地观测卫星i当前所认为的针对第o批应急观测任务中第j个任务的胜出者,具体取值与胜出者编号相对应。当waioj=φ时,对地观测卫星i认为当前对于该任务无胜出者。
中标出价列表WBio@{wbio1,...,wbiou},长度为u,其中wbioj表示对应的胜出者所给出的出价,当取值为0时,表示该任务当前无胜出者。
一致性束算法CBBA包括两个阶段的迭代:束构建阶段和一致性构建阶段,其中前者对应于每个对地观测卫星以贪婪搜索的形式生成排序后的任务束,后者对应于识别任务分配冲突,并通过与相邻对地观测卫星s之间的局部通信进行冲突消解。这两个阶段不断地迭代重复,直至达到收敛。
在一致性束算法CBBA的基础上,改进一致性束算法m-CBBA在原有的两个阶段之前增加一个新的阶段,同步通信环路预测阶段。
由于实际应用中存在严格的通信约束,则对地观测卫星s之间的强连通性是不可能达到的,因此该通信环路预测可在提升***全局性能和保证应急观测任务无冲突调度之间取得均衡,而不是实现在胜出者列表上的全局收敛。
阶段1:同步通信环路预测阶段
同步通信环路预测阶段用于预测同步通信环路的存在。同步通信环路是用于两个或多个对地观测卫星(卫星)的双向通信链路,以向其他对地观测卫星共享信息向量并从环路中接收来自其他对地观测卫星的反馈信息。一个同步通信环路构建于发起者和一个或多个响应者的通信时间窗口之间均存在重叠时间的基础上。同步通信环路的描述形式是一个多元组<ID,I,R,ES,LF>SCL:
—ID是一个标识符;
—I是通信发起者编号,即第一个对其信息进行分享;
—R是接收并分享信息的通信响应者编号;
—ES是通信环路的最早开始时间;
—LF是通信环路的最晚结束时间。
其中通信发起者I和响应者R都是对地观测卫星,并且ES和LF是由发起者和所有响应者的通信时间窗口上的重叠时间区间确定的。
每次同步共享反馈交互的都是在一可用的同步通信环路上触发形成的,其中分享上行时间戳与分享下行时间戳相同,而反馈上行时间戳与反馈下行时间戳相同。所对应的描述形式<id,i,i′,tUS-DS,tUF-DF>源于同步通信环路<ID,I,R,ES,LF>SCL,其中i∈I,i′∈R,tUS-DS,tUF-DF∈[ES,LF],tUS-DS<tUF-DF。
对于每个对地观测卫星,当一批应急观测任务到达时,都需要预测自身作为发起者的可用同步通信环路。此外,当一组共享信息向量到达时,需要对自身作为响应者的该同步通信环路进行记录。通过上述两种情况,每个卫星都会从自身的角度建立和更新同步通信环路。对于每个卫星,对所有的同步通信环路进行更新并按照最早开始时间的先后进行排序和执行。
为了实现本地局部收敛,在环路可用的时间区间内,每当环路发起者或响应者所构造的任务束改变时,则在环路上触发一次同步共享反馈交互,同时生成最后一次信息更新的时间戳。
同步通信环路中的一个共享反馈交互认定为一次通信,则在m-CBBA算法的在线协同调度中最多有un1(n1-1)次通信。在每次共享反馈交互中,进行以下两个阶段。
阶段2:束构建
相比于枚举所有可能的任务束,在m-CBBA中,每个对地观测卫星都构建了各自相应的任务束,并随着调度过程的进行而更新。在算法的这一阶段中,每个对地观测卫星以顺序性贪婪方式不断将任务添加到其任务束中,直到不能再添加任何任务为止。任务束中的任务是根据依次添加的顺序排列的,而执行计划中的任务是按照其对应的实际开始时间进行排序的。
束构建过程是基于面向单星的在线调度算法。对于当前不在任务束中的每个可用任务,对地观测卫星将其得分值与当前中标向量中的任务得分进行比较,如果更大,则保留为新的任务中标得分。对于未调度的任务集合,对地观测卫星从中选取具有最高得分的任务,并将其添加到任务束中,并对任务束,执行计划,时间戳,胜出者和中标向量进行更新,以包含新添加的任务。
束构建的递归过程一直持续到任务束已达到其容量限额,或者不再有任务能够添加到任务束中,即该卫星对于剩余任务的执行不再能够胜过其他卫星。
所述的单星在线调度算法在LEO对该单星的任务进行调度。单星在线调度算法调度该星的常规任务和协同分配算法分配给该星的应急观测任务。所述单星在线调度算法包括:
(1)在T-驱动的调度时刻点,采用渐进式方法中的完全重调度策略,生成下一个周期时间区间内的新任务计划,T-驱动的调度时刻点是根据给定的时间间隔T来确定特定的调度时间点lT,0≤l≤L,LT≤H<(L+1)T,每到达一个调度时间点lT,则计算生成后一调度区间[lT,(l+1)T]的任务计划,其中l为正整数,T为给定的时间间隔,L为最大T-驱动调度次数,H为总调度区间,以及
(2)在C*-驱动的重调度时刻点,采用修订式方法中的调度计划修复策略,当卫星运行在给定的调度区间内时,若在某一时刻t(0<t<H),星上的应急观测任务累积量Ct超过给定的阈值C*时,则执行重调度计算,其中阈值C*为应急观测任务的临界累积数,
除上述两种调度时刻点之外,不在任何其他时刻点进行调度。
更具体地,在所述单星在线调度算法中,在T-驱动的调度时刻点的具体调度算法如下:
输入:
输出:
具体步骤如下:
步骤13按照设定的启发式规则,对整合后的观测任务集合中的任务进行排序;
在C*-驱动的重调度时刻点的调度算法如下:
输入:
输出:
具体步骤如下:
阶段3:一致性构建
一旦对地观测卫星s构建起自身的任务束,则需要与其他对地观测卫星进行通信,以消解分配冲突。在接收到相邻对地观测卫星s传递过来的胜出者和相应的中标信息后,每个对地观测卫星可以确定其任务束中的任何任务是否已被其他对地观测卫星胜出并中标。
在一致性构建阶段中,每对相邻对地观测卫星s会同步共享以下信息向量:胜出者列表WAio,中标出价列表WBio和时间戳向量TSio,对应表示其他对地观测卫星s所接收到最新信息更新的时间戳。
对于在发送方i和响应方i′之间传递的每条消息,对地观测卫星i′会根据所接收到的信息,执行一系列操作,从而对其自身的信息向量进行更新。这些具体操作包括将其自身信息WAi′o,WBi′o和TSi′o,与对地观测卫星i所对应的信息进行比较,以确定对于每个任务哪些是最新的信息。对地观测卫星i′可对每个任务j采取三种可能的操作:
1.更新:wai′oj=waioj,wbi′oj=wbioj;
2.重置:wai′oj=φ,wbi′oj=0;
3.舍弃:wai′oj=wai′oj,wbi′oj=wbi′oj;.
对对地观测卫星i发送来的消息进行响应时,针对应急观测任务oj,对地观测卫星i′所采用的操作决策规则,如表4.1。
表.1同步通信下当接收到消息时的操作决策规则
如果在通信过程中胜出者列表WAio或中标出价列表WBio有所变动,则对地观测卫星检查所更新或重置的任务是否在自身的任务束中。如果存在,则任务束中这些任务及其后续的任务都需要清空释放,然后当是任务束中第一个被胜出的任务位置序号,则对于任务束中后续的所有任务及其任务脚标,会进行如下更新:
从而任务束被裁剪,以移除后续的任务:
同时,相对应的任务条目也从执行计划和时间戳向量中移除。
如果后续已调度任务不释放,则结果过于保守,会导致算法收益降低,则释放后续的已调度任务,重新进行束构建,是由任务调度安排的前后时间关联性所决定的。至此,算法重新进入第二阶段,将新的任务添加到任务束中。m-CBBA算法则在后两个阶段中进行迭代直至信息变量不再编队。在实际应用中,同一同步通信环路上的同步分享反馈次数会受到通信成本代价的限制。
基于异步通信的改进异步一致性束算法
异步通信应用于算法中各个模块之间相互独立运行的情况。因此,异步通信在分散式算法中广泛应用,即能够随时利用可用的信息,并不按照严格的调度节点。
在m-CBBA算法中,束构建阶段可以独立而异步地执行的,因为在此阶段对地观测卫星s只需要本地局部信息。然而,m-CBBA算法的一致性构建阶段采用同步通信策略会影响整个***的响应能力和性能表现,因为冲突消解策略要求每个对地观测卫星在每次通信迭代中知道每个相邻对地观测卫星的最新信息。可见,这是采用了全局一致性的方法,其中每个对地观测卫星分享自己的当前状态,并等待接收所有相邻对地观测卫星s的反馈消息(或者等待一定的时间),然后进入下一个算法阶段。
对于实时***而言,这样强制的算法延迟是不切合实际的,影响到算法的性能表现和收敛速率,因此本节提出m-ACBBA算法,采用了异步通信环路和异步冲突消解规则。
与m-CBBA算法相比,m-ACBBA不仅将中继节点视为通信节点,而且还采用中继节点发挥协同器的作用。因此,m-CBBA和m-ACBBA之间的不同之处在于m-ACBBA使用异步通信环路预测和异步局部决策规则。
阶段1:异步通信环路预测
当对地观测卫星在对新一批应急观测任务进行调度之后立即将其信息向量分享给中继节点,中继节点作为协同器,负责对异步通信环路进行预测。
异步子通信环路(ASCLs,Asynchronous sub-communication loops)是由单星自身的调度周期,应急观测任务的截止日期和中继节点与对地观测卫星之间的通信时间窗口共同决定的。具体的描述形式是一多元组<ID,I,R,ES,LF>ASCL,与同步通信环路的描述形式相同,但不同之处在于I或R为中继节点,当一对地观测卫星确定为I时,R则设定为中继节点。否则,当I为中继节点时,则R设定为对地观测卫星。可见,通信时间窗口ctwigq等同于异步子通信环路<id,i,g,cesigq,clfigq>或<id,g,i,cesigq,clfigq>。
异步通信环路(ACL)<ID,I,R,ESUS,LFUS,ESDS-UF,LFDS-UF,ESDF,LFDF>ACL是由两个或三个独立的异步子通信环路构成的。
图5(a)表示两个异步子通信环路<ID1,I1,R1,ES1,LF1>ASCL和<ID2,I2,R2,ES2,LF2>ASCL形成上述异步通信环路。在这两个异步子通信环路中,存在五种相关关系,具体如下:
(1)I和R均是对地观测卫星。ACL发起者I是前一ASCL的发起者I1,同时ACL响应者R是后一ASCL的响应者R2。
I=I1
R=R2
(2)前一ASCL的响应者R1和后一ASCL的发起者I2均是通信中继节点。当通信中继节点相互连通时,则两ASCL对应的中继节点可以不相同。否则,对应的中继节点是同一个。
(3)ESUS和LFUS表示ACL的分享上行通信时间窗口。
ESUS=ES1
LFUS=min{LF1,LF2}
ESUS≤LFUS
(4)ESDS-UF和LFDS-UF表示ACL的分享下行-反馈上行通信时间窗口。
ESDS-UF=max{ES1,ES2}
LFDS-UF=LF2
ESDS-UF≤LFDS-UF
(5)ESDF和LFDF表示ACL的反馈下行通信时间窗口。
ESDF=max{ES1,ES2}
LFDF=LF1
ESDF≤LFDF
图5(b)表示三个异步子通信环路ASCLs<ID1,I1,R1,ES1,LF1>ASCL,<ID2,I2,R2,ES2,LF2>ASCL和<ID3,I3,R3,ES3,LF3>ASCL,按照最早开始时间的先后顺序,组成一异步通信环路ACL。同样,这三个ASCLs中存在五个相关关系,具体如下:
(1)I和R均是对地观测卫星。ACL的发起者I是第一个ASCL的发起者I1,同时,响应者R是最后一个ASCL的响应者R3。
I=I1
R=R3
(2)第一个ASCL的发起者I1与最后一个ASCL的响应者R3相同,以构成一完整通信环路。
I1=R3
(3)第一个ASCL的响应者R1和后两个ASCL的发起者I2和I3均是通信中继节点。当中继节点之间相互连通时,则相对应的中继节点可以不同,否则,是同一个中继节点。
(4)ESUS和LFUS表示ACL的分享上行通信时间窗口。
ESUS=ES1
LFUS=min{LF1,LF2}
ESUS≤LFUS
(5)ESDS-UF和LFDS-UF表示ACL的分享下行-反馈上行通信时间窗口。
ESDS-UF=max{ES1,ES2}
LFDS-UF=LF2
ESDS-UF≤LFDS-UF
(6)ESDF和LFDF表示ACL的反馈下行通信时间窗口。
ESDF=max{ES2,ES3}
LFDF=LF3
ESDF≤LFDF
异步分享反馈交互是在可用的异步通信环路上触发生成的。在该交互中,共有四个时间戳,包括分享上行,分享下行,反馈上行和反馈下行。其描述形式<id,i,i′,tUS,tDS,tUF,tDF>源于异步通信环路<ID,I,R,ESUS,LFUS,ESDS-UF,LFDS-UF,ESDF,LFDF>ACL,其中i∈I,i′∈R,tUS∈[ESUS,LFUS],tDS,tUF∈[ESDS-UF,LFDS-UF],tDF∈[ESDF,LFDF],tUS≤tDS,tDS<tUF,tUF≤tDF。
对于每个对地观测卫星,在ACL的可用时间区间内,当发起者或响应者的任务束发生变动时,则异步分享反馈交互在环路上执行。由于通信延迟和调度计算所占用的时间相对较短,则将分享下行与反馈上行的时间戳等同,tDS≈tUF。
ACL上的一次异步分享反馈交互共包括三次通信,即分享上行,分享下行-反馈上行和反馈下行。对于一批次的应急观测任务,采用m-ACBBA算法进行在线协同调度,最多需要3un1(n1-1)次通信。
在每次异步分享反馈交互中,束构建和一致性构建均执行。与m-CBBA算法相比,不同之处在于m-ACBBA算法在一致性构建阶段,采用异步局部决策规则。
异步通信协议中消息的到达是乱序的,比如较早时间生成的消息会比较晚时间生成的信息晚到。因此,m-CBBA算法中对地观测卫星基于消息接收时间的时间戳更新机制不适合这种异步情况。则m-ACBBA算法的配置中,用消息上的时间戳表示胜出对地观测卫星s中标的实际时刻,即TSio@{tsio1,...,tsiou},长度为u,其中tsioj表示针对第o批次中应急观测任务j,对地观测卫星i上这一胜出标的生成时间。而不是m-CBBA算法中所采用的对地观测卫星s消息更新时间,即长度为n1,其中tsioi′表示针对第o批次应急观测任务,对地观测卫星i进行最新信息更新的时间戳。
在m-ACBBA算法中,采用一系列局部消解规则,而不需要全局状态信息。同时,该消解协议能够对乱序的消息和冗余信息进行一致化处理。
mACBBA算法中的冲突消解规则,不仅是接收者对胜出者列表和中标出价列表进行更新,而是确定需要进行再次分享的消息。
进行再次分享的原因,在于:一是减少网络中的通信负担,避免冗余信息的再次分享;二是处理异步***中时间戳的混乱。
接收者可采用以下五种可能的行动:
1.更新&再分享:接收者根据发送者过来的信息,更新胜出者列表,中标出价列表以及胜出时间向量,然后再分享这次更新;
2.舍弃&再分享:接收者不改变其信息状态,然后对自身所认定的胜出者信息进行再分享;
3.舍弃&不再分享:接收者不改变自身的信息状态,也不对其再分享,应用于所接收到的消息与现有的信息相一致的情况。
4.重置&再分享:接收者对胜出者和胜出值进行清零,然后对原接收到的消息进行再分享,则消除混乱;
5.更新时间&再分享:接收者是胜出者,当接收到消息时,对其胜出生成时间戳更新为当前时间,以确认自己仍然认为自己是胜出者。
对对地观测卫星i发送来的消息进行响应时,针对应急观测任务oj,对地观测卫星i′所采用的操作决策规则,如表2。
表2异步通信下当接收到消息时的局部决策规则
对地观测任务的目标是监测发现,详查和跟踪森林火灾或火山爆发活动。更准确地说,森林火灾和火山爆发需要自动监测发现,定位和识别。在这些应急事件发生的情况下,我们着重提高***的快速响应能力,并通过在线协同提高***的全局收益。
我们利用地球静止轨道(GEO)通信卫星作为通信中继节点,在一定程度上缓解了通信短缺的情况,但是通信条件仍然受限,约束仍然苛刻。首先,GEO和LEO之间的通信时间窗口比地面站和对地观测卫星之间的通信时间窗口要大得多;其次,GEO与其覆盖范围内的地面站时时相连,从而保证任务需求上传和数据下传的实时可用性,所以以这种通信方式,多星在线协同变为现实,可支持更高水平上的灾难管理。图6给出了空间段和地面段之间可用的通信情况,具体如下,
(1)当LEO通过GEO的覆盖区域时,对地观测卫星与地球静止轨道通信卫星之间的通信可用。该通信链路用于通过GEO,实现LEO和地面站之间的应急观测任务双向传输,并将低速率数据从LEO下传到地面站(GS)。
(2)GEO与GS之间的通信交互是可以时时进行的,因为GS位于GEO的覆盖区域。该通信链路用于应急观测任务上传和低速率数据下载。
(3)GEO之间的双向通信交互是可行的,实现在GEO之间传输与任务有关的信息。
(4)对地观测卫星之间无法进行直接通信。
应用场景中,在同一轨道上配置3颗对地观测卫星,在赤道上方配置3颗地球静止轨道通信卫星,具体卫星参数见表4.3。仿真的时间长度设定为6h。
表.3三颗对地观测卫星和三颗地球静止轨道通信卫星的轨道参数设置
该配置场景的特点如下:
1.当同一轨道上的三颗对地观测卫星对地面上同一目标进行观测成像时,LEO_A和LEO_B之间的时间间隔约为12分钟,LEO_B与LEO_C之间的时间间隔也是12分钟。
2.地球静止轨道通信卫星的通信覆盖区域为一锥角为7.5°的圆锥形区域。
3.对地观测卫星的观测成像区域为一矩形。若在同一轨上对同一目标进行观测成像,则LEO_A和LEO_B之间成像角度的变化值约为20°,同样LEO_B与LEO_C之间成像角度的变化值也约为20°。
对当前分布式卫星***的通信情况进行分析,包括面向集中-分布式架构的高低轨卫星之间通信情况,如表4,和面向分散式架构的低轨之间通信情况,如表5。
表4当前分布式卫星***中高低轨卫星之间的通信约束情况
表5当前分布式卫星***中对地观测卫星之间的通信约束情况
我们对比了面向分布式卫星***的四种在线协同调度算法,包括单项任务下的合同网协议算法SI-CNP,批次任务下的合同网协议算法BA-CNP,改进的一致性束算法m-CBBA和改进的异步一致性束算法m-ACBBA。
在m-CBBA和m-ACBBA算法中,每次共享反馈交互的分享上行,分享下行,反馈上行和反馈下行时间戳分别是从通信环路对应的最早开始时间和最晚结束时间之间的时间区间内随机生成。在每个完全通信环路上,最多进行一次共享反馈交互。
我们采用三项性能指标,分别是总收益,已调度成功的应急观测任务比例和相应的通信总次数。这三项指标的具体描述如下:
(1)***总收益,指分布式卫星***在整个调度区间内的所有已调度成功的任务收益之和。
(2)已调度成功的应急观测任务比例,指分布式卫星***中已调度成功的应急观测任务在进入调度求解的应急观测任务中所占的百分比。
(3)通信总次数。通信总次数是是三种单向通信情况的通信次数之和,分别是LEO向GEO发起通信,GEO向LEO发起通信以及GEO之间单向通信。
为了评估四种在线协同调度算法的性能,我们重点关注分布式卫星***对应急观测任务到达的响应性,因此两个与任务相关的参数量,应急观测任务的到达速率和相机指向角度,以及一个与调度机制相关的参数量,滚动调度周期,对于测试算例的生成是十分重要的。可见,这两个与任务相关的参数量分别代表应急观测任务的时间分布特性和空间分布特性,而滚动调度周期则决定了在线调度计算的规模大小。此外,其他所有参数量均由给定的均匀分布生成。与卫星、任务和调度机制相关的参量设置情况,如表6所示。
表6相关参量设置
图7给出了不同滚动调度周期下分别针对地球静止轨道通信卫星之间不互联与互联情况的m-CBBA算法结果。与地球静止轨道通信卫星之间不互联的情况相比,互联互通的地球静止轨道通信卫星网络有助于实现更高的总收,完成更高比例的应急观测任务。因为与不互联的情况相比,互联互通使得通信时间窗口数量增加或时间区间延长。同时,互联互通需要更多的通信,特别是GEO之间的通信交互。
表7显示随着应急观测任务到达速率增加,每个算法所取得的总收益也在增加。在应急观测任务到达速率相同的情况下,应急观测任务的空间分布范围越广,所取得的总收益越低,调度成功的应急观测任务比例越低。
当滚动调度周期T为6min时,在应急观测任务的空间和时间分布相同时,m-ACBBA和m-CBBA都比SI-CNP和BA-CNP获得更高的总收益,而m-ACBBA的总收益通常高于m-CBBA。而当滚动调度周期T为12min时,m-ACBBA在四种算法中得到最高的总收益,而SI-CNP算法所得到的收益高于BA-CNP和m-CBBA算法。可见,对于BA-CNP和m-CBBA算法,滚动调度周期越长,总收益越低,因为BA-CNP算法针对批次任务进行调度,而每个调度周期对应于一个应急观测任务批次,同时m-CBBA算法在每个滚动调度周期内的同步通信环路上最多进行一次分享反馈交互,则滚动调度周期越长,分享反馈次数越少,从而影响应急观测任务的及时响应调度。
对于SI-CNP和m-ACBBA算法,滚动调度周期越长,总收益越高,特别是m-CBBA算法在通信次数减少的情况,总收益仍有所增长,说明异步通信策略表现优异。
表7各算法在不同参数条件下所取得的总收益
表8表明在相同的参数量下,m-ACBBA在四种算法中调度成功的应急观测任务比例最高。整体上按调度成功的应急观测任务比例由高到低的方式,对各算法进行排序,依次是:m-ACBBA>m-CBBA>BA-CNP>SI-CNP。
表8各算法在不同参数条件下所取得的应急观测任务调度成功比例(%)
表9显示随着应急观测任务的到达速率增加,通信数量有所增加,而SI-CNP的通信量增长幅度最大,因为其依次对单个任务进行协同调度。在这四种算法中,m-CBBA的通信次数最少,SI-CNP的通信次数最多,因为同步通信环路的构建条件严格,则环路数量少,进而影响通信交互次数。随着滚动调度周期增长,SI-CNP对应的通信量大大增加,而其他三种算法对应的通信量均有所减少。
表9各算法在不同参数条件下的通信次数
表10表明三种算法相比SI-CNP在总收益上均有提高,其中m-ACBBA的收益增加幅度最大。同样,表11表明三种算法相比SI-CNP算法在应急观测任务调度成功的比例上均有提高,其中m-ACBBA的提高幅度最大。
表12显示随着应急观测任务的到达速率增加,平均计算时间增加。在相同的参数条件下,这四种算法所需的计算时间几乎相同,且均与星上实际计算能力相匹配。
表10当T为6min时与SI-CNP相比,各算法的收益增加(%)
表11当T为6min时与SI-CNP相比,各算法的应急观测任务调度成功增量(%)
表12一颗对地观测卫星上一调度周期内调度平均计算时间(/s)
综上所述,当***中的通信成本代价高时,m-CBBA算法可在***总收益和通信次数之间取得平衡,而当***的通信成本低时,m-ACBBA算法是获得高***总收益和高应急观测任务调度成功比例的最佳选择。
本发明对严格通信约束下面向简单任务的同构分布式对地观测卫星***在线协同调度问题进行研究。面向简单任务的多星协同调度问题,均属于ID[ST-SR-TA]问题。因为单星上可以分配并调度多个任务,而且由于时间窗口约束和资源(电量,固存)限制,针对某个任务的单星调度,会受到其他任务的影响,则存在调度计划内部的依赖关系。首先针对通信时间窗口和观测时间窗口,定义应急观测任务的批次可用性和时间有效性,构建每批应急观测任务到达时的子问题MILP模型;然后,以第三章提出的单星在线调度机制为基础,分别面向集中-分布式架构和分散式架构,各提出两种在线协同任务调度算法。具体而言,面向集中-分布式架构的两种算法是单项任务下的合同网协议算法SI-CNP和批次任务下的合同网协议算法BA-CNP;面向分散式架构的两种算法是改进一致性束算法m-CBBA和改进异步一致性束算法m-ACBBA,两者分别是基于同步通信和异步通信。计算实验结果表明,相比不互联互通,中继节点的互联互通会使得通信时间窗口有所增加和扩展,有助于***总收益和应急观测任务调度成功比例的增加。当***中的通信成本高时,m-CBBA算法在平衡利润和所需通信数量上占优。当***中的通信成本代价高时,m-CBBA算法可在***总收益和通信次数之间取得平衡,而当***的通信成本低时,m-ACBBA算法是获得高***总收益和高应急观测任务调度成功比例的最佳选择。
最后需要指出的是:以上实施例仅用以说明本发明的技术方案,而非对其限制。本领域的普通技术人员应当理解:可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (6)
1.一种同构多星在线协同方法,其特征在于,对地观测卫星作为资源观测卫星,地球静止轨道通信卫星作为通信中继节点,对地观测卫星与地球静止轨道通信卫星通信,对地观测卫星之间不能直接进行卫星通信,所述同构多星在线协同方法采用集中-分布式协同架构,协同决策功能仅配置在地球静止轨道通信卫星上,而对地观测卫星进行分布式计算,所述同构多星在线协同方法采用基于拍卖机制的任务协同分配算法,拍卖由作为拍卖商的地球静止轨道通信卫星组织进行,每个对地观测卫星根据自身的本地信息来计算与应急观测任务相对应的投标出价值,并以此进行投标,每个对地观测卫星将其出价通知拍卖商,然后地球静止轨道通信卫星决定每个应急观测任务的胜出者,如果有多个通信中继节点相互关联,则拍卖商的角色能够在这些通信中继节点之间进行转移,
所述任务协同分配算法包括两个阶段的迭代:束构建阶段和一致性构建阶段,其中前者对应于每个对地观测卫星以贪婪搜索的形式生成排序后的任务束,后者对应于识别任务分配冲突,并通过与相邻对地观测卫星之间的局部通信进行冲突消解,上述的两个阶段不断地迭代重复,直至达到收敛,
对地观测卫星与地球静止轨道通信卫星之间的通信是间歇性的,所述任务协同分配算法进一步包括一个同步通信环路预测阶段,其中,同步通信环路是用于两个或多个对地观测卫星的双向通信链路,以向其他对地观测卫星共享信息向量并从环路中接收来自其他对地观测卫星的反馈信息,一个同步通信环路构建于发起者和一个或多个响应者的通信时间窗口之间均存在重叠时间的基础上,同步通信环路的描述形式是一个多元组<ID,I,R,ES,LF>SCL:
—ID是一个标识符;
—I是通信发起者编号;
—R是接收并分享信息的通信响应者编号;
—ES是通信环路的最早开始时间;
—LF是通信环路的最晚结束时间,
其中通信发起者I和响应者R都是对地观测卫星,并且ES和LF是由发起者和所有响应者的通信时间窗口上的重叠时间区间确定的,
同步通信环路预测阶段用于预测同步通信环路的存在,且提供所述同步通信环路的多元组<ID,I,R,ES,LF>SCL,同步通信环路预测阶段是束构建阶段的前提。
2.如权利要求1所述的同构多星在线协同方法,其特征在于,所述地球静止轨道通信卫星采用下述的中标决策方式:
1).如果有多个对地观测卫星能够在截止时间之前完成应急观测任务,则选择使***总收益增量最大的对地观测卫星作为中标者;如果有多个卫星的投标出价值对应于最高的总收益增量,则以实际开始观测时间最早的对地观测卫星具有中标优先权;
2).如果只有一个对地观测卫星可以在截止日期之前完成应急观测任务,同时***总收益的增量为正值,则选择该卫星作为中标者;
3).如果当对地观测卫星执行该应急观测任务产生的***总收益增量均是负值,或者没有对地观测卫星可以在截止时间之前完成应急观测任务,则地球静止轨道通信卫星舍弃该应急观测任务。
3.如权利要求1所述的同构多星在线协同方法,其特征在于,以下式计算应急观测任务的时间可用性wioj,计算如下:
其中,olsioj为第o批次中第j个应急观测任务的最晚观测开始时间,且olsioj=max{olfioj-poj-(2·maxθ)/sproll,oesioj},
其中,
olfioj为针对对地观测卫星i,第o批次中第j个应急观测任务的最晚观测结束时间,
poj为第o批次中第j个应急观测任务的成像时长,
maxθ为对地观测卫星的最大侧摆角,
sproll为对地观测卫星的侧摆角速度,
oesioj为针对对地观测卫星i,第o批次中第j个应急观测任务的最早观测开始时间
trio为第o批次应急观测任务的时间到达对地观测卫星i的时间
wio为第o批次应急观测任务的时间可用性wio。
4.如权利要求1所述的同构多星在线协同方法,其特征在于,每颗对地观测卫星具有六个信息向量,具体如下:
1.)任务束,1≤n≤|BUio|,该任务束表示对地观测卫星i从第o批应急观测任务中选取并调度成功的任务集合,同时按照各任务添加到束中的时间先后进行排序,当前任务束的长度为|BUio|,且要小于一批次内的应急观测任务数量u,|BUio|≤u;当任务束为空时,用BUio=φ和|BUio|=0表示,
4.)胜出者列表WAio@{waio1,...,waiou},长度为u,其中waioj表示对地观测卫星i当前所认为的针对第o批应急观测任务中第j个任务的胜出者,具体取值与胜出者编号相对应,当waioj=φ时,对地观测卫星i认为当前对于该任务无胜出者,
5.)中标出价列表WBio@{wbio1,...,wbiou},长度为u,其中wbioj表示对应的胜出者所给出的出价,当取值为0时,表示该任务当前无胜出者,
5.如权利要求1所述的同构多星在线协同方法,其特征在于,对地观测卫星上执行单星在线调度算法,调度该对地观测卫星的常规任务和上述协同分配算法分配给该对地观测卫星的应急观测任务,所述单星在线调度算法包括:
(1)在T-驱动的调度时刻点,采用渐进式方法中的完全重调度策略,生成下一个周期时间区间内的新任务计划,T-驱动的调度时刻点是根据给定的时间间隔T来确定特定的调度时间点lT,0≤l≤L,LT≤H<(L+1)T,每到达一个调度时间点lT,则计算生成后一调度区间[lT,(l+1)T]的任务计划,其中l为正整数,T为给定的时间间隔,L为最大T-驱动调度次数,H为总调度区间,以及
(2)在C*-驱动的重调度时刻点,采用修订式方法中的调度计划修复策略,当卫星运行在给定的调度区间内时,若在某一时刻t(0<t<H),星上的应急观测任务累积量Ct超过给定的阈值C*时,则执行重调度计算,其中阈值C*为应急观测任务的临界累积数,
除上述两种调度时刻点之外,不在任何其他时刻点进行调度。
6.如权利要求5所述的同构多星在线协同方法,其特征在于,在所述单星在线调度算法中,在T-驱动的调度时刻点的具体调度算法如下:
输入:
输出:
具体步骤如下:
步骤13 按照设定的启发式规则,对整合后的观测任务集合中的任务进行排序;
在C*-驱动的重调度时刻点的调度算法如下:
输入:
输出:
具体步骤如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711279412.3A CN108023637B (zh) | 2017-12-06 | 2017-12-06 | 一种同构多星在线协同方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711279412.3A CN108023637B (zh) | 2017-12-06 | 2017-12-06 | 一种同构多星在线协同方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108023637A CN108023637A (zh) | 2018-05-11 |
CN108023637B true CN108023637B (zh) | 2020-06-23 |
Family
ID=62078722
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711279412.3A Expired - Fee Related CN108023637B (zh) | 2017-12-06 | 2017-12-06 | 一种同构多星在线协同方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108023637B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115173926A (zh) * | 2022-07-08 | 2022-10-11 | 中国科学院微小卫星创新研究院 | 基于拍卖机制的星地融合中继网络的通信方法和通信*** |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109670678B (zh) * | 2018-11-29 | 2021-05-18 | 合肥工业大学 | 虚拟星座动态环境下的任务分配方法和*** |
CN110422341A (zh) * | 2019-07-31 | 2019-11-08 | 南京航空航天大学 | 一种用于火星探测的火星车系留飞艇***及其工作方法 |
CN111077770B (zh) * | 2019-10-29 | 2022-05-24 | 合肥工业大学 | 多星协同资源的配置方法和*** |
CN110851926B (zh) * | 2019-10-31 | 2023-09-01 | 武汉科技大学 | 一种基于遍历环路的行星轮系同构判定方法、***和介质 |
CN111124637B (zh) * | 2019-12-09 | 2023-08-29 | 中国空间技术研究院 | 一种分布式卫星***动态任务规划方法 |
CN111898098B (zh) * | 2020-08-18 | 2021-07-13 | 哈尔滨工业大学 | 面向多方协作的利益相关者服务价值冲突发现和消解方法 |
CN114906353A (zh) * | 2020-09-29 | 2022-08-16 | 中国科学院微小卫星创新研究院 | 导航卫星空间态势感知方法 |
CN115081884B (zh) * | 2022-06-23 | 2023-05-02 | 哈尔滨工业大学 | 一种分布式星上在线多对多任务规划方法 |
CN115616902B (zh) * | 2022-11-07 | 2023-03-07 | 中国人民解放军国防科技大学 | 基于鲁棒拍卖算法的集群航天器任务分配方法及装置 |
CN115535297B (zh) * | 2022-11-30 | 2023-03-28 | 南京航空航天大学 | 一种卫星星座在轨分布式协同调度方法 |
CN115913340B (zh) * | 2023-01-06 | 2023-05-02 | 武汉大学 | 无人机辅助通信方法、装置、设备及存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102780523A (zh) * | 2012-07-31 | 2012-11-14 | 中国电子科技集团公司第五十四研究所 | 一种多卫星协同观测业务调度方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020077099A1 (en) * | 2000-12-18 | 2002-06-20 | Space Systems/Loral, Inc. | Method and system for providing satellite communications |
-
2017
- 2017-12-06 CN CN201711279412.3A patent/CN108023637B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102780523A (zh) * | 2012-07-31 | 2012-11-14 | 中国电子科技集团公司第五十四研究所 | 一种多卫星协同观测业务调度方法 |
Non-Patent Citations (3)
Title |
---|
分布式卫星自主构型重构技术研究;张健;《中国优秀博硕士学位论文全文数据库》;20071115;全文 * |
基于Agent的对地观测卫星分布式协同任务规划研究;王冲;《中国博士学位论文全文数据库》;20140415;正文第3-4页、第19-45页,附图1.1、附图2.2 * |
基于合同网的任务协同研究;李志翔;《通讯世界》;20161130;第10-11页 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115173926A (zh) * | 2022-07-08 | 2022-10-11 | 中国科学院微小卫星创新研究院 | 基于拍卖机制的星地融合中继网络的通信方法和通信*** |
CN115173926B (zh) * | 2022-07-08 | 2023-07-07 | 中国科学院微小卫星创新研究院 | 基于拍卖机制的星地融合中继网络的通信方法和通信*** |
Also Published As
Publication number | Publication date |
---|---|
CN108023637A (zh) | 2018-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108023637B (zh) | 一种同构多星在线协同方法 | |
CN109088667B (zh) | 一种面向简单任务的同构多星在线协同方法 | |
CN111912412B (zh) | 面向应用的异构星座天地一体化任务规划方法及装置 | |
Wang et al. | Multi-resource coordinate scheduling for earth observation in space information networks | |
CN108510145B (zh) | 一种通信约束下面向复合任务的异构多星在线协同方法 | |
Du et al. | A new multi-satellite autonomous mission allocation and planning method | |
CN112580906A (zh) | 卫星遥感任务规划与地面资源调度联合求解方法 | |
Augenstein et al. | Optimal scheduling of a constellation of earth-imaging satellites, for maximal data throughput and efficient human management | |
WO2019127946A1 (zh) | 一种基于学习型遗传算法的多任务多资源滚动分配方法 | |
CN113612525A (zh) | 基于约束满足的低轨互联网星座星地链路规划方法 | |
CN112183928B (zh) | 一种适用于遥感卫星成像的调度*** | |
Zheng et al. | Onboard mission allocation for multi-satellite system in limited communication environment | |
Li | Online scheduling of distributed Earth observation satellite system under rigid communication constraints | |
Yang et al. | Onboard coordination and scheduling of multiple autonomous satellites in an uncertain environment | |
CN109347536A (zh) | 一种基于态势知识的空间网络资源状态监控*** | |
Picard et al. | Autonomous agents and multiagent systems challenges in earth observation satellite constellations | |
CN108377163B (zh) | 一种基于异步通信的多星在线协同方法 | |
Schmidt | Ground station networks for efficient operation of distributed small satellite systems | |
CN115545413A (zh) | 一种基于约束合同网的多星自主任务规划方法 | |
CN112423346A (zh) | 一种中继多址接入调度方法、装置、电子设备及存储介质 | |
Cahoy et al. | Initial results from ACCESS: an autonomous cubesat constellation scheduling system for earth observation | |
Kennedy | Planning and scheduling for earth-observing small satellite constellations | |
CN109711702A (zh) | 面向多管理平台的对地观测任务分布式协同规划方法 | |
CN116862167A (zh) | 基于多背包模型的低轨遥感卫星星座应急任务规划方法 | |
Krupke et al. | Automated data retrieval from large-scale distributed satellite systems |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200623 Termination date: 20201206 |