CN115607126A - A non-contact blood pressure measurement method based on pulsed ultra-wideband radar - Google Patents
A non-contact blood pressure measurement method based on pulsed ultra-wideband radar Download PDFInfo
- Publication number
- CN115607126A CN115607126A CN202211194354.5A CN202211194354A CN115607126A CN 115607126 A CN115607126 A CN 115607126A CN 202211194354 A CN202211194354 A CN 202211194354A CN 115607126 A CN115607126 A CN 115607126A
- Authority
- CN
- China
- Prior art keywords
- signal
- blood pressure
- pulse wave
- radar
- layer
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
- A61B5/021—Measuring pressure in heart or blood vessels
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- Theoretical Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Medical Informatics (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Physiology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- General Engineering & Computer Science (AREA)
- Cardiology (AREA)
- Software Systems (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Fuzzy Systems (AREA)
- Algebra (AREA)
- Vascular Medicine (AREA)
- Computational Linguistics (AREA)
- Computing Systems (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
Description
技术领域technical field
本发明属于非接触生命体征检测领域,具体涉及一种基于脉冲超宽带雷达的非接触式血压测量方法。The invention belongs to the field of non-contact vital sign detection, and in particular relates to a non-contact blood pressure measurement method based on pulse ultra-wideband radar.
背景技术Background technique
高血压是最常见的慢性病之一,据统计,中国约有四分之一的人都患有高血压,它会诱发很多心脑血管疾病。所以,血压监测对于疾病的预防具有重大意义。目前常用的血压计大多是袖带式,需要紧绑在被测者的胳膊上,这对于胳膊受伤不适宜进行皮肤接触的患者来说并不适用,袖带也给测量增加了不便。Hypertension is one of the most common chronic diseases. According to statistics, about a quarter of the people in China suffer from hypertension, which can induce many cardiovascular and cerebrovascular diseases. Therefore, blood pressure monitoring is of great significance for disease prevention. Most of the commonly used sphygmomanometers are cuff-type and need to be tightly tied to the arm of the subject. This is not suitable for patients whose arms are injured and are not suitable for skin-to-skin contact. The cuff also adds inconvenience to the measurement.
非接触式测量为长期血压监测提供了便利,超宽带雷达具有功耗低、抗干扰能力强以及穿透力强等优势,已经存在许多使用超宽带雷达进行非接触生命体征监测的研究,随着技术的发展,基于雷达的非接触式血压测量方法逐渐出现,具有良好的发展前景。Non-contact measurement provides convenience for long-term blood pressure monitoring. Ultra-wideband radar has the advantages of low power consumption, strong anti-interference ability, and strong penetration. There have been many researches on non-contact vital sign monitoring using ultra-wideband radar. With the development of technology, the radar-based non-contact blood pressure measurement method has gradually appeared and has a good development prospect.
目前基于雷达的血压测量方法有使用心电图或光电容积图与雷达结合的方式,这类方法同样需要传感器附着身体,由此带来不便。仅使用雷达的方法一般是:通过获取心脏或颈动脉和腕部桡动脉的跳动波形,根据波形的某些特征,例如收缩周期、幅值、斜率等参数,计算脉搏波传输时间,根据线性回归拟合血压的表达式。但这些方法通常需要被试者在测量时平躺,保持完全静止甚至屏住呼吸,然而当伴随身体微动时,噪声干扰较大,脉搏波会存在失真,难以直接提取出血压相关的特征,拟合度较差。Currently, radar-based blood pressure measurement methods combine electrocardiography or photoplethysmography with radar. This type of method also requires the sensor to be attached to the body, which brings inconvenience. The method of only using radar is generally: by obtaining the beating waveform of the heart or carotid artery and wrist radial artery, according to certain characteristics of the waveform, such as systolic period, amplitude, slope and other parameters, calculate the pulse wave transit time, according to linear regression Fit an expression for blood pressure. However, these methods usually require the subject to lie flat, remain completely still or even hold his breath during the measurement. However, when the body moves slightly, the noise interference is large, the pulse wave will be distorted, and it is difficult to directly extract blood pressure-related features. The fit is poor.
现有技术中,如文献1:CN202111244372.5提供了一种基于毫米波雷达信号的人体血压检测方法,在距离人体手腕处预设高度的位置采集毫米波雷达信号,通过三部分神经网络生成预测血压。该方法需要对三部分神经网络分别进行训练,操作较为繁琐,不利于此后长期实时化的监测,且神经网络的泛化能力较弱,全部依赖于神经网络进行数据处理与预测,在真实场景中的鲁棒性也可能会有所不足。文献2:CN202210010289.X也提供了一种基于毫米波雷达的血压检测方法,利用小波包分解重构脉搏波信号,并对其进行特征点检测,根据特征点数据进行回归分析,获得血压检测结果,但在实际场景中往往波形存在畸变,特征点数据会模糊甚至异常,无法良好地拟合与血压的关系。In the prior art, document 1: CN202111244372.5 provides a human blood pressure detection method based on millimeter-wave radar signals, collecting millimeter-wave radar signals at a preset height from the human wrist, and generating predictions through a three-part neural network blood pressure. This method needs to train the three parts of the neural network separately, and the operation is cumbersome, which is not conducive to long-term real-time monitoring in the future, and the generalization ability of the neural network is weak, and all rely on the neural network for data processing and prediction. The robustness may also be insufficient. Document 2: CN202210010289.X also provides a blood pressure detection method based on millimeter-wave radar, using wavelet packet decomposition to reconstruct the pulse wave signal, and performing feature point detection on it, and performing regression analysis based on the feature point data to obtain blood pressure detection results , but in the actual scene, the waveform is often distorted, and the feature point data will be blurred or even abnormal, which cannot fit the relationship with blood pressure well.
发明内容Contents of the invention
有鉴于此,本发明提出一种基于脉冲超宽带雷达的非接触血压测量方法,利用神经网络,能够自动捕捉脉搏波信号的深层时序特征,准确且便捷地进行血压监测。In view of this, the present invention proposes a non-contact blood pressure measurement method based on pulse ultra-wideband radar, which can automatically capture the deep temporal characteristics of the pulse wave signal by using a neural network, and accurately and conveniently monitor blood pressure.
所述基于脉冲超宽带雷达的非接触血压测量方法,具体步骤如下:The non-contact blood pressure measurement method based on pulsed ultra-wideband radar, the specific steps are as follows:
步骤一、超宽带雷达持续发射脉冲信号,经有效探测区域内人体、墙体以及地面反射后,被接收天线接收,逐行堆叠累积,形成二维雷达信号矩阵M1;
二维雷达回波信号矩阵记为: The two-dimensional radar echo signal matrix is recorded as:
其中行向量表示快时间维度,与探测距离呈正相关;列向量表示慢时间维度,与数据累计时间呈正相关,矩阵各元素xij代表雷达信号采样值。The row vector represents the fast time dimension, which is positively correlated with the detection distance; the column vector represents the slow time dimension, which is positively correlated with the data accumulation time, and each element x ij of the matrix represents the radar signal sampling value.
步骤二、对雷达信号矩阵M1进行预处理操作,去除直流分量,进行带通滤波和去杂波,得到预处理后的雷达信号矩阵M。Step 2: Preprocessing the radar signal matrix M1, removing the DC component, performing band-pass filtering and removing clutter, to obtain the preprocessed radar signal matrix M.
步骤三、从雷达信号矩阵M中按列计算各列的能量值,选择最大的一列信号,视为生命体征信号S1;
步骤四、通过变分模态分解算法对生命体征信号S1消除运动和呼吸干扰,获得心跳信号S2和脉搏波信号S3。Step 4: Eliminate movement and respiration interference on the vital sign signal S1 through a variational mode decomposition algorithm to obtain a heartbeat signal S2 and a pulse wave signal S3.
首先,对生命体征信号S1,通过变分模态分解获得具有不同中心频率的模态信号,进行快速傅里叶变换计算其频率,选择标准心率值;First, for the vital sign signal S1, obtain modal signals with different center frequencies through variational mode decomposition, perform fast Fourier transform to calculate its frequency, and select a standard heart rate value;
变分模态分解的公式为:The formula for variational modal decomposition is:
其中{μi}和{ωi}分别对应分解后的第i个模态分量和中心频率,δ(t)为狄拉克函数,β为二次惩罚因子,γ为拉格朗日乘数算子,f(t)为原始信号。Where {μ i } and {ω i } correspond to the i-th modal component and center frequency after decomposition respectively, δ(t) is the Dirac function, β is the quadratic penalty factor, γ is the Lagrangian multiplier sub, f(t) is the original signal.
然后,选取中心频率为标准心率值的模态信号作为心跳信号S2。Then, a modal signal whose center frequency is a standard heart rate value is selected as the heartbeat signal S2.
最后,除了能量最小的模态信号外,选取心跳信号及能量小于心跳信号的模态信号叠加,视为脉搏波信号S3。Finally, except for the modal signal with the smallest energy, the heartbeat signal and the modal signals with energy smaller than the heartbeat signal are selected to be superimposed, and regarded as the pulse wave signal S3.
步骤五、根据心跳信号S2和脉搏波信号S3的周期及相关性,提取单拍脉搏波信号S4;
单拍脉搏波提取包括:Single-beat pulse wave extraction includes:
首先,以心跳信号极小值点为分割点,切分单拍心跳信号。Firstly, the single-beat heartbeat signal is segmented by taking the minimum value point of the heartbeat signal as the segmentation point.
然后,以最接近心跳信号极小值点的脉搏波信号极小值点为分割点,切分单拍脉搏波信号。Then, the single-beat pulse wave signal is segmented by taking the minimum value point of the pulse wave signal closest to the minimum value point of the heartbeat signal as the segmentation point.
最后,计算单拍心跳信号与单拍脉搏波信号之间的相关系数,选取相关系数大于0.6的单拍脉搏波信号S4;Finally, calculate the correlation coefficient between the single-beat heartbeat signal and the single-beat pulse wave signal, and select the single-beat pulse wave signal S4 with a correlation coefficient greater than 0.6;
相关系数计算公式为: The formula for calculating the correlation coefficient is:
其中gp为单拍脉搏波信号,gh为单拍心跳信号,Cov()为协方差公式,σ2()为方差公式。Where g p is a single-beat pulse wave signal, g h is a single-beat heartbeat signal, Cov() is a covariance formula, and σ 2 () is a variance formula.
步骤六、设置血压预测神经网络,将单拍脉搏波信号S4标准化后输入血压预测神经网络进行训练。Step 6, setting up the blood pressure prediction neural network, and inputting the blood pressure prediction neural network into the blood pressure prediction neural network after normalizing the single-shot pulse wave signal S4 for training.
步骤七、针对新的受试者,重复上述步骤一到五,得到单拍脉搏波信号S4,标准化后输入进训练好的血压预测神经网络,直接得到血压结果。
本发明的优点在于:The advantages of the present invention are:
本发明一种基于脉冲超宽带雷达的非接触式血压测量方法,无需接触皮肤,能够方便且准确地进行血压测量。与基于信号处理的方法不同,本发明不仅使用超宽带雷达提取出与血压相关的脉搏波信号,还通过神经网络自适应地分析信号的时序特征,根据特征映射为血压结果进行输出。The invention provides a non-contact blood pressure measurement method based on pulse ultra-wideband radar, which can conveniently and accurately measure blood pressure without touching the skin. Different from methods based on signal processing, the present invention not only uses ultra-wideband radar to extract pulse wave signals related to blood pressure, but also adaptively analyzes the timing characteristics of signals through neural networks, and outputs blood pressure results according to feature mapping.
附图说明Description of drawings
图1为本发明一种基于脉冲超宽带雷达的非接触式血压测量方法的流程图。FIG. 1 is a flowchart of a non-contact blood pressure measurement method based on pulsed ultra-wideband radar according to the present invention.
图2为本发明经预处理后的雷达信号矩阵示意图。Fig. 2 is a schematic diagram of the radar signal matrix after preprocessing according to the present invention.
图3为本发明提取的生命体征信号示意图。Fig. 3 is a schematic diagram of vital sign signals extracted by the present invention.
图4为本发明提取的脉搏波信号示意图。Fig. 4 is a schematic diagram of the pulse wave signal extracted by the present invention.
图5为本发明提取的单拍脉搏波示意图。Fig. 5 is a schematic diagram of a single beat pulse wave extracted by the present invention.
图6为本发明血压预测神经网络结构示意图。Fig. 6 is a schematic diagram of the structure of the blood pressure prediction neural network of the present invention.
具体实施方式detailed description
下面结合实施例和附图对本发明做进一步详细解释说明。下面描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。The present invention will be further explained in detail below in conjunction with the embodiments and the accompanying drawings. The embodiments described below are exemplary only for explaining the present invention and should not be construed as limiting the present invention.
本发明一种基于脉冲超宽带雷达的非接触式血压测量方法,利用神经网络,能够自动捕捉脉搏波信号的深层时序特征,准确且便捷地进行血压监测;具体分为以下步骤:A non-contact blood pressure measurement method based on pulsed ultra-wideband radar in the present invention can automatically capture the deep time sequence characteristics of pulse wave signals by using neural network, and accurately and conveniently monitor blood pressure; it is specifically divided into the following steps:
步骤1:雷达信号获取。超宽带雷达持续发射脉冲信号,信号经有效探测区域内物体、墙体以及地面反射后,被接收天线接收,逐行堆叠累积,形成二维雷达信号矩阵M1;Step 1: Radar signal acquisition. The ultra-wideband radar continuously transmits pulse signals. After the signals are reflected by objects, walls and the ground in the effective detection area, they are received by the receiving antenna and stacked and accumulated row by row to form a two-dimensional radar signal matrix M1;
步骤2:雷达信号预处理。去除直流分量,进行带通滤波和去杂波。Step 2: Radar signal preprocessing. Remove the DC component, perform bandpass filtering and clutter removal.
步骤3:生命体征信号提取。对预处理后的雷达信号矩阵M,在快时间维度挑选能量最大的一列信号,视为生命体征信号S1,每列信号的能量可以通过计算其方差的形式来表示。Step 3: Vital sign signal extraction. For the preprocessed radar signal matrix M, select a column of signals with the largest energy in the fast time dimension as the vital sign signal S1, and the energy of each column of signals can be expressed by calculating its variance.
步骤4:运动和呼吸干扰消除。通过变分模态分解(Variational ModeDecomposition,VMD),对生命体征信号S1消除运动和呼吸干扰,获得心跳信号S2和脉搏波信号S3。Step 4: Movement and Breathing Interference Elimination. Through variational mode decomposition (Variational ModeDecomposition, VMD), the vital sign signal S1 is eliminated from movement and breathing interference, and a heartbeat signal S2 and a pulse wave signal S3 are obtained.
步骤5:单拍脉搏波提取。根据心跳信号S2和脉搏波信号S3的周期及相关性,提取单拍脉搏波信号S4,标准化后作为血压预测神经网络的输入。Step 5: Single beat pulse wave extraction. According to the period and correlation between the heartbeat signal S2 and the pulse wave signal S3, the single-beat pulse wave signal S4 is extracted, and standardized as the input of the neural network for blood pressure prediction.
步骤6:设置血压预测神经网络。Step 6: Set up the blood pressure prediction neural network.
步骤7:划分训练集、验证集和测试集,对神经网络进行训练。Step 7: Divide the training set, verification set and test set, and train the neural network.
步骤8:通过训练后的血压预测神经网络对受试者进行血压预测。Step 8: Predict the subject's blood pressure through the trained blood pressure prediction neural network.
如图1所示,具体步骤如下:As shown in Figure 1, the specific steps are as follows:
步骤一、超宽带雷达持续发射脉冲信号,经有效探测区域内人体、墙体以及地面反射后,被接收天线接收,逐行堆叠累积,形成二维雷达信号矩阵M1;
二维雷达回波信号矩阵记为: The two-dimensional radar echo signal matrix is recorded as:
其中行向量表示快时间维度,与探测距离呈正相关;列向量表示慢时间维度,与数据累计时间呈正相关,矩阵各元素xij代表雷达信号采样值。The row vector represents the fast time dimension, which is positively correlated with the detection distance; the column vector represents the slow time dimension, which is positively correlated with the data accumulation time, and each element x ij of the matrix represents the radar signal sampling value.
步骤二、对雷达信号矩阵M1进行预处理操作,去除直流分量,进行带通滤波和去杂波,得到预处理后的雷达信号矩阵M。Step 2: Preprocessing the radar signal matrix M1, removing the DC component, performing band-pass filtering and removing clutter, to obtain the preprocessed radar signal matrix M.
预处理操作包括:首先,对雷达信号矩阵M1的每行根据特定应用场景进行分段,每段取平均值,各段减去其对应的平均值,去除直流分量。The preprocessing operation includes: firstly, segmenting each row of the radar signal matrix M1 according to a specific application scenario, taking an average value for each segment, subtracting the corresponding average value from each segment, and removing the DC component.
本实施例选择每行以156列数据分段;In this embodiment, each row is selected to be segmented with 156 columns of data;
然后,对去除直流分量的雷达信号矩阵,使用通带为6.5GHz-8GHz的带通滤波器对每一行数据进行滤波。Then, for the radar signal matrix from which the DC component has been removed, a bandpass filter with a passband of 6.5GHz-8GHz is used to filter each row of data.
最后,对带通滤波后的雷达信号矩阵,通过滑动平均算法去除静态杂波,得到雷达信号矩阵M;Finally, for the radar signal matrix after bandpass filtering, the static clutter is removed by the moving average algorithm, and the radar signal matrix M is obtained;
每行的静态杂波为前一行静态杂波与本行数据的加权和,表示为:The static clutter of each row is the weighted sum of the static clutter of the previous row and the data of this row, expressed as:
C(t,τ)=a·C(t-1,τ)+(1-a)·x(t,τ)C(t,τ)=a C(t-1,τ)+(1-a)x(t,τ)
其中C(t,τ)为本行的静态杂波,a为权重,在本实施例中设置为0.9,C(t-1,τ)为前一行静态杂波,x(t,τ)为本行数据,初始化时第一行的静态杂波为0.1·x(1,τ)。Among them, C(t,τ) is the static clutter of this row, a is the weight, which is set to 0.9 in this embodiment, C(t-1,τ) is the static clutter of the previous row, and x(t,τ) is For the data in this row, the static clutter of the first row is 0.1·x(1,τ) during initialization.
步骤三、从雷达信号矩阵M中按列计算各列的能量值,选择最大的一列信号,视为生命体征信号S1;
每列信号的能量通过计算其平方和的形式来表示,第j列信号的能量计算公式为:The energy of each column signal is expressed by calculating the sum of its squares, and the energy calculation formula of the jth column signal is:
其中m为信号长度。Where m is the signal length.
步骤四、通过变分模态分解算法对生命体征信号S1消除运动和呼吸干扰,获得心跳信号S2和脉搏波信号S3。Step 4: Eliminate movement and respiration interference on the vital sign signal S1 through a variational mode decomposition algorithm to obtain a heartbeat signal S2 and a pulse wave signal S3.
首先,对生命体征信号S1,通过变分模态分解获得具有不同中心频率的模态信号,进行快速傅里叶变换计算其频率,选择标准心率值;First, for the vital sign signal S1, obtain modal signals with different center frequencies through variational mode decomposition, perform fast Fourier transform to calculate its frequency, and select a standard heart rate value;
正常人心率范围为0.8-2Hz,即每分钟50-120次,所以选取落在0.8-2Hz的频率为心率值。The normal human heart rate range is 0.8-2Hz, that is, 50-120 beats per minute, so the frequency falling within 0.8-2Hz is selected as the heart rate value.
变分模态分解的公式为:The formula for variational modal decomposition is:
其中{μi}和{ωi}分别对应分解后的第i个模态分量和中心频率,δ(t)为狄拉克函数,β为二次惩罚因子,γ为拉格朗日乘数算子,f(t)为原始信号。Where {μ i } and {ω i } correspond to the i-th modal component and center frequency after decomposition respectively, δ(t) is the Dirac function, β is the quadratic penalty factor, γ is the Lagrangian multiplier sub, f(t) is the original signal.
然后,选取中心频率为标准心率值的模态信号作为心跳信号S2。Then, a modal signal whose center frequency is a standard heart rate value is selected as the heartbeat signal S2.
最后,除了能量最小的模态信号外,选取心跳信号及能量小于心跳信号的模态信号叠加,视为脉搏波信号S3。Finally, except for the modal signal with the smallest energy, the heartbeat signal and the modal signals with energy smaller than the heartbeat signal are selected to be superimposed, and regarded as the pulse wave signal S3.
步骤五、根据心跳信号S2和脉搏波信号S3的周期及相关性,提取单拍脉搏波信号S4;
单拍脉搏波提取包括:Single-beat pulse wave extraction includes:
首先,以心跳信号极小值点为分割点,切分单拍心跳信号。Firstly, the single-beat heartbeat signal is segmented by taking the minimum value point of the heartbeat signal as the segmentation point.
然后,以最接近心跳信号极小值点的脉搏波信号极小值点为分割点,切分单拍脉搏波信号。Then, the single-beat pulse wave signal is segmented by taking the minimum value point of the pulse wave signal closest to the minimum value point of the heartbeat signal as the segmentation point.
最后,计算单拍心跳信号与单拍脉搏波信号之间的相关系数,选取相关系数大于0.6的单拍脉搏波信号S4;Finally, calculate the correlation coefficient between the single-beat heartbeat signal and the single-beat pulse wave signal, and select the single-beat pulse wave signal S4 with a correlation coefficient greater than 0.6;
相关系数计算公式为: The formula for calculating the correlation coefficient is:
其中gp为单拍脉搏波信号,gh为单拍心跳信号,Cov()为协方差公式,σ2()为方差公式。Where g p is a single-beat pulse wave signal, g h is a single-beat heartbeat signal, Cov() is a covariance formula, and σ 2 () is a variance formula.
步骤六、设置血压预测神经网络,将单拍脉搏波信号S4标准化后输入血压预测神经网络进行训练。Step 6, setting up the blood pressure prediction neural network, and inputting the blood pressure prediction neural network into the blood pressure prediction neural network after normalizing the single-shot pulse wave signal S4 for training.
血压预测神经网络包括以下结构:The blood pressure prediction neural network includes the following structures:
首先,将标准化后的单拍脉搏波信号S4作为样本,输入两个分支进行不同尺度的卷积;First, the standardized single-beat pulse wave signal S4 is used as a sample, and input to two branches for convolution of different scales;
第一分支的结构依次为5×1一维卷积层,BN层,最大池化层,5×1一维卷积层,BN层,5×1一维卷积层,BN层,1×1一维卷积层,最大池化层;第二分支的结构依次为3×1一维卷积层,BN层,最大池化层,3×1一维卷积层,BN层,3×1一维卷积层,BN层,1×1一维卷积层,最大池化层。The structure of the first branch is 5×1 one-dimensional convolution layer, BN layer, maximum pooling layer, 5×1 one-dimensional convolution layer, BN layer, 5×1 one-dimensional convolution layer, BN layer, 1× 1 one-dimensional convolution layer, maximum pooling layer; the structure of the second branch is 3×1 one-dimensional convolution layer, BN layer, maximum pooling layer, 3×1 one-dimensional convolution layer, BN layer, 3× 1 one-dimensional convolutional layer, BN layer, 1×1 one-dimensional convolutional layer, maximum pooling layer.
然后,将两个分支卷积得到的特征图拼接,依次经过Dropout层,两层门控循环单元(Gated Recurrent Unit,GRU),多头注意力模块,Dropout层和两层全连接层,输出预测结果;Then, the feature maps obtained by the convolution of the two branches are spliced, and then pass through the Dropout layer, two layers of gated recurrent units (Gated Recurrent Unit, GRU), multi-head attention module, dropout layer and two layers of fully connected layers, and output the prediction results ;
其中Dropout层用于防止过拟合,门控循环单元和多头注意力模块能捕获信号的深层时序特征,全连接层能将特征映射为输出的血压值。Among them, the Dropout layer is used to prevent overfitting, the gated recurrent unit and the multi-head attention module can capture the deep temporal features of the signal, and the fully connected layer can map the features to the output blood pressure value.
步骤七、针对新的受试者,重复上述步骤一到五,得到单拍脉搏波信号S4,标准化后输入进训练好的血压预测神经网络,直接得到血压结果。
实施例:Example:
本实例中将脉冲超宽带雷达安装在距离受试者胸部60厘米左右的支架上,被检测的人员为一名成年男性,年龄为23岁,体重为65kg,身高为180cm。In this example, the pulse ultra-wideband radar is installed on a bracket about 60 cm away from the subject's chest. The person to be detected is an adult male, aged 23, weighing 65 kg, and 180 cm tall.
处理流程的步骤如下:The steps of the processing flow are as follows:
步骤1:使用脉冲超宽带雷达获取雷达信号,脉冲超宽带雷达对测试场地持续地发送脉冲信号,经测试场地内的待测试人体以及地面的反射后,被雷达接收天线接收,接收信号以二维矩阵的数据形式进行储存和进一步分析;Step 1: Use the pulse ultra-wideband radar to obtain radar signals. The pulse ultra-wideband radar continuously sends pulse signals to the test site. After being reflected by the human body to be tested and the ground in the test site, it is received by the radar receiving antenna. matrix data form for storage and further analysis;
其中,该二维矩阵的行向量表示脉冲超宽带雷达快时间维度,与探测距离呈正相关,探测距离在本实例中设置为3米,即本实例中二维矩阵的每一个行向量拥有437列的固定长度;列向量表示慢时间维度,与数据累计时间呈正相关,雷达采样帧率被设置为20帧每秒,即本实例中二维矩阵的列向量长度每秒钟增加20行。每组采集10秒数据,即形成200*437的二维矩阵M1,代表实数域。Among them, the row vector of the two-dimensional matrix represents the fast time dimension of the pulse ultra-wideband radar, which is positively correlated with the detection distance. The detection distance is set to 3 meters in this example, that is, each row vector of the two-dimensional matrix in this example has 437 columns The fixed length of ; the column vector represents the slow time dimension, which is positively correlated with the data accumulation time. The radar sampling frame rate is set to 20 frames per second, that is, the length of the column vector of the two-dimensional matrix in this example increases by 20 rows per second. Each group collects data for 10 seconds to form a 200*437 two-dimensional matrix M1, represents the field of real numbers.
步骤2:对雷达信号矩阵M1进行预处理操作,去除直流分量、进行带通滤波和去杂波;得到预处理后的雷达信号矩阵MStep 2: Preprocess the radar signal matrix M1, remove the DC component, perform band-pass filtering and remove clutter; obtain the preprocessed radar signal matrix M
预处理后的雷达信号矩阵如图2所示;The preprocessed radar signal matrix is shown in Figure 2;
步骤3:生命体征信号提取。在静止环境下,由胸膛振动引起的起伏占据最大能量,所以对预处理后的雷达信号矩阵在快时间维度挑选能量最大的一列信号,视为生命体征信号S1,每列信号的能量可以通过计算其方差的形式来表示,生命体征信号如图3所示。Step 3: Vital sign signal extraction. In a static environment, the ups and downs caused by the chest vibration occupy the maximum energy, so the preprocessed radar signal matrix Select a series of signals with the largest energy in the fast time dimension as the vital sign signal S1, The energy of each column signal can be expressed by calculating its variance, and the vital sign signal is shown in Figure 3.
步骤4:运动和呼吸干扰消除。通过变分模态分解,对生命体征信号S1消除运动和呼吸干扰,获得心跳信号和脉搏波信号脉搏波信号如图4所示;Step 4: Movement and Breathing Interference Elimination. Through variational mode decomposition, the vital sign signal S1 eliminates motion and respiration interference, and obtains the heartbeat signal and pulse wave signal The pulse wave signal is shown in Figure 4;
具体步骤包括:Specific steps include:
步骤4-1:对生命体征信号S1,通过变分模态分解获得具有不同中心频率的模态信号,进行快速傅里叶变换(Fast Fourier Transform,FFT)计算其频率,正常人心率范围为0.8-2Hz,所以选取落在0.8-2Hz的频率为心率值。Step 4-1: For the vital sign signal S1, obtain modal signals with different center frequencies through variational mode decomposition, and perform Fast Fourier Transform (FFT) to calculate its frequency. The heart rate range of normal people is 0.8 -2Hz, so select the frequency that falls within 0.8-2Hz as the heart rate value.
步骤4-2:以心率值为标准,选取中心频率最接近心率的模态信号作为心跳信号S2。Step 4-2: Based on the heart rate value as a standard, select the modal signal whose center frequency is closest to the heart rate as the heartbeat signal S2.
步骤4-3:因为环境噪声通常占据最小的能量,为了恢复出具有较多血压相关信息的波形,所以除了能量最小的模态信号外,选取心跳信号及能量小于心跳信号的模态信号叠加,视为脉搏波信号S3。Step 4-3: Because environmental noise usually occupies the smallest energy, in order to restore the waveform with more blood pressure related information, in addition to the modal signal with the smallest energy, select the heartbeat signal and the modal signal with energy less than the heartbeat signal to superimpose, It is regarded as the pulse wave signal S3.
步骤5:单拍脉搏波提取。根据心跳信号S2和脉搏波信号S3的周期及相关性,提取单拍脉搏波信号标准化后作为神经网络的输入,如图5所示,具体步骤包括:Step 5: Single beat pulse wave extraction. According to the cycle and correlation of the heartbeat signal S2 and the pulse wave signal S3, extract the single beat pulse wave signal After normalization, it is used as the input of the neural network, as shown in Figure 5, and the specific steps include:
步骤5-1:以心跳信号极小值点为分割点切分单拍心跳信号。Step 5-1: Segment the single-beat heartbeat signal with the minimum point of the heartbeat signal as the segmentation point.
步骤5-2:以最接近心跳信号极小值点的脉搏波信号极小值点为分割点切分单拍脉搏波信号。Step 5-2: Segment the single-beat pulse wave signal with the minimum value point of the pulse wave signal closest to the minimum value point of the heartbeat signal as the segmentation point.
步骤5-3:计算单拍心跳信号与单拍脉搏波信号之间的相关系数,选取相关系数大于0.6的单拍脉搏波信号,因为脉搏波信号的周期一般不大于1.5秒,所以对单拍脉搏波以1.5秒即30个采样点进行截取或补零操作,标准化后作为神经网络的输入。Step 5-3: Calculate the correlation coefficient between the single-beat heartbeat signal and the single-beat pulse wave signal, and select the single-beat pulse wave signal with a correlation coefficient greater than 0.6, because the period of the pulse wave signal is generally not greater than 1.5 seconds, so the single-beat The pulse wave is intercepted or zero-filled at 30 sampling points in 1.5 seconds, and then standardized as the input of the neural network.
步骤6:设置血压预测神经网络,网络结构如图6所示,包括:Step 6: Set up the neural network for blood pressure prediction. The network structure is shown in Figure 6, including:
首先通过两个分支进行不同尺度的卷积,结构依次为一维卷积层,BN层,最大池化层,一维卷积层,BN层,一维卷积层,BN层,1×1一维卷积层,最大池化层,其中两个分支的前三个卷积层的卷积核大小不同,分别为5和3,其余层相同。Firstly, convolution of different scales is performed through two branches, and the structure is sequentially one-dimensional convolutional layer, BN layer, maximum pooling layer, one-dimensional convolutional layer, BN layer, one-dimensional convolutional layer, BN layer, 1×1 One-dimensional convolution layer, maximum pooling layer, the convolution kernel sizes of the first three convolution layers of the two branches are different, 5 and 3 respectively, and the rest of the layers are the same.
不同尺度的卷积能提取波形的不同特征,1×1卷积核用于升维,BN层用于防止梯度爆炸,最大池化层用于减少计算量,增加鲁棒性。Convolutions of different scales can extract different features of the waveform. The 1×1 convolution kernel is used to increase the dimension, the BN layer is used to prevent gradient explosion, and the maximum pooling layer is used to reduce the amount of calculation and increase robustness.
接着,将两个分支得到的特征图拼接,经过Dropout层,两层门控循环单元,多头注意力模块,Dropout层,两层全连接层,输出预测结果。Dropout层用于防止过拟合,GRU网络能获取序列的时序特征,多头注意力模块能捕捉序列的深层时序关系。Then, the feature maps obtained by the two branches are spliced, and the prediction results are output through the dropout layer, two layers of gated recurrent units, multi-head attention module, dropout layer, and two layers of fully connected layers. The dropout layer is used to prevent over-fitting, the GRU network can obtain the sequence characteristics of the sequence, and the multi-head attention module can capture the deep sequence relationship of the sequence.
步骤7:基于单拍脉搏波信号和其对应标签,对神经网络进行训练,其中标签为血压计测得的真实血压值。损失函数采用均方根误差(Mean Square Error,MSE),激活函数为Relu,优化器为Adam,批量大小为32,共进行50轮训练。Step 7: Train the neural network based on the single-beat pulse wave signal and its corresponding label, where the label is the real blood pressure value measured by the sphygmomanometer. The loss function uses Mean Square Error (Mean Square Error, MSE), the activation function is Relu, the optimizer is Adam, the batch size is 32, and a total of 50 rounds of training are performed.
步骤8:通过训练后的血压预测神经网络对被试者进行血压预测。Step 8: Predict the blood pressure of the subjects through the trained blood pressure prediction neural network.
显然,本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的精神和范围。这样,倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围之内,则本申请也意图包含这些改动和变型在内。Obviously, those skilled in the art can make various changes and modifications to the application without departing from the spirit and scope of the application. In this way, if these modifications and variations of the present application fall within the scope of the claims of the present application and their equivalent technologies, the present application is also intended to include these modifications and variations.
Claims (6)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211194354.5A CN115607126B (en) | 2022-09-28 | 2022-09-28 | Non-contact blood pressure measurement method based on pulse ultra-wideband radar |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211194354.5A CN115607126B (en) | 2022-09-28 | 2022-09-28 | Non-contact blood pressure measurement method based on pulse ultra-wideband radar |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115607126A true CN115607126A (en) | 2023-01-17 |
| CN115607126B CN115607126B (en) | 2023-06-02 |
Family
ID=84861092
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211194354.5A Active CN115607126B (en) | 2022-09-28 | 2022-09-28 | Non-contact blood pressure measurement method based on pulse ultra-wideband radar |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115607126B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116616726A (en) * | 2023-04-21 | 2023-08-22 | 中国人民解放军空军军医大学 | Method and system for suppressing respiratory harmonic wave in UWB impact biological radar heartbeat signal |
| CN117582203A (en) * | 2023-10-19 | 2024-02-23 | 南京理工大学 | Non-contact continuous blood pressure monitoring devices, computer equipment and storage media for nighttime sleep |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150112220A1 (en) * | 2013-10-23 | 2015-04-23 | King Abdullah University Of Science And Technology | Apparatus and method for wireless monitoring using ultra-wideband frequencies |
| WO2017035579A1 (en) * | 2015-08-28 | 2017-03-09 | Monash University | Cuffless blood pressure monitoring system |
| CN108478203A (en) * | 2018-02-08 | 2018-09-04 | 南京理工大学 | A kind of blood pressure measuring method monitoring radar based on single vital sign |
| US20210169361A1 (en) * | 2019-12-10 | 2021-06-10 | Hill-Rom Services, Inc. | Micro-movement and gesture detection using radar |
| CN113303776A (en) * | 2021-05-28 | 2021-08-27 | 西安理工大学 | Non-contact blood pressure measuring method based on cyclic neural network |
| CN114305364A (en) * | 2022-01-05 | 2022-04-12 | 北京科技大学 | Blood pressure detection method, system and equipment based on millimeter wave radar |
-
2022
- 2022-09-28 CN CN202211194354.5A patent/CN115607126B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150112220A1 (en) * | 2013-10-23 | 2015-04-23 | King Abdullah University Of Science And Technology | Apparatus and method for wireless monitoring using ultra-wideband frequencies |
| WO2017035579A1 (en) * | 2015-08-28 | 2017-03-09 | Monash University | Cuffless blood pressure monitoring system |
| CN108478203A (en) * | 2018-02-08 | 2018-09-04 | 南京理工大学 | A kind of blood pressure measuring method monitoring radar based on single vital sign |
| US20210169361A1 (en) * | 2019-12-10 | 2021-06-10 | Hill-Rom Services, Inc. | Micro-movement and gesture detection using radar |
| CN113303776A (en) * | 2021-05-28 | 2021-08-27 | 西安理工大学 | Non-contact blood pressure measuring method based on cyclic neural network |
| CN114305364A (en) * | 2022-01-05 | 2022-04-12 | 北京科技大学 | Blood pressure detection method, system and equipment based on millimeter wave radar |
Non-Patent Citations (1)
| Title |
|---|
| 田泽懿 等: "基于射频信号的非接触式血压监测系统", 《光学仪器》, pages 81 - 87 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116616726A (en) * | 2023-04-21 | 2023-08-22 | 中国人民解放军空军军医大学 | Method and system for suppressing respiratory harmonic wave in UWB impact biological radar heartbeat signal |
| CN116616726B (en) * | 2023-04-21 | 2025-08-15 | 中国人民解放军空军军医大学 | Method and system for suppressing respiratory harmonic wave in UWB impact biological radar heartbeat signal |
| CN117582203A (en) * | 2023-10-19 | 2024-02-23 | 南京理工大学 | Non-contact continuous blood pressure monitoring devices, computer equipment and storage media for nighttime sleep |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115607126B (en) | 2023-06-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Khosrow-Khavar et al. | Automatic and robust delineation of the fiducial points of the seismocardiogram signal for noninvasive estimation of cardiac time intervals | |
| CN106037694B (en) | A kind of continuous blood pressure measurer based on pulse wave | |
| CN108416367B (en) | Sleep staging method based on multi-sensor data decision-level fusion | |
| CN109907752A (en) | A kind of cardiac diagnosis and monitoring method and system of the interference of removal motion artifacts and ecg characteristics detection | |
| CN110537910A (en) | Non-invasive screening system for coronary heart disease based on joint analysis of ECG and heart sound signals | |
| CN115089139B (en) | Personalized physiological parameter measurement method combined with biometric recognition | |
| CN107137071A (en) | It is a kind of to analyze the method that heart impact signal is used for calculating short-term heart beat value | |
| CN110420019A (en) | A kind of depth recurrence heart rate estimation method of ballistocardiography signal | |
| CN115607126B (en) | Non-contact blood pressure measurement method based on pulse ultra-wideband radar | |
| CN109833034A (en) | The method and device of blood pressure data is extracted in a kind of pulse wave signal | |
| CN111528814A (en) | Method for monitoring blood pressure through machine learning based on LSTM neural network | |
| CN111839488A (en) | Non-invasive continuous blood pressure measurement device and method based on pulse wave | |
| CN113116321A (en) | Non-invasive continuous blood pressure measuring system based on PSO-GRNN neural network | |
| Saxena et al. | QRS detection using new wavelets | |
| CN114098679B (en) | Vital sign monitoring waveform recovery method based on deep learning and radio frequency sensing | |
| KR20210066332A (en) | Method and apparatus for determining biometric information of target | |
| CN116451110A (en) | Method for constructing blood glucose prediction model based on signal energy characteristics and pulse cycle | |
| CN114469041A (en) | A method for analyzing the characteristics of heart rate change data during exercise | |
| Quanyu | Pulse signal analysis based on deep learning network | |
| CN114145725B (en) | A PPG sampling rate estimation method based on non-invasive continuous blood pressure measurement | |
| CN113057617B (en) | Non-invasive monitoring system for cardiac output | |
| CN114680848A (en) | Millimeter wave vital sign detection method introducing compressed sensing theory | |
| CN114983412A (en) | Non-contact type brain cognitive load objective detection method based on distributed radar | |
| CN117582203B (en) | Non-contact continuous blood pressure monitoring device, computer equipment and storage medium for nighttime sleep | |
| Lu et al. | Accurate heart beat detection with Doppler radar using bidirectional GRU network |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |