CN113242568A - 一种不确定网络环境中的任务卸载和资源分配方法 - Google Patents

Abstract

本发明涉及无线通信技术领域，尤其涉及一种不确定网络环境中的任务卸载和资源分配方法，包括将任务卸载过程建模为两阶段卸载模型，并将该模型优化为基于两阶段随机规划的任务卸载和资源分配问题，利用随机模拟方法将其化为样本均值近似问题，将该问题解耦为本地计算资源分配子问题、传输功率和边缘计算资源联合分配子问题以及卸载决策子问题；采用标准拉格朗日乘子法、遗传算法以及分析本地计算和边缘计算的时延估计和能耗预算求解三个子问题；用户根据求解三个子问题获得最优分配策略进行任务卸载；本发明能够在时延不确定的网络中满足任务计算时延的要求，同时保证系统能耗最小化。

Description

一种不确定网络环境中的任务卸载和资源分配方法

技术领域

本发明涉及无线通信技术领域，尤其涉及一种不确定网络环境中的任务卸载和资源分配方法。

背景技术

随着物联网技术的快速发展和5G/6G新型应用的普及，时延敏感性应用已经被广泛设想，如虚拟现实、无人驾驶和人脸识别等正以前所未有的速度发展。移动边缘计算(Mobile Edge Computing，MEC)通过将任务卸载网络边缘节点，如基站和无线接入点，为时延敏感型应用提供了固有的低时延优势。

移动边缘计算中任务卸载和资源分配策略往往是影响用户卸载时延和能耗的关键因素，为提高用户的服务质量(Quality of Experience，QoE)，同时最小化计算卸载能耗，需要根据任务的计算能耗预算、计算时延，以及MEC服务器的可用资源等性能约束条件来制定任务卸载和资源分配策略。国内外研究人员针对此问题进行了大量深入研究，一些主要的成果有：

(1)移动边缘计算中的任务分配和计算频率任务算法(参考文献：Dinh T Q,TangJ,La Q D,et al.Offloading in Mobile Edge Computing:Task Allocation andComputational Frequency Scaling[J].IEEE Transactions on Communications,2017,65(8):3571-3584.)：该算法考虑了单个用户将任务卸载到多个MEC服务器的场景，通过优化任务的卸载决策和用户的中央处理单元(Central Process Unit，CPU)频率值来实现任务计算时延和用户能耗的最小化。

(2)联合任务卸载和无线资源分配算法(参考文献：Chen H,D Zhao,Chen Q,etal.Joint Computation Offloading and Radio Resource Allocations in Small-CellWireless Cellular Networks[J].IEEE Transactions on Green Communications andNetworking,2020,4(3):745-758.)：该算法考虑在多用户多MEC服务器场景下，对任务卸载决策、传输功率和信道分配进行了联合优化，目的是在满足任务计算时延约束条件下，最小化所有用户的能耗。

任务计算时延是保证用户QoE的重要指标之一，上述研究工作均考虑了任务的计算时延，然而都忽略了MEC服务器任务队列的排队等待时延。在实际的MEC网络环境中，由于MEC服务器的计算、存储等资源非常有限，通常无法快速响应海量突发的计算请求，因此，任务在MEC服务器端的排队等待时间是不可忽略的。由于任务到达的随机性，以及MEC服务器任务队列的时变性，欲得到任务排队等待时间的准确预测值是不现实的，这种不确定性因素对传统任务卸载和资源分配造成了严峻挑战。此外，大部分研究工作都只针对用户端能耗，而计算资源有限的MEC服务器端将面临着大量计算密集型和时延敏感型用户的访问，随之而来的是计算资源匮乏，以及高能耗等问题。因此，如何在任务计算时延不确定的MEC网络环境中制定高效、绿色的任务卸载和资源分配策略具有重要研究价值。

发明内容

为解上述问题，本发明考虑任务在MEC服务器的随机排队等待时间导致的计算时延不确定问题，提出了一种不确定网络环境中的任务卸载和资源分配方法，具体包括以下步骤：

利用随机规划理论将任务在MEC服务器的排队等待时间建模为一组随机参数，并将任务卸载过程建模为两阶段卸载模型；

基于两阶段卸载模型，在时延约束条件下，以最小化系统总能耗为目标，将优化问题建模为基于两阶段随机规划的任务卸载和资源分配问题；

利用随机模拟方法将基于两阶段随机规划的任务卸载和资源分配问题转化为样本均值近似问题；

将样本均值近似问题解耦为本地计算资源分配子问题、传输功率和边缘计算资源联合分配子问题，以及卸载决策子问题；

