CN111277307A - 一种有限反馈欠秩信道时mu-mimo系统的资源分配方法 - Google Patents

一种有限反馈欠秩信道时mu-mimo系统的资源分配方法 Download PDF

Info

Publication number
CN111277307A
CN111277307A CN202010069963.2A CN202010069963A CN111277307A CN 111277307 A CN111277307 A CN 111277307A CN 202010069963 A CN202010069963 A CN 202010069963A CN 111277307 A CN111277307 A CN 111277307A
Authority
CN
China
Prior art keywords
user
optimization problem
rank
dual
channel
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.)
Granted
Application number
CN202010069963.2A
Other languages
English (en)
Other versions
CN111277307B (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.)
Nanjing University of Posts and Telecommunications
Original Assignee
Nanjing University of Posts and Telecommunications
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 Nanjing University of Posts and Telecommunications filed Critical Nanjing University of Posts and Telecommunications
Priority to CN202010069963.2A priority Critical patent/CN111277307B/zh
Publication of CN111277307A publication Critical patent/CN111277307A/zh
Application granted granted Critical
Publication of CN111277307B publication Critical patent/CN111277307B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems
    • H04B7/0426Power distribution
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems
    • H04B7/0452Multi-user MIMO systems
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems
    • H04B7/0456Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
    • H04B7/0486Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting taking channel rank into account
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W52/00Power management, e.g. TPC [Transmission Power Control], power saving or power classes
    • H04W52/04TPC
    • H04W52/18TPC being performed according to specific parameters
    • H04W52/26TPC being performed according to specific parameters using transmission rate or quality of service QoS [Quality of Service]
    • H04W52/265TPC being performed according to specific parameters using transmission rate or quality of service QoS [Quality of Service] taking into account the quality of service QoS
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W52/00Power management, e.g. TPC [Transmission Power Control], power saving or power classes
    • H04W52/04TPC
    • H04W52/18TPC being performed according to specific parameters
    • H04W52/26TPC being performed according to specific parameters using transmission rate or quality of service QoS [Quality of Service]
    • H04W52/267TPC being performed according to specific parameters using transmission rate or quality of service QoS [Quality of Service] taking into account the information rate
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W72/00Local resource management
    • H04W72/04Wireless resource allocation
    • H04W72/044Wireless resource allocation based on the type of the allocated resource
    • H04W72/0473Wireless resource allocation based on the type of the allocated resource the resource being transmission power
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W72/00Local resource management
    • H04W72/50Allocation or scheduling criteria for wireless resources
    • H04W72/54Allocation or scheduling criteria for wireless resources based on quality criteria
    • H04W72/542Allocation or scheduling criteria for wireless resources based on quality criteria using measured or perceived quality
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W72/00Local resource management
    • H04W72/50Allocation or scheduling criteria for wireless resources
    • H04W72/54Allocation or scheduling criteria for wireless resources based on quality criteria
    • H04W72/543Allocation or scheduling criteria for wireless resources based on quality criteria based on requested quality, e.g. QoS

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Quality & Reliability (AREA)
  • Power Engineering (AREA)
  • Mobile Radio Communication Systems (AREA)
  • Radio Transmission System (AREA)

Abstract

本发明揭示了一种有限反馈欠秩信道时MU‑MIMO系统的资源分配方法,包括如下步骤:S1、建立MU‑MIMO系统模型,分析处于富散射环境和非富散射环境下的用户容量上限和等效信道状况;S2、推导有限反馈欠秩情况下的预编码矩阵,推导用户容量与SINR之和间的关系;S3、以同时实现最大化能量效率和最大化有效吞吐量为目标、建立双目标优化问题模型;S4、求解双目标优化问题,通过功率分配和用户选择实现系统能量效率和有效吞吐量的双目标优化。本发明的方法能够在保证精细QoS要求的前提下,同时实现对系统吞吐量和能量效率的优化。

Description

