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

41 % Modified 23.07.09 SW
107
42 % ========================================================================
43 xdot = HovorkaModel2pancreas(0,x,u,ug,d,par);
Listing D.5: HovorkaModelSimulation2Pancreas.m
1 function [Tx,G,I,X]=HovorkaModelSimulation2Pancreas(T,x0,U,ugvec,D,par)
2 % HOVORKAMODELSIMULATION Simulation using the Hovorka model
3 %
4 % Syntax: [Tx,G,I,X]=HovorkaModelSimulation2Pancreas(T,x0,U,ugvec,D,par)
5
6 nx = length(x0);
7 N = length(T);
8 Tx(1) = T(1);
9 X = x0';
10 for k=1:N−1
11 x = X(end,:)';
12 [Tk,Xk]=ode45(@HovorkaModel2pancreas,[T(k) T(k+1)],x,[],U(:,k),...
13 ugvec(:,k),D(:,k),par);
14 X = [X; Xk];
15 Tx = [Tx; Tk];
16 end
17
18 G = X(:,3)/par.VG;
19 I = X(:,5);
Listing D.6: experiment simulation pan con 5.m
1 clear all; close all;
2
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);
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 = 5; % Target value of blood glucose [mmol/L].