采用标准拉格朗日乘子法获得所述本地计算资源分配子问题的最优分配策略；

采用遗传算法获得所述传输功率和边缘计算资源联合分配子问题的最优分配策略；

通过分析本地计算和边缘计算的时延估计和能耗预算获得所述卸载决策子问题的最优策略；

用户基于所述本地计算资源分配子问题的最优分配策略、传输功率和边缘计算资源联合分配子问题的最优分配策略，以及卸载决策子问题的最优策略进行任务卸载。

进一步的，在两阶段卸载模型中，将卸载决策过程分为两个阶段，包括：

第一阶段的决策变量为任务的传输功率，即用户在没有观察到不确定性MEC服务器排队等待时间的情况下，考虑未来所有可能排队时间的影响，做出第一阶段的传输功率分配策略p_i；

第二阶段的决策变量为MEC服务器的CPU频率资源，即当任务上传至MEC服务器，排队等待时间实现已知，在获得排队等待时间和第一阶段的传输功率分配策略p_i的条件下，MEC服务器将根据约束条件采取追索动作，对第一阶段策略进行补偿。

进一步的，若第i个用户的任务卸载决策为1，即第i个用户将任务卸载到服务器计算，则通过遗传算法计算最优的卸载策略，具体包括：

将每一个可行的传输功率进行浮点向量编码，每一个浮点向量表示一个染色体，浮点向量维数与解向量维数一致；

从用户i传输功率的可行域中随机产生一个点，检验其是否满足用户传输时延小于最大传输时延要求，即

如果满足则作为一个染色体，否则，重新产生一个随机点，直到满足约束条件，重复以上过程M次，产生M个染色体作为初始种群，记为p_i,1,p_i,2,...,p_i,M；

对于每个染色体p_i,m，求解其适应度函数，并根据适应度求其评价函数，且对于每个染色体评价函数值越大，被选择产生后代的概率越大，其中适应度为样本均值近似问题函数；

对于染色体p_i,m，计算前m个染色的累积概率

从区间(0,q_i,M)中生成一个随机数r，若满足q_i,m-1＜r＜q_i,m，q₀＝0，则选择染色体p_i,m作为父本繁殖下一代，依此规则选择M个染色体；

定义P_c为交叉概率，则种群中有期望值为P_c×M个染色体进行交叉操作，即假设p_i,1和p_i,2为要交叉的两个染色体，则从区间(0,1)中生成一个随机数c，直到p'_i,1＝cp_i,1+(1-c)p_i,2和p'_i,2＝(1-c)p_i,1+cp_i,2满足用户传输时延小于最大传输时延要求，即

则p'_i,1和p'_i,2将代替原染色体p_i,1和p_i,2成为两个新的染色体；

定义P_u为变异概率，则种群中有期望值为P_u×M个染色体进行变异操作，即假设p_i,1为要变异的染色体，随机选择一个变异方向d，从

中生成一个随机数α，直到p'_i,1＝p_i,1+αd满足约束条件，则p'_i,1将取代原染色体p_i,1成为新的染色体；

经过选择、交叉和变异操作，可以得到一个新的种群，并准备进行下一代进化，如果上述步骤达到了给定的循环次数，则遗传算法终止；算法终止后，从最后一代选择一个适应度最高的染色体，即得到优化问题的全局最优解。

本发明针对移动边缘计算中任务在边缘服务器的随机排队等待时间导致的计算时延不确定问题，提出了一种不确定网络环境中的任务卸载和资源分配算法。首先，在计算时延约束条件下，以最小化系统总能耗为目标，将优化问题建模为两阶段随机规划问题。其次，为了降低两阶段随机规划问题的计算复杂度，利用随机模拟方法将原问题转化为基于样本均值近似的混合整数非线性规划问题，并将该问题解耦为本地计算资源分配、传输功率和边缘计算资源联合分配，以及卸载决策3个子问题。接着，采用拉格朗日乘子法获得本地计算资源最优分配策略，同时采用遗传算法获得传输功率和边缘计算资源最优联合分配策略。最后，通过分析本地和边缘计算的时延估计和能耗预算获得最优任务卸载决策。通过仿真实验，与传统算法相比，所提算法能够在时延不确定的网络中满足任务计算时延的要求，同时保证系统能耗最小化。

附图说明

图1为本发明中一种不确定网络环境中的任务卸载和资源分配算法流程图；

图2为本发明两阶段卸载模型图；

图3为不同算法下系统总能耗随着任务量的变化过程图；

图4为不同算法下系统总能耗随着时延约束的变化过程图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明提出一种不确定网络环境中的任务卸载和资源分配方法，如图1，包括以下步骤：

