CN115363547A - 一种基于超宽带雷达相参积累的人体生命体征检测方法 - Google Patents

Abstract

本发明公开了一种基于超宽带雷达相参积累的人体生命体征检测方法，属于检测生命体征信息领域，方法包括：超宽带雷达向人体发射短脉冲信号，采集回波信号得到快时间、慢时间二维回波数据矩阵；首先进行去除直流分量、静止杂波以及快时间维滤波等预处理；其次根据最大标准差方法锁定人体基准距离门，提取前后若干个距离门系列一维信号，完成相位补偿后进行时域相参积累处理，获得相参积累信号；最后根据变分模态分解算法分别分离出呼吸模态和心跳模态，通过傅里叶变换，完成对人体呼吸频率和心跳频率的检测。本发明可有效抑制环境的各类噪声和干扰，实现呼吸和心跳信号从雷达回波信号中准确的提取与分离。

Description

一种基于超宽带雷达相参积累的人体生命体征检测方法

技术领域

本发明涉及检测生命体征信息技术领域，涉及一种基于超宽带雷达相参积累的人体生命体征检测方法。

背景技术

生命体征作为衡量健康情况的一项重要指标，可用于辨别是否存在生命以及当前生命体的生理机能。生命体征主要指呼吸、心率、体温、脉搏、血压等生命信息，其中，呼吸和心率最能直接反映当前人体的生理状态。目前生命体征检测方式主要分为两大类，接触式检测和非接触式检测。医院常见的心电监护仪和市场上的一些便携式或可穿戴式的体征检测产品，如指甲式脉搏血氧计、穿戴式心电记录仪，智能手环等大部分是采用接触式生命体征检测技术来实现的。尽管接触式生命体征检测能获得较高的准确度，但在一些特定场合，由于传染性病毒传播能力极强，接近患者进行检测是极其危险的；在一些难以使用复杂的有线连接的情况下，如对婴儿、皮肤大面积烧伤受害者进行心电图（ECG）监测，或在建筑物倒塌、地震期间人类被埋的救援任务中，非接触式生命体征检测发挥着至关重要的作用。基于电磁波检测的超宽带雷达技术作为一种新兴的非接触式生命体征检测手段，具有可以穿透非金属障碍物；不易受到天气、温度、光照等环境因素的影响；保护隐私等优势，在非接触式生命体征检测方面具有独特的优势。但是由于人体的呼吸心跳信号强度微弱，心跳信号强度更是只有呼吸信号强度的十分之一左右，存在心跳信号容易受到环境噪声和呼吸信号谐波干扰的问题，利用常规信号处理方式并无法准确有效的检测生命体征信息。

有鉴于此，特提出本发明。

发明内容

本发明的目的是提供了一种基于超宽带雷达相参积累的人体生命体征检测方法，能从待检测人体的超宽带雷达回波中准确地提取并分离呼吸信号和心跳信号，实现非接触式准确检测人体的呼吸和心跳，进而解决现有技术中存在的上述技术问题。

本发明的目的是通过以下技术方案实现的：

一种基于超宽带雷达相参积累的人体生命体征检测方法，包括：

步骤S1，通过超宽带雷达向待检测人体发射超宽带短脉冲信号；

步骤S2，接收经过待检测人体反射后形成的回波数据，对回波数据采样后得到原始二维回波矩阵；

步骤S3，对原始二维回波矩阵通过去除直流分量、去除静止杂波以及抑制噪声进行预处理得到预处理后的二维回波矩阵；

步骤S4，从预处理后的二维回波矩阵中根据最大标准差方法锁定人体基准距离门；

步骤S5，提取人体基准距离门前后若干个距离门的一维信号作为系列一维信号，对系列一维信号进行相位补偿后，采用时域相参积累方法进行处理，获得相参积累信号；

步骤S6，采用变分模态分解算法从所述相参积累信号分离得出呼吸模态和心跳模态，对分离得到的呼吸模态和心跳模态进行频谱变换得出频谱信息，通过所述频谱信息得到待检测人体的呼吸频率和心跳频率。

与现有技术相比，本发明所提供的基于超宽带雷达相参积累的人体生命体征检测方法，其有益效果包括：

