26.07.2013 Views

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 ...

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!