CN114095062B - 毫米波大规模相控阵系统的宽带修正模拟波束成形方法 - Google Patents

毫米波大规模相控阵系统的宽带修正模拟波束成形方法 Download PDF

Info

Publication number
CN114095062B
CN114095062B CN202111490919.XA CN202111490919A CN114095062B CN 114095062 B CN114095062 B CN 114095062B CN 202111490919 A CN202111490919 A CN 202111490919A CN 114095062 B CN114095062 B CN 114095062B
Authority
CN
China
Prior art keywords
vector
beam forming
array
phased array
millimeter wave
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.)
Active
Application number
CN202111490919.XA
Other languages
English (en)
Other versions
CN114095062A (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.)
Chongqing Aerospace Launch Vehicle Electronic Technology Co ltd
Original Assignee
Chongqing Aerospace Launch Vehicle Electronic Technology Co ltd
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 Chongqing Aerospace Launch Vehicle Electronic Technology Co ltd filed Critical Chongqing Aerospace Launch Vehicle Electronic Technology Co ltd
Priority to CN202111490919.XA priority Critical patent/CN114095062B/zh
Publication of CN114095062A publication Critical patent/CN114095062A/zh
Application granted granted Critical
Publication of CN114095062B publication Critical patent/CN114095062B/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
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/10Monitoring; Testing of transmitters
    • H04B17/11Monitoring; Testing of transmitters for calibration
    • H04B17/12Monitoring; Testing of transmitters for calibration of transmit antennas, e.g. of the amplitude or phase
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/20Monitoring; Testing of receivers
    • H04B17/21Monitoring; Testing of receivers for calibration; for correcting measurements
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE 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/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)
  • Radio Transmission System (AREA)

Abstract

本发明涉及一种毫米波大规模相控阵系统的宽带修正模拟波束成形方法,属于无线通信技术领域。该方法包括:S1:根据信号入射角度,构造标准模拟波束成形矢量;S2:采用CAZAC序列中的Zadoff‑Chu序列构造修正矢量;S3:根据标准模拟波束成形矢量和修正矢量,构造修正模拟波束成形矢量。本发明应用于实际毫米波远距离高容量大规模相控阵通信系统接,在不增加复杂度的情况下,仅仅对标准模拟波束成形矢量进行修正。使用本发明方法所得到的空间等效信道频域响应不存在深衰落点,从而消除了波束倾斜引起的码间干扰对系统误比特率性能的影响。

Description

