CN109116392A - 复杂环境条件下多星座gnss融合高精度动态定位方法 - Google Patents
复杂环境条件下多星座gnss融合高精度动态定位方法 Download PDFInfo
- Publication number
- CN109116392A CN109116392A CN201810824778.2A CN201810824778A CN109116392A CN 109116392 A CN109116392 A CN 109116392A CN 201810824778 A CN201810824778 A CN 201810824778A CN 109116392 A CN109116392 A CN 109116392A
- Authority
- CN
- China
- Prior art keywords
- fuzziness
- satellite
- gnss
- subset
- integer ambiguity
- 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
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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining 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/42—Determining position
- G01S19/421—Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system
- G01S19/423—Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system by combining or switching between position solutions derived from different satellite radio beacon positioning systems
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
本发明公开了复杂环境条件下多星座GNSS融合高精度动态定位方法,包括以下步骤:步骤1:获取浮点模糊度解集合及其方差‑协方差阵对和进行分解得到整周模糊度值,若整周模糊度值满足设定要求则转入步骤4,否则转入步骤2;步骤2:从中剔除第i颗卫星的模糊度,然后进行模糊度分解,得到Ratio值ratioi和浮点模糊度与最优整周模糊度间的距离sni;步骤3:将ratioi和sni根据设定阈值分组后确定最优模糊度子集,根据最优模糊度子集进行整周模糊度解算得到整周模糊度值;步骤4:根据整周模糊度值进行GNSS固定解的解算,得到GNSS动态定位结果;本发明能有效避免偏差较大的浮点模糊度对模糊度分解造成的影响,实现模糊度的快速固定,确保GNSS定位结果的正确性和可靠性。
Description
技术领域
本发明涉及高精度全球卫星导航定位技术领域,具体涉及一种适用于复杂环境条件下的多星座GNSS融合高精度动态定位方法。
背景技术
实现GNSS快速高精度动态定位的先决条件是模糊度快速、正确分解。随着多星座(包括:北斗、GPS、GLONASS、Galileo等)卫星导航系统的发展,可视卫星数越来越多,多星座GNSS融合定位成为重要发展方向。多系统融合在提高定位结果精度和可靠性的同时,也改善了模糊度分解的条件,有利于提高模糊度分解的成功率。由于GNSS信号遮挡是影响GNSS定位的至关重要因素,尤其是在复杂环境下(高楼、桥梁、深山峡谷等强信号反射环境),极易造成卫星信号被遮挡或多源反射,多路径误差更为严重,从而严重影响模糊度的快速正确分解,影响GNSS进行快速高精度动态定位。在复杂环境条件下利用GNSS进行定位时,由于受到来自观测值中噪声及各种未改正偏差的影响,浮点模糊度解中会含有偏差。而少数含有较大偏差的模糊度会严重影响模糊度的分解,甚至造成模糊度分解失败。与采用常规的全模糊度分解方法相比,部分模糊度分解方法通过选择最优模糊度子集进行分解。该模糊度子集中剔除了含偏差较大的模数度,具有较高的精度和可靠性,有利于提高模糊度分解成功的概率,实现GNSS快速高精度动态定位。因此,如何构建选择模糊度子集的准则,将直接影响到所选择模糊度子集的可靠性,也直接影响到模糊度分解的成功与否。
在选择模糊度子集的过程中,根据构建的准则剔除或保留某颗卫星的模糊度。当选择的模糊度子集满足模糊度分解成功的各项条件时,则该模糊度子集将被用于卫星导航定位固定解的解算。在部分模糊度分解方法中,常用的选择模糊度子集的准则有卫星高度角、浮点模糊度方法、模糊度分解成功率等。将卫星高度角作为选择模糊度子集的准则是剔除高度角最小的卫星,并对剩余卫星的模糊度进行分解。其原理是建立在卫星高度角越低,观测值精度也越低的假设之上。但是,在实际情况下,由于受到观测条件的影响,使得观测值精度并不一定与卫星高度角一致。这是由于卫星高度角仅仅是从理论的角度表征观测值精度的,而未考虑实际观测条件的影响。根据浮点模糊度方差的大小选择模糊度子集,该方法只考虑了单个模糊度的精度,而忽略了该模糊度与其它模糊度之间的相关关系。以模糊度分解的成功率作为准则选择模糊度子集,虽然考虑了模糊度之间的相关关系,但是却没有考虑浮点模糊度的因素。因此,在复杂的城市观测环境下进行GNSS快速高精度动态定位时,采用卫星高度角、浮点模糊度解的方差大小、模糊度分解成功率作为模糊度子集的选择准则并不能确保快速地选择一组可靠的模糊度子集。
发明内容
本发明提供一种有效避免偏差较大的浮点模糊度对模糊度分解造成的影响,实现模糊度的快速固定,确保GNSS快速高精度动态定位结果的正确性和可靠性的适用于复杂环境条件下的多星座GNSS融合高精度动态定位方法。
本发明采用的技术方案是:复杂环境条件下多星座GNSS融合高精度动态定位方法,包括以下步骤:
步骤1:获取浮点模糊度解集合及其方差-协方差阵对和进行分解得到整周模糊度值,若整周模糊度值满足设定要求则转入步骤4,否则转入步骤2;
步骤2:从中剔除第i颗卫星的模糊度,然后进行模糊度分解,得到Ratio值ratioi和浮点模糊度与最优整周模糊度间的Mahalanobis距离sni;
步骤3:将步骤2得到的ratioi和sni根据设定阈值分组后确定最优模糊度子集,根据最优模糊度子集进行整周模糊度解算得到整周模糊度值;
步骤4:根据整周模糊度值进行GNSS固定解的解算,得到GNSS动态定位结果。
进一步的,所述步骤3中确定最优模糊度子集的过程如下:
将ratioi大于设定阈值的模糊度子集作为候选模糊度子集;
若候选模糊度子集个数不为零,则选择sni值最小的候选模糊度子集作为最优模糊度子集;
若候选模糊度子集个数为零,则将ratioi中最大值对应的浮点模糊度子集取代重复步骤2~3;直至剩余卫星个数小于设定值为止。
进一步的,所述步骤2中ratioi的计算过程如下:
式中:为剔除第i颗卫星的模糊度后剩余卫星的浮点模糊度解集合;为剔除第i颗卫星的模糊度后剩余卫星的方差-协方差阵,为的逆矩阵,为模糊度分解后获得的整周模糊度最优值,为分解后获得的整周模糊度次优值,q为卫星数。
进一步的,所述步骤2中sni的计算过程如下:
式中:为剔除第i颗卫星的模糊度后剩余卫星的浮点模糊度解集合;为剔除第i颗卫星的模糊度后剩余卫星的方差-协方差阵,为的逆矩阵,为模糊度分解后获得的整周模糊度最优值,q为卫星数。
进一步的,所述步骤1中通过卡尔曼滤波方法对观测数据进行处理得到浮点模糊度解集合及其方差-协方差阵
进一步的,所述步骤4中通过将整周模糊度值代入观测方程,利用卡尔曼滤波方法进行GNSS固定解的解算。
本发明的有益效果是:
(1)本发明在顾及模糊度残差及模糊度间协方差的基础上,以模糊度分解的Ratio值及浮点模糊度与最优整周模糊度间的Mahalanobis距离为指标构建模糊度子集,确保快速地选择一组可靠的模糊度子集;
(2)本发明不仅能够提高模糊度正确固定率,而且可以降低被剔除卫星的个数,快速实现GNSS高精度动态定位;
(3)本发明能够有效避免偏差较大的浮点模糊度对模糊度分解造成的影响,实现模糊度的快速固定,确保GNSS快速高精度动态定位结果的正确性和可靠性。
附图说明
图1为观测数据I本发明方法和传统方法模糊度正确固定率对比图。
图2为观测数据II本发明方法和传统方法模糊度正确固定率对比图。
图3为观测数据I和观测数据II采用本发明和传统方法剔除卫星个数的统计信息对比图。
图4为观测数据I本发明动态定位结果图。
图5为观测数据II本发明动态定位结果图。
具体实施方式
下面结合附图和具体实施例对本发明进一步说明。
复杂环境条件下多星座GNSS融合高精度动态定位方法,包括以下步骤:
步骤1:获取浮点模糊度解集合及其方差-协方差阵对和进行分解得到整周模糊度值,若整周模糊度值满足设定要求则转入步骤4,否则转入步骤2;
通过扩展卡尔曼滤波方法对观测数据进行处理得到浮点模糊度解集合及其方差-协方差阵并对其进行分解。若模糊度分解结果满足要求,则转入步骤4;若模糊度分解结果不能满足要求,则转入步骤2。
步骤2:从中剔除第i颗卫星的模糊度,然后对剩余卫星的模糊度进行分解,得到Ratio值ratioi和浮点模糊度与最优整周模糊度间的Mahalanobis距离sni;
ratioi的计算过程如下:
sni的计算过程如下:
式中:为剔除第i颗卫星的模糊度后剩余卫星的浮点模糊度解集合;为剔除第i颗卫星的模糊度后剩余卫星的方差-协方差阵,为的逆矩阵,为模糊度分解后获得的整周模糊度最优值,为分解后获得的整周模糊度次优值,q为卫星数。
步骤3:将步骤2得到的ratioi和sni根据设定阈值分组后确定最优模糊度子集,根据最优模糊度子集进行整周模糊度解算得到整周模糊度值;
确定最优模糊度子集的过程如下:
将ratioi大于设定阈值的模糊度子集作为候选模糊度子集;
根据设置的Ratio值阈值ratio_threshold,将ratioi分为两组,并将Ratio值大于阈值ratio_threshold的模糊度子集作为候选模糊度子集;
ratio′={ratiok|ratiok≥ratio_threshold,k=1,…,q}
ratio″={ratiok|ratiok<ratio_threshold,k=1,…,q}
对应的,sni也分为两组:
sn′={snk|ratiok≥ratio_threshold,k=1,…,q}
sn″={snk|ratiok<ratio_threshold,k=1,…,q}
若候选模糊度子集个数不为零,则选择sn′中最小的Mahalanobis距离对应的候选模糊度子集作为最优模糊度子集,使用该模糊度子集进行整周模糊度的解算;
若候选模糊度子集个数为零,则将ratio″中最大的Ratio值对应的浮点模糊度子集取代重复步骤2~3,直至剩余卫星个数小于设定值为止。当剩余卫星个数不足以计算固定解时,该循环结束,此时所有模糊度均无法固定,只能得到基线浮点定位结果。
步骤4:根据整周模糊度值进行GNSS固定解的解算,得到GNSS动态定位结果;
在采用全模糊度(步骤1)或部分模糊度(步骤2~3)分解方式分解模糊度成功后,将固定的整周模糊度值代入观测方程,利用卡尔曼滤波方法进行GNSS固定解的解算,得到GNSS动态定位结果。
为了验证本发明的有效性,对七天的观测数据I和观测数据II进行处理,并与根据卫星高度角大小选择模糊度子集的方法进行了对比,结果如图1和图2所示。
图1和图2分别是观测数据I和观测数据II采用本发明方法和根据卫星高度角选择模糊度子集方法的模糊度正确固定率对比结果图。由图1和图2可知,与根据卫星高度角选择模糊度子集的方法相比,采用本发明时,观测数据I和观测数据II模糊度正确固定率分别提高了13.3%和6.8%。说明采用本发明方法能更加可靠地选择一组模糊度子集,快速进行模糊度的正确分解,从而在城市环境下,实现利用GNSS进行快速高精度定位。
图3是观测数据I和观测数据II采用本发明方法和根据卫星高度角方法在选择模糊度子集的过程中剔除的卫星个数统计信息。从图中可以看出,本发明方法剔除1~2颗卫星的历元约占历元总数的95.9%,优于根据卫星高度角选择模糊度子集方法的86.6%。说明本发明方法能更加有效地剔除低精度的模糊度,避免其对模糊度分解造成的影响。
图4和图5分别是观测数据I和观测数据II采用本发明时第一天观测数据的动态定位结果。从图中可以看出,观测数据I和观测数据II在North和East方向动态定位结果达到0.04m内,在Up方向的定位结果在0.1m内。基线定位结果达到了厘米级的定位精度,满足GNSS高精度动态定位的精度要求。
本发明综合考虑了模糊度残差及模糊度间协方差的影响,以模糊度分解的Ratio值及浮点模糊度与最优整周模糊度间的Mahalanobis距离作为选择模糊度子集的准则。由于综合考虑了各个因素的影响,该准则能快速地选择一组可靠的模糊度子集,有效避免偏差较大的浮点模糊度对模糊度分解造成的影响,实现模糊度的快速固定,确保GNSS快速高精度动态定位结果的正确性和可靠性。
Claims (6)
1.一种复杂环境条件下多星座GNSS融合高精度动态定位方法,其特征在于,包括以下步骤:
步骤1:获取浮点模糊度解集合及其方差-协方差阵对和进行分解得到整周模糊度值,若整周模糊度值满足设定要求则转入步骤4,否则转入步骤2;
步骤2:从中剔除第i颗卫星的模糊度,然后进行模糊度分解,得到Ratio值ratioi和浮点模糊度与最优整周模糊度间的距离sni;
步骤3:将步骤2得到的ratioi和sni根据设定阈值分组后确定最优模糊度子集,根据最优模糊度子集进行整周模糊度解算得到整周模糊度值;
步骤4:根据整周模糊度值进行GNSS固定解的解算,得到GNSS动态定位结果。
2.根据权利要求1所述的一种复杂环境条件下多星座GNSS融合高精度动态定位方法,其特征在于,所述步骤3中确定最优模糊度子集的过程如下:
将ratioi大于设定阈值的模糊度子集作为候选模糊度子集;
若候选模糊度子集个数不为零,则选择sni值最小的候选模糊度子集作为最优模糊度子集;
若候选模糊度子集个数为零,则将ratioi中最大值对应的浮点模糊度子集取代重复步骤2~3;直至剩余卫星个数小于设定值为止。
3.根据权利要求1所述的一种复杂环境条件下多星座GNSS融合高精度动态定位方法,其特征在于,所述步骤2中ratioi的计算过程如下:
式中:为剔除第i颗卫星的模糊度后剩余卫星的浮点模糊度解集合;为剔除第i颗卫星的模糊度后剩余卫星的方差-协方差阵,为的逆矩阵,为模糊度分解后获得的整周模糊度最优值,为分解后获得的整周模糊度次优值,q为卫星数。
4.根据权利要求1所述的一种复杂环境条件下多星座GNSS融合高精度动态定位方法,其特征在于,所述步骤2中sni的计算过程如下:
式中:为剔除第i颗卫星的模糊度后剩余卫星的浮点模糊度解集合;为剔除第i颗卫星的模糊度后剩余卫星的方差-协方差阵,为的逆矩阵,为模糊度分解后获得的整周模糊度最优值,q为卫星数。
5.根据权利要求1所述的一种复杂环境条件下多星座GNSS融合高精度动态定位方法,其特征在于,所述步骤1中通过卡尔曼滤波方法对观测数据进行处理得到浮点模糊度解集合及其方差-协方差阵
6.根据权利要求1所述的一种复杂环境条件下多星座GNSS融合高精度动态定位方法,其特征在于,所述步骤4中通过将固定整周模糊度值代入观测方程,利用卡尔曼滤波方法进行GNSS固定解的解算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810824778.2A CN109116392B (zh) | 2018-07-25 | 2018-07-25 | 复杂环境条件下多星座gnss融合高精度动态定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810824778.2A CN109116392B (zh) | 2018-07-25 | 2018-07-25 | 复杂环境条件下多星座gnss融合高精度动态定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109116392A true CN109116392A (zh) | 2019-01-01 |
CN109116392B CN109116392B (zh) | 2023-01-10 |
Family
ID=64863363
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810824778.2A Active CN109116392B (zh) | 2018-07-25 | 2018-07-25 | 复杂环境条件下多星座gnss融合高精度动态定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109116392B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110376629A (zh) * | 2019-07-16 | 2019-10-25 | 东南大学 | 基于ratio值极大原则的卫星差分系统间偏差确定方法 |
CN112083464A (zh) * | 2019-06-14 | 2020-12-15 | 北京合众思壮科技股份有限公司 | 一种部分模糊度的固定方法及装置 |
CN112230261A (zh) * | 2020-09-30 | 2021-01-15 | 桂林电子科技大学 | 一种部分模糊度固定方法、装置、计算机设备和存储介质 |
CN115480278A (zh) * | 2022-08-01 | 2022-12-16 | 北方雷科(安徽)科技有限公司 | 一种适应于复杂多径环境下的双天线定向抗差算法 |
CN117452462A (zh) * | 2023-12-22 | 2024-01-26 | 华芯拓远(天津)科技有限公司 | 一种模型和数据相结合的部分模糊度固定方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010021656A2 (en) * | 2008-08-19 | 2010-02-25 | Trimble Navigation Limited | Gnss signal processing methods and apparatus with ambiguity selection |
CN103576175A (zh) * | 2013-11-06 | 2014-02-12 | 西南交通大学 | 一种双频多星座gnss整周模糊度otf解算方法 |
CN104459745A (zh) * | 2014-12-25 | 2015-03-25 | 东南大学 | 一种多星座长基线网络rtk部分模糊度快速解算方法 |
-
2018
- 2018-07-25 CN CN201810824778.2A patent/CN109116392B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010021656A2 (en) * | 2008-08-19 | 2010-02-25 | Trimble Navigation Limited | Gnss signal processing methods and apparatus with ambiguity selection |
CN103576175A (zh) * | 2013-11-06 | 2014-02-12 | 西南交通大学 | 一种双频多星座gnss整周模糊度otf解算方法 |
CN104459745A (zh) * | 2014-12-25 | 2015-03-25 | 东南大学 | 一种多星座长基线网络rtk部分模糊度快速解算方法 |
Non-Patent Citations (5)
Title |
---|
冯威等: "FirCAR算法的OTF快速定位方法", 《测绘学报》 * |
舒宝等: "基于BDS/GPS组合定位的部分模糊度固定效果分析", 《武汉大学学报(信息科学版)》 * |
赵丽等: "基于QR分解的快速解算初始整周模糊度方法的研究", 《南京航空航天大学学报》 * |
赵兴旺等: "精密单点定位部分模糊度固定方法", 《中国空间科学技术》 * |
邵晓东等: "分步组合部分模糊度固定法在中短基线中的应用", 《测绘科学》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112083464A (zh) * | 2019-06-14 | 2020-12-15 | 北京合众思壮科技股份有限公司 | 一种部分模糊度的固定方法及装置 |
CN112083464B (zh) * | 2019-06-14 | 2023-12-26 | 北京合众思壮科技股份有限公司 | 一种部分模糊度的固定方法及装置 |
CN110376629A (zh) * | 2019-07-16 | 2019-10-25 | 东南大学 | 基于ratio值极大原则的卫星差分系统间偏差确定方法 |
CN110376629B (zh) * | 2019-07-16 | 2021-03-19 | 东南大学 | 基于ratio值最大原则的卫星差分系统间偏差确定方法 |
CN112230261A (zh) * | 2020-09-30 | 2021-01-15 | 桂林电子科技大学 | 一种部分模糊度固定方法、装置、计算机设备和存储介质 |
CN112230261B (zh) * | 2020-09-30 | 2024-05-10 | 桂林电子科技大学 | 一种部分模糊度固定方法、装置、计算机设备和存储介质 |
CN115480278A (zh) * | 2022-08-01 | 2022-12-16 | 北方雷科(安徽)科技有限公司 | 一种适应于复杂多径环境下的双天线定向抗差算法 |
CN117452462A (zh) * | 2023-12-22 | 2024-01-26 | 华芯拓远(天津)科技有限公司 | 一种模型和数据相结合的部分模糊度固定方法 |
CN117452462B (zh) * | 2023-12-22 | 2024-03-29 | 华芯拓远(天津)科技有限公司 | 一种模型和数据相结合的部分模糊度固定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109116392B (zh) | 2023-01-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109116392A (zh) | 复杂环境条件下多星座gnss融合高精度动态定位方法 | |
CN106168672B (zh) | 一种gnss多模单频rtk周跳探测方法及装置 | |
CN106873017B (zh) | 模糊度固定方法 | |
CN112230261B (zh) | 一种部分模糊度固定方法、装置、计算机设备和存储介质 | |
Groves et al. | Combining inertially-aided extended coherent integration (supercorrelation) with 3D-mapping-aided GNSS | |
CN108919321A (zh) | 一种基于尝试法的gnss定位粗差探测方法 | |
CN102998681A (zh) | 一种卫星导航系统的高频钟差估计方法 | |
CN108254773A (zh) | 一种多gnss的实时钟差解算方法 | |
CN112146663B (zh) | 基于偏置星载加速度计的gnss卫星精密定轨增强方法及系统 | |
CN104076381B (zh) | 实时精密单点定位方法 | |
CN109633703B (zh) | 一种应对遮挡场景的北斗导航无源定位方法 | |
Seepersad | Reduction of initial convergence period in GPS PPP data processing | |
CN112130174B (zh) | 一种改进的gnss-ir提取积雪深度方法 | |
Dai et al. | Innovative algorithms to improve long range RTK reliability and availability | |
CN112230254B (zh) | 一种gps载波相位多径误差的校正方法及装置 | |
CN105785416A (zh) | 基线约束下单频单历元gnss快速定向方法 | |
CN108254774A (zh) | 基于gnss多频信号的单基站长距离实时定位方法 | |
CN113419262A (zh) | 一种全系统rtk快速选星方法 | |
WO2023236643A1 (zh) | 定位方法、装置、设备及存储介质 | |
CN109738912B (zh) | 一种基于gnss卫星信号实现固定点授时的方法 | |
CN115586555A (zh) | 一种基于ost与adop耦合的gnss部分模糊度固定方法 | |
Zhang et al. | Sparse or dense: Challenges of Australian network RTK | |
CN115616637A (zh) | 一种基于三维格网多径建模的城市复杂环境导航定位方法 | |
CN103837130B (zh) | 用于机载激光扫描系统的数据处理方法及装置 | |
CN108957487A (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 |