心率识别程序求解释
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////心率计算模块
-*--------------------------------------------------------------*/
#include "main.h"
#include "rhythmeheart.h"
/*************************************************
Function: rhythmcount
Description:计算心率值
Calls: --
Called By:real_sample()
函数入口:数据区首地址
Return: rhythm 心率值
*************************************************/
unsigned char rhythmcount(unsigned char *fifo)
{
signed char i;
signed char diff; //求导数数据区
signed char dif_max=0; //导数最大值
signed char thr_max; //阈值最大值
unsigned char j=0;
signed char k; //最大值数组
unsigned char rhythm;
float DR_R; //R-R间距
*diff=0; //数据区首两个数据导数置0
*(diff+1)=0;
for(i=2;i<max_data-2;i++) //按DIFF(i)=f(i+1)-f(i-1)+2*f(i+2)-2*f(i-2)公式计算导数
{
*(diff+i)=(signed char)(*(fifo+i+1))+2*(signed char)(*(fifo+i+2))-(signed char)(*(fifo+i-1))-2*(signed char)(*(fifo+i-2));
}
*(diff+max_data-1)=0; //数据区末两个数据导数置0
*(diff+max_data)=0;
for(i=0;i<max_data;i++) //求导数最大值
{
if((dif_max)>*(diff+i))
dif_max=*(diff+i);
}
thr_max=(dif_max>>1)-(dif_max>>3); //设定阈值,为导数最大值的0.375倍
for(i=1;i<max_data-1;i++) //将满足f(i)>thr_max和f(i)*f(i+1)<0这两个条件的值找到,即为R波最高点
{
if((*(diff+i)>thr_max)&&(((*(diff+i))^(*(diff+i+1)))>>7))
k=i;
}
DR_R=(k-k)/(j-1); //计算R-R间距
rhythm=(unsigned char)(12000/DR_R); //通过公式200*60/DR-R求得心率值
return rhythm;
}
/////////////////////////////////////////////////////////////////////////////////////////////////
心率计算模块头文件
-*--------------------------------------------------------------*/
#ifndef _RHYTHMEHEART_H
#define _RHYTHMEHEART_H
#define max_data 512 //数组最大值,调试时可往下修正
#define default_rhythm 75 //初始心率值
unsigned char rhythmcount(unsigned char *fifo);
#endif
/*-----------------------------------------------------------------*-
--------------------END OF FILE------------------------------------
-*-----------------------------------------------------------------*/
/////////////////////////////////////////////////////////////////////////////////////////////////// 这个还是自己理解吧 请问,这个心率模块是什么型号的 来看看 mark留个记号 没看懂 DIFF(i)=f(i+1)-f(i-1)+2*f(i+2)-2*f(i-2) 怎么算出导数的(高数重修N次才过关的)。
不过整段代码的大意是找一段数据里的极大值(波峰顶点),然后和上一个极大值所在的时间点相减得出两个波峰之间的时间差,就是RR了(当然,他代码里是把几个波峰的总时间差做平均)。然后最后返回值哪里又不懂了,60是60秒的意思,200比较纳闷,估计采样频率是200HZ吧。
所以关键不在这段代码上,而是在传入的FIFO是串啥数据,是怎么处理后得到的。
感觉不是胸率带(只检测R峰)就是脉博测量的。 wingerchen 发表于 2014-10-24 20:02
没看懂 DIFF(i)=f(i+1)-f(i-1)+2*f(i+2)-2*f(i-2) 怎么算出导数的(高数重修N次才过关的)。
不过整段代码 ...
不是,它不仅是检测极大值,还有斜率。差分阀值法 请问,这个心率模块是啥名字,我也想买来学习一下 maimaige 发表于 2014-10-24 23:22
请问,这个心率模块是啥名字,我也想买来学习一下
不知道,这代码是csdn上,你可以去看看 心率模块,这也有啊 求导是计算斜率的,就是通过找电压变化最快的点来辅助定位R波位置。有时T波的幅度会大于QRS波群,但斜率一定是Q-->R最大。所以要先用斜率找到R波的上升沿,然后再搜索幅度最大值处的R波位置。 但是这种算法并不可取,因为每次需要计算512组心电数据 好专业,过来学习~ 算法了解不多,过来学习学习! 学习学习~ 心率测试现在很火啊 我这边有心率测量模块 不过没研究过 这方面的算法 Name_006 发表于 2014-10-30 13:08
我这边有心率测量模块 不过没研究过 这方面的算法 ...
啊,你那模块还能自己计算心率吗 机器人天空 发表于 2014-10-30 13:11
啊,你那模块还能自己计算心率吗
模块上 有个运放做的电压比较器心率波形 输出电压高出 设定的电压值 比较器后边DO就会输出高电平 低于设定电压值 就输出低电平, 这样就转化成了 脉冲输出 可以用单片机来计数 Name_006 发表于 2014-10-30 13:31
模块上 有个运放做的电压比较器心率波形 输出电压高出 设定的电压值 比较器后边DO就会输出高电平 低于 ...
懂了,以前的监护仪测心率也是的这样的
页:
[1]