Front. Physiol. Frontiers in Physiology Front. Physiol. 1664-042X Frontiers Media S.A. 1206471 10.3389/fphys.2023.1206471 Physiology Original Research A novel target state detection method for accurate cardiopulmonary signal extraction based on FMCW radar signals Zhang et al. 10.3389/fphys.2023.1206471 Zhang Xiaozheng 1 † Yang Chenxi 1 † Xiao Zhengyu 2 Lu Binbin 2 Zhang Ji 2 Li Jianqing 1 * Liu Chengyu 1 * 1 The State Key Laboratory of Digital Medical Engineering, School of Instrument Science and Engineering, Southeast University, Nanjing, China 2 Chuhang Technology Co. Ltd., Nanjing, China

Edited by: Hao Gao, University of Glasgow, United Kingdom

Reviewed by: Rishi Raj Sharma, Defence Institute of Advanced Technology (DIAT), India

Ewaryst Janusz Tkacz, Silesian University of Technology, Poland

 *Correspondence: Jianqing Li, ljq@seu.edu.cn; Chengyu Liu, chengyu@seu.edu.cn

These authors have contributed equally to this work and share first authorship

 26 06 2023 2023 14 1206471 15 04 2023 08 06 2023 Copyright © 2023 Zhang, Yang, Xiao, Lu, Zhang, Li and Liu. 2023 Zhang, Yang, Xiao, Lu, Zhang, Li and Liu

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Frequency-modulated continuous wave radar is capable of constant, real-time detection of human presence and monitoring of cardiopulmonary signals such as respiration and heartbeat. In highly cluttered environments or when the human body moves randomly, noise signals may be relatively large in some range bins, making it crucial to accurately select the range bin containing the target cardiopulmonary signal. In this paper, we propose a target range bin selection algorithm based on a mixed-modal information threshold. We introduce a confidence value in the frequency domain to determine the state of the human target and employ the range bin variance in the time domain to determine the range bin change status of the target. The proposed method accurately detects the state of the target and effectively selects the range bin containing the cardiopulmonary signal with a high signal-to-noise ratio. Experimental results demonstrate that the proposed method achieves better accuracy in cardiopulmonary signal rate estimation. Moreover, the proposed algorithm is lightweight in data processing and has good real-time performance.

 FMCW radar target state detection range bin cardiopulmonary signal respiration rate heart rate section-at-acceptance Computational Physiology and Medicine 1 Introduction

Respiration and heartbeat can indicate the basic physiological health status of the human body and are also adopted for clinical diagnosis (Gu, 2016; Weenk et al., 2017). However, traditional vital sign detection methods, such as electrocardiography or photoplethysmography, are contact-based and require electrodes to be in contact with the body surface, which can cause discomfort with long-term use (Castaneda et al., 2018; Kebe et al., 2020). Non-contact radar for vital sign detection is a promising approach that offers increased flexibility in usage. Radars can be placed anywhere in a given space and provide continuous monitoring without the patient’s perception (Morgan and Zierdt, 2009; Obeid et al., 2012; Will et al., 2016; Yang et al., 2020). It overcomes the discomfort of the long-term use of electrodes in wearable sensors and is of great value to newborn babies, burned patients, and infectious disease patients in hospitals (Zhao et al., 2013). The remote monitoring of vital signs of the radar system is more and more widely used. Radar can detect human presence in indoor environments (Mercuri et al., 2013; Gennarelli et al., 2022). Radar can track the periodic movement of body parts caused by breathing, which can be used for system diagnosis of breathing problems such as apnea (Sun, 2019; Islam et al., 2021; Ishrak et al., 2023). Radar sensing technology can also perform gesture recognition and classification for human-computer interaction (Sharma et al., 2023). The latest research shows that radar systems can also accurately extract cardiac physiological parameters for estimating heartbeat intervals and cardiac timing (Xia et al., 2021). Its application scenarios also incorporate the detection of diseases such as sleep (Hong et al., 2018); monitoring of the driver’s health status (Park et al., 2019; Wang et al., 2021); earthquake rescue detection (Liu et al., 2011), etc.

The radars currently used for vital signs detection include frequency-modulated continuous wave (FMCW) radar (He et al., 2017; Ahmad et al., 2018), ultra-wideband (UWB) radar (Deng et al., 2019; Yang et al., 2021; Sharma et al., 2023), continuous wave (CW) radar (Hong et al., 2018; Gennarelli et al., 2022), etc. The FMCW radar has a wide range of applications for continuous non-contact vital sign monitoring in a variety of scenarios. It provides high-ranging accuracy and can detect micro-displacement of the abdomen and chest caused by cardiopulmonary activity. Low-cost low-power radar sensors have advantages in continuous time cardiopulmonary signal monitoring (Suijker et al., 2014). The radar intermediate frequency (IF) signal undergoes Fast Fourier Transform (FFT) to obtain micro-displacement information in each range bin. For cardiopulmonary signal monitoring, it is necessary to select the range bin where the target is located to extract the micro-displacement signal as the cardiopulmonary signal (He et al., 2017; Lee et al., 2019). Then, the respiration rate (RR) and heart rate (HR) are extracted separately from the cardiopulmonary signal. Therefore, the accuracy of respiration rate and heart rate estimation can be improved from two aspects. One is the precise selection of range bins to obtain cardiopulmonary signals with high signal-to-noise ratio (SNR). The other is the separation algorithm of heartbeat and respiration, which eliminates more noise signals.

Firstly, the accurate selection of the range bin is crucial to extract the cardiopulmonary signal with high SNR. The traditional range bin selection method is based on the maximum amplitude or maximum phase (Li et al., 2013; Van Loon et al., 2016; Alizadeh et al., 2019; Sacco et al., 2020), which cannot avoid the large noise caused by the random motion of the human body or other small displacements, so the accuracy of the cardiopulmonary signal may be affected. In (Choi et al., 2020), the selection method based on the magnitude-phase coherence index improves the range bin selection accuracy, but it requires a lot of computing power and lacks real-time evaluation performance.

Then, various signal processing algorithms are employed for heartbeat and respiration signal separation, each with its advantages and limitations (Singh et al., 2020). The most straightforward way to extract RR and HR is to use an FFT on a segment of the phase signal and find the peaks in the spectrum. However, the leakage problem due to the limited data length of FFT leads to degraded detection performance. Both the FFT algorithm and the autocorrelation function are applied to the received signal to extract the cardiopulmonary signal (Chioukh et al., 2011). In (Lee et al., 2016), the authors addressed the FFT smearing and leakage problem by using the multiple signal classification algorithm. In (Lv et al., 2021), an algorithm combining complete ensemble empirical mode decomposition with adaptive noise and fast independent component analysis was employed to investigate short-term heart rate measurement techniques. In order to obtain higher estimation accuracy, these methods increase the complexity of the algorithm and require more computing power.

There are currently few research algorithms for cardiorespiratory monitoring at the low sampling rate of low-power radars. To meet the real-time and lightweight requirements of the algorithm, this paper proposes a target state detection method based on mixed-modal thresholds that accurately selects the target range bin and extracts cardiopulmonary signals with high SNR. In the frequency domain, the target energy confidence value is compared to a threshold to precisely identify the target state. In the time domain, the range bin variance is compared to a threshold to select the optimal range bin in each data frame. For respiration rate and heart rate estimation, this paper utilizes the FFT algorithm, the most commonly used spectral analysis method. While assuming a static human body relative to the radar, this approach provides acceptable accuracy and requires minimal computing power, making it suitable for embedded systems with low-power radars. Furthermore, the proposed method demonstrates robust performance even with lower sampling rate data. The main contributions of this work are summarized as the following:

The proposed method performs well for human target state detection and can select range bins containing cardiopulmonary signals with high SNR. We estimate respiration rate and heart rate with the FFT method. The results demonstrate that the proposed method is more accurate than traditional range bin localization methods.

The proposed method only needs lightweight calculations and has certain adaptability and real-time performance under low sampling rate data, which is extremely valuable for realizing low-cost low-power radar chip embedded systems.

The rest of this paper is organized as follows. Methods are introduced in Section 2. The experimental results are described in detail in Section 3. In Section 4, the conclusion is drawn.

 2 Methods

The process of radar signal processing is presented in Figure 1, which consists of three parts: signal pre-processing, target state detection, and cardiopulmonary signal processing. The innovative part of our proposed algorithm is mainly the target state detection part.

The overall signal processing flow. (A) Pre-processing. (B) Target state detection flow. (C) Cardiopulmonary signal processing.

 2.1 FMCW radar signal model in cardiopulmonary measurements

Any object in the radar path will reflect the echo signal, and the radar system can detect the distance, speed, and angle of the target by capturing and processing the reflected signal. The experimental setup for radar cardiopulmonary signal monitoring is presented in Figure 2A. Assume that in a chirp period T c , the frequency of the radar transmit signal f t and the transmit wave signal s T t are denoted as: f t = f c + γ t , 0 ≤ t < T c s T t = A T ⁡ exp ⁡ ⁡ j 2 π f c t + γ t 2 / 2 , 0 ≤ t < T c where f c is the carrier frequency, γ is the slope of the chirp and indicates the change rate of the modulation frequency over time, t is the sweep time, T c is the chirp period, A T is the transmit signal power, and j is an imaginary unit. The radar received signal s R t corresponds to a delay in the transmit signal of time τ can be described as (Alizadeh et al., 2019): s R t = A R ⁡ exp j 2 π f c t − τ + γ t − τ 2 / 2 , τ ≤ t < T c τ = 2 R 0 + x t / c where A R represents the received signal power, τ is the received time delay, R 0 stands for the initial distance when the human body remains relatively stationary from the radar, x t is the micro-displacement of the body surface caused by respiration and heartbeat, and c is the speed of light. The transmit signal and the receive signal overlap on τ , T c , and the IF signal s I F t is obtained by inputting these two signals into the mixer. The definition of s I F t is expressed as: s I F t = s T t s R * t ≈ A T A R ⁡ exp j 2 π γ τ t + f c τ = A T A R ⁡ exp j 2 π f I F t + φ t , τ < t ≤ T c where the term τ 2 ( τ ≪ 1 ) is neglected, f I F is the frequency of the IF signal, and φ t is the phase of the IF signal, as drawn in Equation (6) (Alizadeh et al., 2019). φ t = 2 π f c τ = 4 π R 0 + x t / c

Radar cardiopulmonary monitoring. (A) Experimental setup. (B) FMCW radar signal model.

An FFT calculation is performed on each chirp sample of the radar-acquired IF signal, this FFT reveals range information, called the Range FFT. Each chirp (time axis) is sampled N times (range axis), and a distance block (range bin) is generated for each sample. Each Range FFT bin represents a specific distance, and this process is repeated for each chirp, so it will produce a Range FFT block of N × M dimensions, as displayed in Figure 2B where N is the number of FFT points and M is the number of chirps per frame. We need to locate the range bin where the target is located, and then extract the cardiopulmonary signals from that range bin.

In Figure 1A, chirp integration is performed for each Frame. The average value of n chirps in each Frame is used as the Chirp selected by this Frame. This operation assists to eliminate phase noise. Only one chirp is selected per frame, which reduces the sampling rate. It also reduces the amount of data processing calculations. We propose a lightweight algorithm for subsequent processing.

The proposed target state detection algorithm.

The proposed target state detection algorithm based on mixed-modal information is demonstrated in Figure 1B, which consists of three procedures: static clutter removal, confidence value calculation, and range bin variance. Static clutter removal is used to remove static background noise in the signal. The confidence value calculation is used to determine whether the target exists. The range bin variance is used to determine whether the target is stable or moving after detecting the presence of the target. The proposed algorithm can accurately detect three human target states.

1) Stable state is a stable human target state. It indicates that the range bin containing cardiopulmonary signal information with high SNR can be accurately selected;

2) Motion state is an intermediate state. It demonstrates that a human target can be detected. However, the human body has relatively large movements, with a large noise signal, and it is difficult to accurately extract cardiopulmonary signal information;

3) Unmanned state means that no human target is detected.

When the human target is in a stable state, the range bin is selected, and the cardiopulmonary signal with a higher SNR is extracted. Then we employ a simple FFT spectral analysis method to estimate the cardiorespiratory rate.

2.1.1 Static clutter removal

In real environments, radar signals regularly contain a lot of noise. Radar receives reflected signals not only from the target but also from the environment and unwanted targets. Before locating potential live targets and extracting weak cardiopulmonary signal signals, static environmental noise needs to be removed from the data.

Then perform static clutter removal on the Range FFT information to remove more static points from the background. The model is expressed in Equation 8. X m ′ k = α X m k + 1 − α X m − 1 ′ k Y m k = X m k − X m ′ k where m is the frame index, k is the range bin index, α is the filter coefficient, the larger the α value the stronger the clutter filtering effect, X is the Range FFT value, which is a complex number, X ′ is the calculated filtering parameter, and Y is the filtered Range FFT value.

Calculate the Range FFT amplitude of each frame after clutter removal, and obtain the range profile matrix containing the instantaneous distance information. Figure 3B manifests the range profile of the slow time dimension. The human target is located near the 25th range bin. The distance between the human and the radar is about 1.5 m.

Target range bin localization. (A) Range FFT for a single frame. Some static clutter is suppressed after clutter removal. (B) Range profile across slow time.

 2.1.2 Confidence value calculation

Based on the set range bin scan range, the target range bin is searched from the range bin start index to the range bin end index. Constant False Alarm Rate (CFAR) (Kronauge and Rohling, 2013) is a technique for radar systems to detect the presence of target signals. We refer to this method to design a human target detection method based on confidence value. The range bin where the human target is located has much larger energy than other bins, and the range bin with the peak energy in the range profile is selected as the detection cell. The target confidence value is the energy ratio of detection cells to training units, as in Equation 9. V m = P D m / P T m P T = ∑ i = D − G − T D − G P i / T + ∑ j = D + G D + G + T P j / T where V is the confidence value, m is the frame index, P D is the peak energy of detection cells, P T is the average energy of training cells, D is the range bin index of peak energy, G is the number of single-side guard cells, and T is the number of single-side training cells.

Each window consists of cells as expressed in Figure 4, and the detection cells are the range bins of peak energy. Guard cells are on both sides of the detection cells and are used to prevent target signal energy from leaking into the training cells. When the target has a large random motion, some noise signal energy will enter the guard cells. Training cells are background noise energy cells. The number of guard cells and training cells is determined according to the environmental noise. Confidence calculations are performed on the range profile of each frame to detect targets. The edge range bin is processed in a circular supplementary manner, which can keep the same number of guard cells and training cells on both sides of the detection cells. Copy the range profile three times and merge it to generate a new matrix, and do calculations on the middle range profile. The cells that are not enough on the left are supplemented with the cells at the end of the first range profile, and the cells on the right are supplemented with the cells at the front of the third range profile.

Confidence value calculation window.

There is an error in the determination of a single frame due to the influence of background noise and random body movements. A method of doing a cumulative calculation of V value data of v frames window is used, which can reduce the error, as in Equation 11. V s u m m = ∑ i = m m + v V i

Figure 5 presents the effect of confidence value accumulation on human presence detection. The false detection rate decreases after accumulation.

Effect of confidence value accumulation on human presence detection. (A) before the accumulation. (B) after the accumulation.

In our work, we set a sliding window of 2 s (40 frames) and a step size of 1 frame for human presence detection, which provides real-time detection. The parameter of V s u m greater than the set threshold means that the human target is detected. Otherwise, noise signals are detected and no human target signals are detected. The threshold is trained according to the background scene without people, which has a certain scene adaptive performance. To avoid target signals with small amplitudes being classified as noise, set the threshold lower.

 2.1.3 Range bin variance

In our proposed algorithm, the target should be in a stable state to extract the cardiopulmonary signal accurately. We introduce the variance s 2 to analyze the dispersion of the range bin in the time domain. After detecting a human target, it can be used to determine whether the target is in stable state or motion state. s 2 = ∑ i = 1 n r i − r 2 / n − 1 where s 2 is the range bin variance, n is the number of frames, r i represents the index of the range bin of the i^th frame, and r stands for the mean value of the range bin of the n frames.

Figure 6B shows that human objects are detected in the 60-s time. we set a sliding window of 2 s range bin variance buffer updated with 1 frame step size, which stores the state of the target every frame. After a human target is detected at each frame, the target state is determined by comparing the variance with a set threshold. The variance is less than the threshold, indicating that the human target is in stable state. If the variance is greater than the threshold, the human target is in motion state. Ideally, the human body should not move very much during the 1-frame period, so the range bin position and the current state will not change. When the proportion of any state in the buffer exceeds 80%, it is determined that the frame is in that state.

(A) Range profile across slow time. (B) Confidence value calculation results. (C) Range bin variance. (D) Target state detection.

The time in the first 20 s of Figure 6C indicates that the range bins of the human body were detected to be cluttered, which means that the person is in a state of relatively large motion. The human body may have a large movement with excessive noise signals, which is not suitable for cardiopulmonary signal extraction. In the last 40 s, the human body is in a stable sitting state, which is stable state. We perform median smoothing of the 2 s window on the target range bin of stable state. Then select the range bin where the target is located, and extract the cardiopulmonary signal information. Our algorithm can make a target state decision within 2 s. The selection of range bins also has the same real-time performance.

 2.2 Cardiopulmonary signal processing

Cardiopulmonary signals are processed using a sliding window of 20 s and a stepping update of 1 s to estimate respiration rate and heart rate. Cardiopulmonary signal processing consists of the following parts.

1) Range bin selection: After using the proposed target state detection algorithm, select the range bin where the target is located. And extract the cardiopulmonary signal phase and amplitude information of the target.

2) Phase signal processing: Phase unwrapping recovers continuous phase changes. Then do the phase difference processing, which helps to eliminate the phase drift caused by the noise signal.

3) Separation of respiration and heartbeat: Respiration and heartbeat signals occupy different frequency ranges, which can be separated with appropriate band-pass filters. A fourth-order IIR cascaded biquad filter limits the frequency range to 0.1–0.6 Hz for respiration detection and 0.8–4 Hz for heartbeat detection. Figure 7 exhibits the comparison between the extracted heartbeat signal and the ECG collected by the wearable sensor.

4) Respiration rate estimation: FFT spectrum analysis is performed on the separated respiration signal. The respiration rate is selected by the largest peak within the spectrum. It can also be estimated by calculating the distance between the peaks of the time-domain waveform.

5) Heart rate estimation: FFT spectrum analysis is performed on the separated heartbeat signal. Respiratory harmonics will overlap with the spectrum of the heartbeat signal, which will affect the heart rate estimation, as shown in Figure 8. Where f r is the respiratory rate and f h is the heart rate. Extract N peaks (from large to small) in the frequency spectrum of the heartbeat signal, and eliminate the peaks corresponding to the respiratory harmonics. Then, the heart rate is selected by the largest peak within the spectrum.

Comparison of extracted heartbeat signal by radar with ECG.

Extracted cardiopulmonary signal (after unwrapping and difference) and its spectrum.

The reference data collected by the wearable sensor is processed in the time domain, and the respiration rate and heart rate are estimated through a sliding window of the same size.

 3 Experimental results 3.1 Experimental setup

The measured experimental scene is presented in Figure 2A. The FMCW radar has 1 Rx antenna and 2 Tx antennas. It has two receiving channel data, and we only select one of them for processing. The operating frequency of the radar is 23.5GHz–27.5 GHz. Each frame period is 50 ms, and each frame has 32 chips. The reference value data was collected using the Shimmer3 wearable wireless sensor, which is considered a reliable sensor for measuring vital rates, with a sampling frequency of 250 Hz.

There were 20 subjects tested. Each subject was fully informed of the experimental procedure and purpose and agreed to share measurement data anonymously. The patient experimental protocols of the First Affiliated Hospital of Nanjing Medical University were approved by the Ethics Committee (Reference code: 2020-SRFA-183). Each subject had a normal vital rate and rested for at least 10 min before the recording procedure. The radar was fixed in front of the subject at the same chest height to capture chest displacement. The subjects sat on a chair about 1.0 m away from the radar, and each subject was measured for 3 min.

 3.2 Experimental results

To evaluate the work of the proposed method, we tested three states of the human target. The method can make an accurate decision on the target state within 2 s, as displayed in Figure 9A. In order to accurately extract the cardiopulmonary signal of the target, we test the person in a resting state and let the target leave the test area after a certain period. Figure 9B exhibits the results obtained before and after using the proposed method. Other objects and random background noise can interfere with range bin extraction. Before the proposed algorithm, we use the range bin selection method based on the maximum amplitude as a comparison method. Other objects and random background noise interfere with range bin extraction, so range bin extraction is accompanied by random jitter. After applying the proposed method, the smooth range bin information can be accurately obtained.

Human target localization. (A) Range bin selection before and after the proposed algorithm. (B) Target state detection.

We performed spectral analysis using the FFT method on the extracted cardiopulmonary signal with a sliding window of 20 s and a step size of 1 s. The reported beats per minute (BPM) for respiration and heartbeat are averaged over 1-min segments. The mean absolute error (MAE) and root mean square error (RMSE) for all 5 subjects who participated in the experiments are compared with their related reference signals. They are represented as: M A E = 1 n ∑ i = 1 n y i − x i = 1 n ∑ i = 1 n e i R M S E = 1 n ∑ i = 1 n y i − x i 2 where y i is the radar-measured value, x i is the reference value and e i is their absolute error.

One of the comparison results of heart rate estimation is shown in Figure 10. The average heart rate was 75.0 BPM estimated by the wearable sensor. Before applying the proposed method, the estimated average heart rate was 67.5 BPM, the overall MAE was 7.4 BPM, and the RMSE was 7.4 BPM. After applying the proposed method, the estimated average heart rate is 75.6 BPM, the overall MAE is 1.5 BPM, and the RMSE is 0.6 BPM. Due to the influence of small body random motion and background noise, the range bin selection will have a large error, which leads to an increase in the cardiopulmonary signal estimation error. After applying the proposed algorithm, the accuracy of range bin selection is improved, and the signal-to-noise ratio of the extracted cardiopulmonary signal is improved. It can be seen from Figure 10 that the estimated heart rate values are in good agreement with the reference heart rate after applying the proposed algorithm. The results show that the range bins selected by this method contain cardiopulmonary signals with higher SNR.

Comparison of heart rate estimates before and after applying the algorithm and the reference value.

The results of the heart rate and respiration rate estimation errors before and after applying the proposed method are displayed in Table 1. The results were measured with 5 subjects with different situations in this experiment. The results give the average MAE and average RMSE of 1.0 BPM and 0.5 BPM for respiration rate estimation, and 1.9 BPM and 1.3 BPM for heart rate estimation after applying the proposed method. The results are all better than the comparison method. Before applying the algorithm, the heart rate estimation errors of targets A and D are both large. Because the background noise of the test environment and the random movement of the body are large, it interferes with the range bin selection of the human target. Therefore, the SNR of the extracted cardiopulmonary signal is low, and the estimation error of the cardiopulmonary rate is large.

Measurement error results before and after applying the proposed algorithm. (MAE ± RMSE, BPM).

Subject Respiration Heartbeat
Before After Before After
A 3.6 ± 3.1 1.5 ± 0.9 10.3 ± 10.3 2.8 ± 2.3
B 0.8 ± 0.3 0.8 ± 0.3 2.8 ± 1.1 2.7 ± 1.0
C 0.2 ± 0.1 0.2 ± 0.1 1.3 ± 1.0 1.2 ± 1.0
D 1.2 ± 0.9 1.1 ± 0.8 7.4 ± 7.4 1.5 ± 0.6
E 1.9 ± 1.6 1.3 ± 0.3 1.8 ± 1.7 1.5 ± 1.4
Average 1.5 ± 1.2 1.0 ± 0.5 4.7 ± 4.3 1.9 ± 1.3

Respiration rate accuracy was defined as the proportion of instances where the radar estimated respiration rate was within 2 BPM of the reference sensor, and heart rate accuracy was defined as within 5 BPM. Our experiments demonstrate an accuracy of 90.8% for respiration rate and 94.1% for heart rate. This result is better than the respiration rate accuracy of 83.7% and the heart rate accuracy of 70.7% reported by (Choi et al., 2020). Therefore, the proposed algorithm can accurately select the range bins containing cardiopulmonary signals with high SNR, and the accuracy performance is better in the respiration rate and heart rate estimation.

 4 Conclusion

In this paper, we propose a target state detection algorithm based on a mixed-modal threshold for target state detection and cardiopulmonary signal range bin selection. First, the algorithm can accurately detect the target state in different scenarios in real time. It can avoid range bins with larger noise signals and select range bins containing cardiopulmonary signals with high SNR, which improves the accuracy of respiration rate and heart rate estimation. In the experimental results, we analyze the radar recordings of 5 subjects in different scenarios and different motion states. We select range bins and extract cardiopulmonary signals by the proposed algorithm, and estimate respiration rate and heart rate by the FFT method. The results demonstrate that the proposed algorithm has better performance than traditional range bin selection methods in terms of cardiopulmonary signal estimation. Furthermore, it performs well on low sampling rate data. it requires less computing power and high performance for real-time evaluation, which is of great value for forming a lightweight embedded system on a low-cost low-power radar chip.

Nevertheless, this work also has some limitations. During the experiment, the subjects were required to remain relatively still. Because the vibration noise of the environment with a frequency close to the heart rate will interfere with the heart rate estimation. In future research, we intend to conduct embedded system testing experiments in various practical situations such as home care, medical monitoring, driver monitoring systems, etc.

 Data availability statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding authors.

 Ethics statement

The studies involving human participants were reviewed and approved by The patient experimental protocols of the First Affiliated Hospital of Nanjing Medical University were approved by the Ethics Committee (Reference code: 2020-SRFA-183). The patients/participants provided their written informed consent to participate in this study.

 Author contributions

XZ, CY, JL, and CL contributed to the study design, data acquisition, and data analysis. ZX, BL, and JZ provided critical assessments of data and provided resources. XZ drafted the article under the supervision of CY, JL, and CL, and All authors contributed to the article and approved the submitted version.

 Funding

This work was supported by the National Natural Science Foundation of China (62101120, 62171123, 62071241) and the Natural Science Foundation of Jiangsu Province (BK20210208, BE2022160). This work was partly supported by Young Elite Scientists Sponsorship Program by CAST (2022QNRC001) and the Fundamental Research Funds for the Central Universities (2242022R40026).

 Conflict of interest