利用随机规划理论将任务在MEC服务器的排队等待时间建模为一组随机参数，并将任务卸载过程建模为两阶段卸载模型；

基于两阶段卸载模型，在时延约束条件下，以最小化系统总能耗为目标，将优化问题建模为基于两阶段随机规划的任务卸载和资源分配问题；

利用随机模拟方法将基于两阶段随机规划的任务卸载和资源分配问题转化为样本均值近似问题；

将样本均值近似问题解耦为本地计算资源分配子问题、传输功率和边缘计算资源联合分配子问题，以及卸载决策子问题；

采用标准拉格朗日乘子法获得所述本地计算资源分配子问题的最优分配策略；

采用遗传算法获得所述传输功率和边缘计算资源联合分配子问题的最优分配策略；

通过分析本地计算和边缘计算的时延估计和能耗预算获得所述卸载决策子问题的最优策略；

用户基于所述本地计算资源分配子问题的最优分配策略、传输功率和边缘计算资源联合分配子问题的最优分配策略，以及卸载决策子问题的最优策略进行任务卸载。

在本实施例中，考虑一个MEC系统模型，该模型包含一个配置了服务器的基站，N个请求任务处理的用户，集合表示为

令UEi表示第i个用户，其中

UEs可以通过无线网络与MEC服务器进行通信。定义用户需要处理的任务为一个三元组A_i＝{D_i,L_i,τ_i}，其中，D_i表示需要计算的任务量大小；L_i表示单位bit任务所需要的CPU周期数，单位为cycles/bit；τ_i表示计算完成该任务的最大时延要求；假设任务不可分割，定义π_i∈{0,1}表示UEi的卸载决策，若选择本地计算任务，则π_i＝0，若选择将任务卸载到MEC服务器处理，则π_i＝1。

每个用户的处理器均支持动态电压频率调整(Dynamic Voltage FrequencyScaling，DVFS)技术，DVFS技术可以动态调整CPU频率值，从而达到节省功耗的目的本地计算的时延为：

本地计算时UEi产生的计算能耗为：

其中，f_i ^l表示本地计算时UEi的CPU频率值，κ_i表示与UEi芯片架构相关的有效能量系数。

用户选择将任务卸载到MEC服务器处理主要经历四个过程：任务上传、MEC服务器排队等待、MEC服务器计算，以及计算结果返回。由于返回结果时的任务量较小，为便于分析，本发明忽略计算结果返回的通信时延。UEi将任务卸载到MEC服务器可获得的传输速率为：

其中，B₀表示通信带宽，h_i表示UEi和MEC服务器之间的信道增益，N₀为高斯白噪声功率谱，p_i为UEi卸载任务时的传输功率。因此，当UEi选择将任务卸载到MEC服务器处理时，任务上传的通信时延为：

当任务上传至MEC服务器后，由于MEC服务器计算资源的有限性，通常无法快速响应突发的计算请求，因此，MEC服务器排队等待时间是不可忽略的，定义T_i ^wait表示UEi在MEC服务器的排队等待时间；任务经过排队等待，MEC服务器将为其提供计算服务，定义f_i ^e表示MEC服务器处理任务时的CPU频率值，则UEi在MEC服务器的计算时延为：

进一步，边缘计算时MEC服务器为UEi提供计算服务所产生的计算能耗为：

其中，γ表示与MEC服务器芯片架构相关的有效能量系数。

实施例1

本实施例提供基于两阶段随机规划的任务卸载及资源分配问题建模的具体实施。

由于任务到达的随机性，以及MEC服务器任务队列的时变性，欲得到任务排队等待时间的准确预测值是不现实的。本实施例首先对排队等待时间进行不确定性分析，然后在满足用户计算时延要求的条件下，以最小化系统总能耗为目标，将优化问题建模为一种基于两阶段随机规划的任务卸载和资源分配问题，具体包括：

(一)、对排队等待时间进行不确定性分析

为了处理MEC服务器排队等待时间的不确定性问题，本发明利用随机规划理论将不确定的排队等待时间建模为一组由概率分布描述的随机参数。不失一般性，假设MEC服务器排队等待时间服从指数分布，定义Ω_i表示UEi卸载任务到MEC服务器时，所有可能的排队等待时间的集合，称为场景，令T_i ^wait∈Ω_i为该场景中的一个实现。进一步，考虑所有的UEs，定义组合场景Ω表示所有UEs在MEC服务器的排队等待时间的集合，可表示为笛卡尔乘积

令

表示组合场景中的一个组合实现。

(二)、构建基于两阶段随机规划的任务卸载和资源分配问题

