# A Method of UWB Radar Vital Detection Based onP Time Extraction of Strong Vital Signs

Authors: Zhen Yang; Jiming Cheng; Qingjie Qi; Xin Li; Yuning Wang
Journal: Journal of Sensors
 (2021)

---

## Abstract

The vital sign information in the echo signal of the UWB radar is weak, because of the interference of complex noise. In this paper, a method namedP times extraction of strong vital signs for processing echo signals of UWB radars is proposed. Different noises can be distinguished by the cumulative probability distribution of the echo signal and using different methods for processing according to corresponding characteristics. The vital sign information which most clearly represents the trapped person is selected using P times extraction of strong vital signs; then, the respiration and heartbeat rates are extracted. At 5 different distances, multiple sets of tests were carried out on static trapped persons and micromovement trapped persons and using a computer to extract vital signs from the obtained data. Experimental data shows that the algorithm proposed in this paper can extract the respiration and heartbeat rates of trapped persons, with small relative errors and variances, and has a certain reference value for UWB radar signal processing.


---

## Body

## 1. Introduction

The vital detection using a radar has broad application prospects in both military and civilian use, such as search and rescue of the wounded on the battlefield after the war, the search for survivors in the rubble after an earthquake or landslide, and patient monitoring during medical procedures [1–3]. Currently, two main types of radars are used for noncontact vital signal detection: continuous wave radar and UWB radar [4]. Continuous wave radar has been used for this purpose [5, 6], which has the advantages of simple structure, strong frequency domain recognition ability, and simple signal processing [7]. But it can only obtain the Doppler frequency information of the target, not distance information. Compared with the Doppler radar, the UWB radar can not only obtain the distance and frequency information of the trapped person at the same time but also has a great improvement in penetration ability, anti-interference ability, resolution, and power consumption [8, 9]. These advantages ensure the accuracy and sensitivity of the radar echo signal. Therefore, the UWB radar has gradually become the main application technology in the field of vital detection and monitoring and its echo signal preprocessing methods and vital sign extraction methods have gradually become research hotspots.When performing vital sign detection, the signal-to-noise ratio (SNR) of the echo signal is low [10–12], because the signal received by the radar contains not only the target reflection echo but also a large amount of complex noise. References [13, 14] used their own algorithms to process the echo signals and achievements have been made in improving the SNR. Reference [15] used the vital sign monitoring method of Nth peak capture to filter out clutter and noise, extract the respiration rate, and suppress its higher harmonics. Reference [16] proposed an UWB radar signal processing method with periodic sampling and superposition, which transfers the frequency information to the baseband and avoids the information loss caused by the dispersion of the spectrum. Reference [17] used singular value decomposition (SVD) to eliminate the noise in the vital signal and used FFT and Hilbert-Huang transform to extract the vital signs. Reference [18] used empirical mode decomposition (EMD) to decompose the radar echo signal, selecting intrinsic mode function (IMF) for reconstruction to obtain the respiration and heartbeat rates. References [19–21] made improvements on the basis of the EMD and achieved certain results of processing the UWB radar echo signals.In the algorithm proposed above, the basis of vital signs detected is one or more slow-time slices containing vital signs. However, in the case of complex noise interference, the vital signs contained in a single slow-time slice are weak, which is vulnerable to noise interference with poor stability and accuracy. What is worse is that external reasons cause the trapped person to move slightly, so the vital information contained in the selected slow-time slice is incomplete, and it is impossible to accurately obtain the trapped person’s vital information. If we want to obtain complete information, the combined superposition method of multiple slices will be used. However, this method will be disturbed by other noise slices and weaken the vital signs of the trapped person and the quality of the signs obtained cannot be guaranteed.In order to solve the above problems, an UWB radar vital information detection algorithm based onP times strong vital sign extraction is proposed. This method extracts P times vital signs from the regional echo signal that contains vital information, which obtains the most representative P group information of the trapped person. So, the method can improve the quality of the vital signs and ensure its accuracy and stability. Then, P group information of the trapped people will be parsed in the spectrum and the respiration and heartbeat rates were extracted according to the peak distribution of P group results. The accuracy and stability of the results were further guaranteed by P times different analysis.


## 2. Basic Principles of UWB Radar Vital Detection

