A spectrum estimation approach for accurate heartbeat detection using Doppler radar based on combination of FTPR and TWV

Non-contact heartbeat detection using Doppler radar is extremely valuable for remotely monitoring and medical diagnosis on special occasions. Nevertheless, fast and accurate heart rate (HR) detection endures several challenges due to influential respiration interference and insufficient frequency spectrum resolution. In this paper, a novel heartbeat detection method with a compact Doppler radar is employed to accurately estimate some indicators of HR and heart rate variability. Firstly, a multiresolution analysis approach based on maximal overlap discrete wavelet transform is introduced to accomplish the preliminary separation of respiration and heartbeat. Subsequently, a template matched filter algorithm is further implemented to maximize the enhancement of the concealed heartbeat component and retrieve the heartbeat signal. Eventually, a novel spectrum estimation method combining frequency–time phase regression with time-window-variation technology is proposed to evaluate the real-time HR. It solves serious dominant frequency estimation deviation and insufficient frequency spectrum resolution in short-period time windows. The accuracy and timeliness of our proposed method are validated by 6 sets of experimental data sampled at the actual office. As a result, the HR detection accuracy is up to 99.70% in different-period time windows of 10 s and 92.09% in 3s. In addition, the mean relative error of extracted beat-to-beat intervals in 3 s ranges from 0.76 to 1.02%.

As a kind of radar sensor, frequency-modulated-continuous-wave radar has been widely concerned with remotely monitoring vital signs such as respiration and heartbeat during the last decades [10][11][12]. So far, many quadrature demodulation methods have been presented one after another due to the widespread application of quadrature receiver architecture with I/Q channel demodulation [13]. The small-angle approximation [14] performs well only by presetting the target position. The complex signal demodulation [15] introduces harmonic interference and seriously reduces HR detection accuracy. Besides, the arctangent demodulation [16] is insensitive to the target position and harmonic interference but remains strong nonlinearity and discontinuity. Furthermore, the extended differentiate and cross-multiply [17] algorithm overcomes the shortcomings of phase wrapping and harmonic interference. Nevertheless, an error accumulation problem due to the discrete integration needs to be solved.
One of the challenges for heartbeat detection is accurately estimating HR under the strong interference of respiration and its harmonics. This is because the chest-wall displacements due to heartbeat (0.1-0.5 mm) are far smaller than those caused by respiration (1-12 mm). A succession of representative methods has been executed for cardiopulmonary signal separation and extraction to eliminate respiration interference, such as continuous wavelet transform (CWT) [18,19], ensemble empirical mode decomposition (EEMD) [19,20], singular spectrum analysis [21], and variational mode decomposition [22]. Although the effectiveness of these methods has been demonstrated for the separation and recovery of respiration and heartbeat, there are still some shortcomings. On the one hand, the above methods are remarkably susceptible to the signal-to-noise ratio. On the other hand, their computational cost is greatly increased. Recently, the capability and robustness of recovering the heartbeat signal based on the characteristic that matched filters can maximize the output signal-to-noise ratio have been substantiated in the presence of large-scale random body movements [23]. Furthermore, multiresolution analysis (MRA) proposed in [24] has been widely applied for fault detection [25], radar systems [26], and medical diagnosis [27]. In contrast, the MRA performance in vital signs research using Doppler radar is looking forward to being further validated.
Another challenge is the insufficient frequency spectrum resolution in short-period time windows. Sufficient frequency spectrum resolution can be achieved using longperiod time windows of more than 10 s [20,28,29]. Nevertheless, the calculated HRV results will be the long-period averages in long-period time windows. In other words, the longer the time windows are, the more the losses of heartbeat details will be. Therefore, for HRV detection in frequency domain analysis, radar data processing in shortperiod time windows is required for fast HR acquisition. A time-window-variation (TWV) technique [30][31][32] is introduced to acquire the best time window and improve the HR accuracy in short-period time windows. Furthermore, the frequency-time phase regression (FTPR) [33,34] can assess the HR robustly in low signal-to-noise ratio conditions and significantly advance the frequency spectrum resolution. Yet, a solution demands to be advised in response to serious dominant frequency estimation deviation in short-period time windows.
In this article, we present a fast and accurate approach for non-contact heartbeat detection using Doppler radar. We first demodulate the quadrature radar echo signal to calculate the chest-wall displacement signal. Then, the heartbeat signal is accurately separated and retrieved from the chest-wall displacement signal. Furthermore, an effective and robust spectrum estimation method is employed to estimate the real-time HR. Finally, human subject experiments in an actual office environment will further verify the consistency of HR and BBIs between radar and ECG.

