Loading [MathJax]/jax/element/mml/optable/BasicLatin.js
高级搜索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于带缓存的云接入网络最优能效设计

孙远 李春国 黄永明 杨绿溪

孙远, 李春国, 黄永明, 杨绿溪. 基于带缓存的云接入网络最优能效设计[J]. 电子与信息学报, 2019, 41(7): 1525-1532. doi: 10.11999/JEIT180722
引用本文: 孙远, 李春国, 黄永明, 杨绿溪. 基于带缓存的云接入网络最优能效设计[J]. 电子与信息学报, 2019, 41(7): 1525-1532. doi: 10.11999/JEIT180722
Yuan SUN, Chunguo LI, Yongming HUANG, Lüxi YANG. Optimal Energy-efficient Design for Cache-based Cloud Radio Access Network[J]. Journal of Electronics & Information Technology, 2019, 41(7): 1525-1532. doi: 10.11999/JEIT180722
Citation: Yuan SUN, Chunguo LI, Yongming HUANG, Lüxi YANG. Optimal Energy-efficient Design for Cache-based Cloud Radio Access Network[J]. Journal of Electronics & Information Technology, 2019, 41(7): 1525-1532. doi: 10.11999/JEIT180722

基于带缓存的云接入网络最优能效设计

doi: 10.11999/JEIT180722
基金项目: 国家863计划(2015AA01A703),国家自然科学基金(61372101, 61671144)
详细信息
    作者简介:

    孙远:男,1988年生,博士生,研究方向为云接入网络资源分配

    李春国:男,1983年生,教授,博士生导师,研究方向为新一代无线通信理论与技术,基于人工智能的视频信息处理,水下无线通信

    黄永明:男,1977年生,教授,博士生导师,研究方向为下一代移动通信技术,毫米波MIMO

    杨绿溪:男,1964年生,教授,博士生导师,研究方向为下一代移动通信技术,信号与信息处理

    通讯作者:

    杨绿溪 lxyang@seu.edu.cn

  • 中图分类号: TN92

Optimal Energy-efficient Design for Cache-based Cloud Radio Access Network

