Abstract
Breathing and heartbeat are critical vital signs which reflect the health status of human beings. Aiming to accurately measure the vital sign in short time window, a novel signal processing method for Doppler radar vital sign detection is proposed. Firstly, a two-step I/Q mismatch correction method which, respectively, estimates the time invariant phase imbalance and gain ratio of I/Q channels in the calibration step and the direct-current offsets during normal operation has been proposed. By decreasing the number of estimation parameters from 5 to 2, the parameters can be effectively estimated with data distributed over shorter arc lengths. Then, to solve the discontinuity occurred in arctangent demodulation, the displacement information of chest movement is extracted from the calibrated I/Q signals by extended differentiate and cross multiply algorithm. Finally, instead of Fourier transform-based methods which require long time windows to guarantee sufficient frequency resolution, the optimal parameters of respiration and heartbeat are found by the intelligent search of the differential evolution algorithm. The experimental results show that the proposed method can accurately measure respiratory rate and heartbeat rate with a short time window. For the 8 s time window, the mean absolute errors of respiration and heartbeat were 0.52 bpm and 0.79 bpm, respectively, demonstrating its promise in real-time applications.
1. Introduction
Breathing and heartbeat are important vital signs as well as key indicators of health condition, sleep quality [1–3], etc. Most physiological diseases and physical abnormalities usually lead to changes in respiration and heartbeat. Therefore, continuous and accurate monitoring of respiration and heartbeat in daily life is very critical for early diagnosis of disease.
Polysomnography (PSG) [4] is the primary clinical tool for monitoring respiration and heartbeat. However, PSG is usually limited to clinical usage due to the professional installation of electrodes and expensive monitors. Sensors such as thermistor [5], magnetometer [6], accelerometer [7, 8], and gyroscope [9] have also been utilized to provide ubiquitous low-cost vital sign monitoring by detecting the temperature change in the nasal cavity [5] or the chest wall motion [6–9]. However, these sensors need to be in close contact with the human body, which will inevitably cause discomfort. In order to avoid the discomfort, some noncontact approaches using the monocular camera [10], the depth camera [11], the received signal strength of ZigBee [12], the channel state information of Wi-Fi [13], the ultrawideband (UWB) radar [14], and the continuous wave (CW) Doppler radar [15–17] have been proposed. Among them, the continuous wave Doppler radar working at higher frequency has attracted the attention of researchers due to its merits of high sensitivity, antijamming performance, and simple radio frequency circuit structure.
The signal processing of Doppler radar-based vital sign detection is regarded as a baseband phase demodulation problem [18]. Small angle approximation method [19] can be used to extract vital signs when the amplitude of chest movement is much smaller than the carrier wavelength of radar. However, for high frequency radar with short wavelength, the problems of null detection [19] and small-angle limitation [20] are inevitable. The complex signal demodulation technique [21, 22] was proposed for solving the null detection problem reliably with a quadrature receiver. However, the measurement of heartbeat can be easily affected by the interference of respiration. To reduce the interference of respiratory harmonics, a long time window is required to obtain a sufficient frequency spectrum solution [23]. To overcome the limitation of small angle approximation, the arctangent demodulation [24, 25] was proposed. However, it relies on accurate direct-current offset compensation and quadrature channel imbalance elimination. Besides, the extended differentiate and cross multiply algorithm was proposed [26] to solve the phase entanglement problem in arctangent demodulation. The abovementioned methods usually use fast Fourier transform (FFT) to extract the frequency of breathing and heartbeat. However, as the frequency resolution of FFT is inversely proportional to the time window, the time window must be increased to improve the measurement accuracy, which leads to the degradation of real-time performance.
In order to avoid the frequency resolution problem of FFT and improve the real-time performance, we propose a novel signal processing method for Doppler radar vital sign detection. The main contribution of this paper is as follows. (1)To effectively estimate and compensate the radar imbalance from data distributed over short arc length, a two-step correction method has been proposed. Since the phase imbalance and the gain ratio of I/Q channels are time invariant, we only estimate them in the calibration step and treat them as constants in normal operation. As a result, only direct-current offsets of I/Q channels need to be estimated in normal operation, which effectively reduce the requirement of arc length(2)To accurately measure vital signs with data in short time window, the phase demodulation is regarded as a parameter optimization problem, and the vital sign parameters which best fit the displacement of chest movement extracted from the calibrated I/Q signals can be estimated by the differential evolution algorithm. Therefore, the frequency resolution problem of FFT-based methods [21, 27] can be avoided, and the respiratory and heartbeat parameters can be accurately estimated with limited data of short time window(3)We conducted extensive experiments in real environments and compared them with state-of-the-arts to evaluate the performance of the proposed method. The experimental results show that the proposed approach can achieve accurate and robust vital sign measurement with a limited data length from 24 GHz Doppler radar
The rest of the paper is organized as follows. Section 2 describes the I/Q mismatch correction and displacement signal extraction method in detail. In Section 3, the respiratory and heartbeat parameter extraction method based on the differential evolution algorithm is presented. The experiments and results are given in Section 4, and Section 5 concludes the work.
2. Signal Preprocessing
2.1. I/Q Mismatch Correction
Ideally, the locus of error-free measurements from the I/Q channels of Doppler radar is a circle centered at the origin. However, the shape of an ideal circle can be easily altered by errors in practice [28]. where is a constant phase shift determined by the initial phase shift 0 and the distance between the radar and the target object. and are the amplitude of two orthogonal output signals. is the chest wall displacement due to respiration and heartbeat. is the carrier wavelength. and are the direct-current offsets in I/Q channels. is the total residual phase error.
In order to observe the variation of radar error parameters, the hardware-in-the-loop simulation setup shown in Figure 1 is developed to simulate chest wall movement and collect I/Q signals of Doppler radar. Table 1 describes the simulation setup’s components. A carton with a size of cm is installed on the sliding table to provide a reflective surface. The sliding table is driven by a stepper motor that is controlled by a microcontroller simulating the chest wall movement. To improve the accuracy of parameter estimation, the displacement of the sliding table is set to ±6 mm to make the collected data distributed on the whole ellipse. Different obstructions such as box, glass, and cloth are placed between the radar and the measured object. A set of radar parameters are estimated by the least square-based ellipse fitting algorithm [29–31] from collected data of I/Q signals under different distances and obstructions. As can be seen from Table 2, the direct-current offsets and the amplitudes vary greatly with distances and obstructions, while the total residual phase error and the ratio of amplitudes remain nearly unchanged. In order to reduce the number of parameters to be estimated, a two-step I/Q mismatch correction method, which estimates the invariant parameters in the calibration step and the variant parameters during normal operation, is proposed.

The calibration step only needs to be performed once before the normal operation. The least-based ellipse fitting algorithm [29] is used to estimate the 5 parameters with sufficient data. Then, the time invariant parameters and / can be obtained and finally used for I/Q mismatch correction during the normal operation.
Equation (1) can be expressed as
Substituting the initial values , of the I/Q signals collected by radar in Equation (2), we yield
By subtracting (3) from (2), we obtain
Note that and / are the invariant parameters estimated in the calibration step. According to Equation (4), the standard recursive least square algorithm can be applied to estimate parameter from I/Q signals, and then, and can be calculated. Finally, the radar signals can be calibrated by substituting , , /, and in Equation (1)
Figure 2 shows the raw and calibrated I/Q signals of radar. As can be seen from Figure 2, the proposed two-step correction method can effectively correct the distorted ellipse into a standard ellipse with center at the origin.

(a)

(b)
2.2. Displacement Signal Extraction
To solve the discontinuity occurred in the arctangent function when the demodulation exceeds the native codomain range of (-/2, /2), the extended differentiate and cross multiply algorithm [26] is adopted for the demodulation of absolute displacement. The derivative of the displacement information can be expressed as where and are the time derivative of and .
In the digital field, the integration is replaced with accumulation. Hence, the displacement information can be rewritten in the discrete form
3. Respiratory and Heartbeat Parameter Extraction
The chest wall displacement caused by respiration and heartbeat can be modeled as sinusoids with amplitudes and , frequencies and , and initial phases and [32]. where subscripts and indicate the corresponding parameters of respiration and heartbeat, respectively. Instead of FFT [33], we directly estimate the frequencies by a swarm intelligence algorithm to reduce the size of time windows.
As the chest wall displacement can be computed by Equation (7) from the calibrated I/Q signals, we can find the most matched parameters , , , , , and by solving where represents the displacement calculated in the th I/Q signals.
We choose the differential evolution algorithm [34] to solve this nonlinear optimization problem. The differential evolution algorithm shown in Figure 3 is a very simple yet fairly powerful stochastic global optimizer for continuous search domain. Practically, the parameters of respiration and heartbeat can be set up within a reasonable range for the purpose of improving the speed optimization and reducing the calculation amount. In this paper, the amplitude, frequency, and phase range of respiration and heartbeat are set between [0 mm 0.2 Hz 0 mm 1 Hz ] and [6 mm 0.8 Hz 0.8 mm 2 Hz ]. As discussed in Section 4.3, the optimal population size and the number of iterations are set as 80 and 100, respectively.

