如何从 HR（心率）获取 HRV（心率变异性）

Question

我正在开发一个应用程序，我必须读取直到一分钟的心跳并从中计算心率变异性。

为了阅读心跳，我使用的是表带（Polar）

Heart beat reading code I have used from link

为了计算 HRV，我浏览了以下链接，但没有任何帮助：

1. Research article on SDNN/RMSSD

2. RMSSD for HRV Analysis

3. Comparison between HR and HRV

请提供我可以从 HR 中得到 HRV(RMSSD) 的公式，得到 HR 的持续时间是一分钟。

任何帮助将不胜感激..

编辑：

我通过以下代码获得了 RR 值：

- (void) updateWithHRMData:(NSData *)datas {

const uint8_t *reportData = [datas bytes];

uint16_t bpm = 0;
uint16_t bpm2 = 0;

if ((reportData[0] & 0x04) == 0)
{
    NSLog(@"%@", @"Data are not present");
}
else
{

    bpm = CFSwapInt16LittleToHost(*(uint16_t *)(&reportData[2]));

    bpm2 = CFSwapInt16LittleToHost(*(uint16_t *)(&reportData[4]));

    if (bpm != 0 || bpm2 != 0) {

            NSLog(@"%u", bpm);

            if (bpm2 != 0) {
                NSLog(@"%u", bpm2);
            }
    }

}

}

我的问题是我得到的 RR 值如下：666,636,645 .... 等 但是当我使用HRV + 应用程序并通过电子邮件导出 RR 值时，它会显示 0.785、0.734、0.724 等值。如果我使用以下公式计算我的 RR 值：

RMSSD =

它给了我完全错误的结果。

请帮忙。

编辑：

 //==================RR
if ((reportData[0] & 0x04) == 0)
{
    NSLog(@"%@", @"Data are not present");
}
else
{
    uint16_t rr2 = CFSwapInt16LittleToHost(*(uint16_t *)(&reportData[2]));
    RR = CFSwapInt16LittleToHost(*(uint16_t *)(&reportData[4]));

    //sometimes two values of RR may found so there are two RR: RR and rr2 which added in array and all calculation goes after one minute

    if (rr2 > 0 || RR > 0) {

        NSLog(@"RR2 %u", rr2);
        if (rr2>0) {
            [RRArray addObject:[NSNumber numberWithInt:rr2]];
        }
        if (RR != 0) {
            NSLog(@"RR %u", RR);
            [RRArray addObject:[NSNumber numberWithInt:RR]];

            }

        }

    }

}

Answer 1

您需要执行您的第一篇链接文章第 2.2.1 节中描述的统计计算。

您无法根据心率计算 HRV - 您需要使用 RR 间隔。

您需要检查传感器返回的数据以查看标志字节的第 4 位是否已设置 - 这表示存在 RR 数据。然后，您可以将 RR 数据读入一个数组，将 16 位样本数据转换为 NSNumbers 。您可以修改教程中的代码，将 16 位心跳值转换为 int。

一旦您收集了足够的样本（我不确定 1 分钟的样本在统计上是否有效，但我不是统计学家），然后您就可以执行分析了。

THB 将是数组中 RR 间隔的数量

MRR 或 I(BAR) 你可以计算如下 -

float rrTotal=0;

for (int i=1;i<[rrIntervals count]; i++) {
   rrTotal+=[[rrIntervals objectAtIndex:i] intValue];
}

float mrr = rrTotal/([rrIntervals count]-1);

SDNN（即您的 HRV）计算如下 -

float sdnnTotal = 0;

for (int i=1;i<[rrIntervals count]; i++) {
   sdnTotal+=pow( [[rrIntervals objectAtIndex:i] intValue] - mrr,2) ;
}

float sdnn = sqrt(sdnTotal/([rrIntervals count]-1));

您可以在向数组添加更多数据时不断重新计算 HRV。