The UWB radar vital detection signal is a micromotion signal of the trapped person. When the trapped person is alive, respiration and heartbeat will drive the chest to produce two different regular movements. It can be seen in Figure1, when the radar system and the trapped person are both fixed, the distance d changes with the movement of the thoracic cavity Δd. By continuously measuring the regular of d, the vital sign information of the trapped person can be obtained.Figure 1
Schematic diagram of UWB radar vital detection. Whered0 is the reference distance between the radar system and the trapped person, d is the actual distance, and Δd is the movement distance of the chest cavity.Research shows [22], the frequency band of respiration is between 0.2 Hz–0.5 Hz, the amplitude of chest cavity movement is 5 mm–15 mm, moving area is about 50 cm2; heartbeat frequency band is between 0.8 Hz–2.5 Hz, the amplitude change of chest cavity is about 2 mm-3 mm, moving area is about 10 cm2. With a complex application background, the heart rate is more difficult to extract than the respiration rate.In the process of respiration and heartbeat, the contraction and relaxation of muscles cause the cyclical fluctuations of the chest cavity, which are approximately sinusoidal. According to the basic principles of the UWB radar, the instantaneous distance between the radar antenna and the trapped person chest is [23]
(1)dt=d0+Δd=d0+Arsin2πfrt+Ahsin2πfht,where Ar and Ah are the amplitude changes of the thoracic cavity caused by respiration and heartbeat, respectively; fr is the trapped person respiration rate, fh is the heartbeat rate, and t is the slow time. According to equation (1), the time delay between the transmitted pulse and the received pulse can be obtained as
(2)τvt=2×dtv=τ0+τrsin2πfrt+τhsin2πfht,where v is the speed of light in air, τr=2Ar/v, and τh=2Ah/v. With the assumption that the environment is static and the thoracic cavity of the trapped person moves periodically, this movement is manifested in a time-varying channel impulse response:
(3)hτ,t=avδτ−τvt+∑iaiδτ−τi,where av and ai are the reflection coefficients of the trapped person and the ith environmental factor to radar waves; τv and τi are the corresponding time delays of the trapped person and the ith environmental factor. The signal measured by the receive antenna can be written as the convolution of the transmit pulse and the channel impulse response:
(4)rτ,t=sτ∗hτ,t=avsδτ−τvt+∑iaisδτ−τi+nτ,t,where sτ is the first-order Gaussian pulse signal as the propagation signal, ∗ is the convolution operation, and nτ,t is the random noise. Then, discrete the received signals and store it in the matrix Rm,n:
(5)Rm,n=rmδT,nTs=avsmδT−τvnTs+∑iaisτ−τi+nmδT,nTs,where t=nTsn=1,2,⋯,N is the slow time, Ts is the sampling interval in slow time, τ=mδT,m=1,2,⋯,M is the fast time, and δT is the sampling interval in fast time.


## 3. Received Signal Preprocessing

### 3.1. Noise Analysis

In order to facilitate subsequent signal processing and improve the SNR of the signal, the received signal needs to be preprocessed before the vital signal extraction. In addition to the trapped person vital information, the received signal of the UWB radar system also includes numerous noise interferences, including static clutter, linear trend, and random noise. In order to facilitate the understanding of the distribution of various noises, the analysis based on the cumulative probability distribution of slow time slices is carried out. Choosing a slow time slice arbitrarily, then, valuate the probability ofX=xii=1,2,⋯,I as pxi, where xi is the value appeared in the slice. The condition of amplitude distribution of all slices satisfies equation (6) as follows:
(6)∑i=1Ipxi=1,0≤pxi≤1.For anyX in this slice, the corresponding cumulative probability distribution function is
(7)FX=PX≤xi,where 0≤FX≤1 and FX is a nondecreasing function of X.Select the most stable static clutter slow time slice of the received signal, calculate its cumulative probability distribution, and draw the schematic diagram as shown in Figure2.Figure 2
Schematic diagram of the most stable static clutter slow time slice. (a) Time domain of slow time. (b) Cumulative probability distribution chart. The red straight line in (b) is the reference normal distribution curve of all amplitudes in the slow time slice.
(a)(b)Combined with equation (5), the environment of the trapped person does not change and the law of reflection and refraction of the radar signal does not change. It can be obtained that the static clutter does not change with time, which can be considered as a fixed value. Therefore, as shown in Figure 2(b), the static clutter reference value is around 55.3 and the slope of this position is larger, which means that the data is relatively concentrated. Due to the linear trend and random noise interference, the cumulative probability distribution curve of area B is consistent with the reference normal distribution curve. The prominent peak in Figure 2(a) corresponds to areas A and C in Figure 2(b), which are linear errors caused by radar instability. The error will affect the entire fast time direction, the amplitude is higher, but the frequency is lower. Random noise is mainly caused by electromagnetic interference, which is characterized by wide distribution and high frequency.And the schematic diagram of the slow time slice containing the vital sign of the trapped person is shown in Figure3.Figure 3
Schematic diagram of the slow time slice containing the vital sign of the trapped person. (a) Time domain of slow time. (b) Cumulative probability distribution chart.
(a)(b)The slopes of the curves inA and C regions in Figure 3(b) are larger, because A and C regions include not only linear errors but also the peak and trough values of the curve shown in Figure 3(a).


### 3.2. Preprocessing Methods

In view of the above noise characteristics, the following preprocessing methods are adopted.


#### 3.2.1. Linear Error Correction

The variance of the slow-time slice can reflect the degree of deviation between the slice and its mean, and the variance of the ideal static clutter slow time slice is 0. Calculate and set the slow time slice with the minimum variance as the radar stability reference slice. By calculating the occurrence probability of the amplitude in the reference slice, it can be considered that the amplitude with the probability of occurrence less than 0.5% is caused by radar instability. Then, the fast time slice which contains unstable values is corrected according to the trend of the reference slow time slice.