To sum up, the proposed vital sign extraction method can be implemented with the following four steps. (1)Collect the orthogonal output signals and detected by Doppler radar(2)Use the proposed correction method shown in Section 2.1 to remove phase imbalance, direct-current offsets, and gain error of I/Q signals(3)Obtain the chest wall displacement information by the extended differentiate and cross multiply method shown in Section 2.2(4)Estimate the parameters of respiration and heartbeat by the differential evolution algorithm
Figure 4 shows the comparison of the displacement extracted from radar outputs and the displacement reconstructed using the vital sign parameters estimated by the differential evolution algorithm. As can be seen from Figure 4, the reconstructed displacement is in good agreement with the reference, demonstrating the effectiveness of the proposed algorithm.

4. Experiments and Results
4.1. Experimental Setup and Methodology
To verify the effectiveness of the proposed algorithm, extensive experiments have been conducted in both controlled and real environments. Figure 5 shows the structure of the experimental setup. The amplified signals of radar are sampled by a microcontroller controlled 24 bit analog-to-digital converter with a sampling frequency of 50 Hz. The negative temperature coefficient (NTC) resistor installing in the breathing valve of mask provides the reference frequency of respiration by sensing the temperature change caused by exhaling hot airflow and inhaling cool airflow during respiration, and the pulse sensor provides the reference frequency of heartbeat by sensing the light transmittance change caused by heartbeat. Besides, the microcontroller also controls the movement of sliding table by using the stepper motor driver. All relevant data is sent to the PC via USB interface and processed by MATLAB R2018a.

We compare the proposed method with the complex signal demodulation method [21], the extended differentiate and cross multiply FFT method [27], and the parameterized demodulation method [32]. The complex signal demodulation method [21] combines the signals of I/Q channels into complex signals and performs spectrum analysis through FFT to obtain the frequencies of respiration and heartbeat. The extended differentiate and cross multiply FFT method [27] extracts displacement information from radar outputs and estimates frequencies using FFT. The parameterized demodulation method [32] defines a demodulation operator to obtain a high spectrum energy concentration and addresses the phase demodulation problem as a parameter optimization problem. Therefore, the frequency resolution problem can be avoided. In the following figures, we use CSD, DACM, and PD as abbreviations of the above methods.
4.2. Verification of I/Q Mismatch Correction
To verify the proposed correction method, we compare the estimated displacement from Equation (7) with the reference displacement of the object on the sliding table in the simulation setup. The amplitude and frequency of the slide movement are set to 0.5-3 mm and 0.3 Hz, respectively. Table 3 shows the mean absolute errors (MAE) of estimated displacement under different arc lengths. As the arc length of the ellipse decreases with the amplitude of movements, the estimation of ellipse parameters becomes difficult due to insufficient data incentive. As shown in Table 3, the performance of the least square-based ellipse fitting algorithm [29] is seriously affected by the reduction of arc length, while the proposed two-step correction algorithm is slightly affected. As shown in Figure 6, the proposed algorithm can accurately estimate the displacement of a moving object, even if the amplitude is 0.5 mm. However, the least square-based ellipse fitting algorithm [29] fails to estimate the displacement due to incorrect parameter estimation of an ellipse. The above results show that the proposed algorithm can accurately estimate the parameters of an ellipse with very short arc length and therefore is more suitable for the detection of slight breathing.

4.3. Determination of Differential Evolution Parameters
The population size and number of iterations are the key parameters affecting the accuracy and computational complexity of a differential evolution algorithm. A larger population size and a larger number of iterations help to reduce error but also increase computational complexity. Figure 7 shows the mean absolute errors (MAE) of estimated respiratory and heartbeat frequencies with different population sizes and numbers of iterations. It can be clearly seen that the errors decrease with the increase of population size and the number of iterations. Finally, the population size and the number of iterations are determined as 80 and 100, respectively, after which the downward gradients of errors have become very small.

(a)

