CN105223592B - 一种高动态环境下卫星导航接收机载波环跟踪装置 - Google Patents
一种高动态环境下卫星导航接收机载波环跟踪装置 Download PDFInfo
- Publication number
- CN105223592B CN105223592B CN201510634326.4A CN201510634326A CN105223592B CN 105223592 B CN105223592 B CN 105223592B CN 201510634326 A CN201510634326 A CN 201510634326A CN 105223592 B CN105223592 B CN 105223592B
- Authority
- CN
- China
- Prior art keywords
- carrier wave
- output
- doppler frequency
- reception signal
- integrator
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/29—Acquisition or tracking or demodulation of signals transmitted by the system carrier including Doppler, related
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种高动态环境下卫星导航接收机载波环跟踪装置,包括鉴别器:用于生成剥离伪码后的接收信号与本地复现载波信号的相位差或频率差,并将所述相位或频率差送入环路滤波器;环路滤波器:根据所述鉴别器的输出生成对接收信号速度和加速度的估计值;数控振荡器:根据所述环路滤波器的输出得到本地复现载波信号的频率,并生成本地复现载波信号输入鉴别器;开环加速度估计器:用于估计接收信号的多普勒频率变化率,并将所述多普勒频率变化率赋值给环路滤波器。本发明环路跟踪高动态信号的能力强,能有效提高环路性能,且当接收信号存在几十g、甚至上百g的加速度,或接收信号动态变化很大时,也能跟踪接收信号。
Description
技术领域
本发明涉及一种高动态环境下卫星导航接收机载波环跟踪装置。
背景技术
目前的高动态接收机方案中,主要通过优化环路参数和改变鉴别器结构两种途径实现高动态信号的跟踪。优化环路参数的方法一般通过使用较大的环路带宽提高环路的动态性能;改变鉴别器结构的方法一般通过增大鉴别器的牵引范围、采用鉴频器代替鉴相器等方法提高环路动态性能。上述两类方法可以在一定程度上提高环路跟踪高动态信号的能力,但当接收信号存在几十g、甚至上百g的加速度,或接收信号动态变化很大(如,短时间内加速度突然增大或减小)时,很难跟踪接收信号。
发明内容
本发明所要解决的技术问题是,针对现有技术不足,提供一种高动态环境下卫星导航接收机载波环跟踪装置。
为解决上述技术问题,本发明所采用的技术方案是:一种高动态环境下卫星导航接收机载波环跟踪装置,包括载波跟踪环和开环加速度估计器;
所述载波环包括:
鉴别器:用于生成剥离伪码后的接收信号与本地复现载波信号的相位差或频率差,并将所述相位或频率差送入环路滤波器;
环路滤波器:根据所述鉴别器的输出生成对接收信号速度和加速度的估计值;
数控振荡器:根据所述环路滤波器的输出得到本地复现载波信号的频率,并生成本地复现载波信号输入鉴别器;
开环加速度估计器:用于估计接收信号延时T1后的多普勒频率和接收信号的多普勒频率变化率,并将所述多普勒频率和多普勒频率变化率赋值给载波环的环路滤波器。
所述环路滤波器包括加速度估计积分器、速度估计积分器和比例支路;所述鉴别器的输出乘以系数K3,作为加速度估计积分器的输入;所述鉴别器的输出乘以系数K1后作为比例支路的输入;鉴别器输出值乘以系数K2,然后与加速度估计积分器的输出相加,得到的和值作为速度估计积分器的输入;速度估计积分器的输出和比例支路的输出相加后的和值,作为数控振荡器的输入;数控振荡器的输出送入所述鉴别器。
所述开环加速度估计器包括上侧支路和下侧支路;所述上侧支路和下侧支路结构相同;所述上侧支路计算当前接收信号的多普勒频率f2,下侧支路计算t1-T1时刻接收信号的多普勒频率f1。
所述上侧支路具体计算流程包括:
1)时长Ts的积分清零器对接收信号进行长度为Ts的积分操作,完成一次积分后积分器清零并开始下一次积分运算,每完成一次积分输出一个积分值;
2)时长Ts的积分清零器的输出作为N点存储器的输入,N点存储器存储完N点积分值后,将所述N点积分值并行输入L点FFT运算模块,L的值大于或等于N,且L的值为2的幂次;
3)L点FFT运算模块对N点积分值进行傅里叶变换,当N小于L时,在N点积分值后面补L-N个零,以使得L点FFT运算模块的输入为L点;
4)完成傅里叶变换后,对傅里叶变换结果取包络并搜索包络的最大值及该最大值的位置;
5)根据包络最大值的位置计算多普勒频率f2。
所述下侧支路具体计算流程包括:
1)将接收信号延时T1,时长Ts的积分清零器对延时T1后的接收信号进行长度为Ts的积分操作,完成一次积分后积分器清零并开始下一次积分运算,每完成一次积分输出一个积分值;
2)时长Ts的积分清零器的输出作为N点存储器的输入,N点存储器存储完N点积分值后,将所述N点积分值并行输入L点FFT运算模块,L的值大于或等于N,且L的值为2的幂次;
3)L点FFT运算模块对N点积分值进行傅里叶变换,当N小于L时,在N点积分值后面补L-N个零,以使得L点FFT运算模块的输入为L点;
4)完成傅里叶变换后,对傅里叶变换结果取包络并搜索包络的最大值及该最大值的位置;
5)根据包络最大值的位置计算多普勒频率f1。
接收信号的多普勒频率变化率a计算公式为:a=(f2-f1)/T1。
所述多普勒频率变化率a赋值给所述加速度估计积分器,所述多普勒频率f2赋值给所述速度估计积分器。
在引导阶段,所述载波环不工作,开环加速度估计器独立工作;在跟踪阶段,所述载波环和开环加速度估计器分别独立工作,每经过时间间隔T2后,分别提取开环加速度估计器得到的多普勒频率变化率a和载波环加速度估计积分器估计的频率变化率a1,并对两者做差,若差值大于门限Th,则用多普勒频率变化率a和多普勒频率f2分别初始化载波环的加速度估计积分器和速度估计积分器;若差值小于等于门限,则载波环继续独立运行。
与现有技术相比,本发明所具有的有益效果为:本发明捕获和引导阶段,通过FFT估计接收信号加速度,实现高动态下的GNSS信号捕获与引导;在跟踪阶段,通过FFT开环估计接收信号加速度,当载波环估计的加速度与FFT开环估计的加速度之差超过门限时,将FFT开环估计的加速度累加至环路的加速度估计器,本发明环路跟踪高动态信号的能力强,能有效提高环路性能,且当接收信号存在几十g、甚至上百g的的加速度,或接收信号动态变化很大时,也能跟踪接收信号。
附图说明
图1为本发明一实施例引导阶段环路原理图;
图2为本发明一实施例开环加速度估计器结构图;
图3为本发明一实施例跟踪阶段载波环原理图。
具体实施方式
本发明的目的是提供一种能够在复杂动态条件下捕获、跟踪卫星导航信号的基带处理装置,主要用于完成导航信号捕获之后的引导和跟踪阶段。捕获成功后,环路首先进入引导阶段。该阶段主要进行多普勒频率精捕获和环路的初始化(包括环路速度积分器的初始化和加速度积分器的初始化)。频率的精捕获通过对剥离伪码和粗捕获频率后的接收信号做FFT实现。为了得到频率变化率的估计(加速度),需要对相隔一段时间的两段信号分别进行频率精捕获,通过两次精捕获的频率差估计接收信号的频率变化率。得出的频率变化率对载波环的加速度估计积分器进行初始化,第二次频率精捕获结果对载波环的速度估计积分器进行初始化。
本发明提供的引导阶段环路结构如图1所示,包括两部分,即载波跟踪环和开环加速度估计器。引导阶段载波环不工作,只启动开环加速度估计器。图1中的载波跟踪环包括鉴别器、环路滤波器和数控振荡器(NCO)三部分。鉴别器生成接收信号(此处及以后所指的接收信号,均为剥离伪码后的信号)与本地复现载波信号的相位或频率差;鉴别器的输出送入环路滤波器,环路滤波器根据鉴别器的输出结果生成对接收信号速度和加速度的估计;NCO根据环路滤波器的输出得到本地复现载波信号的频率,并生成本地复现载波信号输入鉴别器。环路滤波器包括加速度估计积分器、速度估计积分器和比例支路三部分,环路滤波器中的系数T为载波环的预检测积分时间,K1、K2、K3为环路滤波器系数,上述系数根据环路带宽及其他环路相关指标确定(如环路带宽Bn为30Hz时,K1=2.4×Bn/0.7845=91.78,K2=1.1×Bn×Bn/0.61544=1608.6,K3=Bn×Bn×Bn/0.482813=55922.3)。鉴别器的输出乘以系数K3和K1后,分别作为加速度估计积分器和比例支路的输入;鉴别器输出值乘以K2和加速度估计积分器的输出相加,得到的和值作为速度估计积分器的输入;速度估计积分器的输出和比例支路的输出相加后的和值,为数控振荡器的输入;数控振荡器的输出送入鉴别器。
开环加速度估计器的结构如图2所示,包含两个支路,上侧支路计算当前信号的多普勒频率f2(为了便于描述,当前时刻记为t1),下侧支路计算t1-T1时刻信号的多普勒频率,两个支路的结构相同。下面只给出上侧支路详细运算流程如下:
1)时长Ts的积分清零器对接收信号进行长度为Ts的积分操作,完成一次积分后积分器清零并开始下一次积分运算,每完成一次积分输出一个积分值;
2)时长Ts的积分清零器的输出作为N点存储器的输入,N点存储器存储完N点积分值后并行输入L点FFT运算模块,L的值大于等于N且为2的幂次;
3)L点FFT计算N点积分值的傅里叶变换,当N小于L时,在N点积分值后面补L-N个零,以使得FFT的输入为L点;
4)完成FFT运算后,由搜索包络最大值模块对FFT结果取包络并搜索包络的最大值及其位置;
5)根据包络最大值的位置计算多普勒频率f2。
相关积分的长度Ts根据接收机捕获指标确定(动态指标要求越高,Ts值越小),两次FFT使用的输入数据的起始时间间隔为T1(T1的值需要在加速度估计精度和加速度估计范围之间折中,T1越小加速度估计范围越大,加速度估计精度越差),每次FFT的输入点数为N(N的值需要在加速度估计精度和加速度估计范围之间折中,N越小加速度估计范围越大,加速度估计精度越差),输出点数为L,L大于等于N且为2的幂次,L大于N时需要在输入的N点数据后面补零,使输入点数为L。两次FFT估计出的两个多普勒频率分别记为f1和f2,多普勒频率变化率的计算表达式为(f2-f1)/T1。开环加速度估计器得到的多普勒频率变化率赋值给图1中载波环的加速度估计积分器,得到的多普勒频率f2赋值给图1中载波环的速度估计积分器。
引导阶段结束后,环路进入跟踪阶段,载波环在这个阶段开始工作,与传统环路相比,本发明设计的跟踪环路增加了一个开环加速度估计器,其原理如图3所示(载波环各个模块的功能与图1相同;图3与图1的不同之处在于:图1为引导阶段,主要任务是初始化载波跟踪环,因此图1中只启动开环加速度估计器,载波跟踪环不工作;图3为跟踪阶段,载波跟踪环和开环加速度估计器分别独立工作),载波环和开环加速度估计器分别独立工作,每经过时间间隔T2后,分别提取开环加速度估计器得到的频率变化率a和载波环加速度估计积分器估计的频率变化率a1,并对两者做差,若差值大于门限Th,则说明当前的接收信号动态超过了载波跟踪环的动态容忍能力,需要用开环加速度估计器输出的a和f2分别初始化载波环的加速度估计积分器和速度估计积分器;若差值小于等于门限,则说明当前的接收信号动态在载波跟踪环动态容忍能力范围内,载波跟踪环继续独立运行。门限Th的值根据环路能够容忍的加速度误差确定(如接收信号在短时间内加速度发生突变,导致接收信号的实际加速度与载波环估计的加速度相差2g(g为重力加速度)时,环路能够自行收敛,当加速度误差大于2g时,环路容易失锁,则门限Th可设置为2g加速度对应的多普勒频率变化率)。T2的值由接收机的动态指标确定,T2的值越小,载波环适应复杂信号动态的能力越强,但开环加速度估计器估计的精度也越低。
Claims (5)
1.一种高动态环境下卫星导航接收机载波环跟踪装置,其特征在于,包括载波跟踪环和开环加速度估计器;
所述载波环包括:
鉴别器:用于生成剥离伪码后的接收信号与本地复现载波信号的相位差或频率差,并将所述相位差或频率差送入环路滤波器;
环路滤波器:根据所述鉴别器的输出生成对接收信号速度和加速度的估计值;所述环路滤波器包括加速度估计积分器、速度估计积分器和比例支路;所述鉴别器的输出乘以系数K3,作为加速度估计积分器的输入;所述鉴别器的输出乘以系数K1后作为比例支路的输入;鉴别器输出值乘以系数K2,然后与加速度估计积分器的输出相加,得到的和值作为速度估计积分器的输入;速度估计积分器的输出和比例支路的输出相加后的和值,作为数控振荡器的输入;数控振荡器的输出送入所述鉴别器;
数控振荡器:根据所述环路滤波器的输出得到本地复现载波信号的频率,并生成本地复现载波信号输入鉴别器;
开环加速度估计器:用于估计接收信号延时T1后的多普勒频率和接收信号的多普勒频率变化率,并将当前接收信号的多普勒频率f2和多普勒频率变化率赋值给载波环的环路滤波器;所述开环加速度估计器包括上侧支路和下侧支路;所述上侧支路和下侧支路结构相同;所述上侧支路计算当前接收信号的多普勒频率f2,下侧支路计算t1-T1时刻接收信号的多普勒频率f1;所述上侧支路具体计算流程包括:
1)时长Ts的积分清零器对接收信号进行长度为Ts的积分操作,完成一次积分后积分器清零并开始下一次积分运算,每完成一次积分输出一个积分值;
2)时长Ts的积分清零器的输出作为N点存储器的输入,N点存储器存储完N点积分值后,将所述N点积分值并行输入L点FFT运算模块,L的值大于或等于N,且L的值为2的幂次;
3)L点FFT运算模块对N点积分值进行傅里叶变换,当N小于L时,在N点积分值后面补L-N个零,以使得L点FFT运算模块的输入为L点;
4)完成傅里叶变换后,对傅里叶变换结果取包络并搜索包络的最大值及该最大值的位置;
5)根据包络最大值的位置计算多普勒频率f2。
2.根据权利要求1所述的高动态环境下卫星导航接收机载波环跟踪装置,其特征在于,所述下侧支路具体计算流程包括:
1)将接收信号延时T1,时长Ts的积分清零器对延时T1后的接收信号进行长度为Ts的积分操作,完成一次积分后积分器清零并开始下一次积分运算,每完成一次积分输出一个积分值;
2)时长Ts的积分清零器的输出作为N点存储器的输入,N点存储器存储完N点积分值后,将所述N点积分值并行输入L点FFT运算模块,L的值大于或等于N,且L的值为2的幂次;
3)L点FFT运算模块对N点积分值进行傅里叶变换,当N小于L时,在N点积分值后面补L-N个零,以使得L点FFT运算模块的输入为L点;
4)完成傅里叶变换后,对傅里叶变换结果取包络并搜索包络的最大值及该最大值的位置;
5)根据包络最大值的位置计算多普勒频率f1。
3.根据权利要求2所述的高动态环境下卫星导航接收机载波环跟踪装置,其特征在于,接收信号的多普勒频率变化率a计算公式为:a=(f2-f1)/T1。
4.根据权利要求3所述的高动态环境下卫星导航接收机载波环跟踪装置,其特征在于,所述多普勒频率变化率a赋值给所述加速度估计积分器,所述多普勒频率f2赋值给所述速度估计积分器。
5.根据权利要求4所述的高动态环境下卫星导航接收机载波环跟踪装置,其特征在于,在引导阶段,所述载波环不工作,开环加速度估计器独立工作;在跟踪阶段,所述载波环和开环加速度估计器分别独立工作,每经过时间间隔T2后,分别提取开环加速度估计器得到的多普勒频率变化率a和载波环加速度估计积分器估计的频率变化率a1,并对两者做差,若差值大于门限Th,则用多普勒频率变化率a和多普勒频率f2分别初始化载波环的加速度估计积分器和速度估计积分器;若差值小于等于门限,则载波环继续独立运行。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510634326.4A CN105223592B (zh) | 2015-09-30 | 2015-09-30 | 一种高动态环境下卫星导航接收机载波环跟踪装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510634326.4A CN105223592B (zh) | 2015-09-30 | 2015-09-30 | 一种高动态环境下卫星导航接收机载波环跟踪装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105223592A CN105223592A (zh) | 2016-01-06 |
CN105223592B true CN105223592B (zh) | 2017-11-17 |
Family
ID=54992648
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510634326.4A Active CN105223592B (zh) | 2015-09-30 | 2015-09-30 | 一种高动态环境下卫星导航接收机载波环跟踪装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105223592B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105699992B (zh) * | 2016-01-13 | 2017-12-15 | 上海交通大学 | 高动态gnss载波信号跟踪方法及*** |
CN107340528A (zh) * | 2016-04-28 | 2017-11-10 | 大唐半导体设计有限公司 | 一种实现高动态导航跟踪的方法和装置 |
CN107329155B (zh) * | 2016-12-23 | 2020-08-14 | 湖南北云科技有限公司 | 一种基于悲观计数器的卫星导航环路参数设置装置及方法 |
CN110068842B (zh) * | 2019-05-06 | 2021-01-29 | 西安开阳微电子有限公司 | 一种卫星信号高精度捕获方法 |
CN110784279B (zh) * | 2019-10-29 | 2022-08-26 | 北京无线电计量测试研究所 | 一种基于卫星双向的远程时间复现方法 |
CN113965435B (zh) * | 2021-10-14 | 2024-04-16 | 中国电子科技集团公司第五十四研究所 | 一种适用于低速信号的高动态频率跟踪方法 |
CN113965436B (zh) * | 2021-10-20 | 2023-06-27 | 北京临近空间飞行器***工程研究所 | 一种多普勒频偏及多普勒频偏加速度联合估计方法 |
CN114252873B (zh) * | 2021-12-09 | 2022-09-30 | 珠海安自达科技有限公司 | 一种毫米波雷达目标加速度精确实时估计方法及其*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5363112A (en) * | 1989-07-05 | 1994-11-08 | The Boeing Company | Noise suppression processor for a carrier tracking loop |
CN102176034A (zh) * | 2011-03-09 | 2011-09-07 | 东南大学 | Gps接收机卫星失锁判决方法 |
CN103592662A (zh) * | 2013-11-29 | 2014-02-19 | 中国航天科工信息技术研究院 | 一种gps信号接收机的载波跟踪方法及环路 |
-
2015
- 2015-09-30 CN CN201510634326.4A patent/CN105223592B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5363112A (en) * | 1989-07-05 | 1994-11-08 | The Boeing Company | Noise suppression processor for a carrier tracking loop |
CN102176034A (zh) * | 2011-03-09 | 2011-09-07 | 东南大学 | Gps接收机卫星失锁判决方法 |
CN103592662A (zh) * | 2013-11-29 | 2014-02-19 | 中国航天科工信息技术研究院 | 一种gps信号接收机的载波跟踪方法及环路 |
Non-Patent Citations (3)
Title |
---|
GNSS软件接收机高动态载波跟踪环路关键技术研究;张晶泊;《中国博士学位论文全文数据库 信息科技辑》;20120815;I136-28 * |
High Dynamic Carrier Tracking Using Kalman Filter Aided Phase-Lock Loop;Weibin Li, etal;《Wireless Communications, Networking and Mobile Computing, 2007》;20071008;673-676 * |
一种高动态BOC扩频信号跟踪算法设计与实现;倪盼盼;《中国优秀硕士学位论文全文数据库 信息科技辑》;20140115;I136-1230 * |
Also Published As
Publication number | Publication date |
---|---|
CN105223592A (zh) | 2016-01-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105223592B (zh) | 一种高动态环境下卫星导航接收机载波环跟踪装置 | |
CN101082664B (zh) | 一种用于高动态卫星导航接收机中信号快捕的装置及其方法 | |
CN102621563B (zh) | 一种基于fpga的gps软件接收机信号跟踪方法及其*** | |
CN109633711B (zh) | 一种超大动态、高灵敏度的扩频测控基带接收方法及装置 | |
CN103885053B (zh) | 一种基于追踪滤波器的雷达数据动目标检测处理方法 | |
CN102116866B (zh) | 整周载波的gps p和/或y码信号的跟踪方法和装置 | |
CN101839987B (zh) | 一种自适应gps软件接收机的实现方法 | |
CN103439697B (zh) | 基于动态规划的目标检测方法 | |
CN103116038B (zh) | 一种利用卫星接收机载波跟踪环路测量加速度的方法 | |
CN106772455A (zh) | 一种基于惯性信息辅助与参数估计的gnss反欺骗环路跟踪方法 | |
CN105807264B (zh) | 雷达脉冲重复频率检测与初始脉冲到达时间的估计方法 | |
CN106443726A (zh) | 一种基于预滤波的gnss矢量跟踪环路及其实现方法 | |
CN105093242B (zh) | 一种卫星接收机失锁后快速捕获和定位方法 | |
CN105607088A (zh) | 一种卫星导航多频接收机信号快速引导跟踪装置 | |
CN104880720A (zh) | 一种弱信号卫星导航***的位同步方法 | |
CN106707271A (zh) | 一种基于数字锁相环的自适应角度跟踪方法 | |
CN103439718A (zh) | 一种高阶boc调制信号的无模糊跟踪单元 | |
CN103344968A (zh) | 一种ins辅助的北斗信号捕获方法 | |
CN110780320B (zh) | 一种软硬件一体化的卫星导航信号处理方法 | |
CN202583467U (zh) | 一种基于fpga的gps软件接收机信号跟踪*** | |
CN105391539B (zh) | Qpsk载波同步锁定检测器 | |
CN104931940B (zh) | 一种基于锁相跟踪环路的穿墙探测墙体抑制方法 | |
CN105954775A (zh) | 惯导加速度辅助接收机跟踪环路算法 | |
CN105572701B (zh) | 一种用于北斗二代接收机弱信号跟踪的码环鉴相方法 | |
CN109613568A (zh) | 一种卫星导航接收机鉴频器测量噪声估算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |