Research Collection
Master Thesis
Modelling and control of the human cardiovascular system
Author(s): Gisler, Stefan
Publication Date: 2011
Permanent Link: https://doi.org/10.3929/ethz-a-007207574
Rights / License: In Copyright - Non-Commercial Use Permitted
This page was generated automatically upon download from the ETH Zurich Research Collection. For more information please consult the Terms of use.
ETH Library Master Thesis
Modelling and control of the human cardiovascular system
Stefan Gisler
Advisers Martin Wieser and Dr. Heike Vallery and Prof. Dr. Robert Riener Sensory Motor Systems Lab (SMS) Swiss Federal Institute of Technology Zurich (ETH)
Submission: April 2011 Contents
1 Introduction 1
2 Human cardiovascular system 5 2.1 Hemodynamic system ...... 5 2.2 Blood pressure regulation ...... 7 2.3 Orthostatic reaction and muscle pump ...... 9 2.4 Cardiovascular pathology ...... 10 2.5 Literature review ...... 11 2.5.1 Cardiovascular responses to passive tilting ...... 11 2.5.2 Cardiovascular modelling ...... 13
3 Cardiovascular model 15 3.1 Hemodynamic system ...... 16 3.2 Blood pressure regulation ...... 18 3.3 Influence of gravity ...... 21 3.4 Influence of stepping ...... 22 3.5 Model simulations ...... 24 3.5.1 Fast tilt-up and tilt-down ...... 24 3.5.2 Stepping ...... 25 3.5.3 Quasi-static ...... 26 3.6 Model validation ...... 33
4 Control design 37 4.1 Model predictive control (MPC) design ...... 37 4.2 Simulation ...... 46
5 Methods 49 5.1 Healthy subjects ...... 49 5.1.1 Implementation ...... 49 5.1.2 Blood pressure recording ...... 49 5.1.3 Experimental design ...... 50
i 5.2 Patients ...... 53 5.2.1 Implementation ...... 53 5.2.2 Blood pressure recording ...... 53 5.2.3 Experimental design ...... 54
6 Results 55 6.1 Healthy subjects ...... 55 6.1.1 Heart rate control ...... 55 6.1.2 Blood pressure control ...... 55 6.1.3 Combined heart rate and blood pressure control . . . . 56 6.2 Patients ...... 57 6.3 Controller performance ...... 60
7 Discussion 61
8 Conclusion and Outlook 65
A Model summary 73 A.1 List of variables ...... 75 A.2 List of parameters ...... 76 A.3 Model equations in non-linear state-space form ...... 78 A.4 Steady-state equations in non-linear state-space form . . . . . 80 A.5 Parameter identification ...... 82 A.6 Model constraints ...... 85
B Summarised results 87
References 95
ii Abstract
Bed-rest leads to cardiovascular deconditioning and may induce a decline in stroke volume, cardiac output and oxygen uptake. Further, it increases the risk of orthostatic intolerance. In an early phase of rehabilitation, it is there- fore important to prevent the development of cardiovascular deconditioning which can be done by verticalisation and mobilisation. In the future, the enhanced ERIGO tilt-table will be able to control physiological signals and hence, stabilise the patient’s cardiovascular system. This thesis focuses on the control of heart rate and blood pressure by means of verticalisation (tilting) and mobilisation (stepping). In a first step, a car- diovascular non-linear model with two inputs (tilting and stepping) and three outputs (heart rate, systolic and diastolic blood pressure) is developed based on physiological principles and existing work. The model is then used for designing a model predictive controller which was found well suited for the given control problem. Five healthy subjects have been tested with three different configurations: isolated heart rate control, isolated blood pressure control and combined control. One patient has been tested with blood pressure control which yiel- ded promising results.
Keywords– Orthostatic intolerance, cardiovascular modelling, model pre- dictive control
iii Acknowledgements
First, I want to thank Prof. Dr. Riener for being accepted to do this thesis at the Sensory-Motor Systems Lab. Then I want to thank my advisers Martin Wieser and Dr. Heike Vallery for their valuable support during the work. Special thanks go to Martin Wieser for his great efforts while testing and debugging the system. This thesis would not have been possible without the probands and patients. A big thanks goes to all the probands, the “Z¨urcher H¨ohenklinik”in Wald, and all the patients that participated in this study. At this point, I also want to thank Rafael R¨ustand Lilith B¨utlerfor their support during the patient measurements in Wald. Last but not least, I want to thank all the students in the student room for the nice and inspiring atmosphere.
iv Chapter 1
Introduction
One major problem with neurological patients suffering from stroke, trau- matic brain injury or paraplegia is the long bed rest after the accident. It leads to deconditioning of the patients’ cardiovascular system and evokes secondary complications such as orthostatic intolerance. Further complica- tions can include venous thrombosis, muscle atrophy, joint contractures and osteoporosis [1], [2]. Therefore, early mobilisation of the patient is crucial as it can reduce the risk of cardiovascular deconditiong and improves the state of health.
This thesis focuses on the cardiovascular aspects of bed-ridden patients, i.e. how the cardiovascular system can be prevented from deconditioning and be- coming unstable. Prolonged bed rest leads to a decrease in circulating blood volume, a decrease in stroke volume and pulse pressure, and an increased heart rate. A direct result of these indications is the inability of the patient’s cardiovascular system to regulate blood pressure when standing up (ortho- static intolerance). In the upright position, the patient suddenly starts to feel dizzy or even faints due to excessive blood pooling in the lower extre- mities and reduced blood perfusion of the upper body. However, orthostatic intolerance is not only caused by prolonged bedrest but can also be a conse- quence of an impaired vegetative nervous system. In paraplegia patients, the sympathetic effector nerves to the heart and the smooth musculature are disrupted or even broken. This leads to a malfunction of the baroreflex which is responsible for regulating arterial blood pressure (see chapters 2.2, 2.4). As a consequence, the sudden decrease in arterial blood pressure cannot be regulated and the patient faints.
A tilt-table therapy is aimed at reconditioning the patient’s cardiovascular system by verticalising to an angle of about 80 degrees. Additional leg mo-
1 vements which can include stepping or cycling movements increase venous return due to the effects of the muscle pump and improve orthostatic tole- rance. The ERIGO device which has been used at the institute since the beginning of the AwaCon project combines these therapies and allows for an optimal treatment of patients with neurological disorders (Figure 1.1). More information about the ERIGO device can be found on the homepage of HO- COMA AG 1. On the ERIGO, physiological signals such as blood pressure,
Figure 1.1: Left: Schematic representation of the ERIGO device with the three inputs. Right: ERIGO during therapy session. heart rate, respiration frequency, skin conductance, oxygen saturation, EEG and EMG can be recorded. However, for this thesis only blood pressure and heart rate need to be recorded, where EMG recordings may be helpful to analyse muscle activity during mobilisation.
The goal of the project is to control and stabilise the cardiovascular system of patients with neurological disorders by verticalisation, mobilisation and cyclic loading of the lower limbs (Figure 1.1). This will help to improve the cardiovascular status of these patients and will have the potential to reduce medication, enhance physiotherapy and shorten the duration of early reha- bilitation [3]. Furthermore, the risk of deconditioning of the cardiovascular system, and complications resulting from this, can be decreased. Additional project information is available on the homepage of the SMS Lab 2. In earlier projects at the SMS, isolated control of heart rate and diastolic
1http://www.hocoma.com/en/products/erigo/ 2http://www.sms.mavt.ethz.ch/research/projects/awacon
2 blood pressure with the inclination angle α as the only control input has been done [4], [5]. In a next step, combined control of heart rate and dias- tolic blood pressure has been succesfully tested with healthy subjects [6]. This latest version also contained another technical innovation: the idea was to not only use α as a control input, but also the stepping frequency fstep which enables the controller to operate over an enlarged bandwith. For this project, the described line of innovation is continued: the goal of this thesis is to control heart rate, systolic and diastolic blood pressure with the two control inputs α and fstep. It is a fact, that the control strategy from the isolated control problem, which consisted of an ordinary PI controller can not be adopted for the new more complex control problem. The challenge is that with an increasing number of inputs and outputs, there are more couplings inside the system and PI control is not suitable anymore. For a multi-input multi-output (MIMO) system, other control strategies have to be applied. The first step consists of developing a cardiovascular model which is the topic of chapter 3 which directly follows after the subsequent chapter about human cardiovascular physiology (chapter 2). Chapter 4 continues with the control design, followed by the results, the discussion and the conclusion (chapters 6, 7 and 8).
3 4 Chapter 2
Human cardiovascular system
This chapter will give a short introduction to physiology and pathophysio- logy of the human cardiovascular system and summarises some results from literature.
2.1 Hemodynamic system
Head
Right lung Left lung
Right Left heart heart
Splanchnic & renal circulation
Legs
Figure 2.1: Schematic representation of the human circulatory system. Adap- ted from: http://www.ionwave.ca
The major task of the hemodynamic system is to supply every single cell
5 of the organsim with oxygen and nutrients and carry away carbon dioxide (CO2) as well as metabolic waste products. In the circulation, the heart acts as a pump which produces a pressure gradient between arterial and venous circulation. Driven by this pressure gradient, deoxygenated blood from the venous circulation flows back to the right heart where it is pumped through the lung. In the lung the blood is enriched with oxygen and reenters systemic circulation when pumped into the aorta by the left ventricle. The arterial tree then supplies the whole body with oxygen and nutrients. From the peripheral regions, where the oxygen and the nutrients are used, the blood returns to the right heart and the circulation is closed (Figure 2.1). Flows and pressures within the human hemodynamic system are characteri- sed by the following list of hemodynamic variables:
• Stroke volume (SV ) defines the amount of blood pumped into the aorta within one beat.
• Cardiac output (CO) is calculated as the product of stroke volume and heart rate (HR) CO = SV · HR
• Systolic blood pressure (sBP ) is the maximal blood pressure that oc- curs during the contracting heart phase (systole).
• Diastolic blood pressure (dBP ) is the minimal blood pressure that occurs during the filling period of the heart, when the ventricles are relaxed (diastole).
• Mean arterial pressure (MAP ) is defined as the integrated blood pres- sure over one heart period divided by the time of one heart period.
R t+tRR BP (t)dt MAP = t tRR
where BP (t) is the continuous blood pressure and tRR is the time of one heart period (R-R interval). A common approximation is given as
1 2 MAP = · sBP + · dBP 3 3
• Central venous pressure (CVP ) is the pressure in the intrathoracic veins and the right atrium. Normal values range from 2 to 4 mmHG [7].
6 • Total peripheral resistance (TPR) is a rather hypothetic measure of vessel resistance in the systemic circulation. In duality to Ohm’s law U = R · I, total peripheral resistance is defined as MAP − CVP TPR = CO
MAP CVP is usually neglected in this calculation and we get TPR = CO .
2.2 Blood pressure regulation
Regulation mechanisms in the cardiovascular system are responsible for adap- ting the hemodynamic variables such as blood pressure and cardiac output according to the body needs. During exercise for example, cardiac output is strongly increased in order to cover the high oxygen need in the skeletal muscles. Another situation where these regulation mechanisms are active is when the body adapts to changes in environmental conditions such as tem- perature differences. And last, these mechanisms are also active in response to orthostatic stress what will be of interest for blood pressure and heart rate control on the ERIGO. Hemodynamic variables can be influenced in several ways. Natural control mechanisms include neurogenic (over the vegetative nervous sytem), hormo- nal (over circulating hormones), humoral (with locally formed substances) or myogenic regulation (vasoconstriction with smooth musculature). For short-term regulation the neurogenic mechanisms which include baroreflex, cardiopulmonary reflex and chemoreceptor reflex are most important. These three types of neurogenic blood pressure regulation will now be described in more detail: The baroreflex plays a central role in short-term blood pressure regulation. The baroreceptors which are located in the aortic arch and the carotid sinus are the sensors in this reflex mechanism. They transmit neural signals to the central nervous system or more precisely to the cardiovascular centre in the medulla oblongata. The impulse frequency of the afferent neurons is determi- ned by the course of the arterial blood pressure: Low arterial blood pressure leads to a high impulse frequency. However, impulse frequency is not only determined by absolute value of the arterial blood pressure but also by its time rate of change. This proportional-derivative (PD) sensor characteristics enable the baroreceptors to send all relevant information about heart func- tion to the central nervous system. In the medulla oblongata the information from the baroreceptors is transmitted to the efferent vegetative nervous sys- tem which determines heart rate, heart contractility and vasoconstriction of
7 peripheral blood vessels, closing the reflex arch. It has to be added that an inhibitory interneuron in the medulla provokes negative feedback which is essential for regulating and stabilising arterial blood pressure.
Head Baroreceptors Cardiovascular centre
Right lung Left lung
Regulate heart rate & Baro- cardiac contractility reflex Right heart Left heart
Splanchnic & renal circulation
Regulate peripheral resistance
Legs
Figure 2.2: Blood pressure regulation with the baroreflex loop.
The cardiopulmonary reflex is another blood pressure regulating mechanism that works synergistically with the baroreflex. The cardiopulmonary recep- tors are located in the venous system, more precisely in the atria and A. pulmonalis. However, cardiopulmonary receptors are not only responsible for blood pressure regulation but also for volume regulation. Stimulation of the receptors by dilated atria leads to an inhibited production of the anti- diuretic hormone (ADH). As a consequence, urine secretion is increased and the circulating blood volume can be reduced. Furthermore, activation of the cardiopulmonary receptors decreases sympathetic activity and inhibits Renin production in the kidneys. Renin promotes the formation of Angiotensin II which has a direct vasoconstrictive effect on the smooth musculature in the vessels. Moreover, Angiotensin II stimulates the production of Aldosterone in the kidneys which increases reabsorption of sodium and water. In the long term, this leads to a higher blood volume and an increased blood pressure. Hence, the Renin-Angiotensin-Aldosterone system (RAAS) is capable of in- creasing arterial blood pressure by the vasoconstrictive effect of Angiotensin II and the volume retention caused by Aldosterone. Note that volume re- gulation is a long-term regulation because it includes hormonal mechanisms
8 and because it takes some time until body fluids have diffused through the capillary walls. The chemoreceptor reflex is mainly responsible for respiration control, but can also influence cardiovascular regulation if partial pressure of oxygen in the blood decreases [7]. The reflex mechanism particularly becomes active if blood pressure falls below 80 mmHg and once active, it acts in the same feedback structure as the baroreflex. As a result, arterial blood pressure is increased. In detail, myogenic regulation is also a kind of neural regulation mechanism if we consider the sympathetic effected vasoconstriction in the peripheral ar- terioles. However, there is also a mechanism called autoregulation that is attributed to myogenic regulation. Autoregulation is the ability of a blood vessel to keep the blood flow constant under changing perfusion pressures. When perfusion pressures are increased, the smooth musculature is activa- ted and prohibits further expansion of the vessel walls (myogenic reaction: Bayliss effect).
2.3 Orthostatic reaction and muscle pump
Everybody knows the dizzy feeling after standing up too fast in the morning. The body’s internal regulation mechanisms are strongly challenged in such si- tuations. Normally, the neural regulation mechanisms as discussed above are able to maintain homeostasis quite fast. Nevertheless, there might be situa- tions where the regulation is incapable of keeping arterial blood pressure and cerebral perfusion at a safe level. Low blood volume or high temperatures for example are conditions that increase the risk of defective homeostasis. This can lead to a syncope which can be rather dangerous when the fainting person falls down on the floor or hits a hard object. One cause of such a syncope is the venous blood pooling in the legs. In a heal- thy person, up to half a litre of blood is shifted from the upper body to the lower extremities [7]. Arterial blood pressure falls immediately and the reflex mechanisms are activated. However, peripheral vasoconstriction caused by sympathetic regulation is usually too weak in order to lower venous blood pooling effectively. Fortunately, there is another mechanism besides the neu- ral regulation which is capable of stabilising the cardiovascular system. The principle behind this mechansim is that the contraction of the skeletal leg muscles efficiently compresses the venous compartments, decreases venous pooling and increases venous return to the heart. Because the venous valves are closed, backflow is not possible and the blood is forced to return back to the heart (figure 2.3). This “muscle pump” is always active when the skeletal
9 leg musculature is active, for example during walking.
Figure 2.3: The muscle pump mechanism stabilises the cardiovascular system by efficiently reducing venous blood pooling and increasing venous return by repeated contractions of the skeletal leg musculature. Source: University of Minnesota
2.4 Cardiovascular pathology
Cardiovascular instability and orthostatic hypotension are common deficits in bed-ridden patients [1], [2]. In spinal cord injury (SCI) patients, for example, one reason for these deficits are the disrupted efferent sympathetic pathways regulating heart rate, heart contractility and peripheral vasocons- triction. Therefore, neural regulation mechanisms can not work properly and blood pressure often drops dramatically in reaction to orthostatic stress. The disturbed balance between sympathicus and parasympathicus leads to an exaggerated increase in heart rate as a compensatory reaction to the blood pressure decrease. This happens because parasympathetic heart rate regu- lation is still intact in SCI patients as efferent parasympathetic nerves are connected to the Vagus nerve and not to the spine. Naturally, sympathetic nervous system disfunction is not the only reason for orthostatic intolerance in neurological patients. Claydon et al. [8] summarise these factors for SCI patients as follows:
10 • Sympathetic nervous system disfunction
• Altered baroreceptor sensitivity
• Lack of skeletal muscle pump
• Cardiovascular deconditioning
• Altered salt and water balance
Baroreceptor sensitivity which is typically reduced in SCI patients is in tight connection with the sympathetic nervous system disfunction. As explained above, baroreflex regulation is severly damaged because of an impaired sym- pathetic nervous system. The lack of the skeletal muscle pump together with immobilisation and pro- longed bed-rest are the reason for cardiovascular deconditioning which in turn negatively affects the overall recovery. Lastly, Claydon et al. report evidence that SCI patients have a decreased plasma volume as a result of an impaired salt and water balance. This leads to problems in volume regula- tion, i.e. hypovolemia and low resting blood pressure with a predisposition to orthostatic intolerance.
2.5 Literature review
2.5.1 Cardiovascular responses to passive tilting Passive tilting leads to an immediate increase of blood volume in the leg veins of about half a liter [7]. Venous return is decreased and because of the Frank-Starling mechanism stroke volume and pulse pressure are diminished as well. To counter the blood pressure drop, neural reflexes are instantly activated and sympathetic action is increased. This has two consequences: Firstly, heart rate rises by approximately 20 % [7] and secondly, diastolic blood pressure rises because of increased peripheral resistance. In contrast, systolic blood pressure is normally rather constant [7], [9]. The above description is considered the healthy cardiovascular response to passive tilting according to standard physiological work of reference such as [7]. Table 2.1 lists the outcome of several studies about cardiovascular responses to passive tilting involving healthy subjects. Note that most of these experimental results conform with the “standard” physiological res- ponse. As the aim of the thesis and the whole project is to enhance therapy of neu- rological patients, a quick survey of typical pathophysiological cardiovascular
11 year HR sBP dBP MAP Hainsworth and Al-Shamma [10] 1988 ↑ ≈ ↑ ↑ Mukai et al. [11] 1995 ↑ ≈ ≈ ≈ Tanaka et al. [12] 1996 ↑ ↓ ≈ ↓ Cooke et al. [13] 1999 ↑ ↑ ↑ ↑ Yokoi and Aoki [14] 1999 ↑ ≈ ↑ ↑ Petersen et al. [15] 2000 ↑ ≈ ↑ ↑ Tulppo et al. [16] 2000 ↑ ≈ ↑ ↑ Toska and Walloe [17] 2002 ↑ n/a n/a ↑ Heldt et al. [18, 19] 2003/04 ↑ ≈ ↑ ↑ Colombo et al. [20] 2005 n/a ↑ ↑ ↑ Masuki et al. [21, 22] 2007 ↑ ↑ ↑ ↑ Chi et al. [23] 2008 ↑ ↑ ↑ ↑ Ramirez et al. [24] 2008 ↑ ≈ ↑ ↑
Table 2.1: Literature summary. ≈ means no significant change, ↑ means significant increase, ↓ means significant decrease. (Adapted and completed with HR from [5]) responses will be done. Table 2.2 presents standard cardiovascular responses of SCI patients. All of these studies conform with the normal pathological reaction to orthostatic stress in SCI patients as described in section 2.4. In addition, on the basis of the work of Houtman [25] and Legramante [26] it can be stated that the higher the lesion the bigger are the implications on the cardiovascular system and the cardiovascular regulation.
year aetiology HR sBP dBP MAP Corbett et al. [27] 1971 Tetrapl. ↑ ↓ ↓ ↓ Houtman et al. [25] 2000 Normal ↑ n/a n/a ↑ Parapl. ↑ n/a n/a ↑ Tetrapl. ↑ n/a n/a ↓ Legramante et al. [26] 2001 Normal ↑ ≈ ≈ ≈ Parapl. ↑ ≈ ≈ ≈ Tetrapl. ↑ ↓ ↓ ↓
Table 2.2: Literature summary. ≈ means no significant change, ↑ means significant increase and ↓ means significant decrease
12 2.5.2 Cardiovascular modelling Computational models of the human cardiovascular system have been deve- loped for many different purposes. An elaborate cardiovascular model can be used to identify aetiologies of cardiovascular diseases such as orthostatic intolerance (OI). Heldt et al. [28] have presented a complex mathematical mo- del which reproduces cardiovascular responses to orthostatic stress. In their study the model was used to investigate the mechanisms that cause postspa- ceflight OI. Leaning et al. [29] formulated a detailed model intended to study and predict the overall effects of an injected drug. However, a cardiovascular model can also be used to examine specific aspects of the cardiovascular sys- tem such as blood-pressure fluctuations and heart-rate variability [30], [31]. Most of these models are aimed at explaining a certain cardiovascular phe- nomenon and are usually rather complex with a high model order. They are normally based on a large number of compartments representing the different parts of the circulation (heart chambers, ventricles, venous and arterial seg- ments). Each compartment or reservoir has a certain pressure Pj and volume Vj Vj − Vj0 Pj = (2.1) Cj where Cj is the compliance and Vj0 the unstressed or zero-pressure volume. Most models that describe the overall cardiovascular system incorporate some elements of nervous system regulation. The baroreflex plays an essential role because it governs the short-term dynamics of blood pressure and heart rate. Long-term dynamics are most often less important than short-term effects and can be neglected in the model description. Therefore, blood pressure regulation mechanisms such as RAAS do not need to be modelled. There are hardly any cardiovascular models in literature which incorporate an orthostatic component and are kept simple. One exception is in the work of Akkerman [32] who presented a mathematical beat-to-beat model designed for tilt-table experiments. He analysed the dynamics of cardiovascular signals after fast tilt-up and tilt-down. The model forms the basis of the whole controller design and will be explained in detail in the following section.
13 14 Chapter 3
Cardiovascular model1
In order to control physiological quantities such as blood pressure and heart rate, an appropriate model of the cardiovascular system is needed. This model should have two inputs, namely the inclination angle α of the ERIGO device and the stepping frequency fstep. Based on these inputs the model should output heart rate, systolic and diastolic blood pressure (Figure 3.1).
HR α = u = y P S fstep PD
Cardiovascular model
Figure 3.1: Inputs and outputs of the cardiovascular model
In order to develop a mathematical model of the human cardiovascular sys- tem for blood pressure regulation, the following two assumptions are made: • For the internal blood pressure regulation, only the baroreflex is taken into account. Other mechanisms such as the cardiopulmonary reflex and the RAAS system are not needed to explain the main blood pres- sure characteristics in tilt-up and tilt-down because they govern the long-term dynamics. 1The material in this chapter is closely related to Akkerman’s work [32]
15 • The blood volume is constant and fluid movements through the capil- lary walls are not considered.
So the model only contains the most important elements that are needed to simulate orthostatic reactions, namely a closed hemodynamic system, a kind of internal blood pressure regulation system and the influence of gravity and stepping (Figure 3.2). These parts will be explained in detail in the following sections.
3.1 Hemodynamic system
The hemodynamic system as explained in section 2.1 can be modelled as a connected system of pipes representing blood vessels. The heart acts as a pump, maintains systemic blood pressure and transports oxygen-poor blood to the lung. The microcirculation in the peripheral parts of the body is the bottleneck in the pipe system and is therefore also called the peripheral re- sistance. The fact that blood vessels are not stiff tubes but compliant vessels is accounted for by introducing a venous and an arterial reservoir which is common engineering practice. In fact, the flattening effect that arterial com- pliance has on the systolic blood pressure peaks is called the “Windkessel” effect which is in accordance with the above mentioned engineering principle of introducing reservoirs for the modelling of compliant tubes. In Akkerman’s model, only the lung, the arteries and the veins are modelled as proper compartments as defined by Equation 2.1. The volumes of these compartments are denoted by VP , VA and VV respectively. In addition, each of these compartments is attributed a compliance (CP , CA and CV ) and a zero-pressure volume (VP 0, VA0 and VV 0). According to Equation 2.1 the compartment pressures at heart beat k + 1 can then be expressed as:
VV (k) − VV 0 PR(k + 1) = (3.1) CV VP (k) − VP 0 PL(k + 1) = (3.2) CP VA(k) − VA0 PD(k + 1) = (3.3) CA where PR(k) is the right atrial pressure, PL(k) the left atrial pressure and PD(k) the diastolic blood pressure which directly depends on the arterial blood volume. The flow between these three reservoirs is characterised by the following set of equations where VPP (k) describes the volume in the pulmonary pipeline
16 Head
Cardiovascular Lung Baroreceptors centre B VP, C P
Windkessel
Right VA, C A P R heart QR
I Baro- P Q reflex L, L P, P Left heart S D
VV, P V , C V
Venous reservoir Q W R
Peripheral circulation
Figure 3.2: Simplified representation of the human cardiovascular system used for model synthesis. Adapted from [32]
which is needed to model the delay between right and left atrium. ξP denotes the number of right stroke volumes that are in the pulmonary pipeline.
VP (k) = VP (k − 1) + QR(k − ξP ) − QL(k) (3.4)
VA(k) = VA(k − 1) + QL(k) − QW (k) (3.5)
VV (k) = VV (k − 1) + QW (k) − QR(k) (3.6)
VPP (k) = VPP (k − 1) + QR(k) − QR(k − ξP ) (3.7)
Based on the Frank-Starling law and the restitution properties of ventricular myocardium, the left and right stroke volumes QL(k) and QR(k) depend on the preload and the length of the previous R-R interval I(k − 1). Akkerman
17 adapted these findings from [33]:
QR(k) = γRPR(k)I(k − 1) (3.8)
QL(k) = γLPL(k)I(k − 1) (3.9) where γR and γL are constant factors called “Starling” factors. The periphe- ral flow QW (k) depends on the peripheral resistance R(k) and the pressure difference between the arterial and the venous segment.
−I(k) QW (k) = CA (PS(k) − PV (k)) 1 − exp (3.10) R(k)CA where PV (k) denotes venous pressure evaluated just after systole when the right stroke volume has been ejected into the pulmonary pipeline:
VV (k − 1) − VV 0 − QR(k) PV (k) = (3.11) CV
The equations for pulse pressure PP (k) and systolic blood pressure PS(k) complete the hemodynamic system:
QL(k) PP (k) = (3.12) CA PS(k) = PD(k) + PP (k) (3.13)
All the introduced variables are beat-to-beat variables which means that they are updated at each heart beat. It is not clear, however, at which instant of the heart beat these variables are refreshed. The systolic blood pressure PS(k) for example is updated during the systole when the continuous blood pressure curve peaks at its maximum value. In contrast, the diastolic blood pressure PD(k) is updated at the end of the diastole. Each hemodynamic variable has its natural physiological sampling instant. Another example are the right and the left stroke volumes QR(k) and QL(k). These variables are updated at the beginning of the systole when the stroke volumes are ejected into the pulmonary pipeline and the aorta respectively. Figure 3.3 graphically summarises the sampling instants of the introduced hemodynamic variables.
3.2 Blood pressure regulation
Blood pressure regulation is performed by different body mechanisms. There are short-term regulations (minutes, hours) and long-term regulations (days,
18 Figure 3.3: Hemodynamic timetable describing at which moment of the heart beat each hemodynamic variable is evaluated. Source: Akkerman [32] weeks 2) as mentioned above. For the purpose of blood pressure and heart rate control, only the short-term regulations have to be considered. There- fore, the modeling will focus on the baroreflex mechanism. Katona et al. [34] have developed a baroreflex model which is composed of a sympathetic and a parasympathetic branch (Figure 3.4) which is widely used in computatio- nal modelling of the human cardiovascular system. Based on a hypothetic barosignal which is a function of arterial blood pressure and pulse pressure, the model outputs the heart period. The model is split in two parts because sympathetic and parasympathetic dynamics are rather different. Parasym- pathetic activity leads to a fast decrease of heart rate which can be shown by electrical stimulation of the Vagus nerve. In contrast, the sympathetic contri- bution on heart rate is slower. In Katona’s model there is a fixed boundary between sympathetic and parasympathetic regulation. Of course, in reality there is a smooth transition between these two types of blood pressure regu- lation. However, it is not needed to map this behaviour to the model and this intuitive simplification is very well applicable. It is even the case, that for applications where these subtle dynamics are of minor importance, Katona’s baroreflex model can be further simplified. In this thesis, the two branches are merged to one neglecting the different dynamics of sympathetic and pa- rasympathetic regulation. More important is the extension of the model by
2Time specifications: http://homepages.uel.ac.uk/M.S.Meah/bs250page4clec3. htm
19 Figure 3.4: Katona’s baroreflex model for heart rate regulation [34]. The neural input signal f(t) is divided in a sympathetic (bottom) and parasym- pathetic part (top) where µ defines the borderline between sympathetic and parasympathetic regulation.
a branch for regulation of the peripheral resistance as proposed by Akker- man [32]. This regulation is based on the sympathetic part of the hypothetic neural barosignal and the implementation is straightforward (Figure 3.5). Put in equations, the baroreflex model can be stated as follows:
B scaling PT1 + I β Bc +
out scaling BS PT1 - ρ + R in B + Bc c RP +
Figure 3.5: Simplified baroreflex model (based on Katona [34] and Akker- man [32]): the two branches regulating heart rate have been merged to one, a second branch has been added for regulation of peripheral resistance.
20 B(k) = PB(k) + kP PP (k) − kB (3.14) 1 2 = P (k) + P (k) − σ sin α(k) + k P (k) − k (3.15) 3 S 3 D B P P B BS(k) = min(B(k),Bc) (3.16) −1 −1 τ τ JBI (k) = e BI JBI (k − 1) + 1 − e BI B(k) (3.17) −1 −1 τ τ JBR(k) = e BR JBR(k − 1) + 1 − e BR BS(k) (3.18)
I(k) = (JBI (k) + Bc)β (3.19)
R(k) = (Bc − JBR(k))ρ + RP (3.20) Please note, that this baroreflex model is a simplification of Akkerman’s model. Please refer to Akkerman [32] for the original work.
3.3 Influence of gravity
The modelling of the orthostatic component describes how the angle α of the tilt-table influences the cardiovascular system and the physiological variables. We are only interested in the gravity component Fg along the body axis which is Fg = g · sin α. (3.21)
Based on that, the first model input u1(k) can be stated as follows:
u1(k) = sin α (3.22) Gravity acts on every single blood vessel in the cardiovascular system and creates rather large hydrostatic pressure differences in a standing human. Arterial pressure is decreased by 25 mmHG at head level and increased by 95 mmHG at leg level [7]. The question arises how gravitational forces can be integrated into the existing hemodynamic model. It is chosen to let gravity directly affect the right and left atrial pressures PR(k) and PL(k) which is a mathematically convenient alternative to modelling the whole hydrostatic column [32]. The atrial pressures then depend on the gravity factors ζR(k) and ζL(k):
(VV (k − 1) − VV 0)ζR(k) PR(k) = (3.23) CV (VP (k − 1) − VP 0)ζL(k) PL(k) = (3.24) CP ζR(k) = 1 − σR sin α(k) (3.25)
ζL(k) = 1 + σL sin α(k) (3.26)
21 Besides the atrial pressures, also the mean arterial pressure at the level of the baroreceptors PB(k) has to be corrected for the gravity influence. The reason is the height difference between the baroreceptors and the heart, where arterial pressure is evaluated. 1 2 P (k) = P (k) + P (k) − σ sin α(k) (3.27) B 3 S 3 D B 3.4 Influence of stepping
Akkerman’s model does not contain a component which describes the effects of the muscle pump when the stepping mechanism is activated. Therefore, these effects were analysed and subsequently added to the model.
The stepping mechanism acts on the cardiovascular system by activating the muscle pump through continuous leg movements. This has the following three immediate effects:
• Compression of the venous leg compartments leads to an increase of peripheral resistance.
• The contracting skeletal muscles decrease expandability of the venous vessels and hence, venous compliance is decreased.
• The muscle pump alters the functionality of the baroreflex mechanism. Similar to the situation of exercise, a “resetting” takes place and the hypothetic pressure level at which neural regulation is switched from parasympathetic to sympathetic action is increased.
Although the stepping mechanism moves the legs passively and we can only speak of a “passive” muscle pump, the stabilising effects on the cardiovascular system are still present, although diminished. Czell et al. [35] have conclu- ded after their pilot study with healthy adults, that passive leg movements stabilises blood circulation and prevents from syncopes. So fortunately, the stabilising effects on the cardiovascular system are still there and can be ex- ploited in the early rehabilitation process of neurological patients. The above listed effects are transformed to mathematical equations so that they can take influence on the existing cardiovascular model of Akkerman. As stepping is the second input after the inclination angle, u2(k) will be the expression for the normalised stepping frequency:
fstep(k) u2(k) = (3.28) fstep,max
22 steps where fstep,max is normally 48 min . As it takes some time for the cardiovas- cular system to adapt to the stepping movements, u2(k) has to be modelled as a first-order system with the time constant τstep which is usually chosen around 40 beats. In addition, the stepping influence at supine position has experimentally been found to be very low (figure 3.11, first 10 minutes). Thus, it is easiest to make u2(k) linearly dependent on u1(k). The adapted stepping input is denoted as κ(k):
−1 −1 τ τ κ(k + 1) = e step κ(k) + 1 − e step u2(k)u1(k) (3.29)
κ(k) now operates in an additive nature on peripheral resistance, venous compliance and the neural barosignal:
R(k) = (Bc − BS(k))ρ + RP + kSRκ(k) (3.30)
CV (k) = CV + kSC κ(k) (3.31)
B(k) = PB(k) + kP PP (k) − kB − kSBκ(k) 1 2 = P (k) + P (k) − σ sin α(k) + k P (k) − k − k κ(k) 3 S 3 D B P P B SB (3.32)
23 3.5 Model simulations
The cardiovascular model can now be used to simulate and analyse heart rate and blood pressure in response to various inputs. In addition, it is possible to investigate other cardiovascular signals such as stroke volume, peripheral resistance or cardiac output. In order to get an idea for what happens in the body during a tilt manoeuvre, a standard fast tilt-up and tilt-down should be examined first. Simulations have been done with standard steady-state values as given in table 3.1. These values were used in combination with a set of fixed parameters (table A.2) for identification of the unknown parameters (appendix A.5).
Table 3.1: Standard steady-state values used for the model simulations: “−” ◦ stands for supine position (α = 0 , fstep = 0); “+” stands for tilted position ◦ ◦ (α = 76 , fstep = 0); “s” stands for stepping (α = 76 , fstep = fstep,max) Steady-state value HR− 65 HR+ 80 − PS 120 + PS 125 − PD 80 + PD 95 HRs 75 s PS 130 s PD 95
3.5.1 Fast tilt-up and tilt-down The adjective “fast” refers to the fact that the tilt-table angle α changes from the minimal angle of zero degrees to the maximal angle of 76 degrees in two or three heart beats (vice versa for tilt-down). Of course, this is not feasible in reality where a full tilt may take up to 30 seconds. However, it is a good way to analyse the dynamics of such a fast tilt, which probably would not be that pronounced when tilting at a slower rate.
Fast tilt-up Model responses with the most important physiological variables are depicted in Figure 3.6. It can be seen that these responses are in accordance with the “standard” physiological response of tilt-up. Details about the dynamic cha-
24 racteristics will be explained in the following paragraph about fast tilt-down simulation. The reason is that tilt-down responses are usually much faster than tilt-up responses and that the dynamic features are easier to identify and explain.
Fast tilt-down When a person is tilted from the initial upright position back to the supine position, blood is shifted in the body under the influence of gravity. This has two immediate effects:
• Blood in the pulmonal pathways is shifted into the lung reservoir cau- sing a lack of blood in the left atrium.
• Blood from the venous reservoir is forced back to the right atrium and venous return is increasing rapidly.
The first effect leads to a fast decrease in arterial blood pressure and left stroke volume. As the blood supply in the left atrium is abruptly dimi- nished, the left stroke volume is immediately decreased according to the Frank-Starling law. This process is visible in the simulated model responses as the initial negative peak in blood pressure, left stroke volume and left atrial pressure. The second effect causes an immediate rise of right stroke volume in response to the increased venous return. After some time, this extra blood volume has made its way through the pulmonal pathways and ends up in the left atrium. This in turn causes the left stroke volume to rise again and leads to the positive blood pressure peak 7 to 8 seconds after the start of the tilt manoeuvre. In a third phase the phyisological signals settle to their steady-state values which is the case after approximately 20 seconds. Left and right stroke volume are balanced and the above description nicely shows how the Frank- Starling mechanism enables the adjustment of left and right stroke volume according to respective ventricle load.
3.5.2 Stepping As the stepping influence on the cardiovascular system is biggest when the table is fully tilted, only the simulation results for α = 76◦ are shown (Fi- gure 3.8). The heart rate shows the expected non-minimum phase behaviour as described by [6], the diastolic blood pressure is hardly influenced and the
25 systolic blood pressure rises, as described by [5]. The barosignal shows the in- verse behaviour of the heart rate, which makes perfect sense as the barosignal directly determines heart rate. Peripheral resistance is decreased when step- ping is activated which can be compared to the adaptation of the peripheral resistance to exercise. The increase of the stroke volumes and the pulse pres- sure point out the stabilising effect of stepping on the cardiovascular system.
3.5.3 Quasi-static The reason for a quasi-static simulation of the cardiovascular model is the analysis of the steady-state behaviour of heart rate and blood pressure at all angles α in the admitted range. Only the angle input is considered for this simulation because the stepping acts smoothly on the outputs whereas the system is expected to show rather different behaviour in the sympathetic and the parasympathetic region respectively. Remember that although the baroreflex regulation on heart rate is active over the whole range, periphe- ral resistance is only influenced by sympathetic regulation (see section 3.2). Figure 3.9 depicts the dependencies of the relevant cardiovascular variables on the inclination angle α. The following observations can be made: • Heart rate strictly increases with α and shows an “S” shape: The heart rate characteristics directly follow from the baroreceptor signal which is based on the arterial pressure at the level of the baro- receptors. • Systolic blood pressure both increases and decreases at lower angles, and strictly increases at higher angles: Systolic blood pressure PS is calculated as the sum of diastolic blood pressure PD and pulse pressure PP . At small angles, stroke volumes don’t change much, but PD is increased. This leads to the increase in PS at small angles. However, as soon as the stroke volumes and subsequently the pulse pressure is decreased, PS is decreased as well. In the sympathetic regulation domain, PS strictly increases because PD grows faster than PP declines. • Diastolic blood pressure strictly increases with α, but at a lower rate at lower angles: The reason is that at higher angles peripheral resistance is increased by the baroreflex which leads to higher arterial pressures. • Peripheral resistance by design only increases at higher angles, when sympathetic regulation becomes active.
26 • Stroke volumes are diminished when tilting.
Figure 3.10 compares the results from the quasi-static simulation with results from Hainsworth [10], Matalon [36], Heldt [28], Fisler [37] and Nguyen [5]. It can be deduced that the accordance of the model results with literature studies and previous work at SMS is satisfying.
27 90 80 60 80 40 [deg] [bpm] Angle 70 20 Heart rate 0 60 0 50 100 0 50 100 Heart beats Heart beats 140 90
120 80 P S P 70 100 D [mmHG] [mmHG]
Barosignal 60
Blood pressure 80 50 0 50 100 0 50 100 Heart beats Heart beats 5 P 80 Q L L 4 60 P [ml] R [mmHG] 3 Q R
Stroke volumes 40 Atrial pressures 2 0 50 100 0 50 100 Heart beats Heart beats 5.5 1600 5 1400 4.5 [l/min] 1200 4 Cardiac output [mmHG ms/ml]
Peripheral resistance 1000 3.5 0 50 100 0 50 100 Heart beats Heart beats
Figure 3.6: Simulation of a fast tilt-up without stepping
28 90 80 60 80 40 [deg] [bpm] Angle 70 20 Heart rate 0 60 0 50 100 0 50 100 Heart beats Heart beats 140 90
120 80 P S 70 100 [mmHG] [mmHG] P
D Barosignal 60
Blood pressure 80 50 0 50 100 0 50 100 Heart beats Heart beats 5 P 80 Q L R 4 60 [ml] Q L [mmHG] 3 P R 40 Stroke volumes Atrial pressures 2 0 50 100 0 50 100 Heart beats Heart beats 5.5 1600 5 1400 4.5 [l/min] 1200 4 Cardiac output [mmHG ms/ml]
Peripheral resistance 1000 3.5 0 50 100 0 50 100 Heart beats Heart beats
Figure 3.7: Simulation of a fast tilt-down without stepping
29 82
40 80 78
20 [bpm] Stepping Heart rate [steps/min] 76
0 74 0 50 100 150 0 50 100 150 Heart beats Heart beats 132 65
130 60 128
[mmHG] [mmHG] 55
126 Barosignal
124 50
Systolic blood pressure 0 50 100 150 0 50 100 150 Heart beats Heart beats 100 65
60 95 [ml]
[mmHG] 55 Stroke volumes 90 50
Diastolic blood pressure 0 50 100 150 0 50 100 150 Heart beats Heart beats 1550 4.6
1500 4.4
1450 [l/min] 4.2 Cardiac output [mmHG ms/ml]
Peripheral resistance 1400 4 0 50 100 150 0 50 100 150 Heart beats Heart beats
Figure 3.8: Simulation of an activation of the stepping mechanism (α = 76◦)
30 80 80
60 75 40 [deg] [bpm] Angle 70 20 Heart rate 65 0 0 500 1000 1500 0 500 1000 1500 Heart beats Heart beats 100 126 124 80 122
[mmHG] [mmHG] 60 120 Barosignal 118 40
Systolic blood pressure 0 500 1000 1500 0 500 1000 1500 Heart beats Heart beats 100 80
90 70 [ml]
[mmHG] 80 60 Stroke volumes 70 50
Diastolic blood pressure 0 500 1000 1500 0 500 1000 1500 Heart beats Heart beats 4.8 1500 4.6 1400 4.4 [l/min] 1300 4.2 Cardiac output [mmHG ms/ml]
Peripheral resistance 1200 4 0 500 1000 1500 0 500 1000 1500 Heart beats Heart beats
Figure 3.9: Quasi static simulation without stepping
31 20 20 Hainsworth Hainsworth Matalon Smith 18 18 Fisler Model Model 16 16
14 14
12 12
10 10
HR [bpm] 8 8 ∆ dBP [mmHG] ∆ 6 6
4 4
2 2
0 0
−2 −2 0 20 40 60 80 0 20 40 60 80 Angle [deg] Angle [deg]
Figure 3.10: Comparison of steady-state behaviour. Left: HR as a function of α. Right: Diastolic BP as a function of α
32 3.6 Model validation
The step of model validation will be performed using measurements from three healthy subjects (see chapter 5.1.3). Evaluation will be done in a qualitative way analysing each measurement separately. Although averaging over all subjects would probably yield better agreement between the model simulation and the measurement, interesting details from the individual cases would be lost. The measurement was divided into an identification and a validation part. The according measurement protocol is illustrated in the lowermost plot of figure 3.11. Note that between the identification and the validation there was a recalibration of the blood pressure measurement device. This can introduce offsets in some cases whereas diastolic blood pressure seems to be affected the most. In figure 3.11 for example, this offset amounted to about 4 mmHg and has been corrected accordingly. Validation results for the first subject (MW) are satisfying and it demons- trates that it is possible to simulate or predict heart rate and blood pressure dynamics.
80
60 HR [bpm]
0 10 20 30 40 50 60 70 Time [min] 120
100
sBP [mmHG] 80 0 10 20 30 40 50 60 70 Time [min]
80
60
dBP [mmHG] 40 0 10 20 30 40 50 60 70 Time [min]
f α step 50 [deg] / [steps/min] α
step 0 f 0 10 20 30 40 50 60 70 Time [min]
Figure 3.11: Model validation with subject MW
33 However, the identified models for the other two subjects deviate more from the measured signal than it was the case for the first subject. This should be analysed in more detail: For the second subject (figure 3.12) it can be said that heart rate and diastolic blood pressure were well reproduced by the mo- del. Systolic blood pressure, however, did not show clear trends: During the identification phase, systolic blood pressure stayed constant when tilting but increased in the end of the experiment during the slow ramp of the inclina- tion angle. Another issue are the calibration offsets, that have already been mentioned above. For the second subject, diastolic blood pressure jumped by about 10 mmHg which has been corrected for in the modelled diastolic blood pressure curve. Already the low values of about 40 mmHg after 17 minutes are unrealistic compared to the baseline values at the beginning of the experiment which were around 55 mmHg. The worse thing however is that after the recalibration the value is not set back to 55 mmHg but is even increased to about 65 mmHg. These huge jumps in the measured signals are physiologically improbable in such a short timespan and it unveils the weaknesses of the blood pressure measurement device.
90 80 70
HR [bpm] 60 50 0 10 20 30 40 50 60 70 Time [min] 120
100
sBP [mmHG] 80 0 10 20 30 40 50 60 70 Time [min] 80
60
dBP [mmHG] 40 0 10 20 30 40 50 60 70 Time [min]
f α step 50 [deg] / [steps/min] α
step 0 f 0 10 20 30 40 50 60 70 Time [min]
Figure 3.12: Model validation with subject MSW
34 The validation measurement for the third subject emphasises the above men- tioned problems: First, systolic blood pressure is hard to reproduce or model. Second, diastolic blood pressure measurement is tampered with calibration offsets. However, it has to be added that the last measurement is an extreme example for what can happen with physiological signals.
100
80
60 HR [bpm]
0 10 20 30 40 50 60 70 Time [min] 180 160 140 120 100 sBP [mmHG] 80 0 10 20 30 40 50 60 70 Time [min] 80
60 dBP [mmHG] 0 10 20 30 40 50 60 70 Time [min]
f α step 50 [deg] / [steps/min] α
step 0 f 0 10 20 30 40 50 60 70 Time [min]
Figure 3.13: Model validation with subject DH
Generally, it can be concluded that the model reproduces heart rate and diastolic blood pressure with a satisfying accuracy. Problems occur, if repro- ducibility is not given, i.e. if the subject responds differently for the same inputs. However, this is a more general issue as any deterministic model would struggle with low reproducibility. In contrast to heart rate and diastolic blood pressure, systolic blood pressure is more difficult to predict and model for healthy subjects. Still, we decided to go on with the strategy of controlling all three variables, because in the end, the system will be used with patients. Patients usually react much bet- ter with systolic blood pressure to verticalisation because neural regulation is impaired. To counter the problem with calibration offsets, it will be important not to
35 do a calibration between identification and control or to reidentify the model when a calibration is necessary. If the offset is only small, it may not be nee- ded to reidentify the model because for control only the relative responses are important. However, if the offset is too high and the new values are out of the identified range, it gets impossible to start the control experiment without adapting the model.
36 Chapter 4
Control design
Starting from the non-linear MIMO system derived in the previous chapter, an appropriate controller will now be developed. This controller has to be able to keep heart rate and blood pressure within reasonable bounds and minimise fluctuations by adjusting the inclination angle and the stepping frequency. For the given MIMO system which is apparently strongly coupled, a SISO approach trying to control each output in an isolated manner seems infeasible. The fact that the system has less inputs than outputs makes it even harder to do so. It is therefore advisable to choose a control strategy which can handle these issues. A linear optimal control approach is frequently used in advanced control applications, and has been chosen for this thesis as well. Two controllers have been developed, implemented and tested: a Linear Quadratic Regulator (LQR) and a Model Predictive Controller (MPC). The LQ controller which was augmented by an integral part in order to eliminate steady-state control errors was experimentally found to be very hard to tune. The reason is that the system is likely to operate near or on the constraints boundaries for the control inputs. As a result, control inputs are often saturated and anti-windup strategies are therefore necessary. However, for a true MIMO system with strong couplings as has been developed in the previous chapter it is difficult to apply standard anti-windup techniques. It is therefore desirable to have a technique which intrinsically accounts for input constraints. This explains why Model Predictive Control (MPC) is suited for the application at hand and why it is to prefer to a common LQ regulator.
4.1 Model predictive control (MPC) design
The advantages of Model Predictive Control are manifold. Two of the most important features are that MPC takes account of actuator limitations and
37 that it is suited for multivariable control problems. The principle behind MPC is as follows: based on the system model the controller predicts future outputs and finds the optimal control inputs by minimising a certain cost function. An intuitive analogon for MPC is driving a car [38]. Imagine that the reference path is the lane, the plant is the car and the controller is represented by the driver. The control objective is to keep the car on the lane, while steering as little as possible, keeping a certain distance to the kurbs, obey speed limitations and so on. The driver now has an internal belief or model of how the car reacts to his inputs. He uses this knowledge to predict future behaviour of the car and give according control inputs in order to stay on the reference path, minimise steering effort and meeting all given constraints. This control problem can generally be formulated with a cost function and according constraints [39].
T T min J = (rs − yp) Q(rs − yp) + ∆u R∆u (4.1) ∆u M∆u ≤ γ (4.2) where ∆u(mNc × 1) is the control input, yp(pNp × 1) the predicted output, rs(pNp × 1) the reference, Q(pNp × pNp) the output weighting matrix and R(mNc × mNc) the control input weighting matrix. The matrix M(4mNc × mNc) and the vector γ(4mNc × 1) define the constraints. The scalars n, m and p are the number of states, the number of inputs and the number of outputs of the MIMO system. The control horizon is denoted as Nc. Table 4.1 summarises these notations. As we will only need constraints on the control inputs, state and output constraints are neglected in this formulation. Note that the vectors rs and yp contain Np discrete samples over the prediction horizon and the vector ∆u contains Nc discrete samples over the control horizon: