CN109031328A - 含时序增益功能激光测距机的探测概率曲线簇测量方法 - Google Patents
含时序增益功能激光测距机的探测概率曲线簇测量方法 Download PDFInfo
- Publication number
- CN109031328A CN109031328A CN201810602572.5A CN201810602572A CN109031328A CN 109031328 A CN109031328 A CN 109031328A CN 201810602572 A CN201810602572 A CN 201810602572A CN 109031328 A CN109031328 A CN 109031328A
- Authority
- CN
- China
- Prior art keywords
- laser
- detection probability
- range finder
- probability
- laser range
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 92
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000012360 testing method Methods 0.000 claims abstract description 28
- 238000005259 measurement Methods 0.000 claims abstract description 15
- 238000002356 laser light scattering Methods 0.000 claims description 22
- 238000000691 measurement method Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000002474 experimental method Methods 0.000 claims description 8
- 239000000523 sample Substances 0.000 claims description 6
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims 1
- 230000035945 sensitivity Effects 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000000149 argon plasma sintering Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/06—Systems determining position data of a target
- G01S17/08—Systems determining position data of a target for measuring distance only
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/497—Means for monitoring or calibrating
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Optical Radar Systems And Details Thereof (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
本发明公开了含时序增益功能激光测距机的探测概率曲线簇测量方法,包括以下几个步骤;第一步:测量激光接收系统不同衰减系数下的探测概率;第二步:确定出激光回波功率的概率分布曲线;第三步:采用“一种激光接收系统探测概率曲线测量解算方法”中的解算方法,解算得到此距离时对应的激光接收系统探测概率的实际响应曲线;第四步:改变距离,重复上述过程,即可得到不同距离的接收系统探测概率的实际响应曲线;直至随着距离的增大,激光接收系统探测概率的实际响应曲线不再改变,测试过程结束;本方法在保证了测试准确性的同时,得到了不同时序增益情况下探测概率与接收功率的完整变化曲线簇。
Description
技术领域
本发明涉及激光设备性能测试领域,具体是含时序增益功能激光测距机的探测概率曲线簇测量方法。
背景技术
激光测距机为有效抑制大气后向散射造成的误测,其接收系统中均设计有时序增益控制功能,在激光发射的瞬间,放大倍率最小,灵敏度阈值较高;而后随时间增加,其放大倍率呈指数形式增加,相应的灵敏度阈值不断降低;在超过一定时间后,其放大倍率达到最大值,相应的灵敏度阈值达到最低,并在之后的时间内保持不变。激光接收系统探测概率与接收功率的变化曲线是描述激光测距机极限工作能力的一种最有效描述手段,最为全面、直观。但现有测试方法都在实验室进行,要求使用稳定性极高的激光发射光源,且因为不能精确控制激光光源发射时间,故只能屏蔽掉激光测距机的时序增益控制功能,而只能测量放大倍率最大(相应的灵敏度阈值最低)时,激光接收系统探测概率与接收功率的变化曲线,而不能测量激光测距机在时序增益过程中的探测概率与接收功率的变化曲线。
发明内容
本发明的目的在于提供一种适用于含时序增益功能激光测距机的探测概率曲线簇测量方法,以解决上述背景技术中提出的问题。
为实现上述目的,本发明提供如下技术方案:
含时序增益功能激光测距机的探测概率曲线簇测量方法,包括以下几个步骤。
第一步:测量激光接收系统不同衰减系数下的探测概率
1-1、将可调衰减器的衰减系数调节设置到最大,保证激光发射系统发出的激光经漫反射靶板散射和衰减器衰减后不能被待测激光测距机的接收系统响应。
1-2、减小可调衰减器的衰减系数,直至待测激光测距机的接收系统可部分响应输出,如待测激光测距机发射N(≥100)次激光,激光接收系统探测响应输出M次,则记此衰减系数下待测激光激光测距机接收系统的探测概率为(100×M/N)%。
1-3、继续减小可调衰减器的衰减系数,重复1-2环节,得到其它衰减系数下待测激光测距机接收系统的探测概率;直至激光接收系统能够完全探测响应到测试光源发出的光(M=N)。记有效测量到的探测概率组数为K,每组对应的衰减系数记为 。
注意:可调衰减器可采用多个定值衰减片组合,也可采用连续可调衰减装置。如采用连续可调衰减装置,则在1-3环节尽量保证探测概率为0至100%之间均匀分布;如采用多个定值衰减片组合方式,则尽量使探测概率接近均匀分布。
第二步:确定出激光回波功率的概率分布曲线
2-1、采用激光散射测量仪同步记录第一步中K组有效测量过程中的K×N次激光散射回波的功率数据。采用现代数学分析软件(如MatLAB)确定这K×N个激光散射回波功率数据的概率分布图。
2-2、在确定的概率分布图上取一定数量的特征值,用这些值在数学分析软件(如MatLAB或DataFit)中计算概率密度函数的结构参数,得到概率密度函数的具体表达式。
第三步:采用“一种激光接收系统探测概率曲线测量解算方法”中的解算方法,即可解算得到此距离时对应的激光接收系统探测概率的实际响应曲线。
第四步:改变距离R的大小,重复上述过程,即可得到不同距离(对应不同激光回波时间)的接收系统探测概率的实际响应曲线;直至随着距离R的增大,激光接收系统探测概率的实际响应曲线不再改变,说明此距离对应的时间后,激光测距机的时序增益功能已达到最大值,测试过程结束。
含时序增益功能激光测距机的接收系统探测概率测量系统,由激光测距机、可调衰减器、漫反射靶板、激光散射测量仪组成,在激光测距机的接收系统前放置可调衰减器,在一定距离R处放置超过激光光斑大小的满反射靶板,在待测激光测距机同一位置放置激光散射测量仪。
与现有技术相比,本发明的有益效果是:通过测试激光测距机对不同距离上的标准靶板的探测概率,同步监测激光接收系统位置处激光回波的功率密度分布,可有效测量激光测距机在时序增益时间段的探测概率与接收功率的变化曲线,在保证测试准确性的同时,可得到含时序增益功能激光测距机完整的探测概率与接收功率的变化曲线簇。相比较现有要求采用高稳定性激光光源,且仅能检测激光接收系统最大增益情况下的探测概率与接收功率的变化曲线,本方法在保证了测试准确性的同时,得到了不同时序增益情况下探测概率与接收功率的完整变化曲线簇。
附图说明
图1为含时序增益功能激光测距机的接收系统探测概率测量系统的测试系统框图。
图2为含时序增益功能激光测距机的探测概率曲线簇测量方法中在350m距离上激光接收系统探测概率的测量曲线图。
图3为含时序增益功能激光测距机的探测概率曲线簇测量方法中在350m距离上测试时激光散射回波所有有效数据的概率分布图。
图4为含时序增益功能激光测距机的探测概率曲线簇测量方法中在350m距离上不同衰减系数下,激光散射回波的概率分布图。
图5为含时序增益功能激光测距机的探测概率曲线簇测量方法中在350m距离上激光接收系统探测概率实际响应曲线的傅立叶逆变换直接得到的数据图。
图6为含时序增益功能激光测距机的探测概率曲线簇测量方法中在350m距离上激光接收系统探测概率的实际响应曲线图。
图7为含时序增益功能激光测距机的探测概率曲线簇测量方法中激光接收系统在不同距离上探测概率的实际响应曲线簇图。
具体实施方式
下面结合具体实施方式对本专利的技术方案作进一步详细地说明。
含时序增益功能激光测距机的探测概率曲线簇测量方法,包括以下几个步骤。
第一步:测量激光接收系统不同衰减系数下的探测概率
1-1、将可调衰减器的衰减系数调节设置到最大,保证激光发射系统发出的激光经漫反射靶板散射和衰减器衰减后不能被待测激光测距机的接收系统响应。
1-2、减小可调衰减器的衰减系数,直至待测激光测距机的接收系统可部分响应输出,如待测激光测距机发射N(≥100)次激光,激光接收系统探测响应输出M次,则记此衰减系数下待测激光激光测距机接收系统的探测概率为(100×M/N)%。
1-3、继续减小可调衰减器的衰减系数,重复1-2环节,得到其它衰减系数下待测激光测距机接收系统的探测概率;直至激光接收系统能够完全探测响应到激光发射系统发出的光(M=N)。记有效测量到的探测概率组数为K,每组对应的衰减系数记为。
注意:可调衰减器可采用多个定值衰减片组合,也可采用连续可调衰减装置。如采用连续可调衰减装置,则在1-3环节尽量保证探测概率为0至100%之间均匀分布;如采用多个定值衰减片组合方式,则尽量使探测概率接近均匀分布。
第二步:确定出激光回波功率的概率分布曲线
2-1、采用激光散射测量仪同步记录第一步中K组有效测量过程中的K×N次激光散射回波的功率数据。采用现代数学分析软件(如MatLAB)确定这K×N个激光散射回波功率数据的概率分布图。
2-2、在确定的概率分布图上取一定数量的特征值,用这些值在数学分析软件(如MatLAB或DataFit)中计算概率密度函数的结构参数,得到概率密度函数的具体表达式。
第三步:采用“一种激光接收系统探测概率曲线测量解算方法”中的解算方法,即可解算得到此距离时对应的激光接收系统探测概率的实际响应曲线。
第四步:改变距离R的大小,重复上述过程,即可得到不同距离(对应不同激光回波时间)的接收系统探测概率的实际响应曲线;直至随着距离R的增大,激光接收系统探测概率的实际响应曲线不再改变,说明此距离对应的时间后,激光测距机的时序增益功能已达到最大值,测试过程结束。
“一种激光接收系统探测概率曲线测量解算方法”中的解算方法,其特征在于包括以下步骤:
3-1、第一步中K组有效测量过程中的K×N次激光散射回波的功率的平均值记为,则第i组探测概率时对应的激光散射回波平均功率可按照下式计算:
则可得到K组激光接收系统探测概率对应的激光散射回波平均功率。
3-2、对3-1环节得到的K组探测概率与激光散射回波平均功率数据进行曲线拟合,即可得到此距离下激光接收系统探测概率的当前测量曲线。记为。
3-3、采用现代数学分析软件(如MatLAB),对激光散射回波功率的概率分布曲线和激光接收系统探测概率的当前测量曲线分别进行傅立叶变换得到和。
3-4、在这种测量方式下,激光接收系统探测概率的当前测量曲线在数学上等效为实际响应曲线与激光散射回波功率的概率分布曲线的相关。故可按照下式计算,得到探测概率实际响应曲线的傅立叶变换结果。
3-5、对3-4环节的计算结果进行傅立叶逆变换,即可得到激光接收系统探测概率的实际响应曲线。
。
如图1所示,含时序增益功能激光测距机的接收系统探测概率测量系统,由激光测距机、可调衰减器、漫反射靶板、激光散射测量仪组成,在激光测距机的接收系统前放置可调衰减器,在一定距离R处放置超过激光光斑大小的满反射靶板,在待测激光测距机同一位置放置激光散射测量仪。
对于某型激光测距机的探测概率响应曲线进行测量,漫反射靶板放置在350m时,不同衰减系数对应的探测概率如表1所示。
表1不同衰减系数对应的探测概率测量值(350m)
衰减系数 | 10 | 10.5 | 10.9 | 11.1 | 11.4 | 11.6 | 12 | 12.6 | 13 |
探测概率 | 1.0 | 1.0 | 0.95 | 0.9 | 0.85 | 0.5 | 0.3 | 0.1 | 0 |
在350m距离上,激光接收系统探测概率的测量曲线如图2所示。
在350m距离上,测试过程激光散射回波所有有效数据的概率分布如图3所示。
在350m距离上,不同衰减系数下,激光散射回波的概率分布如图4所示。
在350m距离上,傅立叶逆变换直接得到的数据如图5所示。
对图5进行去噪和曲线拟合,得到350m处,激光接收系统探测概率的实际响应曲线如图6所示。
不同距离处,激光接收系统探测概率的实际响应曲线簇如图7所示。
上面对本专利的较佳实施方式作了详细说明,但是本专利并不限于上述实施方式,在本领域的普通技术人员所具备的知识范围内,还可以在不脱离本专利宗旨的前提下做出各种变化。
Claims (3)
1.含时序增益功能激光测距机的探测概率曲线簇测量方法,其特征在于,包括以下几个步骤:
第一步:测量激光接收系统不同衰减系数下的探测概率
1-1、将可调衰减器的衰减系数调节设置到最大,保证激光发射系统发出的激光经漫反射靶板散射和衰减器衰减后不能被待测激光测距机的接收系统响应;
1-2、减小可调衰减器的衰减系数,直至待测激光测距机的接收系统可部分响应输出,如待测激光测距机发射N(≥100)次激光,激光接收系统探测响应输出M次,则记此衰减系数下待测激光激光测距机接收系统的探测概率为(100×M/N)%;
1-3、继续减小可调衰减器的衰减系数,重复1-2环节,得到其它衰减系数下待测激光测距机接收系统的探测概率;直至激光接收系统能够完全探测响应到激光发射系统发出的光(M=N);
记有效测量到的探测概率组数为K,每组对应的衰减系数记为;
注意:可调衰减器可采用多个定值衰减片组合,也可采用连续可调衰减装置;
如采用连续可调衰减装置,则在1-3环节尽量保证探测概率为0至100%之间均匀分布;如采用多个定值衰减片组合方式,则尽量使探测概率接近均匀分布;
第二步:确定出激光回波功率的概率分布曲线
2-1、采用激光散射测量仪同步记录第一步中K组有效测量过程中的K×N次激光散射回波的功率数据;
采用现代数学分析软件(如MatLAB)确定这K×N个激光散射回波功率数据的概率分布图;
2-2、在确定的概率分布图上取一定数量的特征值,用这些值在数学分析软件(如MatLAB或DataFit)中计算概率密度函数的结构参数,得到概率密度函数的具体表达式;
第三步:采用“一种激光接收系统探测概率曲线测量解算方法”中的解算方法,即可解算得到此距离时对应的激光接收系统探测概率的实际响应曲线;
第四步:改变距离R的大小,重复上述过程,即可得到不同距离(对应不同激光回波时间)的接收系统探测概率的实际响应曲线;直至随着距离R的增大,激光接收系统探测概率的实际响应曲线不再改变,说明此距离对应的时间后,激光测距机的时序增益功能已达到最大值,测试过程结束。
2.根据权利要求1用于含时序增益功能激光测距机的探测概率曲线测量方法的含时序增益功能激光测距机的接收系统探测概率测量系统,其特征在于,由激光测距机、可调衰减器、漫反射靶板、激光散射测量仪组成,在激光测距机的接收系统前放置可调衰减器,在一定距离R处放置超过激光光斑大小的满反射靶板,在待测激光测距机同一位置放置激光散射测量仪。
3.根据权利要求1含时序增益功能激光测距机的探测概率曲线测量方法中“一种激光接收系统探测概率曲线测量解算方法”中的解算方法,其特征在于,包括以下步骤:
3-1、第一步中K组有效测量过程中的K×N次激光散射回波的功率的平均值记为,则第i组探测概率时对应的激光散射回波平均功率可按照下式计算:
则可得到K组激光接收系统探测概率对应的激光散射回波平均功率;
3-2、对3-1环节得到的K组探测概率与激光散射回波平均功率数据进行曲线拟合,即可得到此距离下激光接收系统探测概率的当前测量曲线;
记为;
3-3、采用现代数学分析软件(如MatLAB),对激光散射回波功率的概率分布曲线和激光接收系统探测概率的当前测量曲线分别进行傅立叶变换得到和;
3-4、在这种测量方式下,激光接收系统探测概率的当前测量曲线在数学上等效为实际响应曲线与激光散射回波功率的概率分布曲线的相关;
故可按照下式计算,得到探测概率实际响应曲线的傅立叶变换结果;
3-5、对3-4环节的计算结果进行傅立叶逆变换,即可得到激光接收系统探测概率的实际响应曲线;
。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810602572.5A CN109031328B (zh) | 2018-06-12 | 2018-06-12 | 含时序增益功能激光测距机的探测概率曲线簇测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810602572.5A CN109031328B (zh) | 2018-06-12 | 2018-06-12 | 含时序增益功能激光测距机的探测概率曲线簇测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109031328A true CN109031328A (zh) | 2018-12-18 |
CN109031328B CN109031328B (zh) | 2022-07-26 |
Family
ID=64612836
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810602572.5A Active CN109031328B (zh) | 2018-06-12 | 2018-06-12 | 含时序增益功能激光测距机的探测概率曲线簇测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109031328B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112888959A (zh) * | 2019-09-30 | 2021-06-01 | 深圳市大疆创新科技有限公司 | 激光测距机的量程的测量方法及系统和存储介质 |
WO2021134410A1 (zh) * | 2019-12-31 | 2021-07-08 | 深圳市大疆创新科技有限公司 | 激光测距机的量程的测量方法及系统和存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102324982A (zh) * | 2011-07-21 | 2012-01-18 | 中国科学院上海技术物理研究所 | 空间量子通信单光子远场分布检测装置 |
KR101487411B1 (ko) * | 2013-09-02 | 2015-01-29 | 주식회사 엘지실트론 | 라이너 및 이를 포함하는 에피텍셜 반응기 |
CN105137450A (zh) * | 2015-08-10 | 2015-12-09 | 哈尔滨工业大学 | 低虚警双Gm-APD探测器光子计数激光雷达 |
CN107918119A (zh) * | 2017-10-23 | 2018-04-17 | 孝感华中精密仪器有限公司 | 一种激光报警器灵敏度测试系统 |
-
2018
- 2018-06-12 CN CN201810602572.5A patent/CN109031328B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102324982A (zh) * | 2011-07-21 | 2012-01-18 | 中国科学院上海技术物理研究所 | 空间量子通信单光子远场分布检测装置 |
KR101487411B1 (ko) * | 2013-09-02 | 2015-01-29 | 주식회사 엘지실트론 | 라이너 및 이를 포함하는 에피텍셜 반응기 |
CN105137450A (zh) * | 2015-08-10 | 2015-12-09 | 哈尔滨工业大学 | 低虚警双Gm-APD探测器光子计数激光雷达 |
CN107918119A (zh) * | 2017-10-23 | 2018-04-17 | 孝感华中精密仪器有限公司 | 一种激光报警器灵敏度测试系统 |
Non-Patent Citations (1)
Title |
---|
马娜等: "基于探测概率曲线的探测灵敏度测试方法", 《激光与红外》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112888959A (zh) * | 2019-09-30 | 2021-06-01 | 深圳市大疆创新科技有限公司 | 激光测距机的量程的测量方法及系统和存储介质 |
WO2021134410A1 (zh) * | 2019-12-31 | 2021-07-08 | 深圳市大疆创新科技有限公司 | 激光测距机的量程的测量方法及系统和存储介质 |
CN113330295A (zh) * | 2019-12-31 | 2021-08-31 | 深圳市大疆创新科技有限公司 | 激光测距机的量程的测量方法及系统和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109031328B (zh) | 2022-07-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11592553B2 (en) | Distance measurement system and method using lidar waveform matching | |
CN106226761B (zh) | 一种高性能相干高频雷达多频探测方法 | |
CN104199026B (zh) | 基于线迹扫描二维近场成像的反向散射截面测量方法 | |
CN109238436B (zh) | 混响水池中瞬态声源特性的测量方法 | |
CN102819008A (zh) | 一种基于非线性最小二乘的非合作雷达辐射源定位方法 | |
RU2529441C1 (ru) | Способ обработки гидролокационной информации | |
US20080307025A1 (en) | Location determination in sensor networks using time-window algorithm | |
CN109031328A (zh) | 含时序增益功能激光测距机的探测概率曲线簇测量方法 | |
CN109917344A (zh) | 一种紧缩场内角反射器rcs单天线测量系统及测量方法 | |
CN108152836A (zh) | 一种基于fda天线抗三星频差定位体制的欺骗方法 | |
CN101907708A (zh) | 目标回波亮点测量方法 | |
RU2550576C1 (ru) | Способ измерения дистанции до шумящего объекта | |
Scafidi et al. | Robust picking and accurate location with RSNI‐Picker2: Real‐time automatic monitoring of earthquakes and nontectonic events | |
RU2690223C1 (ru) | Способ определения координат морской шумящей цели | |
RU2334244C1 (ru) | Способ определения местоположения источника радиоизлучения | |
CN108801596A (zh) | 一种激光接收系统探测概率曲线测量解算方法 | |
RU2193782C2 (ru) | Способ оценки характеристик радиолокационной станции при действии активных шумовых помех | |
CN106707251A (zh) | 应答机功率校准方法及装置 | |
CN113406580B (zh) | Mimo雷达低截获性能评估方法、装置和存储介质 | |
RU2768011C1 (ru) | Способ одноэтапного адаптивного определения координат источников радиоизлучений | |
RU2559310C2 (ru) | Способ оценки дистанции до шумящего в море объекта | |
KR20140070723A (ko) | Em 존의 다중 잡음 환경을 고려한 확률적 간섭 평가 장치 및 방법 | |
Li et al. | Pulse train detection algorithm based on edge enhancement | |
RU2378660C1 (ru) | Способ обнаружения и определения координат источника радиоизлучения | |
CN113820675B (zh) | 一种雷达测试方法、装置、设备及存储介质 |
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 |