CN102302361A - 一种心率及身体状态监测装置及方法 - Google Patents
一种心率及身体状态监测装置及方法 Download PDFInfo
- Publication number
- CN102302361A CN102302361A CN201110174039A CN201110174039A CN102302361A CN 102302361 A CN102302361 A CN 102302361A CN 201110174039 A CN201110174039 A CN 201110174039A CN 201110174039 A CN201110174039 A CN 201110174039A CN 102302361 A CN102302361 A CN 102302361A
- Authority
- CN
- China
- Prior art keywords
- circuit
- heart rate
- axle
- fpga
- 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.)
- Granted
Links
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明提供一种心率及身体状态监测装置,传感器安装在腰带上,本装置包括数据采集模块、FPGA控制模块和DSP数据处理模块。受试者心率和身体状态改变时,传感器采集由此产生的体动信号;其中心脏跳动引起的体动信号经模拟放大滤波电路处理后,和身体状态变化引起的体动信号分别经A/D转换电路输出至FPGA;FPGA选择一路或多路传感器输出传至DSP进行体动信号处理,计算心率并判别身体状态,得到的数据反馈到FPGA;FPGA根据DSP反馈的信息,将数据通过串口外接通信模块发送给监护中心工作站。本发明能完成多路信号的实时、快速、精确采集,实现人体生理信息获取,且与穿戴技术结合,实现生理信息的长期动态监测与辨识。
Description
技术领域
本发明涉及医疗器械领域,具体涉及一种心率及身体状态监测装置。
背景技术
21世纪被称为“银发世纪”,世界范围内的老龄化浪潮滚滚而来。根据中国人口信息研究中心的调查统计,2000年中国60岁以上人口比例为10.31%,65岁以上人口占总人口的比例为7.17%,按照国际标准来衡量,中国已经步入老龄化社会,而且在今后几十年其老龄化趋势将愈加明显。另外,随着社会经济的发展,居住方式的变化,家庭结构的小型化,以及人口流动的加速,子女数的减少,其老人家庭空巢率正在不断的加大。因此,解决“空巢家庭”中老年人的家庭监护问题日趋突出。
随着年龄的增加,人体解剖组织结构和生理代谢发生一系列变化,机体功能衰退,应变能力减退,骨骼也变得较为松脆,这些生理或其他原因引发的变化常常可以通过人体的心律和身体状态表现出来。当人生病时,特别是心脏病发作时,心率都会发生明显的改变。另外,由于自身疾病如心脑血管疾病或外界影响等因素,人的身体状态也常常发生改变。根据美国国家安全委员会的报告指出,在65岁以上的人口中,跌倒所造成的死亡居所有意外死亡原因的第一位,占此年龄段意外死亡的33%。近年来,我国心脑血管疾病发病率持续上升,每年有54.4万人心脏性猝死。面对越来越严峻的现实,我们应该做好相对的预防措施。首先,心脑血管疾病的发生是有先兆的,如果刚出现病灶时就立刻救治,很多人是能够缓解过来的。现在有很多心脏性猝死的病人是由于发病时无人知晓,没有得到及时的救治,才导致严重的后果。因此建立一套完善的心率和身体状态监测与通信系统是相当重要的。如果能实时地监测人的心率和身体状态,对及时发现心脑血管疾病的病情是有很大帮助的。另外,心率与身体状态也能反映了人的生理信息,如疲劳、醉酒等也会引起心率与身体状态的改变。所以,实时地监测心率与身体状态不仅对于病人,对于健康人也是有用的。但是,临床上使用的心、肺功能检测方法需要医护人员使用专业设备进行操作,且其中大多数方法需要在人体指定部位贴附传感器,人体会感觉不适,不便用于心、肺功能的日常监护。因此有必要研究一种无感觉、且可长时间使用的非临床心肺功能实时监测方法,用于长期监测心肺工作情况并在其出现异常时发出警示。
发明内容
针对现有技术之不足,本发明提供一种心率及身体状态监测装置,将传感器安装在腰带上,本装置包括数据采集模块、FPGA控制模块和DSP数据处理模块;
所述数据采集模块包括传感器电路、模拟放大滤波电路和A/D转换电路,其中传感器电路包括一个固定腰带里的电阻应变式传感器和一个三轴加速度传感器,固定腰带里的电阻应变式传感器用于测量心脏跳动引起的体动信号,主要用来监测心率;三轴加速度传感器用来监测身体状态引起的体动信号;模拟放大滤波电路包括前置放大电路、隔直电路、陷波电路、低通滤波电路及三级放大电路,前置放大电路的输出端接至隔直电路的输入端,第一级放大电路的输入端与隔直电路输出端相连接,第一级放大电路输出端接至陷波电路输入端,陷波电路输出端接至第二级放大电路输入端,第二级放大电路输出端接至低通滤波电路输入端,低通滤波电路输出端再接至陷波电路输入端,陷波电路输出端接至后级放大电路输入端。
数据采集模块主要通过电阻应变式传感器和三轴加速度传感器采集体动信号并转换成电信号,送到后级模拟放大滤波电路中进行信号调理,经A/D转换电路转换为数字信号送往FPGA内部进行预处理。
所述FPGA控制模块,主要进行数据采集、滤波、传输、存储、显示及任务管理。FPGA控制模块包括电源电路,无线通信模块,SD存储接口,LCD,按键、开关和LED显示。FPGA采集A/D转换电路输出的体动信号数字量,并对这些数字量进行滤波,将滤波以后数据传递给DSP进行数据处理,同时数据存储在SD卡中,并在LCD上实时显示体动信号波形,FPGA根据DSP反馈的指令将采集回来体动信号数据发送给具有串口接口通信模块。
所述DSP数据处理模块包括JTAG接口电路、SRAM存储电路和DSP与FLASH连接电路,DSP对FPGA输出的体动信号进行特征提取,监测心率和身体状态。DSP外设SRAM存储器来扩展程序运行期间的存储空间,DSP外设FLASH存储器,实现DSP程序的装载。
所述传感器电路中的电阻应变式传感器输出端连接至模拟放大滤波电路输入端,模拟放大滤波电路输出端连接至一个A/D转换电路输入端,三轴加速度传感器的三个输出端分别连接至三个A/D转换电路输入端,四个A/D转换电路的输出端连至FPGA控制模块的GPIO,FPGA控制模块通过GPIO与DSP数据处理模块I/O输入端相连,FPGA通过串口外接通信模块,FPGA与SD卡通过SD存储接口连接。
所述传感器电路中电阻应变式传感器采集心脏跳动引起的体动信号,微弱的体动信号输出至放大电路,放大电路的输出经由低通滤波电路和陷波电路进行滤波;滤波后输出的信号输入A/D转换电路,产生的数字信号经电平转换并行输出至FPGA控制模块的GPIO;三轴加速度传感器采集身体状态变化引起的体动信号,输出引脚XYZ分别接至三个相同的A/D转换电路的串行输入端,产生的数字信号经电平转换并行输出至FPGA控制模块的GPIO,FPGA控制模块根据输入的键控信息,选择一路或多路传感器输出产生的数字信号传递至DSP数据处理模块I/O输入端进行算法处理。
采用本发明装置监测心率及身体状态方法,具体按以下步骤执行:
步骤一:采集信号。首先为受试者带上心率及身体状态监测装置,当受试者心率和身体状态发生改变时,传感器采集由此产生的信号,并转换为交流电压输出;
步骤二:转换信号。经电阻应变式传感器输出的信号经过模拟放大滤波电路处理后,同三轴加速度传感器输出的信号分别输出至四个A/D转换电路转换为数字信号输出至FPGA中;
步骤三:信号预处理。FPGA采集A/D转换电路输出的体动信号,对采集到的信号进行滤波,将滤波后的数据输出至DSP;
步骤四:心率及身体状态监测。DSP模块进行体动信号的处理,计算心率并判别身体状态,所得到的数据反馈到FPGA,FPGA将数据存储在SD卡中,并在LCD上实时显示心率信号波形,FPGA根据DSP反馈的信息,将数据通过串口外接通信模块发送给监护中心工作站。
所述步骤四中,DSP模块对体动信号进行滤波、J波群检测、体动周期计算,从而计算心率。具体方法如下:
步骤1:选择合适小波基,确定分解尺度,对采集的原始体动信号进行尺度分解;
步骤2:确定阈值函数,对体动信号进行小波去噪;
阈值函数如式(1)所示,
步骤3:利用小波模极大值提取局部极值点,局部极值点位置记为J1;
步骤4:选取单调上升幅度阈值,具体地:任选t为幅度的分割阈值,所取的上升沿幅度个数占所有上升沿个数的比例为w0,平均幅值为μ0;去除的幅度个数占所有上升沿个数的比例为w1,平均幅值为μ1。信号的总平均上升沿幅度的平均值为:μ=w0×μ0+w1×μ1。从最小上升幅度值到最大上升幅度值遍历t,当t=t1,使g=w0×(μ0-μ)2+w1×(μ1-μ)2最大时,t1即为最佳的幅度分割阈值。
步骤5:判断J1中的点的幅值是否大于阈值t1,将大于阈值的点保存到J2中。
步骤6:确定J波中心点的位置。由于成人平静时的心率不超过400ms,系统采样率为500Hz,设定200个点为标准,在J2中去除相邻点之间的距离小于200个点的部分,剩下的点构成J。
步骤7:计算心率。利用这些J波的中心点位置计算一个平均周期即心动周期,心动周期的倒数就是所求心率。J为J波中心点对应的采样序数。
所述步骤四中,判别身体状态的方法是通过身体所处的方向(水平或垂直)、加速度能量变化和(是否在指定阈值区间内)来判断处于平卧、直立、行走、跑步等状态。
所述身体状态包括以下四种情况:
身体所处方向水平,三轴加速度能量变化和小于阈值,则身体属于平卧状态;
身体所处方向垂直,三轴加速度能量变化和为零,则身体处于直立状态;
身体所处方向垂直,三轴加速度能量变化和在阈值区间里,则身体处于行走状态;
身体所处方向垂直,三轴加速度能量变化和大于阈值,则身体处于跑步状态;
具体判别过程按如下步骤进行:
步骤1:计算加速度能量和E,并判定身体所处方向。
身体处于一个三维立体空间中,坐标轴x、y、z满足右手定则,z轴与x轴、y轴相互正交于原点。身体所在平面平行于z轴方向定义为垂直方向;垂直于z轴方向定义为水平方向。
三轴加速度能量变化和E的定义:
其中x′ac(n)、y′ac(n)、z′ac(n)表示三轴加速度x轴、y轴、z轴在采样点n时输出变化量,T为求和的时间长度,fs表示三轴加速度输出采样率。
判定身体所处方向(处于垂直或水平状态)的方法如下:当人在站立或行走时在z轴的加速度值约为g,x轴、y轴的加速度值接近于0,全局加速度方向平行于z轴;当身体在平卧时,在z轴方向上的加速度值将会减小接近于0,而x轴和y轴中的一个轴的加速度值约为g。由于重心偏离,全局加速度方向发生改变,因此计算全局加速度与z轴夹角φ来判别身体状态。如果φ<20°,或φ>160°,判定身体处于垂直状态,否则20°<φ<160°时,判定身体处于水平状态。
其中φ(n)第n时刻全局加速度与z轴夹角,x′dc(n)、y′dc(n)、z′dc(n)表示三轴加速度x轴、y轴、z轴第n时刻输出的直流分量。
步骤2:根据判定的身体所处方向和三轴加速度能量和E,确定身体状态。
有益效果
本发明装置同时具备DSP与FPGA的优势,支持更高的计算处理能力,可以很好地完成多路信号、大量信号的实时、快速、精确的采集,适用于恶劣复杂的环境,且开发成本低,具有较高的使用价值,有很好的应用前景。将传感器安装在腰带上,这与人们正常使用腰带的情况是相同的,佩戴者不会感觉到有监测装置存在,较为舒适,由于人们佩戴腰带的时间长,监测的时间也相应较长。可使在家庭,办公等环境下,长时间实时监护,还可以辅助医生在病人完全放松无心理负担的情况下更准确地了解心肺功能,预测和诊断心肺疾病,评估生理或心理测试结果,以及监测睡眠质量等。实现人体生理信号获取和处理平台,并且与穿戴技术相结合,实现生理信息的连续动态监测与辨识。本发明具备无创、可移动操作、使用简便、可长时间连续使用特点,实现远程监控,减少佩戴者的心理和生理负担。
附图说明
图1本发明总体结构图;
图2本发明数据采集模块框图;
图3本发明实施例全桥测量电路;
图4本发明实施例三轴加速度传感器电路原理图;
图5本发明实施例前置放大电路原理图;
图6本发明实施例第一二级放大电路原理图;
图7本发明实施例陷波电路原理图;
图8本发明实施例的低通滤波电路原理图;
图9本发明实施例的后级放大电路原理图;
图10本发明实施例FPGA与A/D芯片连接原理图;
图11本发明实施例5V供电稳压电路原理图;
图12本发明实施例3.3V供电稳压电路原理图;
图13本发明实施例2.5V供电稳压电路原理图;
图14本发明实施例1.2V供电稳压电路原理图;
图15本发明实施例FPGA与UART接口连接原理图;
图16本发明实施例FPGA与SD接口连接原理图;
图17本发明实施例FPGA与LCD触摸屏连接原理图;
图18本发明实施例DSP与FPGA连接原理图;
图19本发明实施例JTAG接口电路原理图;
图20本发明实施例SRAM存储电路原理图;
图21本发明实施例FLASH和DSP的连接电路原理图;
图22本发明装置工作流程图;
图23本发明心率计算流程图;
图24本发明实施例标定的J波位置信号图;
图25本发明判别身体状态流程图。
具体实施方式
本发明提供一种心率及身体状态监测装置,将传感器安装在腰带上,本装置包括数据采集模块、FPGA控制模块和DSP数据处理模块,总体结构如图1所示,数据采集模块如图2所示。首先电阻应变式传感器和三轴加速度传感器采集体动信号,电阻应变式传感器一路输出经模拟放大滤波电路处理送入A/D转换电路,三轴加速度传感器三路输出直接送入A/D转换电路,四路传感器的输出经A/D转换电路转换输出至FPGA处理。FPGA内部的NoisII处理器对挂载在Avalon总线上的存储、显示和通信设备进行调用。由DSP模块完成体动信号的处理,提取心率和身体状态的信息,进行心率计算和身体状态判别,通过LCD实时显示心率和预警信息;FPGA通过串口外接通信模块完成与上位机的通信;通过SD接口,将大量的体动信号数据存入SD卡中,为病情的诊断提供充足的数据来源。另外还有按键、开关和LED对装置进行设置和工作状态提示。
受试者配带本发明装置后,心脏跳动引起的体动信号通过传感器转换为交流电压输出,经过模拟放大滤波电路处理后,信号进入A/D转换电路转换成数字信号;身体状态引起的体动信号通过三轴加速度传感器输出至A/D转换电路转换成数字信号。A/D转换电路输出的信号以中断方式、连续方式或定时方式输出至FPGA控制模块,FPGA控制模块通过串口连接通信模块将监控对象的心率和身体状态传送到监护中心。当受试者的心率或身体状态发生异常,医院的中心工作站终端显示界面进行异常提示,本装置向中心工作站的个人急救中心发出报警信息。本装置设有输入输出设备,输出设备可以是LED指示灯,指示灯常亮提示装置处于工作状态,不断闪烁指示装置电力不足;输入设备可以是按键或开关,控制装置的操作和信息收发。
所述数据采集模块包括传感器电路、模拟放大滤波电路和A/D转换电路,传感器电路包括一个固定腰带里的电阻应变式传感器和一个三轴加速度传感器,经电阻应变式传感器输出的信号,受外界的作用力作用,使传感器的弹性梁发生形变,使贴在弹性梁部位的应变片阻值发生变化,在激励电压的作用下,输出的电压信号也发生正比例变化,电阻应变式传感器全桥测量电路如图3所示,三轴加速度传感器型号为MMA7260Q,MMA7260Q的g_Select1和g_Select2两引脚分别接FPGA芯片两个GPIO,输出为高电平,此时设定灵敏度200mV/g,输出引脚X,Y,Z分别接RC低通滤波电路,R阻值1kΩ,C容值0.1μf,电路如图4所示。
模拟放大滤波电路中,传感器电路输出的体动信号,首先经过高稳定性、高共模抑制比的差分前置放大电路消除共模信号,保留差模信号;然后经过隔直电路滤除仪器等引起的直流噪声;所得信号经由第一级放大电路传入50Hz陷波器消除50Hz工频干扰;再经过第二级放大电路,送入低通滤波器中去除高频成分,保留有用的低频信号;再次送入50Hz陷波器;再进行隔直滤波去除干扰信号;最后送入后级放大倍数可调的放大电路根据实际需要进行放大得到体动信号。
模拟放大滤波电路中,前置放大电路型号为AD620AN,电路如图5所示,图中,隔直电路采用两个电解电容正极对正极串联,第一二级放大电路型号为ICL7650,电路如图6所示,陷波电路滤波芯片选用UAF42AP,电路如图7所示,其中,2脚输入,6脚输出。10,9脚分别为UAF42AP的正负电源。为了去除高频干扰,在电源端需连接0.1uf的瓷片电容到地。第一、二级放大电路放大倍数为40倍,后级放大电路的放大倍数根据需要可以调节。低通滤波电路采用两个2阶压控低通滤波器级联,电路如图8所示,其中有源器件型号为ICL7650,滤波电容的选择可根据实际调整。后级放大电路型号为ICL7650,电路图如图9所示,通过可变电阻实现增益可调。
AD转换电路选用型号为ADC0832AP芯片,分辨率8位,其中CS为片选信号,DI为读数据选通信号,CLK为时钟信号,DO为转换后的数据信号,与FPGA的GPIO相连,数据输入通道CH0、CH1与采集电路的输出相连。FPGA与A/D芯片连接原理如图10所示。
FPGA控制模块电路包括电源电路,无线通信模块,SD存储接口,LCD,按键、开关和LED显示。
FPGA电源电路需要使用5V、3.3V、2.5V和1.2V四组直流电源。5V由三端稳压芯片LM2940得到,电路如图11所示;3.3V采用TPS75833,电路如图12所示;2.5V采用TPS76825Q,如图13所示;1.2V部分采用TPS76801Q,电路如图14所示。
FPGA与无线通信模块采用MAX232芯片来完成RS-232-C串行接口的电平转换。MAX3232的11脚T1IN与FPGA的一个GPIO相连,用网络标号TXD表示,12脚R1OUT与FPGA的另一个GPIO相连,用网络标号RXD实现相互连接。13脚R1IN与DB9型连接器的3号引脚相连,14脚T1OUT与DB9型连接器的2号引脚相连,电路如图15所示。MAX3232芯片的电容C3、C4、C5、C6分别接在1和3、4和5、2和6、6和15之间,C7是芯片电源引脚滤波电容。
FPGA与SD接口电路如图16所示,7脚D0、5脚CLK、2脚CMD、1脚CD/D3分别与FPGA的GPIO相连,分别用网络标号SD DAT、SD CLK、SD CMD和SD DATA3来实现连接。
FPGA与LCD接口电路如图17所示,LCD模块负责显示体动信号波形、心率等信息,其中V0为LCD显示对比度调整端,接正电源对比度最弱,接地对比度最高,通过10kΩ电位器调整。寄存器选择引脚RS、读写信号引脚RW、芯片使能引脚E、双向数据线引脚D7~D0均与FPGA控制模块的GPIO连接,分别用网络标号LCD_RS、LCD_RW、LCD_E和LCD[7..0]来实现连接。
FPGA与DSP连接,主要完成对ADC0832AP控制、数据的传输以及相应的逻辑控制,电路如图18所示。装置采用基于DSP与FPGA协同处理,DSP芯片型号为TMS320VC5410,FPGA芯片型号为EP2C6Q208。
DSP芯片TMS320VC5410主要进行数据的计算及传输,通过向FPGA发送准备好信号,使FPGA通过A/D转换电路采集数据,并将数据传输给TMS320VC5410,完成输入数据的计算、打包等处理。
DSP模块电路包括JTAG接口电路、SRAM存储电路和DSP与FLASH连接电路。
JTAG接口电路用于仿真和测试,JTAG接口共有14针,其中测试数据输入引脚TDI、测试数据输出引脚TDO、测试模式选择引脚TMS和测试时钟输入引脚TCK分别同FPGA的GPIO相连。JTAG接口电路如图19所示。
DSP的SRAM芯片型号为IS61LV51216-12T,用来存储DSP在工作过程中所使用到的数据。电路原理如图20所示,芯片11脚和13脚接电源端,12脚和34脚接地端;6脚为芯片使能端,41脚为读使能端,17脚为写使能端,39脚和40脚接地;A0~A18为19位地址线;I/O0~I/O15为16位数据线,这些引脚均与FPGA的GPIO相连。
本发明软件部分主要包括数据采集、数据处理及传输。其中,由于采用串口外接无线通信模块的协议已经固化在模块中,因此,只需要在DSP和FPGA中编写数据采集、处理和收发程序。本发明软件流程如图22所示。DSP和FPGA主要是进行初始化、逻辑控制、前端数据采集和数据的处理传输。系统上电复位后,首先完成初始化,包括FPGA、DSP的初始化;然后FPGA等待接收DSP准备好信号,FPGA根据键入信息控制A/D转换器完成数据的采集;再由FPGA向DSP发出写指令,FPGA把读取数据并发送给DSP,由DSP完成数据的计算、识别等处理,进行异常判断,如有异常,DSP将数据发送给FPGA,FPGA通过无线通信模块发送给上位机。
采用本发明装置监测心率及身体状态方法,具体按以下步骤执行:
步骤一:采集信号。首先为受试者带上心率及身体状态监测装置,当受试者心率和身体状态发生改变时,传感器采集由此产生的信号,并转换为交流电压输出;
步骤二:转换信号。经电阻应变式传感器输出的信号经过模拟放大滤波电路处理后,同三轴加速度传感器输出的信号分别输出至四个A/D转换电路转换为数字信号输出至FPGA中;
步骤三:信号预处理。FPGA采集A/D转换电路输出的体动信号,对采集到的信号进行滤波,将滤波后的数据输出至DSP;
步骤四:心率及身体状态监测。DSP模块进行体动信号的处理,计算心率并判别身体状态,所得到的数据反馈到FPGA,FPGA将数据存储在SD卡中,并在LCD上实时显示心率信号波形,FPGA根据DSP反馈的信息,将数据通过串口外接通信模块发送给医院的监护中心工作站。
所述步骤四中,DSP模块对体动信号进行滤波、J波群检测、体动周期计算,从而计算心率。具体流程如图23所示,具体步骤如下:
步骤1:选择合适小波基,确定分解尺度,对采集的原始体动信号进行尺度分解;
本具体实施方式中选择了尺度函数比较接近体动信号波形且小波函数对称性较好的sym8小波基。对含噪的体动信号进行小波分解,分解尺度j需满足公式(5)。
其中fc为待分析信号x(t)中心频率,fs为体动信号的采样频率,N为信号长度。
步骤2:确定阈值函数,对体动信号进行小波去噪,
阈值函数如式(1)所示,
步骤3:利用小波模极大值提取局部极值点,局部极值点位置记为J1;
对于通过小波分解得到去噪后的小波系数,根据极大值处的二阶导数大于或小于零来得到,满足(f(j)-f(j-1))*(f(j)-f(j+1))>0,则点j即为局部极大值点,并将局部极大值点j位置记为J1。
步骤4:选取单调上升幅度阈值,具体地:任选t为幅度的分割阈值,所取的上升沿幅度个数占所有上升沿个数的比例为w0,平均幅值为μ0;去除的幅度个数占所有上升沿个数的比例为w1,平均幅值为μ1。信号的总平均上升沿幅度的平均值为:μ=w0×μ0+w1×μ1。从最小上升幅度值到最大上升幅度值遍历t,当t=t1,使g=w0×(μ0-μ)2+w1×(μ1-μ)2最大时,t1即为最佳的幅度分割阈值。
步骤5:判断J1中的点的幅值是否大于阈值t1,将大于阈值的点保存到J2中。
步骤6:确定J波中心点的位置。由于成人平静时的心率不超过400ms,系统采样率为500Hz,设定200个点为标准,在J2中去除相邻点之间的距离小于200个点的部分,剩下的点构成J,标定的J波位置如图24所示。
步骤7:计算心率。利用这些J波的位置计算一个平均周期即心动周期,心动周期的倒数就是所求心率。J为J波中心点对应的采样点数。
心率监测对心律失常(期前收缩、心动过速、扑动与颤动、传导阻滞等)的监测具有重要意义。
所述步骤四中,判别身体状态的方法是通过身体所处的方向(水平或垂直)、加速度能量变化和(是否在指定阈值区间内)来判断处于平卧、站立、行走、跑步等状态。身体状态的监测通过三轴加速度传感器来实现,即把加速度值分解到x、y、z轴上,根据基态的值,计算出加速度和x、y、z轴之间的夹角,从而判断身体所处状态。其中,x轴的正方向指向人前进的方向,y轴的正方向指向身体的正左方,z轴的正方向指向身体的正上方然后送入A/D转换电路中。
判别身体状态具体过程按如下步骤进行:
步骤1:计算加速度能量和E,并判定身体所处方向。
身体处于一个三维立体空间中,坐标轴x、y、z满足右手定则,z轴与x轴、y轴相互正交于原点。身体所在平面平行于z轴方向定义为垂直方向;垂直于z轴方向定义为水平方向。
三轴加速度能量变化和E的定义:
其中x′ac(n)、y′ac(n)、z′ac(n)表示三轴加速度x轴、y轴、z轴在采样点n时输出变化量,T为求和的时间长度,fs表示三轴加速度输出采样率。
判定身体所处方向(处于垂直或水平状体)的方法如下:当人在站立或行走时全局加速度在z轴的加速度值约为g,x轴、y轴的加速度值接近于0,全局加速度方向平行于z轴;当人体在平卧时,在z轴方向上的加速度值将会减小接近于0,而x轴和y轴中的一个轴的加速度值约为g。由于重心偏离,全局加速度方向发生改变,因此计算全局加速度与z轴夹角φ来判别身体状态,如果φ<20°,或φ>160°,判定身体处于垂直状态,否则20°<φ<160°时,判定身体处于水平状态。
其中φ(n)第n时刻全局加速度与z轴夹角,x′dc(n)、y′dc(n)、z′dc(n)表示三轴加速度x轴、y轴、z轴第n时刻输出的直流分量。
步骤2:根据判定的身体所处方向和三轴加速度能量和E,确定身体状态。
根据监测三轴加速度输出,判定身体所处方向和三轴加速度能量和E,确定身体状态包括以下几种情况:
身体所处方向水平,三轴加速度能量变化和小于阈值E1,则身体属于平卧状态;
身体所处方向垂直,三轴加速度能量变化和为零,则身体处于直立状态;
身体所处方向垂直,三轴加速度能量变化和在阈值区间{E2,E3}里,则身体处于行走状态;
身体所处方向垂直,三轴加速度能量变化和大于阈值E4,则身体处于跑步状态。
本实施例以判定身体跌倒为例,判定身体处于跌倒状态的方法如下:
首先判断身体是否从垂直状态转变接近水平状态(条件1),再判定全局加速度与z轴夹角变化率是否超过阈值(条件2),从而确定身体是否处于跌倒状态。
初始时身体处于直立状态,跌倒后大多数情况下身体处于水平或接近水平状态。因此先判定身体是否从垂直状态转变接近水平状态,当人在站立或行走时z轴的加速度值约为g,x轴、y轴的加速度值接近于0,当身体发生跌倒时,由于重心偏离,z轴方向上的加速度值将会减小,当身体由垂直变为水平状态时,z轴方向上的加速度值变为0,而x轴和y轴中的一个轴或两个轴的加速度值也会有相应变化。因此计算全局加速度与z轴夹角φ,φ如果从垂直状态取值φ1(t1)(φ<20°,或φ>160°)变化到水平状态取值φ2(t2)(20°<φ<160°),则认为满足了跌倒发生的条件1;根据身体发生跌倒的形态会发现,发生跌倒时身体重心的改变速度会比一般日常活动中的重心改变的速度快,且幅度较大。因此,对全局加速度与z轴夹角进行差分运算求夹角变化率判定结果是否超过设定的阈值λ,超过满足条件2判定为跌倒。判定身体跌倒的流程如图25所示。
本发明装置同时具备DSP与FPGA的优势,支持更高的计算处理能力,可以很好地完成多路信号、大量信号的实时、快速、精确的采集,适用于恶劣复杂的环境,且开发成本低,具有较高的使用价值。
Claims (4)
1.一种心率及身体状态监测装置,将传感器安装在腰带上,其特征在于:该监测装置包括:数据采集模块、FPGA控制模块和DSP数据处理模块;
所述数据采集模块包括传感器电路、模拟放大滤波电路和A/D转换电路,其中传感器电路包括一个固定腰带里的电阻应变式传感器和一个三轴加速度传感器,固定腰带里的电阻应变式传感器用于监测心脏跳动引起的体动信号,主要用来监测心率;三轴加速度传感器用来监测身体状态变化引起的体动信号;模拟放大滤波电路包括前置放大电路、隔直电路、陷波电路、低通滤波电路及三级放大电路,前置放大电路的输出端接至隔直电路的输入端,第一级放大电路的输入端与隔直电路输出端相连接,第一级放大电路输出接至陷波电路输入端,陷波电路输出端接至第二级放大电路输入端,第二级放大电路输出接至低通滤波电路输入端,低通滤波电路输出端再接至陷波电路输入端,陷波电路输出端接至后级放大电路输入端;
所述FPGA控制模块包括电源电路,无线通信模块,SD存储接口,LCD,按键、开关和LED显示;
所述DSP数据处理模块包括JTAG接口电路、SRAM存储电路和DSP与FLASH连接电路;
所述传感器电路中的电阻应变式传感器输出端连接至模拟放大滤波电路输入端,模拟放大滤波电路输出端连接至一个A/D转换电路输入端,三轴加速度传感器的三个输出端分别连接至三个A/D转换电路输入端,四个A/D转换电路的输出端连至FPGA控制模块的GPIO,FPGA控制模块通过GPIO与DSP数据处理模块I/O输入端相连,FPGA通过串口外接通信模块,FPGA与SD卡通过SD存储接口连接。
2.采用权利要求1所述的心率及身体状态监测装置监测心率及身体状态的方法,其特征在于按如下步骤进行:
步骤一:采集信号,首先为受试者配带上心率及身体状态监测装置,当受试者心率和身体状态发生改变时,传感器采集由此产生的体动信号,并转换为交流电压输出;
步骤二:转换信号,经电阻应变式传感器输出的信号经过模拟放大滤波电路处理后,同三轴加速度传感器输出的信号分别输出至四个A/D转换电路转换为数字信号输出至FPGA中;
步骤三:信号预处理,FPGA采集A/D转换电路输出的体动信号,对采集到的信号进行滤波,将滤波后的信号输出至DSP;
步骤四:心率及身体状态监测,DSP模块进行体动信号的处理,计算心率并判别身体状态,所得到的数据反馈到FPGA,FPGA将数据存储在SD卡中,并在LCD上实时显示体动信号波形,FPGA根据DSP反馈的信息,将数据通过串口外接通信模块发送给监护中心工作站。
3.根据权利要求2所述的心率及身体状态监测方法,其特征在于:所述步骤四中,DSP模块对体动信号处理,计算心率是DSP模块对体动信号进行滤波、J波群检测、心动周期计算,从而计算心率,计算心率具体方法如下:
步骤1:选择合适小波基,确定分解尺度,对采集的体动信号进行尺度分解;
步骤2:确定阈值函数,对体动信号进行小波去噪;
步骤3:利用小波模极大值提取局部极值点,局部极值点位置记为J1;
步骤4:选取单调上升幅度阈值,具体地:任选t为幅度的分割阈值,所取的上升沿幅度个数占所有上升沿个数的比例为w0,平均幅值为μ0;去除的幅度个数占所有上升沿个数的比例为w1,平均幅值为μ1;信号的总平均上升沿幅度的平均值为:μ=w0×μ0+w1×μ1,从最小上升幅度值到最大上升幅度值遍历t,当t=t1,使g=w0×(μ0-μ)2+w1×(μ1-μ)2最大时,t1即为最佳的幅度分割阈值;
步骤5:判断J1中的点的幅值是否大于阈值t1,将大于阈值的点保存到J2中;
步骤6:确定J波中心点的位置;
步骤7:计算心率,利用这些J波的中心点位置计算一个平均周期即心动周期,心动周期的倒数就是所求心率,J为J波中心点对应的采样序数,
4.根据权利要求2所述的心率及身体状态监测方法,其特征在于:所述步骤四中,DSP模块进行体动信号处理,判别身体状态的方法按如下步骤进行:
步骤1:计算加速度能量和E,并判定身体所处方向;
身体处于一个三维立体空间中,坐标轴x、y、z满足右手定则,z轴与x轴、y轴相互正交于原点,身体所在平面平行于z轴方向定义为垂直方向,垂直于z轴方向定义为水平方向,
三轴加速度能量变化和E的定义:
其中x′ac(n)、y′ac(n)、z′ac(n)表示三轴加速度x轴、y轴、z轴在采样点n时输出变化量,T为求和的时间长度,fs表示三轴加速度输出采样率;
判定身体所处方向(处于垂直或水平状体)的方法如下:当人在站立或行走时在z轴的加速度值约为g,x轴、y轴的加速度值接近于0,全局加速度方向平行于z轴;当身体在平卧时,在z轴方向上的加速度值将会减小接近于0,而x轴和y轴中的一个轴的加速度值约为g,由于重心偏离,全局加速度方向发生改变,因此计算全局加速度与z轴夹角φ来判别身体状态,如果φ<20°,或φ>160°,判定身体处于垂直状态,否则20°<φ<160°时,判定身体处于水平状态;
其中φ(n)第n时刻全局加速度与z轴夹角,x′dc(n)、y′dc(n)、z′dc(n)表示三轴加速度x轴、y轴、z轴第n时刻输出的直流分量;
步骤2:根据判定的身体所处方向和三轴加速度能量和E,确定身体状态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110174039 CN102302361B (zh) | 2011-06-24 | 2011-06-24 | 一种心率及身体状态监测装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110174039 CN102302361B (zh) | 2011-06-24 | 2011-06-24 | 一种心率及身体状态监测装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102302361A true CN102302361A (zh) | 2012-01-04 |
CN102302361B CN102302361B (zh) | 2013-02-13 |
Family
ID=45376356
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110174039 Expired - Fee Related CN102302361B (zh) | 2011-06-24 | 2011-06-24 | 一种心率及身体状态监测装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102302361B (zh) |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102662686A (zh) * | 2012-03-09 | 2012-09-12 | 中国科学院微电子研究所 | Fpga的加载方法及其装置 |
CN102743154A (zh) * | 2012-07-02 | 2012-10-24 | 深圳市长江力伟股份有限公司 | 一种便携式人体体征监视系统和装置 |
CN103605304A (zh) * | 2013-11-09 | 2014-02-26 | 电子科技大学 | 一种多路微弱信号采集处理电路 |
CN103837696A (zh) * | 2012-11-21 | 2014-06-04 | 通用汽车环球科技运作有限责任公司 | 传感器装置和通过陷波滤波器进行信号处理的方法 |
WO2014086261A1 (zh) * | 2012-12-07 | 2014-06-12 | 河南华南医电科技有限公司 | 一种心脏电刺激信号分配电路 |
CN104490401A (zh) * | 2015-01-15 | 2015-04-08 | 张�杰 | 基于加速度数据的智能终端进行人体行为自动识别的方法 |
CN104510460A (zh) * | 2014-12-22 | 2015-04-15 | 哈尔滨师范大学 | 一种基于物联网的独居老人自动监测系统 |
CN104678745A (zh) * | 2015-02-10 | 2015-06-03 | 柳州市金旭节能科技有限公司 | 一种具有跌倒检测功能的手表 |
CN105223424A (zh) * | 2015-11-13 | 2016-01-06 | 张希梅 | 变电站在线监测装置 |
WO2016019894A1 (en) * | 2014-08-07 | 2016-02-11 | Beijing Zhigu Tech Co., Ltd. | Dominant limb identification method and device |
CN105982659A (zh) * | 2016-03-17 | 2016-10-05 | 厦门鑫奥力电器有限公司 | 一种用于测量心率的耳夹及其控制方法 |
CN106031636A (zh) * | 2015-03-11 | 2016-10-19 | 济宁智城中科信息技术有限公司 | 一种跌倒检测方法及相应的可穿戴监护装置及监护方法 |
CN106344006A (zh) * | 2016-11-03 | 2017-01-25 | 太原理工大学 | 基于极点对称模态分解和支持向量机的j波检测方法 |
CN106725897A (zh) * | 2017-01-24 | 2017-05-31 | 新博医疗技术有限公司 | 一种呼吸跟踪监测通讯系统 |
CN106725394A (zh) * | 2016-12-30 | 2017-05-31 | 北京工业大学 | 一种血压数据采集装置及方法 |
CN107088062A (zh) * | 2017-05-08 | 2017-08-25 | 歌尔科技有限公司 | 一种心率模组、采集心率的电子装置和心率采集方法 |
CN109269395A (zh) * | 2018-11-01 | 2019-01-25 | 重庆地质矿产研究院 | 一种岩土工程结构远程交互式应变监测装置及监测方法 |
CN109350050A (zh) * | 2018-11-20 | 2019-02-19 | 锐捷网络股份有限公司 | 心跳间期检测方法及装置 |
CN109387279A (zh) * | 2018-12-19 | 2019-02-26 | 南京工程学院 | 一种对汽轮机发电机末级的振动信号的检测电路 |
CN109635355A (zh) * | 2018-11-19 | 2019-04-16 | 北京时代民芯科技有限公司 | 一种面向gpio的频率可调的滤波电路 |
CN110031653A (zh) * | 2019-05-09 | 2019-07-19 | 保定开拓精密仪器制造有限责任公司 | 三轴石英挠性加速度计小型信号采样装置及其加工工艺 |
CN110650685A (zh) * | 2017-03-24 | 2020-01-03 | 爱尔西斯有限责任公司 | 评估人的心理生理状态的方法 |
CN112870555A (zh) * | 2021-02-07 | 2021-06-01 | 丹源医学科技(杭州)有限公司 | 一种休息起搏频率的自动调节装置和心脏起搏装置 |
CN112971720A (zh) * | 2021-02-07 | 2021-06-18 | 中国人民解放军总医院 | 检测入睡点的方法 |
CN113616163A (zh) * | 2021-10-14 | 2021-11-09 | 慕思健康睡眠股份有限公司 | 睡眠监测电路和睡眠监测系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006320734A (ja) * | 2006-06-30 | 2006-11-30 | Matsushita Electric Works Ltd | 睡眠検査装置および睡眠時無呼吸検査装置 |
CN101474066A (zh) * | 2009-01-07 | 2009-07-08 | 上海大学 | 带有加速度传感器的血压实时监测远程适时服务的方法与系统 |
-
2011
- 2011-06-24 CN CN 201110174039 patent/CN102302361B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006320734A (ja) * | 2006-06-30 | 2006-11-30 | Matsushita Electric Works Ltd | 睡眠検査装置および睡眠時無呼吸検査装置 |
CN101474066A (zh) * | 2009-01-07 | 2009-07-08 | 上海大学 | 带有加速度传感器的血压实时监测远程适时服务的方法与系统 |
Non-Patent Citations (2)
Title |
---|
王子洪: "新型穿戴式人体状态监测技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》, no. 4, 15 April 2011 (2011-04-15), pages 140 - 630 * |
秦晓华: "一种老年人移动健康监护系统的研究", 《中国医学物理学杂志》, vol. 28, no. 1, 31 January 2011 (2011-01-31), pages 2407 - 2410 * |
Cited By (35)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102662686B (zh) * | 2012-03-09 | 2015-05-13 | 中国科学院微电子研究所 | Fpga的加载方法及其装置 |
CN102662686A (zh) * | 2012-03-09 | 2012-09-12 | 中国科学院微电子研究所 | Fpga的加载方法及其装置 |
CN102743154A (zh) * | 2012-07-02 | 2012-10-24 | 深圳市长江力伟股份有限公司 | 一种便携式人体体征监视系统和装置 |
CN103837696A (zh) * | 2012-11-21 | 2014-06-04 | 通用汽车环球科技运作有限责任公司 | 传感器装置和通过陷波滤波器进行信号处理的方法 |
WO2014086261A1 (zh) * | 2012-12-07 | 2014-06-12 | 河南华南医电科技有限公司 | 一种心脏电刺激信号分配电路 |
CN103605304A (zh) * | 2013-11-09 | 2014-02-26 | 电子科技大学 | 一种多路微弱信号采集处理电路 |
CN103605304B (zh) * | 2013-11-09 | 2015-11-18 | 电子科技大学 | 一种多路微弱信号采集处理电路 |
WO2016019894A1 (en) * | 2014-08-07 | 2016-02-11 | Beijing Zhigu Tech Co., Ltd. | Dominant limb identification method and device |
CN104510460A (zh) * | 2014-12-22 | 2015-04-15 | 哈尔滨师范大学 | 一种基于物联网的独居老人自动监测系统 |
CN104490401A (zh) * | 2015-01-15 | 2015-04-08 | 张�杰 | 基于加速度数据的智能终端进行人体行为自动识别的方法 |
CN104678745A (zh) * | 2015-02-10 | 2015-06-03 | 柳州市金旭节能科技有限公司 | 一种具有跌倒检测功能的手表 |
CN106031636A (zh) * | 2015-03-11 | 2016-10-19 | 济宁智城中科信息技术有限公司 | 一种跌倒检测方法及相应的可穿戴监护装置及监护方法 |
CN105223424A (zh) * | 2015-11-13 | 2016-01-06 | 张希梅 | 变电站在线监测装置 |
CN105467224A (zh) * | 2015-11-13 | 2016-04-06 | 张希梅 | 变电站在线监测装置 |
CN105982659A (zh) * | 2016-03-17 | 2016-10-05 | 厦门鑫奥力电器有限公司 | 一种用于测量心率的耳夹及其控制方法 |
CN106344006A (zh) * | 2016-11-03 | 2017-01-25 | 太原理工大学 | 基于极点对称模态分解和支持向量机的j波检测方法 |
CN106725394A (zh) * | 2016-12-30 | 2017-05-31 | 北京工业大学 | 一种血压数据采集装置及方法 |
CN106725897A (zh) * | 2017-01-24 | 2017-05-31 | 新博医疗技术有限公司 | 一种呼吸跟踪监测通讯系统 |
CN110650685A (zh) * | 2017-03-24 | 2020-01-03 | 爱尔西斯有限责任公司 | 评估人的心理生理状态的方法 |
CN110650685B (zh) * | 2017-03-24 | 2024-02-20 | 爱尔西斯有限责任公司 | 评估人的心理生理状态的方法 |
CN107088062A (zh) * | 2017-05-08 | 2017-08-25 | 歌尔科技有限公司 | 一种心率模组、采集心率的电子装置和心率采集方法 |
CN107088062B (zh) * | 2017-05-08 | 2023-09-01 | 歌尔科技有限公司 | 一种心率模组、采集心率的电子装置和心率采集方法 |
CN109269395A (zh) * | 2018-11-01 | 2019-01-25 | 重庆地质矿产研究院 | 一种岩土工程结构远程交互式应变监测装置及监测方法 |
CN109635355B (zh) * | 2018-11-19 | 2023-04-18 | 北京时代民芯科技有限公司 | 一种面向gpio的频率可调的滤波电路 |
CN109635355A (zh) * | 2018-11-19 | 2019-04-16 | 北京时代民芯科技有限公司 | 一种面向gpio的频率可调的滤波电路 |
CN109350050A (zh) * | 2018-11-20 | 2019-02-19 | 锐捷网络股份有限公司 | 心跳间期检测方法及装置 |
CN109387279A (zh) * | 2018-12-19 | 2019-02-26 | 南京工程学院 | 一种对汽轮机发电机末级的振动信号的检测电路 |
CN109387279B (zh) * | 2018-12-19 | 2020-10-16 | 南京工程学院 | 一种对汽轮机发电机末级振动信号进行检测的检测电路 |
CN110031653A (zh) * | 2019-05-09 | 2019-07-19 | 保定开拓精密仪器制造有限责任公司 | 三轴石英挠性加速度计小型信号采样装置及其加工工艺 |
CN110031653B (zh) * | 2019-05-09 | 2021-07-27 | 保定开拓精密仪器制造有限责任公司 | 三轴石英挠性加速度计小型信号采样装置及其加工工艺 |
CN112870555A (zh) * | 2021-02-07 | 2021-06-01 | 丹源医学科技(杭州)有限公司 | 一种休息起搏频率的自动调节装置和心脏起搏装置 |
CN112971720B (zh) * | 2021-02-07 | 2023-02-03 | 中国人民解放军总医院 | 检测入睡点的方法 |
CN112870555B (zh) * | 2021-02-07 | 2021-11-09 | 丹源医学科技(杭州)有限公司 | 一种休息起搏频率的自动调节装置和心脏起搏装置 |
CN112971720A (zh) * | 2021-02-07 | 2021-06-18 | 中国人民解放军总医院 | 检测入睡点的方法 |
CN113616163A (zh) * | 2021-10-14 | 2021-11-09 | 慕思健康睡眠股份有限公司 | 睡眠监测电路和睡眠监测系统 |
Also Published As
Publication number | Publication date |
---|---|
CN102302361B (zh) | 2013-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102302361B (zh) | 一种心率及身体状态监测装置及方法 | |
CN106971059B (zh) | 一种基于神经网络自适应健康监测的可穿戴设备 | |
CN101357062B (zh) | 基于容积脉搏波检测的生命体征参数测评装置 | |
CN102138789B (zh) | 一种动态心电和运动记录与分析系统 | |
CN105769173B (zh) | 一种具有心电信号去噪功能的心电监测系统 | |
EP3692900A1 (en) | Artificial intelligence self-learning-based static electrocardiography analysis method and apparatus | |
CN101732040B (zh) | 一种无创多路脉搏波检测方法 | |
CN100515327C (zh) | 一种心血管系统联合变异性指标的检测方法和装置 | |
CN101721206B (zh) | 具有人体运动状态及姿态检测功能的动态心电图装置 | |
CN104523263A (zh) | 基于移动互联网的孕产妇健康监护系统 | |
CN111084618A (zh) | 一种可穿戴式多功能呼吸循环检测系统及方法 | |
US20060287605A1 (en) | Heart rate variability analyzing device | |
CN203662751U (zh) | 一种便携式心电仪 | |
CN103027675A (zh) | 新型便携式三导联实时无线心电监测系统及分析方法 | |
CN202027572U (zh) | 一种动态心电和运动记录与分析装置 | |
CN102462494A (zh) | 一种新型智能心电测试保健器 | |
CN101612044A (zh) | 穿戴式人体手臂震颤检测装置及其信号处理方法 | |
CN110148466A (zh) | 一种基于迁移学习的心冲击信号房颤计算机辅助诊断方法 | |
CN202920160U (zh) | 中医脉象采集系统 | |
CN202960481U (zh) | 中医脉象采集装置 | |
CN201227272Y (zh) | 基于Zigbee的穿戴式心电和呼吸率监护仪 | |
CN113160921A (zh) | 一种基于血流动力学的数字人体心血管系统的构建方法和应用 | |
CN101297753A (zh) | 基于无线传感器网络的标准三导联心电监护系统 | |
CN103211591A (zh) | 心电信号采集盒 | |
CN113288168A (zh) | 一种智能可穿戴疲劳监测及预警系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130213 Termination date: 20180624 |