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

84 MATLAB kode til kantfindingsprogram
24 C2 = load plot count(Data2,2);
25 C3 = load plot count(Data3,3);
26 C4 = load plot count(Data4,4);
27 C5 = load plot count(Data5,5);
28 C6 = load plot count(Data6,6);
29 C7 = load plot count(Data7,7);
30 C8 = load plot count(Data8,8);
31 C9 = load plot count(Data9,9);
32 C10 = load plot count(Data10,10);
33 C11 = load plot count(Data11,11);
34 C12 = load plot count(Data12,12);
35 C13 = load plot count(Data13,13);
36 C14 = load plot count(Data14,14);
37 C15 = load plot count(Data15,15);
38
39 % Creating vectors containing number of counts for each isotope
40 % together with a time vector.
41 Time = [C1(1);C2(1);C3(1);C4(1);C5(1);C6(1);C7(1);C8(1);C9(1);C10(1);...
42 C11(1);C12(1);C13(1);C14(1);C15(1)] − C1(1);
43 CTc99m = [C1(2);C2(2);C3(2);C4(2);C5(2);C6(2);C7(2);C8(2);C9(2);C10(2);...
44 C11(2);C12(2);C13(2);C14(2);C15(2)];
45 CIn111 = [C1(3);C2(3);C3(3);C4(3);C5(3);C6(3);C7(3);C8(3);C9(3);C10(3);...
46 C11(3);C12(3);C13(3);C14(3);C15(3)];
47
48 % Correction of phyical decay
49 kTc = log(2)/360; % Defining transformation constant for 99mTc.
50 DTc99m = CTc99m.*exp(kTc*Time);
51 kIn = log(2)/(2.83*24*60); % Defining transformation constant for 111In.
52 DIn111 = CIn111.*exp(kIn*Time);
53
54 % Plotting number of counts as a function of time for each isotope.
55 h=figure(16)
56 subplot(1,2,1)
57 plot(Time, CTc99m,'−b.','MarkerSize',15)
58 hold on
59 plot(Time, DTc99m,'−r.','MarkerSize',15)
60 xlabel('Time [min]','fontsize',14);ylabel('Counts','fontsize',14);
61 title('ˆ9ˆ9ˆmTc','fontsize',14)
62 axis([0 Time(end) 0 max(CTc99m)]);
63 subplot(1,2,2)
64 plot(Time, CIn111,'−b.','MarkerSize',15)
65 hold on
66 plot(Time, DIn111,'−r.','MarkerSize',15)
67 xlabel('Time [min]','fontsize',14);ylabel('Counts','fontsize',14);
68 title('ˆ1ˆ1ˆ1In','fontsize',14)
69 axis([0 Time(end) 0 max(CIn111)]);
70 set(get(h,'CurrentAxes'),'fontsize',14)
71 legend('Uncorrected data','Corrected data');
72
73 % Correction of downscatter
74 STc99m = DTc99m − 0.35*DIn111;
75 h2=figure(17)
76 plot(Time, CTc99m,'−b.','MarkerSize',15)
77 hold on
78 plot(Time, STc99m,'−g.','MarkerSize',15)