CN102100553A - 一种基于rls自适应滤波器的胃电慢波信号检测方法 - Google Patents
一种基于rls自适应滤波器的胃电慢波信号检测方法 Download PDFInfo
- Publication number
- CN102100553A CN102100553A CN2009102485415A CN200910248541A CN102100553A CN 102100553 A CN102100553 A CN 102100553A CN 2009102485415 A CN2009102485415 A CN 2009102485415A CN 200910248541 A CN200910248541 A CN 200910248541A CN 102100553 A CN102100553 A CN 102100553A
- Authority
- CN
- China
- Prior art keywords
- filter
- data
- signal
- digital
- input
- 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
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开一种基于RLS自适应滤波器的胃电慢波信号检测方法,包括以下步骤:a.信号采集;b.信号放大;c.模拟信号转换成数字信号;d.信号处理;e.数据预处理;f.噪声估计;g.基于RLS自适应滤波器处理;本发明具有以下有益效果及优点:1.面向胃电慢波信号的自适应处理滤波器设计,即把心电、呼吸信号、抖动,以及信号放大装置本身的噪声的都作为胃电慢波信号的噪声,而对其进行自适应滤波器的设计,以降低滤波器的复杂程度。2.测量信号通过自适应滤波后的信号复现,即可以完成对胃电慢波信号的精确复现。
Description
技术领域
本发明涉及一种胃电信号处理的方法,具体的说是为获得可以进行医学诊断的胃电慢波信号处理方法及装置。
背景技术
胃病是一种常见、多发、难根治性疾病。随着医学诊断技术的进展,胃病诊断已经不是难事。从钡餐、内窥镜、胃液分析到活体组织检查再到C14色谱诊断,器质性胃病容易被发现。但是对于功能性胃病,却无济于事,上述方法会给病人带来一定的痛苦、损伤,并有传染疾病的危险。作为一种、安全、无创、无痛苦的胃肠动力诊断方法,即通过体表贴测量电极是消化道生理检测、诊断的无创检测方法。
胃电信号十分微弱,变化十分缓慢,因此通过体表电极的观测胃电难度很大,有另一方面,对胃电图的临床意义认识不足,目前无论在胃电检测手段、分析方法及有效性方面都需要很大的提高。
胃电信号是一种典型的微弱信号,体表胃电幅度为几百微伏到几毫伏量级,正常人胃电信号慢波成分频率极低,带宽为0.02-0.3Hz;体表采集到的信号又包含了各种干扰,如心电、呼吸信号、抖动等,部分干扰的幅度还可能远大于胃电信号。测量过程中会引入大量的强背景噪声,其中又以心电、呼吸伪迹最难去除。
自适应滤波器可以不必事先给定信号及噪声的自相关函数,它可以利用前一时刻已获得的滤波器参数自动地调节现时刻的滤波器参数使得滤波器输出和未知的输入之间的均方误差最小化或者(加权)最小二乘误差最小化,从而它可以实现最优滤波。
自适应滤波器的算法有很多,有RLS(递归最小二乘法)和LMS(最小均方算法)等。
发明内容
本发明针对胃电中的慢波信号,通过自适应信号处理方法,完成对慢波信号的观测。
本发明的目的是提供一种能够对体表测量电极获取信号的自适应滤波处理方法,去除心电、呼吸信号、抖动等带来的测量影响,完成对胃电慢波信号的观测。
为解决上述技术问题,本发明公开一种基于递归最小二乘RLS(Recursiveleast squares)自适应滤波器的胃电慢波信号检测方法,包括以下步骤:
a.信号采集;
b.信号放大;
c.模拟信号转换成数字信号;
d.信号处理;
e.数据预处理;
f.噪声估计;
g.基于RLS自适应滤波器处理;
h.数据输出。
所述步骤b中,信号放大500-1000倍;
所述步骤e包括数据重采样步骤和数字带通/低通滤波步骤,所述数据重采样步骤即对输入的数字信号依次进行重新采样,数据输入端输入为采样率Fs的数据,经过重采样运算,数据输出端输出为采样率fs的数据,其中有条件Fs>fs,并且fs为1-10Hz。
数字带通/低通滤波步骤所述即对0.01-0.5Hz数字带通或截止频率为0.5Hz数字低通滤波处理;。
所述步骤f包括数据主频查找步骤和数字带阻滤波步骤,所述数据主频查找步骤为数字信号输入端输入一串数字信号,经过功率谱分析,主峰位置查找后,从输出端输出主峰对应主频值;
所述数字带阻滤波步骤为:首先从主频输入端输入主频数值,经过阻带设置部分得到滤波器阻带,该滤波器阻带经过滤波器参数动态生成部分得到一组滤波器参数,利用这组参数动态生成一个带阻滤波器,然后,从数字信号输入端输入数字信号,经过该带阻滤波器,在数字信号输出端处输出滤波后的数字信号。
本发明具有以下有益效果及优点:
1.面向胃电慢波信号的自适应处理滤波器设计,即把心电、呼吸信号、抖动,以及信号放大装置本身的噪声的都作为胃电慢波信号的噪声,而对其进行自适应滤波器的设计,以降低滤波器的复杂程度。
2.测量信号通过自适应滤波后的信号复现,即可以完成对胃电慢波信号的精确复现。
附图说明
图1为本发明的结构框图;
图2为数据重采样模块示意图;
图3为数字带通/低通滤波模块示意图;
图4为数据主频查找模块示意图;
图5为数字带阻滤波模块示意图;
图6为胃电慢波信号复现示意图。
具体实施方式
下面结合实施例及附图详细说明本发明。
本发明的方法可由下列装置实现:
如图1所示,一种基于RLS自适应滤波器的胃电慢波信号检测装置,包括硬件部分1和软件处理部分2,所述硬件部分1包括体表测量电极、模拟信号放大器11、A/D转换卡12和数字运算设备13;所述体表电极用于信号采集,所述体表电极通过导联和模拟信号放大器11相连,所述模拟信号放大器11将体表微弱的胃电信号放大;所述模拟信号放大器11通过导线和A/D转换卡12相连,所述A/D转换卡12将模拟信号转换为数字信号;所述A/D转换卡12通过数据线和数字运算设备13相连,所述数字运算设备13将数字信号进行信号处理;
所述软件部分2包括数字预处理模块21、噪声估计模块22、RLS自适应滤波模块23及显示、存储等设备,其中所述数据预处理模块21包括数据重采样模块211和数字带通/低通滤波模块212;所述噪声估计模块22包括数据主频查找模块221和数字带阻滤波模块222;所述数据重采样模块211和数字带通/低通滤波模块212将数字运算设备13输出的数字信号进行预处理;所述的预处理结果一部分依次通过数据主频查找模块221和数字带阻滤波模块222,得到噪声估计;最后所述预处理结果以及所述噪声估计输入RLS自适应滤波模块,由RLS自适应滤波模块输出数据到显示、存储等设备。
如图2所示,所述数据重采样模块211包括数据输入端,重采样运算,数据输出端三个部分,这三个部分依次连接。数据输入端输入为采样率Fs的数据,经过重采样运算,数据输出端输出为采样率fs的数据,其中有条件Fs>fs,并且fs约为1-10Hz。
如图3所示,所述数字带通/低通滤波模块212包括数字信号输入端,数字滤波,数字信号输出端三个部分,这三个部分依次连接。数字信号输入端输入一串数字信号,经过通带约0.01-0.5Hz数字带通或者截止频率约0.5Hz数字低通滤波,数字信号输出端输出为滤波后的数字信号。
如图4所示,所述数据主频查找模块221包括数字信号输入端,功率谱分析,主峰位置查找,输出端四个部分,这四个部分依次连接。数字信号输入端输入一串数字信号,经过功率谱分析,主峰位置查找后,从输出端输出主峰对应主频值。
如图5所示,所述数字带阻滤波模块222包括主频输入端,数字信号输入端,阻带设置,滤波器参数动态生成,带阻滤波器(动态生成),数字信号输出端六个部分。其中主频输入端连接阻带设置部分,阻带设置部分连接滤波器参数动态生成部分,滤波器参数动态生成部分连接带阻滤波器部分,数字信号输入端与数字信号输出端分别连接带阻滤波器部分。首先从主频输入端输入主频数值,经过阻带设置部分得到滤波器阻带,该滤波器阻带经过滤波器参数动态生成部分得到一组滤波器参数,利用这组参数动态生成一个带阻滤波器,然后,从数字信号输入端输入数字信号,经过该带阻滤波器,在数字信号输出端处输出滤波后的数字信号。
图6为胃电慢波信号复现示意图,由此可见,本发明可以完成对胃电慢波信号的精确复现。
如上所述为本发明的一个典型的实施方式,由于实施方式较多,不能一一列举,其他任何在本发明的整体框架下进行组合、简化,替代等,都应在本发明的保护范围之内。
Claims (4)
1.一种基于RLS自适应滤波器的胃电慢波信号检测方法,其特征在于包括以下步骤:
a.信号采集;
b.信号放大;
c.模拟信号转换成数字信号;
d.信号处理;
e.数据预处理;
f.噪声估计;
g.基于RLS自适应滤波器处理;
h.数据输出。
2.按权利要求1所述的检测方法,其特征在于:所述步骤b中,信号放大500-1000倍。
3.按权利要求1所述的检测方法,其特征在于:所述步骤e包括数据重采样步骤和数字带通/低通滤波步骤,所述数据重采样步骤即对输入的数字信号依次进行重新采样,数据输入端输入为采样率Fs的数据,经过重采样运算,数据输出端输出为采样率fs的数据,其中有条件Fs>fs,并且fs为1-10Hz。
数字带通/低通滤波步骤所述即对0.01-0.5Hz数字带通或截止频率为0.5Hz数字低通滤波处理;。
4.按权利要求1所述的检测方法,其特征在于:所述步骤f包括数据主频查找步骤和数字带阻滤波步骤,所述数据主频查找步骤为数字信号输入端输入一串数字信号,经过功率谱分析,主峰位置查找后,从输出端输出主峰对应主频值;
所述数字带阻滤波步骤为:首先从主频输入端输入主频数值,经过阻带设置部分得到滤波器阻带,该滤波器阻带经过滤波器参数动态生成部分得到一组滤波器参数,利用这组参数动态生成一个带阻滤波器,然后,从数字信号输入端输入数字信号,经过该带阻滤波器,在数字信号输出端处输出滤波后的数字信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102485415A CN102100553A (zh) | 2009-12-18 | 2009-12-18 | 一种基于rls自适应滤波器的胃电慢波信号检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102485415A CN102100553A (zh) | 2009-12-18 | 2009-12-18 | 一种基于rls自适应滤波器的胃电慢波信号检测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102100553A true CN102100553A (zh) | 2011-06-22 |
Family
ID=44153780
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009102485415A Pending CN102100553A (zh) | 2009-12-18 | 2009-12-18 | 一种基于rls自适应滤波器的胃电慢波信号检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102100553A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105078442A (zh) * | 2014-05-16 | 2015-11-25 | 中国科学院沈阳自动化研究所 | 基于特征谱的胃磁慢波信号频率检测方法 |
WO2018018570A1 (zh) * | 2016-07-29 | 2018-02-01 | 华为技术有限公司 | 用于心电测量的装置和方法 |
-
2009
- 2009-12-18 CN CN2009102485415A patent/CN102100553A/zh active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105078442A (zh) * | 2014-05-16 | 2015-11-25 | 中国科学院沈阳自动化研究所 | 基于特征谱的胃磁慢波信号频率检测方法 |
WO2018018570A1 (zh) * | 2016-07-29 | 2018-02-01 | 华为技术有限公司 | 用于心电测量的装置和方法 |
CN108601544A (zh) * | 2016-07-29 | 2018-09-28 | 华为技术有限公司 | 用于心电测量的装置和方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109907752B (zh) | 一种去除运动伪影干扰与心电特征检测的心电诊断与监护系统 | |
Nayak et al. | Filtering techniques for ECG signal processing | |
Obeid et al. | Two multichannel integrated circuits for neural recording and signal processing | |
JP6980011B2 (ja) | 心臓信号をフィルタリングするためのシステム | |
CN113397555A (zh) | 生理参数监测用c-lstm的心率失常分类算法 | |
US11147517B2 (en) | Physiological measurement device with common mode interference suppression | |
CN101032396A (zh) | 基于心率变异信号的脑缺血缺氧监护系统 | |
US20070073184A1 (en) | Sensor Device for Detecting LEEG Signals and Detecting Method Thereof | |
CN203555751U (zh) | 一种基于时频分析的腹部肠鸣音分析仪系统 | |
CN104473631A (zh) | 一种基于非负盲分离胎儿心电瞬时心率识别方法及系统 | |
CN111419208A (zh) | 一种基于加速度传感器的无束缚实时心率监测方法与系统 | |
CN102008316A (zh) | 基于听诊音频谱分析的多功能数字式听诊器电路 | |
CN101449970A (zh) | 生物电放大器 | |
CN201641998U (zh) | 一种基于rls自适应滤波器的胃电慢波信号检测装置 | |
CN102100553A (zh) | 一种基于rls自适应滤波器的胃电慢波信号检测方法 | |
Ali et al. | A novel technique for extraction foetal electrocardiogram using adaptive filtering and simple genetic algorithm | |
JP2009039540A (ja) | ペースメーカパルス検出のための方法及びシステム | |
CN107928636B (zh) | 一种具有温度补偿功能的脉诊仪 | |
Irianto et al. | A low-cost electro-cardiograph machine equipped with sensitivity and paper speed option | |
CN102100554A (zh) | 一种基于rls自适应滤波器的胃电慢波信号检测装置 | |
CN113100735B (zh) | 一种基于静脉血流图的血栓检测装置及方法 | |
CN212853501U (zh) | 心音信号采集分析系统 | |
CN103654840A (zh) | 信号处理系统、信号处理装置和存储介质 | |
CN109498055B (zh) | 心音信号采集分析系统及信号分析方法 | |
Lee et al. | Intelligent Stethoscope System and Diagnosis Platform With Synchronized Heart Sound and Electrocardiogram Signals |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20110622 |