CN113176451A - 一种基于查表法的脉冲信号检测方法 - Google Patents

一种基于查表法的脉冲信号检测方法 Download PDF

Info

Publication number
CN113176451A
CN113176451A CN202110436780.4A CN202110436780A CN113176451A CN 113176451 A CN113176451 A CN 113176451A CN 202110436780 A CN202110436780 A CN 202110436780A CN 113176451 A CN113176451 A CN 113176451A
Authority
CN
China
Prior art keywords
value
pulse
pulse signal
code
preset
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
CN202110436780.4A
Other languages
English (en)
Other versions
CN113176451B (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.)
Zhejiang Viewshine Intelligent Meter Co Ltd
Original Assignee
Zhejiang Viewshine Intelligent Meter 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 Zhejiang Viewshine Intelligent Meter Co Ltd filed Critical Zhejiang Viewshine Intelligent Meter Co Ltd
Priority to CN202110436780.4A priority Critical patent/CN113176451B/zh
Publication of CN113176451A publication Critical patent/CN113176451A/zh
Application granted granted Critical
Publication of CN113176451B publication Critical patent/CN113176451B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/02Measuring characteristics of individual pulses, e.g. deviation from pulse flatness, rise time or duration
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/02Measuring characteristics of individual pulses, e.g. deviation from pulse flatness, rise time or duration
    • G01R29/023Measuring pulse width

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Analogue/Digital Conversion (AREA)

Abstract

本发明公开了一种基于查表法的脉冲信号检测方法,包括:S1、基于要采集的脉冲信号制作编码表;S2、获取预设宽度的窗口内的脉冲编码值;S3、若编码表中查找到获取的脉冲编码值,则对应的脉冲信号有效;采用预置编码表和检测的实时脉冲信号的编码值进行比对,减少了脉冲信号有效性判断的复杂性和运算量,并可以有效的滤除因干扰引起的编码值,在获取脉冲信号编码值过程中只用到了移位操作,无需复杂的逻辑运算,还可以根据脉冲周期和所需检测精度调整采集脉冲编码值的窗口宽度。

Description

一种基于查表法的脉冲信号检测方法
技术领域
本发明涉及信号处理技术领域,尤其涉及一种基于查表法的脉冲信号检测方法。
背景技术
脉冲信号的检测通常需要捕捉脉冲的起始和结束、还需要计算脉冲的周期和占空比数值,才能确定波形、滤除干扰。采用上述方法运行时细节步骤多、逻辑较为复杂,计算量大,后期修改维护不方便。
发明内容
本发明为了克服以上技术的不足,提供了一种基于查表法的脉冲信号检测方法,通过预置编码表来对脉冲信号编码进行比对,减少了脉冲信号有效性判断的复杂性和运算量。
本发明提供了一种基于查表法的脉冲信号检测方法,包括S1、基于要采集的脉冲信号制作编码表;S2、获取预设宽度的窗口内的脉冲编码值;S3、若编码表中查找到获取的脉冲编码值,则对应的脉冲信号有效。
进一步的,步骤S1具体包括:S11,基于要采集的脉冲信号预设周期为T,波形长度为n×T的脉冲波形,其中n大于等于2;S12,基于预设宽度的窗口和预设采样间隔对脉冲波形取编码值;S13,基于预设滑动步进值滑动窗口,每滑动一次取一次新编码值;S14,直至窗口的最低位滑动至脉冲波形长度的第n个周期结束位置,完成所有的编码值的取值,从而制作编码表。
取出所有可能的编码值后,用这些编码值来组成编码表,不同的脉冲波形对应的编码表不同,有效的滤除干扰。
进一步的,步骤S12还包括:若编码值中任意一个bit的采样位置位于预设脉冲波形的上升沿或下降沿,则对应bit的值为0或1。
进一步的,步骤S2具体包括:S21、基于预设宽度和编码表预设置编码值;S22、基于预设采样间隔采集实时脉冲信号,并得到一个电平值;S23,丢弃预设值编码值的最高位,将其余的编码值左移一位,将检测的电平值放入最低位,得到新的编码值。
编码值每间隔一个预设采样间隔更新一次。
进一步的,步骤S3具体包括:在编码表中查找新的编码值,若编码表中查找到新的编码值,则认为新的编码值有效,对应的脉冲为有效脉冲。
进一步的,还包括S4,重复执行步骤S22-S23,得到新的编码值后执行步骤S3。
每次得到一个新的编码值都去编码表中进行查找。
进一步的,所述预设宽度基于脉冲周期和检测精度确定。
进一步的,所述预设宽度为10s。
脉冲周期越长或要求的精度越高、窗口的宽度越大,窗口宽度在选定后就不再更改。
本发明的有益效果是:
1、采用预置编码表和检测的实时脉冲信号的编码值进行比对,减少了脉冲信号有效性判断的复杂性和运算量;
2、因干扰引起的编码值不会出现在预置的编码表中,可有效的滤除干扰;
3、在获取脉冲信号编码值过程中只用到了移位操作,无需复杂的逻辑运算;4、根据脉冲周期和所需检测精度调整采集脉冲编码值的窗口宽度。
附图说明
图1为本发明实施例方法流程示意图;
图2为本发明实施例中编码表中的第1个编码值的获取;
图3为本发明实施例中编码表中的第2个编码值的获取。
具体实施方式
为了进一步理解本发明,下面结合实施例对本发明优选实施方案进行描述,但是应当理解,这些描述只是为进一步说明本发明的特征和优点,而不是对本发明权利要求的限制。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,为本发明实施例的流程示意图,本实施例应用的系统包括一个智能燃气表和一个119报警器组成,其中、119报警器输出的信号有两种:平时输出低电平、在检测到燃气泄漏报警时输出脉冲信号,其中输出的脉冲信号周期5秒、占空比50%,该信号通过线缆与智能燃气表相连接,智能燃气表通过检测线缆上的脉冲信号来判断是否发生燃气泄漏报警。
下面结合附图1和具体实施例对本发明做进一步详细说明,本实施例所述的一种基于查表法的脉冲信号检测方法具体包括如下:
S1、基于要采集的脉冲信号制作编码表。首先根据119报警器输出脉冲的波形制作编码表。
S11,基于要采集的脉冲信号预设周期为T,波形长度为n×T的脉冲波形;
根据119报警器输出的脉冲波形模拟相同的脉冲波形,从而做出预设周期为5s,占空比为50%,波形长度为3个周期的脉冲波形。其中,脉冲波形的波形长度至少为2个周期。
S12,基于预设窗口宽度和预设采样间隔对脉冲波形取编码值。
使用窗口宽度为10s的窗口对脉冲波形进行滑动取值,预设窗口宽度为10s的窗口对预设脉冲波形采样,预设采样间隔为1s。
在每1秒的间隔上取1个值,高电平为1、低电平为0,一共10个值,组成一个10bit的编码值。如图2中可取出第1个编码值为0bX1100X1100、其中最高位的X对应第10位和第5位采样,采集脉冲波形的上升沿,所以X对应的可选值为0或1。因此如图2所示,0b0110001100,0b0110011100,0b1110001100,0b1110011100四个值都用来制作编码表。
S13,基于预设滑动步进值t1滑动窗口,每滑动一次取一次新编码值;
预设滑动取值步进为1s,将窗口向右移动1s宽度。如图3所示,取第二个编码值为0b1100X1100X,对应的0b1100011000,0b1100011001,0b1100111000,0b1100111001四个值用来制作编码表。
S14,直至窗口的最低位滑动至脉冲波形长度的第3个周期结束位置,完成所有的编码值的取值,从而制作编码表。
S2、获取宽度t1的窗口内的脉冲编码值。具体包括:
S21、基于预设宽度和编码表预设置编码值。
对应预设宽度设置一个宽度为10bit的编码值,在本发明的一个实施例中。预设编码值为0b0110001100。
S22、基于预设采样间隔采集实时脉冲信号,并得到一个电平值。
每1s检测1次实时脉冲信号的高低电平、得到1个电平值,其中高电平为1、低电平为0。
S23,丢弃预设值编码值的最高位,将其余的编码值左移一位,将检测的电平值放入最低位,得到新的编码值。
先丢弃编码值的bit9,再将剩余9bit左移1位、以空出bit0,最后将得到的电平值放入bit0,得到新的编码值。
S3,若编码表中查找到获取的脉冲编码值,则对应的脉冲信号有效。
将得到的新的编码值与编码表中的值进行对照,若该编码值在编码表中,则认为该编码是报警器发出的有效脉冲。
S4、重复执行步骤S22-S23,得到新的编码值后执行步骤S3。
通过以上步骤,每1秒检测1次脉冲信号,编码值也就会每1秒更新1次。每次得到1个新的编码值之后,将得到的编码值与编码表中的值进行对照,若该编码值在编码表中,则认为该编码是报警器发出的有效脉冲。因干扰引起的编码值不会出现在预置的编码表中,可有效的滤除干扰。编码值判定为有效脉冲后、燃气表判定为报警器报警、其所处环境发生了燃气泄漏事件,同时燃气表向数据中心发送燃气泄漏报警信息。
以上仅此外,描述了本发明的基本原理和优选实施方式,本领域人员可以根据上述描述做出许多变化和改进,这些变化和改进应该属于本发明的保护范围。

Claims (8)

1.一种基于查表法的脉冲信号检测方法,其特征在于,包括:
S1、基于要采集的脉冲信号制作编码表;
S2、获取预设宽度的窗口内的脉冲编码值;
S3、若编码表中查找到获取的脉冲编码值,则对应的脉冲信号有效。
2.根据权利要求1所述的一种基于查表法的脉冲信号检测方法,其特征在于,步骤S1具体包括:
S11,基于要采集的脉冲信号预设周期为T,波形长度为n×T的脉冲波形,其中n大于等于2;
S12,基于预设宽度的窗口和预设采样间隔对脉冲波形取编码值;
S13,基于预设滑动步进值滑动窗口,每滑动一次取一次新编码值;
S14,直至窗口的最低位滑动至脉冲波形长度的第n个周期结束位置,完成所有的编码值的取值,从而制作编码表。
3.根据权利要求2所述的一种基于查表法的脉冲信号检测方法,其特征在于,步骤S12还包括:若编码值中任意一个bit的采样位置位于预设脉冲波形的上升沿或下降沿,则对应bit的值为0或1。
4.根据权利要求1所述的一种基于查表法的脉冲信号检测方法,其特征在于,步骤S2具体包括:
S21、基于预设宽度和编码表预设置编码值;
S22、基于预设采样间隔采集实时脉冲信号,并得到一个电平值;
S23,丢弃预设值编码值的最高位,将其余的编码值左移一位,将检测的电平值放入最低位,得到新的编码值。
5.根据权利要求4所述的一种基于查表法的脉冲信号检测方法,其特征在于,步骤S3具体包括:在编码表中查找新的编码值,若编码表中查找到新的编码值,则认为新的编码值有效,对应的脉冲为有效脉冲。
6.根据权利要求5所述的一种基于查表法的脉冲信号检测方法,其特征在于,还包括S4,重复执行步骤S22-S23,得到新的编码值后执行步骤S3。
7.根据权利要求1-6任一所述的一种基于查表法的脉冲信号检测方法,其特征在于,所述预设宽度基于脉冲周期和检测精度确定。
8.根据权利要求7所述的一种基于查表法的脉冲信号检测方法,其特征在于,所述预设宽度为10s。
CN202110436780.4A 2021-04-22 2021-04-22 一种基于查表法的脉冲信号检测方法 Active CN113176451B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110436780.4A CN113176451B (zh) 2021-04-22 2021-04-22 一种基于查表法的脉冲信号检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110436780.4A CN113176451B (zh) 2021-04-22 2021-04-22 一种基于查表法的脉冲信号检测方法

Publications (2)

Publication Number Publication Date
CN113176451A true CN113176451A (zh) 2021-07-27
CN113176451B CN113176451B (zh) 2023-08-15

Family

ID=76924739

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110436780.4A Active CN113176451B (zh) 2021-04-22 2021-04-22 一种基于查表法的脉冲信号检测方法

Country Status (1)

Country Link
CN (1) CN113176451B (zh)

Citations (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB849891A (en) * 1958-04-15 1960-09-28 Standard Telephones Cables Ltd Improvements in or relating to electric pulse code modulation systems of communication
JPH0355190U (zh) * 1989-09-29 1991-05-28
CN1062986A (zh) * 1990-12-28 1992-07-22 黑龙江省公安厅 无线治安防盗微机报警系统的编码方法及中继方式
US5818878A (en) * 1993-12-31 1998-10-06 Nec Corporation Method and apparatus for decoding a diphase-coded digital signal
CN1369088A (zh) * 1999-08-16 2002-09-11 关卡系统股份有限公司 采用可变时移的电子物品保安系统
CN1657740A (zh) * 2004-02-16 2005-08-24 中国石油勘探开发研究院钻井工艺研究所 一种接收和检测泥浆压力脉冲信号的方法及装置
US20060013086A1 (en) * 2004-07-16 2006-01-19 Ali Corporation Method for generating write pulse in optical disk drives
CN101321015A (zh) * 2008-06-20 2008-12-10 中国科学院上海光学精密机械研究所 全数字化的脉冲位置调制信号接收机
JP2008309756A (ja) * 2007-06-18 2008-12-25 Kawasaki Microelectronics Kk パルス幅測定方法および回路
CN101499880A (zh) * 2009-03-03 2009-08-05 四川九洲电器集团有限责任公司 基于广播式自动相关监视的信号检测方法及装置
CN102054186A (zh) * 2010-12-17 2011-05-11 深圳市远望谷信息技术股份有限公司 电子标签的信息编码方法及其装置
CN103424566A (zh) * 2013-08-20 2013-12-04 杭州日鼎控制技术有限公司 一种电机转速的检测装置及其检测方法
CN104166356A (zh) * 2014-07-30 2014-11-26 廊坊市新思维科技有限公司 一种管道爬行器的控制编码方法
CN104579352A (zh) * 2015-02-12 2015-04-29 中国科学技术大学 基于fpga的温度计码到二进制码的编码转换装置和方法
CN104597748A (zh) * 2015-02-12 2015-05-06 中国科学技术大学 一种基于fpga的时间数字变换器
CN104635126A (zh) * 2015-01-27 2015-05-20 国家电网公司 一种基于滑动窗口的局部放电单脉冲提取方法
CN104935311A (zh) * 2014-11-19 2015-09-23 中国科学院自动化研究所 一种数字信号隔离器及相应的脉宽编解码方法
CN105305755A (zh) * 2015-11-06 2016-02-03 中国科学院电工研究所 一种永磁同步电动机转子位置的检测方法
CN107153174A (zh) * 2017-05-18 2017-09-12 陕西长岭电子科技有限责任公司 Gjb661‑89标准信标机信号的识别方法
CN107231192A (zh) * 2017-02-28 2017-10-03 华南师范大学 双电平脉冲宽度编码的数据传输方法及led可见光通信系统
CN107247962A (zh) * 2017-05-23 2017-10-13 华中科技大学 一种基于滑动窗口的实时电器识别方法和系统
US20180291594A1 (en) * 2014-08-14 2018-10-11 Reliance Worldwide Corporation Methods and apparatus for fluid flow monitoring and leak detection
CN109407500A (zh) * 2018-11-22 2019-03-01 深圳天眼激光科技有限公司 一种基于fpga的时间间隔测量方法
CN109751043A (zh) * 2017-11-01 2019-05-14 中国石油化工股份有限公司 用于地层压力随钻测量工具的压力脉冲编解码系统和方法
CN111103490A (zh) * 2019-12-26 2020-05-05 广东电网有限责任公司电力科学研究院 一种负荷辨识事件检测方法、装置及设备
CN111161097A (zh) * 2019-12-31 2020-05-15 四川大学 基于假设检验的事件检测算法检测开关事件的方法和装置
CN111693785A (zh) * 2020-05-14 2020-09-22 湖南毂梁微电子有限公司 一种数字脉冲信号宽度测量电路及测量方法
CN112230178A (zh) * 2020-05-06 2021-01-15 青岛鼎信通讯股份有限公司 一种电能表需量计量方法

Patent Citations (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB849891A (en) * 1958-04-15 1960-09-28 Standard Telephones Cables Ltd Improvements in or relating to electric pulse code modulation systems of communication
JPH0355190U (zh) * 1989-09-29 1991-05-28
CN1062986A (zh) * 1990-12-28 1992-07-22 黑龙江省公安厅 无线治安防盗微机报警系统的编码方法及中继方式
US5818878A (en) * 1993-12-31 1998-10-06 Nec Corporation Method and apparatus for decoding a diphase-coded digital signal
CN1369088A (zh) * 1999-08-16 2002-09-11 关卡系统股份有限公司 采用可变时移的电子物品保安系统
CN1657740A (zh) * 2004-02-16 2005-08-24 中国石油勘探开发研究院钻井工艺研究所 一种接收和检测泥浆压力脉冲信号的方法及装置
US20060013086A1 (en) * 2004-07-16 2006-01-19 Ali Corporation Method for generating write pulse in optical disk drives
JP2008309756A (ja) * 2007-06-18 2008-12-25 Kawasaki Microelectronics Kk パルス幅測定方法および回路
CN101321015A (zh) * 2008-06-20 2008-12-10 中国科学院上海光学精密机械研究所 全数字化的脉冲位置调制信号接收机
CN101499880A (zh) * 2009-03-03 2009-08-05 四川九洲电器集团有限责任公司 基于广播式自动相关监视的信号检测方法及装置
CN102054186A (zh) * 2010-12-17 2011-05-11 深圳市远望谷信息技术股份有限公司 电子标签的信息编码方法及其装置
CN103424566A (zh) * 2013-08-20 2013-12-04 杭州日鼎控制技术有限公司 一种电机转速的检测装置及其检测方法
CN104166356A (zh) * 2014-07-30 2014-11-26 廊坊市新思维科技有限公司 一种管道爬行器的控制编码方法
US20180291594A1 (en) * 2014-08-14 2018-10-11 Reliance Worldwide Corporation Methods and apparatus for fluid flow monitoring and leak detection
CN104935311A (zh) * 2014-11-19 2015-09-23 中国科学院自动化研究所 一种数字信号隔离器及相应的脉宽编解码方法
CN104635126A (zh) * 2015-01-27 2015-05-20 国家电网公司 一种基于滑动窗口的局部放电单脉冲提取方法
CN104579352A (zh) * 2015-02-12 2015-04-29 中国科学技术大学 基于fpga的温度计码到二进制码的编码转换装置和方法
CN104597748A (zh) * 2015-02-12 2015-05-06 中国科学技术大学 一种基于fpga的时间数字变换器
CN105305755A (zh) * 2015-11-06 2016-02-03 中国科学院电工研究所 一种永磁同步电动机转子位置的检测方法
CN107231192A (zh) * 2017-02-28 2017-10-03 华南师范大学 双电平脉冲宽度编码的数据传输方法及led可见光通信系统
CN107153174A (zh) * 2017-05-18 2017-09-12 陕西长岭电子科技有限责任公司 Gjb661‑89标准信标机信号的识别方法
CN107247962A (zh) * 2017-05-23 2017-10-13 华中科技大学 一种基于滑动窗口的实时电器识别方法和系统
CN109751043A (zh) * 2017-11-01 2019-05-14 中国石油化工股份有限公司 用于地层压力随钻测量工具的压力脉冲编解码系统和方法
CN109407500A (zh) * 2018-11-22 2019-03-01 深圳天眼激光科技有限公司 一种基于fpga的时间间隔测量方法
CN111103490A (zh) * 2019-12-26 2020-05-05 广东电网有限责任公司电力科学研究院 一种负荷辨识事件检测方法、装置及设备
CN111161097A (zh) * 2019-12-31 2020-05-15 四川大学 基于假设检验的事件检测算法检测开关事件的方法和装置
CN112230178A (zh) * 2020-05-06 2021-01-15 青岛鼎信通讯股份有限公司 一种电能表需量计量方法
CN111693785A (zh) * 2020-05-14 2020-09-22 湖南毂梁微电子有限公司 一种数字脉冲信号宽度测量电路及测量方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ZHENYUAN JI等: "Identification of radio signal pulse modulation characteristics", 2012 IEEE 11TH INTERNATIONAL CONFERENCE ON SIGNAL PROCESSING, pages 1910 - 1913 *
刘昌华 主编: "《数字逻辑原理与EPGA设计 第3版》", 28 February 2021, 北京:北京航空航天大学出版社, pages: 1 - 2 *
唐勇,裴先登: "光记录的调制、编码与再生信号检测", 计算机应用研究, no. 06, pages 8 - 10 *
马芮: "二总线火灾报警系统编码与通信的实现", 河北工业大学学报, vol. 46, no. 3, pages 18 - 22 *

Also Published As

Publication number Publication date
CN113176451B (zh) 2023-08-15

Similar Documents

Publication Publication Date Title
CN1314215C (zh) 基于频谱分析的串行数据抖动测量方法
US7492846B2 (en) Standard wave receiver and time code decoding method
EP1795870B1 (de) Schaltungsanordnung und Verfahren zur Feststellung von Zählfehlern in einer Positionsmesseinrichtung
CN1118965C (zh) 数字电信站的信道估算装置
US7957464B2 (en) Digital signal receiver
CN101741626B (zh) 一种告警信息的处理方法和装置
CN101938711A (zh) 用于垃圾消息检测的方法和设备
CN109412703B (zh) 利用短脉冲短时能量谱定时的时延差编码方法
DE69729016T2 (de) Verfahren und Einrichtung zur Taktrückgewinnung
CN109088621B (zh) 一种信号滤波方法及装置
CN112559417A (zh) 空调通信的控制方法、装置、通信系统和可读存储介质
CN113176451B (zh) 一种基于查表法的脉冲信号检测方法
CN111786865A (zh) 一种数据处理方法及设备
CN109818713B (zh) 副载波调制的米勒码解码方法、装置和设备及存储介质
CN104243095A (zh) 一种卷积码与线性分组码的码字类型盲识别方法
CN201426113Y (zh) 用于智能家居电子设备的数字信号滤波整形电路
CN113375503B (zh) 激光制导抗有源干扰的方法和系统
CN108734944A (zh) 一种无线遥控解码方法及其装置
CN114253879A (zh) 一种基于软件模拟的多通道串口通信方法及系统
CN111262652B (zh) 一种基于边沿检测的fm0编码进行解码的方法及系统
CN1625833A (zh) 输入突发信号中包含的加性直流分量检测
CN112950921A (zh) 一种空调红外遥控信号匹配方法
CN109412901B (zh) 基于时域处理的采集数据连续性检测方法及检测系统
DE4418622C2 (de) Verfahren und Anordnung zur Bestimmung der Übertragungsrate in einem Bussystem
CN106849955B (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