一种基于毫米波雷达的呼吸心跳监测系统技术方案

本发明专利技术公开了一种基于毫米波雷达的呼吸心跳监测系统，包括FMCW毫米波雷达，采用所述FMCW毫米波雷达对人体发射电磁波信号，接收携带生命体征的反射信号，对接收的雷达回波信号进行解调，并从中分离提取出呼吸与心跳信息，数据采集以及解析：雷达接收到回波数据后，将原始回波数据通过网络电缆传输到PC端并进行解析从而获取Range

【技术实现步骤摘要】
一种基于毫米波雷达的呼吸心跳监测系统


[0001]本专利技术涉及信号处理
，尤其涉及一种基于毫米波雷达的呼吸心跳监测系统。

技术介绍

[0002]现代社会生活水平的提高，促使我国人口平均寿命不断增长，但人口老龄化的问题也日趋严重。目前，我国65岁及以上老人所占比重已经接近10％，根据报告预测，到2030年时，我国65岁以上的人口占全国总人口的比重将提高到18.2％。日常生活中老年人群体经常遭受各种慢性健康问题的困扰，这些疾病不易被察觉，且持续时间长。通常，呼吸、心跳等生命体征信号一定程度上反映了个体的健康状况，医护人员根据这些体征参数是否异常判断病人的健康状况。因此，长期对老年人的呼吸和心跳进行监测，可及时发现各种慢性疾病的先兆，达到有效预防和提前治疗的目的。
[0003]经过呼吸、心跳周期运动调制的雷达回波信号是一个非平稳信号，当前从雷达回波从提取分离生命体征信号的解决方案主要有短时傅里叶变换(Short
‑
Time Fourier Transform，STFT)、小波变换、自适应滤波和模态分解等。STFT信号处理方法是在傅里叶变化的基础上发展起来处理非平稳信号的。STFT的主要思想是把长的非平稳信号看成是一系列等时长的短时平稳信号的叠加,信号做傅里叶变换之前乘一个时间有限的窗函数，并假定非平稳信号在分析窗的短时间隔内是平稳的，通过窗函数在时间轴上的移动，对信号进行逐段分析得到信号的一组局部“频谱”。STFT过程中，窗的长度决定频谱图的时间分辨率和频率分辨率，窗长越长，截取的信号越长，信号越长，傅里叶变换后频率分辨率越高，时间分辨率越差，即STFT中，时间分辨率和频率分辨率两者不可兼得，必须根据具体需求进行取舍；小波变换可以较好的分离呼吸、心跳信息的信号分量，但分离的效果取决于小波基函数和小波尺度的选择。
[0004]不同的场景对于小波变换的参数要求不定相同，而且小波算法十分复杂，在代码实现上较为困难。自适应滤波的方法，按照噪声自适应抵消的原理，将实验中提取的呼吸信号作为噪声参考信号，体动信号作为原始信号，进行自适应处理，检测出心跳信号。虽然该算法结构简单，但是对噪声信号相关性要求比较高，在环境噪声影响较大且信号非平稳时，算法处理效果不理想。模态分解可以自适应将信号依据频率从高到低分离出来，可以将复合心肺信号分解为若干个子信号，从而达到心跳呼吸分离的目的。
[0005]模态分解的包括基于经验模态分解和变分模态分解方法，前者通过自适应迭代筛选分解信号，对于生命体征信号的分解效果并不理想；后者通过构造复杂的约束变分表达式以迭代求解的方式非递归地进行信号处理，这一方法算法复杂，难以在硬件设备上实现，而且存在模态混叠的问题。

技术实现思路

[0006]本专利技术公开一种基于毫米波雷达的呼吸心跳监测系统，旨在解决
技术介绍
中提出
的长期使用接触式的生命体征监测设备可能会使人们感到不舒服，比如穿戴心率带会使胸部遭到挤压导致呼吸不畅；长时间使用鼻探头这一侵入式的设备，心理上和生理上都会使人产生不适感，而这种不适有可能对呼吸产生影响的技术问题。
[0007]为了实现上述目的，本专利技术采用了如下技术方案：
[0008]一种基于毫米波雷达的呼吸心跳监测系统，包括FMCW毫米波雷达，采用所述FMCW毫米波雷达对人体发射电磁波信号，接收携带生命体征的反射信号，对接收的雷达回波信号进行解调，并从中分离提取出呼吸与心跳信息，具体包括以下步骤：
[0009]步骤一：数据采集以及解析：雷达接收到回波数据后，将原始回波数据通过网络电缆传输到PC端并进行解析从而获取Range
‑
FFT数据；
[0010]步骤二：有无目标判定：通过单帧Range
‑
FFT数据计算当前帧回波信号能量，并将其与预先设定阈值进行比较，从而判断当前帧是否包含有目标信息，若小于阈值，则当前帧无目标，若大于等于阈值，则当前帧有目标；
[0011]步骤三：目标检测：当确定当前帧包含有目标信息后，对当前帧进行目标检测，通过Range
‑
FFT数据中最大峰值点确定人体目标距离雷达位置，并统计连续10帧的最大峰值点，选取出现次数最多的最大峰值点作为人体目标的索引；
[0012]步骤四：确定人体目标的索引后，根据当前帧当前索引的点计算当前帧的相位信息，并将其循环存储到Phase_buffer中，Phase_buffer为1*1000的实数矩阵；而后通过0.2～0.5Hz，2阶带通滤波器对Phase_buffer进行滤波，取得滤波后的实数相位矩阵breath_phase_buffer，通过0.8～1.6Hz，2阶带通滤波器对Phase_buffer进行滤波，取得滤波后的实数相位矩阵heart_phase_buffer；
[0013]步骤五：对上述取得的两个实数相位矩阵进行频率估计，分别取得相位变化频率，进而转换为呼吸及心跳频率；
[0014]步骤六：对每一帧呼吸及心跳数据进行循环存储，分别统计其中出现次数最多的作为当前帧呼吸及心跳结果。
[0015]本专利技术采用已在市场上商用的宽带毫米波FMCW雷达作为测量平台，通过一系列信号处理方法，实现较高精度的实时呼吸率和心率的监测功能，适用于养老院、医院临床监护场所和其他需要呼吸、心跳监测的场所。
[0016]在一个优选的方案中，所述FMCW毫米波雷达的数据采样频率为20Hz，每帧间隔为50ms，雷达接收到目标返回中频信号后，进行FFT得到Range
‑
FFT维数据，假定该数据长度为N，则该数据为包含目标与雷达的距离信息和目标呼吸心跳引起的相位变化信息的N个复数，所述步骤二中，判断当前检测区域内是否有人体目标：由Range
‑
FFT数据计算N个复数的模值之和作为感应强度，将其与第一阈值进行比较，若其大于第一阈值，则判定当前检测区域内存在人体，若否，则输出“当前无明显生命体征”信息。
[0017]在一个优选的方案中，所述FMCW毫米波雷达监测时需要对数据进行处理，具体包括有慢多普勒维数据累积、静态杂波滤除、目标检测、目标筛选、目标点相位提取、带通滤波及保存呼吸心跳波形和频谱估计法，所述慢多普勒维数据累积：将该帧所有chirp的N点Range
‑
FFT数据对应相加，连续保存M帧数据获得M*N大小的二维矩阵，定义为Range
‑
cube；当Range
‑
cube中的数据存满以后，将最新帧的Range
‑
FFT数据存入Range
‑
cube最后位置，将Range
‑
cube第一帧数据移除，其它数据依次前移，所述静态杂波滤除：Range
‑
cube存满之
后，首先对其进行静态杂波滤除，由于静态杂波往往与雷达的距离值固定，随着时间变化，雷达回波信号在该距离点处的值基本保持不变，通过计算Range
‑
cube中M帧本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于毫米波雷达的呼吸心跳监测系统，包括FMCW毫米波雷达，其特征在于，采用所述FMCW毫米波雷达对人体发射电磁波信号，接收携带生命体征的反射信号，对接收的雷达回波信号进行解调，并从中分离提取出呼吸与心跳信息，具体包括以下步骤：步骤一：数据采集以及解析：雷达接收到回波数据后，将原始回波数据通过网络电缆传输到PC端并进行解析从而获取Range
‑
FFT数据；步骤二：有无目标判定：通过单帧Range
‑
FFT数据计算当前帧回波信号能量，并将其与预先设定阈值进行比较，从而判断当前帧是否包含有目标信息，若小于阈值，则当前帧无目标，若大于等于阈值，则当前帧有目标；步骤三：目标检测：当确定当前帧包含有目标信息后，对当前帧进行目标检测，通过Range
‑
FFT数据中最大峰值点确定人体目标距离雷达位置，并统计连续10帧的最大峰值点，选取出现次数最多的最大峰值点作为人体目标的索引；步骤四：确定人体目标的索引后，根据当前帧当前索引的点计算当前帧的相位信息，并将其循环存储到Phase_buffer中，Phase_buffer为1*1000的实数矩阵；而后通过0.2～0.5Hz，2阶带通滤波器对Phase_buffer进行滤波，取得滤波后的实数相位矩阵breath_phase_buffer，通过0.8～1.6Hz，2阶带通滤波器对Phase_buffer进行滤波，取得滤波后的实数相位矩阵heart_phase_buffer；步骤五：对上述取得的两个实数相位矩阵进行频率估计，分别取得相位变化频率，进而转换为呼吸及心跳频率；步骤六：对每一帧呼吸及心跳数据进行循环存储，分别统计其中出现次数最多的作为当前帧呼吸及心跳结果。2.根据权利要求1所述的一种基于毫米波雷达的呼吸心跳监测系统，其特征在于，所述FMCW毫米波雷达的数据采样频率为20Hz，每帧间隔为50ms，雷达接收到目标返回中频信号后，进行FFT得到Range
‑
FFT维数据，假定该数据长度为N，则该数据为包含目标与雷达的距离信息和目标呼吸心跳引起的相位变化信息的N个复数。3.根据权利要求2所述的一种基于毫米波雷达的呼吸心跳监测系统，其特征在于，所述步骤二中，判断当前检测区域内是否有人体目标：由Range
‑
FFT数据计算N个复数的模值之和作为感应强度，将其与第一阈值进行比较，若其大于第一阈值，则判定当前检测区域内存在人体，若否，则输出“当前无明显生命体征”信息。4.根据权利要求3所述的一种基于毫米波雷达的呼吸心跳监测系统，其特征在于，所述FMCW毫米波雷达监测时需要对数据进行处理，具体包括有慢多普勒维数据累积、静态杂波滤除、目标检测、目标筛选、目标点相位提取、带通滤波及保存呼吸心跳波形和频谱估计法，所述慢多普勒维数据累积：将该帧所有chirp的N点Range
‑
FFT数据对应相加，连续保存M帧数据获得M*N大小的二维矩阵，定义为Range
‑
cube；当Range
‑
cube中的数据存满以后，将最新帧的Range
‑
FFT数据存入Range
‑
cube最后位置，将Range
‑
cube第一帧数据移除，其它数据依次前移。5.根据权利要求4所述的一种基于毫米波雷达的呼吸心跳监测系统，其特征在于，所述静态杂波滤除：Range
‑
cube存满之后，首先对其进行静...

【专利技术属性】
技术研发人员：杨绍分，袁文忠，
申请(专利权)人：德心智能科技常州有限公司，
类型：发明
国别省市：