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

110 MATLAB kode til simulering af forsøgsscenarie
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.5;
45 elseif Gx < BS
46 iv glu bolus = iv glu bolus + 0.5;
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)],[10 10],'r−',Tx,G,'b−','linewidth',2); grid;
61 xlabel('time [min]','fontsize',14); ylabel('G (mmol/L)','fontsize',14)
62 ylim([0 12]); xlim([Tx(1) Tx(end)]);
63 text(60,10.5,'\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','SouthEast')
71
72 % Plot of insulin concentration in plasma.
73 subplot(222)
74 plot(Tx,I,'b−','linewidth',2); grid;
75 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;
81 subplot(223)
82 stairs(T,D,'b−','linewidth',2)
83 hold on
84 stairs(T,IV glu,'m−','linewidth',2); grid;
85 xlabel('time [hr]','fontsize',14); ylabel('Glucose (g/min)','fontsize',14);
86 xlim([T(1) T(end)]); ylim([0 32]);
87 text(2.9,11,['Total glucose = ' num2str(sum(IV glu)) 'g'],'fontsize',12)
88 set(get(h,'CurrentAxes'),'fontsize',14)
89 legend('Meal','IV glucose','Location','SouthEast')