CN111207776A - 一种适用于火星探测的星敏感器与陀螺联合标定方法 - Google Patents

一种适用于火星探测的星敏感器与陀螺联合标定方法 Download PDF

Info

Publication number
CN111207776A
CN111207776A CN202010116766.1A CN202010116766A CN111207776A CN 111207776 A CN111207776 A CN 111207776A CN 202010116766 A CN202010116766 A CN 202010116766A CN 111207776 A CN111207776 A CN 111207776A
Authority
CN
China
Prior art keywords
star sensor
attitude
gyroscope
quaternion
gyro
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
CN202010116766.1A
Other languages
English (en)
Other versions
CN111207776B (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.)
Shanghai Aerospace Control Technology Institute
Original Assignee
Shanghai Aerospace Control Technology Institute
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 Shanghai Aerospace Control Technology Institute filed Critical Shanghai Aerospace Control Technology Institute
Priority to CN202010116766.1A priority Critical patent/CN111207776B/zh
Publication of CN111207776A publication Critical patent/CN111207776A/zh
Application granted granted Critical
Publication of CN111207776B publication Critical patent/CN111207776B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
    • G01C25/005Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices

Landscapes

  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Gyroscopes (AREA)

Abstract

本发明公开了一种适用于火星探测的星敏感器与陀螺联合标定方法,使火星探测器在全生命周期任意姿态基准下,都能够完成星敏感器和陀螺的联合标定。火星探测器在稳态飞行过程中,根据星敏感器及陀螺的测量模型在惯性基准下构造标定算法。与现有技术相比,其有益效果是,一种适用于火星探测的星敏感器与陀螺联合标定方法,使火星探测器可以在惯性空间任意姿态下进行星敏感器和陀螺的联合标定,提高了标定的弧段,增加了星敏感器和陀螺联合标定的可靠性。

Description

