- Open Access
Modelling and regulating of cardio-respiratory response for the enhancement of interval training
© Haddad et al.; licensee BioMed Central Ltd. 2014
Received: 4 November 2013
Accepted: 22 January 2014
Published: 6 February 2014
The interval training method has been a well known exercise protocol which helps strengthen and improve one’s cardiovascular fitness.
To develop an effective training protocol to improve cardiovascular fitness based on modelling and analysis of Heart Rate (HR) and Oxygen Uptake (VO2) dynamics.
In order to model the cardiorespiratory response to the onset and offset exercises, the (K 4 b 2 , Cosmed) gas analyzer was used to monitor and record the heart rate and oxygen uptake for ten healthy male subjects. An interval training protocol was developed for young health users and was simulated using a proposed RC switching model which was presented to accommodate the variations of the cardiorespiratory dynamics to running exercises. A hybrid system model was presented to describe the adaptation process and a multi-loop PI control scheme was designed for the tuning of interval training regime.
By observing the original data for each subject, we can clearly identify that all subjects have similar HR and VO2 profiles. The proposed model is capable to simulate the exercise responses during onset and offset exercises; it ensures the continuity of the outputs within the interval training protocol. Under some mild assumptions, a hybrid system model can describe the adaption process and accordingly a multi-loop PI controller can be designed for the tuning of interval training protocol. The self-adaption feature of the proposed controller gives the exerciser the opportunity to reach his desired setpoints after a certain number of training sessions.
The established interval training protocol targets a range of 70-80% of HRmax which is mainly a training zone for the purpose of cardiovascular system development and improvement. Furthermore, the proposed multi-loop feedback controller has the potential to tune the interval training protocol according to the feedback from an individual exerciser.
Heart rate (HR) and oxygen uptake (VO2) are key indicators of functional health status; their measurements can aid early detection of cardiac diseases [1, 2]. Furthermore, these cardio respiratory endurances have long been recognized as one of the fundamental components of physical fitness.
To build effective training protocols, VO2 and HR are often applied to assess the exercise intensity for achieving different training goals, such as fat burning or cardiovascular system improvement .
One of the most effective training regimens is called interval training which was first described by Reindell and Roskamm  and was popularized in the 1950s by the Olympic champion, Emil Zatopek. Swedish physiologist Per OlØf Astrand provided the first scientific study  on interval training in 1960. Since then, it has been the basis for athletic training programs for many years.
Nowadays, interval training method becomes a well known exercise protocol which helps strengthen and improve cardiovascular fitness. Previous researchers in this area [6, 7] have investigated the effectiveness of interval training and its significance in improving factors associated with O2 transport along with muscle uptake. Interval training has been shown to raise skeletal muscle enzyme activity [8, 9], improve vascular health , and decrease the risk of cardiovascular-disease in obese adolescents . Additionally, interval training increases proteins that transport fatty acids across the mitochondrial membrane . This study aims to develop an effective interval training protocol for an individual exerciser to improve cardiovascular fitness based on modelling and analysis of HR and VO2 dynamics.
Interval training interleaves high intensity exercises with recovery or rest periods . The idea behind this approach is breaking up a set amount of work into smaller segments rather than performing a greater volume of work at a higher intensity . In addition, it is well known that exercise at varying intensities helps improve exerciser’s aerobic capacity to exercise longer . A well designed interval training protocol can guide the exerciser to timely shift from low to high training intensities and vice versa within the exercise zone to achieve desired exercise targets with safe.
Any interval training protocol has at least three different periods, warm-up, exercise (switching between high intensity period and recovery period), and cool-down. The most important characters for the design of an appropriate interval training protocol are time (or distance), intensity level (speed level in our case at onset segment), time of each recovery period (offset), and number of repetitions within the exercise period. For example, both the time intervals of the onset and offset segments should be properly selected based on their transient properties to avoid reaching a steady state value during interval training.
To achieve better interval training effect, this study first investigates the dynamic responses of VO2 and HR by building a model, which can depict both the onset and offset dynamics. Although existing literatures [15, 16] have studied the HR and/or VO2 responses to the onset and offset exercises, we need to re-explore the onset and offset dynamics with a certain maximum intensity (70 percent of VO2max) for the purpose of improving cardiovascular fitness. We noticed the distinct differences of two dynamic characteristics (time constants and steady state gains) for the onset and offset of running exercises.
A switching RC circuit model for healthy male subjects has then been developed, which can well accommodate these observed differences.
Based on this model, an interval training protocol has been proposed for this group of subjects. In order to adapt to different individual exerciser, a multi-loop integral control scheme has been proposed and achieved desired simulation results. It should be pointed out that in moderate intensity exercise the relationship between HR and VO2 is nearly linear . As the proposed interval training protocol is particularly designed in moderate intensity range, HR can be used as an indicator of exercise strength during the proposed training.
This paper is organized as follows. Section ‘Experimental equipment and exercise protocol’ describes the experimental equipment and exercise protocol setup for model estimation. Section ‘Results and discussion’ presents the main results and discussions including the proposed switching model and interval training protocol. Section ‘Individualized adaptation for the proposed interval training protocol’ describes the adaption framework as well as the multiloop PI controller design. Section ‘Control system verification’ provides a practical method to examine and verify the functionality as well as performance of the proposed controller. Finally, Section ‘Conclusion’ gives the conclusions.
Experimental equipment and exercise protocol
Since both HR and VO2 can vary with the time of the day, emotional state, food and caffeine intake, previous activity, fatigue, temperature, humidity, altitude and dehydration , the participants were asked to have a light meal at least two hours before the experiment and not to engage in intense or prolonged exercise for 24 hours prior to each experiment. Environmental conditions were the same for all participants. The temperature of the laboratory was set at 25 degrees Celsius and the humidity was set at approximately 50%. All subjects were requested to be familiar with an automated treadmill before experiments.
All laboratory analyses were performed using a portable gas analyzer (K 4 b 2 , Cosmed). The analyzer comes with a compatible HR probe which consists of two parts: the elastic belt containing the transmitter, and the receiver.
Using chest electrodes to monitor HR signals is considered to be valid, accurate and reliable during physical activity monitoring .
The K4b2 gas analyzer was used because it has previously been reported to be valid, accurate and reliable .
The experiments have been done on healthy however untrained subjects. A submaximal running protocol, that is practical and easy to apply, was implemented in the setup. It is required that the subjects where to perform at less than 80% of their maximum heart rate (HRmax), this is also a recommendation for unfit people and for those with respiratory or cardiac risks .
The relationship between HR and VO2 becomes non-linear during light and very highly intense activity . Therefore, and for the values of HR and VO2 to remain located in the linear part, a moderate speed of 9 km/h was examined and set to be the highest level in the training protocol; each participant, at varying times, reached a reading of 80% of his HRmax at this speed.
Each subject completed two running exercises on separate occasions. In each session and as a warming up interval, the subjects were requested to walk on a treadmill at 5 km/h for 4 minutes to avoid any kind of stiffness or fatiguing effect may occur during the running interval. Then, they were asked to run at 9 km/h for 6 minutes, before a recovery time or cooling down period of 5 minutes at 5 km/h speed.
During experiments, HR was monitored to ensure the safety of the subjects. Whenever noticeable change in heart rhythm is detected the experiment will be terminated.
Results and discussion
It is known that many internal and environmental factors may have drastic effects on the readings during the experiment , even when all of the laboratory conditions have been standardized and set to be the same for all participants. Some environmental and undetectable factors such as dehydration, heat loss mechanisms as well as altitude may still affect the results. To reduce the influences of these uncontrollable factors, the proposed training protocol was repeated twice and the data was interpolated, averaged and filtered.
Where k and τ are the steady state gain and the time constant respectively.
Estimated time constants and normalized steady state gains of HR response
SS Gain (k1)
Time Constant (τ 1)(sec)
SS Gain (k2)
Time Constant (τ 2)(sec)
Estimated time constants and normalized steady state gains of VO 2 response
SS Gain (k1)
Time constant (τ 1)(sec)
SS Gain (k2)
Time constant (τ 2)(sec)
By observing the original data for each subject we can clearly identify that all subjects have similar HR and VO2 profiles. For each individual subject, steady state gain of offset is smaller than steady state gain of onset for both HR and VO2.
Averaged values of time constants and steady state gains
Based on the results listed in Tables 2, 3 and 4, it can be seen that the time constant of offset is larger than that of onset for both HR and VO2. Paper  has shown the same result for HR response during walking exercises; nevertheless, it did not investigate the dynamic characteristics of a running protocol.
The proposed interval training protocol
In this subsection, a sensible interval training protocol will be built based on the experimental results. To build a training protocol, targets have to be set first.
For instance, if the aim of a training protocol is to develop basic endurance and aerobic capacity, all running intervals should be completed at a maximum of 70% of HRmax while developing cardiovascular system can be done by training at a range of 70% to 80% of HRmax.
Relationship between HR max and VO 2max
Exercise intensity levels that coincide with HR max and VO 2max
Fat burning (low)
Anaerobic (very vigorous)
As mentioned before, a linear relationship between VO2 and HR exists during moderate intense activity and that is why our training protocol will target the aerobic zone (70-80% of HRmax or 50-70% VO2max) which aims to develop the exerciser’s cardiovascular system. Based on equation (2), for an averaged model HRmax = 205.8-0.685 × (average age) = 205.8-0.685 × (30.63) = 184.82 ≈ 185 bpm. Thus, building an interval training which reaches a 148 bpm (80% of HRmax) in each interval will guarantee achieving ∼70% of VO2max.
The experiment results show that for an average model, the exerciser exceeds an 80% of his HRmax at the 297th second which means 57 seconds after the start of the running period, also 1 minute at speed 5km/h is considered to be the recovery period after each running interval.
In order to verify the new proposed interval training protocol shown in Figure 3, more subjects need to be recruited, and this will be the next step of this study.
The parameters of the model are tuned based on Table 4. ROn multiplied by C (ROn ×C) represents the time constant for the onset response, while ROff multiplied by C (ROff ×C) represents the time constant for the offset response.
The gain value can be controlled by manipulation of the voltage sources VOn and VOff. It should be noted that when the steady state gains of the onset and offset dynamics are approximately the same, VOn will be constant and VOff is constant zero. As another special case, if both the onset and offset exercises reach steady state, then both VOn and VOff can be constant.
The switch is utilized to replicate the switching behavior between the two responses based on the training protocol segments. National Instruments Multisim has been used to simulate the proposed interval training protocol.
The experimental results approximately match the simulation results. However, in contrast with HR, the experimental results for VO2 were skewed due to the presence of extensive noise. As the intensity of this protocol is anticipated in the moderate range for the purpose of improving cardiovascular fitness, the relationship between HR and VO2 is nearly linear . Furthermore, HR is easier to be measured and has wider frequency bandwidth. Therefore, HR is suitable and therefore recommended to be used as an indicator of exercise strength for the proposed interval training protocol.
Individualized adaptation for the proposed interval training protocol
The proposed interval training protocol is based on the established average model for eight healthy young male subjects. However, for an individual exerciser (who is young and healthy) the proposed protocol may need to be adjusted due to the differences of the intra- and inter-subjects.
In this section, we present a hybrid system model (a dynamic system that exhibits both continuous and discrete dynamic behaviors) to describe the adaptation process and propose a multi-loop PI (Proportional and Integral) control approach for the tuning of interval training protocol.
The adaptation framework
Specifically, each single experiment can be considered as a discrete event and the dynamics of cardio-respiratory responses to exercise can be depicted by a continuous linear model. It will be shown under modest assumptions, we can simplify this special hybrid system as a simple discrete time system.
The intermission between the interval training experiments can be hours or days. However, in order to simplify the process, we first assume that the intermission is constant (similar as the sampling time of a discrete time system) given that the subject has similar physiological conditions before each training (i.e., the intermission between the interval training is long enough to recover the exercisers).
For each training, the cardio respiratory responses (HR/VO2) are actually the outputs of continuous systems. In order to simplify the system as a static discrete time system, we first pick up some key characteristics from the continuous processes, which are associated with exercise effects, such as mean value and standard deviation. The proposed training protocol (see Figure 3) contains three-period square wave. In this study we pick up the lowest point and the highest point of the third period HR response y(t4) and y(t5) respectively as the reflections of exercise effects (see Figure 6).
We can then treat every single interval training experiment as a simple static system, which has two inputs (the duty cycle (Δ t1) and the period (Δ t1 + Δ t2)) and two outputs (y(t4) and y(t5)). Under the above assumptions and simplification, the overall adaptation process can now be simply treated as a two-input two-output (2I2O) discrete time system.
In order to determine the model of the 2I2O static system, we investigate HR response for the proposed interval training protocol. Figure 6 shows the protocol and its response.
Controller design and simulation study
As the adaptation can be depicted as a discrete time system, we can design a discrete time multi-loop controller (Proportional and Integral controller) to adjust the two inputs (the duty cycle (Δ t1) and the period (Δ t1 + Δ t2)) to regulate the two outputs (y(t4) and y(t5)) to reach the desired reference values.
Up to now, we can see even though our aim is to regulate a continuous process; we proposed a simple discrete time control technique to fulfill the goal.
During interval training of an individual subject, at the first iteration the exerciser is asked to run under the predefined training protocol. HR response at t4 and t5 (as shown in Figure 6) will then be measured and feedbacked to the multi-loop control system to update controller outputs’ inputs u1 and u2 (the duty cycle (Δ t1) and the period (Δ t1 + Δ t2)).
This modification will try to adjust the measured output values (y(t4) and y(t5)) to approach the desired setpoints for the next training exercise. This recursive process will be repeated until we finally reach the desired set points.
In simulation, from Figure 8, we conclude that after at least 12 iterations (12 training exercises) the system will reach the desired setpoints.
Since our controller is based on discrete events regulating a continuous process, it can perform one time between two training experiments. Therefore, it is very easy to be implemented in low cost portable devices which have limited computation power.
Control system verification
This section provides a practical method to examine and verify the functionality as well as performance of the proposed controller. The new interval training protocol combined with the designed controller have been implemented and tested on stair climbing exercise as a free environment activity. This can show how the proposed interval training protocol that has been created earlier can be implemented using different exercising techniques, such as stairs climbing, cycling and even swimming.
The setup, concept and structure of the interval training remain the same. However, the change would be within the type of exercise during the onset and offset stages. Subjects in the new implementation will be requested to ascend a staircase of steps without any pause; the exercisers will attempt to achieve the training protocol without holding the banisters or the handrail.
The experiments were carried out at free speed and the participants were given the instruction to walk at their normal pace. However, each exerciser was instructed to stay within 4-5 km/h climbing speed. This will represent the onset stage in our training protocol. On the other hand, offset stages can be represented as walking aside on the same step the exercisers reach when the onset period finishes; this continuous walking on the same step is the perfect match of the offset stage in the treadmill interval training exercise. All participants were instructed to use their left legs for the first step, to only place one foot on each step (foot-over-foot ascent) and to continue walking in straight line.
Based on the simulation results in the previous section, the exerciser needs 12 iterations to reach his desired setpoints. Practically, the experiments results for a new male subject (Subject No.10, Age=31 year) show that the participant has almost reached his desired setpoints after three iterations only.
Control system parameters and their corresponding experimental results for subject no.10 after three iterations
Duty cycle (%)
Referring to the results from Table 7 and Figure 10, we can conclude that the exerciser has almost reached his desired setpoints after only 3 iterations, when the simulation results showed that it can be done after 12 iterations. This shows how our controller has accomplished its duty with minimum number of iterations.
It is worth mentioning that for a proposed training regime based on an averaged model; we expect each subject to respond differently to a particular exercise or training protocol, of interest was Subject 10 who reached his desired set points after only three iterations. We may reach the desired set points or a steady state for another subject after 10, 11 or even 12 iterations, it depends totally on the intra- and inter- model uncertainties.
During the study, it was considered that after a certain number of training sessions, the training capacity of the subject will improve and this may affect the result of his HR response. In the next step of the study, we will recruit more subjects in order to verify our results and to even improve the performance of our controller.
A designed square-wave exercise protocol was applied to investigate the dynamic characteristics of HR and VO2 responses for both the onset and offset exercises. The K4b2 portable device was used to measure breath-by-breath VO2 and beat-by-beat HR. According to the experimental results, it was concluded that the time constants and steady state gains of the onset exercise for both HR and VO2 are distinctively different with those of the offset exercise. Based on the identified dynamic characteristics for both HR and VO2, we proposed an interval training protocol to improve cardiovascular fitness. A switching RC circuit model was presented to simulate HR and VO2 responses for the proposed interval training protocol. In order to adapt to individual users, a multi-loop integral control scheme has been proposed to regulate the HR response for interval running exercise under free living conditions. In the next step of the study, we will recruit more subjects for the validation of the proposed models, interval training protocol and the feedback control based adaption scheme.
The authors are thankful for the supports from the Centre for Health Technologies, the University of Technology, Sydney (UTS), Australia, and the school of human movement studies, the Charles Sturt University (CSU), Australia. Special thanks go to Dr. Cheyne Donges (the school of human movement studies, CSU), Dr. Rob Duffield, and Professor Aaron Coutts (the Centre for Health Technologies, UTS) for their help and guidance.
- Acharya R, Kumar A, Bhat IP, et al.: Classification of cardiac abnormalities using heart rate signals. Med Boil Eng Comput 2004,42(3):288–93. 10.1007/BF02344702View ArticleGoogle Scholar
- Achten J, Jeukendrup AE: Heart rate monitoring applications and limitations. Sports Med 2003 33(7):517–38.Google Scholar
- Ekelund U, Poortvliet E, Yngve A, Hurtig-Wennlöv A, Nilsson A, Sjöström M: Heart rate as an indicator of the intensity of physical activity in human adolescents. Eur J Appl Physiol 2001, 85: 244–249. 10.1007/s004210100436View ArticleGoogle Scholar
- Billat LV: Interval training for performance: A scientific and empirical practice. Sports Med 2001,31(1):13–19. 10.2165/00007256-200131010-00002View ArticleGoogle Scholar
- Per-Olof A, Rohdahl K: Work Physiology. New York: McGraw-Hill Book Company; 1970.Google Scholar
- Duffield R, Edge J, Bishop D: Effects of high-intensity interval training on the VO 2 response during severe exercise. J Sci Med Sport 2006, 9: 249–55. 10.1016/j.jsams.2006.03.014View ArticleGoogle Scholar
- Helgerud J, HØydal K, Wang E, Karlsen T, Berg P, Bjerkaas M, Simonsen T, Helgesen C, Hjorth N, Bach R, Hoff J: Aerobic high-intensity intervals improve VO 2 more than moderate training. Medicine and science in sports and exercise. Med Sci Sports Exerc 2007,39(4):665–71.View ArticleGoogle Scholar
- Gibala MJ, et al.: Short-term sprint interval versus traditional endurance training: similar initial adaptations in human skeletal muscle and exercise performance. J Physiol 2006,575(3):901–911. 10.1113/jphysiol.2006.112094View ArticleGoogle Scholar
- Talanian JL, et al.: Two weeks of high-intensity aerobic interval training increases the capacity for fat oxidation during exercise in women. J Appl Physiol 2007, 102: 1439–47.View ArticleGoogle Scholar
- Rakobowchuk M, et al.: Sprint interval and traditional endurance training induce similar improvements in peripheral arterial stiffness and flow-mediated dilation in health humans. Am J Physiol-Regulatory Integr Comp Physiol 2008,295(1):R236–42. 10.1152/ajpregu.00069.2008View ArticleGoogle Scholar
- TjØnna AE, et al.: Aerobic interval training reduces cardiovascular risk factors more than a multi-treatment approach in overweight adolescents. Clin Sci 2009,116(4):317–26. 10.1042/CS20080249View ArticleGoogle Scholar
- Talanian JL, et al.: Exercise training increases sarcolemmal and mitochondrial fatty acid transport proteins in human skeletal muscle. Am J Physiol Endocrinol Metab 2010,299(2):E180–88.Google Scholar
- Suh M-k, Rofouei M, Nahapetian A, Kaiser WJ, Sarrafzadeh M: Optimizing interval training protocols using data mining decision trees. Wearable and Implantable Body Sensor Networks 2009, BSN 2009: 318–323.View ArticleGoogle Scholar
- Haskell WL, Montoye HJ, Orenstein D: Physical activity and exercise to achieve health-related physical fitness components. Public Health Rep 1985,100(2):202–212.Google Scholar
- Koga S, Shiojiri T, Shibasaki M, Fukuba Y, Fukuoka Y, Kondo N: Kinetics of oxygen uptake and cardiac output at onset of arm exercise. Respir Physiol 1996,103(2):195–202. 10.1016/0034-5687(95)00082-8View ArticleGoogle Scholar
- Su SW, Chen W, Liu D, Fang Y, Kuang W, Yu X, Guo T, Celler BG, Nguyen HT: Dynamic modeling of heart rate response under different exercise intensity. Open Med Inform J 2010, 4: 81–85.View ArticleGoogle Scholar
- Bouchard DR, Trudeau F: Reliability of the assessment of the oxygen/heart rate relationship during a workday. Appl Ergon 2007,38(5):491–497. 10.1016/j.apergo.2006.10.002View ArticleGoogle Scholar
- Mc.Crory MA, Mole PA, Nommsen-Rivers LA, Dewey KG: Between-day and within-day variability in the relation between heart rate and oxygen consumption: effect on the estimation of energy expenditure by heart-rate monitoring. Am J Clin Nutr 1997, 66: 18–25.Google Scholar
- Duffield R, Dawson B, Pinnington HC, Wong P: Accuracy and reliability of a cosmed K4b 2 portable gas analysis system. J Sci Med Sport 2004,7(1):11–22. 10.1016/S1440-2440(04)80039-2View ArticleGoogle Scholar
- Balderrama C, Ibarra G, De La Riva J, Lopez S: Evaluation of three methodologies to estimate the VO 2max in people of different ages. Appl Ergon 2010, 42: 162–68. 10.1016/j.apergo.2010.06.017View ArticleGoogle Scholar
- Su SW, Wang L, Celler BG, Savkin AV, Guo Y: Identification and control for heart rate regulation during treadmill exercise. IEEE Trans Biomed Eng 2007,54(7):1238–46.View ArticleGoogle Scholar
- Su SW, Huang S, Wang L, Celler BG, Savkin AV, Guo Y, Cheng T: Nonparametric Hammerstein model based model predictive control for heart rate regulation. Conf Proc IEEE Eng Med Biol Soc 2007, 2007: 2984–7.Google Scholar
- Stanula A, Gabrys T, Szmatlan-Gabrys U, Roczniok R, Maszczyk A, Pietraszewski P: Calculating lactate anaerobic thresholds in sports involving different endurance preparation. J Exerc Sci Fit 2013,11(1):12–18. 10.1016/j.jesf.2012.12.001View ArticleGoogle Scholar
- Robergs RA, Landwehr R: The surprising history of the ” HRmax=220-age ” equation. J Exerc Physiol Online 2002,5(2):1–10.Google Scholar
- Inbar O, Oten A, Scheinowitz M, Rotstein A, Dlin R, Casaburi R: Normal cardiopulmonary responses during incremental exercise in 20–70-yr-old men. Med Sci Sport Exerc 1994,26(5):538–546.Google Scholar
- Goebel R, Sanfelice R, Teel A: Hybrid dynamical systems. IEEE Control Syst Mag 2009,29(2):28–93.MathSciNetView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.