心率识别程序求解释

2020-01-21 21:47发布

///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
心率计算模块
-*--------------------------------------------------------------*/
#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[max_data];             //求导数数据区
        signed char dif_max=0;                  //导数最大值
        signed char thr_max;                    //阈值最大值
        unsigned char j=0;
        signed char k[8];                       //最大值数组
        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[j++]=i;
        }
        DR_R=(k[j-1]-k[0])/(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------------------------------------
-*-----------------------------------------------------------------*/
///////////////////////////////////////////////////////////////////////////////////////////////////
友情提示: 此问题已得到解决,问题已经关闭,关闭后问题禁止继续编辑,回答。
该问题目前已经被作者或者管理员关闭, 无法添加新回复
19条回答
maimaige
1楼-- · 2020-01-22 14:53
 精彩回答 2  元偷偷看……
机器人天空
2楼-- · 2020-01-22 16:27
maimaige 发表于 2014-10-24 23:22
请问,这个心率模块是啥名字,我也想买来学习一下

不知道,这代码是csdn上,你可以去看看
guoliezhi
3楼-- · 2020-01-22 18:50
心率模块,这也有啊
机器人天空
4楼-- · 2020-01-23 00:12
求导是计算斜率的,就是通过找电压变化最快的点来辅助定位R波位置。有时T波的幅度会大于QRS波群,但斜率一定是Q-->R最大。所以要先用斜率找到R波的上升沿,然后再搜索幅度最大值处的R波位置。
机器人天空
5楼-- · 2020-01-23 02:57
但是这种算法并不可取,因为每次需要计算512组心电数据
nnmnnm
6楼-- · 2020-01-23 06:03
 精彩回答 2  元偷偷看……

一周热门 更多>