CN111024218B - 基于自动跟踪的超声波相关检测方法 - Google Patents
基于自动跟踪的超声波相关检测方法 Download PDFInfo
- Publication number
- CN111024218B CN111024218B CN201911383600.XA CN201911383600A CN111024218B CN 111024218 B CN111024218 B CN 111024218B CN 201911383600 A CN201911383600 A CN 201911383600A CN 111024218 B CN111024218 B CN 111024218B
- Authority
- CN
- China
- Prior art keywords
- time
- waveform
- value
- real
- max
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于自动跟踪的超声波相关检测方法,该方法包括以下步骤:包括以下步骤:S1:计算参考波形的所有的极值,并提取以最大值为中心对称的2m+1个极值点建立参考波形极值序列XVm;S2:计算实时接收波形的所有的极值,并提取以最大值为中心对称的2n+1个极值点建立实时波形极值序列XSn;S3:将参考波形极值序列XVm和实时波形极值序列XSn进行滑动相关计算得到参考波形极值序列XVm和实时波形极值序列XSn的最大相关值τmax以及本次超声波渡越时间t;S4:将最大相关值τmax与预设预置有效阈值进行比较当,最大相关值τmax大于等于有效阈值时,本次测量结果有效,输出测量结果t;否则,本次测量结果无效并建立激励重发机制。
Description
技术领域
本发明涉及一种基于自动跟踪的超声波相关检测方法。
背景技术
如图1所示,在超声波换能器能量转换和声波空气传输过程中,由于机械弹性(超声波换能器压电陶瓷片和空气的弹性)的存在,在接收到的超声波脉冲串的头部和尾部会产生超声波波形振幅衰减,从起振开始到最后停止至少需要10多个周期,包络呈现纺锤体形状。超声波到达时的强度(振幅)存在随机变化的成分,采用常规的阈值检测方法判断超声波到达的时间点是不准确的,不能准确地检测超声波渡越时间。一般仅用于机器人避障、倒车防撞装置等要求不高的场合。为了准确地判断超声波到达的时间点,准确检测超声波渡越时间,人们进行了很多探索和改进,目前应用较多的主要有以下4种:
(1)阈值检测法。接收判识简单,但是接收波形幅度变化直接影响测量精度。
(2)相位差法。该方法短距离应用非常准确,已应用于超声波风速仪,自动气象站的大气测温等。但该方法只能用于超声波(机械能波)传输距离很短的场合。
(3)调制波形标记法。对声波进行调制、编码、扩频等均属于此类。应用较多的是双频频移键控法,已用于轿车空调测温和超声波测距仪。
(4)相关检测法或包络检测法。这种方法能有效地利用波形的整体信息,弱化突发噪声、白噪声对测量结果的影响。但是参考波形的合理选择和计算量大一定程度上限制了它的应用。
发明内容
本发明的目的是提供一种基于自动跟踪的超声波相关检测方法,利用参考波形和实时接收波形进行极值相关运算,针对运算结果实时评审参考波形的合理性,差异较大且在有效范围内,直接用当前接收波形替换参考波形,实现参考波形自动跟踪接收波形,周而复始进行超声波接收极值相关检测。
为解决上述技术问题,本发明提供一种基于自动跟踪的超声波相关检测方法,该方法包括以下步骤:
S1:计算参考波形的所有的极值,并提取以最大值为中心对称的2m+1个极值点建立参考波形极值序列XVm;
S2:计算实时接收波形的所有的极值,并提取以最大值为中心对称的2n+1个极值点建立实时波形极值序列XSn;
S3:将参考波形极值序列XVm和实时波形极值序列XSn进行滑动相关计算得到参考波形极值序列XVm和实时波形极值序列XSn的最大相关值τmax以及本次超声波渡越时间t;
S4:将最大相关值τmax与预设预置有效阈值进行比较当,最大相关值τmax大于等于有效阈值时,本次测量结果有效,输出测量结果t;否则,本次测量结果无效并建立激励重发机制。
进一步地,所述步骤S1具体包括:
S11:从参考波形中提取信号数据序列Vm,计算参考波形的所有的极值,同时获取参考波形的最大峰值PV和参考波形的起始时刻TV0;
S12:以PV为中心对称提取2m+1个极值点建立参考极值序列XVm,根据取出的2m+1个极值点的时刻点构建参考时间序列TVm。
进一步地,所述步骤S2具体包括:
S21:从实时接收波形中提取信号数据序列Sn,计算实时接收波形的所有的极值,同时获取实时接收波形的最大峰值PS和实时接收波形的起始时刻TS0;
S22:以PS为中心对称提取2n+1个极值点建立实时接收极值序列XSn,根据取出的2n+1个极值点的时刻点构建实时接收时间序列TSn;其中,n≥m。
进一步地,所述步骤S3具体包括:
S31:将参考波形极值序列XVm和实时波形极值序列XSn进行滑动相关计算所采用的具体公式如下:
其中,序数i的取值范围0~m;ΔN为采样点序数,取值范围为0~(n-m);
S32:获取τ(△N)的最大值作为最大相关值τmax,最大相关值τmax对应最大相关值对应序数为ΔNmax,则本次超声波渡越时间t所采用的具体公式为:
Δtmax=ΔNmax/fs t=TS0-TV0+Δtmax (2)
其中,△tmax为最大相关时间差,且Δtmax=ΔNmax/fs,fs为采样速率。
进一步地,该方法还包括步骤S5:
当最大相关值τmax大于等于有效阈值时,则计算该实时接收波形的可信度值,当计算得到的所述可信度值在有效阈值区间(ρth,ρtv)时,则自动更新参考波形,实现自动跟踪。
进一步地,计算所述实时接收波形的可信度值所采用的具体公式为:
当ρ<ρth时:本次结果判为无效,波形重发;
当ρth≤ρ<ρtv时:本次结果有效,超声波渡越时间为t=TS0-TV0+Δtmax,本次实时接收波形替换当前参考波形,实现实时接收波形自动跟踪;
当ρ≥ρtv时:本次结果有效,超声波渡越时间为t=TS0-TV0+Δtmax,参考波形继续使用不替换。
本发明的有益效果为:通过采用超声波实时接收信号极值相关检测法判断实时接收波形有效性,该方法大大减少了参与乘加运算的数据量,比普通整体计算量小;并且,波形重发机制消除了无效接收波形的影响,提高了测量结果的准确性。
附图说明
此处所说明的附图用来提供对本申请的进一步理解,构成本申请的一部分,在这些附图中使用相同的参考标号来表示相同或相似的部分,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1为超声波接收信号极值图。
具体实施方式
一种基于自动跟踪的超声波相关检测方法,该方法包括以下步骤:
S1:计算参考波形的所有的极值,并提取以最大值为中心对称的2m+1个极值点建立参考波形极值序列XVm;步骤S1具体包括:
S11:从参考波形中提取信号数据序列Vm,计算参考波形的所有的极值,同时获取参考波形的最大峰值PV和参考波形的起始时刻TV0;
S12:以PV为中心对称提取2m+1个极值点建立参考极值序列XVm,根据取出的2m+1个极值点的时刻点构建参考时间序列TVm;m取值范围为5~20。
S2:计算实时接收波形的所有的极值,并提取以最大值为中心对称的2n+1个极值点建立实时波形极值序列XSn;步骤S2具体包括:
S21:从实时接收波形中提取信号数据序列Sn,计算实时接收波形的所有的极值,同时获取实时接收波形的最大峰值PS和实时接收波形的起始时刻TS0;
S22:以PS为中心对称提取2n+1个极值点建立实时接收极值序列XSn,根据取出的2n+1个极值点的时刻点构建实时接收时间序列TSn;其中,n≥m,且n取值范围为10~30。
S3:将参考波形极值序列XVm和实时波形极值序列XSn进行滑动相关计算得到参考波形极值序列XVm和实时波形极值序列XSn的最大相关值τmax以及本次超声波渡越时间t;步骤S3具体包括:
S31:将参考波形极值序列XVm和实时波形极值序列XSn进行滑动相关计算所采用的具体公式如下:
其中,序数i的取值范围0~m;ΔN为采样点序数,取值范围为0~(n-m);
S32:获取τ(△N)的最大值作为最大相关值τmax,最大相关值τmax对应最大相关值对应序数为ΔNmax,则本次超声波渡越时间t所采用的具体公式为:
Δtmax=ΔNmax/fs t=TS0-TV0+Δtmax (2)
其中,△tmax为最大相关时间差,且Δtmax=ΔNmax/fs,fs为采样速率。
S4:将最大相关值τmax与预设预置有效阈值进行比较当,最大相关值τmax大于等于有效阈值时,本次测量结果有效,输出测量结果t;否则,本次测量结果无效并建立激励重发机制。
该方法通过采用超声波实时接收信号极值相关检测法判断实时接收波形有效性,该方法大大减少了参与乘加运算的数据量,比普通整体计算量小;并且,波形重发机制消除了无效接收波形的影响,提高了测量结果的准确性。
S5:当最大相关值τmax大于等于有效阈值时,则计算该实时接收波形的可信度值,当计算得到的所述可信度值在有效阈值区间(ρth,ρtv)时,则自动更新参考波形,实现自动跟踪;计算所述实时接收波形的可信度值所采用的具体公式为:
当ρ<ρth时:本次结果判为无效,波形重发;
当ρth≤ρ<ρtv时:本次结果有效,超声波渡越时间为t=TS0-TV0+Δtmax,本次实时接收波形替换当前参考波形,实现实时接收波形自动跟踪;
当ρ≥ρtv时:本次结果有效,超声波渡越时间为t=TS0-TV0+Δtmax,参考波形继续使用不替换。
该方法通过采用在有效范围内判定为本次接收波形有效并给出可信度值,根据可信度值判断接收波形数据是否需要替换参考波形,实现参考波形自动跟踪实时接收波形,避免参考波形一成不变,提高了环境适应性。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (4)
1.一种基于自动跟踪的超声波相关检测方法,其特征在于,包括以下步骤:
S1:计算参考波形的所有的极值,并提取以最大值为中心对称的2m+1个极值点建立参考波形极值序列XVm;
S2:计算实时接收波形的所有的极值,并提取以最大值为中心对称的2n+1个极值点建立实时波形极值序列XSn;
S3:将参考波形极值序列XVm和实时波形极值序列XSn进行滑动相关计算得到参考波形极值序列XVm和实时波形极值序列XSn的最大相关值τmax以及本次超声波渡越时间t;
S4:将最大相关值τmax与预设预置有效阈值进行比较,当最大相关值τmax大于等于有效阈值时,本次测量结果有效,输出测量结果t;否则,本次测量结果无效并建立激励重发机制;
S5:当最大相关值τmax大于等于有效阈值时,则计算该实时接收波形的可信度值,当计算得到的所述可信度值在有效阈值区间(ρth,ρtv)时,则自动更新参考波形,实现自动跟踪;计算所述实时接收波形的可信度值所采用的具体公式为:
当ρ<ρth时:本次结果判为无效,波形重发;
当ρth≤ρ<ρtv时:本次结果有效,超声波渡越时间为t=TS0-TV0+Δtmax,本次实时接收波形替换当前参考波形,实现实时接收波形自动跟踪;
当ρ≥ρtv时:本次结果有效,超声波渡越时间为t=TS0-TV0+Δtmax,参考波形继续使用不替换。
2.根据权利要求1所述的基于自动跟踪的超声波相关检测方法,其特征在于,所述步骤S1具体包括:
S11:从参考波形中提取信号数据序列Vm,计算参考波形的所有的极值,同时获取参考波形的最大峰值PV和参考波形的起始时刻TV0;
S12:以PV为中心对称提取2m+1个极值点建立参考极值序列XVm,根据取出的2m+1个极值点的时刻点构建参考时间序列TVm。
3.根据权利要求2所述的基于自动跟踪的超声波相关检测方法,其特征在于,所述步骤S2具体包括:
S21:从实时接收波形中提取信号数据序列Sn,计算实时接收波形的所有的极值,同时获取实时接收波形的最大峰值PS和实时接收波形的起始时刻TS0;
S22:以PS为中心对称提取2n+1个极值点建立实时接收极值序列XSn,根据取出的2n+1个极值点的时刻点构建实时接收时间序列TSn;其中,n≥m。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911383600.XA CN111024218B (zh) | 2019-12-27 | 2019-12-27 | 基于自动跟踪的超声波相关检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911383600.XA CN111024218B (zh) | 2019-12-27 | 2019-12-27 | 基于自动跟踪的超声波相关检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111024218A CN111024218A (zh) | 2020-04-17 |
CN111024218B true CN111024218B (zh) | 2021-11-26 |
Family
ID=70194906
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911383600.XA Active CN111024218B (zh) | 2019-12-27 | 2019-12-27 | 基于自动跟踪的超声波相关检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111024218B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111693732B (zh) * | 2020-06-24 | 2021-12-24 | 中煤科工集团重庆研究院有限公司 | 基于滑动参考波形的超声波渡越时间互相关计算方法 |
Citations (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08160060A (ja) * | 1994-12-02 | 1996-06-21 | Kaijo Corp | 超音波利用測定装置の受信波利得制御方法 |
JP2001324318A (ja) * | 2000-05-15 | 2001-11-22 | Non-Destructive Inspection Co Ltd | 配管等構造体における膜厚測定方法 |
CN1804611A (zh) * | 2005-12-27 | 2006-07-19 | 哈尔滨工业大学 | 基于合成孔径聚焦的超声渡越时间检测方法 |
CN102238915A (zh) * | 2008-12-02 | 2011-11-09 | 株式会社东芝 | 超声波诊断装置、多普勒测量装置和多普勒测量方法 |
CN102645552A (zh) * | 2012-04-19 | 2012-08-22 | 北京航空航天大学 | 一种用于渡越时间测量的自适应相关算法 |
CN103464359A (zh) * | 2013-08-21 | 2013-12-25 | 惠州市诺丹富超声波设备有限公司 | 一种超声波发生器的自适应方法 |
CN105181045A (zh) * | 2015-05-22 | 2015-12-23 | 重庆川仪自动化股份有限公司 | 超声波流量计消除温度影响声速的方法 |
CN105606170A (zh) * | 2015-12-31 | 2016-05-25 | 浙江威星智能仪表股份有限公司 | 一种带自学习模板的超声波气体计量装置及方法 |
CN106248157A (zh) * | 2016-06-27 | 2016-12-21 | 浙江大学 | 互相关时差法气体超声波流量计参考波形的确定方法 |
CN106643939A (zh) * | 2017-02-20 | 2017-05-10 | 重庆川仪自动化股份有限公司 | 用于超声波流量计计算超声波传播时间的方法 |
CN106643843A (zh) * | 2016-12-14 | 2017-05-10 | 广州大学 | 超声波检测装置中信号电路延迟时间的检测装置和方法 |
CN106768109A (zh) * | 2017-02-21 | 2017-05-31 | 合肥工业大学 | 基于回波上升段峰值拟合和基于回波能量点定位的气体超声波流量计信号处理方法 |
CN106768297A (zh) * | 2016-12-28 | 2017-05-31 | 清华大学苏州汽车研究院(吴江) | 一种全自动超声波指向性测试系统 |
CN106813602A (zh) * | 2015-11-30 | 2017-06-09 | 江苏中海达海洋信息技术有限公司 | 一种基于跳频信号的全自动超声波测深仪信号处理方法 |
CN106871980A (zh) * | 2017-02-13 | 2017-06-20 | 合肥工业大学 | 基于相邻峰值最大差值的气体超声流量计间歇式激励和信号处理方法和系统 |
CN108052352A (zh) * | 2017-11-30 | 2018-05-18 | 努比亚技术有限公司 | 屏幕开关状态的控制方法、终端及计算机可读存储介质 |
JP2018138919A (ja) * | 2018-03-27 | 2018-09-06 | 辻野 次郎丸 | 超音波複合振動加工装置の振動検出器 |
CN109029598A (zh) * | 2018-06-15 | 2018-12-18 | 浙江大学 | 一种互相关法气体超声波流量计中基于主次包络结合处的参考波形选取方法 |
CN109297550A (zh) * | 2018-09-21 | 2019-02-01 | 浙江大学 | 基于互相关函数包络特征点的超声波渡越时间计算方法 |
-
2019
- 2019-12-27 CN CN201911383600.XA patent/CN111024218B/zh active Active
Patent Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08160060A (ja) * | 1994-12-02 | 1996-06-21 | Kaijo Corp | 超音波利用測定装置の受信波利得制御方法 |
JP2001324318A (ja) * | 2000-05-15 | 2001-11-22 | Non-Destructive Inspection Co Ltd | 配管等構造体における膜厚測定方法 |
CN1804611A (zh) * | 2005-12-27 | 2006-07-19 | 哈尔滨工业大学 | 基于合成孔径聚焦的超声渡越时间检测方法 |
CN102238915A (zh) * | 2008-12-02 | 2011-11-09 | 株式会社东芝 | 超声波诊断装置、多普勒测量装置和多普勒测量方法 |
CN103431876A (zh) * | 2008-12-02 | 2013-12-11 | 株式会社东芝 | 超声波诊断装置、多普勒测量装置和多普勒测量方法 |
CN102645552A (zh) * | 2012-04-19 | 2012-08-22 | 北京航空航天大学 | 一种用于渡越时间测量的自适应相关算法 |
CN103464359A (zh) * | 2013-08-21 | 2013-12-25 | 惠州市诺丹富超声波设备有限公司 | 一种超声波发生器的自适应方法 |
CN105181045A (zh) * | 2015-05-22 | 2015-12-23 | 重庆川仪自动化股份有限公司 | 超声波流量计消除温度影响声速的方法 |
CN106813602A (zh) * | 2015-11-30 | 2017-06-09 | 江苏中海达海洋信息技术有限公司 | 一种基于跳频信号的全自动超声波测深仪信号处理方法 |
CN105606170A (zh) * | 2015-12-31 | 2016-05-25 | 浙江威星智能仪表股份有限公司 | 一种带自学习模板的超声波气体计量装置及方法 |
CN106248157A (zh) * | 2016-06-27 | 2016-12-21 | 浙江大学 | 互相关时差法气体超声波流量计参考波形的确定方法 |
CN106643843A (zh) * | 2016-12-14 | 2017-05-10 | 广州大学 | 超声波检测装置中信号电路延迟时间的检测装置和方法 |
CN106768297A (zh) * | 2016-12-28 | 2017-05-31 | 清华大学苏州汽车研究院(吴江) | 一种全自动超声波指向性测试系统 |
CN106871980A (zh) * | 2017-02-13 | 2017-06-20 | 合肥工业大学 | 基于相邻峰值最大差值的气体超声流量计间歇式激励和信号处理方法和系统 |
CN106643939A (zh) * | 2017-02-20 | 2017-05-10 | 重庆川仪自动化股份有限公司 | 用于超声波流量计计算超声波传播时间的方法 |
CN106768109A (zh) * | 2017-02-21 | 2017-05-31 | 合肥工业大学 | 基于回波上升段峰值拟合和基于回波能量点定位的气体超声波流量计信号处理方法 |
CN108052352A (zh) * | 2017-11-30 | 2018-05-18 | 努比亚技术有限公司 | 屏幕开关状态的控制方法、终端及计算机可读存储介质 |
JP2018138919A (ja) * | 2018-03-27 | 2018-09-06 | 辻野 次郎丸 | 超音波複合振動加工装置の振動検出器 |
CN109029598A (zh) * | 2018-06-15 | 2018-12-18 | 浙江大学 | 一种互相关法气体超声波流量计中基于主次包络结合处的参考波形选取方法 |
CN109297550A (zh) * | 2018-09-21 | 2019-02-01 | 浙江大学 | 基于互相关函数包络特征点的超声波渡越时间计算方法 |
Non-Patent Citations (4)
Title |
---|
《Automatic Ultrasonic Robotic Array》;Gordon Dobie;《2013 IEEE International Ultrasonics Symposium (IUS)》;20130725;全文 * |
《Multipath ultrasonic gas flow-meter based on multiple reference waves》;Hongliang Zhou;《Ultrasonics》;20170718;全文 * |
《基于温度自适应的超声波渡越时间测量方法研究》;李勇;《自动化与仪表》;20190115;全文 * |
《时差法互相关气体超声波流量计低流速计算及参考波形研究》;汤士桢;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20190815;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111024218A (zh) | 2020-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107015230B (zh) | 一种超声波测距方法 | |
US8427905B2 (en) | Optimum pseudo random sequence determining method, position detection system, position detection method, transmission device and reception device | |
US8363511B2 (en) | Method and device for detection of surroundings | |
Fortuna et al. | Chaotic pulse position modulation to improve the efficiency of sonar sensors | |
CN102073040B (zh) | 超音波感测器的控制方法 | |
US11231493B2 (en) | Methods and apparatus for narrowband ranging systems using coarse and fine delay estimation | |
CN101594185B (zh) | 移动水声通信信号的多普勒估计与同步方法 | |
US9146307B2 (en) | Ultrasonic detection device | |
JP2021524910A (ja) | コード化された超音波送信に対する飛行時間及びコードシグネチャ検出 | |
KR102311559B1 (ko) | 초음파 트랜스듀서의 에코 신호 처리 방법 | |
US20090242310A1 (en) | Object detecting apparatus, opening and closing control system for vehicle using object detecting apparatus, and method of detecting upstroke of envelope | |
EP0102074A2 (en) | Method and apparatus for detecting position of object with ultrasonic wave | |
CN100510775C (zh) | 一种低信噪比下超声波测距的方法及其测距仪 | |
CN111024218B (zh) | 基于自动跟踪的超声波相关检测方法 | |
CN101458332A (zh) | 一种超声波测距方法及其系统 | |
CN105911308A (zh) | 一种风速风向的测量方法 | |
CN108562906B (zh) | 声波测距方法及装置 | |
CN110233677A (zh) | 一种基于光正交码的激光通信测距装置及方法 | |
Huang et al. | High-precision ultrasonic ranging system platform based on peak-detected self-interference technique | |
CN105116371B (zh) | 一种基于连续发射调频信号的目标定位方法与装置 | |
US9103764B2 (en) | Method and system for determining the time-of-flight of an acoustic signal | |
CN110988853B (zh) | 基于有效峰值的超声波到达时间计算方法 | |
JP5103793B2 (ja) | 障害物検出装置および位置特定方法 | |
CN103969648B (zh) | 超声波测距方法 | |
JP2007315892A (ja) | 障害物検出装置および受信時刻推測方法 |
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 |