本发明旨在满足用户处理时延要求的条件下，通过优化本地和MEC服务器的CPU频率资源、任务传输功率，以及任务卸载决策来最小化系统的总能耗。如图2所示，用户首先在没有观察到不确定性MEC服务器排队等待时间的情况下，考虑未来所有可能排队等待时间的影响，并通过衡量本地和边缘计算时的能耗预算和时延估计得到卸载决策π_i，当用户选择将任务卸载到MEC服务器处理时，本发明基于两阶段随机规划将卸载决策过程分为两个阶段：

1)第一阶段的决策变量为任务的传输功率。同样地，用户在没有观察到不确定性MEC服务器排队等待时间的情况下，考虑未来所有可能排队时间的影响，做出第一阶段的传输功率分配策略p_i；

2)第二阶段的决策变量为MEC服务器的CPU频率资源。当任务上传至MEC服务器，排队等待时间实现已知，在获得排队等待时间的实现T_i ^wait和第一阶段的策略p_i的条件下，MEC服务器将根据约束条件(比如卸载时延)采取追索动作f_i ^e(p_i,T_i ^wait)，以弥补第一阶段策略的不准确预测。

结合以上分析，可得到基于两阶段随机规划的任务卸载和资源分配问题的期望值模型，表示为：

P1:

目标函数：

f_i ^l，min≤f_i ^l≤f_i ^l，max (8-b)

f^e,min≤f_i ^l≤f^e,max (8-c)

T_i ^loca(f_i ^l)≤τ_i (8-d)

π_i∈{0,1} (8-f)

其中，π＝{π₁,π₂,...,π_N}表示任务卸载决策集合，

表示本地CPU频率资源分配策略集合；

表示所有组合场景下MEC服务器CPU频率资源分配策略集合，其中

p＝{p₁,p₂,...p_N}表示UEs任务传输功率分配策略集合；

表示组合场景的期望；

和

分别表示传输功率的最小值和最大值；f_i ^l,min和f_i ^l,max分别表示本地CPU频率的最小值和最大值；f^e,min和f^e,max分别表示MEC服务器CPU频率的最小值和最大值；T_i ^wait(ω)表示组合实现为ω时UEi的排队等待时间。约束条件式(8-d)和式(8-e)分别表示在本地和MEC服务器执行任务的时间须在其时延要求范围内。

实施例2

本实施例为了降低实施例1提出的两阶段随机规划问题的计算复杂度，将问题P1中的期望值模型转化为基于样本均值近似的MINLP问题。接着，将MINLP问题解耦为本地计算资源分配、传输功率和边缘计算资源联合分配，以及卸载决策3个子问题进行求解，具体包括：

(一)、将问题P1中的期望值模型转化为基于样本均值近似的MINLP问题

本发明考虑了MEC服务器排队等待时间的不确定环境下的任务卸载和资源分配优化问题，并将优化问题建模为一个两阶段随机规划问题。然而，求解两阶段随机规划问题通常面临着“维数灾难”的挑战，这将导致较高的计算复杂度。例如，当用户的场景空间Ω_i的数量为1000时，组合场景Ω的数量将达到1000^N，求解场景数量如此庞大的两阶段随机规划问题是困难且不现实的。

为了解决两阶段随机规划中的“维数灾难”问题，本发明首先采用随机模拟方法从场景空间Ω_i中抽取K个独立同分布的样本，并组成样本场景

其中，

表示该样本中的第k个实现。进一步考虑所有用户，可得到基于样本空间的样本组合场景Ω'＝{ω₁,ω₂,...,ω_S}，其中，S＝K^N表示该样本组合场景的数量，ω_s表示该样本组合场景中的第s个组合实现。接着，基于样本均值近似方法将问题P1中的期望值模型近似为样本均值模型，表示为：

P2：

目标条件：式(8-a)-式(8-f)

在优化问题P2中，π_i为{0,1}二元决策变量，边缘计算第一阶段目标函数是关于p_i的非线性函数，因此，优化问题P2为MINLP问题，且为NP难问题。由于本地计算资源分配变量f^l、传输功率分配变量p和边缘计算资源分配变量f^e，以及卸载决策变量π之间完全解耦，因此，为解决MINLP问题，本发明将优化问题P2解耦为本地计算资源分配子问题、传输功率和边缘计算资源联合分配子问题，以及卸载决策子问题进行求解。为便于分析，本发明接下来只针对一个用户求解其最优策略，由于每个用户之间相互独立，最优策略容易扩展到多用户场景。

(二)、本地计算资源分配

假设卸载决策变量π_i＝0，即UEi选择将任务在本地处理，则优化问题P2可写为本地计算资源分配子问题，表示为：

