Energy Consumption Optimization Algorithm for Full-Duplex Relay-Assisted Mobile Edge Computing Systems
-
摘要: 为缓解终端设备处理大数据量、低时延业务的压力，该文提出一种基于全双工中继的移动边缘计算网络资源分配算法。首先，在满足计算任务时延约束、用户最大计算能力、用户和中继的最大发射功率约束条件下，考虑中继选择与子载波分配因子、用户任务卸载系数、用户与中继的传输功率的联合优化，建立了系统总能耗最小化问题。其次，利用交替迭代和变量代换的方法，将原非凸问题分解为两个凸优化子问题，并利用内点法和拉格朗日对偶原理分别进行求解。仿真结果表明，所提算法具有较低的能量消耗。Abstract: In order to alleviate the pressure of terminal devices to deal with the big-data and low-delay services, a resource allocation algorithm is studied for mobile edge computing networks with full-duplex relays. Firstly, the constraints of the maximum task latency, the maximum computing ability of users, and the maximum transmit power of users and relays are considered for achieving total energy consumption minimization by jointly optimizing the relay selection and subcarrier allocation factor, user task offloading coefficient, and the transmission power of users and relays. Secondly, based on the alternating iteration method and the variable-substitution approach, the originally non-convex problem is decomposed into two convex subproblems, which are solved by using the interior-point method and Lagrange dual theory, respectively. Simulation results show that the proposed algorithm has low energy consumption.
-
Key words:
- Full-duplex based relays /
- Mobile edge computing /
- Resource allocation
-
1. 引言
雷达信号识别是雷达威胁告警(Radar Warning Receiver, RWR)、雷达支援侦察(Electronic Support Measurement, ESM)、雷达情报侦察(ELetronic INTelligence, ELINT)的重要内容，是雷达对抗首先要解决的关键性问题。正确识别雷达信号可以帮助己方对目标雷达类型、载体、用途、威胁等级的识别，进而影响进一步的作战决策[1]。
在过去十几年里，雷达信号识别受到了非常广泛的研究，纵观这些研究主要集中在筛选(提取)能有效表征辐射源信号特征，构建新的识别特征参数体系上。这些方法如果按处理对象的不同，可以粗略分成两条研究路线：一条是实数序列处理，另一条是符号序列的处理。尽管雷达信号的实数序列处理仍然是最主要的分析手段，并且得到非常多的研究成果[2–8]。但是，随着雷达信号的日益复杂，一些尝试从符号序列分析的角度进来挖掘雷达数据隐含信息的方法也受到关注[9–11]，如文献[9]是在雷达信号符号化的基础上，提取了反映序列随机程度的信息熵进行雷达信号识别，并获得较好的性能；文献[11]进一步研究了对于相同雷达信号，不同调制参数的可分性。这些方法的都是在实数序列符号化的基础上进行特征提取，但是，所用的符号化方法过于简单，容易损失信号的大量信息导致提取的特征存在偏差。所以，选择一种重构误差较小的符号化方法是非常重要的。
纵观当前关于符号化方法的研究中，符号化聚合近似(Symbolic Aggregate approXimation, SAX)方法被认为是一种简单且非常有效的实数序列符号化的方法[12,13]，它可以有效地降低序列的维度，同时能保留原始数据的大部分信息(重构误差小)，并且具有一定的降噪能力，该方法已经在许多时间序列分析上得到应用。因此，在本文所提雷达信号识别框架中，首先使用SAX方法将雷达信号序列转化为符号化序列，再对符号化序列进行统计学分析。在对符号化序列的统计分析研究上，本文提出了一种基于信息理论的多尺度信息熵(Multi-Scale Information Entropy, MSIE)特征，该特征是多个字符集尺度下提取的信息熵集合，它反映了符号化序列在不同尺度上的复杂度特性。具体的提取方法将在第3节进行详细阐述。
本文的结构安排如下：第2节详细介绍SAX符号化方法的理论与方法，在此基础上，第3节给出多尺度信息熵的定义和提取方法，并且分析其抗噪性能；然后，在第4节对文中提出的雷达信号识别模型进行详细介绍。最后，在第5节进行仿真实验的设计和分析。
2. 实数序列符号化
面对越来越复杂的雷达信号数据，一些传统的实数序列处理方法受到了一定的挑战，因此，一些学者尝试从另一种途径来处理这些复杂数据的分析问题，即实数序列符号化。实数序列符号化具有下列两个优势[14]：
(1) 序列的符号化是一个对原始数据“粗粒化”的过程，它可以去除数据无关紧要的非本质特征，留下数据关键的特征信息；
(2) 符号化序列在计算速度上具有很大的优势。
实数序列符号化的关键问题是如何选取合适的符号化策略，但是目前并没有一个指导性的结论。目前，普遍认为SAX方法是一种非常有效的符号化方法，因此，本文选用该方法对雷达信号序列进行符号化处理。下面对SAX方进行详细描述。
2.1 分段聚集近似算法
SAX方法是实数序列符号化的一个经典算法，是由文献[15]在分段聚集近似(Piecewise Aggregate Approximation, PAA)算法的基础上提出的。PAA算法是一种数据的重表达方法，它能基本保留数据原有形状，并且算法时间复杂度低，因此，PAA算法被广泛应用在时间序列分析问题上。经典的PAA算法可以描述如下：
将一个长度为n的序列
Q={qi,i=1,2,···,n} 映射到m维空间P={pi,i=1,2,···,m} ，序列P的第i个元素可以通过式(1)计算pi=mnn/mi∑j=n/m(i−1)+1qj,1≤i≤m (1) 式(1)可以简单描述为，序列Q被分成m个大小相等的“段”，然后计算每一个“段”的均值来构建新的序列P。从式(1)可以发现，n需要满足m的整数倍。然而，现实情况很难满足数据可分的情况，所以，当
n/m 不是整数时可以进行如下操作[12]：当
1<i<m 时pi=mn((1−ci−1)qt(i−1)+t(i−1)+b(i)∑j=t(i−1)+1qj+ciqt(i−1)+b(i)+1) (2) 假设
F(∙) 为对数值“∙ ”取整数部分的函数，则式(2)中，t(i)=t(i−1)+b(i−1)+1 ，并且令t(1)=a+1 ，其中a=F(n/m) ，另外，b(i)=F(n/m+ci−1−1) 。所以，当i=1 或i=m 时pi={mn(a∑j=1qj+c1qa+1),i=1mn(cmqn−a+n∑j=n−a+1qj),i=m (3) 其中，序列
{ci,i=1,2,···,m} 是数据点的权值，并且c1=cm=G(n/m) ,G(∙) 是对数值“∙ ”取小数部分的函数。同时，式(2)和式(3)满足每一个“段”的数据权值之和等于n/m 。为了更形象地描述上面的公式，下面列举一个简单的示例：假设要将序列{k,k=1,2,···,12} 划分成包含5个数值的新序列，记为{si,i=1,2,···,5} ，因为G(12/5)=2/5 ，所以，当i=1 ,i=5 时，s1=(1+2+(2/5)3) (5/12)=1.75 ,s5=((2/5)10+11+12)(5/12)= 11.25。当1<i<5 时，s2=((3/5)3+4+(4/5)5) (5/12)≈4.08 ,s3=((1/5)5+6+7+(1/5)8)(5/12) =6.5 ,s4=((4/5)8+9+(3/5)10)(5/12)≈8.92 。通过上面的计算，可以获得一个新序列
{1.754.086.58.9211.25} 。新序列的每一个元素都包含原始序列12/5=2.4 个数据点，这样就解决了序列不可分的问题。2.2 SAX符号化方法
直接根据实数序列的数值特点进行符号划分是最常见的符号化方法，最常用的是静态划分方法，它需要先选定一个大小为k的字符集，假设为
{12···k} ，然后根据式(4)可以将时间序列Q={qi,i=1,2,···,n} 符号化为符号序列S={si,i=1,2,···,w} ，其中w≪n ，从而实现数据的降维。si={1,−∞<qi<β12,β1≤qi<β2⋮k,βk<qi<+∞ (4) 其中，
{β1β2···βk} 是对原始序列Q进行空间划分的断点集合。实数序列符号化的关键在于如何确定这一断点集合。通行的方法有等区间法和等概率法，SAX采用等概率法获得不同的分区，从信息论的角度来说，等概率法有利于保留数据的有效信息。SAX具体的算法可以简单描述为：在假设时间序列满足正态分布的条件下，(1)对时间序列标准化，将一般正态分布的序列转变为标准正态分布；(2)对标准化的时间序列用PAA算法进行分段划分，这里需要设置参数w大小，表示划分后序列的维数；(3)根据设定的符号集大小(参数a)，查询等概率断点表(表1给出了参数a在3～8之间的断点划分)确定断点集合，最后完成时间序列的符号化。
表 1 参数a从3～8的等概率断点查询表[12]断点(βi) 字符集大小(a) 3 4 5 6 7 8 β1 0.43 0.67 0.84 0.97 1.07 1.15 β2 0.43 0 0.25 0.43 0.57 0.67 β3 – 0.67 0.25 0 0.18 0.32 β4 – – 0.84 0.43 0.18 0 β5 – – – 0.97 0.57 0.32 β6 – – – – 1.07 0.67 β7 – – – – – 1.15 图1是6类常见雷达信号(包括具有恒定频率的常规信号(CP)、线性调频信号(LFM)、二相编码信号(BPSK)、二进制频率编码信号(2FSK)、非线性调制信号(NLFM)和COSTAS编码信号)的SAX符号化示例。其中，字符集大小a=5；符号化后的维度w=25。从图1可以看出，不同调制类型的雷达信号，它的符号化序列有明显区别，因此，在这些符号化序列的基础上提取的统计特征具有良好的类别可分性。
3. 多尺度信息熵的特征提取
统计不同字符在符号序列中出现的相对频率是分析符号化序列常见方法，一种被广泛应用的频率统计方法是修正的信息熵，它定义为[9]
H=−1lg(ns)ns∑i=1pilg(pi) (5) 其中，
pi 是第i个字符出现的概率，实际中可以用频率来近似，即si/N (si 是第i个字符出现的个数；N 是符号化序列的长度)；ns 是符号化序列中出现的不同字符集大小。修正的信息熵刻画了离散数据的随机程度，对于完全随机的序列，修正信息熵为1，而对于不完全随机的序列，它的修正信息熵在0～1之间。在修正信息熵的基础上，本文提出了一种多尺度信息熵(Multi-Scale Information Entropy, MSIE)的提取方法，用于描述符号化序列在不同字符集尺度上的无规则程度。因此，字符集大小为j的修正信息熵定义为
Hj=−1lg(j)nj∑i=1εj,ilg(εj,i),j=a1,a2,···,am (6) 式(6)中，
εj,i 是在该尺度下第i个字符出现的概率；a1 和am 分别是字符集大小的下界和上界，m为字符集大小区间。令T为多尺度信息熵向量，选定一个字符集大小区间，假设为
{a1a2···am} ，则T 的计算方法如下：(1) 对于一个实数序列
Q={qi,i=1,2,···,n} ，采用SAX对其进行符号化，可以得到m个符号化序列，记为{St,t=1,2,···,m} ；(2) 对每一个符号化序列
St 依据式(6)计算其在第t个字符集尺度下的修正信息熵，得到m个不同字符集尺度下的修正信息熵，级联这些值所组成的特征向量称为多尺度信息熵，记为T={Ha1,Ha2,···Ham} 。对于SAX符号化方法，如果对超参数a和w设置合适的值，符号化后的序列可以很好地保留原序列的结构信息，同时摒弃一些细微的信息。因此，SAX方法本身存在一定的降噪能力，在此基础上，提取的统计量参数(修正的信息熵)也就具一定的抗噪性能。
为了更好地说明多尺度信息熵的抗噪性能，下面用文献[9]中的仿真信号来分析不同字符集大小下的修正信息熵的抗噪性能，该仿真信号为
S1(t)=2sin(30πt)+sin(150πt)+5cos(21.6πt) (7) 图2信号S1在不同字符集尺度(a在3～6之间)下的信息熵与信噪比(Signal to Noise Ratio, SNR)的关系图。其中，SNR在–10～20 dB范围内。从图2中可看出，随着SNR的升高，信息熵的值趋于稳定，在SNR>0 dB时，信息熵的值基本保持不变，所以，信息熵具有较好的抗噪性能。因此，多尺度信息熵也就具有较强的抗噪性能。
4. 雷达信号识别模型
综合上面分析和介绍，本文提出了如图3的雷达信号识别方法。其中，数据预处理阶段主要对信号进行消噪处理。考虑到不同调制类型的雷达信号在频率域有不同的幅度谱，同时，在频率域进行特征提取可以削弱载频变化的影响[16]，本文以雷达信号的频率域序列作为处理对象。在频率域，由于有用信息的能量集中在较窄的频带范围，而噪声能量分布在整个频带上，所以可以按式(8)对信号序列进行消噪处理。令原始雷达信号序列为
{ri,i=1,2,···,n} ，并假设消噪后的信号序列为{fi,i=1,2,···,n} ，则移动平均滤波的计算公式为fi=ri+ri+1+···+ri+d−1d (8) 其中，d表示移动窗口大小。由于雷达信号的幅度谱是对称的，所以可以只取左边的幅度谱。
在选择分类算法的问题上，考虑到k-NN算法具有实现简单，易于理解，且对于特征维数较低的多分类问题通常具有较好的性能。因此，选用k-NN算法作为信号的分类器，该算法的具体实现步骤如下：
步骤 1 给定一个未知的实数序列Q，计算Q和数据库(训练集)中所有样本的相似度，通常采用欧式距离作为相似性度量方法；
步骤 2 寻找k个和Q最为相似的样本，作为未知对象Q的k个近邻；
步骤 3 在这k个近邻中出现次数最多的类别即为对象Q的类别。
但是，k-NN算法也面临k值选择的问题，如果k值过大，对于不平衡数据集的分类性能较差；而当k值过小时，算法对噪声较为敏感。根据以往的实验经验，在接下来的仿真实验中，设置超参数k=1。
另外，需要注意的是，使用SAX符号化方法需要确定符号集的大小a和符号化后的维数w的值，然而，目前没有这方面的理论研究，尽管文献[12]给出了参数a的经验性选择范围，但是实际应用中还需要根据不同的数据集进行考虑。一般来说，两个超参数都可以通过交叉验证来选择。
5. 仿真实验
通过仿真6类常见的辐射源信号来验证本文方法的有效性，包括CP, LFM, BPSK, BFSK, NLFM和COSTAS编码信号。BPSK和BFSK采用13位Barker码，BFSK信号的两个频点分别为10 MHz和12 MHz; LFM的频率变化范围是10～15 MHz; NLFM采用正弦波调制。并且，各辐射源信号的载频为10 MHz，采样频率为60 MHz，脉宽为10 μs。同时，考虑到实际情况中存在的测量误差，在仿真信号时将载频的均方根误差设置为1 MHz；将脉宽的均方根误差设置为1 μs。
仿真下列一组数据进行测试：在20～25 dB的信噪比范围内，每种信号样式随机产生20个脉冲信号，总计120个脉冲信号作为模型的训练样本。同时，分别在5 dB, 10 dB, 15 dB, 20 dB信噪比下各仿真1200个脉冲信号作为测试样本，每一信号样式各200个脉冲。实验中的信号附加的噪声为高斯白噪声。本次实验是在Matlab R2010b平台下完成的。
根据第3节介绍的方法对上面仿真的雷达信号提取多尺度信息熵(设置a的范围为3～9；w=140)，该特征是一个7维的向量，为了便于观察多维数据集的分布情况，这里采用文献[17]提出的t-分布式随机邻域嵌入(t-distributed Stochastic Neighbor Embedding, t-SNE)方法进行数据可视化，将7维的数据映射到2维空间中。图4是在信噪比为5～20 dB范围内6类雷达信号数据的分布图，从图中可以看出，除了BPSK信号和BFSK信号出现小部分重叠无法分开外，多尺度信息熵的类间分离性较好，具有较好的类别可分性，因此，采用k-NN(k =1)分类器可以获得良好的性能。
表2是本文所提方法的识别结果，从表中可以看出，在信噪比为15 dB时，6类信号的识别正确率都达到90%以上，说明了本文方法的有效性。另外，从图中还可以看出，CP, LFM, BPSK, NLFM和COSTAS信号在信噪比为5 dB时，仍保持98%以上的正确识别率，这是因为这5类雷达信号的幅度谱明显不同，符号化后提取的多尺度信息熵就具有较好的类别可分性。
表 2 不同SNR下6种雷达信号的识别率雷达信号 信噪比SNR (dB) 20 15 10 5 LFM 1.000 1.000 1.000 0.985 CP 1.000 1.000 1.000 1.000 BPSK 0.975 0.990 0.990 1.000 BFSK 0.930 0.910 0.800 0.700 NLFM 1.000 1.000 1.000 1.000 COSTAS 1.000 1.000 1.000 1.000 为了圈定本文所提方法的适用范围，需要进一步增加实验进行对比。因此，为上述6类雷达信号增加2种雷达信号，包括四相编码信号(QPSK)和Chirp子脉冲步进频率信号(CSF)。其中，QPSK采用16位Frank码调制。8类信号在20 dB信噪比下仿真800个样本，每类信号包含100个脉冲作为测试集，在20～25 dB的信噪比范围内，每种信号样式随机产生20个脉冲信号，总计180个脉冲信号作为模型的训练样本，实验中用交叉验证选择字符集大小a=11，以及符号化后的维度w=160。考虑到QPSK信号和BPSK信号的幅度谱较为相似，因此，理论上这两种信号提取的多尺度信息熵存在部分重叠，分类器将难以区分这两类信号，从图5中的混淆矩阵也可以明显看出，QPSK和BPSK的区分度较差。
传统的方法可以较好地区分BPSK与QPSK两种信号，但是一般需要将雷达信号转换到时-频域上提取特征，而大多数基于时-频变换方法或多或少的存在问题[16]，并且受噪声的影响严重，容易导致时-频信号失真，而傅里叶变换能获得稳定的频率域序列。另外，基于时-频变换方法的时间复杂度通常要高于傅里叶变换算法。这里选择一种经典的时-频域特征—小波脊频级联特征(Wavelet Ridge Frequency Cascade-Connection Feature, WRFCCF)[16]进行实验对比，该特征是一个7维的特征向量，其中每个特征是雷达信号在时-频域上的统计量。首先，对提取WRFCCF和MSIE所耗费的时间进行对比，实验采用上述的100个训练样本进行实验，结果重复10次取平均(如表3所示)。从表中可以看出，提取WRFCCF所耗费的时间要远远高于提取多尺度信息熵的时间。
表 3 提取两种特征耗费的时间对比特征向量 耗费时间(s) WRFCCF 135.102 MSIE 1.704 使用小波脊频级联特征，并选择k-NN分类器，可以得到如图6的混淆矩阵。从图中可以看出，WRFCCF能有效的区分BPSK, QPSK两种信号，但是对于其它调制类型的雷达信号，如图6中的CP、BFSK的区分度较差，而从图5可知，本文所提方法除了BPSK, QPSK外，其它信号的识别率基本上达到了100%(如表4，本文方法获得更好的总体识别率)，这是因为多尺度信息熵是描述雷达信号的幅度谱的结构形态，当雷达信号的幅度谱存在较弱相似性的情况下，可以获得好的识别性能。
表 4 两种方法的总体识别正确率比较(%)识别方法 总体识别率 WRFCCF+k-NN 92.13 MSIE+k-NN 95.63 多尺度信息熵是一种典型的频率域特征，这里选择两种经典的频率域特征进行识别性能的比较，包括文献[2]提出的盒维数和稀疏性，以及文献[9]提出的信息熵特征。盒维数和稀疏性这两个复杂度特征(complexity characteristics)与信息熵类似，是从不同层面来描述信号的结构。同样采用t-SNE方法对数据进行可视化，图7是在信噪比为5～20 dB范围内6类基于复杂度的雷数据分布图，显然，复杂度特征的类间分离性较差，分布图中有许多数据出现重叠。
在测试特征的识别性能之前，首先对特征提取算法的计算量进行对比。提取盒维数特征时，需要计算序列的绝对值差分，所以算法的执行次数
fb(n)=N+(N+1)/2 ，其中，N是序列的长度；提取稀疏性特征时，首先需要先对序列进行二进制符号化，将序列划分为L个长度为N0 的子序列，统计窗内序列中“1”元素的占比，最后求和。所以，该算法的执行次数fs(n)=N+LN0 。假设样本集的大小为M，则复杂度特征提取算法的时间复杂度为T(n)=O(M(fb(n)+fs(n)))=O(M(N+(N+1)/2+N+LN0))=O(MN)。 (9) SAX是在PAA的基础上进行符号化，而且容易知道PAA算法的时间复杂度为
O(L1) ，其中L1 是算法划分的段数。另外，修正信息熵是统计不同字符出现的相对频率，因此，它的时间复杂度也为O(L1) 。所以，多尺度信息熵的时间复杂度为T(n)=O(cML1) ，其中，c是字符集尺度的取值范围，实际情况中通常较小，可以认为cL1≈N ，则多尺度信息熵的时间复杂度约等于O(MN) 。另外，由文献[9]中关于信息熵的提取算法，可知其时间复杂度为O(MN) 。通过上面的分析可知，3种特征具有相同的时间复杂度，能在合理的时间内提取特征。表5分别是3种方法在5 dB, 10 dB, 15 dB, 20 dB信噪比下的总体识别正确率。为了便于描述，这里将本文所提方法表述为“MSIE+k-NN”；将文献[2]提出的方法表述为“CC+k-NN”；将文献[9]提出的方法表述为“SIE+k-NN”。3种方法的总体识别率如表所示，从表中可以容易看出，本文方法相比两种经典的频率域特征具有更好的识别正确率。
表 5 3种方法的总体识别正确率比较(%)识别方法 信噪比SNR(dB) 20 15 10 5 MSIE+k-NN 98.42 97.25 94.25 91.25 CC+k-NN 80.25 73.08 54.33 <50 SIE+k-NN 81.42 79.08 71.25 62.92 6. 结论
本文提出了一种基于多尺度信息熵的雷达信号识别方法，能够较好地对雷达信号进行识别。本文方法将信号从实数域映射到符号域，去除信号中的细微特征的同时保留其结构特征，从而提取具有较好抗噪性能的特征。通过仿真实验表明，本文方法在雷达信号的幅度谱存在较弱相似性的情况下，具有很好的识别率，并且优于传统的基于复杂度特征的识别方法，说明了本文方法具有一定的参考价值。
1)基于自干扰消除技术，本文假设干扰消除残留项可以建模为噪声功率。 -
算法1 基于交替迭代的资源分配算法 1.初始化系统参数：pmaxn,Pmaxm, Tmaxn, Fmaxn, κ, B, K，hkn,m, σ2, gkm, fMn, Sn, Cn, φ, σ2SI；定义交替迭代算法收敛精度ℓ；初始化交替迭
代次数t；定义外层最大迭代次数Tmax；定义Dinkelbach迭代算法收敛精度ζ以及相应最大迭代次数Lmax；初始化梯度迭代次数l，初始化
q=0；
2.给定αkn,m(t)和ˉpkn,m(t)，利用内点法求解问题式(18)，得到xkn,m(t+1)；
3.给定xkn,m(t+1)，给定q，求解问题式(20)，得到当前最优值ˉpkn,m(l)，αkn,m(l)；
4.当|xkn,m(t+1)Snˉpkn,m(l)−qˉRkn,m(l)|≥ζ，或者l≤Lmax；
5.令Flag=0，更新l=l+1；
6.将q更新为q=xkn,m(t+1)Snˉpkn,m(l)/xkn,m(t+1)Snˉpkn,m(l)ˉRkn,m(l)ˉRkn,m(l)，结束并执行步骤3；
7.当|xkn,m(t+1)Snˉpkn,m(l)−qˉRkn,m(l)|≤ζ，或者l=Lmax;
8.令Flag=1，更新并输出αn,m(t+1) = αln,m(l)，ˉpkn,m(t+1)=ˉpkn,m(l)；
9.当|{ELn(t+1)+αkn,m(t+1)EUTn(t+1)}−{ELn(t)+αkn,m(t)EUTn(t)}|≥ℓ，或者t≤Tmax；
10.更新t=t+1，执行步骤2；
11.结束并输出。
下载: 导出CSV
表 1 仿真参数
参数 值 参数 值 ϕ/(cycles⋅bit−1) 40 κ 10−24 Fmaxn/cycles 109 Pmaxm/W 5 fMn/cycles 1010 pmaxn/W 1 K 5 χ 3 σ2/mW 10−6 B/Hz 106
下载: 导出CSV
-
[1] LIU Yaqiong, PENG Mugen, SHOU Guochu, et al. Toward edge intelligence: Multiaccess edge computing for 5G and internet of things[J]. IEEE Internet of Things Journal, 2020, 7(8): 6722–6747. doi: 10.1109/JIOT.2020.3004500 [2] MAO Yuyi, YOU Changsheng, ZHANG Jun, et al. A survey on mobile edge computing: The communication perspective[J]. IEEE Communications Surveys & Tutorials, 2017, 19(4): 2322–2358. doi: 10.1109/COMST.2017.2745201 [3] 谢人超, 廉晓飞, 贾庆民, 等. 移动边缘计算卸载技术综述[J]. 通信学报, 2018, 39(11): 138–155. doi: 10.11959/j.issn.1000-436x.2018215XIE Renchao, LIAN Xiaofei, JIA Qingmin, et al. Survey on computation offloading in mobile edge computing[J]. Journal on Communications, 2018, 39(11): 138–155. doi: 10.11959/j.issn.1000-436x.2018215 [4] 李国权, 徐勇军, 陈前斌. 基于干扰效率多蜂窝异构无线网络最优基站选择及功率分配算法[J]. 电子与信息学报, 2020, 41(2): 957–964. doi: 10.11999/JEIT190419LI Guoquan, XU Yongjun, and CHEN Qianbin. Interference efficiency-based base station selection and power allocation algorithm for multi-cell heterogeneous wireless networks[J]. Journal of Electronics &Information Technology, 2020, 41(2): 957–964. doi: 10.11999/JEIT190419 [5] ZHANG Jing, XIA Weiwei, YAN Feng, et al. Joint computation offloading and resource allocation optimization in heterogeneous networks with mobile edge computing[J]. IEEE Access, 2018, 6: 19324–19337. doi: 10.1109/ACCESS.2018.2819690 [6] WANG Jingrong, LIU Kaiyang, and PAN Jianping. Online UAV-mounted edge server dispatching for mobile-to-mobile edge computing[J]. IEEE Internet of Things Journal, 2020, 7(2): 1375–1386. doi: 10.1109/JIOT.2019.2954798 [7] QIAN Liping, FENG Anqi, HUANG Yupin, et al. Optimal SIC ordering and computation resource allocation in MEC-aware NOMA NB-IoT networks[J]. IEEE Internet of Things Journal, 2019, 6(2): 2806–2816. doi: 10.1109/JIOT.2018.2875046 [8] SUN Haijian, ZHOU Fuhui, and HU R Q. Joint offloading and computation energy efficiency maximization in a mobile edge computing system[J]. IEEE Transactions on Vehicular Technology, 2019, 68(3): 3052–3056. doi: 10.1109/tvt.2019.2893094 [9] SONG Zhengyu, LIU Yuanwei, and SUN Xin. Joint radio and computational resource allocation for NOMA-based mobile edge computing in heterogeneous networks[J]. IEEE Communications Letters, 2018, 22(12): 2559–2562. doi: 10.1109/LCOMM.2018.2875984 [10] REN Jinke, YU Guanding, HE Yinghui, et al. Collaborative cloud and edge computing for latency minimization[J]. IEEE Transactions on Vehicular Technology, 2019, 68(5): 5031–5044. doi: 10.1109/TVT.2019.2904244 [11] CHEN Xihan, CAI Yunlong, SHI Qingjiang, et al. Efficient resource allocation for relay-assisted computation offloading in mobile-edge computing[J]. IEEE Internet of Things Journal, 2020, 7(3): 2452–2468. doi: 10.1109/JIOT.2019.2957728 [12] XIE Biyuan, ZHANG Qi, and QIN Jiayin. Joint optimization of cooperative communication and computation in two-way relay MEC systems[J]. IEEE Transactions on Vehicular Technology, 2020, 69(4): 4596–4600. doi: 10.1109/TVT.2020.2973292 [13] LAN Yanwen, WANG Xiaoxiang, WANG Dongyu, et al. Mobile-edge computation offloading and resource allocation in heterogeneous wireless networks[C]. Proceedings of 2019 IEEE Wireless Communications and Networking Conference, Marrakesh, Morocco, 2019: 1–6. doi: 10.1109/WCNC.2019.8885793. [14] ALAM S, MARK J W, and SHEN X S. Relay selection and resource allocation for multi-user cooperative OFDMA networks[J]. IEEE Transactions on Wireless Communications, 2013, 12(5): 2193–2205. doi: 10.1109/TWC.2013.032113.120652 [15] WEN Dingzhu, YU Guanding, LI Rongpeng, et al. Results on energy- and spectral-efficiency tradeoff in cellular networks with full-duplex enabled base stations[J]. IEEE Transactions on Wireless Communications, 2017, 16(3): 1494–1507. doi: 10.1109/TWC.2016.2647593 [16] CHEN Lei, YU F R, JI Hong, et al. Distributed virtual resource allocation in small-cell networks with full-duplex self-backhauls and virtualization[J]. IEEE Transactions on Vehicular Technology, 2016, 65(7): 5410–5423. doi: 10.1109/TVT.2015.2469149 [17] DINKELBACH W. On nonlinear fractional programming[J]. Management Science, 1967, 13(7): 492–498. doi: 10.1287/mnsc.13.7.492 [18] BOYD S and VANDENBERGHE L. Convex Optimization[M]. Cambridge: Cambridge University Press, 2004. [19] HUANG Kaizhi and ZHANG Bo. Robust secure transmission for wireless information and power transfer in heterogeneous networks[J]. IET Communications, 2019, 13(4): 411–422. doi: 10.1049/iet-com.2018.5529 [20] CROUZEIX J P and FERLAND J A. Algorithms for generalized fractional programming[J]. Mathematical Programming, 1991, 52(1/3): 191–207. doi: 10.1007/BF01582887 -
下载:
计量
- 文章访问数: 851
- HTML全文浏览量: 470
- PDF下载量: 76
- 被引次数: 0
下载:
