<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">OJAppS</journal-id><journal-title-group><journal-title>Open Journal of Applied Sciences</journal-title></journal-title-group><issn pub-type="epub">2165-3917</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ojapps.2013.31009</article-id><article-id pub-id-type="publisher-id">OJAppS-29447</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Biomedical&amp;Life Sciences</subject><subject> Chemistry&amp;Materials Science</subject><subject> Computer Science&amp;Communications</subject><subject> Engineering</subject><subject> Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  An Optimal Control Problem for Hypoxemic Hypoxia Tissue-Blood Carbon Dioxide Exchange during Physical Activity
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>ean</surname><given-names>Marie Ntaganda</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Benjamin</surname><given-names>Mampassi</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Department of Mathematics and Computer Science, Dakar University, Dakar, Senegal</addr-line></aff><aff id="aff1"><addr-line>National University of Rwanda, Butare, Rwanda</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>jmnta@yahoo.fr(EMN)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>29</day><month>03</month><year>2013</year></pub-date><volume>03</volume><issue>01</issue><fpage>56</fpage><lpage>61</lpage><history><date date-type="received"><day>November</day>	<month>2,</month>	<year>2012</year></date><date date-type="rev-recd"><day>December</day>	<month>3,</month>	<year>2012</year>	</date><date date-type="accepted"><day>December</day>	<month>10,</month>	<year>2012</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
   This paper aims at solving an optimal control problem for determining the response of hypoxia to heart rate and alveolar ventilation that are cardiovascular and respiratory control respectively during a physical activity. A two nonlinear coupled ordinary differential equations is presented. The cost function of optimal control problem is discretized using the linear B-splines functions defined on a regular grid. The results show the determinant parameters stabilized at normal value. 
 