P2-1：

约束条件：式(8-b)-式(8-d)

在优化问题P2-1中，目标函数是关于f_i ^l的仿射函数，约束条件式(8-c)是关于f_i ^l的凸约束，因此，优化问题P2-1为凸规划问题。采用拉格朗日乘子法求解优化问题P2-1，定义拉格朗日函数为：

其中，μ_i≥0、λ_i≥0和ν_i≥0均表示拉格朗日乘子。进一步可得问题式(X)的KKT(Karush-Kuhn-Tucker)条件如下：

通过求解KKT条件，可得到任务在本地计算时最优CPU频率分配策略的闭合表达式为：

(三)、传输功率和边缘计算资源联合分配

假设卸载决策变量π_i＝1，即UEi选择将任务卸载到MEC服务器计算，则优化问题P2可写为传输功率和边缘计算资源分配子问题，表示为：

P2-2：

约束条件：式(8-a)、式(8-c)和式(8-e)

在优化问题(12)中，由于边缘计算第一阶段目标函数是关于p_i的非凸函数，为此，本发明接下来采用遗传算法求取P2-2问题的全局最优解，遗传算法是受到根据生物进化论的启发而得出的一种全局优化算法，并且特别适用于高维决策问题。

实施例3

本实施例提出采用遗传算法求取P2-2问题的全局最优解的具体实施，具体包括以下步骤：

1)编码

在两阶段卸载模型中，主要目标是确定第一阶段任务传输功率，因此，本发明将每一个可行的传输功率进行浮点向量编码，每一个浮点向量表示一个染色体，浮点向量维数与解向量维数一致。

2)初始种群

定义M表示种群大小，初始化过程随机产生M个染色体。从用户传输功率的可行域中随机产生一个点，并检验其是否满足约束条件，如果满足则作为一个染色体，否则，重新产生一个随机点，直到满足约束条件。重复以上过程M次，产生M个染色体作为初始种群，记为p_i,1,p_i,2,...,p_i,M。

3)评价函数

评价函数用来对种群中的每个染色体设定一个概率，以使该染色体被选择的可能性与种群中其他染色体的适应度成比例，适应度强的染色体被选择产生后代的几率较大。本发明利用优化问题P2-2的最优值作为染色体的适应度，利用适应度的大小决定一个序，并按照此序确定染色体被选择的概率。最优值越小，适应度越强，即作为父本繁殖下一代的概率越大。具体操作为对于每个染色体p_i,m，m＝1,2,...,M，首先求解如下最小化问题的目标函数最优值作为适应度值，即：

S.t.：式(8-c)和式(8-e)

在优化问题式(15)中，边缘计算第一阶段的目标值为常数，边缘计算第二阶段目标函数是关于

的仿射函数，约束条件式(8-d)是关于

的凸约束，因此，优化问题式(15)为凸规划问题。同样地，采用拉格朗日乘子法和KKT条件可以求得每个场景下MEC服务器最优CPU频率分配策略的闭合表达式为：

将所有场景下最优CPU频率值带入式(15)，可求得每个染色体相应的适应度值，并按照适应度值从小到大对染色体进行排序，据此定义评价函数如下

eval(p_i,m)＝a(1-a)^m-1,m＝1,2,...,M (17)

其中，a∈(0,1)。对于染色体p_i,m，评价函数值越大，被选择产生后代的概率越大。

4)选择

对于每一个染色体p_i,m，计算前m个染色体的累积概率，表示为：

从区间(0,q_i,M)中生成一个随机数r，若满足q_i,m-1＜r＜q_i,m，q_i,0＝0，则选择染色体p_i,m作为父本繁殖下一代，依此规则选择M个染色体。

5)交叉

定义P_c为交叉概率，则种群中有期望值为P_c×M个染色体进行交叉操作。假设p_i,1和p_i,2为要交叉的两个染色体，首先从区间(0,1)中生成一个随机数c，直到p'_i,1＝cp_i,1+(1-c)p_i,2和p'_i,2＝(1-c)p_i,1+cp_i,2满足约束条件，则p'_i,1和p'_i,2将代替原染色体p_i,1和p_i,2成为两个新的染色体。

6)变异

定义P_u为变异概率，则种群中有期望值为P_u×M个染色体进行变异操作。假设p_i,1为要变异的染色体，随机选择一个变异方向d，然后从

中生成一个随机数α，直到p'_i,1＝p_i,1+αd满足约束条件，则p'_i,1将取代原染色体p_i,1成为新的染色体。

7)终止条件