毫米波大规模相控阵系统的宽带修正模拟波束成形方法
技术领域
本发明属于无线通信技术领域,涉及一种毫米波大规模相控阵系统的宽带修正模拟波束成形方法。
背景技术
毫米波频段的范围为30GHz~300GHz,具有众多可用频谱,高吞吐量毫米波通信已被视为第五代(5G)蜂窝网络及以后时代的关键启用技术,可以支持超高速率无线传输。然而,毫米波频段信号的严重衰减限制了有效通信距离。幸运的是,由于天线的尺寸更小,频率更高,毫米波系统允许在一个小区域内封装大量天线,以提供显著的波束形成增益,以补偿严重的信号衰减。尽管利用大规模阵列天线进行毫米波通信有诸多好处,但同时考虑到系统性能和复杂性,要充分发挥大规模阵列天线的优势还存在诸多挑战。
考虑远程点对高吞吐量毫米波通信场景,毫米波收发机可以采用大规模相控阵进行模拟波束成形来增加通信距离。在相控阵中,所有阵元通过移相器网络连接到单个射频链路上。然而在天线规模数和信号带宽较大时,电磁波在大规模阵列孔径上的物理传播延迟可以与符号周期相当。因此各个阵元对基带信号的时延不可忽略,阵列导向矢量将与频率相关。需要注意的是,移相器的移相值与频率是不相关的。如果仍采用传统窄带标准模拟波束成形方法,会导致大规模相控阵出现宽带波束倾斜现象,即产生的空间波束在不同频率下波束指向不同。从通信角度而言,大规模相控阵宽带波束倾斜现象不仅仅会导致阵列增益损失,还会导致物理无线信道和移相器网络组成的空间等效信道变为频率选择性信道。此外,该频率选择性信道存在多个深衰落点。信道的频率选择性会导致接收机出现码间干扰,从而降低了系统误比特率性能。
为了消除宽带波束倾斜现象,常采用实时延迟线(True Time Delay,TTD)取代移相器。TTD阵列可以由波导、同轴电缆或光纤组成,但存在体积大、功耗大、成本高和受温度等环境影响等缺点。此外,全数字阵列结构每个阵元通道都单独配置一路RF链路,从而使得所有阵列数据可以在数字基带进行处理。TTD阵列和全数字阵列能够产生与频率相关的波束成形矢量,因此能够消除波束倾斜现象。但与传统相控阵相比,会导致毫米波大规模天线系统的复杂度和功耗急剧增加。基于大规模天线阵列宽带系统可实现角度,需设计低复杂度的宽带模拟波束方法以降低相控阵宽带波束倾斜现象引起的性能损失。
发明内容
有鉴于此,本发明的目的在于提供一种宽带修正模拟波束成形方法,应用于实际毫米波远距离高容量大规模相控阵通信系统接,在不增加复杂度的情况下,仅仅对标准模拟波束成形矢量进行修正。使用本发明方法所得到的空间等效信道频域响应不存在深衰落点,从而消除了波束倾斜引起的码间干扰对系统误比特率性能的影响。
为达到上述目的,本发明提供如下技术方案:
一种毫米波大规模相控阵系统的宽带修正模拟波束成形方法,基于恒包络零自相关(CAZAC)序列设计一种修正的宽带模拟波束成形方法,从而消除物理无线信道和移相器网络组成的空间等效信道中由波束倾斜引起的深衰落点。本质上,空间等效信道频域响应相当于模拟波束成形矢量关于天线空间采样的离散傅里叶变换。为使空间等效信道尽可能平坦,可以寻找一种在频域具有平坦特性的序列来对模拟波束成形矢量进行修正。CAZAC序列即是满足该要求的序列。本发明使用CAZAC序列中的Zadoff-Chu序列对标准模拟波束成形矢量进行修正。参阅图1,该方法具体包括以下步骤:
S1:根据信号入射角度θ,构造标准模拟波束成形矢量wstd
S2:采用CAZAC序列中的Zadoff-Chu序列构造修正矢量wcor
S3:根据标准模拟波束成形矢量wstd和修正矢量wcor,构造修正模拟波束成形矢量w。
进一步,步骤S1中,构造的标准模拟波束成形矢量wstd为:
Figure BDA0003399315310000021
Figure BDA0003399315310000022
m=1,…,M
其中,fc表示载波频率,d表示阵列单元间距,θ0表示期望波束指向,c表示光速,M表示天线数目,φstd,m为第m个天线对应的移相器的移相值。
进一步,步骤S2中,采用Zadoff-Chu序列构造的修正矢量wcor为:
Figure BDA0003399315310000023
Figure BDA0003399315310000024
m=1,…,M
其中,B表示信号频带带宽,φcor,m为第m个天线对应的移相器的修正相位。
进一步,步骤S3中,构造的修正模拟波束成形矢量W为:
w=wstd⊙wcor
其中,⊙为向量Hadamard积操作。
本发明的有益效果在于:本发明主要应用于远距离高容量毫米波大规模相控阵通信系统,其优势在于:
1)与传统的相控阵标准模拟波束成形方法相比,本发明能够消除标准模拟波束成形方法下的空间等效信道中由波束倾斜引起的深衰落点,从而在不增加复杂度的情况下改善毫米波通信系统的误比特率性能。
2)与TTD阵列和全数字阵列宽带波束成形方法相比,本发明降低了系统的成本、功耗和复杂度。本发明方法仍然使用模拟相控阵结构,并不需要在每个阵元通道使用实时延迟线和RF链路,而是直接修正移相器网络各移相器的移相值消除宽带影响。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。
附图说明
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作优选的详细描述,其中:
图1为本发明宽带修正模拟波束成形方法的流程图;
图2为本发明修正模拟波束成形方法与传统相控阵模拟波束成形方法(标准模拟波束成形方法)在60GHz毫米波通信系统下的空间等效信道频域响应对比效果图;
图3为60GHz毫米波单载波频域均衡(SC-FDE)通信系统中大规模相控阵使用本发明修正模拟波束成形方法、传统相控阵模拟波束成形方法(标准模拟波束成形方法)与理想TTD阵列或全数字阵列(理想曲线)的误比特率性能对比效果图。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需要说明的是,以下实施例中所提供的图示仅以示意方式说明本发明的基本构想,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
为方便理解,首先对本发明实施例考虑的宽带毫米波视距单输入多输出(LoS-SIMO)系统进行简要介绍:
本实施例考虑一个宽带毫米波LoS-SIMO系统。接收机部署大规模相控阵形成单个空间波束指向固定方位单天线发射机。空间无线信道假设为单径视距(LoS)信道。接收机天线阵列为均匀分布直线阵列(ULA),天线数目为M,阵列单元间距为d。假设发射机发送信号从远场入射到接收机相控阵列,入射角度为θ∈[-π/2,π/2]。在不考虑信号衰减,并假定所有阵元为理想全向天线的情况下,阵列接收机接收的频域信号矢量r(f)可以表示为:
r(f)=a(θ,f)s(f)+n(f)
Figure BDA0003399315310000041
其中,s(f)为频域等效基带远场信号,f∈[-B/2,B/2],B表示基带信号带宽;矢量n(f)为接收机的高斯白噪声,矢量a(θ,f)为宽带频域阵列导向矢量,由信号频率和方向共同决定。参数τ=dsinθ/c表示相邻阵元间时延差,c为光速。
接收机使用移相器网络进行模拟波束成形,将阵列接收信号进行合并。移相器网络中每个天线通道对应的移相器移相值与频率无关。在标准模拟波束成形方法中,移相器的移相值根据载波频率fc设计。如果期望波束指向为θ0,则移相器网络的模拟波束成形矢量为:
Figure BDA0003399315310000042
Figure BDA0003399315310000043
m=1,…,M
接收机使用模拟波束成形矢量w得到的频域基带合并信号可以表示为:
Figure BDA0003399315310000044
从通信系统角度考虑,可以将模拟波束成形和物理无线信道看做为一个空间等效信道。则该空间等效信道的频域响应为H(f)=wHa(θ,f)。假设接收机知道信号入射角度θ,并设置θ0=θ。如果使用标准模拟波束成形矢量wstd,则得到的空间等效信道的频域响应为:
Figure BDA0003399315310000045
当天线数目M和信号带宽B较大时,使用标准模拟波束成形方法得到的空间等效信道频域响应将出现频率选择性,并存在多个深衰落点。从而引起接收机码间干扰,使得系统误比特率性能下降。
其中,修正宽带模拟波束成形方法,主要包括以下内容:
1)根据信号入射角度θ,构造标准模拟波束成形矢量wstd
2)对于第m个天线对应的移相器,在标准波束成形方法下的移相值φstd,m基础上添加修正相位φcor,m,得到对应的修正波束成形加权值为:
Figure BDA0003399315310000051
m=1,…,M
相应地,修正后的模拟波束成形矢量为:
w=wstd⊙wcor
Figure BDA0003399315310000052
Figure BDA0003399315310000053
m=1,…,M
其中,⊙为向量Hadamard积操作,wcor表示修正矢量。使用修正模拟波束成形矢量得到的空间等效信道频域响应可以表示为:
Figure BDA0003399315310000054
3)修正模拟波束成形矢量下得到的空间等效信道频域响应H(f)相当于修正矢量关于天线空间采样的离散傅里叶变换。因此,可以通过设计
Figure BDA0003399315310000055
m=1,…,M使得等效信道频域响应尽可能平坦。在本发明中,采用CAZAC序列中的Zadoff-Chu序列设计修正矢量。
定义s=fdsinθ/c,则有s∈[-Bdsinθ/2c,-Bdsinθ/2c]。对函数H(s)进行M点均匀采样,得到
Figure BDA0003399315310000056
m=0,…,M-1的离散傅里叶变换表达式,即
Figure BDA0003399315310000057
k=0,1,…,M-1
严格来讲,应该讨论变量s在范围[-Bdsinθ/2c,-Bdsinθ/2c]内的连续变化情况。不过当采样数目M(天线数目)很大时,离散傅里叶变换能够很好地近似连续傅里叶变换变化。因此,Hk,k=0,1,…,M-1近似函数H(s)是可信的。接下来,使用Zadoff-Chu序列设计矢量
Figure BDA0003399315310000058
m=0,…,M-1,即
Figure BDA0003399315310000061
本发明主要考虑天线数目M为偶数的情况。对变量
Figure BDA0003399315310000062
的相位进行求导,可以得到其频域范围为s∈[rq/M,r(M-1+q)/M]。已知期望频率范围为s∈[-Bdsinθ/2c,-Bdsinθ/2c]。则可以得到参数r和q别为:
Figure BDA0003399315310000063
根据参数r和q可以得到修正相位的具体表达式为:
Figure BDA0003399315310000064
m=0,…,M-1
最后,使用该修正相位对标准模拟波束成形矢量进行修正得到修正后的模拟波束成形矢量。
对比实验:
图2为本发明修正模拟波束成形方法与传统相控阵模拟波束成形方法(标准模拟波束成形方法)在60GHz毫米波通信系统下的空间等效信道频域响应对比效果图,其中天线数目M=256,信号入射角度θ=30°。图3为60GHz毫米波单载波频域均衡(SC-FDE)通信系统中大规模相控阵使用本发明修正模拟波束成形方法、传统相控阵模拟波束成形方法(标准模拟波束成形方法)与理想TTD阵列或全数字阵列(理想曲线)的误比特率性能对比;其中天线数目M=256,信号入射角度θ=30°,接收机均衡算法采用迫零(ZF)算法或最小均方误差(MMSE)算法。从图2可以看出,使用本发明所提出的修正宽带模拟波束成形方法所得到的空间等效信道频域响应不存在深衰落点,消除了波束倾斜引起的码间干扰对系统误比特率性能的影响。从图3可以看出,本发明方法可应用于实际毫米波远距离高容量大规模相控阵通信系统接。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (1)

