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

27
90 MATLAB kode til fit af data
28 % Plot of Cobelli Model and data
29 h1 = figure(1);
30 plot(T,Qsto,'LineWidth',2)
31 hold on
32 plot(Time,Data,'.r','MarkerSize',20)
33 xlabel('Time [min]','fontsize',14);
34 ylabel('Relative amount of glucose in stomach','fontsize',14);
35 legend('Cobelli','Data')
36 set(get(h1,'CurrentAxes'),'fontsize',14)
37 axis([0 400 0 1])
38 print('−depsc', '−tiff', ['cobelli1 stomach fit'])
39 print('−dpng', '−loose', ['cobelli1 stomach fit'])
40
41 %%%%%%%%%%%%%%%%%%%%%%%% Elashoff %%%%%%%%%%%%%%%%%%%%%%%%
42
43 % Guess of initial value
44 x0 = [0.015; 1.03]; % [k; beta]
45 XE = fmincon(@Efun,x0,A,B);
46
47 D = 75; % [g]
48 kabs = 0.231; % rate constant of intestinal absorption [1/min].
49 k = XE(1); % rate of emptying [1/min].
50 beta = XE(2); % shape factor [−].
51
52 % Løsning af differentialligning
53 [T,Q] = ode15s(@Elashoff,[0 400],[0;0],[],D,kabs,k,beta);
54 qduo = Q(:,1);
55 qgut = Q(:,2);
56 qsto = D − qduo;
57 Qsto = qsto/max(qsto);
58
59 % Plot of Elashoff Model and data
60 h2 = figure(2);
61 plot(T,Qsto,'LineWidth',2)
62 hold on
63 plot(Time,Data,'.r','MarkerSize',20)
64 xlabel('Time [min]','fontsize',14);
65 ylabel('Relative amount of glucose in stomach','fontsize',14);
66 legend('Elashoff','Data')
67 set(get(h2,'CurrentAxes'),'fontsize',14)
68 axis([0 400 0 1])
69 print('−depsc', '−tiff', ['elashoff stomach fit'])
70 print('−dpng', '−loose', ['elashoff stomach fit'])
71
72 %%%%%%%%%%%%%%%%%%%%%% Goetze %%%%%%%%%%%%%%%%%%%%%%%%%%
73
74 % Guess of initial value
75 x0 = [1.14; 16.23]; % [k; tempt]
76 XG = fmincon(@Gfun,x0,A,B);
77
78 K = XG(1); % Lag phase [−]
79 tempt = XG(2); % Emptying time constant [min]
80
81 D = 75; % Amount of glucose [g]