一种有限反馈欠秩信道时MU-MIMO系统的资源分配方法
技术领域
本发明涉及一种针对多目标的资源分配方法,具体涉及一种基于有限反馈下的信道欠秩时MU-MIMO系统的多目标优化资源分配方法,属于无线通信技术领域。
背景技术
近年来,随着无线通信技术的不断发展,各种相关技术层出不穷。以多用户多输入多输出(Multi-User Multiple Input Multiple Output,MU-MIMO)技术为例,该技术在用户端和基站端使用多根天线,利用空间资源,可极大地提高通信系统的频谱效率。在其应用过程中,由于多用户在同一时隙内使用相同频率,所以需要利用信道状态信息(ChannelState Information,CSI)在基站端进行预编码,从而消除用户间的共信道干扰(Co-Channel Interference)。
在传统的时分双工(TDD)系统中,用户端可以通过信道估计直接获知下行信道的CSI,基站端则可以通过信道互易原则获知下行信道CSI。但由于TDD系统要求全网同步、存在一定的技术难度,因此在实际的无线网络中,频分双工(FDD)系统的应用更为普遍。
在FDD系统中,基站端可以通过用户端的信息反馈间接获知下行信道CSI,但是在一些存在有大量用户和大量收发天线的场景中,信道矩阵的规模很大,有限的上行信道频谱资源负担不起庞大的反馈量。为了降低反馈量,可在用户端和基站端预存相同的码本,码本反映的是信道信息所有可能的取值,每种具体的取值称为码字,用户端只需向基站端反馈所选择的码字索引即可。但在实际的应用中,由于码字和真实的信道信息间存在一定的误差,会导致残留多用户干扰,这些干扰将会影响用户容量,进而对用户选择造成影响。
MU-MIMO技术可以在不消耗额外频谱带宽的情况下有效地利用空间资源去提高无线通信系统的吞吐量,因此MU-MIMO技术已经成为5G网络的关键技术之一。随着人们对无线通信需求以及社会对节能环保要求的不断提高,优化系统吞吐量以及能量效率(每单位能量所能传输的数据速率)是研究MU-MIMO资源分配的两个重要目标。资源分配算法通过用户选择和功率分配来优化吞吐量和能量,用户选择和功率分配与被选择用户的MIMO信道矩阵的秩相关,信道的秩决定了空分信道的个数,从而影响一定功率下用户可获得速率。
在现有的一些研究及相关文献中,进行资源分配时仅考虑了信道满秩的情况、即用户信道矩阵和所有用户的联合信道矩阵都是满秩的。然而在实际的应用环境中,由于散射情况和物理天线的分布,很多情况下满秩信道和欠秩信道是同时存在的,例如基站天线和用户终端天线存在视线路径(LOS)时。而在欠秩信道存在的情况下,用户可获得速率和系统吞吐量会随着秩的不同而变化、系统总功率也随之变化,而且系统能选择的最大用户数会随着被选择的用户集的不同而变化。因此满秩情况下所设计的资源分配算法并不能实现欠秩情况下系统吞吐量和能量效率的优化。
此外,虽然优化系统吞吐量以及优化能量效率一直是MU-MIMO系统资源分配的研究热点,但是,目前的相关研究和文献都只考虑了单一目标优化而没有考虑到两者的联合优化、即限定系统最低吞吐量去实现系统功率最小化或者限定系统最大功率去实现系统吞吐量最大化,仅仅能够实现对能量效率的优化。这是因为系统吞吐量和能量效率相互耦合,传统的贪婪算法无法同时优化这两个目标,而穷举法又因为计算量太大而无法实现。
在进行资源分配时,必须考虑到业务的QoS要求,包括时延要求和速率要求。由于不同种类的业务有不同的QoS要求,但在研究资源分配时,QoS要求过多可能会导致约束条件过多或者优化问题非凸的情况,因此大多数现有研究都没有对业务的QoS要求提供精细的保证。这些现有研究大多只考虑业务QoS速率要求的上限而没有考虑下限。实际上,用户业务可以分为实时业务和非实时业务,大多数实时业务不需要过大的速率,例如语音业务,64kbps就可使得话音信号质量几乎达到上限、MOS(Mean Opinion Score,平均意见值)4.4分。所以对实时业务来说,超过上限的速率对业务质量来说是无意义的,也就是说是无效的。
综上所述,针对上述各种问题,本发明提出了一种在真实环境下(存在欠秩信道),基于MU-MIMO有限反馈系统,充分利用所有等效信道资源,精细保证QoS要求(包括速率上下限和时延)并同时优化系统有效吞吐量和能量效率的资源分配算法。
发明内容
鉴于现有技术存在上述缺陷,本发明的目的是提出一种基于有限反馈下的信道欠秩时MU-MIMO系统的多目标优化资源分配方法,具体如下。
一种有限反馈欠秩信道时MU-MIMO系统的资源分配方法,包括如下步骤:
S1、建立MU-MIMO系统模型,分析处于富散射环境和非富散射环境下的用户容量上限和等效信道状况;
S2、推导有限反馈欠秩情况下的预编码矩阵,推导用户容量与SINR之间的关系;
S3、以同时实现最大化能量效率和最大化有效吞吐量为目标、建立双目标优化问题模型;
S4、求解双目标优化问题,通过功率分配和用户选择实现系统能量效率和有效吞吐量的双目标优化。
优选地,所述S1包括如下步骤:
建立有限反馈欠秩情况下MU-MIMO系统的下行链路模型,确定分别处于富散射环境和非富散射环境下的用户信道,利用块对角化技术进行预编码处理时的用户容量限制,并得到用户选择过程中新用户加入对被服务用户的等效信道矩阵维度以及秩的影响。
优选地,所述S1包括如下步骤:
S11、在MU-MIMO系统中,设系统基站有N根发射天线,系统内总的用户数为M0,第k个用户终端的接收天线数为nk,1≤k≤M0,令M为基站能同时服务的最大用户数;
S12、建立MU-MIMO系统模型,表达式为
Figure BDA0002376874510000041
其中,k=1,2,3…,M0且k≠i,
Figure BDA0002376874510000042
表示用户k的信道矩阵,
Figure BDA0002376874510000043
表示用户k的预编码矩阵,
Figure BDA0002376874510000044
为用户k的加性高斯白噪声、其均值为零、方差为σ2,Lk表示Hk的秩;
S13、令
Figure BDA0002376874510000045
其中,
Figure BDA0002376874510000046
Figure BDA0002376874510000047
表示
Figure BDA0002376874510000048
的秩,
Figure BDA0002376874510000049
进行奇异值分解为
Figure BDA00023768745100000410
其中,
Figure BDA00023768745100000411
表示
Figure BDA00023768745100000412
的对角矩阵,
Figure BDA00023768745100000413
表示
Figure BDA00023768745100000414
的酉矩阵,
Figure BDA0002376874510000051
用户的预编码矩阵为
Figure BDA0002376874510000052
其中,
Figure BDA0002376874510000053
用于消除用户间干扰,Bk用于实现用户速率最大化,最大用户数M的约束式为
Figure BDA0002376874510000054
其中,
Figure BDA0002376874510000055
表示
Figure BDA0002376874510000056
的秩,随被选择用户数目的增加而增大,在满秩情况下,
Figure BDA0002376874510000057
在欠秩的情况下,
Figure BDA0002376874510000058
优选地,所述S2包括如下步骤:
MU-MIMO系统模型内的用户端从码本中选取出量化误差最小的码字并将其索引反馈给基站端,基站端依据索引从码本中取出对应码字,基站端利用对应码字设计用户的预编码矩阵推导用户信道矩阵奇异值分解以后所有等效信道的SINR,并推导用户容量与SINR之间的关系。
优选地,所述S2包括如下步骤:
S21、MU-MIMO系统模型内的用户端从码本中选取出量化误差最小的码字并将其索引反馈给基站端,所依据的量化误差最小的准则为
Figure BDA0002376874510000059
其中,Hk′为从码本中选取的用户k的码字,V′k为H′k通过SVD分解的右酉矩阵;
S22、设从码本中选取的用户k的码字为H′k
Figure BDA00023768745100000510
Figure BDA00023768745100000511
中,对
Figure BDA00023768745100000512
进行SVD分解,得到用户k的预编码矩阵
Figure BDA00023768745100000513
S23、基站端使用
Figure BDA00023768745100000514
对用户信号做预编码,得到
Figure BDA0002376874510000061
接收端利用MMSE接收机检测用户k接收的信号,根据MMSE接收机的原理推导可得用户k信道矩阵奇异值分解以后第j个等效信道上有用信号功率和干扰及噪声功率比值
Figure BDA0002376874510000062
Figure BDA0002376874510000063
其中,pk,j表示为用户k的第j个平行信道上分配的功率,[·]H为共轭变换矩阵,
Figure BDA0002376874510000064
表示第k个用户所取码字的右酉矩阵的第j列,
Figure BDA0002376874510000065
为用户k信道矩阵的第j个奇异值,
用户k的信道容量为
Figure BDA0002376874510000066
其中,Rk为用户k的速率,rk为H′k的秩。
优选地,所述S3包括如下步骤:
根据同时最大化能量效率和最大化有效吞吐量的双优化目标以及保证系统要求吞吐量下限、天线功率上限及被选择用户速率上下限的优化条件建立优化问题模型,该优化问题的解是由用户平行信道功率和已选择用户集合这两个变量组成的解向量。
优选地,所述S3包括如下步骤:
令min E(pk,j),
max C(pk,j),
`s.t.C(pk,j)≥C0
Figure BDA0002376874510000067
Rk0<Rk<Rk1,k=1,2,…,M,
Figure BDA0002376874510000071
其中,
Figure BDA0002376874510000072
Figure BDA0002376874510000073
Rk1与Rk0分别表示用户k的速率上下限,Rk为用户k的速率,用户调度选中的用户集合为Ψ,基站能同时服务的最大用户数为M;
系统有效吞吐量为所有被选择用户的速率之和表示为C;
Figure BDA0002376874510000074
表示天线i的发射功率,i=1,2.3…,N;
Figure BDA0002376874510000075
P0表示天线功率的上限,E为基站总功率;
e表示功率放大器漏极效率的倒数,Pc代表电路功耗。
优选地,所述S4包括如下步骤:
建立S3中优化问题的拉格朗日对偶凸优化问题,将拉格朗日对偶凸优化问题进行展开化简得到一个与对偶问题具有相同最优解向量的简化凸优化问题,求解简化凸优化问题、求得对偶问题的最优解向量,通过凸优化技术解向量中的用户平行信道功率的最优解关于已选择用户集合的函数,将求得的用户平行信道功率关于已选择用户集合的最优解代入简化凸优化问题目标函数,通过用户选择来求得最优解。
优选地,所述S4包括如下步骤:
S41、建立S3中优化问题的拉格朗日对偶凸优化问题,将拉格朗日对偶凸优化问题进行展开化简得到一个与对偶问题具有相同最优解向量的简化凸优化问题,即
采用主要目标法,得到最终优化问题为
min E(pk,j),
C(pk,j)≥C0
Figure BDA0002376874510000081
Rk0<Rk<Rk1,k=1,2,…,M,
Figure BDA0002376874510000082
上述问题的拉格朗日对偶问题为
max h(u,α,β,γ),pk,j,α,β,γ,μ>0,
Figure BDA0002376874510000083
i,j,
其中,
Figure BDA0002376874510000084
Figure BDA0002376874510000085
Figure BDA0002376874510000086
α=(α12,…,αM),β=(β12,…,βM),γ=(γ12,…,γN)为拉格朗日乘子向量,μ为拉格朗日乘子;
S42、求得对偶问题的最优解向量,即
根据拉格朗日对偶算法的性质,不论原问题是否为凸优化问题,其对偶问题均为凸优化问题,对于凸优化问题,KKT点就是最优点;
假设最优拉格朗日乘子为u****,最优化点
Figure BDA0002376874510000087
满足以下必要条件
Figure BDA0002376874510000088
Figure BDA0002376874510000089
S43、通过用户选择来求得最优值,即
根据业务的时延要求,业务z的时延要求为dz,则该业务的最大等待时隙数
Figure BDA00023768745100000810
tti代表时隙长度,令Wk,z代表使用业务z的用户k的等待时隙数;在进行用户选择时,首先选择等待时隙数达到最大时延的用户,其次以最大化吞吐量为标准进行用户选择;且乘子α,β,γ,μ的最优值通过二分算法来计算,随后将用户选择带入二分算法。
优选地,所述S43中所述二分算法,包括如下步骤:
S431、初始化,μmin=0,αmin=βmin=γmin=0,μmax>>0,αmaxmaxmax,>>0;
S432、
Figure BDA0002376874510000091
Figure BDA0002376874510000092
Figure BDA0002376874510000093
Figure BDA0002376874510000094
S433、用改进的贪婪算法进行用户选择;
S434、将本次迭代得到的u,α,β,γ代入最优解中,计算pk,j
S435、利用所得到
Figure BDA0002376874510000095
的和pk,j计算出∑k∈ΨRk
若∑k∈ΨRk≥C0,则令μmax=μ,否则令μmin=μ;
S436、重复步骤S432~S435,直到μmaxmin≤δ,
其中,δ为控制精度,δ越小,算法精度越高,计算量越大。
与现有技术相比,本发明的优点主要体现在以下几个方面:
本发明提出了一种基于有限反馈下的信道欠秩时MU-MIMO系统的多目标优化资源分配方法,充分地利用了MIMO空间信道资源,并给出了反映用户容量的综合SINR值。本发明的方法能够在保证精细QoS要求的前提下,同时实现对系统吞吐量和能量效率的优化。
此外,本发明的方法还为同领域内的其他相关问题提供了参考,具有十分广阔的应用前景,可以以此为依据进行拓展延伸,将其应用于无线通信技术和频谱资源分配领域内。
以下便结合实施例附图,对本发明的具体实施方式作进一步的详述,以使本发明技术方案更易于理解、掌握。
附图说明
图1为本发明的方法流程示意图;
图2为本发明有限反馈欠秩情况下的MU-MIMO系统下行信道模型示意图;
图3为四种类型业务的时延和速率要求情况表。
具体实施方式
本发明希望设计出一种在有限反馈欠秩信道情况下、联合优化系统有效吞吐量和能量效率的资源分配算法。本发明的方案在设计双优化算法时,根据同时最大化能量效率和最大化有效吞吐量的双优化目标以及保证系统要求吞吐量下限、天线功率上限及被选择用户速率上下限的优化条件建立优化问题模型并求解,该优化问题的解是由用户平行信道功率和已选择用户集合这两个变量组成的解向量。此优化问题为双目标优化且约束条件众多,并不是一个标准的凸优化问题,不能直接用凸优化技术来解决,因此本发明通过拉格朗日对偶算法来解决该优化问题。
如图1~图3所示,本发明提出了一种信道欠秩时MU-MIMO系统的多目标优化资源分配方法,具体包括如下步骤。
S1、建立MU-MIMO系统模型,分析处于富散射环境和非富散射环境下的用户容量上限和等效信道状况;
S2、推导有限反馈欠秩情况下的预编码矩阵,推导用户容量与SINR之间的关系;
S3、以同时实现最大化能量效率和最大化有效吞吐量为目标、建立双目标优化问题模型;
S4、求解双目标优化问题,通过功率分配和用户选择实现系统能量效率和有效吞吐量的双目标优化。
具体而言,所述S1包括如下步骤:
S11、建立有限反馈欠秩情况下MU-MIMO系统的下行链路模型,在MU-MIMO系统中,设系统基站(BS)有N根发射天线,系统内总的用户数为M0,第k个用户终端的接收天线数为nk,1≤k≤M0,令M为基站能同时服务的最大用户数,具体公式为
Figure BDA0002376874510000111
其中公式中的各参数含义为:k=1,2.3…,M0;k≠i,
Figure BDA0002376874510000112
表示用户k的信道矩阵,
Figure BDA0002376874510000113
表示用户k的预编码矩阵,
Figure BDA0002376874510000114
为用户k的加性高斯白噪声,其均值为零,方差为σ2。用Lk来表示Hk的秩。
S12、依据MU-MIMO系统模型,分析分别处于富散射环境和非富散射环境下的用户信道,利用块对角化技术进行预编码处理时的用户容量限制,并得到用户选择过程中新用户加入对被服务用户的等效信道矩阵维度以及秩的影响,具体为;
Figure BDA0002376874510000115
Figure BDA0002376874510000116
Figure BDA0002376874510000117
表示
Figure BDA0002376874510000118
的秩。对
Figure BDA0002376874510000119
进行奇异值分解:
Figure BDA00023768745100001110
其中,
Figure BDA00023768745100001111
表示
Figure BDA00023768745100001112
的对角矩阵,
Figure BDA00023768745100001113
表示
Figure BDA00023768745100001114
的酉矩阵。那么根据奇异值分解的性质,可以得出:
Figure BDA00023768745100001115
用户的预编码矩阵可以设计为:
Figure BDA0002376874510000121
其中,
Figure BDA0002376874510000122
用来消除用户间干扰,Bk用来实现用户速率最大化。
为了使(4)中
Figure BDA0002376874510000123
有非零解,那么等式的个数应该小于变量的个数,因此:
Figure BDA0002376874510000124
其中
Figure BDA0002376874510000125
表示
Figure BDA0002376874510000126
的秩,它随着被选择用户数目的增加而增大,那么(5)就是最大用户数M的约束式。
在满秩情况下,
Figure BDA0002376874510000127
而在欠秩的情况下,
Figure BDA0002376874510000128
因此相比于满秩情况,欠秩情况下,最大用户数M会增大,即系统可以选择更多的用户同时服务。
具体而言,所述S2包括如下步骤:
S21、MU-MIMO系统模型内的用户端从码本中选取出量化误差最小的码字将其索引反馈给基站端,具体为,与信道矩阵弦距(chordal distance)最小的码字量化误差最小,而码字和信道矩阵的弦距可以用其乘积的模值衡量,模值越大弦距越小,依据量化误差最小的准则:
Figure BDA0002376874510000129
S22、基站端依据索引从码本中取出对应码字,基站端利用对应码字设计用户的预编码矩阵推导用户信道矩阵奇异值分解以后所有等效信道的SINR,并推导用户容量与SINR之和间的关系,分别为:
假设从码本中选取的用户k的码字为H′k,基站端获知所有用户的码字以后,将除了用户k码字H′k的其他所有码字带入
Figure BDA00023768745100001210
Figure BDA00023768745100001211
中,再对
Figure BDA00023768745100001212
进行SVD分解,就可以求出用户k的预编码矩阵:
Figure BDA00023768745100001213
基站端利用
Figure BDA0002376874510000131
对用户信号做预编码,将
Figure BDA0002376874510000132
带入(1)式可得:
Figure BDA0002376874510000133
但由于Vi和V′i存在误差,导致求出的预编码矩阵
Figure BDA0002376874510000134
无法完全消除用户的相互干扰,故存在干扰残留
Figure BDA0002376874510000135
在接收端,为了降低多用户干扰以及信道估计误差对信号检测的影响,可以利用MMSE接收机检测用户k接收的信号,检测到的估计信号可以表示为:
Figure BDA0002376874510000136
MMSE接收机的原理是使估计信号和原始信号误差的平方根最小,即
minE{||s′k-sk||2}=E{(s′k-sk)H(s′k-sk)}, (10)
其中,s′k为估计信号,sk为原始信号。
将s′k=gkyk带入上式并进一步变形可得:
Figure BDA0002376874510000137
其中,tr为矩阵的迹。
再对gk求偏导并令其等于零可得:
Figure BDA0002376874510000138
Figure BDA0002376874510000139
其中,[·]H为共轭变换矩阵,[·]-1为逆矩阵,IN为N×N的单位矩阵,
Figure BDA00023768745100001310
再将
Figure BDA00023768745100001311
带入上式可得:
Figure BDA00023768745100001312
Figure BDA0002376874510000141
其中
Figure BDA0002376874510000142
为用户k信道矩阵的第j个奇异值,rk为H′k的秩。
再把(9)式带入(5)式可得:
Figure BDA0002376874510000143
由上式推得用户k信道矩阵奇异值分解以后第j个等效信道上有用信号功率和干扰及噪声功率比值
Figure BDA0002376874510000144
Figure BDA0002376874510000145
用户k的信道容量为:
Figure BDA0002376874510000146
具体而言,所述S3包括如下步骤:
S31、根据同时最大化能量效率和最大化有效吞吐量的双优化目标以及保证系统要求吞吐量下限、天线功率上限及被选择用户速率上下限的优化条件建立优化问题模型,具体为
min E(pk,j), (18)
max C(pk,j), (19)
`s.t.C(pk,j)≥C0
Figure BDA0002376874510000151
Rk0<Rk<Rk1,k=1,2,…,M
Figure BDA0002376874510000152
其中,
Figure BDA0002376874510000153
Figure BDA0002376874510000154
与Rk0分别表示用户k的速率上下限,Rk为用户k的速率,用户调度选中的用户集合为Ψ,基站能同时服务的最大用户数为M。系统有效吞吐量为所有被选择用户的速率之和表示为C。
Figure BDA0002376874510000155
表示天线i的发射功率,i=1,2.3…,N:
Figure BDA0002376874510000156
P0表示天线功率的上限,E为基站总功率。e表示功率放大器漏极效率的倒数,Pc代表电路功耗。
具体而言,所述S4包括如下步骤:
S41、建立S3中优化问题的拉格朗日对偶凸优化问题,将拉格朗日对偶凸优化问题进行展开化简得到一个与对偶问题具有相同最优解向量的简化凸优化问题,分别为:
采用主要目标法,最终优化问题为:
min E(pk,j), (20)
C(pk,j)≥C0
Figure BDA0002376874510000157
Rk0<Rk<Rk1,k=1,2,…,M
Figure BDA0002376874510000161
针对优化问题(20),本发明通过拉格朗日对偶算法来解决该优化问题。首先构建上述优化问题的拉格朗日函数:
Figure BDA0002376874510000162
其中α=(α12,…,αM),β=(β12,…,βM),γ=(γ12,…,γN)为拉格朗日乘子向量,μ为拉格朗日乘子,而且pk,j,α,β,γ,μ>0,
Figure BDA0002376874510000163
i,j。
相应的拉格朗日对偶函数为:
h(u,α,β,γ)=min L(pk,j,α,β,γ,μ), (22)
原问题的拉格朗日对偶问题为:
max h(u,α,β,γ), (23)
Figure BDA0002376874510000164
S42、求得对偶问题的最优解向量,具体为,根据拉格朗日对偶算法的性质,不论原问题是否为凸优化问题,其对偶问题均为凸优化问题。对于凸优化问题,KKT点就是最优点。假设最优拉格朗日乘子为u****,最优化点
Figure BDA0002376874510000165
满足以下必要条件:
Figure BDA0002376874510000166
Figure BDA0002376874510000167
S43、通过用户选择来求得最优值,具体为,根据图3中业务的时延要求,业务z的时延要求为dz,则该业务的最大等待时隙数
Figure BDA0002376874510000168
tti代表时隙长度,令Wk,z代表使用业务z的用户k的等待时隙数。那么在进行用户选择时,为了保证业务的时延要求,首先选择等待时隙数达到最大时延的用户,其次以最大化吞吐量为标准进行用户选择。而且(23)中的乘子α,β,γ,μ的最优值可以通过二分法来计算。
所述二分算法的具体步骤如下:
S431、μmin=0,αmin=βmin=γmin=0,μmax>>0,αmaxmaxmax,>>0。
S432、
Figure BDA0002376874510000171
Figure BDA0002376874510000172
Figure BDA0002376874510000173
Figure BDA0002376874510000174
S433、用改进的贪婪算法进行用户选择。
S434、将本次迭代得到的u,α,β,γ代入公式(24),计算pk,j
S435、将S432、S433中所得到
Figure BDA0002376874510000179
的和pk,j代入公式(20),计算出∑k∈ΨRk。如果∑k∈ΨRk≥C0,则令μmax=μ,否则令μmin=μ。
S436、重复S432~S435,直到μmaxmin≤δ。
将用户选择带入二分算法的第二步,用户选择的计算方法如下:
初始化:未调用用户集Ω={1,2,…,M0},已调用用户集
Figure BDA0002376874510000175
步骤一:计算Ω中所有用户的Wk,z,如果Wk,z=nz,调用用户k。更新Ψ={k:Wk,z≥nz-1},Ω=Ω+Ψ。
步骤二:计算用户k,k∈Ψ的
Figure BDA0002376874510000176
判决
Figure BDA0002376874510000177
k∈Ψ是否成立,如果不成立,结束调度。若成立,首先调度用户容量最大的用户作为k1,满足
Figure BDA0002376874510000178
更新Ψ=Ψ+{k1},Ω=Ω-{k1}。
步骤三:计算用户k,k∈Ψ的
Figure BDA0002376874510000181
判决
Figure BDA0002376874510000182
k∈Ψ是否成立,如果不成立,结束调度。如果成立,对
Figure BDA0002376874510000183
令Ψs=Ψ+{s},并计算集合Ψs的容量:
Figure BDA0002376874510000184
在集合Ω中,找到一个用户使得
Figure BDA0002376874510000185
即调度用户s满足:
Figure BDA0002376874510000186
如果对于
Figure BDA0002376874510000187
均有
Figure BDA0002376874510000188
则结束算法。更新系统容量
Figure BDA0002376874510000189
选中用户集合Ψ=Ψ+{s},Ω=Ω-{s}。
步骤四:重复步骤三。
最终,调度用户集为Ψ。
本发明推导了一般真实环境下、即欠秩信道存在的情况下,基于有限反馈的MU-MIMO系统能同时服务的最大用户数、用户可获得速率以及系统总功率,证明了欠秩情况下系统能选择的用户数大于满秩情况下系统能选择的用户数、即欠秩情况下系统可以选择更多的用户同时服务。
综上所述,本发明提出了一种基于有限反馈下的信道欠秩时MU-MIMO系统的多目标优化资源分配方法,充分地利用了MIMO空间信道资源,并给出了反映用户容量的综合SINR值。本发明的方法能够在保证精细QoS要求的前提下,同时实现对系统吞吐量和能量效率的优化。
此外,本发明的方法还为同领域内的其他相关问题提供了参考,具有十分广阔的应用前景,可以以此为依据进行拓展延伸,将其应用于无线通信技术和频谱资源分配领域内。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神和基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内,不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。

