CN105447306A - 心冲击信号周期的计算方法和装置 - Google Patents
心冲击信号周期的计算方法和装置 Download PDFInfo
- Publication number
- CN105447306A CN105447306A CN201510779087.1A CN201510779087A CN105447306A CN 105447306 A CN105447306 A CN 105447306A CN 201510779087 A CN201510779087 A CN 201510779087A CN 105447306 A CN105447306 A CN 105447306A
- Authority
- CN
- China
- Prior art keywords
- cycle
- energy
- signal
- heart impact
- impact signal
- 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.)
- Pending
Links
Classifications
-
- G06F19/321—
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7271—Specific aspects of physiological measurement analysis
Abstract
本发明揭示了一种心冲击信号周期的计算方法和装置,其中方法包括:获取心冲击微动传感器产生的心冲击信号,并将心冲击信号能量化,形成能量波;计算所述能量波变化的周期,其中,能量波变化的周期与所述心冲击信号周期相同。本发明的心冲击信号周期的计算方法和装置,将微动传感器采集的心冲击信号能量化,因为能量不分正负,所以微动传感器受到心冲击和阻尼的相互作用时,每一次心冲击的能量都是由最大值逐渐变小,直到下一次心冲击再次变大,而将心冲击信号能量化后,可以准确的找到能量波的波峰和波谷,通过简单的运算获取到能量波变化的周期,因为能量波变化的周期与所述心冲击信号周期相同,所以可以准确的得到心冲击信号的周期。
Description
技术领域
本发明涉及到心冲击信号的周期计算领域,特别是涉及到一种心冲击信号周期的计算方法和装置。
背景技术
微动传感器可以传感到呼吸信号、心冲击信号等信号,心冲击信号表现为周期性的包络,每个包络意味着一次心跳,我们可以通过假设每个包络都有一个维持一段较长时间不被超越的极大值和极小值,而且极值都出现在包络的相同相位,就此我们可以通过非常简单的程序获取心跳周期,但是存在两个缺点:第一是不被超越的时间常数决定了此方法测量的范围,假如设置为0.5秒,则测量范围低于120次/分,而且此常数不能设置太小,否则无法锁定周期内真正的极值;第二是无法保证极值总是在周期的相同相位出现,这也就无法保证每次周期测量的准确性。
发明内容
本发明的主要目的为提供一种可以准确测量心冲击信号周期的计算方法和装置。
为了实现上述发明目的,本发明首先提出一种心冲击信号周期的计算方法,包括:
获取心冲击微动传感器产生的心冲击信号,并将心冲击信号能量化,形成能量波;
计算所述能量波变化的周期,其中,能量波变化的周期与所述心冲击信号周期相同。
进一步地,所述将心冲击信号能量化,形成能量波的步骤包括:
按照预设的转化速率将心冲击信号数字化,获取对应的电压信号和电压变化率信号;
根据电压信号与电压变化率信号计算微动传感器的动能和弹性势能,其中微动传感器的弹性势能和动能,与电压信号和电压变化率信号的平方成线性关系。
进一步地,所述预设的转化速率为每秒转化100次~1000次。
进一步地,所述计算所述能量波变化的周期,其中,能量波变化的周期与所述心冲击信号周期相同的步骤之前,包括:
按照指定条件进行滤波处理。
进一步地,所述计算所述能量波变化的周期,其中,能量波变化的周期与所述心冲击信号周期相同的步骤包括:
根据能量波相邻的两个波峰或两个波谷计算能量波变化的周期。
本发明还提供一种心冲击信号周期的计算装置,包括:
能量转化单元,用于获取心冲击微动传感器产生的心冲击信号,并将心冲击信号能量化,形成能量波;
计算单元,用于计算所述能量波变化的周期,其中,能量波变化的周期与所述心冲击信号周期相同。
进一步地,所述能量转化单元包括:
数字转化模块,用于按照预设的转化速率将心冲击信号数字化,获取对应的电压信号和电压变化率信号;
计算模块,用于根据电压信号与电压变化率信号计算微动传感器的动能和弹性势能,其中微动传感器的弹性势能和动能,与电压信号和电压变化率信号的平方成线性关系。
进一步地,所述预设的转化速率为每秒转化100次~1000次。
进一步地,所述心冲击信号周期的计算装置还包括滤波单元,用于按照指定条件进行滤波处理。
进一步地,所述计算单元包括:
极值计算模块,用于根据能量波相邻的两个波峰或两个波谷计算能量波变化的周期。
本发明的心冲击信号周期的计算方法和装置,将微动传感器采集的心冲击信号能量化,可以得到对应的能量波,因为能量不分正负,所以微动传感器受到心冲击和阻尼的相互作用时,每一次心冲击的能量都是由最大值逐渐变小,直到下一次心冲击再次变大,而将心冲击信号能量化后,可以准确的找到能量波的极值,即能量波的波峰和波谷,从而通过简单的运算获取到能量波变化的周期,又因为能量波变化的周期与所述心冲击信号周期相同,所以可以准确的得到心冲击信号的周期。
附图说明
图1为本发明一实施例的心冲击信号周期的计算方法的流程图;
图2为本发明一实施例的心冲击信号能量化并形成能量波的流程图;
图3为本发明另一实施例的心冲击信号周期的计算方法的流程图;
图4为本发明一实施例的心冲击信号周期的计算装置的结构框图;
图5为本发明一实施例的能量转化单元的结构框图;
图6为本发明另一实施例的心冲击信号周期的计算装置的结构框图;
图7为本发明一实施例的计算单元的结构框图。
本发明目的的实现、功能特点及优点将结合实施例,参照附图做进一步说明。
具体实施方式
应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
参照图1,本发明实施例中提供一种心冲击信号周期的计算方法,包括步骤:
S1、获取心冲击微动传感器产生的心冲击信号,并将心冲击信号能量化,形成能量波;
S2、计算所述能量波变化的周期,其中,能量波变化的周期与所述心冲击信号周期相同。
上述步骤S1中所述,一个封闭的系统,其能量是守恒的,当外部对其做功时,可以改变该封闭系统的能量。本实施例中,微动传感器或者微动传感器及其包裹微动传感器的包裹物,可以形成一个近似于封闭的系统,当微动传感器或者微动传感器及其包裹微动传感器的包裹物贴附于人体心脏或周边时,心脏跳动会冲击微动传感器所在的近似封闭的系统,微动传感器会产生相应的心冲击信号,即微动传感器采集的微动信号,所述心冲击信号是一种包络信号,每一个包络信号表示一次心跳,该心冲击信号包含的能量一般为微动传感器的动能和势能,心脏跳动冲击微动传感器,微动传感器接受到的能量最大时,一般是微动传感器采集到的信号最大值时,即微动传感器弹性势能最大时,然后随着阻尼和冲击振荡的过程,能量会逐渐衰减,直到下一次心脏跳动冲击再升高到最大值,如此反复进行。所以将微动传感器采集到的心冲击信号能量化,产生的能量波可以反映心冲击信号的周期。能量波可以为三角波和锯齿波,或近似于三角波和锯齿波等。
如上述步骤S2所述,因为能量波的波形会规律的出现极大值和极小值,所以可以准确的计算出能量波的周期,而能量波的变化周期即为心冲击信号的周期。
本实施例中,上述心冲击信号周期的计算方法,将微动传感器采集的心冲击信号能量化,可以得到对应的能量波,因为能量不分正负,所以微动传感器受到心冲击和阻尼的相互作用时,每一次心冲击的能量都是由最大值逐渐变小,直到下一次心冲击再次变大,而将心冲击信号能量化后,可以准确的找到能量波的极值,即能量波的波峰和波谷,从而通过简单的运算获取到能量波变化的周期,又因为能量波变化的周期与所述心冲击信号周期相同,所以可以准确的得到心冲击信号的周期。
参照图2,本实施例中,上述步骤S1中将心冲击信号能量化,形成能量波的步骤包括:
S11、按照预设的转化速率将心冲击信号数字化,获取对应的电压信号和电压变化率信号;
S12、根据电压信号与电压变化率信号计算微动传感器的动能和弹性势能,其中微动传感器的弹性势能和动能,与电压信号和电压变化率信号的平方成线性关系。
如上述步骤S11所述,首先按照一定的转换速率将心冲击信号数字化,其主要体现在将微动传感器产生的电压信号由模拟信号转化为数字信号,方便进行数字化处理,然后通过相邻的两个电压值计算出电压变化率信号。
如上述步骤S12所述,上述的电压信号相当于微动传感器的弹性势能,电压变化率信号相当于微动传感器的动能,微动传感器的弹性势能和动能,与电压信号和电压变化率信号的平方成线性关系,通过这个关系可以计算出微动传感器的能量。
本实施例中,上述预设的转化速率为每秒转化100次~1000次,可以满足能量转化的速度,不会遗漏关键的能量点。
本实施例中,上述计算所述能量波变化的周期,其中,能量波变化的周期与所述心冲击信号周期相同的步骤S2之前,包括步骤:
S1`、按照指定条件进行滤波处理。
如上述步骤S1`所述,将指定条件的滤波是指对一些工频干扰、电器干扰等产生的干扰波过滤掉,提高心冲击信号周期的计算准确性。
参照图3,本实施例中,上述计算所述能量波变化的周期,其中,能量波变化的周期与所述心冲击信号周期相同的步骤S2包括:
S21、根据能量波相邻的两个波峰或两个波谷计算能量波变化的周期。
如上述步骤S21所述,即是顶点标记法,把最大能量值(波峰)作为一个周期的开始和上一周期的结束,或者把最小能量值(波谷)作为一个周期的开始和上一周期的结束就可以通过简单的运算获取准确的周期。
在一具体实施例中,获取微动传感器采集到的心冲击信号,然后按照每秒转化100次~1000次的速率将心冲击信号进行数字化,得到数字化的电压信号和电压变化率信号,然后通过电压信号和电压信号变化率与微动传感器的弹性势能和动能的线性关系,得到心冲击信号所具有的能量,形成能量波,然后对能量波进行过滤处理,该能量根据数字化转化率的大小,可能形成锯齿波或三角波等可以清晰得到能量波的多个波峰值和波谷值,然后将工频干扰、电器辐射等信号过滤掉,从而可以根据波峰或波谷计算心冲击信号的周期,比如,波峰可以作为一个周期的开始,上一个周期的结束等,计算简单准确。
参照图4,本发明实施例中,还提供一种心冲击信号周期的计算装置,包括:
能量转化单元10,用于获取心冲击微动传感器产生的心冲击信号,并将心冲击信号能量化,形成能量波;
计算单元20,用于计算所述能量波变化的周期,其中,能量波变化的周期与所述心冲击信号周期相同。
如上述能量转化单元10,一个封闭的系统,其能量是守恒的,当外部对其做功时,可以改变该封闭系统的能量。本实施例中,微动传感器或者微动传感器及其包裹微动传感器的包裹物,可以形成一个近似于封闭的系统,当微动传感器或者微动传感器及其包裹微动传感器的包裹物贴附于人体心脏或周边时,心脏跳动会冲击微动传感器所在的近似封闭的系统,微动传感器会产生相应的心冲击信号,即微动传感器采集的微动信号,所述心冲击信号是一种包络信号,每一个包络信号表示一次心跳,该心冲击信号包含的能量一般为微动传感器的动能和势能,心脏跳动冲击微动传感器,微动传感器接受到的能量最大时,一般是微动传感器采集到的信号最大值时,即微动传感器弹性势能最大时,然后随着阻尼和冲击振荡的过程,能量会逐渐衰减,直到下一次心脏跳动冲击再升高到最大值,如此反复进行。所以将微动传感器采集到的心冲击信号能量化,产生的能量波可以反映心冲击信号的周期。能量波可以为三角波和锯齿波,或近似于三角波和锯齿波等。
如上述计算单元20,因为能量波的波形会规律的出现极大值和极小值,所以可以准确的计算出能量波的周期,而能量波的变化周期即为心冲击信号的周期。
本实施例中,上述心冲击信号周期的计算装置,将微动传感器采集的心冲击信号能量化,可以得到对应的能量波,因为能量不分正负,所以微动传感器受到心冲击和阻尼的相互作用时,每一次心冲击的能量都是由最大值逐渐变小,直到下一次心冲击再次变大,而将心冲击信号能量化后,可以准确的找到能量波的极值,即能量波的波峰和波谷,从而通过简单的运算获取到能量波变化的周期,又因为能量波变化的周期与所述心冲击信号周期相同,所以可以准确的得到心冲击信号的周期。
参照图5,本实施例中,上述能量转化单元10包括:
数字转化模块11,用于按照预设的转化速率将心冲击信号数字化,获取对应的电压信号和电压变化率信号;
计算模块12,用于根据电压信号与电压变化率信号计算微动传感器的动能和弹性势能,其中微动传感器的弹性势能和动能,与电压信号和电压变化率信号的平方成线性关系。
如上述数字转化模块11,首先按照一定的转换速率将心冲击信号数字化,其主要体现在将微动传感器产生的电压信号由模拟信号转化为数字信号,方便进行数字化处理,然后通过相邻的两个电压值计算出电压变化率信号。
如上述计算模块12,上述的电压信号相当于微动传感器的弹性势能,电压变化率信号相当于微动传感器的动能,微动传感器的弹性势能和动能,与电压信号和电压变化率信号的平方成线性关系,通过这个关系可以计算出微动传感器的能量。
本实施例中,述预设的转化速率为每秒转化100次~1000次,可以满足能量转化的速度,不会遗漏关键的能量点。
参照图6,本实施例中,上述心冲击信号周期的计算装置,还包括滤波单元30,用于按照指定条件进行滤波处理。
如上述的滤波单元30,将指定条件的滤波是指对一些工频干扰、电器干扰等产生的干扰波过滤掉,提高心冲击信号周期的计算准确性。
参照图7,本实施例中,上述计算单元20包括:
极值计算模块21计算模块12,用于根据能量波相邻的两个波峰或两个波谷计算能量波变化的周期。即是顶点标记法,把最大能量值(波峰)作为一个周期的开始和上一周期的结束,或者把最小能量值(波谷)作为一个周期的开始和上一周期的结束就可以通过简单的运算获取准确的周期。
在一具体实施例中,心冲击信号周期的计算装置,首先获取微动传感器采集到的心冲击信号,然后按照每秒转化100次~1000次的速率将心冲击信号进行数字化,得到数字化的电压信号和电压变化率信号,然后通过电压信号和电压信号变化率与微动传感器的弹性势能和动能的线性关系,得到心冲击信号所具有的能量,形成能量波,然后对能量波进行过滤处理,该能量根据数字化转化率的大小,可能形成锯齿波或三角波等可以清晰得到能量波的多个波峰值和波谷值,然后将工频干扰、电器辐射等信号过滤掉,从而可以根据波峰或波谷计算心冲击信号的周期,比如,波峰可以作为一个周期的开始,上一个周期的结束等,计算简单准确。
以上所述仅为本发明的优选实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (10)
1.一种心冲击信号周期的计算方法,其特征在于,包括:
获取心冲击微动传感器产生的心冲击信号,并将心冲击信号能量化,形成能量波;
计算所述能量波变化的周期,其中,能量波变化的周期与所述心冲击信号周期相同。
2.根据权利要求1所述的心冲击信号周期的计算方法,其特征在于,所述将心冲击信号能量化,形成能量波的步骤包括:
按照预设的转化速率将心冲击信号数字化,获取对应的电压信号和电压变化率信号;
根据电压信号与电压变化率信号计算微动传感器的动能和弹性势能,其中微动传感器的弹性势能和动能,与电压信号和电压变化率信号的平方成线性关系。
3.根据权利要求2所述的根据权利要求1所述的心冲击信号周期的计算方法,其特征在于,所述预设的转化速率为每秒转化100次~1000次。
4.根据权利要求1所述的根据权利要求1所述的心冲击信号周期的计算方法,其特征在于,所述计算所述能量波变化的周期,其中,能量波变化的周期与所述心冲击信号周期相同的步骤之前,包括:
按照指定条件进行滤波处理。
5.根据权利要求1所述的根据权利要求1所述的心冲击信号周期的计算方法,所述计算所述能量波变化的周期,其中,能量波变化的周期与所述心冲击信号周期相同的步骤包括:
根据能量波相邻的两个波峰或两个波谷计算能量波变化的周期。
6.一种心冲击信号周期的计算装置,其特征在于,包括:
能量转化单元,用于获取心冲击微动传感器产生的心冲击信号,并将心冲击信号能量化,形成能量波;
计算单元,用于计算所述能量波变化的周期,其中,能量波变化的周期与所述心冲击信号周期相同。
7.根据权利要求6所述的心冲击信号周期的计算装置,其特征在于,所述能量转化单元包括:
数字转化模块,用于按照预设的转化速率将心冲击信号数字化,获取对应的电压信号和电压变化率信号;
计算模块,用于根据电压信号与电压变化率信号计算微动传感器的动能和弹性势能,其中微动传感器的弹性势能和动能,与电压信号和电压变化率信号的平方成线性关系。
8.根据权利要求7所述的心冲击信号周期的计算装置,其特征在于,所述预设的转化速率为每秒转化100次~1000次。
9.根据权利要求6所述的心冲击信号周期的计算装置,其特征在于,还包括滤波单元,用于按照指定条件进行滤波处理。
10.根据权利要求6所述的心冲击信号周期的计算装置,其特征在于,所述计算单元包括:
极值计算模块,用于根据能量波相邻的两个波峰或两个波谷计算能量波变化的周期。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510779087.1A CN105447306A (zh) | 2015-11-12 | 2015-11-12 | 心冲击信号周期的计算方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510779087.1A CN105447306A (zh) | 2015-11-12 | 2015-11-12 | 心冲击信号周期的计算方法和装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105447306A true CN105447306A (zh) | 2016-03-30 |
Family
ID=55557476
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510779087.1A Pending CN105447306A (zh) | 2015-11-12 | 2015-11-12 | 心冲击信号周期的计算方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105447306A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108836299A (zh) * | 2018-04-23 | 2018-11-20 | 深圳市友宏科技有限公司 | 一种bcg心率提取方法、存储介质以及装置 |
CN110115574A (zh) * | 2018-02-07 | 2019-08-13 | 普天信息技术有限公司 | 心率监测的方法和装置 |
CN114246581A (zh) * | 2021-12-20 | 2022-03-29 | 杭州慧光健康科技有限公司 | 基于心冲击信号短时能量的床垫传感心率识别系统及方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101689219A (zh) * | 2006-11-01 | 2010-03-31 | 必安康医疗有限公司 | 用于监测心肺参数的系统和方法 |
WO2010067297A1 (en) * | 2008-12-11 | 2010-06-17 | Koninklijke Philips Electronics N.V. | Method and apparatus for the analysis of ballistocardiogram signals |
CN104523266A (zh) * | 2015-01-07 | 2015-04-22 | 河北大学 | 一种心电信号自动分类方法 |
-
2015
- 2015-11-12 CN CN201510779087.1A patent/CN105447306A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101689219A (zh) * | 2006-11-01 | 2010-03-31 | 必安康医疗有限公司 | 用于监测心肺参数的系统和方法 |
WO2010067297A1 (en) * | 2008-12-11 | 2010-06-17 | Koninklijke Philips Electronics N.V. | Method and apparatus for the analysis of ballistocardiogram signals |
CN104523266A (zh) * | 2015-01-07 | 2015-04-22 | 河北大学 | 一种心电信号自动分类方法 |
Non-Patent Citations (3)
Title |
---|
JAKUB PARAK: "Heart Rate Detection from Ballistocardiogram", 《POSTER 2012 COMPUTER SCIENCE》 * |
JAKUB PARAK: "Heart Rate Detection from Ballistocardiogram", 《POSTER 2012》 * |
王春武等: "基于心冲击信号的心率提取算法", 《东北大学学报(自然科学版)》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110115574A (zh) * | 2018-02-07 | 2019-08-13 | 普天信息技术有限公司 | 心率监测的方法和装置 |
CN108836299A (zh) * | 2018-04-23 | 2018-11-20 | 深圳市友宏科技有限公司 | 一种bcg心率提取方法、存储介质以及装置 |
CN108836299B (zh) * | 2018-04-23 | 2021-05-14 | 深圳市友宏科技有限公司 | 一种bcg心率提取方法、存储介质以及装置 |
CN114246581A (zh) * | 2021-12-20 | 2022-03-29 | 杭州慧光健康科技有限公司 | 基于心冲击信号短时能量的床垫传感心率识别系统及方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103727959B (zh) | 计步方法及装置 | |
CN103954295B (zh) | 一种基于加速度传感器的计步方法 | |
CN105447306A (zh) | 心冲击信号周期的计算方法和装置 | |
CN104406604B (zh) | 一种计步方法 | |
CN102670207B (zh) | 一种基于足底压力的步态分析方法 | |
CN105509763B (zh) | 计步干扰去除方法以及装置 | |
CN101548885B (zh) | 去除电生理信号中工频干扰信号的方法 | |
CN105180959B (zh) | 一种适用于腕式计步器的抗干扰计步方法 | |
WO2013179018A1 (en) | Respiration rate extraction from cardiac signals | |
CN102626307A (zh) | 动态脉搏信号实时检测系统及检测方法 | |
RU2015125907A (ru) | Способ и устройство для идентификации переходов между положениями сидя и стоя | |
CA2807517C (en) | Method for measuring heart rate variability | |
JP2014510545A5 (zh) | ||
CN105651303A (zh) | 一种基于三轴加速度传感器的计步系统及方法 | |
CN106037743A (zh) | 一种提取呼吸率的方法及设备 | |
JP4952162B2 (ja) | データ処理装置,データ処理方法及びデータ処理プログラム | |
RU2015154989A (ru) | Устройство для воздействия на кожу с адаптивной обратной связью по движению | |
CN103211603A (zh) | 一种心理压力检测跟踪反馈系统 | |
CN102018507A (zh) | 一种基于示波法测量血压的数据分析处理方法 | |
CN103810393A (zh) | 一种基于心电向量的心电波形特征点定位方法 | |
CN105997043A (zh) | 一种基于腕式可穿戴设备的脉率提取方法 | |
CN106175695A (zh) | 一种睡眠呼吸暂停综合征的检测系统 | |
CN104354698A (zh) | 基于心电信号及地理信息的驾驶员应激判别系统和方法 | |
CN104545863A (zh) | 基于模糊模式识别的bcg心率提取方法及系统 | |
CN108042107A (zh) | 一种ppg信号伪差校正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160330 |
|
RJ01 | Rejection of invention patent application after publication |