经过选择、交叉和变异操作，可以得到一个新的种群，并准备进行下一代进化。如果上述步骤达到了给定的循环次数，则遗传算法终止。算法终止后，从最后一代选择一个适应度最高的染色体，即得到优化问题的全局最优解，该解即为用户i最优发射功率。

实施例4

本实施例给出卸载决策的求解方法。

用户通过对本地计算和边缘计算的能耗预算和时延估计得到卸载决策，在获得本地计算资源、传输功率和边缘计算资源最优分配策略的情况下，优化问题P2可写为任务卸载决策子问题：

P2-3：

约束条件：式(8-f)

在优化问题P2-3中，唯一的变量是卸载决策π_i∈{0,1}，因此，优化问题P2-3为0-1规划问题。本发明采用枚举法分别将0和1带入问题P2-3，再比较目标函数值的大小，较小的目标函数值对应的π_i即为最优卸载决策，闭合表达式如下：

实施例5

本实施例为了验证本发明的有效性，进行了相关仿真实验。

本发明利用MATLAB工具对本发明提出的一种不确定网络环境下任务卸载和资源分配算法(SS_2SSP)进行仿真验证，为了验证所提算法的性能，对比算法包括：

(1)基于最大排队等待时间的静态卸载策略(Maximum Waiting Time basedStatic Offloading，MWT_SO)：用户基于任务在MEC服务器最大排队时间的静态网络中，对传输功率和边缘计算资源进行优化；

(2)贪婪式卸载策略(Greedy)：用户基于任务在MEC服务器最大排队时间的静态网络中，以及采用MEC服务器最大CPU频率的贪婪式卸载策略对传输功率进优化；

(3)本地计算策略(OnlyLocal)：用户将任务只在本地处理。

本发明考虑单个MEC服务器场景，MEC服务器排队等待时间服从指数分布，平均排队等待时间为0.1～4s；若无特别说明，排队等待时间样本数量K均为100个；计算密度L_i为700cycles/bit；本地CPU频率范围为100～2500MHz；MEC服务器CPU频率范围为500～5000MHz；本地有效能量系数k_i为10^-7；MEC服务器有效能量系数γ为10^-8；信道增益h_i为10^-2；传输功率范围为5～33dBm；传输带宽B₀为30MHz；白噪声功率N₀为-50dBm；种群大小M为140；迭代次数为200；交叉概率P_c为0.6；变异概率P_u为0.01；独立重复试验次数为10次。

图3所示为不同算法下系统总能耗随着任务量的变化过程图，其中，每个任务计算时延要求为1s/Mbit，MEC服务器最大排队等待时间为5s，平均排队等待时间为2s。从图3中可以看出，随着任务量的增加，所有算法的系统总能耗随之增加，而本发明所提SS_2SSP算法具有最低的系统总能耗，这是由于SS_2SSP算法在第一阶段制定传输功率策略时考虑了未来所有可能排队等待时间的影响，待任务上传至MEC服务器后，获得了排队等待时间的实现，并结合第一阶段传输功率分配策略，MEC服务器根据用户的计算时延约束条件采取第二阶段的追索策略，即CPU频率分配策略，这将弥补第一阶段策略的不准确预测。而MWT_SO算法是基于MEC服务器最大排队等待时间且没有追索动作的静态卸载策略，当MEC服务器平均排队等待时间较小时，MWT_SO算法将会以更高的能耗去满足时延需求。Greedy算法基于最大的排队等待时间，且贪婪式地使用MEC服务器最大CPU频率资源，相比MWT_SO算法和SS_2SSP算法具有更高的系统能耗。OnlyLocal算法由于将任务全在本地处理，较其他算法具有最高的系统能耗。

图4所示为不同算法下系统总能耗随着时延约束的变化过程图，其中，用户数1个，任务量大小为20Mbit，MEC服务器最大排队等待时间为2.5s，平均排队等待时间为1s。从图4可以看出，随着时延约束的增加，系统总能耗随之减小，本发明所提SS_2SSP算法具有最低的系统总能耗，而MWT_SO算法基于最大排队等待时间具有较高的系统总能耗，与MWT_SO算法对比，Greedy算法没有对MEC服务器CPU频率进行优化，因此，具有更高的系统能耗，OnlyLocal算法将任务只在本地处理，由于本地计算性能和计算资源的限制，与其他算法相比，具有最高的系统总能耗。

尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。

Claims (9)

1.一种不确定网络环境中的任务卸载和资源分配方法，其特征在于，包括以下步骤：

利用随机规划理论将任务在MEC服务器的排队等待时间建模为一组随机参数，并将任务卸载过程建模为两阶段卸载模型；

基于两阶段卸载模型，在时延约束条件下，以最小化系统总能耗为目标，将优化问题建模为基于两阶段随机规划的任务卸载和资源分配问题；

利用随机模拟方法将基于两阶段随机规划的任务卸载和资源分配问题转化为样本均值近似问题；

将样本均值近似问题解耦为本地计算资源分配子问题、传输功率和边缘计算资源联合分配子问题，以及卸载决策子问题；

采用标准拉格朗日乘子法获得所述本地计算资源分配子问题的最优分配策略；

采用遗传算法获得所述传输功率和边缘计算资源联合分配子问题的最优分配策略；

通过分析本地计算和边缘计算的时延估计和能耗预算获得所述卸载决策子问题的最优策略；

用户基于所述本地计算资源分配子问题的最优分配策略、传输功率和边缘计算资源联合分配子问题的最优分配策略，以及卸载决策子问题的最优策略进行任务卸载。

2.根据权利要求1所述的一种不确定网络环境中的任务卸载和资源分配方法，其特征在于，在两阶段卸载模型中，将卸载决策过程分为两个阶段，包括：

第一阶段的决策变量为任务的传输功率，即用户在没有观察到不确定性MEC服务器排队等待时间的情况下，考虑未来所有可能排队时间的影响，做出第一阶段的传输功率分配策略p_i；

第二阶段的决策变量为MEC服务器的CPU频率资源，即当任务上传至MEC服务器，排队等待时间实现已知，在获得排队等待时间和第一阶段的传输功率分配策略p_i的条件下，MEC服务器将根据约束条件采取追索动作，对第一阶段策略进行补偿。

3.根据权利要求1所述的一种不确定网络环境中的任务卸载和资源分配方法，其特征在于，基于两阶段随机规划的任务卸载和资源分配问题表示为：

约束条件：

f_i ^l，min≤f_i ^l≤f_i ^l，max；

f^e，min≤f_i ^l≤f^e，max；

T_i ^loca(f_i ^l)≤τ_i；

π_i∈{0，1}；

其中，π＝{π₁，π₂，...，π_N}表示任务卸载决策集合，N为请求任务处理的用户数；

表示本地CPU频率资源分配策略集合，f_i ^l表示第i个用户的本地CPU频率资源分配策略；

表示所有组合场景下MEC服务器CPU频率资源分配策略集合，其中

表示第i个用户在MEC服务器的排队等待时间的集合；p＝{p₁，p₂，...p_N}表示用户任务传输功率分配策略集合，p_N表示第N个用户任务传输功率分配策略；π_i表示第i个用户的卸载决策；κ_i为第i个用户芯片构架相关的有效能量系数；D_i为待计算任务的任务量大小；L_i表示单位bit任务所需要的CPU周期数；R_i为第i个用户将任务卸载到MEC服务器的传输速率；γ为与MEC服务器芯片架构相关的有效能量系数；τ_i表示计算完成该任务的最大时延要求；

表示组合场景的期望；

和

分别表示传输功率的最小值和最大值；f_i ^l，min和f_i ^l，max分别表示本地CPU频率的最小

值和最大值；f^e，min和f^e，max分别表示MEC服务器CPU频率的最小值和最大值；T_i ^wait(ω)表示组合实现为ω时第i个用户的排队等待时间；T_i ^loca(f_i ^l)为第i个用户本地计算的时延；T_i ^tran(p_i)为第i个用户选择将任务卸载到MEC服务器处理时任务上传的通信时延；

为第i个用户在MEC服务器的计算时延。

4.根据权利要求1所述的一种不确定网络环境中的任务卸载和资源分配方法，其特征在于，利用随机模拟方法将基于两阶段随机规划的任务卸载和资源分配问题转化为样本均值近似问题包括：

约束条件：

f_i ^l，min≤f_i ^l≤f_i ^l，max；

f^e，min≤f_i ^l≤f^e，max；

T_i ^loca(f_i ^l)≤τ_i；

π_i∈{0，1}；

其中，π＝{π₁，π₂，...，π_N}表示任务卸载决策集合，N为请求任务处理的用户数；

表示本地CPU频率资源分配策略集合，f_i ^l表示第i个用户的本地CPU频率资源分配策略；

表示所有组合场景下MEC服务器CPU频率资源分配策略集合，其中

表示第i个用户在MEC服务器的排队等待时间的集合；p＝{p₁，p₂，...p_N}表示用户任务传输功率分配策略集合，p_N表示第N个用户任务传输功率分配策略；π_i表示第i个用户的卸载决策；κ_i为第i个用户芯片构架相关的有效能量系数；D_i为待计算任务的任务量大小；L_i表示单位bit任务所需要的CPU周期数；R_i为第i个用户将任务卸载到MEC服务器的传输速率；γ为与MEC服务器芯片架构相关的有效能量系数；τ_i表示计算完成该任务的最大时延要求；

