CN107296598B - 基于光电传感器的心率测定方法及装置 - Google Patents

基于光电传感器的心率测定方法及装置 Download PDF

Info

Publication number
CN107296598B
CN107296598B CN201710480340.2A CN201710480340A CN107296598B CN 107296598 B CN107296598 B CN 107296598B CN 201710480340 A CN201710480340 A CN 201710480340A CN 107296598 B CN107296598 B CN 107296598B
Authority
CN
China
Prior art keywords
sample data
pulse signal
group
waveform
photoelectric
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
CN201710480340.2A
Other languages
English (en)
Other versions
CN107296598A (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.)
Wuxi Etek Microelectronics Co ltd
Original Assignee
Wuxi Etek Microelectronics Co ltd
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 Wuxi Etek Microelectronics Co ltd filed Critical Wuxi Etek Microelectronics Co ltd
Priority to CN201710480340.2A priority Critical patent/CN107296598B/zh
Publication of CN107296598A publication Critical patent/CN107296598A/zh
Application granted granted Critical
Publication of CN107296598B publication Critical patent/CN107296598B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02416Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
    • A61B5/02427Details of sensor
    • A61B5/02433Details of sensor for infrared radiation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Physiology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Cardiology (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

本发明公开了一种心率测定方法和装置,所述方法包括:光电传感器感应皮肤中血管的脉动得到光电感应脉冲信号;每预定时长采样所述光电感应脉冲信号得到一组样本数据;根据采集的各组样本数据得到所述光电感应脉冲信号的波形的斜率,直到得到所述光电感应脉冲信号的波形的最高点和/或最低点;根据所述光电感应脉冲信号的波形的最高点和/或最低点得到所述心跳的心率。这样不需要采集电路的精准度,在实际使用中能达到90%以上的准确率。

Description

基于光电传感器的心率测定方法及装置
【技术领域】
本发明涉及在心率测定领域,尤其涉及一种基于光电传感器的心率测定方法及装置。
【背景技术】
心率测定装置,在医疗领域应用已经非常的广泛,但是这种装置非常专业,只能限于在医院或医疗机构使用。
目前,智能电子设备,比如智能手机、智能手环等,被赋予了越来越多的额外的功能,比如测量心率。这种测量心率的方式最终的结果可能不是非常的精确,但是其可以给用户非常好的体验,也能够粗略的了解到心跳的情况。
然而,传统的光电心率算法都是采用包络检测的方法,包络检测方法虽然具有稳定度高,误码率低的优点,但对芯片的机能要求非常高,需要能有快速傅里叶的能力。这是大部分传统单片机不能胜任的。还有一种方法是利用判断心电脉冲的峰值点来检测,但这种方法虽然快速,但光电信号脉冲在不同的客户的感应范围不稳定,对采样电路的要求也比较高,所以准确度并不是很高。
因此,有必要提出一种解决方案来解决上述问题。
【发明内容】
本发明要解决的技术问题在于提供一种基于光电传感器的心率测定方法及装置,其利用判断脉冲的上升下降沿的斜率变化来实现光电传感器的简易快速算法。
为了解决上述问题,根据本发明的一个方面,本发明提供了一种一种心率测定方法,其特征在于,其包括:光电传感器感应皮肤中血管的脉动得到光电感应脉冲信号;每预定时长采样所述光电感应脉冲信号得到一组样本数据;根据采集的各组样本数据得到所述光电感应脉冲信号的波形的斜率,直到得到所述光电感应脉冲信号的波形的最高点和/或最低点;根据所述光电感应脉冲信号的波形的最高点和/或最低点得到所述心跳的心率。
进一步的,采样周期为1ms-10ms,一组样本数据包括有16-256个样本数据。
进一步的,所述光电感应脉冲信号的波形的相邻的两个最高点或最低点之间的时间间隔为一个最小心跳周期,基于最小心跳周期得到所述心率。
进一步的,对于每组样本数据,确定这组样本数据的整体斜率是上升还是下降,从相邻的整体斜率由上升转为下降的两组样本数据中选择最大的一个样本数据作为所述光电感应脉冲信号的波形的最高点,从相邻的整体斜率有下降转为上升的两组样本数据中选择最小的一个样本数据作为所述光电感应脉冲信号的波形的最低点。
进一步的,如果一组样本数据中,有超过预定数目的样本数据相对于此前的样本数据是增大的,则认为该组样本数据的整体斜率是上升,如果一组样本数据中,有超过预定数目的样本数据相对于此前的样本数据是减小的,则认为该组样本数据的整体斜率是下升,如果一组样本数据中,其整体斜率既不是上升、也不是下降,并且其前一组样本数据的整体斜率为上升,下一组样本数据的整体斜率为下降,则从该组样本数据中选择最大的一个样本数据作为所述光电感应脉冲信号的波形的最高点,如果一组样本数据中,其整体斜率既不是上升、也不是下降,并且其前一组样本数据的整体斜率为下降,下一组样本数据的整体斜率为上升,则从该组样本数据中选择最小的一个样本数据作为所述光电感应脉冲信号的波形的最低点。
根据本发明的另一个方面,本发明提供一种心率测定装置,其包括:光电传感器,其用于感应皮肤中血管的脉动得到光电感应脉冲信号;采样模块,用于每预定时长采样所述光电感应脉冲信号得到一组样本数据;处理器,根据采集的各组样本数据得到所述光电感应脉冲信号的波形的斜率,直到得到所述光电感应脉冲信号的波形的最高点和/或最低点,根据所述光电感应脉冲信号的波形的最高点和/或最低点得到所述心跳的心率。
进一步的,采样周期为1ms-10ms,一组样本数据包括有16-256个样本数据。
进一步的,所述光电感应脉冲信号的波形的相邻的两个最高点或最低点之间的时间间隔为一个最小心跳周期,基于最小心跳周期得到所述心率。
与现有技术相比,本发明采用了判断脉冲斜率的方式来实现光电脉冲的快速判断,并不需要采集电路的精准度,在实际使用中能达到90%以上的准确率。
关于本发明的其他目的,特征以及优点,下面将结合附图在具体实施方式中详细描述。
【附图说明】
结合参考附图及接下来的详细描述,本发明将更容易理解,其中同样的附图标记对应同样的结构部件,其中:
图1(a)为光电传感器得到的光电感应脉冲信号的第一种波形;
图1(b)为光电传感器得到的光电感应脉冲信号的第二种波形;
图2为本发明中的心率测定方法在一个实施例中的流程示意图;和
图3为本发明中的心率测定装置在一个实施例中的结构示意图。
【具体实施方式】
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
此处所称的“一个实施例”或“实施例”是指与所述实施例相关的特定特征、结构或特性至少可包含于本发明至少一个实现方式中。在本说明书中不同地方出现的“在一个实施例中”并非必须都指同一个实施例,也不必须是与其他实施例互相排斥的单独或选择实施例。本发明中的“多个”、“若干”表示两个或两个以上。本发明中的“和/或”表示“和”或者“或”。
本发明提出一种基于光电传感器的心率测定方法及装置,其利用判断脉冲的上升下降沿的斜率变化来实现光电传感器的简易快速算法,并不需要采集电路的精准度,在实际使用中能达到90%的准确率。
图2为本发明中的心率测定方法200在一个实施例中的流程示意图。如图2所示的,所述心率测定方法200包括如下步骤。
步骤210,光电传感器感应皮肤中血管的脉动得到光电感应脉冲信号。
在光电感应脉冲方式中存在着两种脉冲波形,第一种是正常的波形,如图1(a),另一种是放大过渡的肖波波形,如图1(b)。由于每个人体的皮肤和佩戴方式不同,会存在两种波形共存。该心率测定方法200通常被应用于佩戴电子装置上,比如手环等。
步骤220,每预定时长采样所述光电感应脉冲信号得到一组样本数据。
每个采样周期会样所述光电感应脉冲信号会得到1个样本数据,采样周期为1ms-10ms。每组样本数据可以包括有16至256个样本数据,比如16、32、64等。
步骤230,根据采集的各组样本数据得到所述光电感应脉冲信号的波形的斜率,直到得到所述光电感应脉冲信号的波形的最高点和/或最低点。
对于每组样本数据,确定这组样本数据的整体斜率是上升还是下降。从相邻的整体斜率由上升转为下降的两组样本数据中选择最大的一个样本数据作为所述光电感应脉冲信号的波形的最高点,从相邻的整体斜率有下降转为上升的两组样本数据中选择最小的一个样本数据作为所述光电感应脉冲信号的波形的最低点。
在一个实施例中,如果一组样本数据中,有超过预定数目的样本数据相对于此前的样本数据是增大的,则认为该组样本数据的整体斜率是上升,如果一组样本数据中,有超过预定数目的样本数据相对于此前的样本数据是减小的,则认为该组样本数据的整体斜率是下升。比如所述预定数目为一组样本数据中样本数据的数目的90%。
如果一组样本数据中,其整体斜率既不是上升、也不是下降,并且其前一组样本数据的整体斜率为上升,下一组样本数据的整体斜率为下降,则从该组样本数据中选择最大的一个样本数据作为所述光电感应脉冲信号的波形的最高点。如果一组样本数据中,其整体斜率既不是上升、也不是下降,并且其前一组样本数据的整体斜率为下降,下一组样本数据的整体斜率为上升,则从该组样本数据中选择最小的一个样本数据作为所述光电感应脉冲信号的波形的最低点。
对数据精度有要求的可以采集最低点,一般只采集最高点就够了,有的信号最低点是一条平的直线较不容易采集到最低点。
步骤240,根据所述光电感应脉冲信号的波形的最高点和/或最低点得到所述心跳的心率。
所述光电感应脉冲信号的波形的相邻的两个最高点或最低点之间的时间间隔为一个最小心跳周期,基于最小心跳周期得到所述心率,比如最小心跳周期是800ms,那么心率是60/0.8=75次/秒。
通常,至少要采集4秒钟的样本数据才能够得到心跳的心率。
根据本发明的另一个方面,本发明提出一种心率测定装置。图3为本发明中的心率测定装置在一个实施例中的结构示意图。
如图3所示的,所述心率测定装置包括光电传感器310、采样模块320、处理器330。
所述光电传感器310用于感应皮肤中血管的脉动得到光电感应脉冲信号。
所述采样模块320用于每预定时长采样所述光电感应脉冲信号得到一组样本数据。采样模块320每个采样周期会样所述光电感应脉冲信号会得到1个样本数据,采样周期可以为1ms-10ms。每组样本数据可以包括有16至256个样本数据,比如16、32、64等。
所述处理器330根据采集的各组样本数据得到所述光电感应脉冲信号的波形的斜率,直到得到所述光电感应脉冲信号的波形的最高点和/或最低点,根据所述光电感应脉冲信号的波形的最高点和/或最低点得到所述心跳的心率。
具体的,对于每组样本数据,确定这组样本数据的整体斜率是上升还是下降。从相邻的整体斜率由上升转为下降的两组样本数据中选择最大的一个样本数据作为所述光电感应脉冲信号的波形的最高点,从相邻的整体斜率有下降转为上升的两组样本数据中选择最小的一个样本数据作为所述光电感应脉冲信号的波形的最低点。
在一个实施例中,如果一组样本数据中,有超过预定数目的样本数据相对于此前的样本数据是增大的,则认为该组样本数据的整体斜率是上升,如果一组样本数据中,有超过预定数目的样本数据相对于此前的样本数据是减小的,则认为该组样本数据的整体斜率是下升。比如所述预定数目为一组样本数据中样本数据的数目的90%。
如果一组样本数据中,其整体斜率既不是上升、也不是下降,并且其前一组样本数据的整体斜率为上升,下一组样本数据的整体斜率为下降,则从该组样本数据中选择最大的一个样本数据作为所述光电感应脉冲信号的波形的最高点。如果一组样本数据中,其整体斜率既不是上升、也不是下降,并且其前一组样本数据的整体斜率为下降,下一组样本数据的整体斜率为上升,则从该组样本数据中选择最小的一个样本数据作为所述光电感应脉冲信号的波形的最低点。对数据精度有要求的可以采集最低点,一般只采集最高点就够了,有的信号最低点是一条平的直线较不容易采集到最低点。
所述光电感应脉冲信号的波形的相邻的两个最高点或最低点之间的时间间隔为一个最小心跳周期,基于最小心跳周期得到所述心率,比如最小心跳周期是800ms,那么心率是60/0.8=75次/秒。
实际使用中,可以根据情况多取几个最大值或最小值,这样就能完成心率的计算。通常,至少要采集4秒钟的样本数据才能够得到心跳的心率。
上述说明已经充分揭露了本发明的具体实施方式。需要指出的是,熟悉该领域的技术人员对本发明的具体实施方式所做的任何改动均不脱离本发明的权利要求书的范围。相应地,本发明的权利要求的范围也并不仅仅局限于前述具体实施方式。

Claims (4)

1.一种心率测定方法,其特征在于,其包括:
光电传感器感应皮肤中血管的脉动得到光电感应脉冲信号;
每预定时长采样所述光电感应脉冲信号得到一组样本数据;
根据采集的各组样本数据得到所述光电感应脉冲信号的波形的趋势,直到得到所述光电感应脉冲信号的波形的多个最高点和/或多个最低点,
对于每组样本数据,确定这组样本数据的整体趋势是上升还是下降,
从相邻的整体趋势由上升转为下降的两组样本数据中选择最大的一个样本数据作为所述光电感应脉冲信号的波形的最高点,从相邻的整体趋势由下降转为上升的两组样本数据中选择最小的一个样本数据作为所述光电感应脉冲信号的波形的最低点,
如果一组样本数据中,有超过预定数目的样本数据相对于前面相邻的一个样本数据是增大的,则认为该组样本数据的整体趋势是上升,如果一组样本数据中,有超过预定数目的样本数据相对于前面相邻的一个样本数据是减小的,则认为该组样本数据的整体趋势是下升,
如果一组样本数据中,其整体趋势既不是上升、也不是下降,并且其前一组样本数据的整体趋势为上升,下一组样本数据的整体趋势为下降,则从该组样本数据中选择最大的一个样本数据作为所述光电感应脉冲信号的波形的最高点,
如果一组样本数据中,其整体趋势既不是上升、也不是下降,并且其前一组样本数据的整体趋势为下降,下一组样本数据的整体趋势为上升,则从该组样本数据中选择最小的一个样本数据作为所述光电感应脉冲信号的波形的最低点;
根据所述光电感应脉冲信号的波形的多个最高点和/或多个最低点得到心跳的心率,所述光电感应脉冲信号的波形的相邻的两个最高点或最低点之间的时间间隔为一个最小心跳周期,基于最小心跳周期得到所述心率。
2.根据权利要求1所述的心率测定方法,其特征在于,
采样周期为1ms-10ms,一组样本数据包括有16-256个样本数据。
3.一种心率测定装置,其特征在于,其包括:
光电传感器,其用于感应皮肤中血管的脉动得到光电感应脉冲信号;
采样模块,用于每预定时长采样所述光电感应脉冲信号得到一组样本数据;
处理器,根据采集的各组样本数据得到所述光电感应脉冲信号的波形的趋势,直到得到所述光电感应脉冲信号的波形的多个最高点和/或多个最低点,对于每组样本数据,确定这组样本数据的整体趋势是上升还是下降,
从相邻的整体趋势由上升转为下降的两组样本数据中选择最大的一个样本数据作为所述光电感应脉冲信号的波形的最高点,从相邻的整体趋势由下降转为上升的两组样本数据中选择最小的一个样本数据作为所述光电感应脉冲信号的波形的最低点,
如果一组样本数据中,有超过预定数目的样本数据相对于前面相邻的一个样本数据是增大的,则认为该组样本数据的整体趋势是上升,如果一组样本数据中,有超过预定数目的样本数据相对于前面相邻的一个样本数据是减小的,则认为该组样本数据的整体趋势是下升,
如果一组样本数据中,其整体趋势既不是上升、也不是下降,并且其前一组样本数据的整体趋势为上升,下一组样本数据的整体趋势为下降,则从该组样本数据中选择最大的一个样本数据作为所述光电感应脉冲信号的波形的最高点,
如果一组样本数据中,其整体趋势既不是上升、也不是下降,并且其前一组样本数据的整体趋势为下降,下一组样本数据的整体趋势为上升,则从该组样本数据中选择最小的一个样本数据作为所述光电感应脉冲信号的波形的最低点,
根据所述光电感应脉冲信号的波形的多个最高点和/或多个最低点得到心跳的心率,所述光电感应脉冲信号的波形的相邻的两个最高点或最低点之间的时间间隔为一个最小心跳周期,基于最小心跳周期得到所述心率。
4.根据权利要求3所述的心率测定装置,其特征在于,
采样周期为1ms-10ms,一组样本数据包括有16-256个样本数据。
CN201710480340.2A 2017-06-22 2017-06-22 基于光电传感器的心率测定方法及装置 Active CN107296598B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710480340.2A CN107296598B (zh) 2017-06-22 2017-06-22 基于光电传感器的心率测定方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710480340.2A CN107296598B (zh) 2017-06-22 2017-06-22 基于光电传感器的心率测定方法及装置

Publications (2)

Publication Number Publication Date
CN107296598A CN107296598A (zh) 2017-10-27
CN107296598B true CN107296598B (zh) 2020-12-25

Family

ID=60135065

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710480340.2A Active CN107296598B (zh) 2017-06-22 2017-06-22 基于光电传感器的心率测定方法及装置

Country Status (1)

Country Link
CN (1) CN107296598B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103549945A (zh) * 2013-10-31 2014-02-05 广州视源电子科技股份有限公司 通过心脏收缩过程特征识别脉率和血氧饱和度的方法
CN105078438A (zh) * 2015-06-19 2015-11-25 京东方科技集团股份有限公司 脉搏周期检测设备和方法和可穿戴电子设备
CN106175741A (zh) * 2016-07-13 2016-12-07 北京贞正物联网技术有限公司 一种用于足部的心率检测方法、检测系统及其鞋垫、智能鞋

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101427917B (zh) * 2008-09-08 2011-02-09 电子科技大学 基于固有趋势子序列模式分解的出现新的心脏活动趋势的ecg信号获取方法
CN101766484B (zh) * 2010-01-18 2011-09-07 董军 用于心电图识别与分类的方法和设备
US10531801B2 (en) * 2013-09-09 2020-01-14 Koninklijke Philips N.V. Fetal heart rate extraction from maternal abdominal ECG recordings
JP6315633B2 (ja) * 2014-09-09 2018-04-25 日本電信電話株式会社 心拍検出方法および心拍検出装置
CN105962920B (zh) * 2016-04-20 2019-06-11 广州视源电子科技股份有限公司 血压脉率检测方法及其系统
CN106373111A (zh) * 2016-08-30 2017-02-01 苏州品诺维新医疗科技有限公司 一种心率检测装置及获取心率变化波形的方法
CN106580301B (zh) * 2016-12-21 2019-09-24 广州心与潮信息科技有限公司 一种生理参数的监测方法、装置和手持设备

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103549945A (zh) * 2013-10-31 2014-02-05 广州视源电子科技股份有限公司 通过心脏收缩过程特征识别脉率和血氧饱和度的方法
CN105078438A (zh) * 2015-06-19 2015-11-25 京东方科技集团股份有限公司 脉搏周期检测设备和方法和可穿戴电子设备
CN106175741A (zh) * 2016-07-13 2016-12-07 北京贞正物联网技术有限公司 一种用于足部的心率检测方法、检测系统及其鞋垫、智能鞋

Also Published As

Publication number Publication date
CN107296598A (zh) 2017-10-27

Similar Documents

Publication Publication Date Title
EP3478166B1 (en) On-demand heart rate estimation based on optical measurements
US9936888B2 (en) Pulse period calculation device and biosensor equipped with the same
US11006885B2 (en) Apparatus for determining blood pressure
CN109872820B (zh) 一种无袖带血压测量方法、装置、设备及存储介质
US20200170577A1 (en) System and Methods for Adaptive Noise Quantification in Dynamic Biosignal Analysis
EP3387990A8 (en) System and method for heart rate detection with motion artifact reduction
EP2752154B1 (en) Method and system for obtaining physiological signal period
US9788737B2 (en) Vital signs measurement system, detecting method of the vital signs measurement system, and vital signs measurement earphone
EP3417774A1 (en) System and method for monitoring a heart rate
CN104000581B (zh) 心电信号处理方法及装置
CN105748056A (zh) 血压检测方法及其系统
CN107898443A (zh) 重搏波检测方法、装置及计算机存储介质
CN104873184A (zh) 一种测量心率的方法及系统
CN107296598B (zh) 基于光电传感器的心率测定方法及装置
TW201811261A (zh) 訊號偵測方法
CN106214143B (zh) 瞬时运动干扰识别方法及心率置信度计算方法
US20060167366A1 (en) Method and apparatus for extracting biological signal such as heartbeat or respiration
CN116369888B (zh) 一种非接触式心率变异性数据获取方法和装置
CN108814642A (zh) 一种电子听诊器的心音定位及心率计算方法
WO2018152713A1 (zh) 一种用于血压测量装置的数据处理方法
Tanaka et al. Wearable health monitoring system by using fuzzy logic heart-rate extraction
CN107747950B (zh) 记步方法与装置
US20200069201A1 (en) Heart rate detection system and wearable device using the same
CN105726001B (zh) 一种血压测量方法及装置
CN204813858U (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