心率变换信号分析.pdfVIP

振 动 与 冲 击 第26卷第 1O期 JOURNAL OFVIBRATIONAND SHOCK 基于 Hilbert-Huang变换的心率变异信号分析 宋立新 ， 王玉静 ， 王 祁 (1．哈尔滨理工大学电气与电子工程学院，哈尔滨 150040；2．哈尔滨工业大学 自动化测试与控制系，哈尔滨 150001) 摘 要 ：提出一种基于Hilbert—Huang变换(HHT)的心率变异信号分析的新方法。心率变异分析被广泛应用于评 估心脏 自律功能以及疾病诊断领域。为获得更多心率信号内在特征，首先利用经验模态分解 (EmpiricalModeDecomposi— tion，EMD)方法将信号分解为一组固有模态函数(IntrinsicModeFunction，IMF)，运用 Hilbert变换计算并分析各层 IMF的 瞬时频率和瞬时幅值 ，从而获取信号所包含的内在信息、心率变异突发时刻和变化趋势。结合积分脉冲频率调制 (Inte— gralPulseFrequencyModulation，IPFM)模型模拟出的心率信号以及真实的心率信号，利用小波分析方法以及HHT方法对 心率信号进行对比分析，实验结果证明了HHT方法的可行性 ，显示了该方法相对于小波分析方法的优势。 关键词：Hilbert—Huang变换；经验模态分解；IPFM模型；心率变异 中图分类号：TN911．72 文献标识码：A 心电信息是人体最重要的生物电信息之一，它反 得的IMF必须满足下列 2个条件 ： 映了心脏活动的基本规律，在提高心血管疾病诊断效 1)其极点和过零点的数 目应该相等或至多差 率和准确性方面有着十分重要的意义。心率变异分析 一 点； 近年来也越来越多地被医学研究者所关注。小波分析 2)在任意时刻，由局部极大值点定义的上包络线 方法因其具有多分辨率分析的特点 J，在时频两域都 和极小值点定义的下包络线的平均值为零。各层 IMF 具有表征信号局部特征的能力而成为非平稳信号分析 的筛选过程如下： 的有力工具，但是对信号作小波变换时小波基的选择 根据信号 (t)的局部极大值与局部极小值确定上 问题以及小波分析方法的非 自适应性也使其具有一定 下包络线的平均值 a(t)，得 (t)与 a(t)的差 h 的局限性 J。 h= (t)一a(t) (1) Hilbe~一Huang变换 (HHT) 是 NordenE．Huang 将h视为新的 (t)，重复以上操作直到 h满足 等人在 1998年提出的一种用于非线性和非平稳时间 IMF所需条件 ，记为 序列分析的新方法。该方法具有 自适应提取信号特征 C =h (2) 的能力 J，将其应用到心率信号上，不仅可以获得心 C视为第一层 IMF，作 率信号所包含的主要频率成分及幅值特征，还能够准 X(t)一C =r (3) 确地检测到信号各频率成分和幅值的瞬态变化情况以 将 r视为新的 (t)，重复以上过程，依次得到第二 及突发状况，方便医学工作者了解病人心脏活动的瞬 层 IMFC：，第三层 IMFC……，直到c或 r满足给定的 时生物现象。 终止条件 (如果分解出的IMF或残余函数 r足够小，或 本文基于Hilbert—Huang变换方法，结合积分脉冲 r成为单调函数)为止。于是 (t)可分解为 频率调制(IPFM)模型_6 模拟出的两组心跳发生时间 序列和一组真实心率数据进行科学分析，体现出了该