Next Article in Journal
Raman Spectroscopy and Spectral Signatures of AlScN/Al2O3
Previous Article in Journal
Performance Analysis of Linearly Arranged Concentric Circular Antenna Array with Low Sidelobe Level and Beamwidth Using Robust Tapering Technique
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Real-Time Heart Rate Detection Method Based on 77 GHz FMCW Radar

1
College of Artificial Intelligence, North China University of Science and Technology, Tangshan 063210, China
2
Hebei Key Laboratory of Industrial Intelligent Perception, Tangshan 063210, China
3
College of Management, North China University of Science and Technology, Tangshan 063210, China
*
Author to whom correspondence should be addressed.
Micromachines 2022, 13(11), 1960; https://doi.org/10.3390/mi13111960
Submission received: 15 September 2022 / Revised: 31 October 2022 / Accepted: 7 November 2022 / Published: 11 November 2022

Abstract

:
This paper proposes a real-time heart rate detection method based on 77 GHz FMCW radar. Firstly, the method establishes a new motion model according to respiratory and heartbeat rules, and extracts the motion signals of the chest and the abdomen; then, the random body motion (RBM) signal is eliminated by a combination of polynomial fitting and recursive least squares (RLS) adaptive filtering; lastly, multi-detection-point adaptive harmonics cancellation (AHC) is used to eliminate respiratory harmonics. In addition, the method introduces a spectrum analysis algorithm based on linear predictive coding (LPC). The experimental results show that the method can effectively eliminate the RBM signal and respiratory harmonics, and that the average real-time heart rate detection error rate is 2.925%.

1. Introduction

Heart rate is one of the important vital signs parameters: it can not only reflect the health status of the human body, but can also provide a reference for clinical diagnosis [1,2]. Furthermore, heart rate variability can reflect the activity, balance and related pathological states of the cardiac autonomic nervous system [3]. As a non-contact detection instrument, Doppler radar can detect heartbeat signals from micro-motion signals on the surface of the human body, and it has attracted more and more attention in the medical field [4,5,6,7].
However, there are two problems for Doppler radar, in detecting heart rate in near real-time: the first is how to eliminate the RBM signal, which seriously affects the accuracy of heart rate extraction [8]; the second is how to eliminate respiratory harmonics. Respiratory harmonics and heartbeat are relatively close in frequency; therefore, respiratory harmonics may be mistaken for heartbeat signals [9].
In references [10,11,12], the RBM signal was eliminated using two systems; however, it was difficult for the two systems to operate synchronously without interfering with one other. Gu et al. [13] proposed a method to remove the RBM signal using deep neural networks, which successfully extracted the respiration rate; however, it was difficult to extract the heart rate due to lack of training data. Yang et al. [14] proposed a scheme combining adaptive noise cancellation with polynomial fitting, which would effectively eliminate the RBM signal; however, the step size factor was fixed in the adaptive noise cancellation algorithm, which made it difficult to balance steady-state error and convergence speed.
In references [15,16], to eliminate respiratory harmonics, the respiratory harmonics cancellation method and the adaptive harmonics comb notch digital filter method were proposed; however, both these methods required the respiratory rate to eliminate the harmonics, which increased the dependence of the heart rate on the respiratory rate. In reference [17], to distinguish between respiratory harmonics and heartbeat signals, the heartbeat signals were amplified; however, the respiratory harmonics were also amplified, which reduced the accuracy of the heart rate extraction. Yang et al. [18] proposed a wave-form-driven matched filtering method based on polynomial fitting; however, it was difficult to extract the heart rate quickly.
In addition to the above problems, real-time heart rate detection also needs to solve the problem of insufficient frequency resolution. To this end, Lee et al. [19] proposed an algorithm based on redistribution of the joint time–frequency transform. Hu et al. [20] applied a continuous wavelet filter and ensemble empirical mode decomposition to recover and separate cardiopulmonary signals; they used the peak-to-peak interval to detect the frequency. References [5,21] proposed the time-window-variation technique combined with the fast Fourier transform and wavelet transform, respectively. Park et al. [22] proposed a novel polyphase-basis discrete cosine transform. Ye et al. [23] proposed a stochastic gradient algorithm based on the time-window-variation technique. Lv et al. [24] proposed a non-contact short-term heart rate detection system based on a 120 GHz narrow beam FMCW millimeter wave radar. Gao et al. [25] proposed a method for extracting the life activity spectrum based on millimeter wave radar; although this method was able to extract the heart rate in near real-time, the acquisition time was long, and the accuracy of the heart rate left room for improvement.
In order to solve the above problems, this paper proposes a real-time heart rate detection method based on 77 GHz FMCW radar, and conducts experimental verification. The method establishes a new motion model according to respiratory and heartbeat rules, eliminates the RBM signal by a combination of polynomial fitting and RLS adaptive filtering, and eliminates respiratory harmonics by multi-detection-point AHC. In addition, a spectrum analysis algorithm based on LPC has been introduced, to improve detection accuracy.

2. Theory

Figure 1 shows the block diagram of a typical FMCW radar structure. The transmitted signal can be expressed as:
S ( t ) = cos ( 2 π f c t + π B T c t 2 )
where f c is the carrier frequency, B is the total bandwidth and T c is the duration. Assume the initial distance between the radar and the body is d 0 ; then, the received signal can be expressed as:
R ( t ) = cos 2 π f c t 2 d 0 + x ( t ) c + π B T c t 2 d 0 + x ( t ) c 2
where c is the speed of light, and x ( t ) is the motion of the human body; then, the intermediate frequency signal can be expressed as:
R I F ( t ) = exp j 2 π 2 B d 0 c T c t + 4 π d 0 + x ( t ) λ
where λ is the maximum wavelength of the signal. The sampled signal can be expressed as:
R I F ( m T s + n T f ) = exp j 2 π 2 B d 0 c T C n T f + 4 π d 0 + x ( m T s ) λ
where m and n represent the index number of chirps and the index number of the sampling point in each chirp, respectively; T s is the slow time sampling period, and T f is the fast time sampling period.

3. Methods

3.1. New Motion Model

Figure 2 shows the structure of the chest and diaphragm. When you inhale, the chest and diaphragm expand outward. When you exhale, the chest and diaphragm contract inward. The motion of the chest can be expressed as:
x 1 ( t ) = a r sin ( 2 π f r t ) + a h sin ( 2 π f h t ) + n t ( t )
where a r is the thoracic amplitude due to respiration, f r is the respiration rate, a h is the thoracic amplitude due to heartbeat, f h is the heart rate, and n t ( t ) is the noise signal; n t ( t ) contains the RBM signal. The motion of the abdomen can be expressed as:
x 2 ( t ) = a a sin ( 2 π f r t ) + n a ( t )
where a a is the abdominal amplitude due to respiration, and n a ( t ) is the noise signal; n a ( t ) contains the RBM signal.
Assume that the initial distance from the radar to the chest and abdomen are d 1 and d 2 , respectively. The chest and abdomen belong to the far field targets. After introducing the multi-receiving antenna model, Formula (4) can be rewritten as:
R I F ( m T s + n T f , l ) = k = 1 2 exp j 2 π 2 B d k c T C n T f + 4 π d k + x k ( m T s ) λ + 2 π ( l 1 ) d sin θ k λ
where l represents the index number of receiving antennas, d is the distance between adjacent receiving antennas, and θ 1 and θ 2 are the respective azimuth angles from the chest and from the abdomen to the radar.

3.2. Acquisition of Vital Signs

This paper, to obtain the distance-Doppler-angle dimension data, performed fast time dimension FFT, slow time dimension FFT and angle calculation on R I F ( m T s + n T f , l ) . The Doppler dimensional effective spectral lines were calculated according to the rules of respiration and heartbeat. The formula for calculating the lower and upper spectral lines can be expressed as:
p 1 = min 4 M T a a r t r max λ , 4 M T a a a t r max λ
p 2 = max 4 M T a a r t r min λ , 4 M T a a a t r min λ
where M is the number of chirps, T a is the time interval, t r max is the maximum respiratory cycle, and t r min is the minimum respiratory cycle. The energy of the Doppler slices was accumulated between the Doppler dimensions [ p 1 p 2 ]; the energy value of each point can be expressed as:
E r , a = p = p 1 p 2 M r , p , a r , p , a
where M r , p , a is distance, Doppler slice, and angle dimension data. Finally, the data after energy accumulation was used for distance and angle dimension CFAR detection.

3.3. Elimination of the RBM Signal

Figure 3 shows the block diagram of the RBM signal elimination. A polynomial fitting algorithm was used to fit the RBM signal from x i ( n ) , denoted x ^ i ( n ) . Assuming that the tap-weight vector of the digital filter is w i ( n 1 ) , then the estimation error vector can be expressed as:
e i ( n ) = x i ( n ) y ^ i ( n ) = x i ( n ) w i ( n 1 ) T x ^ i ( n )
where T denotes transposition. The error-weighted sum of the squares of e i ( n ) can be expressed as:
s i n = k = 1 n σ n k e i ( n ) 2
where σ is the forgetting factor. At this point, the tap-weight vector can be expressed as:
w i ( n ) = w i ( n 1 ) + R i 1 ( n ) x ^ i ( n ) e i ( n )
where R i 1 ( n ) is a non-fixed step factor. When the value of s i n is minimum, the cycle will end and the estimation error vector e i ( n ) is output.

3.4. Elimination of Respiratory Harmonics

Figure 4 shows the block diagram of respiratory harmonics elimination. The motion signal of the chest after removing the RBM signal is denoted by e 1 ( n ) , which contains respiratory signals, respiratory harmonics and heartbeat signals. The motion signal of the abdomen after removing the RBM signal is denoted e 2 ( n ) , which contains respiratory signal and respiratory harmonics. The second, third and fourth harmonics of respiration are easily confused with the heartbeat signal; therefore, e 1 ( n ) and e 2 ( n ) are decomposed into five IMF (i.e., I M F 1 1 , I M F 1 2 , I M F 1 3 , I M F 1 4 and I M F 1 5 ) and four IMF (i.e., I M F 2 1 , I M F 2 2 , I M F 2 3 and I M F 2 4 ), respectively. Signal components with frequencies between [0.6–2.5 Hz] are reconstructed separately and denoted h t , r ( n ) + h t , h ( n ) and h a , r ( n ) . When performing harmonic elimination processing, h t , r ( n ) + h t , h ( n ) and h a , r ( n ) · a c are used as the original signal and the reference signal, respectively (the adaptive filtering process based on RLS is detailed in Section 3.3); a c is the energy weight, which can be expressed as:
a c = max E h t , r ( n ) + h t , h ( n ) / max E h a , r ( n )
where E   represents the energy of each signal component; when the cycle ends, an approximate heartbeat signal can be obtained. Using LPC [26] to extend h t , h ( n ) , it can be expressed as:
h ^ t , h ( n ) = l = 1 L a l h t , h ( n l )
where L is the model order, and a l is the prediction coefficient. Finally, analyzing the spectrum of h ^ t , h ( n ) can obtain the heart rate.

4. Results and Discussion

We conducted the experimental validation using the commercial Texas Instruments IWR1443BOOST mmWave radar sensor (Texas Instruments, Dallas, TX, USA), which can form an FMCW radar in a MIMO system. The system parameters were: f c = 77   GHz ; B = 3.98   G H z ; T c = 52   μ s ; λ = 3.9   m m ; N = 512 ; T f = 0.07   μ s ; M = 256 ; T s = 50   m s . Figure 5 shows the experimental scene. During data collection, each subject sat directly in front of the radar, typing on the keyboard with his/her left hand while moving the mouse with his/her right hand. The subject’s left and right arms performed irregular movements, and the subject’s chest and abdomen slowly shook back and forth, which simulated the motion state of office workers in a real environment as far as possible. At the same time, there was a medical finger-pressure pulse sensor YX306 (Yuwell, Suzhou, China) on the finger, which displayed the number of heartbeats. When the value of the figure-pressure pulse sensor was stable, the radar started to collect data and record the value displayed by the finger-pressure pulse sensor every second (reference heart rate). As the data collected by the radar and the numerical value recorded by the figure-pressure pulse sensor needed to be synchronized, this may also have caused small errors. Additionally, we measured the distance with a tape measure, and calculated the angle (expected result). In order to verify the real-time heart rate extraction performance, we collected data for 12.8 s each time, and implemented the method in MATLAB. To quantify the performance of the method, the error rate and the average error rate were introduced into the performance analysis, which can be expressed as:
e r r o r = H R m e a H R r e a H R r e a × 100 %
error ¯ = ( l = 1 L e r r o r l ) / L
where H R m e a is the measured heart rate, H R r e a is the average heart rate recorded by the figure-pressure pulse sensor, and L is the number of measurement groups. It should be noted that the radar board was placed vertically, and clockwise was positive.
Figure 6 shows a set of experimental results. The chest and abdomen belong to the far field targets. Figure 6a shows the localization results of the chest and abdomen, the distance and angle of which were: 2 1 ° , 45.91   c m , 4 4 ° and 28.69   c m . These results were in fair agreement with the expected values. Figure 6b,c show the motion signals of the chest and abdomen, respectively, after removing the RBM signal. As can be seen, the RBM signal has been removed and the baseline is effectively back on the central axis. To better demonstrate the performance of real-time heart rate extraction, the 12.8 s signal was equally divided into four segments (i.e., segment 1, segment 2, segment 3, segment 4). Figure 6d–g show the results of the four-segment elimination of the respiratory harmonics. As can be seen from the figure, the motion signals of the chest and abdomen contained respiratory harmonics of the same frequency. Observing the results, it can be seen that the respiratory harmonics were effectively eliminated. The signal after harmonic elimination was linearly predicted based on LPC, and the heart rates obtained by spectrum analysis were: 1.527 Hz, 1.464 Hz, 1.412 Hz and 1.465 Hz. The error rates were: 3.74%, 1.74%, 0.07% and 3.83%, respectively. The average error rate was 2.35% and the computer calculation time was 3.01 s (MagicBook 16 Pro).
We validated the proposed approach by conducting 15 experiments on five subjects, differing in height (172–188 cm) and in age (26–49 years). During the experiment, each subject sat directly in front of the radar and used his/her computer. We collected data for 12.8 s each time. Table 1 shows the results of the experimental validation. The average error rates of the four segments were: 3.10%, 3.37%, 2.68% and 2.55%. The average error rate for all experiments was 2.925%. At the same time, the computer running time was less than 3.01 s. This shows that the proposed method can effectively extract the real-time heart rate.
Table 2 lists the comparison results with other works. We set the distance between 0.28–0.70 m to simulate an office scene realistically. When detecting the real-time heart rate, the time window length was generally required to be less than 5 s. In order to better verify the performance of the proposed method in detecting the real-time heart rate, we set a 3.2 s time window for detecting the heart rate. In general, the shorter the time window, the lower the accuracy of the heart rate detection. As can be seen from the table, the biggest advantage of the proposed method is that the required detection time is short, and accurate heart rate detection can be performed at any time.

5. Conclusions

In order to improve the application of non-contact detection technology in the medical field, this paper proposes a real-time heart rate detection method based on 77 GHz FMCW radar. The method eliminates the RBM signal by a combination of polynomial fitting and RLS adaptive filtering, and eliminates respiratory harmonics by a multi-detection-point AHC. The heart rate is then obtained using an LPC-based spectral analysis algorithm. Through experiments, it was found that the proposed method effectively eliminated the RBM signal and respiratory harmonics, and detected the real-time heart rate of office workers. The proposed method has great potential for the real-time heart rate detection of key personnel, such as office workers, drivers, etc.

Author Contributions

Conceptualization, Z.J.; methodology, Z.J.; software, Z.J.; validation, X.H. and Z.J.; formal analysis, X.H. and Z.J.; investigation, X.H. and Z.J.; resources, Z.J.; data curation, Z.J.; writing—original draft preparation, Z.J.; writing—review and editing, X.H. and Z.J.; visualization, Z.J.; supervision, R.Z.; project administration, X.H. and R.Z.; funding acquisition, R.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Science and Technology project of Hebei Education Department (Project No. ZD2020152), the graduate teaching quality engineering project of Hebei Education Department (Project No. SFK202201).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the plots within this paper, and other findings of this study, are available from the second author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cavagnaro, M.; Pittella, E.; Pisa, S. UWB pulse propagation into human tissues. Phys. Med. Biol. 2013, 58, 8689–8707. [Google Scholar] [CrossRef]
  2. Pimentel, M.A.; Redfern, O.C.; Hatch, R.; Young, J.D.; Tarassenko, L.; Watkinson, P.J. Trajectories of vital signs in patients with COVID-19. Resuscitation 2020, 156, 99–106. [Google Scholar] [CrossRef] [PubMed]
  3. Khan, A.A.; Lip, G.Y.; Shantsila, A. Heart rate variability in atrial fibrillation: The balance between sympathetic and parasympathetic nervous system. Eur. J. Clin. Investig. 2019, 49, e13174. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Sun, L.; Hong, H.; Li, Y.; Gu, C.; Xi, F.; Li, C.; Zhu, X. Noncontact Vital Sign Detection based on Stepwise Atomic Norm Minimization. IEEE Signal Process. Lett. 2015, 22, 2479–2483. [Google Scholar] [CrossRef]
  5. Li, M.; Lin, J. Wavelet-Transform-Based Data-Length-Variation Technique for Fast Heart Rate Detection Using 5.8-GHz CW Doppler Radar. IEEE Trans. Microw. Theory Tech. 2018, 66, 568–576. [Google Scholar] [CrossRef]
  6. Nosrati, M.; Tavassolian, N. High-Accuracy Heart Rate Variability Monitoring Using Doppler Radar Based on Gaussian Pulse Train Modeling and FTPR Algorithm. IEEE Trans. Microw. Theory Tech. 2018, 66, 556–567. [Google Scholar] [CrossRef]
  7. Yamamoto, K.; Toyoda, K.; Ohtsuki, T. Spectrogram-Based Non-Contact RRI Estimation by Accurate Peak Detection Algorithm. IEEE Access 2018, 6, 60369–60379. [Google Scholar] [CrossRef]
  8. Li, C.; Lin, J. Random Body Movement Cancellation in Doppler Radar Vital Sign Detection. IEEE Trans. Microw. Theory Tech. 2008, 56, 3143–3152. [Google Scholar]
  9. Massagram, W.; Lubecke, V.M.; Host-Madsen, A.; Boric-Lubecke, O. Assessment of Heart Rate Variability and Respiratory Sinus Arrhythmia via Doppler Radar. IEEE Trans. Microw. Theory Tech. 2009, 57, 2542–2549. [Google Scholar] [CrossRef]
  10. Wang, K.F.; Horng, T.S.; Peng, K.C.; Jau, J.K.; Li, J.Y.; Chen, C.C. Single-Antenna Doppler Radars Using Self and Mutual Injection Locking for Vital Sign Detection With Random Body Movement Cancellation. IEEE Trans. Microw. Theory Tech. 2011, 59, 3577–3587. [Google Scholar] [CrossRef]
  11. Gu, C.; Wang, G.; Li, Y.; Lnoue, T.; Li, C. A Hybrid Radar-Camera Sensing System With Phase Compensation for Random Body Movement Cancellation in Doppler Vital Sign Detection. IEEE Trans. Microw. Theory Tech. 2013, 61, 4678–4688. [Google Scholar] [CrossRef]
  12. Tang, M.C.; Kuo, C.Y.; Wun, D.C.; Wang, F.K.; Horng, T.S. A Self- and Mutually Injection-Locked Radar System for Monitoring Vital Signs in Real Time With Random Body Movement Cancellation. IEEE Trans. Microw. Theory Tech. 2016, 64, 4812–4822. [Google Scholar] [CrossRef]
  13. Gu, C.; Wang, J.; Lien, J. Deep Neural Network based Body Movement Cancellation for Doppler Radar Vital Sign Detection. In Proceedings of the 2019 IEEE MTT-S International Wireless Symposium (IWS), Guangzhou, China, 19–22 May 2019; pp. 1–3. [Google Scholar]
  14. Yang, Z.K.; Shi, H.; Zhao, S.; Huang, X.D. Vital sign detection during large-scale and fast body movements based on an adaptive noise cancellation algorithm using a single Doppler radar sensor. Sensors 2020, 20, 4183. [Google Scholar] [CrossRef] [PubMed]
  15. Ti, J.; Lin, J. Respiration harmonics cancellation for Accurate Heart Rate measurement in non-contact vital sign detection. In Proceedings of the 2013 IEEE MTT-S International Microwave Symposium Digest (MTT), Seattle, WA, USA, 2–7 June 2013; pp. 1–3. [Google Scholar]
  16. Huang, T.Y.; Hayward, L.F.; Lin, J. Adaptive harmonics comb notch digital filter for measuring heart rate of laboratory rat using a 60-GHz radar. In Proceedings of the 2016 IEEE MTT-S International Microwave Symposium (IMS), San Francisco, CA, USA, 22–27 May 2016; pp. 1–4. [Google Scholar]
  17. Nosrati, M.; Tavassolian, N. Accurate Doppler Radar-Based Cardiopulmonary Sensing Using Chest-Wall Acceleration. IEEE J. Electromagn. RF Microw. Med. Biol. 2019, 3, 41–47. [Google Scholar] [CrossRef]
  18. Yang, Z.K.; Zhao, S.; Huang, X.D.; Lv, W. Accurate Doppler radar-based heart rate measurement using matched filter. IEICE Electron. Express 2020, 17, 20200062. [Google Scholar] [CrossRef] [Green Version]
  19. Lee, C.; Yoon, C.; Kong, H.J.; Kim, H.C.; Kim, Y. Heart Rate Tracking Using a Doppler Radar With the Reassigned Joint Time-Frequency Transform. IEEE Antennas Wirel. Propag. Lett. 2011, 10, 1096–1099. [Google Scholar]
  20. Hu, W.; Zhao, Z.; Wang, Y.; Zhang, H.; Lin, F. Noncontact Accurate Measurement of Cardiopulmonary Activity Using a Compact Quadrature Doppler Radar Sensor. IEEE Trans. Biomed. Eng. 2014, 61, 725–735. [Google Scholar] [CrossRef]
  21. Tu, J.; Lin, J. Fast Acquisition of Heart Rate in Noncontact Vital Sign Radar Measurement Using Time-Window-Variation Technique. IEEE Trans. Instrum. Meas. 2015, 65, 112–122. [Google Scholar] [CrossRef]
  22. Park, J.; Ham, J.W.; Park, S.; Kim, D.H.; Park, S.J.; Kang, H.; Park, S.O. Polyphase-Basis Discrete Cosine Transform for Real-Time Measurement of Heart Rate With CW Doppler Radar. IEEE Trans. Microw. Theory Tech. 2018, 66, 1644–1659. [Google Scholar] [CrossRef]
  23. Ye, C.; Toyoda, K.; Ohtsuki, T. A stochastic gradient approach for robust heartbeat detection with Doppler radar using time-window-variation technique. IEEE Trans. Biomed. Eng. 2019, 66, 1730–1741. [Google Scholar] [CrossRef]
  24. Lv, W.; Zhao, Y.; Zhang, W.; Liu, W.; Hu, A.; Miao, J. Remote Measurement of Short-Term Heart Rate With Narrow Beam Millimeter Wave Radar. IEEE Access 2021, 9, 165049–165058. [Google Scholar] [CrossRef]
  25. Gao, Z.; Ali, L.; Wang, C.; Liu, R.; Wang, C.; Qian, C.; Sung, H.; Meng, F. Real-Time Non-Contact Millimeter Wave Radar-Based Vital Sign Detection. Sensors 2022, 22, 7560. [Google Scholar] [CrossRef] [PubMed]
  26. Ding, Y.; Yu, X.; Lei, C.; Sun, Y.; Xu, x.; Zhang, J. A Novel Real-Time Human Heart Rate Estimation Method for Noncontact Vital Sign Radar Detection. IEEE Access 2020, 8, 88689–88699. [Google Scholar] [CrossRef]