通过对超宽带雷达采集的原始二维回波数据进行一系列预处理，首先去除原始二维回波数据的直流分量，并采用自适应背景抵消对去除直流分量的二维回波数据矩阵的慢时间维去除静止杂波，进一步采用IIR带通滤波器对去除静止杂波二维回波数据矩阵快时间维进行噪声抑制，增强微弱的心跳呼吸信号；不拘泥于传统的常规方法只选取单一的距离门信号进行生命体征检测，而是选取人体基准距离门前后若干个距离门得到的系列一维信号，对该系列一维信号进行相位补偿，进一步采用时域相参积累方法进行处理来提高信噪比，对时域相参积累信号采用VMD算法进行处理，完成待检测人体的呼吸频率和心跳频率的检测；该方法解决了采用超宽带雷达检测生命体征信息时，心跳信号容易受到噪声和呼吸信号谐波干扰的问题，能够准确的将人体呼吸和心跳信号从回波信号中提取并分离出来，得到准确的检测结果。

附图说明

为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域的普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他附图。

图1是本发明实施例提供的基于超宽带雷达相参积累的人体生命体征检测方法流程图。

图2是本发明实施例提供的检测方法中超宽带雷达采集的原始回波数据矩阵图。

图3是本发明实施例提供的检测方法中对原始回波数据矩阵进行预处理的结果图。

图4是本发明实施例提供的检测方法中预处理后矩阵不同快时间距离门的标准误差函数图。

图5是传统方法只选取一个距离门的信号波形图。

图6是本发明实施例提供的检测方法中时域相参积累处理得到的相参积累信号波形图。

图7是本发明实施例提供的基于超宽带雷达相参积累的人体生命体征检测方法检测到呼吸信号的时域波形图。

图8是本发明实施例提供的基于超宽带雷达相参积累的人体生命体征检测方法检测到呼吸信号的频域波形图。

图9是本发明实施例提供的基于超宽带雷达相参积累的人体生命体征检测方法检测到心跳信号的时域波形图。

图10是本发明实施例提供的基于超宽带雷达相参积累的人体生命体征检测方法检测到心跳信号的频域波形图。

图11是本发明实施例提供的基于超宽带雷达相参积累的人体生命体征检测方法和传统单一距离门法以及心电仪测得的呼吸率对比图。

图12是本发明实施例提供的基于超宽带雷达相参积累的人体生命体征检测方法和传统单一距离门法以及心电仪测得的心率对比图。

具体实施方式

下面结合本发明的具体内容，对本发明实施例中的技术方案进行清楚、完整地描述；显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例，这并不构成对本发明的限制。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明的保护范围。

首先对本文中可能使用的术语进行如下说明：

术语“和/或”是表示两者任一或两者同时均可实现，例如，X和/或Y表示既包括“X”或“Y”的情况也包括“X和Y”的三种情况。

术语“包括”、“包含”、“含有”、“具有”或其它类似语义的描述，应被解释为非排它性的包括。例如：包括某技术特征要素（如原料、组分、成分、载体、剂型、材料、尺寸、零件、部件、机构、装置、步骤、工序、方法、反应条件、加工条件、参数、算法、信号、数据、产品或制品等），应被解释为不仅包括明确列出的某技术特征要素，还可以包括未明确列出的本领域公知的其它技术特征要素。

术语“由……组成”表示排除任何未明确列出的技术特征要素。若将该术语用于权利要求中，则该术语将使权利要求成为封闭式，使其不包含除明确列出的技术特征要素以外的技术特征要素，但与其相关的常规杂质除外。如果该术语只是出现在权利要求的某子句中，那么其仅限定在该子句中明确列出的要素，其他子句中所记载的要素并不被排除在整体权利要求之外。

除另有明确的规定或限定外，术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解，例如：可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本文中的具体含义。

术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”、“顺时针”、“逆时针”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述和简化描述，而不是明示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本文的限制。

下面对本发明所提供的基于超宽带雷达相参积累的人体生命体征检测方法进行详细描述。本发明实施例中未作详细描述的内容属于本领域专业技术人员公知的现有技术。本发明实施例中未注明具体条件者，按照本领域常规条件或制造商建议的条件进行。本发明实施例中所用试剂或仪器未注明生产厂商者，均为可以通过市售购买获得的常规产品。

如图1所示，本发明实施例提供一种基于超宽带雷达相参积累的人体生命体征检测方法，包括：

步骤S1，通过超宽带雷达向待检测人体发射超宽带短脉冲信号；

步骤S2，接收经过待检测人体反射后形成的回波数据，对回波数据采样后得到原始二维回波矩阵；

步骤S3，对原始二维回波矩阵通过去除直流分量、去除静止杂波以及抑制噪声进行预处理得到预处理后的二维回波矩阵；

步骤S4，从预处理后的二维回波矩阵中根据最大标准差方法锁定人体基准距离门；

步骤S5，提取人体基准距离门前后若干个距离门的一维信号作为系列一维信号，对系列一维信号进行相位补偿后，采用时域相参积累方法进行处理，获得相参积累信号；

步骤S6，采用变分模态分解算法(VMD)从所述相参积累信号分离得出呼吸模态和心跳模态，对分离得到的呼吸模态和心跳模态进行频谱变换得出频谱信息，通过所述频谱信息得到待检测人体的呼吸频率和心跳频率。

上述检测方法的步骤S1中，发射的所述超宽带短脉冲信号为调制高斯脉冲信号，第m个周期的超宽带短脉冲信号

表示为：

其中，A为超宽带雷达发射超宽带脉冲信号的幅值；

为控制脉冲宽度的标准差；τ为快时间；T _s 为超宽带雷达发射信号的重复周期；f _c 为超宽带雷达发射信号的中心频率。

上述检测方法的步骤S2中，发射M个超宽带短脉冲信号，采样量化形成的所述原始二维回波矩阵R[m,n]表征为：

上式中，α_l、τ_l分别是周围第l个静止目标冲激响应的幅度和延时；

、

分别是待检测人体胸腔冲激响应的幅度和延时；T _f 表示快时间维的采样周期；n=1，2，…，N，其中N是快时间总采样点数；m=1，2，…，M，其中M是发射超宽带短脉冲信号总数。

上述检测方法的步骤S3中，按以下方式对原始二维回波矩阵去除直流分量、静止杂波以及快时间维滤波进行预处理，包括：

步骤S31，去除直流分量，采用均值相减法去除原始二维回波矩阵R[m,n]的直流分量，得到的去除直流分量的二维回波矩阵O[m,n]为：

其中，c[n]是原始二维回波矩阵R[m,n]每列的均值；

步骤S32，去除静止杂波，采用自适应背景抵消方法对去除直流分量的二维回波矩阵O[m,n]进行去除静止杂波，得到的去除静止杂波的二维回波矩阵V[m,n]为：

其中,静止杂波的

为：

其中，η∈[0,1],一般取值为0.9, Q[m,1]=O[m,1]；

步骤S33，抑制噪声，采用IIR带通滤波器对去除静止杂波的二维回波矩阵V[m,n]快时间维进行噪声抑制，得到预处理后的二维回波矩阵T[m,n]，IIR带通滤波器的差分公式如下：

其中，a0，a1，···a6；b1，b2，···b6为滤波器系数，T[m,0]， T[m,1]，…，T[m, 6]均为0；滤波之后得到预处理后的二维回波矩阵T[m,n]。

上述检测方法的步骤S4中，人体基准距离门是预处理后的二维回波矩阵T[m,n]标准差最大的快时间距离门

，为：

其中，SD _n 为预处理后的二维回波矩阵T[m,n]不同快时间距离门的标准差；

为预处理后的二维回波矩阵T[m,n]不同快时间距离门的均值。

上述检测方法的步骤S5中，按以下方式提取人体基准距离门前后若干个距离门的一维信号作为系列一维信号，对系列一维信号进行相位补偿后，采用时域相参积累方法进行处理，获得相参积累信号，包括：

步骤S51，在预处理后的二维回波矩阵T[m,n]中，根据待检测人体的胸腔厚度，选取人体基准距离门前后各5个快时间距离门，获取系列一维信号S _i (m)表征为：

，

对S _i (m)进行[-1,1]归一化处理得到：

其中，min[ ]表示最小值操作符，max[ ]表示最大值操作符；

步骤S52，按以下方式对系列一维信号进行相位补偿，以第一个信号为基准，对与其反相的信号进行反向处理得到的相位补偿后的信号

为：

上式中, ΔΦ表示为

和

的相位差；

步骤S53，采用时域积累方法获得相参积累信号s(t _m )，该相参积累信号s(t _m )表征为：

其中，t _m 为离散慢时间，t _m =mT _s ；m=1，2，…，M，其中M是发射脉冲信号总数；T _s 是脉冲信号的重复周期。

上述检测方法的步骤S6中，按以下方式采用变分模态分解算法从所述相参积累信号分离得出呼吸模态和心跳模态，对分离得到的呼吸模态和心跳模态进行频谱变换得出频谱信息，通过所述频谱信息得到待检测人体的呼吸频率和心跳频率，包括：

步骤S61，采用变分模态分解算法对相参积累信号s(t _m )进行模态分解，提取呼吸模态与心跳模态，提取模态分量的参数模型表征为：

其中，K为需要分解的模态个数；

分别是分解后第K个模态分量和其对应的中心频率；

分别是所有模态分量和中心频率的集合，呼吸模态分量和心跳模态分量包含于

之中；

是模态分量

的希尔伯特变换，

是冲激函数，

表示梯度运算符，

表示平方范数运算符，表示所有

模态分量之和；

步骤S62，求解模态分量

和中心频率

：引入拉格朗日乘法算子λ，将约束参数模型转变为非约束参数模型，得到增广拉格朗日表达式：

上式中，α为二次惩罚因子，作用是降低高斯噪声的干扰；

为约束函数；利用交替方向乘子算法，求解各模态分量

和对应的中心频率

；

步骤S63，检测呼吸频率和心跳频率：对每个模态分量

进行傅里叶变换，得到

的频谱：

，搜寻每个模态分量

频谱的峰值频率

，如果

∈ [0.2,0.5]Hz，即为当前时间段待检测人体的呼吸频率；如果

∈[0.9,1.7]Hz，即为当前时间段待检测人体的心跳频率，完成对待检测人体呼吸频率和心跳频率的检测。

综上可见，与现有技术相比，本发明实施例的检测方法通过对超宽带雷达采集的原始二维回波数据进行一系列预处理，首先去除原始二维回波数据的直流分量，并采用自适应背景抵消对去除直流分量的二维回波数据矩阵的慢时间维去除静止杂波，进一步采用IIR带通滤波器对去除静止杂波二维回波数据矩阵快时间维进行噪声抑制，增强微弱的心跳呼吸信号；不拘泥于传统的常规方法只选取单一的距离门信号进行生命体征检测，而是选取人体基准距离门前后若干个距离门得到的系列一维信号，对该系列一维信号进行相位补偿，进一步采用时域相参积累方法进行处理来提高信噪比，对时域相参积累信号采用VMD算法进行处理，完成待检测人体的呼吸频率和心跳频率的检测；该方法解决了采用超宽带雷达检测生命体征信息时，心跳信号容易受到噪声和呼吸信号谐波干扰的问题，能够准确的将人体呼吸和心跳信号从回波信号中提取并分离出来，得到准确的检测结果。

为了更加清晰地展现出本发明所提供的技术方案及所产生的技术效果，下面以具体实施例对本发明实施例所提供的基于超宽带雷达相参积累的人体生命体征检测方法进行详细描述。

实施例1

参见图1，本发明实施例一种基于超宽带雷达相参积累的人体生命体征检测方法，包括步骤如下：

步骤S1：产生和发射超宽带信号；

通过超宽带雷达发射的超宽带脉冲信号源产生中心频率为4.3GHz、带宽为2GHz的调制高斯脉冲信号，并由超宽带雷达的发射天线向待检测人体发射，该超宽带脉冲信号

表示为：

（1）

其中，A为超宽带雷达发射超宽带脉冲信号的幅值；

为控制脉冲宽度的标准差；τ为快时间；T _s 为超宽带雷达发射信号的重复周期；f _c 为超宽带雷达发射信号的中心频率；

步骤S2：接收待检测人体的回波数据；

超宽带雷达发射的调制高斯脉冲信号经过待检测人体的反射，在超宽带雷达的接收天线处被接收，通过ADC采样量化回波数据，形成一个大小M×N的原始二维回波矩阵R[m, n]，如图2所示。n=1，2，…，N，其中，N为快时间维，即快时间总采样点数，是图2中的横坐标维；m=1，2，…，M，其中，M为慢时间维，即发射超宽带短脉冲信号总数，是图2中的纵坐标维；快时间为超宽带雷达发射的脉冲信号传播的时间，其单位为ns；慢时间为超宽带雷达对待检测人体的探测时间，其单位为s；

步骤S3：对原始二维回波矩阵R[m,n]进行去除直流分量、去除静止杂波和快时间维滤波的预处理；

