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

108 MATLAB kode til simulering af forsøgsscenarie
26 Tx = [];
27 G = [];
28 I = [];
29 X = [];
30 IV ins = repmat(12,1,N); % Constant rate of insulin infusion [mU/min].
31 IV glu = zeros(1,N);
32 iv glu bolus = 0.0; % Glucosebolus at the beginning [g/min].
33
34 % Control algorithm for iv glucoseadministration.
35 for t = 0:5:N−2
36 if t == 0
37 Xx = xss;
38 end
39 IV glu(t+1:t+5) = IV glu(t+1:t+5)+iv glu bolus;
40 [tx,Gx,Ix,Xx] = HovorkaModelSimulation2Pancreas...
41 (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);
42 Xx = Xx(end,:)';
43 if Gx > BS
44 iv glu bolus = iv glu bolus − 0.25;
45 elseif Gx < BS
46 iv glu bolus = iv glu bolus + 0.25;
47 end
48 if iv glu bolus < 0
49 iv glu bolus = 0;
50 end
51 Tx = [Tx; tx(1)];
52 G = [G; Gx(1)];
53 I = [I; Ix(1)];
54 X = [X; Xx(1)];
55 end
56
57 % Plot of glucose concentration in plasma.
58 h = figure(1);
59 subplot(221)
60 plot([Tx(1) Tx(end)],[5 5],'r−',Tx,G,'b−','linewidth',2); grid;
61 xlabel('time [min]','fontsize',14); ylabel('G (mmol/L)','fontsize',14);
62 ylim([0 10]); xlim([Tx(1) Tx(end)]);
63 text(60,6,'\downarrow','fontsize',14)
64 hold on;
65 fill([Tx(1) Tx(end) Tx(end) Tx(1)],...
66 [70/par.MwG*10 70/par.MwG*10 140/par.MwG*10 140/par.MwG*10],...
67 'g','facealpha',0.2)
68 hold off
69 set(get(h,'CurrentAxes'),'fontsize',14)
70 legend('Target value of BS','Actual value of BS','Location','NorthEast')
71
72 % Plot of insulin concentration in plasma.
73 subplot(222)
74 plot(Tx,I,'b−','linewidth',2);
75 grid; xlabel('time [min]','fontsize',14); ylabel('I (mU/L)','fontsize',14);
76 xlim([Tx(1) Tx(end)])
77 set(get(h,'CurrentAxes'),'fontsize',14)
78
79 % Plot of oral glucose load and iv glucose infusion.
80 T=T/60;