Figure 1. Block diagram of the non-contact FMCW radar. SG: Signal Generator; VCO: Voltage Controlled Oscillator; PA: Power Amplifier; Tx: Transmitting antenna; Rx: Receiving antenna; LNA: Low Noise Amplifier; MIXER: Mixer; ADC: Analog Digital Converter; DSP: Digital Signal Processing.
Figure 1. Block diagram of the non-contact FMCW radar. SG: Signal Generator; VCO: Voltage Controlled Oscillator; PA: Power Amplifier; Tx: Transmitting antenna; Rx: Receiving antenna; LNA: Low Noise Amplifier; MIXER: Mixer; ADC: Analog Digital Converter; DSP: Digital Signal Processing.
Micromachines 13 01960 g001
Figure 2. Structure of the chest and diaphragm.
Figure 2. Structure of the chest and diaphragm.
Micromachines 13 01960 g002
Figure 3. Block diagram of the RBM signal elimination: x i ( n ) is the motion signal of the human body; RLS: recursive least squares; x 1 ( n ) is the motion signal of the chest; x 2 ( n ) is the motion signal of the abdomen.
Figure 3. Block diagram of the RBM signal elimination: x i ( n ) is the motion signal of the human body; RLS: recursive least squares; x 1 ( n ) is the motion signal of the chest; x 2 ( n ) is the motion signal of the abdomen.
Micromachines 13 01960 g003
Figure 4. Block diagram of respiratory harmonics signal elimination: e 1 ( n ) is the motion signal of the chest after removing the RBM signal; e 2 ( n ) is the motion signal of the abdomen after removing the RBM signal. VMD: Variational Mode Decomposition; IMF: Intrinsic mode function; AHC: adaptive harmonics cancellation.
Figure 4. Block diagram of respiratory harmonics signal elimination: e 1 ( n ) is the motion signal of the chest after removing the RBM signal; e 2 ( n ) is the motion signal of the abdomen after removing the RBM signal. VMD: Variational Mode Decomposition; IMF: Intrinsic mode function; AHC: adaptive harmonics cancellation.
Micromachines 13 01960 g004
Figure 5. Experimental environment.
Figure 5. Experimental environment.
Micromachines 13 01960 g005
Figure 6. Experimental results: (a) localization results of the chest and abdomen; (b) the motion signal of the chest after removing the RBM signal; (c) the motion signal of the abdomen after removing the RBM signal; (d) the result of the segment 1 elimination of respiratory harmonics; (e) the result of the segment 2 elimination of respiratory harmonics; (f) the result of the segment 3 elimination of respiratory harmonics; (g) the result of the segment 4 elimination of respiratory harmonics.
Figure 6. Experimental results: (a) localization results of the chest and abdomen; (b) the motion signal of the chest after removing the RBM signal; (c) the motion signal of the abdomen after removing the RBM signal; (d) the result of the segment 1 elimination of respiratory harmonics; (e) the result of the segment 2 elimination of respiratory harmonics; (f) the result of the segment 3 elimination of respiratory harmonics; (g) the result of the segment 4 elimination of respiratory harmonics.
Micromachines 13 01960 g006
Table 1. Results of the experimental validation.
Table 1. Results of the experimental validation.
Segment 1 (3.2 s)Segment 2 (3.2 s)Segment 3 (3.2 s)Segment 4 (3.2 s)
H R m e a / Hz H R r e a / Hz e r r o r / % H R m e a / Hz H R r e a / Hz e r r o r / % H R m e a / Hz H R r e a / Hz e r r o r / % H R m e a / Hz H R r e a / Hz e r r o r / %
11.51241.57223.801.66251.58335.001.48081.52222.721.45841.45830.01
21.44521.47782.211.44611.45000.271.43871.41671.561.42931.41670.89
31.69981.63334.071.64371.59443.091.56251.52222.651.43231.42500.51
41.41761.35005.011.37151.32223.291.41681.32786.701.36971.36110.61
51.52761.46674.151.46381.43891.731.41231.41110.081.46501.40563.82
61.33851.32221.231.44461.32229.361.43051.31678.651.37531.35001.87
71.42521.44171.141.42351.41670.481.39391.40000.441.40581.37222.45
81.43661.38333.851.33471.38333.521.40481.39440.741.34101.38333.06
91.41531.41110.291.41591.45002.351.40471.37921.851.33291.31251.55
101.66601.57785.591.63421.56674.311.46571.49171.741.39951.44443.11
111.37051.45005.481.43121.41111.421.35891.36670.571.34811.37221.76
121.40061.40560.351.41041.44442.801.37331.38891.281.41981.34005.96
131.42151.46673.081.41671.45002.301.36471.38891.741.38641.36671.44
141.36491.30004.991.38541.30006.571.39661.30566.971.36111.29445.15
151.41241.39441.291.39301.33894.041.36681.33332.511.38971.32225.10
error ¯ /%3.10 error ¯ /%3.37 error ¯ /%2.68 error ¯ /%2.55
Table 2. Comparison of the proposed method with the other existing real-time heart rate detection methods.
Table 2. Comparison of the proposed method with the other existing real-time heart rate detection methods.
Ref. No.Detection Distance/mAcquisition Time/sTime Window Length/s error ¯ / %
[19]0.57555
[20]0.524015<3
[21]Not mentioned302–53.4
[5]Not mentioned603–53.5
[22]1–1.180–903, 2, 1.53.75, 5.58, 7.58
[23]0.812084.2
[24]110034.38
[25]<0.5Not mentionedNot mentioned3.65
This work0.28–0.7012.83.22.925
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Huang, X.; Ju, Z.; Zhang, R. Real-Time Heart Rate Detection Method Based on 77 GHz FMCW Radar. Micromachines 2022, 13, 1960. https://doi.org/10.3390/mi13111960

AMA Style

Huang X, Ju Z, Zhang R. Real-Time Heart Rate Detection Method Based on 77 GHz FMCW Radar. Micromachines. 2022; 13(11):1960. https://doi.org/10.3390/mi13111960

Chicago/Turabian Style

Huang, Xiaohong, Zedong Ju, and Rundong Zhang. 2022. "Real-Time Heart Rate Detection Method Based on 77 GHz FMCW Radar" Micromachines 13, no. 11: 1960. https://doi.org/10.3390/mi13111960

APA Style

Huang, X., Ju, Z., & Zhang, R. (2022). Real-Time Heart Rate Detection Method Based on 77 GHz FMCW Radar. Micromachines, 13(11), 1960. https://doi.org/10.3390/mi13111960

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