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

留言板

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

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

弹性光网络中时延感知的降级恢复路由与频谱分配算法

于存谦 张黎 何荣希 李靖宇

黄盛. 两用户非正交多址接入的最优时延均衡和功率控制方法[J]. 电子与信息学报, 2019, 41(8): 1902-1908. doi: 10.11999/JEIT180818
引用本文: 于存谦, 张黎, 何荣希, 李靖宇. 弹性光网络中时延感知的降级恢复路由与频谱分配算法[J]. 电子与信息学报, 2020, 42(10): 2420-2428. doi: 10.11999/JEIT190759
Sheng HUANG. Optimal Delay Equilibrium and Power Control Method for Two-user Non-orthogonal Multiple Access[J]. Journal of Electronics & Information Technology, 2019, 41(8): 1902-1908. doi: 10.11999/JEIT180818
Citation: Cunqian YU, Li ZHANG, Rongxi HE, Jingyu LI. Delay-aware Degradation-recovery Routing and Spectrum Allocation Algorithm in Elastic Optical Networks[J]. Journal of Electronics & Information Technology, 2020, 42(10): 2420-2428. doi: 10.11999/JEIT190759

弹性光网络中时延感知的降级恢复路由与频谱分配算法

doi: 10.11999/JEIT190759
基金项目: 国家自然科学基金(61371091, 61801074),中国博士后科学面上基金(2019M661074),辽宁省自然科学基金(2019-BS-021),中央高校基本科研业务费(3132020205, 3132019221)
详细信息
    作者简介:

    于存谦:男,1983年生,博士,副教授,研究方向为光数据中心网络

    张黎:女,1992年生,硕士,研究方向为弹性光网络

    何荣希:男,1971年生,博士,教授,研究方向为光网络和无线网络技术

    李靖宇:女,1995年生,硕士生,研究方向光数据中心网络

    通讯作者:

    何荣希 hrx@dlmu.edu.cn

  • 中图分类号: TN929.11

Delay-aware Degradation-recovery Routing and Spectrum Allocation Algorithm in Elastic Optical Networks