(b)
4.4. Experiments in Controlled Environments
During the experiments in a controlled environment, the sliding table simulates the thoracic motion according to the movement mode given by Equation (8) and the parameters shown in Table 4. Figure 8 shows the estimated respiratory rate and heartbeat rate from data 3. The real respiratory and heartbeat frequencies are fixed at 18 bpm and 72 bpm, respectively. We then use a short time window of 16 s to increase the sensitivity. Since the frequency resolution of FFT is inversely proportional to the time window ( bpm), the respiratory and heartbeat frequencies estimated by the complex signal demodulation and the extended differentiate and cross multiply methods are 18.75 () bpm and 75 () bpm, respectively. Therefore, the estimation results of these two methods are far away from the real values. To improve measurement accuracy, the time window needs to be increased. But this reduces the real-time performance. While the frequency resolution of the parameterized demodulation method [32] and our proposed differential evolution (DE) method do not depend on the time window size, both methods can accurately extract the frequency information of respiration and heartbeat with a limited time window.

(a)

(b)
Figure 9 shows the influence of time window size on the mean absolute errors (MAE) of each method for data 1-6. Due to the large amplitude of chest motion caused by respiration, the estimated errors of respiratory rate are smaller than those of heartbeat rate. With the decrease of time window, the accuracy of FFT-based methods [21, 27] decreases dramatically. The parameterized demodulation method (PD) [32] and the proposed differential evolution method (DE) are only slightly affected. Our proposed method performs better than the parameterized method on all time windows. Figure 10 shows the influence of amplitudes ( and ) on the MAE of each method. The time window is 8 s, corresponding to the frequency resolution of 7.5 bpm for FFT. Hence, for respiratory rate of 30 bpm, the respiratory errors of the FFT-based method [21, 27] are 0. As can be seen from Figure 10, the proposed differential evolution method is not affected by the amplitudes of respiration and heartbeat. Its MAEs of respiratory rate and heartbeat rate for and (data 1) are only 0.04 bpm and 0.32 bpm. As a comparison, the MAEs of parameterized demodulation method are 0.46 bpm and 0.48 bpm, respectively. Therefore, the proposed differential evolution method is more suitable for detecting slight breathing.

(a)

(b)

(a)

(b)
4.5. Experiments in Real Environments
As shown in Figure 11, during the experiments in real environments, a subject seated still in a chair and breathed naturally with her chest facing the antennas. At the same time, the subject wore the mask with NTC resistor and the finger pulse sensor shown in Figure 5 to measure the reference signals of respiration and heartbeat, respectively. Figure 12 shows the comparison between the signals measured by the NTC resistor, the pulse sensor, and the displacement demodulated from the 24 GHz Doppler radar. By performing peak-seeking operations on the above signals, the reference frequencies of respiration and heartbeat can be obtained.


(a)

(b)
Figure 13 shows the respiratory rate and the heartbeat rate obtained from the above vital sign extraction algorithms. The obtained parameters are estimated with data from a certain time window. A smaller time window can better reflect the changes in respiratory and heartbeat parameters and improve the real-time performance of the algorithm. In this paper, an overlapping sliding window with a window size of 16 s and a shift interval of 1 s was adopted to estimate the parameters. The total time length of the estimated parameters is 60 s. The reference is the average value of respiratory rate and heartbeat rate obtained in the time window of 16 s. It should be noted that the frequency resolution of FFT corresponding to the time window of 16 s is 3.75 bpm, which makes the results of the complex signal demodulation [21] and the extended differentiate and cross multiply methods [27] far away from the reference and fluctuate greatly. The parameterized demodulation algorithm [32] and the proposed differential evolution algorithm are not affected by the problem. The errors of the proposed differential evolution method ( bpm and 0.43 bpm for respiration and heartbeat) are smaller than those of parameterized demodulation method ( bpm and 0.68 bpm for respiration and heartbeat). As shown in Figure 14, the accuracy of the proposed differential evolution algorithm is also better than other algorithms in different time windows, especially in the case of short time window. These results show that the proposed differential evolution method is more suitable for real-time monitoring of vital signs in short time window.

(a)

(b)

(a)

(b)
5. Conclusions
In this paper, we proposed a novel vital sign signal extraction algorithm to accurately extract the respiratory and heartbeat parameters in short time window. The proposed method firstly corrects the I/Q signals according to the proposed two-step I/Q mismatch correction method. Then, the chest wall displacement can be obtained by the extended differentiate and cross multiply method. Finally, the differential evolution method is used to estimate the parameters of respiration and heartbeat. The comprehensive experimental results demonstrate that the proposed method has higher accuracy than other existing vital sign extraction algorithms in Doppler radar vital sign detection, especially in the case of short time window. Our future work will focus on sleep stage detection based on these measured vital signs and machine learning.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the Zhejiang Provincial Natural Science Foundation of China under Grant No. LY19F010010.