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

104 MATLAB kode til simulering af forsøgsscenarie
68 VI = par.VI; % Insulin distribution volume [L]
69 VG = par.VG; % Glucose distribution volume [L]
70 EGP0 = par.EGP0; % Liver glucose production at zero insulin [mmol/min]
71 F01 = par.F01; % Insulin independent glucose consumption [mmol/min]
72 MwG = par.MwG; % Molecular weight of glucose [g/mol]
73
74 gamma = par.gamma; % [kg*mU/min*pmol]
75 K = par.K; % [mU/mmol]
76 alpha = par.alpha; % [1/min]
77 beta = par.beta; % [L*mU/min*mmol]
78 Sb = par.Sb; % [mU/min]
79 Gss = par.Gss; % [mmol/L]
80
81
82 %% Model
83 % =========================================================================
84 % CHO absorption
85 % =========================================================================
86 D = (1000/MwG)*d; % CHO input rate in glucose equivalents [mmol/min]
87 UG1 = D1/tauD; % Glucose flux from compartment 1 to 2 [mmol/min]
88 UG = D2/tauD; % Glucose flux from stomach/gut to blood [mmol/min]
89 D1dot = AG*D−UG1; % Stomach/gut compartment 1 balance [mmol/min]
90 D2dot = UG1−UG; % Stomach/gut compartment 2 balance [mmol/min]
91
92 % =========================================================================
93 % Insulin infusion
94 % =========================================================================
95 UI = u;
96
97 % =========================================================================
98 % Glucose subsystem
99 % =========================================================================
100 G = Q1/VG; % Plasma glucose concentration [mmol/L]
101
102 % Non−insulin dependent glucose consumption [mmol/min]
103 if G ≥ 4.5
104 F01c = F01;
105 else
106 F01c = F01*G/4.5;
107 end
108
109 % Renal glucose excretion [mmol/min]
110 if G ≥ 9.0
111 FR = 0.003*(G−9.0)*VG;
112 else
113 FR = 0.0;
114 end
115
116 % Mass balances for the two glucose compartments
117 Q12 = x1*Q1; % Glucose transport/distrubtion [mmol/min]
118 Q21 = k12*Q2; % Glucose transport distribution [mmol/min]
119 Q2out = x2*Q2; % Glucose disposal in adipose tissue [mmol/min]
120 EGP = EGP0*(1.0−x3); % Endogenous Glucose Production [mmol/min]
121
122 Q1dot = UG + ug − F01c − FR − Q12 + Q21 + EGP; % Compartment 1 [mmol/min]