CN107607102B - 基于干扰观测器的mems陀螺滑模抖振抑制方法 - Google Patents

基于干扰观测器的mems陀螺滑模抖振抑制方法 Download PDF

Info

Publication number
CN107607102B
CN107607102B CN201711073626.5A CN201711073626A CN107607102B CN 107607102 B CN107607102 B CN 107607102B CN 201711073626 A CN201711073626 A CN 201711073626A CN 107607102 B CN107607102 B CN 107607102B
Authority
CN
China
Prior art keywords
neural network
formula
mems
model
gyro
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
CN201711073626.5A
Other languages
English (en)
Other versions
CN107607102A (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.)
Northwestern Polytechnical University
Shanghai Aerospace Control Technology Institute
Shenzhen Institute of Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Shanghai Aerospace Control Technology Institute
Shenzhen Institute of Northwestern Polytechnical University
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 Northwestern Polytechnical University, Shanghai Aerospace Control Technology Institute, Shenzhen Institute of Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201711073626.5A priority Critical patent/CN107607102B/zh
Publication of CN107607102A publication Critical patent/CN107607102A/zh
Application granted granted Critical
Publication of CN107607102B publication Critical patent/CN107607102B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
  • Feedback Control In General (AREA)
  • Gyroscopes (AREA)

Abstract

本发明公开了一种基于干扰观测器的MEMS陀螺滑模抖振抑制方法,用于解决现有MEMS陀螺仪模态控制方法实用性差的技术问题。技术方案是首先设计干扰观测器,在滑模控制中对干扰进行估计与补偿,从而降低抖振;同时根据神经网络预测误差和跟踪误差,设计神经网络权值的复合自适应法则律,修正神经网络的权重系数,实现未知动力学的有效动态估计。本发明通过设计神经网络权值的复合自适应律,修正神经网络的权重系数,实现未知动力学的有效动态估计。结合滑模控制理论,实现对MEMS陀螺未知动力学的前馈补偿,进一步提高MEMS陀螺仪的控制精度。设计干扰观测器,对外部干扰进行估计与补偿,有效降低了滑模抖振,实用性好。

Description

基于干扰观测器的MEMS陀螺滑模抖振抑制方法
技术领域
本发明涉及一种MEMS陀螺仪模态控制方法,特别涉及一种基于干扰观测器的MEMS陀螺滑模抖振抑制方法。
背景技术
随着非线性控制技术的发展,Park S等人将先进的智能学习和非线性控制理论引入MEMS陀螺仪模态控制过程中,对提高系统鲁棒性,改善MEMS陀螺仪性能做出了重要贡献。考虑MEMS陀螺系统中未知且动态变化的不确定及干扰,如何实现未知动力学的有效学习和滑模控制的前馈补偿,是提高陀螺性能的关键。
《Robust adaptive sliding mode control of MEMS gyroscope using T-Sfuzzy model》(Shitao Wang and Juntao Fei,《Nonlinear Dynamics》,2014年第77卷第1–2期)一文中,费俊涛等人采用T-S模糊逻辑系统学习MEMS陀螺动力学的不确定项和干扰,再利用滑模控制器对不确定及干扰进行补偿。这种方法虽然实现了不确定未知情况下的MEMS陀螺控制,但一方面由于违背了模糊逻辑逼近不确定的本意,很难实现未知动力学的有效动态估计,另一方面为消除不确定和干扰,需要很大的切换增益,带来了严重的滑模抖振。
发明内容
为了克服现有MEMS陀螺仪模态控制方法实用性差的不足,本发明提供一种基于干扰观测器的MEMS陀螺滑模抖振抑制方法。该方法首先设计干扰观测器,在滑模控制中对干扰进行估计与补偿,从而降低抖振;同时根据神经网络预测误差和跟踪误差,设计神经网络权值的复合自适应法则律,修正神经网络的权重系数,实现未知动力学的有效动态估计。本发明根据平行估计模型和动力学模型构建神经网络预测误差,结合跟踪误差,设计神经网络权值的复合自适应律,修正神经网络的权重系数,实现未知动力学的有效动态估计。结合滑模控制理论,实现对MEMS陀螺未知动力学的前馈补偿,进一步提高MEMS陀螺仪的控制精度。设计干扰观测器,对外部干扰进行估计与补偿,可有效降低滑模抖振,实用性好。
本发明解决其技术问题所采用的技术方案:一种基于干扰观测器的MEMS陀螺滑模抖振抑制方法,其特点是包括以下步骤:
(a)考虑正交误差的MEMS陀螺仪的动力学模型为:
其中,m为检测质量块的质量;Ωz为陀螺输入角速度;为静电驱动力; x*分别是MEMS陀螺仪检测质量块沿驱动轴的加速度,速度和位移;y*分别是检测质量块沿检测轴的加速度,速度和位移;dxx,dyy是阻尼系数;kxx,kyy是刚度系数;dxy是阻尼耦合系数,kxy是刚度耦合系数。
为提高机理分析准确度,对MEMS陀螺动力学模型进行无量纲化处理。取无量纲化时间t*=ωot,然后在式(1)两边同时除以参考频率的平方参考长度q0和检测质量块质量m,得到MEMS陀螺的无量纲化模型为
其中,
重新定义相关系统参数为
则MEMS陀螺的无量纲化模型化简为
令A=2S-D,B=Ω2-K,考虑环境因素和未建模因素造成的参数波动以及外部干扰的影响,则式(4)表示为
所述的无量纲化模型由状态变量q=[xy]T和控制输入u=[ux uy]T组成。其中,x,y分别为无量纲化后检测质量块沿驱动轴和检测轴的运动位移;ux uy分别表示无量纲化后施加在驱动轴和检测轴的力;A、B、C是模型的参数,且其值与陀螺仪的结构参数和动力学特性有关;P为模型参数不确定带来的未知动力学,且ΔA,ΔB为环境因素和未建模因素造成的未知的参数波动;d(t)为外部干扰。
(b)构造神经网络逼近
其中,Xin是神经网络的输入向量,且为神经网络的权值矩阵;θ(Xin)为M维基向量。基向量的第i个元素为
其中,Xmi,σi分别是该高斯函数的中心和标准差,且
定义最优估计参数w*
其中,ψ是w的集合。
因此,动力学模型的不确定项表示为
其中,ε为神经网络的逼近误差。
且不确定项的估计误差为
其中,
(c)定义神经网络预测误差为
其中,的估计值。
对式(11)求一阶导数,有
由于式(5)的平行估计模型设计为
其中,为外部干扰d(t)的估计值;Kz为正定矩阵。
定义辅助变量
z=d-Kdξnn (14)
其中,Kd为正定矩阵。
考虑式(12)和式(13),式(14)的一阶导数为
其中,
设计的估计值为
其中,Knn为正定矩阵。
则干扰观测器为
(d)建立MEMS陀螺的动力学参考模型为
其中,qd为参考振动位移信号,为qd的二阶导数;Ax,Ay分别为检测质量块沿驱动轴和检测轴振动的参考振幅;ωx,ωy分别为检测质量块沿驱动轴和检测轴振动的参考角频率。
构建跟踪误差为
e=q-qd (19)
定义滑模面
其中,β满足Hurwitz条件。则
滑模控制器设计为
其中,K0为正定矩阵。
将控制器式(22)代入式(21),有
考虑神经网络预测误差式(11)和滑模函数式(20),设计神经网络权值的复合学习律为
其中,r1,r2,r3,δ为正常数。
(e)根据得到的控制器式(22)和复合学习权重更新律式(24),返回到MEMS陀螺的动力学模型式(5),对陀螺检测质量块的振动位移和速度进行跟踪控制。
本发明的有益效果是:该方法首先设计干扰观测器,在滑模控制中对干扰进行估计与补偿,从而降低抖振;同时根据神经网络预测误差和跟踪误差,设计神经网络权值的复合自适应法则律,修正神经网络的权重系数,实现未知动力学的有效动态估计。本发明根据平行估计模型和动力学模型构建神经网络预测误差,结合跟踪误差,设计神经网络权值的复合自适应律,修正神经网络的权重系数,实现未知动力学的有效动态估计。结合滑模控制理论,实现对MEMS陀螺未知动力学的前馈补偿,进一步提高MEMS陀螺仪的控制精度。设计干扰观测器,对外部干扰进行估计与补偿,有效降低了滑模抖振,实用性好。
下面结合附图和具体实施方式对本发明作详细说明。
附图说明
图1是本发明基于干扰观测器的MEMS陀螺滑模抖振抑制方法的流程图。
具体实施方式
参照图1。本发明基于干扰观测器的MEMS陀螺滑模抖振抑制方法具体步骤如下:
(a)考虑正交误差的MEMS陀螺仪的动力学模型为:
其中,m为检测质量块的质量;Ωz为陀螺输入角速度;为静电驱动力; x*分别是MEMS陀螺仪检测质量块沿驱动轴的加速度,速度和位移;y*分别是检测质量块沿检测轴的加速度,速度和位移;dxx,dyy是阻尼系数;kxx,kyy是刚度系数;dxy是阻尼耦合系数,kxy是刚度耦合系数。
为提高机理分析准确度,对MEMS陀螺动力学模型进行无量纲化处理。取无量纲化时间然后在式(1)两边同时除以参考频率的平方参考长度q0和检测质量块质量m,可以得到MEMS陀螺的无量纲化模型为
其中,
重新定义相关系统参数为
则MEMS陀螺的无量纲化模型可化简为
令A=2S-D,B=Ω2-K,考虑环境因素和未建模因素造成的参数波动以及外部干扰的影响,则式(4)可表示为
该模型由状态变量q=[x y]T和控制输入u=[ux uy]T组成。其中,x,y分别为无量纲化后检测质量块沿驱动轴和检测轴的运动位移;ux uy分别表示无量纲化后施加在驱动轴和检测轴的力;A、B、C是模型的参数,且其值与陀螺仪的结构参数和动力学特性有关;P为模型参数不确定带来的未知动力学,且ΔA,ΔB为环境因素和未建模因素造成的未知的参数波动;d(t)为外部干扰。
根据某型号的振动式硅微机械陀螺,选取陀螺各参数为m=0.57×10-7kg,q0=[10-6 10-6]Tm,ω0=1kHz,Ωz=5.0rad/s,kxx=80.98N/m,kyy=71.62N/m,kxy=0.05N/m,dxx=0.429×10-6Ns/m,dyy=0.0429×10-6Ns/m,dxy=0.0429×10-6Ns/m,则可计算得到选取外部干扰
(b)利用神经网络动态估计模型参数不确定带来的未知动力学。
构造神经网络逼近
其中,Xin是神经网络的输入向量,且为神经网络的权值矩阵;θ(Xin)为M维基向量,M为神经网络节点数,选取M=5×5×3×3=225。基向量的第i个元素为
其中,Xmi,σi分别是该高斯函数的中心和标准差,且其值在[-2020]×[-0.240.24]×[-1010]×[-0.120.12]之间任意选取,另外σi=1。
定义最优估计参数w*
其中,ψ是w的集合。
因此,动力学模型的不确定项可表示为
其中,ε为神经网络的逼近误差。
且不确定项的估计误差为
其中,
(c)设计干扰观测器估计并补偿外部干扰。
定义神经网络预测误差为
其中,的估计值。
对式(11)求一阶导数,有
由于式(5)的平行估计模型可设计为
其中,为外部干扰d(t)的估计值;Kz为正定矩阵,取值为
定义辅助变量
z=d-Kdξnn (14)
其中,Kd为正定矩阵,取值为
考虑式(12)和式(13),式(14)的一阶导数为
其中,
设计的估计值为
其中,Knn为正定矩阵,取值为
则干扰观测器为
(d)引入滑模控制,实现未知动力学的前馈补偿,并给出神经网络权值的复合学习律。
建立MEMS陀螺的动力学参考模型为
其中,qd为参考振动位移信号,为qd的二阶导数;Ax,Ay分别为检测质量块沿驱动轴和检测轴振动的参考振幅,且Ax=10μm,Ay=0.12μm;ωx,ωy分别为检测质量块沿驱动轴和检测轴振动的参考角频率,且ωx=2000rad/s,ωy=2000rad/s。
构建跟踪误差为
e=q-qd (19)
定义滑模面
其中,β满足Hurwitz条件,取值为
滑模控制器可设计为
其中,K0为正定矩阵,取值为
将控制器式(22)代入式(21),有
考虑神经网络预测误差式(11)和滑模函数式(20),设计神经网络权值的复合学习律为
其中,r1,r2,r3,δ为正常数,且分别取值r1=0.2,r2=5,r3=2,δ=15。
(e)根据得到的控制器式(22)和复合学习权重更新律式(24),返回到MEMS陀螺的动力学模型式(5),对陀螺检测质量块的振动位移和速度进行跟踪控制。
本发明未详细说明部分属于领域技术人员公知常识。