#### 3.2.2. Linear Trend Suppression

After completing the correction of the linear error, the slow time slice still has a linear trend and this problem will be solved by linear trend suppression in the slow time direction. After estimating the linear trend of the slow time slice, it is subtracted from the original slice to obtain a matrixW which does not contain static clutter and linear trend [24].
(8)WT=ΩT−zzTz−1zTΩT,where z=z1,z2,z1=1,2,⋯,NT,z2=1,1,⋯,1T.


#### 3.2.3. Signal Enhancement

In practical applications, the signal received by the UWB radar system is weak due to electromagnetic interference and multipath effects. Therefore, the received signal needs to be enhanced to improve the SNR. For preserving the vital signal, the Butterworth filtering is adopted to filter out low-frequency and high-frequency interference and its bandpass frequency is 0.1 Hz–3 Hz. For the matrix obtained after the above processing, it mainly contains vital information and zero mean noise. The vital information of the trapped person has a certain periodicity, but zero-mean noise does not which can be reduced to 0 through autocorrelation [10].
(9)Rn′=Exmn1xmn2,where 0≤n1,n2≤N−1,n’=n1–n2. E is for the average.The received signals before and after preprocessing are shown in Figure4 to verify the effectiveness of the preprocessing methods.Figure 4
Preprocessing effect diagram. (a) Grayscale of the original received signal. (b) Grayscale of the received signal after preprocessing.
(a)(b)The schematic diagram of the static clutter slow time slice after preprocessing is shown in Figure5, and the slow time slice containing the vital sign of the trapped person after preprocessing is shown in Figure 6.Figure 5
The schematic diagram of the static clutter slow time slice after preprocessing. (a) Time domain of slow time. (b) Cumulative probability distribution chart.
(a)(b)Figure 6
The schematic diagram of the slow time slice containing the vital sign of the trapped person after preprocessing. (a) Time domain of slow time. (b) Cumulative probability distribution chart.
(a)(b)The vital sign information in the original signal has been drowned in noise (see Figure4(a)). In Figure 4(b), various clutters are filtered effectively by preprocessing and the vital information of the trapped person is displayed. It can be seen in Figure 5 that in the preprocessed slow time slice, linear error, static clutter, and linear trend are all filtered out effectively and the remaining part is low-amplitude random noise with a reference value of 0. In Figure 6, with the vital signs being intact, the noise interference is effectively suppressed and the SNR is significantly improved.


## 3.1. Noise Analysis

In order to facilitate subsequent signal processing and improve the SNR of the signal, the received signal needs to be preprocessed before the vital signal extraction. In addition to the trapped person vital information, the received signal of the UWB radar system also includes numerous noise interferences, including static clutter, linear trend, and random noise. In order to facilitate the understanding of the distribution of various noises, the analysis based on the cumulative probability distribution of slow time slices is carried out. Choosing a slow time slice arbitrarily, then, valuate the probability ofX=xii=1,2,⋯,I as pxi, where xi is the value appeared in the slice. The condition of amplitude distribution of all slices satisfies equation (6) as follows:
(6)∑i=1Ipxi=1,0≤pxi≤1.For anyX in this slice, the corresponding cumulative probability distribution function is
(7)FX=PX≤xi,where 0≤FX≤1 and FX is a nondecreasing function of X.Select the most stable static clutter slow time slice of the received signal, calculate its cumulative probability distribution, and draw the schematic diagram as shown in Figure2.Figure 2
Schematic diagram of the most stable static clutter slow time slice. (a) Time domain of slow time. (b) Cumulative probability distribution chart. The red straight line in (b) is the reference normal distribution curve of all amplitudes in the slow time slice.
(a)(b)Combined with equation (5), the environment of the trapped person does not change and the law of reflection and refraction of the radar signal does not change. It can be obtained that the static clutter does not change with time, which can be considered as a fixed value. Therefore, as shown in Figure 2(b), the static clutter reference value is around 55.3 and the slope of this position is larger, which means that the data is relatively concentrated. Due to the linear trend and random noise interference, the cumulative probability distribution curve of area B is consistent with the reference normal distribution curve. The prominent peak in Figure 2(a) corresponds to areas A and C in Figure 2(b), which are linear errors caused by radar instability. The error will affect the entire fast time direction, the amplitude is higher, but the frequency is lower. Random noise is mainly caused by electromagnetic interference, which is characterized by wide distribution and high frequency.And the schematic diagram of the slow time slice containing the vital sign of the trapped person is shown in Figure3.Figure 3
Schematic diagram of the slow time slice containing the vital sign of the trapped person. (a) Time domain of slow time. (b) Cumulative probability distribution chart.
(a)(b)The slopes of the curves inA and C regions in Figure 3(b) are larger, because A and C regions include not only linear errors but also the peak and trough values of the curve shown in Figure 3(a).


## 3.2. Preprocessing Methods

In view of the above noise characteristics, the following preprocessing methods are adopted.


### 3.2.1. Linear Error Correction

The variance of the slow-time slice can reflect the degree of deviation between the slice and its mean, and the variance of the ideal static clutter slow time slice is 0. Calculate and set the slow time slice with the minimum variance as the radar stability reference slice. By calculating the occurrence probability of the amplitude in the reference slice, it can be considered that the amplitude with the probability of occurrence less than 0.5% is caused by radar instability. Then, the fast time slice which contains unstable values is corrected according to the trend of the reference slow time slice.


### 3.2.2. Linear Trend Suppression

After completing the correction of the linear error, the slow time slice still has a linear trend and this problem will be solved by linear trend suppression in the slow time direction. After estimating the linear trend of the slow time slice, it is subtracted from the original slice to obtain a matrixW which does not contain static clutter and linear trend [24].
(8)WT=ΩT−zzTz−1zTΩT,where z=z1,z2,z1=1,2,⋯,NT,z2=1,1,⋯,1T.


### 3.2.3. Signal Enhancement

In practical applications, the signal received by the UWB radar system is weak due to electromagnetic interference and multipath effects. Therefore, the received signal needs to be enhanced to improve the SNR. For preserving the vital signal, the Butterworth filtering is adopted to filter out low-frequency and high-frequency interference and its bandpass frequency is 0.1 Hz–3 Hz. For the matrix obtained after the above processing, it mainly contains vital information and zero mean noise. The vital information of the trapped person has a certain periodicity, but zero-mean noise does not which can be reduced to 0 through autocorrelation [10].
(9)Rn′=Exmn1xmn2,where 0≤n1,n2≤N−1,n’=n1–n2. E is for the average.The received signals before and after preprocessing are shown in Figure4 to verify the effectiveness of the preprocessing methods.Figure 4
Preprocessing effect diagram. (a) Grayscale of the original received signal. (b) Grayscale of the received signal after preprocessing.
(a)(b)The schematic diagram of the static clutter slow time slice after preprocessing is shown in Figure5, and the slow time slice containing the vital sign of the trapped person after preprocessing is shown in Figure 6.Figure 5
The schematic diagram of the static clutter slow time slice after preprocessing. (a) Time domain of slow time. (b) Cumulative probability distribution chart.
(a)(b)Figure 6
The schematic diagram of the slow time slice containing the vital sign of the trapped person after preprocessing. (a) Time domain of slow time. (b) Cumulative probability distribution chart.
(a)(b)The vital sign information in the original signal has been drowned in noise (see Figure4(a)). In Figure 4(b), various clutters are filtered effectively by preprocessing and the vital information of the trapped person is displayed. It can be seen in Figure 5 that in the preprocessed slow time slice, linear error, static clutter, and linear trend are all filtered out effectively and the remaining part is low-amplitude random noise with a reference value of 0. In Figure 6, with the vital signs being intact, the noise interference is effectively suppressed and the SNR is significantly improved.


## 3.2.1. Linear Error Correction

The variance of the slow-time slice can reflect the degree of deviation between the slice and its mean, and the variance of the ideal static clutter slow time slice is 0. Calculate and set the slow time slice with the minimum variance as the radar stability reference slice. By calculating the occurrence probability of the amplitude in the reference slice, it can be considered that the amplitude with the probability of occurrence less than 0.5% is caused by radar instability. Then, the fast time slice which contains unstable values is corrected according to the trend of the reference slow time slice.


## 3.2.2. Linear Trend Suppression

After completing the correction of the linear error, the slow time slice still has a linear trend and this problem will be solved by linear trend suppression in the slow time direction. After estimating the linear trend of the slow time slice, it is subtracted from the original slice to obtain a matrixW which does not contain static clutter and linear trend [24].
(8)WT=ΩT−zzTz−1zTΩT,where z=z1,z2,z1=1,2,⋯,NT,z2=1,1,⋯,1T.


## 3.2.3. Signal Enhancement

In practical applications, the signal received by the UWB radar system is weak due to electromagnetic interference and multipath effects. Therefore, the received signal needs to be enhanced to improve the SNR. For preserving the vital signal, the Butterworth filtering is adopted to filter out low-frequency and high-frequency interference and its bandpass frequency is 0.1 Hz–3 Hz. For the matrix obtained after the above processing, it mainly contains vital information and zero mean noise. The vital information of the trapped person has a certain periodicity, but zero-mean noise does not which can be reduced to 0 through autocorrelation [10].
(9)Rn′=Exmn1xmn2,where 0≤n1,n2≤N−1,n’=n1–n2. E is for the average.The received signals before and after preprocessing are shown in Figure4 to verify the effectiveness of the preprocessing methods.Figure 4
Preprocessing effect diagram. (a) Grayscale of the original received signal. (b) Grayscale of the received signal after preprocessing.
(a)(b)The schematic diagram of the static clutter slow time slice after preprocessing is shown in Figure5, and the slow time slice containing the vital sign of the trapped person after preprocessing is shown in Figure 6.Figure 5
The schematic diagram of the static clutter slow time slice after preprocessing. (a) Time domain of slow time. (b) Cumulative probability distribution chart.
(a)(b)Figure 6
The schematic diagram of the slow time slice containing the vital sign of the trapped person after preprocessing. (a) Time domain of slow time. (b) Cumulative probability distribution chart.
(a)(b)The vital sign information in the original signal has been drowned in noise (see Figure4(a)). In Figure 4(b), various clutters are filtered effectively by preprocessing and the vital information of the trapped person is displayed. It can be seen in Figure 5 that in the preprocessed slow time slice, linear error, static clutter, and linear trend are all filtered out effectively and the remaining part is low-amplitude random noise with a reference value of 0. In Figure 6, with the vital signs being intact, the noise interference is effectively suppressed and the SNR is significantly improved.