1.一种毫米波大规模相控阵系统的宽带修正模拟波束成形方法,其特征在于,该方法基于一个宽带毫米波LoS-SIMO系统,接收机部署大规模相控阵形成单个空间波束指向固定方位单天线发射机;空间无线信道假设为单径视距信道;接收机天线阵列为均匀分布直线阵列,天线数目为M,阵列单元间距为d;假设发射机发送信号从远场入射到接收机相控阵列,入射角度为θ∈[-π/2,π/2];接收机使用移相器网络进行模拟波束成形,将阵列接收信号进行合并;该方法具体包括以下步骤:
S1:根据信号入射角度θ,构造标准模拟波束成形矢量wstd
Figure FDA0004067751630000011
Figure FDA0004067751630000012
其中,fc表示载波频率,d表示阵列单元间距,θ0表示期望波束指向,c表示光速,M表示天线数目,φstd,m为第m个天线对应的移相器的移相值;
S2:采用CAZAC序列中的Zadoff-Chu序列构造修正矢量wcor
Figure FDA0004067751630000013
Figure FDA0004067751630000014
其中,B表示信号频带带宽,φcor,m为第m个天线对应的移相器的修正相位;
S3:根据标准模拟波束成形矢量wstd和修正矢量wcor,构造修正模拟波束成形矢量W;
w=wstd⊙wcor
其中,⊙为向量Hadamard积操作。
CN202111490919.XA 2021-12-08 2021-12-08 毫米波大规模相控阵系统的宽带修正模拟波束成形方法 Active CN114095062B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111490919.XA CN114095062B (zh) 2021-12-08 2021-12-08 毫米波大规模相控阵系统的宽带修正模拟波束成形方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111490919.XA CN114095062B (zh) 2021-12-08 2021-12-08 毫米波大规模相控阵系统的宽带修正模拟波束成形方法

