Next Article in Journal
A Novel Adaptively-Robust Strategy Based on the Mahalanobis Distance for GPS/INS Integrated Navigation Systems
Previous Article in Journal
An Unsupervised Deep Hyperspectral Anomaly Detector
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Probabilities of False Alarm for Vital Sign Detection on the Basis of a Doppler Radar System

by
Nguyen Thi Phuoc Van
1,
Liqiong Tang
1,
Subhas Chandra Mukhopadhyay
2,*,
Duc Minh Nguyen
3 and
Faraz Hasan
1
1
School of Engineering and Advanced Technology, Massey University, Manawatu Private Bag 11 222, Palmerston North 4442, New Zealand
2
School of Engineering, Macquarie University, Sydney 2109, Australia
3
School of Electronics and Telecommunications, Hanoi University of Science and Technology, Hanoi 100000, Vietnam
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(3), 694; https://doi.org/10.3390/s18030694
Submission received: 17 January 2018 / Revised: 21 February 2018 / Accepted: 22 February 2018 / Published: 26 February 2018
(This article belongs to the Section Remote Sensors)

Abstract

:
Vital detection on the basis of Doppler radars has drawn a great deal of attention from researchers because of its high potential for applications in biomedicine, surveillance, and finding people alive under debris during natural hazards. In this research, the signal-to-noise ratio (SNR) of the remote vital-sign detection system is investigated. On the basis of different types of noise, such as phase noise, Gaussian noise, leakage noise between the transmitting and receiving antennae, and so on, the SNR of the system has first been examined. Then the research has focused on the investigation of the detection and false alarm probabilities of the system when the transmission link between the human and the radar sensor system took the Nakagami-m channel model. The analytical model for the false alarm and the detection probabilities of the system have been derived. The proposed theoretical models for the SNR and detection probability match with the simulation and measurement results. These theoretical models have the potential to be used as good references for the hardware development of the vital-sign detection radar sensor system.

1. Introduction

The non-contact monitoring of vital signals (heart beat and breathing rate) has attracted much attention recently. In this vital-sign monitoring system, the transceivers send a radio microwave signal toward the human position and collect the reflected signal from the person in order to obtain useful information such as vital signs and human location. There are many different approaches. The Doppler radar system has attracted a great deal of interest. The first Doppler radar system was developed for medical applications in the 1970s [1]. This radar system worked at a 10 GHz frequency, and it alarmed when apnoea occurred. This system did not give the exact respiration of the patient, but it did give an alarm when the patient stopped breathing for 30 s.
The application of microwave radar systems in health care systems has attracted a great deal of research attention [2,3,4,5,6,7,8]. In [3], the authors discussed the potential of using an ultra-wideband (UWB) technique for the radar sensor system in biological platform applications. Cianca et al. considered frequency modulation ultra-wideband (FM-UWB) in monitoring heart beats and breathing rates. In their work, the important features for medical applications, such as penetration capability, precision ranging, and electromagnetic radiation, were regarded carefully. The FM-UWB radar showed significant advantages in detecting vital signals [4]. Similarly, in 2011, Domenico Zito et al. [5] investigated UWB system-on-chip radar sensors in 90 nm CMOS technology. This system could detect the heart beat and respiration of adults and babies and allowed for the continuous observation of a baby’s breathing rate. This system could detect the breathing rate at a distance of 45 cm for the respiratory disorder diagnosis purpose and could sense a driver falling asleep. In [7], CMOS technology was also employed with the novel scheme design bondwire-interconnection between the antenna and the Doppler radar sensor system. This system could detect from a distance of up to 105 cm. A new patent of mm-wave radar sensors for detecting driver lethargy was investigated by Brankovic et al. in 2018 [8]. Their system could calculate the fatigue probability of a driver on the basis of the heart beat and respiratory rate. The probability was then compared with the threshold to give the driver an alarm signal. In [6], Stefano et al. reviewed a different type of micro-radar system for medical applications. In their research, the different types of radar sensor systems were examined carefully from continuous-wave (CW) to FM, UWB or hybrid radar systems.
Besides biomedical applications, a signal processing method for the vital-sign-detecting radar sensor system is critical [9,10,11]. In [9], Choi et al. showed that the multi-human detection algorithm could detect multiple people sequentially by detecting the first single cluster and repeating the process for the adjacent cluster. The signal-to-noise ratio (SNR) is one of the most important parameters of the radar system. The theoretical calculation of the SNR in the Doppler radar system should be considered and compared with the measured result to validate the theoretical considerations [10]. A denoising method for through-wall vital signs was proposed by Xiaolin et al. [11]. Several techniques were employed in this research, including short-time Fourier transform (STFT) and empirical mode decomposition (EEMD). This processing method can reduce clutter and is straightforward to implement.
Finding survivors buried under rubble is a crucial application of a Doppler radar sensor system [12,13]. In such an application, the microwave beam must penetrate deeply through the debris or rubble to find living victims. Therefore, L and S bands seem to be good choices. Kun-Mu Chen et al. [12] built systems and conducted experiments at two frequencies, 450 and 1150 MHz, and found that the 1150 MHz system operated better when the debris had metallic wires. In 2016, F. JalaliBidgoli et al. [13] proposed a similar system operating at a 1150 MHz frequency with the focus on minimizing clutter and noise [13]. This system is portable and can detect a victim below a depth of 1.5 m. Contrary to the earlier work, the transmission link between the radar system and the human body is a Friis channel model. This model implies a line-of-sight (LOS) connection between the radar sensor system and the subject. For the purpose of finding survivors, the system has to operate through different materials. Clearly, the LOS channel model is not adequate to describe the environment in which the radar system is likely to be used.
This paper presents research that considers the previous investigations and the drawbacks. In this research work, the SNR is carefully taken into account with respect to the Nakagami-m channel model and its special case (Rayleigh fading channel model). This model was chosen because it can more adequately describe realistic situations [14]. The Nakagami-m fading channel model can be used for a wide range of radio links, from non-fading channels to land-mobile or indoor-mobile prorogations. This channel can also serve as a LOS/non-line-of-sight (NLOS) channel model [15,16]. In comparison with previous studies, in terms of a channel radio link between the vital-sign-detecting radar system and the human, this research considers a more realistic and common channel model, the Nakagami-m, for the vital-sign-detecting radar sensor system. To the best of the authors’ knowledge, this is the first work to use the SNR at the base band under the Nakagami-m radio channel link between a living person and the radar sensor system. Moreover, the detection/false alarm probabilities of the vital-sign-detecting radar sensor system have been established, simulated and analyzed.

