CN112859067A - 一种基于毫米波雷达的动态目标的监测方法 - Google Patents
一种基于毫米波雷达的动态目标的监测方法 Download PDFInfo
- Publication number
- CN112859067A CN112859067A CN202011644479.4A CN202011644479A CN112859067A CN 112859067 A CN112859067 A CN 112859067A CN 202011644479 A CN202011644479 A CN 202011644479A CN 112859067 A CN112859067 A CN 112859067A
- Authority
- CN
- China
- Prior art keywords
- signal
- millimeter wave
- target
- wave radar
- frequency
- 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
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000012544 monitoring process Methods 0.000 title claims abstract description 21
- 230000003071 parasitic effect Effects 0.000 claims abstract description 10
- 238000007781 pre-processing Methods 0.000 claims abstract description 10
- 230000006855 networking Effects 0.000 claims abstract description 9
- 238000005259 measurement Methods 0.000 claims abstract description 8
- 238000005070 sampling Methods 0.000 claims description 6
- 230000003595 spectral effect Effects 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 5
- 238000005286 illumination Methods 0.000 claims description 4
- 238000003775 Density Functional Theory Methods 0.000 abstract description 9
- 230000003340 mental effect Effects 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 208000017667 Chronic Disease Diseases 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000000474 nursing effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/87—Combinations of radar systems, e.g. primary radar and secondary radar
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/08—Systems for measuring distance only
- G01S13/32—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
- G01S13/34—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of continuous, frequency-modulated waves while heterodyning the received signal, or a signal derived therefrom, with a locally-generated signal related to the contemporaneously transmitted signal
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/46—Indirect determination of position data
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/886—Radar or analogous systems specially adapted for specific applications for alarm systems
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/46—Indirect determination of position data
- G01S2013/462—Indirect determination of position data using multipath signals
-
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于毫米波雷达的动态目标的监测方法,包括如下步骤:S1利用毫米波雷达发射线性调频信号,经过回波信号与发射信号混频后得到基带差拍信号;S2对得到的信号进行预处理去除寄生调幅干扰、正交混频镜像干扰等影响信号质量的问题;S3利用DFT对含有目标距离信息的差拍信号进行频率的估计,继而转化成距离的估计而实现目标的距离参数测量;S4多个毫米波雷达分布式组网实现目标的三维空间定位。本发明提供的方法能够用于动态目标的监测,通过毫米波雷达获取的原始数据经过预处理、DFT计算和分布式组网实现动态目标的监测方法可用于老人、婴幼儿、精神病患者、大面积烧伤等特殊人群进行体征动态监护。
Description
技术领域
本发明涉及一种基于毫米波雷达的动态目标的监测方法,具体涉及雷达信号地发送、收集以及预处理,多个毫米波雷达分布式组网三维定位等。
背景技术
近些年,家庭空巢化趋势与独居老人增多,能够实现远程实时居家监控的智能医疗设备需求量不断增加。年轻一代也正在成为慢性病大军,具有患病时间长、服务需求量大的特点,如何合理有效地进行中年慢性疾病的管理,成为目前社会广泛关注的话题。
毫米波雷达是一种具有高频带、窄波束、高抗干扰的雷达***。在室内定位领域,不同于电红外、蓝牙、WIFI等传感器结合TOF和AOA等技术在准确性、误报和环境变化(如黑暗、亮度和烟雾)方面的局限性,基于毫米波雷达的室内人员定位和跌倒检测方法具有较强的穿透能力、非介入式传感、良好的隐私保护性,不易受环境因素(如天气、温度、光照等)影响等优势,可以在完全黑暗或明亮的白天工作。与可穿戴医疗设备相比,非接触式的毫米波雷达有效避免了因接触身体对人造成的不适,相比运动手环而言,毫米波雷达对检测到的数据运用更充分,具有较强的医疗价值,相比可穿戴医疗设备而言,毫米波雷达的灵活性更高,体验感更加舒适便捷,适用人群更广。在市场上也存在同样运用毫米波技术实现生命体征检测的产品,但是仍处于发展阶段,技术不够成熟,场景不够多元。
该方法涉及到利用毫米波雷达进行动态目标的监测,因而需要对采集到的数据进行计算分析。为了获取完整目标的方位信息,需要多个毫米波雷达分布式组网完成三维空间的定位。由于对三维定位的实时性要求较高,故寻求一种精度足够高,实时性较强的方法,能够快速完成雷达数据的采集计算分析,该方法为相关的人体定位的应用提供支持,以便后续很快地检测到老人摔倒或者出现其他异常情况时能够及时报警。
发明内容
为了更加高效便捷地采集毫米波雷达数据,同时也能为动态目标的定位实时提供原始数据,因而有必要提供了一种基于毫米波雷达的动态目标的监测方法。该方法的使用,能够较快完成差拍信号的获取、预处理、组网定位,具有较高精度、实时监测、抗干扰等优点。本发明的具体技术方案如下:
一种基于毫米波雷达的动态目标的监测方法,包括以下步骤:
S1:利用毫米波雷达发射线性调频信号f(t),经过回波信号R(t)与发射信号T(t)混频后得到基带的差拍信号ΔS(t);
S2:对得到的信号进行预处理去除寄生调幅干扰、正交混频镜像干扰;
S3:利用DFT对含有目标距离信息的差拍信号进行频率的估计,继而转化成距离的估计而实现目标的距离参数测量;
S4:多个毫米波雷达分布式组网实现目标的三维空间定位。
进一步的,S1中,所述利用毫米波雷达发射线性调频信号f(t),经过回波信号R(t)与发射信号T(t)混频后得到基带的差拍信号ΔS(t),包括:
毫米波雷达发射线性调频信号f(t);
发射信号T(t)和回波信号R(t)混频后得到基带的差拍信号ΔS(t)。
进一步的,S2,所述对得到的信号进行预处理,包括:
利用寄生调幅干扰波形相对稳定的特点,获取雷达空照时的差拍信号可以近似为寄生调幅信号;
采用主瓣能量大,旁瓣能量小的哈布斯窗函数抑制旁瓣对领近谱峰主瓣的影响,即正交混频镜像干扰。
进一步的,S3,所述对预处理后的差拍信号ΔS(t)进行N点DFT得到距离信息,其方式包括:
对ΔS(t)进行采样率fs的N点采样得到序列ΔS(n);
对序列ΔS(n)进行N点DFT求出最高谱线处的频率实现对ΔS(t)频率的估计,继而转化成距离的估计而实现目标的距离参数测量。
进一步的,S4,所述多个毫米波雷达分布式组网实现目标的定位,其中包括:
多个雷达相互独立测距,建立球面方程实现三维定位;
计算球面方程的最小二乘解作为目标所在位置。
本发明的基于毫米波雷达的动态目标的监测方法,能够用于动态目标的监测,通过毫米波雷达获取的原始数据经过预处理、DFT计算和分布式组网实现动态目标的监测方法可用于老人、婴幼儿、精神病患者、大面积烧伤等特殊人群进行体征动态监护。
附图说明
图1为实施例提供的基于毫米波雷达的动态目标的监测方法流程图。
图2为实施例中毫米波雷达三维组网示意图。
具体实施方式
为了更清楚地说明本申请实施例的技术方案,下面对实施例中所需要的附图作简单的介绍。
下面结合附图对本发明作进一步说明。
本发明实例提供的毫米波雷达动态目标定位的方法,是用于实时采集毫米波雷达数据,计算处理得到目标的三维定位,可为养老院、医院、家庭等的应用提供支持。以下结合附图对本发明实施例的具体实施方式进行详细说明。应当理解的是,此处所描述的具体实施方式仅用于说明和解释本发明实施例,并不用于限制本发明实施例。
如图1是为申请实施例提供的基于毫米波雷达的动态目标的监测方法流程图。该方法包括S1利用毫米波雷达发射线性调频信号,经过回波信号与发射信号混频后得到基带差拍信号;S2对得到的信号进行预处理去除寄生调幅干扰、正交混频镜像干扰;S3利用DFT对含有目标距离信息的差拍信号进行频率的估计,继而转化成距离的估计实现目标的参数测量;S4多个毫米波雷达分布式组网实现目标的三维空间定位。
发射后经过t1=2R/c时间得到的回波信号为:
混频后得到基带差拍信号为:
ΔS(t)=exp{j(T(t)-R(t))}
步骤S2中,对于寄生调幅干扰s(t),会对差拍信号ΔS(t)造成影响和干扰,即会在调频周期内差拍信号上叠加一个调制信号:
ΔS*(t)=ΔS(t)+s(t)
此时利用寄生调幅干扰相对稳定的特点,先获取雷达空照是对应的差拍信号ΔS1(t),这个可以近似为寄生调幅信号,即ΔS1(t)≈s(t),因此可以求出目标真正的差拍信号,从而抑制寄生调幅干扰的影响。
对于正交混频镜像干扰,可以利用窗函数来抑制旁瓣对领近谱峰主瓣的影响,这里可以选用三阶哈布斯窗:
其中N为采样点数,k=3,a0=0.35875,a1=0.48829,a2=0.14128,a3=0.01168。步骤S3中,要求的目标距离R就在进行过预处理过的ΔS(t)中,将t1=2R/c代入并对其进行采样率fs的N点采样得到序列ΔS(n),接着对其求N点DFT求出信号频谱X(k)最高谱线处的频率fm,当增加DFT点数实现频谱的细化,即频谱谱线距离更加密集,最高谱线更加接近于真实频率值,以此实现对ΔS(t)频率的估计,可以提高估计的准确度,目标距离R为:
步骤S4中,多个毫米波雷达相互独立测距,可以建立球面方程实现三维定位,只需要每个毫米波雷达测量距离,不需要角度的测量,不需要计算相位差,因此计算量适度、定位精度较高。
如图2所示,设三个固定雷达的位置互相垂直,坐标分别为(xi,yi,zi),i=1,2,3,目标坐标为(x,y,z),与雷达距离为Ri,
由于实际环境中存在误差,测量值和真实值存在一定的误差,因此设
即求解该方程组的最小二乘解,该解是当前测量条件下最接近真实值的解,就是目标的坐标,Newton法求解非线性方程组的具体过程如下:
忽略余项Ri得到:
f(x0)+J(x0)(x-x0)=0
求其最小二乘解x1
x1=x0-J-1(x0)(x0)
其中J-1应理解成Jacobi阵的广义逆矩阵,因此每次迭代的值:
xk+1=xk-J-1(xk)(xk)
随着迭代误差逐渐减小小于门限值或达到设定得最大迭代次数时,此时可得到最小二乘解。
并且由于在很小时间内动态目标的运动位移很小,可以视为静态目标,因此此方法可以实现较高精度,较高实时性的动态目标监测。
Claims (5)
1.一种基于毫米波雷达的动态目标的监测方法,其特征在于包括以下步骤:
S1:利用毫米波雷达发射线性调频信号f(t),经过回波信号R(t)与发射信号T(t)混频后得到基带的差拍信号ΔS(t);
S2:对得到的信号进行预处理去除寄生调幅干扰、正交混频镜像干扰;
S3:利用DFT对含有目标距离信息的差拍信号进行频率的估计,继而转化成距离的估计而实现目标的距离参数测量;
S4:多个毫米波雷达分布式组网实现目标的三维空间定位。
2.根据权利要求1所述的一种基于毫米波雷达的动态目标的监测方法,其特征在于S1中,所述利用毫米波雷达发射线性调频信号f(t),经过回波信号R(t)与发射信号T(t)混频后得到基带的差拍信号ΔS(t),包括:
毫米波雷达发射线性调频信号f(t);
发射信号T(t)和回波信号R(t)混频后得到基带的差拍信号ΔS(t)。
3.根据权利要求1所述的一种基于毫米波雷达的动态目标的监测方法,其特征在于S2,所述对得到的信号进行预处理,包括:
利用寄生调幅干扰波形相对稳定的特点,获取雷达空照时的差拍信号可以近似为寄生调幅信号;
采用主瓣能量大,旁瓣能量小的哈布斯窗函数抑制旁瓣对领近谱峰主瓣的影响,即正交混频镜像干扰。
4.根据权利要求1所述的一种基于毫米波雷达的动态目标的监测方法,其特征在于S3,所述对预处理后的差拍信号ΔS(t)进行N点DFT得到距离信息,其方式包括:
对ΔS(t)进行采样率fs的N点采样得到序列ΔS(n);
对序列ΔS(n)进行N点DFT求出最高谱线处的频率实现对ΔS(t)频率的估计,继而转化成距离的估计而实现目标的距离参数测量。
5.根据权利要求1所述的一种基于毫米波雷达的动态目标的监测方法,其特征在于S4,所述多个毫米波雷达分布式组网实现目标的定位,其中包括:
多个雷达相互独立测距,建立球面方程实现三维定位;
计算球面方程的最小二乘解作为目标所在位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011644479.4A CN112859067B (zh) | 2020-12-31 | 2020-12-31 | 一种基于毫米波雷达的动态目标的监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011644479.4A CN112859067B (zh) | 2020-12-31 | 2020-12-31 | 一种基于毫米波雷达的动态目标的监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112859067A true CN112859067A (zh) | 2021-05-28 |
CN112859067B CN112859067B (zh) | 2024-01-16 |
Family
ID=76000956
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011644479.4A Active CN112859067B (zh) | 2020-12-31 | 2020-12-31 | 一种基于毫米波雷达的动态目标的监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112859067B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102170318A (zh) * | 2011-01-11 | 2011-08-31 | 南京国睿安泰信科技股份有限公司 | 用于接收机性能测试的频谱分析算法 |
CN102621550A (zh) * | 2011-01-28 | 2012-08-01 | 上海无线电设备研究所 | 一种w波段前车防撞雷达探测***的实现方法 |
CN107907878A (zh) * | 2017-11-08 | 2018-04-13 | 零八电子集团有限公司 | 高精度获取fmcw雷达距离测量值的方法 |
CN109959930A (zh) * | 2017-12-22 | 2019-07-02 | 英飞凌科技股份有限公司 | 使用毫米波雷达传感器监测结构对象的***和方法 |
CN110118966A (zh) * | 2019-05-28 | 2019-08-13 | 长沙莫之比智能科技有限公司 | 基于毫米波雷达的人员检测与计数*** |
-
2020
- 2020-12-31 CN CN202011644479.4A patent/CN112859067B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102170318A (zh) * | 2011-01-11 | 2011-08-31 | 南京国睿安泰信科技股份有限公司 | 用于接收机性能测试的频谱分析算法 |
CN102621550A (zh) * | 2011-01-28 | 2012-08-01 | 上海无线电设备研究所 | 一种w波段前车防撞雷达探测***的实现方法 |
CN107907878A (zh) * | 2017-11-08 | 2018-04-13 | 零八电子集团有限公司 | 高精度获取fmcw雷达距离测量值的方法 |
CN109959930A (zh) * | 2017-12-22 | 2019-07-02 | 英飞凌科技股份有限公司 | 使用毫米波雷达传感器监测结构对象的***和方法 |
CN110118966A (zh) * | 2019-05-28 | 2019-08-13 | 长沙莫之比智能科技有限公司 | 基于毫米波雷达的人员检测与计数*** |
Also Published As
Publication number | Publication date |
---|---|
CN112859067B (zh) | 2024-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | mmHRV: Contactless heart rate variability monitoring using millimeter-wave radio | |
Ahmad et al. | Vital signs monitoring of multiple people using a FMCW millimeter-wave sensor | |
Mercuri et al. | SFCW microwave radar for in-door fall detection | |
Yang et al. | Multi-breath: Separate respiration monitoring for multiple persons with UWB radar | |
CN113261925A (zh) | 一种基于毫米波雷达的人体感知方法及*** | |
CN111481184A (zh) | 基于毫米波雷达技术的多目标呼吸心率监测方法与*** | |
Khan et al. | Wireless health monitoring using passive WiFi sensing | |
WO2021248472A1 (zh) | 基于超宽带雷达的目标跟踪方法、装置、设备及存储介质 | |
Baboli et al. | A new algorithm for detection of heart and respiration rate with UWB signals | |
Li et al. | Opportunistic physical activity monitoring via passive WiFi radar | |
CN114366052A (zh) | 一种基于毫米波雷达的养老院智能监护***及方法 | |
Liu et al. | mmRH: Noncontact vital sign detection with an FMCW mm-wave radar | |
CN115736886A (zh) | 一种基于毫米波雷达的多目标呼吸率参数估计方法 | |
Upadhyay et al. | Vital sign detection via angular and range measurements with mmWave MIMO radars: Algorithms and trials | |
CN115469301A (zh) | 基于超宽带雷达的老年人活动能力评估方法及*** | |
CN116942142A (zh) | 一种基于毫米波雷达的人员智能化监测方法及相关设备 | |
Raheel et al. | Breathing and heartrate monitoring system using IR-UWB radar | |
Wang et al. | Noncontact extraction of biomechanical parameters in gait analysis using a multi-input and multi-output radar sensor | |
Zhang et al. | Lt-fall: The design and implementation of a life-threatening fall detection and alarming system | |
KR20140106795A (ko) | Uwb 기반 비접촉식 생체신호진단기 | |
Schroth et al. | Emergency response person localization and vital sign estimation using a semi-autonomous robot mounted SFCW radar | |
CN112859067A (zh) | 一种基于毫米波雷达的动态目标的监测方法 | |
KR20180088019A (ko) | 레이더를 이용한 생체 신호 측정 방법 및 장치 | |
Almudevar et al. | Home monitoring using wearable radio frequency transmitters | |
Xue et al. | Accurate multi-target vital signs detection method for FMCW radar |
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 |