CN106357318A - 收敛速率可调的大规模mimo迭代检测方法 - Google Patents
收敛速率可调的大规模mimo迭代检测方法 Download PDFInfo
- Publication number
- CN106357318A CN106357318A CN201610926372.6A CN201610926372A CN106357318A CN 106357318 A CN106357318 A CN 106357318A CN 201610926372 A CN201610926372 A CN 201610926372A CN 106357318 A CN106357318 A CN 106357318A
- Authority
- CN
- China
- Prior art keywords
- matrix
- detection method
- iterative
- iteration
- mmse
- 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
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
- H04B7/0837—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
- H04B7/0842—Weighted combining
- H04B7/0848—Joint weighting
- H04B7/0854—Joint weighting using error minimizing algorithms, e.g. minimum mean squared error [MMSE], "cross-correlation" or matrix inversion
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Radio Transmission System (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明公开了一种收敛速率可调的大规模MIMO迭代检测方法(MMSE‑ANSE),将信道矩阵和接收信号向量经过预处理后得到匹配滤波输出向量和MMSE滤波矩阵;通过调节输入收敛速率控制参数,所提出的迭代算法的迭代矩阵将发生变化并始终和MMSE滤波矩阵的下三角阵元素有关;设置初始迭代解为零向量,收敛速率控制参数越大,迭代次数越多,所提出的迭代检测算法性能越好;所述检测算法的算法复杂度O(LM2)。本发明在复杂性和性能方面优于常规的基于Neumann级数展开的检测方法,特别是在传播环境较差的情况下。
Description
技术领域
本发明涉及计算机通信领域,尤其涉及一种收敛速率可调的大规模MIMO上行链路信号迭代检测方法。
背景技术
大规模多输入多输出(Multiple-Input Multiple-Output,MIMO)被认为是第5代(5G)无线***的关键技术之一[1]。通过装备大量天线(例如,在基站处安装数百个并且用户端安装数十个),该技术可以提供更高的频谱效率,更快的峰值数据速率以及比小规模MIMO***上的更好的能量效率[2]。然而,在大规模MIMO上行链路中,随着天线数目的显着增加,诸如最大似然(ML)检测和最大后验(MAP)检测的最佳检测方法在计算复杂度方面变得难以承受[3]。因此研究人员将目光转向近似MAP算法(例如消息传递检测器[4])和近似ML检测算法(例如多分支[5]和可能性上升搜索(LAS)[6][7]检测器)。同时,线性检测方法,如迫零(ZF)和最小均方误差(MMSE),在大规模MIMO***中因其次优性能和低复杂性特性备受关注[1][2]。
不幸的是,ZF和MMSE检测方法都不可避免地涉及计算复杂度很高的矩阵求逆操作。基于Cholesky分解[8]的M×M矩阵的精确求逆算法对于大规模MIMO***来说需要O(M3)的高计算复杂度,其中M表示单天线用户的数量。因此,有人提出了一种基于Neumann系列扩展(NSE)的检测方法[9],以降低矩阵求逆的复杂度。这种基于NSE的检测方法具有易于硬件实现的优点,因为其算法具有内在并行性。然而,它在较差的传播环境,例如天线比N/M较小时(N是基站处的天线的数量)或考虑相关信道的情况下收敛十分缓慢或甚至不能收敛。更糟的是,如果所采用的Neumann级数的项数大于2,传统的基于NSE的方法仍然花费O(M3)的复杂度[9]。为了保持O(M2)的复杂度,在[10]中提出了一种基于NSE的迭代检测算法。然而,该方法只能实现和基于NSE的传统检测算法相同的BER性能,并没有解决在较差的传播环境下的收敛问题。
发明内容
发明目的:本发明针对现有技术存在的问题,提供一种收敛速率可调的大规模MIMO上行链路信号迭代检测方法(MMSE-ANSE),其通过引入一个收敛速率控制参数来调整迭代矩阵形式,使得所提出的方法的收敛速度变得可调。根据数值模拟结果,所提出的方法在复杂性和性能方面优于常规的基于NSE的方法,特别是在传播环境较差的情况下。
技术方案:本发明所述的收敛速率可调的大规模MIMO迭代检测方法,包括:
(1)将信道矩阵H和接收信号y通过匹配滤波器进行预处理,得到输出信号yMF和MMSE滤波矩阵W;其中,yMF=HHy,W=G+N0ΙM,G=HHH为Gram矩阵,N0为噪声方差,ΙM为M维单位阵,(·)H为共轭转置操作;
(2)根据MMSE滤波矩阵W和预设收敛速率控制参数k计算得到迭代矩阵X和Y;其中,X=tril(W)-tril(W,k),Y=X-W,tril(W)表示取W的下三角阵,tril(W,k)表示取W的第k个对角线以下的元素(为下三角形式),则X由矩阵W的第k对角线至主对角线元素构成;
(3)设置迭代初始解为s0=0;
(4)按照公式sl=X-1(Ysl-1+yMF)进行迭代计算;其中,sl表示第l次迭代时s的值,X-1表示X的逆,由于X元素分布稀疏切都在下三角部分,该求逆计算复杂度小于等于O(M2);
(5)达到预设迭代次数后,迭代结束,所计算得到s的值为待检测信号的估计结果。
进一步的,步骤(2)中的控制参数k的取值范围为1~M次,M发射天线数。
进一步的,步骤(4)步骤4中的迭代次数l为1~L,L为设置的最大迭代次数。因此检测算法的所需的复数乘法次数为(LM2+M),L通常较小,因此算法复杂度为O(LM2)。
采用上述步骤计算的原理具体为:
考虑到大规模MIMO***上行链路MMSE检测中滤波矩阵W为Hermitian正定阵且主对角线占优,因此有以下定理:
对于大规模MIMO***上行链路,迭代方法sl=X-1(Ysl-1+yMF)对所有初始解s0收敛。
证明过程如下:令λ为迭代矩阵X-1Y的主特征值。令η为λ对应的特征向量。令m为η的最大模值的下标。缩放η使得|ηm|=1,且|ηi|≤1,对于i≠m。将X-1Yη=λη的第m行改写为
满足不等式
因为在大规模MIMO***中W主对角线占优,所以上面式子的最后一项可以表示为σ2/(d-σ1)且d,σ1,σ2都是非负的并满足d-σ1-σ2>0。因此,谱半径满足
这就保证了迭代算法的收敛性。
有益效果:本发明与现有技术相比,其显著优点是:本发明重点考虑了计算复杂度和算法性能;同时,对于较差的传播环境,本发明可以通过调整迭代矩阵灵活调节算法的迭代特性,保证收敛性和收敛速率,从而为性能需求不同的场合提供了更好的灵活性。且本发明为迭代算法,在软件编程方面可以节省内存消耗,在硬件实现方面具有节省面积的优点。
附图说明
图1是本发明提供的收敛速率可调的大规模MIMO迭代检测算法和传统基于Neumann级数展开的MIMO检测算法的计算复杂度对比图;
图2是发射天线(用户数)数为16,接收天线数为128时,采用本发明信号检测算法和其他传统检测算法的误码率对比图;
图3是发射天线数(用户数)为16,接收天线数为64时,采用本发明信号检测算法和其他传统检测算法的误码率对比图;
图4是发射天线数(用户数)为16,接收天线数为96时,考虑空间相关性的情况下(接收端相关因子为0.4),采用本发明的收敛速率可调的迭代检测算法和其他传统检测算法的误码率对比图。
具体实施方式
本实施例中建立一个大规模MIMO信道模型进行模拟操作。在大规模MIMO***中,一般有N>>M(基站天线数N远大于发射天线数,即用户数M)。首先M个不同用户产生的并行传输比特流分别通过信道编码进行编码,然后映射到星座符号,并采取星座图集合能量归一化。让s=[s1,s2,......,sM]T表示信号向量,s中包含了分别从M个用户产生的传输符号,采用64-QAM方式映射。H表示维度是N×M信道矩阵,故上行链路基站端的接收信号向量y可以表示为
y=Hs+n
其中y的维度为N×1,n为N×1维的加性白噪声向量,其元素服从零均值方差为No的高斯分布。上行链路多用户信号检测任务就是从接收机接收向量y=[y1,y2,......,yN]T估计传输信号符号s。假设H已知,其元素服从均值为0方差为1的独立同分布,采用最小均方误差(MMSE)线性检测理论,对传输信号向量的估计表示为
该估计过程等效为求解线性方程组
第l次迭代所得的检测结果为
sl=X-1(Ysl-1+yMF)
显然迭代矩阵X和Y对本发明中所提出的迭代算法收敛速率有很大影响,为了使得在较差的传播条件下(如天线比N/M较小或考虑天线相关性)保证算法依然收敛和提高收敛速度,我们采用收敛速率控制系数k调整迭代矩阵X和Y,从而得到
X=tril(W)-tril(W,k),Y=X-W
基于上述分析,大规模MIMO迭代检测方法包括以下步骤:
S1、将信道矩阵H和接收信号y通过匹配滤波器进行预处理,得到输出信号yMF和MMSE滤波矩阵W。
其中,yMF=HHy,W=G+N0ΙM,G=HHH为Gram矩阵,N0为噪声方差,ΙM为M维单位阵,(·)H为共轭转置操作。
S2、根据MMSE滤波矩阵W和预设收敛速率控制参数k计算得到迭代矩阵X和Y。
其中,X=tril(W)-tril(W,k),Y=X-W。
S3、设置迭代初始解为s0=0。
S4、按照公式sl=X-1(Ysl-1+yMF)进行迭代计算。
S5、达到预设迭代次数L后,迭代结束,所计算得到s的值为待检测信号的估计结果。
为了验证效果,分别在在迭代次数分别为2,3和4时,其计算复杂度(复数乘法数量)对比结果见图1。对于天线配置为128×16和64×16的大规模MIMO***,采用1/2速率卷积码以及64-QAM映射,本实施例的仿真结果见图2、图3;对于天线配置为96×16的大规模MIMO***,接收端天线相关因子为0.4,发射端相关因子分别为0.3和0.6时,本实施例的仿真结果见图4。通过以上复杂度对比结果和误比特率对比结果可以看出,本实施例在计算复杂度和检测性能方面均优于传统基于Neumann级数展开的大规模MIMO检测算法。
参考文献
[1].T.L.Marzetta,“Noncooperative Cellular Wireless with UnlimitedNumbers of Base Station Antennas,”IEEE Trans.Wireless Commun.,2010.
[2].F.Rusek,D.Persson,B.K.Lau,E.G.Larsson,T.L.Marzetta,O.Edfors,andF.Tufvesson,“Scaling Up MIMO:Opportunities and Challenges with Very LargeArrays,”IEEE Signal Process.Mag.,2013.
[3].E.G.Larsson,O.Edfors,F.Tufvesson,and T.L.Marzetta,“Massive MIMOfor Next GenerationWireless Systems,”IEEE Commun.Mag.,2014.
[4].S.Wu,L.Kuang,Z.Ni,J.Lu,D.Huang,and Q.Guo,“Low-ComplexityIterative Detection for Large-Scale Multiuser MIMO-OFDM Systems UsingApproximate Message Passing,”IEEE J.Sel.Topics in Sig.Proc.,2014.
[5].R.C.de Lamare,“Adaptive and Iterative Multi-Branch MMSE DecisionFeedback Detection Algorithms for Multi-Antenna Systems,”IEEE Trans.WirelessCommun.,2013.
[6].K.V.Vardhan,S.K.Mohammed,A.Chockalingam,and B.S.Rajan,“A Low-Complexity Detector for Large MIMO Systems and Multicarrier CDMA Systems,”IEEE J.Sel.Areas in Commun.,Apr.2008.
[7].P.Li and R.D.Murch,“Multiple Output Selection-LAS Algorithm inLarge MIMO Systems,”IEEE Commun.Lett.
[8].C.Studer,S.Fateh,and D.Seethaler,“ASIC Implementation of Soft-Input Soft-Output MIMO Detection Using MMSE Parallel InterferenceCancellation,”IEEE J.Solid-State Circuits,2011.
[9].M.Wu,B.Yin,G.Wang,C.Dick,J.R.Cavallaro,and C.Studer,“Large-ScaleMIMO Detection for 3GPP LTE:Algorithms and FPGA Implementations,”IEEEJ.Sel.Topics in Sig.Proc.,2014.
[10].F.Wang,C.Zhang,X.Liang,Z.Wu,S.Xu,and X.You,“Efficient IterativeSoft Detection Based on Polynomial Approximation for Massive MIMO,”inProc.IEEE WCSP,2015.
[11].L.Dai,X.Gao,X.Su,S.Han,C.L.I,and Z.Wang,“Low-Complexity Soft-Output Signal Detection Based on Gauss-Seidel Method for Uplink MultiuserLarge-Scale MIMO Systems,”IEEE Trans.Veh.Technol.,2015.
[12].I.B.Collings,M.R.G.Butler,and M.McKay,“Low Complexity ReceiverDesign for MIMO Bit-Interleaved Coded Modulation,”in Proc.IEEE ISSSTA,2004.
[13].Y.Saad,Iterative Methods for Sparse Linear Systems.plus 0.5emminus 0.4emSiam,2003.
[14].G.Stewart,Matrix Algorithms:Volume 1:Basic Decompositions.plus0.5em minus 0.4emSociety for Industrial and Applied Mathematics,1998.
[15].B.E.Godana and T.Ekman,“Parametrization Based Limited FeedbackDesign for Correlated MIMO Channels Using New Statistical Models,”IEEETrans.Wireless Commun.,2013.
Claims (3)
1.一种收敛速率可调的大规模MIMO迭代检测方法,其特征在于该方法包括:
(1)将信道矩阵H和接收信号y通过匹配滤波器进行预处理,得到输出信号yMF和MMSE滤波矩阵W;其中,yMF=HHy,W=G+N0ΙM,G=HHH为Gram矩阵,N0为噪声方差,ΙM为M维单位阵,(·)H为共轭转置操作;
(2)根据MMSE滤波矩阵W和预设收敛速率控制参数k计算得到迭代矩阵X和Y;其中,X=tril(W)-tril(W,k),Y=X-W,tril(W)表示取W的下三角阵,tril(W,k)表示取W的第k个对角线以下的元素;
(3)设置迭代初始解为s0=0;
(4)按照公式sl=X-1(Ysl-1+yMF)进行迭代计算;其中,sl表示第l次迭代时s的值,X-1表示X的逆;
(5)达到预设迭代次数后,迭代结束,所计算得到s的值为待检测信号的估计结果。
2.根据权利要求1所述的收敛速率可调的大规模MIMO迭代检测方法,其特征在于:步骤(2)中的控制参数k的取值范围为1~M次,M为发射天线数。
3.根据权利要求1所述的收敛速率可调的大规模MIMO迭代检测方法,其特征在于:步骤(4)步骤4中的迭代次数l为1~L,L为设置的最大迭代次数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610926372.6A CN106357318B (zh) | 2016-10-31 | 2016-10-31 | 收敛速率可调的大规模mimo迭代检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610926372.6A CN106357318B (zh) | 2016-10-31 | 2016-10-31 | 收敛速率可调的大规模mimo迭代检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106357318A true CN106357318A (zh) | 2017-01-25 |
CN106357318B CN106357318B (zh) | 2019-10-11 |
Family
ID=57863692
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610926372.6A Active CN106357318B (zh) | 2016-10-31 | 2016-10-31 | 收敛速率可调的大规模mimo迭代检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106357318B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107094043A (zh) * | 2017-05-23 | 2017-08-25 | 西安电子科技大学 | 基于块迭代法的改进后的mmse低复杂度信号检测方法 |
CN111478749A (zh) * | 2020-02-16 | 2020-07-31 | 西安电子科技大学 | 基于优化初值快收敛mimo迭代检测方法、***及应用 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1901388A (zh) * | 2005-07-20 | 2007-01-24 | 电子科技大学中山学院 | 一种低复杂度的线性迭代多用户检测装置与方法 |
CN1998160A (zh) * | 2004-04-22 | 2007-07-11 | 高通股份有限公司 | Mimo通信***的迭代本征向量计算 |
WO2008069467A1 (en) * | 2006-12-05 | 2008-06-12 | Electronics And Telecommunications Research Institute | Iterative receiver and method for detecting transmitted signal thereof |
CN101534177A (zh) * | 2009-04-03 | 2009-09-16 | 西安交通大学 | Mimo***中基于rumswf的低复杂度降秩均衡方法 |
CN101582743A (zh) * | 2009-06-15 | 2009-11-18 | 北京邮电大学 | 一种用于迭代接收机的mimo检测方法及*** |
US20110145671A1 (en) * | 2009-12-14 | 2011-06-16 | Young Jin Moon | Method and apparatus for iterative determination of mimo iterative receiver |
CN105680984A (zh) * | 2015-05-11 | 2016-06-15 | 南京理工大学 | 一种宽带多用户mimo-ofdm***的分组交替迭代检测方法 |
CN105978609A (zh) * | 2016-04-25 | 2016-09-28 | 东南大学 | 一种相关信道下大规模mimo线性检测硬件构架及方法 |
CN106059645A (zh) * | 2016-05-18 | 2016-10-26 | 郑州大学 | 大规模mimo***中低复杂度信道估计方法 |
-
2016
- 2016-10-31 CN CN201610926372.6A patent/CN106357318B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1998160A (zh) * | 2004-04-22 | 2007-07-11 | 高通股份有限公司 | Mimo通信***的迭代本征向量计算 |
CN1901388A (zh) * | 2005-07-20 | 2007-01-24 | 电子科技大学中山学院 | 一种低复杂度的线性迭代多用户检测装置与方法 |
WO2008069467A1 (en) * | 2006-12-05 | 2008-06-12 | Electronics And Telecommunications Research Institute | Iterative receiver and method for detecting transmitted signal thereof |
CN101534177A (zh) * | 2009-04-03 | 2009-09-16 | 西安交通大学 | Mimo***中基于rumswf的低复杂度降秩均衡方法 |
CN101582743A (zh) * | 2009-06-15 | 2009-11-18 | 北京邮电大学 | 一种用于迭代接收机的mimo检测方法及*** |
US20110145671A1 (en) * | 2009-12-14 | 2011-06-16 | Young Jin Moon | Method and apparatus for iterative determination of mimo iterative receiver |
CN105680984A (zh) * | 2015-05-11 | 2016-06-15 | 南京理工大学 | 一种宽带多用户mimo-ofdm***的分组交替迭代检测方法 |
CN105978609A (zh) * | 2016-04-25 | 2016-09-28 | 东南大学 | 一种相关信道下大规模mimo线性检测硬件构架及方法 |
CN106059645A (zh) * | 2016-05-18 | 2016-10-26 | 郑州大学 | 大规模mimo***中低复杂度信道估计方法 |
Non-Patent Citations (1)
Title |
---|
ANLAN YU, CHUAN ZHANG, SHUNQING ZHANG, AND XIAOHU YOU: "Efficient SOR-based detection and architecture for 《large-scale MIMO uplink》", 《IEEE ASIA PACIFIC CONFERENCE ON CIRCUITS AND SYSTEMS (APCCAS)》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107094043A (zh) * | 2017-05-23 | 2017-08-25 | 西安电子科技大学 | 基于块迭代法的改进后的mmse低复杂度信号检测方法 |
CN107094043B (zh) * | 2017-05-23 | 2021-07-09 | 西安电子科技大学 | 基于块迭代法的改进后的mmse低复杂度信号检测方法 |
CN111478749A (zh) * | 2020-02-16 | 2020-07-31 | 西安电子科技大学 | 基于优化初值快收敛mimo迭代检测方法、***及应用 |
CN111478749B (zh) * | 2020-02-16 | 2021-08-31 | 西安电子科技大学 | 基于优化初值快收敛mimo迭代检测方法、***及应用 |
Also Published As
Publication number | Publication date |
---|---|
CN106357318B (zh) | 2019-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101378277B (zh) | 多用户预编码及调度方法和实现该方法的基站 | |
CN101015160B (zh) | 具有减少的信道状态信息反馈的闭环多输入多输出*** | |
CN105071843B (zh) | 大规模mimo***低复杂度多项式展开矩阵求逆方法及应用 | |
CN101345592B (zh) | 应用于mimo的自适应信号的检测器及检测方法 | |
CN103685090A (zh) | 多输入多输出信道表现的预测方法 | |
Aref et al. | Deep learning-aided successive interference cancellation for MIMO-NOMA | |
CN101453444A (zh) | 星座图简化装置和方法、接收机 | |
CN106341169B (zh) | 一种多用户大规模mimo***上行链路的天线选择方法 | |
CN103188703A (zh) | 幸存星座点选择方法和qrm-mld信号检测方法 | |
Dong et al. | Improved joint antenna selection and user scheduling for massive MIMO systems | |
CN101026435A (zh) | 通信***中低复杂度的极大似然检测方法及装置 | |
CN106357309A (zh) | 基于非理想信道下大规模mimo线性迭代检测方法 | |
Chataut et al. | Huber fitting based ADMM detection for uplink 5G massive MIMO systems | |
US8107546B2 (en) | Detection method of space domain maximum posteriori probability in a wireless communication system | |
CN107919895A (zh) | 一种大规模多用户mimo***的分布式检测方法 | |
Isukapalli et al. | Packet error probability of a transmit beamforming system with imperfect feedback | |
CN101227254A (zh) | 一种在多入多出***中v-blast的检测方法 | |
Huttunen et al. | DeepTx: Deep learning beamforming with channel prediction | |
CN101964667B (zh) | 用于长期演进方案的高效多天线检测方法 | |
CN106357318B (zh) | 收敛速率可调的大规模mimo迭代检测方法 | |
CN107222246A (zh) | 一种近似mmse性能的高效大规模mimo检测方法及*** | |
Garth et al. | Exact symbol error probabilities for SVD transmission of BPSK data over fading channels | |
CN101616427B (zh) | 一种基站间的多用户mimo检测方法 | |
Mansoor et al. | Superimposed training based estimation of sparse MIMO channels for emerging wireless networks | |
CN110868244B (zh) | 一种基于信道穿刺的低复杂度通信信号检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |