CN110113084A - Mimo闭环传输***的信道预测方法 - Google Patents
Mimo闭环传输***的信道预测方法 Download PDFInfo
- Publication number
- CN110113084A CN110113084A CN201910491468.8A CN201910491468A CN110113084A CN 110113084 A CN110113084 A CN 110113084A CN 201910491468 A CN201910491468 A CN 201910491468A CN 110113084 A CN110113084 A CN 110113084A
- Authority
- CN
- China
- Prior art keywords
- matrix
- channel
- closed loop
- transmission system
- singular value
- 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
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/373—Predicting channel quality or other radio frequency [RF] parameters
-
- 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
-
- 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
- H04B7/0456—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Quality & Reliability (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radio Transmission System (AREA)
Abstract
本发明公开了一种MIMO闭环传输***的信道预测方法,属于数字通信领域。该方法根据前两个时刻的信道状态预测下一时刻的信道状态,然后反馈给发送端作预编码。包括:MIMO***通过一段时间的测量得到信道的统计信息,包括多普勒频偏、发送端和接收端的相关矩阵;接收机将当前时刻和前一时刻的信道矩阵进行奇异值分解,把右奇异矩阵的若干列提取出构建两个新矩阵,建模为Grassmannian流形的两个相邻点;基于Grassmannian流形的测地线理论,经过系列矩阵变换构建一条测地线,对下一时刻的信道状态做出预测,最终反馈给发送端。与传统预测方法相比,本发明方法与下一时刻真实的信道状态更为接近,弦距离误差性能更好。
Description
技术领域
本发明属于无线通信技术领域,特别涉及MIMO闭环传输***的信道预测方法。
背景技术
多天线通信***在收发端均配备多根天线,可以充分利用空间资源,具有较高的通信速率和可靠的通信质量,近年来受到了国内外广泛的关注。由于该***优点突出,目前已经作为关键技术被4G和5G通信***采用。
信道状态信息(channel state information,CSI)对于通信***十分重要,可以用来解码接收数据,或者反馈到发送端做预处理用以提高***性能。CSI通常采用信道估计获取,需要发送端每隔一定的时间***导频作辅助。但是,当信道变化较快时,为了准确的估计信道,不得不频繁的***导频,这无疑会造成导频的比例增加,占用了较多***资源。因此,需要对信道进行跟踪和预测。目前,多数信道跟踪算法是基于线性的信号空间,即欧几里得空间。另一方面,微分流形作为现代几何学的重要研究成果,推广了三维欧式空间中曲线和曲面概念,是拓扑学和几何学中一类重要的空间。以微分流形为基础的现代几何研究成果已经在模式识别、图像处理等工程领域中获得了一些应用。英国《应用概率论进展》(“Bayesian and geometric subspace tracking”,Advances in Applied Probability,2004,36(1):43-56)将时变阵列信号看成Grassmannian流形的点,基于每个点的切空间和贝叶斯理论提出了一种时变信号空间的跟踪方法。但是该方法十分复杂,并且需要数据的部分先验信息。美国《国际电气与电子工程师协会通信学报》(“Transmission subspacetracking for MIMO systems with low-rate feedback”,IEEE Transactions onCommunications,2007,55(8):1629-1639)提出了Grassmannian信道子空间跟踪法,其反馈链路速率很低,仅1比特。该算法整体而言性能良好,但是存在收敛速度较慢的问题。“Grassmannian subspace prediction for precoded spatial multiplexing MIMO withdelayed feedback”,IEEE Signal Processing Letters,2011,18(10):555-558)提出了一种Grassmannian子空间预测技术,但其弦误差性能不理想,预测精度有待提高。
发明内容
针对现有技术存在的上述问题,本发明的目的在于为MIMO闭环传输***的信道提供一种几何跟踪方法,不仅可降低导频信号的占比,减小***负担,而且提高预测性能和精度。
为了解决上述问题,本发明所采用的技术方案如下:
一种MIMO闭环传输***的信道预测方法,包括以下步骤:
(1)通过测量得到信道的统计信息,包括多普勒频偏、发送端和接收端的相关矩阵θT和θR;
(2)分别对当前时刻和前一时刻的信道状态矩阵H(t)和H(t-1)做奇异值分解,得到相应的右奇异值矩阵V(t)和V(t-1),分别提取前D列创建新的矩阵Vd(t)和Vd(t-1);
(3)构造一个特殊矩阵E,对E作奇异值分解,然后选择左奇异值矩阵中非零奇异值对应的列,建立一个新的矩阵
(4)对矩阵[Vd(t-1)]HVd(t)做奇异值分解,构建矩阵U2(t-1)和B(t-1);
(5)根据Vd(t)和Vd(t-1)建立基于Grassmannian流形的测地线方程,搜索找到最优预测步长,进而得到下一时刻的预编码矩阵反馈给发送端。
所述步骤2中的奇异值分解表示为:
H(t)=U(t)Λ(t)VH(t)
H(t-1)=U(t-1)Λ(t-1)VH(t-1)
其中上标H表示矩阵的Hermitian操作。
所述步骤4对矩阵[Vd(t-1)]HVd(t)做奇异值分解可以表述为:
[Vd(t-1)]HVd(t)=U1(t-1)C(t-1)[V1(t-1)]H
所述步骤4构建的矩阵U2(t-1)和B(t-1)如下:
其中,
其中,A(t-1)=U2(t-1)Φ(t-1)[U1(t-1)]H,Φ(t-1)=sin-1(S(t-1)),函数sin-1(.)表示对矩阵逐点求反正弦值。
所述步骤5建立的基于Grassmannian流形的测地线方程如下:
其中,NT×D的矩阵是从NT×NT的单位矩阵中选择前D列组成,函数EXP(.)代表矩阵的指数运算;s表示预测步长。
所述步骤5根据下式找到最优预测步长:
其中,函数f(.)是Grassmannian流形两个点的弦距离,表示如下:
||.||F代表矩阵的Frobenius范数,D是发送数据子流数,E(.)代表求期望。
相比于现有技术,本发明的有益效果为:
本发明将Grassmannia流形的测地线理论用到信道跟踪和预测中,根据前两个时刻的信道状态预测下一时刻的信道状态,然后反馈给发送端作预编码。与传统的预测方法相比,本发明方法与下一时刻真实的信道状态较为接近,弦距离误差性能更好。
附图说明
图1是本发明的MIMO闭环传输***的信道预测方法流程图;
图2是图1中方法在4发4收***中与传统方法的性能对比图。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明。应当了解,以下提供的实施例仅是为了详尽地且完全地公开本发明,并且向所属技术领域的技术人员充分传达本发明的技术构思,本发明还可以用许多不同的形式来实施,并且不局限于此处描述的实施例。对于表示在附图中的示例性实施方式中的术语并不是对本发明的限定。
本发明适用于点对点的MIMO闭环传输***,其中发送端和接收端均配备多根天线,***有一条反馈链路反馈下一时刻(时间段)的信道信息。MIMO***的输入输出关系可以表示为:
y(t)=H(t)W(t)x(t)+n(t) (1)
其中,H(t)是在时刻(时间片段)t的NR×NT维矩阵,W(t)是NT×D的预编码矩阵,该矩阵每列的Frobenius范数是1并且不同列是正交的;n(t)是高斯白噪声,x(t)是发送数据,其自相关矩阵满足E[x(t)xH(t)]=Id,Id是d阶单位矩阵,符号E(.)代表期望运算,上标H代表Hermitian转置。
MIMO信道建模为:
其中,NT×NT维矩阵θT和NR×NR维矩阵θR分别代表发送端和接收端的信道相关矩阵;Hω(t)是在时刻(时间片段)t的NR×NT维随机矩阵,矩阵每个元素都服从均值零、方差1的高斯分布,并且元素之间相互独立。令矩阵Hω(t)=(hij(t)),每个元素之间的时间相关性可以用Jake’s模型描述:
E[hij(t1)(hij(t2))*]=J0(2πfd(t2-t1)) (3)
其中fd是多普勒频偏,函数J0(.)是第一类零阶贝塞尔函数。
如图1所示,对MIMO闭环传输******的信道预测方法包括以下步骤:
步骤1,***通过一段时间的测量得到信道的统计信息,包括多普勒频偏fd、发送端和接收端的相关矩阵θT和θR。
在一个实施例中,MIMO***的收发端均配备4根天线,发送数据子流数目D=2。信道的相关矩阵θT和θR均采用指数相关模型,系数分别设置为0.2和0.3;归一化的多普勒频偏fdTs设为0.05,其中Ts是相邻时刻(时间片段)的间隔。
将统计信息代入公式(2)-(3)描述的信道模型,注意公式(3)的(t2-t1)应换成相邻时刻间隔Ts,得到信道矩阵样本。
步骤2,分别对当前时刻和前一时刻的4×4维信道矩阵H(t)和H(t-1)做奇异值分解,得到相应的右奇异值矩阵V(t)和V(t-1),分别提取前2列创建两个4×2矩阵Vd(t)和Vd(t-1),t=1,2…Ns,其中Ns是样本个数。
这里,H(t)和H(t-1)的奇异值分解分别表示为:
H(t)=U(t)Λ(t)VH(t) (4)
H(t-1)=U(t-1)Λ(t-1)VH(t-1) (5)
其中上标H表示矩阵的Hermitian操作。
通过步骤2提取出前一时刻和当前时刻的信道信息Vd(t-1)和Vd(t),建模为Grassmannian流形的两个点,以便构造这两个点的测地线方程。
步骤3,构造一个特殊矩阵E=I4-Vd(t-1)[Vd(t-1)]H,其中I4是4×4的单位矩阵;对E作奇异值分解,然后选择左奇异值矩阵中非零奇异值对应的列,建立一个新的4×2矩阵构造E是为了求出vd(t-1)的正交补矩阵容易验证E与vd(t-1)正交,但是E还不是正交矩阵,所以需对E作奇异值分解。
步骤4,对矩阵[Vd(t-1)]HVd(t)做奇异值分解,构建矩阵U2(t-1)和B(t-1)。这里,
[Vd(t-1)]HVd(t)=U1(t-1)C(t-1)[V1(t-1)]H (6)
其中,2×2矩阵4×4矩阵B(t-1)表示为:
其中,2×2矩阵A(t-1)=U2(t-1)Φ(t-1)[U1(t-1)]H,2×2矩阵Φ(t-1)=sin-1(S(t-1)),函数sin-1(.)表示对矩阵逐点求反正弦值。
步骤3和4所做的都是为了求出步骤5中的测地线方程。有了测地线方程,才能根据前两个点对未来的预编码做预测。
步骤5,根据Vd(t)和Vd(t-1)建立基于Grassmannian流形的测地线方程,搜索找到最优预测步长,得到下一时刻的预编码矩阵反馈给发送端。包括如下步骤:
5-1)根据Vd(t)和Vd(t-1)建立测地线方程:
其中,4×2的矩阵I4,2是从4×4的单位矩阵中选择前2列组成,函数EXP(.)代表矩阵的指数运算;4×4的矩阵
5-2)找到最优预测步长SOPT。
其中,函数f(.)代表Grassmannian流形两个点的弦距离,表示如下:
||.||F代表矩阵的Frobenius范数,E(.)代表求期望。式(11)表示给定预测步长s,预测的信道信息与真实的信道信息之间的弦距离误差。
当信道统计信息不变时,最优步长也不变,因此,此步骤只需要计算一次。
5-3)得到下一时刻的预编码矩阵。即代入最优预测步长sOPT到测地线方程,得到
5-4)反馈给发送端。
图2是在4发4收***下本发明方法与传统方法的弦距离误差性能对比,传统方法见《国际电气与电子工程师协会通信快报》(“Grassmannian subspace prediction forprecoded spatial multiplexing MIMO with delayed feedback”,IEEE SignalProcessing Letters,2011,18(10):555-558),这里弦距离误差由公式计算得到。可以看到,在预测步长s的整个变化范围内(0~1),提出方法的弦距离误差均比传统方法小。传统方法的弦距离误差在sOPT=0.7时达到最小值-8.80dB;提出方法的弦距离误差在sOPT=0.8时达到最小值-9.96dB,优于传统方法1.16dB。
Claims (7)
1.一种MIMO闭环传输***的信道预测方法,其特征在于,所述方法包括以下步骤:
(1)通过测量得到信道的统计信息,包括多普勒频偏、发送端和接收端的相关矩阵θT和θR,根据信道模型,得到信道状态矩阵;
(2)分别对当前时刻和前一时刻的信道状态矩阵H(t)和H(t-1)做奇异值分解,得到相应的右奇异值矩阵V(t)和V(t-1),分别提取前D列创建新的矩阵Vd(t)和Vd(t-1),D是发送数据子流数;
(3)构造一个特殊矩阵E,其中是NT×NT的单位矩阵,NT是发送端的天线数,对E作奇异值分解,选择左奇异值矩阵中非零奇异值对应的列,建立一个新的矩阵
(4)对矩阵[Vd(t-1)]HVd(t)做奇异值分解,构建矩阵U2(t-1)和B(t-1);
(5)根据Vd(t)和Vd(t-1)建立基于Grassmannian流形的测地线方程,搜索找到最优预测步长,进而得到下一时刻的预编码矩阵反馈给发送端。
2.根据权利要求1所述的MIMO闭环传输***的信道预测方法,其特征在于,所述步骤1中H(t)和H(t-1)的奇异值分解分别表示为:
H(t)=U(t)Λ(t)VH(t)
H(t-1)=U(t-1)Λ(t-1)VH(t-1)
其中上标H表示矩阵的Hermitian操作。
3.根据权利要求1所述的MIMO闭环传输***的信道预测方法,其特征在于,所述步骤4中对矩阵[Vd(t-1)]HVd(t)做奇异值分解表示为:
[Vd(t-1)]HVd(t)=U1(t-1)C(t-1)[V1(t-1)]H。
4.根据权利要求3所述的MIMO闭环传输***的信道预测方法,其特征在于,所述步骤4构建矩阵U2(t-1)和B(t-1)分别为:
其中,
A(t-1)=U2(t-1)Φ(t-1)[U1(t-1)]H,Φ(t-1)=sin-1(S(t-1)),函数sin-1(.)表示对矩阵逐点求反正弦值。
5.根据权利要求4所述的MIMO闭环传输***的信道预测方法,其特征在于,所述步骤5建立的基于Grassmannian流形的测地线方程为:
其中,是从NT×NT的单位矩阵中选择前D列组成的NT×D矩阵,函数EXP(.)代表矩阵的指数运算;s表示预测步长。
6.根据权利要求5所述的MIMO闭环传输***的信道预测方法,其特征在于,所述步骤5根据下式找到最优预测步长SOPT:
其中,E(.)代表求期望,函数f(.)是Grassmannian流形两个点的弦距离,表示形式为:
||.||F代表矩阵的Frobenius范数,D是发送数据子流数。
7.根据权利要求6所述的MIMO闭环传输***的信道预测方法,其特征在于,所述步骤5得到下一时刻的预编码矩阵方法如下:代入最优预测步长sOPT到测地线方程,得到
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910491468.8A CN110113084B (zh) | 2019-06-06 | 2019-06-06 | Mimo闭环传输***的信道预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910491468.8A CN110113084B (zh) | 2019-06-06 | 2019-06-06 | Mimo闭环传输***的信道预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110113084A true CN110113084A (zh) | 2019-08-09 |
CN110113084B CN110113084B (zh) | 2022-02-01 |
Family
ID=67494154
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910491468.8A Active CN110113084B (zh) | 2019-06-06 | 2019-06-06 | Mimo闭环传输***的信道预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110113084B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112311706A (zh) * | 2020-11-03 | 2021-02-02 | 南京大学 | 基于贝叶斯神经网络信道预测的频分双工下行传输方法 |
US11121788B1 (en) | 2020-06-08 | 2021-09-14 | Wuhan University | Channel prediction method and system for MIMO wireless communication system |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110013603A1 (en) * | 2009-07-20 | 2011-01-20 | Qinghua Li | Techniques for MIMO beamforming for frequency selective channels in wireless communication systems |
US20120082198A1 (en) * | 2010-09-30 | 2012-04-05 | Nec (China) Co., Ltd. | Method and apparatus for obtaining channel state information required for beamforming |
CN103299570A (zh) * | 2011-02-28 | 2013-09-11 | 日电(中国)有限公司 | 用于预测mimo***中的预编码矩阵的方法和装置 |
US20130272438A1 (en) * | 2011-03-16 | 2013-10-17 | Nec (China) Co., Ltd. | Method, transmitter and receiver for beamforming |
US8711961B2 (en) * | 2010-07-15 | 2014-04-29 | The Board Of Regents Of The University Of Texas System | Communicating channel state information using predictive vector quantization |
US9071301B2 (en) * | 2006-02-14 | 2015-06-30 | Nec Laboratories America, Inc. | Precoding with a codebook for a wireless system |
CN105515622A (zh) * | 2015-11-25 | 2016-04-20 | 南京邮电大学 | 基于高分辨率动态聚焦码本的有限反馈预测预编码方法 |
-
2019
- 2019-06-06 CN CN201910491468.8A patent/CN110113084B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9071301B2 (en) * | 2006-02-14 | 2015-06-30 | Nec Laboratories America, Inc. | Precoding with a codebook for a wireless system |
US20110013603A1 (en) * | 2009-07-20 | 2011-01-20 | Qinghua Li | Techniques for MIMO beamforming for frequency selective channels in wireless communication systems |
US8711961B2 (en) * | 2010-07-15 | 2014-04-29 | The Board Of Regents Of The University Of Texas System | Communicating channel state information using predictive vector quantization |
US20120082198A1 (en) * | 2010-09-30 | 2012-04-05 | Nec (China) Co., Ltd. | Method and apparatus for obtaining channel state information required for beamforming |
CN103299570A (zh) * | 2011-02-28 | 2013-09-11 | 日电(中国)有限公司 | 用于预测mimo***中的预编码矩阵的方法和装置 |
US20130272438A1 (en) * | 2011-03-16 | 2013-10-17 | Nec (China) Co., Ltd. | Method, transmitter and receiver for beamforming |
CN105515622A (zh) * | 2015-11-25 | 2016-04-20 | 南京邮电大学 | 基于高分辨率动态聚焦码本的有限反馈预测预编码方法 |
Non-Patent Citations (5)
Title |
---|
DALIN ZHU AND MING LEI: "On the Quantization and Prediction for Precoded MIMO with Delayed Limited Feedback", 《2012 IEEE 75TH VEHICULAR TECHNOLOGY CONFERENCE》 * |
DALIN ZHU ET AL: "Grassmannian Subspace Prediction for Precoded Spatial Multiplexing MIMO With Delayed Feedback", 《IEEE SIGNAL PROCESSING LETTERS》 * |
YU ZHANG ET AL: "Predictive Unitary Precoding for Spatial Multiplexing Systems in Temporally Correlated MIMO Channels with Delayed Limited Feedback", 《2011 IEEE GLOBAL TELECOMMUNICATIONS CONFERENCE》 * |
李汀: "Grassmannian流形上基于高分辨率动态聚焦码本的的MIMO信道预测", 《信号处理》 * |
杨子琪,郁进明: "基于格拉斯曼流形的无线信道预测", 《信息通信》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11121788B1 (en) | 2020-06-08 | 2021-09-14 | Wuhan University | Channel prediction method and system for MIMO wireless communication system |
CN112311706A (zh) * | 2020-11-03 | 2021-02-02 | 南京大学 | 基于贝叶斯神经网络信道预测的频分双工下行传输方法 |
CN112311706B (zh) * | 2020-11-03 | 2021-11-16 | 南京大学 | 基于贝叶斯神经网络信道预测的频分双工下行传输方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110113084B (zh) | 2022-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113395221B (zh) | 一种基于正交时频空联合信道估计与符号检测方法 | |
Kang et al. | Deep-learning-based channel estimation for wireless energy transfer | |
CN107332598B (zh) | 一种基于深度学习的mimo***联合预编码和天线选择方法 | |
CN103347268B (zh) | 认知传感器网络中基于能量有效性观测的自适应压缩重构方法 | |
CN110912598A (zh) | 基于长短时注意力机制的大规模mimo***csi反馈方法 | |
US20100054355A1 (en) | Wireless communication system | |
CN108512787B (zh) | 大规模mimo***的超精细信道估计方法 | |
CN113972939B (zh) | 基于双时间尺度和深度学习的天线***预编码方法和装置 | |
TWI591973B (zh) | A signal detection method and device | |
CN104702390A (zh) | 分布式压缩感知信道估计中的导频分配方法 | |
CN107483091A (zh) | 一种fdd大规模mimo‑ofdm***下的信道信息反馈算法 | |
CN110808764B (zh) | 一种大规模mimo中继***联合信息估计方法 | |
Marseet et al. | Application of complex-valued convolutional neural network for next generation wireless networks | |
CN110113084A (zh) | Mimo闭环传输***的信道预测方法 | |
Fan et al. | Fully convolutional neural network-based CSI limited feedback for FDD massive MIMO systems | |
CN115941001A (zh) | 一种基于mimo***的信息传输收发装置、***和方法 | |
CN106788641B (zh) | 一种信息与能量联合传输***的预编码传输方法 | |
CN102710392B (zh) | 基于功率约束的连续梯度搜索垂直分层空时码检测方法 | |
CN108667498A (zh) | 反馈受限下多天线传输的有效容量优化方法 | |
CN105207959A (zh) | 一种miso-ofdm***的信道估计方法 | |
Chu et al. | Wireless channel prediction for multi-user physical layer with deep reinforcement learning | |
CN114006639B (zh) | 一种适用于上行多用户mimo***的分布式预编码方法 | |
CN103475603A (zh) | 基于序参量非正交变换的通信***盲信道估计方法 | |
CN115378768B (zh) | 一种基于空间调制***的卫星通信信道估计方法及*** | |
CN114389730B (zh) | 一种基于深度学习和脏纸编码的miso***波束形成设计方法 |
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 |