CN112034491A - 一种基于误差核心分布的完好性保护级计算方法 - Google Patents

一种基于误差核心分布的完好性保护级计算方法 Download PDF

Info

Publication number
CN112034491A
CN112034491A CN202010892642.2A CN202010892642A CN112034491A CN 112034491 A CN112034491 A CN 112034491A CN 202010892642 A CN202010892642 A CN 202010892642A CN 112034491 A CN112034491 A CN 112034491A
Authority
CN
China
Prior art keywords
protection level
integrity
positioning
probability
satellites
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
CN202010892642.2A
Other languages
English (en)
Other versions
CN112034491B (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.)
CETC 54 Research Institute
Original Assignee
CETC 54 Research Institute
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 CETC 54 Research Institute filed Critical CETC 54 Research Institute
Priority to CN202010892642.2A priority Critical patent/CN112034491B/zh
Publication of CN112034491A publication Critical patent/CN112034491A/zh
Application granted granted Critical
Publication of CN112034491B publication Critical patent/CN112034491B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/10Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing dedicated supplementary positioning signals
    • G01S19/11Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing dedicated supplementary positioning signals wherein the cooperating elements are pseudolites or satellite radio beacon positioning system signal repeaters
    • G01S19/115Airborne or satellite based pseudolites or repeaters
    • 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
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/07Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
    • G01S19/071DGPS corrections
    • 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
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/08Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing integrity information, e.g. health of satellites or quality of ephemeris data
    • 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
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/396Determining accuracy or reliability of position or pseudorange measurements
    • 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
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/40Correcting position, velocity or attitude
    • G01S19/41Differential correction, e.g. DGPS [differential GPS]

Abstract

本发明公开了一种基于误差核心分布的完好性保护级计算方法,属于卫星导航中的地基增强系统技术领域。本方法包括以下步骤:统计差分定位误差分布;计算每个卫星数所对应的真实完好性概率,得到真实完好性概率大于完好性风险的最大卫星数N;遍历所有N颗卫星,得到最大的定位范围;遍历所有的概率值P,计算出不同概率值P所对应的保护级和定位结果;按照不同的用户需求,得到最后输出的保护级和定位结果。本方法构造了一种需要少量先验数据支撑的机载完好性保护级计算方法,其通过计算可覆盖完好性的卫星数进行保护级和定位结果的解算,对于GBAS领域具有重要的研究意义和应用价值。

Description

