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

Bilag B
MATLAB kode til fit af data
1 clear all; close all;
2
3 % Constrains for fmincon
4 A = [ −1 0 ; 0 −1 ];
5 B = [0; 0];
6
7 % Loading experimental data
8 load('eks cor.mat')
9 Data = STc99m/max(STc99m);
10
Listing B.1: fitdata.m
11 %%%%%%%%%%%%%%%%%%%%%% Cobelli %%%%%%%%%%%%%%%%%%%%%%%%
12 % Guess of initial value
13 x0 = [0.18; 0.18]; % [k21; kempt]
14 XC = fmincon(@Cfun,x0,A,B);
15
16 k21 = XC(1); % rate of grinding [1/min]
17 kempt = XC(2); % rate of gastric emptying [1/min]
18 kabs = 0.013; % rate constant of intestinal absorption [1/min]
19
20 q0 = [1; 0; 0]; % Initial condition
21 [T,Q] = ode15s(@Model1,[0 400],q0,[],k21,kempt,kabs);
22
23 qsto1 = Q(:,1); % Mass of glucose in stommach in solid phase [g]
24 qsto2 = Q(:,2); % Mass of glucose in stommach in liquid phase [g]
25 qsto = qsto1 + qsto2;
26 Qsto = qsto/max(qsto);