一种适用于火星探测的星敏感器与陀螺联合标定方法
技术领域
本发明涉及一种适用于火星探测的星敏感器与陀螺联合标定方法,用于火星探测器在飞行全过程星敏感器与陀螺联合标定,适应火星探测器在任意惯性姿态下的星敏感器与陀螺联合标定需求。
背景技术
火星探测器在轨利用星敏感器直接确定探测器的惯性姿态,利用陀螺直接测量探测器相对惯性系的角速度。但星敏感器测量存在随机误差,陀螺存在常值漂移,因此,为保证姿态测量及控制精度,一方面需要对星敏感器测量的随机误差进行抑制,另一方面需要补偿陀螺的常值漂移。
近地卫星的工作模式一般为三轴稳定对地,其星敏感器与陀螺联合标定在对地轨道坐标系下进行,其实施一般要求卫星姿态相对于轨道坐标系为小姿态。而对于火星探测任务,姿态基准不固定,无法像近地卫星一样进行星敏感器和陀螺联合标定。
火星探测任务对地通讯、对日巡航以及特定姿态保持等任务需求要求探测器需具备在任意惯性姿态下进行星敏感器和陀螺联合标定的能力,因此,可以研究一种适用于火星探测的星敏感器与陀螺联合标定方法。
发明内容
本发明的目的在于:克服现有技术的不足,提出了一种一种适用于火星探测的星敏感器与陀螺联合标定方法,用于火星探测器在飞行全过程星敏感器与陀螺联合标定。
为了解决上述技术问题,本发明通过以下的技术方案实现:
一种适用于火星探测的星敏感器与陀螺联合标定方法,包括如下步骤:
(1)判断每拍星敏感器数据是否有效,如果星敏感器数据无效,当拍不再进行标定计算;如果星敏感器数据有效,则进入步骤(2);
(2)判断当前拍是否为首拍,如果是首拍,进入步骤(4);如果不是首拍,则进入步骤(3);
(3)按照星敏感器输出姿态四元数前后拍绝对值最大的分量符号决定是否对星敏感器当前拍的姿态四元数进行取反操作;
(4)进行陀螺漂移扣除计算;
(5)利用扣除常值漂移后的陀螺角速度,积分惯性姿态;
(6)计算星敏感器输出惯性姿态与陀螺积分惯性姿态之间的偏差;
(7)根据观测构型,预测星敏感器与陀螺积分之间的姿态偏差,进而计算陀螺常值漂移的估计值及姿态偏差的估计值;
(8)估计环绕器的惯性姿态。
进一步的,所述步骤(1)中,星敏感器数据有效的条件是星敏感器的模式处于星跟踪且数据有效标志为有效。
进一步的,所述步骤(3)中,按照星敏感器输出姿态四元数前后拍绝对值最大的分量符号决定是否对星敏感器当前拍的姿态四元数进行取反操作,具体为:判断当前星敏感器输出的姿态四元数绝对值最大的分量与前一拍的姿态四元数绝对值最大的分量符号是否相同,如果相同,则进入步骤(4);如果不同,将当前拍星敏感器输出的姿态四元数各分量取反后进入步骤(4)。
进一步的,所述步骤(4)进行陀螺漂移扣除计算,具体为:从惯性角速度中扣除陀螺漂移估计值,得到扣除常值漂移后的陀螺角速度。
进一步的,步骤(5)中,对扣除的陀螺漂移进行限幅,限幅范围为[-0.2,0.2]rad/s,陀螺积分惯性四元数采用按欧拉轴计算。
进一步的,所述步骤(6)中计算星敏感器输出惯性姿态与陀螺积分惯性姿态之间的偏差,计算方法为先对陀螺积分姿态求四元数的逆,然后按照四元数圈乘方法右乘星敏感器输出四元数,并取矢量部分作为姿态偏差的观测值。
进一步的,所述步骤(7)中,偏差姿态的估计值计算方法为上一拍姿态偏差值减去0.1倍的上一拍陀螺常值漂移估计值,偏差姿态的估计值初值为[0 0 0]。
进一步的,陀螺常值漂移计算所用的系数分别为K1=K2=K3=0.00361;K4=K5=K6=-0.0001305。
进一步的,所述步骤(9)中,环绕器惯性姿态的计算方法为将偏差姿态按照四元数构型归一化后,左边圈乘陀螺积分惯性姿态四元数。
本发明采用的方法与现有技术相比,其优点和有益效果是:
(1)本发明通过提出一种适用于火星探测的星敏感器与陀螺联合标定方法,解决了火星环绕器全生命周期下星敏感器与陀螺联合标定的问题。
(2)本发明标定方法可以使火星探测器在惯性空间任意姿态下进行星敏感器和陀螺的联合标定,提高了标定的弧段,增加了星敏感器和陀螺联合标定的可靠性。
(3)本发明对星敏感器输出姿态四元数的最大分量进行前后拍符号一致性判断,避免星敏感器输出姿态四元数标量过0时出现标定量跳变情况。
附图说明
图1为本发明的方法流程图。
具体实施方式
本发明公开了一种适用于火星探测的星敏感器与陀螺联合标定方法,使火星探测器在全生命周期任意姿态基准下,都能够完成星敏感器和陀螺的联合标定。火星探测器在稳态飞行过程中,根据星敏感器及陀螺的测量模型在惯性基准下构造标定算法。
如图1所示,本发明提出的一种适用于火星探测的星敏感器与陀螺联合标定方法,包括如下步骤:
(1)判断每拍星敏感器数据是否有效,如果星敏感器数据无效,当拍不再进行标定计算;如果星敏感器数据有效,则进入步骤(2);
星敏感器数据有效的条件是星敏感器的模式处于星跟踪且数据有效标志为有效。
(2)判断当前拍是否为首拍,如果是首拍,进入步骤(4);如果不是首拍,则进入步骤(3);
(3)按照星敏感器输出姿态四元数前后拍绝对值最大的分量符号决定是否对星敏感器当前拍的姿态四元数进行取反操作;
具体为:判断当前星敏感器输出的姿态四元数绝对值最大的分量与前一拍的姿态四元数绝对值最大的分量符号是否相同,如果相同,则进入步骤(4);如果不同,将当前拍星敏感器输出的姿态四元数各分量取反后进入步骤(4)。
(4)进行陀螺漂移扣除计算;
具体为:从惯性角速度中扣除陀螺漂移估计值,得到扣除常值漂移后的陀螺角速度。
(5)利用扣除常值漂移后的陀螺角速度,积分惯性姿态;
积分惯性姿态过程中,对扣除的陀螺漂移进行限幅,限幅范围为[-0.2,0.2]rad/s,陀螺积分惯性四元数采用按欧拉轴计算。
(6)计算星敏感器输出惯性姿态与陀螺积分惯性姿态之间的偏差;
计算方法为先对陀螺积分姿态求四元数的逆,然后按照四元数圈乘方法右乘星敏感器输出四元数,并取矢量部分作为姿态偏差的观测值。
(7)根据观测构型,预测星敏感器与陀螺积分之间的姿态偏差,进而计算陀螺常值漂移的估计值及姿态偏差的估计值;
偏差姿态的估计值计算方法为上一拍姿态偏差值减去0.1倍的上一拍陀螺常值漂移估计值,偏差姿态的估计值初值为[0 0 0]。陀螺常值漂移计算为已知计算方法,本发明仅需确定陀螺常值漂移计算所用的滤波增益系数分别为K1=K2=K3=0.00361;K4=K5=K6=-0.0001305。
(8)估计环绕器的惯性姿态。环绕器惯性姿态的计算方法为将偏差姿态按照四元数构型归一化后,左边圈乘陀螺积分惯性姿态四元数。
给出本发明实施例:
标定方法流程具体描述如下:
输入:惯性角速度ωxbi、ωybi、ωzbi,星敏计算的惯性系到本体系的姿态四元数qbi,陀螺惯性四元数qgbi(第一拍初值为qbi赋值),计算周期T;
输出:陀螺漂移估计值bx、by、bz
处理过程:
1.陀螺积分一步预测四元数
1)扣除陀螺漂移固定值
xbi=ωxbi-bx
ybi=ωybi-by
zbi=ωzbi-bz
2)调用陀螺积分惯性四元数算法,计算qgbi
注:qgbi的初值为算法启动时刻器上星敏有效值qbi
3)修正qgbi符号
qm=′max(qbi)
注:qm表示取qbi中绝对值最大的量。
如果
qm(n)·qm(n+1)<0
qgbi=-qgbi
2.测量偏差四元数求取
δqgi0=qgbi0·qbi0+qgbi1·qbi1+qgbi2·qbi2+qgbi3·qbi3
δqgi1=-qgbi1·qbi0+qgbi0·qbi1+qgbi3·qbi2-qgbi2·qbi3
δqgbi2=-qgbi2·qbi0-qgbi3·qbi1+qgbi0·qbi2+qgbi1·qbi3
δqgbi3=-qgbi3·qbi0+qgbi2·qbi1-qgbi1·qbi2+qgbi0·qbi3
3.陀螺漂移估计(初始值:δqi1=δqi2=δqi3=0,Δbx=Δby=Δbz=0,陀螺组合常值漂移的初值取bx=by=bz=0;(滤波增益:K1=K2=K3=0.00361;K4=K5=K6=-0.0001305)。
1)偏差四元数预测(仅取矢量部分)
δqif1=δqi1+0.5`ωzbi·δqi2-0.5`ωybiδqi3-0.25Δbx
δqif2=-0.5`ωzbi·δqi1+δqi2+0.5`ωxbiδqi3-0.25Δby
δqif3=0.5`ωybi·δqi1-0.5`ωxbi·δqi2+δqi3-0.25Δbz
2)偏差四元数及常值漂移增量估计
δqi1=δqif1+K1(δqgi1-δqif1)
δqi2=δqif2+K2(δqgi2-δqif2)
δqi3=δqif3+K3(δqgi3-δqif3)
Δbx=K4(δqgi1-δqif1)
Δby=K5(δqgi2-δqif2)
Δbz=K6(δqgi3-δqif3)
3)常值漂移估计
Figure BDA0002391730000000061
4)更新姿态
Figure BDA0002391730000000062
上述方法使火星探测器可以在惯性空间任意姿态下进行星敏感器和陀螺的联合标定,提高了标定的弧段,增加了星敏感器和陀螺联合标定的可靠性。