Authors ZX, BL, JZ were employed by Chuhang Technology Co. Ltd.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 References Ahmad A. Roh J. C. Wang D. Dubey A. (2018). “Vital signs monitoring of multiple people using a FMCW millimeter-wave sensor,” in 2018 IEEE Radar Conference (RadarConf18), Oklahoma City, OK, USA, 23-27 April 2018 (IEEE), 1450–1455. 10.1109/RADAR.2018.8378778 Alizadeh M. Shaker G. De Almeida J. C. M. Morita P. P. Safavi-Naeini S. (2019). Remote monitoring of human vital signs using mm-wave FMCW radar. IEEE Access 7, 54958–54968. 10.1109/ACCESS.2019.2912956 Castaneda D. Esparza A. Ghamari M. Soltanpur C. Nazeran H. (2018). A review on wearable photoplethysmography sensors and their potential future applications in health care. Int. J. Biosens. 4, 195–202. 10.15406/ijbsbe.2018.04.00125 Chioukh L. Boutayeb H. Wu K. Deslandes D. (2011). “Monitoring vital signs using remote harmonic radar concept,” in 2011 41st European Microwave Conference, Manchester, UK, 12-14 October 2011 (IEEE), 1269–1272. 10.23919/EuMC.2011.6101869 Choi H. I. Song H. Shin H. C. (2020). Target range selection of FMCW radar for accurate vital information extraction. IEEE Access 9, 1261–1270. 10.1109/ACCESS.2020.3043013 Deng Q. Le J. Barbat S. Tian R. Chen Y. (2019). “Efficient living subject localization and weak vital-sign signal enhancement using impulse radio based UWB radar,” in 2019 IEEE Intelligent Vehicles Symposium (IV), Paris, France, 09-12 June 2019 (IEEE), 777–782. 10.1109/IVS.2019.8814004 Gennarelli G. Colonna V. E. Noviello C. Perna S. Soldovieri F. Catapano I. (2022). CW Doppler radar as occupancy sensor: A comparison of different detection strategies. Front. Signal Process. 2, 12. 10.3389/frsip.2022.847980 Gu C. (2016). Short-range noncontact sensors for healthcare and other emerging applications: A review. Sensors (Basel) 16, 1169. 10.3390/s16081169 He M. Nian Y. Gong Y. (2017). Novel signal processing method for vital sign monitoring using FMCW radar. Biomed. Signal Process. Control 33, 335–345. 10.1016/j.bspc.2016.12.008 Hong H. Zhang L. Gu C. Li Y. Zhou G. Zhu X. (2018). Noncontact sleep stage estimation using a CW Doppler radar. IEEE J. Emerg. Sel. Top. Circuits Syst. 8, 260–270. 10.1109/JETCAS.2017.2789278 Ishrak M. S. Cai F. Islam S. M. M. Borić-Lubecke O. Wu T. Lubecke V. M. (2023). Doppler radar remote sensing of respiratory function. Front. Physiol. 14, 1130478. 10.3389/fphys.2023.1130478 Islam S. M. Fioranelli F. Lubecke V. M. (2021). Can radar remote life sensing technology help combat COVID-19? Front. Commun. Netw. 2, 648181. 10.3389/frcmn.2021.648181 Kebe M. Gadhafi R. Mohammad B. Sanduleanu M. Saleh H. Al-Qutayri M. (2020). Human vital signs detection methods and potential using radars: A review. Sensors 20, 1454. 10.3390/s20051454 Kronauge M. Rohling H. (2013). Fast two-dimensional CFAR procedure. IEEE Trans. Aerosp. Electron. Syst. 49, 1817–1823. 10.1109/TAES.2013.6558022 Lee H. Kim B. H. Park J. K. Yook J. G. (2019). A novel vital-sign sensing algorithm for multiple subjects based on 24-GHz FMCW Doppler radar. Remote Sens. 11, 1237. 10.3390/rs11101237 Lee K. J. Park C. Lee B. (2016). Tracking driver's heart rate by continuous-wave Doppler radar. Annu. Int. Conf. IEEE Eng. Med. Biol. Soc. 2016, 5417–5420. 10.1109/EMBC.2016.7591952 Li C. Lubecke V. M. Boric-Lubecke O. Lin J. (2013). A review on recent advances in Doppler radar sensors for noncontact healthcare monitoring. IEEE Trans. Microw. theory Tech. 61, 2046–2060. 10.1109/TMTT.2013.2256924 Liu L. Xie H. Liu Z. Bagtzoglou A. Barrowes B. Vincitore A. M. (2011). “Numerical simulation of UWB impulse radar vital signal detection at an earthquake disaster site,” in Proceedings of the 1st international conference on wireless technologies for humanitarian relief (ACWR), Kerala, India, December 18–21, 2011, 278–281. 10.1145/2185216.2185291 Lv W. He W. Lin X. Miao J. (2021). Non-Contact monitoring of human vital signs using FMCW millimeter wave radar in the 120 GHz band. Sensors 21, 2732. 10.3390/s21082732 Mercuri M. Soh P. J. Pandey G. Karsmakers P. Vandenbosch G. A. Leroux P. (2013). Analysis of an indoor biomedical radar-based system for health monitoring. IEEE Trans. Microw. Theory Tech. 61, 2061–2068. 10.1109/TMTT.2013.2247619 Morgan D. R. Zierdt M. G. (2009). Novel signal processing techniques for Doppler radar cardiopulmonary sensing. Signal process. 89, 45–66. 10.1016/j.sigpro.2008.07.008 Obeid D. Zaharia G. Sadek S. El Zein G. (2012). Microwave Doppler radar for heartbeat detection vs electrocardiogram. Microw. Opt. Technol. Lett. 54, 2610–2617. 10.1002/mop.27152 Park J. K. Hong Y. Lee H. Jang C. Yun G. H. Lee H. J. (2019). Noncontact RF vital sign sensor for continuous monitoring of driver status. IEEE Trans. Biomed. circuits Syst. 13, 493–502. 10.1109/TBCAS.2019.2908198 Sacco G. Piuzzi E. Pittella E. Pisa S. (2020). “Vital signs monitoring for different chest orientations using an FMCW radar,” in 2020 XXXIIIrd General Assembly and Scientific Symposium of the International Union of Radio Science, Rome, Italy, 29 August 2020 - 05 September 2020 (IEEE), 1–4. 10.23919/URSIGASS49373.2020.9232333 Sharma R. R. Kumar K. A. Cho S. H. (2023). Novel time-distance parameters based hand gesture recognition system using multi-UWB radars. IEEE Sens. Lett. 7, 1–4. 10.1109/LSENS.2023.3268065 Singh A. Rehman S. U. Yongchareon S. Chong P. H. J. (2020). Multi-resident non-contact vital sign monitoring using radar: A review. IEEE Sens. J. 21, 4061–4084. 10.1109/JSEN.2020.3036039 Suijker E. Bolt R. Van Wanum M. Van Heijningen M. Maas A. Van Vliet F. (2014). “Low cost low power 24 GHz FMCW radar transceiver for indoor presence detection,” in 2014 44th European Microwave Conference, Rome, Italy, 08-10 October 2014 (IEEE), 1758–1761. 10.1109/EuMC.2014.6986797 Sun G. (2019). “Non-contact vital sign measurement with medical radar and its clinical applications,” in 2019 3rd International Conference on Recent Advances in Signal Processing, Telecommunications and Computing (SigTelCom), Hanoi, Vietnam, 21-22 March 2019 (IEEE), 90–92. 10.1109/SIGTELCOM.2019.8696113 Van Loon K. Breteler M. Van Wolfwinkel L. Rheineck Leyssius A. Kossen S. Kalkman C. (2016). Wireless non-invasive continuous respiratory monitoring with FMCW radar: A clinical validation study. J. Clin. Monit. Comput. 30, 797–805. 10.1007/s10877-015-9777-5 Wang F. Zeng X. Wu C. Wang B. Liu K. R. (2021). Driver vital signs monitoring using millimeter wave radio. IEEE Internet Things J. 9, 11283–11298. 10.1109/JIOT.2021.3128548 Weenk M. Van Goor H. Frietman B. Engelen L. J. Van Laarhoven C. J. Smit J. (2017). Continuous monitoring of vital signs using wearable devices on the general ward: Pilot study. JMIR mHealth uHealth 5, e91. 10.2196/mhealth.7208 Will C. Shi K. Lurz F. Weigel R. Koelpin A. (2016). “Instantaneous heartbeat detection using a cross-correlation based template matching for continuous wave radar systems,” in 2016 IEEE Topical Conference on Wireless Sensors and Sensor Networks (WiSNet), Austin, TX, USA, 24-27 January 2016 (IEEE), 31–34. 10.1109/WISNET.2016.7444314 Xia W. Li Y. Dong S. Du M. M. Li X. H. Yang T. (2021). Antihypertensive effects of exercise involve reshaping of gut microbiota and improvement of gut-brain axis in spontaneously hypertensive rat. IEEE Trans. Instrum. Meas. 70, 1–24. 10.1080/19490976.2020.1854642 Yang C. Bruce B. Liu X. Gholami B. Tavassolian N. (2020). A hybrid radar-camera respiratory monitoring system based on an impulse-radio ultrawideband radar. Annu. Int. Conf. IEEE Eng. Med. Biol. Soc. 2020, 2646–2649. 10.1109/EMBC44109.2020.9175267 Yang D. Zhu Z. Zhang J. Liang B. (2021). The overview of human localization and vital sign signal measurement using handheld IR-UWB through-wall radar. Sensors 21, 402. 10.3390/s21020402 Zhao F. Li M. Qian Y. Tsien J. Z. (2013). Remote measurements of heart and respiration rates for telemedicine. PLoS One 8, 71384. 10.1371/journal.pone.0071384