CN110855339A - 一种低复杂度的空时键控离散矩阵生成方法 - Google Patents

一种低复杂度的空时键控离散矩阵生成方法 Download PDF

Info

Publication number
CN110855339A
CN110855339A CN201911051440.9A CN201911051440A CN110855339A CN 110855339 A CN110855339 A CN 110855339A CN 201911051440 A CN201911051440 A CN 201911051440A CN 110855339 A CN110855339 A CN 110855339A
Authority
CN
China
Prior art keywords
matrix
discrete
data
equal
making
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
CN201911051440.9A
Other languages
English (en)
Other versions
CN110855339B (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.)
Donghua University
National Dong Hwa University
Original Assignee
Donghua University
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 Donghua University filed Critical Donghua University
Priority to CN201911051440.9A priority Critical patent/CN110855339B/zh
Publication of CN110855339A publication Critical patent/CN110855339A/zh
Application granted granted Critical
Publication of CN110855339B publication Critical patent/CN110855339B/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/0456Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Radio Transmission System (AREA)

Abstract

本发明提供了一种低复杂度的空时键控离散矩阵生成方法。首先对离散矩阵每个元素的取值区间进行离散化,从而把连续优化问题转化为离散优化问题,搜索区间相应的由无穷多数据集转化为有限数据集。然后提出了一种交替最优方法来搜寻最优离散矩阵集,该方法每次仅优化离散矩阵集中的一个数据,在此基础上再优化下一个数据,直到所有的数据都优化完毕。相比于传统空时键控中通过大量随机搜索寻找最优矩阵的方式,本发明能够大大降低矩阵搜寻过程中的复杂度,并能够保持优异的性能。

Description

一种低复杂度的空时键控离散矩阵生成方法
技术领域
本发明涉及一种低复杂度的空时键控离散矩阵生成方法,属于通信技术领域。
背景技术
随着无线通信行业的迅猛发展,无线通信技术对数据传输量的需求也骤然增 加,人们对一种高效可靠的带宽传输技术提出了更高的要求。多输入多输出 (Multiple InputMultiple Output,MIMO)传输技术就是一种新兴的无线通信技术, 它在发射端和接收端装配多根天线,能够显著提高频谱效率和增强系统可靠性, 经过近些年的发展,MIMO技术已经成为学术界的研究热点。
空间调制(Spatial Modulation,SM)技术和空移键控技术(Spatial-ShiftKeying,SSK)是在MIMO技术基础上发展而来的,其核心思想是在每一个符号 间隔内仅激活发射天线中的一个,将发射天线的索引作为一种额外携带信息的方 式,通过建立不同输入比特与天线序号的映射关系,实现充分利用空间资源的多 维调制。
为了进一步挖掘空间调制传输中的空间和时间资源,在SM和SSK的基础 上,近些年出现了一种新的技术—空时键控(Spatial Time Shift Keying,STSK), 它不再以天线的索引作为信息携带的载体,而是采用了预先设计的离散矩阵的激 活来携带信息,类似于SM的思想,输入比特被分为两个部分,一部分用来选取 离散矩阵的索引,其余部分用于映射AMP星座点。由于STSK可以按照不同的 准则来设计离散矩阵,如容量最大化和行列式准则,与SM相比,它能够灵活的 权衡MIMO系统中的分集和复用增益。
对于STSK来说,空时离散矩阵集对其性能的影响非常大,因此其优化问题 备受关注。通常来说,对STSK离散矩阵的研究通常采用随机搜索方法,即在满 足能量约束的条件下,空时离散矩阵集是使用满足高斯分布的数随机生成的。然 后根据特定的准则(如行列式准则)计算相应的行列式,当随机搜索过程重复一 定的次数后,选择具有对应行列式最小值最大的离散矩阵集作为设计矩阵。
通过多次搜索,可以达到较好的性能,但是随着搜索空间的增加,随机搜索 方法的复杂度会非常的高。
发明内容
本发明的目的是:提供一种比传统随机搜索方法复杂度更低的搜索方法。
为了达到上述目的,本发明的技术方案是提供了
本发明提供的一种低复杂度的空时键控离散矩阵生成方法在保持传统搜索 算法性能的基础上显著降低了搜索的复杂度,具有如下优点:
(1)对离散矩阵每个元素的取值区间进行离散化,从而把连续优化问题转 化为离散优化问题,搜索区间相应的由无穷多数据集转化为有限数据集。
(2)提出了一种交替最优方法来搜寻最优离散矩阵集,该方法摆脱了传统 的高复杂度的大量搜索方法,能够大大降低系统的复杂度。具体的:
A、一般情况下,利用传统随机搜索方法一共需要进行(105~107)N'次行列 式计算,其中N'代表STSK系统使用一次行列式准则需要进行的行列式计算次 数;
B、本发明提供方法的方法一共需要四层循环,从外到内分别需要进行 Q,M,T,22k次循环,故总的循环次数为22kQMT。则采用本发明提供的方法一共需 要进行22kQMTN'次行列式计算。其中Q是STSK系统需要的离散矩阵个数,M 是STSK系统中发射天线的数目,T代表每个时空块的符号数,22k是生成的离 散复数个数。
C、在一般情况下22kQMTN'<<(105~107)N',故本发明能够大大降低系统的 复杂度,相比于传统的随机搜索法有明显的优势。
附图说明
图1为STSK系统框图;
图2为本发明提出的方法与传统搜索方法比较的系统性能图。
具体实施方式
下面结合具体实施例,进一步阐述本发明。应理解,这些实施例仅用于说明 本发明而不用于限制本发明的范围。此外应理解,在阅读了本发明讲授的内容之 后,本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样落于本 申请所附权利要求书所限定的范围。
本发明提供的一种低复杂度的空时键控离散矩阵生成方法首先生成一组离 散数据,然后提出了一种交替最优方法来搜寻最优离散矩阵集,该方法通过遍历 生成的离散数据,每次仅优化离散矩阵集中的一个数据,在此基础上再优化下一 个数据,直到所有的数据都优化完毕,具体包括以下步骤:
第一步:产生离散数据集。
在区间-Am~+Am上等间隔产生2k个实数离散数据,数据间隔为则生成的数据集为:
根据集合R,生成22k个复数数据集P,生成方法是实部和虚部分别遍历集 合R。
第二步、用交替最优方法找出一组最优离散矩阵
步骤2-1、随机生成一组矩阵Bq∈CM×T(q=1,…,Q),其每一个元素都服从均 值为0,方差为1的正态分布。其中Q是STSK系统需要的离散矩阵个数,M是 STSK系统中发射天线的数目,T代表每个时空块的符号数,C代表复数集。
步骤2-2、对矩阵组Bq∈CM×T(q=1,…,Q)进行归一化操作,以便矩阵组能够 满足能量约束tr[Bq HBq]=T(q=1,…,Q),其中tr[·]表示求迹的操作。
步骤2-3、令最优的离散矩阵组DMq∈CM×T(q=1,…,Q)的初始值和矩阵组 Bq∈CM×T(q=1,…,Q)相等。
步骤2-4、用集合P中的第p个数据来替换Bq(m,t),其中Bq(m,t)代表矩阵 组Bq的第q个矩阵的第m行第t列中的数据。在首次执行步骤2-4时, p、t、m、q的初始值都为1。
步骤2-5、对替换后的矩阵Bq进行归一化操作。
步骤2-6、在QL(QL-1)/2种(l,q)和(l',q')(l≠l',q≠q')的组合中计算
Figure BDA0002255418020000033
的行列式的值,并在所有可能的结果中保存一个最小值Dp。其中,(l,q)和(l',q') 分别对应系统的发射符号Sl,q=slBq和Sl',q'=sl'Bq',sl是选择的星座点;L是调 制星座点的总个数;
Figure BDA0002255418020000041
Figure BDA0002255418020000042
的秩的值满足
Figure BDA0002255418020000043
Figure BDA0002255418020000044
代表
Figure BDA0002255418020000045
的共轭转置,H是共轭转置符号。
步骤2-7、当p<22k时,令p=p+1,重复执行步骤2-2~2-7;否则令p=1, 结束循环,进入下一步骤。
步骤2-8、在Dp(p=1,…,22k)找出一个最大值,用这个最大值对应的Bq的归 一化形式来替换DMq
步骤2-9、更新矩阵Bq的值,令Bq=DMq
步骤2-10、当t<T时,令t=t+1,重复执行步骤2-2~2-10;否则令t=1, 结束循环,进入下一步骤。
步骤2-11、当m<M时,令m=m+1,重复执行步骤2-2~2-11;否则令m=1, 结束循环,进入下一步骤。
步骤2-12、当q<Q时,令q=q+1,重复执行步骤2-2~2-12;否则结束 循环,进入下一步。
步骤2-13、返回矩阵组DMq(q=1,…,Q)作为最优矩阵组。
在本发明中选取了行列式准则作为离散矩阵的优化准则(本发明与优化准则 无关),该准则行列式计算部分如步骤2-6所示。由于(l,q)和(l',q')(l≠l',q≠q')共 有
Figure BDA0002255418020000047
种组合,故每次进行步骤2-6操作需要进行 N'=QL(QL-1)/2次的行列式计算。
为方便对比复杂度,下面将在此准则下使用随机搜索方法来搜寻离散矩阵的 步骤列出:
步骤1、根据系统要求,随机产生一组元素服从高斯分布的离散矩阵 Aq∈CM×T(q=1,...,Q),并进行归一化操作以满足能量约束:
tr[Aq HAq]=T(q=1,...,Q)
步骤2、在QL(QL-1)/2种(l,q)和(l',q')(l≠l',q≠q')的组合中计算
Figure BDA0002255418020000046
的行列式的值,并在所有可能的结果中保存一个最小值。其中(l,q)和(l',q')对应系 统发射的空时信号Sl,q=slAq和Sl',q'=sl'Aq';sl是选择的星座点,L代表调制星 座点的总个数;并且
Figure BDA0002255418020000052
的秩的值满足
Figure BDA0002255418020000053
Figure BDA0002255418020000054
代表的共轭转置,H是共轭转置符号。
步骤3、重复步骤1和步骤2,找出所有行列式最小值中最大的一组离散矩 阵作为最优的离散矩阵,一般需要重复的次数为105~107
图1为STSK系统框图,本例STSK系统有M=3根发射天线和N=2根接 收天线,需要Q=4个离散矩阵,每个时空块的符号数T=2,即 STSK(M,N,T,Q)=STSK(3,2,2,4);两种方法调制方式均为QPSK,调制星座 点的总个数是L=4;传统方法搜索的次数106次,选择的数据的区间为-1~+1, 产生22个实数数据。
STSK系统调制信号和离散矩阵的映射表如下表所示:
输入比特 离散矩阵 星座点 输入比特 离散矩阵 星座点
0000 A<sub>1</sub> s<sub>1</sub> 1000 A<sub>3</sub> s<sub>1</sub>
0001 A<sub>1</sub> s<sub>2</sub> 1001 A<sub>3</sub> s<sub>2</sub>
0010 A<sub>1</sub> s<sub>3</sub> 1010 A<sub>3</sub> s<sub>3</sub>
0011 A<sub>1</sub> s<sub>4</sub> 1011 A<sub>3</sub> s<sub>4</sub>
0100 A<sub>2</sub> s<sub>1</sub> 1100 A<sub>4</sub> s<sub>1</sub>
0101 A<sub>2</sub> s<sub>2</sub> 1101 A<sub>4</sub> s<sub>2</sub>
0110 A<sub>2</sub> s<sub>3</sub> 1110 A<sub>4</sub> s<sub>3</sub>
0111 A<sub>2</sub> s<sub>4</sub> 1111 A<sub>4</sub> s<sub>4</sub>
下面结合具体实例,进一步阐述本发明。
第一步:产生离散数据集。
在区间-1~+1上等间隔产生22个实数离散数据,数据间隔为Δ=0.5,则生 成的数据集为:R={-0.75,-0.25,0.25,0.75}。
根据集合R,生成24个复数数据集P,生成方法是,实部和虚部分别遍历 集合R。
第二步、用交替最优方法找出一组最优离散矩阵
步骤2-1、随机生成一组矩阵Bq∈C3×2(q=1,…,4),其每一个元素都服从均 值为0,方差为1的正态分布,其中C代表复数集。
步骤2-2、对矩阵组Bq∈C3×2(q=1,…,4)进行归一化操作,以便矩阵组能够 满足能量约束tr[Bq HBq]=2(q=1,…,4),其中tr[·]表示求迹的操作。
步骤2-3、令最优的离散矩阵组DMq∈C3×2(q=1,…,4)的初始值和矩阵组 Bq∈C3×2(q=1,…,4)相等。
步骤2-4、用集合P中的第p个数据来替换Bq(m,t),其中Bq(m,t)代表矩阵 组Bq的第q个矩阵的第m行第t列中的数据。
步骤2-5、对替换后的矩阵Bq进行归一化操作。
步骤2-6、在120种(l,q)和(l',q')(l≠l',q≠q')的组合中计算
Figure BDA0002255418020000061
的行列式的值,并在所有可能的结果中保存一个最小值Dp。其中,(l,q)和(l',q')分别对 应系统的发射符号Sl,q=slBq和Sl',q'=sl'Bq',sl是选择的星座点;L是调制星座 点的总个数;
Figure BDA0002255418020000062
Figure BDA0002255418020000063
的秩的值满足
Figure BDA0002255418020000064
Figure BDA0002255418020000065
代表的共轭转置,H是共轭转置符号。
步骤2-7、当p<22k时,令p=p+1,重复执行步骤2-2~2-7;否则令p=1, 结束循环,进入下一步骤。
步骤2-8、在Dp(p=1,…,24)找出一个最大值,用这个最大值对应的Bq的归 一化形式来替换DMq
步骤2-9、更新矩阵Bq的值,令Bq=DMq
步骤2-10、当t<2时,令t=t+1,重复执行步骤2-2~2-10;否则令t=1, 结束循环,进入下一步骤。
步骤2-11、当m<3时,令m=m+1,重复执行步骤2-2~2-11;否则令m=1, 结束循环,进入下一步骤。
步骤2-12、当q<4时,令q=q+1,重复执行步骤2-2~2-12;否则结束循 环,进入下一步。
步骤2-13、返回矩阵组DMq(q=1,…,4)作为最优矩阵组。
在STSK(3,2,2,4)情况下,N'=120,则传统的搜索方式需要进行 106N'=1.2×108次行列式的计算,而本发明提出的方法仅需要384N'=46080次 行列式计算,大大降低了复杂度。由图2可知,在相同的情况下,本发明提出的 方法在降低复杂度的同时也能够获得比传统方式更加优异的误码率性能。

Claims (3)

1.一种低复杂度的空时键控离散矩阵生成方法,其特征在于,包括以下步骤:
第一步:产生离散数据集
在区间-Am~+Am上等间隔产生2k个实数离散数据,数据间隔为
Figure FDA0002255418010000011
则生成的数据集为:
Figure FDA0002255418010000012
根据数据集R,生成22k个复数数据集P;
第二步、用交替最优方法找出一组最优离散矩阵,包括以下步骤
步骤2-1、随机生成一组矩阵Bq∈CM×T(q=1,…,Q),其每一个元素都服从均值为0,方差为1的正态分布,其中Q是STSK系统需要的离散矩阵个数,M是STSK系统中发射天线的数目,T代表每个时空块的符号数,C代表复数集;
步骤2-2、对矩阵组Bq∈CM×T(q=1,…,Q)进行归一化操作,以便矩阵组能够满足能量约束tr[Bq HBq]=T(q=1,…,Q),其中tr[·]表示求迹的操作;
步骤2-3、令最优的离散矩阵组DMq∈CM×T(q=1,…,Q)的初始值和矩阵组Bq∈CM×T(q=1,…,Q)相等;
步骤2-4、用复数数据集P中的第p个数据来替换Bq(m,t),其中Bq(m,t)代表矩阵组Bq的第q个矩阵的第m行第t列中的数据;
步骤2-5、对替换后的矩阵Bq进行归一化操作;
步骤2-6、在QL(QL-1)/2种(l,q)和(l',q')(l≠l',q≠q')的组合中计算
Figure FDA0002255418010000013
的行列式的值,并在所有可能的结果中保存一个最小值Dp。其中,(l,q)和(l',q')分别对应系统的发射符号Sl,q=slBq和Sl',q'=sl'Bq',sl是选择的星座点;L是调制星座点的总个数;
Figure FDA0002255418010000014
Figure FDA0002255418010000015
的秩的值满足
Figure FDA0002255418010000016
Figure FDA0002255418010000017
代表
Figure FDA0002255418010000018
的共轭转置,H是共轭转置符号。
步骤2-7、当p<22k时,令p=p+1,重复执行步骤2-2~2-7;否则令p=1,结束循环,进入下一步骤;
步骤2-8、在Dp(p=1,…,22k)找出一个最大值,用这个最大值对应的Bq的归一化形式来替换DMq
步骤2-9、更新矩阵Bq的值,令Bq=DMq
步骤2-10、当t<T时,令t=t+1,重复执行步骤2-2~2-10;否则令t=1,结束循环,进入下一步骤;
步骤2-11、当m<M时,令m=m+1,重复执行步骤2-2~2-11;否则令m=1,结束循环,进入下一步骤;
步骤2-12、当q<Q时,令q=q+1,重复执行步骤2-2~2-12;否则结束循环,进入下一步;
步骤2-13、返回矩阵组DMq(q=1,…,Q)作为最优矩阵组。
2.如权利要求1所述的一种低复杂度的空时键控离散矩阵生成方法,其特征在于,生成22k个所述复数数据集P的方法是实部和虚部分别遍历所述数据集R。
3.如权利要求1所述的一种低复杂度的空时键控离散矩阵生成方法,其特征在于,在首次执行步骤2-4时,p、t、m、q的初始值都为1。
CN201911051440.9A 2019-10-31 2019-10-31 一种低复杂度的空时键控离散矩阵生成方法 Active CN110855339B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911051440.9A CN110855339B (zh) 2019-10-31 2019-10-31 一种低复杂度的空时键控离散矩阵生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911051440.9A CN110855339B (zh) 2019-10-31 2019-10-31 一种低复杂度的空时键控离散矩阵生成方法