一种基于误差核心分布的完好性保护级计算方法
技术领域
本发明涉及一种基于误差核心分布的完好性保护级计算方法,属于卫星导航中的地基增强系统(GBAS)技术领域。
背景技术
全球卫星导航系统(GNSS)是基于性能的导航(PBN)发展的核心技术,是未来导航服务能力的提升的基础。为改善GNSS信号的性能,实现民航精密进近和着陆引导在精度、完好性、可用性等方面的要求,国际民航组织定义了GNSS增强系统,包括ABAS(空基增强系统)、SBAS(海基增强系统)、GBAS(地基增强系统)等。其中,GBAS是一个可以用于航空飞行器进行精密进近的本地差分卫星导航系统。主要是为机场附近的GNSS信号精度、完好性和可用性提供增强服务。
完好性是导航系统由于故障而不能提供导航信息或者保护级超出告警门限时,系统能够在给定时间内及时向用户发出告警的能力。完好性风险时指系统危险误导信息漏检的概率。完好性是衡量卫星导航系统性能的重要指标,完好性评估技术对于GBAS建设以及评估具有重要的意义。
GBAS完好性方面主要基于完好性通道(GIC)与保护级的技术体制。核心问题在完好性风险与误差包络技术两方面。
完好性风险分配是完好性研究的前提。现有卫星导航增强系统的完好性风险分配是只考虑单接收机故障的情况。在多种导航系统并存的情况下,目前认为两个及两个以上接收机同时故障发生概率为0的假设将不成立。完好性风险是基于故障进行分配的,某个特定故障的风险值为该故障发生的先验概率与故障发生时监测系统对该故障的漏检概率的乘积。但是,由于样本的缺乏和技术手段的不成熟,很多故障的漏检概率和先验概率都难以确定。目前采用的方法主要有两种,一种是基于专家经验或根据陆基导航系统的经验数据直接估计某一类故障的完好性风险值,该方法缺乏理论依据,对于某些故障的完好性风险分配值不是很合理;另一种方法为了保证系统的可用性,对事件发生的频数进行一定的扩大,从而得到一个保守的先验概率,同时用高斯分布包络测试统计量的不确定度得到漏检概率,该方法得到的结果比较保守,实际的风险值和其相比有一定的冗余度。
误差真实误差可能来自具有不同标准差的总体分布,而在进行数据处理时通常需要对数据进行一定的混合,因此导致真实误差呈现厚尾的特性,现有的实际观测数据表明伪距校正误差概率分布核近似为高斯分布,其尾部分布未知,但概率分布远高于高斯。因此必须在计算广播标准差时考虑这种厚尾特性。
完好性评估过程中,由于样本数量有限,只能基于模型外推的方法进行评估,常见的外推模型有高斯尾部包络模型、拉普拉斯尾部包络模型等,这些模型都需要确定尾部阈值参数,但是现有的阈值选取方法是根据经验设定的,缺乏合理性。
因此,构造一种需要少量先验数据支撑的机载完好性保护级计算算法具有重要的研究意义和应用价值。
发明内容
本发明的目的是为了克服了现有GBAS机载完好性保护级计算所需的先验尾部概率分布需要长时间数据统计、保护级结果不能够充分体现定位的偏差误差等问题,提出一种基于误差核心分布的完好性保护级计算方法。
本发明是通过下述技术方案实现的:
一种基于误差核心分布的完好性保护级计算方法,其包括如下步骤:
(101)统计差分定位误差分布,得到不同伪距误差值ΔY对应的概率值P;
(102)设定完好性风险,当前共有M颗卫星,从卫星数为M开始遍历,一直遍历到卫星数为4,计算每个卫星数所对应的真实完好性概率,得到真实完好性概率大于完好性风险的最大卫星数N;
(103)利用当前所有卫星中任意N颗卫星进行定位,对每N颗卫星进行计算得到一个用户的定位范围;遍历所有N颗卫星,取所有定位范围的最大值和最小值,得到最大的定位范围;
(104)根据步骤(103)得到的定位范围,计算得到保护级和定位结果;
(105)按照以上方式,遍历所有的概率值P,计算出不同概率值P所对应的保护级和定位结果;
(106)按照以下两种不同的用户需求,得到最后输出的保护级和定位结果:
i)取垂直保护级最小的概率值P所对应的保护级和定位结果作为最后输出的保护级和定位结果;
ii)取保护级平方根最小的概率值P所对应的保护级和定位结果作为最后输出的保护级和定位结果。
本发明所取得的有益效果在于:
1、本发明方法仅需要少量先验数据支撑,可有效完成机载完好性保护级计算。
2、本发明方法通过计算可覆盖完好性的卫星数进行保护级和定位结果的解算,方法新颖,计算方式简单且易于实现。
具体实施方式
一种基于误差核心分布的完好性保护级计算方法,其包括如下步骤:
(101)统计差分定位误差分布f(ΔY)=P。根据伪距误差ΔY统计概率P,得到不同伪距误差值对应的概率值。
(102)设定完好性风险为Pinter,当前共有M颗卫星。从卫星数为M开始遍历,一直遍历到卫星数4,计算每个卫星数所对应的的真实完好性概率,得到真实完好性概率大于完好性风险的最大卫星数N。
(103)利用当前所有卫星中任意N颗卫星进行定位,每N颗卫星的计算可以得到一个用户的定位范围。遍历所有N颗卫星,取所有定位范围的最大值和最小值,得到最大的定位范围。
(104)根据以上得到的定位范围,计算得到保护级(即PL)和定位结果。
(105)按照以上方式,遍历所有的P,计算出不同概率P所对应的保护级和定位结果。
(106)按照以下两种不同的用户需求,得到最后输出的保护级和定位结果。
i)取VPL(垂直保护级)最小的概率P所对应的保护级和定位结果作为最后输出的保护级和定位结果。
ii)取PL平方根最小的概率P所对应的保护级和定位结果作为最后输出的保护级和定位结果。
以下为一个更具体的例子:
一种基于误差核心分布的完好性保护级计算方法,其包括如下步骤:
(101)每个卫星收取10000个差分伪距误差作为采样数据.根据伪距误差ΔY统计概率P,得到不同伪距误差值对应的概率值。
ΔY 1 0.8 0.5 0.3
P 0.999 0.996 0.993 0.99
(102)设定定位精度完好性风险为Pinter=10-7。假设共10颗卫星。选定P=0.999。
N颗卫星都参与定位,仍满足完好性风险时,下式应该成立:
Figure BDA0002657362640000051
下面通过遍历的方式,计算N颗卫星参与定位的真实完好性概率:
(i)N=10,10颗卫星都参与定位,真实完好性概率为:
1-0.99910>10-7
所以10颗卫星都参与定位不满足完好性风险。
(ii)N=9,9颗卫星都参与定位,真实完好性概率为:
Figure BDA0002657362640000061
所以9颗卫星都参与定位不满足完好性风险。
(iii)N=8,8颗卫星都参与定位,真实完好性概率为:
Figure BDA0002657362640000062
所以8颗卫星都参与定位不满足完好性风险。
(iv)N=7,7颗卫星都参与定位,真实完好性概率为:
Figure BDA0002657362640000063
所以7颗卫星都参与定位满足完好性风险。
(103)任取10颗卫星中的7颗卫星进行定位,通过以下公式计算差分位置误差:
ΔX=S×ΔY
其中:
S=(GTWG)-1GTW
Figure BDA0002657362640000064
Figure BDA0002657362640000071
其中,θi和Azi是第i颗卫星的仰角和方位角,σi 2是第i颗卫星对应的用户测距误差方差。以上数据均可由卫星星历和观测数据计算得到。
i)取1,2,3,4,5,6,7号星,首先使用最小二乘法或者卡尔曼滤波法计算定位结果。然后按照下列公式计算差分位置误差ΔX1,通过差分位置误差修正定位结果得到最终的定位结果范围
Figure BDA0002657362640000072
ii)取1,2,3,4,5,6,8号星,首先使用最小二乘法或者卡尔曼滤波法计算定位结果。然后按照下列公式计算差分位置误差ΔX2,通过差分位置误差修正定位结果得到最终的定位结果范围
Figure BDA0002657362640000073
以此类推,将所有情况遍历一遍,得到120个定位结果范围,从中得到最大的范围
Figure BDA0002657362640000074
(104)根据以上定位结果范围,得到保护级PL0.999和定位结果X0.999
Figure BDA0002657362640000081
Figure BDA0002657362640000082
(105)选定P=0.996,P=0.993,P=0.99,完好性风险为Pinter=10-7。按照上述过程,分别计算保护级和定位结果,得到:
Figure BDA0002657362640000083
Figure BDA0002657362640000084
Figure BDA0002657362640000085
Figure BDA0002657362640000086
Figure BDA0002657362640000087
Figure BDA0002657362640000088
(106)根据用户需求,选择不同的约束条件:
i)约束条件1:对比得到的保护级结果,取最小的VPL对应的保护级和定位结果,作为最后输出的保护级结果。
假设,最大的VPL对应的是P=0.993,即:
VPL0.993=MIN(VPL0.999,VPL0.993,VPL0.996,VPL0.99)
则最后输出的保护级结果PL和定位结果X是:
Figure BDA0002657362640000091
Figure BDA0002657362640000092
ii)约束条件2:按照以下公式对比各保护级结果:
Figure BDA0002657362640000093
取最小的上式结果对应的保护级和定位结果,作为最后输出的保护级结果。
假设:
Figure BDA0002657362640000094
则最后输出的保护级结果PL和定位结果X是:
Figure BDA0002657362640000095
Figure BDA0002657362640000096
总之,本发明构造了一种需要少量先验数据支撑的机载完好性保护级计算方法,其通过计算可覆盖完好性的卫星数进行保护级和定位结果的解算,对于GBAS领域具有重要的研究意义和应用价值。

