CN114879229A - 一种适用于gbas***的卫星导航信号畸变检测方法 - Google Patents
一种适用于gbas***的卫星导航信号畸变检测方法 Download PDFInfo
- Publication number
- CN114879229A CN114879229A CN202210580017.3A CN202210580017A CN114879229A CN 114879229 A CN114879229 A CN 114879229A CN 202210580017 A CN202210580017 A CN 202210580017A CN 114879229 A CN114879229 A CN 114879229A
- Authority
- CN
- China
- Prior art keywords
- signal
- intermediate frequency
- zero intermediate
- frequency baseband
- sequence
- 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
Links
Images
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/23—Testing, monitoring, correcting or calibrating of receiver elements
-
- 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/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE 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/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种适用于GBAS***的卫星导航信号畸变检测方法,针对机场GBAS***监测卫星导航信号的场景,首先利用GBAS地面参考接收机完成卫星导航信号的接收,将卫星导航数字化信号进行载波多普勒剥离,利用本地伪码判决得到每个积分周期的导航电文符号,实现对多个积分周期的剥离电文后零中频基带信号累加,利用累加后的零中频基带信号完成相关曲线计算,再将相关曲线进行不同采样长度延迟的循环移位处理,与未移位的相关曲线相减,提取相减结果中相关峰值所在区域的序列并进行拟合求根,最终通过比较判决完成监测信号畸变的检测告警。本发明可快速检测影响定位服务的卫星导航信号畸变,提升了机场GBAS***的应用安全,保证了高精度高可信陆基增强引导。
Description
技术领域
本发明涉及卫星导航信号质量监测评估技术领域,具体为一种适用于GBAS***的卫星导航信号畸变检测方法。
背景技术
GBAS(Ground Based Augmentation System,地基增强***),最早由国际民航组织提出,从现有陆基导航***过渡到星基导航***,通过卫星导航差分定位提升定位导航精度,并增加了一系列完好性监视算法,提高***完好性、可用性、连续性的指标,使机场覆盖空域范围内的配置相应机载设备的飞机获得到达I类精密进近(CAT-I)甚至更高标准的精密进近、着陆引导服务。
GBAS***中最重要的设施就是地面参考站,主要完成对所在区域可见的导航卫星信号的监测接收,并形成差分信息播发给卫星。然而由于导航卫星信号在生成过程中,可能会由于链路器件受环境变化和元器件非线性影响,将产生频率色散、时延缓变、频率/通道一致性失真,导致真实信号的多维度属性分析结果模糊化,严重影响地面参考站生成差分信息的质量,给进近引导定位带来很大的风险。
因此本专利发明一种适用于GBAS***的卫星导航信号畸变检测方法,可以有效的快速监测导航卫星信号畸变并及时告警,保障机载用户的定位精度和完好性。
发明内容
本发明的目的在于提供一种适用于GBAS***的卫星导航信号畸变检测方法。本方法主要针对机场GBAS***监测卫星导航信号的场景,快速检测影响定位服务的卫星导航信号畸变,提升了机场GBAS***的应用安全,保证了高精度高可信陆基增强引导。具体方案如下:
一种适用于GBAS***的卫星导航信号畸变检测方法,包括以下步骤:
步骤1:利用GBAS地面参考接收机完成卫星导航信号的接收,并将通过AD芯片的数字化信号进行载波多普勒剥离,得到零中频基带信号;
步骤2:将每个积分周期的零中频基带信号和本地伪码作相乘求和,判决得到每个积分周期的导航电文符号,并根据符号正负剥离零中频基带信号中的导航电文,并对剥离电文后多个积分周期的零中频基带信号进行累加;
步骤3:将累加后的零中频基带信号与本地伪码作相关运算,得到相关曲线,并对相关曲线进行内插扩充处理;
步骤4:将内插扩充处理发相关曲线进行不同采样长度延迟的循环移位处理,并与移位前的相关曲线相减,提取相减结果中相关峰值所在区域的序列,并对该区域进行拟合,计算拟合方程的根;
步骤5:选取所有不同采样长度延迟下的拟合方程根结果最大值,如果最大值大于设定阈值,则GBAS***完成该时刻监测信号畸变的检测告警。
进一步的,所述步骤1中:
GBAS地面参考接收机AD芯片原始采样的数字化信号为A(n),剥离载波后的零中频基带信号为:
更进一步的,所述步骤2具体包括步骤:
步骤2.1:选取剥离载波后的零中频基带信号B(n)的一个积分周期和本地伪码C(n)作相乘求和,并计算判决得到每个积分周期的导航电文符号:
其中,N为一个积分周期的信号序列长度,M为当前接收的零中频基带信号的第M个积分周期;B(i)为信号周期数为i时的剥离载波后零中频基带信号,i为信号周期数的序号,C(mod(i,N))表示信号周期数为mod(i,N)时的本地伪码,|·|表示绝对值,mod表示取模值;
步骤2.2:选取计算得到的导航电文符号D(M),与剥离载波后的零中频基带信号B(n)实现对多个积分周期的剥离电文后零中频基带信号累加,得到累加后的零中频基带信号E(n),E(n)序列中每个数值的计算方法为:
其中,j为零中频基带信号的周期数序号。
更进一步的,所述步骤3具体包括步骤:
步骤3.1:将步骤2.2得到的累加后的零中频基带信号E(n)与本地伪码C(n)作相关运算,得到相关曲线F(n),计算方法为:
步骤3.2:将所述相关曲线F(n)进行内插扩充处理,得到扩充后的相关曲线H(n),计算方法为:
其中,G(n)表示相关曲线的傅里叶变换结果的扩充序列。
更进一步的,所述步骤4具体包括步骤:
步骤4.1:选取不同的采样长度延迟对相关曲线H(n)依据不同的采样长度延迟进行循环移位处理,得到循环移位序列:
其中,采样长度延迟Delay(k)=N×(k-6)/L,k=1,2,…10,共11个延迟,其中L为本地伪码的码长;
将不同的循环移位序列与未移位的相关曲线相减,得到相减结果:
Qk(n)=Pk(n)-H(n)
步骤4.2:在所述相减结果Qk(n)中提取相关峰所在序列:
Rk(x)=Qk(x),x=1,2,…N/L
构造相关峰码片间隔序列:
S(x)=x-N/L/2,x=1,2,…N/L
其中,x表示相关峰所在序列的值的序号;
实现对相关峰码片间隔序列S(x)与相关峰所在序列Rk(x)的线性拟合,得到的不同的采样长度延迟Delay(k)对应的拟合系数Uk、Vk计算结果为:
其中,Rk(i)表示相关峰所在序列的第i个值;S(i)表示相关峰码片间隔序列的第i个值;
步骤4.3:利用所述拟合系数Uk、Vk得到拟合方程的根Wk=-Vk/Uk。
更进一步的,所述步骤5中,计算拟合方程的根Wk中的最大值,并与0.5比较,如果大于0.5,则该时刻检测到威胁到定位服务的卫星导航信号畸变,从而实现信号质量监测告警。
本发明与现有技术相比的有益效果为:
实现在GBAS***内的导航信号畸变快速检测,有效识别因卫星环境变化和元器件非线性引入的频率色散、时延缓变、频率/通道一致性失真等对信号质量的影响,实现对异常卫星信号定位计算排除,保证高精度高可信的GBAS***进近引导服务。
附图说明
图1是一种适用于GBAS***的卫星导航信号畸变检测方法的流程图。
具体实施方式
以下结合附图和实施例对本发明做进一步的说明。
图1为一种适用于GBAS***的卫星导航信号畸变检测方法的实现流程图。本发明的方法是这样实现的:本方法首先利用GBAS地面参考接收机完成卫星导航信号的接收,将卫星导航数字化信号进行载波多普勒剥离,利用本地伪码判决得到每个积分周期的导航电文符号,实现对多个积分周期的剥离电文后零中频基带信号累加,利用累加后的零中频基带信号完成相关曲线计算,然后将相关曲线进行不同采样长度延迟的循环移位处理,与未移位的相关曲线相减,提取相减结果中相关峰值所在区域的序列并进行拟合求根,最终通过比较判决的方式完成监测信号畸变的检测告警。具体包括以下步骤:
(1)利用GBAS地面参考接收机完成卫星导航信号的接收,并将通过AD芯片的数字化信号进行载波多普勒剥离,得到零中频基带信号。
GBAS地面参考接收机通过卫星导航天线获取卫星导航无线信号,信号通过接收机的射频前端后进入AD芯片,由AD芯片将模拟信号转化为原始采样的数字化信号为A(n)。在导航信号的接收过程中,跟踪环路完成多普勒频率、相位的估计,利用估计值重构载波后,剥离载波后得到零中频基带信号其中fc为信号中心频率,fd为接收机估计信号多普勒频率,为接收机估计载波相位,fs为信号采样率。
(2)将每个积分周期的零中频基带信号和本地伪码作相乘求和,判决得到每个积分周期的导航电文符号,并根据符号正负剥离零中频基带信号中的导航电文,并实现对多个积分周期的剥离电文后零中频基带信号累加。
选取上一步得到的零中频基带信号B(n)的一个积分周期和本地伪码C(n)作相乘求和,并计算判决得到每个积分周期的导航电文符号:
其中N为一个积分周期的信号序列长度,M为当前接收的零中频基带信号的第M个积分周期。
选取上一步计算得到的电文符号D(M),与零中频基带信号B(n)进实现对多个积分周期的剥离电文后零中频基带信号累加,得到累加后的零中频基带信号E(n),E(n)序列中每个数值的计算方法为:
(3)将累加后的零中频基带信号与本地伪码作相关运算,得到相关曲线,并对相关曲线进行内插扩充处理。
将上一步得到的累加后的零中频基带信号E(n)与本地伪码C(n)作相关运算,得到相关曲线F(n),计算方法为:
将上一步得到的相关曲线进行内插扩充处理,得到扩充后的相关曲线H(n),计算方法为:
其中,
(4)将相关曲线进行不同采样长度延迟的循环移位处理,并与未移位的相关曲线相减,提取相减结果中相关峰值所在区域的序列,并对该区域进行拟合,计算拟合方程的根。
选取不同的采样长度延迟Delay(k)=N×(k-6)/L,k=1,2,…10,共11个延迟,其中L为本地伪码的码长,对相关曲线H(n)依据不同的采样长度延迟进行循环移位处理,得到循环移位序列:
将不同的循环移位序列与未移位的相关曲线相减,得到相减结果:
Qk(n)=Pk(n)-H(n)
在上一步计算的相减结果Qk(n)中提取相关峰所在序列:
Rk(x)=Qk(x),x=1,2,…N/L
构造相关峰码片间隔序列:
S(x)=x-N/L/2,x=1,2,…N/L
实现对S(x)与Rk(x)序列的线性拟合,得到的不同的采样长度延迟Delay(k)对应的拟合系数Uk、Vk计算结果为:
利用在上一步计算的拟合系数Uk、Vk可以得到拟合方程的根Wk=-Vk/Uk;
(5)选取所有不同采样长度延迟下的拟合方程根结果最大值,如果最大值大于0.5,则GBAS***完成该时刻监测信号畸变的检测告警。
计算Wk中的最大值,将最大值与0.5比较,如果大于0.5,则证明相关峰两侧的不对称性存在对导航定位产生较大偏差的影响,则该时刻检测到威胁到定位服务的信号畸变,此时需要在GBAS***内实现对该颗卫星的告警,在GBAS***生成差分信息时会排除到此颗卫星,从而实现信号质量监测告警和定位结果可靠性的提升。
Claims (6)
1.一种适用于GBAS***的卫星导航信号畸变检测方法,其特征在于,包括以下步骤:
步骤1:利用GBAS地面参考接收机完成卫星导航信号的接收,并将通过AD芯片的数字化信号进行载波多普勒剥离,得到零中频基带信号;
步骤2:将每个积分周期的零中频基带信号和本地伪码作相乘求和,判决得到每个积分周期的导航电文符号,并根据符号正负剥离零中频基带信号中的导航电文,并对剥离电文后多个积分周期的零中频基带信号进行累加;
步骤3:将累加后的零中频基带信号与本地伪码作相关运算,得到相关曲线,并对相关曲线进行内插扩充处理;
步骤4:将内插扩充处理发相关曲线进行不同采样长度延迟的循环移位处理,并与移位前的相关曲线相减,提取相减结果中相关峰值所在区域的序列,并对该区域进行拟合,计算拟合方程的根;
步骤5:选取所有不同采样长度延迟下的拟合方程根结果最大值,如果最大值大于设定阈值,则GBAS***完成该时刻监测信号畸变的检测告警。
3.根据权利要求2所述的适用于GBAS***的卫星导航信号畸变检测方法,其特征在于,所述步骤2具体包括步骤:
步骤2.1:选取剥离载波后的零中频基带信号B(n)的一个积分周期和本地伪码C(n)作相乘求和,并计算判决得到每个积分周期的导航电文符号:
其中,N为一个积分周期的信号序列长度,M为当前接收的零中频基带信号的第M个积分周期;B(i)为信号周期数为i时的剥离载波后零中频基带信号,i为信号周期数的序号,C(mod(i,N))表示信号周期数为mod(i,N)时的本地伪码,|·|表示绝对值,mod表示取模值;
步骤2.2:选取计算得到的导航电文符号D(M),与剥离载波后的零中频基带信号B(n)实现对多个积分周期的剥离电文后零中频基带信号累加,得到累加后的零中频基带信号E(n),E(n)序列中每个数值的计算方法为:
其中,j为零中频基带信号的周期数序号。
5.根据权利要求4所述的适用于GBAS***的卫星导航信号畸变检测方法,其特征在于,所述步骤4具体包括步骤:
步骤4.1:选取不同的采样长度延迟对相关曲线H(n)依据不同的采样长度延迟进行循环移位处理,得到循环移位序列:
其中,采样长度延迟Delay(k)=N×(k-6)/L,k=1,2,…10,共11个延迟,其中L为本地伪码的码长;
将不同的循环移位序列与未移位的相关曲线相减,得到相减结果:
Qk(n)=Pk(n)-H(n)
步骤4.2:在所述相减结果Qk(n)中提取相关峰所在序列:
Rk(x)=Qk(x),x=1,2,…N/L
构造相关峰码片间隔序列:
S(x)=x-N/L/2,x=1,2,…N/L
其中,x表示相关峰所在序列的值的序号;
实现对相关峰码片间隔序列S(x)与相关峰所在序列Rk(x)的线性拟合,得到的不同的采样长度延迟Delay(k)对应的拟合系数Uk、Vk计算结果为:
其中,Rk(i)表示相关峰所在序列的第i个值;S(i)表示相关峰码片间隔序列的第i个值;步骤4.3:利用所述拟合系数Uk、Vk得到拟合方程的根Wk=-Vk/Uk。
6.根据权利要求1所述的适用于GBAS***的卫星导航信号畸变检测方法,其特征在于,所述步骤5中,计算拟合方程的根Wk中的最大值,并与0.5比较,如果大于0.5,则该时刻检测到威胁到定位服务的卫星导航信号畸变,从而实现信号质量监测告警。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210580017.3A CN114879229B (zh) | 2022-05-25 | 2022-05-25 | 一种适用于gbas***的卫星导航信号畸变检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210580017.3A CN114879229B (zh) | 2022-05-25 | 2022-05-25 | 一种适用于gbas***的卫星导航信号畸变检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114879229A true CN114879229A (zh) | 2022-08-09 |
CN114879229B CN114879229B (zh) | 2023-07-18 |
Family
ID=82676774
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210580017.3A Active CN114879229B (zh) | 2022-05-25 | 2022-05-25 | 一种适用于gbas***的卫星导航信号畸变检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114879229B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6987820B1 (en) * | 2000-10-18 | 2006-01-17 | Honeywell International, Inc. | Apparatus for navigation satellite signal quality monitoring |
CN102436000A (zh) * | 2011-11-24 | 2012-05-02 | 北京航空航天大学 | 信号质量监测方法和装置 |
US20200049831A1 (en) * | 2018-08-08 | 2020-02-13 | Thales Management & Services Deutschland Gmbh | Method for operating a plurality of gnss receivers for detecting satellite signal deformation |
CN111190198A (zh) * | 2020-01-14 | 2020-05-22 | 中国民用航空总局第二研究所 | Gbas测试设备的卫星接收机及其伪距生成方法与*** |
CN113126131A (zh) * | 2021-03-25 | 2021-07-16 | 中国电子科技集团公司第五十四研究所 | 一种超低失真导航信号采集及无混叠分离方法 |
CN113281786A (zh) * | 2021-06-02 | 2021-08-20 | 中国电子科技集团公司第五十四研究所 | 全空域多级触发式gnss空间信号质量监测评估方法 |
-
2022
- 2022-05-25 CN CN202210580017.3A patent/CN114879229B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6987820B1 (en) * | 2000-10-18 | 2006-01-17 | Honeywell International, Inc. | Apparatus for navigation satellite signal quality monitoring |
CN102436000A (zh) * | 2011-11-24 | 2012-05-02 | 北京航空航天大学 | 信号质量监测方法和装置 |
US20200049831A1 (en) * | 2018-08-08 | 2020-02-13 | Thales Management & Services Deutschland Gmbh | Method for operating a plurality of gnss receivers for detecting satellite signal deformation |
CN110824506A (zh) * | 2018-08-08 | 2020-02-21 | 塔莱斯管理与服务德国有限责任公司 | 操作用于检测卫星信号变形的多个gnss接收器的方法 |
CN111190198A (zh) * | 2020-01-14 | 2020-05-22 | 中国民用航空总局第二研究所 | Gbas测试设备的卫星接收机及其伪距生成方法与*** |
CN113126131A (zh) * | 2021-03-25 | 2021-07-16 | 中国电子科技集团公司第五十四研究所 | 一种超低失真导航信号采集及无混叠分离方法 |
CN113281786A (zh) * | 2021-06-02 | 2021-08-20 | 中国电子科技集团公司第五十四研究所 | 全空域多级触发式gnss空间信号质量监测评估方法 |
Non-Patent Citations (3)
Title |
---|
刘建成;范建军;冯晓超;桑怀胜;: "卫星导航信号畸变导致的测距偏差的估计方法", 宇航学报, no. 11 * |
刘瑞: "卫星接收机信号处理关键技术的研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
郝小丽;李锐;: "卫星自主完好性相关峰监视", 航空计算技术, no. 06 * |
Also Published As
Publication number | Publication date |
---|---|
CN114879229B (zh) | 2023-07-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5247056B2 (ja) | 伝搬遅延時間測定装置およびレーダ装置 | |
RU2669700C1 (ru) | Радиоприемник для определения местоположения источника сигнала | |
EP1972962A2 (en) | Transmitter independent techniques to extend the performance of passive coherent location | |
US5874916A (en) | Frequency selective TDOA/FDOA cross-correlation | |
US6370207B1 (en) | Method for mitigating multipath effects in radio systems | |
Lopez-Risueno et al. | Atomic decomposition-based radar complex signal interception | |
US7920093B2 (en) | Methods for improving computational efficiency in a global positioning satellite receiver | |
CN110673168B (zh) | 异步多用户联合欺骗信号检测方法及装置 | |
Barkat | Instantaneous frequency estimation of nonlinear frequency-modulated signals in the presence of multiplicative and additive noise | |
KR101674747B1 (ko) | I/q 불균형을 보상하기 위한 레이더 신호 처리 장치 및 방법 | |
CN115327581B (zh) | 识别gnss信号中欺骗信号的方法、装置和接收机 | |
CN104007421A (zh) | 基于全变差和压缩感知的罗兰c被动雷达toa估计方法 | |
CN111865865B (zh) | 适用于高灵敏度星载ads-b接收机的频偏及相偏估计方法 | |
EP1916540B1 (en) | Gps positioning method and gps position device | |
KR20110032982A (ko) | 도착 시간 차이 및 도착 주파수 차이를 추정하는 방법 및 장치 | |
JP2004340855A (ja) | 衛星測位システム及び衛星測位方法 | |
CN117060985B (zh) | 一种船载双天线pcma***信号重捕方法及装置 | |
CN109669157A (zh) | 基于压缩传输的tdoa估计方法 | |
CN114879229B (zh) | 一种适用于gbas***的卫星导航信号畸变检测方法 | |
CN110749905B (zh) | 单星低复杂度卫星导航欺骗信号检测识别方法及装置 | |
Mao et al. | New code delay compensation algorithm for weak GPS signal acquisition | |
KR101644560B1 (ko) | 통신 신호에 대한 tdoa/fdoa 정보 추정 장치 및 방법 | |
Nezhadshahbodaghi et al. | A new method to extend differential integration for weak GPS signal acquisition | |
RU2338261C1 (ru) | Способ регистрации и обработки телеметрического сигнала с временным разделением каналов | |
Darvin et al. | Analysis of autocorrelation based frequency measurement algorithm for IFM receivers |
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 |