Claims (1)

1.一种基于干扰观测器的MEMS陀螺滑模抖振抑制方法,其特征在于包括以下步骤:
(a)考虑正交误差的MEMS陀螺仪的动力学模型为:
其中,m为检测质量块的质量;Ωz为陀螺输入角速度;为静电驱动力; x*分别是MEMS陀螺仪检测质量块沿驱动轴的加速度,速度和位移;y*分别是检测质量块沿检测轴的加速度,速度和位移;dxx,dyy是阻尼系数;kxx,kyy是刚度系数;dxy是阻尼耦合系数,kxy是刚度耦合系数;
为提高机理分析准确度,对MEMS陀螺动力学模型进行无量纲化处理;取无量纲化时间t*=ωot,然后在式(1)两边同时除以参考频率的平方参考长度q0和检测质量块质量m,得到MEMS陀螺的无量纲化模型为
其中,
重新定义相关系统参数为
则MEMS陀螺的无量纲化模型化简为
令A=2S-D,B=Ω2-K,考虑环境因素和未建模因素造成的参数波动以及外部干扰的影响,则式(4)表示为
所述的无量纲化模型由状态变量q=[x y]T和控制输入u=[ux uy]T组成;其中,x,y分别为无量纲化后检测质量块沿驱动轴和检测轴的运动位移;ux uy分别表示无量纲化后施加在驱动轴和检测轴的力;A、B、C是模型的参数,且其值与陀螺仪的结构参数和动力学特性有关;P为模型参数不确定带来的未知动力学,且ΔA,ΔB为环境因素和未建模因素造成的未知的参数波动;d(t)为外部干扰;
(b)构造神经网络逼近
其中,Xin是神经网络的输入向量,且 为神经网络的权值矩阵;θ(Xin)为M维基向量;基向量的第i个元素为
其中,Xmi,σi分别是该高斯函数的中心和标准差,且
定义最优估计参数w*
其中,ψ是w的集合;
因此,动力学模型的不确定项表示为
其中,ε为神经网络的逼近误差;
且不确定项的估计误差为
其中,
(c)定义神经网络预测误差为
其中,的估计值;
对式(11)求一阶导数,有
由于式(5)的平行估计模型设计为
其中,为外部干扰d(t)的估计值;Kz为正定矩阵;
定义辅助变量
z=d-Kdξnn (14)
其中,Kd为正定矩阵;
考虑式(12)和式(13),式(14)的一阶导数为
其中,
设计的估计值为
其中,Knn为正定矩阵;
则干扰观测器为
(d)建立MEMS陀螺的动力学参考模型为
其中,qd为参考振动位移信号,为qd的二阶导数;Ax,Ay分别为检测质量块沿驱动轴和检测轴振动的参考振幅;ωx,ωy分别为检测质量块沿驱动轴和检测轴振动的参考角频率;
构建跟踪误差为
e=q-qd (19)
定义滑模面
其中,β满足Hurwitz条件;则
滑模控制器设计为
其中,K0为正定矩阵;
将控制器式(22)代入式(21),有
考虑神经网络预测误差式(11)和滑模函数式(20),设计神经网络权值的复合学习律为
其中,r1,r2,r3,δ为正常数;
(e)根据得到的控制器式(22)和复合学习权重更新律式(24),返回到MEMS陀螺的动力学模型式(5),对陀螺检测质量块的振动位移和速度进行跟踪控制。
CN201711073626.5A 2017-11-05 2017-11-05 基于干扰观测器的mems陀螺滑模抖振抑制方法 Active CN107607102B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711073626.5A CN107607102B (zh) 2017-11-05 2017-11-05 基于干扰观测器的mems陀螺滑模抖振抑制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711073626.5A CN107607102B (zh) 2017-11-05 2017-11-05 基于干扰观测器的mems陀螺滑模抖振抑制方法

Publications (2)

Publication Number Publication Date
CN107607102A CN107607102A (zh) 2018-01-19
CN107607102B true CN107607102B (zh) 2019-08-09

Family

ID=61085293

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711073626.5A Active CN107607102B (zh) 2017-11-05 2017-11-05 基于干扰观测器的mems陀螺滑模抖振抑制方法

Country Status (1)

Country Link
CN (1) CN107607102B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110456640B (zh) * 2019-07-18 2022-03-29 西北工业大学 基于非奇异终端滑模设计的mems陀螺仪参数辨识神经网络控制方法
CN110389527B (zh) * 2019-07-18 2022-04-01 西北工业大学 基于异类估计的mems陀螺仪滑模控制方法
CN110471293B (zh) * 2019-09-23 2022-02-25 南通大学 一种估计时变角速度的z轴陀螺仪滑模控制方法
CN111025915B (zh) * 2019-12-31 2022-11-25 南通大学 一种基于干扰观测器的z轴陀螺仪神经网络滑模控制方法
CN113418518A (zh) * 2020-10-21 2021-09-21 河海大学常州校区 一种微陀螺仪系统的控制方法
CN113848721B (zh) * 2021-10-09 2023-05-16 九江学院 基于高增益观测器滑模控制的冷原子重力仪主动隔振方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1821722A (zh) * 2006-03-27 2006-08-23 北京航空航天大学 一种陀螺仪标度因数和输入轴失准角的解耦测试设备
CN101158588A (zh) * 2007-11-16 2008-04-09 北京航空航天大学 一种基于集成神经网络的微小卫星用mems陀螺仪误差补偿方法
CN103616818A (zh) * 2013-11-14 2014-03-05 河海大学常州校区 微陀螺仪的自适应模糊神经全局快速终端滑模控制方法
CN103900610A (zh) * 2014-03-28 2014-07-02 哈尔滨工程大学 基于灰色小波神经网络的mems陀螺随机误差预测方法
CN104281056A (zh) * 2014-09-18 2015-01-14 河海大学常州校区 基于神经网络上界学习的微陀螺仪鲁棒自适应控制方法
CN105045097A (zh) * 2015-05-26 2015-11-11 河海大学常州校区 一种基于神经网络的微陀螺仪反演全局滑模模糊控制方法
CN107289969A (zh) * 2016-04-01 2017-10-24 南京理工大学 一种mems惯性传感器自动批量标定方法及系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8884710B2 (en) * 2011-10-25 2014-11-11 Invensense, Inc. Gyroscope with phase and duty-cycle locked loop
US10451418B2 (en) * 2015-12-09 2019-10-22 Invensense, Inc. MEMS gyroscope amplitude control via quadrature

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1821722A (zh) * 2006-03-27 2006-08-23 北京航空航天大学 一种陀螺仪标度因数和输入轴失准角的解耦测试设备
CN101158588A (zh) * 2007-11-16 2008-04-09 北京航空航天大学 一种基于集成神经网络的微小卫星用mems陀螺仪误差补偿方法
CN103616818A (zh) * 2013-11-14 2014-03-05 河海大学常州校区 微陀螺仪的自适应模糊神经全局快速终端滑模控制方法
CN103900610A (zh) * 2014-03-28 2014-07-02 哈尔滨工程大学 基于灰色小波神经网络的mems陀螺随机误差预测方法
CN104281056A (zh) * 2014-09-18 2015-01-14 河海大学常州校区 基于神经网络上界学习的微陀螺仪鲁棒自适应控制方法
CN105045097A (zh) * 2015-05-26 2015-11-11 河海大学常州校区 一种基于神经网络的微陀螺仪反演全局滑模模糊控制方法
CN107289969A (zh) * 2016-04-01 2017-10-24 南京理工大学 一种mems惯性传感器自动批量标定方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《A Nonsingular Terminal Sliding Mode Approach Using Adaptive Disturbance Observer for Finite-Time Trajectory Tracking of MEMS Triaxial Vibratory Gyroscope》;WeiWang;《Mathematical Problems in Engineering》;20151231;全文 *
《三轴微机电系统陀螺仪自适应干扰补偿方法》;王伟;《控制理论与应用》;20140430;第31卷(第4期);全文 *

