Matematisk Model for Mavesækkens Tømning - Danmarks Tekniske ...

90
91 % Plot of iv insulin infusion.
92 subplot(224)
93 stairs(T,IV ins,'b−','linewidth',2); grid;
94 xlabel('time [hr]','fontsize',14);
95 ylabel('Insulin infusion (mU/min)','fontsize',14);
96 xlim([T(1) T(end)]); ylim([0 uss+1]);
97 set(get(h,'CurrentAxes'),'fontsize',14)
1 clear all; close all;
2
Listing D.8: experiment simulation pan con 25.m
3 % Parameters for the Steady State.
4 BodyWeight = 70.0;
5 par = HovorkaParametersPancreas(BodyWeight);
6 uss = 0.0; % Basal rate of insulin.
7 dss = 0.0; % Basal rate of glucose.
8
9 % Compute Steady State.
10 Ipo0 = par.Sb;
11 xinit = [0 0 55.9483 23.3399 5.7677 0.0295 0.0047 0.2997 Ipo0 0]';
12 [Gss,Iss,xss]=HovorkaBasalState2Pancreas(uss,0,dss,par,xinit);
13 gss = Gss*par.MwG/10;
14
15 Ts = 10; % Time to eat the glucose [min].
16 dmeal = 40/Ts; % Amount of oral glucose [g/min].
17
18 T = 0:1:5*60; % Time of simulation [min].
19 N = length(T);
20
21 D = repmat(dss,1,N);
111
22 idd = round([60 61 62 63 64 65 66 67 68 69] + 1); % Time for meal [min].
23 D(1,idd) = D(1,idd)+dmeal; % Meal vector.
24
25 BS = 2.5; % Target value of blood glucose [mmol/L].
26 Tx = [];
27 G = [];
28 I = [];
29 X = [];
30 IV ins = repmat(0,1,N);
31 IV ins(0+1:60+1) = IV ins(0+1:60+1)+40; % Insulin infusion first
32 % hour [mU/min].
33 IV ins(61+1:300+1) = IV ins(61+1:300+1)+20; % Insulin infusion last four
34 % hours [mU/min].
35 IV glu = zeros(1,N);
36 iv glu bolus = 0; % Glucosebolus at the beginning [g/min].
37
38 % Control algorithm for iv glucoseadministration.
39 for t = 0:5:N−2
40 if t == 0
41 Xx = xss;
42 end
43 IV glu(t+1:t+5) = IV glu(t+1:t+5)+iv glu bolus;