步骤S31：采用均值相减法去除原始二维回波矩阵R[m,n]的直流分量，得到去除直流分量的二维回波矩阵O[m,n]，处理方法如下：

计算原始二维回波矩阵R[m,n]慢时间维均值c[n]：

（2）

通过以下均值相减法公式计算得出去除直流分量的二维回波矩阵O[m,n]：

（3）

得到去直流分量的二维回波矩阵O[m,n]；

步骤S32：对去除直流分量的二维回波矩阵O[m,n]，采用自适应背景抵消对去除静止杂波得到去除静止杂波后的二维回波矩阵V[m,n]，公式为：

（4）

其中，η∈[0,1],一般取值为0.9, Q[m,1]=O[m,1]；

超宽带雷达接收回波数据通常包括待检测人体和周围静止目标（如桌椅、墙壁等）的反射回波，静止杂波的幅度远大于待检测人体回波，从而使得待检测人体回波淹没于静止杂波中，因此采用自适应背景抵消对去除静止杂波；

步骤S33：采用IIR带通滤波器对去除静止杂波的二维回波矩阵V[m,n]快时间维进行噪声抑制：

（5）

其中，a0，a1，···a6；b1，b2，···b6为滤波器系数，T[m,0]， T[m,1]，…，T[m, 6]均为0；滤波之后得到预处理后的二维回波矩阵T[m,n]，如图3所示。

步骤S4：采用标准差最大法锁定人体基准距离门

；

（6）

其中，SD _n 为预处理后的二维回波矩阵T[m,n]不同快时间距离门的标准差；

为预处理后的二维回波矩阵不同快时间距离门的均值；

步骤S5：获取相参积累信号；

步骤S51：根据待检测人体胸腔厚度，选取人体基准距离门前后各5个快时间距离门，则系列一维信号S _i (m)可以表示为：

（7）

对选取的所有一维慢时间信号进行[-1,1]归一化处理：

；

步骤S52：以其中第一个信号为基准，对与其反相的信号进行相位补偿，得到相位补偿后的信号

为：

（8）

上式中,ΔΦ表示为

和

的相位差；

在相位补偿处理后再进行时域相参积累得到相参积累信号s(t _m )，如图6所示，图5则是传统单一距离门法得到的慢时间信号；

（9）

其中，t _m 为离散慢时间，t _m = mT _s， m是慢时间采样点数，m = 1，2，…，M；T _s 是脉冲信号的重复周期，T _s =80ms；

步骤S6：变分模态分解（VMD）算法分解相参积累信号s(t _m )，并对提取的呼吸模态、心跳模态分别做FFT变换：

步骤S61：假设相参积累信号s(t _m )被分解为K个分量，保证分解信号为具有中心频率的有限带宽的模态分量，同时各模态的估计带宽之和最小，约束条件为所有模态之和与原始信号相等，则相应提取呼吸、心跳模态的参数优化模型可以表征为：

（10）

其中，K为需要分解的模态个数；

分别是分解后第K个模态分量和其对应的中心频率；

分别是所有模态分量和中心频率的集合，呼吸模态分量和心跳模态分量包含于

之中；

是模态分量

的希尔伯特变换，

是冲激函数，

表示梯度运算符，

表示平方范数运算符，表示所有

模态分量之和；

步骤S62：然后引入拉格朗日乘法算子λ，将约束参数模型转变为非约束参数模型，得到增广拉格朗日表达式为：

（11）

上式中，α为二次惩罚因子，作用是降低高斯噪声的干扰；

为约束函数；利用交替方向乘子算法，求解各模态分量

和对应的中心频率

；模态分量

和中心频率

的参数模型主要求解过程如下：

（12.1）

（12.2）

（12.3）

式中，γ为噪声容忍度，满足信号分解的保真度要求，

和

分别对应

和

的傅里叶变换。

初始化

和最大迭代次数

；

利用公式 (12.1) 和 (12.2) 更新

和

；

利用公式 (12.3) 更新；

精度收敛判据ε＞0，若不满足

且 n＜N，则返回步骤S62，否则完成迭代，输出最终的

和

；

步骤S63：如图7和图9所示对提取的呼吸模态、心跳模态分别做FFT变换；在正常呼吸频段0.2-0.5Hz内的频谱峰值频率即为当前时间段的呼吸频率，在正常心跳频段0.9-1.7Hz内的频谱峰值频率即为当前心跳频率，如图8所示，测得目标的呼吸率为0.333Hz，约20次/分钟（BPM）；如图10所示，测得的心率为1.167Hz，约70次/分钟（BPM）。

下面通过一组实测数据对本发明检测方法的效果作进一步说明。

实验采集距离雷达1.7m的人体回波数据，一共2分钟，雷达采集同时使用接触式心电仪检测呼吸率和心率，并进行数据进行对比。对2分钟的数据进行滑窗处理，窗长为24s，滑动步长为4s，分别采用传统单一距离门法和本发明相参积累方法进行处理，并且都采用VMD算法提取呼吸模态和心跳模态。图11为本发明的基于超宽带雷达相参积累的人体生命体征检测方法和传统单一距离门法以及心电仪测得的呼吸率对比图；图12为本发明的基于超宽带雷达相参积累的生命体征检测方法和传统单一距离门法以及心电仪测得的心率对比图。通过比较可以看出，本发明的相参积累处理之后心跳信号的检测准确率明显提高，而传统单一距离门法检测的结果，其检测的呼吸率准确率较高，但是检测的心率和心电仪参考数据相比有着较大误差。

综上可见，本发明实施例的检测方法，可有效抑制环境的各类噪声和干扰，能在一定区域内、隔着某些介质，如衣物、木板等，且在不接触被测人体目标的情况下从雷达回波信号中准确的提取与分离呼吸信号和心跳信号，实现对人体体征的准确检测，很好的解决了由于人体的呼吸心跳信号强度微弱，心跳信号强度更是只有呼吸信号强度的十分之一左右，现有方法利用超宽带雷达系统检测生命体征时，容易受到周围环境杂波和噪声干扰，加大了在低信噪比的回波信号中提取呼吸、心跳信号的难度。

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明披露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求书的保护范围为准。本文背景技术部分公开的信息仅仅旨在加深对本发明的总体背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。

Claims (7)

1.一种基于超宽带雷达相参积累的人体生命体征检测方法，其特征在于，包括：

步骤S1，通过超宽带雷达向待检测人体发射超宽带短脉冲信号；

步骤S2，接收经过待检测人体反射后形成的回波数据，对回波数据采样后得到原始二维回波矩阵；

步骤S3，对原始二维回波矩阵通过去除直流分量、去除静止杂波以及抑制噪声进行预处理得到预处理后的二维回波矩阵；

步骤S4，从预处理后的二维回波矩阵中根据最大标准差方法锁定人体基准距离门；

步骤S5，提取人体基准距离门前后若干个距离门的一维信号作为系列一维信号，对系列一维信号进行相位补偿后，采用时域相参积累方法进行处理，获得相参积累信号；

步骤S6，采用变分模态分解算法从所述相参积累信号分离得出呼吸模态和心跳模态，对分离得到的呼吸模态和心跳模态进行频谱变换得出频谱信息，通过所述频谱信息得到待检测人体的呼吸频率和心跳频率。

2.根据权利要求1所述的基于超宽带雷达相参积累的人体生命体征检测方法，其特征在于，所述步骤S1中，发射的所述超宽带短脉冲信号为调制高斯脉冲信号，第m个周期的超宽带短脉冲信号

表示为：

其中，A为超宽带雷达发射超宽带脉冲信号的幅值；

为控制脉冲宽度的标准差；τ为快时间；T _s 为超宽带雷达发射信号的重复周期；f _c 为超宽带雷达发射信号的中心频率。

3.根据权利要求1所述的基于超宽带雷达相参积累的人体生命体征检测方法，其特征在于，所述步骤S2中，发射M个超宽带短脉冲信号，采样量化形成的所述原始二维回波矩阵R[m,n]表征为：

上式中，α_l、τ_l分别是周围第l个静止目标冲激响应的幅度和延时；

、

分别是待检测人体胸腔冲激响应的幅度和延时；T _f 表示快时间维的采样周期；n=1，2，…，N，其中N是快时间总采样点数；m=1，2，…，M，其中M是发射超宽带短脉冲信号总数。

4.根据权利要求1-3任一项所述的基于超宽带雷达相参积累的人体生命体征检测方法，其特征在于，所述步骤S3中，按以下方式对原始二维回波矩阵去除直流分量、静止杂波以及快时间维滤波进行预处理，包括：

步骤S31，去除直流分量，采用均值相减法去除原始二维回波矩阵R[m,n]的直流分量，得到的去除直流分量的二维回波矩阵O[m,n]为：

其中，c[n]是原始二维回波矩阵R[m,n]每列的均值:c[n]

步骤S32，去除静止杂波，采用自适应背景抵消方法对去除直流分量的二维回波矩阵O[m,n]进行去除静止杂波，得到的去除静止杂波的二维回波矩阵V[m,n]为：

其中,静止杂波的

为：

其中,η∈[0,1],一般取值为0.9, Q[m,1]=O[m,1]；

步骤S33，抑制噪声，采用IIR带通滤波器对去除静止杂波的二维回波矩阵V[m,n]快时间维进行噪声抑制，得到预处理后的二维回波矩阵T[m,n]，IIR带通滤波器的差分公式如下：

其中，a0，a1，···a6；b1，b2，···b6为滤波器系数，T[m,0]， T[m,1]，…，T[m,6]均为0；滤波之后得到预处理后的二维回波矩阵T[m,n]。

5.根据权利要求4所述的基于超宽带雷达相参积累的人体生命体征检测方法，其特征在于，所述步骤S4中，人体基准距离门是预处理后的二维回波矩阵T[m,n]标准差最大的快时间距离门

，为：

其中，SD _n 为预处理后的二维回波矩阵T[m,n]不同快时间距离门的标准差；

为预处理后的二维回波矩阵T[m,n]不同快时间距离门的均值。

6.根据权利要求5所述的基于超宽带雷达相参积累的人体生命体征检测方法，其特征在于，所述步骤S5中，按以下方式提取人体基准距离门前后若干个距离门的一维信号作为系列一维信号，对系列一维信号进行相位补偿后，采用时域相参积累方法进行处理，获得相参积累信号，包括：

步骤S51，在预处理后的二维回波矩阵T[m,n]中，根据待检测人体的胸腔厚度，选取人体基准距离门前后各5个快时间距离门，获取系列一维信号S _i (m)表征为：

对S _i (m)进行[-1,1]归一化处理得到：

其中，min[ ]表示最小值操作符，max[ ]表示最大值操作符；

步骤S52，按以下方式对系列一维信号进行相位补偿，以第一个信号为基准，对与其反相的信号进行反向处理得出相位补偿后的信号

为：

上式中,ΔΦ表示为

和

的相位差；

在相位补偿处理后再进行时域相参积累得到相参积累信号s(t _m )，该相参积累信号s(t _m )表征为：

其中，t _m 为离散慢时间，t _m =mT _s ；m=1，2，…，M，其中M是发射脉冲信号总数；T _s 是脉冲信号的重复周期。

7.根据权利要求6所述的基于超宽带雷达相参积累的人体生命体征检测方法，其特征在于，所述步骤S6中，按以下方式采用变分模态分解算法从所述相参积累信号分离得出呼吸模态和心跳模态，对分离得到的呼吸模态和心跳模态进行频谱变换得出频谱信息，通过所述频谱信息得到待检测人体的呼吸频率和心跳频率，包括：

步骤S61，采用变分模态分解算法对相参积累信号s(t _m )进行模态分解，提取呼吸模态与心跳模态，提取模态分量的参数模型表征为：

其中，K为需要分解的模态个数；

分别是分解后第K个模态分量和其对应的中心频率；

分别是所有模态分量和中心频率的集合，呼吸模态分量和心跳模态分量包含于

之中；

是模态分量

的希尔伯特变换，

是冲激函数，

表示梯度运算符，

表示平方范数运算符，表示所有

模态分量之和；

步骤S62，求解模态分量

和中心频率

：引入拉格朗日乘法算子λ，将约束参数模型转变为非约束参数模型，得到增广拉格朗日表达式：

上式中，α为二次惩罚因子，作用是降低高斯噪声的干扰；

为约束函数；利用交替方向乘子算法，求解各模态分量

和对应的中心频率

；

步骤S63，检测呼吸频率和心跳频率：对每个模态分量

进行傅里叶变换，得到

的频谱：

，搜寻每个模态分量

频谱的峰值频率

，如果

∈ [0.2,0.5]Hz，即为当前时间段待检测人体的呼吸频率；如果

∈[0.9,1.7]Hz，即为当前时间段待检测人体的心跳频率，完成对待检测人体呼吸频率和心跳频率的检测。

Images