CN111591472A - 一种调整卫星姿态的方法和相关装置 - Google Patents
一种调整卫星姿态的方法和相关装置 Download PDFInfo
- Publication number
- CN111591472A CN111591472A CN202010413058.4A CN202010413058A CN111591472A CN 111591472 A CN111591472 A CN 111591472A CN 202010413058 A CN202010413058 A CN 202010413058A CN 111591472 A CN111591472 A CN 111591472A
- Authority
- CN
- China
- Prior art keywords
- orientation
- vector
- satellite model
- satellite
- adjusting
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 64
- 230000007547 defect Effects 0.000 abstract description 6
- 230000009471 action Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000005096 rolling process Methods 0.000 description 4
- 230000000007 visual effect Effects 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000002349 favourable effect Effects 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
-
- 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)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Aviation & Aerospace Engineering (AREA)
- Astronomy & Astrophysics (AREA)
- General Physics & Mathematics (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本申请公开了一种调整卫星姿态的方法和相关装置,该方法包括:依据欧拉角和J2000坐标系下X轴的朝向向量计算卫星模型的头朝向目标向量;按照卫星模型的头朝向当前向量和头朝向目标向量调整卫星模型的头朝向;依据欧拉角和J2000坐标系下Y轴的朝向向量计算卫星模型的帆板朝向目标向量;按照卫星模型的帆板朝向当前向量和帆板朝向目标向量调整卫星模型的帆板朝向。计算卫星模型的头朝向目标向量和帆板朝向目标向量,在卫星模型的头朝向当前向量和帆板朝向当前向量基础上,采用向量法调整卫星模型的头朝向和帆板朝向,能够屏蔽基于欧拉角和/或四元数调整卫星姿态的缺点;先调整卫星模型的头朝向再调整帆板朝向,有效避免万向节锁死。
Description
技术领域
本申请涉及卫星数据处理技术领域,尤其涉及一种调整卫星姿态的方法和相关装置。
背景技术
卫星姿态是指卫星星体在轨道上运行所处的空间指向状态。直角坐标系的原点位于卫星星体上,指向地面的Z轴反映偏航方向,Y轴反映俯仰方向,X轴反映滚转方向。卫星星体在轨道上运行过程中,可以切换切向模式、对地模式、对日模式、对月模式等多种模式,在卫星星体对应的卫星模型切换上述多种模式之后,可能会产生对X轴、Y轴和Z轴三轴的偏移,为了能够保持切换模式后卫星模型的卫星姿态是正确的,需要基于对X轴、Y轴和Z轴三轴的偏移调整卫星模型的卫星姿态。
一般的,调整物体姿态使用欧拉角和/或四元数,即,利用欧拉角和/或四元数调整卫星姿态。但是,在使用欧拉角需要转换左右手系时,对于卫星模型的不同轴向来说转换公式不一样,则利用欧拉角调整卫星姿态的方式较为复杂繁琐,且容易出现万向节锁死的问题,无法轻易解决;四元数表示四维向量,对于卫星模型的三维向量而言,利用四元数调整卫星姿态不够直观,不利于想象和观察。
发明内容
有鉴于此,本申请实施例提供一种调整卫星姿态的方法和相关装置,能够屏蔽基于欧拉角和/或四元数调整卫星姿态的缺点,不仅调整方式直观形象,而且有效避免了万向节锁死的情况。
第一方面,本申请实施例提供了一种调整卫星姿态的方法,所述方法包括:
基于欧拉角和J2000坐标系下X轴的朝向向量,获得卫星模型的头朝向目标向量;所述欧拉角包括所述卫星模型的滚转角、俯仰角和偏航角;
基于所述卫星模型的头朝向当前向量和所述头朝向目标向量,调整所述卫星模型的头朝向;
基于所述欧拉角和所述J2000坐标系下Y轴的朝向向量,获得所述卫星模型的帆板朝向目标向量;
基于所述卫星模型的帆板朝向当前向量和所述帆板朝向目标向量,调整所述卫星模型的帆板朝向。
可选的,所述基于欧拉角和J2000坐标系下X轴的朝向向量,获得卫星模型的头朝向目标向量,具体为:
将所述J2000坐标系下X轴的朝向向量绕所述欧拉角旋转后的向量作为所述卫星模型的头朝向目标向量。
可选的,所述基于所述卫星模型的头朝向当前向量和所述头朝向目标向量,调整所述卫星模型的头朝向,包括:
基于所述头朝向当前向量和所述头朝向目标向量,获得所述头朝向当前向量与所述头朝向目标向量之间的第一夹角角度;
根据所述第一夹角角度直接调整所述卫星模型的头朝向。
可选的,所述基于所述欧拉角和所述J2000坐标系下Y轴的朝向向量,获得所述卫星模型的帆板朝向目标向量,具体为:
将所述J2000坐标系下Y轴的朝向向量绕所述欧拉角旋转后的向量作为所述卫星模型的帆板朝向目标向量。
可选的,所述基于所述卫星模型的帆板朝向当前向量和所述帆板朝向目标向量,调整所述卫星模型的帆板朝向,包括:
基于所述帆板朝向当前向量和所述帆板朝向目标向量,获得所述帆板朝向当前向量与所述帆板朝向目标向量之间的第二夹角角度;
基于所述卫星模型的Y轴位置值和所述第二夹角角度,调整所述卫星模型的帆板朝向。
可选的,所述基于所述卫星模型的Y轴位置值和所述第二夹角角度,调整所述卫星模型的帆板朝向,具体为:
若所述卫星模型的Y轴位置值大于0,根据180°与所述第二夹角角度的差值角度调整所述卫星模型的帆板朝向。
可选的,所述基于所述卫星模型的Y轴位置值和所述第二夹角角度,调整所述卫星模型的帆板朝向,具体为:
若所述卫星模型的Y轴位置值小于等于0,根据180°与所述第二夹角角度的加和角度调整所述卫星模型的帆板朝向。
第二方面,本申请实施例提供了一种调整卫星姿态的装置,所述装置包括:
第一获得单元,用于基于欧拉角和J2000坐标系下X轴的朝向向量,获得卫星模型的头朝向目标向量;所述欧拉角包括所述卫星模型的滚转角、俯仰角和偏航角;
第一调整单元,用于基于所述卫星模型的头朝向当前向量和所述头朝向目标向量,调整所述卫星模型的头朝向;
第二获得单元,用于基于所述欧拉角和所述J2000坐标系下Y轴的朝向向量,获得所述卫星模型的帆板朝向目标向量;
第二调整单元,用于基于所述卫星模型的帆板朝向当前向量和所述帆板朝向目标向量,调整所述卫星模型的帆板朝向。
第三方面,本申请实施例提供了一种终端设备,所述终端设备包括处理器以及存储器:
所述存储器用于存储程序代码,并将所述程序代码传输给所述处理器;
所述处理器用于根据所述程序代码中的指令执行上述第一方面任一项所述的调整卫星姿态的方法。
第四方面,本申请实施例提供了一种计算机可读存储介质,所述计算机可读存储介质用于存储程序代码,所述程序代码用于执行上述第一方面任一项所述的调整卫星姿态的方法。
与现有技术相比,本申请至少具有以下优点:
采用本申请实施例的技术方案,首先,依据欧拉角和J2000坐标系下X轴的朝向向量计算卫星模型的头朝向目标向量;欧拉角包括卫星模型的滚转角、俯仰角和偏航角;其次,按照卫星模型的头朝向当前向量和头朝向目标向量调整卫星模型的头朝向;然后,依据欧拉角和J2000坐标系下Y轴的朝向向量计算卫星模型的帆板朝向目标向量;最后,按照卫星模型的帆板朝向当前向量和帆板朝向目标向量调整卫星模型的帆板朝向。由此可见,先计算卫星模型的头朝向目标向量,在卫星模型的头朝向当前向量的基础上,采用向量法调整卫星模型的头朝向,再计算卫星模型的帆板朝向目标向量,在卫星模型的帆板朝向当前向量的基础上,采用向量法调整卫星模型的帆板朝向,向量法能够屏蔽基于欧拉角和/或四元数调整卫星姿态的缺点;先调整卫星模型的头朝向再调整卫星模型的帆板朝向,不仅调整方式直观形象,而且有效避免了万向节锁死的情况。
附图说明
为了更清楚地说明本申请实施例的技术方案,下面将对本申请实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1为本申请实施例中一种应用场景所涉及的系统框架示意图;
图2为本申请实施例提供的一种调整卫星姿态的方法的流程示意图;
图3为本申请实施例提供的一种卫星姿态调整前卫星模型示意图;
图4为本申请实施例提供的一种调整卫星姿态的装置的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
目前,在卫星星体对应的卫星模型切换切向模式、对地模式、对日模式、对月模式等多种模式之后,为了能够保持切换模式后卫星模型的卫星姿态是正确的,利用欧拉角和/或四元数调整卫星模型的卫星姿态。但是,在使用欧拉角需要转换左右手系时,对于卫星模型的不同轴向来说转换公式不一样,则利用欧拉角调整卫星姿态的方式较为复杂繁琐,且容易出现万向节锁死的问题,无法轻易解决;四元数表示四维向量,对于卫星模型的三维向量而言,利用四元数调整卫星姿态不够直观,不利于想象和观察。
为了解决这一问题,在本申请实施例中,依据欧拉角和J2000坐标系下X轴的朝向向量计算卫星模型的头朝向目标向量;欧拉角包括卫星模型的滚转角、俯仰角和偏航角;按照卫星模型的头朝向当前向量和头朝向目标向量调整卫星模型的头朝向;依据欧拉角和J2000坐标系下Y轴的朝向向量计算卫星模型的帆板朝向目标向量;按照卫星模型的帆板朝向当前向量和帆板朝向目标向量调整卫星模型的帆板朝向。可见,先计算卫星模型的头朝向目标向量,在卫星模型的头朝向当前向量的基础上,采用向量法调整卫星模型的头朝向,再计算卫星模型的帆板朝向目标向量,在卫星模型的帆板朝向当前向量的基础上,采用向量法调整卫星模型的帆板朝向,向量法能够屏蔽基于欧拉角和/或四元数调整卫星姿态的缺点;先调整卫星模型的头朝向再调整卫星模型的帆板朝向,不仅调整方式直观形象,而且有效避免了万向节锁死的情况。
举例来说,本申请实施例的场景之一,可以是应用到如图1所示的场景中,该场景包括卫星模型101和处理器102,卫星模型101与处理器102进行交互,处理器102获取欧拉角后,采用本申请实施例的实施方式调整卫星模型101的头朝向和帆板朝向,以实现卫星模型101的卫星姿态调整。
可以理解的是,在上述应用场景中,虽然将本申请实施方式的动作描述由处理器102执行,但是,本申请在执行主体方面不受限制,只要执行了本申请实施方式所公开的动作即可。
可以理解的是,上述场景仅是本申请实施例提供的一个场景示例,本申请实施例并不限于此场景。
下面结合附图,通过实施例来详细说明本申请实施例中调整卫星姿态的方法和相关装置的具体实现方式。
示例性方法
参见图2,示出了本申请实施例中一种调整卫星姿态的方法的流程示意图。在本实施例中,所述方法例如可以包括以下步骤:
步骤201:基于欧拉角和J2000坐标系下X轴的朝向向量,获得卫星模型的头朝向目标向量;所述欧拉角包括所述卫星模型的滚转角、俯仰角和偏航角。
由于现有技术中利用欧拉角和/或四元数调整卫星模型的卫星姿态,在使用欧拉角需要转换左右手系时,对于卫星模型的不同轴向来说转换公式不一样,则利用欧拉角调整卫星姿态的方式较为复杂繁琐,且容易出现万向节锁死的问题,无法轻易解决;四元数表示四维向量,对于卫星模型的三维向量而言,利用四元数调整卫星姿态不够直观,不利于想象和观察。而向量法在综合利用欧拉角和四元数调整卫星姿态的优点的基础上,能够屏蔽利用欧拉角和四元数调整卫星姿态的缺点,且针对卫星模型的三个轴向的特性,在对其中两个轴向进行调整后,剩余一个轴向不需要进行调整,使得调整卫星模型的卫星姿态的误差较小;因此,在本申请实施例中,考虑采用向量法调整卫星模型的卫星姿态。采用向量法调整卫星模型的卫星姿态是指采用向量法调整卫星模型的头朝向和帆板朝向,在调整顺序上需要先采用向量法调整卫星模型的头朝向,再采用向量法调整卫星模型的帆板朝向。
本申请实施例中,采用向量法调整卫星模型的头朝向是指在已知卫星模型的头朝向当前向量的基础上,先确定卫星模型的头朝向目标向量,再进行卫星模型的头朝向调整。
由于卫星模型处于J2000坐标系中,卫星模型的头朝向与J2000坐标系下X轴的朝向相同,因此,依据包括卫星模型的滚转角、俯仰角和偏航角的欧拉角,以及J2000坐标系下X轴的朝向向量,即可计算得到卫星模型的头朝向目标向量。具体地,卫星模型的头朝向目标向量实际上是指将J2000坐标系下X轴的朝向向量绕欧拉角旋转所得到的向量。即,在本申请实施例一种可选的实施方式中,所述步骤201例如具体可以为:将所述J2000坐标系下X轴的朝向向量绕所述欧拉角旋转后的向量作为所述卫星模型的头朝向目标向量。
步骤202:基于所述卫星模型的头朝向当前向量和所述头朝向目标向量,调整所述卫星模型的头朝向。
在步骤201获得卫星模型的头朝向目标向量之后,依据已知的卫星模型的头朝向当前向量和卫星模型的头朝向目标向量,将卫星模型的头朝向当前向量转换到卫星模型的头朝向目标向量,即可实现卫星模型的头朝向调整。
具体地,针对卫星模型的头朝向当前向量和卫星模型的头朝向目标向量,可以计算得到头朝向当前向量与头朝向目标向量之间夹角角度作为第一夹角角度;由于卫星模型的头朝向一直朝向地球的球心,第一夹角角度的正负不会引起卫星模型的头朝向乱转问题,因此按照该第一夹角角度直接将卫星模型的头朝向当前向量转换到卫星模型的头朝向目标向量,即可实现卫星模型的头朝向调整。即,在本申请实施例一种可选的实施方式中,所述步骤202例如可以包括以下步骤:
步骤A:基于所述头朝向当前向量和所述头朝向目标向量,获得所述头朝向当前向量与所述头朝向目标向量之间的第一夹角角度;
步骤B:根据所述第一夹角角度直接调整所述卫星模型的头朝向。
步骤203:基于所述欧拉角和所述J2000坐标系下Y轴的朝向向量,获得所述卫星模型的帆板朝向目标向量。
本申请实施例中,同上述说明可知,采用向量法调整卫星模型的帆板朝向是指在已知卫星模型的帆板朝向当前向量的基础上,先确定卫星模型的帆板朝向目标向量,再进行卫星模型的帆板朝向调整。
由于卫星模型处于J2000坐标系中,卫星模型的帆板朝向与J2000坐标系下Y轴的朝向相同,因此,依据包括卫星模型的滚转角、俯仰角和偏航角的欧拉角,以及J2000坐标系下Y轴的朝向向量,即可计算得到卫星模型的帆板朝向目标向量。具体地,卫星模型的帆板朝向目标向量实际上是指将J2000坐标系下Y轴的朝向向量绕欧拉角旋转所得到的向量。即,在本申请实施例一种可选的实施方式中,所述步骤203例如具体可以为:将所述J2000坐标系下Y轴的朝向向量绕所述欧拉角旋转后的向量作为所述卫星模型的帆板朝向目标向量。
步骤204:基于所述卫星模型的帆板朝向当前向量和所述帆板朝向目标向量,调整所述卫星模型的帆板朝向。
在步骤203获得卫星模型的帆板朝向目标向量之后,依据已知的卫星模型的帆板朝向当前向量和卫星模型的帆板朝向目标向量,将卫星模型的帆板朝向当前向量转换到卫星模型的帆板朝向目标向量,即可实现卫星模型的帆板朝向调整。
具体地,针对卫星模型的帆板朝向当前向量和卫星模型的帆板朝向目标向量,可以计算得到帆板朝向当前向量与帆板朝向目标向量之间夹角角度作为第二夹角角度;卫星模型可能处于赤道上方,也可能处于赤道下方,第二夹角角度的正负容易引起卫星模型的帆板朝向乱转问题,为了解决该问题,需要基于表示卫星模型处于赤道上方还是赤道下方的卫星模型的Y轴位置值结合第二夹角角度,将卫星模型的帆板朝向当前向量转换到卫星模型的帆板朝向目标向量,才能实现卫星模型的帆板朝向调整。即,在本申请实施例一种可选的实施方式中,所述步骤204例如可以包括以下步骤包括:
步骤C:基于所述帆板朝向当前向量和所述帆板朝向目标向量,获得所述帆板朝向当前向量与所述帆板朝向目标向量之间的第二夹角角度;
步骤D:基于所述卫星模型的Y轴位置值和所述第二夹角角度,调整所述卫星模型的帆板朝向。
当卫星模型的Y轴位置值大于0时,表示卫星模型处于赤道上方,此时需要利用180°减去第二夹角角度,得到180°与第二夹角角度的差值角度,按照该差值角度将卫星模型的帆板朝向当前向量转换到卫星模型的帆板朝向目标向量,才能实现卫星模型的帆板朝向调整。因此,在本申请实施例一种可选的实施方式中,所述步骤D例如具体可以为:若所述卫星模型的Y轴位置值大于0,根据180°与所述第二夹角角度的差值角度调整所述卫星模型的帆板朝向。
当卫星模型的Y轴位置值小于等于0时,表示卫星模型处于赤道上方,此时需要利用180°加上第二夹角角度,得到180°与第二夹角角度的加和角度,按照该加和角度将卫星模型的帆板朝向当前向量转换到卫星模型的帆板朝向目标向量,才能实现卫星模型的帆板朝向调整。因此,在本申请实施例一种可选的实施方式中,所述步骤D例如具体可以为:若所述卫星模型的Y轴位置值小于等于0,根据180°与所述第二夹角角度的加和角度调整所述卫星模型的帆板朝向。
作为一种示例,如图3所示的卫星姿态调整前卫星模型示意图;该图中卫星模型处于J2000坐标系中,卫星模型的头朝向与J2000坐标系下X轴的朝向相同,卫星模型的帆板朝向与J2000坐标系下Y轴的朝向相同,在已知欧拉角、卫星模型的头朝向当前向量和卫星模型的帆板朝向当前向量的基础上,执行本申请实施例的步骤201-步骤204,先调整卫星模型的头朝向,再调整卫星模型的帆板朝向,无需调整卫星模型的其他朝向,即可实现卫星模型的卫星姿态调整,且使得调整卫星模型的卫星姿态的误差较小。
通过本实施例提供的各种实施方式,首先,依据欧拉角和J2000坐标系下X轴的朝向向量计算卫星模型的头朝向目标向量;欧拉角包括卫星模型的滚转角、俯仰角和偏航角;其次,按照卫星模型的头朝向当前向量和头朝向目标向量调整卫星模型的头朝向;然后,依据欧拉角和J2000坐标系下Y轴的朝向向量计算卫星模型的帆板朝向目标向量;最后,按照卫星模型的帆板朝向当前向量和帆板朝向目标向量调整卫星模型的帆板朝向。由此可见,先计算卫星模型的头朝向目标向量,在卫星模型的头朝向当前向量的基础上,采用向量法调整卫星模型的头朝向,再计算卫星模型的帆板朝向目标向量,在卫星模型的帆板朝向当前向量的基础上,采用向量法调整卫星模型的帆板朝向,向量法能够屏蔽基于欧拉角和/或四元数调整卫星姿态的缺点;先调整卫星模型的头朝向再调整卫星模型的帆板朝向,不仅调整方式直观形象,而且有效避免了万向节锁死的情况。
示例性装置
参见图4,示出了本申请实施例中一种调整卫星姿态的装置的结构示意图。在本实施例中,所述装置例如具体可以包括:
第一获得单元401,用于基于欧拉角和J2000坐标系下X轴的朝向向量,获得卫星模型的头朝向目标向量;所述欧拉角包括所述卫星模型的滚转角、俯仰角和偏航角;
第一调整单元402,用于基于所述卫星模型的头朝向当前向量和所述头朝向目标向量,调整所述卫星模型的头朝向;
第二获得单元403,用于基于所述欧拉角和所述J2000坐标系下Y轴的朝向向量,获得所述卫星模型的帆板朝向目标向量;
第二调整单元404,用于基于所述卫星模型的帆板朝向当前向量和所述帆板朝向目标向量,调整所述卫星模型的帆板朝向。
在本申请实施例一种可选的实施方式中,所述第一获得单元401具体用于:
将所述J2000坐标系下X轴的朝向向量绕所述欧拉角旋转后的向量作为所述卫星模型的头朝向目标向量。
在本申请实施例一种可选的实施方式中,所述第一调整单元402包括:
第一获得子单元,用于基于所述头朝向当前向量和所述头朝向目标向量,获得所述头朝向当前向量与所述头朝向目标向量之间的第一夹角角度;
第一调整子单元,用于根据所述第一夹角角度直接调整所述卫星模型的头朝向。
在本申请实施例一种可选的实施方式中,所述第二获得单元403具体用于:
将所述J2000坐标系下Y轴的朝向向量绕所述欧拉角旋转后的向量作为所述卫星模型的帆板朝向目标向量。
在本申请实施例一种可选的实施方式中,所述第二调整单元404包括:
第二获得子单元,用于基于所述帆板朝向当前向量和所述帆板朝向目标向量,获得所述帆板朝向当前向量与所述帆板朝向目标向量之间的第二夹角角度;
第二调整子单元,用于基于所述卫星模型的Y轴位置值和所述第二夹角角度,调整所述卫星模型的帆板朝向。
在本申请实施例一种可选的实施方式中,所述第二调整子单元具体用于:
若所述卫星模型的Y轴位置值大于0,根据180°与所述第二夹角角度的差值角度调整所述卫星模型的帆板朝向。
在本申请实施例一种可选的实施方式中,所述第二调整子单元具体用于:
若所述卫星模型的Y轴位置值小于等于0,根据180°与所述第二夹角角度的加和角度调整所述卫星模型的帆板朝向。
通过本实施例提供的各种实施方式,首先,依据欧拉角和J2000坐标系下X轴的朝向向量计算卫星模型的头朝向目标向量;欧拉角包括卫星模型的滚转角、俯仰角和偏航角;其次,按照卫星模型的头朝向当前向量和头朝向目标向量调整卫星模型的头朝向;然后,依据欧拉角和J2000坐标系下Y轴的朝向向量计算卫星模型的帆板朝向目标向量;最后,按照卫星模型的帆板朝向当前向量和帆板朝向目标向量调整卫星模型的帆板朝向。由此可见,先计算卫星模型的头朝向目标向量,在卫星模型的头朝向当前向量的基础上,采用向量法调整卫星模型的头朝向,再计算卫星模型的帆板朝向目标向量,在卫星模型的帆板朝向当前向量的基础上,采用向量法调整卫星模型的帆板朝向,向量法能够屏蔽基于欧拉角和/或四元数调整卫星姿态的缺点;先调整卫星模型的头朝向再调整卫星模型的帆板朝向,不仅调整方式直观形象,而且有效避免了万向节锁死的情况。
此外,本申请实施例还提供了一种终端设备,所述终端设备包括处理器以及存储器:
所述存储器用于存储程序代码,并将所述程序代码传输给所述处理器;
所述处理器用于根据所述程序代码中的指令执行上述方法实施例所述的调整卫星姿态的方法。
本申请实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质用于存储程序代码,所述程序代码用于执行上述方法实施例所述的调整卫星姿态的方法。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
专业人员还可以进一步意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本申请的范围。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述,仅是本申请的较佳实施例而已,并非对本申请作任何形式上的限制。虽然本申请已以较佳实施例揭露如上,然而并非用以限定本申请。任何熟悉本领域的技术人员,在不脱离本申请技术方案范围情况下,都可利用上述揭示的方法和技术内容对本申请技术方案做出许多可能的变动和修饰,或修改为等同变化的等效实施例。因此,凡是未脱离本申请技术方案的内容,依据本申请的技术实质对以上实施例所做的任何简单修改、等同变化及修饰,均仍属于本申请技术方案保护的范围内。
Claims (10)
1.一种调整卫星姿态的方法,其特征在于,包括:
基于欧拉角和J2000坐标系下X轴的朝向向量,获得卫星模型的头朝向目标向量;所述欧拉角包括所述卫星模型的滚转角、俯仰角和偏航角;
基于所述卫星模型的头朝向当前向量和所述头朝向目标向量,调整所述卫星模型的头朝向;
基于所述欧拉角和所述J2000坐标系下Y轴的朝向向量,获得所述卫星模型的帆板朝向目标向量;
基于所述卫星模型的帆板朝向当前向量和所述帆板朝向目标向量,调整所述卫星模型的帆板朝向。
2.根据权利要求1所述的方法,其特征在于,所述基于欧拉角和J2000坐标系下X轴的朝向向量,获得卫星模型的头朝向目标向量,具体为:
将所述J2000坐标系下X轴的朝向向量绕所述欧拉角旋转后的向量作为所述卫星模型的头朝向目标向量。
3.根据权利要求1所述的方法,其特征在于,所述基于所述卫星模型的头朝向当前向量和所述头朝向目标向量,调整所述卫星模型的头朝向,包括:
基于所述头朝向当前向量和所述头朝向目标向量,获得所述头朝向当前向量与所述头朝向目标向量之间的第一夹角角度;
根据所述第一夹角角度直接调整所述卫星模型的头朝向。
4.根据权利要求1所述的方法,其特征在于,所述基于所述欧拉角和所述J2000坐标系下Y轴的朝向向量,获得所述卫星模型的帆板朝向目标向量,具体为:
将所述J2000坐标系下Y轴的朝向向量绕所述欧拉角旋转后的向量作为所述卫星模型的帆板朝向目标向量。
5.根据权利要求1所述的方法,其特征在于,所述基于所述卫星模型的帆板朝向当前向量和所述帆板朝向目标向量,调整所述卫星模型的帆板朝向,包括:
基于所述帆板朝向当前向量和所述帆板朝向目标向量,获得所述帆板朝向当前向量与所述帆板朝向目标向量之间的第二夹角角度;
基于所述卫星模型的Y轴位置值和所述第二夹角角度,调整所述卫星模型的帆板朝向。
6.根据权利要求5所述的方法,其特征在于,所述基于所述卫星模型的Y轴位置值和所述第二夹角角度,调整所述卫星模型的帆板朝向,具体为:
若所述卫星模型的Y轴位置值大于0,根据180°与所述第二夹角角度的差值角度调整所述卫星模型的帆板朝向。
7.根据权利要求6所述的方法,其特征在于,所述基于所述卫星模型的Y轴位置值和所述第二夹角角度,调整所述卫星模型的帆板朝向,具体为:
若所述卫星模型的Y轴位置值小于等于0,根据180°与所述第二夹角角度的加和角度调整所述卫星模型的帆板朝向。
8.一种调整卫星姿态的装置,其特征在于,包括:
第一获得单元,用于基于欧拉角和J2000坐标系下X轴的朝向向量,获得卫星模型的头朝向目标向量;所述欧拉角包括所述卫星模型的滚转角、俯仰角和偏航角;
第一调整单元,用于基于所述卫星模型的头朝向当前向量和所述头朝向目标向量,调整所述卫星模型的头朝向;
第二获得单元,用于基于所述欧拉角和所述J2000坐标系下Y轴的朝向向量,获得所述卫星模型的帆板朝向目标向量;
第二调整单元,用于基于所述卫星模型的帆板朝向当前向量和所述帆板朝向目标向量,调整所述卫星模型的帆板朝向。
9.一种终端设备,其特征在于,所述终端设备包括处理器以及存储器:
所述存储器用于存储程序代码,并将所述程序代码传输给所述处理器;
所述处理器用于根据所述程序代码中的指令执行权利要求1-7任一项所述的调整卫星姿态的方法。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质用于存储程序代码,所述程序代码用于执行权利要求1-7任一项所述的调整卫星姿态的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010413058.4A CN111591472B (zh) | 2020-05-15 | 2020-05-15 | 一种调整卫星姿态的方法和相关装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010413058.4A CN111591472B (zh) | 2020-05-15 | 2020-05-15 | 一种调整卫星姿态的方法和相关装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111591472A true CN111591472A (zh) | 2020-08-28 |
CN111591472B CN111591472B (zh) | 2021-12-10 |
Family
ID=72191065
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010413058.4A Active CN111591472B (zh) | 2020-05-15 | 2020-05-15 | 一种调整卫星姿态的方法和相关装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111591472B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112255599A (zh) * | 2020-10-16 | 2021-01-22 | 北京世冠金洋科技发展有限公司 | 实现卫星雷达波范围在地图上展示的方法、装置及设备 |
Citations (37)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH02262500A (ja) * | 1988-12-09 | 1990-10-25 | Hughes Aircraft Co | 衛星制御システム |
US6135392A (en) * | 1998-09-29 | 2000-10-24 | Hughes Electronics Corporation | Spacecraft attitude control actuator and method |
US20020116993A1 (en) * | 1998-11-06 | 2002-08-29 | Christian Patouraux | Method and apparatus for determining the yaw angle of a satellite |
US6470270B1 (en) * | 2000-11-07 | 2002-10-22 | The Boeing Company | Lost-in-space stellar attitude acquisition using multiple star trackers |
US20040140401A1 (en) * | 2002-08-30 | 2004-07-22 | Nec Corporation | System and method for controlling the attitude of a flying object |
US20040144910A1 (en) * | 2003-01-28 | 2004-07-29 | Peck Mason A | Attitude determination system and method |
US20080315039A1 (en) * | 2007-06-21 | 2008-12-25 | Lael Rudd | System and methods for space vehicle torque balancing |
CN101402398A (zh) * | 2008-11-18 | 2009-04-08 | 航天东方红卫星有限公司 | 一种卫星姿态快速挽救方法 |
CN101825475A (zh) * | 2010-05-17 | 2010-09-08 | 哈尔滨工业大学 | 空间光学遥感传感器的像移补偿方法 |
US20100235132A1 (en) * | 2002-10-04 | 2010-09-16 | Havens Steven W | Method and apparatus for acquiring and processing transducer data |
CN101858969A (zh) * | 2010-03-26 | 2010-10-13 | 航天东方红卫星有限公司 | 一种基于优化虚拟转动的卫星目标姿态预先确定方法 |
CN101858747A (zh) * | 2010-03-26 | 2010-10-13 | 航天东方红卫星有限公司 | 一种有效利用地球辐照能的卫星帆板对日定向目标姿态的解析确定方法 |
CN102099250A (zh) * | 2008-06-09 | 2011-06-15 | 阿斯特里姆简易股份公司 | 用于控制卫星姿态的方法和姿态受控卫星 |
CN102346033A (zh) * | 2010-08-06 | 2012-02-08 | 清华大学 | 基于卫星观测角误差估计的直接定位方法及系统 |
CN102944241A (zh) * | 2012-11-15 | 2013-02-27 | 北京理工大学 | 基于多胞型线性微分包含的航天器相对姿态确定方法 |
US20130092785A1 (en) * | 2008-07-11 | 2013-04-18 | Davidson Technologies, Inc. | System and method for guiding and controlling a missile using high order sliding mode control |
CN103279127A (zh) * | 2013-05-22 | 2013-09-04 | 上海新跃仪表厂 | 一种仅用角度信息的geo轨道卫星自主控制方法 |
CN103344872A (zh) * | 2013-06-28 | 2013-10-09 | 上海新跃仪表厂 | 一种星敏安装极性的测试方法 |
CN103941593A (zh) * | 2014-04-04 | 2014-07-23 | 国家测绘地理信息局卫星测绘应用中心 | 低轨卫星姿态仿真方法 |
US20140361124A1 (en) * | 2013-06-07 | 2014-12-11 | Thales | Propulsion system with four modules for satellite orbit control and attitude control |
CN104332707A (zh) * | 2014-10-27 | 2015-02-04 | 西安空间无线电技术研究所 | 一种用于低轨星载天线跟踪地面站的方法 |
CN104570742A (zh) * | 2015-01-29 | 2015-04-29 | 哈尔滨工业大学 | 基于前馈pid控制的异面交叉快变轨道快速高精度相对指向控制方法 |
CN104848860A (zh) * | 2015-05-19 | 2015-08-19 | 北京空间飞行器总体设计部 | 一种敏捷卫星成像过程姿态机动规划方法 |
CN104881036A (zh) * | 2015-06-19 | 2015-09-02 | 哈尔滨工业大学 | 基于代数Lyapunov方程的控制受限小卫星三轴磁力矩姿态控制方法 |
CN105035371A (zh) * | 2015-08-04 | 2015-11-11 | 北京控制工程研究所 | 一种基于osg三维引擎的经典轨道三维空间关系构建方法 |
CN105139406A (zh) * | 2015-09-08 | 2015-12-09 | 哈尔滨工业大学 | 一种基于序列图像的跟踪精度反演方法 |
CN105160125A (zh) * | 2015-09-24 | 2015-12-16 | 航天东方红卫星有限公司 | 一种星敏感器四元数的仿真分析方法 |
CN105184002A (zh) * | 2015-09-17 | 2015-12-23 | 航天东方红卫星有限公司 | 一种数传天线指向角度的仿真分析方法 |
CN105446346A (zh) * | 2015-11-26 | 2016-03-30 | 航天东方红卫星有限公司 | 遥感卫星对月相对定标姿态调整方法 |
CN105807781A (zh) * | 2014-12-31 | 2016-07-27 | 上海新跃仪表厂 | 一种基于比例导引的空间近距飞越末制导方法 |
CN106096148A (zh) * | 2016-06-14 | 2016-11-09 | 中国空间技术研究院 | 一种简单姿态控制下的大倾角轨道卫星太阳帆板指向方法 |
CN106197425A (zh) * | 2016-06-30 | 2016-12-07 | 中国电子科技集团公司第五十四研究所 | 基于卫星姿态角的地面目标点位置的计算方法 |
CN106292677A (zh) * | 2016-07-15 | 2017-01-04 | 上海航天控制技术研究所 | 基于恒星时角的姿态控制方法和系统 |
US9643740B2 (en) * | 2013-02-15 | 2017-05-09 | The Boeing Company | Gyroless three-axis sun acquisition using sun sensor and unscented kalman filter |
US20170269610A1 (en) * | 2016-03-17 | 2017-09-21 | Mitsubishi Electric Research Laboratories, Inc. | Concurrent Station Keeping, Attitude Control, and Momentum Management of Spacecraft |
CN107885546A (zh) * | 2017-11-07 | 2018-04-06 | 郑州师范学院 | 一种面向全空间信息系统的坐标系转换方法 |
US20190049998A1 (en) * | 2017-08-10 | 2019-02-14 | Mitsubishi Electric Research Laboratories, Inc. | Model predictive control of spacecraft |
-
2020
- 2020-05-15 CN CN202010413058.4A patent/CN111591472B/zh active Active
Patent Citations (37)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH02262500A (ja) * | 1988-12-09 | 1990-10-25 | Hughes Aircraft Co | 衛星制御システム |
US6135392A (en) * | 1998-09-29 | 2000-10-24 | Hughes Electronics Corporation | Spacecraft attitude control actuator and method |
US20020116993A1 (en) * | 1998-11-06 | 2002-08-29 | Christian Patouraux | Method and apparatus for determining the yaw angle of a satellite |
US6470270B1 (en) * | 2000-11-07 | 2002-10-22 | The Boeing Company | Lost-in-space stellar attitude acquisition using multiple star trackers |
US20040140401A1 (en) * | 2002-08-30 | 2004-07-22 | Nec Corporation | System and method for controlling the attitude of a flying object |
US20100235132A1 (en) * | 2002-10-04 | 2010-09-16 | Havens Steven W | Method and apparatus for acquiring and processing transducer data |
US20040144910A1 (en) * | 2003-01-28 | 2004-07-29 | Peck Mason A | Attitude determination system and method |
US20080315039A1 (en) * | 2007-06-21 | 2008-12-25 | Lael Rudd | System and methods for space vehicle torque balancing |
CN102099250A (zh) * | 2008-06-09 | 2011-06-15 | 阿斯特里姆简易股份公司 | 用于控制卫星姿态的方法和姿态受控卫星 |
US20130092785A1 (en) * | 2008-07-11 | 2013-04-18 | Davidson Technologies, Inc. | System and method for guiding and controlling a missile using high order sliding mode control |
CN101402398A (zh) * | 2008-11-18 | 2009-04-08 | 航天东方红卫星有限公司 | 一种卫星姿态快速挽救方法 |
CN101858969A (zh) * | 2010-03-26 | 2010-10-13 | 航天东方红卫星有限公司 | 一种基于优化虚拟转动的卫星目标姿态预先确定方法 |
CN101858747A (zh) * | 2010-03-26 | 2010-10-13 | 航天东方红卫星有限公司 | 一种有效利用地球辐照能的卫星帆板对日定向目标姿态的解析确定方法 |
CN101825475A (zh) * | 2010-05-17 | 2010-09-08 | 哈尔滨工业大学 | 空间光学遥感传感器的像移补偿方法 |
CN102346033A (zh) * | 2010-08-06 | 2012-02-08 | 清华大学 | 基于卫星观测角误差估计的直接定位方法及系统 |
CN102944241A (zh) * | 2012-11-15 | 2013-02-27 | 北京理工大学 | 基于多胞型线性微分包含的航天器相对姿态确定方法 |
US9643740B2 (en) * | 2013-02-15 | 2017-05-09 | The Boeing Company | Gyroless three-axis sun acquisition using sun sensor and unscented kalman filter |
CN103279127A (zh) * | 2013-05-22 | 2013-09-04 | 上海新跃仪表厂 | 一种仅用角度信息的geo轨道卫星自主控制方法 |
US20140361124A1 (en) * | 2013-06-07 | 2014-12-11 | Thales | Propulsion system with four modules for satellite orbit control and attitude control |
CN103344872A (zh) * | 2013-06-28 | 2013-10-09 | 上海新跃仪表厂 | 一种星敏安装极性的测试方法 |
CN103941593A (zh) * | 2014-04-04 | 2014-07-23 | 国家测绘地理信息局卫星测绘应用中心 | 低轨卫星姿态仿真方法 |
CN104332707A (zh) * | 2014-10-27 | 2015-02-04 | 西安空间无线电技术研究所 | 一种用于低轨星载天线跟踪地面站的方法 |
CN105807781A (zh) * | 2014-12-31 | 2016-07-27 | 上海新跃仪表厂 | 一种基于比例导引的空间近距飞越末制导方法 |
CN104570742A (zh) * | 2015-01-29 | 2015-04-29 | 哈尔滨工业大学 | 基于前馈pid控制的异面交叉快变轨道快速高精度相对指向控制方法 |
CN104848860A (zh) * | 2015-05-19 | 2015-08-19 | 北京空间飞行器总体设计部 | 一种敏捷卫星成像过程姿态机动规划方法 |
CN104881036A (zh) * | 2015-06-19 | 2015-09-02 | 哈尔滨工业大学 | 基于代数Lyapunov方程的控制受限小卫星三轴磁力矩姿态控制方法 |
CN105035371A (zh) * | 2015-08-04 | 2015-11-11 | 北京控制工程研究所 | 一种基于osg三维引擎的经典轨道三维空间关系构建方法 |
CN105139406A (zh) * | 2015-09-08 | 2015-12-09 | 哈尔滨工业大学 | 一种基于序列图像的跟踪精度反演方法 |
CN105184002A (zh) * | 2015-09-17 | 2015-12-23 | 航天东方红卫星有限公司 | 一种数传天线指向角度的仿真分析方法 |
CN105160125A (zh) * | 2015-09-24 | 2015-12-16 | 航天东方红卫星有限公司 | 一种星敏感器四元数的仿真分析方法 |
CN105446346A (zh) * | 2015-11-26 | 2016-03-30 | 航天东方红卫星有限公司 | 遥感卫星对月相对定标姿态调整方法 |
US20170269610A1 (en) * | 2016-03-17 | 2017-09-21 | Mitsubishi Electric Research Laboratories, Inc. | Concurrent Station Keeping, Attitude Control, and Momentum Management of Spacecraft |
CN106096148A (zh) * | 2016-06-14 | 2016-11-09 | 中国空间技术研究院 | 一种简单姿态控制下的大倾角轨道卫星太阳帆板指向方法 |
CN106197425A (zh) * | 2016-06-30 | 2016-12-07 | 中国电子科技集团公司第五十四研究所 | 基于卫星姿态角的地面目标点位置的计算方法 |
CN106292677A (zh) * | 2016-07-15 | 2017-01-04 | 上海航天控制技术研究所 | 基于恒星时角的姿态控制方法和系统 |
US20190049998A1 (en) * | 2017-08-10 | 2019-02-14 | Mitsubishi Electric Research Laboratories, Inc. | Model predictive control of spacecraft |
CN107885546A (zh) * | 2017-11-07 | 2018-04-06 | 郑州师范学院 | 一种面向全空间信息系统的坐标系转换方法 |
Non-Patent Citations (2)
Title |
---|
FAN WU: "Dynamics and control of spacecraft with a large misaligned rotational component", 《AEROSPACE SCIENCE AND TECHNOLOGY》 * |
淡鹏: "三轴稳定卫星姿态可视化绘制与建模", 《电子科技》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112255599A (zh) * | 2020-10-16 | 2021-01-22 | 北京世冠金洋科技发展有限公司 | 实现卫星雷达波范围在地图上展示的方法、装置及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN111591472B (zh) | 2021-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hasan et al. | Analysis on Euler angles rotation of a rigid body in three-axis attitude based on RazakSAT data | |
US9073648B2 (en) | Star tracker rate estimation with kalman filter enhancement | |
CN107450582B (zh) | 一种基于星上实时规划的相控阵数传引导控制方法 | |
CN108427427B (zh) | 一种航天器对地表定向目标姿态角计算方法 | |
CN104459608B (zh) | 一种舰船纵横摇综合矢量坐标变换天线波束指向修正方法 | |
CN107329485B (zh) | 一种快速的多约束航天器姿态路径递归规划方法 | |
CN116009596A (zh) | 一种姿态角的计算方法、装置、介质、以及电子设备 | |
US11787569B2 (en) | System and method for optimizing a low-thrust trajectory of a spacecraft trajectory | |
CN111591472B (zh) | 一种调整卫星姿态的方法和相关装置 | |
CN115416877A (zh) | 一种基于小推力的卫星碰撞高度规避方法及智能化系统 | |
CN106569696B (zh) | 一种渲染输出全景影像的方法、系统及便携式终端 | |
CN109460049B (zh) | 基于惯性指向模式的地球同步轨道卫星远地点变轨方法 | |
CN116182781A (zh) | 一种星载三轴转台对地指向旋转角度的计算方法 | |
CN117022679A (zh) | 一种卫星的姿态控制方法、装置及卫星 | |
CN113815903B (zh) | 一种用于遥感卫星的飞轮过零规避方法 | |
CN114879709A (zh) | 一种面向运动目标跟踪观测的卫星姿态控制方法及装置 | |
CN114840812A (zh) | 欧拉角的计算方法、装置、设备及介质 | |
CN108225306A (zh) | 基于遥感卫星凝视姿态的星敏感器安装布局方法 | |
CN113968362B (zh) | 一种卫星在轨自主三轴快速机动控制方法 | |
CN105371851B (zh) | 一种基于频域分析的卫星姿态模型构建方法 | |
CN113830330A (zh) | 一种基于中继星测控的卫星姿态指向方法和系统 | |
CN112498745B (zh) | 一种低轨光学卫星立体成像自主规划方法 | |
CN112960145B (zh) | 一种遥感卫星对地姿态机动扫描的轨迹规划方法及系统 | |
CN113218249A (zh) | 跟随式遥操作战车及控制方法 | |
CN112379680A (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 |