Claims (10)

1.一种有限反馈欠秩信道时MU-MIMO系统的资源分配方法,其特征在于,包括如下步骤:
S1、建立MU-MIMO系统模型,分析处于富散射环境和非富散射环境下的用户容量上限和等效信道状况;
S2、推导有限反馈欠秩情况下的预编码矩阵,推导用户容量与SINR之间的关系;
S3、以同时实现最大化能量效率和最大化有效吞吐量为目标、建立双目标优化问题模型;
S4、求解双目标优化问题,通过功率分配和用户选择实现系统能量效率和有效吞吐量的双目标优化。
2.根据权利要求1所述的一种有限反馈欠秩信道时MU-MIMO系统的资源分配方法,其特征在于,所述S1包括如下步骤:
建立有限反馈欠秩情况下MU-MIMO系统的下行链路模型,确定分别处于富散射环境和非富散射环境下的用户信道,利用块对角化技术进行预编码处理时的用户容量限制,并得到用户选择过程中新用户加入对被服务用户的等效信道矩阵维度以及秩的影响。
3.根据权利要求1所述的一种有限反馈欠秩信道时MU-MIMO系统的资源分配方法,其特征在于,所述S1包括如下步骤:
S11、在MU-MIMO系统中,设系统基站有N根发射天线,系统内总的用户数为M0,第k个用户终端的接收天线数为nk,1≤k≤M0,令M为基站能同时服务的最大用户数;
S12、建立MU-MIMO系统模型,表达式为
Figure FDA0002376874500000021
其中,k=1,2,3...,M0且k≠i,
Figure FDA0002376874500000022
表示用户k的信道矩阵,
Figure FDA0002376874500000023
表示用户k的预编码矩阵,
Figure FDA0002376874500000024
为用户k的加性高斯白噪声、其均值为零、方差为σ2,Lk表示Hk的秩;
S13、令
Figure FDA0002376874500000025
其中,
Figure FDA0002376874500000026
Figure FDA0002376874500000027
表示
Figure FDA0002376874500000028
的秩,
Figure FDA0002376874500000029
进行奇异值分解为
Figure FDA00023768745000000210
其中,
Figure FDA00023768745000000211
表示
Figure FDA00023768745000000212
的对角矩阵,
Figure FDA00023768745000000213
表示
Figure FDA00023768745000000214
的酉矩阵,
Figure FDA00023768745000000215
用户的预编码矩阵为
Figure FDA00023768745000000216
其中,
Figure FDA00023768745000000217
用于消除用户间干扰,Bk用于实现用户速率最大化,最大用户数M的约束式为
Figure FDA00023768745000000218
其中,
Figure FDA00023768745000000219
表示
Figure FDA00023768745000000220
的秩,随被选择用户数目的增加而增大,在满秩情况下,
Figure FDA00023768745000000221
在欠秩的情况下,
Figure FDA00023768745000000222
4.根据权利要求1所述的一种有限反馈欠秩信道时MU-MIMO系统的资源分配方法,其特征在于,所述S2包括如下步骤:
MU-MIMO系统模型内的用户端从码本中选取出量化误差最小的码字并将其索引反馈给基站端,基站端依据索引从码本中取出对应码字,基站端利用对应码字设计用户的预编码矩阵推导用户信道矩阵奇异值分解以后所有等效信道的SINR,并推导用户容量与SINR之间的关系。
5.根据权利要求3所述的一种有限反馈欠秩信道时MU-MIMO系统的资源分配方法,其特征在于,所述S2包括如下步骤:
S21、MU-MIMO系统模型内的用户端从码本中选取出量化误差最小的码字并将其索引反馈给基站端,所依据的量化误差最小的准则为
Figure FDA0002376874500000031
其中,Hk为从码本中选取的用户k的码字,Vk为Hk通过SVD分解的右酉矩阵;
S22、设从码本中选取的用户k的码字为H′k
Figure FDA0002376874500000032
Figure FDA0002376874500000033
中,对
Figure FDA0002376874500000034
进行SVD分解,得到用户k的预编码矩阵
Figure FDA0002376874500000035
S23、基站端使用
Figure FDA0002376874500000036
对用户信号做预编码,得到
Figure FDA0002376874500000037
接收端利用MMSE接收机检测用户k接收的信号,根据MMSE接收机的原理推导可得用户k信道矩阵奇异值分解以后第j个等效信道上有用信号功率和干扰及噪声功率比值
Figure FDA00023768745000000311
Figure FDA0002376874500000038
其中,pk,j表示为用户k的第j个平行信道上分配的功率,[·]H为共轭变换矩阵,
Figure FDA0002376874500000039
表示第k个用户所取码字的右酉矩阵的第j列,
Figure FDA00023768745000000312
为用户k信道矩阵的第j个奇异值,
用户k的信道容量为
Figure FDA00023768745000000310
其中,Rk为用户k的速率,rk为H′k的秩。
6.根据权利要求1所述的一种有限反馈欠秩信道时MU-MIMO系统的资源分配方法,其特征在于,所述S3包括如下步骤:
根据同时最大化能量效率和最大化有效吞吐量的双优化目标以及保证系统要求吞吐量下限、天线功率上限及被选择用户速率上下限的优化条件建立优化问题模型,该优化问题的解是由用户平行信道功率和已选择用户集合这两个变量组成的解向量。
7.根据权利要求5所述的一种有限反馈欠秩信道时MU-MIMO系统的资源分配方法,其特征在于,所述S3包括如下步骤:
令min E(pk,j),
max C(pk,j),
`s.t.C(pk,j)≥C0
Figure FDA0002376874500000041
Rk0<Rk<Rk1,k=1,2,…,M,
Figure FDA0002376874500000042
其中,
Figure FDA0002376874500000043
Figure FDA0002376874500000044
Rk1与Rk0分别表示用户k的速率上下限,Rk为用户k的速率,用户调度选中的用户集合为Ψ,基站能同时服务的最大用户数为M;
系统有效吞吐量为所有被选择用户的速率之和表示为C;
Figure FDA0002376874500000045
表示天线i的发射功率,i=1,2.3…,N;
Figure FDA0002376874500000046
P0表示天线功率的上限,E为基站总功率;
e表示功率放大器漏极效率的倒数,Pc代表电路功耗。
8.根据权利要求1所述的一种有限反馈欠秩信道时MU-MIMO系统的资源分配方法,其特征在于,所述S4包括如下步骤:
建立S3中优化问题的拉格朗日对偶凸优化问题,将拉格朗日对偶凸优化问题进行展开化简得到一个与对偶问题具有相同最优解向量的简化凸优化问题,求解简化凸优化问题、求得对偶问题的最优解向量,通过凸优化技术解向量中的用户平行信道功率的最优解关于已选择用户集合的函数,将求得的用户平行信道功率关于已选择用户集合的最优解代入简化凸优化问题目标函数,通过用户选择来求得最优解。
9.根据权利要求7所述的一种有限反馈欠秩信道时MU-MIMO系统的资源分配方法,其特征在于,所述S4包括如下步骤:
S41、建立S3中优化问题的拉格朗日对偶凸优化问题,将拉格朗日对偶凸优化问题进行展开化简得到一个与对偶问题具有相同最优解向量的简化凸优化问题,即
采用主要目标法,得到最终优化问题为
min E(pk,j),
C(pk,j)≥C0
Figure FDA0002376874500000051
Rk0<Rk<Rk1,k=1,2,…,M,
Figure FDA0002376874500000052
上述问题的拉格朗日对偶问题为
Figure FDA0002376874500000053
其中,
Figure FDA0002376874500000054
Figure FDA0002376874500000061
Figure FDA0002376874500000062
α=(α1,α2,…,αM),β=(β1,β2,…,βM),γ=(γ1,γ2,…,γN)为拉格朗日乘子向量,μ为拉格朗日乘子;
S42、求得对偶问题的最优解向量,即
根据拉格朗日对偶算法的性质,不论原问题是否为凸优化问题,其对偶问题均为凸优化问题,对于凸优化问题,KKT点就是最优点;
假设最优拉格朗日乘子为u*,α*,β*,γ*,最优化点
Figure FDA0002376874500000063
满足以下必要条件
Figure FDA0002376874500000064
Figure FDA0002376874500000065
S43、通过用户选择来求得最优值,即
根据业务的时延要求,业务z的时延要求为dz,则该业务的最大等待时隙数
Figure FDA0002376874500000066
tti代表时隙长度,令Wk,z代表使用业务z的用户k的等待时隙数;在进行用户选择时,首先选择等待时隙数达到最大时延的用户,其次以最大化吞吐量为标准进行用户选择;且乘子α,β,γ,μ的最优值通过二分算法来计算,随后将用户选择带入二分算法。
10.根据权利要求9所述的一种有限反馈欠秩信道时MU-MIMO系统的资源分配方法,其特征在于,所述S43中所述二分算法,包括如下步骤:
S431、初始化,μmin=0,αmin=βmin=γmin=0,μmax>>0,αmax,βmax,γmax,>>0;
S432、
Figure FDA0002376874500000067
Figure FDA0002376874500000068
Figure FDA0002376874500000071
Figure FDA0002376874500000072
S433、用改进的贪婪算法进行用户选择;
S434、将本次迭代得到的u,α,β,γ代入最优解中,计算pk,j
S435、利用所得到
Figure FDA0002376874500000073
的和pk,j计算出∑k∈ΨRk
若∑k∈ΨRk≥C0,则令μmax=μ,否则令μmin=μ;
S436、重复步骤S432~S435,直到μmaxmin≤δ,
其中,δ为控制精度,δ越小,算法精度越高,计算量越大。
CN202010069963.2A 2020-01-21 2020-01-21 一种有限反馈欠秩信道时mu-mimo系统的资源分配方法 Active CN111277307B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010069963.2A CN111277307B (zh) 2020-01-21 2020-01-21 一种有限反馈欠秩信道时mu-mimo系统的资源分配方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010069963.2A CN111277307B (zh) 2020-01-21 2020-01-21 一种有限反馈欠秩信道时mu-mimo系统的资源分配方法

Publications (2)

Publication Number Publication Date
CN111277307A true CN111277307A (zh) 2020-06-12
CN111277307B CN111277307B (zh) 2022-04-08

Family

ID=71001967

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010069963.2A Active CN111277307B (zh) 2020-01-21 2020-01-21 一种有限反馈欠秩信道时mu-mimo系统的资源分配方法

Country Status (1)

Country Link
CN (1) CN111277307B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115987340A (zh) * 2023-03-21 2023-04-18 南京邮电大学 一种5g物联网信道相干且有限反馈条件下的用户调度方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106209191A (zh) * 2016-07-20 2016-12-07 南京邮电大学 一种mu‑mimo系统真实环境低复杂度用户选择方法
CN108093411A (zh) * 2018-01-10 2018-05-29 重庆邮电大学 D2d通信网络中基于信道签名的资源调度优化方法
WO2018205259A1 (en) * 2017-05-12 2018-11-15 Qualcomm Incorporated Precoder resource group allocation methods for mimo communication
CN108832979A (zh) * 2018-06-11 2018-11-16 南京邮电大学 一种信道欠秩时mu-mimo系统多目标优化资源分配算法
CN109617583A (zh) * 2018-12-14 2019-04-12 华南理工大学 一种预编码与功率分配联合优化方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106209191A (zh) * 2016-07-20 2016-12-07 南京邮电大学 一种mu‑mimo系统真实环境低复杂度用户选择方法
WO2018205259A1 (en) * 2017-05-12 2018-11-15 Qualcomm Incorporated Precoder resource group allocation methods for mimo communication
CN110612694A (zh) * 2017-05-12 2019-12-24 高通股份有限公司 用于mimo通信的预编码器资源组分配方法
CN108093411A (zh) * 2018-01-10 2018-05-29 重庆邮电大学 D2d通信网络中基于信道签名的资源调度优化方法
CN108832979A (zh) * 2018-06-11 2018-11-16 南京邮电大学 一种信道欠秩时mu-mimo系统多目标优化资源分配算法
CN109617583A (zh) * 2018-12-14 2019-04-12 华南理工大学 一种预编码与功率分配联合优化方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
SU PAN; YAN YAN; KUSI ANKRAH BONSU; WEIWEI ZHOU: "Resource Allocation Algorithm for MU-MIMO Systems With Double-Objective Optimization Under the Existence of the Rank Deficient Channel Matrix", 《 IEEE ACCESS ( VOLUME: 7)》 *
潘甦,张哓波: "SCMA系统中基于等效频谱带宽的无线资源分配", 《南京邮电大学学报(自然科学版)》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115987340A (zh) * 2023-03-21 2023-04-18 南京邮电大学 一种5g物联网信道相干且有限反馈条件下的用户调度方法

Also Published As

Publication number Publication date
CN111277307B (zh) 2022-04-08

Similar Documents

Publication Publication Date Title
US8464105B2 (en) Method of multiuser precoding and scheduling and base station for implementing the same
US20130163517A1 (en) Methods for optimal collaborative mimo-sdma
CN106060950B (zh) 一种基于机会干扰对齐的蜂窝下行信道中数据传输方法
JP5244172B2 (ja) 無線通信ネットワークにおける複数アンテナ送信を制御する方法及び装置
CN110492915A (zh) 一种基于mimo-noma短包传输的功率分配方法
CN106452525B (zh) 一种大规模MIMO中能量效率、频谱效率和QoS联合优化方法
EP3185434B1 (en) Method and device for beamforming
CN108832979B (zh) 一种信道欠秩时mu-mimo系统多目标优化资源分配算法
CN111970033B (zh) 能效谱效联合优化的大规模mimo多播功率分配方法
CN103117787A (zh) 一种协同多天线系统中的自适应比特分配方法及装置
CN102882570A (zh) 移动通信网络下设备间通信的最优收发联合处理方法
CN102811491A (zh) 多点协作系统中有限反馈比特数联合分配方法
CN102104451A (zh) 多输入多输出系统中多用户收发联合预编码的方法及装置
CN114337976A (zh) 一种联合ap选择与导频分配的传输方法
CN107425891B (zh) Mu-mimo无线能量/数据传输中的能量效率联合优化方法
CN111277307B (zh) 一种有限反馈欠秩信道时mu-mimo系统的资源分配方法
CN108370263A (zh) 用于在mu-mimo系统中增强用户选择的方法和装置
Dong et al. Energy-efficient sparse beamforming for multiuser MIMO systems with nonideal power amplifiers
Li et al. IRS-Based MEC for Delay-Constrained QoS Over RF-Powered 6G Mobile Wireless Networks
CN108064070B (zh) 一种用于大规模mimo多小区网络的用户接入方法
Almasaoodi et al. New Quantum Strategy for MIMO System Optimization.
CN114710187A (zh) 一种用户数动态变换场景下多小区大规模mimo智能通信的功率分配方法
CN115065392A (zh) 脏纸编码条件下实现miso下行总和速率最大化的波束形成设计方法
CN108418615B (zh) 一种基于mu-mimo有限反馈系统的用户调度方法
CN102970106B (zh) 具有码字约束的码字选择以及相应的信道质量估计方法

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