CN110907887A - 一种基于均匀圆阵的超分辨测向装置及处理系统 - Google Patents

一种基于均匀圆阵的超分辨测向装置及处理系统 Download PDF

Info

Publication number
CN110907887A
CN110907887A CN201910998574.5A CN201910998574A CN110907887A CN 110907887 A CN110907887 A CN 110907887A CN 201910998574 A CN201910998574 A CN 201910998574A CN 110907887 A CN110907887 A CN 110907887A
Authority
CN
China
Prior art keywords
super
array
uniform circular
circular array
finding
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
CN201910998574.5A
Other languages
English (en)
Other versions
CN110907887B (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.)
Chengdu Grand Duke Bo Chuan Information Technology Co Ltd
Original Assignee
Chengdu Grand Duke Bo Chuan Information 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 Chengdu Grand Duke Bo Chuan Information Technology Co Ltd filed Critical Chengdu Grand Duke Bo Chuan Information Technology Co Ltd
Priority to CN201910998574.5A priority Critical patent/CN110907887B/zh
Publication of CN110907887A publication Critical patent/CN110907887A/zh
Application granted granted Critical
Publication of CN110907887B publication Critical patent/CN110907887B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S3/00Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/14Systems for determining direction or deviation from predetermined direction

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radio Transmission System (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本发明公开了一种基于均匀圆阵的超分辨测向装置及处理系统,包括计算机,所述计算机分别连接有射频开关驱动器、双通道接收机和校正信号源,所述双通道接收机连接有开关矩阵,所述开关矩阵分别与所述校正信号源、所述射频开关驱动器和九元双通道均匀圆阵连接。本发明通过合适的加权和切换开关的顺序,使得双通道法可以利用两个通道接收到的数据恢复阵列的协方差矩阵,大大减低系统复杂度,减小系统成本与体积,提高系统的测向精度,测试结果更加稳健有效。

Description

一种基于均匀圆阵的超分辨测向装置及处理系统
技术领域
本发明涉及超分辨测向技术领域,尤其涉及一种基于均匀圆阵的超分辨测向装置及处理系统。
背景技术
市面上所有的多通道超分辨测向方法,均是一个通道对应一个接收机,然而这种方法使得接收机数量随着阵元的增多随之增加,造成硬件成本与体积也相应增加,同时通道的数量的增加,也会增加通道间幅相不一致、误差增加,使得高分辨、超分辨测向估计算法性能恶化,导致测向系统稳定性与测向系统精度一直达不到国外先进水平。
上述缺陷,值得改进。
发明内容
为了克服现有技术的不足,本发明提供一种基于均匀圆阵的超分辨测向装置及处理系统。
本发明技术方案如下:
一种基于均匀圆阵的超分辨测向装置,包括计算机,所述计算机分别连接有射频开关驱动器、双通道接收机和校正信号源,所述双通道接收机连接有开关矩阵,所述开关矩阵分别与所述校正信号源、所述射频开关驱动器和九元双通道均匀圆阵连接。
在一些实施例中,所述计算机为工业计算机。
在一些实施例中,所述工业计算机的工作频段为20MHz~6000MHz。
本发明还提出一种基于均匀圆阵的超分辨测向处理系统,包括:
步骤S1:通过接收到的校正数据,计算出从辅助校正源之后到接收机部分的幅相误差值;
步骤S2:用被测信号数据,重构出阵列的协方差矩阵R0;
步骤S3:用矢量奇异值算法CESVD对R0进行解相干处理,得到解相干后的协方差矩阵R;
步骤S4:将协方差矩阵R进行特征分解,得到特征值和对应的特征矢量;
步骤S5:由特征值通过修正的信息论准则估计出信源个数N;
步骤S6:根据估计出的信源个数N,以及小特征值对应的噪声子空间Un,利用MUSIC算法得到谱估计;
步骤S7:通过谱峰搜索,得到N个信源的DOA。
在一些实施例中,所述信源的DOA通过软件界面展现给用户。
根据上述方案的本发明,其有益效果在于,通过合适的加权和切换开关的顺序,使得双通道法可以利用两个通道接收到的数据恢复阵列的协方差矩阵,大大减低系统复杂度,减小系统成本与体积,提高系统的测向精度,测试结果更加稳健有效。
附图说明
图1是本发明装置的结构示意图。
图2是本发明侧向抗干扰算法的流程图。
图3是本发明处理系统的流程图。
在图中,10、计算机;20、射频开关驱动器;30、双通道接收机;40、开关矩阵;50、九元双通道均匀圆阵;60、校正信号源。
具体实施方式
下面结合附图以及实施方式对本发明进行进一步的描述:
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要说明的是,术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
如图1所示,一种基于均匀圆阵的超分辨测向装置,包括计算机10,所计算机10分别连接有射频开关驱动器20、双通道接收机30和校正信号源60,双通道接收机30连接有开关矩阵40,开关矩阵40分别与校正信号源60、射频开关驱动器20和九元双通道均匀圆阵50连接。
在一些实施例中,计算机10为工业计算机10。
在一些实施例中,工业计算机10的工作频段为20MHz~6000MHz。
圆阵:相对于线阵,圆阵的估计性能是360°全方位无阵列模糊的,因此本实验系统采用九元的均匀圆阵,九元圆阵与小于九元的圆阵相比,在信号方向上的分辨率更优,为提高分辨率,采用双圈均匀圆阵,在结构上,这两圈均匀圆阵做成一体装在天线罩里。
校正信号源60为了解决通道不一致的问题,在矩阵开关前灌入一个校准信号源,类似于计算系统的传递函数,通过校准信号便可测量出校准信号到接收通道之间的幅相误差,也就是说,除了天线误差无法校正之外,其他的射频部分均得到校正,保证了开关矩阵40的加入不会给系统造成幅相不一致。
开关矩阵40,校正信号与被测信号必须分时接到接收机上,也就是说,校正信号源60与天线不能同时接入,这就需要开关矩阵40控制接入的顺序,且在接入天线阵列以后,需要开关矩阵40控制双通道接收机30分时接到不同的天线阵元上,以采集完整的数据恢复出协方差矩阵。
射频开关驱动器20,由计算机10给出的指令通过射频开关驱动器20来控制矩阵开关的工作状态,当矩阵开关开始工作时,测向的工作就开始了,先使校正源的信号接入到接收机,得到校正数据,再使来自天线阵列的信号加到接收机。
双通道接收机30,本系统采用双通道接收机30将分时得到校正数据与不同的两个天线阵元接收的被测信号数据提供给计算机10进行处理。
计算机10实现人机对话,控制测向的开始与结束,并且将采集到的数据通过软件进行处理,计算出信源的DOA后将结果通过软件界面直观的展现给用户,且一旦启动测向,当被测信源位置移动时,计算机10会实时地把信源方位展现给用户。
计算机10实现对整个系统控制和调度。
工作原理,首先通过射频开关驱动器20控制九元双通道均匀圆阵50开关,捕获环境中特定频段的无线电波信号,射频开关驱动器20发送握手协议信号与开关矩阵40取得通信,对需要打通的天线进行编码和驻留时间参数设置,然后通过双通道接收机30,把截获到的无线电信号转换数据基带信号,对每组天线采集测向数据进行测向抗干扰算法分析处理。其次,判断编码的天线是否都打通完成,若打通完成,表明轮询测向此轮结束,否则,继续打通未打通天线,最后控制校正信号源60开关,测量出接收通道间误差,把采集的37组测向数据送入到计算机10,完成超分辨阵列测向算法,输出测向示向度。
图2为测向抗干扰算法工作流程。数字中频信号通过数字下变频变为数据基带IQ后,首先发送握手协议信号与天线控制板取得通信,然后对需要打通的天线进行编码和驻留时间参数进行设置,其次对每组天线采集测向数据进行测向抗干扰处理,最后判断编码的天线是否都打通完成。测向抗干扰算法的具体步骤如下:
步骤一,发送握手协议信号与天线控制板取得通信,若没有成功,启动延时重发机制,若连续10次重发都没成功,向上位机发送错误类型码,告知人员检查电路连线和系统板子是否工作正常。
步骤二,对需要打通的天线进行编码(0~N)和驻留时间参数设置。
步骤三,启动天线控制触发信号。
步骤四,查询当前组天线i是否打通成功,若没有打通,重发触发信号,若连续重发5次触发信号天线无法打通,向上位机发送一个错误标志码,告知操作人员检修此组天线高低端子是否有损坏。
步骤五,若此组天线i打通成功,进行延时去抖,这样保证天线处于正常打通状态,采集测向数据没有抖动。
步骤六,延时去抖后,按照上位机指令采集一段测向数据。
步骤七,对采集到此组测向数据进行排序,剔除异常测向数据值,这样就大大提高测向准确度,起到测向抗干扰目的。
步骤八,抗干扰处理完成后,输出天线i的测向数据到存储器中。
步骤九,按照驻留时间间隔判断所有天线组是否打通完成,若没有,启动下一组天线打通流程。若完成打通所有天线组,结束轮询测向,等待下次测向命令。
图3所示,本发明还提出一种基于均匀圆阵的超分辨测向处理系统,包括:
步骤S1:通过接收到的校正数据,计算出从辅助校正源之后到接收机部分的幅相误差值,校正幅相的不一致。
步骤S2:用被测信号数据,重构出阵列的协方差矩阵R0;
步骤S3:用矢量奇异值算法CESVD对R0进行解相干处理,得到解相干后的协方差矩阵R;
步骤S4:将协方差矩阵R进行特征分解,得到特征值和对应的特征矢量;
步骤S5:由特征值通过修正的信息论准则估计出信源个数N;
步骤S6:根据估计出的信源个数N,以及小特征值对应的噪声子空间Un,利用MUSIC算法得到谱估计;
步骤S7:通过谱峰搜索,得到N个信源的DOA。
在一些实施例中,所述信源的DOA通过软件界面展现给用户。
需要说明的是,得到被测信源的DOA,是这里的校正源只能校正与方位无关的通道幅相误差,并且这种误差是时变的,需要经常校正,本系统是每进行一次测量时均校正一次幅相误差,这样保证超分辨测向具有优良的性能。
相对于现有技术的有益效果是,采用上述方案,本发明通过合适的加权和切换开关的顺序,使得双通道法可以利用两个通道接收到的数据恢复阵列的协方差矩阵,大大减低系统复杂度,减小系统成本与体积,提高系统的测向精度,使得测试结果更加稳健有效。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
上面对本发明专利进行了示例性的描述,显然本发明专利的实现并不受上述方式的限制,只要采用了本发明专利的方法构思和技术方案进行的各种改进,或未经改进将本发明专利的构思和技术方案直接应用于其它场合的,均在本发明的保护范围内。

Claims (5)

1.一种基于均匀圆阵的超分辨测向装置,其特征在于,包括计算机,所述计算机分别连接有射频开关驱动器、双通道接收机和校正信号源,所述双通道接收机连接有开关矩阵,所述开关矩阵分别与所述校正信号源、所述射频开关驱动器和九元双通道均匀圆阵连接。
2.根据权利要求1所述的一种基于均匀圆阵的超分辨测向装置,其特征在于,所述计算机为工业计算机。
3.根据权利要求2所述的一种基于均匀圆阵的超分辨测向装置,其特征在于,所述工业计算机的工作频段为20MHz~6000MHz。
4.一种基于均匀圆阵的超分辨测向处理系统,其特征在于,包括:
步骤S1:通过接收到的校正数据,计算出从辅助校正源之后到接收机部分的幅相误差值;
步骤S2:用被测信号数据,重构出阵列的协方差矩阵R0;
步骤S3:用矢量奇异值算法CESVD对R0进行解相干处理,得到解相干后的协方差矩阵R;
步骤S4:将协方差矩阵R进行特征分解,得到特征值和对应的特征矢量;
步骤S5:由特征值通过修正的信息论准则估计出信源个数N;
步骤S6:根据估计出的信源个数N,以及小特征值对应的噪声子空间Un,利用MUSIC算法得到谱估计;
步骤S7:通过谱峰搜索,得到N个信源的DOA。
5.根据权利要求4所述的一种基于均匀圆阵的超分辨测向处理系统,其特征在于,所述信源的DOA通过软件界面展现给用户。
CN201910998574.5A 2019-10-21 2019-10-21 一种基于均匀圆阵的超分辨测向处理系统 Active CN110907887B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910998574.5A CN110907887B (zh) 2019-10-21 2019-10-21 一种基于均匀圆阵的超分辨测向处理系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910998574.5A CN110907887B (zh) 2019-10-21 2019-10-21 一种基于均匀圆阵的超分辨测向处理系统

Publications (2)

Publication Number Publication Date
CN110907887A true CN110907887A (zh) 2020-03-24
CN110907887B CN110907887B (zh) 2021-02-23

Family

ID=69815939

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910998574.5A Active CN110907887B (zh) 2019-10-21 2019-10-21 一种基于均匀圆阵的超分辨测向处理系统

Country Status (1)

Country Link
CN (1) CN110907887B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113253190A (zh) * 2021-04-22 2021-08-13 中国电子科技集团公司第二十九研究所 一种分布式系统的射频通道全频段幅度校正的方法
CN114189300A (zh) * 2021-11-26 2022-03-15 成都大公博创信息技术有限公司 一种短波频段信号分选方法
CN116593956A (zh) * 2023-07-19 2023-08-15 成都大公博创信息技术有限公司 多天线阵与主机一体化设计的全频段空间谱监测测向系统

Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101149429A (zh) * 2006-09-18 2008-03-26 中国人民解放军空军工程大学导弹学院 适用于均匀圆阵的阵列互耦校正与信源测向方法
CN201360050Y (zh) * 2008-12-17 2009-12-09 成都华日通讯技术有限公司 九天线元的双通道相关干涉仪测向天线
CN102445679A (zh) * 2011-09-22 2012-05-09 成都中安频谱科技有限公司 三信道空间谱估计测向系统的测向方法
CN102608565A (zh) * 2012-03-23 2012-07-25 哈尔滨工程大学 一种基于均匀圆阵列的波达方向估计方法
US20140152504A1 (en) * 2012-12-02 2014-06-05 Khalifa University of Science, Technology & Research (KUSTAR) Method and system for measuring direction of arrival of wireless signal using circular array displacement
CN203774463U (zh) * 2013-12-31 2014-08-13 成都华日通讯技术有限公司 九天线元的空间谱估计测向天线
CN104049234A (zh) * 2014-03-18 2014-09-17 电子科技大学 采用均匀圆阵快速测定空间谱的方法
CN104506253A (zh) * 2015-01-13 2015-04-08 重庆大学 一种相控阵天线发射通道幅相误差校正系统及方法
CN104717026A (zh) * 2015-03-05 2015-06-17 成都大公博创信息技术有限公司 一种抗干扰单通道测向处理装置及方法
CN105099577A (zh) * 2015-06-11 2015-11-25 中国电子科技集团公司第三十六研究所 一种无线电监测设备及无线监测方法、装置
CN106932753A (zh) * 2017-03-14 2017-07-07 湖南雷神电子科技有限公司 反无人机被动探测和测向定位系统
CN109633578A (zh) * 2018-12-03 2019-04-16 上海无线电设备研究所 一种双通道高精度相位标校系统及方法
CN109633525A (zh) * 2018-11-29 2019-04-16 上海无线电设备研究所 一种快速的均匀圆阵空间谱测向实现方法
CN109660305A (zh) * 2018-11-16 2019-04-19 上海无线电设备研究所 一种对多通道接收机幅度、相位校正及补偿方法
CN109814063A (zh) * 2019-02-27 2019-05-28 中国电子科技集团公司第三十六研究所 一种干涉仪测向方法和装置
CN109959892A (zh) * 2019-04-11 2019-07-02 中国电子科技集团公司第三十六研究所 一种均匀圆阵双通道干涉仪测向方法、装置和系统
CN110082709A (zh) * 2019-05-14 2019-08-02 中国电子科技集团公司第三十六研究所 一种均匀五元圆阵测向方法

Patent Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101149429A (zh) * 2006-09-18 2008-03-26 中国人民解放军空军工程大学导弹学院 适用于均匀圆阵的阵列互耦校正与信源测向方法
CN201360050Y (zh) * 2008-12-17 2009-12-09 成都华日通讯技术有限公司 九天线元的双通道相关干涉仪测向天线
CN102445679A (zh) * 2011-09-22 2012-05-09 成都中安频谱科技有限公司 三信道空间谱估计测向系统的测向方法
CN102608565A (zh) * 2012-03-23 2012-07-25 哈尔滨工程大学 一种基于均匀圆阵列的波达方向估计方法
US20140152504A1 (en) * 2012-12-02 2014-06-05 Khalifa University of Science, Technology & Research (KUSTAR) Method and system for measuring direction of arrival of wireless signal using circular array displacement
CN203774463U (zh) * 2013-12-31 2014-08-13 成都华日通讯技术有限公司 九天线元的空间谱估计测向天线
CN104049234A (zh) * 2014-03-18 2014-09-17 电子科技大学 采用均匀圆阵快速测定空间谱的方法
CN104506253A (zh) * 2015-01-13 2015-04-08 重庆大学 一种相控阵天线发射通道幅相误差校正系统及方法
CN104717026A (zh) * 2015-03-05 2015-06-17 成都大公博创信息技术有限公司 一种抗干扰单通道测向处理装置及方法
CN105099577A (zh) * 2015-06-11 2015-11-25 中国电子科技集团公司第三十六研究所 一种无线电监测设备及无线监测方法、装置
CN106932753A (zh) * 2017-03-14 2017-07-07 湖南雷神电子科技有限公司 反无人机被动探测和测向定位系统
CN109660305A (zh) * 2018-11-16 2019-04-19 上海无线电设备研究所 一种对多通道接收机幅度、相位校正及补偿方法
CN109633525A (zh) * 2018-11-29 2019-04-16 上海无线电设备研究所 一种快速的均匀圆阵空间谱测向实现方法
CN109633578A (zh) * 2018-12-03 2019-04-16 上海无线电设备研究所 一种双通道高精度相位标校系统及方法
CN109814063A (zh) * 2019-02-27 2019-05-28 中国电子科技集团公司第三十六研究所 一种干涉仪测向方法和装置
CN109959892A (zh) * 2019-04-11 2019-07-02 中国电子科技集团公司第三十六研究所 一种均匀圆阵双通道干涉仪测向方法、装置和系统
CN110082709A (zh) * 2019-05-14 2019-08-02 中国电子科技集团公司第三十六研究所 一种均匀五元圆阵测向方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113253190A (zh) * 2021-04-22 2021-08-13 中国电子科技集团公司第二十九研究所 一种分布式系统的射频通道全频段幅度校正的方法
CN113253190B (zh) * 2021-04-22 2022-07-08 中国电子科技集团公司第二十九研究所 一种分布式系统的射频通道全频段幅度校正的方法
CN114189300A (zh) * 2021-11-26 2022-03-15 成都大公博创信息技术有限公司 一种短波频段信号分选方法
CN116593956A (zh) * 2023-07-19 2023-08-15 成都大公博创信息技术有限公司 多天线阵与主机一体化设计的全频段空间谱监测测向系统
CN116593956B (zh) * 2023-07-19 2023-09-29 成都大公博创信息技术有限公司 多天线阵与主机一体化设计的全频段空间谱监测测向系统

Also Published As

Publication number Publication date
CN110907887B (zh) 2021-02-23

Similar Documents

Publication Publication Date Title
CN110907887B (zh) 一种基于均匀圆阵的超分辨测向处理系统
CN102298685B (zh) 基于虚拟仪器的汽车电子射频识别参数检测系统
CN104717026B (zh) 一种抗干扰单通道测向处理装置及方法
KR102001394B1 (ko) 로그-영역 안테나 어레이 보간에 기반한 수신신호의 도래각 추정 방법과 이를 위한 장치
CN109683151A (zh) 非均匀噪声环境下基于矩阵补全的酉求根music角度估计方法
CN111352063A (zh) 一种均匀面阵中基于多项式求根的二维测向估计方法
CN110417490B (zh) 一种基于fpga的阵列通道校正方法
CN110208737B (zh) 一种超短波双通道宽带测向系统及门限判定测向方法
CN101436890A (zh) 一种天线阵列误差参数估计方法及其设备
CN111133691B (zh) 一种校准阵列天线的方法、装置和系统
CN113655435B (zh) 到达角度的确定方法、装置、信号接收设备、系统和介质
CN111698014B (zh) 一种天线阵列系统
CN116756697A (zh) 一种基于特征融合和maclnn的低慢小无人机识别方法
CN111190135A (zh) 一种适用于任意阵列的校准方法
CN111988740B (zh) 一种AoA估计方法、装置、设备及介质
CN113923087A (zh) 基于特定天线阵列切换顺序的载波频率偏移误差消除方法及系统
JP4679075B2 (ja) 電波到来角推定装置
CN111491368A (zh) 一种适用于aoa算法定位基站的校正方法及校正装置
CN113114315B (zh) 目标对象分组方法、装置、智能终端及存储介质
CN111175690A (zh) 联合对角化l型mimo雷达圆与非圆混合测向方法
CN113766424B (zh) 一种无人机室内定位方法、装置、计算机设备及存储介质
CN114785426B (zh) 多天线抗干扰方法、装置、设备和计算机可读存储介质
CN113382471B (zh) 一种基于离开角估计的定位系统和定位方法
CN110662160B (zh) 波束赋形的方法和装置及计算机可读存储介质
CN113514808B (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