Claims (1)

1.一种基于误差核心分布的完好性保护级计算方法,其特征在于,包括如下步骤:
(101)统计差分定位误差分布,得到不同伪距误差值ΔY对应的概率值P;
(102)设定完好性风险,当前共有M颗卫星,从卫星数为M开始遍历,一直遍历到卫星数为4,计算每个卫星数所对应的真实完好性概率,得到真实完好性概率大于完好性风险的最大卫星数N;
(103)利用当前所有卫星中任意N颗卫星进行定位,对每N颗卫星进行计算得到一个用户的定位范围;遍历所有N颗卫星,取所有定位范围的最大值和最小值,得到最大的定位范围;
(104)根据步骤(103)得到的定位范围,计算得到保护级和定位结果;
(105)按照以上方式,遍历所有的概率值P,计算出不同概率值P所对应的保护级和定位结果;
(106)按照以下两种不同的用户需求,得到最后输出的保护级和定位结果:
i)取垂直保护级最小的概率值P所对应的保护级和定位结果作为最后输出的保护级和定位结果;
ii)取保护级平方根最小的概率值P所对应的保护级和定位结果作为最后输出的保护级和定位结果。
CN202010892642.2A 2020-08-31 2020-08-31 一种基于误差核心分布的完好性保护级计算方法 Active CN112034491B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010892642.2A CN112034491B (zh) 2020-08-31 2020-08-31 一种基于误差核心分布的完好性保护级计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010892642.2A CN112034491B (zh) 2020-08-31 2020-08-31 一种基于误差核心分布的完好性保护级计算方法

Publications (2)

Publication Number Publication Date
CN112034491A true CN112034491A (zh) 2020-12-04
CN112034491B CN112034491B (zh) 2024-04-09

Family

ID=73587540

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010892642.2A Active CN112034491B (zh) 2020-08-31 2020-08-31 一种基于误差核心分布的完好性保护级计算方法

Country Status (1)

Country Link
CN (1) CN112034491B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115408652A (zh) * 2022-10-31 2022-11-29 北京航空航天大学 一种基于截断误差估计的ekf组合导航保护级确定方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102096075A (zh) * 2010-12-23 2011-06-15 中国航天科工信息技术研究院 基于Galileo系统完好性概念的多模用户完好性评估方法
CN104267410A (zh) * 2014-10-10 2015-01-07 北京航空航天大学 机载完好性监测中多故障的排除方法和装置
CN104732085A (zh) * 2015-03-23 2015-06-24 北京航空航天大学 卫星导航星基增强系统可用性预测方法
CN105068088A (zh) * 2015-06-29 2015-11-18 北京航空航天大学 双频卫星导航星基增强系统可用性预测方法
CN108761497A (zh) * 2018-03-06 2018-11-06 北京交通大学 一种列车卫星定位局域完好性监测方法及系统
US20190187295A1 (en) * 2017-12-18 2019-06-20 Korea Advanced Institute Of Science And Technology (Kaist) Method and system for local-area differential gnss for uav navigation, and for generating optimal protection level and geometry screening therefor

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102096075A (zh) * 2010-12-23 2011-06-15 中国航天科工信息技术研究院 基于Galileo系统完好性概念的多模用户完好性评估方法
CN104267410A (zh) * 2014-10-10 2015-01-07 北京航空航天大学 机载完好性监测中多故障的排除方法和装置
CN104732085A (zh) * 2015-03-23 2015-06-24 北京航空航天大学 卫星导航星基增强系统可用性预测方法
CN105068088A (zh) * 2015-06-29 2015-11-18 北京航空航天大学 双频卫星导航星基增强系统可用性预测方法
US20190187295A1 (en) * 2017-12-18 2019-06-20 Korea Advanced Institute Of Science And Technology (Kaist) Method and system for local-area differential gnss for uav navigation, and for generating optimal protection level and geometry screening therefor
CN108761497A (zh) * 2018-03-06 2018-11-06 北京交通大学 一种列车卫星定位局域完好性监测方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
于耕,王寒,赵龙: "基于卡尔曼滤波的卫星地基增强系统 位置域完好性监测分析", 科学技术与工程, vol. 18, no. 26, 30 September 2018 (2018-09-30), pages 243 - 247 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115408652A (zh) * 2022-10-31 2022-11-29 北京航空航天大学 一种基于截断误差估计的ekf组合导航保护级确定方法

Also Published As

Publication number Publication date
CN112034491B (zh) 2024-04-09

Similar Documents

Publication Publication Date Title
CN110007317B (zh) 一种选星优化的高级接收机自主完好性监测方法
CN109061683B (zh) 一种优化水平保护级的h-araim系统
CN110196434B (zh) 一种高级接收机自主完好性监测的星座动态选择方法
EP3206048B1 (en) Use of wide area reference receiver network data to mitigate local area error sources
CN104318089A (zh) 一种用于局域增强系统完好性监测的门限值确定方法
CN104950320A (zh) 一种监视地基增强系统对流层修正参数的方法及系统
Marini-Pereira et al. Ground-based augmentation systems operation in low latitudes-part 1: Challenges, mitigations, and future prospects
Lee et al. GPS and Galileo with RAIM or WAAS for Vertically Guided Approaches
CN110174683B (zh) 基于贝叶斯决策的gbas保护级完好性风险分配方法
CN112034491B (zh) 一种基于误差核心分布的完好性保护级计算方法
CN112198533B (zh) 一种多假设下的地基增强系统完好性评估系统及方法
CN109581430A (zh) 一种基于伪卫星的gbas电离层空间梯度的监测方法
CN111948675B (zh) 一种地基增强卫星导航系统的伪卫星布设方法
Kuna et al. Preliminary analysis of standalone Galileo and NavIC in the context of positioning performance for low latitude region
Caamano et al. Performance Evaluation of the Ionospheric Threat Mitigation Strategies in Dual-Frequency Multi-Constellation GBAS
Ammana et al. Advantage of irnss s-band signal for gbas applications in adverse ionospheric storm conditions
Jochems et al. GBAS use cases beyond what was envisioned–drone navigation
Ma et al. Progress of global ARAIM availability of BDS-2/BDS-3 with TGD and ISB
Madonna et al. NIORAIM algorithm applied to a multiconstellation GNSS: Analysis of integrity monitoring performances in various phases of flight
Bang et al. Preliminary availability assessment to support single-frequency SBAS development in the Korean region
Hu et al. Multiple reference consistency check algorithm in GBAS based on S-values auxiliary
Pandya et al. Using Kriging to optimize WAAS performance over the entire solar cycle
Nietlispach et al. Analysis of the Occurrence of Critical Satellites in GAST-D+ Processing
Saito et al. Impact assessment of ionospheric scintillation associated with plasma bubbles on GAST-D ground integrity monitors
Misra et al. Augmentation of GPS/LAAS with GLONASS: Performance assessment

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