Publications (2)

Publication Number Publication Date
CN114095062A CN114095062A (zh) 2022-02-25
CN114095062B true CN114095062B (zh) 2023-03-14

Family

ID=80306850

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111490919.XA Active CN114095062B (zh) 2021-12-08 2021-12-08 毫米波大规模相控阵系统的宽带修正模拟波束成形方法

Country Status (1)

Country Link
CN (1) CN114095062B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016062190A1 (en) * 2014-10-23 2016-04-28 Huawei Technologies Co., Ltd. System and method for beam alignment
CN107493125A (zh) * 2017-09-21 2017-12-19 北京航空航天大学 一种任意波束宽度的模拟波束赋形方法
WO2018120339A1 (zh) * 2016-12-28 2018-07-05 西安交通大学 一种用于实际宽带大规模mimo系统的混合预编码设计方法
CN111478724A (zh) * 2020-04-15 2020-07-31 南京航空航天大学 一种面向无人机毫米波平台的三维波束搜索方法
CN112636773A (zh) * 2020-12-18 2021-04-09 电子科技大学 基于数字频域补偿的宽带时域波束成形方法
CN113258293A (zh) * 2021-05-12 2021-08-13 哈尔滨工业大学(深圳) 基于组合分辨率移相器的毫米波大规模阵列波束成形系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110365388B (zh) * 2019-07-31 2022-07-29 东南大学 一种低复杂度毫米波多播波束成形方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016062190A1 (en) * 2014-10-23 2016-04-28 Huawei Technologies Co., Ltd. System and method for beam alignment
WO2018120339A1 (zh) * 2016-12-28 2018-07-05 西安交通大学 一种用于实际宽带大规模mimo系统的混合预编码设计方法
CN107493125A (zh) * 2017-09-21 2017-12-19 北京航空航天大学 一种任意波束宽度的模拟波束赋形方法
CN111478724A (zh) * 2020-04-15 2020-07-31 南京航空航天大学 一种面向无人机毫米波平台的三维波束搜索方法
CN112636773A (zh) * 2020-12-18 2021-04-09 电子科技大学 基于数字频域补偿的宽带时域波束成形方法
CN113258293A (zh) * 2021-05-12 2021-08-13 哈尔滨工业大学(深圳) 基于组合分辨率移相器的毫米波大规模阵列波束成形系统