Publications (2)

Publication Number Publication Date
CN110855339A true CN110855339A (zh) 2020-02-28
CN110855339B CN110855339B (zh) 2021-11-12

Family

ID=69599442

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911051440.9A Active CN110855339B (zh) 2019-10-31 2019-10-31 一种低复杂度的空时键控离散矩阵生成方法

Country Status (1)

Country Link
CN (1) CN110855339B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105099625A (zh) * 2015-08-14 2015-11-25 广东顺德中山大学卡内基梅隆大学国际联合研究院 一种用于多维调制的空时键控最优坐标组合搜索方法
CN105119869A (zh) * 2015-06-30 2015-12-02 深圳清华大学研究院 基于空时矩阵星座图的空移键控通信方法
CN105207751A (zh) * 2015-08-14 2015-12-30 广东顺德中山大学卡内基梅隆大学国际联合研究院 基于多维调制的空时键控方法及其联合优化方法
US9813278B1 (en) * 2013-10-31 2017-11-07 Sensor Networks And Cellular System Center, University Of Tabuk Quadrature spatial modulation system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9813278B1 (en) * 2013-10-31 2017-11-07 Sensor Networks And Cellular System Center, University Of Tabuk Quadrature spatial modulation system
CN105119869A (zh) * 2015-06-30 2015-12-02 深圳清华大学研究院 基于空时矩阵星座图的空移键控通信方法
CN105099625A (zh) * 2015-08-14 2015-11-25 广东顺德中山大学卡内基梅隆大学国际联合研究院 一种用于多维调制的空时键控最优坐标组合搜索方法
CN105207751A (zh) * 2015-08-14 2015-12-30 广东顺德中山大学卡内基梅隆大学国际联合研究院 基于多维调制的空时键控方法及其联合优化方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
GUANGTAO ZHENG等: "Multi-Dimensional Space–Time Shift Keying for Wireless Communications", 《IEEE ACCESS》 *
S. SUGIURA等: "Space-Time Shift Keying: A Unified MIMO Architecture", 《2010 IEEE GLOBAL TELECOMMUNICATIONS CONFERENCE GLOBECOM 2010》 *
SHINYA SUGIURA: "Dispersion Matrix Optimization for Space-Time Shift Keying", 《IEEE COMMUNICATIONS LETTERS》 *
SHINYA SUGIURA等: "Coherent and Differential Space-Time Shift Keying: A Dispersion Matrix Approach", 《IEEE TRANSACTIONS ON COMMUNICATIONS》 *