Claims (9)

1.一种适用于火星探测的星敏感器与陀螺联合标定方法,其特征在于包括如下步骤:
(1)判断每拍星敏感器数据是否有效,如果星敏感器数据无效,当拍不再进行标定计算;如果星敏感器数据有效,则进入步骤(2);
(2)判断当前拍是否为首拍,如果是首拍,进入步骤(4);如果不是首拍,则进入步骤(3);
(3)按照星敏感器输出姿态四元数前后拍绝对值最大的分量符号决定是否对星敏感器当前拍的姿态四元数进行取反操作;
(4)进行陀螺漂移扣除计算;
(5)利用扣除常值漂移后的陀螺角速度,积分惯性姿态;
(6)计算星敏感器输出惯性姿态与陀螺积分惯性姿态之间的偏差;
(7)根据观测构型,预测星敏感器与陀螺积分之间的姿态偏差,进而计算陀螺常值漂移的估计值及姿态偏差的估计值;
(8)估计环绕器的惯性姿态。
2.根据权利要求1所示的一种适用于火星探测的星敏感器与陀螺联合标定方法,其特征在于:所述步骤(1)中,星敏感器数据有效的条件是星敏感器的模式处于星跟踪且数据有效标志为有效。
3.根据权利要求1所述的一种适用于火星探测的星敏感器与陀螺联合标定方法,其特征在于:所述步骤(3)中,按照星敏感器输出姿态四元数前后拍绝对值最大的分量符号决定是否对星敏感器当前拍的姿态四元数进行取反操作,具体为:判断当前星敏感器输出的姿态四元数绝对值最大的分量与前一拍的姿态四元数绝对值最大的分量符号是否相同,如果相同,则进入步骤(4);如果不同,将当前拍星敏感器输出的姿态四元数各分量取反后进入步骤(4)。
4.根据权利要求1所述的一种适用于火星探测的星敏感器与陀螺联合标定方法,其特征在于:所述步骤(4)进行陀螺漂移扣除计算,具体为:从惯性角速度中扣除陀螺漂移估计值,得到扣除常值漂移后的陀螺角速度。
5.根据权利要求1所述的一种适用于火星探测的星敏感器与陀螺联合标定方法,其特征在于:步骤(5)中,对扣除的陀螺漂移进行限幅,限幅范围为[-0.2,0.2]rad/s,陀螺积分惯性四元数采用按欧拉轴计算。
6.根据权利要求1所述的一种适用于火星探测的星敏感器与陀螺联合标定方法,其特征在于:所述步骤(6)中计算星敏感器输出惯性姿态与陀螺积分惯性姿态之间的偏差,计算方法为先对陀螺积分姿态求四元数的逆,然后按照四元数圈乘方法右乘星敏感器输出四元数,并取矢量部分作为姿态偏差的观测值。
7.根据权利要求1所述的一种适用于火星探测的星敏感器与陀螺联合标定方法,其特征在于:所述步骤(7)中,偏差姿态的估计值计算方法为上一拍姿态偏差值减去0.1倍的上一拍陀螺常值漂移估计值,偏差姿态的估计值初值为[0 0 0]。
8.根据权利要求7所述的一种适用于火星探测的星敏感器与陀螺联合标定方法,其特征在于:陀螺常值漂移计算所用的系数分别为K1=K2=K3=0.00361;K4=K5=K6=-0.0001305。
9.根据权利要求1所述的一种适用于火星探测的星敏感器与陀螺联合标定方法,其特征在于:所述步骤(9)中,环绕器惯性姿态的计算方法为将偏差姿态按照四元数构型归一化后,左边圈乘陀螺积分惯性姿态四元数。
CN202010116766.1A 2020-02-25 2020-02-25 一种适用于火星探测的星敏感器与陀螺联合标定方法 Active CN111207776B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010116766.1A CN111207776B (zh) 2020-02-25 2020-02-25 一种适用于火星探测的星敏感器与陀螺联合标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010116766.1A CN111207776B (zh) 2020-02-25 2020-02-25 一种适用于火星探测的星敏感器与陀螺联合标定方法