</p></abstract><kwd-group><kwd>Arterial Oxygen Pressure; Arterial Dioxide Pressure; Cardiovascular/Respiratory System; Optimal Control Problem; Numerical Simulation</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Hypoxia, or hypoxiation, is defined as a pathological condition related to adequate oxygen supply in human body. It is in two main types: the generalized hypoxia that is characterized by the deprived adequate oxygen supply in whole body and tissue hypoxia which happens in its region. It differs from hypoxemia called also hypoxaemia in that within the arterial blood the oxygen concentration is abnormally low. Hypoxemia was originally defined as a deficiency of oxygen in arterial blood but standard manuals take this to mean an abnormally low partial pressure of oxygen, content of oxygen or percent saturation of hemoglobin with oxygen, either found singly or in combination. The serious cases of the hypoxemia happen when the decreased partial pressure of oxygen in blood is less than 60 mmHg.</p><p>In addition, the generalized hypoxia occurs in healthy people when they ascend to high altitude, where it causes altitude sickness leading to point constitutes the beginning of the steep portion of the hemoglobin dissociation curve, where a small decrease in the partial pressure of oxygen results in a large decrease in the oxygen content of the blood or when hemoglobin oxygen saturation is less than 90%. The reason of this is this potentially fatal complications including high altitude pulmonary edema (HAPE) and high altitude cerebral edema (HACE) [<xref ref-type="bibr" rid="scirp.29447-ref1">1</xref>]. It also occurs in healthy individuals when breathing mixtures of gases with low oxygen content.</p><p>Hypoxic hypoxia is a result of insufficient oxygen available to the lungs. The examples of how lungs can be deprived of oxygen are a blocked airway, a drowning or a reduction in partial pressure (high altitude above 10,000 feet). Hypoxia is also a serious consequence of pre-term birth in the neonate.</p><p>The main cause for this is that the lungs of the human fetus are among the last organs to develop during pregnancy.</p><p>To assist the lungs to distribute oxygenated blood throughout the body, infants at risk of hypoxia are often placed inside an incubator capable of providing continuous positive airway pressure (also known as a humidicrib). The insufficient delivery of oxygen (low<img src="9-2310102\f4f85dbf-ebdb-442b-9128-c65e7e1ad870.jpg" />) or inability to utilize oxygen (normal<img src="9-2310102\94df10a2-14e8-4416-9f4e-f268d04b8b16.jpg" />) causes also the hypoxia where we assist to oxygen deficiency at the mitochondrial sites. This phenomeno accurs when <img src="9-2310102\3b642546-d4b9-429b-8d28-0564924ceef4.jpg" /> less than 7.3 kPa (55 mmHg). Below this threshold the ventilation starts to stimulate carotid body activity. The hyperventilation reduces <img src="9-2310102\a5c0000f-1908-4a02-a9ec-a9dcf983c842.jpg" /> and<img src="9-2310102\84c534f6-5e57-4ae4-899b-9506fb1d1afc.jpg" />, which limits the initial rise in ventilation, because it decreases the carotid body and central chemoreceptor stimuli. In fact, in humans, hypoxia is detected by chemoreceptors in the carotid body. This response does not control ventilation rate at normal<img src="9-2310102\d9864fb2-b530-48d0-9d09-62683aabe83d.jpg" />, but below normal the activity of neurons innervating these receptors increases dramatically, so much so to override the signals from central chemoreceptors in the hypothalamus, increasing<img src="9-2310102\4a9ebfed-5877-4167-97e4-621c34b8f33a.jpg" /> despite a falling<img src="9-2310102\5532fb9e-9f91-46a5-a4db-5cfd5a85769a.jpg" />.</p><p>Any physical activity obviously causes the body to demand more oxygen for normal functioning. The muscles rob the brain of the marginal amounts of oxygen available in the blood and the time of onset of hypoxic symptoms is shortened. However, the improvement of performance of athlete in high altitude results in a mild and non-damaging intermittent hypoxia used intentionally during training to develop an athletic performance adaptation at both the systemic and cellular level. Mathematical models quantifying hypoxic hypoxia have been proposed [2-4]. The optimal control problem based on the responses of cardiovascular respiratory system parameters to its controls, heart rate and alveolar ventilation, during physical activity has not been considered in the situation of hypoxic hypoxia. This work focuses on this issue where the mathematical model is modified to include the controls of cardiovascular respiratory system.</p><p>The remainder of this paper is structured as follows. In Section 2, we present a mathematical model, an optimal control problem and its descretization. The results of numerical simulation are discussed in Section 3. In Section 4 deals with concluding remarks.</p></sec><sec id="s2"><title>2. Setting of an Optimal Control Problem</title><p>The model we present in this paper involves modifying of model equations as developed by Guillermo Gutierrez [<xref ref-type="bibr" rid="scirp.29447-ref2">2</xref>] in order to include the role of physical activity. The diagram for a two compartmental model is illustrated in the <xref ref-type="fig" rid="fig1">Figure 1</xref> where mass transport model of tissue CO<sub>2</sub> exchange is developed to examine the relative contributions of blood flow and cellular hypoxia (dysoxia) to increases in tissue and venous blood CO<sub>2</sub> concentration.</p><p>From this compartmental diagram the model equations are as follows.</p><disp-formula id="scirp.29447-formula149801"><label>(1)</label><graphic position="anchor" xlink:href="9-2310102\863d0304-18ad-4da1-9648-9a8382b2d89d.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.29447-formula149802"><label>(2)</label><graphic position="anchor" xlink:href="9-2310102\e3ffbfb3-7b5b-467d-870b-6bf540502824.jpg"  xlink:type="simple"/></disp-formula><p>The integral role of physical activity results in the influence of the demand of the tissues for oxygen during hypoxia in altitude. Just as resting ventilation increases dramatically at high altitude, so does ventilation during physical activity. In fact, at moderate levels of physical activity, there is little or no change in arterial <img src="9-2310102\41a1c135-c33d-411e-b47a-d2d41156089f.jpg" />between rest and exercise. Since carbon dioxide production for a given work level is essentially independent of altitude, this means that measured ventilation is independent of altitude at a given work level. At work levels approaching maximal values at any altitude, alveolar and arterial <img src="9-2310102\1c2fa0bb-0761-4694-8de1-3eae7a50dd5e.jpg" />fall compared with the resting level and physical activity ventilation measured at correspondingly rises. Furthermore, during exercise, increases in alveolar ventilation must parallel the increased tissue oxygen consumption and carbon dioxide production by the exercising muscles, both of which rise in direct proportion to the increase in power output. These relationships are governed by the following equation</p><p><img src="9-2310102\add7d8ff-c06a-43bb-bb3d-2405ad8c6999.jpg" /></p><p>where <img src="9-2310102\bb04c4aa-54ad-4a0d-b102-d103134b6434.jpg" /> is alveolar ventilation and <img src="9-2310102\669f7fd5-6394-4cf1-a4c3-855855dcf968.jpg" /> denotes a constant [<xref ref-type="bibr" rid="scirp.29447-ref5">5</xref>], so that we have</p><p><img src="9-2310102\67fb0c00-b821-4e56-bcc0-65082cf9ca21.jpg" /></p><p>where <img src="9-2310102\7ff4089a-9ca1-47da-94ba-2669a80fb1ef.jpg" />the slope of the physiological <img src="9-2310102\77bb0fa5-b755-4af6-9ac2-2104c7bc6c6f.jpg" />dissociation curve and <img src="9-2310102\aad5ff6e-e5e8-4bcf-ac4d-d8e450e1d20f.jpg" />constant for the physiological <img src="9-2310102\d37e6d7b-9b12-4723-91b0-d6fe72075160.jpg" /><sub> </sub>dissociation curve [<xref ref-type="bibr" rid="scirp.29447-ref6">6</xref>].</p><p>Taking RQ as respiratory exchange ration and <img src="9-2310102\daf7fef7-2d60-4f6b-8ccc-66a9833da276.jpg" />as consummation rate of O<sub>2</sub>, the relationship between <img src="9-2310102\5fb0205e-f4b1-484f-bfeb-7e3bd68ea7e0.jpg" /> and <img src="9-2310102\c65b2358-73fd-4659-8e8c-443d4c3122b3.jpg" /> satisfies the following relation [<xref ref-type="bibr" rid="scirp.29447-ref2">2</xref>].</p><disp-formula id="scirp.29447-formula149803"><label>(3)</label><graphic position="anchor" xlink:href="9-2310102\6d62eec6-ff87-4bce-a1c0-ff24fed44fc0.jpg"  xlink:type="simple"/></disp-formula><p>Using Fick’s principle [<xref ref-type="bibr" rid="scirp.29447-ref2">2</xref>] applied to the relation (3) allows to get</p><p><img src="9-2310102\e7fe86ba-06ad-495a-a75b-5e763b868349.jpg" /></p><p>where <img src="9-2310102\66440777-85a9-4e98-bea4-7a8e95ecd033.jpg" /> (resp.<img src="9-2310102\f1bdb78c-9a4d-4147-8a32-b3a591129dd9.jpg" />) is arterial (venous) concentration of<img src="9-2310102\45675a53-2292-4bc6-a8d4-5a351a6509de.jpg" />, <img src="9-2310102\c8813c64-5349-4aa7-8c5e-641ba4e6c507.jpg" />denotes heart rate and <img src="9-2310102\b9e0ae89-8913-42a9-98e6-1dd2c8338d3a.jpg" /> represents stroke volume.</p><p>In addition, it is known that the human respiratory control system varies the ventilation rate <img src="9-2310102\b08ca541-0b34-4e06-b61f-fd05c81f1dfd.jpg" /><sub> </sub>in response to the levels of <img src="9-2310102\127831f4-e9fa-4118-8391-553a898a9c22.jpg" /> and <img src="9-2310102\9fda09cd-c5a2-4e9b-a7aa-43bb86b6fbfc.jpg" /> in the body and the control mechanisms of cardiovascular system influences global control in the blood vessels as well as well as heart rate <img src="9-2310102\6dbca9b9-1d3c-4ba5-aea3-4349eeabae0f.jpg" /> for impacting blood flow <img src="9-2310102\05a6ab1a-a67a-463c-9ca2-a2b62c638c0b.jpg" /> [<xref ref-type="bibr" rid="scirp.29447-ref7">7</xref>]. Generally, during physical activity in altitude and particular in the hypoxia case, the control mechanism of these two systems plays a crucial role.</p><p>Finally, we are interested in the following model equations</p><disp-formula id="scirp.29447-formula149804"><label>(4)</label><graphic position="anchor" xlink:href="9-2310102\7b19186e-b68d-404f-a32a-487a898e9207.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.29447-formula149805"><label>(5)</label><graphic position="anchor" xlink:href="9-2310102\c45d827e-8ef8-41d7-899f-cf7075d1574e.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.29447-formula149806"><label>(6)</label><graphic position="anchor" xlink:href="9-2310102\375c1d7b-4a21-42eb-a16a-47ab8a369405.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.29447-formula149807"><label>(7)</label><graphic position="anchor" xlink:href="9-2310102\142e9973-5619-4f38-9294-a9de8cf8f24a.jpg"  xlink:type="simple"/></disp-formula><p>where two last ordinary differential equations describe the control of cardiovascular and respiratory system is described respectively with <img src="9-2310102\a5906954-a00e-4832-8c8a-ebf440fc1572.jpg" />and <img src="9-2310102\61e9f24f-7316-4239-a538-e5f295b95cbd.jpg" /> the functions to be determined by an optimality criterion.</p><p>The alveolar gas equation allows the calculation of the alveolar partial pressure of oxygen as follows [<xref ref-type="bibr" rid="scirp.29447-ref8">8</xref>]</p><p><img src="9-2310102\4c2427c7-96cb-4dac-a0cc-23b85f88c60e.jpg" /></p><p>where <img src="9-2310102\734070c7-647e-4679-bcd7-761eb41ff0ae.jpg" />is the alveolar partial pressure of oxygen, <img src="9-2310102\12cf77a7-c61c-4af4-bc8a-e25c848158b2.jpg" />denotes the prevailing atmospheric pressure, <img src="9-2310102\09fd7711-d0fc-493d-a6d6-54feb11a6cfe.jpg" />represents the saturated vapor pressure of water at body temperature and the prevailing atmospheric pressure and <img src="9-2310102\089bd69b-12f8-4879-946d-237586fa8de5.jpg" />is the fraction of inspired gas that is oxygen (expressed as a decimal). In addition, the relation between alveolar partial and arterial pressure of oxygen is given by</p><p><img src="9-2310102\bef99d82-de7e-4ac5-9298-2d628c467595.jpg" /></p><p>because the<img src="9-2310102\d9de9572-a57d-4fbf-822c-3610e4bb7a57.jpg" />to <img src="9-2310102\a8c1d9f9-37f7-44fa-b240-edc3f2e480d6.jpg" />gradient is normally close to and is written as follows:</p><p><img src="9-2310102\a7d078b8-5f89-4a77-ac13-b6f1e75755f4.jpg" /></p><p>Similarly, blood <img src="9-2310102\a94d10a7-6778-4093-b6b7-42d2a35a1c38.jpg" />is calculated on the basis of the Henderson-Hasselbach equation [<xref ref-type="bibr" rid="scirp.29447-ref9">9</xref>] as follows.</p><p><img src="9-2310102\8d0b3bff-cb01-45d9-a40e-867465fdff94.jpg" /></p><p>where <img src="9-2310102\bd5e342c-0896-4a25-9711-b7b2bb0b05a5.jpg" /> is <img src="9-2310102\4540e9c7-bbc5-4eb2-9b7c-16552d0e3280.jpg" />content of plasma defined by Douglas [<xref ref-type="bibr" rid="scirp.29447-ref10">10</xref>] as</p><p><img src="9-2310102\d916720d-7639-45dc-b45d-8cd95d239d17.jpg" /></p><p>It results that the alveolar carbon dioxide pressure of oxygen <img src="9-2310102\a39e54e1-520b-4f82-acda-5723717a68e6.jpg" /> is equivalent to <img src="9-2310102\c4d1a4eb-9f28-474c-b7ff-d5eddffdd03a.jpg" /><sub> </sub>(there is no gradient).</p><p>Furthermore, it appears that a main goal of respiratory control is to keep <img src="9-2310102\d5e6df03-03fe-477d-ad98-a9d4a9547f2e.jpg" /><sub> </sub>venous partial pressure as close as possible to an equilibrium value denoted by <img src="9-2310102\ec5b65f2-212d-43d9-8e7b-2afd12656b76.jpg" /> and, to a lesser extent, control <img src="9-2310102\7234f943-ff37-4f85-9fdf-0f4d89706252.jpg" /><sub> </sub>to the equilibrium such that the cost functional can be formulated in the following way.</p><p>Find <img src="9-2310102\a7fdc924-5ba6-4e42-aa85-18a2b7bc5277.jpg" />and <img src="9-2310102\fb14ad97-9223-4eac-be45-78ccdaa2def4.jpg" />solution of</p><disp-formula id="scirp.29447-formula149808"><label>(8)</label><graphic position="anchor" xlink:href="9-2310102\fc16c42a-64f5-4634-b032-9ce2308f1b3a.jpg"  xlink:type="simple"/></disp-formula><p>subject to the systems (4 )-(7).</p><p>In the relation (8), the positive scalar coefficients <img src="9-2310102\4ee25692-e51f-4bb1-8cdc-6d1fee5327bf.jpg" /> and <img src="9-2310102\74b5990e-cfef-40c0-a428-812bcfa2ac05.jpg" /> determine how much weight is attached to each cost component term in the integrand while <img src="9-2310102\79b3feda-f696-4128-a9df-d8d6f2d7280d.jpg" /><sub> </sub>denotes the maximum time that the physical activity can take.</p><p>Let us consider <img src="9-2310102\0ca210d5-9b67-412c-98a6-33183bbfb9c7.jpg" />the vector space that is span of</p><disp-formula id="scirp.29447-formula149809"><label>(9)</label><graphic position="anchor" xlink:href="9-2310102\f3642c90-f356-408a-a9c6-8eb101a7dbc3.jpg"  xlink:type="simple"/></disp-formula><p>a base of linear B-splines functions on a regular grid</p><disp-formula id="scirp.29447-formula149810"><label>(10)</label><graphic position="anchor" xlink:href="9-2310102\70429e29-cc94-49aa-b7b8-3cb62963ea5e.jpg"  xlink:type="simple"/></disp-formula><p>The functions <img src="9-2310102\dc16b9dd-b5e7-46df-94ff-6f6ff797dc2e.jpg" /> verify the following relation <img src="9-2310102\76ccf5a6-1a1c-4492-bf90-979a814563cd.jpg" /> where <img src="9-2310102\42a28fca-92a0-4a5c-910e-832cdcd14eb8.jpg" /> denotes Kronecker symbol. The descretization of the optimal problem (8) is done by setting the state vector</p><p><img src="9-2310102\0d4edaa6-59f2-4726-8e6f-c2839031b592.jpg" /></p><p>and the desired final vector</p><p><img src="9-2310102\537d494f-3ab6-49cb-a411-3a6a82f84094.jpg" /></p><p>such that it can be written as follows.</p><disp-formula id="scirp.29447-formula149811"><label>(11)</label><graphic position="anchor" xlink:href="9-2310102\85921cb2-7482-48f3-85b5-95b834fcfa65.jpg"  xlink:type="simple"/></disp-formula><p>where</p><p><img src="9-2310102\2c5c0ad7-f07a-4a6b-911a-4a4741669e0d.jpg" /></p><p>with <img src="9-2310102\d7868628-d832-4f1f-99ee-d1336d6d72b6.jpg" /> and <img src="9-2310102\bc2c96e2-ec88-4e3c-a6c9-7f612b4af580.jpg" /> respectively the <img src="9-2310102\647b85ce-a026-42a3-a830-ee5379a915ee.jpg" />component of the vectors <img src="9-2310102\952981c1-015c-4579-a5ba-6e7536f0a8d4.jpg" /> and<img src="9-2310102\d3665d00-9f07-48fc-bbd2-9c41eead049d.jpg" />.</p><p>We are looking for <img src="9-2310102\aff2b98a-b588-45c8-8ecd-a03ca22176b6.jpg" />a approximated solution of (11) in the set <img src="9-2310102\589fb328-9df7-48e2-ad50-84c3195a2eb1.jpg" />such that</p><disp-formula id="scirp.29447-formula149812"><label>(12)</label><graphic position="anchor" xlink:href="9-2310102\2172e887-fa17-46b4-930c-08c0fb9e9671.jpg"  xlink:type="simple"/></disp-formula><p>Therefore the cost function (11) becomes</p><disp-formula id="scirp.29447-formula149813"><label>(13)</label><graphic position="anchor" xlink:href="9-2310102\9bd84328-c9ea-4533-8129-bafe722272a3.jpg"  xlink:type="simple"/></disp-formula><p>with <img src="9-2310102\fb8e54b8-9d05-4a0f-a1af-085a355cce4b.jpg" /> and finally the discrete formulation of optimal problem (8) subject to (4)-(7) is written as follows</p><disp-formula id="scirp.29447-formula149814"><label>(14)</label><graphic position="anchor" xlink:href="9-2310102\0ee1b718-8f28-4f13-b474-da902b607273.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="9-2310102\4eaa231b-e92e-4e8a-96c5-c9cf5613903d.jpg" /> is a matrix <img src="9-2310102\ccf9e4e1-9800-41ca-ad1d-0916335a8626.jpg" /> such that the components <img src="9-2310102\17973537-28b3-4bae-815a-976e5a3d97e7.jpg" /> are components of the function <img src="9-2310102\18996c2e-0426-4c2f-b541-eb80e64be604.jpg" /> in the set <img src="9-2310102\2fc596bc-3864-41d2-90fb-19fbf4c1bd1f.jpg" /> and <img src="9-2310102\9e560d80-522a-420f-a1f8-12b19303ca99.jpg" /> represents the matrix with <img src="9-2310102\49dc1d49-5274-4688-b3ac-92142341b7c2.jpg" /> component is</p><p><img src="9-2310102\23ff2eca-fb47-425b-9f87-5be5859acb3f.jpg" /></p><p>where <img src="9-2310102\ae2b0bd3-9b9f-42b2-b703-5f43414fa27e.jpg" /> denotes the solution of the systems (4)-(5) associated to<img src="9-2310102\000b2ed1-34f7-4e5e-b20c-930573595fcc.jpg" />, <img src="9-2310102\78a0b328-4f3c-4555-af36-bc00eecdaa41.jpg" />and <img src="9-2310102\6117d27b-1f6c-4e1f-a085-5e2d7de921c5.jpg" /> are matrix defined by</p><disp-formula id="scirp.29447-formula149815"><label>(15)</label><graphic position="anchor" xlink:href="9-2310102\be5c2352-f1d6-4982-9327-16a12d31be9d.jpg"  xlink:type="simple"/></disp-formula></sec><sec id="s3"><title>3. Test Results</title><p>For solving the optimal control problem (13) subject to the constraints (4)-(7) we consider the parameters presented in the <xref ref-type="table" rid="table1">Table 1</xref>.</p><p>We take <img src="9-2310102\671d5f2a-3cd5-4963-bcbc-7f06a6a3bd2f.jpg" /> and <img src="9-2310102\754b75e6-d2bf-47e1-af82-4a591ed1cd63.jpg" /> as normal value for a healthy individual that <img src="9-2310102\739a8bea-1530-40c7-8987-90a9bf489b4c.jpg" /> and <img src="9-2310102\147e011d-053b-4ab8-8a43-c7c460a5aa9c.jpg" /></p><p>= 95 mmHg. The variation curve of carbon dioxide in the tissues and carbon dioxide in vascular is illustrated in the <xref ref-type="fig" rid="fig2">Figure 2</xref>. The control of cardiovascular and respiratory systems, <img src="9-2310102\9f8a8b4d-2010-4b07-9064-a800ecd6ed6b.jpg" />and <img src="9-2310102\8c0c5b81-6ef5-4ffa-a791-72e77d842363.jpg" /> against the time are given in the <xref ref-type="fig" rid="fig3">Figure 3</xref> whereas the <xref ref-type="fig" rid="fig4">Figure 4</xref> shows the responses of these controls.</p><p><xref ref-type="table" rid="table1">Table 1</xref>. Values of parameters used in numerical simulation.</p><p><img src="9-2310102\e56872ef-00b5-4bcf-a4b7-23ad8a6246ae.jpg" /></p><p>The <xref ref-type="fig" rid="fig2">Figure 2</xref> shows a decrease of arterial and venous dioxide carbon concentration against the time. They are maintained at a level where the variation is small. This is due to the effect of ventilation during physical activity. In fact, ventilation increases abruptly in the initial stages of exercise and is then followed by a more gradual increase. This mechanism of increase results of to motor centre activity and afferent impulses from proprioceptors of the limbs, joints and muscles. Since peripheral chemoreceptors are responsible for increasing ventilation, Central chemoreceptors may be readjusted to increase ventilation to maintain carbon dioxide concentrations.</p><p>At the onset of physical activity, the heart rate and alveolar ventilation increase. Generally, heart rate increases to about 90% of its maximum values during strenuous physical activity. Furthermore, the ventilation increases with increases in work rate at submaximal physical activity intensities. These physiological effects of physical activity on cardiovascular-respiratory system are justified by the variation of its controls in the <xref ref-type="fig" rid="fig3">Figure 3</xref> where they reach a value and they are stabilized with small oscillations.</p><p>The <xref ref-type="fig" rid="fig4">Figure 4</xref> shows the arterial carbon dioxide (resp. oxygen) decreases (increases) at the onset of physical activity to be stabilized at normal value. The changes which occur in<img src="9-2310102\e881e226-c0bc-49e6-99e1-d15378dc6871.jpg" /> and <img src="9-2310102\51ff2da1-d608-4182-8092-735d517ee961.jpg" />values during exercise are usually small. But since hypoxia is characterized by the adequate oxygen supply deprived in whole (generalized hypoxia) or a region of the body (tissue hypoxia) the ventilation process during physical activity plays an integral role. During physical activity, when sufficient oxygen for flux through the is not available, the increased reliance on glycolysis results in increased accumulation of lactic acid, which initially leads to an increase in<img src="9-2310102\c993e878-50b7-4fec-bfdd-5726b3e84b7e.jpg" />. However, this is counteracted by the stimulation of ventilation and as a result <img src="9-2310102\b880751c-1d86-4eee-a2e2-3fb070883790.jpg" />is decreased. In addition, in</p><p>humans, hypoxia is detected by chemoreceptor in the carotid body. This response does not control ventilation rate at normal<img src="9-2310102\7f787041-2c5c-47ae-b7dd-2301674c5d8a.jpg" />, but below normal the activity of neurons innervating these receptors increases dramatically, so much so to override the signals from central chemoreceptors in the hypothalamus, increasing <img src="9-2310102\ebabaa51-35b7-4ecc-93cb-0fe30b4ea179.jpg" /> despite a falling<img src="9-2310102\ddf21055-0bec-472d-9e01-3a66d0db96e1.jpg" />.</p></sec><sec id="s4"><title>4. Conclusion</title><p>In this work we have investigated a role of controls of cardiovascular-respiratory system during physical activity to hypoxia. The heart rate and both minute ventilation and alveolar ventilation increase; in this way the lungs transfer more oxygen and carbon dioxide and keep pace with metabolic demands. In this increase of the controls result the increase of arterial and vascular carbon dioxide. In addition, the partial arterial pressure of carbon dioxide and oxygen are maintained at normal value.</p></sec><sec id="s5"><title>REFERENCES</title></sec><sec id="s6"><title>NOTES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.29447-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">A. Cymerman and P. B. Rock, “Medical Problems in High Mountain Environments,” A Handbook for Medical Officers, USARIEM-TN94-2, US Army Research Institute of Environmental Medicine Thermal and Mountain Medicine Division Technical Report, 2009.</mixed-citation></ref><ref id="scirp.29447-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">G. Gutierrez, “A Mathematical Model of Tissue-Blood Carbon Dioxide Exchange during Hypoxia,” American Journal of Respiratory Critical Care Medecine, Vol. 169. No. 4, 2004, pp. 525-533.  
doi:10.1164/rccm.200305-702OC</mixed-citation></ref><ref id="scirp.29447-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">J. J. Tyson and B. Novak, “Regulation of the Eukariotic Cell-Cycle: Molecular Antagonism, Hysteresis, and Irreversible Transitions,” Journal of Theoretical Biology, Vol. 210, No. 2, 2001, pp. 249-263.  
doi:10.1006/jtbi.2001.2293</mixed-citation></ref><ref id="scirp.29447-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">D. Gammack, H. M. Byrne and C. E. Lewis, “Estimating the Selective Advantage of Mutant p53 Tumour Cells to Repeated Rounds of Hypoxia,” Bulletin of Mathematical Biology, Vol. 63, No. 1, 2001, pp. 135-166.  
doi:10.1006/bulm.2000.0210</mixed-citation></ref><ref id="scirp.29447-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">J. A. Dempsey, H. V. Forster and D. M. Ainsworth, “Regulation of Hyperpnea, Hyperventilation, and Respiratory Muscle Recruitment during Exercise,” In: J. A. Dempsey and A. I. Pack, Eds., Regulation of Breathing, Marcel Dekker, New York, 1995, pp. 1065-1133.</mixed-citation></ref><ref id="scirp.29447-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">S. Timischl, “A Global Model for the Cardiovascular and Respiratory System,” PhD Thesis, Karl Franzens Univer sit of Graz, 1998.</mixed-citation></ref><ref id="scirp.29447-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">S. Timischl, J. J. Batzel and F. Kappel, “Modeling the Human Cardiovascular Respiratory Control System: An Optimal Control Application to the Transition to Non REM Sleep,” Spezialforschungsbereich F-003 Technical Report 190, Karl Franzens Universitat, 2000.</mixed-citation></ref><ref id="scirp.29447-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">D. Curran-Everett, “A Classic Learning Opportunity from Fenn, Rahn, and Otis (1946): The Alveolar Gas Equation,” Advances in Physiology Education, Vol. 30, No. 2, 2006, pp. 58-62. doi:10.1152/advan.00076.2005</mixed-citation></ref><ref id="scirp.29447-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">I. Giovannini, C. Chiarla, G. Boldrini and M. Castagneto, “Calculation of Venoarterial CO2 Concentration Difference,” Journal of Applied Physiology, Vol. 74, No. 2, 1993, pp. 959-964.</mixed-citation></ref><ref id="scirp.29447-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">A. R. Douglas, N. L Jones and J. W. Reed, “Calculation of Whole Blood CO2 Content,” Journal of Applied Physiology, Vol. 65, No. 1, 1988, pp. 473-477.</mixed-citation></ref></ref-list></back></article>