Funds: The National 863 Project of China (2015AA01A703), The National Natural Science Foundation of China (61372101, 61671144)
  • 摘要: 在云接入网络(Cloud-RAN)中,现有工作大多假定射频拉远头(RRH)不具备缓存功能。然而下一代通信网络具有以内容为中心的特性,因此在Cloud-RAN中考虑带缓存的RRHs也变得有必要。该文考虑在Cloud-RAN中有效设计缓存方案,并通过资源分配有效减轻前程链路负担。假设系统采用正交频分多址接入(OFDMA)技术,通过联合优化子载波(SC)分配,RRH选择与传输功率,最小化系统下行总功耗,并通过拉格朗日对偶分解转化非凸问题,获得最优分配方案。仿真结果表明,比起其它缓存方案,该文提出的优化算法可以有效地提升系统能效,满足未来通信需求。
  • 云接入网络(Cloud Radio Access Network, Cloud-RAN),作为中国移动首先提出的技术,有希望是未来5G标准的候补[1]。在Cloud-RAN中,基带处理部分被聚集并且共享在一个虚拟的基带单元(Base Band Unit, BBU)池。作为软中继的射频拉远头(Remote Radio Heads, RRHs),可以通过有线/无线前程连接,将移动用户设备(User Equipment, UE)接收来的信号压缩并转发至中央BBU池。与分布式天线架构[2]和大规模天线系统[3]不同,在Cloud-RAN系统中,分布在RRHs与BBU池中间的前程链路容量是有限的。为了在前程容量约束的条件下提升Cloud-RAN性能,已有相关文献进行信号量化/压缩相关技术的研究[4,5]。但减轻前程链路负担的途径不止一条,通过设计带缓存功能RRHs,同样可以在前程链路容量受限的条件下,提升系统性能。尽管Cloud-RAN中的绿色通信技术也有一些相关的工作[6,7],但大多数已有的工作集中于研究RRHs不带缓存功能场景。在Cloud-RAN研究领域,在“RRHs带缓存”场景下,如何有效设计资源分配方案,提升系统性能,现有的工作十分有限。

    在现有带缓存的Cloud-RAN文献中,文献[8]首先从性能分析角度推导了无线信道的有效容量,并建立了嵌套联合形式博弈模型,进行资源的联合优化。文献[9]将稀疏多播波束成型技术引入带缓存的Cloud-RAN,建立了混合整数非线性规划问题,比较了不同传播、缓存方式的优劣。华为的美国研发团队与法国研发团队分别利用YouTube与优酷网的数据,从协作分层缓存[10]与主动缓存[11]两个不同的角度,设计了优化缓存方式,提升相关性能指标。文献[12]从更普遍的移动多媒体服务需求出发,研究了在线资源分配、内容放置与请求路径的联合优化问题。文献[13]的研究兴趣集中在文件传输,通过联合优化传输波束成型与文件传输策略,最小化用户下载文件的时延。文献[14]将文献[9]的工作进一步拓展,在多播前程链路中,研究了最优缓存大小的策略。文献[15]采用遗传算法作为优化工具,提出了权衡内容缓存策略以减少小区平均中断概率,减轻前程链路负担。

    在以上相关工作中,文献[8]与文献[9]最早发表在2016年,文献[11]发表在2017年,文献[10], 文献[1215]均是2018年的工作。本文工作延续了以上文献的研究思路,但希望解决已有工作的一些缺陷:(1)干扰管理技术过于复杂。大多数文献采用复杂的预编码波束成型技术解决干扰问题,为了有效降低过于复杂的干扰管理,并与现有的4G网络兼容,考虑正交频分多址接入(OFDMA)技术引入带缓存的Cloud-RAN;(2)未能考虑联合优化多种因素。以上文献大多侧重某一方面变量优化,未曾联合考虑“缓存状态”“用户最低需求”“RRHs选择”等多种基本要素。为了使优化问题更加实际,需要进行多种因素的联合优化;(3)RRHs的节能需求未能突出。在Cloud-RAN中,RRHs作为功能弱化的处理转发节点,其续航能力也受到限制。因此,RRHs的能耗需求也成为一个需要关注的问题,而以上的文献并没有特别关注这一方面。

    因此,本文采用OFDMA的基于缓存的Cloud-RAN中下行传输过程中,在考虑用户需求与RRHs缓存状态的情形下,通过联合优化用户子载波分配、RRH选择、子载波传输功率分配,最小化RRHs的传输功率。由于问题是非凸的,因此提出拉格朗日对偶算法,有效地转化问题为凸问题,并给出了优化的资源分配方案。仿真表明,比起其他缓存方案,本文提出的优化算法可以有效地提升系统能效,性能有明显提升。由于篇幅限制,暂不考虑有关RRH调度的传输时延问题。

    图1所示,考虑Cloud-RAN单个聚类中下行传输中,存在数量为M单天线的RRHs, 定义为M={1,2,···,M}K个单天线用户,定义为K={1,2,···,K};所有内容数目为T,定义为T={1,2,···,T},每个RRH的存储数量L不能超过上界T,即LT。当内容t存储到第m-th RRH时,定义im,t=1,否则im,t=0。对于用户k请求内容t,也有类似的定义指示θk,t=1,否则θk,t=0。设定云中心已知所有RRHs的缓存内容集合{im,t},与所有用户的请求内容集合{θk,t}。假定每个用户在同一时刻只能请求一份内容,因此满足Tt=1θk,t=1,tT,但同一份内容可以被多个用户请求,tkT为用户k所请求的内容。

    图 1  带缓存Cloud-RAN下行传输结构图

    在RRH中缓存的业务信息受欢迎的程度满足文献[15]中所提及的Zipf分布,并通过OFDMA方式,根据请求内容具体的大小,在单个或者多个间隔内传输。本文假定在单个间隔内完成资源分配优化。定义φk,n指示子载波SCn是否分配给用户k,当子载波SCn被分配给用户k时,φk,n=1,否则φk,n=0。同时定义φn[φ1,n,φ2,n,···,φK,n]T为在子载波SCn上的用户分配。根据OFDMA的定义,设定任何子载波SCn同时最多只能分配给单个用户,即满足1Tφn1,nN。分配给用户k的子载波集合定义为NkN,当φk,n=1时,Nk={n},并满足Nk1Nk2=, k1k2,k1, k2K

    前程链路容量受限时,RRH一般仅会向中央云请求用户的未缓存内容,并通过OFDMA下行传输给用户。因此定义RRH调用子载波传输的指示因子为χm,n,当RRHm调用子载波n进行传输时,χm,n=1,否则χm,n=0。定义通过子载波n传输的RRHs子集合为Xn={mM|χm,n=1}, nN。因此,Xn中的RRHs可以在子载波n上协作传输内容给用户k,即φk,n=1

    定义hk,m,n为从RRHm到用户k在子载波n上的无线信道增益系数,而云端已知所有信道状态信息。在Xn中所有RRHs同时传输,参考文献[8],分配到子载波n的用户k上的接收信噪比(SNR)可以表示为

    γk,n(χn,pn)=|hTk,nxn|2/σ2=(Mm=1|hk,m,n|χm,npm,n)2/σ2 (1)

    其中,xn代表SCn向用户k传输信号的矢量,σ2是接收端加性高斯白噪声的功率,所有用户端噪声相同。在SCnNk上的可达速率可表示为

    Rk,n(χn,pn)=Blog2(1+γk,n(χn,pn))/N (2)

    χn固定时,Rk,n(χn,pn)对变量{pm,n}是凹函数,由于篇幅限制,证明略。

    图1列举了几种用户请求内容的典型模式,具体如下所示:

    (1) User1模式:User1请求的内容Data1(D1)仅在某个RRH的缓存中,比如RRH1。此时不需要通过前程链路从BBU池调度数据,由RRH1直接下行传输内容。

    (2) User2模式:User2请求的内容Data2(D2)同时存储在多个RRH的缓存中,比如RRH2, RRH3。此时不需要通过前程链路从BBU池调度数据,由RRH2与RRH3联合协同下行传输内容。

    (3) User3模式:User3请求的内容Data3(D3)只有部分存储在多个RRH的缓存中,比如RRH2, RRH4。需要由离User3最近的RRH3通过前程链路从BBU池调度缺失数据,再由RRH2, RRH3, RRH4联合协同下行传输内容。在多个RRHs协同发送数据给User3时,可以参照现有文献[1618]设计一种联合的时频同步算法,保证传输信息不出现乱序。本文默认传输同步已经完成。

    (4) User4模式:User4请求的内容Data4(D4)没有存储在任何RRH的缓存中。需要由离User4最近的RRH4通过前程链路从BBU池调度完整数据,由RRH4直接下行传输内容。

    在实际中,由于多个用户可能都请求最受欢迎的某内容,RRHs将受欢迎的内容按照比例放在缓存中,并按照实际请求从BBU池调度一些独特的内容。无论如何,由于前程链路容量是受限的。因此RRHm通过OFDMA方式,向 N SCs上所有用户传输非缓存数据的速率(即从BBU池调度下行的数据),不应当超过前程链路容量Ftlm(tl是total的简写),即需满足限制条件

    Tt=1(1im,t)max (3)

    基于第2节系统模型的定义,本节将建立具体的优化模型。优化目标为最小化所有RRHs在SCs上的传输功率,优化变量为用户-子载波分配{\left\{ {{{\text{φ}}_n}} \right\}_{n \in \cal{N}}}, RRHs选择{\left\{ {{{\text{χ}}_n}} \right\}_{n \in \cal{N}}},功率分配{\left\{ {{{\text{p}}_n}} \right\}_{n \in \cal{N}}},约束条件为前程链路容量限制与用户最小速率限制。优化问题(P1)为

    {\rm P}1\,\mathop {\min }\limits_{{{\left\{ {{{\text{p}}_n},{{\text{χ}}_n},{{\text{φ}}_n}} \right\}}_{n \in {\cal N}}}} P_{{\rm{RRHs}}}^{{\rm{tl}}} \hspace{140pt} (4)

    s.t.

    \begin{aligned} {\rm{}}& \sum\limits_{t = 1}^{{T}} {\left( {1 - {i_{m,t}}} \right)} \mathop {\max }\limits_{k \in \cal{K}}\Bigg\{ {\theta _{k,t}}\sum\limits_{n = 1}^N {\chi _{m,n}}{\varphi _{k,n}}\\ {\rm{}}&\quad \cdot{R_{k,n}}\left( {{{\text{χ}}_n},{{\text{p}}_n}} \right) \Bigg\} \le F_m^{{\rm{tl}}} \end{aligned} \hspace{50pt} (5)
    \sum\limits_{n = 1}^N {{\varphi _{k,n}}} \frac{B}{N}{\log _2}\left( {1 \!+\! {\gamma _{k,n}}\left( {{{\text{χ}}_n},{{\text{p}}_n}} \right)} \right) \ge R_k^{\min },\forall k \in \cal{K} \hspace{10pt} (6)
    {\chi _{m,n}} \in \left\{ {0,1} \right\},\forall m \in \cal{M},\forall n \in \cal{N} \hspace{70pt} (7)
    {\text{1}^{\rm{T}}}{{\text{φ}}_n} \le 1,{\varphi _{k,n}} \in \left\{ {0,1} \right\},\forall k \in {\cal K},n \in {\cal N} \hspace{40pt} (8)
    P_{{\rm{RRHs}}}^{{\rm{tl}}} = \sum\limits_{m = 1}^M {\sum\limits_{n = 1}^N {{p_{m,n}}} } + {P_c},\forall m \in \cal{M},n \in \cal{N} \hspace{35pt} (9)

    {P_c}为传输固定功率损耗。需要特别指出的是,式(9)中的变量{p_{m,n}}的物理意义是{\rm RRH}m在子载波n上所分配的功率,是一个非负值,而式(10)中的优化变量{{\text{p}}_n} = {\left[ {{p_{1,n}},{p_{2,n}} ·\!·\!·\,{p_{M,n}}} \right]^{\rm{T}}} \in \mathbb{R}_ + ^{M \times 1}是一个矩阵。

    式(5)是第2节所提及的前程链路容量约束,具体可见式(3);式(6)是满足用户需求的最小速率,即传统网络所提及的QoS(Quality of Service)限制;式(7)是{\rm RRH}m是否在{\rm SC}n传输的指示因子,{\chi _{m,n}} = 1代表{\rm RRH}m选择在{\rm SC}n上进行传输;式(8)是OFDMA设定的SC分配方式,{\varphi _{k,n}} = 1代表载波{\rm{SC}}n被分配给用户k, {\text{1}^{\rm{T}}}{{\text{φ}}_n} \le 1代表任何子载波{\rm{SC}}n同时最多只能分配给一个用户。

    为了简化约束条件式(5),考虑引入辅助变量

    {\xi _{m,t}} = {{\mathop {\max }\limits_{k \in \cal{K}} \!\left\{\! {{\theta _{k,t}}\sum\limits_{n = 1}^N {{\chi _{m,n}}{\varphi _{k,n}}{R_{k,n}}\left( {{{\text{χ}}_n},{{\text{p}}_n}} \right)} } \!\right\}}\!\! \Bigg/ \!\!{F_m^{{\rm{tl}}}}} \hspace{18pt} (10)

    优化问题(P1)可以等效转化为问题(P2)

    {\rm P}2 \mathop {\min }\limits_{{{\left\{ {{{\text{p}}_n},{{\text{χ}}_n},{{\text{φ}}_n}} \right\}}_{n \in {\cal N}}},\left\{ {{\xi _{m,t}}} \right\}} P_{{\rm{RRHs}}}^{{\rm{tl}}} \hspace{60pt} (11)

    s.t.

    \begin{aligned} {\rm{}}& \frac{1}{{F_m^{{\rm{tl}}}}}{\theta _{k,t}}\left( {1 \!-\! {i_{m,t}}} \right)\!\!\sum\limits_{n = 1}^N {\chi _{m,n}}{\varphi _{k,n}}{R_{k,n}}\left( {{{\text{χ}}_n},{{\text{p}}_n}} \right)\\ {\rm{}}& \quad \le {\theta _{k,t}}\left( {1 - {i_{m,t}}} \right){\xi _{m,t}} \end{aligned} (12)
    \sum\limits_{t = 1}^{{T}} {\left( {1 - {i_{m,t}}} \right){\xi _{m,t}}} \le 1,\forall m \in \cal{M} \hspace{30pt} (13)

    式(6),式(7),式(8),式(9)

    当RRHs已经缓存用户请求的内容,即{i_{m,t}} = 1时,约束条件式(12), 式(13)无意义。由于式(2)中定义的{R_{k,n}}\left( {{{\text{χ}}_n},{{\text{p}}_n}} \right)对于\left\{ {{p_{m,n}}} \right\}是凹函数,且(P2)中约束条件是非凸的,可以得知问题(P2)整体上是非凸的。

    文献[19]分析了非凸优化问题的对偶间隙,当目标与约束在SCs上可分离,SCs数目趋向无限时,对偶间隙趋向至0。在优化问题(P2)中,目标式(11)与约束式(12)在SCs上可分离,约束式(13)在SCs上虽然不可分离,但不影响整体问题(P2)的凹凸性。因此,当N趋向于无穷时,问题(P2)的对偶界限趋向于0,因此可以利用拉格朗日对偶方法[20]解决问题(P2)。

    设定约束式(12)的对偶因子为{\alpha _{m,k,t}} \ge 0, m \in \cal{M}, k \in \cal{K}, t \in \cal{T}\;;约束式(6)的对偶因子为{\beta _k} \ge 0, k \in \cal{K}。如果将\cal{M}_{{t_k}}^i \triangleq \left\{ {m|{i_{m,{t_k}}} = 0} \right\}, k \in \cal{K}定义为缓存中没有被用户k请求内容{t_k}的RRHs合集,则在约束式(12)中存在I \triangleq \displaystyle\sum\nolimits_{k = 1}^K {\left| {\cal{M}_{{t_k}}^I} \right|}

    定义向量{\text{α}} \in \mathbb{R}_ + ^{I \times 1}, {\text{β}} \in \mathbb{R}_ + ^{K \times 1}为对偶变量的合集,仅有约束式(12)与约束式(6)的问题(P2)的拉格朗日变换为

    \begin{aligned} {\rm{}}& L\left( {{{\left\{ {{{\text{p}}_n},{{\text{χ}}_n},{{\text{φ}}_n}} \right\}}_{n \in {\cal N}}},\left\{ {{\xi _{m,t}}} \right\},{\text{α}},{\text{β}}} \right) \\ {\rm{}}& \quad= \sum\limits_{n = 1}^N {{L_n}\left( {{{\text{p}}_n},{{\text{χ}}_n},{{\text{φ}}_n},{\text{α}},{\text{β}}} \right)} - \sum\limits_{m = 1}^M \sum\limits_{k = 1}^K \sum\limits_{t = 1}^T {\alpha _{m,k,t}} \\ {\rm{}}&\quad\quad \cdot{\theta _{k,t}}\left( {1 - {i_{m,t}}} \right){\xi _{m,t}} + \sum\limits_{k = 1}^K {{\beta _k}} \end{aligned} (14)

    其中,子函数{L_n}\left( {{{\text{p}}_n},{{\text{χ}}_n},{{\text{φ}}_n},{\text{α}},{\text{β}}} \right)定义为

    \begin{aligned} {\rm{}}& {L_n}\left( {{{\text{p}}_n},{{\text{χ}}_n},{{\text{φ}}_n},{\text{α}},{\text{β}}} \right) \\ {\rm{}}& \triangleq \sum\limits_{m = 1}^M {{p_{m,n}}} + \sum\limits_{m = 1}^M {\sum\limits_{k = 1}^K {\sum\limits_{t = 1}^T {{\chi _{m,n}}} } } {\varphi _{k,n}}{\theta _{k,t}}\left(\! {1 \!-\! {i_{m,t}}} \!\right)\!\frac{{{\alpha _{m,k,t}}}}{{F_m^{{\rm{tl}}}}}\\ {\rm{}}& \quad\cdot{R_{k,n}}\!\left( {{{\text{χ}}_n},{{\text{p}}_n}} \right)- \sum\limits_{k = 1}^K {\frac{{{\beta _k}}}{{R_k^{\min }}}{\varphi _{k,n}}{R_{k,n}}\left( {{{\text{χ}}_n},{{\text{p}}_n}} \right)} \end{aligned} (15)

    对应的拉格朗日对偶函数可以表示为

    \begin{aligned} {\rm{}}& g\left( {{\text{α}},{\text{β}}} \right) =\mathop {\min }\limits_{{{\left\{ {{{\text{p}}_n},{{\text{χ}}_n},{{\text{φ}}_n}} \right\}}_{n \in {\cal N}}},\left\{ {{\xi _{m,t}}} \right\}} L\left( {{\left\{ {{{\text{p}}_n},{{\text{χ}}_n},{{\text{φ}}_n}} \right\}}_{n \in {\cal N}}},\right.\\ {\rm{}}& \qquad\qquad\ \left\{ {{\xi _{m,t}}} \right\},{\text{α}},{\text{β}} \big),\\ {\rm{}}& {\rm s.t.}\;\; \text{式}(13), \text{式}(7)-\text{式}(9) \end{aligned} (16)

    参照约束式(12),约束式(6)的对偶函数可以表示为

    g\left( {{\text{α}},{\text{β}}} \right) = {g_1}\left( {{\text{α}},{\text{β}}} \right) + {g_2}\left( {{\text{α}},{\text{β}}} \right) + \sum\limits_{k = 1}^K {{\beta _k}} \hspace{56pt} (17)
    \begin{aligned} {\rm{}}&\! {g_1}\left( {{\text{α}},{\text{β}}} \right) = \mathop {\min }\limits_{{{\left\{ {{{\text{p}}_n},{{\text{χ}}_n},{{\text{φ}}_n}} \right\}}_{n \in {\cal N}}}} \sum\limits_{n = 1}^N {{L_n}} \left( {{{\text{p}}_n},{{\text{χ}}_n},{{\text{φ}}_n},{\text{α}},{\text{β}}} \right),\hspace{23pt}\\ {\rm{}}& \!{\rm s.t.} \;\;\text{式}(7)-\text{式}(9) \end{aligned} (18)
    \begin{aligned} {\rm{}}&\! {g_2}\left( {{\text{α}},{\text{β}}} \right) = \mathop {\min }\limits_{\left\{ {{\xi _{m,t}}} \right\}} - \sum\limits_{m = 1}^M {\sum\limits_{t = 1}^T {\sum\limits_{k = 1}^K {{\theta _{k,t}}\left( {1 \!-\! {i_{m,t}}} \right){\alpha _{m,k,t}}{\xi _{m,t}}} } }, \\ {\rm{}}& \! {\rm s.t.} \;\;\text{式}(13) \end{aligned} (19)

    式(22)中的最小化问题可以分解成N个平行的子问题,每个子问题对应一个单独的{\rm SC}n \in \cal{N}。所有子问题可以用相似的优化问题表示。

    \left.{\begin{aligned} {\rm{}}& {\mathop {\min }\limits_{{{\left\{ {{{\text{p}}_n},{{\text{χ}}_n},{{\text{φ}}_n}} \right\}}_{n \in {\cal N}}}} \sum\limits_{n = 1}^N {{L_n}} \left( {{{\text{p}}_n},{{\text{χ}}_n},{{\text{φ}}_n},{\text{α}},{\text{β}}} \right)}\\ {\rm{}}& {\rm s.t.} \;\;{{\text{p}}_n} \succeq {\text{0}},{{\text{χ}}_n} \in {{\left\{ {0,1} \right\}}^{M \times 1}},\\ {\rm{}}&\quad\ \ {{\text{1}}^{\rm{T}}}{{\text{φ}}_n} \le 1,{{\text{φ}}_n} \in {{\left\{ {0,1} \right\}}^{K \times 1}} \end{aligned}}\right\} (20)

    式(20)中,{L_n}\left( {{{\text{p}}_n},{{\text{χ}}_n},{{\text{φ}}_n},{\text{α}},{\text{β}}} \right)与式(15)中的定义一致。下面的内容将研究如何解决优化问题式(20)。

    定义{\rm SC}n上用户关联策略固定{\varphi _n} = {{\mathop \varphi \limits^{\smile}} _n},当{{\mathop \varphi \limits^{\smile}} _n} = 0时,{\rm SC}n不被任何用户占用。由于{{\text{p}}_n} \succeq {\text{0}},式(20)中的目标函数仅在{{\text{p}}_n} = {\text{0}}时实现最小值,与RRHs选择变量{{\text{χ}}_n}无关。因此,如果没有用户被分配{\rm SC}n,所有RRHs的功率分配{{\text{p}}_n} = {\text{0}}时,可以不失一般性地假设{{\text{χ}}_n} = {\text{0}}。定义{{\mathop k \limits^{\smile}} _n} \in \cal{K}是分配{\rm SC}n的用户,因此{{\mathop \varphi \limits^{\smile}} _{{{{\mathop k \limits^{\smile}} }_n},n}} = 1, {{\mathop \varphi \limits^{\smile}} _{k,n}} = 0;定义{t_{{{{\mathop k \limits^{\smile}} }_n}}} \in \cal{T}\;是用户{{\mathop k \limits^{\smile}} _n}请求内容,满足{\theta _{{{{\mathop k \limits^{\smile}} }_n},{t_{{{{\mathop k \limits^{\smile}} }\!_n}}}}} = 1, {\theta _{{{{\mathop k \limits^{\smile}} }_n},t}} = 0, \forall t \ne {t_{{{{\mathop k \limits^{\smile}} }\!\!_n}}}。式(20)可以转化为

    {\begin{aligned} {\rm{}} & \! \begin{array}{*{20}{c}} {\mathop {\min }\limits_{{p_n},{\chi _n}} }{ -\! \left(\! {{{{\beta _{{{{\mathop k \limits^{\smile}} }_n}}}} / {R_{{{{\mathop k \limits^{\smile}} }_n}}^{\min }}} \!-\! \!\sum\limits_{m = 1}^M \!\!{\left(\!\! {1 \!-\! {i_{m,{{{\mathop k \limits^{\smile}} }_n}}}} \right)\!{\chi _{m,n}}{{{\alpha _{m,{{{\mathop k \limits^{\smile}} }_n},{t_{{{{\mathop k \limits^{\smile}} }\!\!_n}}}}}} \Bigg/\!\! {F_m^{{\rm{tl}}}}}} } \!\!\right)} \end{array}\\ {\rm{}}& \qquad\cdot {R_{{{{\mathop k \limits^{\smile}} }_n},n}}\left( {{{\text{χ}}_n},{{\text{p}}_n}} \right) + \sum\limits_{m = 1}^M {{p_{m,n}}} , \\ {\rm{}} & \;{\rm s.t.}\;\;{{{\text{p}}_n} \succeq {\text{0}},{{\text{χ}}_n} \in {{\left\{ {0,1} \right\}}^{M \times 1}}} \end{aligned}} (21)

    对于给定的RRH选择方案{{\mathop {\text{χ}} \limits^{\smile}} _n},问题式(21)的最优功率分配方案{{\mathop {\text{p}} \limits^{\frown}} _n}由以下推论给出。

    推论1 假定{{\text{χ}}_n} = {{\mathop {\text{χ}} \limits^{\smile}} _n},问题式(21)的最优功率分配方案{{\mathop {\text{p}} \limits^{\frown}} _n}如式(26)

    {{\mathop {p} \limits^{\frown}} _{m,n}} \!= \!{{\Pi \left( {{{{\mathop {\text{χ}} \limits^{\smile}} }_n}} \right){{\!\left| {{h_{{{{\mathop k \limits^{\smile}} }\!_n},m,n}}}\!\right|}^2}{{{\mathop {\chi } \limits^{\smile}} }_{m,n}}}\!\! \bigg/\bigg [\!\!{{{\left(\! {{\varTheta _{_{{{{\mathop k \limits^{\smile}} }\!\!_n},n}}}\!\left(\! {{{{\mathop {\text{χ}} \limits^{\smile}} }\!_n}} \right)} \!\right)}^2}\!{\sigma ^2}}}\!\bigg] \hspace{13pt} (22)

    其中,

    \Pi \!\left( {{{{\mathop {\text{χ}} \limits^{\smile}} }_n}} \right) \!=\! {\left[ {\frac{{B{\varOmega _{_{{{{\mathop k \limits^{\smile}} }\!\!_n},n}}}\left( {{{{\mathop {\text{χ}} \limits^{\smile}} }_n}} \right){\varTheta _{_{{k_n},n}}}\!\left(\! {{{{\mathop {\text{χ}} \limits^{\smile}} }_n}} \right)}}{{N\ln 2}} \!-\! 1} \right]^ + }\hspace{60pt} (22a)
    {\varOmega _{{{{\mathop k \limits^{\smile}} }_n},n}}\!\left(\! {{{{\mathop {\text{χ}} \limits^{\smile}} }_n}} \!\right)\! =\! \frac{{{\beta _{{{{\mathop k \limits^{\smile}} }_n}}}}}{{R_{{{{\mathop k \limits^{\smile}} }_n}}^{\min }}}\! -\!\! \displaystyle\sum\limits_{m = 1}^M \!{\frac{{\left(\! {1 \!-\! {i_{m,{t_{{{{\mathop k \limits^{\smile}} }_n}}}}}}\! \right){\chi _{m,n}}{\alpha _{m,{{{\mathop k \limits^{\smile}} }_n},{t_{{{{\mathop k \limits^{\smile}} }_n}}}}}}}{{F_m^{{\rm{tl}}}}}} \hspace{24pt} (22b)
    {\varTheta _{_{{{{\mathop k \limits^{\smile}} }\!\!_n},n}}}\left( {{{{\mathop {\text{χ}} \limits^{\smile}} }_n}} \right)= \displaystyle\sum\limits_{m = 1}^M {\frac{{{\chi _{m,n}}{{\left| {{h_{{{{\mathop k \limits^{\smile}} }_n},m,n}}} \right|}^2}}}{{{\sigma ^2}}}} \hspace{88pt} (22c)

    推论1说明,对于给定用户关联策略{{\mathop k \limits^{\smile}} _n}与RRH选择策略{{\mathop {\text{χ}} \limits^{\smile}} _n},最优功率分配策略是阈值结构。当{\varOmega _{{{{\mathop k \limits^{\smile}} }\!_n},n}}\left( {{{{\mathop {\text{χ}} \limits^{\smile}} }_n}} \right){\varTheta _{_{{{{\mathop k \limits^{\smile}} }\!_n},n}}}\left( {{{{\mathop {\text{χ}} \limits^{\smile}} }_n}} \right) \le {{\left( {N\ln 2} \right)} / B},所有RRHs在SC n上分配零功率;当{\varOmega _{{{{\mathop k \limits^{\smile}} }\!_n},n}}\left( {{{{\mathop {\text{χ}} \limits^{\smile}} }_n}} \right){\varTheta _{_{{{{\mathop k \limits^{\smile}} }\!_n},n}}}\left( {{{{\mathop {\text{χ}} \limits^{\smile}} }_n}} \right) > {{\left( {N\ln 2} \right)} / B}, {\rm RRH }m上分配的功率取决于无线信道增益\left| {{h_{{{{\mathop k \limits^{\smile}} }_n},m,n}}} \right|, {\beta _{{{{\mathop k \limits^{\smile}} }_n}}}等变量。

    在式(22)中,如果所有RRHs均缓存用户请求的内容{t_{{{{\mathop k \limits^{\smile}} }_n}}},即{i_{m,{t_{{{{\mathop k \limits^{\smile}} }\!\!_n}}}}} = 1,功率分配表达式可简化为

    \begin{aligned} {{\mathop {p} \limits^{\frown}} _{m,n}} =& {{{{\mathop {\text{χ}} \limits^{\smile}} }_{m,n}}{{\left| {{h_{{{{\mathop k \limits^{\smile}} }_\!n},m,n}}} \right|}^2}{{\!\!\left[ \!{{{B{\beta _{{{{\mathop k \limits^{\smile}} }\!_n}}}{\varTheta _{_{{{{\mathop k \limits^{\smile}} }\!\!_n},n}}}\!\!\left(\! {{{{\mathop {\text{χ}} \limits^{\smile}} }_n}}\! \right)} \!\!/\! {R_{{{{\mathop k \limits^{\smile}} }_n}}^{\min }\!N\ln 2 \!-\! 1}}} \!\right]}^ + }} \\ {\rm{}}& \bigg/\left[{{{{\left( {{\varTheta _{{{{\mathop k \limits^{\smile}} }\!\!_n},n}}\left( {{{{\mathop {\text{χ}} \limits^{\smile}} }_n}} \right)} \right)}^2}{\sigma ^2}} }\right] \end{aligned} (23)

    如果{{\mathop {\text{χ}} \limits^{\smile}} _n} = {\text{0}},即没有RRH选择,{{\mathop {\text{p}} \limits^{\frown}} _n} = {\text{0}}。对于聚类仅有一个RRH的特例,功率分配表达式可简化为

    \begin{aligned} {{\mathop {p} \limits^{\frown}}\! _n} =& \left[ {B\left( {{{{\beta _{{{{\mathop k \limits^{\smile}} }\!_n}}}} / {R_{{{{\mathop k \limits^{\smile}} }\!\!_n}}^{\min } - {{\left( {1 - {i_{{t_{{{{\mathop k \limits^{\smile}} }\!\!_n}}}}}} \right){\alpha _{{{{\mathop k \limits^{\smile}} }\!_n},{t_{{{{\mathop k \limits^{\smile}} }\!\!_n}}}}}} / {{F^{\rm tl}}}}}}} \right)} \right.\\ {\rm{}}& \cdot\left. { {(N\ln 2)} - {{{\sigma ^2}}\Big / {{{\left| {{h_{_{{{{\mathop k \limits^{\smile}} }\!\!_n}},n}}} \right|}^2}}}} \right]^ + \end{aligned} (24)

    式(24)与已知的注水算法形式相似,主要区别在于在不同SCn \in \cal{N}上的水位不同。如果\left( {{{{\beta _{{{{\mathop k \limits^{\smile}} }\!\!_n}}}} / {R_{{{{\mathop k \limits^{\smile}} }\!\!_n}}^{\min }}}} \right) \le {{\left( {1 - {i_{{t_{{{{\mathop k \limits^{\smile}} }\!_n}}}}}} \right){\alpha _{{{{\mathop k \limits^{\smile}} }\!\!_n},{t_{{{{\mathop k \limits^{\smile}} }\!\!\!_n}}}}}} / {{F^{\rm tl}}}}, {\rm{SC}}n上不分配任何功率。

    对于任何给定对偶变量{\text{α}}, {\text{β}},式(21)可以按以下步骤,通过推论1得到最优解

    步骤 1 对{{\mathop k \limits^{\smile}} _n} \in \cal{K},固定SC n上的用户;

    步骤 2 对于{2^M}种RRH选择可能,利用式(22)计算最优功率分配策略{{\mathop {\text{p}} \limits^{\frown}} _n},并将式(22)代入式(21),通过使式(21)目标函数最小化,得到最优RRH选择策略{{\mathop {\text{χ}} \limits^{\frown}} _n}

    步骤 3 基于最优功率分配策略{{\mathop {\text{p}} \limits^{\frown}} _n},最优RRH选择策略{{\mathop {\text{χ}} \limits^{\frown}} _n},通过使式(20)目标函数最小化,得到最优用户SC分配策略{{\mathop {\text{φ}} \limits^{\frown}} _n}

    推论2 问题式(19)的最优解如式(25)所示

    {{\mathop {\xi} \limits^{\frown}} _{m,t}} = \!\!\left\{\!\! {\begin{aligned} {\rm{}}& {1,}\; t = \mathop {\arg \max }\limits_{l \in {\cal T}}\Bigg\{ \sum\limits_{k = 1}^K{\theta _{k,t}} \\ {\rm{}}& \quad\cdot {\left(\! {1 \!-\! {i_{m,l}}} \right){\alpha _{m,k,l}}} \!\Bigg\}\\ {\rm{}}& {0,}\; {{\text{其它}}} \end{aligned}} \right. (25)

    推论2说明式(19)可以通过搜寻\displaystyle\sum\nolimits_{k = 1}^K {{\theta _{k,t}}\left( {1 - {i_{m,l}}} \right){\alpha _{m,k,l}}} 的最大值内容,并将其与用户请求的内容比较,如果两者一致,则{{\mathop {\xi} \limits^{\frown}} _{m,t}} = 1。通过分析,可知寻找式(19)最优解的最大复杂度为\cal{O}\left( {MT} \,\right),发生在至少T\,个用户请求不同内容,且所有内容都没有在任何RRHs端缓存。

    此时,问题(P1)的对偶问题可表示为

    \mathop {\max }\limits_{\alpha \succeq {\text{0}},\beta \succeq {\text{0}}} g\left( {{\text{α}},{\text{β}}} \right) (26)

    式(30)是一个凸问题,可以通过经典凸优化方法,如椭球法得到最优对偶变量{\mathop {\text{α}} \limits^{\frown}} , {\mathop {\text{β}} \limits^{\frown}} 。基于推论1与推论2得到的最优解{{\mathop {\text{p}} \limits^{\frown}} _n}, {{\mathop {\text{χ}} \limits^{\frown}} _n}, {{\mathop {\text{φ}} \limits^{\frown}} _n}, {{\mathop {\xi} \limits^{\frown}} _{m,t}},与最优对偶变量{\mathop {\text{α}} \limits^{\frown}} ,{\mathop {\text{β}} \limits^{\frown}} ,本文提出解决式(4)中问题(P1)基于功耗最小的联合优化算法步骤如下:

    步骤 1 相关参数初始化,设置{\text{α}}, {\text{β}},最大迭代次数{I_{\max }},收敛条件\varepsilon 等;

    步骤 2 给定分配{\rm SC}n的用户{{\mathop k \limits^{\smile}} _n}, RRHs选择SCs策略{{\mathop {\text{χ}} \limits^{\frown}} _n},基于推论1得出最优功率分配{{\mathop {\text{p}} \limits^{\frown}} _n}

    步骤 3 将{{\mathop {\text{p}} \limits^{\frown}} _n}代入式(21), 通过使式(21)目标函数最小化,得到最优RRH选择SCs策略{{\mathop {\text{χ}} \limits^{\frown}} _n}

    步骤 4 将{{\mathop {\text{p}} \limits^{\frown}} _n}, {{\mathop {\text{χ}} \limits^{\frown}} _n},通过使式(20)目标函数最小化,得到最优用户SC分配策略{{\mathop {\text{φ}} \limits^{\frown}} _n}

    步骤 5 基于推论2,解决式(19)中的优化问题,得到最优解{{\mathop {\xi} \limits^{\frown}} _{m,t}}

    步骤 6 利用椭球法,更新对偶变量{\mathop {\text{α}} \limits^{\frown}} ,{\mathop {\text{β}} \limits^{\frown}}

    步骤 7 重复步骤2—步骤6,直至算法收敛。

    通过分析,可知解决式(4)中优化问题(P1)的计算复杂度是\cal{O}\left( {NK{2^M}} \right);式(20)的计算复杂度是\cal{O}\left( {K{2^M}} \right);式(16)的计算复杂度最小值\cal{O}\left( M \right),最大值\cal{O}\left( {MT} \;\right);算法1的计算复杂度最小值为\cal{O}\left( {2NK + 1} \right),最大值\; \cal{O}\left( {NK{2^M} + MT} \;\right)

    仿真中,设定一个小区聚类中无线带宽为B=20 MHz,均分为32个SCs。假设用户所经历的是独立同分布6路径瑞利衰落信道,背景噪声的功率谱密度设为–169 dBm/Hz。RRHs个数M = 4,在半径为100 m范围内均匀分布,用户个数K = 8,在半径为100 m范围内随机分布。满足用户最低传输需求R_k^{\min }=20 Mbps,用户请求所有内容数目T = 40,满足Zipf分布[9],用户请求内容t的概率为{\rho _t} = {t^{ - \lambda }}/\displaystyle\sum\nolimits_{t = 1}^{{T}} {{t^{ - \lambda }}} ,t \in \cal{T}\;, \lambda 为分布参数,设定\lambda = 0.9。每个RRH前程容量约束F_m^{{\rm{tl}}}=60 Mbps,最大缓存内容S = 4

    在仿真中,除了本文提出的算法,还涉及以下对比方案:

    方案1  最大概率缓存。每个RRH根据内容的流行概率进行缓存,直到缓存空间全部占用。

    方案2  无缓存。RRHs没有缓存功能,用户请求任何数据,均需要通过前程链路向BBU池请求。

    方案3  平均功率分配。RRHs端在所有子载波上平均分配功率,不进行任何优化。

    方案4  单RRH选择。聚类中仅有1个RRH,不存在RRH选择的问题。

    图2给出了随着RRH总数变大,不同参数下本文算法的最大复杂度变化趋势。当SCs总数N与需求内容总数T\;一致时,用户数目K越大,算法最大复杂度越大。同样的趋势,对参数SCs总数N与需求内容总数T\;也同时成立。但必须指出,RRHs个数M对复杂度影响最大,因为运算中M位于指数位置,其它变量均为线性乘积。

    图 2  RRH个数变化,本文算法最大复杂度变化趋势

    图3给出了随着用户总数变大,不同参数下本文算法的最大复杂度变化趋势。当用户数目K线性增长,本文算法最大复杂度也线性增长。在图2中,当RRHs个数M线性增长,本文算法最大复杂度呈指数增长。这说明不同参数对复杂度影响不同,RRHs个数M占主导作用。

    图 3  用户个数变化,本文算法最大复杂度变化趋势

    图4给出了随着前程容量门限变化时,多种方案功耗对比。比起方案1—方案4,本文算法可以使单个RRH平均功率损耗最低,即保持最优性能。本文算法性能优于基于最大概率缓存的方案1,因为在大多数场景存在多个RRHs各有用户请求一部分内容,并联合发送给用户的状况,方案1更适应RRHs独立为用户缓存内容的场景。当前程容量门限增大时,本文算法的优势逐渐变小。F_m^{{\rm{tl}}}=90 Mbps时,除了不做任何缓存的方案2,本文算法与方案1, 3, 4性能几乎一致。因为随着前程容量门限增大,说明前程链路的负担已经被大大减轻,用户可以直接从BBU池请求发送下行数据而不必顾虑前程约束,RRHs缓存作用则受到削弱。这说明,缓存优化在前程链路严重受限场景更有必要。

    图 4  前程容量门限变化,多种方案功耗对比

    图5给出了随着单个RRH最大缓存数目变化时,多种方案功耗对比。与图4结论类似,本文算法可以一直实现单个RRH最低功率损耗,但算法优势随着单个RRH最大缓存数目变大而变小。因为当单个RRH最大缓存数目变大时,用户更倾向于从更少数目的RRHs请求内容,避免了向所有RRHs或者BBU池发送内容请求,从整体上节约能耗。这说明,RRHs增加缓存功能具有必要性,在节约能耗上有着非常明显的作用。

    图 5  单个RRH最大缓存数目变化,多种方案功耗对比

    本文针对带缓存的Cloud-RAN,以最小化RRHs能耗为目标,对SCs分配、RRHs选择与RRHs下行功率等变量进行了优化设计。通过拉格朗日对偶分解,原始的非凸问题转化为几个可以解决的独立子问题,基于子问题的最优解,本文提出了联合设计优化算法。仿真结果验证了所提算法可以有效地降低RRHs功耗,提升系统能效,符合未来绿色通信的需求。本文并没有涉及RRHs端前程容量与缓存大小动态分配的问题,在未来工作中需要进一步研究。

  • 图  1  带缓存Cloud-RAN下行传输结构图

    图  2  RRH个数变化,本文算法最大复杂度变化趋势

    图  3  用户个数变化,本文算法最大复杂度变化趋势

    图  4  前程容量门限变化,多种方案功耗对比

    图  5  单个RRH最大缓存数目变化,多种方案功耗对比

  • CHECKO A, CHRISTIANSEN H L, YAN Yan, et al. Cloud RAN for mobile networks—A technology overview[J]. IEEE Communications Surveys & Tutorials, 2015, 17(1): 405–426. doi: 10.1109/COMST.2014.2355255
    WANG Dongming, WANG Jiangzhou, YOU Xiaohu, et al. Spectral efficiency of distributed MIMO systems[J]. IEEE Journal on Selected Areas in Communications, 2013, 31(10): 2112–2127. doi: 10.1109/JSAC.2013.131012
    WANG Yi, LI Chunguo, HANG Yongming, et al. Energy-efficient optimization for downlink massive MIMO FDD systems with transmit-side channel correlation[J]. IEEE Transactions on Vehicular Technology, 2016, 65(9): 7228–7243. doi: 10.1109/TVT.2015.2483519
    LIU Liang, BI Suzhi, and ZHANG Rui. Joint power control and fronthaul rate allocation for throughput maximization in OFDMA-based cloud radio access network[J]. IEEE Transactions on Communications, 2015, 63(11): 4097–4110. doi: 10.1109/TCOMM.2015.2469687
    ZHOU Yuhan and YU Wei. Optimized backhaul compression for uplink cloud radio access network[J]. IEEE Journal on Selected Areas in Communications, 2014, 32(6): 1295–1307. doi: 10.1109/JSAC.2014.2328133
    LI Chunguo, SONG Kang, WANG Dongming, et al. Optimal remote radio head selection for cloud radio access networks[J]. Science China Information Sciences, 2016, 59(10): 102315. doi: 10.1007/s11432-016-0060-y
    LUO Shixin, ZHANG Rui, and LIM T J. Downlink and uplink energy minimization through user association and beamforming in C-RAN[J]. IEEE Transactions on Wireless Communications, 2015, 14(1): 494–508. doi: 10.1109/TWC.2014.2352619
    ZHAO Zhongyuan, PENG Mugen, DING Zhiguo, et al. Cluster content caching: An energy-efficient approach to improve quality of service in cloud radio access networks[J]. IEEE Journal on Selected Areas in Communications, 2016, 34(5): 1207–1221. doi: 10.1109/JSAC.2016.2545384
    TAO Meixia, CHEN Erkai, ZHOU Hao, et al. Content-centric sparse multicast beamforming for cache-enabled Cloud RAN[J]. IEEE Transactions on Wireless Communications, 2016, 15(9): 6118–6131. doi: 10.1109/TWC.2016.2578922
    TRAN T X, LE D V, YUE Guosen, et al. Cooperative hierarchical caching and request scheduling in a cloud radio access network[J]. IEEE Transactions on Mobile Computing, 2018, 17(12): 2729–2743. doi: 10.1109/TMC.2018.2818723
    CHEN Mingzhe, SAAD W, YIN Changchuan, et al. Echo state networks for proactive caching in cloud-based radio access networks with mobile users[J]. IEEE Transactions on Wireless Communications, 2017, 16(6): 3520–3535. doi: 10.1109/TWC.2017.2683482
    PU Lingjun, JIAO Lei, CHEN Xu, et al. Online resource allocation, content placement and request routing for cost-efficient edge caching in cloud radio access networks[J]. IEEE Journal on Selected Areas in Communications, 2018, 36(8): 1751–1767. doi: 10.1109/JSAC.2018.2844624
    YANG Xiaolong, ZHENG Jianchao, FEI Zesong, et al. Optimal file dissemination and beamforming for cache-enabled C-RANs[J]. IEEE Access, 2018, 6: 6390–6399. doi: 10.1109/ACCESS.2017.2775198
    DAI Binbin, LIU Yafeng, and YU Wei. Optimized base-station cache allocation for cloud radio access network with multicast backhaul[J]. IEEE Journal on Selected Areas in Communications, 2018, 36(8): 1737–1750. doi: 10.1109/JSAC.2018.2844638
    YE Zhun, PAN Cunhua, ZHU Huiling, et al. Tradeoff caching strategy of the outage probability and fronthaul usage in a cloud-RAN[J]. IEEE Transactions on Vehicular Technology, 2018, 67(7): 6383–6397. doi: 10.1109/TVT.2018.2797957
    GOLNARI A, SHABANY M, NEZAMALHOSSEINI A, et al. Design and implementation of time and frequency synchronization in LTE[J]. IEEE Transactions on Very Large Scale Integration (VLSI) Systems, 2015, 23(12): 2970–2982. doi: 10.1109/TVLSI.2014.2387861
    ABDZADEH-ZIABARI H and SHAYESTEH M G. Robust timing and frequency synchronization for OFDM systems[J]. IEEE Transactions on Vehicular Technology, 2011, 60(8): 3646–3656. doi: 10.1109/TVT.2011.2163194
    KRÄMER J and SEEGER B. Semantics and implementation of continuous sliding window queries over data streams[J]. ACM Transactions on Database Systems, 2009, 34(1): Ariticle No.4. doi: 10.1145/1508857.1508861
    YU Wei and LUI R. Dual methods for nonconvex spectrum optimization of multicarrier systems[J]. IEEE Transactions on Communications, 2006, 54(7): 1310–1322. doi: 10.1109/TCOMM.2006.877962
    BOYD S and VANDENBERGHE L. Convex Optimization[M]. Cambridge: Cambridge University Press, 2004: 32–40.
  • 期刊类型引用(2)

    1. 田晨景,谢钧,曹浩彤,骆西建,刘亚群. 5G网络切片研究进展. 计算机科学. 2023(11): 282-295 . 百度学术
    2. 万剑锋. 通信光传送网子载波比特分配算法研究. 计算机仿真. 2020(12): 131-134+203 . 百度学术

    其他类型引用(2)

  • 加载中
图(5)
计量
  • 文章访问数:  2193
  • HTML全文浏览量:  901
  • PDF下载量:  95
  • 被引次数: 4
出版历程
  • 收稿日期:  2018-07-18
  • 修回日期:  2019-01-28
  • 网络出版日期:  2019-02-15
  • 刊出日期:  2019-07-01

目录

/

返回文章
返回