2. System Model

The working principle of the proposed vital-sign-detecting radar sensor system is presented using the block diagram shown in Figure 1. The transmitting antenna sends continuous waves to the human and waits for the reflected signals from the human body. At the human’s chest, the RF signal is modulated by the periodic displacement of the chest before being reflected back to the receiving antenna [17]. At the receiver, the signal is amplified by a low-noise amplifier (LNA) and is then sent to the phase detector. In this diagram, the in-phase (I) and quadrature (Q) phase detectors are used to reduce the null-point problem.
The transmitted signal from the vital-sign-detecting radar sensor system can be described by
Y T = 2 P T cos 2 π f t + ϕ t
where P T is the power of the transmitter, f is the carrier frequency, t is the transpiring time, and ϕ t is the phase noise. The initial phase of the transmitting signal is assumed to be zero. The signal Y T is then reflected by the human’s chest. The movement of the chest can be considered as the periodic signal x t . The distance from the human’s location to the antennae of the radar sensor system is d (m). The received signal Y R can be expressed as follows [18]:
Y R = A R cos 2 π f t 4 π d λ 4 π x t λ + ϕ t 2 d c
where c is the light velocity, λ is the wavelength of the signal, and A R is the return signal amplitude that depends on te transmitted power, the gain of the antennae, the radar cross-section, and reflection loss from the human body and environment (channel links from the radar to the human body and from the human body to the radar).
In fact, in practice, the transmitting antenna is quite close to the receiving antenna, and thus the leakage noise is unavoidable. Moreover, the radar system also receives the reflected signal from the stationary part of the human and other clusters [17]. Therefore, the stable reflected signal should be considered at the receiving antenna. The receiving signal can be rewritten as
Y R = A R cos 2 π f t 4 π d λ 4 π x t λ + ϕ t 2 d c + A l e a k cos 2 π f t t l + ϕ t t l + A c cos 2 π f t 4 π d 0 λ + ϕ t 2 d c
where A l e a k is the amplitude of the leakage signal, t l is the time delay between the transmitting leakage signal and the local oscillator (LO) signal, and A c is the amplitude of the stationary cluster signal.
In previous studies, channel links between the radar sensor system and the human were considered as in free space [10]. This model implied that there were LOS links between the radar system and the human. However, the vital-sign-detecting radar system has to operate through the NLOS channel. Because of this reason, in this study, the Nakagami-m channel model was chosen as the channel model between the radar sensor system and the human. This channel model was considered to be close to the practical situation and could model both LOS and NLOS environments [14]. On the basis of the theory of radar systems [19], the amplitude of the receiving signal A R V o l t can be formulated as follows [20]:
A R = 2 P T G T G R σ a ρ r ρ t r λ 2 | g a | 2 | g b | 2 4 π 3 d 2 α
where G T = 10 G T d B i 10 and G R = 10 G R d B i 10 are the gains of the transmitting antenna and receiving antenna, respectively; λ (m) is the wavelength of the signal; ρ r and ρ t are the radiated efficiencies of the receiving and transmitting antennae; r is the reflection coefficient of the human body; σ a (m2) is the radar cross-section (RCS) of the human’s chest; and α is the path loss. The channel coefficients for the links from the radar to the human body and from the human body to the radar are g a and g b , respectively.
The leakage amplitude can be estimated as follows [21]:
A l e a k = 2 η l e a k P T
and
A c = 2 P T G T G R σ c ρ r ρ t r λ 2 | g a | 2 | g b | 2 4 π 3 d 2 α
where G r x is the gain of the receiver, and σ c is the RCS of other parts of the human body.
The receiving signal Y R is then mixed with the LO signal to obtain the intermediate frequency (IF) signal. As can be seen from Figure 1, I & Q demodulation is applied to reduce the null-point problem [22]. The output signals are then fed to the bandpass filter (BPF). The low- and high-end frequencies of the BPF are expressed as f L and f H , respectively.
The base band can be written as follows [21]:
B ( t ) = A R G R X cos θ 0 + 4 π x t λ + ϕ s t + A c G R X cos θ 0 + ϕ c t + A l e a k G R X cos θ l + ϕ l t + n ( t ) + n 1 / f ( t )
where θ 0 is a constant phase shift due to the distance from the human body to the radar system; θ l is a constant phase shift that depends on the delay time t l ; ϕ s t , ϕ c t , and ϕ l t are the residual phase noise of the reflected vital signals (respiratory signal and heart signal), clutter signal, and antennae leakage signal, respectively; n ( t ) and n 1 / f ( ) are the additive white Gaussian noise (AWGN) and base-band 1 / f noise. The residual phase noise in Equation (7) can be calculated as follows [21]:
ϕ s t ϕ c t = ϕ t ϕ t 2 d c
ϕ l t = ϕ t ϕ t t l
In the following sections, the SNR is considered carefully, and the comparison between the theoretical calculation and the measurement is examined.

3. SNR Analysis

3.1. Signal Power

In the case in which the optimal demodulation is achieved, the receiving power at the output of the mixer can be given as follows [10,21,22]:
S B = A R 2 G R X 4 π x t λ 2 = 2 P T G T G R σ a ρ r ρ t r λ 2 | g a | 2 | g b | 2 4 π 3 d 2 α G R X 16 π 2 x 2 t ¯ λ 2 = P T G T G R σ a ρ r ρ t r | g a | 2 | g b | 2 G R X x 2 t ¯ 2 π d 2 α
In this work, the main sources of noise have been considered carefully: residual phase noise, AWGN, and the base-band 1 / f noise.

3.2. Residual Phase Noise

Clearly, phase noise degrades the performance of the Doppler radar system. The effect of this phenomenon is to spread energy of the clutter into the target frequency range. The spreading of clutter obscures the target signal and reduces the signal-to-clutter ratio [23]. In this section, the residual phase noise power is calculated. The residual phase noise consists of phase noise by the antennae leakage signal and the clutter. As defined in [10,21,23], the base-band noise spectral density ( S ϕ ( f 0 ) ) is derived from the RF phase noise spectral density S θ ( f 0 ) :
S ϕ ( f 0 ) = S ϕ ( f 0 ) 4 s i n 2 2 π d f 0 c S ϕ ( f 0 ) 16 π 2 d 2 f 0 2 c 2
S θ f 0 can be derived from the phase noise at 1 Hz, S θ 1 , if the close-in phase noise has a 30 dB/decade slope, as follows [21]:
S ϕ s f 0 S ϕ 1 1 H z 3 f 0 3
On the basis of Equations (11) and (12), S ϕ c ( f 0 ) can be given:
S ϕ s f 0 S ϕ 1 1 H z 3 f 0 3 16 π 2 d 2 f 0 2 c 2
The mean square of the residual phase noise in the time domain is calculated through the integral of the spectrum over the frequency range [ f L t o f H ] of the BPF [21]:
ϕ s t 2 ¯ = f L f H S ϕ s f 0 d f 0 = 16 π 2 ( 1 H z ) 3 S ϕ 1 d 2 c 2 l n f H f L
Finally, the residual phase noise power ( N ϕ s ) of the received signal can be calculated [21]:
N ϕ s A R 2 G R X ϕ s t 2 ¯ 2 P T G T G R σ a ρ r ρ t r λ 2 | g a | 2 | g b | 2 4 π 3 d 2 α G R X 16 π 2 ( 1 H z ) 3 S ϕ 1 d 2 c 2 l n f H f L P T G T G R σ a ρ r ρ t r λ 2 | g a | 2 | g b | 2 2 π d 2 α 1 f 2 G R X S ϕ 1 l n f H f L
By following similar steps to Equations (11)–(15), the residual phase noise power ( N ϕ c ) of the clutter signal can be calculated:
N ϕ c P T G T G R σ c ρ r ρ t r λ 2 | g a | 2 | g b | 2 2 π d 2 α 1 f 2 G R X S ϕ 1 l n f H f L
This then turns to antennae leakage phase noise; the base-band spectral density of antennae leakage phase noise S ϕ l f 0 is given as
S ϕ l f 0 S ϕ 1 1 H z 3 f 0 3 16 π 2 t l 2 f 0 2 4
The power of the antennae leakage phase noise ( N ϕ l ) can be expressed as
N ϕ l A l e a k 2 G R X ϕ l t 2 ¯ A l e a k 2 G R X f L f H S ϕ l f 0 d f 0 A l e a k 2 G R X 16 π 2 ( 1 H z ) 3 S ϕ 1 t l 2 2 2 l n f H f L 32 η l e a k P T G R X π 2 S ϕ 1 t l 2 2 2 l n f H f L

3.3. Additive White Gaussian Noise

