CN107479071A - 伪卫星接收机和伪卫星信号接收方法 - Google Patents

伪卫星接收机和伪卫星信号接收方法 Download PDF

Info

Publication number
CN107479071A
CN107479071A CN201610395814.9A CN201610395814A CN107479071A CN 107479071 A CN107479071 A CN 107479071A CN 201610395814 A CN201610395814 A CN 201610395814A CN 107479071 A CN107479071 A CN 107479071A
Authority
CN
China
Prior art keywords
code
doppler
estimation
carrier
renewal
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
CN201610395814.9A
Other languages
English (en)
Other versions
CN107479071B (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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Priority to CN201610395814.9A priority Critical patent/CN107479071B/zh
Publication of CN107479071A publication Critical patent/CN107479071A/zh
Application granted granted Critical
Publication of CN107479071B publication Critical patent/CN107479071B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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/30—Acquisition or tracking or demodulation of signals transmitted by the system code related

Abstract

本申请公开了一种伪卫星接收机和伪卫星信号接收方法，在捕获处理后，通过码环迭代更新码多普勒估计，并基于码多普勒和载波多普勒之间的关系，更新载波多普勒估计，提高跟踪处理的初始载波多普勒估计的精度。这样，使得后续跟踪处理中载波环路能够收敛，以实现伪卫星信号的跟踪。

Description

伪卫星接收机和伪卫星信号接收方法
技术领域
本申请涉及伪卫星接收机和伪卫星信号接收方法。
背景技术
伪卫星系统作为一种可选的地面增强系统，可以应用于增强GNSS系统的可靠性和完好性，在恶劣的接收环境下，如城市、室内等环境里，能够作为对GNSS系统的补充，甚至由伪卫星系统独立提供定位导航功能。
为应对远近效应，伪卫星信号的扩频调制方式与实际卫星导航信号有所不同，常见的伪卫星信号在码分复用的基础上还采用了脉冲调制。图1示出了伪卫星信号脉冲调制的示意图。如图所示，脉冲调制的扩频码的播发周期为1ms，但有效的扩频码的播发时长只有100us，每隔1ms播发一次，利用不同的时隙和伪码序列区分伪卫星。
但这一调制方式对传统的信号接收处理方式提出了挑战，由于积分时间缩短，导致了自相关函数滚降对频率偏差不敏感，使得基带捕获算法得到的载波多普勒值具有较大的偏差，造成载波环路无法收敛而无法实现伪卫星信号的跟踪。
发明内容
本申请的目的是提供一种伪卫星接收机和伪卫星信号接收方法。
根据本申请的一个方面，公开了一种伪卫星接收机，包括捕获单元和跟踪单元，接收机的捕获单元进行捕获处理后，接收机的跟踪单元通过码环迭代更新码多普勒估计，并基于码多普勒和载波多普勒之间的关系，更新载波多普勒估计，提高跟踪处理的初始载波多普勒估计的精度。
根据本申请的另一个方面，公开了一种伪卫星信号接收方法，在捕获处理后，通过码环迭代更新码多普勒估计，并基于码多普勒和载波多普勒之间的关系，更新载波多普勒估计，提高跟踪处理的初始载波多普勒估计的精度。
附图说明
图1示出了伪卫星信号脉冲调制的示意图。
图2示出了根据本申请一种实施方式的伪卫星信号接收方法的示意图。
图3示出了根据本申请一种实施方式的伪卫星接收机的示意图。
图4示出了根据本申请一种实施方式的伪卫星接收机的预估计器的示意图。
图5示出了根据本申请另一种实施方式的伪卫星接收机的预估计器的环路滤波器的示意图。
具体实施方式
下面参照附图对本申请公开的伪卫星接收机和伪卫星信号接收方法进行详细说明。为简明起见，本申请各实施例的说明中，相同或类似的装置使用了相同或相似的附图标记。
根据本申请的实施方式，提供了一种伪卫星接收机和伪卫星信号接收方法，在捕获处理后，接收机通过码环迭代更新码多普勒估计，并基于码多普勒和载波多普勒之间的关系，更新载波多普勒估计，提高接收机跟踪处理前的初始载波多普勒估计的精度。
码多普勒和载波多普勒的对应关系可以从卫星的信号模型中得出，卫星信号的动态离散模型为：
其中：
分别表示采样时刻的码相位、载波相位和载波多普勒频率以及多普勒频率变化率；
表示载波多普勒和码多普勒之间的比例关系；
fcode表示伪码速率；fcarr表示载波频率。
在现有的卫星信号接收机中，通常跟踪处理中，载波环与码环同时工作，并且载波环的输出一直对码环提供多普勒辅助，这是因为载波环颤动比码环颤动的噪声要小几个量级，因而更准确。现有技术普遍认为，载波跟踪环从接收机动态应力门限的角度看是脆弱的环节，企图用码环的输出去辅助载波环将是灾难性的(参考《GPS原理与应用》,第五章第5.4节,电子工业出版社，2012年)。
然而，在伪卫星系统中，由于伪卫星信号采用脉冲调制信号，搜索使用的多普勒频格较宽，导致最终得到的载波多普勒估计的误差较大，无法直接应用于初始化跟踪处理中的载波环的多普勒估计，更无法对码环提供辅助。
根据本申请的伪卫星接收机和伪卫星信号接收方法，在捕获处理后，首先通过码环迭代更新码多普勒估计，并以码环迭代更新的码多普勒估计来辅助更新载波多普勒估计，提高跟踪处理的初始载波多普勒估计的精度，使得后续跟踪处理中载波环路能够收敛，以实现伪卫星信号的跟踪。
图2示出了根据本申请一种实施方式的伪卫星信号接收方法的示意图。如图所示，在步骤S1中进行捕获处理。在步骤S2中，进行码环估计，通过码环迭代更新码多普勒估计，并基于码多普勒和载波多普勒之间的关系，更新载波多普勒估计，提高载波多普勒估计的精度。在步骤S3中进行捕获处理。由于跟踪处理的初始载波多普勒估计的精度提高，使得后续载波环路能够收敛，从而实现伪卫星信号的跟踪。
图3示出了根据本申请一种实施方式的伪卫星接收机的示意图。如图所示，伪卫星接收机包括捕获单元100和跟踪单元200。
捕获单元100的捕获处理主要包括搜索各个可见伪卫星信号，并对捕获搜索到的可见伪卫星信号进行粗略码相位估计和粗略载波多普勒估计。
导航信号的捕获是一个二维搜索的过程，通过码相位和载波多普勒的搜索确定信号的初始码相位和载波多普勒偏移。
将接收到的伪卫星信号(通常可以建模为经过下变频的中频信号SIF(n))与本地复现载波进行混频，再进行本地复现码CL(n)相关运算，得到关于码偏τ，多普勒频偏fe的函数：
S(n)＝aR(τ)sinc(fenTs)，
R(τ)＝normalized[∑SIF(n)·CL(n-τ)]，
其中，n表示采样点，Ts表示采样周期，R(τ)表示归一化的码自相关函数。
对上述函数进行傅里叶变换，以并行算出每一个频格fk上的相关值大小：
捕获单元100通过混频、相关以及傅里叶变换运算后，可以并行得到所有待搜索频格上的相关值。经过比较可以得到最大相关幅值所对应的频格fk和其相关值。
针对伪卫星的脉冲调制信号，搜索使用的多普勒频格较宽，导致最终得到的多普勒频偏估计值的误差较大，无法直接应用于跟踪环节。
跟踪单元200包括预估计器210和跟踪模块220。预估计器210根据粗略的载波多普勒估计得到粗略的码多普勒估计，通过码环迭代更新码多普勒估计和码相位估计，并根据更新的码多普勒估计更新载波多普勒估计。跟踪模块220以更新的载波多普勒估计、更新的码多普勒估计和更新的码相位估计作为跟踪环路的初始值，进行伪卫星信号跟踪处理。
图4示出了根据本申请的一种实施方式的预估计器的示意图。如图所示，预估计器210包括转换器211、本地信号复现单元212、相关器213、码相位鉴别器214、环路滤波器215。
转换器211根据码多普勒和载波多普勒之间的关系，实现码多普勒估计和载波多普勒估计之间的转换。例如，转换器211可以将来自捕获单元100的粗略的载波多普勒估计转换为粗略的码多普勒估计，以作为预估计器210的输入值。当预估计器的码环迭代收敛后，转换器211可以将更新的码多普勒估计转换为更新的载波多普勒估计，输出至跟踪模块220。
本地信号复现单元212以粗略的载波多普勒估计、粗略的码多普勒估计和粗略的码相位估计作为初始值，并以环路滤波器迭代更新的码多普勒估计作为反馈值复现本地信号。在预估计器的初始化阶段，本地信号复现单元212以粗略的载波多普勒估计、粗略的码多普勒估计和粗略的码相位估计作为初始值来复现本地信号。在随后的预估计器工作过程中，本地信号复现单元212根据环路滤波器215所反馈的迭代更新的码多普勒估计，复现本地信号。
相关器213将复现的本地信号与接收到的伪卫星信号进行相关，获得相关值。例如，相关器213可以计算复现的本地信号的超前/滞后信号与接收到的伪卫星信号之间的相关值。
码相位鉴别器214根据所述相关值，鉴别复现的本地信号与接收到的卫星信号之间的码相位偏差。
环路滤波器215根据所述码相位偏差，迭代更新码多普勒估计。
图5示出了根据本申请的一种实施方式的预估计器中的环路滤波器的示意图。其中，环路滤波器的表达式为：
x0(k)＝x0(k-1)+p1Δd
x1(k)＝x0(k)+p2Δd，
其中，Δd表示码相位偏差；p1和p2表示放大器的放大倍数；节点x0表示迭代更新的码多普勒估计，节点x0的初始值为粗略的码多普勒估计，当码环迭代更新收敛后得到的节点x0的更新值，即为更新的码多普勒估计；节点x1表示码多普勒估计的反馈值，反馈至本地信号复现单元。
预估计器210的上述迭代过程可以至码环迭代更新收敛时结束。此时，预估计器210的环路滤波器215输出更新的码多普勒估计。其中，码环迭代更新的收敛通过码环迭代时间确定。例如，可以认为，1s后码环迭代已经收敛。此外，码环迭代更新的收敛也可以通过所述偏差的大小确定。例如，可以认为，当码相位偏差鉴别输出结果小于预设码相位偏差门限时，码环迭代已经收敛。
跟踪模块220以来自预估计器210的更新的载波多普勒估计、更新的码多普勒估计和更新的码相位估计作为的初始值，进行伪卫星信号跟踪处理。例如，跟踪模块220可以先执行牵引过程，再执行跟踪过程，以对载波相位的动态变化进行实时跟踪，提取伪距和载波相位，解调出导航电文等。
此外，根据本申请的一种实施方式的伪卫星信号接收方法还可以包括：通过捕获处理获得伪卫星信号的粗略的载波多普勒估计和粗略的码相位估计；根据粗略的载波多普勒估计得到粗略的码多普勒估计，通过码环迭代更新码多普勒估计和码相位估计，并根据更新的码多普勒估计更新载波多普勒估计；以更新的载波多普勒估计、更新的码多普勒估计和更新的码相位估计作为跟踪处理的初始值，进行伪卫星信号跟踪处理。
更进一步地，根据本申请的一种实施方式的伪卫星信号接收方法还可以进一步包括：根据码多普勒和载波多普勒之间的关系，将粗略的载波多普勒估计转换为粗略的码多普勒估计；以粗略的载波多普勒估计、粗略的码多普勒估计和粗略的码相位估计作为初始值，并以码环迭代更新的码多普勒估计作为反馈值复现本地信号；将复现的本地信号与接收到的伪卫星信号进行相关，获得相关值；根据相关值，鉴别复现的本地信号与接收到的卫星信号之间的码相位偏差；根据码相位偏差，迭代更新码多普勒估计，直到码环迭代更新收敛；将更新的码多普勒多普勒估计转换为更新的载波多普勒估计。
根据本申请的伪卫星接收机和伪卫星信号接收方法在跟踪处理前，首先进行码环估计，以码环迭代更新码多普勒估计，获得了精度更高的、更新的码多普勒估计，并进一步更新载波多普勒估计，提高了载波多普勒估计的精度。此外，根据本申请的预估计器中的码环结构与现有跟踪处理中的码环结构类似，因此仅需要在现有接收机的硬件和/或软件结构基础上略作调整即可，实现方式简单，节省成本。
根据本申请的实施方式可以硬件、软件或其组合的形式来实现。本申请的一个方面提供了包括用于实现根据本申请的实施方式的伪卫星接收机和伪卫星信号接收方法可执行指令的计算机程序。此外，此类计算机程序可使用例如光学或磁性可读介质、芯片、ROM、PROM或其它易失性或非易失性设备的任何形式的存储器来存储。根据本申请的一种实施例，提供了存储此类计算机程序的机器可读存储器。
以上参照附图对本申请的示例性的实施方案进行了描述。本领域技术人员应该理解，上述实施方案仅仅是为了说明的目的而所举的示例，而不是用来进行限制，凡在本申请的教导和权利要求保护范围下所作的任何修改、等同替换等，均应包含在本申请要求保护的范围内。

Claims (10)

1.一种伪卫星接收机，包括捕获单元和跟踪单元，其特征在于，接收机的捕获单元进行捕获处理后，接收机的跟踪单元通过码环迭代更新码多普勒估计，并基于码多普勒和载波多普勒之间的关系，更新载波多普勒估计，提高跟踪处理的初始载波多普勒估计的精度。
2.如权利要求1所述的接收机，其中，
所述捕获单元获得伪卫星信号的粗略的载波多普勒估计和粗略的码相位估计；
所述跟踪单元包括预估计器和跟踪模块，
所述预估计器根据所述粗略的载波多普勒估计得到粗略的码多普勒估计，通过码环迭代更新码多普勒估计和码相位估计，并根据更新的码多普勒估计更新载波多普勒估计；
所述跟踪模块以更新的载波多普勒估计、更新的码多普勒估计和更新的码相位估计作为跟踪模块的初始值，进行伪卫星信号跟踪处理。
3.如权利要求2所述的接收机，其中，所述预估计器进一步包括：
转换器，根据码多普勒和载波多普勒之间的关系，实现码多普勒估计和载波多普勒估计之间的转换；
本地信号复现单元，以所述粗略的载波多普勒估计、粗略的码多普勒估计和粗略码相位估计作为初始值，并以码环迭代更新的码多普勒估计作为反馈值复现本地信号；
相关器，将复现的本地信号与接收到的伪卫星信号进行相关，获得相关值；
码相位鉴别器，根据所述相关值，鉴别复现的本地信号与接收到的卫星信号之间的码相位偏差；
环路滤波器，根据所述码相位偏差，迭代更新码多普勒估计，直到码环迭代更新收敛。
4.如权利要求2所述的接收机，其中，码环迭代更新的收敛通过码环迭代时间确定。
5.如权利要求2所述的接收机，其中，码环迭代更新的收敛通过所述码相位偏差的大小确定。
6.一种伪卫星信号接收方法，其特征在于，在捕获处理后，通过码环迭代更新码多普勒估计，并基于码多普勒和载波多普勒之间的关系，更新载波多普勒估计，提高跟踪处理的初始载波多普勒估计的精度。
7.如权利要求6所述的伪卫星信号接收方法，包括：
通过捕获处理获得伪卫星信号的粗略的载波多普勒估计和粗略的码相位估计；
根据所述粗略的载波多普勒估计得到粗略的码多普勒估计，通过码环迭代更新码多普勒估计和码相位估计，并根据更新的码多普勒估计更新载波多普勒估计；
以更新的载波多普勒估计、更新的码多普勒估计和更新的码相位估计作为跟踪处理的初始值，进行伪卫星信号跟踪处理。
8.如权利要求7所述的接收方法，进一步包括：
根据码多普勒和载波多普勒之间的关系，将粗略的载波多普勒估计转换为粗略的码多普勒估计；
以所述粗略的载波多普勒估计、粗略的码多普勒估计和粗略的码相位估计作为初始值，并以码环迭代更新的码多普勒估计作为反馈值复现本地信号；
将复现的本地信号与接收到的伪卫星信号进行相关，获得相关值；
根据所述相关值，鉴别复现的本地信号与接收到的卫星信号之间的码相位偏差；
根据所述码相位偏差，迭代更新码多普勒估计，直到码环迭代更新收敛；
将更新的码多普勒多普勒估计转换为更新的载波多普勒估计。
9.如权利要求7所述的接收方法，其中，码环迭代更新的收敛通过码环迭代时间确定。
10.如权利要求7所述的接收方法，其中，码环迭代更新的收敛通过所述码相位偏差的大小确定。
CN201610395814.9A 2016-06-07 2016-06-07 伪卫星接收机和伪卫星信号接收方法 Active CN107479071B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610395814.9A CN107479071B (zh) 2016-06-07 2016-06-07 伪卫星接收机和伪卫星信号接收方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610395814.9A CN107479071B (zh) 2016-06-07 2016-06-07 伪卫星接收机和伪卫星信号接收方法

Publications (2)

Publication Number Publication Date
CN107479071A true CN107479071A (zh) 2017-12-15
CN107479071B CN107479071B (zh) 2019-09-20

Family

ID=60594208

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610395814.9A Active CN107479071B (zh) 2016-06-07 2016-06-07 伪卫星接收机和伪卫星信号接收方法

Country Status (1)

Country Link
CN (1) CN107479071B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110737003A (zh) * 2018-07-19 2020-01-31 清华大学 跳时信号捕获装置和方法
CN113758373A (zh) * 2021-09-07 2021-12-07 重庆天箭惯性科技股份有限公司 一种提高弹载接收机定位、测速精度的方法、装置及设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101441260A (zh) * 2008-12-18 2009-05-27 中国科学院微电子研究所 一种全球定位接收机跟踪系统及其跟踪方法
CN102116865A (zh) * 2009-12-30 2011-07-06 中国科学院微电子研究所 一种超高灵敏度的gnss信号的载波跟踪方法
CN103293537A (zh) * 2012-03-05 2013-09-11 安凯(广州)微电子技术有限公司 一种高效的gps数字跟踪方法以及gps数字跟踪环

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101441260A (zh) * 2008-12-18 2009-05-27 中国科学院微电子研究所 一种全球定位接收机跟踪系统及其跟踪方法
CN102116865A (zh) * 2009-12-30 2011-07-06 中国科学院微电子研究所 一种超高灵敏度的gnss信号的载波跟踪方法
CN103293537A (zh) * 2012-03-05 2013-09-11 安凯(广州)微电子技术有限公司 一种高效的gps数字跟踪方法以及gps数字跟踪环

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王修益 等: ""一种北斗伪卫星兼容接收机的设计方案"", 《第四届中国卫星导航学术年会电子文集》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110737003A (zh) * 2018-07-19 2020-01-31 清华大学 跳时信号捕获装置和方法
CN110737003B (zh) * 2018-07-19 2022-03-25 清华大学 跳时信号捕获装置和方法
CN113758373A (zh) * 2021-09-07 2021-12-07 重庆天箭惯性科技股份有限公司 一种提高弹载接收机定位、测速精度的方法、装置及设备
CN113758373B (zh) * 2021-09-07 2023-01-10 重庆天箭惯性科技股份有限公司 一种提高弹载接收机定位、测速精度的方法、装置及设备

Also Published As

Publication number Publication date
CN107479071B (zh) 2019-09-20

Similar Documents

Publication Publication Date Title
CN106772455B (zh) 一种基于惯性信息辅助与参数估计的gnss反欺骗环路跟踪方法
CN104375151B (zh) 导航信号接收机和接收方法
JP3806425B2 (ja) 衛星測位方法及び衛星測位システム
WO2016050055A1 (zh) 一种信号捕获方法、装置及计算机存储介质
US20200271791A1 (en) Tdoa-based positioning system using terrestrial wireless signal sources
CN104280746A (zh) 一种惯性辅助gps的深组合半实物仿真系统
CN104536016A (zh) 一种gnss新体制信号捕获装置及方法
CN101520505B (zh) 全球卫星导航系统弱信号下的自适应相干累积联合捕获方法
CN102183770A (zh) 一种抗多径干扰的gps伪随机码跟踪环路及其抗多径干扰方法
CN105553506B (zh) 一种长码扩频信号的快速捕获方法及装置
US20170285173A1 (en) Method and apparatus for gnss signal tracking
CN104765050A (zh) 一种新的北斗信号二次捕获算法
CN105842713A (zh) 一种基于sft的ins辅助bds信号快速捕获方法
CN103176189A (zh) 高灵敏度卫星导航接收机的远近效应抑制器及其方法
CN104215979A (zh) 一种基于分段相关结合fft运算的大频偏gnss信号捕获方法
CN104735713A (zh) 一种适应于宽带无线定位系统中的精确信号传播时延估计方法及实现装置
CN107479071B (zh) 伪卫星接收机和伪卫星信号接收方法
CN103760578B (zh) 一种gnss卫星导航信号的无模糊跟踪方法
KR101838406B1 (ko) Boc 신호를 추적하는 방법 및 boc 신호 추적 장치
CN105204050B (zh) 一种惯性辅助的多通道混合型矢量跟踪方法
CN105527636A (zh) 导航信号捕获转跟踪方法及系统
CN106291612B (zh) 一种导航卫星星间链路无线信号高性能捕获判决方法
CN103339526B (zh) 设备和方法
CN103983989A (zh) 一种用于卫星定位导航接收系统的数字非整数跟踪环路
JP3837419B2 (ja) 衛星測位方法及び衛星測位システム

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