Publications (2)

Publication Number Publication Date
CN111207776A true CN111207776A (zh) 2020-05-29
CN111207776B CN111207776B (zh) 2022-04-12

Family

ID=70784445

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010116766.1A Active CN111207776B (zh) 2020-02-25 2020-02-25 一种适用于火星探测的星敏感器与陀螺联合标定方法

Country Status (1)

Country Link
CN (1) CN111207776B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112665570A (zh) * 2020-11-30 2021-04-16 北京电子工程总体研究所 一种基于星敏感器的mems陀螺零偏在轨简化工程计算方法
CN113008272A (zh) * 2021-03-08 2021-06-22 航天科工空间工程发展有限公司 一种用于微小卫星的mems陀螺在轨常值漂移标定方法和系统
CN114413883A (zh) * 2021-12-23 2022-04-29 上海航天控制技术研究所 卫星姿态确定精度的提升方法、存储介质和电子设备

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080269988A1 (en) * 2003-03-20 2008-10-30 Feller Walter J Combined gnss gyroscope control system and method
CN101846510A (zh) * 2010-05-28 2010-09-29 北京航空航天大学 一种基于星敏感器和陀螺的高精度卫星姿态确定方法
CN103344872A (zh) * 2013-06-28 2013-10-09 上海新跃仪表厂 一种星敏安装极性的测试方法
CN103398713A (zh) * 2013-04-26 2013-11-20 哈尔滨工程大学 一种星敏感器/光纤惯性设备量测数据同步方法
CN104019817A (zh) * 2014-05-30 2014-09-03 哈尔滨工程大学 一种用于卫星姿态估计的范数约束强跟踪容积卡尔曼滤波方法
CN106767767A (zh) * 2016-11-23 2017-05-31 上海航天控制技术研究所 一种微纳多模星敏感器系统及其数据融合方法
CN107228674A (zh) * 2017-06-06 2017-10-03 上海航天控制技术研究所 一种针对星敏感器和陀螺联合滤波的改进方法
CN108827310A (zh) * 2018-07-12 2018-11-16 哈尔滨工程大学 一种船用星敏感器辅助陀螺仪在线标定方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080269988A1 (en) * 2003-03-20 2008-10-30 Feller Walter J Combined gnss gyroscope control system and method
CN101846510A (zh) * 2010-05-28 2010-09-29 北京航空航天大学 一种基于星敏感器和陀螺的高精度卫星姿态确定方法
CN103398713A (zh) * 2013-04-26 2013-11-20 哈尔滨工程大学 一种星敏感器/光纤惯性设备量测数据同步方法
CN103344872A (zh) * 2013-06-28 2013-10-09 上海新跃仪表厂 一种星敏安装极性的测试方法
CN104019817A (zh) * 2014-05-30 2014-09-03 哈尔滨工程大学 一种用于卫星姿态估计的范数约束强跟踪容积卡尔曼滤波方法
CN106767767A (zh) * 2016-11-23 2017-05-31 上海航天控制技术研究所 一种微纳多模星敏感器系统及其数据融合方法
CN107228674A (zh) * 2017-06-06 2017-10-03 上海航天控制技术研究所 一种针对星敏感器和陀螺联合滤波的改进方法
CN108827310A (zh) * 2018-07-12 2018-11-16 哈尔滨工程大学 一种船用星敏感器辅助陀螺仪在线标定方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112665570A (zh) * 2020-11-30 2021-04-16 北京电子工程总体研究所 一种基于星敏感器的mems陀螺零偏在轨简化工程计算方法
CN112665570B (zh) * 2020-11-30 2022-11-22 北京电子工程总体研究所 一种基于星敏感器的mems陀螺零偏在轨简化工程计算方法
CN113008272A (zh) * 2021-03-08 2021-06-22 航天科工空间工程发展有限公司 一种用于微小卫星的mems陀螺在轨常值漂移标定方法和系统
CN113008272B (zh) * 2021-03-08 2022-04-19 航天科工空间工程发展有限公司 一种用于微小卫星的mems陀螺在轨常值漂移标定方法和系统
CN114413883A (zh) * 2021-12-23 2022-04-29 上海航天控制技术研究所 卫星姿态确定精度的提升方法、存储介质和电子设备
CN114413883B (zh) * 2021-12-23 2023-09-05 上海航天控制技术研究所 卫星姿态确定精度的提升方法、存储介质和电子设备

