CN109839894B - 一种双边遥操作系统的控制方法 - Google Patents

一种双边遥操作系统的控制方法 Download PDF

Info

Publication number
CN109839894B
CN109839894B CN201811568035.XA CN201811568035A CN109839894B CN 109839894 B CN109839894 B CN 109839894B CN 201811568035 A CN201811568035 A CN 201811568035A CN 109839894 B CN109839894 B CN 109839894B
Authority
CN
China
Prior art keywords
matrix
representing
master
robot arm
fixed time
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
CN201811568035.XA
Other languages
English (en)
Other versions
CN109839894A (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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201811568035.XA priority Critical patent/CN109839894B/zh
Publication of CN109839894A publication Critical patent/CN109839894A/zh
Application granted granted Critical
Publication of CN109839894B publication Critical patent/CN109839894B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Feedback Control In General (AREA)

Abstract

本发明公开了一种双边遥操作系统的控制方法,包括以下步骤:步骤1:建立双边遥操作系统模型的动态方程;步骤2:设计系统固定时间控制器;步骤3:利用Lyapunov函数证明该控制器下的系统稳定性;本发明克服了双边遥操作系统稳定性和固定时间控制稳定性两者不能同时研究的难题,将固定时间控制应用到双边遥操作系统的研究中,实现了系统收敛时间的估算不依赖于初值,相对于以往的双边遥操作系统的控制方法来说,有着更好的收敛速度和性能。在一些对收敛时间有着严格要求的应用中,比如航天和海洋探索等领域,有着很大的优势,是一种可靠的、有效的、重要的控制方法。

Description

一种双边遥操作系统的控制方法
技术领域
本发明属于遥操作控制领域,具体为一种双边遥操作系统的固定时间控制方法。
背景技术
随着宇宙航空、海底探险等工作的发展,人们对危险的未知领域充满了好奇,遥操作机器人的需求也开始增大。而双边遥操作因为比传统的遥操作系统有更好的同步性能,便成为了科学家的研究热点。国内外出现了很多研究固定时间控制系统和双边遥操作系统相关建模、仿真与分析的工作案例。文献1(Jinming Huang,Zhongcai Zhang.Nonlinearfeedback design for fixed-time tracking of a class of nonlinearsystems.International Journal of Computer Mathematics.2017,94(7):1349-1362.)对一系列非线性系统的非线性反馈进行固定时间跟踪稳定性分析。文献2(Hongfei Li,Chuandong Li,Tingwen Huang,Deqiang Ouyang.Fixed-time stability andstabilization of impulsive dynamical systems.2017,354:8626-8644.)对脉冲动态系统进行固定时间稳定性的详细分析和研究。文献3(袁利,马广富,董经纬,李传江,姜博严.航天器近距离交会的固定时间终端滑模控制.2018,39(2):195-205.)则基于固定时间概念设计了一种含双幂次项的六自由度位姿终端滑模自适应控制器,不仅能全局提高姿态和位置的跟踪速度及精度,还能估计系统稳定所需的时间上界,且该上界与状态初始值无关。文献4(Yana Yang,Changchun Hua,Xinping Guan.Adaptive Fuzzy Finite-TimeCoordination Control for Networked Nonlinear Bilateral TeleoperationSystem.IEEE Transactions on fuzzy systems.2014,22(3):631-641.)研究了一种快速终端滑模下的有限时间双边遥操作系统的稳定性。相比于有限时间控制,固定时间控制估计的系统收敛时间上界不依赖于系统初值,即不管系统初始状态是在接近平衡点处还是远离平衡点处,都能保证系统在固定的时间内稳定。此外,相比于快速有限时间稳定,固定时间控制器中显式含有大于1和小于1的分数幂,能同时提高系统状态在接近平衡点阶段或远离平衡点阶段的收敛速度,即具有更快的收敛速度。
现在的研究双边遥操作系统稳定性和固定时间控制稳定性的成果越来越多,可是同时研究双边遥操作系统和固定时间控制的成果很少。固定时间控制有更高的收敛速度,并且能满足更高的性能指标要求,所以将固定时间控制应用到双边遥操作系统的研究中具有非常重要的意义。
发明内容
本发明目的在于将固定时间控制应用到双边遥操作系统的研究中,提供一种设计合理并具有良好效果的双边遥操作系统的固定时间控制方法。
为达成上述目的,本发明所采用的技术方案如下:
一种双边遥操作系统的固定时间快速终端滑模控制方法,包括以下步骤:
步骤1:建立双边遥操作系统模型的动态方程;
步骤2:设计系统固定时间控制器;
步骤3:利用Lyapunov函数证明该控制器下的系统稳定性。
更进一步的,所述步骤1建立双边遥操作系统模型的动态方程具体为:
建立系统动态方程,根据欧拉拉格朗日方程,主系统模型为:
Figure GDA0002022604150000021
式中,m表示主系统,s表示从系统,n表示系统的自由度;
qm(t)∈Rn表示主系统机械臂的位移矢量,
Figure GDA0002022604150000022
表示主机械臂的速度矢量;
Figure GDA0002022604150000023
表示主机械臂的加速度矢量,
Mm(qm)∈Rn×n表示一个对称的正定矩阵,在对称正定的条件下,存在正常数
Figure GDA0002022604150000024
Figure GDA0002022604150000025
使得
Figure GDA0002022604150000026
Figure GDA0002022604150000027
表示科里奥利力矩阵,Cm矩阵和Mm矩阵存在以下关系:矩阵
Figure GDA0002022604150000028
是反对称矩阵;对任意矩阵A∈Rn,有:
Figure GDA0002022604150000029
τdm∈Rn表示主端干扰矩阵,
gm∈Rn表示重力矩阵;
Fh∈Rn表示作用在主端的人力;
τm∈Rn表示应用矩阵;
建立系统动态方程,根据欧拉拉格朗日方程,从系统模型为:
Figure GDA0002022604150000031
式中,qs(t)∈Rn表示从系统机械臂的位移矢量;
Figure GDA0002022604150000032
表示从机械臂的速度矢量;
Figure GDA0002022604150000033
表示从机械臂的加速度矢量;
Ms(qs)∈Rn×n表示一个对称的正定矩阵;在对称正定的条件下,存在正常数
Figure GDA0002022604150000034
Figure GDA0002022604150000035
使得
Figure GDA0002022604150000036
Figure GDA0002022604150000037
表示科里奥利力矩阵;Cs矩阵和Ms矩阵存在以下关系:矩阵
Figure GDA00020226041500000311
是反对称矩阵,对任意矩阵A∈Rn,有:
Figure GDA0002022604150000038
τds∈Rn表示干扰矩阵;
gs∈Rn表示重力矩阵;
Fe∈Rn表示从端机械臂受到的环境的作用力;
τs∈Rn表示应用矩阵。
更进一步的,所述步骤2设计系统固定时间控制器具体为:
步骤2-1:建立所应用的快速终端滑模控制模型:
Figure GDA0002022604150000039
式中,αm、αs、βm和βs表示正常数,
gm<hm<2gm,gs<hs<2gs
em=qs(t-Ts)-qm表示主机械臂到从机械臂的位置同步误差;
Figure GDA00020226041500000310
表示主机械臂到从机械臂的速度跟踪误差;
es=qm(t-Ts)-qs表示从机械臂到主机械臂的位置同步误差;
Figure GDA0002022604150000041
表示从机械臂到主机械臂的速度跟踪误差;
Tm表示从主端到从端的单向传输时间延迟;
Ts表示从从端到主端的单向传输时间延迟;
Tm和Ts为已知量;
步骤2-2:根据快速终端滑模控制器,得出双边遥操作系统的动态方程:
Figure GDA0002022604150000042
其中:
Figure GDA0002022604150000043
Figure GDA0002022604150000044
式中,m表示主系统,s表示从系统,
qm(t)∈Rn表示主系统机械臂的位移矢量,
Figure GDA0002022604150000045
表示主机械臂的速度矢量;
Mm(qm)∈Rn×n表示一个对称的正定矩阵;
Figure GDA0002022604150000046
表示科里奥利力矩阵;
τdm∈Rn表示主端干扰矩阵,
τm∈Rn表示应用矩阵;
Fh∈Rn表示作用在主端的人力;
Fm(Zm)和Fs(Zs)定义如下:
Figure GDA0002022604150000047
Figure GDA0002022604150000048
步骤2-3:根据步骤2-2得出固定时间控制器:
Figure GDA0002022604150000051
式中,
Figure GDA0002022604150000052
Figure GDA0002022604150000053
表示正定对角矩阵,pm,ps<1,rm,rs>1,ξm和ξs表示正标量,且存在||τdm||≤ξm,||τds||≤ξs
更进一步的,所述步骤3利用Lyapunov函数证明该控制器下的系统稳定性具体为:
步骤3-1:设计李雅普诺夫函数为V,具体表达形式为:
Figure GDA0002022604150000054
因为矩阵Mm和Ms都为正定矩阵,所以:
Figure GDA0002022604150000055
由上式可得,李雅普诺夫函数正定;
步骤3-2:计算李雅普诺夫函数的导数:
李雅普诺夫函数求导的表达式为:
Figure GDA0002022604150000056
根据
Figure GDA0002022604150000057
为反对称矩阵,得出
Figure GDA0002022604150000058
推导出:
Figure GDA0002022604150000059
将李雅普诺夫函数的导数表达式变形为:
Figure GDA00020226041500000510
忽略重力作用,将控制器代入表达式得:
Figure GDA0002022604150000061
根据上式可得,李雅普诺夫函数的导数负定;李雅普诺夫函数正定,所以可以证明得到系统稳定;
步骤3-3:证明系统在固定时间内稳定:
Figure GDA0002022604150000062
其中,
Figure GDA0002022604150000063
Figure GDA0002022604150000064
得出:
Figure GDA0002022604150000065
所以,系统在固定时间稳定,且稳定时间为:
Figure GDA0002022604150000066
与现有技术相比,本发明的有益效果为:
本发明克服了双边遥操作系统稳定性和固定时间控制稳定性两者不能同时研究的难题,将固定时间控制应用到双边遥操作系统的研究中,实现估算不依赖于初值,相对于以往的双边遥操作系统的控制方法来说,有着更好的收敛速度,满足更高的性能指标要求。
附图说明
图1是本发明主端机械臂的人类操作力示意图。
图2是本发明在设计的固定时间控制器下系统主从端机械臂的位置跟踪仿真图。
图3是本发明在设计的固定时间控制器下系统主从端机械臂的位置跟踪误差仿真图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细的描述。
如图1、图2和图3所示,本发明是在一种固定时间控制下,双边遥操作系统在受外力作用时主从端达到跟踪稳定的研究。具体操作步骤如下:
步骤1:建立系统动力学模型:
主端:
Figure GDA0002022604150000071
从端:
Figure GDA0002022604150000072
外界环境力模型:人类操作力矩Fh=Jm T*[0 1]T*F
其中惯性矩阵Mm、Ms,科里奥利力矩阵Cm和Cs以及Jm和Js的表达式如下:
Figure GDA0002022604150000073
Figure GDA0002022604150000074
Figure GDA0002022604150000075
仿真中,各参数选择如下:
m1=0.5kg,m2=1kg,l1=1m,l2=0.8m,g=9.81m/s2
另外有
Figure GDA0002022604150000081
F的示意图如图1所示。
步骤2:建立控制器模型:
Figure GDA0002022604150000082
其中,
Figure GDA0002022604150000083
Figure GDA0002022604150000084
Figure GDA0002022604150000085
参数选择为:
αm=αs=5,βm=βs=8,gm/hm=gs/hs=7/10,pm=ps=2/3,rm=rs=4,
Figure GDA0002022604150000086
ξm=ξs=0.5。初始状态为qm(0)=qs(0)=[0 0]T
最后,对以上仿真结果作出总结,在设计的固定时间控制下,在图1所示的人类操作力的作用下,从端能在固定时间内稳定跟踪主端如图2,且从图3可以看出跟踪误差保持在0.08N.m之内。
虽然本发明已以较佳实施例揭露如上,然其并非用以限定本发明。本发明所属技术领域中具有通常知识者,在不脱离本发明的精神和范围内,当可作各种的更动与润饰。因此,本发明的保护范围当视权利要求书所界定者为准。

Claims (1)

1.一种双边遥操作系统的固定时间快速终端滑模控制方法,其特征在于,包括以下步骤:
步骤1:建立双边遥操作系统模型的动态方程;
所述步骤1建立双边遥操作系统模型的动态方程具体为:
建立系统动态方程,根据欧拉拉格朗日方程,主系统模型为:
Figure FDA0003181533320000011
m表示主系统,s表示从系统,n表示系统的自由度;
qm(t)∈Rn表示主系统机械臂的位移矢量,
Figure FDA0003181533320000012
表示主机械臂的速度矢量;
Figure FDA0003181533320000013
表示主机械臂的加速度矢量,
Mm(qm)∈Rn×n表示一个对称的正定矩阵,在对称正定的条件下,存在正常数
Figure FDA0003181533320000014
Figure FDA0003181533320000015
使得
Figure FDA0003181533320000016
Figure FDA0003181533320000017
表示科里奥利力矩阵,Cm矩阵和Mm矩阵存在以下关系:矩阵
Figure FDA0003181533320000018
是反对称矩阵;对任意矩阵A∈Rn,有:
Figure FDA0003181533320000019
τdm∈Rn表示主端干扰矩阵,
gm∈Rn表示重力矩阵;
Fh∈Rn表示作用在主端的人力;
τm∈Rn表示应用矩阵;
建立系统动态方程,根据欧拉拉格朗日方程,从系统模型为:
Figure FDA00031815333200000110
式中,qs(t)∈Rn表示从系统机械臂的位移矢量;
Figure FDA00031815333200000111
表示从机械臂的速度矢量;
Figure FDA00031815333200000112
表示从机械臂的加速度矢量;
Ms(qs)∈Rn×n表示一个对称的正定矩阵;在对称正定的条件下,存在正常数
Figure FDA00031815333200000113
Figure FDA00031815333200000114
使得
Figure FDA00031815333200000115
Figure FDA0003181533320000021
表示科里奥利力矩阵;Cs矩阵和Ms矩阵存在以下关系:矩阵
Figure FDA0003181533320000022
是反对称矩阵,对任意矩阵A∈Rn,有:
Figure FDA0003181533320000023
τds∈Rn表示干扰矩阵;
gs∈Rn表示重力矩阵;
Fe∈Rn表示从端机械臂受到的环境的作用力;
τs∈Rn表示应用矩阵;
步骤2:设计系统固定时间控制器;
所述步骤2设计系统固定时间控制器具体为:
步骤2-1:建立所应用的快速终端滑模控制模型:
Figure FDA0003181533320000024
式中,αm、αs、βm和βs表示正常数,
gm<hm<2gm,gs<hs<2gs
em=qs(t-Ts)-qm表示主机械臂到从机械臂的位置同步误差;
Figure FDA0003181533320000025
表示主机械臂到从机械臂的速度跟踪误差;
es=qm(t-Ts)-qs表示从机械臂到主机械臂的位置同步误差;
Figure FDA0003181533320000026
表示从机械臂到主机械臂的速度跟踪误差;
Tm表示从主端到从端的单向传输时间延迟;
Ts表示从从端到主端的单向传输时间延迟;
Tm和Ts为已知量;
步骤2-2:根据快速终端滑模控制器,得出双边遥操作系统的动态方程:
Figure FDA0003181533320000031
其中:
Figure FDA0003181533320000032
Figure FDA0003181533320000033
式中,m表示主系统,s表示从系统,
qm(t)∈Rn表示主系统机械臂的位移矢量,
Figure FDA0003181533320000034
表示主机械臂的速度矢量;
Mm(qm)∈Rn×n表示一个对称的正定矩阵;
Figure FDA0003181533320000035
表示科里奥利力矩阵;
τdm∈Rn表示主端干扰矩阵,
τm∈Rn表示应用矩阵;
Fh∈Rn表示作用在主端的人力;
Fm(Zm)和Fs(Zs)定义如下:
Figure FDA0003181533320000036
Figure FDA0003181533320000037
步骤2-3:根据步骤2-2得出固定时间控制器:
Figure FDA0003181533320000038
式中,
Figure FDA0003181533320000039
Figure FDA00031815333200000310
表示正定对角矩阵,pm,ps<1,rm,rs>1,ξm和ξs表示正标量,且存在||τdm||≤ξm,||τds||≤ξs
步骤3:利用Lyapunov函数证明该控制器下的系统稳定性,具体为:
步骤3-1:设计李雅普诺夫函数为V,具体表达形式为:
Figure FDA0003181533320000041
因为矩阵Mm和Ms都为正定矩阵,所以:
Figure FDA0003181533320000042
由上式可得,李雅普诺夫函数正定;
步骤3-2:计算李雅普诺夫函数的导数:
李雅普诺夫函数求导的表达式为:
Figure FDA0003181533320000043
根据
Figure FDA0003181533320000044
为反对称矩阵,得出
Figure FDA0003181533320000045
推导出:
Figure FDA0003181533320000046
将李雅普诺夫函数的导数表达式变形为:
Figure FDA0003181533320000047
忽略重力作用,将控制器代入表达式得:
Figure FDA0003181533320000048
根据上式可得,李雅普诺夫函数的导数负定;李雅普诺夫函数正定,所以可以证明得到系统稳定;步骤3-3:证明系统在固定时间内稳定:
Figure FDA0003181533320000051
其中,
Figure FDA0003181533320000052
Figure FDA0003181533320000053
得出:
Figure FDA0003181533320000054
所以,系统在固定时间稳定,且稳定时间为:
Figure FDA0003181533320000055
CN201811568035.XA 2018-12-21 2018-12-21 一种双边遥操作系统的控制方法 Active CN109839894B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811568035.XA CN109839894B (zh) 2018-12-21 2018-12-21 一种双边遥操作系统的控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811568035.XA CN109839894B (zh) 2018-12-21 2018-12-21 一种双边遥操作系统的控制方法

Publications (2)

Publication Number Publication Date
CN109839894A CN109839894A (zh) 2019-06-04
CN109839894B true CN109839894B (zh) 2021-11-23

Family

ID=66883300

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811568035.XA Active CN109839894B (zh) 2018-12-21 2018-12-21 一种双边遥操作系统的控制方法

Country Status (1)

Country Link
CN (1) CN109839894B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110262256B (zh) * 2019-07-18 2020-09-11 浙江大学 一种非线性遥操作系统的多边自适应滑模控制方法
CN110568759B (zh) * 2019-09-26 2022-06-28 南京理工大学 分数阶混沌系统的鲁棒同步控制方法
CN111427264B (zh) * 2020-03-15 2021-12-14 中国地质大学(武汉) 一种复杂遥操作技术的神经自适应固定时间控制方法
CN112783046B (zh) * 2020-12-31 2022-03-15 西北工业大学 基于模糊策略的双边遥操作末端平滑行为规划控制方法
CN114488791B (zh) * 2021-12-15 2023-07-21 西北工业大学 基于操作者意图理解的遥操作事件触发固定时间控制方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002087382A (ja) * 2000-09-14 2002-03-27 Akiyoshi Ishii 船舶の動揺低減制御の方法および装置
CN105319972A (zh) * 2015-11-27 2016-02-10 燕山大学 基于快速终端滑模的遥操作机器人固定时间控制方法
CN106647281A (zh) * 2017-01-18 2017-05-10 燕山大学 一种基于终端滑模的遥操作系统干扰有限时间补偿方法
CN107577147A (zh) * 2017-09-08 2018-01-12 北京邮电大学 一种基于自适应Smith预估器的遥操作双边PID控制方法
CN108983606A (zh) * 2018-07-09 2018-12-11 南京理工大学 一种机械臂系统的鲁棒滑模自适应控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002087382A (ja) * 2000-09-14 2002-03-27 Akiyoshi Ishii 船舶の動揺低減制御の方法および装置
CN105319972A (zh) * 2015-11-27 2016-02-10 燕山大学 基于快速终端滑模的遥操作机器人固定时间控制方法
CN106647281A (zh) * 2017-01-18 2017-05-10 燕山大学 一种基于终端滑模的遥操作系统干扰有限时间补偿方法
CN107577147A (zh) * 2017-09-08 2018-01-12 北京邮电大学 一种基于自适应Smith预估器的遥操作双边PID控制方法
CN108983606A (zh) * 2018-07-09 2018-12-11 南京理工大学 一种机械臂系统的鲁棒滑模自适应控制方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Adaptive Fuzzy Finite-Time Coordination Control for Networked Nonlinear Bilateral Teleoperation System;Yana Yang;《IEEE TRANSACTIONS ON FUZZY SYSTEMS》;20140630;第641-641页 *
Fixed-time stability and stabilization of impulsive dynamical systems;Hongfei Li;《Journal of the Franklin Institute》;20171102;第8626-8644页 *
含有未知参数的中立型神经网络自适应同步;屈艺等;《南京理工大学学报》;20121230(第06期);第968-972页 *

Also Published As

Publication number Publication date
CN109839894A (zh) 2019-06-04

Similar Documents

Publication Publication Date Title
CN109839894B (zh) 一种双边遥操作系统的控制方法
CN108015763B (zh) 一种抗噪声干扰的冗余度机械臂路径规划方法
CN107045557B (zh) 面向约束的可重构机械臂非奇异终端滑模力位置控制方法
Mohammadi et al. Nonlinear disturbance observers: Design and applications to Euler? Lagrange systems
Guo et al. Zhang neural network, Getz–Marsden dynamic system, and discrete-time algorithms for time-varying matrix inversion with application to robots' kinematic control
Delavari et al. Fractional order controller design for a flexible link manipulator robot
CN109669479B (zh) 一种基于事件触发的移动机器人轨迹跟踪控制方法
CN110262256B (zh) 一种非线性遥操作系统的多边自适应滑模控制方法
CN106681343B (zh) 一种航天器姿态跟踪低复杂度预设性能控制方法
Islam et al. Teleoperation systems with symmetric and unsymmetric time varying communication delay
Zhai et al. Adaptive finite‐time control for nonlinear teleoperation systems with asymmetric time‐varying delays
CN110543184B (zh) 一种刚性飞行器的固定时间神经网络控制方法
CN106584455B (zh) 一种遥操作机械臂系统的时延控制方法
CN108983734B (zh) 一种考虑三角形结构下遥操作系统的有限时间控制方法
CN103389650A (zh) 四通道遥操作系统的双边未知死区自适应控制系统
CN109901603A (zh) 一种输入时延下多空间飞行器姿态协同控制方法
CN108988710B (zh) 考虑长时延的网络化h∞模型参考直流电机调速方法和系统
CN108594656B (zh) 一种双边托举机器人系统高精度抗干扰连续滑模控制方法
Wei et al. Wavelet neural networks robust control of farm transmission line deicing robot manipulators
CN108406766B (zh) 一种基于复合积分滑模的多机械臂系统同步控制方法
Raouf et al. Workspace trajectory tracking control for two-flexible-link manipulator through output redefinition
Huang et al. Forwarding‐based dynamic surface control for antagonistic actuated robots
Crenganis et al. Inverse kinematics of a 7 DOF manipulator using adaptive neuro-fuzzy inference systems
Li et al. Fixed‐time fault‐tolerant control of manipulator systems based on sliding mode observer
CN109397284A (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