一种非接触式生命体征检测系统的心跳频率检测算法
【技术领域】
[0001] 本发明涉及一种基于雷达技术的非接触式生命体征检测系统的信号处理算法，尤 其是心跳频率检测算法。
【背景技术】
[0002] 非接触式生命特征检测技术将雷达技术和生物医学工程技术融为一体，利用电磁 波对人体的呼吸和心跳等生命特征进行非接触式探测。电磁波具有非金属介质穿透特性， 因此非接触式生命体征监测系统能够穿透衣服、被褥检测生命信号，不需要与人体皮肤直 接接触，这一优势使其在婴儿、有睡眠障碍的病人、大面积烧伤病人、精神病人等特殊群体 的生命体征监测中得到应用。此外，经过特殊设计的生物雷达还能穿透砖墙、木板等障碍物 探测到人体的微动信号，这使得其在救灾、反恐、军事等领域也有广泛的应用前景。
[0003] 心跳频率检测的准确度是非接触式生命特征检测系统的最重要的指标之一，由于 心跳信号能量很微弱且具有非平稳特性，所以心跳信号很容易淹没在雷达的噪声和杂波 中。此外，由于人体在监测过程中很有可能进行随机运动，如人体的抖动、翻身等，这些随机 运动产生的扰动干扰信号都具有非平稳特点，其频率很有可能在心跳频率范围内，因此传 统滤波方法不能将这些扰动干扰滤除，而且由于传统FFT方法集中度较低、分辨力较差，因 此采用FFT测频法不能准确测得心跳频率。
【发明内容】
[0004] 为了克服上述现有技术的不足，本发明提供了一种非接触式生命体征检测系统的 心跳频率检测算法，该算法能将心跳信号与人体抖动信号、呼吸信号和环境干扰噪声分离， 而且能准确提取心跳信号的频率信息。
[0005] 本发明采用的技术方案如下：
[0006] 一种用于非接触式生命体征检测系统的心跳频率检测算法，包括如下步骤：
[0007] 步骤1，对非接触式生命体征检测系统的连续波多普勒雷达前端输出的I、Q两路 信号分别进行带通滤波；其中，带通滤波器的低截止频率为〇. 5Hz、高截止频率为5. 2Hz ;
[0008] 步骤2,通过圆心估计算法恢复经过滤波的I、Q两路信号的有用直流分量；
[0009] 步骤3,对步骤2得到的I、Q两路信号进行反正切解调；
[0010] 其特征在于，得到解调信号之后的步骤如下：
[0011] 步骤4,采用完全聚类经验模态分解算法对解调信号进行分解，得到一系列本征模 函数；
[0012] 步骤5,对得到的本征模函数分别做希尔伯特变换，并且求得各本征模函数对应的 希尔伯特边际谱；
[0013] 步骤6,在各本征模函数对应的希尔伯特边际谱上进行峰值检测；
[0014] 步骤7,根据希尔伯特边际谱的谱峰位置和谱峰附近的能量集中度来判断心跳信 号对应的希尔伯特边际谱；
[0015] 步骤8,在心跳信号对应的希尔伯特边际谱上提取心跳信号频率信息。
[0016] 进一步的，步骤4中采用完全聚类经验模态分解算法对解调信号进行分解，具体 步骤如下：
[0017] 步骤4. 1，构造 I个信号χ[η]+ ε，其中x[n]是反正切解调得到的 信号，win]，（i = 1，2，...，1)是分布为N(0，1)的白噪声，利用普通经验模态分解 算法分别求解I个信号Χ[η]+ε<^[η]的第一个本征模分量nffV[n]，并且计算：
[0018] 步骤4. 2,计算第一个剩余分量
[0019] 步骤4. 3,对
进行经验模态分解，得到第一个本 征模分量；其中，算子E,( ·)的功能是利用普通经验模态分解算法求一个信号的第j个模 态，并由此得到：
[0020] 步骤4. 4,根据k = 2, · · ·，K，计算第k个剩余分量：
[0021] 步骤4. 5,根据
进行经验模态分解，得到第一个 本征模分量，并由此得到
[0022] 步骤4. 6,重复步骤4. 4到步骤4. 5步，当剩余本征模分量不再满足可以分解的 要求时，即剩余本征模分量的极值不足两个时，停止运算，最后的剩余本征模分量表示为：
是总模态数；
[0023] 至此，解调信号X [η]被分解为K个本征模函数MFk [η]和一个剩余分量R [η]。进 一步的，步骤5中求解各本征模函数mFk[n]的希尔伯特边际谱，步骤如下：
[0024] 步骤5. 1，令uk(n) = IMFk[n]，首先对uk(n)做离散希尔伯特变换：
[0025] 步骤5. 2,计算
[0026] 步骤5. 3,计算
对上式 进行修正
[0027] 步骤5. 4,求希尔伯特谱
[0028] 步骤5. 5,求得各本征模函数頂Fk[n]对应的希尔伯特边际谱：
[0029] 进一步的，步骤7中根据谱峰位置和谱峰附近的能量集中度来判断心跳信号对应 的希尔伯特边际谱，具体方法如下：
[0030] 步骤7. 1，判断各边际谱hk(f)的谱峰对应的频率fpeak,k是否在0. 5Hz到3Hz之间；
[0031] 步骤7. 2,对于满足步骤7. 1条件的边际谱，计算其在f_k,k-0. 1Hz到f_k,k+0. 1Hz 之间各频率点的幅度的总和Ak，并计算其所有频率点的幅度总和Bk，求能量集中度Ek= A k/ Bk;
[0032] 步骤7.3,当^大于阈值α时，认为hk(f)是心跳信号的边际谱；若有多个边际谱 满足^大于阈值α的条件，则认为峰值幅度最大的边际谱是心跳信号的边际谱。
[0033] 进一步的，步骤8中将心跳信号的希尔伯特边际谱的谱峰位置对应的频率作为心 跳信号的频率。
[0034] 进一步的，所述心跳信号的希尔伯特边际谱峰值对应的是频率为1. 2402Hz。
[0035] 与现有技术相比，本发明可在人体不稳定、环境干扰噪声大的条件下有效提取心 跳信号并且获得准确心跳频率信息。
[0036] 下面结合说明书附图对本发明做进一步描述。
【附图说明】
[0037] 图1为本发明的非接触式生命体征检测系统的心跳频率检测算法流程图。图中的 输入信号是连续波多普勒雷达前端送来I、Q两路信号。
[0038] 图2为对I、Q两路信号进行带通滤波、圆心估计以及反正切解调后得到的解调信 号。此信号中包含心跳信号以及干扰信号，其中13s之间23s有比较大的人体抖动干扰。
[0039] 图3为对解调信号进行完全聚类经验模态分解后得到的一系列本征模函数的波 形。
[0040] 图4为第1到6个本征模函数对应的希尔伯特边际谱。（a)、本征模函数1的边际 谱，（b)、本征模函数2的边际谱，（c)、本征模函数3的边际谱，（d)、本征模函数4的边际谱， (e)、本征模函数5的边际谱，（f)、本征模函数6的边际谱。
【具体实施方式】
[0041] 非接触式生命体征检测系统是基于连续波多普勒雷达探测原理，雷达天线向人体 胸腔位置发射电磁波，呼吸和心跳等引起胸腔的起伏会对电磁波的相位进行调制，胸腔反 射的电磁波被雷达天线接收后经过雷达前端的处理变为I、Q两路信号输出，对I、Q两路信 号进行AD采样后，进行如下处理：
[0042] 采用本发明提供的非接触式生命体征检测系统的心跳频率检测算法对实验采集 到的信号进行处理：
[0043] 步骤1，对非接触式生命体征检测系统的连续波多普勒雷达前端输出的I、Q两路 信号分别进行带通滤波，带通滤波器的低截止频率为〇. 5Hz、高截止频率为5. 2Hz。
[0044] 步骤2,通过圆心估计算法恢复经过滤波的I、Q两路信号的有用直流分量。
[0045] 步骤3,对步骤2得到的两路信号进行反正切解调。图2所示为实验实例中得到的 解调信号波形，由图可见，此信号中包含心跳信号以及干扰信号，其中13s之间23s有比较 大的人体抖动干扰。
[0046] 步骤4,采用完全聚类经验模态分解算法对解调信号进行分解，得到一系列本征模 函数分量，具体计算步骤如下：
[0047] 步骤4. 1，构造 I个信号
其中x[n]是反正切解调得到的 信号，
是分布为N(0，1)的白噪声。利用普通经验模态分 解算法分别求解I个信号
的第一个本征模分量并且计算
[0048] 步骤4. 2,计算第一个剩余分量
[0049] 步骤= 1，···，Ι进行经验模态分解（算子EJ·) 的功能是求一个信号经验模态分解的第j个模态），得到第一个本征模分量，并由此得到
[0050] 步骤4. 4,对于k = 2, · · ·，K计算第k个剩余分量