Also Published As

Publication number Publication date
CN110855339B (zh) 2021-11-12

Similar Documents

Publication Publication Date Title
Swindlehurst et al. Minimum probability-of-error perturbation precoding for the one-bit massive MIMO downlink
CN101427507B (zh) 根据信道相位反馈使用空时格码的方法和设备
CN107359921B (zh) 大规模mimo系统基于标准正交化的混合预编码方法
CN108366036B (zh) 一种面向稀疏码多址接入系统的调制码本设计方法
CN110661555B (zh) 用于大规模mimo的部分连接的移相器网络的混合预编码算法
CN112350814A (zh) 一种高效的上行链路scma码本设计
CN110855328B (zh) 一种基于天线分组的差分空间调制方法、设备及存储介质
CN110336614B (zh) 一种适合于无线光通信的多层空间脉冲调制方法
CN103607262B (zh) 空时分组编码mimo系统中的两阶段预编码方法
CN107425894B (zh) 基于信道范数的广义空间调制系统收发端天线选择方法
CN113708811B (zh) 一种毫米波大规模mimo系统中的混合预编码设计方法
CN110677182A (zh) 基于上行链路分层空时结构scma码本的通信方法
Wen et al. Designing enhanced multidimensional constellations for code-domain NOMA
CN110855339B (zh) 一种低复杂度的空时键控离散矩阵生成方法
CN112003680B (zh) 一种scma系统中低复杂度多用户检测方法
CN106899530B (zh) 一种基于门限值的天线组合和调制阶数联合优化方法
CN116260502B (zh) 一种基于可重构智能表面的双域索引调制通信方法
CN101129037B (zh) 经由特征坐标变换的酉矩阵码本的改进Hochwald构造
CN110190880B (zh) 一种广义空间调制递增发射天线组合选择方法
CN106301495A (zh) 基于最大化信道容量准则的广义空间调制系统天线选择方法
CN107493123B (zh) 基于预编码辅助广义正交空间调制的低复杂度检测方法
CN113067672A (zh) 一种非正交多址接入移动通信低复杂度接收方法
CN109257077B (zh) 一种结合霍夫曼编码和bch编码的新型空间调制方法
CN112543163A (zh) 基于星座衍生的scma母码本设计方法
CN106788630A (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