和

分别表示传输功率的最小值和最大值；f_i ^l,min和f_i ^l,max分别表示本地CPU频率的最小值和最大值；f^e,min和f^e,max分别表示MEC服务器CPU频率的最小值和最大值；T_i ^wait(ω_s)表示组合实现为ω_s时第i个用户的排队等待时间；T_i ^loca(f_i ^l)为第i个用户本地计算的时延；T_i ^tran(p_i)为第i个用户选择将任务卸载到MEC服务器处理时任务上传的通信时延；

为第i个用户在MEC服务器的计算时延；S＝K^N表示该样本组合场景的数量，K为从第i个用户的场景集合Ω_i中抽取的独立同分布样本数量。

5.根据权利要求4所述的一种不确定网络环境中的任务卸载和资源分配方法，其特征在于，若第i个用户的任务卸载决策为0，即第i个用户将任务在本地处理，则任务在本地计算时最优CPU频率分配策略为：

其中，f_i ^l*为任务在本地计算时最优CPU频率分配策略。

6.根据权利要求4所述的一种不确定网络环境中的任务卸载和资源分配方法，其特征在于，若第i个用户的任务卸载决策为1，即第i个用户将任务卸载到服务器计算，则通过遗传算法计算最优的卸载策略，具体包括：

将每一个可行的传输功率进行浮点向量编码，每一个浮点向量表示一个染色体，浮点向量维数与解向量维数一致；

从用户i传输功率的可行域中随机产生一个点，并检验其是否满足用户传输时延小于最大传输时延要求，如果满足则作为一个染色体，否则，重新产生一个随机点，直到满足约束条件，重复以上过程M次，产生M个染色体作为初始种群，记为p_i,1,p_i,2,...,p_i,M；

对于每个染色体p_i,m，求解其适应度函数，并根据适应度求其评价函数，且对于每个染色体评价函数值越大，被选择产生后代的概率越大，其中适应度为样本均值近似问题函数；

对于染色体p_i,m，计算前m个染色的累积概率

从区间(0,q_i,M)中生成一个随机数r，若满足q_i,m-1＜r＜q_i,m，q_i,0＝0，则选择染色体p_i,m作为父本繁殖下一代，依此规则选择M个染色体；

定义P_c为交叉概率，则种群中有期望值为P_c×M个染色体进行交叉操作，即假设p_i,1和p_i,2为要交叉的两个染色体，则从区间(0,1)中生成一个随机数c，直到p'_i,1＝cp_i,1+(1-c)p_i,2和p'_i,2＝(1-c)p_i,1+cp_i,2满足用户传输时延小于最大传输时延要求，则p'_i,1和p'_i,2将代替原染色体p_i,1和p_i,2成为两个新的染色体；

定义P_u为变异概率，则种群中有期望值为P_u×M个染色体进行变异操作，即假设p_i,1为要变异的染色体，随机选择一个变异方向d，从

中生成一个随机数α，直到p'_i,1＝p_i,1+αd满足约束条件，则p'_i,1将取代原染色体p_i,1成为新的染色体；

经过选择、交叉和变异操作，可以得到一个新的种群，并准备进行下一代进化，如果上述步骤达到了给定的循环次数，则遗传算法终止；算法终止后，从最后一代选择一个适应度最高的染色体，即得到优化问题的全局最优解。

7.根据权利要求6所述的一种不确定网络环境中的任务卸载和资源分配方法，其特征在于，每个染色体相应的适应度值，并按照适应度值从小到大对染色体进行排序，则染色体p_i,m的评价函数表示为：

eval(p_i,m)＝a(1-a)^m-1,m＝1,2,...,M；

其中，eval(p_i,m)为染色体p_i,m的评价函数；a为常数，且a∈(0,1)。

8.根据权利要求1所述的一种不确定网络环境中的任务卸载和资源分配方法，其特征在于，最优卸载决策表示为：

其中，

为最优卸载决策；

为任务在本地计算时最优CPU频率分配策略；

为每个场景下MEC服务器最优CPU频率分配策略；

为第i个用户的最优发射功率。

9.根据权利要求7所述的一种不确定网络环境中的任务卸载和资源分配方法，其特征在于，每个场景下MEC服务器最优CPU频率分配策略

表示为：

其中，

为第i个用户选择发射功率为p_i,m时将任务卸载到MEC服务器处理时任务上传的通信时延。

Images