CN111884703B - 一种基于通信卫星间协同计算的业务请求分配方法 - Google Patents

一种基于通信卫星间协同计算的业务请求分配方法 Download PDF

Info

Publication number
CN111884703B
CN111884703B CN202010567261.7A CN202010567261A CN111884703B CN 111884703 B CN111884703 B CN 111884703B CN 202010567261 A CN202010567261 A CN 202010567261A CN 111884703 B CN111884703 B CN 111884703B
Authority
CN
China
Prior art keywords
service request
omega
mapping
service
wteg
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
Application number
CN202010567261.7A
Other languages
English (en)
Other versions
CN111884703A (zh
Inventor
张亚生
任智源
窦志斌
程文驰
何辞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
CETC 54 Research Institute
Original Assignee
Xidian University
CETC 54 Research Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xidian University, CETC 54 Research Institute filed Critical Xidian University
Priority to CN202010567261.7A priority Critical patent/CN111884703B/zh
Publication of CN111884703A publication Critical patent/CN111884703A/zh
Application granted granted Critical
Publication of CN111884703B publication Critical patent/CN111884703B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/18521—Systems of inter linked satellites, i.e. inter satellite service
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/90—Details of database functions independent of the retrieved data types
    • G06F16/901—Indexing; Data structures therefor; Storage structures
    • G06F16/9024—Graphs; Linked lists
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/004—Artificial life, i.e. computing arrangements simulating life
    • G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • 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/18578—Satellite systems for providing broadband data service to individual earth stations
    • H04B7/18584—Arrangements for data networking, i.e. for data packet routing, for congestion control
    • 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/195—Non-synchronous stations
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Astronomy & Astrophysics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Computing Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Molecular Biology (AREA)
  • Computational Linguistics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • Biophysics (AREA)
  • Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Radio Relay Systems (AREA)

Abstract

本发明公开的基于通信卫星间协同计算的业务请求分配方法，涉及通信技术领域，根据预设的基于业务请求的有向无环图DAG及基于低轨通信卫星网络的加权时间扩展图WTEG之间的映射规则，生成业务请求对应的多个映射结果，根据二值粒子群优化算法，从多个映射结果中获取最优的映射结果，根据映射结果，将业务请求分配至各个通信卫星，有效地降低了卫星间协作计算的业务请求处理时延。

Description

一种基于通信卫星间协同计算的业务请求分配方法
技术领域
本发明涉及通信技术领域，具体涉及一种基于通信卫星间协同计算的业务请求分配方法。
背景技术
随着信息技术的不断发展，涌现了如远程手术，实时交通状况监控，无人驾驶等一系列时延敏感的新型业务模式，信息连接从传统的人-人互联，正朝向物-物互联的境界发展，地理的概念在信息层面逐渐变得模糊，用户随时随地享受新型信息业务成为可能。然而，当前地面移动通信网络有限的覆盖现状与和未来无处不在的信息服务需求之间存在巨大的鸿沟，因此，非地面通信网络(Non-terrestrial Networks，NTN)成为当前移动通信领域的研究热点。其中，卫星通信网络具备覆盖范围广、抗毁性强、全天候服务、不受地理环境约束等众多优点，能够有效补齐移动通信网络覆盖范围有限的短板，特别是低轨通信卫星网络(Low earth orbit-communication satellite network,LEO-SCN)凭借星地传输距离短、建造成本低的独特优势获得了广泛的关注，正逐步从概念构想迅速转向实质建设。目前，3GPP已将NTN纳入5G标准讨论稿；波音、空客、Facebook、SpaceX等高科技企业相继投资了OneWeb、Starlink等十余个LEO-SCN，以期实现全球的通信网络无缝覆盖；同期，中国航天科技集团和中国航天科工集团分别提出了“鸿雁”和“虹云”LEO-SCN，希冀引领中国进入LEO卫星通信时代。可以看出，将LEO-SCN作为移动通信网络的补充，实现深山、沙漠、海洋乃至外太空等区域的泛在覆盖已成为主流趋势。但是，随着提供的信息服务类型不断拓展，LEO-SCN在支撑时延敏感的新型信息业务时产生了新的问题。当前已部署和建设中的LEO-SCN均采用弯管式透明转发通信卫星，用户的请求数据统统经过LEO-SCN转发并在地面云中心服务器进行集中处理。这将导致海量用户数据经过LEO-SCN传输到地面站，给LEO-SCN带来巨大的传输压力；同时，频繁的星地链路长距离传输产生较高的业务传输时延，难以满足新型信息业务的低时延需求。因此，如何利用低轨通信卫星进行在轨计算成了研究的热点问题。
为了将业务部署在低轨通信卫星上处理，研究者们提出了一些优化单星服务能力的方案。硬件方面，针对机载处理器更新难度大、成本高的问题，Hofmann等人基于现场可编程门阵列技术提出了一种机载处理平台，能够在轨调整服务功能，以满足不同业务的处理需求；面向业务内在机理，Jacobs等人设计了一个中央处理器和现场可编程门阵列相结合的机载混合处理器架构并给出了高光谱图像监测与分析的案例，通过部署服务功能到最兼容的处理器模块上，提高低轨通信卫星的计算性能。软件优化方面，Wang等人提出了一种面向多终端用户的星载资源动态分配方案，最大化单颗卫星的资源利用率，但是上述方案都受到卫星载荷能力制约，单颗卫星越来越难以适应日益增长的业务请求数量。
发明内容
为解决现有技术的不足，本发明实施例提供了一种基于通信卫星间协同计算的业务请求分配方法，该方法包括以下步骤：
根据预设的基于业务请求的有向无环图DAG及基于低轨通信卫星网络的加权时间扩展图WTEG之间的映射规则，生成业务请求对应的多个映射结果；
根据二值粒子群优化算法，从所述多个映射结果中获取最优的映射结果；
根据所述映射结果，将所述业务请求分配至各个通信卫星。
优选地，所述映射规则的创建过程包括：
采用表达式Φ＝(Ω,Γ)作为有向无环图DAG的数学模型，其中，Ω＝{ω1,...,ωl}为业务请求节点集合，构成Φ的业务请求集群；Γ为边集合，表示业务请求间的依赖关系；Φ为单起点、单终点的数学模型，即ω1为业务请求起点，ωl为业务请求终点，其余为中间业务请求节点；其中，业务请求ωi由元组{εi,ηi,Di}表征，其中，εi为计算复杂度系数，表征执行业务请求ωi耗费计算资源的程度；ηi为数据缩放因子参数，反映执行业务请求ωi过程中对输入数据量放大与缩小的过程；Di表示ωi的输入数据量，定义Θ↑(ωi)＝{ωj|(ωj,ωi)∈Γ}表示业务请求ωi的先驱节点集合，则Di可通过式
Figure BDA0002548295920000031
求得，其中，Dj和ηj属于表征服务功能ωj的元组；
采用表达式N＝(NT,ET,Graph)作为加权时间扩展图WTEG的数学模型，其中，NT为加权时间扩展图WTEG的节点集合，ET为加权时间扩展图WTEG的边集合，Graph为加权时间扩展图WTEG的时延权重组合矩阵，其中，
Figure BDA0002548295920000041
Figure BDA0002548295920000042
表示相邻时隙q与q+1间所有边的权重，
Figure BDA0002548295920000043
表示时隙q内所有边的权重；
制定Ω至NT的映射规则，包括：
根据表达式
Figure BDA0002548295920000044
将业务请求起点映射至业务请求发起通信卫星，将业务请求终点映射至结果接收通信卫星；将中间业务请求节点映射至任意通信卫星，由于业务请求之间不可分割，故每个业务请求都需要分配到单颗低轨通信卫星上计算，但受计算过程影响，存在业务请求被缓存到下一时隙的现象，此时业务请求经过虚拟缓存链路到达下一时隙的副本节点，因此，任意业务请求ωi∈Ω会被映射至NT中同一通信卫星的ρi+1个副本节点，ρi为自然数，ρi为由业务请求ωi计算引起的跨时隙数目；
制定Γ至P的映射规则，包括：
根据表达式Z((ωi,ωj))＝Path(B(ωi)e,B(ωj)s)，将有向无环图DAG任意边
Figure BDA0002548295920000051
映射为加权时间扩展图WTEG节点集合B(ωi)中最后一个节点B(ωi)e到加权时间扩展图WTEG节点集合B(ωj)第一个节点B(ωj)s之间的最短路由用Path(B(ωi)e,B(ωj)s)表示，路由中包含的边序列即为边(ωi,ωj)的映射结果，其中，鉴于低轨通信卫星网络的动态特性，加权时间扩展图WTEG的部分节点之间无连通路由，因此，(ωi,ωj)的映射需满足表达式
Figure BDA0002548295920000052
优选地，根据二值粒子群优化算法，求解表达式T(x)的最小值，得到所述多个映射结果中最优的映射结果包括：
假设粒子群规模为U，同时I表示群体最大迭代次数，则第i次迭代过程中第u个粒子的位置和速度表达式分别如式(1)及式(2)所示：
Figure BDA0002548295920000053
Figure BDA0002548295920000054
式(1)中，
Figure BDA0002548295920000055
式(2)中，
Figure BDA0002548295920000056
当i<I时，粒子u首先通过跟踪其历史最优位置pubest和整个粒子群的历史最优位置gbest更新其速度，如式(3)所示，
Figure BDA0002548295920000057
其中，μ为惯性权重，γ1、γ2为加速因子，β1、β2为均匀分布于区间[0,1]的随机数；随后跟踪速度
Figure BDA0002548295920000058
方向按照式(4)更新其位置，
Figure BDA0002548295920000059
式(4)中，
Figure BDA00025482959200000510
表示数值向上取整运算，
Figure BDA00025482959200000511
表示两个矩阵的元素之间的置换操作,具体过程是找出
Figure BDA0002548295920000061
中数值最大的元素索引作为待置换位置，根据式的约束条件
Figure BDA0002548295920000062
计算
Figure BDA0002548295920000063
中与待置换位置进行置换操作的元素索引，完成置换操作，通过式
Figure BDA0002548295920000064
计算其适应度值；
将适应度值定义为当前粒子映射结果的业务请求处理总时延，将表达式
Figure BDA0002548295920000065
分为最后一项业务请求的计算时延
Figure BDA0002548295920000066
和前项业务请求的累积处理时延
Figure BDA0002548295920000067
运用递归原理，最终计算得到粒子的适应度值，即所述多个映射结果中最优的映射结果，其中，
Figure BDA0002548295920000068
表示ωi在加权时间扩展图WTEG中映射的最后一个节点，
Figure BDA0002548295920000069
表示ωl在加权时间扩展图WTEG中映射的第一个节点，T(X)为业务请求处理总时延。
本发明实施例提供的基于通信卫星间协同计算的业务请求分配方法具有以下有益效果：
利用加权时间扩展图WTEG模型，屏蔽了LEO-SCN的动态特性，从而将基于业务请求的有向无环图DAG驱动的协同计算技术引入到LEO-SCN，达到降低业务请求处理时延的目的。
附图说明
图1为典型的空基云雾网络架构示意图；
图2为本发明实施例提供的基于通信卫星间协同计算的业务请求分配方法流程示意图；
图3a为DAG模型示意图；
图3b为WTEG模型示意图；
图3c为DAG模型映射至WTEG模型过程示意图；
图4为LEO-SCN前3时隙的WTEG模型示意图；
图5为采用地面云处理业务请求时延与星间协同处理业务请求时延对比仿真示意图；
图6为计算复杂度系数对时延的影响仿真示意图；
图7为不同数据缩放因子组合业务时延对比仿真示意图；
图8为采用不同算法的处理时延比较仿真示意图。
具体实施方式
以下结合附图和具体实施例对本发明作具体的介绍。
为了能够利用LEO-SCN机载计算资源实现星间协同计算构想，解决远距离的星地传输带来的高传输时延问题，本发明采用如图1所示的典型的空基云雾网络架构进行研究。业务请求被卸载到传输路径中的卫星上进行分布式计算，不再转发至地面云中心集中计算，在传输过程中完成业务处理。从计算和网络的角度该空基云雾网络架构可以被视为卫星雾层和地面云层两部分。
卫星雾层由具有一定计算能力的低轨通信卫星集群组成，能够为沙漠、海洋、高空等区域的用户提供高速的信息服务。为了保证信息服务质量，当业务请求到达发起卫星后，业务被调度到不同的卫星上协同计算，并将处理结果直接传给目标用户。地面云层由高性能服务器组成，通过地面站与卫星雾层相连，对卫星雾层进行监控和管理，包括网络资源调度、星间链路状态监控、接入用户信息管理等。此外，为了实现星间协同计算，地面云层需要根据卫星雾层网络资源状态、业务图模型等信息，在业务请求到达发起卫星前，预先规划业务调度方案，并将调度指令缓存到相应的通信卫星。
其中，随着业务模式的改变，爆发出灵活多样的用户需求，因此，新型信息业务的处理过程往往复杂多变，难以用统一的单业务表示。为了更好的处理新型信息业务，研究者们通常将其划分为多个可单独执行且不可分割的业务请求，并按照业务请求间的相互依赖关系建立有向无环图(directed acyclic graph,DAG)模型。加权时间扩展图(Weightedtime extension graph，WTEG)模型将LEO-SCN动态拓扑转化为高维稳态图，完整的表征了任意时刻网络的连通情况，能够分析网络的最短路由和最大容量且分别为时隙内和时隙间的边设计时延权重参数，用于星间协同计算方案的时延性能分析。
如图2所示，本发明实施例提供的基于通信卫星间协同计算的业务请求分配方法包括以下步骤：
S101，根据预设的基于业务请求的有向无环图DAG及基于低轨通信卫星网络的加权时间扩展图WTEG之间的映射规则，生成业务请求对应的多个映射结果；
S102，根据二值粒子群优化算法，从多个映射结果中获取最优的映射结果；
S103，根据映射结果，将业务请求分配至各个通信卫星。
可选地，该映射规则的创建过程包括：
采用表达式Φ＝(Ω,Γ)作为有向无环图DAG的数学模型，其中，Ω＝{ω1,...,ωl}为业务请求节点集合，构成数学模型Φ的业务请求集群；Γ为边集合，表示业务请求间的依赖关系；Φ为单起点、单终点的数学模型，即ω1为业务请求起点，ωl为业务请求终点，其余为中间业务请求节点；其中，业务请求ωi由元组{εi,ηi,Di}表征，其中，εi为计算复杂度系数，表征执行业务请求ωi耗费计算资源的程度；ηi为数据缩放因子参数，反映执行业务请求ωi过程中对输入数据量放大与缩小的过程；Di表示ωi的输入数据量，定义Θ↑(ωi)＝{ωj|(ωj,ωi)∈Γ}表示业务请求ωi的先驱节点集合，则Di可通过式
Figure BDA0002548295920000091
求得，其中，Dj和ηj属于表征服务功能ωj的元组；
采用表达式N＝(NT,ET,Graph)作为加权时间扩展图WTEG的数学模型，其中，NT为加权时间扩展图WTEG的节点集合，ET为加权时间扩展图WTEG的边集合，Graph为加权时间扩展图WTEG的时延权重组合矩阵，其中，
Figure BDA0002548295920000092
Figure BDA0002548295920000101
表示相邻时隙q与q+1间所有边的权重，
Figure BDA0002548295920000102
表示时隙q内所有边的权重；
制定Ω至NT的映射规则，包括：
根据表达式
Figure BDA0002548295920000103
将业务请求起点映射至业务请求发起通信卫星，将业务请求终点映射至结果接收通信卫星；将中间业务请求节点映射至任意通信卫星，由于业务请求之间不可分割，故每个业务请求都需要分配到单颗低轨通信卫星上计算，但受计算过程影响，存在业务请求被缓存到下一时隙的现象，此时业务请求经过虚拟缓存链路到达下一时隙的副本节点，因此，任意业务请求ωi∈Ω会被映射至NT中同一通信卫星的ρi+1个副本节点，ρi为自然数，ρi为由业务请求ωi计算引起的跨时隙数目；
制定Γ至P的映射规则，包括：
根据表达式Z((ωi,ωj))＝Path(B(ωi)e,B(ωj)s)，将有向无环图DAG任意边
Figure BDA0002548295920000104
映射为加权时间扩展图WTEG节点集合B(ωi)中最后一个节点B(ωi)e到加权时间扩展图WTEG节点集合B(ωj)第一个节点B(ωj)s之间的最短路由用Path(B(ωi)e,B(ωj)s)表示，路由中包含的边序列即为边(ωi,ωj)的映射结果，其中，鉴于低轨通信卫星网络的动态特性，加权时间扩展图WTEG的部分节点之间无连通路由，因此，(ωi,ωj)的映射需满足表达式
Figure BDA0002548295920000111
作为一个具体的实施例，图3a为DAG模型示意图，图3b为WTEG模型示意图，图3c为DAG模型映射至WTEG模型过程示意图，DAG被映射到WTEG中3个连续的时隙上。其中，节点ω1、ω2、ω3及ω5分别被映射到单个网络节点
Figure BDA0002548295920000112
及
Figure BDA0002548295920000113
而节点ω4受计算过程的影响，在相邻两个时隙的节点进行计算，故被映射到节点集合
Figure BDA0002548295920000114
同时，边映射结果分别为：
Figure BDA0002548295920000115
Figure BDA0002548295920000116
及
Figure BDA0002548295920000117
由此可得，DAG在WTEG上的映射结果与DAG是异构的，DAG节点和边的数目在映射过程中都可能改变。
可选地，根据二值粒子群优化算法，求解表达式T(x)的最小值，得到多个映射结果中最优的映射结果包括：
假设粒子群规模为U，同时I表示群体最大迭代次数，则第i次迭代过程中第u个粒子的位置和速度表达式分别如式(1)及式(2)所示：
Figure BDA0002548295920000118
Figure BDA0002548295920000119
式(1)中，
Figure BDA00025482959200001110
式(2)中，
Figure BDA00025482959200001111
当i<I时，粒子u首先通过跟踪其历史最优位置pubest和整个粒子群的历史最优位置gbest更新其速度，如式(3)所示，
Figure BDA0002548295920000121
其中，μ为惯性权重，γ1、γ2为加速因子，β1、β2为均匀分布于区间[0,1]的随机数；随后跟踪速度
Figure BDA0002548295920000122
方向按照式(4)更新其位置，
Figure BDA0002548295920000123
式(4)中，
Figure BDA0002548295920000124
表示数值向上取整运算，
Figure BDA0002548295920000125
表示两个矩阵的元素之间的置换操作,具体过程是找出
Figure BDA0002548295920000126
中数值最大的元素索引作为待置换位置，根据式的约束条件s.t.
Figure BDA0002548295920000127
计算
Figure BDA0002548295920000128
中与待置换位置进行置换操作的元素索引，完成置换操作，通过式
Figure BDA0002548295920000129
计算其适应度值；
将适应度值定义为当前粒子映射结果的业务请求处理总时延，将表达式
Figure BDA00025482959200001210
分为最后一项业务请求的计算时延
Figure BDA00025482959200001211
和前项业务请求的累积处理时延
Figure BDA00025482959200001212
运用递归原理，最终计算得到粒子的适应度值，即多个映射结果中最优的映射结果，其中，
Figure BDA00025482959200001213
表示ωi在加权时间扩展图WTEG中映射的最后一个节点，
Figure BDA0002548295920000131
表示ωl在加权时间扩展图WTEG中映射的第一个节点，T(X)为业务请求处理总时延。
其中，为了最小化业务处理时延，基于上述过程生成的映射规则，本发明以业务处理时延为优化目标构建理论优化问题，求解时延最优的星间协同计算方案。为了便于描述又不失一般性，本发明忽略业务请求由用户至发起卫星及业务处理结果由接收卫星至用户的过程带来的时延。
对于有向无环图DAG，其进行到节点ωi∈Γ的处理时延可由式(5)表示，
T(ωi)＝Tcomp(ωi)+Taccu(ωi) (5)
其中，Tcomp(ωi)表示ωi的计算时延，如式(6)所示，其中
Figure BDA0002548295920000132
表示节点B(ωi)s的计算能力；
Figure BDA0002548295920000133
Taccu(ωi)表示Φ进行到ωi时的累积时延，为ωi所有前向节点的处理时延加上由前向节点转发所需数据量至ωi的传输时延的最大值，故ωi的累积时延Taccu(ωi)如式(7)所示。
Figure BDA0002548295920000134
其中
Figure BDA0002548295920000135
为最短路由Z((ωj,ωi))传输单位数据量的时延，由式(6)式及(7)可知，式(5)也可以通过式(8)表示。
Figure BDA0002548295920000136
因此，Φ的业务处理时延如式(9)所示。
Figure BDA0002548295920000141
对于相同的DAG与WTEG，存在多个映射结果同时上述映射规则。而不同的映射结果会经过图N中不同的边和节点进行传输和计算，导致不同的业务处理时延，因此，最小化业务请求处理时延就是寻找时延最低的映射结果。对于任意的图Φ与图N，任一个映射结果都可以由一个l行n×d列的决策矩阵X唯一确定，决策矩阵X如式(10)所示。
Figure BDA0002548295920000142
其中，
Figure BDA0002548295920000143
表示节点ωi是否被映射到节点
Figure BDA0002548295920000144
如式(11)所示。当
Figure BDA0002548295920000145
时，节点ωi被映射至节点
Figure BDA0002548295920000146
反之，当
Figure BDA0002548295920000147
时，节点ωi不被映射至节点
Figure BDA0002548295920000148
Figure BDA0002548295920000149
此外，ωi被映射到WTEG中ρi+1个节点，所以ωi在WTEG中的映射还应满足式(12)所示条件。
Figure BDA00025482959200001410
由式(9)至式(12)可得，某个映射结果的业务处理时延可用式(13)表示，其中，
Figure BDA00025482959200001411
表示ωi在图N中映射的最后一个节点，
Figure BDA00025482959200001412
表示ωl在图N中映射的第一个节点。
Figure BDA00025482959200001413
为了验证本发明实施例提供的基于通信卫星间协同计算的业务请求分配方法的优越性，下面提供三组对比仿真实验，验证在LEO-SCN中业务图驱动的协同计算方案的时延性能。首先，对比星间协同计算与地面云计算的业务处理时延，随后，分析计算复杂度系数变化对业务处理时延的影响，最后，比较通信卫星间协同计算对不同数据缩放因子组合业务的承载能力。
仿真平台采用MATLAB仿真软件，实验室计算机CPU为i5-7400，内存为24GB。仿真设计了一个局部的12星LEO-SCN，低轨通信卫星均匀分布在高度大约500km的4个轨道，轨道倾角都为90°，且任意轨道都与其相邻轨道反向。考虑到太空中即将部署数万颗LEO卫星，彼时LEO-SCN拓扑切换将十分频繁，所以时隙时长设为100ms，并且LEO卫星的初始相位被设置在拓扑切换的临界点。不失一般性，假设业务在初始时刻到达接入卫星，此时LEO-SCN前3时隙的WTEG模型如图4所示。LEO-SCN系统参数设置如表1所示，Bsc为星地链路带宽，ρc为地面云服务器计算能力。基于表1给出的参数值，星间设计速率r*设定为1000Mbps。此外，仿真采用业务Φ的DAG如图3a所示，且对于
Figure BDA0002548295920000151
εi＝1900/8cycles/bit且ηi＝1，二值粒子群优化算法的基本参数为：种群规模U＝100，最大迭代次数I＝50，加速因子r1＝r2＝1，惯性权重μ＝1.5。
表1
Figure BDA0002548295920000152
Figure BDA0002548295920000161
下面分别测试星间协同计算与地面云计算的业务处理时延并进行比较：
如图5所示，结果表明，一方面，星间协同计算的时延性能始终优于地面云计算，这是由于地面云计算模式中，星地链路远距离数据传输会带来超高的传输时延，然而星间协同计算将业务调度到多颗卫星上执行，有效地避免了星地远距离传输；另一方面，两者的差距呈现出愈来愈大的趋势，这是因为在数据量较小时，地面云服务器的计算能力可以弥补传输时延造成的损失，但是随着数据量不断增加，传输时延的影响愈发明显。具体而言，当数据量较小，如1Mb时，星间协同计算与地面云计算的业务处理时延分别为18.38ms和44.75ms，此时两者的时延差仅为26.37ms。随着数据量由1Mb增加至5Mb，星间协同计算的业务处理时延平均增速为18.37ms/Mb，计算时延增速稍高于传输时延增速，两者分别为11.87ms/Mb和6.5ms/Mb；相反地，地面云计算的业务处理时延增速一直居高不下，稳定在44.75ms/Mb，其中传输时延增速高达40ms/Mb，占处理总时延的\89.4％，是处理时延增速居高不下的主要原因。当数据量增至5Mb时，两者的处理时延分别为91.87ms和223.75ms，两者的时延差已增至131.88ms。此时，相比地面云计算，星间协同计算将时延性能提高了58.9％。因此，星间协同计算方法更适应新型信息业务的低时延需求。
下面讨论计算复杂度系数对业务处理时延的影响：
以业务请求Φ为例，当其计算复杂度从0.2εi增加到1.8εi时，业务处理时延的变化趋势如图6所示。由图6可知，若数据量一定，随着计算复杂度升高，业务处理时延逐渐增加。而且，业务处理时延的增加在不同数据量约束下呈现两种趋势。当数据量为1Mb时，随着计算复杂度由0.2εi升高到1.8εi，业务总时延由6.4ms平缓增加至27.9ms。这是因为随着计算复杂度升高，业务的计算时延相应升高引起业务处理时延相应增加；然而，当数据量为5Mb时，随着计算复杂度升高，业务总时延由50.9ms跳跃式增加至152.3ms。特别是，业务处理时延在复杂度为1.4εi时由98.8ms跳跃式增加至125.5ms。这是因为当计算复杂度增加到为1.4εi时，由于计算时延增长，Φ不能在原本的时隙完成，数据需要通过虚拟缓存链路到达下一时隙，此缓存过程导致传输时延显著增加，最终，业务处理时延出现跳跃式增长。
下面比较星间协同计算对不同数据缩放因子组合业务的承载能力：
不同服务功能会基于其特性，按照一定的比率放大或缩小所经过的业务数据量，本节采用2、1、0.5三组缩放因子值设计仿真。如图3a所示的节点集合可以被划分为前排节点(ω1,ω2)和后排节点(ω3,ω4,ω5)，则按照不同的数据缩放因子组合，业务可以被分为前缩后放型和前放后缩型。仿真设计了4种业务Φ1,Φ2,Φ3,Φ4与Φ作对比，Φ1,Φ2,Φ3,Φ4的数据缩放因子组合分别设置为{0.5,1,1,1,2}、{2,1,1,1,0.5}、{1,1,0.5,2,1}及{1,1,2,0.5,1}，其中Φ3及Φ4为了讨论在并行业务分支上数据一缩一放的情况对业务处理时延的影响。仿真结果如图7所示。仿真结果表明，当数据量为5Mb时，Φ、Φ1、Φ2、Φ3与Φ4的业务处理时延分别为91.87ms、62.81ms、138.75ms、95.68ms和95.93ms。其中，Φ1的业务处理时延最低，这是因为其在业务起点将数据量压缩为输入数据量的一半，使得后续过程所需处理的数据量减半。同理，由于业务起点将数据量放大为输入数据量的两倍，导致后续过程所需处理的数据量翻倍，所以Φ2的业务处理时延最高。此外，Φ3与Φ4的业务处理时延相当且略高于Φ的总时延。由式(7)可知，ω5的累积时延为ω3与ω4引起的累积时延的最大值，故当输入数据量和计算复杂度相等时，ω5的累积时延由数据缩放因子最大的分支决定。综上所述，星间协同计算方案在处理前缩后放型业务时更能发挥其性能优势，而且数据缩小因子更宜分布在DAG主干节点。
下面分别测试几种经典算法(随机动态算法Pick-KX、加权轮转算法WRR、贪婪负载均衡算法Greedy-LB)的卫星间协同计算时延并与二值粒子群优化算法的时延进行比较，仿真结果如图8所示。结果表明，当数据量小于1Mb时，四种算法的时延差别不明显，随着数据量不断增大，算法间的时延差异逐渐显现。例如，当数据量为5Mb时，随机动态算法、加权轮转算法、贪婪负载均衡算法及二值粒子群优化算法算法的时延分别为181.50ms、155.37ms、141.25ms及91.87ms。这是因为随机动态算法算法在LEO-SCN中随机选取卫星进行映射，所以时延性能最差；加权轮转算法和贪婪负载均衡算法均以卫星的计算能力为依据选取卫星，所以其时延性能相差不大，且较Pick-KX算法时延性能有所提升；BPSO算法将链路传输速率和卫星计算能力联合考虑选取卫星，通过优化业务处理时延搜索最佳映射方案，故相比其他3种经典映射算法能更有效降低星间协作计算的业务处理时延。
本发明实施例提供的基于通信卫星间协同计算的业务请求分配方法，根据预设的基于业务请求的有向无环图DAG及基于低轨通信卫星网络的加权时间扩展图WTEG之间的映射规则，生成业务请求对应的多个映射结果，根据二值粒子群优化算法，从多个映射结果中获取最优的映射结果，根据映射结果，将业务请求分配至各个通信卫星，有效地降低了卫星间协作计算的业务请求处理时延。
在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述的部分，可以参见其他实施例的相关描述。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本申请的实施例可提供为方法、系统或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
以上仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。

Claims (1)

1.一种基于通信卫星间协同计算的业务请求分配方法，其特征在于，包括：
根据预设的基于业务请求的有向无环图DAG及基于低轨通信卫星网络的加权时间扩展图WTEG之间的映射规则，生成业务请求对应的多个映射结果；
根据二值粒子群优化算法，从所述多个映射结果中获取最优的映射结果；
假设粒子群规模为U，同时I表示群体最大迭代次数，则第i次迭代过程中第u个粒子的位置和速度表达式分别如式(1)及式(2)所示：
Figure FDA0003523733640000011
Figure FDA0003523733640000012
式(1)中，
Figure FDA0003523733640000013
式(2)中，
Figure FDA0003523733640000014
当i＜I时，粒子u首先通过跟踪其历史最优位置pubest和整个粒子群的历史最优位置gbest更新其速度，如式(3)所示，
Figure FDA0003523733640000015
其中，μ为惯性权重，γ1、γ2为加速因子，β1、β2为均匀分布于区间[0,1]的随机数；随后跟踪速度
Figure FDA0003523733640000016
方向按照式(4)更新其位置，
Figure FDA0003523733640000017
式(4)中，
Figure FDA0003523733640000018
表示数值向上取整运算，
Figure FDA0003523733640000019
表示两个矩阵的元素之间的置换操作,具体过程是找出
Figure FDA00035237336400000110
中数值最大的元素索引作为待置换位置，根据公式的约束条件
Figure FDA0003523733640000021
计算
Figure FDA0003523733640000022
中与待置换位置进行置换操作的元素索引，完成置换操作，通过式
Figure FDA0003523733640000023
计算其适应度值；
将适应度值定义为当前粒子映射结果的业务请求处理总时延，将表达式
Figure FDA0003523733640000024
分为最后一项业务请求的计算时延
Figure FDA0003523733640000025
和前项业务请求的累积处理时延
Figure FDA0003523733640000026
运用递归原理，最终计算得到粒子的适应度值，即所述多个映射结果中最优的映射结果，其中，
Figure FDA0003523733640000027
表示ωi在加权时间扩展图WTEG中映射的最后一个节点，
Figure FDA0003523733640000028
表示ωl在加权时间扩展图WTEG中映射的第一个节点，T(X)为业务请求处理总时延；
根据所述映射结果，将所述业务请求分配至各个通信卫星；
所述映射规则的创建过程包括：
采用表达式Φ＝(Ω,Γ)作为有向无环图DAG的数学模型，其中，Ω＝{ω1,...,ωl}为业务请求节点集合，构成Φ的业务请求集群；Γ为边集合，表示业务请求间的依赖关系；Φ为单起点、单终点的数学模型，即ω1为业务请求起点，ωl为业务请求终点，其余为中间业务请求节点；其中，业务请求ωi由元组{εi,ηi,Di}表征，其中，εi为计算复杂度系数，表征执行业务请求ωi耗费计算资源的程度；ηi为数据缩放因子参数，反映执行业务请求ωi过程中对输入数据量放大与缩小的过程；Di表示ωi的输入数据量，定义Θ↑(ωi)＝{ωj|(ωj,ωi)∈Γ}表示业务请求ωi的先驱节点集合，则Di可通过式
Figure FDA0003523733640000031
求得，其中，Dj和ηj属于表征服务功能ωj的元组；
采用表达式N＝(NT,ET,Graph)作为加权时间扩展图WTEG的数学模型，其中，NT为加权时间扩展图WTEG的节点集合，ET为加权时间扩展图WTEG的边集合，Graph为加权时间扩展图WTEG的时延权重组合矩阵，其中，
Figure FDA0003523733640000032
Figure FDA0003523733640000033
表示相邻时隙q与q+1间所有边的权重，
Figure FDA0003523733640000034
表示时隙q内所有边的权重；
制定Ω至NT的映射规则，包括：
根据表达式
Figure FDA0003523733640000041
将业务请求起点映射至业务请求发起通信卫星，将业务请求终点映射至结果接收通信卫星；将中间业务请求节点映射至任意通信卫星，由于业务请求之间不可分割，故每个业务请求都需要分配到单颗低轨通信卫星上计算，但受计算过程影响，存在业务请求被缓存到下一时隙的现象，此时业务请求经过虚拟缓存链路到达下一时隙的副本节点，因此，任意业务请求ωi∈Ω会被映射至NT中同一通信卫星的ρi+1个副本节点，ρi为自然数，ρi为由业务请求ωi计算引起的跨时隙数目；
制定Γ至P的映射规则，包括：
根据表达式Z((ωi,ωj))＝Path(B(ωi)e,B(ωj)s)，将有向无环图DAG任意边
Figure FDA0003523733640000042
映射为加权时间扩展图WTEG节点集合B(ωi)中最后一个节点B(ωi)e到加权时间扩展图WTEG节点集合B(ωj)第一个节点B(ωj)s之间的最短路由用Path(B(ωi)e,B(ωj)s)表示，路由中包含的边序列即为边(ωi,ωj)的映射结果，其中，鉴于低轨通信卫星网络的动态特性，加权时间扩展图WTEG的部分节点之间无连通路由，因此，(ωi,ωj)的映射需满足表达式
Figure FDA0003523733640000043
CN202010567261.7A 2020-06-19 2020-06-19 一种基于通信卫星间协同计算的业务请求分配方法 Active CN111884703B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010567261.7A CN111884703B (zh) 2020-06-19 2020-06-19 一种基于通信卫星间协同计算的业务请求分配方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010567261.7A CN111884703B (zh) 2020-06-19 2020-06-19 一种基于通信卫星间协同计算的业务请求分配方法

Publications (2)

Publication Number Publication Date
CN111884703A CN111884703A (zh) 2020-11-03
CN111884703B true CN111884703B (zh) 2022-08-26

Family

ID=73157850

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010567261.7A Active CN111884703B (zh) 2020-06-19 2020-06-19 一种基于通信卫星间协同计算的业务请求分配方法

Country Status (1)

Country Link
CN (1) CN111884703B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112803984B (zh) * 2020-12-29 2023-05-19 国网甘肃省电力公司信息通信公司 一种空天地一体化通信网络中卫星间切换的方法
CN113222468B (zh) * 2021-06-02 2022-04-08 中国电子科技集团公司第五十四研究所 一种基于深度强化学习的成像卫星资源调度方法
CN113986532B (zh) * 2021-10-15 2024-05-03 武汉大学 一种低轨卫星物联网分布式任务协作处理方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6366761B1 (en) * 1998-10-06 2002-04-02 Teledesic Llc Priority-based bandwidth allocation and bandwidth-on-demand in a low-earth-orbit satellite data communication network
CN107343025B (zh) * 2017-06-07 2021-01-26 西安电子科技大学 分布式卫星云雾网络架构及能耗约束下的时延优化方法
US20190074894A1 (en) * 2017-09-05 2019-03-07 Hughes Network Systems, Llc Convergent architectures for multi-orbit satellite communications
CN108804220A (zh) * 2018-01-31 2018-11-13 中国地质大学（武汉） 一种基于并行计算的卫星任务规划算法研究的方法
CN110493131B (zh) * 2019-09-24 2021-06-01 大连大学 一种sdn架构下的空间信息网络路由策略的设计方法
CN110830292B (zh) * 2019-11-01 2021-06-25 西安电子科技大学 面向医疗大数据的云雾混合路径确定方法

Also Published As

Publication number Publication date
CN111884703A (zh) 2020-11-03

Similar Documents

Publication Publication Date Title
CN111884703B (zh) 一种基于通信卫星间协同计算的业务请求分配方法
Liu et al. Deep reinforcement learning based latency minimization for mobile edge computing with virtualization in maritime UAV communication network
CN114362810B (zh) 一种基于迁移深度强化学习的低轨卫星跳波束优化方法
Cui et al. Latency optimization for hybrid GEO–LEO satellite-assisted IoT networks
CN112799784A (zh) 一种基于分散计算的低轨卫星网络最优任务分配方法
CN114665952A (zh) 一种基于星地融合架构下低轨卫星网络跳波束优化方法
Pfandzelter et al. Edge (of the earth) replication: Optimizing content delivery in large leo satellite communication networks
Zheng et al. LEO satellite channel allocation scheme based on reinforcement learning
CN115514769B (zh) 卫星弹性互联网资源调度方法、系统、计算机设备及介质
CN115276754A (zh) 一种基于栅格时延预测的卫星传输优化方法
CN116366134A (zh) 一种基于要素时空扩展图的天基物联网多星协同策略构建方法
CN115276756A (zh) 一种保障服务质量的低轨卫星星座优化设计方法
CN118250750B (zh) 基于深度强化学习的卫星边缘计算任务卸载及资源分配方法
CN113472420A (zh) 一种基于区域用户兴趣感知的卫星网络缓存放置方法
CN115866788A (zh) 面向mec任务主动迁移的天地融合网络3c资源调度方法
Lai et al. Cooperatively constructing cost-effective content distribution networks upon emerging low earth orbit satellites and clouds
Mao et al. Digital twin satellite networks towards 6G: Motivations, challenges, and future perspectives
He et al. Distributed satellite cluster laser networking algorithm with double-layer markov drl architecture
Zhou et al. Intelligent Multi-Objective Routing for Future Ultra-Dense LEO Satellite Networks
Wang et al. Optimal download of dynamically generated data by using ISL offloading in LEO networks
CN115765826B (zh) 一种面向按需服务的无人机网络拓扑重构方法
Wang et al. Task assignment strategy in LEO-muti-access edge computing based on matching game
CN116113050A (zh) 动态波束调度方法及装置
CN116633422A (zh) 面向IoT任务卸载的低轨卫星网络多维资源调度方法
Bao et al. Toward Intelligent Cross-Domain Resource Coordinate Scheduling for Satellite Networks

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