Matematisk Model for Mavesækkens Tømning - Danmarks Tekniske ...
Matematisk Model for Mavesækkens Tømning - Danmarks Tekniske ...
Matematisk Model for Mavesækkens Tømning - Danmarks Tekniske ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
6.3 Simulering af clamp-<strong>for</strong>søg 71<br />
det aktuelle blodsukker. For at gøre controlleren bedre kunne glukoseændringen<br />
ligeledes indg˚a som et argument <strong>for</strong> at ændre glukoseinfusionen, se linje 14 og<br />
18.<br />
1 <strong>for</strong> t = 0:5:N−2<br />
2 % Initial values <strong>for</strong> compartments be<strong>for</strong>e first simulation.<br />
3 if t == 0<br />
4 Xx = xss;<br />
5 end<br />
6 % Giving same glucoseinfusion <strong>for</strong> 5 minutes.<br />
7 IV glu(t+1:t+5) = IV glu(t+1:t+5)+iv glu bolus;<br />
8 % Simulating 5 minutes.<br />
9 [tx,Gx,Ix,Xx] = Hovorka<strong>Model</strong>Simulation2Pancreas...<br />
10 (T(t+1):T(t+5),Xx,IV ins(t+1:t+5),IV glu(t+1:t+5),D(t+1:t+5),par);<br />
11 Xx = Xx(end,:)';<br />
12 % If bloodsugar is higher than target value: reduce glucoseinfusion<br />
13 % <strong>for</strong> next 5 minutes.<br />
14 if Gx > BS<br />
15 iv glu bolus = iv glu bolus − 0.25;<br />
16 % If bloodsugar is lower than target value: enlarge glucoseinfusion<br />
17 % <strong>for</strong> next 5 minutes.<br />
18 elseif Gx < BS<br />
19 iv glu bolus = iv glu bolus + 0.25;<br />
20 end<br />
21 if iv glu bolus < 0<br />
22 iv glu bolus = 0;<br />
23 end<br />
24 % Saving time−vector, plasma−glucose−vector, plasma−insulin−vector.<br />
25 Tx = [Tx; tx(1)];<br />
26 G = [G; Gx(1)];<br />
27 I = [I; Ix(1)];<br />
28 X = [X; Xx(1)];<br />
29 end<br />
6.3.1 Normoglykæmisk clamp<br />
Som det fremg˚ar af figur 6.5, holdes blodsukkeret relativt konstant p˚a 5 mmol/L.<br />
Under hele simuleringen holdes insulininfusionen konstant p˚a 12 mU/min. Dette<br />
er alt <strong>for</strong> lavt i <strong>for</strong>hold til den infusion, som blev givet under <strong>for</strong>søget, dette vil<br />
blive omtalt senere. Den intravenøse glukoseinfusion er fra starten 0 og reguleres<br />
af controlleren med en nøjagtighed p˚a 0.25 g/min. M˚altidet gives 60 minutter<br />
efter simuleringens start, og de 40 g kulhydrat indtages i løbet af 10 minutter.<br />
Ved simuleringen af normoglykæmi er den samlede mængde glukose, som er<br />
givet intravenøst, lidt over 300 g. Dette stemmer godt overens med det virkelige<br />
<strong>for</strong>søg, hvor der ca. blev brugt 1 1<br />
2<br />
pose glukoseopløsning, som hver indeholdt<br />
200 g glukose. Det fremg˚ar af figur 6.5, at glukoseinfusionsraten faldt ca. 30<br />
minutter efter m˚altidet, hvilket stemmer overens med det virkelige <strong>for</strong>søg, se