Funds: The National Natural Science Foundation of China (61371091, 61801074), Chian General Fundation for Postdoctoral Science (2019M661074), The Natural Science Foundation of Liaoning Province (2019-BS-021), The Fundamental Research Funds for The Central Universities (3132020205, 3132019221)
  • 摘要: 移动云计算、人工智能(AI)、5G等新兴技术应用促使弹性光网络(EON)在骨干传输网中发挥更重要的角色,降级服务(DS)技术为降低EON的业务阻塞率、提高频谱利用率提供了新途径。该文首先对现有DS算法的资源分配不公、忽略低等级业务的体验质量(QoE)等问题,建立了以最小化降级频次、降级等级与传输时延损失(TDL)为联合优化目标的混合整数线性规划(MILP)模型,并提出一种时延感知的降级恢复路由与频谱分配(DDR-RSA)算法。为提高降级业务的QoE和运营商收益,在算法的最优DS窗口选择阶段中融入降级恢复策略,在保障传输数据量不变的前提下,将降级业务向空闲频域复原,从而提高频谱效率、减小降级业务TDL和最大化网络收益。最后,通过仿真证明了所提算法在业务阻塞率、网络收益和降级业务成功率等方面的优势。
  • 未来无线通信网络将面临高密度用户、高速率及低时延业务等网络需求的挑战[1]。传统的通信系统常采用正交多址接入(Orthogonal Multiple Access, OMA)技术,如时分多址接入(Time-Division Multiple Access, TDMA)、频分多址接入(Frequency-Division Multiple Access, FDMA)、码分多址接入(Code-Division Multiple Access, CDMA)、正交频分复用(Orthogonal Frequency-Division Multiple Access, OFDMA)等技术。正交多址接入技术将网络资源在时域、频域或码域划分成正交的资源块,从而避免了用户之间的干扰。但是,正交的资源块也限制了单个网络的接入用户数。为了应对高密度用户的通信网络需求,在有限的频谱资源下,亟需更高效的多址接入技术来提升接入用户容量和频谱利用率。因此,非正交多址接入(Non-Orthogonal Multiple Access, NOMA)已成为5G通信系统中的研究热点之一。

    NOMA技术在功率域为不同的用户提供不同的增益,并采用叠加编码(Superposition Coding, SC)实现多路信号在相同时域、频域或码域的资源块上传输[2]。同时,接收端通过串行干扰消除(Successive Interference Cancellation, SIC)技术来区别多路用户信号,并提取本接收端所需的信息[3]。因此,NOMA技术能够成倍地提升通信网络系统的接入用户容量,并显著地提升频谱资源利用率。在前人研究工作中,文献[4]考虑了多用户功率分配、信令负载、SIC错误率和高移动性等网络因素,通过系统级仿真评估了NOMA技术的性能,实验结果表明NOMA技术的优越性。文献[5]分析了两用户非正交多址接入场景下网络容量与功率控制参数的变化关系,并验证了采用NOMA技术的网络总能够为用户提供不差于采用OMA技术时的性能。此外,由于NOMA技术需要在功率域区分不同用户的信号并通过SIC解调用户信号,因此NOMA技术的应用需要配合高效的功率控制算法。在不完美信道状态信息的场景下,文献[6]考虑了信道状态信息误差、用户服务质量需求等因素,以最小化发送功率为目标,研究了联合功率控制、速率控制、用户调度和SIC技术的多维度资源分配算法;文献[7]考虑以最大化网络能效为优化目标,提出了一种次优的联合用户调度与功率控制算法;文献[8]基于统计的信道状态信息提出了联合速率与功率分配的资源管理方法,在满足各用户的吞吐量需求的条件下最小化网络的总传输功率。基于完美信道状态信息的情况下,文献[9]提出了以最大化能量效率为目标的最优功率控制策略;文献[10]通过松弛方法研究了最小化功率消耗的方法;文献[11]提出了以最大化网络吞吐量为目标的联合用户配对与功率分配方法;文献[12]考虑传统电网与可再生能源的融合电能,研究了联合用户接入与功率分配的资源管理方法;文献[13]研究了异构无线网络下以最大化微小区用户的总吞吐量为目标的资源分配方法。但是,上述资源管理方法均需要通过循环迭代的方式来获得功率控制的次优解或最优解,通常具有较高的时间复杂度。因此,本文致力于研究简单可行的资源管理方法,通过推导功率分配的最优解的闭式表达式,使得基站能够在不需要循环迭代的情况下将各个时隙的瞬时信道状态信息与瞬时队列状态信息直接映射为该时隙下的功率控制结果,从而极大地降低了基站的计算开销。

    此外,随着视频分享网络、在线游戏、社交网络等时延敏感性业务的飞速发展,数据流的时延性能逐渐地成为了用户体验的重要指标,单纯最优化吞吐量无法提供良好的时延性能。例如,网络中某个节点具有较差的信道条件和很长的数据队列,仅基于信道条件的资源管理方法倾向于给该节点分配少量的网络资源以提升全网的资源利用率,进一步恶化了该节点的时延性能。因此,高效的资源管理方法应联合信道状态信息和队列状态信息,从而均衡节点排队时延并优化全网吞吐量。基于马尔科夫决策过程和贝尔曼方程,文献[14]设计了时延感知的联合信道分配和功率控制算法来最小化D2D终端的平均时延与平均丢弃率之和。基于李雅普诺夫优化方法(Lyapunov optimization),文献[15]提出了一种在线的功率控制算法来处理能量效率与排队时延之间的均衡关系。文献[16]研究了联合时延均衡和串行干扰消除的多能源融合网络资源管理方法。文献[17]研究了大规模多小区中联合时延均衡和功率控制的分布式资源管理方法。

    本文致力于研究在采用NOMA技术的单小区网络中联合时延均衡和功率控制的资源管理算法。由于功率控制在多条数据流之间的耦合性以及时延均衡和功率控制在长时间维度上的相互影响,联合时延均衡和功率控制的资源管理算法通常具有较高的复杂度。本文的主要贡献如下:

    (1) 考虑信道衰落的随机性、业务到达的动态性以及控制变量的耦合性,本文建模了联合时延均衡和功率控制的随机优化问题,在保障各用户的队列稳定性与基站的功率约束下,最大化网络的总吞吐量。其次,针对目标函数的时间平均特性和约束条件的参数在不同时隙的随机特性,本文采用李雅普诺夫优化方法将随机优化问题等效地分解为单时隙的时延均衡问题与功率控制问题。

    (2) 本文设计了一种低复杂度的最优时延均衡和功率控制的资源管理方法(DEPC-NOMA)。虽然单时隙的功率控制问题仍为非凸的优化问题,但是本文利用该优化问题的KKT(Karush-Kuhn-Tucker)条件的数学特征,通过理论推导获得了最优解在每个时隙的闭式表达式,极大地降低了功率控制算法的复杂度。

    (3) 本文通过仿真对比了DEPC-NOMA方法与其它方法的性能。仿真结果表明了本文提出的DEPC-NOMA方法能够有效地均衡各节点的排队时延并且显著地提升网络总吞吐量。

    本文研究单小区下行链路,一个基站采用NOMA技术在一个信道上同时服务两个用户,用户集合定义为U={1,2}。网络系统在时域上采用时隙结构,其中时隙t{1,2,···,}。定义时隙t的信道状态为G(t)={gu(t)}uU,其中,gu(t)的取值受影响于路径损耗、多径快衰落和阴影慢衰落。在信道的频谱带宽小于信道的相干带宽的情况下,本文考虑信道状态G(t)在不同时隙上是独立同分布的,且G(t)在单个时隙内保持恒定而在时隙边界上可能变化。在每个时隙tg1(t)g2(t)进行排序,并定义gl(t)gs(t)分别为集合G(t)={g1(t),g2(t)}中较大值的元素和较小值的元素,其下标ls对应地分别等于g1(t)g2(t)中较大值的下标和较小值的下标。通过在基站侧采用叠加编码(Superposition Coding, SC),在任意的功率控制分配结果下,基站利用自适应编码调制保证其给用户s传输的数据包能够被用户s成功地解调,则具有较优信道条件的用户l必定能够成功解调基站给用户s传输的数据包,并在用户l采用串行干扰消除技术剥离用户s的数据包带来的干扰[3]。因此,基站在时隙t给用户l和用户s的传输数据量可分别表示为

    Rl(t)=wτlog2(1+Pl(t)gl(t)σ2) (1)
    Rs(t)=wτlog2(1+Ps(t)gs(t)Pl(t)gs(t)+σ2) (2)

    其中,w表示信道带宽,τ表示时隙长度,Pu(t)为基站在时隙t给用户u发送数据的发射功率,σ2为白噪声功率。由于设备能力和功率节约的约束,基站的发射功率存在瞬时功率约束与时间平均的功率约率。基站的瞬时发射功率需大于等于零并且小于等于其允许的最大发射功率Pmax,而基站的时间平均发射功率需小于等于门限值Pmean,即

    Pu(t)0,uU,t (3)
    uUPu(t)Pmax,t (4)
    limsupT1T T1t=0E{uUPu(t)}Pmean (5)

    为了均衡数据排队时延,本文采用接纳控制方法来动态调整数据的接纳过程。考虑动态业务到达模型,将时隙t内到达基站的用户u的到达数据量表示为Au(t),而实际填入数据队列的接入数据量表示为ru(t)。由于接入数据量必定小于到达数据量,可得

    ru(t)[0,Au(t)],uU,t (6)

    基站为每个用户u维护一个缓存接入数据的队列。令Qu(t)表示发送节点u在时隙t的数据队列长度,则节点u的数据队列更新方程可定义为

    Qu(t+1)=max{Qu(t)Ru(t),0}+ru(t),uU (7)

    其中,ru(t)Ru(t)的取值受约束于业务产生速率和调制编码方案,因此两者分别存在上界值rmaxuRmaxu。此外,数据队列的稳定性需要保证时间平均上的数据队列长度为有限值,从而保障接入的业务数据能够在有限时间内离开数据队列,即

    limsupT1T T1t=0E{Qu(t)}<,uU (8)

    用户u的吞吐量是数据队列Qu在单位时间内输出的数据量。在有限排队时延的条件下,数据队列将趋于稳定状态,从而在时间平均上数据队列的输入数据量等于输出数据量。因此,由数据队列更新式(7)可得,在队列稳定状态下,网络u的吞吐量定义为

    limT1T T1t=0E{min{Ru(t),Qu(t)}}=limT1T T1t=0E{ru(t)}ˉru (9)

    综上,在保障基站功率限制和数据队列稳定性的约束下,以最大化网络总吞吐量为目标的联合时延均衡和功率控制的随机优化问题可建模为

    \begin{align} {{\rm P}1:} & \ {\mathop {\max }\limits_{_{\{ r(t),{{P}}(t)\} }} } \quad{\sum\limits_{u \in {\cal U}} {{{\bar r}_u}} }\\ {} & \ {{\rm{s}}{\rm{.t}}{\rm{.}}} \quad{{{式}}(3) {—} {{式}}(6),{{式}}(8)} \end{align} (10)

    其中,{{r}}(t) = {[{r_u}(t)]_{u \in \cal{U}}}{{P}}(t) = {[{P_u}(t)]_{u \in \cal{U}}}为随机优化问题的控制变量。首先,任意节点的在任意时隙的传输数据量不但取决于自身的信道条件,而且受其他节点的功率控制结果的影响。其次,由于网络变量的随机性以及时间平均上的约束条件式(5)和式(8),控制变量\{ {{r}}(t),{{P}}(t)\} 在不同时隙间的耦合性给求解优化问题P1带来了挑战。 虽然本文采用李雅普诺夫优化方法[18]将问题P1转换为单时隙的时延均衡与功率控制问题,然而功率控制优化问题的非凸目标函数使得求解该优化问题具有较高的复杂度。

    本节通过李雅普诺夫优化方法分析随机稳定性的特点,依据瞬时的网络状态进行每个时隙的控制决策,设计了低复杂度的最优联合时延均衡与功率控制方法,并推导了最优解的闭式表达式。

    为了处理时间平均上的功率约束式(5),本节定义虚拟队列Z(t)的队列长度更新方程定义为

    {Z(t + 1)} { = \max \left\{ {Z(t) + \sum\limits_{u \in \cal{U}} {{P_u}} (t) - {P_{\rm{mean}}},0} \right\}} (11)

    依据队列稳定理论(Rate stability theorem)[18],通过保证虚拟队列Z(t)在时间平均上的队列长度是有限的,即可确保虚拟队列Z(t)在时间平均上流入速率小于等于流出速率,从而满足时间平均上的功率约束式(5)。

    {{S}}(t) = [{Q_u}(t),Z(t)]表示时隙t的队列状态信息,依据Lyapunov Drift定理[18],定义李雅普诺夫函数(Lyapunov function)和单时隙的条件李雅普诺夫偏移(One-slot conditional Lyapunov drift)分别为L({{S}}(t)) = \left[\displaystyle\sum\nolimits_{u \in \cal{U}} {Q_u^2(t) + {Z^2}(t)} \right]/2 \Delta ({{S}}(t)) = {\rm E}\{ L({{S}}(t \!+\! 1)) \!-\! L({{S}}(t))|{{S}}(t)\} 。随后,依据Lyapunov优化定理[18],优化问题\rm{P}1的李雅普诺夫优化函数(Lyapunov drift-plus-penalty function)构建为

    \Delta ({{S}}(t)) - V\; {\rm E}\left\{ {\sum\limits_{u \in \cal{U}} {{r_u}} (t)|{{S}}(t)} \right\}, (12)

    其上界可通过定理1给出。

    定理1 假设{{G}}(t)在不同时隙上是独立同分布的,则对于任意的联合时延均衡和功率控制的资源管理方法,所有可能的{{S}}(t)和任意的大于零的控制参数V\;,优化问题P1的李雅普诺夫优化函数式(12)的上界为

    \begin{align} {F_{\rm upper}} =& {\rm E}\left\{ \sum\limits_{u \in \cal{U}} {({Q_u}(} t) - V\,){r_u}(t) - \sum\limits_{u \in \cal{U}} {{Q_u}} (t){R_u}(t)\right. \\ {\rm{}}& \left.+ \sum\limits_{u \in \cal{U}} Z (t){P_u}(t)|{{S}}(t) \right\} \!+ B + Z(t){P_{\rm{mean}}}\quad\quad \end{align} (13)

    其中,

    {B = } {\sum\limits_{u \in \cal{U}} {\frac{{{{(R_u^{\max })}^2} + {{(r_u^{\max })}^2}}}{2}} + \frac{{P_{\max }^2 - P_{\rm{mean}}^2}}{2}} (14)

    证明 首先,对数据队列更新式(7)与虚拟队列更新式(11)分别在等式两边进行取平方、移位、累加等操作可得

    \begin{align} {\rm{}}& \sum\limits_{u \in \cal{U}} {\frac{{Q_u^2(t + 1) - Q_u^2(t)}}{2}} \le \sum\limits_{u \in \cal{U}} {\frac{{R_u^2(t) + r_u^2(t)}}{2}} \\ {\rm{}}&\quad\quad - \sum\limits_{u \in \cal{U}} {{Q_u}} (t)({R_u}(t) - {r_u}(t)) \end{align} (15)
    \begin{align} {\rm{}}& \frac{{{Z^2}(t + 1) - {Z^2}(t)}}{2} \le \frac{{{{\left( {\displaystyle\sum\limits_{u \in \cal{U}} {{P_u}} (t)} \right)}^2} - P_{\rm{mean}}^2}}{2} \\ {\rm{}}& \quad\quad\ + Z(t)\left( {\sum\limits_{u \in \cal{U}} {{P_u}} (t) - {P_{\rm{mean}}}} \right) \end{align} (16)

    其中,不等式符号源自于不等式 {(\max [Q - R,0] + r)^2} \le{Q^2} + {R^2} + {r^2} + 2Q(R - r)。进一步,将不等式(15)和不等式(16)相加可得

    \begin{align} {\rm{}}& L({{S}}(t + 1)) - L({{S}}(t)) \\ {\rm{}}& \le \sum\limits_{u \in {\cal U}} {\frac{{R_u^2(t) + r_u^2(t)}}{2}} + \frac{{{{\left( {\displaystyle\sum\limits_{u \in {\cal U}} {{P_u}} (t)} \right)}^2} - P_{{\rm{mean}}}^2}}{2}\\ {\rm{}}&\quad+ \sum\limits_{u \in {\cal U}} {{Q_u}} (t){r_u}(t) - \left[ \sum\limits_{u \in {\cal U}} {{Q_u}} (t){R_u}(t)\right. \\ {\rm{}} & \quad\left.+ Z(t)\left( {\sum\limits_{u \in {\cal U}} {{P_u}} (t) - {P_{{\rm{mean}}}}} \right) \right] \end{align} (17)

    给定{{S}}(t),对不等式(17)两边同时取条件期望 并减去V\; {\rm E}\left\{ {\sum\limits_{u \in \cal{U}} {{r_u}} (t)|{{S}}(t)} \right\}可得

    \begin{align} \Delta ({{S}}(t)) &- V\; {\rm E}\left\{ {\sum\limits_{u \in \cal{U}} {{r_u}} (t)|{{S}}(t)} \right\} \\ {\rm{}}& \le B + {\rm E}\left\{ {\sum\limits_{u \in \cal{U}} {({Q_u}(} t) - V){r_u}(t)|{{S}}(t)} \right\} \\ {\rm{}}& - {\rm E}\left\{ \left[ \sum\limits_{u \in \cal{U}} {{Q_u}} (t){R_u}(t)- Z(t)\right.\right. \\ {\rm{}}& \left.\left.\left( {\sum\limits_{u \in \cal{U}} {{P_u}} (t) - {P_{\rm{mean}}}} \right) \right]\Biggr|{{S}}(t) \right\} \end{align} (18)

    其中, {B}= {\displaystyle\sum\limits_{u \in \cal{U}} {\dfrac{{{{(R_u^{\max })}^2} + {{(r_u^{\max })}^2}}}{2}} + \dfrac{{P_{\max }^2 - P_{\rm{mean}}^2}}{2}}

    证毕

    依据Min Drift-Plus-Penalty定理[18],对于给定的{{S}}(t),优化问题P1可以转换为最小化上界函数式(13)的优化问题,并受限于瞬时约束式(3)、式(4)和式(6)。研究上界函数式(13)的结构可得:

    (1) 给定{{S}}(t),该上界函数的最后2项在时隙t内是常数;

    (2) 给定{{S}}(t),最小化该上界函数的期望等同于根据当前的信道状态信息和队列状态信息的情况来最小化该期望内部的函数;

    (3) 给定{{S}}(t),该上界函数中期望内部的函数可以分解为2个相互独立的函数并独立地进行最小化操作。具体地,最小化上界函数式(13)中期望内部的函数可独立地分解为:(a)最小化该期望中的第1项对应的时延均衡问题,该时延均衡问题依据各用户的数据队列长度对当前到达的数据进行接纳控制,从而有效地均衡数据传输的时延;(b)最小化该期望中的第2项对应的功率控制问题,该功率控制问题依据各用户的数据队列长度与信道条件建模功率控制对吞吐量的增益并维持关于功率消耗的虚拟队列Z(t)的有限队列长度。

    最小化上界函数式(12)的期望中的第1项可建模为如下优化问题P2

    \left. \begin{aligned} {{{\rm P}}2:} & {\mathop {\min }\limits_{_{r(t)}} } \quad{\:\sum\limits_{u \in {\cal U}} {({Q_u}(} t) - V){r_u}(t)}\\ {} & {{\rm{s}}{\rm{.t}}{\rm{.}}} \quad {{{式}}(6)} \end{aligned}\right\} (19)

    优化问题P2为线性规划问题[19],其在时隙t的最优解{{{r}}^ * }(t)

    {{r}_u^ * (t)} { = {A_u}(t) \cdot {{{1}}_{\{ {Q_u}(t) < V\,\} }},\forall u \in \cal{U}} (20)

    最优解{{{r}}^ * }(t)表明,调整控制参数V\;可以限制数据队列的最大队列长度,从而达到均衡排队时延的目的。

    在每个时隙t下,功率控制的决策是从信道状态信息{{G}}(t)和队列状态信息{{S}}(t)到发射功率向量{{P}}(t)的映射。在两用户的场景下,最小化上界函数式(12)的期望中的第2项可建模为如下优化问题P3

    \begin{align} {\rm{}}& {{{\rm P}}3:} \\ {\rm{}}& {\mathop {\min }\limits_{_{{\hat{{P}}}(t)}} } \ Z(t)({P_l}(t) \!+\! {P_s}(t)) \\ {\rm{}}& \quad\quad- {Q_l}(t)w\tau {{\log }_2}\left( {1 \!+\! \frac{{{P_l}(t){g_l}(t)}}{{{\sigma ^2}}}} \right)\\ {\rm{}}{\rm{}}& \quad\quad { - {Q_s}(t)w\tau {{\log }_2}\left( {1 + \frac{{{P_s}(t){g_s}(t)}}{{{P_l}(t){g_s}(t) + {\sigma ^2}}}} \right)}\\ {\rm{}}& {{\rm{s}}{\rm{.t}}{\rm{.}}} \quad\quad\quad\quad {{{式}}(3),{{式}}(4)} \end{align} (21)

    其中,{\hat{{P}}}(t) = \{ {P_l}(t),{P_s}(t)\} 。需要注意的是,{R_u}(t)的信干噪比中发射功率的相互耦合关系使得优化问题P3的目标函数为非凸函数。如定理2所示,本节分析了优化问题P3的KKT条件[19],并利用KKT条件的函数特性推导了优化问题P3的最优解的闭式表达式,极大地减少了功率控制方法的复杂度。

    定理2:如果优化问题P3取得了最优解{\hat {{P}}^*}(t),则最优解{\hat {{P}}^*}(t)必定为最优解集合\tilde {{P}}(t)中的元素

    \begin{align} \tilde {{P}}(t) =& \left\{ \left(P_l^{(1)}(t),P_s^{(1)}(t)\right),\left(P_l^{(2)}(t),P_s^{(2)}(t)\right),\right.\\ {\rm{}}& \left.\left(P_l^{(3)}(t),P_s^{(3)}(t)\right)\right\} \end{align} (22)

    其中,

    P_l^{(1)}(t) = \frac{{({Q_s}(t){g_s}(t) - {Q_l}(t){g_l}(t)){\sigma ^2}}}{{({Q_l}(t) - {Q_s}(t)){g_l}(t){g_s}(t)}} \hspace{90pt} (23)
    \begin{align} P_s^{(1)}(t) =& \min \biggr\{ \frac{{w\tau {Q_s}(t)}}{{Z(t)}} - \frac{{{\sigma ^2}}}{{{g_s}(t)}} - P_l^{(1)}(t), \\ {\rm{}}& {P_{\max }} - P_l^{(1)}(t) \biggr\}\hspace{160pt} \end{align} (24)
    \begin{align} {\rm{}}&P_l^{(2)}(t) = \max \left\{ {0,\min \left\{ {\frac{{w\tau {Q_l}(t)}}{{Z(t)}} - \frac{{{\sigma ^2}}}{{{g_l}(t)}},{P_{\max }}} \right\}} \right\}\quad\\ {\rm{}}& \hspace{212pt}(25) \end{align}
    P_s^{(2)}(t) = 0 \hspace{183pt} (26)
    P_l^{(3)}(t) = 0 \hspace{183pt} (27)
    \begin{align} {\rm{}}& P_s^{(3)}(t) = \max \left\{ {0,\min \left\{ {\frac{{w\tau {Q_s}(t)}}{{Z(t)}} - \frac{{{\sigma ^2}}}{{{g_s}(t)}},{P_{\max }}} \right\}} \right\} \ \ \\ {\rm{}}& \hspace{211pt}(28) \end{align}

    证明 依据KKT条件定义[19],优化问题P3的最优解满足KKT条件式(29)—式(33)

    \begin{align} {\rm{}}& \frac{{{Q_l}(t)w\tau {g_l}(t)}}{{{P_l}(t){g_l}(t) + {\sigma ^2}}} - \frac{{{Q_s}(t)w\tau {g_s}(t)}}{{{P_l}(t){g_s}(t) + {\sigma ^2}}} \\ {\rm{}}& \quad+ \frac{{{Q_s}(t)w\tau {g_s}(t)}}{{({P_l}(t) + {P_s}(t)){g_s}(t) + {\sigma ^2}}}\! -\! Z(t) \!+\! {\lambda _l} \!-\! \mu = 0\quad\quad\quad\; \end{align} (29)
    \frac{{{Q_s}(t)w\tau {g_s}(t)}}{{({P_l}(t) + {P_s}(t)){g_s}(t) + {\sigma ^2}}} - Z(t) + {\lambda _s} - \mu = 0 \quad\quad\ (30)
    {\lambda _l}{P_l}(t) = 0 \hspace{180pt} (31)
    {\lambda _s}{P_s}(t) = 0 \hspace{180pt} (32)
    \mu ({P_l}(t) + {P_s}(t) - {P_{\max }}) = 0 \hspace{90pt} (33)

    其中,{\lambda _l}{\lambda _s}的取值可分为以下3种情况:

    (1) 当{\lambda _l} = {\lambda _s} = 0时,优化问题P3的最优解满足条件式(34)—式(36)

    \begin{align} {\rm{}}& \frac{{{Q_l}(t)w\tau {g_l}(t)}}{{{P_l}(t){g_l}(t) + {\sigma ^2}}} - \frac{{{Q_s}(t)w\tau {g_s}(t)}}{{{P_l}(t){g_s}(t) + {\sigma ^2}}} \\ {\rm{}}& + \frac{{{Q_s}(t)w\tau {g_s}(t)}}{{({P_l}(t) + {P_s}(t)){g_s}(t) + {\sigma ^2}}} - Z(t) - \mu = 0 \end{align} (34)
    \frac{{{Q_s}(t)w\tau {g_s}(t)}}{{({P_l}(t) + {P_s}(t)){g_s}(t) + {\sigma ^2}}} - Z(t) - \mu = 0 \quad\ (35)
    \mu ({P_l}(t) + {P_s}(t) - {P_{\max }}) = 0 \hspace{65pt} (36)

    求解式(34)—式(36)可得

    P_l^{(1)}(t) = \frac{{({Q_s}(t){g_s}(t) - {Q_l}(t){g_l}(t)){\sigma ^2}}}{{({Q_l}(t) - {Q_s}(t)){g_l}(t){g_s}(t)}} \hspace{25pt} (37)
    \begin{align} P_s^{(1)}(t) =& \min \left\{ \frac{{w\tau {Q_s}(t)}}{{Z(t)}} - \frac{{{\sigma ^2}}}{{{g_s}(t)}} - P_l^{(1)}(t),\right.\\ {\rm{}}& {P_{\max }} - P_l^{(1)}(t) \biggr\} \end{align} (38)

    (2) 当{\lambda _s} \ne 0时,则优化问题P3的最优解满足P_s^{(2)}(t) = 0

    \frac{{{Q_l}(t)w\tau {g_l}(t)}}{{{P_l}(t){g_l}(t) + {\sigma ^2}}} - Z(t) + {\lambda _l} - \mu = 0 (39)
    {\lambda _l}{P_l}(t) = 0 \hspace{105pt} (40)
    \mu ({P_l}(t) - {P_{\max }}) = 0 \hspace{67pt} (41)

    求解式(39)—式(41)可得

    \begin{align} {\rm{}}&P_l^{(2)}(t) = \max \left\{ {0,\min \left\{ {\frac{{w\tau {Q_l}(t)}}{{Z(t)}} - \frac{{{\sigma ^2}}}{{{g_l}(t)}},{P_{\max }}} \right\}} \right\}\\ {\rm{}}& \hspace{212pt}(42) \end{align}

    (3) 当{\lambda _l} \ne 0时,则优化问题P3的最优解满足P_l^{(3)}(t) = 0

    \frac{{{Q_s}(t)w\tau {g_s}(t)}}{{{P_s}(t){g_s}(t) + {\sigma ^2}}} - Z(t) + {\lambda _s} - \mu = 0 (43)
    {\lambda _s}{P_s}(t) = 0 \hspace{107pt} (44)
    \mu ({P_s}(t) - {P_{\max }}) = 0 \hspace{72pt} (45)

    求解式(43)—式(45)可得

    \begin{align} {\rm{}}&P_s^{(3)}(t) = \max \left\{ {0,\min \left\{ {\frac{{w\tau {Q_s}(t)}}{{Z(t)}} - \frac{{{\sigma ^2}}}{{{g_s}(t)}},{P_{\max }}} \right\}} \right\}\\ {\rm{}}& \hspace{212pt}(46) \end{align}

    综上,优化问题{\rm P}3的最优解必定属于最优解集合\tilde {{P}}(t)的元素。 证毕

    依据Min Drift-Plus-Penalty定理[18],数据队列式(7)和虚拟队列式(11)是稳定的,从而满足了时间平均上的有限队列长度约束式(8)和功率约束式(5)。本节通过MATLAB仿真验证本文提出的最优时延均衡和功率控制的非正交多址接入方法(DEPC-NOMA)在网络吞吐量和数据队列长度上的性能。本文采用如下2种方法做性能对比:

    (1) DEEP-NOMA:该方法的时延均衡方法与DEPC-NOMA方法相同。其次,该方法的功率控制采用均分的方式(Equal Power, EP),基站在每个时隙t给两个用户分配相同的发射功率,并保证满足功率约束。

    (2) DEPC-TDMA:该方法的时延均衡方法与DEPC-NOMA方法相同。其次,两用户通过时分多址接入(TDMA)的方式共享频率与功率资源,则该方法的功率控制的最优解{\hat {{P}}^*}(t)属于集合式(47)

    \begin{align} {\rm{}}&{\tilde {{P}}_{\rm TDMA}}(t) \!=\! \left\{ \left(P_l^{(2)}(t),P_s^{(2)}(t)),(P_l^{(3)}(t),P_s^{(3)}(t)\right)\right\} \\ {\rm{}}& \hspace{212pt}(47) \end{align}

    考虑一个半径为250 m的圆形小区,基站位于中心位置,用户1和用户2分别随机分布在距离基站125 m和200 m的圆环上。信道带宽w为10 MHz。每个时隙长度\tau 为1 ms。基站的瞬时功率约束{P_{\max }}与时间平均的功率约束{P_{\rm{mean}}}分别为20 dBm与18 dBm。噪声功率谱密度{\sigma ^2}为–174 dBm/Hz。每个时隙内的平均数据到达量为50 kbit/slot。路径损耗模型采用PL = 128.1 + 37.6{\log _2}d,其中,路径距离d单位为km。多径快衰落服从相互独立的均值为1的指数分布。阴影慢衰落模型采用标准方差为8 dB的独立对数正态分布。在仿真中,数值结果为50次拓扑生成且每个拓扑运行2000个时隙所获得的平均值结果。

    图1描述了参数V的取值范围为[1, 400]时平均数据队列长度与网络总吞吐量的权衡关系。在平均数据队列长度较小时,随着平均数据队列长度逐渐增加,用户有充足的缓存数据来保证网络传输机会得到高效的利用,从而网络总吞吐量迅速地上升。而在平均数据队列长度较大时,网络传输机会的利用趋于饱和状态,从而网络总吞吐量随着平均数据队列长度的增加而缓慢上升。此外,对比DEEP-NOMA方法和DEPC-TDMA方法,本文提出的DEPC-NOMA方法能够显著地提升网络吞吐量并降低排队时延。

    图 1  平均数据队列长度与网络总吞吐量的权衡关系

    图2图3分别说明了各用户的吞吐量和平均数据队列长度与控制参数V\;的关系。各用户的吞吐量随着控制参数V\;的增大而增大,但是各用户的平均数据队列长度及排队时延也随着控制参数V\;的增大而增大。相对于DEPC-TDMA方法,DEPC-NOMA方法能够创造更多的传输机会,从而有效地提升用户的吞吐量性能,同时保持良好的排队时延性能。此外,在DEEP-NOMA方法中,由于其功率控制未考虑时延均衡的需求,虽然信道条件较好的用户1能够获得很高的吞吐量与很低的排队时延,但是信道条件较差的用户2却忍受着极低的吞吐量与极高的排队时延,用户间的公平性极差。在本文提出的DEPC-NOMA方法中,通过联合时延均衡与功率控制,虽然用户1的吞吐量与排队时延性能稍微下降,但是用户2的吞吐量与排队时延能够获得极大地性能提升。

    图 2  各用户的吞吐量与控制参数V 的关系
    图 3  各用户的平均数据队列长度与控制参数V 的关系

    针对采用非正交多址接入技术的两用户单小区网络,考虑信道衰落的随机性、业务到达的动态性以及控制变量的耦合性,本文提出了一种低复杂度的最优时延均衡和功率控制的资源管理方法(DEPC-NOMA),以期最大化网络总吞吐量,并满足数据队列稳定性与基站功率约束。仿真结果表明,本文提出的DEPC-NOMA方法能够显著地提升网络总吞吐量并均衡用户的排队时延。

  • 图  1  降级恢复举例

    图  2  4种算法的数据量阻塞率

    图  3  3种算法的降级服务成功率

    图  4  3种算法的平均延迟时间

    图  5  4种算法的网络收益

    图  6  3种算法的降级等级占比

    表  1  RSA问题符号定义

    变量定义内容
    \overline \omega 正整数,\psi 中的业务优先级上界;
    {w_r}正整数,r所在的起始频谱槽序号;
    f_{u,v}^r二值变量,若r经过光纤链路e(u,v) \in E,则f_{u,v}^r = 1;否则f_{u,v}^r = 0
    {\rho _{i,j}}二值变量,若{r_i}{r_j}经过同一段光纤链路,且{w_i}{w_j}小,则{\rho _{i,j}} = 1;否则{\rho _{i,j}} = 0
    \xi _{s,u}^r二值变量,若r的源节点为u \in N,则\xi _{s,u}^r = 1;否则,\xi _{s,u}^1 = 0
    \xi _{d,v}^r二值变量,若r的目的节点为v \in N,则\xi _{d,v}^r = 1;否则,\xi _{d,v}^r = 0
    {\delta _r}二值变量,若r降级,则{\delta _r} = 1;否则,{\delta _r}{\rm{ = 0}}
    {\chi _r}正整数,r释放的频谱槽数;
    {\beta _r}正整数,r恢复的频谱槽数;
    {v_r}{z_r}正整数,DR后r首/尾频谱槽序号;
    q_k^e正实数,第k个频谱槽可被r用来DR的起始时间;
     t_r^{{\rm{end}}'}正实数,r被降级后的离开时间。
    下载: 导出CSV

    表  2  启发式算法部分的变量

    变量定义内容
    u_{t,c}^{l,k}二值变量,若{p_k}中第l条链路的第c位频谱槽的第t时隙被占用,则u_{t,c}^{l,k} = 1;否则u_{t,c}^{l,k} = 0
    u_{t,c}^p二值变量,若{p_k}的第c位频谱槽的第t时隙被占用,则u_{t,c}^p = 1;否则,u_{t,c}^p = 0
    B_{b,e}^{k,h}{p_k}的空闲频谱窗口,其频谱槽首、末序号为be,时长为h,含频谱槽数为n_{b,e}^{k,h} = e - b + 1
    \tau _{b,e}^{k,h}正整数,B_{b,e}^{k,h}为满足r的带宽尚需的频谱槽数;
    {\chi _{r'}}正实数,降级业务r'释放的频谱槽数;
    \tau _{b,e}^{{\rm{left}},l},\tau _{b,e}^{{\rm{right}},l}正整数,B_{b,e}^{k,h}的每条链路上[b - \tau _{b,e}^{k,h},b)(e,e + \tau _{b,e}^{k,h}]内最少可释放的频谱槽数,l \in {p_k}
    \tau _{b,e}^{{\rm{left}}},\tau _{b,e}^{{\rm{right}}}正整数,B_{b,e}^{k,h}所在路径上[b - \tau _{b,e}^{k,h},b)(e,e + \tau _{b,e}^{k,h}]内可释放的频谱槽数;
    {b_{r'}}正整数,可降级业务r'占用的带宽;
    {\rm{ho}}{{\rm{p}}_{r'}}正整数,r'所在路径的链路数;
    \theta _{r'}^{{\chi _{r'}}}正实数,r'释放的数据量;
    \left[ {s,d} \right]r'占用的频谱,有d - s + 1 = {b_{r'}}
    \theta _t^{r'}正实数,r'在第t时隙可恢复的数据量;
    \theta _{r'}'正实数,r'可恢复数据量之和;
    t_{r'}^{{\rm{end}}'}正实数,r'降级后的离去时间。
    [b - \tau _{b,e}^{k,h},b),(e,e + \tau _{b,e}^{k,h}]B_{b,e}^{k,h}的左/右两侧的降级备选区间;
    [s - {\chi _{r'}},d - {\chi _{r'}}],[s + {\chi _{r'}},d + {\chi _{r'}}]r'左/右两侧分别可恢复的频域;
    下载: 导出CSV

    表  3  DR策略伪码

     输入:\psi , {{G}}\left( {N,E,C} \right).
     输出:t_{r'}^{{\rm{end}}'} and {\chi _{r'}}.
     (1) if {o_{r'}} < {o_r} then
     (2)  t_{r'}^{{\rm{end}}'} \leftarrow t_{r'}^{{\rm{end}}}, {\chi _{r'}} \leftarrow 0; calculate {\chi _{r'}}, \theta _{r'}^{{\chi _{r'}}} in
         [b - \tau _{b,e}^{k,h},b) \cup (e,e + \tau _{b,e}^{k,h}];
     (3)  {\left[ {{{S}}_k^p} \right]_{{\rm{T}} \times \left| {\rm{C}} \right|}} \leftarrow {\left[ {{{U}}_k^l} \right]_{{\rm{T}} \times \left| {\rm{C}} \right|}} in
         \left[ {s + {\chi _{r'}},d + {\chi _{r'}}} \right] \cup \left[ {s - {\chi _{r'}},d - {\chi _{r'}}} \right];
     (4)  while t \ge t_r^a and t \le {\overline \alpha _{r'} } do
     (5)   for u_{t,c}^{{p_{r'}}} = 0 do \theta _t^{r'} \leftarrow Eq.24, c + + ; end for
     (6)   if u_{t,c}^{{p_{r'}}} = 1 then
     (7)    if c \ge c' in u_{t - 1,c'}^{{p_{r'}}} = 1 then calculate
           \theta _{r'}' \leftarrow Eq.25, t++;
     (8)    else then
     (9)     \theta _{r'}' \leftarrow Eqs. 24-25(\left( {c,d - {\chi _{r'}}} \right] or
            \left[ {s + {\chi _{r'}},c} \right)), t++;
     (10)    end if
     (11)   end if
     (12)   if \theta _{r'}' \ge \theta _{r'}^{ {\chi _{r'} } } then return t_{r'}^{{\rm{end}}'} and {\chi _{r'}}; end if
     (13)  end while
     (14)  if \theta _{r'}' < \theta _{r'}^{{\chi _{r'}}} then set{\chi _{r'}} = {\chi _{r'}} - 1; jump to
         Line 1; end if
     (15)  if {\chi _{r'}} = = 0 then return 0; end if
     (16) end if
    下载: 导出CSV
  • ZHU Zuqing, LU Wei, ZHANG Liang, et al. Dynamic service provisioning in elastic optical networks with hybrid single-/multi-path routing[J]. Journal of Lightwave Technology, 2013, 31(1): 15–22. doi: 10.1109/JLT.2012.2227683
    GU Yamei, YUAN Xin, and YOU Shanhong. Blocking performances for fixed/flexible-grid and sub-band conversion in optical networks[C]. 2014 IEEE/CIC International Conference on Communications in China, Shanghai, China, 2014: 116–120. doi: 10.1109/ICCChina.2014.7008254.
    朱丹, 徐威远, 陈文娟, 等. 基于光波分复用网络的分布式多目标定位系统[J]. 雷达学报, 2019, 8(2): 171–177. doi: 10.12000/JR19028

    ZHU Dan, XU Weiyuan, CHEN Wenjuan, et al. Distributed multi-target localization system based on optical wavelength division multiplexing network[J]. Journal of Radars, 2019, 8(2): 171–177. doi: 10.12000/JR19028
    SAVAS S S, HABIB M F, TORNATORE M, et al. Exploiting degraded-service tolerance to improve performance of telecom networks[C]. Optical Fiber Communication Conference, San Francisco, USA, 2014: W2A. 31. doi: 10.1364/OFC.2014.W2A.31.
    MUHAMMAD A, CAVDAR C, WOSINSKA L, et al. Service differentiated provisioning in dynamic WDM networks based on set-up delay tolerance[J]. Journal of Optical Communications and Networking, 2013, 5(11): 1250–1261. doi: 10.1364/JOCN.5.001250
    SAVAS S S, HABIB M F, TORNATORE M, et al. Network adaptability to disaster disruptions by exploiting degraded-service tolerance[J]. IEEE Communications Magazine, 2014, 52(12): 58–65. doi: 10.1109/MCOM.2014.6979953
    VADREVU C S K, WANG Rui, TORNATORE M, et al. Degraded service provisioning in mixed-line-rate WDM backbone networks using multipath routing[J]. IEEE/ACM Transactions on Networking, 2014, 22(3): 840–849. doi: 10.1109/TNET.2013.2259638
    ZHONG Zhizhen, LI Jipu, HUA Nan, et al. On QoS-assured degraded provisioning in service-differentiated multi-layer elastic optical networks[C]. 2016 IEEE Global Communications Conference, Washington, USA, 2016: 1–5. doi: 10.1109/GLOCOM.2016.7842043.
    SANTOS A S, HOROTA A K, ZHONG Zhizhen, et al. An online strategy for service degradation with proportional QoS in elastic optical networks[C]. 2018 IEEE International Conference on Communications, Kansas City, USA, 2018: 1–6. doi: 10.1109/ICC.2018.8422781.
    HOU Weigang, GUO Lei, NING Zhaolong, et al. Appropriate service degradability for virtualized inter-data-center optical networks[C]. 2018 IEEE Global Communications Conference, Abu Dhabi, UAE, 2018: 206–212. doi: 10.1109/GLOCOM.2018.8647198.
    HOU Weigang, NING Zhaolong, GUO Lei, et al. Service degradability supported by forecasting system in optical data center networks[J]. IEEE Systems Journal, 2019, 13(2): 1514–1525. doi: 10.1109/JSYST.2018.2821714
    ROJAS J S, GALLÓN Á R, and CORRALES J C. Personalized service degradation policies on OTT applications based on the consumption behavior of users[C]. The 18th International Conference on Computational Science and Its Applications, Melbourne, Australia, 2018: 543–557. doi: 10.1007/978-3-319-95168-3_37.
    于存谦, 张黎, 何荣希. 弹性光网络基于区分降级服务和自适应调制的动态路由与频谱分配算法[J]. 电子与信息学报, 2019, 41(1): 38–45. doi: 10.11999/JEIT180075

    YU Cunqian, ZHANG Li, and HE Rongxi. Dynamic routing and spectrum assignment algorithm based on differentiated degraded-service and adaptive modulation in elastic optical networks[J]. Journal of Electronics &Information Technology, 2019, 41(1): 38–45. doi: 10.11999/JEIT180075
    YANG Anjia, WENG Jian, CHENG Nan, et al. DeQoS attack: Degrading quality of service in VANETs and its mitigation[J]. IEEE Transactions on Vehicular Technology, 2019, 68(5): 4834–4845. doi: 10.1109/TVT.2019.2905522
    PROTO S, VENTURA F, APILETTI D, et al. PREMISES, a scalable data-driven service to predict alarms in slowly-degrading multi-cycle industrial processes[C]. 2019 IEEE International Congress on Big Data, Milan, Italy, 2019: 139–143. doi: 10.1109/BigDataCongress.2019.00032.
    张海波, 李虎, 陈善学, 等. 超密集网络中基于移动边缘计算的任务卸载和资源优化[J]. 电子与信息学报, 2019, 41(5): 1194–1201. doi: 10.11999/JEIT180592

    ZHANG Haibo, LI Hu, CHEN Shanxue, et al. Computing offloading and resource optimization in ultra-dense networks with mobile edge computation[J]. Journal of Electronics &Information Technology, 2019, 41(5): 1194–1201. doi: 10.11999/JEIT180592
  • 期刊类型引用(18)

    1. 寇子明,李腾宇,吴娟. 矿井提升机天轮智能润滑系统设计及控制策略. 煤炭学报. 2024(02): 1230-1239 . 百度学术
    2. 崔胜胜,李汐,牟颖莹,李振. 间接接入式直流电能表异常计量数据识别算法. 微型电脑应用. 2024(07): 130-133 . 百度学术
    3. 毛冬,张辰,陈又咏,刘永清,焦艳斌. 基于国产CPU环境的国产数据库历史数据迁移技术. 太赫兹科学与电子信息学报. 2024(10): 1154-1160 . 百度学术
    4. 潘海霞,曹宁. 基于SDN的数据中心端口故障检测方法设计. 自动化与仪器仪表. 2023(02): 56-60 . 百度学术
    5. 刘浪,时宏伟. 基于注意力机制的CNN-LSTM的ADS-B异常数据检测. 计算机系统应用. 2023(04): 94-103 . 百度学术
    6. 卢宁. 基于加权广义匹配滤波联合SVDD的FOD检测算法. 现代雷达. 2023(05): 112-119 . 百度学术
    7. 何佳佳,安大庆,叶美丽. 面向学生体质测量的数据异常可视化检测方法. 信息技术. 2023(06): 83-87+93 . 百度学术
    8. 王鹏,董方正,张伟,汪克念,赵长啸. 基于TriLSTM-SVDD的ADS-B欺骗数据检测方法. 电光与控制. 2023(10): 77-81+119 . 百度学术
    9. 李妮锶,林琳,姚成. ADS-B系统安全问题解决策略对比分析. 信息技术. 2022(03): 30-35+41 . 百度学术
    10. 杨章勇,李欢. 工业以太网通信协议关键方法的优化与仿真. 计算机仿真. 2022(03): 377-380+395 . 百度学术
    11. 符士侃,夏元轶,杜钰,石廷川. 基于概率主题模型的大数据平台隐私泄露自动检测方法. 自动化与仪器仪表. 2022(04): 115-118+123 . 百度学术
    12. 曾雄飞. 基于粒子群算法优化BP神经网络的PID控制算法. 电子设计工程. 2022(11): 69-73+78 . 百度学术
    13. 汤双霞. 一种毫米波雷达机场跑道异物检测算法. 激光与红外. 2022(06): 820-826 . 百度学术
    14. 暴佳伟,田小平,刘宇娜,邹长宽,张雨晴. ADS-B安全问题研究综述. 现代防御技术. 2022(05): 28-35 . 百度学术
    15. 李娟,李海瑛,边玲,刘秋红. 基于跨模态深度度量学习的时序数据异常检测. 计算机仿真. 2022(10): 533-537 . 百度学术
    16. 刘明,张弘. AI交互终端异常数据入侵识别与仿真. 计算机仿真. 2022(11): 467-471 . 百度学术
    17. 张维,杜兰. 一种集成式Beta过程最大间隔一类分类方法. 电子与信息学报. 2021(05): 1219-1227 . 本站查看
    18. 郭峰. 远距离多信道光纤通信网络异常数据检测研究. 激光杂志. 2021(09): 98-103 . 百度学术

    其他类型引用(13)

  • 加载中
图(6) / 表(3)
计量
  • 文章访问数:  1934
  • HTML全文浏览量:  578
  • PDF下载量:  51
  • 被引次数: 31
出版历程
  • 收稿日期:  2019-09-13
  • 修回日期:  2020-06-15
  • 网络出版日期:  2020-07-17
  • 刊出日期:  2020-10-13

目录

/

返回文章
返回