CN105539447B - 基于组合波形的汽车变道辅助***的信号处理方法和*** - Google Patents
基于组合波形的汽车变道辅助***的信号处理方法和*** Download PDFInfo
- Publication number
- CN105539447B CN105539447B CN201511009287.5A CN201511009287A CN105539447B CN 105539447 B CN105539447 B CN 105539447B CN 201511009287 A CN201511009287 A CN 201511009287A CN 105539447 B CN105539447 B CN 105539447B
- Authority
- CN
- China
- Prior art keywords
- point
- value
- calculated
- passage
- waveform
- 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
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W40/00—Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Mathematical Physics (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
基于组合波形的汽车变道辅助***的信号处理方法,属于信号处理领域,用于对汽车行驶后方,毫米波雷达波束覆盖范围内的目标,实现相对距离、相对速度以及方向角的解算,技术要点是:S1.对各段波形,A/D采集到的IQ数据,进行FFT计算;S2.将各段波形FFT变换后的复数模值做门限检测,输出过门限对应的复数模值最大的点的位置;S3.计算得到恒频段的多普勒频率值;S4.计算得到锯齿波段的差频频率值;S5.计算相对速度值。
Description
技术领域
本发明属于信号处理领域,尤其涉及一种汽车变道辅助***的信号处理方法。
背景技术
变道辅助***也称盲区监测,作为一项汽车智能安全技术,能够通过安装的电子控制***,在左右两个后视镜内或者其他地方提醒驾驶者后方安全范围内有无障碍物或来车,从而消除视线盲区,提高行车安全。
目前有代表性的变道辅助***有以下三类:
一、奔驰新E级上的变道辅助设备。该装置是通过汽车两侧安装的传感器探知两侧后方是否有车辆,并将信息通过电脑***控制,在左右两个后视镜内或者其他地方提醒驾驶者后方的来车。
二、沃尔沃的变道辅助叫盲点信息***。变道辅助***采用外后视镜根部的摄像头会对距离3米宽,9.5米长的一个扇形盲区进行25帧/秒的图像监控,如果有速度大于10公里/小时,且与车辆本身速度差在20-70公里/小时之间的移动物体(车辆或者行人)进入该盲区,***对比每帧图像,当***认为目标进一步接近时,A柱上的警示灯就会亮起,防止出现事故。
三、奥迪的变道辅助叫侧向辅助***。这套***会在车速超过60公里/小时介入,依靠传感器的帮助,奥迪侧向辅助***可以探测到侧后方最远50米处的车辆,若此时变道有潜在危险,后视镜上就会亮起警示灯。
发明内容
为了对汽车行驶后方,毫米波雷达波束覆盖范围内的目标,实现相对距离、相对速度以及方向角的解算,本发明提出了一种基于组合波形的汽车变道辅助***的信号处理方法。
本发明的技术方案是:一种基于组合波形的汽车变道辅助***的信号处理方法,所述波形包括恒频波CW和锯齿波FMCW,该方法包括以下步骤:
S1.对各段波形,A/D采集到的IQ数据,进行FFT计算;
S2.将各段波形FFT变换后的复数模值做门限检测,输出过门限对应的复数模值最大的点的位置;
S3.计算得到恒频段的多普勒频率值;
S4.计算得到锯齿波段的差频频率值;
进一步的,所述方法还包括步骤:S5.计算相对速度值。
进一步的,所述方法还包括步骤:S6.计算相对距离值。
进一步的,所述方法还包括步骤:S7.计算方向角。
进一步的,在步骤S1中,对通道1中的第一段恒频波CW、第二段锯齿波 FMCW,A/D采集到的IQ数据,选取各段线性度高的256点数据,分别进行 256点FFT,对通道2中的第一段恒频波CW,A/D采集到的IQ数据,选取各段线性度高的256点数据,分别进行256点FFT。
进一步的,步骤S2中,设通道1中,第一段恒频波CW,过门限点的复数模值最大的点的坐标为a,第二段锯齿波FMCW,过门限点的复数模值最大的点的坐标为b,设通道2中,第一段恒频波CW,过门限点的复数模值最大的点的坐标为c;若过门限的位置点等于1,则认为其是直流分量,不作为目标判定,直接剔除该位置点。
进一步的,通道1中,第一段恒频波CW,过门限点的复数模值最大的点的坐标为a,根据如下规则,则该点对应的多普勒频率为fda,通道2中,第一段恒频波CW,过门限点的复数模值最大的点的坐标为c,根据如下规则,则该点对应的多普勒频率为fdc;
该规则为,若点数为1≤x≤128,判断目标靠近,其对应点上的多普勒频率若点数为128<x≤256,判断目标远离,其对应点上的多普勒频率其中fs为采样频率。
进一步的,通道1中,第二段锯齿波FMCW,过门限点的复数模值最大的点的坐标为b,根据如下规则,则该点对应的差频频率值为fob;
该规则为,即若点数为1≤y≤128,其对应点上的差频频率值若点数为128<y≤256,其对应点上的差频频率值其中 fs为采样频率。
进一步的,根据得到的多普勒频率值fda,计算目标的速度v,计算速度公式为其中,c为光速,f为中心频率;
进一步的,根据恒频段计算得到的多普勒频率值fda,以及锯齿波段得到的差频频率值fob,计算目标的距离R,计算距离公式为其中,T 为每一段波形的作用时间,B为调频带宽,fd为多普勒频率值,fo为差频频率值。
进一步的,通道1中,第一段恒频波CW过门限点的位置a,根据如下计算方法计算对应点上的相位值,得到相位矩阵为ψa,通道2中,第一段恒频波CW 过门限点的位置b,根据如下计算方法计算对应点上的相位值,得到相位矩阵为ψb;计算通道1的第一段恒频波CW和通道2的第一段恒频波CW1之间的相位差,根据计算公式得到相位差矩阵为Δψac;根据公式方位角公式,得到目标的方向角,其中,d为天线间距,λ为波长。
本发明还涉及一种汽车变道辅助***,执行上述任一项技术方案的基于组合波形的汽车变道辅助***的信号处理方法。
有益效果:本发明给出了一种基于恒频波与锯齿波组合波形实现汽车变道辅助***的信号处理方法,为相关设计人员提供一种对于波形设计以及目标相关信息解算的新思路;所述方法,可以实现对单目标的相对距离以及相对速度的检测,同时可以实现对目标方向角的检测功能。
附图说明
图1基于组合波形的汽车变道辅助***信号处理流程图;
图2恒频波CW与锯齿波FMCW在一个扫频周期内的频率变化图。
具体实施方式
实施例1:一种基于组合波形的汽车变道辅助***的信号处理方法,所述波形包括恒频波CW和锯齿波FMCW,本实施例给出了在中心频率f为24.128GHz 工作频率下的具体波形图,第一段波形为恒频波CW,工作频率为24.128GHz,第二段波形为上升的锯齿波FMCW,工作频率变化范围为从24.128GHz变化到 24.278GHz,带宽为150MHz,每一段的周期T为5ms。恒频波CW与锯齿波 FMCW在一个扫频周期范围内的频率变化图如图1所示。
本实施例所涉及的方法,主要是对相对速度为80公里/小时,最大探测距离为50m,最小探测距离为1m的相邻车道上的目标车辆进行检测跟踪。根据探测到的目标的相对速度、相对距离以及方向角为变道提供警示依据等。
该方法包括以下步骤:
S1.对各段波形,A/D采集到的IQ数据,进行FFT计算;优选的,在步骤 S1中,对通道1中的第一段恒频波CW、第二段锯齿波FMCW,A/D采集到的 IQ数据,选取各段线性度高的256点数据,分别进行256点FFT,对通道2中的第一段恒频波CW,A/D采集到的IQ数据,选取各段线性度高的256点数据,分别进行256点FFT。其中,线性度高应该用指线性度好来表示会更好,线性度好意思就是采集的数据存在一定的线性关系,线性度不好,就是该段的数据没有呈现出一定的线性关系,表现出的是非线性关系等特点。为了更准确的分析数据,所以在采集到的数据中只选取呈现出线性关系的部分做FFT分析。
S2.将各段波形FFT变换后的复数模值做门限检测,输出过门限对应的复数模值最大的点的位置;优选的,设通道1中,第一段恒频波CW,过门限点的复数模值最大的点的坐标为a,第二段锯齿波FMCW,过门限点的复数模值最大的点的坐标为b,设通道2中,第一段恒频波CW,过门限点的复数模值最大的点的坐标为c;若过门限的位置点等于1,则认为其是直流分量,不作为目标判定,直接剔除该位置点。
S3.计算得到恒频段的多普勒频率值;优选的,通道1中,第一段恒频波CW,过门限点的复数模值最大的点的坐标为a,根据如下规则,则该点对应的多普勒频率为fda,通道2中,第一段恒频波CW,过门限点的复数模值最大的点的坐标为c,根据如下规则,则该点对应的多普勒频率为fdc;
该规则为,若点数为1≤x≤128,判断目标靠近,其对应点上的多普勒频率若点数为128<x≤256,判断目标远离,其对应点上的多普勒频率其中fs为采样频率。
S4.计算得到锯齿波段的差频频率值;优选的,通道1中,第二段锯齿波 FMCW,过门限点的复数模值最大的点的坐标为b,根据如下规则,则该点对应的差频频率值为fob;
该规则为,即若点数为1≤y≤128,其对应点上的差频频率值若点数为128<y≤256,其对应点上的差频频率值其中 fs为采样频率。
本实施例针对中心频率在24GHz或77GHz,基于恒频波调制的CW信号以及锯齿波调制的FMCW信号组合而成的波形,根据该调制波形实现一种汽车变道辅助***信号处理的方法。由该方法设计的汽车变道辅助***可以实现对汽车行驶后方,毫米波雷达波束覆盖范围内的目标,实现相对距离、相对速度以及方向角的解算。
实施例2:具有与实施例1任一项方案相同的技术方案,更为具体的是:所述方法还包括:S5.计算相对速度值。根据得到的多普勒频率值fda,计算目标的速度v,计算速度公式为其中,c为光速,c=3×108,f为中心频率,f=24.128GHz。
实施例3:具有与实施例1或2相同的技术方案,更为具体的是,该方法还包括步骤:S6.计算相对距离值。根据恒频段计算得到的多普勒频率值fda,以及锯齿波段得到的差频频率值fob,计算目标的距离R,计算距离公式为其中,T为每一段波形的作用时间,T=5ms,B为调频带宽, B=150MHz,fd为多普勒频率值,fo为差频频率值。
实施例4:具有与实施例1或2或3相同的技术方案,更为具体的是,该方法还包括步骤:S7.计算方向角。通道1中,第一段恒频波CW过门限点的位置a,根据如下计算方法计算对应点上的相位值,得到相位矩阵为ψa,通道2中,第一段恒频波CW过门限点的位置b,根据如下计算方法计算对应点上的相位值,得到相位矩阵为ψb;计算通道1的第一段恒频波CW和通道2的第一段恒频波 CW1之间的相位差,根据计算公式得到相位差矩阵为Δψac;根据公式方位角公式,得到目标的方向角,其中,d=7.5mm为天线间距,λ=12.4mm为波长。
实施例5:一种基于组合波形的汽车变道辅助***的信号处理方法,与实施例 1-4中装置中的装置相对应的方法,包括:
S1.对各段波形,A/D采集到的IQ数据,进行FFT计算;
S2.将各段波形FFT变换后的复数模值做门限检测,输出过门限对应的复数模值最大的点的位置;
S3.计算得到恒频段的多普勒频率值;
S4.计算得到锯齿波段的差频频率值。
S5.计算相对速度值;
S6.计算相对距离值。
S7.计算方向角。
实施例6:一种汽车变道辅助***,执行实施例1或2或3或4中任一项技术方案的基于组合波形的汽车变道辅助***的信号处理方法。
本实施例所设计恒频波与锯齿波的组合波形,以及波形所涉及的参数选取不局限于本实施例公开的参数,本领域技术人员可以根据具体应用场景,选取不同的设计参数,或是对波形进行改进等;在本发明披露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (9)
1.一种基于组合波形的汽车变道辅助***的信号处理方法,其特征在于,所述波形包括恒频波CW和锯齿波FMCW,该方法包括以下步骤:
S1.对各段波形,A/D采集到的IQ数据,进行FFT计算;
S2.将各段波形FFT变换后的复数模值做门限检测,输出过门限对应的复数模值最大的点的位置;
S3.计算得到恒频段的多普勒频率值;
S4.计算得到锯齿波段的差频频率值;
步骤S3中,通道1中,第一段恒频波CW,过门限点的复数模值最大的点的坐标为a,根据如下规则,则该点对应的多普勒频率为fda,通道2中,第一段恒频波CW,过门限点的复数模值最大的点的坐标为c,根据如下规则,则该点对应的多普勒频率为fdc;
该规则为,若点数为1≤x≤128,判断目标靠近,其对应点上的多普勒频率若点数为128<x≤256,判断目标远离,其对应点上的多普勒频率其中fs为采样频率。
2.如权利要求1所述的基于组合波形的汽车变道辅助***的信号处理方法,其特征在于,还包括步骤:
S5.计算相对速度值;
S6.计算相对距离值。
3.如权利要求1或2所述基于组合波形的汽车变道辅助***的信号处理方法,其特征在于,还包括步骤:S7.计算方向角。
4.如权利要求1所述的基于组合波形的汽车变道辅助***的信号处理方法,其特征在于,在步骤S1中,对通道1中的第一段恒频波CW、第二段锯齿波FMCW,A/D采集到的IQ数据,选取各段线性度高的256点数据,分别进行256点FFT,对通道2中的第一段恒频波CW,A/D采集到的IQ数据,选取各段线性度高的256点数据,分别进行256点FFT;步骤S2中,设通道1中,第一段恒频波CW,过门限点的复数模值最大的点的坐标为a,第二段锯齿波FMCW,过门限点的复数模值最大的点的坐标为b,设通道2中,第一段恒频波CW,过门限点的复数模值最大的点的坐标为c;若过门限的位置点等于1,则认为其是直流分量,不作为目标判定,直接剔除该位置点。
5.如权利要求1所述的方法,其特征在于,步骤S4中,通道1中,第二段锯齿波FMCW,过门限点的复数模值最大的点的坐标为b,根据如下规则,则该点对应的差频频率值为fob;
该规则为,若点数为1≤y≤128,其对应点上的差频频率值若点数为128<y≤256,其对应点上的差频频率值其中fs为采样频率。
6.如权利要求2所述的基于组合波形的汽车变道辅助***的信号处理方法,步骤S5中,根据得到的多普勒频率值fda,计算目标的速度v,计算速度公式为其中,c为光速,f为中心频率。
7.如权利要求2所述的基于组合波形的汽车变道辅助***的信号处理方法,其特征在于,步骤S6中,根据恒频段计算得到的多普勒频率值fda,以及锯齿波段得到的差频频率值fob,计算目标的距离R,计算距离公式为其中,T为每一段波形的作用时间,B为调频带宽,fd为多普勒频率值,fo为差频频率值。
8.如权利要求3所述的基于组合波形的汽车变道辅助***的信号处理方法,其特征在于,步骤S7中,通道1中,第一段恒频波CW过门限点的位置a,根据如下计算方法计算对应点上的相位值,得到相位矩阵为ψa,通道2中,第一段恒频波CW过门限点的位置b,根据如下计算方法计算对应点上的相位值,得到相位矩阵为ψb;计算通道1的第一段恒频波CW和通道2的第一段恒频波CW1之间的相位差,根据计算公式得到相位差矩阵为Δψac;根据公式方位角公式,得到目标的方向角,其中,d为天线间距,λ为波长。
9.一种汽车变道辅助***,执行权利要求1-8中任一项所述的基于组合波形的汽车变道辅助***的信号处理方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511009287.5A CN105539447B (zh) | 2015-12-29 | 2015-12-29 | 基于组合波形的汽车变道辅助***的信号处理方法和*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511009287.5A CN105539447B (zh) | 2015-12-29 | 2015-12-29 | 基于组合波形的汽车变道辅助***的信号处理方法和*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105539447A CN105539447A (zh) | 2016-05-04 |
CN105539447B true CN105539447B (zh) | 2018-03-06 |
Family
ID=55819121
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511009287.5A Active CN105539447B (zh) | 2015-12-29 | 2015-12-29 | 基于组合波形的汽车变道辅助***的信号处理方法和*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105539447B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107783114A (zh) * | 2016-08-25 | 2018-03-09 | 大连楼兰科技股份有限公司 | 旋翼无人机远距离复杂环境防撞毫米波雷达信号处理***及方法 |
CN107783090B (zh) * | 2016-08-25 | 2021-05-11 | 大连楼兰科技股份有限公司 | 基于毫米波雷达的固定翼无人机防撞***雷达信号处理方法 |
CN107783127B (zh) * | 2016-08-25 | 2021-04-09 | 大连楼兰科技股份有限公司 | 旋翼无人机防撞毫米波雷达信号处理方法 |
CN107783131B (zh) * | 2016-08-25 | 2021-05-18 | 大连楼兰科技股份有限公司 | 基于组合波形的自动驾驶汽车防撞雷达***信号处理方法 |
CN107783100B (zh) * | 2016-08-25 | 2021-02-19 | 大连楼兰科技股份有限公司 | 基于组合波形的旋翼无人机短距离防撞***信号处理方法 |
CN107783130B (zh) * | 2016-08-25 | 2021-02-19 | 大连楼兰科技股份有限公司 | 基于组合波形的无人驾驶汽车复杂环境防碰撞***信号处理方法 |
CN107783128B (zh) * | 2016-08-25 | 2021-02-19 | 大连楼兰科技股份有限公司 | 基于毫米波雷达的固定翼无人机多目标防撞*** |
CN107783101B (zh) * | 2016-08-25 | 2021-04-09 | 大连楼兰科技股份有限公司 | 基于组合波形的旋翼无人机短距离防撞***信号处理方法及装置 |
CN107783126B (zh) * | 2016-08-25 | 2021-03-30 | 大连楼兰科技股份有限公司 | 基于组合波形的自动驾驶汽车防撞雷达***信号处理方法及装置 |
CN107783112A (zh) * | 2016-08-25 | 2018-03-09 | 大连楼兰科技股份有限公司 | 基于组合波形的旋翼无人机复杂环境防碰撞***信号处理方法 |
CN109959925A (zh) * | 2017-12-25 | 2019-07-02 | 大连楼兰科技股份有限公司 | 一种电力巡线中信号数据处理方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20110109614A (ko) * | 2010-03-31 | 2011-10-06 | (주)블루포인트 | 차량의 경사각 보정방법 및 그 보정장치 |
CN104442812A (zh) * | 2013-10-16 | 2015-03-25 | 陈武强 | 一种汽车变道辅助探测装置及探测办法 |
CN104554232A (zh) * | 2013-10-17 | 2015-04-29 | 上海汽车集团股份有限公司 | 提高汽车变道安全性的方法和装置 |
CN104554277A (zh) * | 2015-01-26 | 2015-04-29 | 广州通易科技有限公司 | 一种基于重力向量和角速度向量的汽车变道识别方法 |
CN105182341A (zh) * | 2015-09-29 | 2015-12-23 | 西安知几天线技术有限公司 | 基于lfm三角波与恒频波组合波形的车辆防撞雷达多目标频率匹配方法 |
-
2015
- 2015-12-29 CN CN201511009287.5A patent/CN105539447B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20110109614A (ko) * | 2010-03-31 | 2011-10-06 | (주)블루포인트 | 차량의 경사각 보정방법 및 그 보정장치 |
CN104442812A (zh) * | 2013-10-16 | 2015-03-25 | 陈武强 | 一种汽车变道辅助探测装置及探测办法 |
CN104554232A (zh) * | 2013-10-17 | 2015-04-29 | 上海汽车集团股份有限公司 | 提高汽车变道安全性的方法和装置 |
CN104554277A (zh) * | 2015-01-26 | 2015-04-29 | 广州通易科技有限公司 | 一种基于重力向量和角速度向量的汽车变道识别方法 |
CN105182341A (zh) * | 2015-09-29 | 2015-12-23 | 西安知几天线技术有限公司 | 基于lfm三角波与恒频波组合波形的车辆防撞雷达多目标频率匹配方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105539447A (zh) | 2016-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105539447B (zh) | 基于组合波形的汽车变道辅助***的信号处理方法和*** | |
EP3141926B1 (en) | Automated detection of hazardous drifting vehicles by vehicle sensors | |
RU2666010C2 (ru) | Способ и устройство для определения интенсивности движения с помощью электронного контроллера, установленного в транспортном средстве | |
CN105652274A (zh) | 基于组合波形的汽车变道辅助***的信号处理装置 | |
CN106019281B (zh) | 物体检测装置以及物体检测方法 | |
CN105785370B (zh) | 物体检测装置和物体检测方法 | |
CN105109484B (zh) | 目标障碍物体确定方法及装置 | |
CN102288121B (zh) | 一种基于单目视觉的车道偏离距离测量及预警方法 | |
US8791802B2 (en) | Driver assistance system for reducing blind-spot-detection false alerts | |
JP5939357B2 (ja) | 移動軌跡予測装置及び移動軌跡予測方法 | |
KR101752858B1 (ko) | 레이더 기반 고 정밀 돌발상황 검지 시스템 | |
CN107111952A (zh) | 响应于极快接近的车辆的车道辅助*** | |
CN102431495B (zh) | 77GHz毫米波汽车主动防撞雷达弯道虚警抑制*** | |
CN104290753A (zh) | 一种前方车辆运动状态追踪预测装置及其预测方法 | |
US20240053468A1 (en) | Vehicle-mounted bsd millimeter wave radar based method for obstacle recognition at low speed | |
CN106428003B (zh) | 不良天气下高速公路车辆偏离车道的预警装置及方法 | |
KR102022539B1 (ko) | 차선 이탈 경보 시스템 | |
EP1657568B1 (en) | System and method for monitoring the external environment of a motor vehicle | |
JP2017223634A (ja) | レーダ装置及びピーク処理方法 | |
CN106405539B (zh) | 车辆雷达***以及用于去除一非感兴趣目标的方法 | |
KR102106976B1 (ko) | 도플러 정보를 이용한 차량의 후방 또는 사각지대 탐지 장치 및 그 방법 | |
CN107783130A (zh) | 基于组合波形的无人驾驶汽车复杂环境防碰撞***信号处理方法 | |
CN112098969B (zh) | 毫米波大车盲区雷达的目标检测与预警的优化方法 | |
CN114523902A (zh) | 一种前车开门场景的碰撞预警方法及*** | |
EP3816668A1 (en) | Vehicle control device and vehicle control method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |