CN108958281B - 基于微分求积法的无人机悬吊运输稳定性分析和控制方法 - Google Patents
基于微分求积法的无人机悬吊运输稳定性分析和控制方法 Download PDFInfo
- Publication number
- CN108958281B CN108958281B CN201710359396.2A CN201710359396A CN108958281B CN 108958281 B CN108958281 B CN 108958281B CN 201710359396 A CN201710359396 A CN 201710359396A CN 108958281 B CN108958281 B CN 108958281B
- Authority
- CN
- China
- Prior art keywords
- representing
- unmanned aerial
- aerial vehicle
- time
- matrix
- 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 37
- 239000000725 suspension Substances 0.000 title claims abstract description 36
- 239000011159 matrix material Substances 0.000 claims abstract description 43
- 238000001228 spectrum Methods 0.000 claims abstract description 17
- 230000004044 response Effects 0.000 claims description 13
- 230000007704 transition Effects 0.000 claims description 9
- 230000010354 integration Effects 0.000 claims description 7
- 230000001133 acceleration Effects 0.000 claims description 6
- 239000000470 constituent Substances 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 230000000704 physical effect Effects 0.000 claims description 3
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 abstract description 4
- 230000008901 benefit Effects 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000011217 control strategy Methods 0.000 abstract description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 241000282414 Homo sapiens Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提供了一种基于微分求积法的无人机悬吊运输稳定性分析和控制方法,包括:基于无人机***的位置环控制***以及悬吊***的特征参数,首先建立***的动力学方程,然后设计时滞反馈***,并确定该***的状态空间方程,在此基础上,引入微分求积方法,求取相邻两时滞段之间的转移矩阵,该转移矩阵谱半径小于1时,***在时域内渐进收敛。根据这一基本原理,可求取***的稳定边界;然后在此稳定边界内,以极小化谱半径为目标,即可求得最速收敛控制参数。采用上述稳定域判定以及最优控制参数求取方法,可快速确定无人机悬吊运输***的最佳飞行控制策略,从而大幅提升无人机的飞行性能,产生良好的经济效益。
Description
技术领域
本发明涉及无人机技术领域,具体地,涉及基于微分求积法的无人机悬吊运输稳定性分析和控制方法。
背景技术
在无人机的诸多应用当中,载物飞行是其有效辅助人类快速完成远距离搬运作业的一类重要功用。当前无人机单体飞行控制已经取得比较成熟的研究成果,但对于有外接机构的情形,***表现出时滞特征。对于此类***,其稳定性边界更难确定,因此也更难达到最优控制的目标。欲解决上述应用问题,目前通常有两种解决思路。一是引入输入整形,以消除悬吊飞行***固有的输入延迟。二是恰当的配置时滞反馈,以获取最佳的收敛速度。第一种方法虽已被充分讨论,但其配置方法与实施均有别于通常意义上的反馈控制,而第二种方法则与目前普遍采用的反馈控制方法具有相同的结构。配置时滞反馈***的关键在于快速准确地确定***的稳定区域,并搜索该区域获取最速收敛控制参数。鉴于上述考虑,本研究提出基于微分求积方法的无人机悬吊运输控制***稳定域判定方法与最优控制参数搜索策略,以期使无人机***稳定高效的完成相应作业任务,具有重要的理论与现实意义。
本发明依据多旋翼无人机的动力学特性,提供一种针对无人机悬吊运输飞行的稳定域判定以及最优控制参数搜索方法。具体在建立悬吊运输飞行控制***动力学方程的基础上,采用微分求积方法确定***的稳定区域,最后在此区域范围内,求取使***响应最速收敛的控制参数。在提升计算效率的前提下,优化无人机悬吊运输的飞行性能,产生良好的经济效益。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种基于微分求积法的无人机悬吊运输稳定性分析和控制方法。
根据本发明提供的基于微分求积法的无人机悬吊运输稳定性分析和控制方法,包括如下步骤:
步骤1:根据悬吊***的物性特征,确定悬吊***的动力学方程;
步骤2:以一阶时滞***近似描述无人机本体动力学响应特征;
步骤3:设计时滞反馈***,建立闭环***状态空间方程;
步骤4:基于微分求积方法,确定相邻两时滞段之间的转移矩阵;
步骤5:计算转移矩阵谱半径小于1的区域,所述区域即为***稳定域;
步骤6:在稳定域内搜索谱半径极小点,所述谱半径极小点对应控制参数即为所需的最优控制点。
优选地,所述步骤1包括:假设一个质量为ml,悬绳长度为l的悬吊***,其动力学方程可以表述如下:
式中:ml表示重物的质量,xl表示重物在机体坐标系中相对于X轴的坐标,yl表示重物在机体坐标系中相对于Y轴的坐标,g表示重力加速度,b表示无人机机体厚度的一半,θv及φv分别表示无人机横滚角、俯仰角,σ表示空气阻尼系数,表示xl的一阶导数,表示xl的二阶导数,表示yl的一阶导数,表示yl的二阶导数,xv表示无人机X方向位置坐标,yv表示无人机Y方向位置坐标。
优选地,所述步骤2包括:基于悬吊***动力学响应方程以及无人机本体响应特性,设计时滞反馈,可得悬吊运输***闭环响应特性如下:
式中:表示负载在X方向上的状态变量微分,kv表示无人机X方向传递函数增益,x*表示参考轨迹,其中,Kp=K11,K11表示反馈增益矩阵K第一行第一列,由于考察***稳定域时,忽略输入参考项x*,得到如下形式:
其中h=τ+τv,且
其中
式中:Ψ表示区间转移左半矩阵,x(t0)表示t0时刻负载X方向状态,x(tn)表示tn时刻,x(t0-h)表示t0-h时刻负载X方向状态;h表示时滞,Γ表示区间转移右半矩阵,x(tn-h)表示tn-h时刻负载X方向状态,In×n表示n维单位矩阵,D表示矩阵Ψ的构成子矩阵,Hij表示微分求积插值系数,表示插值基函数,τi表示第i个离散时间点,τk表示第k个时间点,τj表示第j个时间点,表示归一化状态,表示时刻负载归一化状态,Aij表示矩阵第i行j列取值。
优选地,所述步骤5包括:时滞区间转移矩阵关系为Φ=Ψ-1Γ,欲使***稳定,则需保证转移矩阵Φ的谱半径小于1,在合理的参数空间范围内均匀选取参数值Kp与反馈时滞τ,求取不同取值条件下的转移矩阵谱半径r(Φ),当r(Φ)<1时,即为***稳定区域。
与现有技术相比,本发明具有如下的有益效果:
本发明提供的基于微分求积法的无人机悬吊运输稳定性分析和控制方法,能够求取无人机***的稳定边界;然后在此稳定边界内,以极小化谱半径为目标,即可求得最速收敛控制参数,采用上述稳定域判定以及最优控制参数求取方法,可快速确定无人机悬吊运输***的最佳飞行控制策略,从而大幅提升无人机的飞行性能。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明实施例中的无人机悬吊运输控制示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
根据本发明提供的基于微分求积法的无人机悬吊运输稳定性分析和控制方法,包括以下步骤:
步骤1:根据悬吊***的物性特征,确定悬吊***的动力学方程;
步骤2:以一阶时滞***近似描述无人机本体动力学响应特征;
步骤3:设计时滞反馈***,建立闭环***状态空间方程;
步骤4:基于微分求积方法,确定相邻两时滞段之间的转移矩阵;
步骤5:计算转移矩阵谱半径小于1的区域,该区域即为***稳定域;
步骤6:在稳定域内搜索谱半径极小点,该点对应控制参数即为所需的最优控制点。
具体地,如图1所示,其中XI-YI-ZI为惯性坐标系,x-y-z为机体坐标系,对于一个质量为ml,悬绳长度为l的悬吊***,其动力学方程可以表述如下:
式中:ml表示重物的质量,xl表示重物在机体坐标系中相对于X轴的坐标,yl表示重物在机体坐标系中相对于Y轴的坐标,g表示重力加速度,b表示无人机机体厚度的一半,θv及φv分别表示无人机横滚角、俯仰角,σ表示空气阻尼系数,表示xl的一阶导数,表示xl的二阶导数,表示yl的一阶导数,表示yl的二阶导数,xv表示无人机X方向坐标,yv表示无人机Y方向坐标。
上述动力学方程可以简写为以下矩阵表达式:
对于无人机本体***,可以采用一阶时滞***描述其响应,其中:Gv表示传函数,Kv表示传递函数增益,τv表示无人机传递函数时滞,s表示拉氏变量。基于上述悬吊***动力学响应方程以及无人机本体响应特性,设计时滞反馈,可得悬吊运输***闭环响应特性如下:
式中:表示负载在X方向上的状态变量微分,kv表示无人机X方向传递函数增益,x*表示参考轨迹,其中,Kp=K11,K11表示反馈增益矩阵K第一行第一列于考察***稳定域时,可忽略输入参考项x*。因此,此飞行悬吊运输***可以写成如下形式:
其中h=τ+τv,且
其中
式中:Ψ表示区间转移左半矩阵,x(t0)表示t0时刻负载X方向状态,x(tn)表示tn时刻,x(t0-h)表示t0-h时刻负载X方向状态;h表示时滞,Γ表示区间转移右半矩阵,x(tn-h)表示tn-h时刻负载X方向状态,In×n表示n维单位矩阵,D表示矩阵Ψ的构成子矩阵,Hij表示微分求积插值系数,表示插值基函数,τi表示第i个离散时间点,τk表示第k个时间点,τj表示第j个时间点,表示归一化状态,表示时刻负载归一化状态,Aij表示矩阵第i行j列取值。
时滞区间转移矩阵关系为Φ=Ψ-1Γ,欲使***稳定,则需保证转移矩阵Φ的谱半径小于1。在合理的参数空间范围内均匀选取参数值Kp与反馈时滞τ,求取不同取值条件下的转移矩阵谱半径r(Φ),当r(Φ)<1时,即为***稳定区域。在该区域内,求取最小谱半径所对应参数取值,即为最优控制参数。
需要说明的是,本发明提供的所述基于微分求积法的无人机悬吊运输稳定性分析和控制方法中的步骤,可以利用基于微分求积法的无人机悬吊运输稳定性分析和控制***中对应的模块、装置、单元等予以实现,本领域技术人员可以参照所述***的技术方案实现所述方法的步骤流程,即,所述***中的实施例可理解为实现所述方法的优选例,在此不予赘述。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (1)
1.一种基于微分求积法的无人机悬吊运输稳定性分析和控制方法,其特征在于,包括如下步骤:
步骤1:根据悬吊***的物性特征,确定悬吊***的动力学方程;
步骤2:以一阶时滞***近似描述无人机本体动力学响应特征;
步骤3:设计时滞反馈***,建立闭环***状态空间方程;
步骤4:基于微分求积方法,确定相邻两时滞段之间的转移矩阵;
步骤5:计算转移矩阵谱半径小于1的区域,所述区域即为***稳定域;
步骤6:在稳定域内搜索谱半径极小点,所述谱半径极小点对应控制参数即为所需的最优控制点;
所述步骤1包括:假设一个质量为ml,悬绳长度为l的悬吊***,其动力学方程表述如下:
式中:ml表示重物的质量,xl表示重物在机体坐标系中相对于X轴的坐标,yl表示重物在机体坐标系中相对于Y轴的坐标,g表示重力加速度,b表示无人机机体厚度的一半,θv及φv分别表示无人机横滚角、俯仰角,σ表示空气阻尼系数,表示xl的一阶导数,表示xl的二阶导数,表示yl的一阶导数,表示yl的二阶导数,xv表示无人机X方向位置坐标,yv表示无人机Y方向位置坐标;
所述步骤2包括:基于悬吊***动力学响应方程以及无人机本体响应特性,设计时滞反馈,可得悬吊运输***闭环响应特性如下:
式中:表示负载在X方向上的状态变量微分,kv表示无人机X方向传递函数增益,x*表示参考轨迹,其中,Kp=K11,K11表示反馈增益矩阵K第一行第一列,由于考察***稳定域时,忽略输入参考项x*,得到如下形式:
其中h=τ+τv,且
其中
式中:Ψ表示区间转移左半矩阵,x(t0)表示t0时刻负载X方向状态,x(tn)表示tn时刻,x(t0-h)表示t0-h时刻负载X方向状态;h表示时滞,Γ表示区间转移右半矩阵,x(tn-h)表示tn-h时刻负载X方向状态,In×n表示n维单位矩阵,D表示矩阵Ψ的构成子矩阵,Hij表示微分求积插值系数,表示插值基函数,τi表示第i个离散时间点,τk表示第k个时间点,τj表示第j个时间点,表示归一化状态,表示时刻负载归一化状态,Aij表示矩阵第i行j列取值;
所述步骤5包括:时滞区间转移矩阵关系为Φ=Ψ-1Γ,欲使***稳定,则需保证转移矩阵Φ的谱半径小于1,在合理的参数空间范围内均匀选取参数值Kp与反馈时滞τ,求取不同取值条件下的转移矩阵谱半径r(Φ),当r(Φ)<1时,即为***稳定区域。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710359396.2A CN108958281B (zh) | 2017-05-19 | 2017-05-19 | 基于微分求积法的无人机悬吊运输稳定性分析和控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710359396.2A CN108958281B (zh) | 2017-05-19 | 2017-05-19 | 基于微分求积法的无人机悬吊运输稳定性分析和控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108958281A CN108958281A (zh) | 2018-12-07 |
CN108958281B true CN108958281B (zh) | 2021-04-02 |
Family
ID=64462104
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710359396.2A Active CN108958281B (zh) | 2017-05-19 | 2017-05-19 | 基于微分求积法的无人机悬吊运输稳定性分析和控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108958281B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109799698B (zh) * | 2019-01-30 | 2020-07-14 | 上海交通大学 | 时滞视觉伺服***最优pi参数优化方法及*** |
CN110032067B (zh) * | 2019-03-07 | 2021-06-11 | 上海交通大学 | 基于时滞反馈的无人机圆周路径悬吊运输控制方法及*** |
CN110045633A (zh) * | 2019-04-28 | 2019-07-23 | 上海交通大学 | 无人机双索运输减振控制方法、***及存储介质 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102749851A (zh) * | 2012-07-24 | 2012-10-24 | 北京航空航天大学 | 一种挠性高超声速飞行器的精细抗干扰跟踪控制器 |
CN102874413A (zh) * | 2012-10-12 | 2013-01-16 | 北京航空航天大学 | 提升直升机吊挂体的飞行稳定性的方法和尾鳍 |
CN102991663A (zh) * | 2012-11-05 | 2013-03-27 | 上海交通大学无锡研究院 | 一种系留旋翼升空平台*** |
US9146557B1 (en) * | 2014-04-23 | 2015-09-29 | King Fahd University Of Petroleum And Minerals | Adaptive control method for unmanned vehicle with slung load |
CN105138006A (zh) * | 2015-07-09 | 2015-12-09 | 哈尔滨工程大学 | 一种时滞非线性多智能体***的协同追踪控制方法 |
CN105460210A (zh) * | 2015-12-04 | 2016-04-06 | 上海浦江桥隧大桥管理有限公司 | 用于桥梁检测的便携式六旋翼飞行器 |
CN105930129A (zh) * | 2016-04-12 | 2016-09-07 | 广东顺德中山大学卡内基梅隆大学国际联合研究院 | 一种应用于直升机吊挂负载摆动角镇定的时变除法方法 |
CN106647781A (zh) * | 2016-10-26 | 2017-05-10 | 广西师范大学 | 基于重复控制补偿神经模糊pid四旋翼飞行器的控制方法 |
CN106647792A (zh) * | 2017-01-25 | 2017-05-10 | 天津大学 | 用于无人机悬挂负载***的抗扰控制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160282861A1 (en) * | 2015-03-24 | 2016-09-29 | Wolfgis, Llc | Method for keeping drones within a designated boundary |
-
2017
- 2017-05-19 CN CN201710359396.2A patent/CN108958281B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102749851A (zh) * | 2012-07-24 | 2012-10-24 | 北京航空航天大学 | 一种挠性高超声速飞行器的精细抗干扰跟踪控制器 |
CN102874413A (zh) * | 2012-10-12 | 2013-01-16 | 北京航空航天大学 | 提升直升机吊挂体的飞行稳定性的方法和尾鳍 |
CN102991663A (zh) * | 2012-11-05 | 2013-03-27 | 上海交通大学无锡研究院 | 一种系留旋翼升空平台*** |
US9146557B1 (en) * | 2014-04-23 | 2015-09-29 | King Fahd University Of Petroleum And Minerals | Adaptive control method for unmanned vehicle with slung load |
CN105138006A (zh) * | 2015-07-09 | 2015-12-09 | 哈尔滨工程大学 | 一种时滞非线性多智能体***的协同追踪控制方法 |
CN105460210A (zh) * | 2015-12-04 | 2016-04-06 | 上海浦江桥隧大桥管理有限公司 | 用于桥梁检测的便携式六旋翼飞行器 |
CN105930129A (zh) * | 2016-04-12 | 2016-09-07 | 广东顺德中山大学卡内基梅隆大学国际联合研究院 | 一种应用于直升机吊挂负载摆动角镇定的时变除法方法 |
CN106647781A (zh) * | 2016-10-26 | 2017-05-10 | 广西师范大学 | 基于重复控制补偿神经模糊pid四旋翼飞行器的控制方法 |
CN106647792A (zh) * | 2017-01-25 | 2017-05-10 | 天津大学 | 用于无人机悬挂负载***的抗扰控制方法 |
Non-Patent Citations (4)
Title |
---|
2015 Asian Conference on Defence Technology (ACDT);Thanan Yomchinda;《2015 Asian Conference on Defence Technology (ACDT)》;20150525;全文 * |
Spot hover control of helicopter and swing control of helicopter sling load by using Zhang-gradient method;Liao Bolin;《2015 34th Chinese Control Conference (CCC)》;20150914;第506-511页 * |
无人机吊挂飞行的非线性控制方法设计;鲜斌;《控制理论与应用》;20160331;第273-279页 * |
直升机吊挂飞行耦合分析;陈元;《南京航空航天大学学报》;20170430;第165-172页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108958281A (zh) | 2018-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | A novel control scheme for quadrotor UAV based upon active disturbance rejection control | |
CN107357166B (zh) | 小型无人直升机的无模型自适应鲁棒控制方法 | |
CN106444799B (zh) | 基于模糊扩张状态观测器和自适应滑模的四旋翼无人机控制方法 | |
He et al. | A Simple Attitude Control of Quadrotor Helicopter Based on Ziegler‐Nichols Rules for Tuning PD Parameters | |
CN107491081B (zh) | 一种抗干扰四旋翼无人机姿态控制方法 | |
Wang et al. | Adaptive dynamic surface control for cooperative path following of marine surface vehicles with input saturation | |
Wang et al. | Adaptive formation control of quadrotor unmanned aerial vehicles with bounded control thrust | |
CN108958281B (zh) | 基于微分求积法的无人机悬吊运输稳定性分析和控制方法 | |
Ding et al. | Chaotic Artificial Bee Colony Algorithm for System Identification of a Small‐Scale Unmanned Helicopter | |
CN110442020B (zh) | 一种新型的基于鲸鱼优化算法的容错控制方法 | |
CN103853157A (zh) | 一种基于自适应滑模的飞行器姿态控制方法 | |
Yu et al. | Attitude tracking control of a quadrotor UAV in the exponential coordinates | |
CN103558857A (zh) | 一种btt飞行器的分布式复合抗干扰姿态控制方法 | |
Roberts et al. | A new position regulation strategy for VTOL UAVs using IMU and GPS measurements | |
CN110794857B (zh) | 考虑外部风干扰的固定翼无人机鲁棒离散分数阶控制方法 | |
CN110865539A (zh) | 一种随机干扰下的无人艇跟踪误差约束控制方法 | |
Liu et al. | Robust command filtered adaptive backstepping control for a quadrotor aircraft | |
CN112631316B (zh) | 变负载四旋翼无人机的有限时间控制方法 | |
CN105116905A (zh) | 一种飞行器姿态控制方法 | |
CN115576341A (zh) | 基于函数微分和自适应变增益的无人机轨迹跟踪控制方法 | |
CN112327926A (zh) | 一种无人机编队的自适应滑模控制方法 | |
Mou et al. | More detailed disturbance measurement and active disturbance rejection altitude control for a flapping wing robot under internal and external disturbances | |
CN106647264A (zh) | 一种基于控制约束的扩展鲁棒h∞的无人机控制方法 | |
CN115129072A (zh) | 固定翼无人机位置跟踪偏差约束下终端滑模控制方法 | |
Wu et al. | Integrated optimization design using improved pigeon-inspired algorithm for a hypersonic vehicle model |
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 |