CN111121788B - 基于双矢量姿态基准的航天器姿态奇异确定方法及系统 - Google Patents
基于双矢量姿态基准的航天器姿态奇异确定方法及系统 Download PDFInfo
- Publication number
- CN111121788B CN111121788B CN201911274943.2A CN201911274943A CN111121788B CN 111121788 B CN111121788 B CN 111121788B CN 201911274943 A CN201911274943 A CN 201911274943A CN 111121788 B CN111121788 B CN 111121788B
- Authority
- CN
- China
- Prior art keywords
- attitude
- spacecraft
- vector
- axis
- dual
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Navigation (AREA)
Abstract
本发明提供了一种基于双矢量姿态基准的航天器姿态奇异确定方法及系统,包括双目标选择、航天器三轴姿态指向确定,航天器相对双目标运动关系计算、航天器三轴姿态变化规律计算和姿态奇异确定五个部分。通过空间中双目标的选择确定航天器双矢量姿态基准,进而确定航天器三轴姿态指向,根据航天器相对双目标运动关系确定航天器三轴姿态变化规律,最终确定航天器是否发生姿态奇异现象。本发明提出的方法简单,适用于高低轨航天器,对航天器在轨三轴姿态指向设计具有指导意义。
Description
技术领域
本发明涉及航天器姿态控制领域,具体地,涉及一种基于双矢量姿态基准的航天器姿态奇异确定方法及系统。
背景技术
航天器姿态控制方式有多种,常用的是三轴稳定姿态控制方式。其方法是限制两轴的指向,利用右手法则确定航天器三轴姿态,保证航天器三轴姿态稳定。大部分在轨航天器一轴指向地心,另一轴指向轨道面的负法线方向,即以地心矢量和轨道面负法线矢量作为姿态基准,从而保持航天器对地指向状态,保证航天器测控能源等需求。
受轨道、测控数传等因素影响,一些航天器并不是根据地心和轨道负法线矢量确定三轴姿态。考虑能源的充分性,可将太阳作为双目标之一,即双矢量之一为太阳矢量。
目前已存在航天器将地心和太阳作为双目标,地心矢量和太阳矢量作为双矢量姿态基准,确定航天器的三轴姿态。因地心矢量和太阳矢量并不是始终垂直,因此由双矢量确定三轴姿态需要进行叉乘运算,使得某一轴与地心矢量垂直,从而确定出三轴姿态。当地心矢量和太阳矢量平行时,叉乘运算无效,当地心矢量和太阳矢量平行处波动时,会引起双矢量确定的三轴姿态指向发生180°翻转现象,即发生姿态奇异现象。
航天器发生姿态奇异的后果严重,一是影响姿态指向精度,奇异严重时造成姿态严重超差,影响卫星姿态和能源安全;二是打乱航天器在轨一系列任务的执行,造成严重的损失。
因此,航天器是否会发生姿态奇异,什么条件下发生姿态奇异,是一个必要的研究课题。
发明专利“大椭圆轨道姿态基准确定方法”(公开号:CN106767811A)提出了采用地心矢量和太阳矢量作为双矢量,但是没有涉及双目标选择、三轴姿态变化规律计算和姿态奇异确定问题。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种基于双矢量姿态基准的航天器姿态奇异确定方法及系统。
根据本发明提供的一种基于双矢量姿态基准的航天器姿态奇异确定方法,包括:
双目标选择步骤:选择空间中的双目标M1,M2,以航天器为矢量起始点,双目标为矢量终止点,确定空间中的双矢量r1=[x1 y1 z1]与r2=[x2 y2 z2]作为姿态基准,其中x1,y1,z1分别为矢量r1在航天器本体系下的三轴坐标,x2,y2,z2分别为矢量r2在航天器本体系下的三轴坐标;
运动关系计算步骤:计算双目标M1,M2相对航天器在惯性空间中的相对运动情况,计算航天器相对双目标的运动关系,体现为双矢量r1=[x1 y1 z1]与r2=[x2 y2 z2]随时间的变化规律r1(t)-[x1(t) y1(t) z1(t)]与r2(t)=[x2(t) y2(t) z2(t)];
姿态奇异确定步骤:根据双矢量随时间的变化规律r1(t),r2(t),确定是否会发生三轴姿态翻转的姿态奇异,若任意某一轴发生姿态翻转,则航天器姿态奇异,反之,航天器在轨无姿态奇异现象。
优选地,所述双矢量是以航天器为矢量起始点,双目标分别为矢量终止点。
优选地,所述运动关系计算步骤利用相对运动学原理计算双目标M1,M2相对航天器在惯性空间中的相对运动情况,计算航天器相对双目标的运动关系。
优选地,所述变化规律计算步骤中,利用规定的双矢量确定三轴姿态法则,确定航天器三轴姿态随时间的变化规律。
根据本发明提供的一种基于双矢量姿态基准的航天器姿态奇异确定系统,包括:
双目标选择模块:选择空间中的双目标M1,M2,以航天器为矢量起始点,双目标为矢量终止点,确定空间中的双矢量r1=[x1 y1 z1]与r2=[x2 y2 z2]作为姿态基准,其中x1,y1,z1分别为矢量r1在航天器本体系下的三轴坐标,x2,y2,z2分别为矢量r2在航天器本体系下的三轴坐标;
运动关系计算模块:计算双目标M1,M2相对航天器在惯性空间中的相对运动情况,计算航天器相对双目标的运动关系,体现为双矢量r1=[x1 y1 z1]与r2=[x2 y2 z2]随时间的变化规律r1(t)=[x1(t) y1(t) z1(t)]与r2(t)=[x2(t) y2(t) z2(t)];
姿态奇异确定模块:根据双矢量随时间的变化规律r1(t),r2(t),确定是否会发生三轴姿态翻转的姿态奇异,若任意某一轴发生姿态翻转,则航天器姿态奇异,反之,航天器在轨无姿态奇异现象。
优选地,所述双矢量是以航天器为矢量起始点,双目标分别为矢量终止点。
优选地,所述运动关系计算模块利用相对运动学原理计算双目标M1,M2相对航天器在惯性空间中的相对运动情况,计算航天器相对双目标的运动关系。
优选地,所述变化规律计算模块中,利用规定的双矢量确定三轴姿态法则,确定航天器三轴姿态随时间的变化规律。
与现有技术相比,本发明具有如下的有益效果:
本发明解决了双矢量姿态基准的航天器姿态奇异确定问题,方法简单,适用于高低轨航天器,对航天器在轨三轴姿态指向设计具有指导意义
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明的工作流程图;
图2为本发明实施例双矢量确定三轴姿态示意图;
图3为姿态奇异示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
下面结合图1说明本发明一种基于双矢量姿态基准的航天器姿态奇异确定方法的具体实施方式。内容包括五个部分:双目标选择、航天器三轴姿态指向确定,航天器相对双目标运动关系计算、航天器三轴姿态变化规律计算和姿态奇异确定。
以地心矢量和太阳矢量作为姿态基准为例:
1、双目标选择
选择地心和太阳作为空间中的双目标M1,M2,以航天器为矢量起始点,地心和太阳为矢量终止点,在地心惯性坐标系下,确定地心矢量Ra1(t)与太阳矢量Ra2(t)作为姿态基准。
2、航天器三轴姿态指向确定
综合考虑航天器帆板安装位置、测控条件、在轨任务特点等因素,利用右手法则,规定地心和太阳矢量姿态基准确定航天器三轴姿态指向的法则,Z轴指向为卫星质心指向地心,即地心矢量方向;Y轴指向为卫星质心指向太阳矢量和地心矢量的叉乘方向r2×r1,X、Y、Z轴构成右手直角坐标系。如附图2所示。
3、航天器相对双目标运动关系计算
设立初始时间,根据航天器的轨道根数和地球绕太阳公转规律,根据轨道根数计算实时航天器和太阳位置,从而得到航天器相对地心和太阳相对于地心位置随时间变化规律。根据航天器轨道六根数(a,e,i,Ω,ω,f),可解算航天器在J2000系下的位置如下:
太阳在J2000历元地心平赤道坐标系中的平均轨道可以采用下面的Kepler根数形式:
其中,AU表示天文单位,1AU=1.49597870*108km,d表示天数,T为J2000起算的儒略世纪数T=(JD-2451545.0)/36525,JD为儒略日数。类似的方法,可计算出太阳相对于地心的位置。
再通过矢量计算法则,确定地心矢量和太阳矢量Ra1(t),Ra2(t),即为r1=[x1 y1z1]与r2=[x2 y2 z2]双矢量随时间的变化规律r1(t)=Ra1(t),r2(t)=Ra2(t)。
4、航天器三轴姿态变化规律计算
5、姿态奇异确定
根据计算航天器相对双目标的运动关系所得到地心和太阳矢量随时间的变化规律的r1(t)=[x1(t) y1(t) z1(t)]与r2(t)=[x2(t) y2(t) z2(t)],确定姿态是否会发生奇异,若出现姿态突然翻转180°的情况,则判断航天器会发生姿态奇异,反之,航天器不会发生姿态奇异。对于确定的航天器轨道,地心矢量和太阳矢量夹角受升交点赤经Ω和太阳赤经Ωs的偏差ΔΩ=Ωs-Ω约束,易知当ΔΩ=0时,太阳矢量在轨道面内,航天器在轨道面内必定存在两点,使得地心矢量和太阳矢量平行,发生姿态奇异现象。
如附图3所示。太阳矢量固定不变,当地心矢量为Z1时,Y轴指向为Y1;当地心矢量为Z2时,Y轴指向为Y2。Y轴指向发生180°翻转,卫星姿态超差,即卫星发生姿态奇异现象。对于一个确定的卫星轨道,一轨内可能两段出现奇异,一段为双矢量夹角接近0°,一段为双矢量夹角接近180°。
本发明在航天器在轨姿态控制方法设计中,支撑解决可能出现的姿态奇异问题,方法简单,工程上易于实现。
在上述一种基于双矢量姿态基准的航天器姿态奇异确定方法的基础上,本发明还提供一种基于双矢量姿态基准的航天器姿态奇异确定系统,包括:
双目标选择模块:选择空间中的双目标M1,M2,以航天器为矢量起始点,双目标为矢量终止点,确定空间中的双矢量r1=[x1 y1 z1]与r2=[x2 y2 z2]作为姿态基准,其中x1,y1,z1分别为矢量r1在航天器本体系下的三轴坐标,x2,y2,z2分别为矢量r2在航天器本体系下的三轴坐标;
运动关系计算模块:计算双目标M1,M2相对航天器在惯性空间中的相对运动情况,计算航天器相对双目标的运动关系,体现为双矢量r1=[x1 y1 z1]与r2=[x2 y2 z2]随时间的变化规律r1(t)=[x1(t) y1(t) z1(t)]与r2(t)=[x2(t) y2(t) z2(t)];
姿态奇异确定模块:根据双矢量随时间的变化规律r1(t),r2(t),确定是否会发生三轴姿态翻转的姿态奇异,若任意某一轴发生姿态翻转,则航天器姿态奇异,反之,航天器在轨无姿态奇异现象。
本领域技术人员知道,除了以纯计算机可读程序代码方式实现本发明提供的系统及其各个装置、模块、单元以外,完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统及其各个装置、模块、单元以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同功能。所以,本发明提供的系统及其各项装置、模块、单元可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置、模块、单元也可以视为硬件部件内的结构;也可以将用于实现各种功能的装置、模块、单元视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (8)
1.一种基于双矢量姿态基准的航天器姿态奇异确定方法,其特征在于,包括:
双目标选择步骤:选择空间中的双目标M1,M2,以航天器为矢量起始点,双目标为矢量终止点,确定空间中的双矢量r1=[x1 y1 z1]与r2=[x2 y2 z2]作为姿态基准,其中x1,y1,z1分别为矢量r1在航天器本体系下的三轴坐标,x2,y2,z2分别为矢量r2在航天器本体系下的三轴坐标;
运动关系计算步骤:计算双目标M1,M2相对航天器在惯性空间中的相对运动情况,计算航天器相对双目标的运动关系,体现为双矢量r1=[x1 y1 z1]与r2=[x2 y2 z2]随时间的变化规律r1(t)=[x1(t) y1(t) z1(t)]与r2(t)=[x2(t) y2(t) z2(t)];
姿态奇异确定步骤:根据双矢量随时间的变化规律r1(t)=Ra1(t),r2(t)=Ra2(t),确定是否会发生三轴姿态翻转的姿态奇异,若任意某一轴发生姿态翻转,则航天器姿态奇异,反之,航天器在轨无姿态奇异现象。
3.根据权利要求1所述的基于双矢量姿态基准的航天器姿态奇异确定方法,其特征在于,所述运动关系计算步骤利用相对运动学原理计算双目标M1,M2相对航天器在惯性空间中的相对运动情况,计算航天器相对双目标的运动关系。
4.根据权利要求1所述的基于双矢量姿态基准的航天器姿态奇异确定方法,其特征在于,所述变化规律计算步骤中,利用规定的双矢量确定三轴姿态法则,确定航天器三轴姿态随时间的变化规律。
5.一种基于双矢量姿态基准的航天器姿态奇异确定系统,其特征在于,包括:
双目标选择模块:选择空间中的双目标M1,M2,以航天器为矢量起始点,双目标为矢量终止点,确定空间中的双矢量r1=[x1 y1 z1]与r2=[x2 y2 z2]作为姿态基准,其中x1,y1,z1分别为矢量r1在航天器本体系下的三轴坐标,x2,y2,z2分别为矢量r2在航天器本体系下的三轴坐标;
运动关系计算模块:计算双目标M1,M2相对航天器在惯性空间中的相对运动情况,计算航天器相对双目标的运动关系,体现为双矢量r1=[x1 y1 z1]与r2=[x2 y2 z2]随时间的变化规律r1(t)=[x1(t) y1(t) z1(t)]与r2(t)=[x2(t) y2(t) z2(t)];
姿态奇异确定模块:根据双矢量随时间的变化规律r1(t),r2(t),确定是否会发生三轴姿态翻转的姿态奇异,若任意某一轴发生姿态翻转,则航天器姿态奇异,反之,航天器在轨无姿态奇异现象。
7.根据权利要求5所述的基于双矢量姿态基准的航天器姿态奇异确定系统,其特征在于,所述运动关系计算模块利用相对运动学原理计算双目标M1,M2相对航天器在惯性空间中的相对运动情况,计算航天器相对双目标的运动关系。
8.根据权利要求5所述的基于双矢量姿态基准的航天器姿态奇异确定系统,其特征在于,所述变化规律计算模块中,利用规定的双矢量确定三轴姿态法则,确定航天器三轴姿态随时间的变化规律。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911274943.2A CN111121788B (zh) | 2019-12-12 | 2019-12-12 | 基于双矢量姿态基准的航天器姿态奇异确定方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911274943.2A CN111121788B (zh) | 2019-12-12 | 2019-12-12 | 基于双矢量姿态基准的航天器姿态奇异确定方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111121788A CN111121788A (zh) | 2020-05-08 |
CN111121788B true CN111121788B (zh) | 2021-07-06 |
Family
ID=70499928
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911274943.2A Active CN111121788B (zh) | 2019-12-12 | 2019-12-12 | 基于双矢量姿态基准的航天器姿态奇异确定方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111121788B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113525720B (zh) * | 2021-07-27 | 2022-01-28 | 北京微纳星空科技有限公司 | 一种航天器目标姿态的确定方法、装置、设备及存储介质 |
CN114526742B (zh) * | 2022-01-25 | 2024-05-07 | 上海卫星工程研究所 | 微纳卫星姿态基准的构件化通用构造方法及系统 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6138061A (en) * | 1997-04-18 | 2000-10-24 | Hughes Electronics Corporation | Onboard orbit propagation using quaternions |
CN1876501A (zh) * | 2006-05-31 | 2006-12-13 | 哈尔滨工业大学 | 基于行为模式的深空三轴稳定姿态定向控制方法 |
CN101236085A (zh) * | 2008-01-31 | 2008-08-06 | 北京控制工程研究所 | 奔月段基于紫外敏感器、星敏感器的自主轨道估计方法 |
EP2472225A2 (en) * | 2010-12-22 | 2012-07-04 | Custom Sensors & Technologies, Inc. | Method and system for initial quaternion and attitude estimation |
CN103646127A (zh) * | 2013-11-20 | 2014-03-19 | 中国空间技术研究院 | 卫星轨道姿态可视化三维显示方法 |
CN104567881A (zh) * | 2014-12-26 | 2015-04-29 | 北京控制工程研究所 | 一种基于太阳、地心及月表重力矢量的巡视器位置姿态确定方法 |
CN104833358A (zh) * | 2015-05-13 | 2015-08-12 | 上海交通大学 | 基于罗德里格斯参数的星敏感器几何线性姿态确定方法 |
CN106323271A (zh) * | 2016-08-01 | 2017-01-11 | 北京理工大学 | 基于特征奇异值的航天器相对姿态测量矢量选取方法 |
CN106767811A (zh) * | 2016-11-22 | 2017-05-31 | 上海卫星工程研究所 | 大椭圆轨道姿态基准确定方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9091552B2 (en) * | 2011-10-25 | 2015-07-28 | The Boeing Company | Combined location and attitude determination system and methods |
-
2019
- 2019-12-12 CN CN201911274943.2A patent/CN111121788B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6138061A (en) * | 1997-04-18 | 2000-10-24 | Hughes Electronics Corporation | Onboard orbit propagation using quaternions |
CN1876501A (zh) * | 2006-05-31 | 2006-12-13 | 哈尔滨工业大学 | 基于行为模式的深空三轴稳定姿态定向控制方法 |
CN101236085A (zh) * | 2008-01-31 | 2008-08-06 | 北京控制工程研究所 | 奔月段基于紫外敏感器、星敏感器的自主轨道估计方法 |
EP2472225A2 (en) * | 2010-12-22 | 2012-07-04 | Custom Sensors & Technologies, Inc. | Method and system for initial quaternion and attitude estimation |
CN103646127A (zh) * | 2013-11-20 | 2014-03-19 | 中国空间技术研究院 | 卫星轨道姿态可视化三维显示方法 |
CN104567881A (zh) * | 2014-12-26 | 2015-04-29 | 北京控制工程研究所 | 一种基于太阳、地心及月表重力矢量的巡视器位置姿态确定方法 |
CN104833358A (zh) * | 2015-05-13 | 2015-08-12 | 上海交通大学 | 基于罗德里格斯参数的星敏感器几何线性姿态确定方法 |
CN106323271A (zh) * | 2016-08-01 | 2017-01-11 | 北京理工大学 | 基于特征奇异值的航天器相对姿态测量矢量选取方法 |
CN106767811A (zh) * | 2016-11-22 | 2017-05-31 | 上海卫星工程研究所 | 大椭圆轨道姿态基准确定方法 |
Non-Patent Citations (1)
Title |
---|
皮卫星姿态确定与控制系统方案设计与实现;蒙涛;《中国优秀博硕士学位论文全文数据库(硕士) 工程科技II辑》;20060915;第1-99页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111121788A (zh) | 2020-05-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Roscoe et al. | Overview and GNC design of the CubeSat Proximity Operations Demonstration (CPOD) mission | |
US9643740B2 (en) | Gyroless three-axis sun acquisition using sun sensor and unscented kalman filter | |
CN109592079A (zh) | 一种限定时间的航天器共面交会变轨策略确定方法 | |
EP3379366B1 (en) | Methods and apparatus to minimize command dynamics of a satellite | |
Zhang et al. | Characteristic analysis for elliptical orbit hovering based on relative dynamics | |
CN111121788B (zh) | 基于双矢量姿态基准的航天器姿态奇异确定方法及系统 | |
Abdelrahman et al. | Sigma-point Kalman filtering for spacecraft attitude and rate estimation using magnetometer measurements | |
Thienel et al. | Accurate state estimation and tracking of a non-cooperative target vehicle | |
CN113091753B (zh) | 用于星敏视场保护的卫星姿态导引方法及其系统 | |
Xie et al. | Guidance, navigation, and control for spacecraft rendezvous and docking: theory and methods | |
Abdollahzadeh et al. | Automatic orbital docking with tumbling target using sliding mode control | |
Wu et al. | Multi-objective output-feedback control for microsatellite attitude control: An LMI approach | |
CN111123961B (zh) | 基于约束分析的双矢量夹角限制范围确定方法及系统 | |
CN116692028B (zh) | 一种小卫星对地快速凝视指向跟踪控制方法及装置 | |
Noumi et al. | Fault-tolerant attitude control systems for satellite equipped with control moment gyros | |
Okasha et al. | Relative motion guidance, navigation and control for autonomous orbital rendezvous | |
CN114802818A (zh) | 晨昏轨道卫星及其对日姿态计算方法、导引方法 | |
Hajiyev et al. | Integration of algebraic method and EKF for attitude determination of small information satellites | |
Chasset et al. | 3-axis magnetic control with multiple attitude profile capabilities in the PRISMA mission | |
Karpenko et al. | Fast attitude maneuvers for the lunar reconnaissance orbiter | |
Okasha et al. | Modeling, dynamics and control of spacecraft relative motion in a perturbed Keplerian orbit | |
CN113247310B (zh) | 一种适用于卫星可连续姿态机动次数的估算方法及系统 | |
Mashtakov et al. | Study of the accuracy provided by small satellite attitude determination & control system | |
Liang et al. | Attitude tracking control based on adaptive sliding mode technique with double closed loop for spacecraft near small body | |
Newman et al. | Investigating solar radiation pressure modeling for operations in near rectilinear halo orbit |
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 |