CN106997476A - 多源无标签数据学习建模的传动系统性能退化评估方法 - Google Patents
多源无标签数据学习建模的传动系统性能退化评估方法 Download PDFInfo
- Publication number
- CN106997476A CN106997476A CN201710119906.9A CN201710119906A CN106997476A CN 106997476 A CN106997476 A CN 106997476A CN 201710119906 A CN201710119906 A CN 201710119906A CN 106997476 A CN106997476 A CN 106997476A
- Authority
- CN
- China
- Prior art keywords
- system performance
- label
- transmission system
- source
- learning model
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Probability & Statistics with Applications (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种多源无标签数据学习建模的传动系统性能退化评估方法。该方法首先将天线传动系统异构信号转化为同构数据,建立同构测试信号与系统退化特征之间的映射关系;然后,建立3层堆叠式深度置信网络拓扑模型,通过对多源无标签信号的编解码重建实现系统退化隐含特征的非监督提取;最后,基于提取的退化隐含特征,计算连续退化特征的局部密度值,并确定无标签条件下系统性能状态的聚类个数,按照距离最大准则进一步计算出各个聚类中心的初始值,利用动态时间弯曲距离计算各个聚类中心的边界,微调聚类中心与聚类边界,完成退化特征向系统性能状态的无标签映射。本发明具有特征提取精度高、状态评估简单有效等特点。
Description
技术领域
本发明涉及机械设备健康评估领域,具体为一种多源无标签数据学习建模的传动系统性能退化评估方法。
背景技术
近十年来机械设备健康评估技术成为长寿命、高可靠机械设备运营管理的关键技术,相比较于传统的设备运营管理技术,健康评估具有可发现早期的系统性能衰退,能给出设备的运行健康程度,可将现行的设备定期维护发展到视情维护,可大幅度降低运营管理费用等优点。而性能退化评估是进行机械设备健康评估的重要前提,因此关于性能退化评估的研究是机械设备健康评估领域的一个关键课题。
当前性能退化评估的方法主要分两类:基于设备健康状态标签的有监督分类和基于监测数据的无监督聚类。基于设备健康状态标签的有监督分类方法,依赖于先验知识构造退化特征,依赖于健康状态标签进行退化状态分类,对于实际中缺乏早期退化知识和退化状态标签的机械设备难以凑效。基于监测数据的无监督聚类方法,由于降低了对退化先验知识和退化状态标签的依赖,在实际中得到了较多的应用,然而现有聚类方法需要事先指定聚类个数,而聚类个数直接决定退化状态个数,聚类个数的先验确定加剧了性能退化状态评估结果的不确定性。如何在多源异构数据下非监督提取退化特征,在无标签条件下确定合理的确定聚类个数并进行聚类是一个有待解决的问题。
发明内容
针对现有技术中存在的问题,本发明提供一种多源无标签数据学习建模的传动系统性能退化评估方法。该方法既克服了退化特征提取对先验知识的依赖,又解决了传统无标签聚类方法需要事先指定聚类个数的问题,具有计算精度高、计算速度快的特点,可实现天线传动系统性能退化状态的有效评估。
为达到以上目的,本发明采取如下技术方案予以实现:
一种多源无标签数据学习建模的传动系统性能退化评估方法,该方法包括以下步骤:
步骤1:对天线传动系统在不同转速、不同采样频率下的异构测试信号通过动态时间规整转化为同构数据;其中异构测试信号包括振动、电流、温度和转角测试信号;
步骤2:利用稀疏约束下的受限玻尔兹曼机模型建立复杂环境条件下多源无标签同构测试信号与系统退化特征之间的映射关系;
步骤3:根据Encoder-Decoder嵌套结构,利用稀疏受限玻尔兹曼机建立系统退化特征提取的6层堆叠式深度置信网络拓扑模型;
步骤4:采用随机梯度下降方法和多源信号进行深度置信网络拓扑模型参数的自适应学习,通过对多源无标签同构测试信号的编解码重建实现系统退化隐含特征的非监督提取;
步骤5:根据步骤4提取的系统退化隐含特征,计算出高维空间下的连续退化特征局部密度值,并确定无标签条件下系统性能状态的聚类个数;
步骤6:根据步骤5确定的系统性能状态聚类个数,按照距离最大准则进一步计算出各个聚类中心的初始值;
步骤7:根据步骤6确定的系统性能状态聚类个数和聚类中心初始值,利用动态时间弯曲距离计算各个聚类中心的边界,将类内距离最小与类间距离最大作为优化目标,微调聚类中心与聚类边界,完成退化特征向系统性能状态的无标签映射。
步骤2中,所述的深度置信网络拓扑模型选择6个受限玻尔兹曼机单元按照堆叠式结构进行组合,每个受限玻尔兹曼机的输出层按照如下式子进行稀疏化约束:
其中Φ为基函数向量,Φ={φ1(x1,x2,x3,x4),φ2(x1,x2,x3,x4),...φn(x1,x2,x3,x4)};x1、x2、x3、x4是振动、电流、温度和转角测试数据;ai是基向量的系数,随着不同的信号而变化。
步骤2中,所述的深度置信网络拓扑模型选择6个受限玻尔兹曼机单元按照堆叠式结构进行组合,其中第一个受限玻尔兹曼机模型输入层为4个单元输出层为10个单元,第二个模型的输入为10个单元输出为5个单元,第三个模型的输入为5个单元输出为1个单元。
步骤3中,建立的3个稀疏受限玻尔兹曼机按照顺序堆叠的方式建立系统退化特征提取的Encoder模型,利用上述3个稀疏受限玻尔兹曼机按照逆序堆叠的方式建立系统退化特征提取的Decoder模型,组合而成的Encoder-Decoder深度置信网络共5个隐含层,Encoder模型和Decoder模型共享一个隐含层,6个受限玻尔兹曼机的输出层均满足稀疏约束条件。
步骤4中,建立的5层深度置信网络拓扑模型共111个网络参数的自适应学习,采用“4-10-5-1-5-10-4”网络对振动、电流、温度、转角等4个测试信号无标签同构测试信号进行编解码重建,从而实现系统退化隐含特征的非监督提取。
步骤5中,计算出的高维空间下的连续退化特征局部密度值,并确定无标签条件下系统性能状态的聚类个数为8。
步骤7中,类内距离最小与类间距离最大按照下式计算:
其中,xi、yi为连续退化特征上任意一点的坐标数据,σi为数据的标准差,p为数据点的个数。
本发明与现有获取超声导波频散曲线方法相比,具有以下优点:
本发明由于对单个受限玻尔兹曼机输出层施加了稀疏约束,相比较于传统特征提取方法,本方法特征提取精度更高,计算速度更快;由于采用4个测试信号进行信息融合,采用6个稀疏受限玻尔兹曼机进行堆叠,与传统深度置信网络相比网络参数求解精度高,尤其在多层网络拓扑结构下更高;由于采用局部密度、距离最大准则和动态时间弯曲距离,本方法可自动确定聚类个数、聚类中心、聚类边界,克服了传统方法聚类中心需要事先指定的问题。该方法利用多层堆叠式稀疏受限玻尔兹曼机和数据密度峰值聚类进行数据深度学习建模,可得到多源异构监测数据下的设备性能退化状态,既克服了退化特征提取对先验知识的依赖,又解决了传统无标签聚类方法(如K均值)需要事先指定聚类个数的问题,具有特征提取精度高、状态评估简单有效等特点,为研究机械设备健康评估中性能退化状态的获取提供一种可行的方法。
附图说明
图1为多源无标签数据深度学习建模的天线传动系统性能退化评估方法流程图;
具体实施方式
为了详细说明本发明的具体实施,以4个信号源(振动、电流、温度、转角)、3种采样频率(5kHz、10kHz、20kHz)、无标退化状态签条件下的船载天线传动系统性能退化状态评估为例进行说明。结合附图及实例是对本发明的解释,并不限制本发明的权利保护范围,如图1所示,具体步骤如下:
1)分别对天线传动系统振动、电流、温度、转角4个测试信号,采用动态时间规整方法,将5kHz、10kHz、20kHz三种频率下的异构信号转化为5kHz条件下的振动、电流、温度、转角同构化数据;
2)利用稀疏约束下的受限玻尔兹曼机模型建立复杂环境条件下4源无标签同构数据与系统退化特征之间的映射关系,其中第一个受限玻尔兹曼机模型输入层为4个单元输出层为10个单元,第二个模型的输入为10个单元输出为5个单元,第三个模型的输入为5个单元输出为1个单元;受限玻尔兹曼机模型的输出层稀疏约束满足下式:
式中Φ为基函数向量,Φ={φ1(x1,x2,x3,x4),φ2(x1,x2,x3,x4),...φn(x1,x2,x3,x4)};x1、x2、x3、x4是振动、电流、温度和转角测试数据;ai是基向量的系数。
3)根据Encoder-Decoder嵌套结构,利用2)中建立的3个稀疏受限玻尔兹曼机按照顺序堆叠的方式建立系统退化特征提取的Encoder模型,利用上述3个稀疏受限玻尔兹曼机按照逆序堆叠的方式建立系统退化特征提取的Decoder模型,组合而成的Encoder-Decoder深度置信网络共5个隐含层,Encoder模型和Decoder模型共享一个隐含层。
4)采用随机梯度下降方法和上述5kHz条件下的振动、电流、温度、转角同构化数据进行5层深度置信网络拓扑模型共111个网络参数的自适应学习,采用建立的“4-10-5-1-5-10-4”网络对4源无标签同构测试信号进行编解码重建,从而实现系统退化隐含特征的非监督提取;
5)根据步骤4提取的系统退化隐含特征,计算出高维空间下的连续退化特征局部密度值,并确定无标签条件下系统性能状态的聚类个数为8;
6)根据步骤5确定的系统性能状态聚类个数8,按照距离最大准则进一步计算出上述8个聚类中心的初始值;
7)根据步骤6确定的8个系统性能状态聚类个数和8个聚类中心初始值,利用动态时间弯曲距离计算各个聚类中心的边界,按照下式将类内距离最小与类间距离最大作为优化目标,并微调上述8个聚类中心初值与聚类边界,完成退化特征向系统性能状态的无标签映射。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明。所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种多源无标签数据学习建模的传动系统性能退化评估方法,其特征在于,该方法包括以下步骤:
步骤1:对天线传动系统在不同转速、不同采样频率下的异构测试信号通过动态时间规整转化为同构数据;其中异构测试信号包括振动、电流、温度和转角测试信号;
步骤2:利用稀疏约束下的受限玻尔兹曼机模型建立复杂环境条件下多源无标签同构测试信号与系统退化特征之间的映射关系;
步骤3:根据Encoder-Decoder嵌套结构,利用稀疏受限玻尔兹曼机建立系统退化特征提取的6层堆叠式深度置信网络拓扑模型;
步骤4:采用随机梯度下降方法和多源信号进行深度置信网络拓扑模型参数的自适应学习,通过对多源无标签同构测试信号的编解码重建实现系统退化隐含特征的非监督提取;
步骤5:根据步骤4提取的系统退化隐含特征,计算出高维空间下的连续退化特征局部密度值,并确定无标签条件下系统性能状态的聚类个数;
步骤6:根据步骤5确定的系统性能状态聚类个数,按照距离最大准则进一步计算出各个聚类中心的初始值;
步骤7:根据步骤6确定的系统性能状态聚类个数和聚类中心初始值,利用动态时间弯曲距离计算各个聚类中心的边界,将类内距离最小与类间距离最大作为优化目标,微调聚类中心与聚类边界,完成退化特征向系统性能状态的无标签映射。
2.根据权利要求1所述的一种多源无标签数据学习建模的传动系统性能退化评估方法,其特征在于,步骤2中,所述的深度置信网络拓扑模型选择6个受限玻尔兹曼机单元按照堆叠式结构进行组合,每个受限玻尔兹曼机的输出层按照如下式子进行稀疏化约束:
其中Φ为基函数向量,Φ={φ1(x1,x2,x3,x4),φ2(x1,x2,x3,x4),...φn(x1,x2,x3,x4)};x1、x2、x3、x4是振动、电流、温度和转角测试数据;ai是基向量的系数,随着不同的信号而变化。
3.根据权利要求2所述的一种多源无标签数据学习建模的传动系统性能退化评估方法,其特征在于,步骤2中,所述的深度置信网络拓扑模型选择6个受限玻尔兹曼机单元按照堆叠式结构进行组合,其中第一个受限玻尔兹曼机模型输入层为4个单元输出层为10个单元,第二个模型的输入为10个单元输出为5个单元,第三个模型的输入为5个单元输出为1个单元。
4.根据权利要求2所述的一种多源无标签数据学习建模的传动系统性能退化评估方法,其特征在于,步骤3中,建立的3个稀疏受限玻尔兹曼机按照顺序堆叠的方式建立系统退化特征提取的Encoder模型,利用上述3个稀疏受限玻尔兹曼机按照逆序堆叠的方式建立系统退化特征提取的Decoder模型,组合而成的Encoder-Decoder深度置信网络共5个隐含层,Encoder模型和Decoder模型共享一个隐含层,6个受限玻尔兹曼机的输出层均满足稀疏约束条件。
5.根据权利要求3所述的一种多源无标签数据学习建模的传动系统性能退化评估方法,其特征在于,步骤4中,建立的5层深度置信网络拓扑模型共111个网络参数的自适应学习,采用“4-10-5-1-5-10-4”网络对振动、电流、温度、转角的4个测试信号无标签同构测试信号进行编解码重建,从而实现系统退化隐含特征的非监督提取。
6.根据权利要求1所述的一种多源无标签数据学习建模的传动系统性能退化评估方法,其特征在于,步骤5中,计算出的高维空间下的连续退化特征局部密度值,并确定无标签条件下系统性能状态的聚类个数为8。
7.根据权利要求1所述的一种多源无标签数据学习建模的传动系统性能退化评估方法,其特征在于,步骤7中,类内距离最小与类间距离最大按照下式计算:
其中,xi、yi为连续退化特征上任意一点的坐标数据,σi为数据的标准差,p为数据点的个数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710119906.9A CN106997476B (zh) | 2017-03-01 | 2017-03-01 | 多源无标签数据学习建模的传动系统性能退化评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710119906.9A CN106997476B (zh) | 2017-03-01 | 2017-03-01 | 多源无标签数据学习建模的传动系统性能退化评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106997476A true CN106997476A (zh) | 2017-08-01 |
CN106997476B CN106997476B (zh) | 2020-04-28 |
Family
ID=59430997
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710119906.9A Active CN106997476B (zh) | 2017-03-01 | 2017-03-01 | 多源无标签数据学习建模的传动系统性能退化评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106997476B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107545112A (zh) * | 2017-09-07 | 2018-01-05 | 西安交通大学 | 多源无标签数据机器学习的复杂装备性能评估与预测方法 |
CN110599416A (zh) * | 2019-09-02 | 2019-12-20 | 太原理工大学 | 一种基于空间目标图像数据库的非合作目标图像盲复原方法 |
CN111090024A (zh) * | 2019-11-14 | 2020-05-01 | 国网上海市电力公司 | 基于外部热、声特征信息的gis状态评估方法及装置 |
CN113588266A (zh) * | 2021-07-19 | 2021-11-02 | 西安交通大学 | 一种故障语义空间内嵌的滚动轴承复合故障诊断方法 |
CN114881917A (zh) * | 2022-03-17 | 2022-08-09 | 深圳大学 | 基于自监督和语义分割的溶栓疗效预测方法及相关装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103413548A (zh) * | 2013-08-16 | 2013-11-27 | 中国科学技术大学 | 一种基于受限玻尔兹曼机的联合频谱建模的声音转换方法 |
CN103729548A (zh) * | 2013-12-18 | 2014-04-16 | 西安交通大学 | 一种机械系统潜在性能退化的间接评估方法 |
CN105206270A (zh) * | 2015-08-20 | 2015-12-30 | 长安大学 | 一种组合pca和rbm的孤立数字语音识别分类系统及方法 |
CN105738107A (zh) * | 2016-01-28 | 2016-07-06 | 西安交通大学 | 基于近似完备可变模式分解的机车轮对轴承故障诊断方法 |
CN105973594A (zh) * | 2016-04-25 | 2016-09-28 | 西北工业大学 | 一种基于连续深度置信网络的滚动轴承故障预测方法 |
-
2017
- 2017-03-01 CN CN201710119906.9A patent/CN106997476B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103413548A (zh) * | 2013-08-16 | 2013-11-27 | 中国科学技术大学 | 一种基于受限玻尔兹曼机的联合频谱建模的声音转换方法 |
CN103729548A (zh) * | 2013-12-18 | 2014-04-16 | 西安交通大学 | 一种机械系统潜在性能退化的间接评估方法 |
CN105206270A (zh) * | 2015-08-20 | 2015-12-30 | 长安大学 | 一种组合pca和rbm的孤立数字语音识别分类系统及方法 |
CN105738107A (zh) * | 2016-01-28 | 2016-07-06 | 西安交通大学 | 基于近似完备可变模式分解的机车轮对轴承故障诊断方法 |
CN105973594A (zh) * | 2016-04-25 | 2016-09-28 | 西北工业大学 | 一种基于连续深度置信网络的滚动轴承故障预测方法 |
Non-Patent Citations (3)
Title |
---|
QING WANG等: "A Universal VAD Based on Jointly Trained Deep Neural Networks", 《网页在线公开:HTTPS://ICCA-SPEECH.ORG/ARCHIVE.INTERSPEECH_2015/I15_2282.HTML》 * |
ZITONG ZHOU等: "Multi-domain description method for bearing fault recognition in varying speed condition", 《IECON2016-42ND ANNUAL CONFERENCE OF THE IEEE INDUSTRIAL ELECTRONICS SOCIETY》 * |
余荣斌: "基于性能退化的光伏组件服役可靠性评估方法研究", 《中国博士学位论文全文数据库 工程科技II辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107545112A (zh) * | 2017-09-07 | 2018-01-05 | 西安交通大学 | 多源无标签数据机器学习的复杂装备性能评估与预测方法 |
CN107545112B (zh) * | 2017-09-07 | 2020-11-10 | 西安交通大学 | 多源无标签数据机器学习的复杂装备性能评估与预测方法 |
CN110599416A (zh) * | 2019-09-02 | 2019-12-20 | 太原理工大学 | 一种基于空间目标图像数据库的非合作目标图像盲复原方法 |
CN110599416B (zh) * | 2019-09-02 | 2022-10-11 | 太原理工大学 | 一种基于空间目标图像数据库的非合作目标图像盲复原方法 |
CN111090024A (zh) * | 2019-11-14 | 2020-05-01 | 国网上海市电力公司 | 基于外部热、声特征信息的gis状态评估方法及装置 |
CN111090024B (zh) * | 2019-11-14 | 2021-11-12 | 国网上海市电力公司 | 基于外部热、声特征信息的gis状态评估方法及装置 |
CN113588266A (zh) * | 2021-07-19 | 2021-11-02 | 西安交通大学 | 一种故障语义空间内嵌的滚动轴承复合故障诊断方法 |
CN114881917A (zh) * | 2022-03-17 | 2022-08-09 | 深圳大学 | 基于自监督和语义分割的溶栓疗效预测方法及相关装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106997476B (zh) | 2020-04-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106997476A (zh) | 多源无标签数据学习建模的传动系统性能退化评估方法 | |
Zhang et al. | [Retracted] Complex System of Vertical Baduanjin Lifting Motion Sensing Recognition under the Background of Big Data | |
Mu et al. | Estimation and inference in spatially varying coefficient models | |
Li et al. | Predicting Short‐Term Electricity Demand by Combining the Advantages of ARMA and XGBoost in Fog Computing Environment | |
Saufi et al. | Differential evolution optimization for resilient stacked sparse autoencoder and its applications on bearing fault diagnosis | |
Shao et al. | Branchy-GNN: A device-edge co-inference framework for efficient point cloud processing | |
CN101832471A (zh) | 一种信号识别分类方法 | |
CN110213003B (zh) | 一种无线信道大尺度衰落建模方法及装置 | |
Marwan et al. | Cross recurrence plots and their applications | |
CN101540047A (zh) | 基于独立高斯混合模型的纹理图像分割方法 | |
CN102842043B (zh) | 基于自动聚类的粒子群优化分类方法 | |
CN104732545A (zh) | 结合稀疏近邻传播和快速谱聚类的纹理图像分割方法 | |
Zaree et al. | Improvement of air pollution prediction in a smart city and its correlation with weather conditions using metrological big data | |
CN106231553B (zh) | 基于无线声传感器网络的多节点信息融合声源定位方法 | |
CN106405683B (zh) | 基于g-l混合噪声特性核岭回归技术的风速预报方法及装置 | |
Guo et al. | Electromagnetic environment portrait based on big data mining | |
CN110838073A (zh) | 一种基于聚类算法的钻井溢流数据异常值处理方法 | |
CN109738852A (zh) | 基于低秩矩阵重建的分布式源二维空间谱估计方法 | |
CN107622476A (zh) | 基于概率生成模型的图像超分辨处理方法 | |
CN103530277A (zh) | 一种基于稀疏感知的多流形邻域点选择方法 | |
CN101465006B (zh) | 由点云构建网格面的方法 | |
Chen et al. | Artificial immune K‐means grid‐density clustering algorithm for real‐time monitoring and analysis of urban traffic | |
Shi et al. | Study on quick judgement of small signal stability using CNN | |
Chaki et al. | A novel framework based on SVDD to classify water saturation from seismic attributes | |
WO2021014495A1 (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 |