System architecture of Doppler radar
The fundamental principle of vital signs detection using continuous wave Doppler radar is to capture the frequency shift caused by the chest-wall displacement due to respiration and heartbeat. The vital signs detection system using 77 GHz millimeter-wave radar is shown in Fig. 1. Except for the PC and the subjects, the rest are all integrated into the radar development hardware platform to collect radar data containing the chest-wall movement information. The radar data are uploaded to the PC through the serial port for subsequent algorithm processing. In this paper, the millimeter-wave radar positioned in front of the subject at a distance of d 0 is employed to extract the chest-wall displacement signal x(t). The radar transmitted signal can be expressed as where f and ϕ(t) are carrier frequency and initial phase, respectively. Firstly, T(t) is transmitted toward the subject from the transmitter (Tx). Then, the phase of T(t) is modulated by the chest-wall displacement signal x(t) and the reflected signal is generated. The reflected signal captured at the radar receiver (Rx) can be depicted as where c denotes the velocity of radio wave, and = c f denotes the wavelength of carrier. Generally, the reflected signal can be approximated as a T(t) after phase modulation, with a time delayed of 2d 0 c . After passing through the low-noise amplifier (LNA) at the front end of Rx, the baseband signal B(t) can be represented as  where θ = 4π d 0 + θ 0 is the constant phase related to the initial phase shift θ 0 and the distance d 0 , and the residual phase noise is simplified as �ϕ(t) = ϕ(t) − ϕ t − 2d 0 c . Subsequently, the baseband signal B(t) is demodulated by the quadrature mixer (I/Q Mixer), and the in-phase and quadrature signals are expressed as, respectively, Finally, the analog baseband signal is converted into a digital signal through A/D. The radar digital signal is processed by the DSP hardware module, and the subjects' chestwall movement data are output and uploaded to PC through the serial port.

Proposed method
To realize accurate heartbeat estimation with Doppler radar in a complex environment, a concrete implementation scheme of radar data processing is fulfilled, mainly consisting of four parts: quadrature signal demodulation, signal preprocessing, cardiopulmonary signal separation and extraction, and heartbeat spectrum estimation. The proposed algorithm flowchart is presented in Fig. 2, where "Radar Raw Signal" means the data uploaded by the radar development platform to the PC through the serial port. The remaining algorithm processing parts in the block diagram are all implemented on the PC side.

Quadrature signal demodulation
In this paper, we introduce the arcsine demodulation algorithm [35] into demodulating the quadrature signal. Due to discarding the arctangent's conventional method, phase wrapping and error accumulation problems can be avoided for the Doppler radar system with a low sampling rate. According to the above principle of vital signs detection, by discretizing the parameters t, the demodulated chest-wall displacement signal can be expressed as where k is the k-th sampled data of radar raw signal, is the wavelength of Doppler radar, [0] is the initial phase, I[k] and Q[k] are the in-phase and quadrature signal of baseband, respectively, and x[n] is the chest-wall displacement signal.

Signal preprocessing
Vital signs are completely overwhelmed by environmental noise. Signal preprocessing is implemented to remove background noise, eliminate polynomial trends, and enhance vital signs.

Clutter suppression
Assuming that the propagation environment within the radar measurement range is static, except for the subject's chest-wall motion, the constant features of stationary scatterers can be captured by averaging the chest-wall displacement data. Then, the background clutter can be removed by subtracting the average. In addition, the amplitude instability of the transmitting unit caused by thermal noise and time drift will result in a polynomial trend of the chest-wall displacement signal [36]. Because of a severe polynomial trend, the lowfrequency components will distort or even drown the dominant frequency. Hence, it is necessary to eliminate the polynomial trend by subtracting the best fit line from the original chest-wall displacement.

Differential enhancement
Differentiating the chest-wall displacement signal can significantly enhance the highfrequency components in the displacement signal, namely the heartbeat and harmonic frequencies.
In [37], the heart-to-respiration ratio increased by more than 16 times after first-order differential processing, validating the effectiveness of differential enhancement for weak heartbeat components. For the chest-wall displacement x[n] after clutter suppression, its first-order differential can be approximately expressed as

Cardiopulmonary signal separation and extraction
The next work is followed by signal preprocessing to accurately separate and extract the heartbeat signal under the influential interference of respiration and its harmonics.

MRA algorithm
Multiresolution analysis (MRA) establishes the connection between wavelet transform and digital filter. Based on this theory, a MRA approach based on MODWT is employed to decompose the preprocessing chest-wall displacement signal in multiple frequency bands and achieve the preliminary separation of heartbeat and respiration signals in this paper. Wavelet transform is essentially a multiscale analysis that is implemented by scale and shift operations of a mother wavelet ψ(t) . For any function f (t) ∈ L 2 (R) , the central frequency ψ(t) of the mother wavelet is recorded as f 0 ; then, f = f 0 /a , and the time-frequency form of wavelet transform can be depicted as follows where a is a scale factor corresponding to frequency information, and b is a shift factor relating to space-time information. By setting the scale parameter a to change the shape of the window, it is possible to provide good time resolution for fast events such as heartbeat and excellent frequency resolution for slower events such as respiration, which is the multiresolution characteristic of the wavelet transform.
However, the continuous wavelet transform basis function ψ a,b (t) has a large computational redundancy for fast HR detection. Transformed into discrete wavelet transform by discretizing the parameters a and b, the calculation efficiency is exceedingly advanced. Therefore, set a = a i 0 , b = jb 0 a i 0 , then the discrete wavelet can then be expressed as Set a 0 = 2 , b 0 = 1 , it can be reduced to binary wavelet, expressed as Unlike traditional Fourier transform, wavelet transform has wavelet basis functions, such as Haar wavelet, Daubechies wavelet, Morlet wavelet, Symlet wavelet, and Coiflet wavelet. In this paper, the Coiflet wavelet is chosen as the mother wavelet according to the characteristics of the heartbeat signal for three reasons. Firstly, the wavelet function and scale function of the Coiflet wavelet have good symmetry, which can effectively avoid phase distortion. Secondly, the Coiflet wavelet function has a suitable tightly supported length and decays rapidly outside the effective supporting area. Finally, the shape of the Coiflet wavelet is similar to that of the heartbeat signal. The scaling function and wavelet function of the Coiflet5 wavelet are shown in Fig. 3, respectively. In this paper, the Coiflet wavelet is selected as the mother wavelet of the discrete wavelet transform, the vanishing moment order is selected as 5 order, and the decomposition layers number is 6 layers. The MODWT is implemented on the preprocessed chest-wall displacement signal x ′ [n] , and the calculated wavelet transform matrix is denoted as WT (L+1)×N , where L is the number of decomposition layers and N is the number of samples of the displacement signal. Then, the Coiflet5 wavelet is used to perform MRA on WT (L+1)×N , and the calculated MRA matrix is denoted as WT mra , whose different rows correspond to the MRA results at different scales.
Since HR typically varies from 0.8 to 2 Hz, it is supposed that the frequency range corresponds to the range from the k-th line to the l-th line of WT mra . Besides, the energy distribution characteristics of wavelet coefficients are closely related to the frequency characteristics of the signal. Therefore, weighted reconstruction of different frequency bands based on the energy ratio can highlight the heartbeat's fundamental frequency component. The energy of each scale can be depicted as where i denotes the i-th layer of the wavelet decomposition, i = 1, 2, . . . , L + 1 . WT (i) mra [n] denotes the n-th wavelet coefficient of the i-th row. E (i) denotes the energy of the i-th layer. Then, the energy-weighted reconstructed heartbeat signal can be expressed as where WT (i) mra represents the wavelet coefficients matrix of the i-th row, and x mra [n] represents the chest-wall displacement signal after preliminary separation and reconstruction using the MRA algorithm.

TMF algorithm
In the field of signal processing, a matched filter is an optimal filter, and its optimal criterion is the maximum signal-to-noise ratio of the output signal [38]. As for radar systems, matched filters have been universally applied in impulse radars, where the transmitted signal is adopted as the template signal [39]. Still, matched filters are rarely used in the field of vital signs detection with Doppler radar. Template matched filters, through (11) Fig. 3 The function of the Coiflet5 wavelet convolution operation between the input signal and the "template signal, " can maximize the recovery of the signal similar to the template signal hidden in the original input signal.
The key to the template matching filter is the selection of the template signal. In [23], the low-speed segment signal is directly selected as the template signal. Then, the body motion in the template signal is removed by polynomial fitting to obtain the heartbeat template signal. The method difficulties lie in the choice of the polynomial fitting order and low extraction accuracy. Therefore, a template selection way is proposed to avoid respiration interference in this paper. During the experiment, the subject is requested to hold his respiration for a few seconds to derive the chest-wall displacement signal caused by only the heartbeat as the template signal. Subsequently, a convolution operation between the energy-reconstructed displacement signal and the "template signal" is carried out in the time domain. Then, the signal after template matching filter can be calculated by convolving x mra [n] with a conjugated, time-reversed version of the template signal h[n], that is, To intuitively illustrate our proposed separation and extraction algorithm, the discrepancy is investigated by examining the chest-wall displacement waveforms' visualization and corresponding spectra, as illustrated in Fig. 4. The chest-wall displacement waveforms of different processing stages are shown in Fig. 4 on the left, and their corresponding spectra are shown in Fig. 4 on the right. Comparing Fig. 4a, b with Fig. 4c, d, it can be easily seen that the background clutter is effectively suppressed by the preprocessing, and the chest-wall displacement signal due to respiration and heartbeat is significantly enhanced. As displayed in Fig. 4e and f, the displacement amplitude decreased approximately 5 times and the respiration's frequency fundamental frequency component is (13)  eliminated, showing that the preliminary separation of respiration and heartbeat is fulfilled by applying the MRA. Furthermore, as presented in Fig. 4g and h, the displacement signal has a powerful periodicity and its dominant spectrum peak is exceedingly significant, demonstrating that the TMF can recover the heartbeat signal concealed in the respiration signal.

Heartbeat spectrum estimation
The HR detection using Doppler radar is essentially a frequency estimation problem, and the frequency resolution determines the HR estimation accuracy. For a short-period time window of less than 5s, it is difficult for the traditional FFT to estimate HR quickly and accurately because of the lack of frequency resolution.

FTPR algorithm
Based on the fact that the phase varies linearly with time when the signal exists a dominant frequency [33], the FTPR algorithm converts the traditional frequency domain peak detection into a linear regression fitting of the phase and time in the time domain, and the slope of the phase is proportional to the dominant frequency of the signal. In [33], taking the time window T = 10 s as an example, the frequency resolution of FFT is 0.1 Hz, but that of FTPR is 0.01 Hz. Owing to the increased frequency resolution, HR detection accuracy is greatly promoted. Additionally, the heartbeat signal clearly has a dominant frequency, which is the heartbeat frequency. Consequently, the FTPR is an ideal choice. The specific steps of the FTPR for HR detection are as follows

FTPR-TWV algorithm
However, the dominant peak of the heartbeat spectrum will seriously deviate from the true heartbeat frequency in short-period time windows. The reasons are as follows: First, there is still residual harmonic noise in the extracted heartbeat signal. Second, the true heartbeat frequency is not an integer multiple of the frequency resolution due to insufficient frequency resolution. Aiming at the above shortcomings, a novel spectrum estimation method combining FTPR with TWV is proposed in this paper. For a time window of Ts, the frequency resolution f is calculated as �f = 1/T . To measure the HR accurately, one of the following conditions or both should be satisfied: (a) f ≪ HR , (b) HR = k f , where k is an integer number. Obviously, the detection accuracy can be satisfied by increasing the length of the time window, but the performance of the real-time detection is greatly reduced. In addition, for a periodic signal with a frequency of f, if the length of the time window is not an integer multiple of the period or in other words, f = m f , the sideband spurs or the harmonics will come out after the FFT, and the fundamental frequency energy will be degraded [30]. The closer the ratio is to an integer multiple, the smaller the fundamental frequency energy leakage is. The specific steps of "Acquire the best time window" are as follows. and draw all the heartbeat spectrograms in the same spectrogram to obtain the spectrum combination. This is the concept of "Spectrum combination. " 4. Based on the fundamental frequency spur theory proposed above, by looking for the maximum value of the spectrum peak in the "Spectrum combination, " the integer multiple of the frequency resolution of the time window corresponding to the maximum value is closest to the HR. This is the concept of "the best time window" proposed in this paper.
According to above steps, the number of time windows is jointly determined by the time interval S and the time variation t and has nothing to do with the original window size T. The more time windows there are, the more frequency resolutions are obtained, the closer the integer multiple of the frequency resolution of the best time window is to the HR, and the higher the detection accuracy. However, the computational cost will increase significantly, and trade-offs are required in the practical applications. In this paper, the time interval as S = 1 f s is selected. An intuitive superiority of the combined spectrum of TWV is visualized in Fig. 5. As an example in the time window of T = 3 s and sampling frequency of f s = 32 Hz, a series of different time windows are generated. As shown in Fig. 5, the spectrum peak of the 3s time window is 0.215, corresponding to a frequency of 1.333 Hz. However, the maximum peak in the combined spectrum is 0.223, corresponding to a time window of 3.179 s and a frequency of 1.255 Hz. Furthermore, because the reference frequency measured by ECG is 1.258 Hz, the time window of 3.179 s is unmistakably the best. Based on the above visual analysis, the severe dominant frequency estimation deviation can be avoided by introducing the TWV. The FTPR-TWV algorithm is illustrated by conducting visual analysis of the specific implementation, as manifested in Fig. 6. This flowchart shows more details of the "Heartbeat Spectrum Estimation" part in Fig. 2. Following the TMF, the samples of the extracted heartbeat signal in the initial time window of T s are selected as input data, and a series of different-period time windows are generated by varying the time window size at a time interval 1/f s from T − t to T + t , where �t < 0.1T . Afterward, a series of the spectra can be acquired by applying FFT to each time window, combining all the spectra into a spectrum. By searching for the largest spectral peak in the combined spectrum, the time window corresponding to this peak is the best. Finally, the HR, closest to the true HR, can be estimated using the FTPR algorithm on the sampled data of the best time window.

Experimental results and discussion
To verify the reliability of the heartbeat detection algorithm proposed in this paper, many experimental datasets are measured under different experimental conditions. Moreover, effective evaluation indicators are adopted to evaluate the performance of the proposed method.

Experimental conditions and parameters
The experiments are completed in an actual office environment, as exhibited in Fig. 7. The HR reference data in this experiment are measured by the traditional contact ECG acquisition module (the lower right corner in Fig. 7). ECG electrodes are attached to the subject's chest near the heart, and the sampling frequency is set to 500 Hz. The subject's ECG data are uploaded to the PC through serial port 2. The subject's ECG data are uploaded to the PC through serial port 2. IWR1443Boost (the upper right corner in Fig. 7) radar platform designed by Texas Instruments is adopted, and the relevant experimental parameters of the Doppler radar are displayed in Table 1. The subjects' chestwall movement data are uploaded to the PC through data serial port 1. The actual office environment makes the experimental data closer to the data in the practical application. As for each experiment, 2 min of data were recorded from the subject 0.5 m away from the Doppler radar, sitting still in front of the radar and breathing normally. The 6 sets of experimental data are the measurement results of two male subjects at different times in the morning, middle and evening. Each set of data contains a 2 min measurement signal. A total of two subjects participate in the vital sign detection experiments. Besides, the time resolution and frequency resolution need to be compromised in practical applications. Generally speaking, the shorter the time window is, the more heartbeat details can be presented, yet the lower the frequency resolution will be, the worse HR detection accuracy will be. Therefore, to prove the universality of the proposed algorithm for different occasions, two different-period time windows, respectively, 3 s and 10 s, are chosen for the HR detection.

Experimental results evaluation
Compared with a band-pass filter (BPF) and EEMD-based algorithm to verify the superiority of the proposed cardiopulmonary signal separation and extraction algorithm, a synchronized data segment of 10 s from the sixth set of data is picked to present the extracted heartbeat signal clearly, respectively. Additionally, the BPF is set as a 10-order Butterworth filter with a passband range from 1 to 3 Hz, and the matched template signal is decided by holding respiration for 5 s in sitting still.
As illustrated in Fig. 8, the blue line represents the reference ECG signal. The red line represents the heartbeat waveform extracted by the incorporation of MRA and TMF. The black line represents the heartbeat waveform extracted by the association of BPF and EEMD. In contrast to the mode mixing of the heartbeat waveform extracted by BPF and EEMD, the heartbeat waveform extracted by MRA and TMF is extremely periodic. Besides, the troughs of the waveform extracted by MRA and TMF basically correspond to the R peak of the ECG signal, indicating the similarity of the real-time HR between radar and ECG.
Based on the superiority of the heartbeat signal extraction algorithm mentioned above, comparing the performances of FFT-TWV, FTPR, and FTPR-TWV algorithms in HR detection, the development of our proposed spectrum estimation method will be evaluated by analyzing the heartbeat parameters of HR and HRV.

