[go: up one dir, main page]

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 PDF

Info

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
Application number
CN202211194354.5A
Other languages
Chinese (zh)
Other versions
CN115607126B (en
Inventor
张金辉
张琳
张静
姜夕康
穆文遥
孟令一
李蕾
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chinese PLA General Hospital
Original Assignee
Chinese PLA General Hospital
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chinese PLA General Hospital filed Critical Chinese PLA General Hospital
Priority to CN202211194354.5A priority Critical patent/CN115607126B/en
Publication of CN115607126A publication Critical patent/CN115607126A/en
Application granted granted Critical
Publication of CN115607126B publication Critical patent/CN115607126B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/021Measuring pressure in heart or blood vessels
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7253Details of waveform analysis characterised by using transforms
    • A61B5/7257Details of waveform analysis characterised by using transforms using Fourier transforms
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information 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

The invention provides a non-contact blood pressure measuring method based on a pulse ultra-wideband radar, belonging to the field of non-contact vital sign detection; the method specifically comprises the following steps: firstly, an ultra-wideband radar continuously transmits pulse signals, and after the pulse signals are reflected by a human body to be detected, a wall body and the ground in an effective detection area, received signals are accumulated line by line to form a two-dimensional radar signal matrix M1; after direct current components are removed, band-pass filtering and noise waves are carried out, a row of signals with the largest energy are selected and used as vital sign signals S1; then, by eliminating the movement and the respiratory interference of the vital sign signal S1, a heartbeat signal S2 and a pulse wave signal S3 are obtained; extracting a single-beat pulse wave signal S4 according to the period and the correlation of the two signals, and training a blood pressure prediction neural network after standardization; and finally, aiming at a new subject, acquiring a single-beat pulse wave signal S4 of the new subject, standardizing the signal, and inputting the standardized signal into a trained blood pressure prediction neural network to predict the blood pressure. The invention does not need to contact the skin and can conveniently and accurately measure the blood pressure.

Description

一种基于脉冲超宽带雷达的非接触血压测量方法A non-contact blood pressure measurement method based on pulsed ultra-wideband radar

技术领域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;Step 1. The ultra-wideband radar continuously transmits pulse signals, which are received by the receiving antenna after being reflected by the human body, wall and ground in the effective detection area, and stacked and accumulated row by row to form a two-dimensional radar signal matrix M1;

二维雷达回波信号矩阵记为:

Figure BDA0003870343730000021
The two-dimensional radar echo signal matrix is recorded as:
Figure BDA0003870343730000021

其中行向量表示快时间维度,与探测距离呈正相关;列向量表示慢时间维度,与数据累计时间呈正相关,矩阵各元素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;Step 3. Calculate the energy value of each column from the radar signal matrix M by column, and select the largest column signal as the vital sign signal 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:

Figure BDA0003870343730000022
Figure BDA0003870343730000022

其中{μ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;Step 5, according to the cycle and correlation of the heartbeat signal S2 and the pulse wave signal S3, extract the single beat pulse wave signal 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;

相关系数计算公式为:

Figure BDA0003870343730000031
The formula for calculating the correlation coefficient is:
Figure BDA0003870343730000031

其中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,标准化后输入进训练好的血压预测神经网络,直接得到血压结果。Step 7. For new subjects, repeat the above steps 1 to 5 to obtain the single-shot pulse wave signal S4, and input it into the trained blood pressure prediction neural network after normalization to directly obtain the blood pressure result.

本发明的优点在于: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;Step 1. The ultra-wideband radar continuously transmits pulse signals, which are received by the receiving antenna after being reflected by the human body, wall and ground in the effective detection area, and stacked and accumulated row by row to form a two-dimensional radar signal matrix M1;

二维雷达回波信号矩阵记为:

Figure BDA0003870343730000041
The two-dimensional radar echo signal matrix is recorded as:
Figure BDA0003870343730000041

其中行向量表示快时间维度,与探测距离呈正相关;列向量表示慢时间维度,与数据累计时间呈正相关,矩阵各元素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;Step 3. Calculate the energy value of each column from the radar signal matrix M by column, and select the largest column signal as the vital sign signal 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:

Figure BDA0003870343730000042
Figure BDA0003870343730000042

其中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:

Figure BDA0003870343730000051
Figure BDA0003870343730000051

其中{μ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;Step 5, according to the cycle and correlation of the heartbeat signal S2 and the pulse wave signal S3, extract the single beat pulse wave signal 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;

相关系数计算公式为:

Figure BDA0003870343730000052
The formula for calculating the correlation coefficient is:
Figure BDA0003870343730000052

其中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,标准化后输入进训练好的血压预测神经网络,直接得到血压结果。Step 7. For new subjects, repeat the above steps 1 to 5 to obtain the single-shot pulse wave signal S4, and input it into the trained blood pressure prediction neural network after normalization to directly obtain the blood pressure result.

实施例: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,

Figure BDA0003870343730000061
代表实数域。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,
Figure BDA0003870343730000061
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:生命体征信号提取。在静止环境下,由胸膛振动引起的起伏占据最大能量,所以对预处理后的雷达信号矩阵

Figure BDA0003870343730000062
在快时间维度挑选能量最大的一列信号,视为生命体征信号S1,
Figure BDA0003870343730000063
每列信号的能量可以通过计算其方差的形式来表示,生命体征信号如图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
Figure BDA0003870343730000062
Select a series of signals with the largest energy in the fast time dimension as the vital sign signal S1,
Figure BDA0003870343730000063
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消除运动和呼吸干扰,获得心跳信号

Figure BDA0003870343730000064
和脉搏波信号
Figure BDA0003870343730000065
脉搏波信号如图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
Figure BDA0003870343730000064
and pulse wave signal
Figure BDA0003870343730000065
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的周期及相关性,提取单拍脉搏波信号

Figure BDA0003870343730000071
标准化后作为神经网络的输入,如图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
Figure BDA0003870343730000071
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)

1. A non-contact blood pressure measuring method based on pulse ultra-wideband radar is characterized by comprising the following specific steps:
firstly, an ultra-wideband radar continuously transmits pulse signals, the pulse signals are reflected by a human body, a wall body and the ground in an effective detection area, and then are received by a receiving antenna, and are stacked and accumulated line by line to form a two-dimensional radar signal matrix M1; removing direct current components, performing band-pass filtering and clutter removal to obtain a preprocessed radar signal matrix M;
then, calculating the energy value of each column from the radar signal matrix M according to the column, selecting the largest column of signals as the vital sign signals S1; eliminating motion and respiration interference on the vital sign signal S1 through a variational modal decomposition algorithm to obtain a heartbeat signal S2 and a pulse wave signal S3; further extracting a single-beat pulse wave signal S4 according to the period and the correlation of the heartbeat signal S2 and the pulse wave signal S3;
the single-beat pulse wave extraction process comprises the following steps:
firstly, taking a heartbeat signal minimum value point as a segmentation point, and segmenting a single-beat heartbeat signal;
then, taking the minimum value point of the pulse wave signal closest to the minimum value point of the heartbeat signal as a segmentation point, and segmenting the single-beat pulse wave signal;
finally, calculating a correlation coefficient between the single-beat heartbeat signal and the single-beat pulse wave signal, and selecting a single-beat pulse wave signal S4 with the correlation coefficient larger than 0.6;
the correlation coefficient calculation formula is:
Figure FDA0003870343720000011
wherein g is p For a single beat of pulse wave signal, g h For a single beat heartbeat signal, cov () is the covariance formula, σ 2 () Is a variance formula;
finally, setting a blood pressure prediction neural network, standardizing the single-beat pulse wave signal S4, and inputting the standardized single-beat pulse wave signal into the blood pressure prediction neural network for training; and inputting a new single-beat pulse wave signal S4 extracted by the subject after the training is finished, and inputting after standardization to directly obtain a blood pressure result.
2. The method for non-contact blood pressure measurement based on the pulse ultra-wideband radar as claimed in claim 1, wherein the two-dimensional radar echo signal matrix
Figure FDA0003870343720000012
Wherein the row vector represents a fast time dimension and is positively correlated with the detection distance; the column vector represents the slow time dimension and is positively correlated with the data accumulation time, and each element x of the matrix ij Representing radar signal sample values.
3. The method for measuring the blood pressure in the non-contact mode based on the pulse ultra-wideband radar as claimed in claim 1, wherein the pre-processing of the two-dimensional radar signal matrix M1 comprises:
firstly, segmenting each row of a radar signal matrix M1 according to a specific application scene, taking an average value of each segment, subtracting the corresponding average value of each segment, and removing a direct current component;
then, filtering each line of data of the radar signal matrix with the direct-current component removed by using a band-pass filter;
finally, removing static clutter from the filtered radar signal matrix through a moving average algorithm to obtain a radar signal matrix M;
the static clutter of each row is the weighted sum of the static clutter of the previous row and the data of the current row, and is expressed as:
C(t,τ)=a·C(t-1,τ)+(1-a)·x(t,τ)
wherein C (t, tau) is the static clutter of the current row, a is the weight, C (t-1, tau) is the static clutter of the previous row, x (t, tau) is the data of the current row, and the static clutter of the first row is 0.1. X (1, tau) during initialization.
4. The method for non-contact blood pressure measurement based on the pulse ultra-wideband radar is characterized in that the energy value of each column is calculated from the radar signal matrix M according to the column, and the energy calculation formula of the jth column signal is as follows:
Figure FDA0003870343720000021
where m is the signal length.
5. The non-contact blood pressure measuring method based on the pulse ultra-wideband radar according to claim 1, wherein the specific process of obtaining the heartbeat signal S2 and the pulse wave signal S3 is as follows:
firstly, obtaining modal signals with different center frequencies from a vital sign signal S1 through variational modal decomposition, carrying out fast Fourier transform to calculate the frequencies of the modal signals, and selecting a standard heart rate value;
the formula of the variation modal decomposition is as follows:
Figure FDA0003870343720000022
wherein [ mu ] i And { omega } and i the i-th modal component and the center frequency after decomposition are respectively corresponding to the f-th modal component and the center frequency, delta (t) is a Dirac function, beta is a secondary penalty factor, gamma is a Lagrangian multiplier operator, and f (t) is an original signal;
then, selecting a modal signal with the center frequency being a standard heart rate value as a heartbeat signal S2;
and finally, except the modal signal with the minimum energy, selecting a heartbeat signal and a modal signal with the energy smaller than the heartbeat signal for superposition, and regarding the superposition as a pulse wave signal S3.
6. The method for non-contact blood pressure measurement based on the pulse ultra-wideband radar is characterized in that the blood pressure prediction neural network comprises the following structures:
firstly, taking a single beat pulse wave signal S4 after standardization as a sample, and inputting two branches to carry out convolution with different scales;
the structure of the first branch is sequentially a 5 multiplied by 1 one-dimensional convolution layer, a BN layer, a maximum pooling layer, a 5 multiplied by 1 one-dimensional convolution layer, a BN layer, a 1 multiplied by 1 one-dimensional convolution layer and a maximum pooling layer; the structure of the second branch is sequentially a 3 multiplied by 1 one-dimensional convolution layer, a BN layer, a maximum pooling layer, a 3 multiplied by 1 one-dimensional convolution layer, a BN layer, a 1 multiplied by 1 one-dimensional convolution layer and a maximum pooling layer;
then, splicing the feature maps obtained by convolution of the two branches, sequentially passing through a Dropot layer, two layers of gating circulation units, a multi-head attention module, a Dropot layer and two layers of full connection layers, and outputting a prediction result;
the Dropout layer is used for preventing overfitting, the gating circulation unit and the multi-head attention module can capture deep timing characteristics of signals, and the full connection layer can map the characteristics to output blood pressure values.
CN202211194354.5A 2022-09-28 2022-09-28 Non-contact blood pressure measurement method based on pulse ultra-wideband radar Active CN115607126B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (6)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
田泽懿 等: "基于射频信号的非接触式血压监测系统", 《光学仪器》, pages 81 - 87 *

Cited By (3)

* Cited by examiner, † Cited by third party
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