AWGN is the main noise at the receiver’s input, and afterwards, this noise is converted to the baseband with the power ( N W ) as follows [10,21]:
N W = 8 G R X k T B N F
where G R X = 10 G R X d B i 10 is the gain of the receiver, B (Hz) is the bandwidth, T (K) is the absolute temperature, k is Boltzman’s constant, and N F = 10 N F d B 10 is the noise figurer of the receiver.

3.4. 1 / f Noise

The mixer is chosen to minimize the 1 / f noise. The power of the 1 / f noise ( N 1 / f ) can be designated by the noise power ( P 1 / f ( 1 ) ) in a 1 Hz bandwidth centered at 1 Hz:
N 1 / f = f L f H P 1 / f 1 f 1 d f = P 1 / f 1 l n f H f L

3.5. SNR

The noise sources are uncorrelated, and the powers of the noise are added together. The SNR is written as
S N R = S B N ϕ s + N ϕ c + N ϕ l + N W + N 1 / f

4. Detection and False Alarm Probabilities

In this section, the analysis and simulation results of the detection probability ( P D ) and false alarm probability ( P f ) of the proposed system are presented. Estimations of the detection/false alarm probabilities are very important for the vital-sign-detecting radar system. The decision hypothesis is detected when the receiving signal and noise are greater than the threshold value P t h . The decision is a false alarm when the power of the noise is greater than the threshold P t h [9,19]. The false alarm/detection probabilities are estimated under two hypotheses, H 0 and H 1 . H 0 corresponds to the case that the vital signal is not presented, and H 1 is the opposite case.

4.1. Detection Probability

The detection probability can be defined as follows:
P D = P r S B + N ϕ s + N ϕ c + N ϕ l + N W + N 1 / f P t h H 1
By using C = P T G T G R σ a ρ r ρ t r G R X x 2 t ¯ 2 π d 2 α and D = P T G T G R ρ r ρ t r λ 2 2 π d 2 α 1 f 2 G R X S ϕ 1 , P D can be rewritten as follows:
P D = P r C | g a | 2 | g b | 2 + D σ a | g a | 2 | g b | 2 + D σ c | g a | 2 | g b | 2 + N ϕ l + N W + N 1 / f P t h = P r C + D σ a + D σ c | g a | 2 | g b | 2 + N ϕ l + N W + N 1 / f P t h = 1 P r | g a | 2 | g b | 2 P t h N ϕ l + N W + N 1 / f C + D σ a + D σ c
The channel model between the radar system and the human body is Nakagami-m, where m is the fading severity parameter. We assume that | g a 2 | and | g b 2 | are independent gamma random variables with the severity fading factor m a and m b ; in the case m = 1 , the Nakagami-m model becomes the Rayleigh fading channel model. By applying the work from [20], the detection probability can be given as
P D = 2 P t h N ϕ l + N W + N 1 / f C + D σ a + D σ c μ a μ b K 1 2 P t h N ϕ l + N W + N 1 / f C + D σ a + D σ c μ a μ b
where μ a and μ b are the mean values of | g a 2 | and | g b 2 | , respectively. K 1 ( . ) is the first-order modified Bessel function of the second kind, K 1 x = 0 e x c o s h t c o s h t d t [24].
When it comes to the normal case of the Nakagami-m channel model, the detection probability of the sensor system can be calculated on the basis of the distribution of the products of two independent variables | g a 2 | and | g b 2 | . Associated with the probability theory in [25], P D can be calculated as
P D = 1 u 2 1 m a m b D m 1 + m 2 1 , m 1 m 2 2 P t h N ϕ l + N W + N 1 / f C + D σ a + D σ c
where u = 2 τ m 1 1 τ m 2 1 and D μ , v ( x ) = 0 x x μ K v ( x ) d x ; τ . and K . are the gamma function and Bessel function, respectively. These functions can be found in [26] and calculated by available software such as Matlab or Wolfram. In the special case m a = m b = I + 1 (where I = 0, 1, 2, ..... ), P D is given as follows [25]:
P D = 2 2 i i ! 1 j = 0 i 2 P t h N ϕ l + N W + N 1 / f C + D σ a + D σ c 2 i + 1 j × K j + 1 2 P t h N ϕ l + N W + N 1 / f C + D σ a + D σ c 2 j / i j !

4.2. False Alarm Probability

Similarly, the false alarm probability is given as
P f = P r N ϕ s + N ϕ c + N ϕ l + N W + N 1 / f P t h H 0 = 1 P r | g a | 2 | g b | 2 P t h N ϕ l + N W + N 1 / f D σ a + D σ c
Under the Rayleigh fading channel model, the false alarm probability can be expressed as
P f = 2 P t h N ϕ l + N W + N 1 / f D σ a + D σ c μ a μ b K 1 2 P t h N ϕ l + N W + N 1 / f D σ a + D σ c μ a μ b
Under the Nakagami-m channel model ( m 1 / 2 ), the false alarm probability can be given as
P f = 1 u 2 1 m a m b D m 1 + m 2 1 , m 1 m 2 2 P t h N ϕ l + N W + N 1 / f D σ a + D σ c
In the case m a = m b = I + 1 (where I = 0, 1, 2, ..... ), P f can be calculated as follows [25]:
P f = 2 2 i i ! 1 j = 0 i 2 P t h N ϕ l + N W + N 1 / f D σ a + D σ c 2 i + 1 j × K j + 1 2 P t h N ϕ l + N W + N 1 / f D σ a + D σ c 2 j / i j !

5. Simulation Result

In this section, the SNR and detection/false alarm probabilities are discussed in terms of simulation and analysis. The values used for simulation and analysis are given in Table 1.

5.1. SNR and Detection/False Alarm Probabilities under Rayleigh Fading Channel

Firstly, a special case of the Nakagami-1 (Rayleigh fading channel) channel model has been considered. This model implied that there was no LOS transmission between the human and the radar system. The power of the receiving signal and different types of noise are shown in Figure 2. The leakage noise was relatively high, around 90 dBm, while the AWGN and 1 / f noise were lower than 120 dBm. The residual phase noise power of the received signal had less effect on the system than the residual phase noise power of the clutter signal; both decreased with the distance increase. It is noted from this simulation that with the transmitting power of 0 dBm, the receiving power was greater than the power of each type of noise when the distance was less than 3 m. This figure gives an important idea to enhance the system. The most dominant noise was leakage noise; this type of noise is defined by the isolation level between the transmitting and receiving antennae. Therefore, this noise could be minimized by hardware development and signal processing techniques.
Besides the power of the signal and noise, the SNR is the most important feature of the radar sensor system. The overall SNR of system is displayed in Figure 3. The first point to note is that the SNR decreased whenever the distance between the human body to the radar system and the antennae leakage increased. At the level of antennae leakage of 20 dB, the SNR reduced from around 30 dB at a distance of 1 m to 20 dB at a distance of 7 m. When the leakage between the two antennae was at 40 dB, the SNR at the same distance increased to 0 dB.
Figure 4 displays the detection probability of the system at three different levels of the threshold power of the receiver’s filter, 70 , 80 , and 90 dBm. The lowest threshold value gave the highest detection probability. When the distance increased, the detection probability decreased. The threshold power of the receiver’s filter defined the steepness of the detection probability line.
To evaluate the accuracy of the vital-sign-detecting radar sensor, the false alarm should be taken into account. Figure 5 shows that when the threshold power of the receiver’s filter was set at 70 dBm, the false alarm probability was high at distances of less than 50 cm. This phenomenon was due to the high level of phase noise when the radar was too close to the human body. When this threshold value reduced to 90 dBm, the false alarm moved to the zero level. The simulation and analysis results are presented in Figure 2, Figure 3, Figure 4 and Figure 5, Figure 6. The performance of the system is discussed under the Rayleigh Fading channel model. This model implied a NLOS link between the human body and the radar system.

5.2. Detection/False Alarm Probabilities of the System under the Nakagami-2 Channel Model

When the quality of the channel increased, for example, when the fading factor m = 2, in comparison with the Rayleigh fading link, the detection probability increased significantly. At a distance of around 7 m, the detection probability was approximately 80% as shown in Figure 6. The false alarm detection probability was relatively high at a close distance (less than 50 cm) because of the high level of phase noise.

6. Measurement Results

This section presents the measurement results of detecting the respiratory rate of humans. The SNR and the detection capability/probability were considered and compared to those of the theoretical models. The measurement was a setup based on the N 5244 A P N A X Microwave Network Analyzer, 43.5 GHz. Two ports of this device were connected to two antennae, as shown in Figure 7. Five people aged from 21 to 25 participated in this experiment. The objects (humans) were required to sit in front of the antennae; one antenna played the role of the transmitting antenna, and the other was the receiving antenna. The operating frequency of the system was set up in L band to obtain higher penetration through different environments [27]. In this study, the frequency was chosen as the optimal operating frequency ( 1.6 GHz) of the microstrip/path antenna. The antennae’s size was 10 cm × 10 cm with a gain of 5 dBi. For each distance (1, 2, 3, and 4 m), objects were measured over 5 min continuously in the normal electronic lab. The internal phase detector of the N 5244 A P N A X Microwave Network Analyzer was utilized to find the phase shift of the receiving signal. This signal was then processed in Matlab to find the breathing rate. At the same time, the object’s breathing rate was measured by the wearable sensor S h i m m e r 3 E C G / E M G with five probe leads connected, as shown in Figure 7; the output signals of the S h i m m e r 3 were sent to the computer by a Bluetooth connection. The measurement results are shown in Table 2. The breathing rates at different distances were close to the reference values. For the rescue detection purpose, the deviation around the beats per minute was acceptable. Regarding the SNR values, at every distance, the measured results were close to the simulation results. The difference was only around 1.5 dBm. The measurement result of the SNR of the second object was the closest to the simulation outcome and the mean values of the five cases. The detailed measured results of object 2 are presented in Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12.
Figure 8 shows the reference signal and the measured signal at a distance of 4 m in the time domain by the radar system. These signals were then processed by the fast Fourier transform (FFT) to find the breathing rate. The signals in the frequency domain can be seen in Figure 9. The breathing rate of the measured person was 16 beats per minute. The window size of the FFT was 30 s and moved forward every second; the number of samples in each second was 256. The SNR and the detection capability of the system were evaluated at each second after the first 30 s.
The SNR of the breathing signal was calculated by the mean value of the signal in the following range: (center frequency − 4 beats) to (center frequency + 4 beats), over the mean value of the remaining signal from 5 to 40 beats/min [28]. The breathing rate of the person was measured at four positions (1, 2, 3, and 4 m). The SNR of the system at each distance was calculated by the mean values of the SNR in each window. Figure 10 shows that the measured result was in line with the theoretical analysis in the case of the Nakagami-2 model. In this situation, the fading parameter was 2, corresponding the LOS connection between the radar system and the human body. Figure 11 presents the detection capability of the system in the time domain.
In Figure 11, the detection capabilities of the system at four distances are presented. When the detection capability equalled 1, the radar system could detect the breathing rate of a human. As can be seen in Figure 11, at most time points, the measured result matched with the simulation outcome. However, when the object was closer to the radar system, the detection capability was higher and the false detection in the measured result was distributed over a period of time. For example, at a distance of 3 m, the false detections occurred at around the 100 th and 150 th seconds, while the simulation showed that the false detection occurred randomly at different time points. This phenomenon could be explained by, for some duration in the laboratory environment, the noise level increasing significantly and negatively affecting the radar system; the simulation result, however, followed the statistical distribution—the errors occurred randomly at different time points.
The detection probability of the system is shown in Figure 12. At distances of 2 and 3 m, the measurement results were close to those of the theoretical model, while at a distance of 1 m, the detection probability was higher than that of the simulation result; the opposite result at a distance of 4 m was found. These differences may have appeared because the measuring time was not long enough. However, the theoretical detection probability model is a good reference to evaluate the performance of the radar sensor system in the long term.