Also Published As

Publication number Publication date
CN107607102A (zh) 2018-01-19

Similar Documents

Publication Publication Date Title
CN107607101B (zh) 基于干扰观测器的mems陀螺滑模控制方法
CN107607102B (zh) 基于干扰观测器的mems陀螺滑模抖振抑制方法
CN107607103B (zh) 基于干扰观测器的mems陀螺仪复合学习控制方法
CN107678282B (zh) 考虑未知动力学和外部干扰的mems陀螺智能控制方法
CN108897226B (zh) 基于干扰观测器的mems陀螺仪预设性能非奇异滑模控制方法
CN103116275B (zh) 基于滑模补偿的微陀螺仪鲁棒神经网络控制系统及方法
CN104281056B (zh) 基于神经网络上界学习的微陀螺仪鲁棒自适应控制方法
CN107608217B (zh) 基于复合学习的mems陀螺仪模糊滑模控制方法
CN102914972B (zh) 基于模型整体逼近的微陀螺仪rbf网络自适应控制方法
CN103324087B (zh) 基于神经网络的微陀螺仪的自适应反演控制系统及方法
CN110421569A (zh) 一种基于有限时间扰动观测器的反演滑模机械臂控制器设计方法
CN104503246B (zh) 微陀螺仪系统的间接自适应神经网络滑模控制方法
CN112181002B (zh) 微陀螺仪双递归扰动模糊神经网络分数阶滑模控制方法
CN107505841A (zh) 一种基于干扰估计器的机械臂姿态鲁棒控制方法
CN107870566B (zh) 基于平行估计复合学习的mems陀螺仪快速启动方法
CN110389528A (zh) 基于扰动观测的数据驱动mems陀螺仪驱动控制方法
CN110703610A (zh) 微陀螺仪的递归模糊神经网络非奇异终端滑模控制方法
CN108406779A (zh) 一种基于未知输入观测器的机械臂运动控制方法
CN105278331A (zh) 一种微陀螺的鲁棒自适应神经网络h无穷控制方法
CN113858218A (zh) 一种机械臂执行器故障诊断方法
CN111427267B (zh) 一种采用力与力矩自适应估计的高速飞行器攻角跟踪方法
Niu et al. Dynamic performance of a slender truss bridge subjected to extreme wind and traffic loads considering 18 flutter derivatives
CN107608216B (zh) 基于平行估计模型的mems陀螺仪复合学习控制方法
CN109062048B (zh) 基于复合学习的mems陀螺仪预设性能非奇异滑模控制方法
CN107861384B (zh) 基于复合学习的mems陀螺仪快速启动方法

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