CN112379599A - 一种针对空中机动目标的三维协同中制导方法 - Google Patents

一种针对空中机动目标的三维协同中制导方法 Download PDF

Info

Publication number
CN112379599A
CN112379599A CN202011314913.2A CN202011314913A CN112379599A CN 112379599 A CN112379599 A CN 112379599A CN 202011314913 A CN202011314913 A CN 202011314913A CN 112379599 A CN112379599 A CN 112379599A
Authority
CN
China
Prior art keywords
interceptor
guidance
target
interceptors
ith
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.)
Granted
Application number
CN202011314913.2A
Other languages
English (en)
Other versions
CN112379599B (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
Original Assignee
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 filed Critical Northwestern Polytechnical University
Priority to CN202011314913.2A priority Critical patent/CN112379599B/zh
Publication of CN112379599A publication Critical patent/CN112379599A/zh
Application granted granted Critical
Publication of CN112379599B publication Critical patent/CN112379599B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive 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/042Adaptive 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)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明公开了一种针对空中机动目标的三维协同中制导方法,将虚拟碰撞点作为拦截器中制导阶段的目标点,同时考虑利用在中制导阶段通过外部传感器直接获得的目标的位置信息、速度信息和加速度信息设计协同中制导方法,并通过添加协同调节项确保所有拦截器的在中制导阶段结束时,几乎同时到达中末交接班区域。本发明结构简单,物理含义清晰,易于工程人员理解,同时具有协同项系数便于调节,制导律鲁棒性强等优点,易于工程实现。

Description

