CN109490925B - 基于双频组合的室内定位方法 - Google Patents

基于双频组合的室内定位方法 Download PDF

Info

Publication number
CN109490925B
CN109490925B CN201811439234.0A CN201811439234A CN109490925B CN 109490925 B CN109490925 B CN 109490925B CN 201811439234 A CN201811439234 A CN 201811439234A CN 109490925 B CN109490925 B CN 109490925B
Authority
CN
China
Prior art keywords
frequency
dual
double
path
phase
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
CN201811439234.0A
Other languages
English (en)
Other versions
CN109490925A (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 CN201811439234.0A priority Critical patent/CN109490925B/zh
Publication of CN109490925A publication Critical patent/CN109490925A/zh
Application granted granted Critical
Publication of CN109490925B publication Critical patent/CN109490925B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/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/42Determining position
    • G01S19/43Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
    • 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

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

本发明提出一种基于双频组合的室内定位方法,主要步骤为发射基站生成多路双频导航信号,并使两路信号相对1pps的偏移量相同;然后计算各频点的实测距离值或组合相位实测距离差值;最后利用最小二乘定位算法进行定位;本发明通过利用双频信号的方式构建观测量,有效扩大室内定位的作用范围,并通过模拟验证可得由此得出的定位精度均能满足1米需求。

Description

基于双频组合的室内定位方法
技术领域
本发明适用范围很广,特别是针对当前的室内外无缝融合定位领域。
背景技术
全球卫星导航系统(global navigation satellite system,GNSS)已广泛应用于人们的日常生活,在海陆空交通,时间同步网络,生命安全,位置服务等国民经济生活方面可谓无处不在,成为公认的现代社会基础设施之一。但是在隧道、城市峡谷、室内等区域,由于信号遮挡无法完成定位服务。当前的伪卫星室内定位系统由于种种缺陷,室内定位精度很难满足定位需求,本文在已有的室内多阵元伪卫星基础上提出基于双频组合的室内定位方法,保障室内定位精度可以达到1米。
发明内容
基于上述原因,本发明解决的技术问题是通过利用双频组合的定位方法有效的解决了室内伪卫星定位精度差难以使用的问题,同时由于与GNSS信号兼容,有效的实现了室内外无缝导航服务能力。
本发明采用的技术方案为:
基于双频组合的室内定位方法,包括以下步骤:
(1)发射基站在同一时钟源的控制下生成多路双频导航信号;
(2)发射基站将每一路双频导航信号分别通过开关滤波器,依次滤掉双频导航信号的其中1路,计算另一路信号与1pps的偏移量,根据偏移量调整零相位点与1PPS上升沿对齐,使两路信号相对1pps的偏移量相同;
(3)接收站接收多路双频导航信号,将每一路双频导航信号通过比较双频小数相位比例关系,计算出各频点信号的载波相位整周数,然后加上周内载波相位值得到各频点的实测距离值;或将每一路双频导航信号进行宽巷组合,得到发射基站每路双频导航信号的双频组合相位值,利用得到的发射基站各路双频导航信号的双频组合相位值进行单差,得到一组组合相位实测距离差值。
(4)根据多组实测距离值或组合相位实测距离差值,利用最小二乘定位算法进行定位;
完成基于双频组合的室内定位。
其中,步骤(3)中每一路双频导航信号通过比较双频小数相位比例关系,计算出各频点信号的载波相位整周数,然后加上周内载波相位值得到各频点的实测距离值,包括以下步骤:
(301)计算每一路双频导航信号的波长差DeltaL
Figure BDA0001884232450000021
其中,DeltaL为双频导航信号的波长差,λL1为双频导航信号中L1频点的波长,λB1为双频导航信号中B1频点的波长;
(302)计算双频实测相位差DeltaPh:
PhL1=CarrL1-floor(CarrL1)
PhB1=CarrB1-floor(CarrB1)
DeltaPh=PhB1-PhL1
其中,CarrL1、CarrB1分别为双频导航信号中频点L1和频点B1的载波相位,floor(·)为向下取整函数,PhL1、PhB1分别为双频导航信号中频点L1和频点B1的周内载波相位,DeltaPh为双频实测相位差值;
(303)计算实测距离中载波相位整周数N
Figure BDA0001884232450000022
(304)计算频点L1和频点B1的实测距离值
Figure BDA0001884232450000023
Figure BDA0001884232450000024
Figure BDA0001884232450000025
其中,接收站距离发射基站的距离满足下式:
Figure BDA0001884232450000026
其中,L为发射基站的发射阵元到最远覆盖范围边界点的距离,λL1为双频导航信号中频点L1的波长,λL2为双频导航信号中频点L2的波长,floor(·)为向下取整函数。
其中,步骤(3)中测算组合相位实测距离差的方式要求各路双频导航信号对应的发射天线阵元间距为组合波长的半波长。
其中,步骤(3)中将每一路双频导航信号进行宽巷组合,得到发射基站每路双频导航信号的双频组合相位值,利用得到的发射基站各路双频导航信号的双频组合相位值进行单差,得到一组组合相位实测距离差值具体计算方式为:
每路双频导航信号中L1、B1两个频点的双频组合相位值为:
Figure BDA0001884232450000031
其中,φL1,φB1分别为L1、B1两个频点的双频组合相位值,λL1、λB1分别为L1、B1两个频点的波长,IL1、IB1为L1、B1两个频点的电离层修正值、T为对流层修正值、g为发射钟差,NL1、NB1为L1、B1两个频点的整周模糊度,εφ,L1,εφ,B1为L1、B1两个频点的噪声误差,r为天线阵元到接收机的的实际距离;
将双频组合相位值做差得到组合相位实测距离差值φw
φw=φL1B1
即:
Figure BDA0001884232450000032
其中,步骤(4)利用组合相位实测距离差值进行定位过程如下:
Figure BDA0001884232450000033
将整周模糊度进行抵消,则组合相位实测距离差值
Figure BDA0001884232450000034
变为
Figure BDA0001884232450000035
则n路双频导航信号时,由上式可得定位方程如下
Figure BDA0001884232450000041
以第一路观测方程为基准,做差得
Figure BDA0001884232450000042
由于多路信号同源发射,则可得g21、g31···gn1为0,将上式转换成距离形式可得,
Figure BDA0001884232450000043
r21=||r2-ru||-||r1-ru||
r31=||r3-ru||-||r1-ru||
rn1=||rn-ru||-||r1-ru||
其中rn为第n路发射天线位置,ru为接收机位置;
上式转换为距离形式可得
Figure BDA0001884232450000044
Figure BDA0001884232450000045
Figure BDA0001884232450000051
最小二乘展开,化简完最后可得:
Figure BDA0001884232450000052
其中x0,y0,z0为初始坐标值,x,y,z为实际结算出的坐标值。
本发明相比现有技术有如下优点:
(1)本发明方法可用性强,避免了周整模糊度的解算过程。
(2)本发明方法测距与定位原理简单,只需用最小二乘即可实现1米精度定位。
(3)本发明成本低,只需借助成熟的导航芯片与软核相结合便可实现。
(4)本发明军民市场均可使用,航空管理、道路交通、室内指引、快速增援、战略引导、防灾减灾等。
附图说明
图1为本发明流程图。
具体实施方式
下面结合附图对本发明作进一步的解释说明。
本发明基于双频组合的室内定位方法,包括以下步骤:
(1)发射基站在同一时钟源的控制下生成多路双频导航信号;
(2)发射基站将每一路双频导航信号分别通过开关滤波器,依次滤掉双频导航信号的其中1路,计算另一路信号与1pps的偏移量,根据偏移量调整零相位点与1PPS上升沿对齐,使两路信号相对1pps的偏移量相同;
(3)接收站接收多路双频导航信号,将每一路双频导航信号通过比较双频小数相位比例关系,计算出各频点信号的载波相位整周数,然后加上周内载波相位值得到各频点的实测距离值;包括以下步骤:
(301)计算每一路双频导航信号的波长差DeltaL
Figure BDA0001884232450000061
其中,DeltaL为双频导航信号的波长差,λL1为双频导航信号中L1频点的波长,λB1为双频导航信号中B1频点的波长;
(302)计算双频实测相位差DeltaPh
PhL1=CarrL1-floor(CarrL1)
PhB1=CarrB1-floor(CarrB1)
DeltaPh=PhB1-PhL1
其中,CarrL1、CarrB1分别为双频导航信号中频点L1和频点B1的载波相位,floor(·)为向下取整函数,PhL1、PhB1分别为双频导航信号中频点L1和频点B1的周内载波相位,DeltaPh为双频实测相位差值;
(303)计算实测距离中载波相位整周数N
Figure BDA0001884232450000062
(304)计算频点L1和频点B1的实测距离值
Figure BDA0001884232450000063
Figure BDA0001884232450000064
Figure BDA0001884232450000065
Figure BDA0001884232450000066
或将每一路双频导航信号进行宽巷组合,得到发射基站每路双频导航信号的双频组合相位值,利用得到的发射基站各路双频导航信号的双频组合相位值进行单差,得到一组组合相位实测距离差值。具体计算方式为:
每路双频导航信号中L1、B1两个频点的双频组合相位值为
Figure BDA0001884232450000071
其中,φL1,φB1分别为L1、B1两个频点的双频组合相位值,λL1、λB1分别为L1、B1两个频点的波长,IL1、IB1为L1、B1两个频点的电离层修正值、T为对流层修正值、g为发射钟差,NL1、NB1为L1、B1两个频点的整周模糊度,εφ,L1,εφ,B1为L1、B1两个频点的噪声误差,r为天线阵元到接收机的的实际距离;
将双频组合相位值做差得到组合相位实测距离差值φw
φw=φL1B1
即:
Figure BDA0001884232450000072
(4)根据多组实测距离值或组合相位实测距离差值,利用最小二乘定位算法进行定位;利用组合相位实测距离差值进行定位过程如下:
Figure BDA0001884232450000073
将整周模糊度进行抵消,则组合相位实测距离差值
Figure BDA0001884232450000074
变为
Figure BDA0001884232450000075
则n路双频导航信号时,由上式可得定位方程如下
Figure BDA0001884232450000076
以第一路观测方程为基准,做差得
Figure BDA0001884232450000081
由于多路信号同源发射,则可得g21、g31···gn1为0,将上式转换成距离形式可得,
Figure BDA0001884232450000082
r21=||r2-ru||-||r1-ru||
r31=||r3-ru||-||r1-ru||
rn1=||rn-ru||-||r1-ru||
其中rn为第n路发射天线位置,ru为接收机位置;
上式转换为距离形式可得
Figure BDA0001884232450000083
Figure BDA0001884232450000084
Figure BDA0001884232450000085
最小二乘展开,化简完最后可得:
Figure BDA0001884232450000091
其中x0,y0,z0为初始坐标值,x,y,z为实际结算出的坐标值。
完成基于双频组合的室内定位。
以上所述,仅为本发明的一具体实例,但本发明的保护范围并不仅限于此,在本发明揭露的技术范围内,可理解想到的变换,都应涵盖在本发明的包含范围内。

Claims (6)

1.基于双频组合的室内定位方法,其特征在于,包括以下步骤:
(1)发射基站在同一时钟源的控制下生成多路双频导航信号;
(2)发射基站将每一路双频导航信号分别通过开关滤波器,依次滤掉双频导航信号的其中1路,计算另一路信号与1pps的偏移量,根据偏移量调整零相位点与1PPS上升沿对齐,使两路信号相对1pps的偏移量相同;
(3)接收站接收多路双频导航信号,将每一路双频导航信号通过比较双频小数相位比例关系,计算出各频点信号的载波相位整周数,然后加上周内载波相位值得到各频点的实测距离值;或将每一路双频导航信号进行宽巷组合,得到发射基站每路双频导航信号的双频组合相位值,利用得到的发射基站各路双频导航信号的双频组合相位值进行单差,得到一组组合相位实测距离差值;
(4)根据多组实测距离值或组合相位实测距离差值,利用最小二乘定位算法进行定位;
完成基于双频组合的室内定位。
2.根据权利要求1所述的基于双频组合的室内定位方法,其特征在于:
步骤(3)中每一路双频导航信号通过比较双频小数相位比例关系,计算出各频点信号的载波相位整周数,然后加上周内载波相位值得到各频点的实测距离值,包括以下步骤:
(301)计算每一路双频导航信号的波长差DeltaL
Figure FDA0002434854840000011
其中,DeltaL为双频导航信号的波长差,λL1为双频导航信号中L1频点的波长,λB1为双频导航信号中B1频点的波长;
(302)计算双频实测相位差DeltaPh
PhL1=CarrL1-floor(CarrL1)
PhB1=CarrB1-floor(CarrB1)
DeltaPh=PhB1-PhL1
其中,CarrL1、CarrB1分别为双频导航信号中频点L1和频点B1的载波相位,floor(·)为向下取整函数,PhL1、PhB1分别为双频导航信号中频点L1和频点B1的周内载波相位,DeltaPh为双频实测相位差值;
(303)计算实测距离中载波相位整周数N
Figure FDA0002434854840000021
(304)计算频点L1和频点B1的实测距离值
Figure FDA0002434854840000022
Figure FDA0002434854840000023
Figure FDA0002434854840000024
Figure FDA0002434854840000025
3.根据权利要求1中所述的基于双频组合的室内定位方法,其特征在于:接收站距离发射基站的距离满足下式:
Figure FDA0002434854840000026
其中,L为发射基站的发射阵元到最远覆盖范围边界点的距离,λL1为双频导航信号中频点L1的波长,λL2为双频导航信号中频点L2的波长,floor(·)为向下取整函数。
4.根据权利要求1所述的基于双频组合的室内定位方法,其特征在于:步骤(3)中测算组合相位实测距离差的方式要求各路双频导航信号对应的发射天线阵元间距为组合波长的半波长。
5.根据权利要求1所述的基于双频组合的室内定位方法,其特征在于:步骤(3)中将每一路双频导航信号进行宽巷组合,得到发射基站每路双频导航信号的双频组合相位值,利用得到的发射基站各路双频导航信号的双频组合相位值进行单差,得到一组组合相位实测距离差值具体计算方式为:
每路双频导航信号中L1、B1两个频点的双频组合相位值为:
Figure FDA0002434854840000027
其中,φL1,φB1分别为L1、B1两个频点的双频组合相位值,λL1、λB1分别为L1、B1两个频点的波长,IL1、IB1为L1、B1两个频点的电离层修正值、T为对流层修正值、g为发射钟差,NL1、NB1为L1、B1两个频点的整周模糊度,εφ,L1,εφ,B1为L1、B1两个频点的噪声误差,r为天线阵元到接收机的的实际距离;
将双频组合相位值做差得到组合相位实测距离差值φw
φw=φL1B1
即:
Figure FDA0002434854840000031
6.根据权利要求5所述的基于双频组合的室内定位方法,其特征在于:步骤(4)利用组合相位实测距离差值进行定位过程如下:
Figure FDA0002434854840000032
将整周模糊度进行抵消,则组合相位实测距离差值
Figure FDA0002434854840000033
变为
Figure FDA0002434854840000034
则n路双频导航信号时,由上式可得定位方程如下:
Figure FDA0002434854840000035
以第一路观测方程为基准,做差得:
Figure FDA0002434854840000036
由于多路信号同源发射,则可得g21、g31···gn1为0,将上式转换成距离形式可得:
Figure FDA0002434854840000041
r21=||r2-ru||-||r1-ru||
r31=||r3-ru||-||r1-ru||
rn1=||rn-ru||-||r1-ru||
其中rn为第n路发射天线位置,ru为接收机位置;
上式转换为距离形式可得:
Figure FDA0002434854840000042
Figure FDA0002434854840000043
Figure FDA0002434854840000044
最小二乘展开,化简完最后可得:
Figure FDA0002434854840000051
其中x0,y0,z0为初始坐标值,x,y,z为实际结算出的坐标值。
CN201811439234.0A 2018-11-29 2018-11-29 基于双频组合的室内定位方法 Active CN109490925B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811439234.0A CN109490925B (zh) 2018-11-29 2018-11-29 基于双频组合的室内定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811439234.0A CN109490925B (zh) 2018-11-29 2018-11-29 基于双频组合的室内定位方法

Publications (2)

Publication Number Publication Date
CN109490925A CN109490925A (zh) 2019-03-19
CN109490925B true CN109490925B (zh) 2020-10-16

Family

ID=65698551

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811439234.0A Active CN109490925B (zh) 2018-11-29 2018-11-29 基于双频组合的室内定位方法

Country Status (1)

Country Link
CN (1) CN109490925B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110703287B (zh) * 2019-10-30 2021-08-31 中国电子科技集团公司第五十四研究所 一种基于阵列双频互备份的多星特征参数识别定位方法
CN111736138B (zh) * 2020-06-23 2023-05-02 西安电子科技大学 基于ofdm信号和三频载波相位测距的室内测距方法
CN112526562B (zh) * 2020-11-26 2022-05-06 中国电子科技集团公司第五十四研究所 一种基于l1/b1双频的pso粒子群测距定位方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5442363A (en) * 1994-08-04 1995-08-15 U.S. Army Corps Of Engineers As Represented By The Secretary Of The Army Kinematic global positioning system of an on-the-fly apparatus for centimeter-level positioning for static or moving applications
CN106646370A (zh) * 2017-02-20 2017-05-10 吕志勤 一种基于载波相位测量的可见光通讯室内定位装置及方法
CN108054522A (zh) * 2017-11-01 2018-05-18 深圳市时空导航科技有限公司 室内gnss天线阵列、定位系统、定位方法及装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104035068B (zh) * 2014-06-26 2016-09-14 桂林电子科技大学 一种基于伪卫星的室内定位系统及方法
CN104898129B (zh) * 2015-06-16 2017-03-29 桂林电子科技大学 通用gps室内外定位系统和方法
CN106595667B (zh) * 2016-12-07 2019-04-26 中国电子科技集团公司第五十四研究所 一种基于互同步的室内伪卫星双点定位系统及方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5442363A (en) * 1994-08-04 1995-08-15 U.S. Army Corps Of Engineers As Represented By The Secretary Of The Army Kinematic global positioning system of an on-the-fly apparatus for centimeter-level positioning for static or moving applications
CN106646370A (zh) * 2017-02-20 2017-05-10 吕志勤 一种基于载波相位测量的可见光通讯室内定位装置及方法
CN108054522A (zh) * 2017-11-01 2018-05-18 深圳市时空导航科技有限公司 室内gnss天线阵列、定位系统、定位方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Indoor positioning based on difference between carrier-phases transmitted from proximately-located antennas of a multi-channel pseudolite;Yoshihiro Sakamoto等;《2012 First International Conference on Innovative Engineering Systems》;20130613;155-156 *

Also Published As

Publication number Publication date
CN109490925A (zh) 2019-03-19

Similar Documents

Publication Publication Date Title
RU2565386C2 (ru) Способ, устройство и система для определения позиции объекта, имеющего приемник глобальной навигационной спутниковой системы, посредством обработки неразностных данных, подобных измерениям фазы несущей, и внешних данных, подобных ионосферным данным
CN105891860B (zh) 一种基于误差分离模式的gnss区域伪距差分增强定位方法
US7990314B2 (en) Method and system for locating a geographical position using broadcast frequency modulation signals
CN109490925B (zh) 基于双频组合的室内定位方法
CN109407126A (zh) 一种多模接收机联合定位解算的方法
CN103529459A (zh) 一种采用单频gps和glonass组合精准定位的方法及其系统
CN105182382A (zh) 一种伪卫星厘米级定位方法
KR101874974B1 (ko) Dgnss 의사거리 보정정보 생성장치 및 그 방법
CN105005019A (zh) 一种基于bds/gps rtk的地下管线测量方法
CN104316943A (zh) 一种伪距离和多普勒组合差分定位系统及方法
CN104730551A (zh) 一种星地双基地差分干涉基线坐标以及形变量测量方法
CN104898144A (zh) 一种北斗系统参考站整周模糊度单历元确定方法
CN110187376A (zh) 一种bds/gps同时钟源的伪卫星多普勒差分测速方法
CN102087363B (zh) 一种中轨道卫星搜救系统定位方法
CN109490828B (zh) 基于同源基线阵列的定位方法
CN105158781A (zh) 一种地球gnss卫星实现外层空间用户定位的方法
CN116359968B (zh) 一种联合北斗二号和北斗三号的三频差分定位方法
JP2000193733A (ja) 測位装置及び測位方法
CN109116396B (zh) 一种多天线gnss差分定位方法
RU2713571C1 (ru) Система высокоточной дифференциальной коррекции для контроля подвижных объектов
CN103869326B (zh) 一种基于伪距指纹匹配的区域快速定位方法
CN103675875B (zh) Bds与gps的b2/l2载波相位星间混频差分方法
CN113156477B (zh) 一种基于Android智能手机的高精度RTK定位方法
Innac et al. Multi-GNSS single frequency precise point positioning
Cao et al. Inter-frequency code bias handling and estimation for multi-frequency BeiDou-3/Galileo uncombined precise point positioning

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