7. Conclusions

In this paper, the SNR and the detection probability of the vital-sign-detecting radar sensor system have been considered and simulated. The measured result of the SNR was close to the theoretical calculation. This theoretical model is a basis for hardware development. The designer could optimize their design by choosing different parameters, such as leakage noise, phase noise, and so on, to obtain a high value of the SNR. In addition, on the basis of different types of noise, the detection probability has also been investigated in terms of simulation and analysis. The simulation followed the analytical models. When the system operated at close distances, the phase noise was relatively high and caused a false alarm problem. The threshold power of the receiver’s filter affected the alarm/detection probabilities. This value should be enhanced to improve the detection probability and reduce the false alarm probability. Moreover, the isolation between the transmitting and receiving antennae is an important factor to boost the SNR and detection probability of the system.

Acknowledgments

This research is funded by the NZ aid program, New Zealand Ministry of Foreign Affairs and Trade, and Faculty for the Future Program- Schlumberger Foundation. Authors would like to express our very great appreciation to Dr. Hung Tran Malardalen University, Vasteras, Sweden for his contribution on development of mathematical model in this paper

Author Contributions

Nguyen Thi Phuoc Van wrote the paper. Duc Minh Nguyen and Nguyen Thi Phuoc Van conceived, designed the experiments. Liqiong Tang, Subhas Chandra Mukhopadhyay and Duc Minh Nguyen analyzed data and revised the paper. Faraz Hasan checked the simulation result and revised the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Caro, C.; Bloice, J. Contactless apnoea detector based on radar. Lancet 1971, 298, 959–961. [Google Scholar] [CrossRef]
  2. McEwan, T.E. Body Monitoring and Imaging Apparatus and Method. US Patent 5,573,012, 12 November 1996. [Google Scholar]
  3. Bilich, C.G. Bio-medical sensing using ultra wideband communications and radar technology: A feasibility study. In Proceedings of the Pervasive Health Conference and Workshops, Innsbruck, Austria, 29 November–1 December 2006; pp. 1–9. [Google Scholar]
  4. Cianca, E.; Gupta, B. FM-UWB for communications and radar in medical applications. Wirel. Pers. Commun. 2009, 51, 793–809. [Google Scholar] [CrossRef]
  5. Zito, D.; Pepe, D.; Mincica, M.; Zito, F.; Tognetti, A.; Lanatà, A.; De Rossi, D. SoC CMOS UWB pulse radar sensor for contactless respiratory rate monitoring. IEEE Trans. Biomed. Circuits Syst. 2011, 5, 503–510. [Google Scholar] [CrossRef] [PubMed]
  6. Pisa, S.; Pittella, E.; Piuzzi, E. A survey of radar systems for medical applications. IEEE Aerosp. Electron. Syst. Mag. 2016, 31, 64–81. [Google Scholar] [CrossRef]
  7. Chan, C.H.; Chou, C.C.; Chuang, H.R. Integrated Packaging Design of Low-Cost Bondwire Interconnection for 60-GHz CMOS Vital-Signs Radar Sensor Chip With Millimeter-Wave Planar Antenna. IEEE Trans. Compon. Packag. Manuf. Technol. 2018, 8, 177–185. [Google Scholar] [CrossRef]
  8. Branković, V.; Jovanović, P.; Mihajlović, V.; Savić, M.; Tasovac, D. MM-Wave Radar Driver Fatigue Sensor Apparatus. US Patent 9,862,271, 9 January 2018. [Google Scholar]
  9. Choi, J.W.; Nam, S.S.; Cho, S.H. Multi-Human Detection Algorithm based on an Impulse Radio Ultra-Wideband Radar System. IEEE Access 2016, 4, 10300–10309. [Google Scholar] [CrossRef]
  10. Droitcour, A.D.; Boric-Lubecke, O.; Kovacs, G.T. Signal-to-noise ratio in Doppler radar system for heart and respiratory rate measurements. IEEE Trans. Microw. Theory Tech. 2009, 57, 2498–2507. [Google Scholar] [CrossRef]
  11. Liang, X.; Zhang, H.; Ye, S.; Fang, G.; Gulliver, T.A. Improved denoising method for through-wall vital sign detection using UWB impulse radar. Digit. Signal Process. 2018, 74, 72–93. [Google Scholar] [CrossRef]
  12. Chen, K.M.; Huang, Y.; Zhang, J.; Norman, A. Microwave life-detection systems for searching human subjects under earthquake rubble or behind barrier. IEEE Trans. Biomed. Eng. 2000, 47, 105–114. [Google Scholar] [CrossRef] [PubMed]
  13. JalaliBidgoli, F.; Moghadami, S.; Ardalan, S. A Compact Portable Microwave Life-Detection Device for Finding Survivors. IEEE Embed. Syst. Lett. 2016, 8, 10–13. [Google Scholar] [CrossRef]
  14. Mishra, M.K.; Sood, N.; Sharma, A.K. Efficient BER Analysis of OFDM System over Nakagami-m Fading Channel. Int. J. Adv. Sci. Technol. 2011, 37, 37–46. [Google Scholar]
  15. Alouini, M.S.; Simon, M.K. Performance of coherent receivers with hybrid SC/MRC over Nakagami-m fading channels. IEEE Trans. Veh. Technol. 1999, 48, 1155–1164. [Google Scholar] [CrossRef]
  16. Suraweera, H.A.; Smith, P.J.; Armstrong, J. Outage probability of cooperative relay networks in Nakagami-m fading channels. IEEE Commun. Lett. 2006, 10. [Google Scholar] [CrossRef]
  17. Kuo, H.C.; Lin, C.C.; Yu, C.H.; Lo, P.H.; Lyu, J.Y.; Chou, C.C.; Chuang, H.R. A Fully Integrated 60-GHz CMOS Direct-Conversion Doppler Radar RF Sensor with Clutter Canceller for Single-Antenna Noncontact Human Vital-Signs Detection. IEEE Trans. Microw. Theory Tech. 2016, 64, 1018–1028. [Google Scholar] [CrossRef]
  18. Girbau, D.; Lazaro, A.; Ramos, A.; Villarino, R. Remote sensing of vital signs using a doppler radar and diversity to overcome null detection. IEEE Sens. J. 2012, 12, 512–518. [Google Scholar] [CrossRef]
  19. Mahafza, B.R. Radar Systems Analysis and Design Using MATLAB Third Edition; CRC Press: Boca Raton, FL, USA, 2013. [Google Scholar]
  20. Van, N.T.P.; Hasan, S.F.; Gui, X.; Mukhopadhyay, S.; Tran, H. Three-Step Two-Way Decode and Forward Relay With Energy Harvesting. IEEE Commun. Lett. 2017, 21, 857–860. [Google Scholar] [CrossRef]
  21. Jang, B.J.; Wi, S.H.; Yook, J.G.; Lee, M.Q.; Lee, K.J. Wireless bio-radar sensor for heartbeat and respiration detection. Prog. Electromagn. Res. C 2008, 5, 149–168. [Google Scholar]
  22. Droitcour, A.D.; Boric-Lubecke, O.; Lubecke, V.M.; Lin, J.; Kovacs, G.T. Range correlation and I/Q performance benefits in single-chip silicon Doppler radars for noncontact cardiopulmonary monitoring. IEEE Trans. Microw. Theory Tech. 2004, 52, 838–848. [Google Scholar] [CrossRef]
  23. Budge, M.; Burt, M. Range correlation effects on phase and amplitude noise. In Proceedings of the Southeastcon’93, Charlotte, NC, USA, 4–7 April 1993. [Google Scholar]
  24. Kreh, M. Bessel Functions: Project for the Penn State—Gttingen Summer School on Number Theory. 2012. Available online: http://www.math.psu.edu/papikian/Kreh.pdf (accessed on 29 July 2012).
  25. Withers, C.S.; Nadarajah, S. On the product of gamma random variables. Qual. Quant. 2013, 47, 545–552. [Google Scholar] [CrossRef]
  26. Jeffrey, A.; Zwillinger, D. Table of Integrals, Series, and Products; Academic Press: Cambridge, MA, USA, 2007. [Google Scholar]
  27. Grandjean, G.; Baghdadi, N.; Paillou, P.; Dreuillet, P.; Dubois, P.; Souyris, J.; Achache, J. Radar penetration in soils: Towards a new system for subsurface Earth observation. In Proceedings of a CEOS SAR Workshop, Toulouse, France, 26–29 October 1999; European Space Agency and CNES: Paris, France, 2000; Volume 450, p. 267. [Google Scholar]
  28. Tariq, A. Vital Signs Monitoring Using Doppler Radar and on-Body Antennas. Ph.D. Thesis, University of Birmingham, Birmingham, UK, 2013. [Google Scholar]
Figure 1. System diagram.
Figure 1. System diagram.
Sensors 18 00694 g001
Figure 2. Power of signal and noise of the system.
Figure 2. Power of signal and noise of the system.
Sensors 18 00694 g002
Figure 3. Signal-to-noise ratio (SNR) versus the distance from the radar system to the human body.
Figure 3. Signal-to-noise ratio (SNR) versus the distance from the radar system to the human body.
Sensors 18 00694 g003
Figure 4. Detection probabilities of the system under different threshold levels of receiver filter.
Figure 4. Detection probabilities of the system under different threshold levels of receiver filter.
Sensors 18 00694 g004
Figure 5. False alarm of the system under different threshold levels of receiver filter.
Figure 5. False alarm of the system under different threshold levels of receiver filter.
Sensors 18 00694 g005
Figure 6. False alarm/detection probabilities of the system under Nakagami-2 channel model.
Figure 6. False alarm/detection probabilities of the system under Nakagami-2 channel model.
Sensors 18 00694 g006
Figure 7. Measurement setup to detect the breathing rate.
Figure 7. Measurement setup to detect the breathing rate.
Sensors 18 00694 g007
Figure 8. Receiving signals in time domain.
Figure 8. Receiving signals in time domain.
Sensors 18 00694 g008
Figure 9. Receiving signals in frequency domain.
Figure 9. Receiving signals in frequency domain.
Sensors 18 00694 g009
Figure 10. Signal-to-noise ratio (SNR) of the system.
Figure 10. Signal-to-noise ratio (SNR) of the system.
Sensors 18 00694 g010
Figure 11. Detection capability of the system in time domain.
Figure 11. Detection capability of the system in time domain.
Sensors 18 00694 g011
Figure 12. Detection probability of the system at different distances.
Figure 12. Detection probability of the system at different distances.
Sensors 18 00694 g012
Table 1. System parameters used for simulation.
Table 1. System parameters used for simulation.
SymbolDescriptionValue
P T Power of transmitter0 dBm
G T Transmitting antenna gain5 dBi
G R Receiving antenna gain5 dBi
G R X Gain of the receiver10 dB
ρ r Radiated efficiency of receiving antenna0.8
ρ t Radiated efficiency of transmitting antenna0.8
rReflection coefficient of the human body0.5
σ a RCS of the human’s chest500 mm 2
N F Received noise figure6 dB
P 1 / f ( 1 ) 1 / f power noise power at 1 Hz 130 dBm/Hz
TAbsolute noise temperature300
f H High frequency5 Hz
f L Low frequency0.1 Hz
S ( 1 ) Phase noise at 1 Hz intercept60 dB/Hz
η l e a k Leakage between transmitting and receiving antennae 20 dB
α Path loss2
Table 2. Respiratory rates and signal-to-noise ratio (SNR) from five objects at different distances.
Table 2. Respiratory rates and signal-to-noise ratio (SNR) from five objects at different distances.
ObjectBreathing RateSNR
Ref1 m2 m3 m4 m1 m2 m3 m4 m
Sim.Meas.Sim.Meas.Sim.Meas.Sim.Meas.
1151617181843.7343.3733.6721.7627.1213.522.3115.77
2161715181643.7343.1333.6733.0327.1227.3222.3121.77
3161818191743.7340.1133.6740.4427.1230.6522.3127.08
4171816171343.7366.2533.6743.1927.1224.5522.3119.55
5161816171643.7334.0933.6744.5627.1227.1222.3118.65
Mean
values
1617.416.417.81643.7345.3933.6736.5927.1224.6222.3120.56

Share and Cite

MDPI and ACS Style

Phuoc Van, N.T.; Tang, L.; Mukhopadhyay, S.C.; Nguyen, D.M.; Hasan, F. Probabilities of False Alarm for Vital Sign Detection on the Basis of a Doppler Radar System. Sensors 2018, 18, 694. https://doi.org/10.3390/s18030694

AMA Style

Phuoc Van NT, Tang L, Mukhopadhyay SC, Nguyen DM, Hasan F. Probabilities of False Alarm for Vital Sign Detection on the Basis of a Doppler Radar System. Sensors. 2018; 18(3):694. https://doi.org/10.3390/s18030694

Chicago/Turabian Style

Phuoc Van, Nguyen Thi, Liqiong Tang, Subhas Chandra Mukhopadhyay, Duc Minh Nguyen, and Faraz Hasan. 2018. "Probabilities of False Alarm for Vital Sign Detection on the Basis of a Doppler Radar System" Sensors 18, no. 3: 694. https://doi.org/10.3390/s18030694

APA Style

Phuoc Van, N. T., Tang, L., Mukhopadhyay, S. C., Nguyen, D. M., & Hasan, F. (2018). Probabilities of False Alarm for Vital Sign Detection on the Basis of a Doppler Radar System. Sensors, 18(3), 694. https://doi.org/10.3390/s18030694

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop