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

Bilag C
MATLAB kode til behandling
af forsøgsdata
1 clear all; close all;
2
3 % Loading data
Listing C.1: curves.m
4 load('normo time.mat') % Time, time vector.
5 load('normo decay.mat') % CTc99m, counts in stomach from Tc99m.
6 load('normo decayI.mat') % CTcIn111, counts in stomach from In111.
7
8 kTc = log(2)/360; % Defining transformation constant for 99mTc.
9 DTc99m = CTc99m.*exp(kTc*Time); % Correction of physical decay.
10
11 kIn = log(2)/(2.83*24*60); % Defining transformation constant for 111In.
12 DIn111 = CIn111.*exp(kIn*Time); % Correction of physical decay.
13
14 STc99m = DTc99m − 0.35*DIn111; % Correction of downscatter.
15 Cdata = STc99m/max(STc99m); % Calculating relative amount in stomach.
16
17 Y = (Cdata(2:end)−Cdata(1:end−1))/15; % Calculating differential
18 % quotient. Y = y(k)−y(k−1)/∆T.
19 % Plot of corrected data.
20 h = figure;
21 plot(Time,Cdata,'.','markersize',15)
22 xlabel('Time [min]','fontsize',14);
23 ylabel('Fraction of food in the stomach','fontsize',14);