CN113099405B - 一种基于发射角估计的音频定位方法 - Google Patents

一种基于发射角估计的音频定位方法 Download PDF

Info

Publication number
CN113099405B
CN113099405B CN202110263517.XA CN202110263517A CN113099405B CN 113099405 B CN113099405 B CN 113099405B CN 202110263517 A CN202110263517 A CN 202110263517A CN 113099405 B CN113099405 B CN 113099405B
Authority
CN
China
Prior art keywords
base station
angle
sound source
receiving end
sound
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
CN202110263517.XA
Other languages
English (en)
Other versions
CN113099405A (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.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
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 Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN202110263517.XA priority Critical patent/CN113099405B/zh
Publication of CN113099405A publication Critical patent/CN113099405A/zh
Application granted granted Critical
Publication of CN113099405B publication Critical patent/CN113099405B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W4/00Services specially adapted for wireless communication networks; Facilities therefor
    • H04W4/30Services specially adapted for particular environments, situations or purposes
    • H04W4/33Services specially adapted for particular environments, situations or purposes for indoor environments, e.g. buildings
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W64/00Locating users or terminals or network equipment for network management purposes, e.g. mobility management

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明属于室内定位技术领域,公开了一种基于发射角估计的音频定位方法,布设至少两个已知位置信息和姿态信息的基站,每个基站均配置有声源发射器阵列,声源发射器阵列包括若干个向不同方向发射声音信号的声源发射器;每个基站通过声源发射器不断发射带有识别码的声音信号,不同的声源发射器发射的声音信号的识别码唯一确定且相互不同;接收端接收来自所有基站的所有声音信号,并对声音信号进行解调得到接收端和基站之间的角度概率分布;根据角度概率分布得到接收端与每个基站对应的最佳空间估计角度;根据基站坐标系下的最佳空间估计角度、基站的位置信息和姿态信息解算得到接收端的空间位置。本发明解决了室内定位的精度较低的问题。

Description

一种基于发射角估计的音频定位方法
技术领域
本发明属于室内定位技术领域,更具体地,涉及一种基于发射角估计的音频定位方法。
背景技术
随着室内外定位技术的发展,基于广覆盖、多层次的位置服务需求巨大。目前,在室外场景下基于GNSS的定位已成熟发展,但在室内、城市峡谷地区、地下等复杂空间由于GNSS难以抵达而无法实现定位,因此已经发展了利用WIFI、蓝牙、超宽带等射频信号进行室内定位研究的技术。但是室内场景复杂多变,电磁场环境不稳定,对定位精度影响较大。此外,较为成熟的超宽带定位方案成本和功耗过高,难以形成普及应用,蓝牙定位方案虽然成本低,但作用距离有限,且在室内等GNSS拒止环境下难以满足高精度的位置服务需求。
发明内容
本发明通过提供一种基于发射角估计的音频定位方法,解决现有技术中室内定位的精度较低的问题。
本发明提供一种基于发射角估计的音频定位方法,包括以下步骤:
步骤1:布设至少两个已知位置信息和姿态信息的基站,每个所述基站均配置有声源发射器阵列,所述声源发射器阵列包括若干个向不同方向发射声音信号的声源发射器;
步骤2:每个所述基站通过所述声源发射器不断发射带有识别码的声音信号,不同的声源发射器发射的声音信号的识别码唯一确定且相互不同;
步骤3:接收端接收来自所有基站的所有声音信号;
步骤4:所述接收端对声音信号进行解调,得到识别码种类和该种类识别码声音信号的强度,结合识别码和声源发射器的编号之间的映射关系、声源发射器的编号和基站之间的归属关系、声源发射器的编号和声音信号的发射角度之间的映射关系,得到所述接收端和所述基站之间的角度概率分布;
步骤5:根据所述角度概率分布得到所述接收端与每个所述基站对应的最佳空间估计角度;
步骤6:根据所述基站坐标系下的最佳空间估计角度、所述基站的位置信息和姿态信息解算得到所述接收端的空间位置。
优选的,所述步骤2中,所述识别码包含该识别码分别与对应的声源发射器、基站之间的所属关系信息;
所述识别码包含该识别码所属声源发射器的编号与该识别码所属声源发射器发射的声音信号的波束方向角度之间的映射关系信息。
优选的,所述步骤3中,所述接收端通过麦克风接收所述声音信号。
优选的,所述步骤4中,得到所述接收端和所述基站之间的角度概率分布的具体实现方式为:
结合已知的识别码和声源发射器的编号之间的映射关系、声源发射器的编号和基站之间的归属关系,对不同的识别码的声音信号强度和基站归属建立的拓扑关系;针对每一个基站,对该基站的所有信号强度分别进行归一化处理,归一化处理后得到的值视为所述接收端位于该基站的角度概率分布。
优选的,所述步骤5中,对每个基站对应的所述角度概率分布进行处理,通过加权求和的处理方法得到所述接收端与该基站对应的最佳空间估计角度。
优选的,所述步骤5中,对每个基站对应的所述角度概率分布进行处理,将最大角度概率对应的波束方向角度作为所述接收端与该基站对应的最佳空间估计角度。
优选的,所述步骤6的具体实现方式为:
所述最佳空间估计角度作为相对角度,针对每个所述基站,结合接收端与该基站对应的最佳空间估计角度以及该基站的所述姿态信息,得到所述接收端和该基站之间的绝对角度;得到所有基站与所述接收端之间的绝对角度后,结合所有基站的所述位置信息,解算得到所述接收端的空间位置。
本发明中提供的一个或多个技术方案,至少具有如下技术效果或优点:
在发明中,首先布设至少两个已知位置信息和姿态信息的基站,每个基站均配置有声源发射器阵列,声源发射器阵列包括若干个向不同方向发射声音信号的声源发射器;然后每个基站通过声源发射器不断发射带有识别码的声音信号,不同的声源发射器发射的声音信号的识别码唯一确定且相互不同;接收端接收来自所有基站的所有声音信号并对声音信号进行解调,得到接收端和基站之间的角度概率分布;之后根据角度概率分布得到接收端与每个基站对应的最佳空间估计角度;最后根据基站坐标系下的最佳空间估计角度、基站的位置信息和姿态信息解算得到接收端的空间位置。本发明基于声音信号进行定位,能够有效避免复杂电磁环境的干扰,具有较高的定位精度;基于声音信号的波长特性,使得信号作用距离比一般射频信号要远,能够实现室内场景的广覆盖。相比于信号接收端的到达角(AOA)估计方法,本发明利用声源发射器能够形成定向信号波束的优势,能够有效减少移动终端安置多天线的成本,广泛适用于能够接收声音信号的移动终端。相比于基于距离估计的方法,本发明对基站的时间同步要求低,能够有效降低基站侧的成本。
附图说明
图1为本发明实施例提供的一种基于发射角估计的音频定位方法的流程图;
图2为本发明实施例提供的一种基于发射角估计的音频定位方法中声源发射器和基站之间的归属关系、声源发射器和信号发射角度之间的映射关系、声源发射器和识别码之间的映射关系示意图。
具体实施方式
由于声音传输速度远低于电磁波信号,易于实现同步,具有低能耗、低成本的优势;声学特征不易受干扰,作用距离可达数十米,能够有效满足室内场景的信号覆盖。此外,基于单声源发射器能够形成定向信号波束,利用该优势,本发明提供一种基于角度估计的音频定位方法,能够实现低成本、高精度、广覆盖的室内定位。
为了更好的理解上述技术方案,下面将结合说明书附图以及具体的实施方式对上述技术方案进行详细的说明。
本实施例提供了一种基于发射角估计的音频定位方法,包括以下步骤:
步骤1:布设至少两个已知位置信息和姿态信息的基站,每个所述基站均配置有声源发射器阵列,所述声源发射器阵列包括若干个向不同方向发射声音信号的声源发射器;
步骤2:每个所述基站通过所述声源发射器不断发射带有识别码的声音信号,不同的声源发射器发射的声音信号的识别码唯一确定且相互不同;
步骤3:接收端接收来自所有基站的所有声音信号;
步骤4:所述接收端对声音信号进行解调,得到识别码种类和该种类识别码声音信号的强度,结合识别码和声源发射器的编号之间的映射关系、声源发射器的编号和基站之间的归属关系、声源发射器的编号和声音信号的发射角度之间的映射关系,得到所述接收端和所述基站之间的角度概率分布;
步骤5:根据所述角度概率分布得到所述接收端与每个所述基站对应的最佳空间估计角度;
步骤6:根据所述基站坐标系下的最佳空间估计角度、所述基站的位置信息和姿态信息解算得到所述接收端的空间位置。
本发明基于声音(机械波)信号进行定位,能够有效避免复杂电磁环境的干扰,具有较高的定位精度,且基于声音信号的波长特性,使得信号作用距离比一般射频信号要远,能够实现室内场景广覆盖。利用单声源发射器能够形成定向信号波束的优势,本发明是一种基于信号发射角(AOD)估计的声音定位方法,相比于信号接收端的到达角(AOA)估计方法,本发明能够有效减少移动终端安置多天线的成本,广泛适用于能够接收声音信号的移动终端。此外,基于角度测量的方法相比于基于距离估计的方法对基站的时间同步要求低,降低基站侧的成本。因此,本发明在有效提高室内定位精度的同时,能够实现低成本和广覆盖。
下面对本发明做进一步的说明。
本发明提供的一种基于发射角估计的音频定位方法,参见图1,包括以下步骤:
步骤1:布设至少两个已知位置信息和姿态信息的基站,每个基站均配置有声源发射器阵列,声源发射器阵列包括若干个向不同方向发射声音信号的声源发射器。
其中,声源发射器和基站之间的归属关系已知;声源发射器对应的发射声音信号的方向角度已知;所有的基站位置信息,姿态信息已知并被存储。已知参考信息参见本发明实例,如图1所示,具体描述如下:
1)布设两个基站,分别记为基站M和基站N,其中,编号为1至5的声源发射器归属基站M、编号为6至10的声源发射器归属基站N,该归属关系已知并存储与接收端本地或云端数据库;
2)编号为1至5的声源发射器对应基站M的θ1至θ5方向,编号为6至10的声源发射器对应基站N的θ6至θ10方向,该映射关系已知并存储与接收端本地或云端数据库;
3)基站M和基站N的位置信息(xM,yM,zM)和(xN,yN,zN)以及姿态信息(αM,βM,γM)和(αN,βN,γN)已知并被存储于接收端本地或云端数据库。
步骤2:每个基站通过声源发射器不断发射带有识别码的声音信号,整个系统中不同编号的声源发射器发射的声音信号的识别码唯一确定且相互不同,识别码和声源发射器的编号之间的映射关系已知。已知参考信息参见本发明实例,参见图1、图2,具体描述如下:
1)基站M和基站N上的声源发射器(编号为1至10)发射的声音信号所带识别码均相互不同(识别码a至j),而识别码a至j归属编号为1至10的声源发射器的映射关系均已知,且被存储于接收端本地或云端数据库。
步骤3:接收端接收来自所有基站的所有声音信号。
即接收端对声音信号进行接收采集,具体的,接收端通过麦克风即可完成声音信号的接收采集。
步骤4:接收端对声音信号进行解调,得到识别码种类和该种类识别码声音信号的强度,结合识别码和声源发射器的编号之间的映射关系、声源发射器的编号和基站之间的归属关系、声源发射器的编号和声音信号的发射角度(即声音信号的波束方向角度)之间的映射关系,经过处理得到接收端和基站之间角度概率分布。
即接收端对波束方向信息进行探测。具体包括:
(1)对步骤3中获取的声音信号进行解调,得到识别码种类和不同识别码声音信号的强度。
(2)结合步骤1和步骤2中已知的识别码和声源发射器的编号之间的映射关系、声源发射器的编号和基站之间的归属关系,对不同的识别码的声音信号强度和基站归属建立的拓扑关系,再对来自同一基站的所有信号强度分别进行归一化处理,处理后的值视为接收端位于该基站的角度概率分布。
该步骤可进一步通过发明实例进行说明,如图1所示,具体描述如下:
接收端对接收信号解调,得到不同识别码种类(编号从a至j)及其信号强度,再结合步骤1、步骤2中已知的归属和映射关系结合步骤1和步骤2中已知的识别码和声源发射器的编号之间的映射关系、声源发射器的编号和基站之间的归属关系,得到接收端位于基站M的不同方向(θ1至θ5)的概率分布和接收端位于基站N的不同方向(θ6至θ10)的信号强度,再对信号强度归一化处理,得到接收端位于该基站的角度概率分布。
步骤5:根据角度概率分布得到接收端与每个基站对应的最佳空间估计角度。
具体的,将不同波束方向角度对应的角度概率进行处理,通过加权求和、取最大值或插值等方法得到接收端与某个基站对应的最佳估计空间角度。例如,参见图1,分别根据基站M和基站N角度概率分布进行处理,得到接收端位于基站M和基站N的、在基站坐标系下的最佳估计角度θM和θN
步骤6:根据接收端基站坐标系下的最佳空间估计角度、基站的位置信息和姿态信息解算得到接收端的空间位置。
在基站坐标系下的最佳估计角度是相对角度,结合已知的基站的姿态信息,可以得到接收端和基站之间的绝对角度,当得到多个基站与接收端的绝对角度信息后,再结合已知的基站的位置信息,解算得到接收端的空间位置。该步骤可进一步通过发明实例进行说明,如图1所示,具体描述如下:
角度θM和θN为接收端相对基站坐标系的相对角度,结合步骤1中已知的基站姿态信息(αM,βM,γM)和(αN,βN,γN),可以得到接收端和基站M、基站N之间的绝对角度关系,再结合步骤1中已知的基站位置信息(xM,yM,zM)和(xN,yN,zN),通过解三角形等办法得到接收端的空间位置。
本发明实施例提供的一种基于发射角估计的音频定位方法至少包括如下技术效果:
(1)本发明基于声音信号进行定位,能够有效避免复杂电磁环境的干扰,具有较高的定位精度。
(2)基于声音信号的波长特性,使得信号作用距离比一般射频信号要远,能够实现室内场景的广覆盖。
(3)本发明提供的是一种基于信号发射角(AOD)估计的声音定位方法,相比于信号接收端的到达角(AOA)估计方法,本发明利用声源发射器能够形成定向信号波束的优势,能够有效减少移动终端安置多天线的成本,广泛适用于能够接收声音信号的移动终端。
(4)相比于基于距离估计的方法,本发明对基站的时间同步要求低,能够降低基站侧的成本。
最后所应说明的是,以上具体实施方式仅用以说明本发明的技术方案而非限制,尽管参照实例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (6)

1.一种基于发射角估计的音频定位方法,其特征在于,包括以下步骤:
步骤1:布设至少两个已知位置信息和姿态信息的基站,每个所述基站均配置有声源发射器阵列,所述声源发射器阵列包括若干个向不同方向发射声音信号的声源发射器;
步骤2:每个所述基站通过所述声源发射器不断发射带有识别码的声音信号,不同的声源发射器发射的声音信号的识别码唯一确定且相互不同;
步骤3:接收端接收来自所有基站的所有声音信号;
步骤4:所述接收端对声音信号进行解调,得到识别码种类和该种类识别码声音信号的强度,结合识别码和声源发射器的编号之间的映射关系、声源发射器的编号和基站之间的归属关系、声源发射器的编号和声音信号的发射角度之间的映射关系,得到所述接收端和所述基站之间的角度概率分布;
得到所述接收端和所述基站之间的角度概率分布的具体实现方式为:结合已知的识别码和声源发射器的编号之间的映射关系、声源发射器的编号和基站之间的归属关系,对不同的识别码的声音信号强度和基站归属建立的拓扑关系;针对每一个基站,对该基站的所有信号强度分别进行归一化处理,归一化处理后得到的值视为所述接收端位于该基站的角度概率分布;
步骤5:根据所述角度概率分布得到所述接收端与每个所述基站对应的最佳空间估计角度;
步骤6:根据所述基站坐标系下的最佳空间估计角度、所述基站的位置信息和姿态信息解算得到所述接收端的空间位置。
2.根据权利要求1所述的基于发射角估计的音频定位方法,其特征在于,所述步骤2中,所述识别码包含该识别码分别与对应的声源发射器、基站之间的所属关系信息;
所述识别码包含该识别码所属声源发射器的编号与该识别码所属声源发射器发射的声音信号的波束方向角度之间的映射关系信息。
3.根据权利要求1所述的基于发射角估计的音频定位方法,其特征在于,所述步骤3中,所述接收端通过麦克风接收所述声音信号。
4.根据权利要求1所述的基于发射角估计的音频定位方法,其特征在于,所述步骤5中,对每个基站对应的所述角度概率分布进行处理,通过加权求和的处理方法得到所述接收端与该基站对应的最佳空间估计角度。
5.根据权利要求1所述的基于发射角估计的音频定位方法,其特征在于,所述步骤5中,对每个基站对应的所述角度概率分布进行处理,将最大角度概率对应的波束方向角度作为所述接收端与该基站对应的最佳空间估计角度。
6.根据权利要求1所述的基于发射角估计的音频定位方法,其特征在于,所述步骤6的具体实现方式为:
所述最佳空间估计角度作为相对角度,针对每个所述基站,结合接收端与该基站对应的最佳空间估计角度以及该基站的所述姿态信息,得到所述接收端和该基站之间的绝对角度;得到所有基站与所述接收端之间的绝对角度后,结合所有基站的所述位置信息,解算得到所述接收端的空间位置。
CN202110263517.XA 2021-03-11 2021-03-11 一种基于发射角估计的音频定位方法 Active CN113099405B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110263517.XA CN113099405B (zh) 2021-03-11 2021-03-11 一种基于发射角估计的音频定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110263517.XA CN113099405B (zh) 2021-03-11 2021-03-11 一种基于发射角估计的音频定位方法