## 4.P Times Extraction of Strong Vital Signs

Vital sign information of the trapped person is contained in multiple slow time slices. Each slice contains the respiration and heartbeat rates of the trapped person, but there are certain differences in the strength of the vital information and whether it is complete. If you only select one slow time data slice as the benchmark for vital sign extraction, it is very easy to cause the wrong estimation of vital sign information about the trapped person. What is worse is that the trapped person shakes or moves in a small distance; only one slow time slice cannot fully reflect the vital signs of the trapped person. The method ofP times extraction of strong vital signs will be used to extract the respiration and heartbeat rates of the trapped person.The energy of the slow time slice is used as the basis for the judgment of the initial search starting point. The signal energy is defined as(10)en=∑m=1Mfm,n2,where fm,n is the amplitude of the mth row and nth column in the matrix and en is the energy to the nth slow time slice. The vital sign information of the trapped person is generally contained in 4–6 slow time slices, and an energy window with a width of 3 is used to prevent the influence of individual slow time slices on the positioning results.
(11)En=∑n−1n+1en,where En is the energy of the nth sampling window of the received signal, n=2,3,⋯,N−1 and the window with maximum energy will be used as the first sampling position.Generally, within 2 s, a person will produce respiration (exhaling or inhaling) movements and heartbeats multiple times. Therefore, choose 2 s as the base time length and step length for extraction ofP times. In the first 2 s, select P starting points equally in the reference window and search backward in the slow-time direction with 2 s as the step length. In every 2 s, select and save the strongest vital sign information, then record the position of this extraction, and use this position as the reference position for the next step, until all of the P times vital signal extractions are completed. The P results obtained contain the strongest vital sign information of multiple slow time slices at different moments.When extracting the respiration and heartbeat rates, perform FFT on theP times vital sign information extracted, find the peak value in the frequency range of respiration, and record the frequency and amplitude. Synthesize the P times results, the frequency with the most number of peaks is the respiration rate. If the number of peaks at multiple rates is the same maximum, the respiration rate will be selected which has the largest average amplitude. If the average amplitude is still the same, the frequency will be averaged. According to the above methods, the respiration rate of the trapped person can be determined; then, the one with the maximum average amplitude of all the alternative frequencies is selected as the final respiratory frequency. The vital sign information flow chart is shown in Figure 7, and the steps of the P times extraction of strong vital sign algorithm are shown in Algorithm 1.Figure 7
Flow chart of vital sign extraction, where step A is the preprocessing. Step B calculates energy from the slow time dimension and determines the maximum energy window position through step C. Step D is the process ofP times extraction of strong vital signs, and finally, obtain the corresponding frequency through step E.Algorithm 1: Steps of the algorithm for P times extraction of strong vital signs.
Name:P times extraction of strong vital signs1) According to the signal energy, determine the initial sampling slow-time slice position and search area;2) In the first 2 s of the area, selectP extraction starting points evenly;3) Start extraction, initializep=1;4) Select thepth starting point to start vital signs extraction;5) In every 2 s, select and save the strongest vital sign information, then record the extraction position, update the next reference position, until the completion of thepth extraction;6) Ifp=P, go to step 7), otherwise p=p+1, go to step 4);7) CompleteP times extraction of strong vital signs


## 5. Results and Discussion

### 5.1. Laboratory Equipment Test of the Stationary Trapped Person

The main experimental equipment is an UWB radar system based on NVA-6100 and a computer. The emission wave of this radar system is the first-order Gaussian pulse wave with a bandwidth of 0.85–9.55 GHz, which conforms to the standard of the UWB radar. The transmitting and receiving antennas of the system are both Vivaldi antennas, the gain reaches 6 dB, and the parallel distance between the two antennas is 10.2 cm. The received signal is connected to the computer through the USB; then, signal preprocessing and vital sign extraction are performed to obtain the corresponding respiration and heart rates.A large number of tests are carried out in MATLAB, differentP are performed on the received signal with a total length of 60 s, the result is obtained, when P>16, and the extracted respiration and heartbeat rate values tend to be stable. Therefore, P=16 is the best parameter of P times extraction of strong vital signs and the accurate vital information can be extracted in the shortest time. In the following experiment, P=16 is used as the extraction standard parameter.


### 5.2. Test of Stationary Trapped Person

In order to verify the effectiveness and accuracy of theP times extraction of strong vital signs to the stationary trapped person, set up the experimental scene as shown in Figure 8.Figure 8
Experimental scene of the stationary trapped person. The radar system is placed on a table 0.85 m from the ground. The volunteer instead of the trapped person.The volunteer was a 26-year-old male with a height of 180 cm and a weight of 80 kg. After multiple measurements, the respiration rate of the volunteer is 0.236 Hz and the heartbeat rate is 1.081 Hz. The volunteer sits in front of the radar system with the chest cavity facing the radar system, staying still and breathing evenly. According to the above method, the volunteer sits at 1 m, 2 m, 3 m, 4 m, and 5 m away from the radar system and collects 3 sets of received signals. The 16 times extraction of strong vital signs is used on the received signals and obtains respiration and heartbeat rates. The slow time slice containing the most vital sign information is used as a comparison, and the obtained respiration rate is shown in Table1, and the heartbeat rate is shown in Table 2.Table 1
Respiration rate of the stationary volunteer at different distances.
Distances (m)16 times (Hz)Slice (Hz)16 times (Hz)Slice (Hz)16 times (Hz)Slice (Hz)10.2260.2260.2440.2440.2120.21220.2150.2420.2030.2030.2070.20730.2390.1790.2350.2090.2570.19440.2680.2680.2770.2770.2430.30450.2620.2620.2260.2900.2350.277Table 2
Heartbeat rate of stationary volunteer at different distances.
Distances (m)16 times (Hz)Slice (Hz)16 times (Hz)Slice (Hz)16 times (Hz)Slice (Hz)11.1011.1011.0850.9951.1211.12121.1031.1571.0920.9911.0581.00731.0761.1651.0431.0431.0570.97141.1030.9831.0261.0261.1231.04251.0481.0481.1291.0321.1091.109Draw the data in Tables1 and 2 into the graphs shown in Figures 9 and 10, where 1–3 experiments are 3 sets of data obtained by volunteers at 1 m and 4–6 experiments are the data at 2 m and so on.Figure 9
The respiration rate of the stationary volunteer.Figure 10
The heartbeat rate of the stationary volunteer.The error at each distance, average error, relative, error and variance are calculated as shown in Table3.Table 3
Error comparison of the stationary volunteer.
ErrorRespiration rateHeartbeat rate16 timesSlice16 timesSliceError at 1 m (Hz)0.01400.01400.02130.0487Error at 2 m (Hz)0.02770.02270.01870.0800Error at 3 m (Hz)0.00830.04200.02230.0773Error at 4 m (Hz)0.02670.04700.03970.0640Error at 5 m (Hz)0.01230.04030.03630.0367Average error (Hz)0.01780.03320.02770.0613Relative error (%)7.5414.072.565.64Variance (×10−4)4.7314.039.7645.18Combining Figure9 and Table 3, it can be seen that for the respiration rate of the stationary volunteer at 1 m and 2 m, a single slice has a higher similarity to the volunteer’s vital sign information extracted by 16 times extraction of strong vital signs. In the range of 3–5 m with weak vital sign information, the extraction effect of 16 times extraction of strong vital signs is significantly better than that of single slice. For the heartbeat rate, although the single slice extraction result can reflect volunteer’s heart rate characteristics but the error is relatively larger and the error of this algorithm is smaller (see Figure 10 and Table 3).


### 5.3. Test of Micromotion Trapped Person

The test of the micromoving volunteer is carried out to verify the algorithm’s ability to extract signs of the micromoving trapped person. During the experiment, the volunteer moves his chest cavity slowly forward while maintaining uniform breathing and the overall moving distance was about 20 cm. According to the above measurement, the volunteer sits at 1 m, 2 m, 3 m, 4 m, and 5 m from the radar system, the respiratory rate obtained is shown in Table4, and the heartbeat rate is shown in Table 5.Table 4
Respiration rate of the micromoving volunteer at different distances.
Distances (m)16 times (Hz)Slice (Hz)16 times (Hz)Slice (Hz)16 times (Hz)Slice (Hz)1 m0.2710.2850.2380.2640.2900.2642 m0.2430.2910.2190.2400.2070.2073 m0.2270.2560.2850.2210.2440.2444 m0.2080.2860.2420.2910.2800.2025 m0.2600.1830.2600.2600.2960.227Table 5
Heartbeat rate of the micromoving volunteer at different distances.
Distances (m)16 times (Hz)Slice (Hz)16 times (Hz)Slice (Hz)16 times (Hz)Slice (Hz)1 m1.0570.9491.0311.0051.0310.8992 m1.0441.0441.0770.9621.0630.9743 m1.0580.9831.0810.8611.0970.9474 m1.0140.8231.0660.9690.9520.8135 m1.3951.1351.1460.9841.0970.947Plot the data in Tables4 and 5 as shown in Figures 11 and 12.Figure 11
The respiration rate of the micromoving volunteer.Figure 12
The heartbeat rate of the micromoving volunteer.Similarly, the error at each distance, average error, relative error, and variance are calculated as shown in Table6.Table 6
Error comparison of the micromoving volunteer.
ErrorRespiration rateHeartbeat rate16 timesSlice16 timesSliceError at 1 m (Hz)0.03030.03500.04130.1300Error at 2 m (Hz)0.01770.02930.01970.0877Error at 3 m (Hz)0.02200.01430.01300.1507Error at 4 m (Hz)0.02600.04630.07030.2127Error at 5 m (Hz)0.03600.02870.13170.0950Average error (Hz)0.02640.03070.05520.1352Relative error (%)11.1913.015.1112.51Variance (×10−4)10.2512.4688.34226.98Combining Figure11 and Table 6, it can be seen that compared with the result of the stationary volunteer, the accuracy of the algorithm for extracting the respiration rate of the micromoving volunteer has decreased but it is still better than the single slice.Combining Figure12 and Table 6, it can be seen that the respiration rate obtained by a single slice is basically below the reference value and the average error is relatively larger. The extraction result of this algorithm is basically stable near the reference value and only has a significant error at the position of 4–5 m.


### 5.4. Performance Comparison

As the results shown in Tables1–6, for stationary trapped persons, the relative error of the respiration rate extracted by this method is 7.54% and the relative error of heartbeat rate is 2.56%, which are better than the result of single slice. Because of the movement of the micromovement trapped person, the results of respiration and heartbeat rates have a larger error. However, compared with the extraction results of a single slice, the relative errors of respiratory and heartbeat rate extraction results of this algorithm are reduced by 1.82% and 7.3%, respectively. And the variance of this algorithm is smaller than that of single slice results in all comparisons, which can reflect that this algorithm has better stability.


## 5.1. Laboratory Equipment Test of the Stationary Trapped Person

The main experimental equipment is an UWB radar system based on NVA-6100 and a computer. The emission wave of this radar system is the first-order Gaussian pulse wave with a bandwidth of 0.85–9.55 GHz, which conforms to the standard of the UWB radar. The transmitting and receiving antennas of the system are both Vivaldi antennas, the gain reaches 6 dB, and the parallel distance between the two antennas is 10.2 cm. The received signal is connected to the computer through the USB; then, signal preprocessing and vital sign extraction are performed to obtain the corresponding respiration and heart rates.A large number of tests are carried out in MATLAB, differentP are performed on the received signal with a total length of 60 s, the result is obtained, when P>16, and the extracted respiration and heartbeat rate values tend to be stable. Therefore, P=16 is the best parameter of P times extraction of strong vital signs and the accurate vital information can be extracted in the shortest time. In the following experiment, P=16 is used as the extraction standard parameter.


## 5.2. Test of Stationary Trapped Person

In order to verify the effectiveness and accuracy of theP times extraction of strong vital signs to the stationary trapped person, set up the experimental scene as shown in Figure 8.Figure 8
Experimental scene of the stationary trapped person. The radar system is placed on a table 0.85 m from the ground. The volunteer instead of the trapped person.The volunteer was a 26-year-old male with a height of 180 cm and a weight of 80 kg. After multiple measurements, the respiration rate of the volunteer is 0.236 Hz and the heartbeat rate is 1.081 Hz. The volunteer sits in front of the radar system with the chest cavity facing the radar system, staying still and breathing evenly. According to the above method, the volunteer sits at 1 m, 2 m, 3 m, 4 m, and 5 m away from the radar system and collects 3 sets of received signals. The 16 times extraction of strong vital signs is used on the received signals and obtains respiration and heartbeat rates. The slow time slice containing the most vital sign information is used as a comparison, and the obtained respiration rate is shown in Table1, and the heartbeat rate is shown in Table 2.Table 1
Respiration rate of the stationary volunteer at different distances.
Distances (m)16 times (Hz)Slice (Hz)16 times (Hz)Slice (Hz)16 times (Hz)Slice (Hz)10.2260.2260.2440.2440.2120.21220.2150.2420.2030.2030.2070.20730.2390.1790.2350.2090.2570.19440.2680.2680.2770.2770.2430.30450.2620.2620.2260.2900.2350.277Table 2
Heartbeat rate of stationary volunteer at different distances.
Distances (m)16 times (Hz)Slice (Hz)16 times (Hz)Slice (Hz)16 times (Hz)Slice (Hz)11.1011.1011.0850.9951.1211.12121.1031.1571.0920.9911.0581.00731.0761.1651.0431.0431.0570.97141.1030.9831.0261.0261.1231.04251.0481.0481.1291.0321.1091.109Draw the data in Tables1 and 2 into the graphs shown in Figures 9 and 10, where 1–3 experiments are 3 sets of data obtained by volunteers at 1 m and 4–6 experiments are the data at 2 m and so on.Figure 9
The respiration rate of the stationary volunteer.Figure 10
The heartbeat rate of the stationary volunteer.The error at each distance, average error, relative, error and variance are calculated as shown in Table3.Table 3
Error comparison of the stationary volunteer.
ErrorRespiration rateHeartbeat rate16 timesSlice16 timesSliceError at 1 m (Hz)0.01400.01400.02130.0487Error at 2 m (Hz)0.02770.02270.01870.0800Error at 3 m (Hz)0.00830.04200.02230.0773Error at 4 m (Hz)0.02670.04700.03970.0640Error at 5 m (Hz)0.01230.04030.03630.0367Average error (Hz)0.01780.03320.02770.0613Relative error (%)7.5414.072.565.64Variance (×10−4)4.7314.039.7645.18Combining Figure9 and Table 3, it can be seen that for the respiration rate of the stationary volunteer at 1 m and 2 m, a single slice has a higher similarity to the volunteer’s vital sign information extracted by 16 times extraction of strong vital signs. In the range of 3–5 m with weak vital sign information, the extraction effect of 16 times extraction of strong vital signs is significantly better than that of single slice. For the heartbeat rate, although the single slice extraction result can reflect volunteer’s heart rate characteristics but the error is relatively larger and the error of this algorithm is smaller (see Figure 10 and Table 3).


## 5.3. Test of Micromotion Trapped Person

The test of the micromoving volunteer is carried out to verify the algorithm’s ability to extract signs of the micromoving trapped person. During the experiment, the volunteer moves his chest cavity slowly forward while maintaining uniform breathing and the overall moving distance was about 20 cm. According to the above measurement, the volunteer sits at 1 m, 2 m, 3 m, 4 m, and 5 m from the radar system, the respiratory rate obtained is shown in Table4, and the heartbeat rate is shown in Table 5.Table 4
Respiration rate of the micromoving volunteer at different distances.
Distances (m)16 times (Hz)Slice (Hz)16 times (Hz)Slice (Hz)16 times (Hz)Slice (Hz)1 m0.2710.2850.2380.2640.2900.2642 m0.2430.2910.2190.2400.2070.2073 m0.2270.2560.2850.2210.2440.2444 m0.2080.2860.2420.2910.2800.2025 m0.2600.1830.2600.2600.2960.227Table 5
Heartbeat rate of the micromoving volunteer at different distances.
Distances (m)16 times (Hz)Slice (Hz)16 times (Hz)Slice (Hz)16 times (Hz)Slice (Hz)1 m1.0570.9491.0311.0051.0310.8992 m1.0441.0441.0770.9621.0630.9743 m1.0580.9831.0810.8611.0970.9474 m1.0140.8231.0660.9690.9520.8135 m1.3951.1351.1460.9841.0970.947Plot the data in Tables4 and 5 as shown in Figures 11 and 12.Figure 11
The respiration rate of the micromoving volunteer.Figure 12
The heartbeat rate of the micromoving volunteer.Similarly, the error at each distance, average error, relative error, and variance are calculated as shown in Table6.Table 6
Error comparison of the micromoving volunteer.
ErrorRespiration rateHeartbeat rate16 timesSlice16 timesSliceError at 1 m (Hz)0.03030.03500.04130.1300Error at 2 m (Hz)0.01770.02930.01970.0877Error at 3 m (Hz)0.02200.01430.01300.1507Error at 4 m (Hz)0.02600.04630.07030.2127Error at 5 m (Hz)0.03600.02870.13170.0950Average error (Hz)0.02640.03070.05520.1352Relative error (%)11.1913.015.1112.51Variance (×10−4)10.2512.4688.34226.98Combining Figure11 and Table 6, it can be seen that compared with the result of the stationary volunteer, the accuracy of the algorithm for extracting the respiration rate of the micromoving volunteer has decreased but it is still better than the single slice.Combining Figure12 and Table 6, it can be seen that the respiration rate obtained by a single slice is basically below the reference value and the average error is relatively larger. The extraction result of this algorithm is basically stable near the reference value and only has a significant error at the position of 4–5 m.


## 5.4. Performance Comparison

As the results shown in Tables1–6, for stationary trapped persons, the relative error of the respiration rate extracted by this method is 7.54% and the relative error of heartbeat rate is 2.56%, which are better than the result of single slice. Because of the movement of the micromovement trapped person, the results of respiration and heartbeat rates have a larger error. However, compared with the extraction results of a single slice, the relative errors of respiratory and heartbeat rate extraction results of this algorithm are reduced by 1.82% and 7.3%, respectively. And the variance of this algorithm is smaller than that of single slice results in all comparisons, which can reflect that this algorithm has better stability.


## 6. Conclusions

In this paper, a new vital sign detection algorithm based on the UWB radar is proposed. This method parses different noise interferences by the cumulative probability density and uses linear error correction, linear trend suppression, and signal enhancement to realize the preprocessing. The region of the trapped person is judged by the energy of the processed signal, andP times of extraction is used to obtain the strongest vital sign. Then, the respiration and heartbeat rates were determined by the distribution of the peak from the P group frequency spectrum. Experimental data show that this method can more accurately determine the respiration and heartbeat rates of the stationary trapped person and has a better performance for the micromoving one.



---

*Source: 7294604-2021-09-20.xml*