一种针对空中机动目标的三维协同中制导方法
技术领域
本发明属于制导技术领域,具体涉及一种三维协同中制导方法。
背景技术
随着反拦截器技术日趋完善,空中机动目标一方面可以通过释放诱饵、蛇形机动或者其他对抗手段对拦截器进行干扰和对抗,使得现有拦截器的拦截效率大大降低;另一方面,空中目标的性能随着人工智能、无人机技术的快速发展而大幅提升,单个拦截器的探测能力和机动能力都难以满足高精度的打击要求,需要采用协同包围攻击技术。因此,多拦截器协同突防或拦截的思路应运而生。多拦截器协同编队协同作战方式是指多个拦截器通过信息共享、功能互补,不仅能够实现对目标从多个方向进行饱和攻击,提高打击效果和命中概率,还能实现单个拦截器不易完成的协同突防等任务。针对军事价值较高的目标,多拦截器协同编队作战方式是实施突发性和高密集的精确打击的重要手段,是未来突发性地区战争的重要打击形式。
对于中远程拦截器来说,要想实现多拦截器编队进行饱和打击需要经历四个飞行阶段,分别为编队集结、编队飞行、编队中末交班和编队末制导阶段。对于中远程拦截器编队来说,要想实现协同打击目标的任务,需要从拦截器发射阶段就开始进行协同,经过协同中制导和中末制导交班,最后进行协同末制导,对目标完成协同聚能饱和打击。若单研究其中一个飞行阶段的协同制导问题是不可行的,从而无法完成协同打击目标的任务,例如:当有两枚拦截器同时指向同一目标,一枚拦截器已经进入末制导阶段,距离目标低于15KM;另一枚拦截器在中制导阶段,距离目标有30KM以上,此时,这两枚拦截器很难完成协同攻击,不符合进行协同制导的基础条件(离目标的距离误差在一定范围内)。所以,对于中远程多拦截器编队协同制导,多拦截器编队整个飞行阶段的协同都至关重要,每一阶段的协同完成效果直接影响下一阶段的协同任务完成度。因此,对于多拦截器编队协同中制导段,需要实现可以在同一时刻或者近似同一时刻到达中末制导交班区,为多拦截器编队进行协同末制导做好准备。
目前,主要的理论研究与工程方法大都集中在协同末制导段,而末制导阶段的协同与中制导阶段的协同有本质上的区别。末制导阶段的目标信息获取由拦截器本身的导引头完成,目标信息精度较高,而且多拦截器编队队形和通信拓扑结构基本不发生变化,拦截器之间距离较小,通信连接情况良好,很容易达到强连通状态。在中制导阶段,目标信息由外部传感器探测和估计所得,拦截器之间距离较大,通信连接情况较差。因此,拦截器在协同中制导段所采用的制导律与在协同末制导段采用的制导律有很大区别,主要体现在:性能指标的要求、约束条件、制导目的以及获取目标信息的方式和目标信息的质量等。
同时,因为末制导时间太短,拦截器可以调整的飞行距离和飞行速度等非常有限,必须要有非常好的协同中制导律,确保在中末制导交班时就能达到时间和空间初步的协同,即给末制导段提供非常好的协同初始条件,才能完成末制导段的协同,从而实现中远程拦截器的时间和空间协同。因此,协同中制导方法是中远程拦截器实现协同的前提和基础。
目前,对协同中制导的理论研究成果鲜有报道,已经严重影响了多拦截器编队协同控制系统的研制。因此,提出一套适合工程应用的多拦截器协同中制导方法就显得非常迫切。
发明内容
为了克服现有技术的不足,本发明提供了一种针对空中机动目标的三维协同中制导方法,将虚拟碰撞点作为拦截器中制导阶段的目标点,同时考虑利用在中制导阶段通过外部传感器直接获得的目标的位置信息、速度信息和加速度信息设计协同中制导方法,并通过添加协同调节项确保所有拦截器的在中制导阶段结束时,几乎同时到达中末交接班区域。本发明结构简单,物理含义清晰,易于工程人员理解,同时具有协同项系数便于调节,制导律鲁棒性强等优点,易于工程实现。
本发明解决其技术问题所采用的技术方案包括以下步骤:
步骤1:建立多拦截器和目标的运动模型:
Figure BDA0002791053000000021
Figure BDA0002791053000000022
Figure BDA0002791053000000031
Figure BDA0002791053000000032
其中,(xi,yi,zi)和(xT,yT,zT)分别代表第i个拦截器和目标在惯性坐标系三个坐标方向上的位置信息,(Vxi,Vyi,Vzi)和(VTx,VTy,VTz)分别代表第i个拦截器和目标在惯性坐标系三个坐标方向上的速度信息,(axi,ayi,azi)和(aTx,aTy,aTz)分别代表第i个拦截器和目标在惯性坐标系三个坐标方向上的加速度信息,i=1,2,...,N,N代表拦截器的个数;
步骤2:计算多拦截器到达目标的剩余时间tgoi
Figure BDA0002791053000000033
其中ri代表第i个拦截器到达目标的剩余距离,计算公式如下:
Figure BDA0002791053000000034
步骤3:计算虚拟碰撞点;
设置虚拟碰撞点作为中制导阶段拦截器的目标点,虚拟碰撞点的计算如下:
Figure BDA0002791053000000035
其中,
Figure BDA0002791053000000036
代表第i个拦截器的虚拟碰撞点;
步骤4:计算三维惯性坐标系下拦截器加速度为:
Figure BDA0002791053000000037
其中,tf代表中制导结束时间,(VTx(tf),VTy(tf),VTz(tf))和(Vxi(tf),Vyi(tf),Vzi(tf))分别代表目标和拦截器在惯性坐标系下三个方向的末速度,(Pxi,Pyi,Pzi)代表协同项调节系数,
Figure BDA0002791053000000041
代表所有拦截器虚拟碰撞点坐标的平均值,计算如下:
Figure BDA0002791053000000042
步骤5:设置协同项调节系数(Pxi,Pyi,Pzi)为大于0的任意正实数,从而三维协同中制导律满足一致性要求。
本发明提出的一种针对空中机动目标的三维协同中制导方法,能够确保中远程拦截器编队几乎同时到达中末交接班区域,为拦截器编队的协同末制导阶段提供良好的初始条件。本发明方法结构简单,物理含义清晰,易于工程人员理解,同时具有协同项系数便于调节,制导律鲁棒性强等优点,易于工程实现。
附图说明
图1是本发明多个拦截器与目标在惯性系下的三维示意图。
图2是本发明实施例各拦截器与目标的运动轨迹示意图。
图3是本发明实施例各拦截器的剩余时间曲线示意图。
图4是本发明实施例各拦截器的剩余距离曲线示意图。
图5是本发明实施例拦截器1的加速度指令曲线示意图。
图6是本发明实施例拦截器2的加速度指令曲线示意图。
图7是本发明实施例拦截器3的加速度指令曲线示意图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
如图1所示,一种针对空中机动目标的三维协同中制导方法,包括如下步骤:
步骤1:建立多拦截器和目标的运动模型:
Figure BDA0002791053000000043
Figure BDA0002791053000000051
Figure BDA0002791053000000052
Figure BDA0002791053000000053
其中,(xi,yi,zi)和(xT,yT,zT)分别代表第i个拦截器和目标在惯性坐标系三个坐标方向上的位置信息,(Vxi,Vyi,Vzi)和(VTx,VTy,VTz)分别代表第i个拦截器和目标在惯性坐标系三个坐标方向上的速度信息,(axi,ayi,azi)和(aTx,aTy,aTz)分别代表第i个拦截器和目标在惯性坐标系三个坐标方向上的加速度信息,i=1,2,...,N,N代表拦截器的个数;
步骤2:计算多拦截器到达目标的剩余时间tgoi
Figure BDA0002791053000000054
其中ri代表第i个拦截器到达目标的剩余距离,计算公式如下:
Figure BDA0002791053000000055
步骤3:计算虚拟碰撞点;
为了提高协同末制导对机动目标的拦截成功率,设置虚拟碰撞点作为中制导阶段拦截器的目标点,虚拟碰撞点的计算如下:
Figure BDA0002791053000000056
其中,
Figure BDA0002791053000000057
代表第i个拦截器的虚拟碰撞点;
步骤4:计算三维惯性坐标系下拦截器加速度为:
Figure BDA0002791053000000061
其中,tf代表中制导结束时间,(VTx(tf),VTy(tf),VTz(tf))和(Vxi(tf),Vyi(tf),Vzi(tf))分别代表目标和拦截器在惯性坐标系下三个方向的末速度,(Pxi,Pyi,Pzi)代表协同项调节系数,
Figure BDA0002791053000000062
代表所有拦截器虚拟碰撞点坐标的平均值,计算如下:
Figure BDA0002791053000000063
步骤5:为了使协同中制导律达到协同一致性,协同项调节系数(Pxi,Pyi,Pzi)为大于0的任意正实数。
以X轴为例进行分析,Y轴和Z轴证明类似,具体证明如下:
由一致性条件
Figure BDA0002791053000000064
即所有拦截器的虚拟碰撞点在协同中制导阶段会趋于一致。因此,通过建立函数
Figure BDA0002791053000000065
将一致性问题转化为当ηi(t+△t)<ηi(t)即ηi(t)单调递减时,则ηi(t)→0,从而
Figure BDA0002791053000000066
接下来分析确保不等式ηi(t+△t)<ηi(t)的具体条件:
由于
Figure BDA0002791053000000067
其中,
Figure BDA00027910530000000612
Bi=MiVTx
Figure BDA0002791053000000069
Figure BDA00027910530000000610
代表所有
Figure BDA00027910530000000611
的平均值。
Figure BDA0002791053000000071
由(10)和(11)式可得:
Figure BDA0002791053000000072
进一步可将axi改写为axi=abi-Pxiaci,其中:
Figure BDA0002791053000000073
由不等式ηi(t+△t)<ηi(t),△t>0可推出:
Figure BDA0002791053000000074
其中,Oi=NAi
Figure BDA0002791053000000075
Wi=Bi(2-Di(Ei-(xT-xi)abi))。
接下来分析不等式(14)右边项是小于0的,从而Pi选择任意大于0的实数可满足ηi(t+△t)<ηi(t)。
由于一般情况下,VTx>0且xT>xi,可得:
Figure BDA0002791053000000076
因为
Figure BDA0002791053000000077
可得:
Figure BDA0002791053000000078
Bi=MiVTx<0 (17)
在拦截器编队飞行过程中,
Figure BDA0002791053000000079
因此可得
Figure BDA00027910530000000710
(Ei-(xT-xi)axi)Di<2 (19)
由于abi是axi的一部分,可以近似处理为axi≈abi,因而获得:
(Ei-(xT-xmi)abi)Di<2 (20)
Wi=Bi(2-(Ei-(xT-xmi)abi)Di)<0 (21)
从(16)、(18)和(21)式可获得
Figure BDA0002791053000000081
因此(Pxi,Pyi,Pzi)可以选取任意大于0的实数,确保ηi(t+△t)<ηi(t)成立,从而三维协同中制导律满足一致性要求。
本发明将虚拟碰撞点作为拦截器中制导阶段的目标点,同时考虑利用在中制导阶段通过外部传感器直接获得的目标的位置信息、速度信息和加速度信息设计协同中制导方法,并通过添加协同调节项确保所有拦截器的在中制导阶段结束时,几乎同时到达中末交接班区域。
具体实施例:
本实施例通过用3个拦截器同时攻击机动目标来验证所提出的三维协同中制导导引律的有效性。实施例中3个拦截器和目标的仿真初始条件:
拦截器1:初始位置(1000m,9000m,-2300m),初始速度(1100m/s,100m/s,100m/s)
拦截器2:初始位置(500m,8500m,-500m),初始速度(1000m/s,100m/s,90m/s)
拦截器3:初始位置(2000m,7600m,3000m),初始速度(950m/s,60m/s,110m/s)
目标:初始位置(100000m,8000m,1000m),初始速度(150m/s,10m/s,15m/s)
选取
Figure BDA0002791053000000082
显然(pxi,pyi,pzi)是大于0的正实数。
图2为在三维平面下多拦截器与目标的轨迹,其中黑色实线为拦截器在协同中制导的轨迹,黑色点实线为拦截器为协同末制导阶段的轨迹。图3为所有拦截器的飞行剩余时间变化曲线。图4为所有拦截器与目标的剩余距离的变化曲线。图5到图7分别表示各拦截器的加速度指令变化曲线。从图2可以看出,本发明方法能确保所有拦截器几乎同时到达中末交班区域,为协同末制导提供了良好的初始条件,也为最终拦截机动目标提供的拦截环境。从图3到图4可以看出,本发明方法能确保所有拦截器的剩余时间与剩余距离在协同中制导阶段趋近收敛,达到一致的状态。从图5到图7可以看出本发明方法的加速度指令平滑,且数值较小,易于工程实现。

Claims (1)

1.一种针对空中机动目标的三维协同中制导方法,其特征在于,包括以下步骤:
步骤1:建立多拦截器和目标的运动模型:
Figure FDA0002791052990000011
Figure FDA0002791052990000012
Figure FDA0002791052990000013
Figure FDA0002791052990000014
其中,(xi,yi,zi)和(xT,yT,zT)分别代表第i个拦截器和目标在惯性坐标系三个坐标方向上的位置信息,(Vxi,Vyi,Vzi)和(VTx,VTy,VTz)分别代表第i个拦截器和目标在惯性坐标系三个坐标方向上的速度信息,(axi,ayi,azi)和(aTx,aTy,aTz)分别代表第i个拦截器和目标在惯性坐标系三个坐标方向上的加速度信息,i=1,2,...,N,N代表拦截器的个数;
步骤2:计算多拦截器到达目标的剩余时间tgoi
Figure FDA0002791052990000015
其中ri代表第i个拦截器到达目标的剩余距离,计算公式如下:
Figure FDA0002791052990000016
步骤3:计算虚拟碰撞点;
设置虚拟碰撞点作为中制导阶段拦截器的目标点,虚拟碰撞点的计算如下:
Figure FDA0002791052990000017
其中,
Figure FDA0002791052990000021
代表第i个拦截器的虚拟碰撞点;
步骤4:计算三维惯性坐标系下拦截器加速度为:
Figure FDA0002791052990000022
其中,tf代表中制导结束时间,(VTx(tf),VTy(tf),VTz(tf))和(Vxi(tf),Vyi(tf),Vzi(tf))分别代表目标和拦截器在惯性坐标系下三个方向的末速度,(Pxi,Pyi,Pzi)代表协同项调节系数,
Figure FDA0002791052990000023
代表所有拦截器虚拟碰撞点坐标的平均值,计算如下:
Figure FDA0002791052990000024
步骤5:设置协同项调节系数(Pxi,Pyi,Pzi)为大于0的任意正实数,从而三维协同中制导律满足一致性要求。
CN202011314913.2A 2020-11-20 2020-11-20 一种针对空中机动目标的三维协同中制导方法 Active CN112379599B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011314913.2A CN112379599B (zh) 2020-11-20 2020-11-20 一种针对空中机动目标的三维协同中制导方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011314913.2A CN112379599B (zh) 2020-11-20 2020-11-20 一种针对空中机动目标的三维协同中制导方法

Publications (2)

Publication Number Publication Date
CN112379599A true CN112379599A (zh) 2021-02-19
CN112379599B CN112379599B (zh) 2022-06-14

Family

ID=74589015

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011314913.2A Active CN112379599B (zh) 2020-11-20 2020-11-20 一种针对空中机动目标的三维协同中制导方法

Country Status (1)

Country Link
CN (1) CN112379599B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113433428A (zh) * 2021-05-10 2021-09-24 广西电网有限责任公司玉林供电局 一种多元混合绝缘气体的协同效应分析方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106843265A (zh) * 2016-12-30 2017-06-13 哈尔滨工业大学 有限时间收敛的三维多导弹协同制导方法与系统
CN108168381A (zh) * 2018-01-04 2018-06-15 北京理工大学 一种多枚导弹协同作战的控制方法
CN110187640A (zh) * 2019-06-29 2019-08-30 东南大学 针对机动目标和允许通信时滞的多导弹协同作战制导律设计方法
CN110412874A (zh) * 2019-07-25 2019-11-05 东南大学 针对机动目标和时延通信的多导弹协同制导律设计方法
CN111026137A (zh) * 2019-12-26 2020-04-17 西北工业大学 一种攻击角约束下同时攻击目标的三维分布式协同导引方法
CN111506101A (zh) * 2019-10-21 2020-08-07 北京理工大学 基于通信网络拓扑结构的飞行器协同制导控制方法及系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106843265A (zh) * 2016-12-30 2017-06-13 哈尔滨工业大学 有限时间收敛的三维多导弹协同制导方法与系统
CN108168381A (zh) * 2018-01-04 2018-06-15 北京理工大学 一种多枚导弹协同作战的控制方法
CN110187640A (zh) * 2019-06-29 2019-08-30 东南大学 针对机动目标和允许通信时滞的多导弹协同作战制导律设计方法
CN110412874A (zh) * 2019-07-25 2019-11-05 东南大学 针对机动目标和时延通信的多导弹协同制导律设计方法
CN111506101A (zh) * 2019-10-21 2020-08-07 北京理工大学 基于通信网络拓扑结构的飞行器协同制导控制方法及系统
CN111026137A (zh) * 2019-12-26 2020-04-17 西北工业大学 一种攻击角约束下同时攻击目标的三维分布式协同导引方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LIU ZHAOHUI 等: "Cooperative Guidance Law Design on Simultaneous Attack for Multiple Missiles under Time-Delayed Communication Topologies", 《 2019 IEEE SYMPOSIUM SERIES ON COMPUTATIONAL INTELLIGENCE (SSCI)》 *
ZIHAO WU,等: "Three-Dimensional Cooperative Mid-Course Guidance Law Against the Maneuvering Target", 《IEEE ACCESS》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113433428A (zh) * 2021-05-10 2021-09-24 广西电网有限责任公司玉林供电局 一种多元混合绝缘气体的协同效应分析方法
CN113433428B (zh) * 2021-05-10 2023-09-01 广西电网有限责任公司玉林供电局 一种多元混合绝缘气体的协同效应分析方法

Also Published As

Publication number Publication date
CN112379599B (zh) 2022-06-14

Similar Documents

Publication Publication Date Title
CN111580556B (zh) 一种时空约束下多无人机协同路径规划与制导方法
CN103884237B (zh) 基于目标概率分布信息的多对一协同制导方法
CN106529073B (zh) 基于拦截几何的高超声速目标拦截弹交接班条件分析方法
CN110412874B (zh) 针对机动目标和时延通信的多导弹协同制导律设计方法
CN108534614B (zh) 一种三维全向实时预测制导方法
CN110764523B (zh) 基于抗饱和平滑变换的比例积分前置导引攻击目标的方法
CN114138000B (zh) 考虑全捷联导引头视场约束的弹群协同制导控制一体化设计方法
CN110345814B (zh) 一种不依赖自身导引头测量信息的末制导算法
CN112648886B (zh) 一种组合制导目标拦截方法及系统
CN114153143B (zh) 一种导弹非奇异固定时间滑模制导律的设计方法
CN112379600A (zh) 一种基于通信时变延迟的分布式协同制导律构建方法
CN115857538A (zh) 三维空间下满足落角约束的多飞行器协同制导方法
CN112379599B (zh) 一种针对空中机动目标的三维协同中制导方法
CN112033235A (zh) 一种采用黄金分割搜索的空战三体对抗最优预测制导方法
CN109543135B (zh) 一种基于信息最大化的多auv协同同时导引方法
CN114489101A (zh) 一种无人机的末端制导控制方法及系统
CN117932894A (zh) 一种多模导引头信息融合目标状态估计方法
CN116774714A (zh) 一种基于事件触发机制的多约束协同制导方法
CN116225048A (zh) 基于落角和时间约束的巡飞弹末端协同打击方法
CN115542939A (zh) 空空导弹分布式协同中制导律分析方法及其导引系统
Tang et al. Research on Target State Estimation and Terminal Guidance Algorithm in the Process of Multi-UAV Cooperative Attack
CN115617063A (zh) 一种带落角约束的飞行器制导控制装置及方法
CN113759954A (zh) 一种针对机动目标的复合制导方法
CN113138604A (zh) 一种反低小慢目标的无人机拦截方法
CN115755955B (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