CN108712196A - 低分辨率毫米波大规模mimo混合预编码系统和方法 - Google Patents
低分辨率毫米波大规模mimo混合预编码系统和方法 Download PDFInfo
- Publication number
- CN108712196A CN108712196A CN201810151431.6A CN201810151431A CN108712196A CN 108712196 A CN108712196 A CN 108712196A CN 201810151431 A CN201810151431 A CN 201810151431A CN 108712196 A CN108712196 A CN 108712196A
- Authority
- CN
- China
- Prior art keywords
- low resolution
- signal
- quantization
- precoding
- digital
- 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
- 238000000034 method Methods 0.000 title claims description 24
- 238000013139 quantization Methods 0.000 claims abstract description 99
- 238000004088 simulation Methods 0.000 claims abstract description 50
- 238000006243 chemical reaction Methods 0.000 claims abstract description 5
- 239000013598 vector Substances 0.000 claims description 46
- 238000004422 calculation algorithm Methods 0.000 claims description 28
- 238000013507 mapping Methods 0.000 claims description 28
- 239000011159 matrix material Substances 0.000 claims description 23
- 230000005540 biological transmission Effects 0.000 claims description 11
- 238000000354 decomposition reaction Methods 0.000 claims description 6
- 230000010363 phase shift Effects 0.000 claims description 4
- 230000017105 transposition Effects 0.000 claims description 2
- 238000004891 communication Methods 0.000 description 16
- 238000005516 engineering process Methods 0.000 description 16
- 238000013461 design Methods 0.000 description 13
- 238000001228 spectrum Methods 0.000 description 11
- 238000012545 processing Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 238000005457 optimization Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 238000010295 mobile communication Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005562 fading Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005538 encapsulation Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
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/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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Radio Transmission System (AREA)
Abstract
本发明公开了一种低分辨率毫米波大规模MIMO混合预编码系统,包括:发射机和接收机,其中所述发射机通过数字预编码器对发射信号进行数字预编码,经发送量化模块转换为低分辨率数字信号,再通过低分辨率DAC芯片转换为模拟信号,经射频链上变频后通过模拟预编码器进行模拟预编码,最后经发射天线发送;所述接收机通过接收天线接收接收信号,经模拟合并模块进行模拟合并,经射频链下变频后通过低分辨率ADC芯片转换为数字信号,通过接收量化模块转换为高分辨率数字信号,最后通过数字合并模块进行数字解调。本发明在保证系统性能损失在合理范围下配置了低分辨率的DAC和ADC芯片,并有效降低系统功耗、成本和复杂度。
Description
技术领域
本发明涉及无线通信领域,特别是涉及一种适用于低分辨率毫米波大规模MIMO混合预编码系统和方法。
背景技术
毫米波大规模MIMO(Multiple-Input Multiple-Output)通信技术是支撑5G(5thGeneration communication system,5G)的通信关键技术。毫米波大规模MIMO系统通常部署成百上千根天线,若再配置价格昂贵的高精度的DAC(Digital to Analog Converter,DAC)和ADC(Analog to Digital Converter,ADC)芯片,并且采用全数字射频链路架构势必导致射频功率开销大,信号处理复杂度高的问题。5G采用大规模MIMO技术。收发信机通过配置大规模的天线,在时域和频域资源内利用大规模天线的空分复用成倍地提升频谱效率,增强信号的网络覆盖范围和系统的容量,大规模MIMO技术通过依托现有的网络,可以成倍地提升网络用户的接入速率以及网络的整体容量。
5G采用毫米波通信技术。毫米波通信技术将用在高低频混合组网和链路回传。毫米波通信的一优势是毫米波频段有大量的、可用的和非授权的频谱而且比传统的低频段如,4G(4th Generation communication system,4G)无线通信频段,具有更大的带宽。毫米波通信的另一优势是毫米波波长处在毫米单位段,使用较小的天线尺寸(波长的一半)和天线间隔(波长的一半),使得几十根天线可被放置在1cm2内,从而使基站和终端侧均可在相对较小的空间里获得较大的波束赋形增益。
5G采用毫米波大规模通信技术。毫米波通信技术和大规模MIMO技术都能分别大幅度提高无线通信的系统容量。在未来5G移动通信系统中,将这两种技术相结合是非常自然的,因为现有的移动通信工作频段主要集中在6GHz以下,使得频谱资源非常拥挤,而在毫米波频段(30GHz-300GHz)可用的频谱资源丰富,能够有效缓解目前低频段信号频谱十分拥挤的现状,可以实现极高速的通信,并支持5G传输速率和容量等方面的需求。此外,毫米波频段的信号具有更短的波长,可以减少大规模天线之间的间距,十分有利用大规模天线的部署,从而提供更显著的波束赋形增益和复用增益。因此,在未来5G移动通信系统中采用毫米波大规模MIMO技术以大幅提高通信系统容量颇有潜力。
虽然理论上已证明毫米波大规模MIMO技术可以将频谱效率提高数个量级,但毫米波大规模MIMO技术在实用化过程中仍面临诸多挑战。首先,射频天线设计的难度明显增大,天线间急剧增加的复杂干扰使得基带信号更加难以处理等等。其次,对于大规模MIMO系统,基站侧部署的大规模天线阵列,通常每根天线都连接一个专门的射频链路,这个射频链路包含模数转换器或数模转换器、调制解调器、混频器、功率放大器等。传统的低频通信系统中,这些器件已发展的较为成熟。但在大规模MIMO通信系统中,这些器件会在功耗、芯片尺寸和信号处理算法等方面带来全新的挑战。虽然目前的高频段芯片封装技术已经极大地降低了电子电路成本,但是较高的功耗仍然是限制毫米波通信的主要瓶颈之一。比如,ADC(或DAC)通过采样、量化将射频信号转化为数字基带信号,其能量消耗通常与其分辨率成指数关系。目前的商用级高速(5G Sample/s)高精度(8-12比特)ADC(或DAC)的功耗大概是500mw。这意味着128根大规模天线系统中每根天线若配备两个ADC(或DAC)(I/Q两路),仅射频链路的功耗就将高达128W。
为了大幅度降低毫米波大规模MIMO系统的功耗和成本,一条可行的技术路线是在收发机采用低分辨率的DAC和ADC芯片。但是,这将不可避免带来信号的非线性量化失真问题,进而导致现有的预编码算法无法直接应用于低分辨率大规模MIMO系统中。首先,对于全数字架构的低分辨率大规模MIMO系统,预编码算法、量化函数以及量化预编码设计之间的关系仍然是未知的;其次,针对模数混合架构的低分辨率大规模MIMO系统,尚未见可行的预编码设计算法;然后,现有的预编码算法研究并未考虑收发信机均低分辨率量化的特点;最后,由于毫米波通信的电磁波趋向于在视距方向上传播,易被物体、建筑物和行人遮挡穿透性性弱,路径损耗高,大规模MIMO和毫米波结合的特殊性不能利用3GHz频谱以下的研究方法来实现,传统的预编码技术不能直接应用于毫米波大规模MIMO系统中。
目前无线通信硬件电路中,DAC和ADC芯片位数通常都大于10位,其中,16位DAC和ADC芯片主要用于基站和手机终端,由于毫米波大规模MIMO系统电路中通常配置128根天线,那么基站一侧就要128*(ADC*2+DAC*2)=512个芯片,然后由于DAC和ADC芯片功耗与其分辨率成指数关系,DAC和ADC芯片价格与其分辨率和采样频率成线性关系,所以,必须降低系统部署的成本和运营维护成本。系统配置低分辨率的DAC和ADC芯片以及混合的射频链路架构是主要的解决方法。
因此,为解决上述大规模MIMO预编码设计遇到的诸多挑战和难题,需要提供一种低分辨率毫米波大规模MIMO混合预编码系统和一种解决低分辨率DAC和ADC条件下毫米波大规模MIMO混合预编码系统容量的求解方法。
发明内容
为了解决上述问题至少之一,本发明第一方面提供一种低分辨率毫米波大规模MIMO混合预编码系统,包括:发射机和接收机,其中所述发射机包括多个数字预编码器、发送量化模块、低分辨率DAC芯片、射频链、模拟预编码器和发射天线;所述接收机包括多个接收天线、模拟合并模块、射频链、低分辨率ADC芯片、接收量化模块和数字合并模块;
所述发射机通过所述数字预编码器对发射信号进行数字预编码,通过所述发送量化模块将所述发射信号的数字信号转换为低分辨率数字信号,通过所述低分辨率DAC芯片将所述发射信号的数字信号转换为模拟信号,经所述射频链上变频后通过所述模拟预编码器进行模拟预编码,最后经所述发射天线发送;
所述接收机通过所述接收天线接收接收信号,通过所述模拟合并模块对所述接收信号进行模拟合并,经所述射频链下变频后通过所述低分辨率ADC芯片将所述接收信号的模拟信号转换为数字信号,通过所述接收量化模块将所述接收信号的数字信号转换为高分辨率数字信号,最后通过所述数字合并模块进行数字解调。
进一步地,所述射频链与所述发射天线或接收天线的连接方式为全连接阵列混合架构,即一个射频链能够驱动多个所述发射天线或接收天线以减少所述射频链的个数。
进一步地,所述模拟预编码器为模拟移相器。
进一步地,所述低分辨率DAC芯片的位数应大于等于所述低分辨率ADC芯片的位数。
本发明第二方面提供一种利用第一方面所述的预编码系统的预编码方法,包括:
S101:根据毫米波信道特点使用Saleh-Valenzuela模型建立毫米波信道矩阵;
S103:根据奇异值分解算法分解所述毫米波信道矩阵,分别计算发射机的理想模拟混合预编码Fopt和接收机的理想模拟混合预编码Wopt;
S105:根据所述Fopt和Wopt,使用正交匹配追踪算法计算发射机的模拟预编码FRF和数字预编码FBB,以及接收机的模拟预编码WRF和数字预编码WBB;
S107:根据加性量化噪声模型分别计算发射机低分辨率DAC条件下的量化函数和接收机低分辨ADC条件下的量化函数;
S109:推导出低分辨率DAC和ADC条件下的毫米波大规模MIMO量化混合预编码系统的表达式。
进一步地,所述发射机的数字基带信号是QPSK信号并且服从高斯分布。
进一步地,S105中使用正交匹配追踪算法计算发射机的模拟预编码FRF和数字预编码FBB,具体包括:
S1051:初始化残差矩阵,设置稀疏度并定义字典;
S1053:根据所述残差和字典求解索引序号,求解索引序号和字典最大相关向量;
S1055:使用新选择的字典求解所述模拟预编码FRF,更新索引模拟预编码FRF;
S1057:求解所述数字预编码FBB;
S1059:更新残差矩阵。
进一步地,S107中根据加性量化噪声模型计算发射机低分辨率DAC条件下的量化函数具体包括:
S1071:根据加性量化噪声原理分解量化映射输出矢量
S1073:求解因低分辨率DAC引起的量化失真噪声的方差
其中,QT(FBBs)表示DAC量化函数,η为低分辨率DAC造成的量化误差因子,s为发射数据流,d为DAC量化失真噪声,diag(.)表示取矩阵的对角元素,表示矩阵FBB的共轭转置运算。
进一步地,S107中根据加性量化噪声模型计算接收机低分辨率ADC条件下的量化函数具体包括:
S2071:根据加性量化噪声原理分解量化映射输出矢量
S2073:求解因低分辨率ADC引起的量化失真噪声的方差
其中,QR(rH)表示ADC量化函数,rH表示量化之前的接收信号,rL表示量化之后的接收信号,μ为低分辨率ADC造成的量化误差因子,q为ADC量化失真噪声,y为接收机天线接收的接收信号。
进一步地,所述低分辨率DAC芯片和所述低分辨率ADC芯片的位数大于等于5。
本发明的有益效果如下:
本发明提供的一种低分辨率毫米波大规模MIMO混合预编码系统和方法针对现有大规模MIMO量化预编码设计的存在的问题,应用SVD和OMP算法推出系统具体公式并分析系统性能,在保证系统性能损失在合理范围下配置了低分辨率的DAC和ADC芯片,并有效降低系统功耗、成本和复杂度,对于大规模MIMO系统预编码的研究具有重要意义。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明。
图1示出本发明的一个实施例所述的低分辨率毫米波大规模MIMO混合预编码系统框图;
图2示出本发明的一个实施例所述的低分辨率毫米波大规模MIMO混合预编码系统结构图;
图3示出本发明的一个实施例所述的低分辨率毫米波大规模MIMO混合预编码方法的流程图;
图4示出本发明的一个实施例所述的使用正交匹配追踪算法计算发射机的模拟预编码FRF和数字预编码FBB的流程图;
图5示出本发明的一个实施例所述的根据加性量化噪声模型计算发射机低分辨率DAC条件下的量化函数的流程图;
图6示出本发明的一个实施例所述的根据加性量化噪声模型计算接收机低分辨率ADC条件下的量化函数的流程图;
图7示出本发明所述混合预编码方法的效果对比图。
具体实施方式
为了更清楚地说明本发明,下面结合优选实施例和附图对本发明做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解,下面所具体描述的内容是说明性的而非限制性的,不应以此限制本发明的保护范围。
本发明的一个实施例提供了一种低分辨率毫米波大规模MIMO混合预编码系统,包括发射机和接收机,其中所述发射机包括多个数字预编码器、发送量化模块、低分辨率DAC芯片、射频链、模拟预编码器和发射天线;所述接收机包括多个接收天线、模拟合并模块、射频链、低分辨率ADC芯片、接收量化模块和数字合并模块;所述发射机通过所述数字预编码器对发射信号进行数字预编码,通过所述发送量化模块将所述发射信号的数字信号转换为低分辨率数字信号,通过所述低分辨率DAC芯片将所述发射信号的数字信号转换为模拟信号,经所述射频链上变频后通过所述模拟预编码器进行模拟预编码,最后经所述发射天线发送;所述接收机通过所述接收天线接收接收信号,通过所述模拟合并模块对所述接收信号进行模拟合并,经所述射频链下变频后通过所述低分辨率ADC芯片将所述接收信号的模拟信号转换为数字信号,通过所述接收量化模块将所述接收信号的数字信号转换为高分辨率数字信号,最后通过所述数字合并模块进行数字解调。
如图1所示,发射机采用低分辨率DAC芯片,接收机采用低分辨率ADC芯片。对于单用户场景,发射机端配置个射频通过全连接阵列混合架构驱动Nt个发射天线。其中,基带处理器和射频链之间采用低分辨率DAC芯片,发射机端把Ns数据流发射出去之前需要对信号做混合预编码预处理,首先在数字基带域通过数字预编码器进行数字预编码FBB处理,通过发送量化模块将经过数字预编码后的发射信号转换为低分辨率数字信号,通过低分辨率DAC芯片将数字信号转换为模拟信号,然后在模拟域经射频链上变频后通过模拟预编码器进行模拟预编码FRF处理。
在一个具体的示例中,在发射机端,发射天线发送信号为基带信号数据流是Ns×1维列向量且基带数字信号为零均值,方差为的QPSK信号,为发射机的发送量化模块的量化映射函数,是对FBBs高分辨率数字信号做量化映射使其变成低分辨率数字信号,从而使系统可以采用低分辨率DAC芯片,FBBs是经变成低分辨率数字信号,发射机端的混合预编码F=FRFFBB是由模拟域的维模拟预编码FRF和数字域的维数字预编码FBB构成。
发射机到接收机之间的传输信道为毫米波信道矩阵H,且传输信号受到的干扰噪声为零均值、单位方差的复高斯加性噪声根据上述表述,假设毫米波大规模MIMO系统一个窄带快衰落毫米波传播信道传输信息,则系统在接收机端的系统模型可以表示为
类似的,接收机端配置个射频链也通过全连接阵列混合架构Nr接收天线,其中,射频链路和基带处理器之间采用低分辨率的ADC芯片。接收机在接收到信号之后,首先,在模拟域通过模拟合并模块进行模拟信号的模拟合并WRF;其次,经过量化映射函数变成低分辨率数字信号;经低分辨率ADC芯片将所述接收信号转换为数字信号,通过接收量化模块将所述接收信号转换为稿分辨率十足自信号。最后,在数字域通过数字合并模块进行WBB数字信号解调。毫米波大规模MIMO系统中为了保证多数据流的传输,限定因此,毫米波大规模MIMO系统中,接收机端的信号处理系统公式为
在接收机端,WRF是模拟域维合并矩阵,为接收机量化映射函数,WBB是数字域维合并矩阵。
综上所述,毫米波大规模MIMO系统,在发射机采用低分辨率DAC芯片和接收机采用低分辨率ADC芯片条件下,收发信机系统模型可以表示为
如图2所示为收发信机结构,在发射机端,首先,做数字预编码FBB成为xH=FBBs,其次,xH经量化映射函数变成然后,xL再做模拟预编码FRF变成发射向量最后,发射向量x通过空中的毫米波信号介质传输到接收机处;接收机端的接收到的信号为y=Hx+n,其中,n为发射信号x在空中受到的高斯噪声这里,接收端的信号y在天线阵列处做模拟合并WRF变成rH,xH通过量化映射函数分解为低分辨率信号矢量rL和量化失真向量q,rL做基带解调WBB恢复
进一步地,在毫米波大规模MIMO系统中,收发信机端的射频链路都是经模拟移相器与天线阵列连接,使用模拟移相器实现模拟预编码器,收发信机中模拟移相器仅可以改变信号的相位而不能改变信号的幅度,因此矩阵FRF和WRF里的所有的元素均匀应满足|(FRF)ij|=|(WRF)ij|。
更进一步地,收发信机的射频链与天线连接的方式为全连接阵列混合架构下,即收发信机只需配置多个射频链,不需要每根天线配置一个射频链,但是,收发信机每一个射频链必须与所有的天线连接。这种架构可以使毫米波大规模MIMO系统减少大量的射频链,由于每根收发信机天线与所有射频链相连,其发送信号是所有射频链信号经过模拟移相网络移相后的叠加。由于收发信机的每条射频链要驱动所有天线,其所需的模拟移相器数等于发射或接收天线与射频链数的乘积。全连接阵列混合架构相比于全数字直接连接架构可以有效降低系统能量消耗和硬件实现复杂度,系统性能高。
更进一步地,在所述低分辨率毫米波大规模MIMO混合预编码系统中,发射机配置DAC分辨率要尽可能的高于接收机ADC的分辨率。
如图3所示,本发明的另一个实施例提供一种利用上述预编码系统的预编码方法,包括:
S101:根据毫米波信道特点使用Saleh-Valenzuela模型建立毫米波信道矩阵;
由于毫米波传播路径中的高自由空间路径损耗和有限散射数,毫米波信道将不再服从常规瑞利衰落。在本发明中,将采用Saleh-Valenzuela模型为毫米波信道分布模型。因此,毫米波信道矩阵可以表示为
其中,Ncl和Nray分别表示散射簇的数量和每个散射簇中的传播路径的数量,表示归一化因子。表示第l条散射簇的第i条传播路径的增益因子,假定其服从复高斯分布,并且为了限定其也满足此外,和分别表示传播信号的到达(Angles of Arrival,AoAs)和离开(Angles of Departure,AoDs)的方位角(Azimuth)和仰角(Elevation)。矢量和分别表示AoAs和AoDs的归一化响应向量。和分别表示AoAs和AoDs的接收和发射天线阵列增益向量。为了便于简单运算,规定和均为单位矩阵。在下列仿真分析,发射天线和接收天线均为线性阵列天线,其相应阵列响应矢量公式可以表示为
其中,d表示天线之间的天线间隔,λ为毫米波频段的信号的波长,N表示线性阵列天线的天线数量,0≤n≤N是天线在2D平面中的y方向的序号。
S103:根据奇异值分解算法分解所述毫米波信道矩阵,分别计算发射机的理想模拟混合预编码Fopt和接收机的理想模拟混合预编码Wopt;
利用奇异值分解(Singular Value Decomposition,SVD)算法分解毫米波信道矩阵H,奇异值分解算法为[U,∑,V]=SVD(H),求得理想模拟混合预编码Fopt和Wopt,其中,Fopt和Wopt分别由Nt×rank(H)维酉矩阵V前Ns个列构成的和Nr×rank(H)维酉矩阵U的前Ns个列构成。
S105:根据所述Fopt和Wopt,使用正交匹配追踪算法计算发射机的模拟预编码FRF和数字预编码FBB,以及接收机的模拟预编码WRF和数字预编码WBB。
根据所述Fopt和Wopt,以最小均方误差准则分解Fopt和Wopt。
由于模拟预编码和组合器都是由模拟移相器实现的,因此FRF和WRF具有相同的模值,因此,通过凸优化的方法解决FRF和WRF求解问题,并找到一个确切的解决方案是不可能。
在本实施例中,发射端利用正交匹配追踪(OMP)算法分别求解模拟预编码FRF和数字预编码FBB。接收端利用OMP算法分别求解模拟预编码WRF和数字预编码WBB。
OMP算法优化公式如下
其中,D是完备的字典,这个字典通常是已经给定的集合,Y是原始信号矩阵,X是待求矩阵。OMP算法思想的本质是:以贪婪迭代的方法选择字典D的列,每一次,让残差向量与已经选出来的所有向量都正交。首先,使得在每次迭代的过程中所选择的字典D的列与当前残差向量的相关性最大,然后,从原始信号向量中减去相关部分并反复迭代这一过程,只到迭代次数满足稀疏度K的要求,最后,停止迭代。
如图4所示,使用正交匹配追踪算法计算发射机的模拟预编码FRF和数字预编码FBB,具体包括:
S1051:初始化残差矩阵,设置稀疏度并定义字典;
初始化残差矩阵Fres=Fopt,稀疏度为Λt为字典。
S1053:根据所述残差和字典求解索引序号,求解索引序号和字典最大相关向量;
求索引序号和字典最大相关向量,
S1055:使用新选择的字典求解所述模拟预编码FRF,更新索引模拟预编码FRF;
更新索引模拟预编码,将新所选择的字典的第i列向量组合到FRF=[FRF|Λt(k)]。
S1057:求解所述数字预编码FBB;
求解数字预编码,
S1059:更新残差矩阵。
更新残差矩阵Fres=(Fopt-FRFFBB)/||Fopt-FRFFBB||F。
本领域技术人员应当理解,接收机预编码WRF和WBB的求解过程和发射机预编码的求解过程类似,这里不再赘述。
然后将上述OMP算法的思想直接应用于预编码求解和设计中,在毫米波大规模MIMO系统中,混合预编码和解码的设计问题的具有类似的解决方案,除了在发射机的预编码求解的过程中存在额外的功率限制,其余过程基本类似。
因此,发射机混合预编码优化问题可以表示为
s.t.FRF∈{at,1,…at,L}
其中,Fopt是最佳的无约束全数字预编码,而,FRF和FBB是待求解的模拟预编码和数字预编码器进行优化,L=NclNray是基矢量at,1,...at,L的个数,基矢量at,1,...at,L可以从发射阵列响应矢量at中矢量集合中选择。{at,1,...at,L}是模拟预编码备选的基矢量集合。是基站处的发射功率限制。
因此,接收机混合预编码优化问题可以表示为
s.t.WRF∈{ar,1,...ar,L} (8)
其中,Wopt是最佳的无约束全数字预编码,而,WRF和WBB是待求解的模拟预编码和数字预编码器进行优化,L=NclNray是基矢量ar,1,...ar,L的个数,基矢量ar,1,...ar,L可以从发射阵列响应矢量ar中矢量集合中选择。{ar,1,...ar,L}是模拟预编码备选的基矢量集合。
S107:根据加性量化噪声模型分别计算发射机低分辨率DAC条件下的量化函数和接收机低分辨ADC条件下的量化函数;
在加性量化噪声模型(AQNM)中,量化映射函数的输入信号复信号z,其通过所述的低分辨率量化函数输出为复信号表示加性量化函数。令量化器输出以使得因此,信号通过加性量化函数模型可以由下列公式表示:
zq=εz+nq (9)
其中,nq是加性量化噪声,并且z和zq是不相关的。
由此可得,加性量化噪声的均值nq为
由此可得,加性量化噪声的方差nq为
上述式子中,ε的计算公式为
ε=1-μ (12)
其中,为失真因子,为加性量化函数输入的方差,加性量化误差eq的方差,eq=z-zq。
如图5所示,根据加性量化噪声模型计算发射机低分辨率DAC条件下的量化函数具体包括:
S1071:根据加性量化噪声原理分解量化映射输出矢量
在本实施例中应用AQNM分解发射机量化映射函数并推导公式。由于应用AQNM要求输入信号必须是高斯分布的,因此,本实施例中发射机输出的信号是服从高斯分布的。
根据上述表述AQNM阐述了高斯信号在通过非线性映射函数之前和之后的关系,它将量化映射信号分解成两个不相关的部分:一个是关于输入量化映射信号矢量的线性函数,另一个是与量化映射的输入矢量不相关的失真项。
发射机的输入向量经过发送量化模块的量化映射函数应用AQNM,则量化映射输出矢量可以分解为量化映射函数:
其中,η为量化失真因子,不同分辨率或位数的DAC芯片,量化失真因子η不同,具体为DAC芯片的位数越高,量化失真因子η就越小,d表示量化失真向量。在发射机端由于量化映射函数的非线性,非量化预编码矢量与量化预编码矢量之间会产生量化误差因此,发射机必须最小化量化误差,以使量化失真向量对系统性能的影响最低化。
S1073:求解因低分辨率DAC引起的量化失真噪声的方差
根据上述AQNM原理,发射机量化失真向量d的方差为
根据上述AQNM原理,量化失真因子η为
因此,根据公式(13),公式(1)可化为:
y=(1-η)HFRFFBBs+HFRFd+n (16)
类似的,如图6所示,根据加性量化噪声模型计算接收机机低分辨率ADC条件下的量化函数具体包括:
S2071:根据加性量化噪声原理分解量化映射输出矢量
在接收机中应用AQNM分解并推导公式。由于发射机使用的输入信号s服从高斯分布,发射机端的高分辨率信号xH=FBBs经过低分辨率量化分解为xL=GFBBs和量化失真向量q,xL经低分辨率DAC芯片转成模拟信号atl再通过射频调制变成att由天线发射出,att受低分辨率DAC芯片影响以及受空中高斯噪声的影响,其波形不再非常平稳光滑,但是只要DAC的位数大于5,由DAC输出的模拟信号仍然平稳光滑,可以认为模拟信号att仍然是服从高斯分布,由此可知接收机接收到的信号y也是服从高斯分布的。因此就满足接收机应用AQNM要求,即接收机收到的接收信号必须是高斯分布的条件。
接收机的输入向量通过接收量化模块的量化映射函数应用AQNM,则量化映射输出矢量rH可以分解为
S2073:求解因低分辨率ADC引起的量化失真噪声的方差
根据上述AQNM原理,接收机量化失真向量q的方差为
根据上述AQNM原理,量化失真因子μ为
S109:推导出低分辨率DAC和ADC条件下的毫米波大规模MIMO量化混合预编码系统的表达式。
发射机和接收机应用AQNM分解接收端量化映射函数从而设计毫米波大规模MIMO量化混合预编码系统。其中,接收机解调信号为,接收机因发射机配置低分辨率DAC芯片引入量化噪声为,接收机因配置低分辨率ADC芯片引入量化噪声为,发射信号x因传播环境干扰引入的噪声为,
收发信机均应用AQNM分解接收端量化映射函数从而设计毫米波大规模MIMO量化混合预编码系统,因此,系统公式可以具体表达为
低分辨率ADC/DA的位数和量化失真因子的关系主要有公式来决定。1位ADC条件下μ1=0.3634,2位ADC条件下μ2=0.1175,当ADC的位数的b大于2时,量化失真因子可以为
毫米波大规模MIMO量化混合预编码系统的因低分辨率DAC引起的量化失真噪声的方差为
毫米波大规模MIMO量化混合预编码系统的因低分辨率ADC引起量化失真向量噪声的方差为
毫米波大规模MIMO量化混合预编码系统因噪声引起的等效高斯噪声其方差为为
因此,毫米波大规模MIMO量化混合预编码系统的系统容量的表达式可以为
R=log2(1+SDNR) (25)
其中,R为传输速率,SDNR为信号噪声比。
如图7所示为本发明所述混合预编码方法的效果对比图,从图中可知,OMP HybridPrecoder是系统未采用低分辨率DAC和ADC芯片,预编码采用OMP算法状态下频谱效率曲线,1-bit DAC and 1-bit ADC、3-bit DAC and 3-bit ADC、5-bit DAC and 5-bit ADC,6-bitDAC and 6-bit ADC仿真曲线分别是系统配置1位、3位、5位、6位DAC和ADC芯片并且利用AQNM分解量化映射函数并且应用OMP算法设计量化混合预编码系统下的频谱效率曲线。1-bit DAC、3-bit DAC、5-bit DAC仿真曲线分别是系统配置1位、3位、5位DAC且利用AQNM分解量化映射函数并且应用OMP算法设计量化混合预编码下的频谱效率曲线。1-bit ADC、3-bitADC、5-bit ADC仿真曲线分别是系统配置1位、3位、5位ADC芯片并且利用AQNM分解量化映射函数并且应用OMP算法设计量化预编码系统下的频谱效率曲线,由图中可知,单独使用位数大于等于5的DAC芯片或单独使用位数大于等于5的ADC芯片时能够满足系统性能的要求。
因此,在本实施例所述的低分辨率毫米波大规模MIMO混合预编码系统中,单独使用发射机的低分辨率DAC芯片的位数为5时,单独使用接收机的低分辨率ADC芯片的位数为5时,同时使用收发信机的低分辨率DAC芯片和低分辨率ADC芯片为6时,能够确保系统的性能在可接受的范围内时实现系统功耗、成本和复杂度的显著降低;当DAC芯片和ADC芯片同时存在收发信机中,发射机配置DAC的分辨率要尽可能的高于接收机ADC的分辨率。
本发明针对现有大规模MIMO量化预编码设计的存在的问题,在研究全连接和低分辨率的DAC和ADC的实际条件下,在理论上根据AQNM和Bussgang原理分析量化映射函数,应用SVD和OMP算法推出系统具体公式并分析系统性能,对于大规模MIMO系统预编码的研究具有重要意义。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。
Claims (10)
1.一种低分辨率毫米波大规模MIMO混合预编码系统,其特征在于,包括:发射机和接收机,其中
所述发射机包括多个数字预编码器、发送量化模块、低分辨率DAC芯片、射频链、模拟预编码器和发射天线;
所述接收机包括多个接收天线、模拟合并模块、射频链、低分辨率ADC芯片、接收量化模块和数字合并模块;
所述发射机通过所述数字预编码器对发射信号进行数字预编码,通过所述发送量化模块将所述发射信号的数字信号转换为低分辨率数字信号,通过所述低分辨率DAC芯片将所述发射信号的数字信号转换为模拟信号,经所述射频链上变频后通过所述模拟预编码器进行模拟预编码,最后经所述发射天线发送;
所述接收机通过所述接收天线接收接收信号,通过所述模拟合并模块对所述接收信号进行模拟合并,经所述射频链下变频后通过所述低分辨率ADC芯片将所述接收信号的模拟信号转换为数字信号,通过所述接收量化模块将所述接收信号的数字信号转换为高分辨率数字信号,最后通过所述数字合并模块进行数字解调。
2.根据权利要求1所述的混合预编码系统,其特征在于,所述射频链与所述发射天线或接收天线的连接方式为全连接阵列混合架构,即一个射频链能够驱动多个所述发射天线或接收天线以减少所述射频链的个数。
3.根据权利要求1所述的混合预编码系统,其特征在于,所述模拟预编码器为模拟移相器。
4.根据权利要求1所述的混合预编码系统,其特征在于,所述低分辨率DAC芯片的位数应大于等于所述低分辨率ADC芯片的位数。
5.一种应用权利要求1-4所述的任一项预编码系统的预编码方法,其特征在于,包括:
S101:根据毫米波信道特点使用Saleh-Valenzuela模型建立毫米波信道矩阵;
S103:根据奇异值分解算法分解所述毫米波信道矩阵,分别计算发射机的理想模拟混合预编码Fopt和接收机的理想模拟混合预编码Wopt;
S105:根据所述Fopt和Wopt,使用正交匹配追踪算法计算发射机的模拟预编码FRF和数字预编码FBB,以及接收机的模拟预编码WRF和数字预编码WBB;
S107:根据加性量化噪声模型分别计算发射机低分辨率DAC条件下的量化函数和接收机低分辨ADC条件下的量化函数;
S109:推导出低分辨率DAC和ADC条件下的毫米波大规模MIMO量化混合预编码系统的表达式。
6.根据权利要求5所述的混合预编码方法,其特征在于,所述发射机的数字基带信号是QPSK信号并且服从高斯分布。
7.根据权利要求5所述的混合预编码方法,其特征在于,S105中使用正交匹配追踪算法计算发射机的模拟预编码FRF和数字预编码FBB,具体包括:
S1051:初始化残差矩阵,设置稀疏度并定义字典;
S1053:根据所述残差和字典求解索引序号,求解索引序号和字典最大相关向量;
S1055:使用新选择的字典求解所述模拟预编码FRF,更新索引模拟预编码FRF;
S1057:求解所述数字预编码FBB;
S1059:更新残差矩阵。
8.根据权利要求5所述的混合预编码方法,其特征在于,S107中根据加性量化噪声模型计算发射机低分辨率DAC条件下的量化函数具体包括:
S1071:根据加性量化噪声原理分解量化映射输出矢量
S1073:求解因低分辨率DAC引起的量化失真噪声的方差
其中,QT(FBBs)表示DAC量化函数,η为低分辨率DAC造成的量化误差因子,s为发射数据流,d为DAC量化失真噪声,diag(.)表示取矩阵的对角元素,表示矩阵FBB的共轭转置运算。
9.根据权利要求5所述的混合预编码方法,其特征在于,S107中根据加性量化噪声模型计算接收机低分辨率ADC条件下的量化函数具体包括:
S2071:根据加性量化噪声原理分解量化映射输出矢量
S2073:求解因低分辨率ADC引起的量化失真噪声的方差
其中,QR(rH)表示ADC量化函数,rH表示量化之前的接收信号,rL表示量化之后的接收信号,μ为低分辨率ADC造成的量化误差因子,q为ADC量化失真噪声,y为接收机天线接收的接收信号。
10.根据权利要求5所述的混合预编码方法,其特征在于,所述低分辨率DAC芯片和所述低分辨率ADC芯片的位数大于等于5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810151431.6A CN108712196B (zh) | 2018-02-14 | 2018-02-14 | 低分辨率毫米波大规模mimo混合预编码系统和方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810151431.6A CN108712196B (zh) | 2018-02-14 | 2018-02-14 | 低分辨率毫米波大规模mimo混合预编码系统和方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108712196A true CN108712196A (zh) | 2018-10-26 |
CN108712196B CN108712196B (zh) | 2021-04-09 |
Family
ID=63866174
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810151431.6A Expired - Fee Related CN108712196B (zh) | 2018-02-14 | 2018-02-14 | 低分辨率毫米波大规模mimo混合预编码系统和方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108712196B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110224730A (zh) * | 2019-06-03 | 2019-09-10 | 东南大学 | 一种毫米波通信的混合预编码结构、混合合并结构及方法 |
CN110350963A (zh) * | 2019-08-01 | 2019-10-18 | 湖南国科锐承电子科技有限公司 | 毫米波mimo通信系统中波束成形的方法和系统 |
CN111342874A (zh) * | 2020-03-02 | 2020-06-26 | 北京大学 | 用于下行大规模宽带毫米波多用户系统的预编码方法 |
CN112073105A (zh) * | 2020-11-11 | 2020-12-11 | 华东交通大学 | 一种低能耗的毫米波mimo通信预编码设计方法 |
WO2021081848A1 (en) * | 2019-10-30 | 2021-05-06 | Nokia Shanghai Bell Co., Ltd. | Quantized precoding in massive mimo system |
CN113014296A (zh) * | 2021-02-18 | 2021-06-22 | 哈尔滨工业大学 | 全双工毫米波双向中继中基于正交匹配追踪的自干扰抑制混合预编码方法 |
CN113179109A (zh) * | 2021-04-26 | 2021-07-27 | 南京盛航海运股份有限公司 | 一种去蜂窝大规模mimo上行频谱效率优化方法 |
CN114039638A (zh) * | 2021-11-04 | 2022-02-11 | 广东工业大学 | 一种混合波束成形器与模数转换器联合设计方法 |
CN115136636A (zh) * | 2019-12-13 | 2022-09-30 | 株式会社Ntt都科摩 | 终端以及无线通信方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011155569A1 (ja) * | 2010-06-09 | 2011-12-15 | 国立大学法人九州大学 | 無線通信システム、送信機、伝搬路特性推定装置、伝搬路特性推定方法及びプログラム |
CN105610477A (zh) * | 2016-01-28 | 2016-05-25 | 西南交通大学 | 基于压缩感知的多发多收系统增强信号复用技术 |
CN105959048A (zh) * | 2016-06-23 | 2016-09-21 | 北京科技大学 | 一种大规模天线的预编码方法 |
CN106603130A (zh) * | 2016-12-20 | 2017-04-26 | 西安电子科技大学 | 一种大规模mimo系统中数模混合预编码方法 |
CN107017927A (zh) * | 2017-02-28 | 2017-08-04 | 东南大学 | 一种大规模mimo系统中基站dac精度配置方法 |
CN107359917A (zh) * | 2017-07-26 | 2017-11-17 | 东南大学 | 一种大规模mimo最优用户调度数目配置方法 |
CN107370493A (zh) * | 2017-06-08 | 2017-11-21 | 东南大学 | 低精度adc与混合预编码结合的毫米波传输方法及通信系统 |
US20180019791A1 (en) * | 2016-07-15 | 2018-01-18 | Futurewei Technologies, Inc. | Digital to analog converter apparatus, system, and method with quantization noise that is independent of an input signal |
-
2018
- 2018-02-14 CN CN201810151431.6A patent/CN108712196B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011155569A1 (ja) * | 2010-06-09 | 2011-12-15 | 国立大学法人九州大学 | 無線通信システム、送信機、伝搬路特性推定装置、伝搬路特性推定方法及びプログラム |
CN105610477A (zh) * | 2016-01-28 | 2016-05-25 | 西南交通大学 | 基于压缩感知的多发多收系统增强信号复用技术 |
CN105959048A (zh) * | 2016-06-23 | 2016-09-21 | 北京科技大学 | 一种大规模天线的预编码方法 |
US20180019791A1 (en) * | 2016-07-15 | 2018-01-18 | Futurewei Technologies, Inc. | Digital to analog converter apparatus, system, and method with quantization noise that is independent of an input signal |
CN106603130A (zh) * | 2016-12-20 | 2017-04-26 | 西安电子科技大学 | 一种大规模mimo系统中数模混合预编码方法 |
CN107017927A (zh) * | 2017-02-28 | 2017-08-04 | 东南大学 | 一种大规模mimo系统中基站dac精度配置方法 |
CN107370493A (zh) * | 2017-06-08 | 2017-11-21 | 东南大学 | 低精度adc与混合预编码结合的毫米波传输方法及通信系统 |
CN107359917A (zh) * | 2017-07-26 | 2017-11-17 | 东南大学 | 一种大规模mimo最优用户调度数目配置方法 |
Non-Patent Citations (1)
Title |
---|
张运良等: "毫米波大规模MIMO系统量化混合预编码", 《科技经济导刊》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110224730A (zh) * | 2019-06-03 | 2019-09-10 | 东南大学 | 一种毫米波通信的混合预编码结构、混合合并结构及方法 |
CN110224730B (zh) * | 2019-06-03 | 2021-06-08 | 东南大学 | 一种毫米波通信的混合预编码结构、混合合并结构及方法 |
CN110350963A (zh) * | 2019-08-01 | 2019-10-18 | 湖南国科锐承电子科技有限公司 | 毫米波mimo通信系统中波束成形的方法和系统 |
WO2021081848A1 (en) * | 2019-10-30 | 2021-05-06 | Nokia Shanghai Bell Co., Ltd. | Quantized precoding in massive mimo system |
US11855719B2 (en) | 2019-10-30 | 2023-12-26 | Nokia Solutions And Networks Oy | Quantized precoding in massive MIMO system |
CN115136636A (zh) * | 2019-12-13 | 2022-09-30 | 株式会社Ntt都科摩 | 终端以及无线通信方法 |
CN111342874A (zh) * | 2020-03-02 | 2020-06-26 | 北京大学 | 用于下行大规模宽带毫米波多用户系统的预编码方法 |
CN111342874B (zh) * | 2020-03-02 | 2021-03-23 | 北京大学 | 用于下行大规模宽带毫米波多用户系统的预编码方法 |
CN112073105A (zh) * | 2020-11-11 | 2020-12-11 | 华东交通大学 | 一种低能耗的毫米波mimo通信预编码设计方法 |
CN113014296B (zh) * | 2021-02-18 | 2022-01-04 | 哈尔滨工业大学 | 全双工毫米波双向中继中基于正交匹配追踪的自干扰抑制混合预编码方法 |
CN113014296A (zh) * | 2021-02-18 | 2021-06-22 | 哈尔滨工业大学 | 全双工毫米波双向中继中基于正交匹配追踪的自干扰抑制混合预编码方法 |
CN113179109B (zh) * | 2021-04-26 | 2022-03-29 | 南京盛航海运股份有限公司 | 一种去蜂窝大规模mimo上行频谱效率优化方法 |
CN113179109A (zh) * | 2021-04-26 | 2021-07-27 | 南京盛航海运股份有限公司 | 一种去蜂窝大规模mimo上行频谱效率优化方法 |
CN114039638A (zh) * | 2021-11-04 | 2022-02-11 | 广东工业大学 | 一种混合波束成形器与模数转换器联合设计方法 |
CN114039638B (zh) * | 2021-11-04 | 2023-05-23 | 广东工业大学 | 一种混合波束成形器与模数转换器联合设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108712196B (zh) | 2021-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108712196A (zh) | 低分辨率毫米波大规模mimo混合预编码系统和方法 | |
Ribeiro et al. | Energy efficiency of mmWave massive MIMO precoding with low-resolution DACs | |
Roth et al. | Achievable rate and energy efficiency of hybrid and digital beamforming receivers with low resolution ADC | |
CN109302224B (zh) | 用于大规模mimo的混合波束赋形算法 | |
CN106603130B (zh) | 一种大规模mimo系统中数模混合预编码方法 | |
CN108494455B (zh) | 采用单比特模数转换器mimo系统的混合预编码设计方法 | |
CN107135024A (zh) | 一种低复杂度的混合波束赋形迭代设计方法 | |
CN108933745A (zh) | 一种基于超分辨率角度和时延估计的宽带信道估计方法 | |
CN108199753A (zh) | 一种毫米波通信中基于迭代最小的预编码方法 | |
CN101170333B (zh) | 分布式多输入多输出无线通信系统真实环境测试平台 | |
CN108023620A (zh) | 应用于毫米波频段的大规模mimo系统混合预编码方法 | |
CN106452535A (zh) | 一种盲信道下的混合波束成形方法及装置 | |
CN109714091B (zh) | 一种在毫米波mimo系统中基于分层设计的迭代混合预编码方法 | |
CN108234101A (zh) | 能效最大化导频信号设计方法及大规模多天线系统 | |
CN106487725A (zh) | 一种多用户mimo系统毫米波信道估计方法 | |
CN108063634A (zh) | 一种低精度量化大规模mimo中最优正则预编码方法 | |
CN109167623A (zh) | 一种应用于毫米波多天线系统的混合波束成形结构及系统 | |
Ding et al. | Performance analysis of mixed-ADC massive MIMO systems over spatially correlated channels | |
Zhao et al. | Energy efficient hybrid beamforming for multi-user millimeter wave communication with low-resolution A/D at transceivers | |
Xiong et al. | Experimental study of plane spiral OAM mode-group based MIMO communications | |
CN108270471A (zh) | 大规模天线阵列系统中低反馈混合预编码与接收合并方法 | |
Lu et al. | Quantized hybrid precoding design for millimeter-wave large-scale MIMO systems | |
Zhang et al. | Analysis of indoor THz communication systems with finite-bit DACs and ADCs | |
Atzeni et al. | Low-resolution massive MIMO under hardware power consumption constraints | |
CN107276657B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210409 |
|
CF01 | Termination of patent right due to non-payment of annual fee |