CN111988108B - 蓝牙接收机中基于接入码的抗干扰高精度信号同步方法 - Google Patents

蓝牙接收机中基于接入码的抗干扰高精度信号同步方法 Download PDF

Info

Publication number
CN111988108B
CN111988108B CN202011138305.0A CN202011138305A CN111988108B CN 111988108 B CN111988108 B CN 111988108B CN 202011138305 A CN202011138305 A CN 202011138305A CN 111988108 B CN111988108 B CN 111988108B
Authority
CN
China
Prior art keywords
value
sampling
access code
correlation
peak value
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
CN202011138305.0A
Other languages
English (en)
Other versions
CN111988108A (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.)
Nanjing Qinheng Microelectronics Co ltd
Original Assignee
Nanjing Qinheng Microelectronics 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 Nanjing Qinheng Microelectronics Co ltd filed Critical Nanjing Qinheng Microelectronics Co ltd
Priority to CN202011138305.0A priority Critical patent/CN111988108B/zh
Publication of CN111988108A publication Critical patent/CN111988108A/zh
Application granted granted Critical
Publication of CN111988108B publication Critical patent/CN111988108B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04JMULTIPLEX COMMUNICATION
    • H04J11/00Orthogonal multiplex systems, e.g. using WALSH codes
    • H04J11/0023Interference mitigation or co-ordination
    • H04J11/0026Interference mitigation or co-ordination of multi-user interference
    • H04J11/0036Interference mitigation or co-ordination of multi-user interference at the receiver
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04JMULTIPLEX COMMUNICATION
    • H04J13/00Code division multiplex systems
    • H04J13/0077Multicode, e.g. multiple codes assigned to one user
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2656Frame synchronisation, e.g. packet synchronisation, time division duplex [TDD] switching point detection or subframe synchronisation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04JMULTIPLEX COMMUNICATION
    • H04J13/00Code division multiplex systems
    • H04J2013/0096Network synchronisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Synchronisation In Digital Transmission Systems (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明公开了一种蓝牙接收机中基于接入码的抗干扰高精度信号同步方法,将双方协定的接入码作为本地训练序列;对接收信号进行采样,不断对间隔一个码元周期的两个采样值进行延时差分运算,根据延时差分运算的结果得到待相关序列,将待相关序列与本地训练序列进行滑动相关运算,寻找相关值的峰值,峰值对应的采样时刻作为接入码的起始位置,实现帧定时同步。本发明利用蓝牙接入码进行信号同步,不容易受到同频段其他信号的干扰,信号检测准确度高,大大减少漏检或虚检问题,从而提高蓝牙接收机的工作效率,保证蓝牙通信准确安全,使通信***更加稳定可靠。

Description

蓝牙接收机中基于接入码的抗干扰高精度信号同步方法
技术领域
本发明属于蓝牙信号传输领域,尤其涉及一种蓝牙接收机中的基于接入码抗干扰高精度信号同步方法。
背景技术
在蓝牙接收机中,由于其帧结构的帧头部分具有很强的规律性(例如蓝牙4.2中帧头为01010101或10101010),因此通常利用帧头来进行信号检测及符号同步。但蓝牙通信所处频段为免授权频段,利用该频段进行通信的通信设备除了蓝牙还有wifi等众多遵循不同协议标准的设备。因此蓝牙通信过程中非常容易遭受同频干扰、邻道干扰等众多干扰的严重影响。此时,如果采用帧头的方式来进行信号检测,很容易出现漏检和虚检等问题。
发明内容
发明目的:为了解决现有技术中蓝牙接收机检测信号易出现漏检和虚检的问题,本发明提供一种蓝牙接收机中基于接入码的抗干扰高精度信号同步方法。
技术方案:一种蓝牙接收机中基于接入码的抗干扰高精度信号同步方法,包括以下步骤:
步骤一:将双方协定的接入码作为本地训练序列;
步骤二:对接收信号进行采样,不断对间隔一个码元周期的两个采样值进行延时差分运算,所述延时差分运算包括相位运算,将延时差分运算得到的相位差作为待相关序列,将待相关序列与本地训练序列进行滑动相关运算,设置第一门限值,在相关值超过第一门限值后出现的最大值作为相关值的峰值,峰值对应的采样时刻作为接入码的起始位置,实现帧定时同步。
进一步地,步骤二中,相关值的峰值的获取方法为:设置第一门限值和第一时间长度,在相关值超过第一门限值后的第一时间长度内出现的最大值,确定为相关值的峰值。
进一步地,步骤二中,相关值的峰值的获取方法为:设置第一门限值;当相关值超过第一门限值后首次出现极值时,确定为相关值的峰值。
进一步地,还包括步骤三:
基于接入码的位置与帧头的位置,根据接入码的起始位置确定帧头的位置,令接收信号r(t)在iTb/M时刻的抽样值为r i ,对帧头部分的采样值进行延时为2Tb的延时差分相关,并进行求和取角度运算
Figure 361111DEST_PATH_IMAGE001
其中v表示参与频偏估计的符号数,M表示过抽样因子,Tb为采样周期;
获得频偏的估计值
Figure 747093DEST_PATH_IMAGE002
为:
Figure 995672DEST_PATH_IMAGE003
一种蓝牙接收机中基于接入码的抗干扰高精度信号同步方法,包括以下步骤:
步骤一:将双方协定的接入码作为本地训练序列;
步骤二:对接收信号进行采样,不断对间隔一个码元周期的两个采样值进行延时差分运算,将延时差分运算的差分结果作为待相关序列,将待相关序列与本地训练序列进行滑动相关运算,寻找相关值的峰值,峰值对应的采样时刻作为接入码的起始位置,实现帧定时同步;
相关值的峰值的获取方法为:设置第二门限值和第二时间长度,在相关值超过第二门限值后的第二时间长度内出现的最大值,确定为相关值的峰值;
设第二门限值为Th1Ek,所述第二门限值的确定方法为:
设采样值为r(t),采样周期为Tb
Figure 859722DEST_PATH_IMAGE005
延时差分运算的差分结果
Figure 142936DEST_PATH_IMAGE006
为:
Figure 434240DEST_PATH_IMAGE007
其中,a0为增益,h为调制指数,I k 为传输的第k个码元,
Figure 904536DEST_PATH_IMAGE008
,g(t)为码元周期长度的矩形脉冲通过高斯滤波器的响应,n(t)为加性高斯白噪声,
Figure 572278DEST_PATH_IMAGE009
为解调频偏;
定义能量值Ek
Figure 208533DEST_PATH_IMAGE010
设第二门限的比例参数为Th1,根据第二门限的比例参数和能量值获得第二门限值Th1Ek的大小。
进一步地,步骤二中,相关值的峰值的获取方法为:设置第二门限值;当相关值超过第二门限值后首次出现极值时,确定为相关值的峰值。
进一步地,还包括步骤三:
基于接入码的位置与帧头的位置,根据接入码的起始位置确定帧头的位置,令接收信号r(t)在iTb/M时刻的抽样值为r i ,对帧头部分的采样值进行延时为2Tb的延时差分相关,并进行求和取角度运算
Figure 670738DEST_PATH_IMAGE001
其中v表示参与频偏估计的符号数,M表示过抽样因子,Tb为采样周期;
获得频偏的估计值
Figure 893909DEST_PATH_IMAGE011
为:
Figure 99763DEST_PATH_IMAGE003
有益效果:本发明提供一种蓝牙接收机中基于接入码的抗干扰高精度同步方法,相比较现有技术,利用蓝牙接入码进行信号同步,不容易受到同频段其他信号的干扰,信号检测准确度高,大大减少漏检或虚检问题,从而提高蓝牙接收机的工作效率,保证蓝牙通信准确安全,使通信***更加稳定可靠。
附图说明
图1为滑动相关运算的相关值的输出图;
图2为帧定时误差的概率分布图;
图3为频偏估计的均方根误差的输出图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步解释说明。
实施例一:
一种蓝牙接收机中基于接入码的抗干扰高精度信号同步方法,包括以下步骤:
步骤一:将双方协定的接入码作为本地训练序列;
步骤二:对接收信号进行采样,不断对间隔一个码元周期的两个采样值进行延时差分运算,延时差分运算包括相位运算,将相位运算的结果,即相位差作为待相关序列,将待相关序列与本地训练序列进行滑动相关运算,寻找相关值的峰值,峰值对应的采样时刻作为接入码的起始位置,实现帧定时同步。
其中,相关值的峰值获取方法为:设置第一门限值和第一时间长度,在相关值超过第一门限值后的第一时间长度内出现的最大值,确定为相关值的峰值。
本实施例,蓝牙BLE采用调制系数为0.5的GFSK作为具体调制方式,其调制形成的频带信号可表示如下:
Figure 91989DEST_PATH_IMAGE012
其中T b 为码元传输周期,I n 为传输的第n个码元,其取值为+1或-1。g(t)为码元周期长度的矩形脉冲通过高斯滤波器的响应,h为调制指数,其取值为0.5,f c为载波调制频率,而g(t)满足
Figure 990675DEST_PATH_IMAGE013
,令
Figure 701142DEST_PATH_IMAGE014
Figure 445107DEST_PATH_IMAGE015
。由于高斯滤波器的能量主要集中在
Figure 55955DEST_PATH_IMAGE016
范围内,可近似认为
Figure 859963DEST_PATH_IMAGE017
假设解调频偏为
Figure 57726DEST_PATH_IMAGE018
,信道为加性高斯白噪声信道,则接收下变频后的基带信号可以表示为:
Figure 870961DEST_PATH_IMAGE019
其中n(t)为加性高斯白噪声,
Figure 837780DEST_PATH_IMAGE020
为收发双方的频偏。
Figure 78269DEST_PATH_IMAGE021
其中
Figure 232170DEST_PATH_IMAGE022
对时间间隔为Tb两个接收信号样值进行差分运算,并进行相位运算,忽略噪声的影响,则
Figure 849096DEST_PATH_IMAGE023
令t0=Tb,则
Figure 168956DEST_PATH_IMAGE024
令发端接入码为cm(0≤m≤31),进行如下相关运算
Figure 580346DEST_PATH_IMAGE025
如果Ik+m=cm(0≤m≤31),则
Figure 487122DEST_PATH_IMAGE026
由于接入码的伪随机性,则
Figure 376581DEST_PATH_IMAGE027
,因此
Figure 317992DEST_PATH_IMAGE028
从上述可知,当接收信号和本地训练序列对齐时,相关值出现峰值,从而实现高精度的帧定时同步。图1显示了在过抽样因子为8的条件下采用本方法进行相关运算的结果,从图中可知,当接收信号与本地序列对齐时将出现明显的相关峰值,从而保证精确的帧定时同步。图2显示了SNR=9dB、过抽样因子为8时定时误差的概率分布。超过60%的概率定时误差为零,定时误差超过两个抽样间隔的概率接近零。
为找寻最佳帧定时位置,需要设置第一门限Th0,当yk>Th0时,则认为接近最佳帧定时位置,在随后的第一时间长度内的若干采样值中找寻使得yk最大的位置,则该位置为最佳帧定时位置。也可以采用不设定第一时间长度的方法,在yk>Th0后第一次达到极值点时认为达到了峰值,方法并不唯一,本实施例为了便于操作,选择前者找寻峰值。
实施例二:
实施例二与实施例一的区别在于在延时差分运算过程中不对相位进行运算,直接使用延时差分值执行后续的相关运算。
当t0=Tb时,忽略噪声的影响,延时差分值为
Figure 900283DEST_PATH_IMAGE007
但此时相关值将受到增益值a0的影响,因为
Figure 294355DEST_PATH_IMAGE029
如果Ik+m=cm(0≤m≤31),则
Figure 987505DEST_PATH_IMAGE030
由于G(0)远小于1/2,所以
Figure 744940DEST_PATH_IMAGE031
显然,此相关值受到了增益a0的影响。为此,需要计算如下能量值
Figure 498132DEST_PATH_IMAGE032
设置第二门限Th1Ek,其中,Th1为第二门限的比例参数,当
Figure 379501DEST_PATH_IMAGE033
时,则接近最佳帧定时位置。在随后的第二时间长度内的采样值中找寻使得
Figure 876341DEST_PATH_IMAGE034
最大的位置,则该位置为最佳帧定时位置。也可以采用不设定第一时间长度的方法,在
Figure 526765DEST_PATH_IMAGE035
后第一次达到极值点时认为达到了峰值,方法并不唯一,本实施例为了便于操作,选择前者找寻峰值。
相对实施例一,实施例二的方法不需要计算相位值,但需要计算能量值。实际中可根据具体需求,灵活选择。
实施例三:
实施例三在实施例一或实施例二的基础上增加了频偏估计的步骤,即步骤三:
在实现帧定时同步后,确定了接入码的起始位置,根据接入码的位置与帧头的位置之间的关系,也就确定了帧头的位置。由于帧头为01010101或10101010序列,因此帧头部分的波形是以周期为2Tb重复出现的。令接收信号r(t)在iTb/M时刻的抽样值为r i ,对帧头采样值进行延时为2Tb的延时差分相关,并进行求和取角度运算
Figure 185280DEST_PATH_IMAGE036
其中v表示参与频偏估计的符号数,M表示过抽样因子。
由此可以获得频偏的估计值
Figure 553944DEST_PATH_IMAGE037
图3显示了采用此估计方法后的均方根频偏估计误差(RMSE)。SNR=9dB时,此误差小于2.2kHz。由于符号速率为1M符号/s,因此归一化误差约为2.2×10-4。对于GFSK解调来说,此频偏估计精度完全满足需求。

Claims (7)

1.一种蓝牙接收机中基于接入码的抗干扰高精度信号同步方法,其特征在于,包括以下步骤:
步骤一:将双方协定的接入码作为本地训练序列;
步骤二:对接收信号进行采样,不断对间隔一个码元周期的两个采样值进行延时差分运算
Figure 469097DEST_PATH_IMAGE001
Figure 979713DEST_PATH_IMAGE002
Figure 228292DEST_PATH_IMAGE003
分别为两个采样值,并计算
Figure 485485DEST_PATH_IMAGE001
的相位得到相位差,将延时差分运算得到的相位差作为待相关序列,将待相关序列与本地训练序列进行滑动相关运算,设置第一门限值,在相关值超过第一门限值后出现的最大值作为相关值的峰值,峰值对应的采样时刻作为接入码的起始位置,实现帧定时同步。
2.根据权利要求1所述的蓝牙接收机中基于接入码的抗干扰高精度信号同步方法,其特征在于,步骤二中,相关值的峰值的获取方法为:设置第一门限值和第一时间长度,在相关值超过第一门限值后的第一时间长度内出现的最大值,确定为相关值的峰值。
3.根据权利要求1所述的蓝牙接收机中基于接入码的抗干扰高精度信号同步方法,其特征在于,步骤二中,相关值的峰值的获取方法为:设置第一门限值;当相关值超过第一门限值后首次出现极值时,确定为相关值的峰值。
4.根据权利要求1至3任一所述的蓝牙接收机中基于接入码的抗干扰高精度信号同步方法,其特征在于,还包括步骤三:
基于接入码的位置与帧头的位置,根据接入码的起始位置确定帧头的位置,令接收信号r(t)在iTb/M时刻的抽样值为r i ,对帧头部分的采样值进行延时为2Tb的延时差分相关,并进行求和取角度运算
Figure 96595DEST_PATH_IMAGE004
其中v表示参与频偏估计的符号数,M表示过抽样因子,Tb为采样周期;
获得频偏的估计值
Figure 653479DEST_PATH_IMAGE005
为:
Figure 717249DEST_PATH_IMAGE006
5.一种蓝牙接收机中基于接入码的抗干扰高精度信号同步方法,其特征在于,包括以下步骤:
步骤一:将双方协定的接入码作为本地训练序列;
步骤二:对接收信号进行采样,不断对间隔一个码元周期的两个采样值进行延时差分运算
Figure 775204DEST_PATH_IMAGE001
Figure 912924DEST_PATH_IMAGE002
Figure 499764DEST_PATH_IMAGE003
分别为两个采样值,将延时差分运算的差分结果作为待相关序列,将待相关序列与本地训练序列进行滑动相关运算,寻找相关值的峰值,峰值对应的采样时刻作为接入码的起始位置,实现帧定时同步;
相关值的峰值的获取方法为:设置第二门限值和第二时间长度,在相关值超过第二门限值后的第二时间长度内出现的最大值,确定为相关值的峰值;
设第二门限值为Th1Ek,所述第二门限值的确定方法为:
设采样值为r(t),采样周期为Tb
Figure 722935DEST_PATH_IMAGE008
延时差分运算的差分结果
Figure 53422DEST_PATH_IMAGE009
为:
Figure 45649DEST_PATH_IMAGE010
其中,a0为增益,h为调制指数,I k 为传输的第k个码元,
Figure 68968DEST_PATH_IMAGE011
,g(t)为码元周期长度的矩形脉冲通过高斯滤波器的响应,n(t)为加性高斯白噪声,
Figure 107331DEST_PATH_IMAGE012
为解调频偏;
定义能量值Ek
Figure 851296DEST_PATH_IMAGE013
设第二门限的比例参数为Th1,根据第二门限的比例参数和能量值获得第二门限值Th1Ek的大小。
6.根据权利要求5所述的蓝牙接收机中基于接入码的抗干扰高精度信号同步方法,其特征在于,步骤二中,相关值的峰值的获取方法为:设置第二门限值;当相关值超过第二门限值后首次出现极值时,确定为相关值的峰值。
7.根据权利要求5或6所述的蓝牙接收机中基于接入码的抗干扰高精度信号同步方法,其特征在于,还包括步骤三:
基于接入码的位置与帧头的位置,根据接入码的起始位置确定帧头的位置,令接收信号r(t)在iTb/M时刻的抽样值为r i ,对帧头部分的采样值进行延时为2Tb的延时差分相关,并进行求和取角度运算
Figure 91172DEST_PATH_IMAGE004
其中v表示参与频偏估计的符号数,M表示过抽样因子,Tb为采样周期;
获得频偏的估计值
Figure 895180DEST_PATH_IMAGE014
为:
Figure 217577DEST_PATH_IMAGE006
CN202011138305.0A 2020-10-22 2020-10-22 蓝牙接收机中基于接入码的抗干扰高精度信号同步方法 Active CN111988108B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011138305.0A CN111988108B (zh) 2020-10-22 2020-10-22 蓝牙接收机中基于接入码的抗干扰高精度信号同步方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011138305.0A CN111988108B (zh) 2020-10-22 2020-10-22 蓝牙接收机中基于接入码的抗干扰高精度信号同步方法

Publications (2)

Publication Number Publication Date
CN111988108A CN111988108A (zh) 2020-11-24
CN111988108B true CN111988108B (zh) 2020-12-29

Family

ID=73451010

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011138305.0A Active CN111988108B (zh) 2020-10-22 2020-10-22 蓝牙接收机中基于接入码的抗干扰高精度信号同步方法

Country Status (1)

Country Link
CN (1) CN111988108B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113497775B (zh) * 2021-09-07 2021-12-03 南京沁恒微电子股份有限公司 一种蓝牙接收机coded模式下的高灵敏度接收方法
CN115226197A (zh) * 2022-07-19 2022-10-21 山东闻远通信技术有限公司 无线通信的帧定时同步方法、装置、电子设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10257042A (ja) * 1997-03-13 1998-09-25 Nippon Telegr & Teleph Corp <Ntt> フレーム同期回路
CN101453238A (zh) * 2007-11-30 2009-06-10 安捷伦科技有限公司 Tdma通信***的差分中导码获取和频偏估计
CN110071888A (zh) * 2019-02-28 2019-07-30 成都坤恒顺维科技股份有限公司 一种高速数传中的快速时间同步方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08237173A (ja) * 1994-12-13 1996-09-13 Canon Inc スペクトラム拡散通信システム
CN103023529B (zh) * 2012-12-05 2014-08-20 天津光电通信技术有限公司 一种基于FPGA实现过采样Golay序列的同步检测的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10257042A (ja) * 1997-03-13 1998-09-25 Nippon Telegr & Teleph Corp <Ntt> フレーム同期回路
CN101453238A (zh) * 2007-11-30 2009-06-10 安捷伦科技有限公司 Tdma通信***的差分中导码获取和频偏估计
CN110071888A (zh) * 2019-02-28 2019-07-30 成都坤恒顺维科技股份有限公司 一种高速数传中的快速时间同步方法

Also Published As

Publication number Publication date
CN111988108A (zh) 2020-11-24

Similar Documents

Publication Publication Date Title
US9954701B2 (en) Bluetooth signal receiving method and device using improved packet detection and symbol timing acquisition
US5131006A (en) Carrier detection for a wireless local area network
CN111988108B (zh) 蓝牙接收机中基于接入码的抗干扰高精度信号同步方法
US20060269019A1 (en) Synchronization signal detector and method
US8254437B2 (en) Transmitting apparatus, receiving apparatus and communication system
US7457378B1 (en) Apparatus and method for RF packet detection and symbol timing recovery
US7245654B2 (en) Carrier sensing, signal quality and link quality in a receiver
CN109379314B (zh) 高速突发数字解调方法和设备
EP1436913B1 (en) A method and apparatus for the detection and classification of signals utilizing known repeated training sequences
US20040151269A1 (en) Receiver sampling in an ultra-wideband communications system
WO2000035159A1 (en) Delay spread estimation for multipath fading channels
CN110249542B (zh) 数字无线电通信
He et al. Adaptive synchronization for non-coherent UWB receivers
WO2015176597A1 (zh) 一种基于ebpsk的通信方法及***
US7532687B2 (en) Dual stage automatic gain control in an ultra wideband receiver
US7349493B2 (en) Demodulation with separate branches for phase and amplitude
CN115174322B (zh) 基于混沌成型滤波器的雷达通信一体化方法
US20050013390A1 (en) Hybrid UWB receiver with matched filters and pulse correlator
US6937672B1 (en) Determining the position of a constant frequency interval in a telecommunication signal
CN110535620B (zh) 一种基于判决反馈的信号检测与同步方法
Olonbayar et al. Synchronisation performance of wireless sensor networks
CN115361037B (zh) 一种慢衰落信道下的超宽带rake接收方法及装置
CN110324272B (zh) 一种基于组合环的快速载波恢复***及方法
CN111007550B (zh) 一种卫星导航理想信号相关功率估计方法
CN109547371B (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