CN106510716A - 一种基于心电信号计算呼吸率的方法及装置 - Google Patents
一种基于心电信号计算呼吸率的方法及装置 Download PDFInfo
- Publication number
- CN106510716A CN106510716A CN201610976312.5A CN201610976312A CN106510716A CN 106510716 A CN106510716 A CN 106510716A CN 201610976312 A CN201610976312 A CN 201610976312A CN 106510716 A CN106510716 A CN 106510716A
- Authority
- CN
- China
- Prior art keywords
- signal
- consequent
- parameter
- sef
- adapting filter
- 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
- 238000000034 method Methods 0.000 title claims abstract description 69
- 230000029058 respiratory gaseous exchange Effects 0.000 title claims abstract description 47
- 238000012549 training Methods 0.000 claims abstract description 20
- 238000012360 testing method Methods 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 11
- 210000002569 neuron Anatomy 0.000 claims description 21
- 230000005284 excitation Effects 0.000 claims description 15
- 238000001914 filtration Methods 0.000 claims description 10
- 239000013598 vector Substances 0.000 claims description 9
- 239000000284 extract Substances 0.000 claims description 6
- 238000013144 data compression Methods 0.000 claims description 4
- 230000017105 transposition Effects 0.000 claims description 3
- 230000000241 respiratory effect Effects 0.000 abstract description 14
- 230000036391 respiratory frequency Effects 0.000 abstract description 6
- 230000003044 adaptive effect Effects 0.000 abstract description 4
- 238000007781 pre-processing Methods 0.000 abstract 1
- 230000036387 respiratory rate Effects 0.000 abstract 1
- 230000006978 adaptation Effects 0.000 description 11
- 238000013528 artificial neural network Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 230000003190 augmentative effect Effects 0.000 description 2
- 238000000205 computational method Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 238000000718 qrs complex Methods 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
- A61B5/0816—Measuring devices for examining respiratory frequency
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- 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/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- 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/7235—Details of waveform analysis
-
- 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/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Physiology (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Pulmonology (AREA)
- Cardiology (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明提供了一种基于心电信号计算呼吸率的方法和装置,该方法包括步骤:获取心电信号并进行预处理,从所述心电信号中提取参数信号,参数信号包括训练用参数信号和测试用参数信号;自适应滤波器对训练用参数信号进行处理得到预估的呼吸信号;计算所述预估的呼吸信号与阻抗法获得的呼吸信号的差值,并调整所述自适应滤波器的参数,使所述差值在预设的范围内;通过调整参数后的自适应滤波器对所述测试用参数信号进行处理得到最终的呼吸信号,并计算呼吸率。本发明通过自适应滤波器的自适应调整,得到满足呼吸信号的纯净度要求的自适应滤波器的模型,通过调整后的自适应滤波器,将心电信号上的呼吸频率信号解调出来,实现对受试者呼吸行为的监测。
Description
技术领域
本发明涉及心电信号处理领域,特别是涉及一种基于心电信号计算呼吸率的方法及装置。
背景技术
目前用于计算呼吸率的方法主要有:阻抗容积法,用高频恒流源测量胸部阻抗的变化来提取呼吸信息;传感器法,使用温度、压力、湿度和气流传感器作为鼻孔传感器;电容法,当呼吸时导致电容值产生相应的变化;呼吸音法,通过拾取呼吸音识别呼吸;超声法,利用超声波产生多谱勒现象,检测出呼吸频率。使用这些方法不但需要增加信号采集部件,而且受到运动和环境的影晌,不适合用于日常监护。
大量临床资料显示,呼吸运动会引起心电图的变化。通过心电图,我们可以观察到在呼吸周期内由胸部运动和心脏位置变化所引起的心电波形峰峰值的改变。这是由于呼吸周期内,描述心脏电波主要传播方向的心脏电轴旋转造成QRS波群形态发生了变化。从心电信号中提取呼吸信号(ECG-DerivedRespiration,EDR)的方法不需要专用传感器和硬件模块检测呼吸信号,只需要用心电监护仪获取心电信号,避免了上述两种检测方法对人体的束缚,使动态呼吸检测成为可能。
然而现有技术中,从心电信号提取呼吸信号的技术并不成熟,仍存在一些问题。例如通过心电信号提取呼吸信号的过程中由于心电、呼吸信号的非线性、随机性和非平稳等特性,造成呼吸信号的丢失和交叉干扰的问题,使计算得到的呼吸率不准确,而且系统的鲁棒性差。
发明内容
鉴于上述状况,有必要针对现有技术中呼吸率计算不准确的问题,提供一种基于心电信号计算呼吸率的方法及装置。
本发明实施例提供了一种基于心电信号计算呼吸率的方法,包括:
获取心电信号并进行预处理,从所述心电信号中提取参数信号,所述参数信号包括训练用参数信号和测试用参数信号;
通过自适应滤波器对所述训练用参数信号进行处理得到预估的呼吸信号;
计算所述预估的呼吸信号与阻抗法获得的呼吸信号的差值,并调整所述自适应滤波器的参数,使所述差值在预设的范围内;
通过调整参数后的所述自适应滤波器对所述测试用参数信号进行处理得到最终的呼吸信号,并根据所述最终的呼吸信号计算呼吸率。
上述方法,其中,所述自适应滤波器对所述参数信号进行处理的步骤包括:
根据预设的模糊规则对所述参数信号进行模糊推理,得到所述模糊规则的第一后件;
通过正交基函数对所述参数信号进行函数扩展,得到所述模糊规则的第二后件;
根据所述第一后件和第二后件,确定所述自适应滤波器的输出信号。
上述方法,其中,所述根据预设的模糊规则对所述参数信号进行模糊推理,得到所述模糊规则的第一后件的步骤包括:
计算所述参数信号在自适应滤波器的每个神经元节点的模糊子集的隶属度;
根据所述隶属度计算每条模糊规则的激励强度;
对所述激励强度进行归一化处理,得到所述模糊规则的第一后件。
上述方法,其中,所述隶属度的计算公式为:
其中,
上述公式中,表示输入信号向量在第j个神经元节点的隶属度,l=1,2,…,r表示输入信号变量的个数,j=1,2…,n表示隶属度函数的个数,μlj(xl)表示第l个输入信号变量xl在第j个神经元节点的隶属度函数,clj表示第l个参考输入信号在第j 个神经元节点的隶属度函数的中心,表示第j个神经元节点的隶属度函数的宽度。
上述方法,其中,所述激励强度的计算公式为:
其中,为第j条模糊规则的激励强度,clj表示第l个参考输入信号在第j个神经元节点的隶属度函数的中心,表示第j个神经元节点的隶属度函数的宽度。
上述方法,其中,所述第一后件的计算公式为:
其中,为第j条模糊规则的激励强度,为第j条模糊规则的第一后件,n为模糊规则的数量。
上述方法,其中,所述正交基函数为:
所述第二后件的计算公式为:
其中,Chm(xl)为输出向量中第l个元素的第m个切比雪夫正交多项式,T为转置操作,M为所述正交基函数的个数,wj为第j条模糊规则的第二后件,α1j,α2j,...,αMj为第j条模糊规则的后件参数集,M为所述正交基函数的个数。
上述方法,其中,所述根据所述第一后件和第二后件,确定所述自适应滤波器输出的信号的步骤包括:
根据如下公式确定所述自适应滤波器输出的信号:
其中,y为所述自适应滤波器输出的信号,Oj为第j条模糊规则的第一后件,wj为第j条模糊规则的第二后件。
上述方法,其中,所述获取心电信号并进行预处理,得到参数信号的步骤包括:
获取心电信号并进行工频滤波,得到滤波信号;
提取所述滤波信号中的RR间期序列信号和RW幅值序列信号。
上述方法,其中,所述根据所述呼吸信号计算呼吸率的步骤包括:
寻找呼吸信号的波峰或者波谷;
计算所述波峰或波谷的周期并换算为呼吸率。
本发明还提供了一种基于心电信号计算呼吸率的装置,包括:
心电信号处理模块,用于获取心电信号并进行预处理,从所述心电信号中提取参数信号,所述参数信号包括训练用参数信号和测试用参数信号;
预估呼吸信号获取模块,用于通过自适应滤波器对所述训练用参数信号对应进行处理得到预估的呼吸信号;
调整模块,计算预估的呼吸信号与参考呼吸信号的差值,并调整所述自适应滤波器的参数,使所述差值在预设的范围内,得到调整后的自适应滤波器;
最终呼吸信号获取模块,用于通过调整后的所述自适应滤波器对所述测试用参数信号进行处理得到最终的呼吸信号。
呼吸率计算模块,用于寻找最终的呼吸信号的波峰或者波谷、计算所述波峰或波谷的周期并换算为呼吸率。
上述装置,其中,所述自适应滤波器包括:
模糊推理模块,用于根据预设的模糊规则对所述参数信号进行模糊推理,得到所述模糊规则的第一后件;
函数扩展模块,用于通过正交基函数对所述参数信号进行函数扩展,得到所述模糊规则的第二后件;
信号确定模块,用于根据所述第一后件和第二后件,确定所述自适应滤波器的输出信号。
本发明以阻抗法获得的呼吸信号作为参考呼吸信号,通过自适应滤波器的自适应调整,使预估的呼吸信号趋近参考呼吸信号,得到满足呼吸信号的纯净度要求的自适应滤波器的模型。将预处理后的心电信号输入调整后的自适应滤波器中,通过调整后的自适应滤波器将调制于心电信号上的呼吸频率信号解调出来,从而实现对受试者呼吸行为的监测。
附图说明
图1为本发明第一实施例提供的基于心电信号计算呼吸率的方法流程图;
图2为心电信号图;
图3为工频滤波后的心电信号图;
图4为自适应滤波器的原理框图;
图5为本发明第二实施例提供的基于心电信号计算呼吸率的方法流程图;
图6为本发明实施例提供的基于心电信号计算呼吸率的装置的结构框图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
参照下面的描述和附图,将清楚本发明的实施例的这些和其他方面。在这些描述和附图中,具体公开了本发明的实施例中的一些特定实施方式,来表示实施本发明的实施例的原理的一些方式,但是应当理解,本发明的实施例的范围不受此限制。相反,本发明的实施例包括落入所附加权利要求书的精神和内涵范围内的所有变化、修改和等同物。
请参阅图1为本发明实施例提供的一种基于心电信号计算呼吸率的方法,包括步骤S10~S13
步骤S10,获取心电信号,并对所述心电信号进行预处理,从所述心电信号中提取参数信号,所述参数信号包括训练用参数信号和测试用参数信号。本发明通过加压单极肢体导联获得原始的心电信号,如图2所示,为原始的心电信号图,由于原始心电信号包含大量的工频干扰,需要进行50Hz工频陷波,滤除工频干扰,工频滤波后的心电信号如图3所示。然后通过阈值法提取工频滤波后的RR间期序列信号和RW幅值序列信号,即为所述参数信号。
步骤S11,通过自适应滤波器对所述训练用参数信号进行处理得到预估的呼吸信号。
步骤S12,计算所述预估的呼吸信号与阻抗法获得的呼吸信号的差值,并调整所述自适应滤波器的参数,使所述差值在预设的范围内,得到调整后的自适应滤波器。
自适应滤波器的原理框图如图4所示,提取预处理的心电信号的RR间期序列信号和RW幅值序列信号n(k)输入到自适应滤波器中。通过自适应滤波器处理所述参数信号得到预估的呼吸信号y(k)。通过计算阻抗法得到的呼吸信号q(k)与预估的呼吸信号的差值,即q(k)-y(k),得到偏差e(k),所述偏差越小越好,偏差越小,表示自适应滤波器输出的呼吸信号越纯净。自适应滤波器不断地重新调整,使偏差e(k)达到预设的范围内,以满足呼吸信号的纯净度要求,从而得到调整后的自适应滤波器。通过阻抗法检测人体的阻抗可获取呼吸信号,以阻抗法获取的呼吸信号作为参考的呼吸信号,调整自适应滤波器的参数。
步骤S13,通过调整后的自适应滤波器对所述测试用参数信号进行处理得到最终的呼吸信号,并根据所述最终的呼吸信号计算呼吸率。呼吸率的计算方法为:寻找呼吸信号的波峰或者波谷;计算所述波峰或波谷的周期并将其换算为呼吸率。
本发明以阻抗法获得的呼吸信号作为参考呼吸信号,自适应滤波器通过自适应调整,使预估的呼吸信号趋近参考呼吸信号,得到满足呼吸信号的纯净度要求的自适应滤波器的模型。选取随呼吸运动发生周期性波动的RR间期序列信号和RW幅值序列信号输入调整后的自适应滤波器中,通过调整后的自适应滤波器将调制于心电信号上的呼吸频率信号解调出来,从而实现对受试者呼吸行为的监测。
请参阅图5,为本发明第二实施例提供的一种基于心电信号计算呼吸率方法的流程图。如图所示,包括步骤S20~S25。
步骤S20,获取心电信号并进行预处理,从所述心电信号中提取参数信号,所述参数信号包括训练用参数信号和测试用参数信号。该步骤中,参考信号为从心电信号中提取的RR间期序列信号和RW幅值序列信号。心电信号进行预处理提取参数信号的过程参照第一实施例。
步骤S21,根据预设的模糊规则对所述训练用参数信号进行模糊推理,得到所述模糊规则的第一后件。得到所述模糊规则的第一后件具体包括步骤S211~S213。
步骤S211,计算所述训练用参数信号在自适应滤波器的每个神经元节点的模糊子集的隶属度。
考虑到信号从产生到电极采集的过程中有一定延迟,对训练用参数信号n(k)进行延时处理后再进行模糊推理。在n(k)后接入抽头延迟线,伪迹信号经过r-1次延迟,得到r维输出向量X(k)=[x1(k),x2(k),...,xr(k)]T。该r维输出向量进入自适应滤波器的第一层,该层的每一个神经元节点为一个隶属度函数,可选取高斯函数为隶属度函数,具体如下:
所述隶属度的计算公式为:
其中,
上述公式中,表示输入信号向量在第j个神经元节点的隶属度,l=1,2,…,r表示输入信号变量的个数,j=1,2…,n表示隶属度函数的个数,μlj(xl)表示第l个输入信号变量x在第j个神经元节点的隶属度函数,clj表示第l个参考输入信号在第j个神经元节点的隶属度函数的中心,表示第j个神经元节点的隶属度函数的宽度。
步骤S212,根据所述隶属度计算每条模糊规则的激励强度。所述激励强度计算公式为:
其中,为第j条模糊规则的激励强度,clj表示第l个输入信号在第j个神经元节点的隶属度函数的中心,表示第j个神经元节点的隶属度函数的宽度。
步骤S213,对所述激励强度进行归一化处理,得到所述模糊规则的第一后件。所述第一后件为:
其中,为第j条模糊规则的激励强度,为第j条模糊规则的第一后件,n为模糊规则的数量。
步骤S22,通过正交基函数对所述训练用参数信号进行函数扩展,得到所述模糊规则的第二后件。
本实施例通过函数链神经网络(FLNN)进行参数信号的函数扩展,函数链神经网络的正交基采用切比雪夫正交多项式(ChebyshevOrthogonalPolynomials,COP):
Ch0(x)=1
Ch1(x)=x
Ch2(x)=2x2-1
……
Chm+1(x)=2xChm(x)-Chm-1(x)
FLNN的基函数T如公式:
其中,Chm(xl)为输出向量中第l个元素的第m个切比雪夫正交多项式,T为转置操作,M为所述正交基函数的个数。在实际情况下,也可以采用其他方法进行函数扩展。经过函数扩展,可将低维扩展到高维空间,实现非线性。
FLNN输出所述模糊规则的第二后件:
其中,wj为第j条模糊规则的第二后件,α1j,α2j,...,αMj为第j条模糊规则的后件参数集,M为所述正交基函数的个数。
本实施例将函数链神经网络应用到自适应滤波器中,通过一组正交基函数将原输入矢量进行维数扩展,将线性参数扩展为非线性,得到模糊规则的第二后件,以增强自适应滤波器的非线性处理能力。
步骤S23,根据所述第一后件和第二后件,确定所述自适应滤波器的输出信号,即预估的呼吸信号。根据所述第一后件和第二后件确定所述自适应滤波器输出的信号的公式为:
其中,y为所述自适应自适应滤波器输出的信号,Oj为第j条模糊规则的第一后件,wj为第j条模糊规则的第二后件。
步骤S24,计算预估的呼吸信号与阻抗法获得的呼吸信号的差值,并调整所述自适应滤波器的参数,使所述差值在预设的范围内,得到调整后的自适应滤波器。
步骤S25,通过调整参数后的自适应滤波器对所述测试用参数信号进行处理得到最终的呼吸信号,并根据所述最终的呼吸信号计算呼吸率。调整参数后的自适应滤波器对所述参数信号进行处理的过程可参考训练用参数信号处理的过程,即步骤S211~S213和步骤S22~S23。
本发明以阻抗法获得的呼吸信号作为参考呼吸信号,自适应滤波器通过自适应调整,使预估的呼吸信号趋近参考呼吸信号,得到满足呼吸信号的纯净度要求的自适应滤波器的模型。通过调整后的自适应滤波器将心电信号上的呼吸频率信号解调出来,实现对受试者呼吸行为的实时监测。本实施例融合函数链神经网络和自适应滤波器的模糊推理,从心电波中提取呼吸信号,充分考虑了心电、呼吸信号的非线性,随机性和非平稳等特性,尽可能地减少有用呼吸信号的丢失,系统的鲁棒性更好。并且通过模糊非线性处理输入自适应滤波器的心电信号,能够解决信号提取过程中的交叉干扰问题。
请参阅图6,本发明实施例还提供了一种基于心电信号计算呼吸率的装置,包括心电信号处理模块、预估呼吸信号获取模块、自适应滤波器、调整模块最终呼吸信号获取模块和呼吸率计算模块。
心电信号处理模块用于获取心电信号并进行预处理,从所述心电信号中提取参数信号,所述参数信号包括训练用参数信号和测试用参数信号。本发明通过加压单极肢体导联获得原始的心电信号,并通过50Hz工频陷波,滤除工频干扰。然后通过阈值法提取工频滤波后的心电信号的RR间期序列信号和RW幅值序列信号,即为所述参数信号。
自适应滤波器用于处理所述参数信号。该自适应滤波器具体包括:
模糊推理模块,用于根据预设的模糊规则对所述参数信号进行模糊推理,得到所述模糊规则的第一后件;
函数扩展模块,用于通过正交基函数对所述参数信号进行函数扩展,得到所述模糊规则的第二后件;
信号确定模块,用于根据所述第一后件和第二后件,确定所述自适应滤波器的输出信号。
预估呼吸信号获取模块用于通过自适应滤波器对所述训练用参数信号对应进行处理得到预估的呼吸信号;
调整模块用于计算预估的呼吸信号与参考呼吸信号的差值,并调整所述自适应滤波器的参数,使所述差值在预设的范围内,得到调整后的自适应滤波器。
最终呼吸信号获取模块,用于通过调整后的所述自适应滤波器对所述测试用参数信号进行处理得到最终的呼吸信号。
通过自适应滤波器处理所述训练用参数信号得到预估的呼吸信号y(k)。通过计算阻抗法得到的呼吸信号q(k)与预估的呼吸信号的差值,即q(k)-y(k),得到偏差e(k),所述偏差越小越好,偏差越小,表示自适应滤波器输出的呼吸信号越纯净。自适应波器不断地重新调整,使偏差e(k)达到预设的范围内,以满足呼吸信号的纯净度要求,从而得到调整后的自适应滤波器。通过阻抗法检测人体的阻抗可获取呼吸信号,以阻抗法获取的呼吸信号作为参考的呼吸信号,调整自适应滤波器的参数。
调整后的自适应滤波器对所述测试用参数信号进行处理得到最终的呼吸信号,并通过呼吸率计算模块计算呼吸率。呼吸率的计算方法为:寻找呼吸信号的波峰或者波谷;计算所述波峰或波谷的周期并将其换算为呼吸率。
本发明以阻抗法获得的呼吸信号作为参考呼吸信号,自适应滤波器通过自适应调整,使预估的呼吸信号趋近参考呼吸信号,得到满足呼吸信号的纯净度要求的自适应滤波器的模型。选取随呼吸运动发生周期性波动的RR间期序列信号和RW幅值序列信号输入调整后的自适应滤波器中,通过调整后的自适应滤波器,将调制于心电信号上的呼吸频率信号解调出来,从而实现对受试者呼吸行为的监测。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由权利要求及其等同物限定。
Claims (12)
1.一种基于心电信号计算呼吸率的方法,其特征在于,包括:
获取心电信号并进行预处理,从所述心电信号中提取参数信号,所述参数信号包括训练用参数信号和测试用参数信号;
通过自适应滤波器对所述训练用参数信号进行处理得到预估的呼吸信号;
计算所述预估的呼吸信号与阻抗法获得的呼吸信号的差值,并调整所述自适应滤波器的参数,使所述差值在预设的范围内;
通过调整参数后的所述自适应滤波器对所述测试用参数信号进行处理得到最终的呼吸信号,并根据所述最终的呼吸信号计算呼吸率。
2.如权利要求1所述的方法,其特征在于,所述自适应滤波器对所述参数信号进行处理的步骤包括:
根据预设的模糊规则对所述参数信号进行模糊推理,得到所述模糊规则的第一后件;
通过正交基函数对所述参数信号进行函数扩展,得到所述模糊规则的第二后件;
根据所述第一后件和第二后件,确定所述自适应滤波器的输出信号。
3.如权利要求2所述的方法,其特征在于,所述根据预设的模糊规则对所述参数信号进行模糊推理,得到所述模糊规则的第一后件的步骤包括:
计算所述参数信号在自适应滤波器的每个神经元节点的模糊子集的隶属度;
根据所述隶属度计算每条模糊规则的激励强度;
对所述激励强度进行归一化处理,得到所述模糊规则的第一后件。
4.如权利要求3所述的方法,其特征在于,所述隶属度的计算公式为:
其中,
上述公式中,表示输入信号向量在第j个神经元节点的隶属度,l=1,2,…,r表示输入信号变量的个数,j=1,2…,n表示隶属度函数的个数,μlj(xl)表示第l个输入信号变量xl在第j个神经元节点的隶属度函数,clj表示第l个参考输入信号在第j个神经元节点的隶属度函数的中心,表示第j个神经元节点的隶属度函数的宽度。
5.如权利要求3所述的方法,其特征在于,所述激励强度的计算公式为:
其中,为第j条模糊规则的激励强度,clj表示第l个参考输入信号在第j个神经元节点的隶属度函数的中心,表示第j个神经元节点的隶属度函数的宽度。
6.如权利要求3所述的方法,其特征在于,所述第一后件的计算公式为:
其中,为第j条模糊规则的激励强度,为第j条模糊规则的第一后件,n为模糊规则的数量。
7.如权利要求2或3所述的方法,其特征在于,所述正交基函数为:
所述第二后件的计算公式为:
其中,Chm(xl)为输出向量中第l个元素的第m个切比雪夫正交多项式,T为转置操作,M为所述正交基函数的个数,wj为第j条模糊规则的第二后件,α1j,α2j,...,αMj为第j条模糊规则的后件参数集,M为所述正交基函数的个数。
8.如权利要求2或3所述的方法,其特征在于,所述根据所述第一后件和第二后件,确定所述自适应滤波器输出的信号的步骤包括:
根据如下公式确定所述自适应滤波器输出的信号:
其中,y为所述自适应滤波器输出的信号,Oj为第j条模糊规则的第一后件,wj为第j条模糊规则的第二后件。
9.如权利要求1所述的方法,其特征在于,所述获取心电信号并进行预处理,得到参数信号的步骤包括:
获取心电信号并进行工频滤波,得到滤波信号;
提取所述滤波信号中的RR间期序列信号和RW幅值序列信号。
10.如权利要求1所述的方法,其特征在于,所述根据所述呼吸信号计算呼吸率的步骤包括:
寻找呼吸信号的波峰或者波谷;
计算所述波峰或波谷的周期并换算为呼吸率。
11.一种基于心电信号计算呼吸率的装置,其特征在于,包括:
心电信号处理模块,用于获取心电信号并进行预处理,从所述心电信号中提取参数信号,所述参数信号包括训练用参数信号和测试用参数信号;
预估呼吸信号获取模块,用于通过自适应滤波器对所述训练用参数信号对应进行处理得到预估的呼吸信号;
调整模块,计算预估的呼吸信号与参考呼吸信号的差值,并调整所述自适应滤波器的参数,使所述差值在预设的范围内,得到调整后的自适应滤波器;
最终呼吸信号获取模块,用于通过调整后的所述自适应滤波器对所述测试用参数信号进行处理得到最终的呼吸信号。
呼吸率计算模块,用于寻找最终的呼吸信号的波峰或者波谷、计算所述波峰或波谷的周期并换算为呼吸率。
12.如权利要求11所述的装置,其特征在于,所述自适应滤波器包括:
模糊推理模块,用于根据预设的模糊规则对所述参数信号进行模糊推理,得到所述模糊规则的第一后件;
函数扩展模块,用于通过正交基函数对所述参数信号进行函数扩展,得到所述模糊规则的第二后件;
信号确定模块,用于根据所述第一后件和第二后件,确定所述自适应滤波器的输出信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610976312.5A CN106510716B (zh) | 2016-11-07 | 2016-11-07 | 一种基于心电信号计算呼吸率的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610976312.5A CN106510716B (zh) | 2016-11-07 | 2016-11-07 | 一种基于心电信号计算呼吸率的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106510716A true CN106510716A (zh) | 2017-03-22 |
CN106510716B CN106510716B (zh) | 2019-09-20 |
Family
ID=58349981
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610976312.5A Active CN106510716B (zh) | 2016-11-07 | 2016-11-07 | 一种基于心电信号计算呼吸率的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106510716B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108704209A (zh) * | 2017-05-19 | 2018-10-26 | 纳智源科技(唐山)有限责任公司 | 呼吸频率监测装置、系统、呼吸机及吸氧机 |
WO2018209993A1 (zh) * | 2017-05-19 | 2018-11-22 | 纳智源科技(唐山)有限责任公司 | 呼吸频率监测装置、系统、呼吸机及吸氧机 |
CN109984742A (zh) * | 2019-04-22 | 2019-07-09 | 深圳大学 | 心阻抗信号处理系统和方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101953685A (zh) * | 2009-04-20 | 2011-01-26 | 迈瑞Ds美国有限责任公司 | 测量心电图和呼吸的装置及方法 |
US20110066037A1 (en) * | 2009-09-14 | 2011-03-17 | Matt Banet | Body-worn monitor for measuring respiration rate |
CN104398252A (zh) * | 2014-11-05 | 2015-03-11 | 深圳先进技术研究院 | 一种心电信号处理方法及装置 |
CN104720808A (zh) * | 2013-12-20 | 2015-06-24 | 中国移动通信集团公司 | 一种检测人体睡眠呼吸的方法和装置 |
CN105722457A (zh) * | 2013-11-20 | 2016-06-29 | 通用电气公司 | 用于确定呼吸率的方法和系统 |
-
2016
- 2016-11-07 CN CN201610976312.5A patent/CN106510716B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101953685A (zh) * | 2009-04-20 | 2011-01-26 | 迈瑞Ds美国有限责任公司 | 测量心电图和呼吸的装置及方法 |
US20110066037A1 (en) * | 2009-09-14 | 2011-03-17 | Matt Banet | Body-worn monitor for measuring respiration rate |
US20140276175A1 (en) * | 2009-09-14 | 2014-09-18 | Sotera Wireless, Inc. | Body-worn monitor for measuring respiration rate |
CN105722457A (zh) * | 2013-11-20 | 2016-06-29 | 通用电气公司 | 用于确定呼吸率的方法和系统 |
CN104720808A (zh) * | 2013-12-20 | 2015-06-24 | 中国移动通信集团公司 | 一种检测人体睡眠呼吸的方法和装置 |
CN104398252A (zh) * | 2014-11-05 | 2015-03-11 | 深圳先进技术研究院 | 一种心电信号处理方法及装置 |
Non-Patent Citations (3)
Title |
---|
VARANINI M: "Adaptive filtering of ECG signal for deriving respiratory activity", 《COMPUTERS IN CARDIOLOGY》 * |
席涛,杨国胜,汤池,焦腾: "基于自适应滤波的心电图中呼吸信号的提取方法", 《第四军医大学学报》 * |
席涛: "USB心电采集系统的研制及从心电图提取呼吸信号的方法探讨", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108704209A (zh) * | 2017-05-19 | 2018-10-26 | 纳智源科技(唐山)有限责任公司 | 呼吸频率监测装置、系统、呼吸机及吸氧机 |
WO2018209993A1 (zh) * | 2017-05-19 | 2018-11-22 | 纳智源科技(唐山)有限责任公司 | 呼吸频率监测装置、系统、呼吸机及吸氧机 |
CN108704209B (zh) * | 2017-05-19 | 2024-04-16 | 纳智源科技(唐山)有限责任公司 | 呼吸频率监测装置、系统、呼吸机及吸氧机 |
CN109984742A (zh) * | 2019-04-22 | 2019-07-09 | 深圳大学 | 心阻抗信号处理系统和方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106510716B (zh) | 2019-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108416367B (zh) | 基于多传感器数据决策级融合的睡眠分期方法 | |
CN110327036B (zh) | 从穿戴式心电中提取呼吸信号和呼吸频率的方法 | |
CN106691474A (zh) | 融合脑电信号与生理信号的疲劳检测系统 | |
CN108388912A (zh) | 基于多传感器特征优化算法的睡眠分期方法 | |
CN114052744B (zh) | 基于脉冲神经网络的心电信号分类方法 | |
CN109907752A (zh) | 一种去除运动伪影干扰与心电特征检测的心电诊断与监护方法及系统 | |
CN110946556A (zh) | 基于可穿戴式体感网的帕金森静息态震颤评估方法 | |
CN104545818A (zh) | 一种基于脉搏与血氧信号的睡眠呼吸暂停综合征检测方法 | |
CN106037743A (zh) | 一种提取呼吸率的方法及设备 | |
CN106510716A (zh) | 一种基于心电信号计算呼吸率的方法及装置 | |
CN106037655A (zh) | 一种睡眠周期信号的分离及周期计算方法及装置 | |
CN108681685A (zh) | 一种基于人体表面肌电信号的躯干动作意图识别方法 | |
CN111248879A (zh) | 一种基于多模态注意力融合的高血压老年人活动分析方法 | |
CN106388825A (zh) | 一种呼吸率提取方法及装置 | |
CN105662454B (zh) | 一种罗音检测的方法和装置 | |
CN107894965A (zh) | 一种用于两组不同类型信号的耦合处理方法 | |
CN106539586A (zh) | 一种呼吸率计算方法及装置 | |
CN106580324A (zh) | 一种呼吸信号提取方法及装置 | |
Ankitha et al. | Literature review on sleep APNEA analysis by machine learning algorithms using ECG signals | |
Chan et al. | Estimating heart rate from seismocardiogram signal using a novel deep dominant frequency regressor and domain adversarial training | |
CN106388824B (zh) | 一种呼吸率提取方法及装置 | |
Xu et al. | ART2 neural network for surface EMG decomposition | |
CN106725487A (zh) | 一种呼吸率提取方法及装置 | |
Zhou et al. | Embedded real-time QRS detection algorithm for pervasive cardiac care system | |
CN113208615B (zh) | 一种用于心肺复苏仪的连续脑电监测及反馈系统与方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |