Open Access

Modelling and regulating of cardio-respiratory response for the enhancement of interval training

  • Azzam Haddad1,
  • Yi Zhang1,
  • Steven Su1Email author,
  • Branko Celler2 and
  • Hung Nguyen1
BioMedical Engineering OnLine201413:9

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.


Interval training Oxygen uptake Heart rate Multi-loop integral controller


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 [3].

One of the most effective training regimens is called interval training which was first described by Reindell and Roskamm [4] and was popularized in the 1950s by the Olympic champion, Emil Zatopek. Swedish physiologist Per OlØf Astrand provided the first scientific study [5] 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 [10], and decrease the risk of cardiovascular-disease in obese adolescents [11]. Additionally, interval training increases proteins that transport fatty acids across the mitochondrial membrane [12]. 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 [13]. 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 [5]. In addition, it is well known that exercise at varying intensities helps improve exerciser’s aerobic capacity to exercise longer [14]. 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 [17]. 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


Eight healthy non-smoking males were invited to join the experiments. They were free from any known cardiac or metabolic disorders, hypertension, and were not under any medication. The UTS Human Research Ethics Committee (UTS HREC 2009000227) approved the study and an informed consent was obtained from all participants before each experiment. The physical characteristics of the participants are presented in Table 1.
Table 1

Physical characteristics


Age (yr)

Height (cm)

Weight (kg)









































SD: Standard Deviation.

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 [18], 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.

Experimental settings

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 [2].

The K4b2 gas analyzer was used because it has previously been reported to be valid, accurate and reliable [19].

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 [20].

The relationship between HR and VO2 becomes non-linear during light and very highly intense activity [2]. 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

Experimental results

It is known that many internal and environmental factors may have drastic effects on the readings during the experiment [2], 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.

Original signals of HR and VO2 for all subjects are shown in Figures 1 and 2 respectively.
Figure 1

Measured HR signals for 8 male subjects running on a treadmill under a predefined protocol.

Figure 2

Measured VO 2 signals for 8 male subjects running on a treadmill under a predefined protocol.

According to literatures [21, 22], heart rate response to exercise can be approximated as a first-order process. Since the response of VO2 has similar shape with HR, we model both HR and VO2 as first-order linear systems for both the onset and offset exercises. The first order system model can be written as follows:
G ( S ) = k τS + 1

Where k and τ are the steady state gain and the time constant respectively.

Matlab System Identification Toolbox has been used to identify the parameters of the model for each step response data. Time constants and steady state gains of HR and VO2 for each subject are shown in Tables 2 and 3.
Table 2

Estimated time constants and normalized steady state gains of HR response


Onset Exercise

Offset Exercise


SS Gain (k1)

Time Constant (τ 1)(sec)

SS Gain (k2)

Time Constant (τ 2)(sec)









































SS: Steady State.

Table 3

Estimated time constants and normalized steady state gains of VO 2 response


Onset exercise

Offset exercise


SS Gain (k1)

Time constant (τ 1)(sec)

SS Gain (k2)

Time constant (τ 2)(sec)









































SS: Steady State

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.

Average values of onset and offset time constants as well as steady state gains for VO2 and HR profiles are shown in Table 4. These values will be used to simulate the VO2 and HR responses during the proposed interval training protocol.
Table 4

Averaged values of time constants and steady state gains





τ 1


τ 2






















SD: Standard Deviation

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 [16] 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.

On the other hand, training in anaerobic zone which is 80% to 90% of HRmax will develop the lactic acid system [23]. Therefore, determining an actual HRmax is the key to constructing a well-designed training program. According to [24], a population specific formula should be used for better estimation of HRmax for an individual subject. However, the most accurate general equation to estimate HRmax is that of Inbar [25]:
HR max = 205.8 0.685 × ( Age )
It is well known that HR is a valid measure of exercise intensity only if it reflects the metabolic rate which can be measured by VO2. Improving aerobic capacity can be achieved by increasing oxygen delivered to the muscles during exercise. Table 5 shows the relationship between VO2max and HRmax [3].
Table 5

Relationship between HR max and VO 2max











Table 6 shows the training zones and their corresponding value for both HRmax and VO2max.
Table 6

Exercise intensity levels that coincide with HR max and VO 2max




Fat burning (low)



Aerobic (moderate)



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.

Figure 3 shows an interval training protocol based on the above calculations.
Figure 3

Proposed interval training protocol targeting an aerobic zone for an averaged model.

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.

Before verifying this protocol by experiment, we tested it by simulation first. Switching RC circuit as shown in Figure 4 was used to simulate both HR and VO2 responses.
Figure 4

Designed RC circuit that can be used in simulating HR and VO 2 responses to interval training protocol.

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.

New male subject (Subject No.9, Age = 30 year, Weight = 84 Kg, and Height = 185 cm) was asked to run on the treadmill under the new interval training protocol. Experimental results for subject 9 are shown in Figure 5.
Figure 5

HR and VO 2 experimental results for subject No.9 under 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 [2]. 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

The final goal of the adaptation of the proposed interval training protocol is tuning the duty cycle (Δ t1) as shown in Figure 6, and the period (Δ t1 + Δ t2) to achieve desired training effects. The adaptation process is a set of successive interval training experiments for a particular user. This process can be described by using a hybrid system model [26].
Figure 6

Controller outputs’ inputs (duty cycle ( Δ t 1 ) and period( Δ t 1 + Δ t 2 )) and HR response during interval training protocol.

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.

Based on the previous model, the outputs y(t4) and y(t5) can be obtained as follows:
y HR ( t 4 ) = y HR ( t 3 ) . ( 1 ( k 2 / k 1 ) ( 1 e Δ t 2 T 2 ) ) y HR ( t 5 ) = y HR ( t 1 ) + y HR ( t 4 ) . e Δ t 1 T 1
Where k1, k2, T1 and T2 are the onset, offset gains and time constants respectively, and
y HR ( t 1 ) = k 1 . ( 1 e Δ t 1 T 1 ) y HR ( t 2 ) = y HR ( t 1 ) k 2 . ( 1 e Δ t 1 T 1 ) . ( 1 e Δ t 2 T 2 ) y HR ( t 3 ) = y HR ( t 1 ) + y HR ( t 2 ) . e Δ t 1 T 1

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.

Figure 7 is the block diagram of the proposed control system. Two discrete PI controllers C1 and C2 have been built and connected to the HR model to control the desired output values y(t4) and y(t5) separately. The controller parameters (the coefficients of P and I actions) of these two simple multi PI controllers are tuned by trial and error. We simulated the control performance in Matlab Simulink.
Figure 7

Multi-loop PI control system to regulate output values y(t 4 ) and y(t 5 ).

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.

The self-adaption feature gives the exerciser the opportunity to reach his desired setpoints after a number of iterations. Figure 8 below shows the simulated controller outputs, y(t4) and y(t5).
Figure 8

Simulation results of PI control system outputs. Offset output signal y(t4) of HR response (a) and onset output signal y(t5) of HR response (b).

In simulation, from Figure 8, we conclude that after at least 12 iterations (12 training exercises) the system will reach the desired setpoints.

In other words, t4 and t5 have to be observed 12 times and re-entered to the controller to achieve the desired outputs. Figure 9 below shows the HR response for the 13th iteration. The third waveform in Figure 9 shows that the output signal finally reached the desired reference values R1 and R2 at t4 and t5 respectively.
Figure 9

Continuous interval training HR response simulation from the controller.

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.

Based on the exerciser’s age, the controller R1 and R2 are set to be 148 bpm(80% of HRmax) and 111 bpm(60% of HRmax) respectively. In the first iteration, 60 seconds is set to be the onset and offset times. At the end of the first training session, the reading of HR at t4 and t5 were observed and feed backed to the multi-loop control system to update the controller output’s inputs (the duty cycle and the period). Accordingly, the controller adjusts the time of onset and offset periods for the next training session. Table 7 shows the controller parameters and their corresponding experimental results for Subject No.10 in three training sessions.
Table 7

Control system parameters and their corresponding experimental results for subject no.10 after three iterations


Onset (sec)

Offset (sec)

Period (sec)

Duty cycle (%)































Figure 10 shows the new subject’s HR response under the stairs climbing activity.
Figure 10

Subject 10’s HR responses under stairs climbing exercise. HR response after the first iteration (a), HR response after the second iteration (b) and HR response after the third iteration (c).

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.

Authors’ Affiliations

Faculty of Engineering and IT, University of Technology, Sydney (UTS)


  1. 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
  2. Achten J, Jeukendrup AE: Heart rate monitoring applications and limitations. Sports Med 2003 33(7):517–38.Google Scholar
  3. 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
  4. 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
  5. Per-Olof A, Rohdahl K: Work Physiology. New York: McGraw-Hill Book Company; 1970.Google Scholar
  6. 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
  7. 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
  8. 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
  9. 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
  10. 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
  11. 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
  12. 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
  13. 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
  14. 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
  15. 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
  16. 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
  17. 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
  18. 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
  19. 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
  20. 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
  21. 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
  22. 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
  23. 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
  24. Robergs RA, Landwehr R: The surprising history of the ” HRmax=220-age ” equation. J Exerc Physiol Online 2002,5(2):1–10.Google Scholar
  25. 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
  26. Goebel R, Sanfelice R, Teel A: Hybrid dynamical systems. IEEE Control Syst Mag 2009,29(2):28–93.MathSciNetView ArticleGoogle Scholar


© Haddad et al.; licensee BioMed Central Ltd. 2014

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 (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.