CN110109361B - 一种高超声速飞行器快速平滑自适应滑模控制方法 - Google Patents
一种高超声速飞行器快速平滑自适应滑模控制方法 Download PDFInfo
- Publication number
- CN110109361B CN110109361B CN201910433547.3A CN201910433547A CN110109361B CN 110109361 B CN110109361 B CN 110109361B CN 201910433547 A CN201910433547 A CN 201910433547A CN 110109361 B CN110109361 B CN 110109361B
- Authority
- CN
- China
- Prior art keywords
- sliding mode
- aircraft
- speed
- above formula
- order
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 11
- 238000009795 derivation Methods 0.000 claims abstract description 4
- 239000013598 vector Substances 0.000 claims description 17
- 239000000446 fuel Substances 0.000 claims description 7
- 125000004432 carbon atom Chemical group C* 0.000 claims description 5
- 230000009194 climbing Effects 0.000 claims description 5
- 238000013461 design Methods 0.000 claims description 4
- 241000764238 Isis Species 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000006978 adaptation Effects 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 238000004088 simulation Methods 0.000 description 16
- 230000003044 adaptive effect Effects 0.000 description 8
- 230000004044 response Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 3
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 239000005433 ionosphere Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035484 reaction time Effects 0.000 description 1
- 239000005437 stratosphere Substances 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 239000005436 troposphere Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种高超声速飞行器快速平滑自适应滑模控制方法。建立高超声速飞行器纵向动力学模型;对飞行速度和飞行高度进行求导,得到控制输入的显式表达式,设计滑模微分器精确跟踪飞行速度和飞行高度的各阶导数,实现精确反馈线性化;选取积分滑模面,确定快速平滑自适应二阶滑模控制器的控制律和自适应律,利用快速平滑自适应二阶滑模控制器对高超声速飞行器进行控制。本发明解决了现有控制算法不平滑的问题。
Description
技术领域
本发明属于高超声速飞行器技术领域,特别涉及了一种高超声速飞行器快速平滑控制方法。
背景技术
高超声速飞行器突防、侦察能力强,反应时间短,飞行速度快,在军事和民用上有着众多需求。但是由于其具有严重非线性、强耦合性、不确定性等特性,给飞控系统的设计增加了难度,特别是在爬升段,飞行器要穿越对流层、平流层、中间层以及部分电离层,不确定性的因素更多。
超螺旋二阶滑模算法因其快速收敛性且对抖振有良好抑制效果,已广泛应用于干扰估计、极值搜索控制、容错控制等相关领域。现有文献通过在算法中加入线性项,加快了算法的收敛速度。但由于其积分项中含有非连续的符号函数,算法本身不平滑。另有文献在假设干扰上界未知的前提下,设计了一种快速自适应二阶滑模控制算法。但其没有考虑自适应参数可能过大估计的情况,并且算法本身不平滑。还有文献通过设计一种可增可减的自适应律来避免自适应参数可能过大估计的情况,但算法仍然不平滑。
发明内容
为了解决上述背景技术提到的技术问题,本发明提出了一种高超声速飞行器快速平滑自适应滑模控制方法。
为了实现上述技术目的,本发明的技术方案为:
一种高超声速飞行器快速平滑自适应滑模控制方法,包括以下步骤:
(1)建立高超声速飞行器纵向动力学模型;
(2)对飞行速度和飞行高度进行求导,得到控制输入的显式表达式,设计滑模微分器精确跟踪飞行速度和飞行高度的各阶导数,实现精确反馈线性化;
(3)选取积分滑模面,确定快速平滑自适应二阶滑模控制器的控制律和自适应律,利用快速平滑自适应二阶滑模控制器对高超声速飞行器进行控制。
进一步地,在步骤(1)中,所述高超声速飞行器纵向动力学模型包括:
(a)爬升段动力学模型:
上式中,V为飞行速度,T为发动机的推力,α为飞行迎角,D为飞行器的阻力,m为飞行器质量,μ为地球重力常数,γ为航迹倾斜角,r为飞行器地心距,r=R+h,R为地球半径,h为飞行高度,L为飞行器的升力,M为俯仰力矩,Iyy为纵向转动惯量,Isp为燃料比冲,g0为海平面重力加速度,dv1、dγ、dq为因参数不确定和外界干扰所引起的复合干扰;字母上方的一点表示该字母的一阶微分;
(b)发动机推力模型:
上式中,β为发动机节流阀调定,βc为发动机节流阀调定的指令信号,ζ为阻尼比,ωn为固有频率;字母上方的两点表示该字母的二阶微分;
推力系数公式:
上式中,CT为推力系数;
燃料比冲公式:
上式中,Ma为马赫数;
(c)空气动力与动力矩模型:
上式中,Sw为机翼浸润面积,cA为平均气动弦长,ρ为空气密度,CL、CD、CM分别为升力系数、阻力系数和俯仰力矩系数。
进一步地,在步骤(2)中,首先进行输入输出线性化:
定义向量x=[V,γ,α,β,h]T,上标T表示矩阵转置,对飞行速度V进行三次微分,对飞行高度h为四次微分,则微分表达式中出现控制输入量:发动机节流阀的调定指令信号βc和升降舵偏角δe:
其中,ce=0.0292;
然后设计滑模微分器:
进一步地,在步骤(2)中,选取积分滑模面s:
上式中,λV、λh为待设计参数;sV和sh分别为关于速度和高度的积分滑模面,其初值均为0;τ为时间积分变量,t为积分时长;eV、eh分别为V、h的跟踪误差,eV=V-Vd,eh=h-hd,Vd、hd分别为飞行器速度、高度的指令信号,d为微分符号。
进一步地,快速平滑自适应二阶滑模控制器的控制律如下:
上式中,bVi、bhi分别为关于速度、高度的控制器结构参数,bVi,bhi>0,aVi、ahi分别为关于速度、高度的自适应增益,aVi,ahi>0,令bi=[bVi bhi]T,ai=[aVi ahi]T,下标i=1,2,n为控制器结构参数,n≥2;
自适应律如下:
上式中,自适应律参数ω1,γ1>0,μ为正常数,μ≤0.1;
快速平滑自适应二阶滑模控制器的趋近算法如下:
采用上述技术方案带来的有益效果:
(1)本发明控制输入平滑,彻底消除了抖振;
(2)本发明在考虑系统存在不确定性以及外界干扰等复合干扰时,仍具有较好的动态性能;
(3)本发明由于加入线性项,具有更快的跟踪速度;
(4)本发明的自适应参数可增可减,不会造成参数的过大估计,且可以大胆设置初值,使收敛速度更快。
附图说明
图1是本发明的基本流程图;
图2是仿真1的速度和高度响应曲线图,其中(a)为速度响应曲线图,(b)为高度响应曲线图;
图3是仿真2的升降舵变化曲线图,其中(a)为FNSSTC算法,(b)为STC和FSSTC算法;
图4是仿真2的燃流量变化曲线图,其中(a)为FNSSTC算法,(b)为STC和FSSTC算法;
图5是仿真2的速度和高度响应曲线图,其中(a)为速度响应曲线图,(b)为高度响应曲线图。
具体实施方式
以下将结合附图,对本发明的技术方案进行详细说明。
本发明设计的一种高超声速飞行器快速平滑自适应滑模控制方法,如图1所示,步骤如下:
步骤1:建立高超声速飞行器纵向动力学模型;
步骤2:对飞行速度和飞行高度进行求导,得到控制输入的显式表达式,设计滑模微分器精确跟踪飞行速度和飞行高度的各阶导数,实现精确反馈线性化;
步骤3:选取积分滑模面,确定快速平滑自适应二阶滑模控制器的控制律和自适应律,利用快速平滑自适应二阶滑模控制器对高超声速飞行器进行控制。
在本实施例中,上述步骤1采用如下优选方案实现:
所述高超声速飞行器纵向动力学模型包括:
(a)爬升段动力学模型:
上式中,V为飞行速度,T为发动机的推力,α为飞行迎角,D为飞行器的阻力,m为飞行器质量,μ为地球重力常数,γ为航迹倾斜角,r为飞行器地心距,r=R+h,R为地球半径,h为飞行高度,L为飞行器的升力,M为俯仰力矩,Iyy为纵向转动惯量,Isp为燃料比冲,g0为海平面重力加速度,dv1、dγ、dq为因参数不确定和外界干扰所引起的复合干扰;字母上方的一点表示该字母的一阶微分;
(b)发动机推力模型:
上式中,β为发动机节流阀调定,βc为发动机节流阀调定的指令信号,ζ为阻尼比,ωn为固有频率;字母上方的两点表示该字母的二阶微分;
推力系数公式:
上式中,CT为推力系数;
燃料比冲公式:
上式中,Ma为马赫数;
(c)空气动力与动力矩模型:
上式中,Sw为机翼浸润面积,cA为平均气动弦长,ρ为空气密度,CL、CD、CM分别为升力系数、阻力系数和俯仰力矩系数。
在本实施例中,上述步骤2采用如下优选方案实现:
(1)输入/输出线性化
定义向量x=[V,γ,α,β,h,m]T,将V、h分别对时间求三次和四次微分,微分表达式中将出现控制输入量:发动机节流阀的调定值βc和升降舵偏角δe,将其作为整体提取出来,可得:
其中,f为Rn上n维的光滑向量场,表示V对向量场f微分3次,表示h对向量场f微分4次,同理,g为Rn上的光滑向量场,表示对向量场g微分一次,表示对向量场g微分一次。dv1、dv2、dh1、dh2表示非匹配不确定性,dv3、dh3表示匹配不确定性,u=[βc δe]T为控制输入。
其中,ce=0.0292为常量。
(2)设计滑模微分器:
输入输出线性化的过程是对V和h求微分,线性化后的模型(8)理应与原模型(1)等价,但由于传统的输入输出线性化没有考虑(6)-(7)中的非匹配不确定性,造成线性化后的模型与愿模型存在偏差,本专利用滑模微分器来逼近 以得到更加精确的值,消除或者减小线性化后的模型(8)与原模型(1)的偏差。
根据滑模微分器的性质,只要选择合适的参数,就可以保证滑模微分器的状态量在有限时间内精确跟踪速度和高度的各阶导数。
在本实施例中,上述步骤3采用如下优选方案实现:
(1)选取积分滑模面:
上式中,λV、λh为待设计参数;sV和sh分别为关于速度和高度的积分滑模面,其初值均为0;τ为时间积分变量,t为积分时长;eV、eh分别为V、h的跟踪误差,eV=V-Vd,eh=h-hd,Vd、hd分别为飞行器速度、高度的指令信号,d为微分符号。
(2)快速平滑自适应二阶滑模控制器的控制律如下:
上式中,bVi、bhi分别为关于速度、高度的控制器结构参数,bVi,bhi>0,aVi、ahi分别为关于速度、高度的自适应增益,aVi,ahi>0,令bi=[bVi bhi]T,ai=[aVi ahi]T,下标i=1,2,上标T表示矩阵转置,n为控制器结构参数,n≥2;
自适应律如下:
上式中,自适应律参数ω1,γ1>0,μ为正常数,μ≤0.1。
根据式(8)、(10)、(11)可计算得出快速平滑趋近算法:
情况1:当n=2时,该算法由于相比传统二阶滑模加入了线性项,可保证|s|快速收敛到μ邻域内,但是该算法并不平滑。
情况2:当n>2时,不仅保证快速收敛,且控制输入式(8)的积分项中不含非连续的符号函数项,故此时的控制输入是平滑的。
仿真实验:
(a)仿真1
爬升过程的初始条件以及指令信号分别设为:V0=590m/s,h0=3528m,Vd=590+10t,hd=3528+75t。仿真时长为300s。
为验证滑模微分器对反馈线性化过程中非匹配不确定性的抑制作用,进行仿真1。滑模微分器的参数设置为:控制器参数设置为:n=4,bV1=bh1=9.51,bV2=bh2=0.01,aV1(0)=ah1(0)=1。仿真条件:对参数Sw,cA,Iy,ρ,CD,CM,CL加入40%的参数摄动,此外加入除参数不确定性的复合干扰:dv1=10m/s2,dγ=0.03rad/s,dq=0.1rad/s2。仿真结果如图2所示。
(b)仿真2
为验证本专利控制算法既快速又平滑的特性,进行仿真2,FSSTC表示本发明中n=4的快速且平滑的算法,FNSSTC表示本发明中n=2的快速但不平滑算法,STC表示一般超螺旋二阶滑模算法。其余参数和仿真条件如仿真1,仿真结果如图2。
从图3、4可以看,传统的二阶滑模的抖振依旧没有消除,控制输入不平滑,呈现锯齿状,而本专利提出的算法,控制输入平滑,基本消除了抖振。从图5可以看出,三种算法中,FSSTC算法的收敛速度最快,FNSSTC算法比传统二阶滑模控制(STC)算法的收敛速度略快一些,虽然其在收敛速度上有些许提升,但由于本专利加入的不确定性及干扰较大,该算法的控制效果差,逐渐振荡至发散。
从仿真结果可以看出,使用滑模微分器估计速度和高度各阶导数的控制算法,输出信号能准确地跟踪上指令信号,而未使用滑模微分器的控制器在接近100s的时候发散了。
综上所述,本发明控制算法中引入了线性项,原理上收敛速度较传统二阶滑模控制算法要快,但是设计自适应律可能减小其收敛速度,又由于本发明设计的自适应律,其参数既可增大又可减小,不会造成参数过大估计,即可选取较大的初始值,保持更快收敛速度。
实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (1)
1.一种高超声速飞行器快速平滑自适应滑模控制方法,其特征在于,包括以下步骤:
(1)建立高超声速飞行器纵向动力学模型;
所述高超声速飞行器纵向动力学模型包括:
(a)爬升段动力学模型:
上式中,V为飞行速度,T为发动机的推力,α为飞行迎角,D为飞行器的阻力,m为飞行器质量,μ为地球重力常数,γ为航迹倾斜角,r为飞行器地心距,r=R+h,R为地球半径,h为飞行高度,L为飞行器的升力,M为俯仰力矩,Iyy为纵向转动惯量,Isp为燃料比冲,g0为海平面重力加速度,dv1、dγ、dq为因参数不确定和外界干扰所引起的复合干扰;字母上方的一点表示该字母的一阶微分;
(b)发动机推力模型:
上式中,β为发动机节流阀调定,βc为发动机节流阀调定的指令信号,ζ为阻尼比,ωn为固有频率;字母上方的两点表示该字母的二阶微分;
推力系数公式:
上式中,CT为推力系数;
燃料比冲公式:
上式中,Ma为马赫数;
(c)空气动力与动力矩模型:
上式中,Sw为机翼浸润面积,cA为平均气动弦长,ρ为空气密度,CL、CD、CM分别为升力系数、阻力系数和俯仰力矩系数;
(2)对飞行速度和飞行高度进行求导,得到控制输入的显式表达式,设计滑模微分器精确跟踪飞行速度和飞行高度的各阶导数,实现精确反馈线性化;
首先进行输入输出线性化:
定义向量x=[V,γ,α,β,h]T,上标T表示矩阵转置,对飞行速度V进行三次微分,对飞行高度h为四次微分,则微分表达式中出现控制输入量:发动机节流阀的调定指令信号βc和升降舵偏角δe:
其中,ce=0.0292;
然后设计滑模微分器:
(3)选取积分滑模面,确定快速平滑自适应二阶滑模控制器的控制律和自适应律,利用快速平滑自适应二阶滑模控制器对高超声速飞行器进行控制;
选取积分滑模面s:
上式中,λV、λh为待设计参数;sV和sh分别为关于速度和高度的积分滑模面,其初值均为0;τ为时间积分变量,t为积分时长;eV、eh分别为V、h的跟踪误差,eV=V-Vd,eh=h-hd,Vd、hd分别为飞行器速度、高度的指令信号,d为微分符号;
快速平滑自适应二阶滑模控制器的控制律如下:
上式中,bVi、bhi分别为关于速度、高度的控制器结构参数,bVi,bhi>0,aVi、ahi分别为关于速度、高度的自适应增益,aVi,ahi>0,令bi=[bVi bhi]T,ai=[aVi ahi]T,下标i=1,2,n为控制器结构参数,n≥2;
自适应律如下:
上式中,自适应律参数ω1,γ1>0,μ为正常数,μ≤0.1;
快速平滑自适应二阶滑模控制器的趋近算法如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910433547.3A CN110109361B (zh) | 2019-05-23 | 2019-05-23 | 一种高超声速飞行器快速平滑自适应滑模控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910433547.3A CN110109361B (zh) | 2019-05-23 | 2019-05-23 | 一种高超声速飞行器快速平滑自适应滑模控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110109361A CN110109361A (zh) | 2019-08-09 |
CN110109361B true CN110109361B (zh) | 2020-09-15 |
Family
ID=67491821
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910433547.3A Active CN110109361B (zh) | 2019-05-23 | 2019-05-23 | 一种高超声速飞行器快速平滑自适应滑模控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110109361B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110794863B (zh) * | 2019-11-20 | 2021-05-28 | 中山大学 | 一种控制性能指标可定制的重型运载火箭姿态控制方法 |
CN112034867A (zh) * | 2019-12-25 | 2020-12-04 | 西北工业大学 | 一种无人机控制算法 |
CN111125971B (zh) * | 2019-12-26 | 2022-03-25 | 北京航空航天大学 | 一种吸气式高超声速飞行器推力不确定性确定方法 |
CN111596686B (zh) * | 2020-05-13 | 2022-06-07 | 中国人民解放军海军工程大学 | 高超声速飞行器纵向系统的预设性能控制方法 |
CN112882388B (zh) * | 2021-01-12 | 2022-02-11 | 北京航空航天大学 | 一种执行器受损下的四旋翼无人机故障检测与诊断方法 |
CN113110028B (zh) * | 2021-04-14 | 2022-06-21 | 福州大学 | 一种飞行环境模拟控制系统的环境信号处理方法 |
CN116736740A (zh) * | 2023-08-09 | 2023-09-12 | 东方空间技术(山东)有限公司 | 一种用于多飞行器的仿真系统及方法 |
-
2019
- 2019-05-23 CN CN201910433547.3A patent/CN110109361B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110109361A (zh) | 2019-08-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110109361B (zh) | 一种高超声速飞行器快速平滑自适应滑模控制方法 | |
CN109358634B (zh) | 一种高超声速飞行器鲁棒自适应控制方法 | |
CN107479384B (zh) | 高超声速飞行器神经网络复合学习非反步控制方法 | |
CN109426146B (zh) | 高超声速飞行器的高阶非奇异Terminal滑模控制方法 | |
CN109144084B (zh) | 一种基于固定时间收敛观测器的垂直起降重复使用运载器姿态跟踪控制方法 | |
CN108828957B (zh) | 基于切换机制的飞行器全局有限时间神经网络控制方法 | |
CN111665857B (zh) | 基于复合智能学习的变体飞行器控制方法 | |
CN108663940B (zh) | 基于集总复合估计的飞行器神经网络学习控制方法 | |
CN107479383B (zh) | 基于鲁棒设计的高超声速飞行器神经网络复合学习控制方法 | |
CN107390531B (zh) | 参数学习有限时间收敛的高超声速飞行器控制方法 | |
CN111367182A (zh) | 考虑输入受限的高超声速飞行器抗干扰反步控制方法 | |
CN106896722B (zh) | 采用状态反馈与神经网络的高超飞行器复合控制方法 | |
CN109062055A (zh) | 一种基于Back-stepping鲁棒自适应动态面的近空间飞行器控制系统 | |
CN110568765A (zh) | 面向攻角跟踪的高超声速飞行器非对称输出受限控制方法 | |
CN108170886B (zh) | 基于预设性能的高超声速飞行器纵向减损控制方法 | |
CN107632518B (zh) | 基于高增益观测器的高超声速飞行器神经网络复合学习控制方法 | |
CN109164708B (zh) | 一种高超声速飞行器神经网络自适应容错控制方法 | |
CN108459505A (zh) | 一种适合控制迭代设计的非常规布局飞行器快速建模方法 | |
CN109703769B (zh) | 一种基于预瞄策略的空中加油对接控制方法 | |
CN110320927A (zh) | 智能变形飞行器的飞行控制方法及系统 | |
CN112631316B (zh) | 变负载四旋翼无人机的有限时间控制方法 | |
CN112327926B (zh) | 一种无人机编队的自适应滑模控制方法 | |
CN108873923B (zh) | 应急指挥控制固定翼无人机舰面紧急起飞控制方法 | |
CN109062234A (zh) | 一种非最小相位飞行器复合学习滑模控制方法 | |
CN107831653B (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 |