CN111060150B - 基于mems加速度传感器和气压传感器的运动载体监测方法 - Google Patents

基于mems加速度传感器和气压传感器的运动载体监测方法 Download PDF

Info

Publication number
CN111060150B
CN111060150B CN201911232727.1A CN201911232727A CN111060150B CN 111060150 B CN111060150 B CN 111060150B CN 201911232727 A CN201911232727 A CN 201911232727A CN 111060150 B CN111060150 B CN 111060150B
Authority
CN
China
Prior art keywords
carrier
acceleration
acceleration sensor
air pressure
speed
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
Application number
CN201911232727.1A
Other languages
English (en)
Other versions
CN111060150A (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201911232727.1A priority Critical patent/CN111060150B/zh
Publication of CN111060150A publication Critical patent/CN111060150A/zh
Application granted granted Critical
Publication of CN111060150B publication Critical patent/CN111060150B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D21/00Measuring or testing not otherwise provided for
    • G01D21/02Measuring two or more variables by means not covered by a single other subclass
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D18/00Testing or calibrating apparatus or arrangements provided for in groups G01D1/00 - G01D15/00

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Fluid Pressure (AREA)

Abstract

本发明公开了一种基于MEMS加速度传感器和气压传感器的运动载体监测方法,将MEMS加速度传感器与气压传感器固连在运动载体上,利用消抖滤波的方法采集加速度传感器信号通过设定运行状态标志位,以此划分载体运行状态,以加速度积分计算速度变化量和高度变化量,基于识别载体运行状态中的制动过程的时间和制动末端的速度信息,修正该运行过程中由于MEMS加速度传感器零偏误差造成的积分误差,最后与气压传感器解算高度进行综合,最终实现基于加速度传感器和气压传感器,对目标载体的速度、位置的精确测量与监测。

Description

基于MEMS加速度传感器和气压传感器的运动载体监测方法
技术领域
本发明涉及MEMS传感器信号处理方法技术领域,尤其涉及一种基于MEMS加速度传感器和气压传感器的运动载体监测方法。
背景技术
近年来,随着微电子技术、集成电路技术和加工工艺的蓬勃发展,MEMS传感器(微型电子机械系统)以体积小、重量轻、功耗低、价格便宜的特点,极大促进了传感器的微型化、智能化。
加速度传感器是一种能够测量加速度的MEMS传感器,通常由质量块、阻尼器、弹性元件、敏感元件和适调电路等部分组成。将加速度传感器构成的惯性测量单元直接固连在运动载体上,可以感受运载体运行过程中的加速度变化并且将其转化为可用的输出信号。但是由于加速度MEMS传感器精度低,干扰大,在解算过程中需要对加速度做二次积分处理,引入了累积误差的问题,此外由于数据采集过程中出现LINUX系统定时器时间不够精准,导致积分时间的减少,都对最终的监测结果具有较大的影响。气压传感器可以提供载体运行时气压变化,通过气压变化可以转化为高度变化,然而气压变化只能解算高度变化,无法判断速度变化,这就限制了载体运行时的速度数据显示,不具备完整的监测功能。因此,针对MEMS加速度传感器对目标载体运行过程进行监测算法研究,具有非常重要的意义。
发明内容
本发明所要解决的技术问题是针对背景技术中所涉及到的缺陷,提供一种基于MEMS加速度传感器和气压传感器的运动载体监测方法,减小加速度传感器易受干扰且精度低所带来的影响。
本发明为解决上述技术问题采用以下技术方案:
基于MEMS加速度传感器和气压传感器的运动载体监测方法,所述载体具有静止、加速启动、匀速运动、减速制动、制动停止或者静止、加速启动、匀速运行、制动停止的周期运动特征,包括以下步骤:
步骤1),在载体上设置用于测量载体加速度的MEMS加速度传感器和用于测量载体所在处大气压的气压传感器,并设定MEMS加速度传感器以及气压传感器数据输出方式,令速度V为0,位移S为0,补偿量X为0,时间参数γ为0;
步骤2),设定MEMS加速度传感器数据积分周期时间T;
步骤3),采集MEMS加速度传感器静止时的三个轴向均值加速度,设定MEMS加速度传感器相应于运动载体运动轴的轴向为主轴向,令MEMS加速度传感器静止时其主轴向的均值加速度为M;
步骤4),从气压传感器获取气压数据P后解算成高度数据H;
步骤5),将MEMS加速度传感器主轴向的加速度与M作差,得到主轴向加速度的变化N,并将包含当前周期在内的前10周期T的加速度变化值N的均值作为当前时刻主轴向加速度的启动值A;
步骤6),根据A判断载体处于周期运动特征中的阶段位置,将A分别和预先设定的加速度阈值α、β比较,判断载体运动状态:
若A大于等于α:此时载体处于周期运动特征中的加速启动阶段;
若A小于α并大于等于β时,此时载体处于周期运动特征中的匀速运行阶段或静止阶段;
若A小于β,此时载体处于周期运动特征中的减速制动阶段;
步骤7),如果载体处于加速启动阶段:
步骤7.1),以周期T对N进行积分得到速度V与当前位移S,并记录载体运行周期中加速的时间T1;
步骤7.2),令速度V等于加速度N以周期T积分后的结果减去0.1X,对当前速度进行补偿;
步骤7.3),更新补偿量X,补偿量X为当前位移S与H的差值;
步骤8),如果载体处于匀速运行或静止阶段,将位移S以速度V与周期T相乘后进行累加,若此时V不为0,记录载体运行周期中匀速运行阶段时间T2;
步骤9),如果载体处于减速阶段,将时间参数γ加1:
步骤9.1),若γ小于预设的制动时间变化量ε:
步骤9.1.1),以周期T对N进行积分得到速度V与当前位移S;
步骤9.1.2),令速度V等于加速度N以周期T积分后的结果减去0.1X,对当前速度进行补偿;
步骤9.1.3),更新补偿量X,补偿量X为当前位移S与H的差值;
步骤9.2),若γ等于ε,此刻时间为载体制动阶段结束时间,记录载体运行周期中减速制动时间T3;
步骤9.3),若γ大于ε:
步骤9.3.1),将γ置零,以周期T对N进行积分得到速度V与当前位移S;
步骤9.3.2),将速度V乘以 (T3+T2+T1)得到位移补偿量s;
步骤9.3.3),将位移S和s作差得载体一次周期运动过程中的最终位移;
步骤9.3.4),将速度V置零;
步骤10),跳转执行步骤4)。
作为本发明一种基于气压传感器辅助的MEMS加速度传感器运动监测方法进一步的优化方案,步骤3)中求取静止时刻加速度传感器三个轴向均值加速度并设定主轴向的步骤如下:
读取MEMS加速度传感器三个轴向的加速度数据,如果其中有两轴的加速度小于α并大于β,读取不少于200个积分周期时间T的MEMS加速度传感器三轴数据,计算MEMS加速度传感器的三个轴向均值加速度;然后将绝对值最接近重力加速度g的轴向均值加速度所在的轴线设定为相应于运动载体的运动轴的轴向、即主轴向。
作为本发明一种基于气压传感器辅助的MEMS加速度传感器运动监测方法进一步的优化方案,所述MEMS加速度传感器采用ADXL362加速度传感器,所述气压传感器采用BMP280气压传感器。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
本发明针对MEMS加速度传感器传感器进行载体运动监测,采用控制载体运行过程加速时间以及运行过程减速制动时间,利用气压高度数据对加速度信号进行反馈控制的同时并利用载体停止运行时的残留速度,输出位置速度信息。本发明很好的减小了加速度传感器易受干扰所带来的噪声影响,在完成采集加速度信息后通过运动状态的识别,控制运动过程中的加速阶段和减速阶段,利用停止时的速度误差进行位移修正,可以大大提高监测数据精度,从而稳定可靠地实现载体运动监测。
附图说明
图1是本发明运动监测方法示意图;
图2是本发明设定MEMS加速度传感器主轴向处理流程图;
图3是本发明载体加速度数据初步处理流程图;
图4是本发明的实例载体示意图;
图5是本发明的实例所提供的在载体中运行实际监测结果。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
本发明可以以许多不同的形式实现,而不应当认为限于这里所述的实施例。相反,提供这些实施例以便使本公开透彻且完整,并且将向本领域技术人员充分表达本发明的范围。在附图中,为了清楚起见放大了组件。
现有的载体运动监测方法主要是通过高精度的传感器分析处理数据监测载体运动,本发明设计了一种基于MEMS加速度传感器和气压传感器的运动监测方法,具体设计思想如图1所示,通过设定系统周期,采集MEMS加速度传感器数据以及气压传感器数据存储,采用消抖滤波的方法对加速度数据处理减小传感器易受干扰的影响,通过识别载体运动状态过程,将气压数据作用于加速度数据形成反馈控制,并识别载体运行状态中的制动过程的时间和制动末端的速度信息,修正该运行过程中由于MEMS加速度传感器零偏误差造成的积分误差。
本发明所述的运动监测方法的具体步骤如下:
步骤1),在载体上设置用于测量载体加速度的MEMS加速度传感器和用于测量载体所在处大气压的气压传感器,并设定MEMS加速度传感器以及气压传感器数据输出方式,令速度V为0,位移S为0,补偿量X为0,时间参数γ为0;
步骤2),设定MEMS加速度传感器数据积分周期时间T;
步骤3),如图2,读取MEMS加速度传感器静止时三个轴向的加速度数据,如果其中有两轴的加速度小于α并大于β,读取不少于200个积分周期时间T的MEMS加速度传感器三轴数据,计算MEMS加速度传感器的三个轴向均值加速度;然后将绝对值最接近重力加速度g的轴向均值加速度所在的轴线设定为相应于运动载体的运动轴的轴向、即主轴向,令MEMS加速度传感器静止时其主轴向的均值加速度为M ;
步骤4),从气压传感器获取气压数据P后解算成高度数据H;
步骤5),如图3,将MEMS加速度传感器主轴向的加速度与M作差,得到主轴向加速度的变化N,并将包含当前周期在内的前10周期T的加速度变化值N的均值作为当前时刻主轴向加速度的启动值A;
步骤6),根据A判断载体处于周期运动特征中的阶段位置,将A分别和预先设定的加速度阈值α、β比较,判断载体运动状态:
若A大于等于α:此时载体处于周期运动特征中的加速启动阶段;
若A小于α并大于等于β时,此时载体处于周期运动特征中的匀速运行阶段或静止阶段;
若A小于β,此时载体处于周期运动特征中的减速制动阶段;
步骤7),如果载体处于加速启动阶段:
步骤7.1),以周期T对N进行积分得到速度V与当前位移S,并记录载体运行周期中加速的时间T1;
步骤7.2),令速度V等于加速度N以周期T积分后的结果减去0.1X,对当前速度进行补偿;
步骤7.3),更新补偿量X,补偿量X为当前位移S与H的差值;
步骤8),如果载体处于匀速运行或静止阶段,将位移S以速度V与周期T相乘后进行累加,若此时V不为0,记录载体运行周期中匀速运行阶段时间T2;
步骤9),如果载体处于减速阶段,将时间参数γ加1:
步骤9.1),若γ小于预设的制动时间变化量ε:
步骤9.1.1),以周期T对N进行积分得到速度V与当前位移S;
步骤9.1.2),令速度V等于加速度N以周期T积分后的结果减去0.1X,对当前速度进行补偿;
步骤9.1.3),更新补偿量X,补偿量X为当前位移S与H的差值;
步骤9.2),若γ等于ε,此刻时间为载体制动阶段结束时间,记录载体运行周期中减速制动时间T3;
步骤9.3),若γ大于ε:
步骤9.3.1),将γ置零,以周期T对N进行积分得到速度V与当前位移S;
步骤9.3.2),将速度V乘以 (T3+T2+T1)得到位移补偿量s;
步骤9.3.3),将位移S和s作差得载体一次周期运动过程中的最终位移;
步骤9.3.4),将速度V置零;
步骤10),跳转执行步骤4)。
本发明的实验效果:
本发明的实例载体对象如图4所示,利用MATLAB对载体对象为电梯时运行的实时数据分析结果如图5所示,从图中可以看出MEMS加速度传感器在上述方法的实现下,能够大幅度提高载体运行监测精度。本发明提供了对一类具有静止——加速启动——匀速运动——减速制动——制动停止(或者静止——加速启动——匀速运行——制动停止)等周期运动特征的运动载体采用MEMS传感器实现高精度运动监测的算法本技术领域技术人员可以理解的是,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (3)

1.基于MEMS加速度传感器和气压传感器的运动载体监测方法,所述载体具有静止、加速启动、匀速运行、减速制动、制动停止的周期运动特征,其特征在于,包括以下步骤:
步骤1),在载体上设置用于测量载体加速度的MEMS加速度传感器和用于测量载体所在处大气压的气压传感器,并设定MEMS加速度传感器以及气压传感器数据输出方式,令速度V为0,位移S为0,补偿量X为0,时间参数γ为0;
步骤2),设定MEMS加速度传感器数据积分周期时间T;
步骤3),采集MEMS加速度传感器静止时的三个轴向均值加速度,设定MEMS加速度传感器相应于运动载体运动轴的轴向为主轴向,令MEMS加速度传感器静止时其主轴向的均值加速度为M;
步骤4),从气压传感器获取气压数据P后解算成高度数据H;
步骤5),将MEMS加速度传感器主轴向的加速度与M作差,得到主轴向加速度变化值N,并将包含当前周期在内的前10周期T的加速度变化值N的均值作为当前时刻主轴向加速度的启动值A;
步骤6),根据A判断载体处于周期运动特征中的阶段位置,将A分别和预先设定的加速度阈值α、β比较,判断载体运动状态:
若A大于等于α:此时载体处于周期运动特征中的加速启动阶段;
若A小于α并大于等于β时,此时载体处于周期运动特征中的匀速运行阶段或静止阶段;
若A小于β,此时载体处于周期运动特征中的减速制动阶段;
步骤7),如果载体处于加速启动阶段:
步骤7.1),以周期T对N进行积分得到速度V与当前位移S,并记录载体运行周期中加速的时间T1;
步骤7.2),令速度V等于加速度N以周期T积分后的结果减去0.1X,对当前速度进行补偿;
步骤7.3),更新补偿量X,补偿量X为当前位移S与H的差值;
步骤8),如果载体处于匀速运行或静止阶段,将位移S以速度V与周期T相乘后进行累加,若此时V不为0,记录载体运行周期中匀速运行阶段时间T2;
步骤9),如果载体处于减速阶段,将时间参数γ加1:
步骤9.1),若γ小于预设的制动时间变化量ε:
步骤9.1.1),以周期T对N进行积分得到速度V与当前位移S;
步骤9.1.2),令速度V等于加速度N以周期T积分后的结果减去0.1X,对当前速度进行补偿;
步骤9.1.3),更新补偿量X,补偿量X为当前位移S与H的差值;
步骤9.2),若γ等于ε,此刻时间为载体减速制动阶段结束时间,记录载体运行周期中减速制动时间T3;
步骤9.3),若γ大于ε:
步骤9.3.1),将γ置零,以周期T对N进行积分得到速度V与当前位移S;
步骤9.3.2),将速度V乘以 (T3+T2+T1)得到位移补偿量s;
步骤9.3.3),将位移S和s作差得载体一次周期运动过程中的最终位移;
步骤9.3.4),将速度V置零;
步骤10),跳转执行步骤4)。
2.根据权利要求1所述的基于MEMS加速度传感器和气压传感器的运动载体监测方法,其特征在于,步骤3)中求取静止时刻加速度传感器三个轴向均值加速度并设定主轴向的步骤如下:
读取MEMS加速度传感器三个轴向的加速度数据,如果其中有两轴的加速度小于α并大于β,读取不少于200个积分周期时间T的MEMS加速度传感器三轴数据,计算MEMS加速度传感器的三个轴向均值加速度;然后将绝对值最接近重力加速度g的轴向均值加速度所在的轴线设定为主轴向。
3.根据权利要求1所述的基于MEMS加速度传感器和气压传感器的运动载体监测方法,其特征在于,所述MEMS加速度传感器采用ADXL362加速度传感器,所述气压传感器采用BMP280气压传感器。
CN201911232727.1A 2019-12-05 2019-12-05 基于mems加速度传感器和气压传感器的运动载体监测方法 Active CN111060150B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911232727.1A CN111060150B (zh) 2019-12-05 2019-12-05 基于mems加速度传感器和气压传感器的运动载体监测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911232727.1A CN111060150B (zh) 2019-12-05 2019-12-05 基于mems加速度传感器和气压传感器的运动载体监测方法

Publications (2)

Publication Number Publication Date
CN111060150A CN111060150A (zh) 2020-04-24
CN111060150B true CN111060150B (zh) 2021-04-20

Family

ID=70299669

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911232727.1A Active CN111060150B (zh) 2019-12-05 2019-12-05 基于mems加速度传感器和气压传感器的运动载体监测方法

Country Status (1)

Country Link
CN (1) CN111060150B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111776904B (zh) * 2020-07-16 2023-03-28 广州独角鲸信息科技有限公司 一种电梯运行数据采集器及其数据采集方法
CN114475718B (zh) * 2020-10-26 2023-11-21 中车株洲电力机车研究所有限公司 列车停车的控制方法、设备、列车和存储介质
CN113900159B (zh) * 2021-08-27 2022-11-04 山东大学 一种基于北斗的智慧城市车载气象采集监控装置和方法
CN119953994A (zh) * 2025-04-09 2025-05-09 杭州海康威视数字技术股份有限公司 基于加速度的楼层检测方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103743992A (zh) * 2014-01-24 2014-04-23 浙江省质量检测科学研究院 电器插座的耐久性测试系统和高效一致性测试方法
CN105712142A (zh) * 2016-03-22 2016-06-29 上海点络信息技术有限公司 一种电梯运行状态的检测系统及检测方法
CN106168790A (zh) * 2016-02-29 2016-11-30 华南理工大学 一种在线改变目标速度和位置的s形加减速控制方法
CN107344688A (zh) * 2016-05-05 2017-11-14 腾讯科技(深圳)有限公司 机器人乘坐电梯时的楼层监测方法和装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06106945A (ja) * 1992-09-24 1994-04-19 Nippondenso Co Ltd 減衰力可変ショックアブソーバ制御装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103743992A (zh) * 2014-01-24 2014-04-23 浙江省质量检测科学研究院 电器插座的耐久性测试系统和高效一致性测试方法
CN106168790A (zh) * 2016-02-29 2016-11-30 华南理工大学 一种在线改变目标速度和位置的s形加减速控制方法
CN105712142A (zh) * 2016-03-22 2016-06-29 上海点络信息技术有限公司 一种电梯运行状态的检测系统及检测方法
CN107344688A (zh) * 2016-05-05 2017-11-14 腾讯科技(深圳)有限公司 机器人乘坐电梯时的楼层监测方法和装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于加速度传感器的油井功图位移测量的算法设计与优化;王永举;《电子设计工程》;20131231;第21卷(第24期);第72-77页 *

Also Published As

Publication number Publication date
CN111060150A (zh) 2020-04-24

Similar Documents

Publication Publication Date Title
CN111060150B (zh) 基于mems加速度传感器和气压传感器的运动载体监测方法
JP5474938B2 (ja) 加速度計および誤差補償方法
CN108680189B (zh) 一种基于卡尔曼滤波的mems陀螺仪z轴零偏动态补偿方法
US7822550B2 (en) Environmental characteristic determination
US7938004B1 (en) Systems and methods for angular rate and position measurement
US20070100550A1 (en) Systems and methods for reducing vibration-induced errors in inertial sensors
EP1593931A4 (en) DIFFERENCE CORRECTION METHOD FOR POSTURE DETERMINATION INSTRUMENT AND DISPLACEMENT MEASUREMENT INSTRUMENT
US20140149061A1 (en) Inertia sensor and method for reducing operation error of the same
JP2011017688A (ja) 角度測定の方法およびそれを実施するための角度測定ジャイロシステム
CN110132271A (zh) 一种自适应卡尔曼滤波姿态估计算法
CN107505845A (zh) 一种提高倾斜镜控制系统扰动抑制能力的控制方法
Zwahlen et al. Open-loop and closed-loop high-end accelerometer platforms for high demanding applications
US11619492B2 (en) Sensor linearization based upon correction of static and frequency-dependent non-linearities
CN115655272A (zh) 基于mems加速度计零偏和标度因数的温度补偿方法及系统
CN106197376A (zh) 基于单轴mems惯性传感器的车身倾角测量方法
CN118275734A (zh) 一种微机电惯性测量单元加速度计温度补偿方法
CN118209136A (zh) 一种基于在线mems陀螺仪偏置估计的温度补偿方法
JP2002131077A (ja) 移動体の停止判定方法及び停止判定装置ならびに停止判定プログラムを記録した記憶媒体
JP2024043518A (ja) Mems加速度計システム、およびアビオニクスシステム
CN113865621B (zh) 任意六位置MEMS陀螺仪及其g值敏感系数标定方法
CN116026328A (zh) 微惯导的零偏滞回效应补偿模型的构建方法和补偿方法
WO2015126910A1 (en) Instrument and method for measuring low indicated air speed
CN113859307B (zh) 车辆运行状态的检测方法、装置及车辆
CN110487529B (zh) 一种利用角速度传感器测量大展弦比翼面弯矩的方法
CN109788200B (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