CN111207775A - 卫星gnc系统半物理仿真闭环测试陀螺零位标定方法 - Google Patents
卫星gnc系统半物理仿真闭环测试陀螺零位标定方法 Download PDFInfo
- Publication number
- CN111207775A CN111207775A CN202010112786.1A CN202010112786A CN111207775A CN 111207775 A CN111207775 A CN 111207775A CN 202010112786 A CN202010112786 A CN 202010112786A CN 111207775 A CN111207775 A CN 111207775A
- Authority
- CN
- China
- Prior art keywords
- satellite
- gyroscope
- gyro
- zero position
- semi
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
本发明涉及一种卫星GNC系统半物理仿真闭环测试陀螺零位标定方法。利用稳态半物理闭环试验中动力学角速度与星上遥测陀螺输出角速度进行实时比较,修正陀螺零位(含地速),以满足系统对高精度陀螺的性能要求。该方法不需要中断半物理仿真测试,可以提高陀螺零位标定的效率,精度。
Description
技术领域
本发明涉及某高轨卫星GNC系统半物理仿真试验中陀螺零位标定方法,本方法可以提高标定的效率和标定的精度,以满足GNC系统对高精度陀螺的使用要求,达到卫星GNC系统仿真测试的要求。
背景技术
某些卫星对陀螺精度要求很高,要利用陀螺输出的角速度进行长时间(4000s)积分,实时得到卫星姿态;陀螺在交付GNC系统时会对陀螺零位进行标定,标定值不包括地速,且很多卫星用陀螺零位开机零位重复性、零位长时间稳定性不高,卫星GNC系统半物理仿真试验有地速影响,在系统对陀螺零位标定精度较高时,GNC系统陀螺零位需要再次进行标定。GNC系统半物理仿真测试陀螺零位标定一般在试验开始前进行标定,待陀螺上电稳定两个小时后,再在无外界干扰环境下,采集陀螺敏感到的地速以及陀螺零位本身输出,持续半小时,对陀螺零位数据取平均值,标定陀螺零位。该方法效率比较低,精度差,一旦标定精度不符合,面临再次标定的困境。
发明内容
针对现有技术存在的不足,本发明要解决的技术问题是提供不中断卫星GNC系统半物理仿真闭环试验,快速、高精度标定陀螺零位的方法,提高GNC系统闭环仿真测试的效率。
本发明解决技术的方案是:卫星GNC系统半物理仿真闭环测试陀螺零位标定方法,通过下述方式实现:
设置卫星GNC系统半物理仿真闭环试验中卫星采用稳态飞轮PD控制算法闭环控制,将卫星动力学模型中干扰力矩设置为0Nm;
在卫星GNC系统半物理仿真闭环试验过程中,卫星GNC系统处于稳定状态后,在每个控制周期即控制节拍根据遥测下传的姿态角偏差,结合PD控制参数,计算出当前节拍陀螺零位,并进行陀螺零位标定。
优选的,利用下式计算当前节拍的陀螺零位:
其中,Tx(k),Ty(k),Tz(k)分别为当前节拍滚动、俯仰、偏航控制力矩,为0Nm;分别为当前节拍遥测下传的滚动、俯仰、偏航姿态角;KPx,KPy,KPz及Kdx,Kdy,Kdz为飞轮PD控制器参数,为常数;即为陀螺的零位。
优选的,卫星GNC系统处于稳定状态后,由飞轮PD控制器计算出当前节拍陀螺零位,并进行陀螺零位标定。
卫星GNC系统半物理仿真闭环测试陀螺零位标定方法,通过下述方式实现:
设置卫星GNC系统半物理仿真闭环试验中卫星采用稳态飞轮PID控制算法闭环控制;
在卫星GNC系统半物理仿真闭环试验过程中,卫星GNC系统处于稳定状态后,将陀螺实时输出角速度与动力学角速度进行比较,得到当前节拍陀螺零位,生成陀螺零位参数注数包,并发给星载计算机,对陀螺零位进行修正即注数修正,进而完成陀螺零位标定。
优选的,卫星GNC系统半物理仿真闭环试验过程中,卫星姿态确定方式先选择星敏+陀螺融合定姿方式;星敏提供初始姿态,然后选择单陀螺,仅由陀螺提供姿态信息进行姿态确定。
优选的,所述的注数修正目标至动力学三轴姿态角速度均值达到预期的0值位。
优选的,动力学三轴姿态角速度均值为不少于1分钟时长内计算的所有节拍动力学三轴姿态角速度的均值。
优选的,所述的0值位为不小于2×10-4°/s。
本发明与现有技术相比的有益效果是:
卫星用陀螺交付前一般会对其零位进行标定,标定值不包括地速,某高轨卫星用半液浮陀螺零位开机零位重复性、零位长时间稳定性不高,卫星GNC系统半物理仿真试验有地速影响,在系统对陀螺零位标定精度较高时,GNC系统陀螺零位需要再次进行标定。GNC系统半物理仿真测试陀螺零位标定一般在试验开始前进行标定,待陀螺上电稳定两个小时后,再在无外界干扰环境下,采集陀螺敏感到的地速以及陀螺零位本身输出,持续半小时,对陀螺零位数据取平均值,标定陀螺零位。该方法效率比较低,精度差,一旦标定精度不符合,面临再次标定的困境。本发明在卫星GNC闭环测试过程中利用稳态PD控制器参数标定陀螺零位(含地速),可以在卫星初态完成陀螺零位标定,提高卫星GNC系统半物理仿真测试效率。利用稳态半物理闭环试验中动力学角速度与星上遥测陀螺输出角速度进行实时比较,修正陀螺零位(含地速),可以更直接修正陀螺零位,以满足系统对高精度陀螺的性能要求。该方法不需要中断半物理仿真测试,可以提高陀螺零位标定的效率,精度。
本发明闭环仿真测试过程中快速标定陀螺零位,选择单陀螺进行卫星姿态测量,避免了星敏感器测量噪声的影响,充分利用陀螺短期内高精度、高稳定度的特性。
本发明将动力学角速度输出值作为陀螺零位标定效果的评判依据,动力学角速度精度高、更新快,陀螺零位修正残差可以快速获得,以便做进一步修正。
附图说明
图1为本发明卫星GNC闭环测试利用稳态PD控制器参数标定陀螺零位流程;
图2为本发明卫星GNC闭环测试利用动力学角速度比较标定陀螺零位流程;
图3为本发明卫星GNC系统半物理闭环仿真系统示意图。
具体实施方式
下面结合实施例对本发明作进一步阐述。
如图3所示,本发明卫星GNC系统半物理闭环仿真试验系统,采用PXI自动化实时仿真设备运行动力学模型,用于模拟卫星在轨运行,动力学将姿态信息发送给星上测量敏感器(星敏+陀螺)地测口或环境模拟设备,星上测量敏感器模拟在轨工作状态,产生姿态信息,星上计算机采集到姿态敏感(星敏+陀螺)信息后,发送控制指令给执行机构进行控制,动力学再采集到执行机构输出信息,形成闭环控制。
实施例1
本例给出一种卫星GNC系统半物理仿真闭环测试陀螺零位标定方法,如图1所示,利用稳态PD控制算法计算来标定陀螺零位:
卫星采用稳态飞轮PD控制算法闭环控制,将卫星动力学模型中干扰力矩设置为0Nm,利用PD控制器特性,此时卫星遥测姿态在陀螺常值漂移影响下,姿态角会有偏差,从而间接计算出陀螺零位值。具体介绍如下:
A)使GNC闭环测试按照流程进入稳态
卫星GNC系统采用轮控PD控制算法闭环仿真测试;卫星GNC系统处于稳定状态。
B)卫星动力学模型去掉干扰力矩影响
卫星动力学模型中设置干扰力矩为0Nm,使卫星不受干扰力矩影响,卫星又处于轮控稳态,无控制力矩输出。
C)利用GNC稳态轮控PD控制器计算陀螺零位
半物理仿真试验不需要中断,直接进行陀螺零位标定;PD控制器由于控制力矩为0,则根据遥测下传的姿态角偏差,结合控制器参数,可以计算出陀螺零位。利用下式进行计算。
Tx(k),Ty(k),Tz(k)分别为当前节拍滚动、俯仰、偏航控制力矩,因为卫星处于飞轮稳态控制,动力学设置干扰力矩为0Nm,故控制力矩为0Nm;分别为当前节拍滚动、俯仰、偏航姿态角,由闭环试验中遥测下传获取。KPx,KPy,KPz及Kdx,Kdy,Kdz为飞轮PD控制器参数,为常数。即为陀螺的零位。
该步骤卫星GNC处于稳定控制模式,用PD控制器计算陀螺零位,实现陀螺零位的标定。
实施例2
本例给出一种卫星GNC系统半物理仿真闭环测试陀螺零位标定方法,如图2所示,与动力学姿态角速度比较来标定陀螺零位:
卫星采用稳态飞轮PID控制算法闭环控制,将陀螺实时输出角速度与动力学角速度进行比较,修正陀螺零位。具体介绍如下:
A)用高精度星敏信息定姿,进行GNC系统闭环测试
卫星姿态确定方式先选择星敏+陀螺融合定姿方式;星敏提供初始姿态,然后选择单陀螺,仅由陀螺提供姿态信息进行姿态确定。
B)使GNC系统按照测试流程进入稳态闭环模式
卫星GNC系统采用轮控PID控制算法闭环仿真测试;卫星GNC控制系统处于稳定状态;经稳定控制后,卫星遥测姿态角小于0.009°,卫星遥测姿态角速度均值小于0.0002°/s,如果陀螺零位标定有误差,则动力学角速度不会到0°/s。
卫星姿态确定方式先选择星敏+陀螺融合定姿方式;星敏提供初始姿态,然后选择单陀螺,仅由陀螺提供姿态信息进行姿态确定。在卫星GNC飞轮稳定控制模式下,采用单陀螺进行姿态确定,使动力学三轴姿态角速度稳定在0°/s,此时卫星遥测角速度即为陀螺零位(包含地速影响)
C)进行多次实时迭代修正
半物理仿真试验不需要中断,直接进行陀螺零位标定;
动力学当前三轴姿态角速度ωx_DLX(k),ωy_DLX(k),ωz_DLX(k),由动力学模型计算得到;与星上遥测角速度ωx(k),ωy(k),ωz(k)进行比较,差值即为陀螺零位,对陀螺零位进行标定,注数修正,直到动力学三轴姿态角速度ωx_DLX(k),ωy_DLX(k),ωz_DLX(k)均值为0°/s。当然实际工程应用过程中不可能要求达到绝对的零值,本例中给出一种可参照的工程实现方式:即可以在不少于1分钟时长内计算的所有节拍动力学三轴姿态角速度的均值不小于2×10-4°/s。
由于星上计算机资源限制,遥测三轴角速度分辨率为0.0001°/s,动力学三轴角速度输出为double类型,远比星上遥测角速度分辨率高,对零位修正有较好的效果。
本发明未详细说明部分属于本领域技术人员的公知常识。
Claims (9)
1.卫星GNC系统半物理仿真闭环测试陀螺零位标定方法,其特征在于通过下述方式实现:
设置卫星GNC系统半物理仿真闭环试验中卫星采用稳态飞轮PD控制算法闭环控制,将卫星动力学模型中干扰力矩设置为0Nm;
在卫星GNC系统半物理仿真闭环试验过程中,卫星GNC系统处于稳定状态后,在每个控制周期即控制节拍根据遥测下传的姿态角偏差,结合PD控制参数,计算出当前节拍陀螺零位,并进行陀螺零位标定。
3.根据权利要求1或2所述的方法,其特征在于:卫星GNC系统处于稳定状态后,由飞轮PD控制器计算出当前节拍陀螺零位,并进行陀螺零位标定。
5.卫星GNC系统半物理仿真闭环测试陀螺零位标定方法,其特征在于通过下述方式实现:
设置卫星GNC系统半物理仿真闭环试验中卫星采用稳态飞轮PID控制算法闭环控制;
在卫星GNC系统半物理仿真闭环试验过程中,卫星GNC系统处于稳定状态后,将陀螺实时输出角速度与动力学角速度进行比较,得到当前节拍陀螺零位,生成陀螺零位参数注数包,并发给星载计算机,对陀螺零位进行修正即注数修正,进而完成陀螺零位标定。
6.根据权利要求5所述的方法,其特征在于:卫星GNC系统半物理仿真闭环试验过程中,卫星姿态确定方式先选择星敏+陀螺融合定姿方式;星敏提供初始姿态,然后选择单陀螺,仅由陀螺提供姿态信息进行姿态确定。
7.根据权利要求5所述的方法,其特征在于:所述的注数修正目标至动力学三轴姿态角速度均值达到预期的0值位。
8.根据权利要求7所述的方法,其特征在于:动力学三轴姿态角速度均值为不少于1分钟时长内计算的所有节拍动力学三轴姿态角速度的均值。
9.根据权利要求7所述的方法,其特征在于:所述的0值位为不小于2×10-4°/s。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010112786.1A CN111207775B (zh) | 2020-02-24 | 2020-02-24 | 卫星gnc系统半物理仿真闭环测试陀螺零位标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010112786.1A CN111207775B (zh) | 2020-02-24 | 2020-02-24 | 卫星gnc系统半物理仿真闭环测试陀螺零位标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111207775A true CN111207775A (zh) | 2020-05-29 |
CN111207775B CN111207775B (zh) | 2023-07-07 |
Family
ID=70789124
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010112786.1A Active CN111207775B (zh) | 2020-02-24 | 2020-02-24 | 卫星gnc系统半物理仿真闭环测试陀螺零位标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111207775B (zh) |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5205518A (en) * | 1991-11-25 | 1993-04-27 | General Electric Co. | Gyroless yaw control system for a three axis stabilized, zero-momentum spacecraft |
CN101214859A (zh) * | 2007-12-26 | 2008-07-09 | 北京控制工程研究所 | 一种变轨期间自主故障检测恢复控制的方法 |
CN101219713A (zh) * | 2007-12-26 | 2008-07-16 | 北京控制工程研究所 | 一种卫星的自主变轨方法 |
CN103017745A (zh) * | 2011-09-27 | 2013-04-03 | 上海航天控制工程研究所 | 一种微机械陀螺(mems)在深空探测中的使用方法 |
WO2013063909A1 (zh) * | 2011-11-03 | 2013-05-10 | 中国科学院自动化研究所 | 一种基于遗传算法的惯性/地磁传感器标定方法 |
RU2011148682A (ru) * | 2011-11-30 | 2013-06-10 | Общество с ограниченной ответственностью "Научно-производственная фирма "Эпсилон" | Способ повышения точности калибровки трехосных лазерных гироскопов с одним общим вибратором |
US20140104445A1 (en) * | 2012-10-12 | 2014-04-17 | Qualcomm Incorporated | Gyroscope conditioning and gyro-camera alignment |
CN104729530A (zh) * | 2014-10-13 | 2015-06-24 | 北京航空航天大学 | 一种惯导/北斗组合导航的半实物仿真方法 |
WO2015200850A2 (en) * | 2014-06-26 | 2015-12-30 | Lumedyne Technologies Incorporated | System and methods for determining rotation from nonlinear periodic signals |
CN105487402A (zh) * | 2014-09-17 | 2016-04-13 | 上海新跃仪表厂 | 一种星敏与陀螺组合定姿全物理仿真试验方法 |
CN106200614A (zh) * | 2016-07-15 | 2016-12-07 | 北京控制工程研究所 | 一种采用控制力矩陀螺真实力矩的航天器姿控测试系统及方法 |
CN106227935A (zh) * | 2016-07-20 | 2016-12-14 | 浙江大学 | 一种卫星姿态控制一体化仿真的系统及实现方法 |
CN107228672A (zh) * | 2017-06-27 | 2017-10-03 | 上海航天控制技术研究所 | 一种适用于姿态机动工况下的星敏和陀螺数据融合方法 |
CN108803373A (zh) * | 2018-04-16 | 2018-11-13 | 上海微小卫星工程中心 | 一种三轴转台的地速消除方法 |
CN110723309A (zh) * | 2019-08-27 | 2020-01-24 | 广东工业大学 | 一种四旋翼无人机转动惯量测量方法 |
-
2020
- 2020-02-24 CN CN202010112786.1A patent/CN111207775B/zh active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5205518A (en) * | 1991-11-25 | 1993-04-27 | General Electric Co. | Gyroless yaw control system for a three axis stabilized, zero-momentum spacecraft |
CN101214859A (zh) * | 2007-12-26 | 2008-07-09 | 北京控制工程研究所 | 一种变轨期间自主故障检测恢复控制的方法 |
CN101219713A (zh) * | 2007-12-26 | 2008-07-16 | 北京控制工程研究所 | 一种卫星的自主变轨方法 |
CN103017745A (zh) * | 2011-09-27 | 2013-04-03 | 上海航天控制工程研究所 | 一种微机械陀螺(mems)在深空探测中的使用方法 |
WO2013063909A1 (zh) * | 2011-11-03 | 2013-05-10 | 中国科学院自动化研究所 | 一种基于遗传算法的惯性/地磁传感器标定方法 |
RU2011148682A (ru) * | 2011-11-30 | 2013-06-10 | Общество с ограниченной ответственностью "Научно-производственная фирма "Эпсилон" | Способ повышения точности калибровки трехосных лазерных гироскопов с одним общим вибратором |
US20140104445A1 (en) * | 2012-10-12 | 2014-04-17 | Qualcomm Incorporated | Gyroscope conditioning and gyro-camera alignment |
WO2015200850A2 (en) * | 2014-06-26 | 2015-12-30 | Lumedyne Technologies Incorporated | System and methods for determining rotation from nonlinear periodic signals |
CN105487402A (zh) * | 2014-09-17 | 2016-04-13 | 上海新跃仪表厂 | 一种星敏与陀螺组合定姿全物理仿真试验方法 |
CN104729530A (zh) * | 2014-10-13 | 2015-06-24 | 北京航空航天大学 | 一种惯导/北斗组合导航的半实物仿真方法 |
CN106200614A (zh) * | 2016-07-15 | 2016-12-07 | 北京控制工程研究所 | 一种采用控制力矩陀螺真实力矩的航天器姿控测试系统及方法 |
CN106227935A (zh) * | 2016-07-20 | 2016-12-14 | 浙江大学 | 一种卫星姿态控制一体化仿真的系统及实现方法 |
CN107228672A (zh) * | 2017-06-27 | 2017-10-03 | 上海航天控制技术研究所 | 一种适用于姿态机动工况下的星敏和陀螺数据融合方法 |
CN108803373A (zh) * | 2018-04-16 | 2018-11-13 | 上海微小卫星工程中心 | 一种三轴转台的地速消除方法 |
CN110723309A (zh) * | 2019-08-27 | 2020-01-24 | 广东工业大学 | 一种四旋翼无人机转动惯量测量方法 |
Non-Patent Citations (13)
Title |
---|
DC WOFFINDEN 等: "Relative angles-only navigation and pose estimation for autonomous orbital rendezvous" * |
JIAO REN 等: "Cascade PID controller for quadrotor" * |
R. WILLIAMSON 等: "Detecting absolute human knee angle and angular velocity using accelerometers and rate gyroscopes" * |
于会 等: "基于RBF神经网络的微机械陀螺仪温度控制系统" * |
何亮: "面向立方体卫星的姿态控制策略与控制算法研究" * |
刘宇 等: "基于振动陀螺的动态零位校准及其误差补偿" * |
刘成阳: "光纤陀螺的数字闭环控制方案的设计与实现" * |
周本川 等: "弹载捷联惯导系统的在线标定方法" * |
孙兆伟 等: "微小卫星高精度三轴稳定控制算法研究" * |
李荣 等: "轴向速度组合的陀螺随钻测量方法研究" * |
胡梦纯 等: "三轴转台误差对光纤捷联惯组陀螺标定精度的影响" * |
袁彦红 等: "利用高精度陀螺对星敏感器在轨标定算法研究" * |
陈杰春 等: "一种标定陀螺仪的新方法" * |
Also Published As
Publication number | Publication date |
---|---|
CN111207775B (zh) | 2023-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109459065B (zh) | 一种基于卫星惯性空间旋转姿态的陀螺安装矩阵标定方法 | |
CN105352527B (zh) | 一种基于双轴转位机构光纤陀螺标定方法 | |
CN103983276A (zh) | 一种基于导航基准系的三框架四轴惯性平台误差标定方法 | |
CN110986928B (zh) | 光电吊舱三轴陀螺仪漂移实时修正方法 | |
CN112461261B (zh) | 一种校正mems陀螺仪温度漂移的装置和方法 | |
CN111623770B (zh) | 一种基于速度误差开环修正的提高惯性制导精度的方法 | |
CN103017765A (zh) | 应用于微机械组合导航系统的偏航角修正方法和修正装置 | |
CN111780753B (zh) | 一种姿态误差反馈修正的提高惯性制导精度的方法 | |
CN101226206A (zh) | 一种对加速度计零位偏差进行在轨标定的方法 | |
CN110895418B (zh) | 补偿舵机动力学滞后的低速旋转飞行器控制方法及系统 | |
CN112797979B (zh) | 一种应用于agv的惯性姿态导航系统 | |
CN110220415B (zh) | 制导弹药外弹道闭环修正模拟平台及模拟方法 | |
CN111780752A (zh) | 一种姿态误差可观测的提高惯性制导精度方法 | |
CN113419431B (zh) | 一种基于事件触发的平流层飞艇轨迹跟踪控制方法及系统 | |
CN111561948A (zh) | 四轴冗余捷联惯导的系统级标定方法 | |
CN112525144B (zh) | 一种非线性姿态检测补偿方法及终端 | |
CN110940357B (zh) | 一种用于旋转惯导单轴自对准的内杆臂标定方法 | |
CN113532481A (zh) | 一种mems惯性导航系统的零偏补偿方法 | |
CN111207775A (zh) | 卫星gnc系统半物理仿真闭环测试陀螺零位标定方法 | |
CN111637883B (zh) | 基于速度和位置误差闭环反馈的提高惯性制导精度的方法 | |
CN115586724B (zh) | 一种齿轮巡检机器人系统自适应分数阶全局滑模控制方法 | |
CN110006455A (zh) | 用于冗余惯导系统中加速度计误差参数的快速标定方法 | |
CN112230540B (zh) | 一种基于过载控制的飞行器侧向并行质心控制方法 | |
CN110285711B (zh) | 基于信息物理系统的制导弹药飞行姿态修正模拟系统及方法 | |
CN111089606B (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 |