CN108332776B - Mems陀螺随机误差组合预测模型的构建方法 - Google Patents

Mems陀螺随机误差组合预测模型的构建方法 Download PDF

Info

Publication number
CN108332776B
CN108332776B CN201810240778.8A CN201810240778A CN108332776B CN 108332776 B CN108332776 B CN 108332776B CN 201810240778 A CN201810240778 A CN 201810240778A CN 108332776 B CN108332776 B CN 108332776B
Authority
CN
China
Prior art keywords
random error
angular rate
gyro
prediction model
output angular
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
Application number
CN201810240778.8A
Other languages
English (en)
Other versions
CN108332776A (zh
Inventor
刘亦男
黄玲玲
周卫华
雷小江
余常恒
王文政
刘杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hubei Sanjiang Aerospace Wanfeng Technology Development Co Ltd
Original Assignee
Hubei Sanjiang Aerospace Wanfeng Technology Development Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hubei Sanjiang Aerospace Wanfeng Technology Development Co Ltd filed Critical Hubei Sanjiang Aerospace Wanfeng Technology Development Co Ltd
Priority to CN201810240778.8A priority Critical patent/CN108332776B/zh
Publication of CN108332776A publication Critical patent/CN108332776A/zh
Application granted granted Critical
Publication of CN108332776B publication Critical patent/CN108332776B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
    • G01C25/005Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices

Landscapes

  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Gyroscopes (AREA)

Abstract

本发明公开了一种MEMS陀螺随机误差组合预测模型的构建方法,包括如下步骤:1、采集MEMS陀螺的输出,进行确定性误差标定,并构建陀螺输出数据集;2、对陀螺输出数据集中的数据进行数据预处理;3、获取陀螺的随机误差,并构建陀螺随机误差数据集;4、建立基于指数平滑算法的随机误差预测模型;5、建立基于灰色理论的随机误差预测模型;6、建立组合预测模型。因此该方法能对MEMS陀螺随机误差性进行有效预测,从而提高陀螺精度。

Description

MEMS陀螺随机误差组合预测模型的构建方法
技术领域
本发明涉及陀螺仪测试技术领域,具体地指一种MEMS(Micro-Electro-Mechanical System,微机电系统)陀螺随机误差组合预测模型的构建方法。
背景技术
近些年,随着微电子加工技术的发展,MEMS惯性传感器特别是陀螺仪的精度已得到大幅提高,MEMS陀螺凭借其在体积、成本、功耗和抗冲击能力等方面的优势,已在惯性导航领域有了广泛的应用,如无人机控制导航技术,机器人运动控制,具有动作识别的可穿戴设备等。基于MEMS陀螺的惯性系统己成为当今惯性技术领域的一个重要的研究热点。
然而低成本的MEMS陀螺虽然瞬态性能较好,但误差漂移会随着时间逐渐累积,对测量精度影响很大。如何提高MEMS陀螺的精度一直以来是研究人员们比较关注的一个重要问题。MEMS陀螺的误差包括确定性误差和随机误差,即陀螺的漂移。对确定性误差,可以比较容易的用代数方程来表示,并通过标定补偿。而随机误差由于制造工艺及使用环境的影响,存在很大的不确定性,是影响其精度的一个重要原因。因此,对MEMS陀螺进行随机误差补偿是提高其精度的一个重要手段。
据不完全统计,现有预测方法已经超过300种,但没有一种方法可以至始至终的适应每一种情况,即使是一开始预测效果较好的方法,随着时间变化,预测效果也会变差,反之,一开始预测效果不好的方法之后预测效果也可能变好。单一预测方法的适应范围有限,难以全面反映MEMS陀螺随机漂移的变化情况。
尽管MEMS陀螺随机漂移数据具有不确定性,但有些时段随机误差总体相对稳定,不确定性程度低,而有些时段随机误差不确定程度相对较高。
发明内容
本发明的目的就是要提供一种MEMS陀螺随机误差组合预测模型的构建方法,该组合预测模型能保证MEMS陀螺随机误差预测的准确性。
为实现此目的,本发明所设计的MEMS陀螺随机误差组合预测模型的构建方法,其特征在于,它包括如下步骤:
步骤1:在静态下,以固定的采样速率采集MEMS陀螺输出的角速率数据,对采集到的陀螺输出角速率数据进行确定性误差标定,根据上述确定性误差标定的结果按照时间序列将陀螺输出角速率数据顺序排列组成陀螺输出角速率数据集;
步骤2:对陀螺输出角速率数据集中的数据依据MEMS陀螺的输出特性,进行奇点去除、趋势项去除和周期项去除的预处理;
步骤3:获取陀螺的随机误差,并构建陀螺输出角速率随机误差数据集;
步骤4:根据陀螺输出角速率随机误差数据集建立基于指数平滑算法的陀螺输出角速率随机误差预测模型;
步骤5:建立基于灰色理论的陀螺输出角速率随机误差预测模型;
步骤6:根据基于指数平滑算法的陀螺输出角速率随机误差预测模型和基于灰色理论的陀螺输出角速率随机误差预测模型的预测结果建立陀螺随机误差组合预测模型。
本发明的有益效果:
本发明提供的方法基于组合的思想,考虑了MEMS陀螺随机误差的特性,分别建立了针对数据随机性小的基于指数平滑算法的随机误差预测模型和针对数据随机性大的基于灰色理论的随机误差预测模型,综合利用两种不同预测方法所提供的信息,以获得一个较好的预测结果。因此该方法能对MEMS陀螺随机误差性进行有效预测,从而提高陀螺精度。
附图说明
图1为本发明的流程图;
图2为本发明基于指数平滑算法的随机误差预测模型的建模方法流程图;
图3为本发明基于灰色理论的随机误差预测模型的建模方法流程图。
具体实施方式
以下结合附图和具体实施例对本发明作进一步的详细说明:
本发明所设计的一种MEMS陀螺随机误差组合预测模型的构建方法,如图1所示,它包括如下步骤:
步骤1:在静态下,以固定的采样速率采集MEMS陀螺输出的角速率数据,对采集到的陀螺输出角速率数据进行确定性误差标定,根据上述确定性误差标定的结果按照时间序列将陀螺输出角速率数据顺序排列组成陀螺输出角速率数据集;
步骤2:对陀螺输出角速率数据集中的数据依据MEMS陀螺的输出特性,进行奇点去除、趋势项去除和周期项去除的预处理;
步骤3:获取陀螺的随机误差,并构建陀螺输出角速率随机误差数据集;
步骤4:根据陀螺输出角速率随机误差数据集建立基于指数平滑算法的陀螺输出角速率随机误差预测模型;
步骤5:建立基于灰色理论的陀螺输出角速率随机误差预测模型;
步骤6:根据基于指数平滑算法的陀螺输出角速率随机误差预测模型和基于灰色理论的陀螺输出角速率随机误差预测模型的预测结果建立陀螺随机误差组合预测模型。
上述技术方案的步骤3中陀螺输出角速率随机误差数据集的构建方法为:首先,计算陀螺输出角速率数据集的均值;然后,取出陀螺输出角速率数据集中的每一个数据,减去均值,即得到陀螺输出角速率的随机误差;最后,将得到的陀螺输出角速率的随机误差构建成陀螺输出角速率随机误差数据集。
上述技术方案的步骤4中根据陀螺输出角速率随机误差数据集建立基于指数平滑算法的陀螺输出角速率随机误差预测模型的具体方法包括以下步骤,如图2所示:
步骤401:建立陀螺输出角速率随机误差数据的一次指数平滑模型,模型公式如下:
Figure GDA0002303511300000041
式中,t为时间序号,xt-i为第t-i期陀螺输出角速率的实际值,i的取值为0-(t-1),
Figure GDA0002303511300000042
为陀螺输出角速率预测初始值,
Figure GDA0002303511300000043
为第t+1期的陀螺输出角速率预测值,α为一次指数平滑模型的平滑系数,0<α<1;,平滑系数α的确定方法为先给出一个大致取值范围,算法中有取值规则为通用规则,从中试选各组平滑系数α,经过反复实验比较预测误差来确定使预测误差最小的平滑系数α;
步骤402:确定陀螺输出角速率预测初始值
Figure GDA0002303511300000044
其中
Figure GDA0002303511300000045
上述技术方案的步骤5中,所述步骤5中选用灰色模型GM(1,1)作为基于灰色理论的陀螺输出角速率随机误差预测模型GM(M,N),其中M表征灰色模型中微分方程的阶数,N表征灰色模型中的变量个数,按照以下步骤建模,具体操作步骤如下:
步骤501:建立GM(1,1)模型,设前k期陀螺随机误差的原始数据序列X(0)(k)为
X(0)(k)={x(0)(1),x(0)(2),…,x(0)(k)},其中,k>0;
步骤502:采用一次累加生成法,生成下列新的数据序列X(1)(k):
X(1)(k)={x(1)(1),x(1)(2),…,x(1)(k)}
式中,
Figure GDA0002303511300000051
步骤503:设有x(0)(n)+ax(1)(n)=u,式中a、u为方程的系数,n=1,2,…,k
步骤504:将步骤503的公式表示为离散形式的一阶线性微分方程:
Figure GDA0002303511300000052
即表示时间序列X(1)(k)中的元素对于时间t'的一阶微分方程是连续的;
步骤505:则对于连续时间序列步骤504中对于微分方程的解
Figure GDA0002303511300000053
为:
Figure GDA0002303511300000054
对于离散值时间序列步骤504中对于微分方程的解
Figure GDA0002303511300000055
为:
Figure GDA0002303511300000056
e为自然常数;
步骤506:按以下公式求解上述方程的系数a、u:
[a,u]T=(BTB)-1BTY;式中
Figure GDA0002303511300000057
z(1)(n)=ωx(1)(n-1)+(1-ω)x(1)(n),n=2,3,…,k,通常情况下,ω为常数,ω=0.5,T为转置符号;
步骤507:按以下公式求解原始数据序列X(0)(k)中的预测值
Figure GDA0002303511300000061
Figure GDA0002303511300000062
步骤508:设残差
Figure GDA0002303511300000063
表示原始数据序列的所有预测结果,n为大于1的自然数,计算每次预测的残差,构成残差序列。
上述技术方案中,采取后验差检验方式验证基于灰色理论的陀螺输出角速率随机误差预测模型GM(M,N)的可靠性,后验差检验方式中后验差比值C和小误差概率P分别按以下公式计算:
C=Se/S0
Figure GDA0002303511300000064
式中S0和Se分别为原始数据序列的均方差和残差序列的均方差,
Figure GDA0002303511300000065
为残差序列的均值。
通常要求后验差比值越小越好,表明预测误差的离散性小,而小误差概率越大越好,表明模型的预测精度较高。
上述技术方案的步骤6中,陀螺随机误差组合预测模型如下,如图3所示:
f=ω1f12f2;式中f为陀螺随机误差组合预测模型预测结果;f1和f2分别为基于指数平滑算法的陀螺输出角速率随机误差预测模型和基于灰色理论的陀螺输出角速率随机误差预测模型的预测结果;ω1、ω2分别为基于指数平滑算法的陀螺输出角速率随机误差预测模型和基于灰色理论的陀螺输出角速率随机误差预测模型的权值。
所述ω1和ω2按以下公式确定:
Figure GDA0002303511300000066
式中,E1为基于指数平滑算法的陀螺输出角速率随机误差预测模型中预测误差绝对值的和,E2为基于灰色理论的陀螺输出角速率随机误差预测模型中预测误差绝对值的和。
为了适应MEMS陀螺随机漂移的变化特点,本发明基于组合的思想,同时考虑随机漂移数据的平稳性和随机性。首先针对陀螺随机漂移数据具有相对稳定的特征,由于指数平滑算法能够反映时间序列的平均趋势及其变化,建立基于指数平滑算法的随机误差预测模型;然后针对陀螺随机漂移数据具有随机变化的特征,由于灰色模型可以弱化原始数据序列的随机性,使其呈现出明显的规律性,建立基于灰色理论的随机误差预测模型;最后通过适当的方式对单一预测方法的结果进行组合,从而综合利用不同预测方法所提供的信息,以获得一个较好的预测结果。
本说明书未作详细描述的内容属于本领域专业技术人员公知的现有技术。

Claims (7)

1.一种MEMS陀螺随机误差组合预测模型的构建方法,其特征在于,它包括如下步骤:
步骤1:在静态下,以固定的采样速率采集MEMS陀螺输出的角速率数据,对采集到的陀螺输出角速率数据进行确定性误差标定,根据上述确定性误差标定的结果按照时间序列将陀螺输出角速率数据顺序排列组成陀螺输出角速率数据集;
步骤2:对陀螺输出角速率数据集中的数据依据MEMS陀螺的输出特性,进行奇点去除、趋势项去除和周期项去除的预处理;
步骤3:获取陀螺的随机误差,并构建陀螺输出角速率随机误差数据集;
步骤4:根据陀螺输出角速率随机误差数据集建立基于指数平滑算法的陀螺输出角速率随机误差预测模型;
步骤5:建立基于灰色理论的陀螺输出角速率随机误差预测模型;
步骤6:根据基于指数平滑算法的陀螺输出角速率随机误差预测模型和基于灰色理论的陀螺输出角速率随机误差预测模型的预测结果建立陀螺随机误差组合预测模型。
2.根据权利要求1所述的MEMS陀螺随机误差组合预测模型的构建方法,其特征在于:所述步骤3中陀螺输出角速率随机误差数据集的构建方法为:首先,计算陀螺输出角速率数据集的均值;然后,取出陀螺输出角速率数据集中的每一个数据,减去均值,即得到陀螺输出角速率的随机误差;最后,将得到的陀螺输出角速率的随机误差构建成陀螺输出角速率随机误差数据集。
3.根据权利要求1所述的MEMS陀螺随机误差组合预测模型的构建方法,其特征在于:所述步骤4中根据陀螺输出角速率随机误差数据集建立基于指数平滑算法的陀螺输出角速率随机误差预测模型的具体方法包括以下步骤:
步骤401:建立陀螺输出角速率随机误差数据的一次指数平滑模型,模型公式如下:
Figure FDA0002303511290000021
式中,t为时间序号,xt-i为第t-i期陀螺输出角速率的实际值,i的取值为0-(t-1),
Figure FDA0002303511290000022
为陀螺输出角速率预测初始值,
Figure FDA0002303511290000023
为第t+1期的陀螺输出角速率预测值,α为一次指数平滑模型的平滑系数,0<α<1;
步骤402:确定陀螺输出角速率预测初始值
Figure FDA0002303511290000024
其中
Figure FDA0002303511290000025
4.根据权利要求1所述的MEMS陀螺随机误差组合预测模型的构建方法,其特征在于:所述步骤5中选用灰色模型GM(1,1)作为基于灰色理论的陀螺输出角速率随机误差预测模型GM(M,N),其中M表征灰色模型中微分方程的阶数,N表征灰色模型中的变量个数,按照以下步骤建模,具体操作步骤如下:
步骤501:建立GM(1,1)模型,设前k期陀螺随机误差的原始数据序列X(0)(k)为
X(0)(k)={x(0)(1),x(0)(2),…,x(0)(k)},其中,k>0;
步骤502:采用一次累加生成法,生成下列新的数据序列X(1)(k):
X(1)(k)={x(1)(1),x(1)(2),…,x(1)(k)}
式中,
Figure FDA0002303511290000026
n=1,2,…,k
步骤503:设有x(0)(n)+ax(1)(n)=u,式中a、u为方程的系数,n=1,2,…,k
步骤504:将步骤503的公式表示为离散形式的一阶线性微分方程:
Figure FDA0002303511290000031
即表示时间序列X(1)(k)中的元素对于时间t'的一阶微分方程是连续的;
步骤505:则对于连续时间序列步骤504中对于微分方程的解
Figure FDA0002303511290000032
为:
Figure FDA0002303511290000033
对于离散值时间序列步骤504中对于微分方程的解
Figure FDA0002303511290000034
为:
Figure FDA0002303511290000035
e为自然常数;
步骤506:按以下公式求解上述方程的系数a、u:
[a,u]T=(BTB)-1BTY;式中
Figure FDA0002303511290000036
z(1)(n)=ωx(1)(n-1)+(1-ω)x(1)(n),n=2,3,…,k,通常情况下,ω为常数,ω=0.5,T为转置符号;
步骤507:按以下公式求解原始数据序列X(0)(k)中的预测值
Figure FDA0002303511290000037
Figure FDA0002303511290000038
步骤508:设残差
Figure FDA0002303511290000039
Figure FDA00023035112900000310
表示原始数据序列的所有预测结果,n为大于1的自然数,计算每次预测的残差,构成残差序列。
5.根据权利要求4所述的MEMS陀螺随机误差组合预测模型的构建方法,其特征在于:采取后验差检验方式验证基于灰色理论的陀螺输出角速率随机误差预测模型GM(M,N)的可靠性,后验差检验方式中后验差比值C和小误差概率P分别按以下公式计算:
C=Se/S0
Figure FDA0002303511290000041
式中S0和Se分别为原始数据序列的均方差和残差序列的均方差,
Figure FDA0002303511290000042
为残差序列的均值。
6.根据权利要求4所述的MEMS陀螺随机误差组合预测模型的构建方法,其特征在于:步骤6中,陀螺随机误差组合预测模型如下:
f=ω1f12f2;式中f为陀螺随机误差组合预测模型预测结果;f1和f2分别为基于指数平滑算法的陀螺输出角速率随机误差预测模型和基于灰色理论的陀螺输出角速率随机误差预测模型的预测结果;ω1、ω2分别为基于指数平滑算法的陀螺输出角速率随机误差预测模型和基于灰色理论的陀螺输出角速率随机误差预测模型的权值。
7.根据权利要求6所述的MEMS陀螺随机误差组合预测模型的构建方法,其特征在于:所述ω1和ω2按以下公式确定:
Figure FDA0002303511290000043
式中,E1为基于指数平滑算法的陀螺输出角速率随机误差预测模型中预测误差绝对值的和,E2为基于灰色理论的陀螺输出角速率随机误差预测模型中预测误差绝对值的和。
CN201810240778.8A 2018-03-22 2018-03-22 Mems陀螺随机误差组合预测模型的构建方法 Active CN108332776B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810240778.8A CN108332776B (zh) 2018-03-22 2018-03-22 Mems陀螺随机误差组合预测模型的构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810240778.8A CN108332776B (zh) 2018-03-22 2018-03-22 Mems陀螺随机误差组合预测模型的构建方法

Publications (2)

Publication Number Publication Date
CN108332776A CN108332776A (zh) 2018-07-27
CN108332776B true CN108332776B (zh) 2020-06-02

Family

ID=62931441

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810240778.8A Active CN108332776B (zh) 2018-03-22 2018-03-22 Mems陀螺随机误差组合预测模型的构建方法

Country Status (1)

Country Link
CN (1) CN108332776B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0685705A1 (en) * 1994-06-02 1995-12-06 Matsushita Electric Industrial Co., Ltd. Offset-drift correcting device for gyro-sensor
CN101666662A (zh) * 2009-09-25 2010-03-10 北京航空航天大学 基于模糊理论的加速退化试验预测方法
CN103900610A (zh) * 2014-03-28 2014-07-02 哈尔滨工程大学 基于灰色小波神经网络的mems陀螺随机误差预测方法
CN105675017A (zh) * 2016-01-12 2016-06-15 山东理工大学 一种应用于光电平台的光纤陀螺随机漂移补偿方法
CN107330149A (zh) * 2017-05-27 2017-11-07 哈尔滨工业大学 基于arma和bpnn组合模型的mimu陀螺随机漂移预测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0685705A1 (en) * 1994-06-02 1995-12-06 Matsushita Electric Industrial Co., Ltd. Offset-drift correcting device for gyro-sensor
CN101666662A (zh) * 2009-09-25 2010-03-10 北京航空航天大学 基于模糊理论的加速退化试验预测方法
CN103900610A (zh) * 2014-03-28 2014-07-02 哈尔滨工程大学 基于灰色小波神经网络的mems陀螺随机误差预测方法
CN105675017A (zh) * 2016-01-12 2016-06-15 山东理工大学 一种应用于光电平台的光纤陀螺随机漂移补偿方法
CN107330149A (zh) * 2017-05-27 2017-11-07 哈尔滨工业大学 基于arma和bpnn组合模型的mimu陀螺随机漂移预测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Hybrid Grey Modification Model with a Fourier Series and Exponential Smoothing Optimization Technique;Yu-Tung Lin et al.;《The Journal of Grey System》;20141231;第26卷(第4期);第108-121页 *
微惯性传感器零点漂移融合滤波算法研究;李雪莲等;《传感器与微系统》;20091231;第28卷(第6期);第50-52页 *

Also Published As

Publication number Publication date
CN108332776A (zh) 2018-07-27

Similar Documents

Publication Publication Date Title
CN109472418B (zh) 基于卡尔曼滤波的机动目标状态预测优化方法
CN111156987A (zh) 基于残差补偿多速率ckf的惯性/天文组合导航方法
Sushchenko et al. Processing of redundant information in airborne electronic systems by means of neural networks
CN106767780B (zh) 基于Chebyshev多项式插值逼近的扩展椭球集员滤波方法
CN101871782B (zh) 基于set2fnn的gps/mems-ins组合导航系统定位误差预测方法
US8386098B2 (en) Method for a deeper search in a time-limited image satellite planning environment
CN103940433B (zh) 一种基于改进的自适应平方根ukf算法的卫星姿态确定方法
CN113091768B (zh) 一种mimu整体动态智能标定补偿方法
CN103776449B (zh) 一种提高鲁棒性的动基座初始对准方法
CN109507706B (zh) 一种gps信号丢失的预测定位方法
CN103162678B (zh) 批量mems陀螺信息融合方法
CN104048676A (zh) 基于改进粒子滤波的mems陀螺随机误差补偿方法
CN111190211A (zh) 一种gps失效位置预测定位方法
CN105675017A (zh) 一种应用于光电平台的光纤陀螺随机漂移补偿方法
CN110727968A (zh) 一种基于强化学习的导航滤波器参数优化方法
CN109000638A (zh) 一种小视场星敏感器量测延时滤波方法
CN109682383B (zh) 一种使用深空三向测量距离和数据的实时滤波定位方法
CN116734864B (zh) 一种常值观测偏差条件下航天器自主相对导航方法
CN105136150A (zh) 一种基于多次星敏感器测量信息融合的姿态确定方法
CN108332776B (zh) Mems陀螺随机误差组合预测模型的构建方法
CN110186482B (zh) 一种提高惯性制导航天器的落点精度的方法
CN109655057B (zh) 一种六推无人机加速器测量值的滤波优化方法及其系统
CN110146084B (zh) 面向卫星故障的多星编队系统分布式协同导航滤波方法
CN110502849B (zh) 一种应用于四维变分同化系统的扰动模式构建方法
CN117029881A (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