Publications (2)

Publication Number Publication Date
CN113099405A CN113099405A (zh) 2021-07-09
CN113099405B true CN113099405B (zh) 2021-12-17

Family

ID=76666998

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110263517.XA Active CN113099405B (zh) 2021-03-11 2021-03-11 一种基于发射角估计的音频定位方法

Country Status (1)

Country Link
CN (1) CN113099405B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105891808A (zh) * 2016-05-27 2016-08-24 复旦大学 一种多目标定位的声波发射装置
CN105988101A (zh) * 2015-03-02 2016-10-05 广州松下空调器有限公司 两设备之间的定位装置及定位方法
CN108226862A (zh) * 2016-12-15 2018-06-29 电信科学技术研究院 一种便携设备、信标及导航系统
CN110568408A (zh) * 2019-08-31 2019-12-13 苏州普息导航技术有限公司 基于单一信号源的音频定位系统及方法
CN110972033A (zh) * 2018-09-28 2020-04-07 硅实验室公司 基于一个或多个射频(rf)信号接收和/或传输特性来修改音频数据信息的系统和方法
CN111971989A (zh) * 2018-04-11 2020-11-20 高通股份有限公司 使用无线电波束的导航和定位系统
WO2021030444A1 (en) * 2019-08-12 2021-02-18 Qualcomm Incorporated Configurable coordinate system for angle reporting for positioning

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6933627B2 (ja) * 2018-10-31 2021-09-08 株式会社東芝 電子装置、到来角推定システム及び信号処理方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105988101A (zh) * 2015-03-02 2016-10-05 广州松下空调器有限公司 两设备之间的定位装置及定位方法
CN105891808A (zh) * 2016-05-27 2016-08-24 复旦大学 一种多目标定位的声波发射装置
CN108226862A (zh) * 2016-12-15 2018-06-29 电信科学技术研究院 一种便携设备、信标及导航系统
CN111971989A (zh) * 2018-04-11 2020-11-20 高通股份有限公司 使用无线电波束的导航和定位系统
CN110972033A (zh) * 2018-09-28 2020-04-07 硅实验室公司 基于一个或多个射频(rf)信号接收和/或传输特性来修改音频数据信息的系统和方法
WO2021030444A1 (en) * 2019-08-12 2021-02-18 Qualcomm Incorporated Configurable coordinate system for angle reporting for positioning
CN110568408A (zh) * 2019-08-31 2019-12-13 苏州普息导航技术有限公司 基于单一信号源的音频定位系统及方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
利用离去角度的无线传感器网络分布式节点定位方法;王文杰 等;《西安交通大学学报》;20100210;全文 *
到达角与出发角方法显著提升蓝牙位置服务性能;Dave Hollander;《电子产品世界》;20190504(第05期);全文 *
基于智能手机的室内定位技术的发展现状和挑战;陈亮 等;《测绘学报》;20171015;全文 *
基于离开角的商用毫米波设备定位方法研究;钱堃 等;《中国科学(信息科学)》;20210131;全文 *

Also Published As

Publication number Publication date
CN113099405A (zh) 2021-07-09

Similar Documents

Publication Publication Date Title
CN102045837B (zh) 移动节点定位方法及装置
CN1852591B (zh) 一种定位移动终端的方法
US20120013508A1 (en) Wireless localization techniques in lighting systems
JP5265631B2 (ja) 1つ以上の基地局に関して無線移動装置の位置を推定すること
JP2010151807A (ja) 無線測位装置及び座標構成方法
CN1979211B (zh) 一种定位系统和定位方法
CN102573050A (zh) 一种室内定位移动终端的系统及方法
CN109541662B (zh) 一种定位方法和装置
CN113347572A (zh) 利用空中基站实现终端定位的方法和系统
CN106888504B (zh) 基于fm与dtmb信号的室内位置指纹定位方法
Li et al. Review on positioning technology of wireless sensor networks
CN113099405B (zh) 一种基于发射角估计的音频定位方法
CN104683949A (zh) 一种应用于无线Mesh网中基于天线阵列的混合自定位方法
CN102035588A (zh) 基于角度信息的发射多播赋形方法及系统
El Abkari et al. Real time positioning over WSN and RFID network integration
CN110456307B (zh) 一种基于室内伪卫星信号载噪比的终端定位方法
CN211019247U (zh) 基于多通信融合的超宽带uwb定位系统
CN109856597B (zh) 一种新体制超视距短波定位系统及定位方法
Zhao et al. High accuracy indoor visible light positioning considering the shapes of illuminators
CN106908757A (zh) 一种单系统实现大范围室内外连续高精度定位的方法
CN100403850C (zh) 可适用于td-scdma网络的圆极化波定位方法
Brida et al. The accuracy of RSS based positioning in GSM networks
Juang et al. Hybrid SADOA/TDOA location estimation scheme for wireless communication systems
CN105282843A (zh) 一种基于方位角电平差值的定位方法及装置
CN114339600B (zh) 一种基于5g信号和声波信号的电子设备室内定位系统和方法

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