一种基于脉冲超宽带雷达的非接触血压测量方法技术

本发明专利技术提出一种基于脉冲超宽带雷达的非接触式血压测量方法，属于非接触生命体征检测领域；具体为：首先，超宽带雷达持续发射脉冲信号，经有效探测区域内待测人体、墙体以及地面反射后，接收信号逐行累积，形成二维雷达信号矩阵M1；经去除直流分量，带通滤波和杂波后，挑选能量最大的一列信号，作为生命体征信号S1；然后，通过消除生命体征信号S1的运动和呼吸干扰，获得心跳信号S2和脉搏波信号S3；根据两信号的周期及相关性，提取单拍脉搏波信号S4，标准化后训练血压预测神经网络；最后针对新的受试者，获取其单拍脉搏波信号S4，标准化后输入训练好的血压预测神经网络，进行血压预测。本发明专利技术无需接触皮肤，能够方便且准确地进行血压测量。测量。测量。

【技术实现步骤摘要】
一种基于脉冲超宽带雷达的非接触血压测量方法


[0001]本专利技术属于非接触生命体征检测领域，具体涉及一种基于脉冲超宽带雷达的非接触式血压测量方法。

技术介绍

[0002]高血压是最常见的慢性病之一，据统计，中国约有四分之一的人都患有高血压，它会诱发很多心脑血管疾病。所以，血压监测对于疾病的预防具有重大意义。目前常用的血压计大多是袖带式，需要紧绑在被测者的胳膊上，这对于胳膊受伤不适宜进行皮肤接触的患者来说并不适用，袖带也给测量增加了不便。
[0003]非接触式测量为长期血压监测提供了便利，超宽带雷达具有功耗低、抗干扰能力强以及穿透力强等优势，已经存在许多使用超宽带雷达进行非接触生命体征监测的研究，随着技术的发展，基于雷达的非接触式血压测量方法逐渐出现，具有良好的发展前景。
[0004]目前基于雷达的血压测量方法有使用心电图或光电容积图与雷达结合的方式，这类方法同样需要传感器附着身体，由此带来不便。仅使用雷达的方法一般是：通过获取心脏或颈动脉和腕部桡动脉的跳动波形，根据波形的某些特征，例如收缩周期、幅值、斜率等参数，计算脉搏波传输时间，根据线性回归拟合血压的表达式。但这些方法通常需要被试者在测量时平躺，保持完全静止甚至屏住呼吸，然而当伴随身体微动时，噪声干扰较大，脉搏波会存在失真，难以直接提取出血压相关的特征，拟合度较差。
[0005]现有技术中，如文献1：CN202111244372.5提供了一种基于毫米波雷达信号的人体血压检测方法，在距离人体手腕处预设高度的位置采集毫米波雷达信号，通过三部分神经网络生成预测血压。该方法需要对三部分神经网络分别进行训练，操作较为繁琐，不利于此后长期实时化的监测，且神经网络的泛化能力较弱，全部依赖于神经网络进行数据处理与预测，在真实场景中的鲁棒性也可能会有所不足。文献2：CN202210010289.X也提供了一种基于毫米波雷达的血压检测方法，利用小波包分解重构脉搏波信号，并对其进行特征点检测，根据特征点数据进行回归分析，获得血压检测结果，但在实际场景中往往波形存在畸变，特征点数据会模糊甚至异常，无法良好地拟合与血压的关系。

技术实现思路

[0006]有鉴于此，本专利技术提出一种基于脉冲超宽带雷达的非接触血压测量方法，利用神经网络，能够自动捕捉脉搏波信号的深层时序特征，准确且便捷地进行血压监测。
[0007]所述基于脉冲超宽带雷达的非接触血压测量方法，具体步骤如下：
[0008]步骤一、超宽带雷达持续发射脉冲信号，经有效探测区域内人体、墙体以及地面反射后，被接收天线接收，逐行堆叠累积，形成二维雷达信号矩阵M1；
[0009]二维雷达回波信号矩阵记为：
[0010]其中行向量表示快时间维度，与探测距离呈正相关；列向量表示慢时间维度，与数据累计时间呈正相关，矩阵各元素x
ij
代表雷达信号采样值。
[0011]步骤二、对雷达信号矩阵M1进行预处理操作，去除直流分量，进行带通滤波和去杂波，得到预处理后的雷达信号矩阵M。
[0012]步骤三、从雷达信号矩阵M中按列计算各列的能量值，选择最大的一列信号，视为生命体征信号S1；
[0013]步骤四、通过变分模态分解算法对生命体征信号S1消除运动和呼吸干扰，获得心跳信号S2和脉搏波信号S3。
[0014]首先，对生命体征信号S1，通过变分模态分解获得具有不同中心频率的模态信号，进行快速傅里叶变换计算其频率，选择标准心率值；
[0015]变分模态分解的公式为：
[0016][0017]其中{μ
i
}和{ω
i
}分别对应分解后的第i个模态分量和中心频率，δ(t)为狄拉克函数，β为二次惩罚因子，γ为拉格朗日乘数算子，f(t)为原始信号。
[0018]然后，选取中心频率为标准心率值的模态信号作为心跳信号S2。
[0019]最后，除了能量最小的模态信号外，选取心跳信号及能量小于心跳信号的模态信号叠加，视为脉搏波信号S3。
[0020]步骤五、根据心跳信号S2和脉搏波信号S3的周期及相关性，提取单拍脉搏波信号S4；
[0021]单拍脉搏波提取包括：
[0022]首先，以心跳信号极小值点为分割点，切分单拍心跳信号。
[0023]然后，以最接近心跳信号极小值点的脉搏波信号极小值点为分割点，切分单拍脉搏波信号。
[0024]最后，计算单拍心跳信号与单拍脉搏波信号之间的相关系数，选取相关系数大于0.6的单拍脉搏波信号S4；
[0025]相关系数计算公式为：
[0026]其中g
p
为单拍脉搏波信号，g
h
为单拍心跳信号，Cov()为协方差公式，σ2()为方差公式。
[0027]步骤六、设置血压预测神经网络，将单拍脉搏波信号S4标准化后输入血压预测神经网络进行训练。
[0028]步骤七、针对新的受试者，重复上述步骤一到五，得到单拍脉搏波信号S4，标准化后输入进训练好的血压预测神经网络，直接得到血压结果。
[0029]本专利技术的优点在于：
[0030]本专利技术一种基于脉冲超宽带雷达的非接触式血压测量方法，无需接触皮肤，能够方便且准确地进行血压测量。与基于信号处理的方法不同，本专利技术不仅使用超宽带雷达提取出与血压相关的脉搏波信号，还通过神经网络自适应地分析信号的时序特征，根据特征映射为血压结果进行输出。
附图说明
[0031]图1为本专利技术一种基于脉冲超宽带雷达的非接触式血压测量方法的流程图。
[0032]图2为本专利技术经预处理后的雷达信号矩阵示意图。
[0033]图3为本专利技术提取的生命体征信号示意图。
[0034]图4为本专利技术提取的脉搏波信号示意图。
[0035]图5为本专利技术提取的单拍脉搏波示意图。
[0036]图6为本专利技术血压预测神经网络结构示意图。
具体实施方式
[0037]下面结合实施例和附图对本专利技术做进一步详细解释说明。下面描述的实施例是示例性的，仅用于解释本专利技术，而不能理解为对本专利技术的限制。
[0038]本专利技术一种基于脉冲超宽带雷达的非接触式血压测量方法，利用神经网络，能够自动捕捉脉搏波信号的深层时序特征，准确且便捷地进行血压监测；具体分为以下步骤：
[0039]步骤1：雷达信号获取。超宽带雷达持续发射脉冲信号，信号经有效探测区域内物体、墙体以及地面反射后，被接收天线接收，逐行堆叠累积，形成二维雷达信号矩阵M1；
[0040]步骤2：雷达信号预处理。去除直流分量，进行带通滤波和去杂波。
[0041]步骤3：生命体征信号提取。对预处理后的雷达信号矩阵M，在快时间维度挑选能量最大的一列信号，视为生命体征信号S1，每列信号的能量可以通过计算其方差的形式来表示。
[0042]步骤4：运动和呼吸干扰消除。通过变分模态分解(Variational Mode Decomposition,VMD)，对生命体征信号S本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于脉冲超宽带雷达的非接触式血压测量方法，其特征在于，具体步骤如下：首先，超宽带雷达持续发射脉冲信号，经有效探测区域内人体、墙体以及地面反射后，被接收天线接收，逐行堆叠累积，形成二维雷达信号矩阵M1；并去除直流分量，进行带通滤波和去杂波，得到预处理后的雷达信号矩阵M；然后、从雷达信号矩阵M中按列计算各列的能量值，选择最大的一列信号，视为生命体征信号S1；通过变分模态分解算法对生命体征信号S1消除运动和呼吸干扰，获得心跳信号S2和脉搏波信号S3；并进一步根据心跳信号S2和脉搏波信号S3的周期及相关性，提取单拍脉搏波信号S4；单拍脉搏波提取过程包括：首先，以心跳信号极小值点为分割点，切分单拍心跳信号；然后，以最接近心跳信号极小值点的脉搏波信号极小值点为分割点，切分单拍脉搏波信号；最后，计算单拍心跳信号与单拍脉搏波信号之间的相关系数，选取相关系数大于0.6的单拍脉搏波信号S4；相关系数计算公式为：其中g
p
为单拍脉搏波信号，g
h
为单拍心跳信号，Cov()为协方差公式，σ2()为方差公式；最后、设置血压预测神经网络，将单拍脉搏波信号S4标准化后输入血压预测神经网络进行训练；直至训练完成后，输入新的受试者提取的单拍脉搏波信号S4，标准化后输入，直接得到血压结果。2.根据权利要求1所述的一种基于脉冲超宽带雷达的非接触式血压测量方法，其特征在于，所述二维雷达回波信号矩阵其中行向量表示快时间维度，与探测距离呈正相关；列向量表示慢时间维度，与数据累计时间呈正相关，矩阵各元素x
ij
代表雷达信号采样值。3.根据权利要求1所述的一种基于脉冲超宽带雷达的非接触式血压测量方法，其特征在于，所述对二维雷达信号矩阵M1预处理的过程包括：首先，对雷达信号矩阵M1的每行根据特定应用场景进行分段，每段取平均值，各段减去其对应的平均值，去除直流分量；然后，对去除直流分量的雷达信号矩阵，使用带通滤波器对每一行数据进行滤波；最后，对滤波后的雷达信号矩阵，通过滑动平均算法去除静态杂波，得到雷达信号矩阵M；每行的静态杂波为前一行静态杂波与本行数据的加权和，表示为：C(t,τ)＝a
·
C(t
‑
1,τ)+(1
‑
a)
·
x(t,τ...

【专利技术属性】
技术研发人员：张金辉，张琳，张静，姜夕康，穆文遥，孟令一，李蕾，
申请(专利权)人：中国人民解放军总医院，
类型：发明
国别省市：