Also Published As

Publication number Publication date
CN111207776B (zh) 2022-04-12

Similar Documents

Publication Publication Date Title
CN111207776B (zh) 一种适用于火星探测的星敏感器与陀螺联合标定方法
Huang et al. Online initialization and automatic camera-IMU extrinsic calibration for monocular visual-inertial SLAM
EP1585939B1 (en) Attitude change kalman filter measurement apparatus and method
KR101294623B1 (ko) 저궤도 영상촬영위성의 향상된 지상 정밀자세결정 방법
US8321076B2 (en) On-line inertia estimation for use in controlling an aerospace vehicle
US20100114517A1 (en) Method and system for orientation sensing
Sabatelli et al. A double stage Kalman filter for sensor fusion and orientation tracking in 9D IMU
JP3165474B2 (ja) 移動体の位置および配向を電磁的に測定するためのシステム
JP2012194175A (ja) 姿勢判定方法、位置算出方法及び姿勢判定装置
CN112798021B (zh) 基于激光多普勒测速仪的惯导系统行进间初始对准方法
JP2015148450A (ja) センサ誤差補正装置、imu較正システム、imu較正方法、および、imu較正プログラム
Chiang et al. Constrained filtering method for attitude determination using GPS and gyro
CN104819717B (zh) 一种基于mems惯性传感器组的多旋翼飞行器姿态检测方法
US20150134294A1 (en) Attitude calculation apparatus and attitude calculation method
JP2013122384A (ja) カルマンフィルタ、及び、状態推定装置
CN110375773B (zh) Mems惯导系统姿态初始化方法
JP2013061309A (ja) カルマンフィルタ、状態推定装置、カルマンフィルタの制御方法、及びカルマンフィルタの制御プログラム
CN105180946B (zh) 基于宽频测量的卫星高精度姿态确定方法及系统
RU2092402C1 (ru) Способ калибровки гироинерциальных измерителей бесплатформенной инерционной навигационной системы ориентации космического аппарата
KR100666160B1 (ko) 궤도기하학을 이용한 인공위성 자기센서의 바이어스보정방법
JP7016101B2 (ja) 情報処理システム、情報処理装置、プログラム、及び情報処理方法
JP2004045385A (ja) 移動体の姿勢検出装置
CN113049005A (zh) Gnss位置法辅助dvl误差标定方法及系统
US20210404811A1 (en) Method for solving attitude of rigid body based on function iterative integration
KR101505866B1 (ko) 자장 센서 바이어스 보정 장치 및 이를 이용한 바이어스 보정 방법

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