HR analysis
Heart rate (HR) refers to the average number of heartbeats per minute, often used as a basic indicator of health monitoring. HR accuracy is defined as the percentage of time where HR detected by Doppler radar is within ± 2% the deviation of the reference HR detected by ECG, which is also used as an HR indicator in this paper.
To visually demonstrate the performance advantages of our proposed algorithm in different-period time windows, the HR estimation results of different algorithms at different time windows from dataset 6 are exhibited in Fig. 9 as a reference from ECG. As shown in Fig. 9a, the HR accuracies of the three algorithms within the deviation are about 96%. However, two severe deviation points in the HR curve of the FTPR remain to be discovered due to the dominant frequency divergence, and the losses of the heartbeat details are grave. As depicted in Fig. 9b, the HR accuracy of the FTPR algorithm is only 43.22%, resulting from the severe dominant frequency estimation deviation in short-period time windows. Thus, the TWV algorithm is introduced into the FTPR algorithm to avoid the problem, and the best time window is determined.
As a result, the HR accuracy of the FTPR-TWV is up to 94.07%. Furthermore, the high agreement can be accomplished for the FTPR-TWV algorithm in tracking the large HR mutation points and presenting the HR variety, supporting the accuracy and timeliness of our proposed algorithm.
In addition, an error evaluation index is utilized: The root-mean-squared error (RMSE) is calculated as Taking the estimated HR results in Table 2 into consideration, we can further assess the achievements of different algorithms. Table 2 is divided into two sub-tables of basic time window T = 10 s and T = 3 s. Each sub-table has 6 groups of experimental data, each group of data corresponds to the estimation results of three different algorithms, and each row corresponds to the results of different indicators of the algorithm. As shown in Table 2(a), the HR accuracies of the three algorithms are all over 92.79%, especially the average HR accuracy of FTPR-TWV is 99.70%, which is in high agreement with the reference value. As illustrated in Table 2(b), the average HR accuracies of FFT-TWV, FTPR, and FTPR-TWV are 82.20%, 40.68%, and 92.09%, respectively. The HR accuracy of FTPR-TWV is significantly higher than the comparison algorithm. In addition, the average RMSE of the FTPR-TWV is 0.90 BPM, which is less than that of the FFT-TWV, further verifying that more details of the shortperiod HR variety can be presented (Table 3). Furthermore, the relative error mean error is used as an indicator to compare the algorithm performances in recent years.

HRV analysis
Heart rate variability (HRV) refers to the time interval between consecutive heartbeats, revealing the regularity of each heartbeat. HRV analysis is of great significance in the diagnosis and medicine of heart diseases. Different from determining the BBIs by detecting the time interval of adjacent heartbeats in the time domain, the HR detection method by searching for the interested frequency in the frequency domain is applied in this paper and [31], where the BBIs can be calculated by the HR estimation value BPM, in units of beats per minute.

The mean relative error (MRE) is defined as
The evaluation indicators of HRV analysis, including the standard deviation of normalto-normal intervals (SDNN) and the root-mean-square successive difference of intervals (RMSSD), are defined as where N bbi is the number of successive heartbeat intervals and represents the time interval between successive heartbeats. In addition, the Bland-Altman analysis [42] is often used in medical clinical research to assess the consistency of different parameters. This method is also adopted to evaluate the agreement of the BBIs from radar and ECG in this paper, where the Bias denotes the average difference, and SD denotes the standard deviation of the difference. The lower and upper 95% limits of agreement (95% LoA) are defined as The results of HRV analysis are visualized in Table 4. As shown in

Discussion
As a result, accurate heartbeat estimation results have been implemented in short-period time windows of 3 s by using our proposed heartbeat detection algorithm. This method solves serious dominant frequency estimation deviation and insufficient frequency spectrum resolution in short-period time window. However, there are still many limitations as well as research topics to be further explored, mainly summarized as follows 1. The experimental test scenario in this paper is that the subject's chest is located directly in front of the radar for monitoring, which is relatively ideal and has limitations for wide applications in different fields. In practical applications, the human body may not be located at a certain angle or even back to the radar, resulting in a lower energy of the measured vital sign. Therefore, enhancement algorithm for vital sign components should be further studied. 2. Research on the methods of the heartbeat signal extraction in the presence of large body movements. The chest-wall motion caused by the heartbeat and respiration is much smaller than the large body movements, so it is difficult to extract vital signs from the radar echo signal in the presence of large body movements. As for heart- beat detection in the presence of large body movements, more advanced algorithms will be required to exclude abnormal parameter values in the radar echo signal. 3. Research on the time domain characteristics of vital signs. The frequency domain analysis of vital signs can only obtain information such as heart rate and respiration rate in a specific time window. Nevertheless, as for clinical medical diagnosis, the time domain information of vital signs at each moment should be observed, which is a major focus of future research work.