Model Identification and Physical Exercise Control using Nonlinear Heart Rate Model and Particle Filter

Physical exercise has been proven to be beneficial for both healthy subjects and cardiac patients. It can improve cardiovascular health and promote recovery from various heart conditions. Heart Rate (HR) is a cardiovascular variable, which can be easily monitored and provides important insights about cardiac functions during and after physical exercise. This study presents a HR-based modeling and control framework to monitor physiological changes during exercise, from which the exercise intensity is optimized to capitalize exercise outcomes. HR models were previously developed to investigate exercise physiology, but efficient model identification has not been extensively discussed in the literature. Most existing HR models are nonlinear state-space models, and traditional optimization techniques may fail to provide accurate model identification results. In this work, we propose to use particle filter (PF) to identify HR model parameters and further optimize the intensity of exercise, e.g., walking or running speed, based on the calibrated model. Specifically, sequential importance sampling and resampling (SISR) and smoothing were chosen to estimate state variables, and particle marginal Metropolis-Hastings method was used to identify model parameters from HR observations. In addition, using predictions calculated from the HR model, treadmill speed was optimized by minimizing the difference between predictions and the target HR. The modeling and control framework is validated with different case studies. The results demonstrate that the proposed method is a useful tool for personalized HR modeling and exercise control, which can benefit both regular exercise training and cardiac rehabilitation.

[1]  M. Pitt,et al.  Filtering via Simulation: Auxiliary Particle Filters , 1999 .

[2]  J. Suykens,et al.  Recurrent least squares support vector machines , 2000 .

[3]  Feng Ding,et al.  Hierarchical gradient-based identification of multivariable discrete-time systems , 2005, Autom..

[4]  A. Doucet,et al.  Monte Carlo Smoothing for Nonlinear Time Series , 2004, Journal of the American Statistical Association.

[5]  Maria Zakynthinaki,et al.  Simulating heart rate kinetics during incremental and interval training , 2016 .

[6]  Andrey V. Savkin,et al.  Nonlinear Modeling and Control of Human Heart Rate Response During Exercise With Various Work Load Intensities , 2008, IEEE Transactions on Biomedical Engineering.

[7]  L. Ljung,et al.  Identification of composite local linear state-space models using a projected gradient search , 2002 .

[8]  Peter J. Schwartz,et al.  Heart-rate profile during exercise as a predictor of sudden death. , 2005, The New England journal of medicine.

[9]  Rainer Rauramaa,et al.  Heart rate response during exercise test and cardiovascular mortality in middle-aged men. , 2006, European heart journal.

[10]  Ignacio Refoyo,et al.  A Model of Heart Rate Kinetics in Response to Exercise , 2008 .

[11]  Yong Zhang,et al.  Unbiased identification of a class of multi-input single-output systems with correlated disturbances using bias compensation methods , 2011, Math. Comput. Model..

[12]  Tien C. Hsia,et al.  System identification: Least-squares methods , 1977 .

[13]  Arnaud Doucet,et al.  On Particle Methods for Parameter Estimation in State-Space Models , 2014, 1412.8695.

[14]  M. Fornage,et al.  Heart Disease and Stroke Statistics—2017 Update: A Report From the American Heart Association , 2017, Circulation.

[15]  Tor Arne Johansen,et al.  On Tikhonov regularization, bias and variance in nonlinear system identification , 1997, Autom..

[16]  Lu Wang,et al.  Identification and Control for Heart Rate Regulation During Treadmill Exercise , 2007, IEEE Transactions on Biomedical Engineering.