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

102 MATLAB kode til simulering af forsøgsscenarie
19 par.ke = 0.138; % Insulin elimination rate [1/min]
20 par.tauD = 40.0; % CHO absorption time constant [min]
21 par.tauS = 55.0; % Insulin absorption time constant [min]
22 par.AG = 0.8; % CHO utilization [−]
23
24 SI1 = 51.2e−4; % Transport insulin sensitivity [L/mU]
25 SI2 = 8.2e−4; % Disposal insulin sensitivity [L/mU]
26 SI3 = 520e−4; % EGP insulin sensitivity [L/mU]
27
28 par.kb1 = par.ka1*SI1; % Activation rate [(L/mU)/min]
29 par.kb2 = par.ka2*SI2; % Activation rate [(L/mU)/min]
30 par.kb3 = par.ka3*SI3; % Activation rate [(L/mU)/min]
31
32 par.VI = 0.12*BodyMass; % Insulin distribution volume [L]
33 par.VG = 0.16*BodyMass; % Glucose distribution volume [L]
34 par.EGP0 = 0.0161*BodyMass; % Liver glucose production at zero insulin
35 % [mmol/min]
36 par.F01 = 0.0097*BodyMass; % Insulin independent glucose consumption
37 % [mmol/min]
38
39 par.MwG = 180.1577; % Molecular weight of glucose [g/mol]
40
41 % Parameters from Cobelli
42 norm = 1/6.00; % [mU/pmol]
43
44 par.gamma = 0.5; % Transfer rate constant between portal vein
45 % and liver [1/min]
46 par.K = (2.3*norm*par.MwG*0.1)*(BodyMass/par.VG); % [mU/mmol]
47 par.alpha = 0.05; % [1/min]
48 par.beta = 0.11*norm*0.1*par.MwG*BodyMass; % [(mU/min)/(mmol/L)]
49
50 HEb = 0.6; % [−]
51 m5 = 0.0304*(1/norm)*(1/BodyMass); % [min/mU]
52 m6 = 0.6471; % [−]
53
54 par.Sb = ((m6−HEb)/m5)*0.36; % [mU/min]
55 % For BW = 70 kg is par.Sb = ((m6−HEb)/m5)*0.36 and I(0) = 5.7677
56 % For BW = 60 kg is par.Sb = ((m6−HEb)/m5)*0.12 and I(0) = 5.7253
57 par.Gss = 4.9656; % [mmol/L] for BW = 70 kg
Listing D.2: HovorkaModel2pancreas.m
1 function xdot = HovorkaModel2pancreas(t,x,u,ug,d,par)
2 %HOVORKAMODEL The Hovorka model for the glucoregulatory system
3 %
4 % The function implements the Hovorka model for the gluco−regulatory system
5 % including absorption of food and iv administration of short−acting
6 % insulin together with Cobelli model for endogeneous insulin production.
7 %
8 % The model is in the form
9 %
10 % xdot(t) = (dx/dt)(t) = f(t,x(t),u,ug,d,par)
11 %
12 % with