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

82 t = 0:1:400;
83 Q = (1 + K*(t/tempt)).*exp(−(t/tempt)); % Amount of glucose in stomach [g]
84
85 % Plot of Goetze Model and data
86 h3 = figure(3);
87 plot(t,Q,'LineWidth',2)
88 hold on
89 plot(Time,Data,'.r','MarkerSize',20)
90 xlabel('Time [min]','fontsize',14);
91 ylabel('Relative amount of glucose in stomach','fontsize',14);
92 legend('Goetze','Data')
93 set(get(h3,'CurrentAxes'),'fontsize',14)
94 axis([0 400 0 1])
95 print('−depsc', '−tiff', ['goetze stomach fit'])
96 print('−dpng', '−loose', ['goetze stomach fit'])
1 function N = Cfun(x)
2
Listing B.2: Cfun.m
3 load('eks cor.mat') % Time, STc99m
4 Data = STc99m/max(STc99m);
5
6 kabs = 0.013;
7 q0 = [1; 0; 0];
8 [T,Q] = ode15s(@Model1,Time',q0,[],x(1),x(2),kabs);
9
10 qsto1 = Q(:,1); % Mass of glucose in stommach in solid phase [g]
11 qsto2 = Q(:,2); % Mass of glucose in stommach in liquid phase [g]
12 qsto = qsto1 + qsto2;
13 Qsto = qsto/max(qsto);
14
15 N = norm(Qsto − Data);
1 function N = Efun(x)
2
Listing B.3: Efun.m
3 load('eks cor.mat') % Time, STc99m
4 Data = STc99m/max(STc99m);
5 N = norm(exp(−(x(1)*Time).ˆx(2)) − Data);
1 function N = Gfun(x)
2
Listing B.4: Gfun.m
3 load('eks cor.mat') % Time, STc99m
4 Data = STc99m/max(STc99m);
5 N = norm((1 + x(1)*Time/x(2)).*exp(−Time/x(2)) − Data);
91