Also Published As

Publication number Publication date
CN114095062A (zh) 2022-02-25

Similar Documents

Publication Publication Date Title
Liu et al. Hybrid beamforming with sub-arrayed MIMO radar: Enabling joint sensing and communication at mmWave band
US8508410B2 (en) Adaptive antenna beamforming
KR101846832B1 (ko) 빔형성 기반 통신 방법 및 장치
CN108365873B (zh) 采用低精度adc毫米波大规模mimo自适应传输方法
Palacios et al. Speeding up mmWave beam training through low-complexity hybrid transceivers
WO2019128924A1 (zh) 一种天线装置及波束状态切换方法
Elbir et al. Terahertz-band integrated sensing and communications: Challenges and opportunities
CN112636773B (zh) 基于数字频域补偿的宽带时域波束成形方法
Alluhaibi et al. Capacity maximisation for hybrid digital-to-analog beamforming mm-wave systems
Nguyen et al. Beam squint effects in THz communications with UPA and ULA: Comparison and hybrid beamforming design
US10505596B2 (en) Wireless communication system
CN114095062B (zh) 毫米波大规模相控阵系统的宽带修正模拟波束成形方法
Blandino et al. Multi-user frequency-selective hybrid MIMO demonstrated using 60 GHz RF modules
Hong et al. Limitations of phased arrays for 5G wireless communications
CN113141203B (zh) 一种太赫兹通信中基于循环延迟的宽带混合预编码方法
Lou et al. Multi-stream signals separation based on space-time-isomeric (SPATIO) array using metasurface antennas
Gao et al. Wideband hybrid analog-digital beamforming massive MIMO systems based on Rotman lens
KR101748814B1 (ko) 하이브리드 빔포밍을 위한 송신기, 수신기 및 신호 송수신 방법
WO2002025776A1 (en) Shaped antenna beam
Torkzaban et al. Codebook design for composite beamforming in next-generation mmwave systems
Baghel et al. MIMO Antennas: A 5G Communication Perspective
CN114039638A (zh) 一种混合波束成形器与模数转换器联合设计方法
Bing A simplified genetic algorithm in multi-beam combination for mmWave communication system
Xiao et al. Multipath grouping for millimeter-wave communications
Bharadwaj et al. MIMO versus phased array antenna systems for 5G mobile communication systems

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