基于二次傅里叶变换的雷达心率监测方法和系统与流程

本发明涉及信号处理和健康监测，尤其是一种基于二次傅里叶变换的雷达心率监测方法和系统。

背景技术：

1、率是人体生命体征的基本指标，也是人体健康状况的重要参考。相比较传统的接触式测量方法(如ecg、红外传感器、振动传感器等)，利用雷达进行无接触心率测量具有使用舒适度高、实时性高等特点，十分适合家庭日常健康监测。

2、从信号角度分析：毫米波雷达通过测量心脏跳动引起体表位移实现心率测量，一般情况下，人体胸腹部有呼吸和心跳两种振动：通常，呼吸引起的胸腹部位移为1-12mm，心跳引起的胸腹部位移为0.1-0.5mm，故呼吸位移一般远大于心跳，呼吸对应的回波信号相位变化也远大于心跳，在利用相位信号进行呼吸率心率分析时，呼吸对应的相位能量远强于心跳；一般情况下，人体呼吸频率一般在0.6hz以下，心跳频率一般在0.8hz至2.5hz之间，而人体的呼吸模式十分复杂，可能存在如类似正弦式、三角形式、矩形式、高斯形式等模式，甚至存在以上两种或多种波形的组合模式，由此可能会产生复杂多变且功率较强的谐波，且该谐波可能会落入心率区间，从而影响对心率的精确测量。

3、从心跳的振动模式分析：心跳振动来自血液对心房壁/心室壁/动脉血管壁等冲击所致，该过程短促而有力，从频谱上看，该过程存在丰富的谐波，而通常检测心跳频率往往关注于心跳基频，而忽视了丰富的谐波信息，实际上，目前的研究缺少合理的手段对心跳谐波进行利用，从而使心率测量更加鲁棒且准确。

4、基于毫米波雷达心率监测的重要技术挑战之一就是，存在呼吸干扰时心率的准确测量。目前常见的处理方法主要有带通滤波、最优频率估计、机器学习等。这些方法能够抑制心率区间之外的直流、呼吸谐波等干扰量，但是对于落在心率区间内的呼吸谐波或调制分量则无法消除，从而影响心率的准确估计。

技术实现思路

1、为了克服上述现有技术中呼吸谐波及调制引起的心率监测干扰问题，本发明提出了一种基于二次傅里叶变换的雷达心率监测方法，提高心率计算的鲁棒性和准确性，且复杂度较低。

2、本发明提出的一种基于二次傅里叶变换的雷达心率监测方法，包括：

3、配置雷达，通过雷达在采样时长内结合采样频率采集胸腔监测数据，并整理为i/q正交数据；对i/q正交数据进行解调得到i/q解调信号；

4、对i/q解调信号中各距离单元的相位数据进行傅里叶变换，获得各距离单元的频谱信息f(i)，i为距离单元序号；

5、对各距离单元i进行频谱处理，令距离单元上的频谱信息f(i)在指定频率区间外置零以形成频谱信息f(i,1)；然后对频谱信息f(i,1)进行滑动窗口归一化，以形成频谱信息f(i,2)；对频谱信息f(i,2)进行傅里叶变换，获取频谱信息f’(i)；

6、确定频谱信息f’(i)的峰值的信噪比(snr)，选择信噪比最大的频谱信息f’(i)作为目标频谱，以目标频谱上两个谱峰之间的频率差值作为心率。

7、优选的，频谱信息f(i,2)的获取方式为：令频谱信息f(i,1)中频率点j的幅值记作u(j)，频谱信息f(i,2)中频率点j的幅值记作u’(j)，频谱信息f(i,1)中频率范围[j，j+f]上的幅度均值记作a_u(j)，频谱信息f(i,1)中频率范围[j，j+f]上的最大幅值记作max_u(j，j+f)；f为设定步长；u’(j)＝a_u(j)/max_u(j，j+f)。

8、优选的，设定步长f的取值范围为[0.1，5]hz。

9、优选的，选择信噪比最大且大于设定信噪比阈值的频谱信息f’(i)作为目标频谱。

10、优选的，采样时长的取值范围为[3s,180s]。

11、优选的，采样频率的取值范围为每秒10-1000次。

12、本发明提出的一种基于二次傅里叶变换的雷达心率监测系统，包括存储器和处理器，存储器中存储有计算机程序，处理器连接存储器，处理器用于执行所述计算机程序，以实现所述的基于二次傅里叶变换的雷达心率监测方法。

13、优选的，还包括雷达，处理器与雷达连接，以获取雷达采集的胸腔监测数据。

14、优选的，雷达采用毫米波雷达。

15、本发明提出的一种存储介质，存储有计算机程序，所述计算机程序被执行时用于实现所述的基于二次傅里叶变换的雷达心率监测方法。

16、本发明的优点在于：

17、(1)本发明提出的基于二次傅里叶变换的雷达心率监测方法，通过对积累的雷达距离数据解析相位，得到各距离单元的相位信息，通过第一次傅里叶变换得到相位变化的频谱，通过对频谱进行频谱带通滤波、滑动归一等操作，得到周期特性良好的心跳频率谐波分布，将此分布进行二次傅里叶变换，可计算心跳频率谐波的周期，该周期即为心跳周期频率。

18、(2)本发明中，通过计算二次傅里叶变换的谱峰的位置和信噪比，可以判断心跳频率，最后综合所有距离单元的计算结果可以判定心率。该方法一方面可以大大减小呼吸干扰对心率测量的影响，另一方面可以极好的利用心跳振动丰富的谐波信息，从而具有较好的鲁棒性；

19、(3)本发明主要处理过程都在频域进行，包括带通滤波、滑动归一、计算周期等，可以省去复杂且耗时的时域滤波、自相关等算法，具有复杂度低、计算量小、适合实时处理等优势，十分适合日常健康监测。

20、(4)本发明能准确且快速地计算出人体的心率，并且能较好地抵抗呼吸等干扰，十分适合心率检测。

技术特征：

1.一种基于二次傅里叶变换的雷达心率监测方法，其特征在于，包括：

2.如权利要求1所述的基于二次傅里叶变换的雷达心率监测方法，其特征在于，频谱信息f(i,2)的获取方式为：令频谱信息f(i,1)中频率点j的幅值记作u(j)，频谱信息f(i,2)中频率点j的幅值记作u’(j)，频谱信息f(i,1)中频率范围[j，j+f]上的幅度均值记作a_u(j)，频谱信息f(i,1)中频率范围[j，j+f]上的最大幅值记作max_u(j，j+f)；f为设定步长；u’(j)＝a_u(j)/max_u(j，j+f)。

3.如权利要求2所述的基于二次傅里叶变换的雷达心率监测方法，其特征在于，设定步长f的取值范围为[0.1，5]hz。

4.如权利要求1所述的基于二次傅里叶变换的雷达心率监测方法，其特征在于，选择信噪比最大且大于设定信噪比阈值的频谱信息f’(i)作为目标频谱。

5.如权利要求1所述的基于二次傅里叶变换的雷达心率监测方法，其特征在于，采样时长的取值范围为[3s,180s]。

6.如权利要求1所述的基于二次傅里叶变换的雷达心率监测方法，其特征在于，采样频率的取值范围为每秒10-1000次。

7.一种基于二次傅里叶变换的雷达心率监测系统，其特征在于，包括存储器和处理器，存储器中存储有计算机程序，处理器连接存储器，处理器用于执行所述计算机程序，以实现如权利要求1-6任一项所述的基于二次傅里叶变换的雷达心率监测方法。

8.如权利要求7所述的基于二次傅里叶变换的雷达心率监测系统，其特征在于，还包括雷达，处理器与雷达连接，以获取雷达采集的胸腔监测数据。

9.如权利要求7所述的基于二次傅里叶变换的雷达心率监测系统，其特征在于，雷达采用毫米波雷达。

10.一种存储介质，其特征在于，存储有计算机程序，所述计算机程序被执行时用于实现如权利要求1-6任一项所述的基于二次傅里叶变换的雷达心率监测方法。

技术总结本发明涉及信号处理和健康监测技术领域，尤其是一种基于二次傅里叶变换的雷达心率监测方法和系统。本发明提出的基于二次傅里叶变换的雷达心率监测方法，通过对积累的雷达距离数据解析相位，得到各距离单元的相位信息，通过第一次傅里叶变换得到相位变化的频谱，通过对频谱进行频谱带通滤波、滑动归一等操作，得到周期特性良好的心跳频率谐波分布，将此分布进行二次傅里叶变换，可计算心跳频率谐波的周期，该周期即为心跳周期频率。本发明方法一方面可以大大减小呼吸干扰对心率测量的影响，另一方面可以极好的利用心跳振动丰富的谐波信息，从而具有较好的鲁棒性。技术研发人员：魏正武,万锦伟受保护的技术使用者：数据空间研究院技术研发日：技术公布日：2024/9/29