[发明专利]一种基于瞬时声音刺激的心脏功能评价方法

权利要求书

1.一种基于瞬时声音刺激的心脑耦合指标计算方法，其特征在于，包括以下步骤：

步骤1)、进行多次听觉刺激实验并记录刺激前及刺激过程中被测试者的心电信号；

步骤2)、对获取的心电信号进行预处理得到逐拍RR间期RR(n)：对采集到的原始心电信号ecg(t)，选择db2小波对信号进行9层分解，去除最后一层的基线漂移以及前三层的高频噪声，得到去噪后的心电信号x(t)；利用小波模极大值的方法对R波峰进行检测，采用支持紧支集且具有一阶消失矩的db2小波，对R波峰进行检测，db2小波的等效滤波器系数如下：

h1＝0.3750,h2＝0.1250,h3＝0.0000

g1＝0.5798.g2＝0.0869,g3＝0.0061

h_k＝h_1-k,g_k＝-g_1-k

若k3,h_k＝g_k＝0

对所采集到的心电信号用Mallat算法作小波分解；

再对检测结果进行检视修改，获得一连串连续R波峰时间点的集合R(N)；对R波峰时间点的集合R(N)，利用后一个波峰时间点R(n+1)减去前一个波峰的时间点R(n)，获得逐拍RR间期RR(n)，其中(1≤nN)；

步骤3)、对步骤2)得到的RR间期RR(n)进行RR间期响应序列计算，得到归一化的RR间期序列RR’(t’,i)；具体的，对得到的RR间期序列RR(t’,i)进行归一化，得到归一化的RR间期序列RR’(t’,i)：

其中，RR_baseline为RR间期序列在-2.5t0段的平均值，将RR′(t′,i)沿着刺激次数求平均，得到被测试者的RR间期声音刺激响应序列RR′(t′)；

步骤4)、根据得到的归一化的RR间期序列RR’(t’,i)，计算被测试者在刺激后内RR间期平均增加量，即心脑耦合指标ΔRRI_in，计算方法如下：

ΔRRI_in＝(RR′(0.5)+RR′(1))/2。

2.根据权利要求1所述的一种基于瞬时声音刺激的心脑耦合指标计算方法，其特征在于，记录声音刺激的时间点t(i)，其中i＝1,2,…,20，将步骤2)中所得的每个刺激时间点t(i)附近的RR间期序列分别进行从t(i)时间点前后等时间s内的2Hz的重采样，得到采样率统一的RR间期序列RR(t,i)，并将新得到RR间期序列的刺激出现点定义为时间点0，数据起始时间点为-s，结束时间点为s，即RR(t’,i)。

3.根据权利要求2所述的一种基于瞬时声音刺激的心脑耦合指标计算方法，其特征在于，具体的，分别进行从t(i)-2.5s至t(i)+2.5s的2Hz的重采样，得到采样率统一的RR间期序列RR(t,i)，并将新得到RR间期序列的刺激出现点定义为时间点0，数据起始时间点为-2.5，结束时间点为2.5，即RR(t’,i)。